数值模拟：第五讲 平面问题(二)——三角形单元分析.ppt

* 5.1 离散化要解决的问题 如何进行离散结构的求解？ 如何得到小单元的刚度特性？ 为什么离散化结构的解能作为原连续问题的近似解？ 研究上述问题就涉及到有限元法的基本原理和基本理论。 5 三节点三角形单元解平面问题 目 标: 掌握平面问题 简单三角形单元位移模式。 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 单元应变和应力 已知节点位移插值形式的单元位移模式： 代入平面问题几何方程（应变~位移关系）得到单元应变： 上式简写为： 为应变矩阵，其一个子块的计算式为： 该式建立了用单元节点位移表达单元上应变分布的关系。 对简单三角形单元，应变矩阵为： 把单元应变代入平面问题物理方程，即得到单元应力： 应力矩阵的分块形式： 对于平面应力问题，应力矩阵的子块为： 由于均质材料的弹性系数为常数，应力矩阵也是常数矩阵，即该单元上应力分布也是常数。 对平面应变问题的应力矩阵，只要把上式中弹性系数作相应变换： 平面问题简单三角形单元应力的讨论 单元上应力和应变均为常数，其大小与单元几何和节点位移有关，一般各单元应力值不相同，应力在单元之间不连续，这是有限元解的近似性的反映。一般把单元上这个常应力值作为单元中心的应力值较合理，有较高的精度。 理论和数值试验均可证明，用该单元求解问题时，误差随单元尺寸的减小和单元数目的增加而减小，这就是有限元解的收敛性。事实上，当单元变得越来越小时，结构在一个小单元区域上的应力趋于均