椭圆形封头与筒体连接区域结构进行了非线性有限元建模和仿真研究及ANSYS有限元分析.doc

符 号 说 明 D0 筒体外径，mm Di 筒体内径，mm δ 回转体壁厚，mm E2 强化模量，MPa a 椭圆形封头长半轴，mm b 椭圆形封头短半轴，mm μ 泊松比 x 边缘应力的作用范围，mm R1 封头中面第一曲率半径，mm R2 封头中面第二曲率半径，mm L 削薄长度，mm δ1 封头厚度，mm δ2 筒体厚度，mm δm 平均厚度，mm σs 材料屈服强度，MPa σb 材料强度极限，MPa σ1 第一主应力，MPa σ2 第二主应力，MPa σ3 第三主应力，MPa σr 径向应力，MPa σθ 周向应力，MPa σφ 经向应力，MPa m 椭圆系数 E1 弹性模量，MPa K 应力增强系数 p 介质压力，MPa 1 前言 1.1论文选题的背景和意义 封头是压力容器重要的受压元件之一是压力容器必不可少的重要组成部分。其种类很多，有椭圆形封头、碟形封头、无折边球面封头、锥形封头、带法兰凸形封头及平盖等。由于封头与筒体连接几何不连续部位，容器受作用时各部分的变形不一致，在连接处附近便产生附加的弯曲变形而这种附加的弯曲变形在局部区域引起的力比内压引起的应力大得多，称这部分应力为不连续应力。一些情况下，由于不连续应力与主要载荷作用下产生的应力为一个数量级的集中峰值应力，会引起设备失效。因而它是疲劳断裂、应力腐蚀破坏的根源，是压力容器中比较薄弱的部位。从而严重影响的承载能力，该部位很有可能成为设备的破坏源[]。无论国内还是国外的著名压力容器规范，都特别强调对此进行详细的应力分析。通过研究各种不连续结构建立的应力分析方法和强度评定对于确保压力容器的经济性和安全性 椭圆形封头是由半个椭球面和短圆筒组成。由于封头的椭球部分经线曲率变化平滑连续，故应力分布比较均匀，且椭圆形封头深度较半球形封头小得多，易于冲压成型，是目前中低压容器中应用较多的封头之一。标准椭圆形封头a/b=2，既便于加工，受力也比较有利[]，在工程中应用广泛。 本文从不同角度对非标准椭圆形封头与筒体连接结构的应力分布状况进行研究论证，从而有助于对这种特定结构强度性能做全面而深入的了解，探讨非标准椭圆形封头应力分析和设计计算方法的可靠性和安全性，以便在实际工程中加以应用。 目前研究压力容器封头与筒体连接结构采用的方法主要有解析计算法、法和实验研究法[]。这几种方法各有优缺点 1.2.1解析法 依据各种分析理论的解析方法，其解析计算结果便于实际工程设计采用。并且通过对各种表达式进行分析，可以得到各种结构参数变化对结构性能的影响，并且来比较各种设计方案和进行结构优化设计。但是解析法往往需要对实际结构进行合理抽象和简化，所建立的计算模型一般与结构的实际情况有差别，在各种假定基础上 1.2.2数值法 目前数值分析方法主要有有限单元法、边界元法、有限差分法等。其中有限元法已成为当今工程问题中应用最广泛的数值计算方法。数值计算可以有针对性地解决结构计算问题，但只能就已知尺寸和工况参数的结构进行数值计算，普遍性较差，而且无法对结构参数的变化进行分析如果要进行此类分析，则设计成本也非常高因此不如解析计算结果的普遍性强。而在计算精度方面，有限元分析中的网格划分、约束和载荷的正确简化处理里以及计算机的容量等因素都极大地影响结果的精度，从而使得设计成本大大提高 1.2.3实验法 采用与实际物体相同或按比例缩小的模型进行实际测量获得实验参数数据，以便归纳出经验公式用于工程设计中去但由于实验要受实验仪器、实验方法、实验装置的制造精度等条件的限制，而且在结构的峰值应力区，现有的实验技术还不能准确地反映实际结构的应力分布情况，费用也比较高。一般来说，适用于特别重要的核心设备或无法采用前两种方法设计的设备。 40多年以来随着理论分析方法的计算机技术的发展科技人员将有限元理论、数值计算技术和计算机辅助设计技术等相结合，开发出一系列通用的大型有限元分析设计软件。采用实验方法进行设计已经逐渐有减少的趋势[]。无论在国际还是在国内，各行各业中越来越多的有限元应用促进了产品设计和水平的提高。CAD）技术相结合。这些软件功能强大、使用方便、结果可靠成为解决涉及机械、土木、冶金、气象、宇航等工程问题强有力和灵活通用的工具其计算结果已经成为各类工业产品设计和性能分析的重要依据。其中比较常用的有：SAP、ADINA、ANSYS、ALGOR、NASTRAN、ABAQUS、COSMOS、MARC等。 对于已有成熟的分析方法的不连续结构，一般还是规范推荐的设计方法，而对于没有成熟的分析方法的不连续结构，主要还是借用常规设计的理论基础进行经验设计或实验设计的方法来解决[]。 另一种更科学更严密的设计规范是JB4732，它采用分析设计方法，要求对压力容器进行应力分析和疲劳分析，由于这种定