弹性力学第6章(1-2-3节).pptVIP

第六章用有限单元法解平面问题;简史;20世纪50年代，平面问题的FEM建立，

并应用于工程问题。

20世纪60年代后，FEM应用于各种力学

问题和非线性问题，并得到迅速发展。

1970年后，FEM被引入我国，并很快地得

到应用和发展。;导出方法;§6-1基本量和基本方程的矩阵表示;面力

位移函数

应变

应力

结点位移列阵

结点力列阵

;物理方程

其中D为弹性矩阵，对于平面应力问题是;——结点虚位移，

——对应的虚应变。

在FEM中，用结点的平衡方程代替平衡微分

方程。;以下来导出FEM。

1.结构离散化——将连续体变换为离散

化结构；;结力研究的对象是离散化结构。如桁架，各单元（杆件）之间除结点铰结外，没有其他联系（图（a））。;将连续体变换为离散化结构（图（c））：即将连续体划分为有限多个、有限大小的单元，并使这些单元仅在一些结点处用绞连结起来，构成所谓‘离散化结构’。;图（c）与图（a）相比，两者都是离散化结构；区别是，桁架的单元是杆件，而图（c）的单元是三角形块体（注意：三角形单元内部仍是连续体）。;2.应用结构力学方法(位移法)进行求解:

;（2)应用插值公式,由单元结点位移，求单元的位移函数;（5）应用虚功方程，由单元的应力，求出

单元的结点力，表示为;——结点对单元的作用力，作用

于单元，称为结点力，以正标向为正。

——单元对结点

的作用力，与数

值相同,方向相反，

作用于结点。;（6）将每一单元中的各种外荷载，按虚功

等效原则移置到结点上，化为结点荷

载，表示为;各单位移置到i结点上的结点荷载

其中表示对围绕i结点的单元求和；;结力法求解;思考题

1.桁架的单元为杆件，而平面体的单元为三角形块体，在三角形内仍是作为连续体来分析的。试考虑后者在用结构力学方法求解时，将会遇到什么困难？

2.在平面问题中，是否也可以考虑其它的单元形状，如四边形单元？;FEM是取结点位移为基本未知数的。但其中每一个单元仍是连续体，所以按弹力公式求应变、应力时，必须首先解决：如何由单元的结点位移来求出单元的位移函数应用插值公式，可由求出位移d。这个插值公式表示了单元中位移的分布形式，因此称为位移模式。;

泰勒级数展开式中，低次幂项是最重要的。∴三角形单元的位移模式，可取为

位移模式在结点应等于结点位移值由此可求出;求解前三个方程，即可得：;其中;其中包含将式按未知数归纳，可表示为

或用矩阵表示为

;N—称为形（态）函数矩阵。;其中;三结点三角形单元的位移模式，略去了二次以上的项，因而其误差量级是且其中只包含了的一次项，所以在单元中的分布如图（a）所示，的分布如图所示。;FEM中以后的一系列工作，都是以位移模式为基础的。

所以当单元趋于很小时，即时，为了使FEM之解逼近于真解，即为了

保证FEM收敛性,位移模式应满足下列

条件：;（1）位移模式必须能反映单元的刚体位移。

（2）位移模式必须能反映单元的常量应变。

因为当单元时，单元中的位移和应

变都趋近于基本量——刚体位移和常量

位移。

;收敛性条件;（3）位移模式应尽可能反映位移的连续性。

即应尽可能反映原连续体的位移连续

性。在三角形单元内部，位移为连续；

在两单元边界ij上，之间均为线

性变化，也