第三章 等参单元.pptVIP

船舶结构有限元分析Finite Element Analysis of Ship Structure;3.1等参单元的引入;3.2四节点矩形双线性单元3.21位移插值函数; 式中所示的位移插值函数的特点：固定x，它 是y的线性函数；固定y，它是x的线性函数。因而， 我们称它为双线性插值函数，把采用这种位移插值 函数的四节点矩形单元称为矩形双线性单元。 矩形双线性单元有着常应变三角形单元所不具 有的反映单元弯曲变形的能力，从而计算精度有望 提高。 由于上式满足常应变准则即完备性条件，又满 足位移连续性条件即协调性条件。因而，矩形双线 性单元是完备的协调单元。;3.2.2形函数的确定;根据条件2，可得 根据矩形单元各个节点的坐标，可将形函数Ni统一表示为 （i=1,2,3,4） 对任意双线性函数f(x , y),都有 为f(x , y)在节点i的值，分别取f(x , y)=1， f(x , y)=x和f(x , y)=y，得： ; 形函数是定义在给定单元上的，所以上式对单元的每一点（x , y）都成立。单元内一点的坐标通过单元节点坐标 的表达式，即坐标插值函数。 位移插值函数和坐标插值函数具有完全相同的构造，都是以节点处的函数值作为参数，并具有相同的形函数。因而矩形双线性单元也是一种等参数单元。;3.2.2矩形双线性单元的完备性和协调性;协调性：只要相邻单元在其公共边界上具有共同节点，相邻单元间公共边界上任一点的位移就必定相同，所以位移连续条件（即协调性条件）得到满足。 在x=1边界上，N1=0，N4=0，有 同理，在x=1的边界有 。可知，矩形双线性单元的边界在手里变形后，仍保持为直线边界。此协调性论证适用于任何等参单元。;3.3四节点四边形等参单元;对于a所示的四边形单元，将各边的等分点用之线连接，并规定它与局部坐标系下单元的相应对边、等分点的连线对应，这样，就建立了图a所示四边形单元和图b所示正方形单元个点间的一一对应关系。;3.3.2自然坐标系下的位移插值函数;3.3.3坐标变换;上式是一条以ξ为参数的直线方程，因而变换式把母单元的12边变换成整体坐标下连接节点的直边。对其他边也有同样结论。 位移插值公式和坐标变换公式具有完全相同的构造，因此这种四节点任意四边形单元也是等参单元，简称为四节点等参数单元，在整体坐标下也必然满足常应变准则。 位移插值函数在自然坐标系下的协调性自动保证了坐标变换式的合理性及位移插值函数在整体坐标系下的协调性 今后，对每个单元，我们把立足点放在自然坐标上，一切计算都转换到自然坐标系下进行;位移插值函数可用矩阵记号记为： 3.3.4应变矩阵B建立 对平面应力问题;上式B是一个三行8列矩阵 （i=1，2，3，4） 利用复合函数求导法则，得 写成矩阵形式，;式中坐标变换矩阵（雅可比矩阵） 由矩阵求逆法则得 式中;3.3.5应力—应变关系;变换后微元面积记为dA，则 2.单元刚度矩阵 把积分区域变换到自然坐标系平面上的相应区域 ;3.等效节点载荷;3.3.7等参变换的条件;3.4等参单元的收敛性;关于单元的安全性，对于 型单元，要求插值函数中包含完全的线性项，本章讨论的所有单元在自然坐标中是满足此要求的。 研究经过等参变换后，笛卡尔坐标下有： 将上式代入，就得到单元内的函数表达式 可以看出，如果插值函数满足 ，说明单元满足完备性要求;某些情况下，坐标插值多项式次数低于位移插值多项式，计算更为便捷，这样的单元称为次参数单元，如图为典型的二维次参数单元。 坐标插值函数的阶数高于位移插值函数的阶数的单元，称为超参数单元，可以更准确的描述壳体的复杂几何形状;3.5数值积分方法;记 则有 插值多项式可写为 于是有 若记 就得到插值求积公式;式中的系数Hi(i=1,2,3…n)称为求积系数或加权系数。它的值与被积函数无关，只依赖于求积节点的个数和位置。 当F(ξ)本身是次数不超过n-1的多项式时， 即具有n个求积节点的插值求积公式至少具有n-1次的代数精确度。 3.5.2高斯求积法 对求积节点的位置予以优选，使插值求积公式对任何次数不超过2n-1次的多项式函数都精确成立。这时求积公式的最高代数精确度为2n-1次，称为高斯求积公式;设F(ξ)为次数不超过2n-1次的任意多项式， 因此有 只需要选取ξi(i=1,2,…，n)，使以下n个等式 成立，则插值求积公式 即为高斯求积公式。积分点的个数n，称为数值积分的阶 ;N次勒让德多项式 在区间[-1，