第4章等参数单元(免费阅读).pptVIP

* * 第4章 等参数单元 4.2 等参数单元概念 4.3 平面等参数单元 4.1 自然坐标和插值函数 第三章有限元分析的思路，十分严谨、明确，但是过于繁琐。对于常应变三角形单元还可行。 如二维8节点四边形，16个待定系数，人工几乎难以完成； 三维20节点单元，20×3个待定系数。 解决途径： 形函数是单元位移的插值函数，可以利用数学中的“插值函数”直接给出形函数。 4.1 自然坐标和插值函数 （局部）自然坐标 取2×2正方形（母单元）的两对边中点作为该正方形的局部坐标系，以Oξη表示，并且规定ξ，η的坐标值有效取值范围在0到±1之间。 坐标值：在中心点（原点）：ξ=η=0； 正方形的各边：ξ=±1，η=±1； 各角点：ξ=1或-1， η=1或-1； 正方形内部的任一点：由局部坐标来定义， ξ、η在±1之间，即P（ξ,η)。 任意四边形单元（子单元）：定义同上。 二者具有相互映射的关系。 自然坐标：这种基于四边形形状而定义的坐标系，是同任意四边形的自然形状相关的，故称之为“自然坐标”。 如果在四边形的各边中点增加节点，即构成更高精度的单元。 曲边8节点四边形单元，相应的自然坐标ξ,η轴也呈曲线形。同样与母单元具有相互映射的关系。 以此类推，更多节点的任意四边形（12节点，16节点等）也具有同样的特性。 高精度四边形单元： 自然坐标的优点： 1.无论四边形的大小和形状如何，其坐标特性是相同的，因此可以用统一的表达式来描述，这有利于推导统一的有限元公式； 2.以自然坐标导出的有限元公式，有利于实现数值积分运算，这就克服了高精度单元的刚度矩阵、等效节点力矩阵等因无法导出显式而必须进行积分所遇到的困难； 3.高精度（8节点及以上）单元，单元边界可以是直边或曲边，就可以更好地逼近实际物体复杂的几何形状。 拉格朗日插值:若函数F(x)在m个点上有已知值，则可以用在这些点上的已知值构造一个不超过n=(m-1)阶的代数多项式来近似。该函数f(x)在所有已知点处均满足f(xi)=F(xi)。 记为： 拉格朗日插值 式中 称为n阶拉格朗日插值函数，也叫“拉格朗日乘子”，可由下式求得： 可简写为 上标n表示函数的幂次，下标i表示对应的已知点（或节点）。 1.一维拉格朗日插值 两点插值： 三点插值： 高次的拉格朗日插值，均可由最简单的线性插值的某种组合来表示。 如采用自然坐标系，取 则可以把拉格朗日插值函数的公式表示如下： 两点插值： 三点插值： 2.二维拉格朗日插值 m,p分别为两个坐标方向上的插值点数目。 n=m-1,q=p-1，为插值函数的幂次。 通常都在两个坐标方向取相同数目的点。 若在ξ,η方向均取两点插值，则： 4.2 等参数单元的概念 平面问题的单元，最简单的是三节点三角形单元，其次是四节点矩形单元。 三角形单元适应性强，较容易进行网格划分，能应用与曲折的几何边界，但精度较低。特别是在应力集中部位容易产生较大误差。 矩形单元精度较高，形状规整，但适应性差，不便用于曲线边界和非正交的直线边界。应用不多。 等参数四边形单元：将整体坐标系（x,y)中的四边形（4节点、8节点，12节点等）单元，变换为局部(自然）坐标系（ξ,η）中的规则正方形，运用插值函数，建立位移模式，进行有限元分析，其坐标变换式和位移模式采用同样的形函数和相同的参数。 等参数的思想，由易到难，由规则单元的特殊情况推广到不规则单元的一般情况。 等参数单元应用最广，至今国际上流行的大型结构分析软件中，几乎都包含有等参数单元库。 应用实践表明，采用等参数单元离散结构，可以达到更高的计算精度，而且结构离散和数据准备工作量相对减少。 等参数三角形单元应用不多。 等参数单元的基本思想：首先导出规则单元（母单元）的形函数，然后进行坐标变换，导出对应的不规则单元的形函数和单元刚度矩阵。 等参数单元的优点： 1.应用范围广。在杆件结构、平面和空间连续体和板壳中都可应用。 2.将不规则的单元变换成规则的单元后，易于构造位移模式。 3.在原结构中可以采用不规则的单元，易于适用边界的形状和改变单元的大小。 4.可以灵活地增减节点，容易构造各种过渡单元。 5.推导过程具有通用性。一维、二维和三维的推导过程基本相同。 4.3 平面等参数单元 四节点四边形是最简单的四边形单元。 通过增加节点（从而也就增加了节点自由度及位移函数的幂次），可以构成8节点、12节点、16节点等高精度单元。 四边形等参数单元应用最为广泛。高精度三角形单元,由于其自身的某些不足,应用不够方便,相对较少。 利用拉格朗日插值方法，直接选定相应的形函数来构造位移模式非常方便。 四边形单元的类型 利用局部自然坐标来建立单元特性，可以避开繁琐的数学推演。 特别是对于高次的单元，由于不可能导出显式，