有限元设计——单元位移模式及插值函数的构造.pptVIP

考虑微段dx，内力 dx的伸长为 x截面上的位移： 根据几何方程求应变，物理方程求应力 应变 应力 对于弹性力学平面问题，单元位移函数可以用多项式表示， 多项式中包含的项数越多，就越接近实际的位移分布，越精确。具体取多项，由单元形式来确定。即以结点位移来确定位移函数中的待定系数。 如图所示的3结点三角形单元，结点I、J、M的坐标分别为(xi,yi)、(xj,yj)、(xm,ym)，结点位移分别为ui、vi、uj、vj、um、vm 。 位移只能确定六个多项式的系数，所以3结点三角形单元的位移函数如下， 将3个结点上的坐标和位移分量代入（1）式就可以将六个待定系数用结点坐标和位移分量表示出来。 （1） 将水平位移分量和结点坐标代入（1）第一式， 写成矩阵形式， 令 则有 A为三角形单元的面积。 [T]的伴随矩阵为， 令 则 同样，将垂直位移分量与结点坐标代入公式（1）中的第二式，可得， 将两式代回（1）整理后可得， 令 （下标i，j，m轮换） 可得 单元内的位移记为 单元的结点位移记为 单元内的位移函数可以简写成， 把[N]称为插值函数矩阵或形态矩阵、形函数矩阵，Ni称为插值函数或形态函数、形函数。 必须注意：插值函数的构成不取决于求解的微分方程式，插值函数构造方法仅取决于：几何图形（单元形状）、 结点数量与位置以及在单元结点处规定的因变量的数量。 选择单元位移函数应满足以下条件： 1）反映单元的刚体位移与常量应变，称为完备性条件。 2）相邻单元在公共边界上的位移连续，单元之间不能重叠，也不能脱离。即位移函数在单元之间连续，称为协调性条件。 单元位移函数满足以上两个条件，就满足收敛性要求。 由（1）可以将单元位移表示成以下的形式， 反映了刚体位移和常应变。 单元位移函数是线性插值函数，因此单元边界上各点的位移可以由两个结点的位移完全确定。两个单元的边界共用两个结点，所以边界上的位移连续。 形态函数Ni具有以下性质： 1）在单元结点上形态函数的值为1或为0。 2）在单元中的任意一点上，三个形态函数之和等于1。 用来计算三角形面积时，要注意单元结点的排列顺序， 当三个结点i，j，m取逆时针顺序时， ； 当三个结点i，j，m取顺时针顺序时， 。 单元类型与分析 三角形单元形态矩阵的计算方法 三 位移函数和插值函数表示方法 二 等参单元简介 四 一 三角形单元的形态函数Ni具有明确的几何含义。 P点的坐标为(x,y)，则矩阵 为三角PIJ面积的两倍，因此形态函数， Ni为三角形面积之比。同样， 例题2：如图所示等腰三角形单元，求其形态矩阵[N]。 解：由 在公式中轮换下标可以计算得 三角形积为 形态函数为 根据形态函数的几何意义，可以直接得到上面的结果。 形态矩阵为 三角形面积的计算公式可得， 如果把三个结点按顺时针方向排列， 即i（a，0），j（0，0），m（0，a） 位移函数和插值函数表示方法 单元类型与分析 等参单元简介 二 三角形单元形态矩阵的计算方法 三 四 一 x y 等参单元使得我们可以处理非矩形、具有曲边的单元，能够使有限元模型更适应形状复杂的物体，等参元技术可应用于平面，空间，板，壳所有单元。 基本概念： 基本单元——形状规则，定义在参考坐标系中 实际单元——形状可能不规则 坐标插值 位移插值 等参单元的要点是：将任意物理空间坐标的单元变换到参考坐标系中，采用的变换函数和描写场函数的单元形函数相同。 如果位移映射的形函数阶数高，则称次参单元，否则称超参单元。 单元内任一点的坐标 x ，y 采用节点的坐标进行插值 例如，对四边形等参单元，单元内任一点的坐标写成矩阵形式为： (2) 八节点四边形单元 1 2 3 4 5 6 7 8 (3) 八节点六边形单元 (1) 四节点四边形等参元 小 结 单元位移函数和插值函数的概念 插值函数的性质和几何含义 单元位移函数和插值函数表示方法 三角形单元形态矩阵的计算方法 等参单元概念及表示方法 试证明，在三角形单元中的任意一点上，所有形态函数之和等于1。 三角形单元的结点坐标如图所示，设单元中一点A的坐标为（0.5，0.2）。已知三角形三结点单元的i结点位移为（2.0,1.0），j结点位移为(2.1,1.1)，m结点位移为(2.15,1.05)。1）写出单元的位移函数；2）求A点的位移分量。 m(0,0) j(1,1) i(1,0) x y 课后思考题 第3章 单元位移模式及插值函数的构造 一、单元类型与分