有限元第7章等参数单元讲义.ppt

例题7-1 ： 求该四变形单元的雅可比矩阵， 并求 取什么值时， 在单元内不变号？ 不变号。 以及用 例题7-2： 如图四节点四边形单元，各节点的坐标为， 求坐标变换公式 雅可比矩阵， 表示 (1)求坐标变换公式 由等参数元的坐标变换公式(7-4)，可得 (2)求 (3)求雅可比矩阵 (4)求 可以由第一步和第二步所得到的结果直接求出， 也可以通过雅可比矩阵求解 例题7-6： 求如图四节点等参数元在重力载荷作用下节点等效外载荷。重力载荷集度为 由题意可知： 各节点等效载荷为 已知图中单元2-3边受到线性分布载荷的作用，求节点等效外载荷。 解：由例题7-2可得到坐标变换关系如下： 由图可知线性分布载荷表示为 只需要考虑节点2，3的等效载荷， 第四节 试凑法来确定形函数 一、四节点正方形单元的形函数 坐标系下，单元的四条边界线的方程分别为 根据形函数应该具有的性质，即 本点为1，它点为0。 二、五节点正方形单元的形函数 五节点正方形单元，单元内任意一点的位移可表示为 满足本点为1，它点为0的特性。 从四节点矩形单元的形函数出发，为了避免混淆。用以下记号表示 在节点5 在1，2，4，5都为零，在点3为1，因此它满足条件，有 我们设法构造 它在1，2，3，4点处为零，以下的表达式能满足上述要求. 又需要满足在点5处为1。代入5点的坐标 是二次曲线，而该边上有三个节点，故能满足连续性要求。 在节点1为1，在点2，3，4为0，在点5 对图7-6所示的五节点正方形单元，我们可得到其形函数表达式： 三、八节点等参数单元的形状函数 对照四节点的情况，八节点正方形单元在局部坐标系下的位移插值函数 利用形状函数的性质，八节点等参数单元的形状函数可由下述两个条件所唯一确定： 1 是形如位移插值函数的双二次函数。 2. 在节点i其值为1，其余节点值为0。 由于形状函数是唯一的，则我们可以采用任意方法来求形状函数。 在节点2-8其值为零。 它在节点1和3-8应为0，注意到直线17、35、57通过这些点，而这三条直线的方程分别是 第五节 八节点曲边四边形等参数单元 任意四边形的四节点等参数单元可以较方便地对求解区域进行分割，但许多情况下仍嫌净度不够理想。一方面是因为位移插值函数是双线形函数，次数仍较低；另一方面因为整体坐标系下的任意四边形是直边四边形，对于具有曲线边界的求解区域的模拟仍有一定误差。 为进步提高精度，在四节点等参数单元的基础上增加节点数目，提高位移插值函数的阶次。使用中采用得最多的是八节点曲边四边形等参数单元。 一、平面八节点等参数单元位移插值函数 局部坐标下，八节点等参数单元仍然是边长为2的正方形。除了原来四节点单元的4个角节点外，又将各边中间点取为节点。节点排列和单元形状如图7-7所示。8个节点的局部坐标 因此这样的位移插值函数是双二次函数，相应的插值称为双二次插值。显然，它比双线性的位移插值函数阶次提高了，计算精度必然也会提高。 那么这一双二次的位移插值函数是否在单元边界上满足变形协调条件呢？ 二次函数，完全可以由该边上三个节点处的函数值所唯一确定，相邻单元的公共边上，三个节点位两相邻单元所共有。所以插值函数在此边上的连续性可以得到保证。 在此局部坐标系下单元变形的协调性条件被满足，再加上等参数坐标变换的相容性，则整体坐标下的变性协调性也将得到满足。 二、八节点等参数单元的形状函数 三、等参数变换极其实现条件 由形状函数以及等参数变换的思想，局部坐标到整体坐标的坐标变换式为 * 第七章 等参数单元 三节点三角形单元具有如下特点： 1. 位移插值函数为线性函数，因此称为三角形线性元。 2. 线性单元的位移在单元内呈线性变化，应力、应变在单元内是一个常量。 3. 应力和应变在求解区域内都不是连续的。 为提高计算精度，实际分析时可以采取的方法： 1. 单元分细； 2. 构造高精度新单元。 将单元分细可提高计算精度，因为有限元法的计算基础就是当单元无限分细时计算结果将收敛于精确解。但是单元分细会增加单元数目和节点数目，从而增加所要求解的方程组，占用和耗费大量的计算机资源。所以，用细分单元的方法来提高精度有时是不经济的。 构造具有较高精度的单元也可以提高计算精度。单元节点数增多，则自由度数目增多，允许采用较高阶次的位移插值函数，从而提高计算精度。例如6节点三角形单元其位移插值函数为完全二次多项式。 单元内位移为二次函数变化，应变和应力呈线性变化。但这种单元的面积小，节点