第五讲空间等参薄板.ppt

6-1等参数单元2(1,-1)3(1,1)4(-1,1)1(-1,-1)1、等参变换 将局部坐标下的规则形状单元转换为总体坐标下几何形状扭曲的单元，以满足任意形状离散的要求。*第31页，共62页，星期日，2025年，2月5日6-1等参数单元局部坐标总体坐标变换函数*第32页，共62页，星期日，2025年，2月5日6-1等参数单元变换实例xyzt3(1,1)4(-1,1)2(1,-1)1(-1,-1)y=1=1=-1=12(x2,y2)1(x1,y1)3(x3,y3)4(x4,y4)uvP(x,y)*第33页，共62页，星期日，2025年，2月5日6-2等参数单元形函数2、形函数的性质同前在矩形单元上的变化如图注意：不是直线*第34页，共62页，星期日，2025年，2月5日6-2等参数单元形函数形函数N1的正确表示直线直线（1,-1）不是平面-1-0.500.51-1-0.500.5100.250.50.751-1-0.500.51*第35页，共62页，星期日，2025年，2月5日6-3等参数单元位移函数3、等参单元位移函数：单元内任意点p的位移函数(2D)： 其中：Ni和坐标变换式的形函数相同。*第36页，共62页，星期日，2025年，2月5日6-4等参数单元刚度矩阵4、等参单元刚阵 1）应变矩阵 注意：应变为位移对x，y的导数，如四节点四边形单元计算式如右：用于二维等参元*第37页，共62页，星期日，2025年，2月5日6-4等参数单元刚度矩阵2）复合求导 记为矩阵(如四节点四边形单元) [J]称为Jacobi矩阵，由坐标变换式确定，当[J]的逆存在时，则形函数对x，y的导数可求，即应变阵可求。*第38页，共62页，星期日，2025年，2月5日6-4等参数单元刚度矩阵应变矩阵*第39页，共62页，星期日，2025年，2月5日6-4等参数单元刚度矩阵3）刚度矩阵 一般而言，等参单元的刚度积分很难有解析式，必须进行数值积分，目前普遍采用高斯数值积分法。(略)*第40页，共62页，星期日，2025年，2月5日6-5空间六面体单元8(x5,y5,z5)1234(x4,y4,z4)5(x5,y5,z5)67xzy23(1,-1,1)48(1,1,-1)657*第41页，共62页，星期日，2025年，2月5日6-5空间六面体单元1)形函数(i=1,2,…8)其中：例：*第42页，共62页，星期日，2025年，2月5日6-6等参数单元说明等参单元的几点说明：1）等参单元为协调元，满足有限元解收敛的充要条件。证明略。2）等参单元存在的充要条件是： 为了保证能进行等参变换(即总体坐标与局部坐标一一对应)，通常要求总体坐标系下的单元为凸，即不能有内角大于或等于或接近180度情况。*第43页，共62页，星期日，2025年，2月5日6-6等参数单元说明3）等参单元的优点是当单元边界呈二次以上的曲线时，容易用很少的单元去逼近曲线边界。4）上述等参单元的理论公式可适应三次以上的曲线型等参元，只是阶次提高，单元自由度相应增加，计算更复杂，积分更困难，实际中，很少超过3次曲线型。5）上述推导要求：保持坐标变换中几何模式阶次与描述单元位移函数中形函数的阶次相同。如取坐标变换的几何模式阶次较单元的位移函数的阶次高，则称此单元为超单元，反之，为亚单元。这两类单元的收敛性也可得到满足。略6）当然，也可取描述单元几何形状的几何模式不是形函数的，如p-element*第44页，共62页，星期日，2025年，2月5日第五讲空间等参薄板*第1页，共62页，星期日，2025年，2月5日空间与轴对称问题等参数单元薄板弯曲问题*第2页，共62页，星期日，2025年，2月5日5-1空间问题简介工程实际中的很多问题难于简化为平面问题，如受任意空间载荷作用的任意形状几何体，受对称于轴线载荷作用的回转体，本章简单介绍两类问题：轴对称问题和空间问题的有限元计算。空间问题的主要困难：（1）离散化不直观；————（网格自动生成）（2）未知量的数目剧增。————（对某些问题简化）——————（轴对称问题）空间分析的优点：精确。*第3页，共62页，星期日，2025年，2月5日5-1空间问题简介某一平面