等参单元概述.pptVIP

等参数单元 若以母单元上12边为例，通过映射可得在平面内任一直线，12边的方程为 η= -1,代入 同理可得 * * 前面介绍的三角形单元和四面体单元，其边界都是直线和平面，对于结构复杂的曲边和曲面外形，只能通过减小单元尺寸，增加单元数量进行逐渐逼近。这样，自由度的数目随之增加，计算时间长，工作量大。另外，这些单元的位移模式是线性模式，是实际位移模式的最低级逼近形式，问题的求解精度受到限制。 　　为了克服以上缺点，人们试图找出这样一种单元：一方面，单元能很好地适应曲线边界和曲面边界，准确地模拟结构形状；另一方面，这种单元要具有较高次的位移模式，能更好地反映结构的复杂应力分布情况，即使单元网格划分比较稀疏，也可得到较好的计算精度。等参数单元（等参元）就具备了以上两条优点，因此，得到广泛应用。 §6-1 等参元的概念 第六章 等参数单元 每一种等参数单元有二种形式：一种形式在局部座标 ξ η 或 ξ η ζ 中，是以边长为2的正方形单元或边长为2的立方体单元用作计算称为母单元；另一种形式在整体座标 x y或x y z中，通过母单元映射到整体座标上的单元用于离散结构物称为子单元。由于母单元上的位移函数与座标变换的函数具有相同的参数。这就是等参数单元名称的来由。 母单元的位移函数: 座标变换的映射关系: 等参元的基本思想：首先导出关于局部坐标系的规整形状的单元（母单元）的高阶位移模式的形函数，然后利用形函数进行坐标变换，得到关于整体坐标系的复杂形状的单元（子单元），这种子单元的位移函数插值公式与位置坐标变换式都用相同的形函数与结点参数进行插值，称其为等参元。 一、形函数 　　对于单元形函数的确定，首先假设单元的位移模式，代入结点的位移和坐标，从而推导出单元的任意一点的位移插值函数，即形函数。实际上，形函数是定义在单元内部的、满足一定条件的、坐标的连续函数。形函数不仅可以用于单元位移函数的插值，还可以用于单元形状的变换。形函数应满足的条件是： (6-1) 1.? 在结点i处 ， 在其他结点处 ； 2. 能保证用它定义的未知量（位移或坐标）在相邻单元之 间的连续性； 3. 应包含任意线性项，以保证用它定义的单元位移可满足常应变条件； 应满足下列等式 以保证用它定义的单元位移能反映刚体位移。 　　如图6-1所示，坐标原点在单位形心上。单元边界是四 条直线： ， 。为保证用形函数定义的未知量在相 邻单元之间的连续性，单元结点数目应与形函数阶次相适应 。因此，对于线性、二次和三次形函数，单元每边的结点数 分别为两个、三个和四个。除四个角点外，其他结点位于各 边的二分点或三分点上。 二维母单元是平面中的2×2正方形 二、母单元 　　首先，根据形函数的定义，在局部坐标中，建立起几何形状简单且规整的单元，称为母单元。 图6-1 二维母单元 (a) 线 性 单 元 (b) 二 次 单 元 1 2 3 4 1 2 3 4 8 7 5 6 1) 线性单元（4结点） (6-3) 以上形函数也可以合并表示为 （i=1, 2, 3, 4） (6-4) 其中 (6-5) 三、坐标变换 　　母单元可以直接用来进行有限元分析，其单元特性可以按照前面几章中讲述的步骤进行。但是这些单元形状规整，难以适应实际工程中出现的各种结构的复杂形状。为了解决这个矛盾，需要用坐标变换的方法，把形状规整的母单元，转换成具有曲线（面）边界的形状复杂的单元。转换后的单元称为子单元。子单元在几何上可以适应各种实际结构的复杂外形。这样，对于一个实际结构，就可以采用各种形状复杂的子单元在整体坐标系中进行划分，来逼近其复杂的曲线或曲面边界。而每个子单元，通过坐标变换，都可以映射成一个局部坐标系下的规整单元，即母单元，计算比较简单。 　　为了进行坐标变换，必须在局部坐标 和整体坐标 之间建立一一对应关系。这种对应关系可以利用形函数建立起来。 1. 平面坐标变换 在整体坐标系中，子单元内任一点的坐标用形函数表示如下 (6-4) 其中，Ni 是用局部坐标表示的形函数，xi，yi 是结点i的整体坐标，上式即为平面坐标变换公式。 返回 通过局部座标与整体座标的映射关系把母单元变换到整体座标上成为一个任意四边形用于离散结构物，它能适合于任意曲边的形状。 把以参数 ξ 代表的 x方程和 y方程消去 ξ ，则得