有限元作业第二次作业.docxVIP

土木工程专业 有限元第二次作业 姓 名： 班 级： 学 号： 指导教师： 二〇一五年6月12日 习 题：平面应力问题的八节点等参元，已给定 8 个节点 的坐标。试查资料并论述： 4 7 3 A 8 6 Y 1 5 2 X 实际单元  4 η = 1 7 3 ˊ η A 8 ξ 6 ξ =﹣1 ξ= 1 1 5 η=﹣1 2 母单元 1、单元中位移函数 u（ξ，η），v（ξ，η）和单元节点位 移 { δe } 的关系式； 2、 [ B ] 矩阵的计算步骤和计算式； 3、单元刚度矩阵 [ k e ] 的一般计算方法和计算步骤； 4、论述相邻单元间公共边界上位移的连续性； ， 5、如果给定母单元中点 A（ ξ， η），怎样求实际单元中与 A， 相对应的点 A（ x，y）；反之，如果给定实际单元中的 ， 点 A（ x，y），怎样求其在母单元中对应点 A（ξ，η）？ 6、如果已经求解得到单元 8 个节点的位移值 { δe } 怎样求单 元中某一点 B（ x，y）的应力？ 解： 1、 此题分两步进行： 单元位移场的表达： 如图 1所示 , 在任意四边形的每 边中间设一附加节点， 则单元边界就变 成二次曲线的了。如果直接在整体坐标 系 x, y 下，像八节点矩形元那样， 构造双二次多项式的位移插值函数， 则因 图 1：在总坐标系中具有二 曲边四边形单元边界是二次曲线， 故边 次曲边的四边形单元 界上的位移是 x(或y) 的五次多项式， 它不能由曲边上三个节点的位移分量唯一地决定， 从而不能保证相邻两个单元在 公共边上位移的协调条件， 所以在整体坐标系 x, y 下构造完全协调的位移插值 函数是很困难的， 利用坐标变换， 可将曲边四 边形单元变换成基本单元，如图 2所示的在自 然坐标 , 下具有边长为 的八节点正方形 2 单元，自然坐标系 , 是外节点坐标值为± 1的局部坐标系。在自然坐标系的单元上构造 图 2：在自然坐标系中的 协调的位移插值函数， 其形状函数是较普通的， 曲边四边形的基本单元 取位移分量为 , 的双二次多项式 , 即： u a1 a2 a3 a4 2 a5a6 2 a7 2 a8 2 （ 1-1 ） v a9 a10 a11 a12 2 a13 2 a15 2 2 a14 a 16 利用 8 个节点的 16 个位移分量可唯一确定 16 个待定常数 a1, a2 ,?, a16 ， 若代入 8 个节点的局部坐标值，得： 1 u1 1 -1-1111-1-1 a1 u5 1 0 -100100 a2 u2 1 1 -1 1 -11 -1 -1 a3 u6 1 1 0 1 0000 a4 u3 1 1 1 11111 a5 （1-2 ） u7 1 0 1 0 0100 a6 u4 1 -1 1 1 -11 1-1 a7 u8 1 -1010000 a8 v1 1 -1 -11 11 -1-1 a9 v5 1 0 -100100 a10 v2 1 1 -1 1 -11 -1-1 a11 v6 1 1 0 10000 a12 v3 1 1 1 1 1111 a13 （1-3 ） v7 1 0 1 00100 a14 v4 1 -1 1 1 -11 1 -1 a15 v8 1 -1010000 a16 将解出的 16 个待定常数 a1 , a2 ,?, a16 代入式（ 1-1 ）即得 : 8 u N1u1 N5 u5 N2 u2 N6 u6 N3u3 N7 u7 N4u4 N8u8 Ni ui i 1 （ 1-4a ） 8 v N1v1 N5 v5 N2 v2 N6 v6 N3 v3 N7 v7 N4 v4 N8v8 Ni vi i 1 也即： u N1I N2 I N8 I e e u N v （1-4b ） 其中 I 为二阶单元矩阵， e N 为形状函数矩阵。 为等参元节点位移列阵， 形状函数的建立： 按等参元思想，在整体坐标系 XY 下, 任何形状歪斜四边形单元都将变换到 局部坐标系 下的正方形单元。 对 8 节点等参元 , 其移模式为： 8 u Ni , ui （1-5 ） i 1 2 式中 , ui 为歪斜单元 8 节点的位移， Ni , 为形状函数。 查阅相关资料，得形函数公式公式为： 8 Fk , Ni , k 1 （1-6） 8 Fki , i k 1 又由形状函数的性质可具体地求出 Ni 的表达式为： N1= 1 1 N2= 1 1 N3= 1 1 N4= 1 1 N5= 1 2 N6= 1 2 N7= 1 2 N8= 1 2  4 4 4 4 1 2 （1-7） 1 2 2 2 2、 根据平面问题的几何方程，单元应变可用节点位移表示如下： x e e y