第4章 轴对称问题和空间问题有限元法课件

1轴对称及空间问题有限元法第四章轴对称和空间问题有限元法上海工程技术大学机械工程学院工程力学部1轴对称及空间问题有限元法第四章轴对称和空间问题有限元法上2第一节轴对称问题有限元法第二节空间问题有限元法2第一节轴对称问题有限元法第二节空间问题有限元法3一.轴对称问题的定义工程中有一类结构，它们的几何形状、约束条件及作用的荷载都对称于某一固定轴(可视为子午面内平面物体绕轴旋转一周的结果)，其力学分析称为轴对称问题。典型例子为烟囱、储液罐等受恒载作用。1.离散化

由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。2.应力和应变

对轴对称问题进行分析一般取柱坐标系，对称轴为Z轴，径向为r轴，环向为θ轴。xyz对称面上任一点p只会在该对称面上发生位移，即所有的应力、应变和位移只是z和r的函数，而与坐标θ无关。那么轴对称问题就可转化为二维平面问题来进行研究。但因与平面问题有区别，常称二维半问题。§4-1

轴对称问题有限元法3一.轴对称问题的定义工程中有一类结构，它们的几何4

如图所示的受均布内压作用的长圆筒，通过z轴的一个纵截面就是对称面。由于对称性，轴对问题共有4个应力分量：其中表示沿半径方向的正应力，称为径向应力；表示沿方向的正应力，称为环向应力或切向应力；表示沿ｚ方向的正应力，称为轴向应力；表示在圆柱面上沿ｚ方向作用的剪应力。4如图所示的受均布内压作用的长圆筒，通过z轴的一个纵截面5

同样，轴对称问题共有4个应变分量：其中表示沿半径方向的正应变，称为径向正应变；表示沿方向的正应变，称为环向正应变或切向正应变；表示沿ｚ方向的正应变，称为轴向正应变；表示沿ｒ和ｚ方向的剪应变。

在轴对称问题中，弹性体内任意一点上，不存在切向位移，只存在径向位移u

和轴向位移w

，两个位移分量表示为，xyz基本方程1.平衡方程5同样，轴对称问题共有4个应变分量：其中表示沿半径62.几何方程62.几何方程72.物理方程72.物理方程8

轴对称单元的特点(与平面三角形单元的区别)轴对称单元为圆环体，单元与单元间为节圆相连接；节点力与节点载荷是施加于节圆上的均布力；单元边界是一回转面；应变分量中出现了，即应变不是常量；且应变矩阵在r→0时，存在奇异点，需特殊处理，通常用该单元的形心坐标替代节点坐标。轴对称结构8轴对称单元的特点(与平面三角形单元的区别)轴对称结构9轴对称问题的有限元法1.离散化

由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。2.单元分析参照平面问题的三角形单元位移函数，轴对称问题的三结点三角形单元位移函数取为，9轴对称问题的有限元法1.离散化由于可视为子午面内平面10其中：形函数：10其中：形函数：11用矩阵表示的单元位移为：单元应变：

将单元位移函数代入几何方程得：11用矩阵表示的单元位移为：单元应变：其中，用几何矩阵表示单元的应变：

由于在是坐标r、z

的函数，分量在单元中不为常量，其它三个应变分量在单元中仍为常量。几何矩阵不再是常数，轴对称三角形单元内的应变也不全为常量。其中，用几何矩阵表示单元的应变：由于在是坐标r、z13单元应力：由弹性矩阵D

和几何矩阵B

可以得到应力矩阵S，并计算出单元内的应力分量，其中：式中：

由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标代替B

矩阵中的变量，将单元中的r和z近似地当作常量，并且分别等于。13单元应力：其中：式中：由于几何矩阵中的元素不是常量，14经过简化，就可以把各个单元近似地当作常应变单元。14经过简化，就可以把各个单元近似地当作常应变单元。153.单元刚度矩阵有了单元应力场和应变场，可以利用虚位移原理或最小势能原理建立单元刚度矩阵单元刚度矩阵的分块矩阵为

由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标代替B

矩阵中的变量。

实践证明采用近似积分也能达到一定的精度，具体对于三角形环单元用形心处坐标代替应变矩阵中的坐标变量。153.单元刚度矩阵有了单元应力场和应变场，可以利用16应变矩阵变成：其中：单元刚度矩阵的近似表达式为：单元刚度矩阵的分块矩阵近似表达式为：16应变矩阵变成：其中：单元刚度矩阵的近似表达式为：单元刚度174.总刚度矩阵集成

求出了每一个三角形单元的刚度矩阵后，按照平面问题介绍的总刚矩阵的集成方法，就可以得到结构的总刚矩阵。5.等效节点载荷计算

计算轴对称问题的等效节点载荷与平面问题的有所不同，因为周对称结构的子午面上的一个节点是关于对称轴中心对称的圆环，故当计算集中力、表面力和体积力时，应在整个环上积分。（1）集中力移置集中力为

--集中力作用点的径向坐标。174.总刚度矩阵集成求出了每一个三角形单元的刚度矩阵后18（2）体积力移置若体积力为重，则单位体积的力为ρ为密度度如果单元离开对称轴较远（）可认为将1/3的自重移置到每个节点上。单元自重移置到节点i，j，m上的等效节点载荷为(i，j，m)面积坐标积分式（在三角形全面积上）面积坐标积分式（在三角形单元边i-j上）18（2）体积力移置若体积力为重，则单位体积的力为ρ为密度度19若体积力为惯性离心力，则单位体积的力为单元离心力移置到节点i，j，m上的等效节点力(i，j，m)设rz平面上单元ijm的ij边有线性分布的径向表面力，在节点i,j的集度分别为，，ij边的边长l。此时有，，，于是节点i的等效节点载荷为(3)分布面力移置19若体积力为惯性离心力，则单位体积的力为单元离心力移置到节20若，则得6.约束处理和求解线性方程组（总体刚度方程）K称为总体刚度矩阵

U称为位移向量

F称为载荷向量

高斯消去法、波前法等20若，则得6.约束处理和求解线性方程组（总体刚度方程）21二.空间问题的有限元法如果物体的形状、尺寸和边界条件不具备某种特殊性，不能简化为平面问题和轴对称问题处理，则必须作为空间问题来分析。

与平面问题和轴对称问题相比，空间问题的计算要复杂的多。

单元为块体形状。常用单元：四面体单元、长方体单元、直边六面体单元、曲边六面体单元、轴对称单元。（2）结点位移3个分量。（3）基本方程比平面问题多。3个平衡方程，6个几何方程，6个物理方程。

§4-2

空间问题有限元法21二.空间问题的有限元法如果物体的形状、尺寸224结点四面体单元：是空间问题最简单的单元，也是常应变、常应力单元，可以类似平面问题三结点三角形单元进行分析。8结点长方体单元：可以类似平面四结点矩形单元进行分析。8结点直边六面体单元：可以类似平面四结点任意四边形等参元分析。20结点曲边六面体单元：等参单元,可以类似平面八结点曲边四边形等参元进行分析。轴对称单元：一平面单元绕一对称轴旋转形成的空间问题。只需在rz平面划分网格，就像平面问题xy平面中的网格一样，这样这类空间问题可以得到简化。（环向位移等于零）224结点四面体单元：是空间问题最简单的单元，也是常应变、常231.四结点四面体单元分析1）位移函数

单元结点位移向量单元位移函数xzyijmn231.四结点四面体单元分析1）位移函数单元位移函数xzy24将节点坐标代入位移函数

由此可解出代定常数α1~α4再代回到式单元位移函数的第1式，可得位移函数u:

24将节点坐标代入位移函数由此可解出代定常数α1~25编号约定：当沿i,j,m的方向转动时，n在大拇指所指的方向25编号约定：当沿i,j,m的方向转动时，n在大拇指所指的方26采用同样的方法，可得单元位移：

26采用同样的方法，可得单元位移：272）单元应变将单元中位移代入上式272）单元应变将单元中位移代入上式283）单元应力弹性矩阵[D]:283）单元应力弹性矩阵[D]:29代入单元应变计算公式，整理后:其中[S]为应力矩阵，且：常量29代入单元应变计算公式，整理后:其中[S]为应力矩阵，且：304)单元刚度矩阵分块矩阵的形式

304)单元刚度矩阵分块矩阵的形式31式中子矩阵[krs]为3×3的矩阵:5)等效结点荷载

体积力与表面力的计算公式与平面三角形单元公式相似，可以采用静力等效原则简化计算。31式中子矩阵[krs]为3×3的矩阵:5)等效结点荷载2六面体类单元的形函数

1）八结点单元12345678ξηζ类似平面问题矩形线性单元，由试凑法可建立形函数如下：2）二十结点单元和平面问题一样，基于试凑法，可以根据上述八结点低阶单元形函数构造各顶点形函数。ξηζ1234567891011121417202六面体类单元的形函数1）八结点单元12345678ξη33轴对称及空间问题有限元法第四章轴对称和空间问题有限元法上海工程技术大学机械工程学院工程力学部1轴对称及空间问题有限元法第四章轴对称和空间问题有限元法上34第一节轴对称问题有限元法第二节空间问题有限元法2第一节轴对称问题有限元法第二节空间问题有限元法35一.轴对称问题的定义工程中有一类结构，它们的几何形状、约束条件及作用的荷载都对称于某一固定轴(可视为子午面内平面物体绕轴旋转一周的结果)，其力学分析称为轴对称问题。典型例子为烟囱、储液罐等受恒载作用。1.离散化

由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。2.应力和应变

对轴对称问题进行分析一般取柱坐标系，对称轴为Z轴，径向为r轴，环向为θ轴。xyz对称面上任一点p只会在该对称面上发生位移，即所有的应力、应变和位移只是z和r的函数，而与坐标θ无关。那么轴对称问题就可转化为二维平面问题来进行研究。但因与平面问题有区别，常称二维半问题。§4-1

轴对称问题有限元法3一.轴对称问题的定义工程中有一类结构，它们的几何36

如图所示的受均布内压作用的长圆筒，通过z轴的一个纵截面就是对称面。由于对称性，轴对问题共有4个应力分量：其中表示沿半径方向的正应力，称为径向应力；表示沿方向的正应力，称为环向应力或切向应力；表示沿ｚ方向的正应力，称为轴向应力；表示在圆柱面上沿ｚ方向作用的剪应力。4如图所示的受均布内压作用的长圆筒，通过z轴的一个纵截面37

同样，轴对称问题共有4个应变分量：其中表示沿半径方向的正应变，称为径向正应变；表示沿方向的正应变，称为环向正应变或切向正应变；表示沿ｚ方向的正应变，称为轴向正应变；表示沿ｒ和ｚ方向的剪应变。

在轴对称问题中，弹性体内任意一点上，不存在切向位移，只存在径向位移u

和轴向位移w

，两个位移分量表示为，xyz基本方程1.平衡方程5同样，轴对称问题共有4个应变分量：其中表示沿半径382.几何方程62.几何方程392.物理方程72.物理方程40

轴对称单元的特点(与平面三角形单元的区别)轴对称单元为圆环体，单元与单元间为节圆相连接；节点力与节点载荷是施加于节圆上的均布力；单元边界是一回转面；应变分量中出现了，即应变不是常量；且应变矩阵在r→0时，存在奇异点，需特殊处理，通常用该单元的形心坐标替代节点坐标。轴对称结构8轴对称单元的特点(与平面三角形单元的区别)轴对称结构41轴对称问题的有限元法1.离散化

由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。2.单元分析参照平面问题的三角形单元位移函数，轴对称问题的三结点三角形单元位移函数取为，9轴对称问题的有限元法1.离散化由于可视为子午面内平面42其中：形函数：10其中：形函数：43用矩阵表示的单元位移为：单元应变：

将单元位移函数代入几何方程得：11用矩阵表示的单元位移为：单元应变：其中，用几何矩阵表示单元的应变：

由于在是坐标r、z

的函数，分量在单元中不为常量，其它三个应变分量在单元中仍为常量。几何矩阵不再是常数，轴对称三角形单元内的应变也不全为常量。其中，用几何矩阵表示单元的应变：由于在是坐标r、z45单元应力：由弹性矩阵D

和几何矩阵B

可以得到应力矩阵S，并计算出单元内的应力分量，其中：式中：

由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标代替B

矩阵中的变量，将单元中的r和z近似地当作常量，并且分别等于。13单元应力：其中：式中：由于几何矩阵中的元素不是常量，46经过简化，就可以把各个单元近似地当作常应变单元。14经过简化，就可以把各个单元近似地当作常应变单元。473.单元刚度矩阵有了单元应力场和应变场，可以利用虚位移原理或最小势能原理建立单元刚度矩阵单元刚度矩阵的分块矩阵为

由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标代替B

矩阵中的变量。

实践证明采用近似积分也能达到一定的精度，具体对于三角形环单元用形心处坐标代替应变矩阵中的坐标变量。153.单元刚度矩阵有了单元应力场和应变场，可以利用48应变矩阵变成：其中：单元刚度矩阵的近似表达式为：单元刚度矩阵的分块矩阵近似表达式为：16应变矩阵变成：其中：单元刚度矩阵的近似表达式为：单元刚度494.总刚度矩阵集成

求出了每一个三角形单元的刚度矩阵后，按照平面问题介绍的总刚矩阵的集成方法，就可以得到结构的总刚矩阵。5.等效节点载荷计算

计算轴对称问题的等效节点载荷与平面问题的有所不同，因为周对称结构的子午面上的一个节点是关于对称轴中心对称的圆环，故当计算集中力、表面力和体积力时，应在整个环上积分。（1）集中力移置集中力为

--集中力作用点的径向坐标。174.总刚度矩阵集成求出了每一个三角形单元的刚度矩阵后50（2）体积力移置若体积力为重，则单位体积的力为ρ为密度度如果单元离开对称轴较远（）可认为将1/3的自重移置到每个节点上。单元自重移置到节点i，j，m上的等效节点载荷为(i，j，m)面积坐标积分式（在三角形全面积上）面积坐标积分式（在三角形单元边i-j上）18（2）体积力移置若体积力为重，则单位体积的力为ρ为密度度51若体积力为惯性离心力，则单位体积的力为单元离心力移置到节点i，j，m上的等效节点力(i，j，m)设rz平面上单元ijm的ij边有线性分布的径向表面力，在节点i,j的集度分别为，，ij边的边长l。此时有，，，于是节点i的等效节点载荷为(3)分布面力移置19若体积力为惯性离心力，则单位体积的力为单元离心力移置到节52若，则得6.约束处理和求解线性方程组（总体刚度方程）K称为总体刚度矩阵

U称为位移向量

F称为载荷向量

高斯消去法、波前法等20若，则得6.约束处理和求解线性方程组（总体刚度方程）53二.空间问题的有限元法如果物体的形状、尺寸和边界条件不具备某种特殊性，不能简化为平面问题和轴对称问题处理，则必须作为空间问题来分析。

与平面问题和轴对称问题相比，空间问题的计算要复杂的多。

单元为块体形状。常用单元：四面体单元、长方体单元、直边六面体单元、曲边六面体单元、轴对称单元。（2）结点位移3个分量。（3）基本方程比平面问题多。3个平衡方程，6个几何方程，6个物理方程。

§4-2

空间问题有限元法21二.空间问题的有限元法如果物体的形状、尺寸544结点四面体单元：是空间问题最简单的单元，也是常应变、常应力单元，可以类似平面问题三结点三角形单元进行分析。8结点长方体单元：可以类似平面四结点矩形单元进行分析。8结点直边六面体单元：可以类似平面四结点任意四边形等参元分析。20结点曲边六面体单元：等参单元,可以类似平面八结点曲边四边形等参元进行分析。轴对称单元：一平面单元绕一对称轴旋转形成的空间问题。只需在rz平面划分网格，就像平面