材料成形计算机模拟第二章

1、材料成形计算机模拟第二章材料成形计算机模拟-第二章主要内容 2-1 弹性力学的基本方程弹性力学的基本方程 2-2 弹性力学平面问题的有限元列式弹性力学平面问题的有限元列式 4节点矩形单元节点矩形单元 2-3 轴对称问题轴对称问题 2-4 三维问题三维问题单单 元元 分分 析析 的的 内内 容容节点位移(1)单元内部各点位移单元应变单元应力(2)(3)节点力(4)位移协调模式几何方程物理方程平衡方程边界条件单元分析单元刚度矩阵单元刚度矩阵 2-2 弹性力学平面问题的弹性力学平面问题的有限元列式有限元列式 八、八、4节点矩形单元节点矩形单元 矩形单元采用比常应变三角形单元次数矩形单元采用比常应变三

2、角形单元次数更高的位移模式，故可以更好地反映弹更高的位移模式，故可以更好地反映弹性体中的位移状态和应力状态。性体中的位移状态和应力状态。 xyyxvxyyxu87654321该单元的位移模式该单元的位移模式这种单元的位移模式是完备的和协调的，满足解的收敛条件，因此4节点矩形单元是协调单元。将4个节点的坐标和位移代入上式可求的18引入一个局部坐标系、，局部坐标的原点取在矩形的形心， 和轴分别与整体坐标轴x和y平行，其坐标变换的关系为00 xxayyb0123401234()/2()/2()/2()/2xxxxxyyyyy式中1 12233441 1223 344uN uN uN uN uvN v

3、N vN vN v用节点位移表示的单元位移模式可在此局部坐标系中表达为1234(1)(1)/4(1)(1)/4(1)(1)/4(1)(1)/4NNNN其中4个角点的自然坐标分别是（-1，-1），（1，-1），（1，1），（-1，1） 0000,(1,2,3,4)1(1)(1)(1,2,3,4)4iiiiNieeeeCCuS uB单元应力为单元应变为eeLNuBu 4节点矩形单元的位移模式比常应变三角节点矩形单元的位移模式比常应变三角形单元中采用的线性模式增添了形单元中采用的线性模式增添了xy项，项，所以矩形单元内的应变分量、应力分量所以矩形单元内的应变分量、应力分量都不是常量，而是沿着都不是常

4、量，而是沿着x及及y方向呈线性方向呈线性变化。因此，在弹性体中采用相同数目变化。因此，在弹性体中采用相同数目节点时，矩形单元的精度要比常应变三节点时，矩形单元的精度要比常应变三角形单元高。角形单元高。 但是，矩形单元亦有明显的缺点，一是不能适应斜交边界和曲线边界，二是不便于对不同部位采用不同大小的单元，因此直接应用受到限制。 2-3 轴对称问题 许多机械零件和结构的几何形状、约束条件以及作用的载荷都对称于某一对称轴，在这种条件下的物体中的位移、应变和应力也对称于此轴。这种问题称为轴对称问题。 应力、应变都与无关，仅是坐标r和z的函数。沿方向的位移为0，因此轴对称问题可作为二维问题处理。 对轴对

5、称问题进行计算时，只需取出一个截面进行网格划分和分析，但应注意到单元是圆环状的，所有的结点载荷都应理解为作用在单元结点所在的圆周上。 本节主要以3节点三角形轴对称环状单元为例进行讨论。它具有与平面三角形单元同样的特点。 在轴对称问题中，通常采用圆柱坐标（r，z）。以对称轴作为z轴，所有应力、应变和位移都与方向无关，只是r和z的函数。任一点的位移只有两个方向的分量，即沿r方向的径向位移u和沿z方向的位移w。由于轴对称，方向的位移v等于零，因此轴对称问题是二维问题。 离散轴对称物体时，采用的单元是一些圆环，各单元在rz平面内形成网格。 一、3节点三角形轴对称单元的插值函数及应力应变矩阵（一）形函数

6、（一）形函数 子午面内的环单元与前面讨论过的平面问题3节点三角形是一样的，它们的形函数也完全一样。形函数是用单元位移分量来描述位移函数的插值函数。形函数是用单元位移分量来描述位移函数的插值函数。i ijjm mi ijjm muNuNuN uvNvN vN v1()2iiiiNabx cyA(, , )i j m其中 （二）单元应变和应力（二）单元应变和应力 将位移插值函数代入轴对称问题的几何方程，得到单元应变：(2)()()TTeerzrzrzijmuuuwuBuBBBurzrzr其中001( , ,)02iiiiiibci j mfAcbB( , ,)iiiiac zfbi j mrr由上

7、式可见，应变分量r、z、rz都是常量，环向应变不是常量。单元应力： ()()TeeeeerzrzijmCC BuSuSSSu101110(1)1(1)(1 2 )101 22(1)eEC对称轴对称问题的弹性矩阵为由于不是常量，所以单元中除切应力rz外其它应力分量也不是常量。 二、3节点轴对称单元的单元刚度矩阵和等效节点载荷eeepTeeTTTcVVSB C BdVuN bdVN pdSN PeTeeVB C BdVKeTebVN bdVP令epTesSN pdSPTeccN PPeeeebscPPPP将轴对称问题的N、B、Ce和dV=rddrdz代入上式，即可求得轴对称单元的单元刚度矩阵和等效

8、节点载荷 为了简化计算和消除在轴对称上r0对积分所带来的麻烦，将积分式中的自变量r、z用单元截面形心处的坐标来近似。)(31)(31mjimjizzzzzrrrrr),(mjirzcbraffiiiii这样(3-56)就近似为 作了这样的近似后，应变矩阵B和应力矩阵S都成了常量阵，简化了计算。 在轴对称问题中，对于单元刚度矩阵和等效节点载荷向量采用上述近似积分方法，就位移和应力而言，其精度是能够满足工程计算要求的。 2-4 三维问题 一、常应变四面体单元 二、六面体单元 一、常应变四面体单元 （一）位移函数（一）位移函数 四面体单元以4个角点为节点。每个节点有3个自由度，一个单元共有12个自由

9、度。()()TeeijmluuvwNuININININ u其中I为三阶单位矩阵。 ()eTiiijjjmmmllluuvwuvwuvwuvw1()6iiiiiNab xc yd zV( , , )i j m l123456789101112uxyzvxyzwxyz位移函数为整理后可得111( , , )111111jjjimmmllljjimmlljjimmlljjimmllxyzaxyzxyzyzbyzyzi j m lyzcyzyzyzdyzyz V是四面体ijml的体积。为了使四面体的体积不为负值，单元节点编号 必须依照一定的顺序。, ,i j m l111161iiijjjmmmlll

10、xyzxyzVxyzxyz设P(x，y，z)为四面体中任一点，记四面体jmlP的体积为Vi，则 111161jjjimmmlllxyzxyzVxyzxyz上式按第4列展开得1()6iiiiiVab xc yd z 定义四面体单元中节点i的体积坐标为 iiVLV1()6iiiiiiiVLab xc yd zNVV与3节点三角形单元的面积坐标相对应,这里形函数,ijmlN NNN即是四面体单元的体积坐标。ijjjjijjjizyxNzyxL),(),(由体积坐标定义可知又因1jimlijmlVVVVLLLLVVVVijmlVVVVV由于位移函数是线性的，相邻单元交界面上的位移由该界面上三个节点位移

11、所决定，因此是连续的，所以常应变四边形单元是协调元。 （二）应变矩阵和应力矩阵（二）应变矩阵和应力矩阵 将三维问题的应变分量写成向量的形式为将三维问题的应变分量写成向量的形式为(222)()()TTxyzxyyzzxeeijmluvwuvvwwuxyzyxzyxzBuBBBB u式中00010006000rrrTrrrrrrrbcdBcbdVdcb应变矩阵Br的元素均为常数,故四面体单元是一种常应变单元。 单元应力为eeeeCC BuSu平面应变问题和轴对称问题的弹性矩阵与三维问题弹性矩阵的对应元素是相同。100011100011000(1)1 200(1)(1 2 )2(1)1 202(1)1 22(1)eEC对称 （三）单元刚度矩阵和等效节点载荷（三）单元刚度矩阵和等效节点载荷 四面体单元的刚度矩阵和等效节点载荷可利用式（3-31）求得，对于三维问题，dV=dxdydz。其中应变矩阵B B为常数矩阵，可提到积分号外。 四面体单元优点：1）适应多种复杂边界形状；2）容易实现网格密度的变化；3）有利于对不规则三维空间进行全自动网格剖分。因此得到广泛应用。 缺点：四面体的拼合较复杂，划分时容易出错，不容易直观地理解。二、六面体单元二、六面体单元 8节点六面体单元，每个节点有3个自由度，一个单元共有24个自由度。插值多项式中包括如下各项： 以