第三讲 有限单元法.doc

第3章 有限单元法 在工程技术领域内，工程师常常运用数学和力学的知识将实际问题抽象成它们应遵循的基本方程（常微分方程或偏微分方程）和相应的边界条件。对于大多数的工程技术问题，由于物体的几何形状和载荷作用方式是很复杂的，除了方程性质比较简单且几何边界相当规则的少数问题之外，试图按经典的弹性力学和塑性力学方法获得解析解是十分困难的，甚至是不可能的。为了克服这种困难，有两条解决途径：一是引入简化假设，将方程和边界条件简化为能够处理的问题，从而得到它在简化状态下的解答。这种方法只在有限的情况下可行，因为过多的简化将可能导致不正确的甚至错误的答案。另一条解决途径就是数值解法，如有限差分法、边界元法、有限单元法和离散元法等。对于非线性问题，有限单元法更为有效，且已经出现了许多通用程序。 有限单元法的主要优点是：①建立于严格理论基础上的可靠性。因为用于建立有限元方程的变分原理或加权余量法在数学上已被证明是微分方程和边界条件的等效积分形式。只要原问题的数学模型是正确的，同时用来求解有限元方程的算法是稳定、可靠的，如果单元满足收敛准则，则近似解最后收敛于原数学模型的精确解；②适应性强，应用范围广，不仅能成功地分析具有复杂边界条件、非线性、非均质材料、动力学等难题，而且还可以推广到解答数学方程中的其它边值问题，如热传导、电磁场、流体力学等问题；③适合计算机实现的高效性。由于有限元分析的各个步骤可以表达成规范化的矩阵形式，最后导致求解方程可以统一为标准的矩阵代数问题，特别适合计算机的编程和执行。已经出现了许多大型结构分析通用程序，如：NASTRAN、ASKA、ADINA、ANSYS、ABAQUS等，可以直接应用。这些优点使有限单元法得到了广泛的应用和发展。 3.1有限单元法分析的基本步骤 在工程或物理问题的数学模型（基本变量、基本方程、求解域和边界条件等）确定以后，有限单元法作为对其进行分析的数值计算方法的基本步骤如下： (1) 离散化 一个复杂的弹性体可以看成是由无限个质点组成的连续体，它具有无限个自由度。将一个受外力作用的连续弹性体离散成一定数量的有限个小单元的集合体，单元之间只在节点上互相联系，亦即只有节点才能传递力。因此，该集合体只具有有限个自由度，这就为解算提供了可能。由无限个质点的连续体转化为有限个单元的集合体的过程，称为离散化。在数学意义上说，就是把微分方程的连续形式转化为代数方程组，以便于进行数值求解。 将求解域离散为有限单元，根据基本场变量与坐标的关系而决定采用一维、二维和三维单元。一维单元用线段表示，二维单元可为三角形元或四边形元，三维单元常用四面体或六面体元。单元划分越密，计算精度越高，但计算工作量也越大。 (2) 单元分析 根据弹性力学的基本方程和变分原理建立单元节点力和节点位移之间的关系，形成单元有限元方程。 1）确定插值函数（形函数） 有限单元法将整个求解域离散为一系列仅靠公共节点联接的单元，而每一单元本身却视为光滑的连续体。单元内任一点的场变量（如位移），可根据其在单元中的假定分布规律（插值函数）Ne 由本单元的节点值插值求得。 N称为单元的形函数矩阵，它与单元节点坐标、节点数目（即单元形状）及插值形式有关。形函数矩阵分量的数目应与单元节点自由度数相等。 所谓“形函数”是一种视单元节点量为“已知量”的插值函数，由单元节点的坐标可直接求出形函数的具体形式。 2）建立单元方程 当问题比较简单时，可以直接根据问题的物理概念建立单元方程。不过，在一般情况下，特别是二维和三维单元，直接法会显得过于繁杂而难以应用。为此，需要采用更为一般的数学方法，如变分法、加权余量法或虚功原理。 (3) 计入边界条件，求解有限元方程 组集后的总体特征矩阵（或称为总刚度矩阵）是奇异的，必须计入边界条件才能求得唯一解。计入边界条件的方法有三种： 1）直接代入法 将自由度的已知量（边界条件）从总体方程组中消去，得到一组阶数降低了的修正方程，其原理是按节点位移已知和待定重新组合方程。由于这种方法改变了方程组的阶数，使程序编制复杂化，故程序中一般不采用。 2）对角线元素置1法 当给定位移值是零位移时，例如无移动的铰支座，可以将系数矩阵K中与零节点位移相对应行列的主对角元素改为1，其它元素改为0；再将载荷列阵中与零节点位移相对应的元素改为0即可。 这种计入边界条件的方法简单，不改变原方程组的阶数和未知量顺序，但只适用于边界条件为零值的情况。 3）对角元素乘大数法 当有节点位移为给定值时，第j个方程作如下修改：对角元素Kjj乘以大数（可取1010左右量级），并将Pj用取代，即： （3-1） 经过修改后的第j个方程为： （3-2） 对于多个给定位移则按序将每个给定位移都作上述修正，得到全部进行修正后的K和P，然后解方程则