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

第二章 有限元分析基本理论 有限元法的基本思路是将一个连续求解区域分割成有限个不重叠且按一定方式相互连接在一起的子域 单元 ，利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数。单元内的场函数通常由未知场函数或其导数在单元各个节点的数值和其插值函数来近似表示。这样，未知场函数或其导数在各个节点上的数值即成为未知量 自由度 。根据单元在边界处相互之间的连续性，将各单元的关系式集合成方程组，求出这些未知量，并通过插值函数计算出各个单元内场函数的近似值，从而得到全求解域上的近似解。 有限元将一个连续的无限自由度问题变成离散的有限自由度问题进行求解。如果将区域划分成很细的网格，也即单元的尺寸变得越来越小，或随着单元自由度的增加及插值函数精度的提高，解的近似程度将不断被改进。如果单元是满足收敛要求的，近似解最后可收敛于精确解。 2.1 有限元分析的基本概念和计算步骤 首先以求解连续梁为例，引出结构有限元分析的一些基本概念和计算步骤。 如图2-1，连续梁承受集中力矩作用。将结构离散为三个节点，两个单元。结构中的节点编号为1、2、3；单元编号为①、②。 图2-1 受集中力矩作用的连续梁 2.1.1单元分析 在有限元分析过程中，第一步是进行结构离散，并对离散单元进行分析，分析的目的是得到单元节点的力与位移的关系。单元分析的方法有直接法和能量法，本节采用直接法。 从连续梁中取出一个典型单元e，左边为节点i，右边为节点j。将节点选择在支承点处，单元两端只产生转角位移、，顺时针转动为正。独立的单元杆端内力为弯矩、，顺时针为正。 记：为单元e的节点位移向量；为单元e的杆端力向量。 根据结构力学位移法可得如下平衡方程： 2-1 式中：，，、分别为单元的抗弯刚度和长度。 i，j 1，2 的物理意义为单元处发生单位转角引起的处的力矩，将式 2-1 写成矩阵形式 2-2 或 2-3 式 2-2 、 2-3 称为梁单元的刚度方程。式中，称为梁单元的刚度矩阵，只要已知梁单元的、就可计算出单元刚度矩阵。 以上分析实现了单元分析的目的，即得到单元刚度方程和单元刚度矩阵。 2.1.2整体分析 有限元分析的第二步要将离散的单元集成整体，组集过程可见图2-2。在组集过程中，必须满足以下条件： 图2-2 离散的单元集成整体 1 变形协调 2-4 2 节点平衡 2-5 式 2-2 代入 2-5 可得： 2-6 式 2-4 代入 2-6 整理可得： 2-7 写成矩阵形式，得 2-8 式 2-8 称为结构刚度方程，它实际上是结构的节点平衡方程，记为 2-9 式中：称为该结构的原始刚度矩阵；称为该结构的位移向量；称为该结构的节点荷载向量。 以上分析实现了整体分析，即得到结构原始刚度矩阵和结构刚度方程。 2.1.3用直接刚度法形成结构刚度矩阵 通过整体分析，建立了节点的平衡方程，即结构的刚度方程，从而得到结构刚度矩阵。但是，要实现电算，不可能对每一具体结构都作一次总体分析，而应该找一种规律，在确定了节点位移和荷载的排序后，使计算机能够直接由单元刚度矩阵集成结构刚度矩阵，从单元刚度方程得到结构的刚度方程，这一方法称为直接刚度法。下面介绍用直接刚度法直接由单元刚度矩阵集成结构刚度矩阵的过程。 确定结构刚度矩阵的阶数。 结构刚度方程中第i行，表示该结构第i个位移分量上力的平衡方程，因此，如果结构有N个独立位移分量，就可列出N个独立平衡方程，结构刚度矩阵就是N×N阶的。本例有3个独立的位移分量，故总刚必然为3×3阶的，写成： 2-10 确定单元刚度矩阵中元素与结构刚度矩阵中元素的关系 若将单元刚度矩阵下标写成位移分量编号的形式。 单元1：， 2-11 单元2：， 2-12 有：，，，，，，，，。 可见，若将单元刚度矩阵中元素下标写成位移分量编号的形式，则结构刚度矩阵中任一刚度元素与单元刚度矩阵中元素有如下关系： 2-13 式中：—单元号，—结构单元总数 因此，用直接刚度法集成总刚，可归纳为以下几步： 结构未知量进行编号，确定各未知量在结构刚度方程中的位置 行号 ； 确定结构刚度矩阵的阶数N； 对单元e进行循环，寻找e单元刚度矩阵中各元素下标对应于整体刚度方程中的未知量编号；并按此编号，根据式 2-13 分别叠加到结构总体刚度矩阵中的对应位置上去。 对单元循环完毕时，结构刚度矩阵就形成了。形成结构刚度矩阵是有限元分析过程中十分重要的环节，为了节约计算机存储空间，加快刚度方程求解速度，我们还必须了解结构刚度矩阵具有如下性质： 结构刚度矩阵是N×N阶的方阵，N为结构的未知量总数。 结构刚度矩阵是对称阵，即，这一性质由力—位移互等定理决定。 处于同一单元上的两个未知量称相关未知量。若两个未知量不相关，则。由式 2-13 可知，两个未知量不相