压力容器封头有限元分析

1、精选优质文档-倾情为你奉上有限元原理及工程计算课程设计说明书指导老师： 专业名称：工程力学学 号： 设计作者： 提交时间：2015.6.25 目 录一、设计任务书1.1题目压力容器封头有限元分析1.2结构图和计算参数1)、构件的几何形状和尺寸；2、材料及性能；3、约束和边界条件；4、载荷）1.3计算内容1） 构件的几何变形计算；2） 构件的各方向的位移计算；3） 构件的应力分布场计算；1.4计算要求1） 计算说明书应包括主要的计算步骤，计算公式、计算简图均应列入，并尽量详细描述计算过程；2） 计算说明书应书写清楚，字体工整，图表清晰规范；3） 在规定时间内完成计算并提交计算说明书。二、课程设计

2、指导书2.1目的要求本课程设计是有限元理论及其应用的重要实践环节之一，是一门了解有限元理论、应用有限元方法对工程问题进行计算的实战性很强的课程。通过本课程的学习，熟练掌握通用有限元软件ANSYS对工程问题的计算全过程，进一步综合和深化对有限元以及力学理论中基本理论和基本概念的理解。同时初步培养学生具备利用有限元方法对工程应用问题进行基本建模、加载、计算和结果分析的能力，使学生在计算、分析解决实践工程力学问题的能力有较大的提高，为学生在今后的毕业设计和进入社会工作后对工程问题的力学分析打下扎实的基础。2.2主要计算步骤1) 明确研究对象，合理简化研究对象的力学分析。a、 对实际构件进行合理简化，

3、确定主要研究对象；b、 对载荷和约束进行合理简化处理；c、 画出简化结构分析图；d、 了解构件的材料特性和几何特性。2) 启动ANSYA分析软件，设定ANSYS分析环境；a、 设置项目工程名；b、 设置计算过程、结果存取工作目录；3) 选择适当分析单元模型和材料模型，输入相关材料模型参数。a、 选择适当的单元模型，设置相关截面参数、实常数等；b、 选择适当的材料模型，设置材料模型的参数。4) 创建有限元模型a、 构思建模步骤，创建模型；b、 应用布尔运算对模型进行整合；c、 控制和合并重合单元、节点等。5) 网格划分a、 确定每部分构件的材料、单元；b、 选择网格划分方式；c、 调整已划分网格

4、；6) 加载和计算a、 设置载荷：集中力、分布力等；b、 设置约束：位移、对称面等；c、 设置计算参数并进行计算；7) 后处理a、 取后处理构件变形对比图；b、 取后处理构件节点位移图；c、 取后处理构件节点应力图，主要包括：第一主应力图、应力强度图、mises等效应力图等。2.3计算结果分析利用ANSYS后处理得到构件的变形图、位移图和应力图，针对这些图，分析最大应力发生点、最大变形发生点，分析该点的危险性，提出一个解决方案并进行计算对比。2.4计算说明书内容计算说明书是学生提交计算结果及其分析的重要报告，其主要内容包括下面几个方面：1) 构件结构图2) 材料和单元选择3) 有限元建模4)

5、网格划分5) 加载6) 计算和后处理7) 结果分析三、计算说明书3.1构件结构图图1 构件结构图3.2材料和单元选择表1 构件材料明细表材料弹性模量/Mpa泊松比屈服强度/Mpa筒体Q245R2.07E50.3245封头Q245R2.07E50.3245开口接管20#2.06E50.3245由上表可知，该模型有两种材料组成，但力学性能相差不大，为简便起见，可简化为同种材料，弹性模量为2e5Mpa，泊松比为0.3，选用常用的solid185单元。3.3实体（几何）建模1、建模的基本思路。采用自底向上的建模，即先建立关键点，通过关键点建立线，由线生成面，再由面旋转180度生成体。2、描述对工程实际

6、对象简化模拟方法，包括构件的几何简化和约束、载荷的处理方式等。因为压力容器为轴对称结构，取模型的一半进行建模，在对称面施加对称约束。3、建模后的模型图。图2 等轴侧立体图图3 左视图图4 俯视图3.4网格划分与有限元建模采用扫虐网格划分，划分后的模型如下：图5 模型的网格划分3.5加载在模型的对称面施加对称约束，上下表面y方向上的自由度为0，内壁施加0.58Mpa的压力。图6 边界条件加载图3.6计算和后处理图7 模型总位移图图8 第一主应力云图图9 应力强度云图图10 Mises等效应力云图应力集中处沿壁厚方向的三向应力：S1:径向应力 S2：周向应力 S3:轴向应力图11 三向应力图应力线

7、性化：如下图定义三个路径pathA，pathB，pathCCAB图12 强度评定路径设定（A,B,C）显示该路径上的应力线性化结果，如下图：PathA14图13 PathA上的应力分布图列表显示结果如下： * MEMBRANE * SX SY SZ SXY SYZ SXZ 1.387 9.972 54.56 -3.902 0.3871E-01 -0.8276 S1 S2 S3 SINT SEQV 54.57 11.48 -0.1316 54.70 49.92 * BENDING * I=INSIDE C=CENTER O=OUTSIDE SX SY SZ SXY SYZ SXZ I -0.81

8、58 -8.840 2.721 0.000 0.000 0.000 C -0.1303 -0.1865E-10 0.5740E-11 0.000 0.000 0.000 O 0.5551 8.840 -2.721 0.000 0.000 0.000 S1 S2 S3 SINT SEQV I 2.721 -0.8158 -8.840 11.56 10.26 C 0.000 0.000 -0.1303 0.1303 0.1303 O 8.840 0.5551 -2.721 11.56 10.32 * MEMBRANE PLUS BENDING * I=INSIDE C=CENTER O=OUTSI

9、DE SX SY SZ SXY SYZ SXZ I 0.5715 1.132 57.28 -3.902 0.3871E-01 -0.8276 C 1.257 9.972 54.56 -3.902 0.3871E-01 -0.8276 O 1.943 18.81 51.84 -3.902 0.3871E-01 -0.8276 S1 S2 S3 SINT SEQV I 57.29 4.757 -3.066 60.36 56.85 C 54.57 11.46 -0.2452 54.82 50.00 O 51.85 19.67 1.071 50.78 44.50 * PEAK * I=INSIDE C

10、=CENTER O=OUTSIDE SX SY SZ SXY SYZ SXZ I 0.3331E-15 0.2312 0.3252 2.450 0.7848E-01 -0.6109E-02 C 0.1915 -0.2060 -0.3578E-02 0.1394 0.1277E-02 -0.1329E-02 O 0.2220E-15 0.5352 -0.5339 -2.885 -0.8444E-01 0.1077E-01 15 S1 S2 S3 SINT SEQV I 2.570 0.3252 -2.339 4.908 4.256 C 0.2355 -0.3571E-02 -0.2501 0.4

11、856 0.4206 O 3.166 -0.5341 -2.631 5.797 5.084 * TOTAL * I=INSIDE C=CENTER O=OUTSIDE SX SY SZ SXY SYZ SXZ I 0.5715 1.363 57.61 -1.452 0.1172 -0.8337 C 1.449 9.766 54.56 -3.762 0.3999E-01 -0.8290 O 1.943 19.35 51.31 -6.787 -0.4573E-01 -0.8169 S1 S2 S3 SINT SEQV TEMP I 57.62 2.465 -0.5435 58.16 56.72 0

12、.000 C 54.57 11.21 -0.1148E-01 54.58 49.92 O 51.32 21.68 -0.4032 51.72 44.95 0.000 PathA图14 PathB上的应力分布图16列表显示结果如下： * MEMBRANE * SX SY SZ SXY SYZ SXZ 0.8554 2.488 56.28 3.027 0.8991 -0.2777 S1 S2 S3 SINT SEQV 56.30 4.801 -1.474 57.77 54.91 * BENDING * I=INSIDE C=CENTER O=OUTSIDE SX SY SZ SXY SYZ SXZ

13、 I -2.135 -3.909 7.396 0.000 0.000 0.000 C -0.5262 -0.8247E-11 0.1561E-10 0.000 0.000 0.000 O 1.082 3.909 -7.396 0.000 0.000 0.000 S1 S2 S3 SINT SEQV I 7.396 -2.135 -3.909 11.30 10.53 C 0.000 0.000 -0.5262 0.5262 0.5262 O 3.909 1.082 -7.396 11.30 10.19 * MEMBRANE PLUS BENDING * I=INSIDE C=CENTER O=O

14、UTSIDE SX SY SZ SXY SYZ SXZ I -1.279 -1.421 63.68 3.027 0.8991 -0.2777 C 0.3292 2.488 56.28 3.027 0.8991 -0.2777 O 1.938 6.396 48.89 3.027 0.8991 -0.2777 S1 S2 S3 SINT SEQV I 63.69 1.675 -4.388 68.08 65.26 C 56.30 4.616 -1.815 58.11 55.18 O 48.91 7.915 0.3988 48.51 45.22 * PEAK * I=INSIDE C=CENTER O

15、=OUTSIDE SX SY SZ SXY SYZ SXZ I -0.4441E-15 0.7497 -0.6557E-01 -2.636 0.6301E-01 0.4904E-01 C 0.7272 -0.5043 -0.6791E-01 -0.7179 0.5062E-02 0.1488E-02 O -0.2220E-15 0.5838 -0.8698 4.659 -0.7442E-01 -0.6157E-01 S1 S2 S3 SINT SEQV17 I 3.037 -0.6289E-01 -2.290 5.327 4.634 C 1.057 -0.6787E-01 -0.8343 1.

16、892 1.648 O 4.962 -0.8714 -4.376 9.338 8.170 * TOTAL * I=INSIDE C=CENTER O=OUTSIDE SX SY SZ SXY SYZ SXZ I -1.279 -0.6714 63.61 0.3913 0.9622 -0.2287 C 1.056 1.983 56.22 2.309 0.9042 -0.2763 O 1.938 6.980 48.02 7.686 0.8247 -0.3393 S1 S2 S3 SINT SEQV TEMP I 63.63 -0.4889 -1.477 65.11 64.62 0.000 C 56

17、.23 3.870 -0.8462 57.08 54.87 O 48.04 12.54 -3.641 51.68 45.78 0.000 PathC图15 PathC上的应力分布图18列表显示结果如下： * MEMBRANE * SX SY SZ SXY SYZ SXZ 2.495 -0.6592 39.24 2.230 -0.6205 -0.6114 S1 S2 S3 SINT SEQV 39.26 3.630 -1.815 41.07 38.64 * BENDING * I=INSIDE C=CENTER O=OUTSIDE SX SY SZ SXY SYZ SXZ I -1.396 -1

18、0.75 0.2959 0.000 0.000 0.000 C -0.9630 -0.2268E-10 0.6243E-12 0.000 0.000 0.000 O -0.5305 10.75 -0.2959 0.000 0.000 0.000 S1 S2 S3 SINT SEQV I 0.2959 -1.396 -10.75 11.04 10.30 C 0.000 0.000 -0.9630 0.9630 0.9630 O 10.75 -0.2959 -0.5305 11.28 11.16 * MEMBRANE PLUS BENDING * I=INSIDE C=CENTER O=OUTSI

19、DE SX SY SZ SXY SYZ SXZ I 1.099 -11.41 39.53 2.230 -0.6205 -0.6114 C 1.532 -0.6592 39.24 2.230 -0.6205 -0.6114 O 1.964 10.09 38.94 2.230 -0.6205 -0.6114 S1 S2 S3 SINT SEQV I 39.55 1.472 -11.80 51.35 46.17 C 39.26 2.901 -2.049 41.31 39.07 O 38.97 10.64 1.387 37.58 33.91 * PEAK * I=INSIDE C=CENTER O=O

20、UTSIDE SX SY SZ SXY SYZ SXZ I 0.6661E-15 1.931 0.4210 -3.336 -0.1954E-01 -0.1752E-01 C 1.444 -0.9774 0.8625E-01 -0.1336 0.5039E-02 0.3980E-02 O 0.6661E-15 1.979 -0.7660 3.871 -0.6152E-03 0.1599E-02 S1 S2 S3 SINT SEQV I 4.439 0.4212 -2.508 6.947 6.041 C 1.452 0.8627E-01 -0.9848 2.437 2.115 O 4.985 -0

21、.7660 -3.006 7.991 7.139 * TOTAL * I=INSIDE C=CENTER O=OUTSIDE SX SY SZ SXY SYZ SXZ I 1.099 -9.477 39.95 -1.106 -0.6401 -0.6289 C 2.976 -1.637 39.32 2.096 -0.6155 -0.6074 O 1.964 12.07 38.18 6.101 -0.6211 -0.6098 S1 S2 S3 SINT SEQV TEMP I 39.97 1.206 -9.602 49.57 45.15 0.000 C 39.34 3.769 -2.450 41.79 39.06 O 38.21 14.91 -0.9068 39.11 34.08 0.000 3.7、结果分析应力强度校核：Q245R的力学性能如下：表2 Q245R的力学性能钢板厚度/mm抗拉强度R/(N/mm2)屈服强度R/(N/mm2)伸长率A/%温度/冲击吸收能量KV2/J180°弯曲试验 弯曲直径（b35m