《压力容器有限元分析》由会员分享，可在线阅读，更多相关《压力容器有限元分析（32页珍藏版）》请在金锄头文库上搜索。
1、精品文档有限元的分析简介 随着科技的进一步发展，传 统的分析方法已不能满足现在社会的 需求，以及更不能满足一些问题的精 确分析， 而有限元的出现和应用 给机电、土木、航天等工业领域带来了历 史性的突破。 ANSYS 是 有限元的应用软件， 主要用于几何和网格划分、 多物理场、结构力学、 流体动力学、非线性结构、 仿真过程及数据管理、 显示动力学等多领 域的应用有限元法是求解工程科学中数学物理问题的一种通用数值方 法。本 书介绍有限元法的基本原理、建模方法及工程应用，强调理论 与实践的结 合。全书包括两篇共 16 章，第 1篇由第 1? 10 章组成， 介绍有限元法的 基本理论和方法，内容包括：
2、有限元法基本理论、平 面问题、轴对称问题 和空间问题、杆梁结构系统、薄板弯曲问题以及 热传导问题、结构动力学 问题、非线性问题的有限元法。有限元主要 介绍有限元建模技术及基于 ANS YS 勺有限元分析工程应用，内容包 括：有限元建模的基本流程、模 型简化技术、网格划分技术、边界条 件处理与模型检查以及基于 ANS YS 勺有限元分析工程应用实例。创新实践课题：压力容器的有限元应力分析与设计一、问题描述1、如图1所示为一台700立式储罐，其手孔的直径为88,材料为16MnR设计压力为13.5Mpa,工作压力为12.3Mpa,弹性模量 为201GPa泊 松比为0.3 ,要求利用有限元分析对此压力
3、容器进行应力分析设计。2、立式储罐用途：主要用于储存气体，如燃气等，因为储罐密封性能好且能承受较高的压力，所以将气体压缩成液体后，方便于储存在储罐内。二、设计基本参数如下表：弹性模量201Gpa工作压力12.3Mpa设计压力13.5Mpa储罐直径700手孔直径88泊松比0.3材料16MnR螺栓力82109N封头直径146壁厚34圆弧面直径18封头厚15立式储縫结构示意简图如下图所示:在压力容器的应力分析中，压力容器部件设计关心的是应力沿壁厚的分布规律及其大小，可采用沿壁厚方向的“校核线”代替校核截面。该容器轴对称，所以只需考虑对储罐上半部分进行分析设计。法 兰上的螺栓力可以转 化为一个集中力
4、F,且F=82109N三、结构壁厚计算1.筒体厚度2KS mpc设计厚度:dC2 C1名义厚度:n 34mm2.椭圆形封头厚度有效厚度:enC2C1标准椭圆封头计算厚度：R 0.0165dC2Ci设计厚度：名义厚度： 有效厚度:3.手孔厚度n 18mmenC2 C i计算厚度:PcDi有限元建模分析本次分析采用 ansys10.0 建立有限元分析和应力设计一、 GUI 操作方式定义工作文件名和工作标题（ 1）定义工作文件名：执行 change jobname ，文件名命名为 wuzu（ 2 ）定义工作标题： 执行 change title 命令，对文件的压力进 行分析（3）关闭三角坐标符号 定
5、义单元类型和材料属性（ 1 ）选择单元类型：在 elementtypes 命令中选择 strucral solid和 quad 8node82（ 2 ）设置单元选项：在 element type option 命令框中选择 k3 为 axisymmertic（ 3 ）设置材料属性：在 material number 命令框中设置 ex 为2.01e11 ， prxy 为 0.3Linear Isotropic Material Properties for Material Number 1EXPRXYT12.0111Add Temper JturDelete TemperatureGraph
6、/ UinedrlQtrQpic Pr-cprties for Material Number 1KP.jric elHelp二、建立几何模型(1)生成矩形面，在by dimensios中设置三个矩形面数据如下：350 ,384，0,28044,146,795.3,846.3,44，(2)生成部分圆环面：执行 partial a nn ulus命令框中设置62,600,742.8 wpx,wpy,theta1,rad,theta2,分别为 0,280,355,8.75,373,63 0,280,355,67,6383,90（3）面叠分操作：执行 booleans下的areas命令A3和A5的面
7、进行叠加（4）删除面操作：执行delete下的below命令，选择A6和A8的面（5）线倒角操作：执行lines下的lines fillet命令，选择要倒角 的两条线在rad文本框中输入20，这则另一倒角的两条线，在RAD文本中输入10（6）线生成面操作：点击by lines命令 选择如下 两组三条线 生 成两 个面（7）面相减操作，在boolsean下的areas选择A10和A3两个面生成如下图（8）面向加操作：在boolsean下的add aaread选择如下四个面生 成如下图 生成关键点：选择 A2面中的L5在line ratio选项框中输入0.348生成关(9) 生成关键点:在rati
8、o下输入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线，生成的图形如图所示。ANV J“吕 2J .益 1KJL11 JL _ 3L-IU1 上 Ti&uk.(12) 进行编号压缩操作(13) 合并关键点：执行 Merge Items命令，弹出图所示。选择Keypoints。(14)保存几何模型三、划分有
9、限元网格连接线过渡圆弧线：执行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) 设置全局单元尺寸：MeshTool命令，在SIZE文本框中输入20。OK C t(1)在线上施加面载荷：执行Pressurel On Lines 命令，拾取L4、L17、L14、L20、L30、L22、L8、L28、L24、L9 和 L7 的线
10、。弹出对话 框，在Load PRES value 文本框中输入13.5。(2)施加集中载何：执行On Nodes命令，拾取编号为708的关键点。弹出如图所示对话框。在然后在 Force/moment value文本框中输入82109Direction of force/mom 中选择 FY,(3)显示单元：执行Plot Element命令，生成结果如图所示(4) 求解运算：执行 Solve-Current LS 命令，浏览后执行 File/Close 命令，在Solution is done对话框时单击Close按钮完成求解运算v(5) 保存分析结果五、浏览运算结果(1)显示变形形状:W ；
11、W1 暑 Mkf?a- AiM.LysX de axcrt of Tfai 容器施加载荷就会产生力，而容器壁受到力的作用将产生形变，变形量随力的大小改变，所以图中所出现的变形不同是因为其受到的力的大小不同，所以容器的制作要根据实际用途设计，其各个部位所 使用的材料的量要合适。变形公式为：F=ps所出结果符合设计要求(2)显示节点位移云图:生成结果如图所示:AN&.flit iKIt-3 LJti a- Jfcse* 丄二 pjt 上 de ax-crt Gf Tfh .14SC-KEjSE-M.UZ11-E4昇J#眶tb舀容器受到力的作用容器壁产生变形就会出现节点的位移变化，变形量不同，节点位
12、移也不同，符合设计要求 显示节点的Von Mises应力：生成结果如图所示1-0$SL.SUItL.W131. ri!&41.511am-fiiS汕 l.曲 Ifcz LJti si AFr-Lj s de axcrc af Tlx s:对设计载荷作用下进行有限元分析，并对分析结果进行应力强度 评定。评定的依据为JB4732-199钢制压力容器一一分析设计标准应力线性化路径的选择原则为：（1）通过应力强度最大节点，并 沿壁厚方向的最短距离设定线性化路径；（2）对于相对高应力强度区，沿壁厚方向设定路径。设计工况下的评定线性化路径如下图所示，线性化结果如下图所示,都符合设计要求。（4 ）在路径上映
13、射数据的图形显示 生成如图所示:ANflhfciNlTI FUJItfSHE MJ7of Trzj :-TI1rWMlBdLBLltt a- Azsff- LA sX deaf Trii zAJtl- s- AL3ff_Ly2- d* jXFt(5)图形显示当前路径上应力线性化结果：以图形显示当前路劲上的所有应力分量、主应力及相关的应力强度BdTJANJATF-2心-Iuht=, . jurrawpm-If s .2-fe -(6)疲劳分析校核最高压力工况与最低压力工况下设备的最大应力强度均出现在接管与圭寸头相贯区的内壁，通过计算，疲劳校核通过。A List Linearized Stress
14、!分析结论附录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 COORDINATES. MEMBRANE *SXSYSZSXYSYZSXZ2.339124.6141.1-0.31360.0000.000S1S2S3SINTSEQV141.1124.6-2.340143.4135.9* BENDING* I=INSIDE C=CENTER O=OUTSIDESXSYSZS
