3章 平面问题的有限单元法.pptVIP

3.1 3结点三角形单元 3.2 3节点三角形单元内任意一点的应变 3.3 3节点三角形单元内任意一点的应力 3.4 3节点三角形单元全部节点位移的求解 3.5 3节点三角形单元刚度矩阵的求解 3.7 结构刚度矩阵的快速建立法 3.7.1 集合的基本原则 3.9 位移函数与解的收敛性 3.10 3节点三角形单元求解公式总结 3.11 3结点常应变单元应用实例 3.12 6结点三角形单元 3.12.1 面积坐标 3.12.2面积坐标的特点 3.12.2面积坐标的微积分计算 3.12.4面积坐标与直角坐标的关系 3.12.5 位移函数 3.12.6 用面积坐标给出的三角形单元的形函数 3.12.7 几何矩阵 3.12.8 应力矩阵 3.12.9 单元刚度矩阵 3.12.10 等效结点载荷 1.集中力的等效结点载荷 若单元内任意一点受有集中力P＝[Px,Py]T，根据载荷移置的普遍公式有： Re=NTP 【接下页】 [补充]：当三角形单元采用面积坐标作为自然坐标时，被积函数可以方便地在三角形单元全域或沿边界积分。如3结点三角形单元载荷移置问题。 1.自重产生的等效结点载荷： 2. 3结点三角形单元ij边作用有均布测压q（如图所示） i j m x y ? 这种计算的简便，在高次单元中尤为明显。 下图为一等厚度的矩形薄板，其一端固定，另一端承受载荷集度为q(kg/cm2)的均布拉力。板长l＝200cm，宽h=100cm，厚度为t。材料的弹性模量为E，泊松比u=1/3。求薄板端角点的位移及板的应力。 y 200cm 100cm q(kg/cm2) x R2x=qht/2 R3x=qht/2 ? ? 1(0,0) 2(200,0) 4(0,100) 3(200,100) (a)结构模型图 (b)有限元模型图 图3－2 等厚矩形薄板 1.结构离散化 为简单起见，将该矩形薄板划分为两个单元、4个结点的单元组合体，单元与结点的编号如图（b)所示。作用于单元?的均布拉力，按能量等效原则移置到结点2与结点3后，有： R2x=qht/2=50qt，R2y=0，R3x=qht/2=50qt，R3y=0 将一端固定的约束条件简化为固定端点的结点1、4处的平面铰。直角坐标系的选取如图(b)所示。 200cm 100cm q(kg/cm2) x R2x=qht/2 R3x=qht/2 ? ? 1(0,0) 2(200,0) 4(0,100) 3(200,100) (a)结构模型图 (b)有限元模型图 2.计算单元刚度矩阵 4.建立结构的结点载荷列阵 因结点载荷列阵R中不必考虑约束反力的作用，故在R1x、R1y、R4x、 R4y处置零，于是得到 3.建立结构刚度矩阵 采用按单元或结点形成结构刚度矩阵的方法，由K?及K?直接形成总体刚度矩阵K： R=[0， 0， 50qt， 0，50qt， 0，0 ，0]T 5.引入位移边界条件，求解线性代数方程组 在结点1和4处，有u1=v1=u4=v4=0。这里采用降阶法引入上述位移边界条件。 [注：事实上，可以通过总体刚度方程求出约束反力，当不考虑约束反力时才可用降阶法] 6.计算单元应力 7.计算结点处的应力 采用绕结点平均法 ?i ?j ?m m(0,0,1) j(0,1,0) i(1,0,0) d di Li=0 Lj=0 Lm=0 P 图5－4 三角形单元的面积坐标 （1)三角形内与结点i的对边j-m平行的直线上的诸点有相同的Li坐标，这个坐标值等于该直线至jm边的距离di与结点i至jm边的距离d的比值。 （2）三个面积坐标并不相互独立，三个面积坐标必然满足： Li+Lj+Lm=1 [注：因为?i+ ?j+ ?m= ?] 从上式可见3个面积坐标中只有2个是独立的。 由于三角形的面积坐标与该三角形的具体形状及其在总体坐标x,y中的位置无关，因此它是三角形的一种自然坐标。 面积坐标表示的函数对直角坐标求导时，采用一般的复合函数求导的法则： 6结点三角形单元的位移函数，可以采用完全的二次多项式，即单元中任意一点的位移分量为： 用面积坐标作为三角形单元的自然坐标时，可以直接得到用面积坐标表示的插值函数而不必经过先求解广义坐标，然后再由广义坐标得到单元插值函数的琐长计算过程。并且在得到插值函数以后，可以利用面积坐标的微积分计算公式方便地计算单元的刚度矩阵和 等效结点载荷列阵。 在讨论二次三角形单元用面积坐标给出的插值函数之前，先讨论线性三角形单元的面积坐标插值函数形式。 i j m 1 2 3 Lm=0 L