压力容器有限元分析报告.docx

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

2、面 问题、轴对称问题和空间问题、杆梁结构系统、薄板弯曲问题以及热 传导问题、结构动力学问题、非线性问题的有限元法。有限元主要介 绍有限元建模技术及基于ANSYS勺有限元分析工程应用，容包括：有 限元建模的基本流程、模型简化技术、网格划分技术、边界条件处理 与模型检查以及基于ANSYS勺有限兀分析工程应用实例创新实践课题：压力容器的有限元应力分析与设计一、 问题描述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图1在压力容器的应力分析中，压力容器部件设计关心的是应力沿壁 厚的分布规律及其大小，可采用沿壁厚方向的“校核线”代替校核截 面。该容器轴对称，所以只需考虑对储罐上半部分进行分析设计。法 兰上的螺栓力可以转化为一个集中力 F,且F=82109N三、结构壁厚计算1.筒体厚度计算厚度

4、：PCDi2KSmPc设计厚度：d C2C1名义厚度：n 34mm有效厚度：e n C2C12.椭圆形封头厚度标准椭圆封头计算厚度：R 0.0165设计厚度：dC2Cl名义厚度：n18mm有效厚度：en C2Cl3.手孔厚度有限元建模分析本次分析采用 ansys10.0 建立有限元分析和应力设计一、GUI操作方式定义工作文件名和工作标题（1）定义工作文件名：执行 change jobname ，文件名命名为 wuzu （ 2）定义工作标题： 执行 change title 命令，对文件的压力进 行分析（3）关闭三角坐标符号 定义单元类型和材料属性（ 1）选择单元类型：在 elementtype

5、s 命令中选择 strucral solid和 quad 8node82（ 2）设置单元选项：在 element type option 命令框中选择 k3 为 axisymmertic（ 3）设置材料属性：在 material number 命令框中设置 ex 为 2.01e11 ，prxy 为 0.3二、建立几何模型(1)生成矩形面，在by dimensios中设置三个矩形面数据如下：350 ,384，0,28044,146,795.3,846.3,44 ，62,600,742.8(2)生成部分圆环面：执行 partial a nn ulus 命令框中设置wpx,wpy,theta1,ra

6、d,theta2, 分别为 0,280,355,8.75,373,63 ，0,280,355,67,6383,90(3)面叠分操作：执行booleans下的areas命令A3和A5的面进行叠加(4)删除面操作：执行delete下的below命令，选择A6和A8的面(5)线倒角操作：执行lines下的lines fillet 命令，选择要倒角 的两条线在rad文本框中输入20，这则另一倒角的两条线，在 RAD 文本中输入10(6)线生成面操作：点击by lines 命令 选择如下 两组三条线 生 成两个面(7)面相减操作，在boolsean下的areas选择A10和A3两个面生 成如下图(8)面

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

8、2J .益1kJL11 JL _ 3L-IU 1 上 Ti&uk.(12)进行编号压缩操作(13)合并关键点：执行 Merge Items命令，弹出图所示。选择Keypoints。(14)保存几何模型三、划分有限元网格(1)连接线过渡圆弧线：执行Lines命令，拾取L22、L8和L28的 线，单击0K按钮。(2)设置单元尺寸：执行 MeshTool命令，弹出划单击Set按钮，拾取 L1、L3、L15、L13、L19、L25、L26、L10、L27 和 L31 的线，弹 出 Element Sixes on Picked Lines 对话框。输入 NDIV=4(3)设置全局单元尺寸：MeshTo

9、ol命令，在SIZE文本框中输入20。(4)戈吩网格：生成的有限元网格如图所示:ANmi iOlfliIIbOS U四、施加载荷并求解(1)在线上施加面载荷：执行 Pressurel On Lines命令，拾取L4、L17、L14、L20、L30、L22、L8、L28、L24、L9 和 L7 的线。弹出对 话框，在Load PRES value文本框中输入13.5。(2)施加集中载荷：执行 On Nodes命令，拾取编号为708的关键点。弹出如图所示对话框。在 Direction of force/mom 中选择FY,然后在Force/moment value文本框中输入82109。(3)显示

10、单元：执行Plot Element命令，生成结果如图所示(4)求解运算：执行Solve-Current LS命令，浏览后执行File/Close 命令，在Solution is done对话框时单击Close按钮完成求解运算(5)保存分析结果五、浏览运算结果(1)显示变形形状:-1wj moa- AiM.LysX de ax-crt of Tfi 容器施加载荷就会产生力，而容器壁受到力的作用将产生形变， 变形量随力的大小改变，所以图中所出现的变形不同是因为其受到的 力的大小不同，所以容器的制作要根据实际用途设计， 其各个部位所 使用的材料的量要合适。变形公式为：F=ps所出结果符合设计要求。(

11、2)显示节点位移云图：生成结果如图所示：& .flit iKI .14SC-KEjSE-M .UZ11-E4 昇 j#眶TB 舀t-3 LJtia- Jfcse*丄二pjt上 de ax-crt Gf Tfh 容器受到力的作用容器壁产生变形就会出现节点的位移变化， 变形量不同，节点位移也不同，符合设计要求(3)显示节点的Von Mises应力：生成结果如图所示1-0$ SL.SU ItL.W 131.41.511 am-fiiS 汕 l.曲 Ifcz LJti si Air-Lj s de ax-crc af Tlx s:对设计载荷作用下进行有限元分析，并对分析结果进行应力强度评定。评定的依据

12、为JB4732-199钢制压力容器一一分析设计标准应力线性化路径的选择原则为：（1）通过应力强度最大节点，并 沿壁厚方向的最短距离设定线性化路径； （2）对于相对高应力强度区，沿壁厚方向设定路径。设计工况下的评定线性化路径如下图所示， 线性化结果如下图所 示，都符合设计要求。（4）在路径上映射数据的图形显示:生成如图所示：.ti-s- d*iXrc of Trzj：iTXl-2 I -1 riHfciNlTI FUJI以图形显示当前路劲上(6)主应力及相关的应力强度(5)图形显示当前路径上应力线性化结果:的所有应力分疲劳分析校核最高压力工况与最低压力工况下设备的最大应力强度均出现在接管与圭寸头

13、相贯区的壁，通过计算，疲劳校核通过。分析结论附录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+14THE FOLLOWING X,Y,Z STRESSES ARE IN SECTION COORD

14、INATES.SXSYSZSXYSYZSXZ2.339124.6141.1-0.31360.0000.000S1S2S3SINTSEQV141.1124.6 -2.340143.4.9* BENDING* I=INSIDE C=CENTER O=OUTSIDESXSYSZSXYSYZSXZ2.4153.624MEMBRANE *I -8.7260.000C -2.5280.000O 3.6690.000S1 S23.624 2.4150.1974E-010.3828E-120.0000.000-2.376-3.6240.0000.0000.0000.0000.1974E-01 0.0003.669 -2.376S3 SINT-8.726 12.35-2.528 2.548-3.624 7.293SEQV11.792.5386.756* MEMBRANE PLUS BENDING * I=INSIDE C=CENTER O=OUTSIDESXI -11.07SYSZ SXY127.0SYZ144.7SXZ-0.31360.0000.000C -4.8