第5讲、空间问题有限元法.pptVIP

四面体常应变单元 位移函数 应变矩阵 单元刚度矩阵 载荷分配 四面体的组合 * * 空间问题(三维)有限元分析 空间三维应力状态 一般的实际物体都是立体的,弹性体受力作用后,其内部各点将沿,X,Y,Z三个坐标的方向发生位移,是三维问题.如各点沿X,Y,Z方向的位移以μ,ν,ω表示,这些位移一般应为各点坐标的函数,即: 弹性体一般变形情况下,有三个方向的线应变 和三对剪应变 由弹性力学可知,应变与位移间的几何关系是: 上式可视为平面问题的推广，而平面问题则为三维问题得特殊情况。 三维弹性体得应变一般有上述六项分量，这些分量可用一列阵{ε}表示，则式（1）可用矩阵形式表示为 弹性体受力作用,内部任一点的应力状态也是三维的,有三项正应力: σx σy σz三对剪切应力: τxy = τyx , τyz = τzy , τzx = τxz 在线弹性范围内,应力应变间的物理关系可用矩阵形式表示为: 对于各向同性的弹性体,在三维应力状态下,弹性系数矩阵[D]的一般形式为: 最简单的空间单元一四面体单元如图所示,i , j , k , m为四个结点,为使单元体积不出现负值,结点的编号按下规定:在右手坐标系中,当右手螺旋按i—j—k转向时,拇指指向m. 单元变形时,各结点都有沿x ,y ,z的三项位移,单元有四个结点,共有12项结点位移,合起来以列阵表示: 单元变形时,单元内各点也有沿x, y, z方向的位移 μ,ν,ω ,一般应为坐标x, y, z,的函数.对于这种简单的四面体单元,其内部位移可假设为坐标的线性函数,全单元有12个自由度,每个单元有三个位移函数,每个位移函数取四项为: 式中a1,a5,a9反映刚体平移,a2,a7,a12,反映常量正应变,其余各项反映绕三坐标轴旋转的刚体转动及常量剪应变. 　　上式含12个待定常数,可以由单元的12项结点位移决定.将四个结点的坐标值代入式(5)中的 μ式:在i ,j ,k ,m四个结点,分别有: 由上式求出a1,a2,a3,a4,代入(5)式可得: 同理 将上面三式用矩阵表示为: 其中[N]为形函数矩阵,可表示为: 式中[I]为三阶单位矩阵,而各结点的形函数可按下式计算得到,即 V为四面体单元的体积,其他 将式(7)代入几何关系式(2),经过微分运算,可以得到单元内应变为: 其中应变矩阵(几何矩阵)是形状函数矩阵[N]经微分算子矩阵[L]作用所得的结果,[B]中任一个子矩阵[Bi]的显式应为: 式中的V及bi, ci, di等按式( 9)决定.可见,这里〔BiJ的每项元素都是由结点坐标决定的常数;因而简单四面体单元内,各点的应变都是一样的,这是一种常应变单元.,因而由位移一阶导数组成的应变,自然就是常值了. 将(4)及(10)代入 即可计算得单元刚度矩阵.由于简单四面体单元为常应变单元.故积分有: 其中Ve为单元体积,按结点分块,此单元刚度矩阵表示为 其中任一子矩阵为: 三维弹性体内如受有均布的体积力(如重力)作用,对于这种简单的四面体单元,可以逐个单元计算出其整个单元的全部体积力,再平均分配到四个结点上,即每个结点分配到1/4的单元体积力.如果单元的某个表面作用有均布的面积力,也可将此面上的全部面积力平均分配到相连的三个结点上,即每个结点分配到三角形上面积力总和的1/3,如果体积力、面积力不是均布的,则不应平均分配,而应按下式计算 其中{PV}e , {PS}e 为e单元内分布体积力和分布面积力分配到单元结点得载荷,[N]为形状函数矩阵，{q}和{p}分别为单元体积力和单元面积力，Ve，Se 则为受有分布力得单元体积和面积。 用四面体单元分析三维结构,单元划分是很灵活的,可以逼近较复杂的几何形状,可以按需要在局部区域随意加密或变稀网格,为了减小计算误差,应尽可能使四面体的六个棱边有相近的长度,而接近正四面体. 　　实际上的三维结构,除在边缘,棱角处为适应边界几何形状,可能需要划分多个四面体单元之外,其余大部分区域都可以划分成一般的六面体或三棱体,采用较为复杂的单元.这样,三维结构可以作为这种复杂单元的组合.也可以按图所示,把六面体,三棱体按一定的规则分解为数个四面体单元的组合.