轴对称问题有限元法分析报告

轴对称问题的有限元模拟分析一、摘要：轴对称问题是弹性空间问题的一个特殊问题，这类问题的特点是物体为某一平面绕其中心轴旋转而成的回转体。

由于一般形状是轴对称物体，用弹性力学的解析方法进行应力计算，很难得到精确解，因此采用有限元法进行应力分析，在工程上十分需要，同时用有限元法得到的数值解，近似程度也比较好。

轴对称问题的有限元分析，可以将要分析的问题由三维转化为二维平面问题来解决。

先是结构离散，然后是单元分析，再进行总纲集成，再进行载荷移置，最后是约束处理和求解线性方程组。

分析完成之后用ABAQUS软件建模以及分析得出结果。

关键字：有限元法轴对称问题ABAQUS软件二、前言:1、有限元法领域介绍：有限单元法是当今工程分析中获得最广发应用的数值计算方法，由于其通用性和有效性，受到工程技术界的高度重视，伴随着计算机科学和技术的快速发展，现在已经成为计算机辅助设计和计算机辅助制造的重要组成部分。

由于有限元法是通过计算机实现的，因此有限元程序的编制以及相关软件的研发就变得尤为重要，从二十世纪五十年代以来，有限元软件的发展按目的和用途可分为专用软件和大型通用商业软件，而且软件往往集成了网络自动划分，结果分析和显示等前后处理功能，而且随着时间的发展，大型通用商业软件的功能由线性扩展到非线性，由结构扩展到非结构等等，这一系列强大功能的实现及运用都要求我们对有限元法的基础理论知识有较为清楚的认识以及对程序编写的基本能力有较好掌握。

2、研究报告目的:我们小组研究的问题是：圆柱体墩粗问题。

毛坯的材料假设为弹塑性，弹性模量210000MPa，泊松比0.3，塑性应力应变为圆柱体毛坯直径d=50mm，高度l=60mm；凸模直径D1=70mm；凹模直径D2=80mm；凸模从接触圆柱体上表面开始向下运动10mm；模具及板材之间的摩擦系数为0.1。

确定圆柱体变形后，凸模所受的反作用力大小.3、研究报告的预期结果：用abaqus建模，通过后期处理计算出应力大小。

第1节基本知识本节的有限元对象为轴对称问题，目的是学习将3D问题转化为2D问题分析的轴对称方法，涉及如何选取轴对称单元、建模规律、载荷的施加方法和后处理技术。

一、轴对称问题的定义轴对称问题是指受力体的几何形状、约束状态，以及其它外在因素都对称于某一根轴（过该轴的任一平面都是对称面）。

轴对称受力体的所有应力、应变和位移均对称于这根轴。

二、用ANSYS解决2D轴对称问题的规定用ANSYS解决2D轴对称问题时，轴对称模型必须在总体坐标系XOY平面的第一象限中创建，并且Y轴为轴旋转的对称轴。

求解时，施加自由约束、压力载荷、温度载荷和Y方向的加速度可以像其它非轴对称模型一样进行施加，但集中载荷有特殊的含义，它表示的是力或力矩在360°范围内的合力，即输入的是整个圆周上的总的载荷大小。

同理，在求解完毕后进行后处理时，轴对称模型输出的反作用力结果也是整个圆周上的合力输出，即力和力矩按总载荷大小输出。

在ANSYS中，X方向是径向，Z方向是环向，受力体承载后的环向位移为零，环向应力和应变不为零。

常用的2D轴对称单元类型和用途见表11-1。

表11-1 2D轴对称常用结构单元列表的高阶单的高阶单在利用ANSYS进行有限元分析时，将这些单元定义为新的单元后，设置单元配置项KEYOPT（3）为Axisymmetric(Shell51和Shell61单元本身就是轴对称单元，不用设置该项)，单元将被指定按轴对称模型进行计算。

后处理时，可观察径向和环向应力，它对应的是SX与SZ应力分量，并且在直角坐标系下观察即可。

可以通过轴对称扩展设置将截面结果扩展成任意扇型区域大小的模型，以便更加真实地观察总体模型的各项结果。

轴对称问题有限元分析实例 2D节2第p=1000 N/mF2y611xO61211-1 圆柱筒壳示意图图——圆柱筒的静力分析一、案例1问题，直0.1 m1000 N/m的压力作用，其厚度为如图11-1所示，圆柱筒材质为A3钢，受，并且圆柱筒壳的下部轴线方向固定，其它方向自由，试计算其变形、mm，高度为16 径12径向应力和轴向应力。

思考题5-1 轴对称问题的定义答：工程中又一类结构，其几何形状、边界条件、所受载荷都对称于某一轴线，这种情况下结构再载荷作用下位移、应变和应力也对称于这个轴线，这种问题成为轴对称问题。

5-2 轴对称问题一般采用的坐标系？作图说明每个坐标分量的物理意义答：在描述轴对称弹性体问题的应力及变形时常采用圆柱坐标r，θ，z。

5-3 轴对称问题中每个点有几个位移分量？各位移分量是那几个自变量的函数？答：位移分量u, w,都只是rz的函数，与θ无关。

5-4 轴对称问题中的每个点有哪几个应力分量？是那几个自变量的函数。

答：4个应力分量；5-5 轴对称问题中的每个点有哪几个应变分量？是那几个自变量的函数答：4个应变分量5-6 轴对称问题是三维问题？二维问题？最简单的轴对称单元是哪种单元？作图说明答：由于轴对称，沿θ方向的环向（周向）位移v等于零。

因此轴对称问题是二维问题；三角形环单元。

（三角形轴对称单元，这些圆环单元与r z平面（子午面）正交的截面是三角形）5-7 写出三角形环单元的位移函数。

满足完备性要求吗？答：满足完备性要求。

5-8 三角形环单元形函数的表达式？指出形函数的性质。

5-9 三角形环单元的应力和应变的特点。

其单元刚度矩阵是几阶的？答：应力分量：剪应力为常量，其他3个正应力分量均随位置变化；应变分量：面内（子五面）3个应变分量为常量，环向应变不是常应变，而是与单元中各点的位置有关。

单元刚度矩阵为六阶。

5-10 有限元方法求解对称问题的基本步骤？1.结构离散化：对整个结构进行离散化，将其分割成若干个单元，单元间彼此通过节点相连；2.求出各单元的刚度矩阵[K](e)：[K](e)是由单元节点位移量{Φ}(e)求单元节点力向量{F}(e)的转移矩阵，其关系式为:{F}(e)= [K](e) {Φ}(e);3.集成总体刚度矩阵[K]并写出总体平衡方程：总体刚度矩阵[K]是由整体节点位移向量{Φ}求整体节点力向量的转移矩阵，其关系式为{F}= [K] {Φ}，此即为总体平衡方程。

第1节基本知识本节的有限元对象为轴对称问题，目的是学习将3D问题转化为2D问题分析的轴对称方法，涉及如何选取轴对称单元、建模规律、载荷的施加方法和后处理技术。

一、轴对称问题的定义轴对称问题是指受力体的几何形状、约束状态，以及其它外在因素都对称于某一根轴（过该轴的任一平面都是对称面）。

轴对称受力体的所有应力、应变和位移均对称于这根轴。

二、用ANSYS解决2D轴对称问题的规定用ANSYS解决2D轴对称问题时，轴对称模型必须在总体坐标系XOY平面的第一象限中创建，并且Y轴为轴旋转的对称轴。

求解时，施加自由约束、压力载荷、温度载荷和Y方向的加速度可以像其它非轴对称模型一样进行施加，但集中载荷有特殊的含义，它表示的是力或力矩在360°范围内的合力，即输入的是整个圆周上的总的载荷大小。

同理，在求解完毕后进行后处理时，轴对称模型输出的反作用力结果也是整个圆周上的合力输出，即力和力矩按总载荷大小输出。

在ANSYS中，X方向是径向，Z方向是环向，受力体承载后的环向位移为零，环向应力和应变不为零。

常用的2D轴对称单元类型和用途见表11-1。

在利用ANSYS进行有限元分析时，将这些单元定义为新的单元后，设置单元配置项KEYOPT（3）为Axisymmetric(Shell51和Shell61单元本身就是轴对称单元，不用设置该项)，单元将被指定按轴对称模型进行计算。

后处理时，可观察径向和环向应力，它对应的是SX与SZ应力分量，并且在直角坐标系下观察即可。

可以通过轴对称扩展设置将截面结果扩展成任意扇型区域大小的模型，以便更加真实地观察总体模型的各项结果。

第2节 2D轴对称问题有限元分析实例图11-1 圆柱筒壳示意图 一、案例1——圆柱筒的静力分析问题 如图11-1所示，圆柱筒材质为A3钢，受1000 N/m 的压力作用，其厚度为0.1 m ，直径12 m ，高度为16 m ，并且圆柱筒壳的下部轴线方向固定，其它方向自由，试计算其变形、径向应力和轴向应力。

第11章轴对称问题的有限元分析第1节基本知识本节的有限元对象为轴对称问题，目的是学习将3D问题转化为2D问题分析的轴对称方法，涉及如何选取轴对称单元、建模规律、载荷的施加方法和后处理技术。

一、轴对称问题的定义轴对称问题是指受力体的几何形状、约束状态，以及其它外在因素都对称于某一根轴（过该轴的任一平面都是对称面）。

轴对称受力体的所有应力、应变和位移均对称于这根轴。

二、用ANSYS解决2D轴对称问题的规定用ANSYS解决2D轴对称问题时，轴对称模型必须在总体坐标系XOY平面的第一象限中创建，并且Y轴为轴旋转的对称轴。

求解时，施加自由约束、压力载荷、温度载荷和Y方向的加速度可以像其它非轴对称模型一样进行施加，但集中载荷有特殊的含义，它表示的是力或力矩在360 °范围内的合力，即输入的是整个圆周上的总的载荷大小。

同理，在求解完毕后进行后处理时，轴对称模型输出的反作用力结果也是整个圆周上的合力输出，即力和力矩按总载荷大小输出。

在ANSYS中, X方向是径向，Z方向是环向，受力体承载后的环向位移为零，环向应力和应变不为零。

常用的2D轴对称单元类型和用途见表11-1。

在利用ANSYS进行有限元分析时，将这些单元定义为新的单元后，设置单元配置项KEYOPT3）为Axisymmetric（Shell51 和Shell61单元本身就是轴对称单元，不用设置该项），单元将被指定按轴对称模型进行计算。

后处理时，可观察径向和环向应力，它对应的是SX与SZ应力分量，并且在直角坐标系下观察即可。

可以通过轴对称扩展设置将截面结果扩展成任意扇型区域大小的模型，以便更加真实地观察总体模型的各项结果。

第2节2D 轴对称问题有限元分析实例、案例1――圆柱筒的静力分析如图11-1所示，圆柱筒材质为A3钢，受1000 N/m 的压力作用，其厚度为 0.1 m ，直径12 m 高度为16 m,并且圆柱筒壳的下部轴线方向固定，其它方向自由，试计算其变形、 径向应力和轴向应力。

3 轴对称问题弹性力学空间问题中的轴对称问题是指，物体的几何形状、约束情况及所受的外力都对称于空间的某一根轴，因此在物体中通过该轴的任何平面都是对称面，所有应力、应变和位移也对称于该轴，这类问题称为轴对称问题。

研究轴对称问题时通常采用圆柱坐标系（r，θ，z），以z轴为对称轴。

轴对称问题实例如图3.1所示的受均布内压作用的长圆筒，通过Z轴的一个纵截面就是对称面图3.1受均布内压作用的长圆筒3.1 三角形截面环单元三结点单元位移函数图4-2 三结点单元轴对称问题分析中所使用的三结点单元，在对称面上是三角形，在整个弹性体中是三棱圆环，各单元中圆环形铰相联接。

三角形截面环单元的结点位移在轴对称问题中，弹性体内任意一点上，不存在切向位移，只存在径向位移u 和轴向位移w ，两个位移分量表示为，⎭⎬⎫⎩⎨⎧=w u f }{[][]Tmm j j i iT mT jT iew u w u w u==δδδδ}{单元结点位移轴对称问题的三结点三角形单元位移函数取为，⎭⎬⎫++=++=z r z r u 654321w αααααα⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡∆=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧m j i m jim j i m j i u u u c c c b b b a a a 21321ααα根据结点位移，可得：⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡∆=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧m j i m jim j i m j i w w w c c c b b b a a a 21654ααα单元形函数jm m j i r z z r a -=mmj ji iz r z r z r 11121=∆mj i z z b -=jm i r r c -=（i ，j ，m ）)(21z c r b a N i i i i ++∆=单元内任一点的位移{}[]{}em jim m j j i i m jim j iN N N w u w u w u N N N N N N w u f δ=⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎦⎤⎢⎣⎡=⎭⎬⎫⎩⎨⎧=00003.2 应变矩阵（几何矩阵）根据几何方程及单元内位移的表达式，可得：⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧∂∂+∂∂∂∂∂∂=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧r w z u z w ru r u zr z r γεεεθ应变矩阵)(21m m j j i i u b u b u b r u ++∆=∂∂)(21m m j j i i u f u f u f r u ++∆=rcz b r a f i i i ++=（下标轮换）)(21m m j j i i w c w c w c z w ++∆=∂∂)(21m m j j i i u c u c u c z u ++∆=∂∂)(21m m j j i i w b w b w b r w ++∆=∂∂应变矩阵[]{}em ji m m mm m jj jj j ii ii i zr z r B B B b c c f b b c c f b b c c f b δγεεεθ=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡∆=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧00000000021),,(00021][m j i b c c f b A B i i i iii ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=3.3 应力矩阵由轴对称问题的物理方程，得到弹性矩阵，⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡---------+-=)1(22100011101110111)21)(1()1(][μμμμμμμμμμμμμμμμμE D应力矩阵11A =-μμ2)1(221A =--μμ3)21)(1(4)1(A E=-+-μμμ令：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-+-=21111110010101)21)(1()1(][A A A A AA A E D μμμ则弹性矩阵为：]][[][B D S =][][m j iS S S S =),,()(2]][[][2211113m j i b A c A c f b A c A f b A c A f b A B D S i ii i i i ii i i i i i ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++∆==由弹性矩阵[D ]和几何矩阵[B ]可以得到应力矩阵[S ]，由应力矩阵可知，除剪应力为常量，其它三个正应力分量都是r 、z 的函数。