第五讲 等参数单元.pdfVIP

维普资讯 TechnicalLesson技术讲座 第五讲 等参数单元 Lesson5 The IsoparametricElements 陈乐生 (福州大学机械工程及 自动化学院) 从平面有限元的分析我们注意到，弹性体有 单元和棱边为直边的空间体单元还会产生较大的 限元求解过程中产生计算误差的地方主要有两 误差，大家比较认可的方法是采用等参数单元。 处：第一是单元计算误差，这是由于单元内假设位 等参数单元的基本思想是 ：先在局部坐标系 移场与物体真实位移场不一致造成的。第二是求 中对简单几何形状的单元 (称之为母单元)按照高 解基本方程时，计算机运算过程的舍人误差造成 阶插值多项式来构造形状函数，形成局部坐标系 的。虽然细分网格、增加单元数可以使假设的位 的单元位移函数；然后通过坐标变换 ，将简单几何 移场逼近真实位移场，但是另一方面，由于单元数 形状的母单元在总体坐标系中映射成实际网格划 增加，整体刚度矩阵扩大，又增加了舍人误差。因 分的曲边或曲面单元。参见图l。虽然经过坐标 此，提高单元精度最有效的办法是提高单元的计 变换之后单元的位移函数在总体坐标系中会变得 算精度，增加单元的节点数是其 中一种有效的办 比较复杂，但是我们可以通过数值积分的办法加 法。但是对某些曲边和曲面结构或者结构周边的 以解决。 曲率变化比较明显的构件，采用如同三角形 、矩形 厂] ． II I二_ J ．．．． ．．．一 图 1 对于二维等参数单元 ，局部坐标系中的母单 点的位移ui和vi(i=l，2，…，ITI)。这里，ITI为单元 元一般采用边长为2的正方形单元和直角边长为 节点总数。 l的等腰直角三角形单元。而空间体单元一般也 将(1)式简记为： 采用边长为2的立方体单元。这里，我们只简单介 f(，r／)= (，叼) ． 绍两种平面等参数单元的母单元的形状函数和坐 { ：‘‘ (2) 标变换，数值积分等问题不作介绍。 【(，r／)=∑ (，r／) 式中， (，叼)(i=l，2，^，m)称为插值函数 1 形状函数的性质 或形状函数。 现假设局部坐标系为考11，在该坐标系中母单 根据上一讲