有限元-第4讲-轴对称问题与空间问题有限元法ppt课件.ppt

BIPT BIPT第4章 轴对称问题与空间问题有限元方法张 洪 伟 BIPT BIPT轴对称问题的有限元法轴对称问题的有限元法空间问题的有限元法空间问题的有限元法12内容提要内容提要 BIPT BIPT3一.轴对称问题的定义 工程中有一类结构，它们的几何形状、约束条件及作用的荷载都对称于某一固定轴(可视为子午面内平面物体绕轴旋转一周的结果)，其力学分析称为轴对称问题。1.离散化离散化 由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。2.应力和应变应力和应变 对轴对称问题进行分析一般取柱坐标系，对称轴为Z轴，径向为r 轴，环向为轴。xyz 对称面上任一点p只在该对称面上发生位移，即所有的应力、应变和位移只是z和r的函数，而与坐标无关。那么轴对称问题就可转化为二维平面问题来进行研究。但因与平面问题有区别，常称二维半问题。轴对称问题的有限元法轴对称问题的有限元法 BIPT BIPT如图所示的受均布内压作用的长圆筒，通过z轴的一个纵截面就是对称面。由于对称性，轴对问题共有4个应力分量：其中 表示沿半径方向的正应力，称为径向应力；表示沿 方向的正应力，称为环向应力或切向应力；表示沿方向的正应力，称为轴向应力；表示在圆柱面上沿方向作用的剪应力。BIPT BIPT同样，轴对称问题共有4个应变分量：其中表示沿半径方向的正应变，称为径向正应变；表示沿方向的正应变，称为环向正应变或切向正应变；表示沿方向的正应变，称为轴向正应变；表示沿和方向的剪应变。在轴对称问题中，弹性体内任意一点上，不存在切向位移，只存在径向位移 u 和轴向位移 w，两个位移分量表示为:xyz基本方程1.平衡方程 BIPT BIPT2.几何方程 BIPT BIPT3.物理方程 BIPT BIPT 轴对称单元的特点(与平面三角形单元的区别)l轴对称单元为圆环体，单元与单元间为节圆相连接；l节点力与节点载荷是施加于节圆上的均布力；l单元边界是一回转面；l应变分量 中出现了 ，即应变不是常量；且应变矩阵在r=0时，存在奇异点，需特殊处理，通常用该单元的形心坐标替代节点坐标。轴对称结构 BIPT BIPT9轴对称问题的有限元法1.离散化 由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。2.单元分析 参照平面问题的三角形单元位移函数，轴对称问题的三结点三角形单元位移函数取为，BIPT BIPT其中：形函数：BIPT BIPT用矩阵表示的单元位移为：单元应变：将单元位移函数带入几何方程得：BIPT BIPT其中，用几何矩阵表示单元的应变：由于是坐标 r、z 的函数，分量在单元中不为常量，其它三个应变分量在单元中仍为常量。应变矩阵 不再是常数，轴对称三角形单元内的应变也不全为常量。BIPT BIPT单元应力：由弹性矩阵 D 和应变矩阵 B 可以得到应力矩阵 S，并计算出单元内的应力分量：其中：式中：由于应变矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标 代替 B 矩阵中的变量，将单元中的r和z近似地当作常量，并且分别等于 。BIPT BIPT经过简化，就可以把各个单元近似地当作常应变单元。BIPT BIPT3.单元刚度矩阵 有了单元应力场和应变场，可以利用虚位移原理或最小势能原理建立单元刚度矩阵，注意体积积分是沿单元的整个圆环求体积积分单元刚度矩阵的分块矩阵为，由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标 代替 B 矩阵中的变量。实践证明采用近似积分也能达到一定的精度，具体对于三角形环单元用形心处坐标代替应变矩阵中的坐标变量。BIPT BIPT应变矩阵变成：其中：单元刚度矩阵的近似表达式为：单元刚度矩阵的分块矩阵近似表达式为：BIPT BIPT4.总刚度矩阵集成求出了每一个三角形单元的刚度矩阵后，按照平面问题介绍的总刚矩阵的集成方法，就可以得到结构的总刚矩阵。5.等效节点载荷计算计算轴对称问题的等效节点载荷与平面问题的有所不同，因为周对称结构的子午面上的一个节点是关于对称轴中心对称的圆环，故当计算集中力、表面力和体积力时，应在整个环上积分。（1）集中力移置 集中力为-集中力作用点的径向坐标。BIPT BIPT（2）体积力移置若体积力为自重，则单位体积的力为为密度度如果单元离开对称轴较远（）可认为将1/3的自重移置到每个节点上。单元自重移置到节点i，j，m上的等效节点载荷为 (i，j，m)BIPT BIPT若体积力为惯性离心力，则单位体积的力为单元离心力移置到节点i，j，m上的等效节点力 (i，j，m)(3)分布面力移置 BIPT BIPT 如果物体的形状、尺寸和边界条件不具备某种特殊性，不能简化为平面问题和轴对称问题处理，则必须作为空间问题来分析。与平面问题和轴对称问题相比，空间问题的计算要复杂的多。（1）单元为块体形状。常用单元：四面体单元、长方体单元、直边六面体单元、曲边六面体单元、轴对称单元。（2）结点位移3个分量。（3）基本方程比平面问题多。3个平衡方程，6个几何方程，6个物理方程。空间问题的有限元法空间问题的有限元法 BIPT BIPT4结点四面体单元：是空间问题最简单的单元，也是常应变、常应力单元，可以类似平面问题三结点三角形单元进行分析。8结点长方体单元：可以类似平面四结点矩形单元进行分析。8结点直边六面体单元：可以类似平面四结点任意四边形等参元分析。20结点曲边六面体单元：等参单元,可以类似平面八结点曲边四边形等参元进行分析。轴对称单元：一平面单元绕一对称轴旋转形成的空间问题。只需在rz平面划分网格，就像平面问题xy平面中的网格一样，这样这类空间问题可以得到简化。（环向位移等于零）BIPT BIPT1.四结点四面体单元分析1）位移函数 单元结点位移向量 单元位移函数xzyijmn BIPT BIPT将节点坐标代入位移函数 由此可解出代定常数1 4再代回到式单元位移函数的第1式，可得位移函数u:BIPT BIPT编号约定：当沿i,j,m的方向转动时，n在大拇指所指的方向 BIPT BIPT采用同样的方法，可得单元位移：BIPT BIPT2）单元应变将单元中位移代入上式 BIPT BIPT3）单元应力弹性矩阵D:BIPT BIPT代入单元应变计算公式，整理后:其中S为应力矩阵，且：常量 BIPT BIPT4)单元刚度矩阵分块矩阵的形式 BIPT BIPT式中子矩阵krs为33的矩阵:5)等效结点荷载 体积力与表面力的计算公式与平面三角形单元公式 相似，可以采用静力等效原则简化计算。崇尚实践崇尚实践 知行并重知行并重