有限元方法 章 平面三角形单元.pptVIP

第五章 平面三角形单元; §5–1 有限元法的基本思想 §5–2 三角形常应变单元 §5–3 形函数的性质 §5–4 刚度矩阵 §5–5 等效节点力载荷列阵 §5–6 有限元分析的实施步骤 §5–7 计算实例;一、有限元法的基本思想;二、经典解与有限元解的区别：;x; 根据题目的要求，可选择适当的单元把结构离散化。对于平面问题可用三角元，四边元等。 ; 　　结构离散化后，要用单元内结点的位移通过插值来获得单元内各点的位移。在有限元法中，通常都是假定单元的位移模式是多项式，一般来说，单元位移多项式的项数应与单元的自由度数相等。它的阶数至少包含常数项和一次项。至于高次项要选取多少项，则应视单元的类型而定。;4. 单元的力学特性分析; 　　最后利用弹性体的虚功方程建立单元结点力阵与结点位移 列阵之间的关系，即形成单元的刚度方程式： 　　; 　　考虑整体结构的约束情况，修改整体刚度方程之后，（5-6）式就变成以结点位移为未知数的代数方程组。解此方程组可求出结点位移。; 求解出整体结构的位移和应力后，可有选择地整理输出某些关键点的位移值和应力值，特别要输出结构的 变形图、应力图、应变图、结构仿真变形过程动画图及整体结构的弯矩、剪力图等等。;一、离散化;? ? 图3-1 弹性体和有限元计算模型;? ? 图3--2 平面三角形单元;二、位 移; 从弹性力学平面问题的解析解法中可知，如果弹性体内的位移分量函数已知，则应变分量和应力分量也就确定了。但是，如果只知道弹性体中某几个点的位移分量的值，那么就不能直接求得应变分量和应力分量。因此，在进行有限元分析时，必须先假定一个位移模式。由于在弹性体内，各点的位移变化情况非常复杂，很难在整个弹性体内选取一个恰当的位移函数来表示位移的复杂变化，但是如果将整个区域分割成许多小单元，那么在每个单元的局部范围内就可以采用比较简单的函数来近似地表示单元的真实位移，将各单元的位移式连接;起来，便可近似地表示整个区域的真实位移函数。这种化繁为简、联合局部逼近整体的思想，正是有限单元法的绝妙之处。;(c);? 图5-2 平面三角形单元;其中; 式中 I是二阶单位矩阵；Ni 、Nj 、Nm 是坐标的函数，它们反映了单元的位移状态，所以一般称之为形状函数，简称形函数。矩阵 [N] 叫做形函数矩阵。三节点三角形单元的形函数是坐标的线性函数。单元中任一条直线发生位移后仍为一条直线，即只要两单元在公共节点处保持位移相等。则公共边线变形后仍为密合。 ;三、应 变;可简写成;四、应 力;其中 [S]叫做应力矩阵，若写成分块形式，有;;;在上节中，提出了形函数的概念，即;;;;;;;;;;一. 单元刚度矩阵;;;;; 继而再引入整个弹性体的载荷列阵{R}2n×1 ，它是移置到节点上的等效节点载荷依节点号码从小到大的顺序排列组成，即;;; 图中有两种编码：一是节点总码：1、2、3、4；二是节点局部码，是每个单元的三个节点按逆时针方向的顺序各自编码为1，2，3。;整体刚度矩阵分块形式为：; 总码 1 2 3 4 1 2 3 4 3 1 2 局部码 ; 总码 1 2 3 4 1 2 3 （5-42） 3 1 2