材料成形计算机模拟 第二章 线性弹性力学问题的有限元法1.ppt

第二章 线性弹性力学问题的有限元法 本章内容提要 本章介绍有限元的基础知识,仅考虑线弹性力学问题。 针对弹性力学平面问题，通过最简单的三角形常应变单元,阐明有限元用于弹性体应力分析的基本原理和方法,然后简要地介绍轴对称问题和三维问题的有限元法。 简要介绍等参单元以及数值积分方法. 本章最后讨论了有限元的一般化，说明有限元是一种求解微分方程边值问题和初值问题的普遍适用的离散化方法。 §2-1 弹性力学的基本方程 一、线弹性静力分析问题的假设 二、基本方程 三、利用虚功原理建立有限元方程 一、线弹性静力分析问题的假设 1 位移梯度是小量，因此应变与位移之间的关系是线性的； 2 物体始终保持为弹性状态，应力与应变之间的关系是线性的； 3 边界条件中不包含接触边界条件。 三、利用虚功原理建立有限元方程 1、虚功的概念 变形体在产生虚位移过程中真实力(外力，应力)所作的功，称为虚功。 虚位移是一种加在变形体上假想的微小的位移。虚位移必须满足变形体的约束条件并保持其连续性。而虚位移的类型(移动，转动)和方向是任意的。所谓微小的就是在产生虚位移过程中，各力的作用线保持不变。 2、虚功原理 如果变形体处于平衡状态，对其给一任意微小虚位移，则外力所作的总虚功必须等于内力(应力合力)所作的总虚功，这就是变形体的虚功原理。 3、虚功原理的应用 在弹性力学中，利用虚功原理可以对某些比较复杂的问题求得近似的解析解，而在有限元法中可对某些更复杂的问题求得近似的数值解。通常用虚功方程代替平衡方程与几何方程。虚功原理是力学中的一个普通原理，它不仅可用于线弹性问题，而且可以用于非线性弹性及弹塑性等材料的非线性问题。 4、有限元方程的建立 利用虚功原理来建立有限元方程。设物体处于平衡状态，给物体内各点以任意的虚位移，它仅是坐标x的函数，同时在Su上满足。根据虚功原理，外力在虚位移上所作的虚功等于因虚位移引起的虚应变能，即 §2-2 弹性力学平面问题的有限元列式 一、平面问题3节点三角形单元的插值函数 二、单元刚度方程 三、单元等效节点载荷列阵 四、整体刚度矩阵和总节点载荷列阵的集成 五、引入边界条件 六、实施步骤与注意事项 一、平面问题3节点三角形单元的插值函数 在有限单元法中，选择节点位移作为基本未知量时称为位移法； 选择节点力作为基本未知量时称为力法； 取一部分节点力和一部分节点位移作为基本未知量时称为混合法。 位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。 当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如位移，应变和应力等由节点位移来表示。 这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。 通常，有限元法我们就将位移表示为坐标变量的简单函数，这种函数称为位移模式或位移函数。 （一）单元的位移模式及插值函数 在一个单元中，每个节点有两个位移分量，节点i的位移用列向量记为 每个单元共有3个节点，6个位移自由度，可表示如下单元节点位移列阵 2、形函数的性质 （1）在节点上形函数的值满足 2、形函数的性质 （2）在单元中任一点有 3、面积坐标与直角坐标的转换关系 二、单元刚度方程 采用虚功原理建立有限元刚度方程。 三、单元等效节点载荷列阵 在有限元分析中，要针对节点自由度建立求解方程，因此单元所受到的面积力和体积力都要由作用于其节点的等效载荷来表示。 几种常见载荷的计算方法： 1、均匀等厚单元的自重 2、均布测压q 四、整体刚度矩阵和总节点载荷列阵的集成 （一）集成方法 假设弹性体被离散成ne个三角形单元和n个节点，则整个弹性体的总虚功应是ne个单元虚功的总和。将扩阶后的全部ne个单元的单元刚度方程叠加起来，则弹性体的整体刚度方程可写成KU=P。 K是所有单元的刚度矩阵之和，称为整体刚度矩阵或总刚度矩阵。P是由作用于弹性体上的载荷（包括体力、面力和集中力）分别移置到有关单元的节点上，并逐个节点加以合成得到的总等效节点力，按照节点编号由小到大排列组成的，所以它称为总节力列阵或载荷列阵。 在实际编程计算的过程中不采用上述先将单元刚度矩阵和单元等效节点载荷列阵扩阶，然后再相加的方法，在计算得到Ke和Pe的各元素后，只需按照单元中各节点的整体自由度编码，“对号入座”地叠加到整体刚度矩阵和整体节点载荷列阵地相应位置上即可。这个过程又被称为组装或集成。 （二）整体刚度矩阵的特点 1）与单元刚度矩阵类似，整体刚度矩阵K中的一列元素kij的物理意义是要使弹性体的节点自由度i发生单位位移，而其他节点自由度都保持零位移的状态下，需要施加的与节点自由度j对应的节点力。 2）整体刚度矩阵K的主元素总是正的。 3）整体刚度矩阵K是对称矩阵。 4）整体刚度矩阵K是一个稀疏矩阵，如果遵守一定的节点编号规则，可使非零元素都