有限元2-弹性力学平面问题有限单元法(21三角形单元,22几个问题讨论).docxVIP

第2章 弹性力学平面问题有限单元法 三角形单元(triangular Element) 三角形单元是有限元分析中的常见单元形式之一，它的优点是: ①对边界形状的适应性较好，②单刚形式及其推导比较简单,故首先介绍之。一、结点位移和结点力列阵 设右图为从某一结构中取出的一典型三角形单元。 在平面应力问题中，单元的每个结点上有沿x、y两个方向的力和位移，单元的结点位移列阵规定为： 相应结点力列阵为: (式2-1-1) u i v i d u i d e j d v j X i Y i F X i F e j F Y j  -1-1) (2 j d j F m X m m m Y m m 二、单元位移函数和形状函数 前已述及，有限单元法是一种近似方法，在单元分析中，首先要求假定（构造）一组在单元内有定义的位移函数作为近似计算的基础。即以结点位移为已知量，假定一个能表示单元内部（包括边界）任意点位移变化规律的函数。 构造位移函数的方法是：以结点(i,j,m)为定点。以位移(ui ,vi ,…um vm ) 为定点上的函数值，利用普通的函数插值法构造出一个单元位移函数。 在平面应力问题中，有u,v两个方向的位移,若假定单元位移函数是线性的, 则可表示成: u u(x,y) 1 2x 3y v v(x,y) 4 5x 6 y (2-1-2)a 式中的6个待定常数α1 ,…, α6 可由已知的6个结点位移分量(3个结点的坐标) 确定。将3个结点坐标(xi,yi ),(xj,yj ),(xm,ym )代入上式得如下两组线性方程: u x y ui 1 2 xi 3 yi j 1 2 j u x 3 j (a) y m 1 2 m 3 m 和 v x y vi 4 5 xi 6 yi j 4 5 j v x 6 j (b) y m 4 5 m 6 m 利用线性代数中解方程组的克来姆法则,由(a)可解出待定常数 、 、 1 2 : A1AA A 1 A A 2 A A 3 A 1 2 3 式中行列式: u x y A ui xi yi 1 u y i i A 1 u 1 j j u x y j 2 y j j 1 u y m m m 1 x u i i A m m 1 x y i i 3 1 xj uj A 1 xj yj 2A 1 xy1 x 1 x y m m m m A为△ 的面积,只要A不为0,则可由上式解出: ijm 1 (a u a u a u ) 1 2A i i j j m m 1 (bu b u b u ) （C） 2 2A i i j j m m 式中: 1 (cu c u c u ) 3 2A i i j j m m a x y x y i yj m y m j a x y x y mi mj m i y i m i m a x y x y m i j j i b i xj m bj y x b y y ii j i  （d） mc c x x m i j j i m cm xj x a为了x书y写方便x，y可将上式记为: i yj m y m i m(i,j,m ) m xxi j x x m i j (i,j,m ) 表示按顺序调换下标,即代表采用i,j,m作轮换的方式便可得到(d) 式。 将(c)式代入2-1-2中,整理后可得: u N (x,y)u N (x,y)u N i i j j m (x,y)u m 同理: v N (x,y)v N (x,y)v N m (x,y)v 式中: i i j j m （2-1-2）b N 1 (a b x c y) (i,j,m ) i 2A i i i (2-1-3) 将三角形单元的位移函数用矩阵表示： u i Nv N f(x,y) u N i 0 N j 0 i m 0 u j  (2-1 -4)a v 0 N i 0 N j 0 N m vi u m v m 或： u { f} v [N ]{d}e (2-1-4)b 三、单元的应变和应力 １、应变──几何矩阵[B] 由弹性力学知,弹性力学平面问题的几何方程： u x x ； y v u v y ； xy y x 用矩阵表示 x 0 x 0 u y y v xy y x 或， H f (2-1-5) [H]称为微分符矩阵,又称为微分算子,“[H]{f}”实际上不是一般的矩阵乘,可以称为微分符矩阵[H]作用在{f}上,其作用规律符合矩阵乘积规定，实际上是按[H]对{f}求导。 将2-1-4式的{f}=[N]{d}代入: {ε}=[H][N]{d}=[B]{d} 2-1-6 x 0 式中: B H N 0 N i 0 N 0 N 0 j