基于ANSYS理论地高压储气罐维静态分析.docxVIP

基于ANSYS的高压储气罐三维静力学分析 罗旭 201220679 动力工程 摘要：本文利用ANSYS软件对储气罐进行了静力学应力、应变分析，得到了储气罐的应力、应变分布图，所得结果与理论计算值能够较好的吻合，为高压储气罐的设计提供依据。 关键词：储气罐；有限元；应力；应变 1工程背景 高压储气罐用以储存压缩气体，保证在正常或应急情况下，提供足量的气体。高压储气罐具有建造结构简单、使用寿命长、对环境污染少、造价低等优点，随着工业生产的发展，储气罐在我国的应用越来越广泛。为提高储气罐的储存效率，人们总是尽可能提高储气罐内气体的压力和容积，因此对高压储气罐的结构进行全面的计算分析，具有重要的理论意义和工程实用价值。 0.060本文以如图1.1所示储气罐作为对象进行分析，储气罐材料选用16MnR低合金钢，并采用正火处理。储气罐筒体直径为2米，高4米，壁厚0.04米；封头为椭圆型封头，封头高度0.4米，壁厚0.04米， 0.060 2.080.0820.44 2.08 0.08 2 0.4 4 图1.1 储气罐结构简图 2 理论计算 2.1 筒体应力计算 传统的高压储气罐设计是利用薄壳理论进行计算的。根据该理论，储气罐器壁中的应力计算公式为： σ 其中：σ1 te pc Di φ——焊缝系数； [σ 表2.1 焊接接头系数φ值 焊缝形式 无损探伤要求 φ值 焊缝形式 无损探伤要求 φ值 双面焊或相当于双面焊的全焊透对接接头 100% 1.00 带垫板的单面焊或对接接头 100% 0.90 局部 0.85 局部 0.80 其中焊缝系数可根据表1选取。假设储液罐的所用焊缝接为双面焊的全焊透对接焊缝，焊缝进行100%探伤，则取焊缝系数φ=1。 由于储液罐材料选用16MnR低合金钢，钢材采用正火处理，查机械设计手册表32.1-12可得，材料的σb=470MPa，σs=285MPa，钢材安全系数?=1.25 代入数据可得，圆筒内壁的计算应力为： σ 故σ1 2.2椭球形封头 椭圆形封头中的应力，包括由内压引起的薄膜应力和封头与筒体连接处的不连续应力。按照Huggenberger对薄膜应力好和Coatses对不连续应力的分析，椭圆型封头中的最大应力对圆筒周向薄膜应力的比值可表示为（ab σ 其中：σ2 te pc Di φ——焊缝系数； [σ K——应力增强系数。 当ab K= 代入数据得： K= 故可得椭圆封头中的应力为： σ 故σ2 3 有限元分析 3.1建立模型 为简化计算，建立模型是不考虑筒体、封头和接管处的焊缝影响，同时去除支座。简化后，建立的储气罐模型如图3.1所示。 图3.1 储气罐三维模型 3.2 单元材料的选取及参数设置 在本次分析中，储气罐的筒体、接管和封头的单元类型均选用solid 185号8节点四面体单元。材料属性均为线弹性材料，材料的弹性模量为2.06×1011Pa，泊松比为0.3，密度为7850Kg/m3。 3.3网格划分 网格划分采用6面体映射网格，网格全局尺寸为0.05.划分网格后的储气罐模型如图3.2所示。 图3.2 网格划分 3.3载荷施加及求解 3.3.1施加约束条件 储气罐底部有裙座支撑，但在实际建模时并没有把裙座剪出来。因此在施加约束时约束裙座与储气罐下封头相连的部分，即约束下封头与筒体的共同边界线。执行Main Menusolutiondefine loadsapply structure displacementon lines命令，拾取底座地面的所有外边界线，并约束3个方向的自由度。 3.3.2载荷的施加及求解 选取储气罐的所有内表面，并在内表面施加5.7MPa的设计压力。模型约束和载荷的施加如图3.3所示。当完成载荷施加后，进入求解器求解。 图3.3 载荷施加 3.4 后处理分析 进入ANSYS的后处理模块可以得到一系列分析结果。后处理模块包括通用后处理模块和时间历程响应后处理模块。通用后处理模块可以用于查看整个模型或选定的部分模型在某一子步或时间步的结果。而时间历程响应后处理模块用于查看特定点在某一时间步的结果。 3.4.1查看变形结果 选择Main Menu/General Postproc/Plot Result/Nodal Solu命令,出现Contour Nodal Solution Date对话框,选择Nodal Solution/DOF Solution/displacement vector sum,得到位移分布图,如图3.4、3.5和3.6所示。根据分析结果可知，储气罐筒体的变形较小，而上封头和下封头的变形较大，最大变形为4.57mm，最大变形位于上封头的中心位置。 3.4 筒体应变分布图 3.5 上封头应变分布图 3.6 下封头应变分布图 3.4.2 查看应力分析