形函数的定义
形函数是定义于单元内部的、坐标的连续函数,它应满足下列要求:
(1)在结点i上,Ni=1;在其它结点上,Ni=0;
(2)能保证用它定义的未知量(u,v或x,y)在相邻单元之间的连续性;
(3)应包含任意线性项,以便用它定义的单元位移可满足常应变条件;
(4)应满足下列等式:∑Ni=1,以便用它定义的单元位移可反映刚体移动;
形函数的阶次越高,单元形状越复杂,单元适应能力也越强,求解应力问题时所需单元数量越少,因而,平衡方程组的阶较低,求解过程所花费的时间也越少。
一维形函数举例:
二维形函数举例:
三维形函数举例:
坐标变化
上述举例的单元的边界都是直线或平面,使用的是局部坐标系(也称为自然坐标系),此类单元称为母单元。形函数是在母单元中用自然坐标定义的。
母单元几何形状简单,便于运算,但难以适应实际工程中出现的各种复杂形状的结构。为了贴合比较复杂的结构外形,可以利用坐标变换,建立在整体坐标系下具有复杂外形的曲边单元,此类单元称为子单元。
经过坐标变换,单元具有双重特性:
(1)子单元的几何特征、载荷等,都来自实际结构,充分反映了实际情况;
(2)大量计算工作是在母单元中进行的,由于它形状简单且规则,运算方便。
为了进行坐标变换,必须在自然坐标系(ξ,η,ζ)和整体坐标系(x,y,z)之间建立一一对应的映射关系,这种对应关系也是通过形函数建立起来的。
在整体坐标系(x,y,z)中,子单元内任一点的坐标用形函数表示:
回忆一下,形函数的另一个作用,反映单元的位移形态:
单元坐标变换公式和单元位移函数中都用到了形函数,它们可以是自然坐标下的一次、二次、三次、甚至更高次的函数。如果单元坐标变换公式和单元位移函数中所用的形函数的阶次相等,那么用于规定单元形状的结点数和用于规定单元位移的结点数应相等,这种单元称为等参数单元。
在等参数单元中,坐标变换和位移函数一般采用相同的结点。
一维单元的坐标变换,如下图所示:
通过形函数,用结点位移定义单元内任一点的位移:
单元应变为:
由于形函数是由局部坐标给出的,根据偏微分法则,可知:
在整体坐标系中,子单元内任一点的坐标用形函数表示:
所以:
由此可求得矩阵[B]:
单元刚度矩阵的普遍公式为:
又因为:
所以:
二维单元的坐标变换,如下图所示:
通过形函数,用结点位移定义单元内任一点的位移:
单元应变为:
由于形函数是由局部坐标给出的,根据偏微分法则,可知:
在整体坐标系中,子单元内任一点的坐标用形函数表示:
由此,可解得矩阵[Bi]:
由坐标变换公式,可直接解得雅克比矩阵:
对其求逆后,得到形函数在整体坐标中的导数如下:
即可得到矩阵[Bi]。
三维单元的坐标变换,如下图所示:
通过形函数,用结点位移定义单元内任一点的位移:
单元应变为:
由于形函数是由局部坐标给出的,根据偏微分法则,可知:
集合起来,得到:
由于形函数是由局部坐标给出的,上式左端可根据Ni(ξ,η,ζ)直接求得。同时,由坐标变换公式求解雅克比矩阵:
根据坐标变换公式,先计算出雅克比矩阵中的矩阵[j],这是各单元共用的。再把结点坐标代入上式,求得雅克比矩阵,并求逆后,得到形函数在整体坐标中的导数如下:
从而计算得到矩阵[Bi],进一步计算得到矩阵[B]。
令:
矩阵分块后,子矩阵为:
又由于:
可得到:
如果子单元的几何形状较为简单,如二维的矩形、三维的正六面体,可直接求出上式右端的积分并得到用显式表示的刚度矩阵。但一般情况下,子单元形状较为复杂,刚度矩阵难以用显式表示,必须依靠数值积分计算[krs]的值。
雅克比矩阵对单元应变有重要的影响。以简单的平面矩形单元为例:
该单元的雅克比矩阵是个常数。
得到矩阵[Bi]:
显然,单元内应变的变化规律取决于形函数的阶次。