摘要：一、对比 对于受均布荷载的简支梁，假设梁的长度L=1，均布荷载大小为1，弹模E为1，惯性矩I为1，那么梁的挠曲线方程的解析解为 $$v(x)=-\left(\frac{x^4}{24}-\frac{x^3}{12}+\frac{x}{24}\right)$$ 梁单元的形函数如下 根据形函数的定义，只要知道了单元两端节点的位移，即可用形函数求出单元上任意位置的挠度。即 $$v(x)=N(x)\c... 阅读全文

摘要：不得不说，Mathematica真是个好东西，以前学习有限元的时候，对于书中的方程推导，看到了就看过去了，从没有想过要自己推导一遍，原因是手工推导太复杂。有了MM，原来很复杂的东西突然变得简单了。1.单元几何描述上图是纯弯梁单元，长度l，弹模E，面积A，惯性矩I。两个节点1和2的位移列阵为\[ q^... 阅读全文

摘要：曾攀老师的《有限元分析基础教程》第三章有二维杆单元的推导，并结合一个例题进行了解析解和基于Matlab的程序求解。但是我感觉书中的MATLAB代码有点罗嗦，而且一些实现方法也比较麻烦，比如已经知道了杆单元的起点和终点坐标，仍然需要另外给出单元局部坐标与整体坐标的夹角，这完全没必要。于是我就用Pyth... 阅读全文