第三讲补充 等参数单元.pptVIP

作业: 简述等参单元的概念及单元分析过程. * 第一节 等参元的概念 前面介绍的三角形单元和四面体单元，其边界都是直线和平面，对于结构复杂的曲边和曲面外形，只能通过减小单元尺寸，增加单元数量进行逐渐逼近。这样，自由度的数目随之增加，计算时间长，工作量大。另外，这些单元的位移模式是线性模式，是实际位移模式的最低级逼近形式，问题的求解精度受到限制。 　　为了克服以上缺点，人们试图找出这样一种单元：一方面，单元能很好地适应曲线边界和曲面边界，准确地模拟结构形状；另一方面，这种单元要具有较高次的位移模式，能更好地反映结构的复杂应力分布情况，即使单元网格划分比较稀疏，也可以得到比较好的计算精度。等参数单元（等参元）就具备了以上两条优点，因此，得到广泛应用。 等参元的基本思想是：首先导出关于局部坐标系的规整形状的单元（母单元）的高阶位移模式的形函数，然后利用形函数进行坐标变换，得到关于整体坐标系的复杂形状的单元（子单元），如果子单元的位移函数插值结点数与其位置坐标变换结点数相等，其位移函数插值公式与位置坐标变换式都用相同的形函数与结点参数进行插值，则称其为等参元。 用形函数单元位移模式 坐标变换式 三节点三角形单元 四节点矩形单元 §4-2 平面四节点四边形等参单元 一、正方形单元 如图，设各节坐标为四点 i (-1,-1) j (1,-1) k(1,1) l (-1,1) 在单元内假定 因为 图（一） 代入节点处的位移连续性条件有: 共有8个方程,可以解出{ }.整理后可表示为 其中 称形函数,形函数有两个重要性质 2、4个之和为1. 对应于 r s 写成矩阵形式有 其中 设单元厚度为t，则以上积分可以表示为 以上积分域边界规则,很易计算 对复杂问题也很容易写成以下数值积分 其中： ——高斯积分点坐标； ——高斯积分系数； 对于任意四边形单元,可以通过一种特定的坐标变换,将其变换到一个矩形区域。这样就可以利用上述过程,计算其单元刚度 二、等参变换: 令 其中： 即正方形单元的形函数。该变换可将一任意四边形区域变换到一个矩形区域,现在看任意四边形该变换将任意四边形区域变换到正方形区域 三、单元分析 在单元内部分 假定： 写成矩阵形式为 其中： 这里Ni （1，2，3，4）是r s坐标下的形函数. 要注意 Ni是r、s的函数。ｒ、ｓ是ｘ、ｙ的隐函数（或反函数）。 因而微分运算要用以下连锁规则 写成矩阵形式有 ——Jacobi 矩阵 所以有 将坐标变换式代入 为计算方便,在r、s坐标下计算以上积分，即利用等参变换公式进行 变量替换， 则有 以上积分表示为 写成高斯积分形式为 单元体力载荷向量可以表示为 写成高斯积分形式为 同理可以得出等效面力载荷向量 应力由 及 表示 注意 不是常矩阵,各元素均为 的函数。在单元内不同点上， 元素的值也不相同.要计算哪点的应力,在单元分析时就要同时给 出该点处的 。要计算多点，就要给出多个 。由于高斯积分点 的 在计算单元刚度时必须计算，为减少计算量，可以直接用来 计算应力。所以很多应用程序一般都输出高斯点的应力。 边界位移连续性讨论 以 边为例 ，在边上有 其中 均退化为的一次式，此时由 可知x,y都均退化为 的一次式，所表示的正是直线的参数方程。 因为该直线必然通过 点,且由量点确定的直线是唯一的确定 实际计算中应注意的问题 1、旋向，指推导时I,J,K,L。 2、单元形状 ，线性应变单元。 §4-3 平面八节点四边形单元 一，等参数变换 其中 令 图（三） 二，单元分析 在单元内部分 假定： 图（三） 成矩阵形式为 其中： 这里，（i=1，2，3，4,5,6,7,8）是r、s的函数，r、s与x、y的关系由 等参变换公式给出。 这里 （1，2，3，…8）是r s坐标下的形函数. 是 的隐函数 因而微分运算要用以下连锁规则 写成矩阵形式有 ——Jacobi 矩阵 所以有