压力容器应力分析的有限元分析.docxVIP

压力容器应力分析的有限元分析 随着电动汽车和化工行业的快速发展，压力容器的设计精度和高可靠性尤为重要。正确分析压力容器的强度是整个设计的关键，它关系到压力容器的工作性能和制造效率等许多问题。压力容器传统设计方法是将简单计算得到的名义薄膜应力或弯曲应力与单一的许用应力对比来设计其壁厚,这样固然简单,但是常常不符合实际受力情况。压力容器在工作条件下,除有一次薄膜应力外,还有弯曲应力、热应力等。由于传统设计方法中没有进行详细应力分析,而只有通过简单的计算和对结构、制造工艺、使用运行等的强制规定以及使用经验的总结来保护设备的安全性,这从某种程度上掩盖了事实的本质,并可能带来不安全因素。目前在工程领域涌现出越来越多的先进的计算分析方法,有限元分析法就是其中之一。本文首先对压力容器应力进行分类,然后以大型有限元分析软件ANSYS为平台,对实例进行计算分析及应力分类评定,重点论述了分析评定过程。 1 壁厚的分布情况 就应力的范围而言,可以分为总体应力和局部应力;按照沿壁厚的分布情况可以分为均匀分布(薄膜应力)、线性分布(弯曲应力)和非线性分布的应力;按其性质可以分为一次应力、二次应力和峰值应力,这些应力往往又相互交叉包含。 1.1 薄膜应力概述 是由外加机械荷载作用而产生的正应力或剪应力,必须满足外载与内力的平衡关系,且无自限性,当其大大超过材料的屈服极限时,将引起容器的过度变形而破坏。一次应力又分为一次总体薄膜应力(Pm)、一次弯曲应力(Pb)和一次局部薄膜应力(P1)。一次总体薄膜应力是指影响范围遍及整个结构的一次薄膜应力,内压在圆筒形或球形壳体中产生的薄膜应力就属于此应力;一次弯曲应力是平衡压力或其它机械荷载所需的沿厚度线性分布的应力;一次局部薄膜应力是由于内压或机械荷载在容器局部范围内所引起的薄膜应力,其应力水平大于一次总体薄膜应力,但影响范围只限于局部区域。 1.2 二次力q 是由于部件的自身约束或相邻部件的约束而产生的正应力或剪应力,必须满足变形协调的要求,其主要特征是具有自限性,热应力就属于此类应力。 1.3 角半径或孔径边缘应力集中 是由荷载或几何形状的突变载接管根部、小的圆角半径或小孔边缘等处所引起的应力集中,其基本特征是在此应力作用下,结构不产生任何形式的变形,仅是疲劳破坏或脆性断裂的可能根源。 2 应力的分类及评定 随着有限元分析软件的迅速发展,它越来越成为实施弹性应力分析的主要手段。如何对有限元计算得到的结果进行应力分类及评定是目前亟待解决的关键问题。而等效线性化处理就是一种解决方法,已为世界各国普遍采用,进入了通用有限元软件。等效线性化处理方法的基本思想来自材料力学和板壳理论中薄膜应力和弯曲应力沿截面均匀分布和线性及非线性分布的峰值应力理论。 2.1 弯曲应力调查的tdt法 由于任意主横截面上的任何部位的几何特征都可由该部位的中面曲率半径R和厚度S来描述,那么沿厚度上的各点就可通过中面法向坐标t来确定(见图1)。 设沿截面厚度方向的实际分布应力为σ,平均应力(薄膜应力)和等效弯曲应力分别为σm和σb=a+bt,这里a和b为待定系数,连同σm一起共3个未知数。根据静力等效、静弯矩等效以及弯曲应力σb沿截面的合力为零的原则得: σm∫h?h(R+t)dt=∫h?hσ(R+t)dt(1)∫h?hσb(R+t)tdt=∫h?h(σ?σm)(R+t)tdt(2)∫h?hσb(R+t)dt=0(3)σm∫-hh(R+t)dt=∫-hhσ(R+t)dt(1)∫-hhσb(R+t)tdt=∫-hh(σ-σm)(R+t)tdt(2)∫-hhσb(R+t)dt=0(3) 式中:h=S/2。 由式(1)可得σm的表达式;将σb=a+bt代入式(2)、式(3)可以求出a和b,并得到: σm=(I1+kI2)/S(4)σb=(?kS+12tS3)[I2+k12?k2S2(12I3?S2I1)](5)I1=∫h?hσdtI2=∫h?hσtdtI3=∫h?hσt2dt(6)σm=(Ι1+kΙ2)/S(4)σb=(-kS+12tS3)[Ι2+k12-k2S2(12Ι3-S2Ι1)](5)Ι1=∫-hhσdtΙ2=∫-hhσtdtΙ3=∫-hhσt2dt(6) 式中:k为中面曲率,k=1/R。 2.2 anasas仿真模型 有限元分析并不能直接得到应力沿路径分布的公式,因此进行线性化时,首先应拟合路径上各点的应力得到应力分布曲线,然后积分得到I1、I2、I3的值,计算出各项应力。也可以通过分段数值积分(ANSYS软件的方法)得到各项应力值。在线性化时,ANSYS软件自动将所定义的路径平均分割为48份(见图2)。用σi(i=1,2,…,49)表示各分段点位置σx、σy、σz、τxy、τyz、τzx6个单项应力中的一个,则沿路径的薄膜应力表达式为: