有限元作业第二次作业

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

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

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

4、a2u211-11-11 -1 -1a3u611010000a4u311111111a5（1-2 ）u710100100a6u41-1 11-11 1-1a7u81-1010000a8v11-1-11 11 -1-1a9v510-100100a10v211-1 1-11 -1-1a11v611010000a12v311111111a13（1-3 ）v710100100a14v41-1 11-11 1 -1a15v81-1010000a16将解出的 16 个待定常数 a1 , a2 , , a16代入式（ 1-1 ）即得 :8u N1u1N5 u5N2 u2N6 u6N3u3N7 u7N4u4N

5、8u8Ni uii 1（ 1-4a）8v N1v1N5 v5N2 v2N6 v6N3 v3N7 v7N4 v4N8v8Ni vii1也即：uN1I N2 IN8 IeeuNv（1-4b ）其中 I 为二阶单元矩阵，eN 为形状函数矩阵。为等参元节点位移列阵，形状函数的建立：按等参元思想，在整体坐标系XY 下,任何形状歪斜四边形单元都将变换到局部坐标系下的正方形单元。对 8 节点等参元 ,其移模式为：8uNi ,ui（1-5 ）i 12式中 ,ui 为歪斜单元8 节点的位移， Ni, 为形状函数。查阅相关资料，得形函数公式公式为：8Fk,Ni,k1（1-6）8Fki, ik1又由形状函数的性质可

6、具体地求出Ni 的表达式为：N1= 11N2= 11N3= 11N4= 11N5= 12N6= 12N7= 12N8= 121 41 41 41 412（1-7）121 21 22、 根据平面问题的几何方程，单元应变可用节点位移表示如下：xeey = BB8=B1 B2（2-1）xy其中：Ni0xBi =0N i（2-2 ）yNiNiyx即要求出矩阵 Bi 中的元素Ni， N i (i1,2,8) 。xy另根据符合函数求导法则，可知：3NixyNiNi=xx（2-3 ）NiyNiJxNiyy其中， J 为二维坐标变化下的 Jacobi 矩阵，即：xyJxy（2-4）其元素计算式为：xy8N i

7、 xi ，xi 18N i yi ，yi 18N i xi ，i 18N i yi（2-5）i 1又根据式（ 2-3 ），有NiNixJ 1（2-6）NiNiy根据公式（ 2-2 ）即可得出 Bi 矩阵，其中 Ni 可由问题 1 方法求出。3、 单元刚度矩阵按普遍公式计算，公式如下：keBT DBhdxdy（3-1 ）BT DBdVee其中e为单元体积域， ke为 16× 16 的方阵（具体形式见下文）， D 为材料的弹性系数矩阵，各向同性材料的弹性系数矩阵为：4101-E 110（3-2 ）D1121-001-2（ 1）2 -上述积分应在局部坐标系内进行，因此面积元素dxdy 需表

8、示成 d d. 如图3 所示为子单元内任一点a x, y 处的微小正方形，它是由局部坐标系中点,处的微元体 d d变换而成的。以 i, j 表示 x, y 轴的单位基矢量，ab, ac 分别由d , d 变换而成，则：ab=x iy j dac =x iy j d（3-3 ）上述 2 个矢量的叉积表示它们所构成的平行四边形面积，故：dVabacJd d（3-4 ）xy其中， J 为矩阵 J 的行列式，即 Jxy图 3：子单元内任一将上式带入式（ 3-1 ），并写成分块形式：点处的微小正方形k11k12k18ek21k22k28k（3-5 ）k81k82k88其中子矩阵的计算公式为：BiT11k

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

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

11、标变换式，即：x1xN10N50 N20 N60N30N70 N40 N8y10（5-1 ）y0 N10N50 N20 N60 N30 N70 N40 N8x8y8式中 x1 , y1 , x2 , y2 , , x8 , y8 为节点坐标，形状函数 N1 , N2 , , N8与前面相同。由上可见，在整体坐标系下的曲边四边形单元和自然坐标系下的正方形单元存在着一一对应的映射关系，只要已知xi , yi (i1,28) 后，由（ 5-1 ）式，利6用自然坐标系下的形状函数，即可完全确定x, y 。即：如果给定母单元中点A*,，通过求出形状函数Ni,(i1,28) ，利用式（ 5-1 ），可求出实际单元中与 A* 相对应的点 A x, y