03-06空间问题的高次单元.doc

§3-6 空间问题的高次单元 和平面问题一样，三维问题也可按上述方式构造高次单元。对四面体单元系列可按某次完全多项式的形式构造单元，对六面体系列单元，可按不完全多项式的形式构造单元。 一、高次四面体单元 三元1、2、3次多项式分别有4、10、20项，所以对应的四面体N次单元取4、10、20个节点，取法应保证每条边上有2、3、4个点。如此构成1、2、3次四面体空间单元。 1. 10节点四面体单元 四面体二次单元构造 用直角坐标(x,y,z)表示的空间完全二次多项式共有十项，位移函数取为 单元的位移函数中有30个待定系数，为了确定这些函数，需要30个参数。今在四面体单元上取10个节点，每个节点有3个参数，正好可以确定位移函数中的系数。 如图取法，四个顶角点，加上六条棱边的中点。每条边上有3个点。保证连续性(如图3-24)， 分析 利用10个节点，每个节点有3个参数，正好可以确定位移函数中的系数。不过对于高次单元，如仍采用笛卡尔直角坐标，计算颇为复杂。 现改用体积坐标，单元法位移函数取为 (3-127) 其中ui，vi，wi，为节点i的位移分量，Ni为用体积坐标表示的二次形函数： 角点： 边点： (3-128) 利用单元分析的普遍公式，不难求出单元的应力、应变、刚度矩阵及节点载荷算式。 由几何方程及物理方程可知，此时单元内的应变及应力已不再是常量，而是坐标的线性函数。 2. 20节点四面体单元 构造 在直角坐标系中空间完全的三次多项式共有20项，四面体单元的位移函数如取为完全三次多项式，则需有20个节点共60个位移参量来确定其中待定系数。 今取四个角点，六条边的十二个三等分点及四个表面形式心为节点(如图3-25)。 分析 此时，用体积坐标表示的形函数为： 角点： 边点： 面形心点： (3-129) 采用单元分析普遍公式，可以求出单元应力、应变，单元刚度矩阵，节点载荷等算式。对20节点四面体单元，单元的应力和应变将是坐标的二次函数。 二、六面体单元 和在平面问题中采用矩形单元一样，空间问题中也可以采用六面体单元。 平面问题取所谓的双线性等，三维问题可取三线性等。 按项数取法对应的常用的六面体单元有8节点24自由度，20节点60自由度单元等。 1. 8节点六面体单元 对于8节点六面体单元(图3-26)每边二点，位移模式取三线性，取如下八项1,x,y,z,xy,yz,zx,xyz。 设单元边长在x,y,z方向分别为2a,2b,2c;选单元的形心为局部坐标系ξηζ的形心，令 则用局部坐标表示的单元形函数为 (i＝1，2，…，8) (3-130) 其中ξi，ηi，ζi为各点的局部坐标。 位移函数可表示为 2. 20节点六面体单元 二次单元：除了8个角点外，再引入12条棱边的中点为节点就得到如图3-27所示的20节点二次六面体单元， 按照上面同样的方法可以构造出用局部坐标表示的形函数： (3-131) 3．N12单元 构造 以上单元构造时均取节点位移为线位移为基本未知量。这种取法可保证相邻节点处位移连续，但并不能保证转角连续。转角相关于某个一阶导数。 基于此，每个节点除三个位移分量外，再引进位移分量的九个一阶导数(相当于应变)作为自由度，采用此种单元可以大大提高计算精度。每个节点自由度达到了12个，对应的力分量也做了相应的变化。 此类的六面体单元有8节点96自由度单元(H96，见图3-28)。四面体单元则有4节点48自由度单元。（T48）。 案例 图3-29表示了采用几种不同的空间单元计算一个悬壁梁弯曲问题的结果。 T１２—4节点12个自由度四面体元；T４８—4节点48自由度四面体单元；H２４—8节点24自由度六面体单元；Ｈ９６—8节点96自由度六面体单元。 该图以结构总自由度数为横坐标，以计算挠度和正确挠度的比值为纵坐标，在每个点旁注明了所用单元数目。 由此图可见： 高次单元的精度优于简单单元，而且六面体单元优于四面体单元。 但必须注意，在结构自由度数相同的条件下，高次单元刚度矩阵计算要花费较多机时。 六面体单元由于形体过于规则，与实际工程结构复杂的外形难于十分吻合，应用上受到限制。 在采用简单的T１２单元分析空间问题时，由于单元应力、应变是常量，不易适应变化的应力场，必须采用密集的网格、大量的单元，才能得到较好的计算精度。例如，对于悬臂梁之类以弯曲为主的结构，在厚度方向如果只布置三四排节点，计算误差将是很大的。 图3-29 空间单元计算悬臂梁结果比较 单元总结： 自由度：由M节点，每个节点L个自由度，计M*L个自由 度； 位移： 单刚LM阶。 总刚LNP阶。 4