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

土木工程专业
有限元第二次作业

姓 名：
班 级：
学 号：
指导教师：

二 〇 一 五年 6月12日
习 题：平面应力问题的八节点等参元，已给定8个节点的坐标。试查资料并论述：

1、单元中位移函数u（ξ，η），v（ξ，η）和单元节点位移{δe }的关系式；
2、[ B ]矩阵的计算步骤和计算式；
3、单元刚度矩阵[ k e ]的一般计算方法和计算步骤；
4、论述相邻单元间公共边界上位移的连续性；
5、如果给定母单元中点A，（ξ，η），怎样求实际单元中与A，相对应的点A（x，y）；反之，如果给定实际单元中的点Ax，y），怎样求其在母单元中对应点A，（ξ，η）？
6、如果已经求解得到单元8个节点的位移值{δe }怎样求单元中某一点B（x，y）的应力？

解：
图1：在总坐标系中具有二次曲边的四边形单元
1、此题分两步进行：
Ø 单元位移场的表达：
图2：在自然坐标系中的曲边四边形的基本单元
如图1所示,在任意四边形的每边中间设一附加节点，则单元边界就变成二次曲线的了。如果直接在整体坐标系下，像八节点矩形元那样，构造双二次多项式的位移插值函数，则因曲边四边形单元边界是二次曲线，故边界上的位移是的五次多项式，它不能由曲边上三个节点的位移分量唯一地决定，从而不能保证相邻两个单元在公共边上位移的协调条件，所以在整体坐标系下构造完全协调的位移插值函数是很困难的，利用坐标变换，可将曲边四边形单元变换成基本单元，如图2所示的在自然坐标下具有边长为2的八节点正方形单元，自然坐标系是外节点坐标值为±1的局部坐标系。在自然坐标系的单元上构造协调的位移插值函数，其形状函数是较普通的，取位移分量为的双二次多项式, 即：
（1-1）
利用8 个节点的16 个位移分量可唯一确定16 个待定常数，若代入8个节点的局部坐标值，得：
（1-2） （1-3）
将解出的16 个待定常数代入式（1-1）即得:
（1-4a）
也即：
（1-4b）
其中为二阶单元矩阵，为等参元节点位移列阵，为形状函数矩阵。
Ø 形状函数的建立：
按等参元思想，在整体坐标系下, 任何形状歪斜四边形单元都将变换到局部坐标系下的正方形单元。
对8节点等参元, 其移模式为：
（1-5）
式中, 为歪斜单元8节点的位移，为形状函数。
查阅相关资料，得形函数公式公式为：
（1-6）
又由形状函数的性质可具体地求出的表达式为：
（17）

2、根据平面问题的几何方程，单元应变可用节点位移表示如下： （2-1）
其中：
（2-2）
即要求出矩阵中的元素，。
另根据符合函数求导法则，可知：
（2-3）
其中，为二维坐标变化下的Jacobi矩阵，即：
（2-4）
其元素计算式为：
，，
， （2-5）
又根据式（2-3），有
（2-6）
根据公式（2-2）即可得出矩阵，其中可由问题1方法求出。

3、单元刚度矩阵按普遍公式计算，公式如下：
（3-1）
其中为单元体积域，为16×16的方阵（具体形式见下文），为材料的弹性系数矩阵，各向同性材料的弹性系数矩阵为：
（3-2）
上述积分应在局部坐标系内进行，因此面积元素需表示成.如图3所示为子单元内任一点处的微小正方形，它是由局部坐标系中点处的微元体变换而成的。以表示轴的单位基矢量，分别由变换而成，则：
（3-3）
图3：子单元内任一点处的微小正方形
上述2个矢量的叉积表示它们所构成的平行四边形面积，故：
（3-4）
其中，为矩阵的行列式，即
将上式带入式（3-1），并写成分块形式： （3-5）
其中子矩阵的计算公式为：
（3-6）
其中是板的厚度。由于被积函数极为复杂,很难得到明显的解析式,必须利用数值积分。程序中采用高斯求积法,对于二维问题的等参元,高斯求积公式为:
（3-7）
式中，，为一维求积点的积分系数，，为沿一维编号的求积积分点的横坐标。对于8节点等参元取三个积分点，即n=3已足够精确。

4、证明：局部坐标系下的单元是规则的正方形，单元边界上的三个节点按线性变化的位移形式，单元变形后这三个节点确定了位移的单元直线边界。所以，局部坐标系下单元是协调的。又由位移插值函数在局部坐标系下的协调性，即可推知坐标变换的协调性（即两个相邻曲边四边形在公共边界上经坐标变换后仍保持连续，不会出现重叠和破缺现象），这也就保证了位移插值函数在整体坐标系下的协调性。即在相邻单元公共边界上位移是连续的。

5、这里，是的函数，在下面的计算中还需知道和的关系，因此必须写出和之间的坐标变换式，这个坐标变换并不难，因为在单元的8个节点上应取值，而单元四条边应为二次曲线，这与的要求完全类同，因此可沿用和位移插值函数完全相同的式子作为坐标变换式，即：
（5-1）
式中为节点坐标，形状函数与前面相同。
由上可见，在整体坐标系下的曲边四边形单元和自然坐标系下的正方形单元存在着一一对应的映射关系，只要已知后，由（5-1）式，利用自然坐标系下的形状函数，即可完全确定。即：如果给定母单元中点，通过求出形状函数，利用式（5-1），可求出实际单元中与相对应的点；同理，如果给定实际单元中的点，利用式（5-1），即可求出的坐标值。

6、根据平面问题的本构方程，单元应变可用节点位移表示如下：
（6-1）
由式（3-2）和式（2-2）可分别求出矩阵和，故由上式，若已知单元8个节点的位移值，可求出单元中某一点B（x，y）的应力。
8