压力容器有限元分析

1、有限元的分析简介随着科技的进一步发展，传统的分析方法已不能满足现在社会的需求，以及更不能满足一些问题的精确 分析，而有限元的出现和应用 给机电、土木、航天等工业领域带来了历史 性的突破。ANSYS是 有限元的应用软件，主要用于几何和网格划分、多 物理场、结构力学、流体动力学、非线性结构、仿真过程及数据管理、显 示动力学等多领域的应用有限元法是求解工程科学中数学物理问题的一种通用数值方法。本书 介绍有限元法的基本原理、建模方法及工程应用，强调理论 与实践的结合。 全书包括两篇共16章，第1篇由第1 10章组成，介绍有限元法的基本理 论和方法，内容包括：有限元法基本理论、平面问题、轴对称问题和空间

2、问 题、杆梁结构系统、薄板弯曲问题以及热传导问题、结构动力学问题、非线 性问题的有限元法。有限元主要介绍有限元建模技术及基于ANS YS勺有限元分析工程应用，内容包括：有限元建模的基本流程、模型简化技术、网格划分技术、边界条件处理与模型检查以及基于ANS YS勺有限元分析工程应用实例。创新实践课题：压力容器的有限元应力分析与设计、问题描述1、如图1所示为一台700立式储罐，其手孔的直径为88,材料为 16MnR设计压力为13.5Mpa,工作压力为12.3Mpa,弹性模量为201GPa泊 松比为0.3 ,要求利用有限元分析对此压力容器进行应力分析设计。2、立式储罐用途：主要用于储存气体，如燃气等

3、，因为储罐密封性能好 且能承受较高的压力，所以将气体压缩成液体后，方便于储存在储罐内。二、设计基本参数如下表：弹性模量201Gpa工作压力12.3Mpa设计压力13.5Mpa储罐直径700手孔直径88泊松比0.3材料16MnR螺栓力82109N封头直径146壁厚34圆弧面直径18封头厚15立式储罐结构示意简图如下图所示:在压力容器的应力分析中，压力容器部件设计关心的是应力沿壁厚的分布规律及其大小，可采用沿壁厚方向的“校核线”代替校核截面。该容器轴对称，所以只需考虑对储罐上半部分进行分析设计。法兰上的螺栓力可以转化为一个集中力F,且F=82109N三、结构壁厚计算1筒体厚度计算厚度：PcD2KS

4、m设计厚度：d名义厚度：34mmn有效厚度：eP cnc2cc21cl2.椭圆形封头厚度标准椭圆封头计算厚度：设R 0.0165计厚度：名义厚度:有效厚度:3.手孔厚度18mmenC2Cj有限元建模分析本次分析采用ansys10.0建立有限元分析和应力设计一、GUI操作方式定义工作文件名和工作标题（1）定义工作文件名：执行change jobname，文件名命名为wuzu（ 2）定义工作标题：执行change title命令，对文件的压力进 行分析（3）关闭三角坐标符号定义单元类型和材料属性（1）选择单元类型：在elementtypes命令中选择strucral solid和 quad 8no

5、de82（2）设置单元选项：在element type option命令框中选择k3为axisymmertic（3）设置材料属性：在material number命令框中设置ex为2.01e11，prxy 为 0.3A Lirr IsctrcicforNumber 1Linear Isotropic Material Properties for Material Number 1EXPRXVAdd Temperature Delete TemperatureGraphKCane 2Help二、建立几何模型(1)生成矩形面，在by dimensios中设置三个矩形面数据如下：350 ,384,

6、 0,28044,146,795.3,846.3,44，62,600,742.8a nnwpx,wpy,theta1,rad,theta2,分别为 0,280,355,8.75,373,630,280,355,67,6383,90AN面叠分操作：执行booleans下的areas命令A3和A5的面进行叠加删除面操作：执行delete下的below命令，选择A6和A8的面线倒角操作：执行lines下的lines fillet命令，选择要倒角的两条线 在rad文本框中输入20，这则另一倒角的两条线，在RAD文本中输入10线生成面操作：点击by lines命令选择如下两组三条线生成两个 面面相减操作

7、，在boolsean下的areas选择A10和A3两个面生成 如下图 面向加操作：在boolsean下的add aaread选择如下四个面生 成如 下图(9)生成关键点：选择A2面中的L5在line ratio选项框中输入0.348生成关生成关键点:在ratio下输入line ratio=0.348,生成关键点21.生成线，拾取 “15,20”“22,19”“14, 17”“11, 21”“3, 13”“4, 16”“5,12”“21,5”“24, 18” 关键点生成九条线。再生成面。面分解操作：在Area by Line下，拾取A2, A3,A10面，拾取 L33,L27,L26线，生成的图

8、形如图所示。莒“ 2J .益 1KJL11 JL_3L-IU1 上 Ti&uk.进行编号压缩操作合并关键点：执行Merge Items命令，弹出图所示。选择KeypointSo保存几何模型三、划分有限元网格连接线过渡圆弧线：执行Lines命令，拾取L22、L8和L28的 线，单击0K按钮。设置单元尺寸：执行MeshTool命令，弹出划单击Set按钮，拾取 L1、L3、L15、L13、L19、L25、L26、L10、L27 和 L31 的线，弹出Element Sixes on Picked Lines对话框。输入 NDIV=4设置全局单元尺寸：MeshTool命令，在SIZE文本框中输入20。

9、戈吩网格：生成的有限元网格如图所示:、施加载荷并求解(1)在线上施加面载荷：执行Pressurel On Lines命令，拾取L4、L17、L14、L20、L30、L22、L8、L28、L24、L9 和 L7 的线。弹出对话 框，在Load PRES value文本框中输入135。 施加集中载荷：执行On Nodes命令，拾取编号为708的关键点。弹出如图所示对话框。在Direction of force/mom中选择FY,然后在Force/moment value文本框中输入82109。显示单元：执行Plot Element命令，生成结果如图所示ANLj . 求解运算：执行Solve-Cur

10、rent LS命令，浏览后执行File/Close命令，在Solution is done对话框时单击Close按钮完成求解运算v保存分析结果五、浏览运算结果AN-1(1)显示变形形状:a- AiMLysX de a-crtof Tfai W； W1 暑 Mkf?容器施加载荷就会产生力，而容器壁受到力的作用将产生形变，变形量随力的大小改变，所以图中所出现的变形不同是因为其受到的力的大小不 同，所以容器的制作要根据实际用途设计，其各个部位所使用的材料的量要 合适。变形公式为：F=ps所出结果符合设计要求。(2)显示节点位移云图：生成结果如图所示：&.flit iKI.14SC-KEjSE-M.U

11、Z11-E4昇 j#眶 TB 舀t-3 LJtia- Jfcse*二 pjt 上 de ax-crt Gf TF围容器受到力的作用容器壁产生变形就会出现节点的位移变化， 变形量不同，节点位移也不同，符合设计要求（3）显示节点的Von Mises应力：生成结果如图所示1-0$SL.SUItL.W131. ri!&41.511am-fiiS汕 l.曲 Ifcz LJti si AiM-Lj s de ax-crc af Tlx s:对设计载荷作用下进行有限元分析，并对分析结果进行应力强度评定。评定的依据为JB4732-199钢制压力容器一一分析设计标准应力线性化路径的选择原则为：（1）通过应力强度

12、最大节点，并沿壁厚方向的最短距离设定线性化路径；（2）对于相对高应力强度区，沿壁厚方向设定路径。设计工况下的评定线性化路径如下图所示，线性化结果如下图所示， 都符合设计要求。（4）在路径上映射数据的图形显示： 生成如图所示：flhfciNLTI FUJIANtfSHE MJ7.AJtl- s- AL3ff_Ly2- d* iXATtof Trzj: (5) 图形显示当前路径上应力线性化结果:以图形显示当前路劲上AN的所有应力分量、主应力及相关的应力强度BdTJJATF-2af Trihz心-Ita- Azsff-LA sX de uht= .jurraw pm-If s.2-fe -疲劳分析校

13、核最高压力工况与最低压力工况下设备的最大应力强度均出现在接管与圭寸头 相贯区的内壁，通过计算，疲劳校核通过。A附录PRINT LINEARIZED STRESSTHROUGA SECTION DEFINED BY PATH= PATH2DSYS= 0* POST1 LINEARIZED STRESS LISTING *INSIDE NODE = 395 OUTSIDE NODE =350LOAD STEP 1 SUBSTEP= 1TIME= 1.0000 LOAD CASE= 0* AXISYMMETRIC OPTION * RHO = 0.28417E+14SXSYSZSXYSYZSXZ2.

14、339124.6141.1-0.31360.0000.000S1S2S3SINTSEQV141.1124.6-2.340143.4135.9* BENDING * I=INSIDE C=CENTER O=OUTSIDESXSYSZSXYSYZSXZI -8.7260.0000.0000.000C -2.5280.1974E-010.3828E-120.0000.0000.000-2.626O 3.6690.0000.0000.000S1S2S3 SINTSEQV3.6242.415-8.72612.3511.790.1974E-01 0.000-2.5282.5482.5383.669-2.3

15、76-3.6247.2936.756THE FOLLOWING X,Y,Z STRESSES ARE IN SECTION COORDINATES. MEMBRANE *SX SYSZ SXYSYZSXZI -11.07127.0144.7-0.31360.0000.000C -4.867124.6141.1-0.31360.0000.000O 1.330122.3137.5-0.31360.0000.000S1 S2S3 SINTSEQVI 144.7 127.0-11.07 155.8147.7C 141.1 124.6-4.868 146.0138.5O 137.5 122.31.330

16、 136.1129.2* PEAK *I=INSIDE C=CENTER O=OUTSIDESX SYSZ SXYSYZSXZI 0.1776E-1411.002.67712.040.0000.000C 4.155-4.801-1.044-1.1390.0000.000O 0.0009.3332.218-7.3230.0000.000S1 S2S3 SINTSEQVI 18.74 2.677-7.742 26.4823.11C 4.297 -1.044-4.944 9.2418.035O 13.35 2.218-4.017 17.3715.24* TOTAL *I=INSIDE C=CENTE

17、R O=OUTSIDESX SYSZ SXYSYZSXZI -11.07138.0147.411.730.0000.000C -0.7129119.8140.0-1.4520.0000.000O 1.330131.6139.7-7.6370.0000.000S1 S2S3 SINTSEQVTEMPI 147.4139.0-11.98159.4155.30.000C 140.0 119.9-0.7304 140.8131.8O 139.7132.00.8841138.8135.10.000附录二PRINT LINEARIZED STRESS THROUGH A SECTION DEFINED B

18、Y PATH= PATH2* POST1 LINEARIZED STRESS LISTING *INSIDE NODE = 395 OUTSIDE NODE = 350LOAD STEP 1 SUBSTEP= 1TIME= 1.0000 LOAD CASE= 0* AXISYMMETRIC OPTION *RHO = 0.28417E+14I=INSIDE C=CENTER O=OUTSIDETHE FOLLOWING X,Y ,Z STRESSES ARE IN SECTION COORDINATES.* MEMBRANE *SX 2.339S1 141.1SY124.6S2124.6SZ1

19、41.1S3 -2.340SXY-0.3136SINT143.40.000S135.9SYZ0.000EQVSXZ* BENDING *I=INSIDE C=CENTER O=OUTSIDESXSYSZSXYSYZSXZI -8.7262.4153.6240.0000.0000.000C -2.5280.1974E-010.3828E-120.0000.0000.000O 3.669-2.376-3.6240.0000.0000.000S1S2S3SINTSEQVI 3.6242.415-8.72612.3511.79C 0.1974E-010.000-2.5282.5482.538O 3.6

20、69-2.376-3.6247.2936.756SXSYSZSXYSYZSXZI -11.07127.0144.7-0.31360.0000.000C -4.867124.6141.1-0.31360.0000.000O 1.330122.3137.5-0.31360.0000.000S1S2S3SINTSEQVI 144.7127.0-11.07155.8147.7C 141.1124.6-4.868146.0138.5O 137.5122.31.330136.1129.2* PEAK *I=INSIDE C=CENTER O=OUTSIDESXSYSZSXYSYZSXZI 0.1776E-

21、1411.002.67712.040.0000.000C 4.155-4.801-1.044-1.1390.0000.000O 0.0009.3332.218-7.3230.0000.000S1S2S3SINTSEQVI 18.742.677-7.74226.4823.11MEMBRANE PLUS BENDING *I=INSIDE C=CENTER O=OUTSIDESXSYSZI-11.07138.0147.4C-0.7129119.8140.0O1.330131.6139.7S1S2S3I147.4139.0-11.98C140.0119.9-0.7304O139.7132.00.88

22、41* TOTAL t TOC o 1-5 h z SXYSYZSXZ11.730.0000.000-1.4520.0000.000-7.6370.0000.000SINTSEQVTEMP159.4155.30.000131.8135.10.000-4.9442233 14i(2479.2418.035C 4.297 -1.044O 13.35附录三PRINT LINEARIZEDSTRESS THROUGHA SECTION DEFINEDBY PATH= PATH3* POST1 LINEARIZED STRESS LISTING *INSIDE NODE = 294 OUTSIDE NO

23、DE = 285LOAD STEP 1 SUBSTEP= 1TIME= 1.0000 LOAD CASE= 0* AXISYMMETRIC OPTION *RHO = 0.81139E+13THE FOLLOWING X,Y ,Z STRESSES ARE IN SECTION COORDINATES.* MEMBRANE *SXSYSZSXYSYZSXZ41.42-1.26199.0617.350.000 0.000S1S2S3SINTSEQV99.0647.58-7.422106.592.24* BENDING *I=INSIDE C=CENTER O=OUTSIDESXSYSZSXYSY

24、ZSXZI -53.125.621-50.890.0000.0000.000C -2.6260.5370-0.8475E-100.0000.0000.000O 47.86-4.54750.890.0000.0000.000S1S2S3SINTSEQVI 5.621-50.89-53.1258.7457.65C 0.53700.000-2.6263.1632.932O 50.8947.86-4.54755.4353.99I=INSIDE C=CENTER O=OUTSIDESX* PEAK tSYI=INSIDE C=CENTER O=OUTSIDESZSXYSYZSXZI -0.1243E-1

25、3-5.756-5.447-19.170.0000.000C 3.3030.86490.59290.18610.0000.000O 0.1421E-13-1.564-1.0847.2120.0000.000S1S2S3SINTSEQVI 16.50-5.447-22.2638.7633.67C 3.3170.85080.59292.7242.605O 6.472-1.084-8.03614.5112.57SX* TOTAL *SYI=INSIDE C=CENTER O=OUTSIDESXZSZSXYSYZI-11.70-1.39642.73-1.8180.0000.000C42.090.141

26、499.6617.530.0000.000O89.28-7.371148.924.560.0000.000S1S2S3SINTSEQVTEMPI42.73-1.084-12.0154.7450.180.000C99.6648.46-6.221105.991.71O148.995.16-13.25162.1143.00.000SXSYI-11.704.360C38.79-0.7235O89.28-5.807S1S2I48.1815.45C99.0645.33O149.992.35SZSXYSYZSXZ48.1817.350.0000.00099.0617.350.0000.000149.917.

27、350.0000.000S3SINTSEQV-22.7970.9661.52-7.258106.392.08-8.873158.8139.3附录四BY PATH= PATH6PRINT LINEARIZED STRESS THROUGH A SECTION DEFINED* POST1 LINEARIZED STRESS LISTING *INSIDE NODE = 360 OUTSIDE NODE = 490LOAD STEP 1 SUBSTEP= 1TIME= 1.0000 LOAD CASE= 0SXSYSZSXYSYZSXZ80.8137.30134.0-59.980.000 0.00

28、0S1S2S3SINTSEQV134.0122.9-4.748138.7133.5* BENDING *I=INSIDE C=CENTER O=OUTSIDESXSYSZSXYSYZSXZI -4.3994.882-3.2570.0000.0000.000C 0.1228-0.4937E-01 -0.4825E-120.0000.0000.000O 4.645-4.9803.2570.0000.0000.000S1S2S3SINTSEQVI 4.882-3.257-4.3999.2818.766C 0.12280.000-0.4937E-010.17220.1536O 4.6453.257-4.9809.6259.012* MEMBRANE PLUS BENDING *I=INSIDE C=CENTER O=OUTSIDESX