第二章有限元分析基本理论

6 第二章 有限元分析基本理论 有限元法的基本思路是将一个连续求解区域分割成有限个不重叠且按一定方式相互连接在一起的子域(单元)，利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数。单元内的场函数通常由未知场函数或其导数在单元各个节点的数值和其插值函数来近似表示。这样，未知场函数或其导数在各个节点上的数值即成为未知量(自由度)。根据单元在边界处相互之间的连续性，将各单元的关系式集合成方程组，求出这些未知量，并通过插值函数计算出各个单元内场函数的近似值，从而得到全求解域上的近似解。 有限元将一个连续的无限自由度问题变成离散的有限自由度问题进行求解。如果将区域划分成很细的网格，也即单元的尺寸变得越来越小，或随着单元自由度的增加及插值函数精度的提高，解的近似程度将不断被改进。如果单元是满足收敛要求的，近似解最后可收敛于精确解。 2.1 有限元分析的基本概念和计算步骤 首先以求解连续梁为例，引出结构有限元分析的一些基本概念和计算步骤。 如图2-1，连续梁承受集中力矩作用。将结构离散为三个节点，两个单元。结构中的节点编号为 1、2、3；单元编号为①、②。 图2-1 受集中力矩作用的连续梁 2.1.1单元分析 在有限元分析过程中，第一步是进行结构离散，并对离散单元进行分析，分析的目的是得到单元节点的力与位移的关系。单元分析的方法有直接法和能量法，本节采用直接法。 从连续梁中取出一个典型单元 e，左边为节点 i，右边为节点 j。将节点选择在支承点处，单元两端只产生转角位移ei 、ej ，顺时针转动为正。独立的单元杆端内力为弯矩im 、jm ，顺时针为正。 记：ejieu 为单元 e 的节点位移向量； ejiemmf为单元 e 的杆端力向量。 根据结构力学位移法可得如下平衡方程： ejeeieejejeeieeikkmkkm22211211 (2-1) 式中：eeeeeeikkikk2412212211，lEIie ， EI 、l 分别为单元e 的抗弯刚度和长度。 7 eijk (i，j=1，2)的物理意义为单元j 处发生单位转角引起的i 处的力矩，将式(2-1)写成矩阵形式 ejieejikkkkmm22211211 (2-2) 或     eeeukf (2-3) 式(2-2)、(2-3)称为梁单元e 的刚度方程。式中， ek称为梁单元e 的刚度矩阵，只要已知梁单元的EI 、l 就可计算出单元刚度矩...