有限元作业第二次作业

1、一 五年 6 月 12 日土木工程专业有限元第二次作业姓名：班级：学号：指导教师：题： 平面应力问题的八节点等参元，已给定8 个节点实际单元4 = 1 73A6= 115 =121、单元中位移函数u（ ， ） ， v（ ， ）和单元节点位移 e 的关系式；2、 B 矩阵的计算步骤和计算式；3、单元刚度矩阵 k e 的一般计算方法和计算步骤；4、论述相邻单元间公共边界上位移的连续性；5、如果给定母单元中点A， （ ， ） ，怎样求实际单元中与A， 相对应的点A（ x， y） ；反之，如果给定实际单元中的点A（ x， y） ， 怎样求其在母单元中对应点A， （ ， ） ？6、如果已经求解得到单元8

2、 个节点的位移值 e 怎样求单元中某一点B（ x， y）的应力？解：图 1： 在总坐标系中具有二次曲边的四边形单元1、 此题分两步进行：单元位移场的表达：如图 1所示 , 在任意四边形的每边中间设一附加节点，则单元边界就变成二次曲线的了。如果直接在整体坐标系 x, y 下， 像八节点矩形元那样，构造双二次多项式的位移插值函数，则因曲边四边形单元边界是二次曲线，故边界上的位移是x（或 y）的五次多项式，它不能由曲边上三个节点的位移分量唯一地决定，从而不能保证相邻两个单元在公共边上位移的协调条件，所以在整体坐标系x,y 下构造完全协调的位移插值2： 在自然坐标系中的 曲边四边形的基本单元函数是很困

3、难的，利用坐标变换，可将曲边四边形单元变换成基本单元，如图2所示的在自然坐标, 下具有边长为2的八节点正方形单元， 自然坐标系, 是外节点坐标值为1的局部坐标系。在自然坐标系的单元上构造协调的位移插值函数，其形状函数是较普通的，取位移分量为, 的双二次多项式, 即：2222ua1a2a3a4a5a6a7a8aaaa2 aa 2a2 a2 （ 1-1 ）v a9a10a11a12 a13a 14a 15 a 16利用 8 个节点的16 个位移分量可唯一确定16 个待定常数a1,a2, , a16 ， 若代入 8 个节点的局部坐标值，得：7u1u5u2u6u3u7u4111111-101110-1

4、-1-1-1011110111011 1 -1 -1a101 00-1 1 -1 -1010-101110101010-1a2a3a4a5a6a7u81-10a8-1-11 1 -1 -1a91-2)v2v6v3v7v40-1-1-10-11001 -1 -10-10-1a10a11a12a13a14a151-3)-1将解出的16 个待定常数也即：u N1I N2I N8I e N1-4b )0a16a1,a2, a16 代入式(1-1 )即得 :8uN1u1N5u5N2u2N6u6N3u3N7u7N4u4N8u8Niuii18vN1v1N5v5N2v2N6v6N3v3N7v7N4v4N8v8

5、Nivii1其中 I 为二阶单元矩阵，e为等参元节点位移列阵，N 为形状函数矩阵。形状函数的建立：按等参元思想，在整体坐标系XY 下 , 任何形状歪斜四边形单元都将变换到局部坐标系下的正方形单元。对 8 节点等参元, 其移模式为：8u Ni , uii1式中 , ui为歪斜单元8节点的位移，Ni, 为形状函数。查阅相关资料，得形函数公式公式为：Ni ,Fk ,k1 81-6)Fk i , i k1又由形状函数的性质可具体地求出Ni 的表达式为：N1= 111 4N2= 111 4N3= 1114N4= 111 4N5= 12 121-7)N6= 1 根据平面问题的几何方程，单元应变可用节点位移

6、表示如下： 12N7= 12 12N8= 1212其中：另根据符合函数求导法则，可知：y = B e = B1B2B8e2-1 )xyBi =x0 Ni2-2)NiyBi 中的元素NixNi(i y1,2, ,8)。x2-3)yyx 8 Ni x 8 Nii xi ，ii1ii12-5)又根据式（2-3） ，有2-6)13根据公式（2-2）即可得出Bi 矩阵，其中Ni 可由问题1 方法求出。3、单元刚度矩阵按普遍公式计算，公式如下：k eBTDBdVBT DBhdxdy（ 3-1 ）ee其中e为单元体积域，k e为 16 16的方阵 （具体形式见下文）， D 为材料1-D E11123-2）1

7、-1-2（2 1- ）上述积分应在局部坐标系内进行，因此面积元素dxdy需表示成d d . 如图3 所示为子单元内任一点a x, y 处的微小正方形，它是由局部坐标系中点,ab, ac分别由处的微元体d d 变换而成的。以i, j表示 x, y轴的单位基矢量，d , d 变换而成，则：xyab= i j dxyac = i j d3-3）上述 2 个矢量的叉积表示它们所构成的平行dV ab ac Jd d （ 3-4）其中，J 为矩阵 J 的行列式，即J3： 子单元内任一 点处的微小正方形将上式带入式（3-1 ） ，并写成分块形式：k12k1k22k283-5）k82k88其中子矩阵的计算公式

8、为：11kijBiTDBjhdVBiTDBjJhd d （ 3-6）11e其中 h 是板的厚度。由于被积函数极为复杂, 很难得到明显的解析式, 必须利用数值积分。程序中采用高斯求积法, 对于二维问题的等参元, 高斯求积公式为:11881 1f , d dHiHjf i, i （ 3-7）2 1i1i1式中， H i ， H j 为一维求积点的积分系数，i ， i 为沿一维编号的求积积分点的横坐标。对于8节点等参元取三个积分点，即n=3已足够精确。4、 证明：局部坐标系下的单元是规则的正方形，单元边界上的三个节点按线性变化的位移形式，单元变形后这三个节点确定了位移的单元直线边界。所以，局部坐标系

9、下单元是协调的。又由位移插值函数在局部坐标系下的协调性，即可推知坐标变换的协调性（即两个相邻曲边四边形在公共边界上经坐标变换后仍保持连续，不会出现重叠和破缺现象），这也就保证了位移插值函数在整体坐标系下的协调性。即在相邻单元公共边界上位移是连续的。5、 这里， u,v是 , 的函数，在下面的计算中还需知道u,v和 x, y的关系，因此必须写出x, y 和 , 之间的坐标变换式，这个坐标变换并不难，因为 x, y 在单元的 8 个节点上应取值xi , yi （i1,2 ， 8）， 而单元四条边应为二次曲线，这与u, v的要求完全类同，因此可沿用和位移插值函数完全相同的式子作为坐标变换式，即：y1

10、（ 5-1）x1xN10N50N20 N60 N3 0 N7 0 N40N80y0N10N50N2 0N6 0 N3 0 N7 0N40N8y8式中x1, y1,x2, y2, , x8, y8 为节点坐标，形状函数N1, N2, N 8与前面相同。由上可见，在整体坐标系下的曲边四边形单元和自然坐标系下的正方形单元存在着一一对应的映射关系，只要已知xi ,yi（i1,2 ， 8）后，由（ 5-1 ）式，利用自然坐标系下的形状函数，即可完全确定x, y 。即：如果给定母单元中点*A , ，通过求出形状函数Ni, （i 1,2，8），利用式（5-1 ） ，可求出实际单元中与A* 相对应的点A x, y ； 同理， 如果给定实际单元中的点A x,