精选数值模拟第五讲平面问题二三角形单元分析讲义

（优选）数值模拟第五讲平面问题二三角形单元分析课件第一页，共四十页。5.2三节点三角形单元的特性分析5.2.1单元作为分析对象按前面结构矩阵位移法分析思想，要求解平面问题的有限元离散结构，需要知道单元（三角形薄片）在节点自由度上受力时的弹性特性或刚度特性。这是一个新问题，一个特殊的弹性力学问题。三角形顶点设为节点，其局部编号为l，m，n（逆时针）。每节点有总体坐标x，y方向两个待求位移分量：u，v。单元共有6个位移分量——6个自由度。下面研究有限元法中特有的求解该特殊弹性力学问题的方法。有限元离散结构受力平衡后，取出一个典型三节点三角形单元e。第二页，共四十页。单元平衡时要在节点处受到节点力（节点对单元的作用力），每节点有2个节点力分量，单元有6个节点力分量。（2）单元节点力列阵下面要研究的问题是该三角形薄片弹性体在保持平衡时所受节点力和节点位移的关系。（1）单元节点位移列阵第三页，共四十页。5.2.2单元位移模式按弹性力学位移法求近似解的思路，位移作为基本未知量时，需要对单元上位移的分布作出假设，即构造含待定参量的简单位移函数——位移模式。为待定系数，称为广义坐标。通常用多项式函数作位移模式，对三节点三角形单元，有6个待定节点位移分量，所以单元上的位移函数只能是含6个待定系数的完全一次多项式：第四页，共四十页。位移多项式写成矩阵形式：代入各节点取值条件后：坐标取节点值坐标取节点值由于有限元法中未知量是节点位移，所以上面单元位移模式需要转换为以节点位移分量为待定参量的形式。过程如下：第五页，共四十页。其中：由第一组方程求解节点坐标行列式分别解出6个待定系数：第六页，共四十页。由第二组方程求解第七页，共四十页。上面求出的待定系数代回位移多项式，得到：至此，单元位移模式已转换为节点位移的插值形式。第八页，共四十页。上式中：位移插值基函数，称为形状函数（形函数）

称为形函数矩阵，是对单元节点位移进行插值得到单元位移分布函数的转换矩阵。用节点位移插值表示单元位移模式是有限元法中除了离散化之外最具代表性，最重要的步骤!!!第九页，共四十页。5.2.3形函数及其性质对于单元位移模式：假设：得到：显然，形函数决定了单元上位移分布的形态。事实上，单元位移模式就是所有形函数的线性组合。一个单元的位移模式决定了该单元描述局部位移场的能力，决定求解的精度、收敛性等，而形函数是最重要的因素。第十页，共四十页。

针对三节点三角形单元，可以导出单元形函数的下列性质。性质1：单元上某节点的形函数在该节点的值为1，在其它节点的值为零。性质2：单元上所有形函数之和等于1。（简单三角形单元的形函数只有2个独立）第十一页，共四十页。性质3（推论）：简单三角形单元的形函数在边界上的性质。某节点的形函数在该点邻边上呈线性分布，取值在0~1之间，在该点对边上值为零。简单三角形单元形函数的几何意义

由形函数表达式和性质1可画出下列形函数几何图形。第十二页，共四十页。根据位移模式表达式及其形函数的性质，可以推断出两个相邻三角形单元上位移分布形状和公共边界上位移的情况：两个单元上位移线性连续分布，各单元在公共边界上位移线性分布，数值相同——边界位移协调！第十三页，共四十页。

由图形几何性质可以推断简单三角形单元形函数的下列结论：1）三角形形心上：2）形函数在边界上的积分：3）形函数在单元上的积分：第十四页，共四十页。5.2.4单元应变和应力已知节点位移插值形式的单元位移模式：代入平面问题几何方程（应变~位移关系）得到单元应变：第十五页，共四十页。上式简写为：

为应变矩阵，其一个子块的计算式为：该式建立了用单元节点位移表达单元上应变分布的关系。对简单三角形单元，应变矩阵为：第十六页，共四十页。把单元应变代入平面问题物理方程，即得到单元应力：应力矩阵的分块形式：第十七页，共四十页。对于平面应力问题，应力矩阵的子块为：由于均质材料的弹性系数为常数，应力矩阵也是常数矩阵，即该单元上应力分布也是常数。对平面应变问题的应力矩阵，只要把上式中弹性系数作相应变换：第十八页，共四十页。平面问题简单三角形单元应力的讨论单元上应力和应变均为常数，其大小与单元几何和节点位移有关，一般各单元应力值不相同，应力在单元之间不连续，这是有限元解的近似性的反映。一般把单元上这个常应力值作为单元中心的应力值较合理，有较高的精度。理论和数值试验均可证明，用该单元求解问题时，误差随单元尺寸的减小和单元数目的增加而减小，这就是有限元解的收敛性。事实上，当单元变得越来越小时，结构在一个小单元区域上的应力趋于均匀，而大量小单元拼接在一起就可能很精确地模拟结构中变化的应力。第十九页，共四十页。5.2.5单元刚度方程与刚度矩阵前面已经得到了用单元节点位移表示的单元内部位移、应变、应力的分布：结构受力平衡时，每个单元在节点力作用下也保持平衡，节点产生相应位移。第二十页，共四十页。下面用弹性体虚位移原理（虚功原理）建立单元的平衡关系。设单元节点虚位移：则相应有单元内部虚位移：相应单元虚应变为：外力虚功计算：对一个单元而言，外力即单元节点力。按虚功原理，虚位移下单元的外力虚功等于虚应变能。第二十一页，共四十页。虚应变能计算单元虚应变能即单元内应力在虚应变上所作虚功。第二十二页，共四十页。由虚功原理得到：根据上面推导，有：外力虚功：虚应变能：考虑到虚位移的任意性，由上式立即得到下列方程：第二十三页，共四十页。上式简写成：其中单元刚度矩阵：该方程建立了单元节点力和节点位移之间的关系，为单元刚度方程。至此，我们已经得到了关于弹性实体小单元的力学特性方程。由于弹性矩阵是对称的，由上式可以看出，单元刚度矩阵是对称方阵。对3节点三角形单元，单元刚度矩阵是6×6方阵。第二十四页，共四十页。上面单元刚度方程和单元刚度矩阵的推导虽然在平面问题的简单三角形单元下得到，但推导的过程和公式的形式具有一般性，推导原理和结果适用一般的连续体力学问题的所有单元。提示：单元刚度方程和单元刚度矩阵的建立是单元分析的核心内容。一般情况下，单元应变矩阵是坐标的函数矩阵，所以单元刚度矩阵的计算需要进行积分运算。所建立的单元刚度矩阵反映了一般弹性体小单元近似的弹性性质，是单元特性的核心。第二十五页，共四十页。单元刚度矩阵的计算弹性力学平面问题的单元刚度矩阵通式：（单元刚度矩阵通式）（平面问题单刚通式）第二十六页，共四十页。由于应变矩阵是常数矩阵，所以由平面问题刚度矩阵通式得到：平面问题简单三角形单元刚度矩阵计算3×3分块形式对平面应力问题，上述刚度矩阵的一个2×2子矩阵为：第二十七页，共四十页。单元刚度方程的分块形式简单三角形单元刚度方程按节点分块形式：展开的单元刚度方程：每行各节点刚度子块反映了节点位移对节点力的贡献。第二十八页，共四十页。5.2.6简单三角形单元分析计算实例(1)

形函数和形函数矩阵前面针对一维杆单元和平面问题简单三角形单元引入的单元形函数，应变、应力矩阵，单元刚度矩阵等是结构有限元分析中最基本、最重要的概念。dd厚度h下面通过分析计算一个具体的单元加以说明。第二十九页，共四十页。按公式计算各形函数代数余子式矩阵：按形函数公式，得到三个形函数：计算节点坐标行列式：dd厚度h容易验证其满足形函数性质第三十页，共四十页。按形函数性质和几何意义，直接写出形函数表达式dd厚度h形函数的图形是一个平面，其方程是：上述两个函数定义在三角形区域的部分就是形函数。同理，直接写出形函数图形的方程为：由形函数的性质，得到第三个形函数：显然，与公式计算得到的结果相同！第三十一页，共四十页。3）形函数矩阵根据计算所得单元形函数，写出单元的形函数矩阵：第三十二页，共四十页。(2)应变矩阵由公式计算单元应变矩阵：第三十三页，共四十页。(3)应力矩阵对平面应力问题：第三十四页，共四十页。(4)单元刚度矩阵计算平面应力问题第三十五页，共四十页。注意：该刚度矩阵每列元素之和等于零矩阵的行列式=0奇异该刚度矩阵元素与单元边长d无关形状相似的单元刚度矩阵相同？平面问题的单元刚度矩阵与单元厚度成正比。第三十六页，共四十页。(5)单元内部力学状态分析单元的形函数、应变矩阵、应力矩阵及刚度矩阵都表征单元的某些特性。对应任何一种单元节点位移状态，通过它们都能给出单元上相应的力学状态。设单元节点位移为：dd厚度h1由单元位移模式：得单元位移：单元位移第三十七页，共四十页。由单元应变方程：得单元应变：dd厚度h1X方向均匀拉伸单元应变第三十八页，共四十页。单