数值分析方法的新进展.docxVIP

数值分析方法的新进展 传统压力计的设计方法是将仅计算的名义膜电压或曲线电压与单调压力的比较，并将壁厚设计为简单的，但往往不符合实际重力的情况。压力容器在工作条件下,除有一次薄膜应力外,还有弯曲应力、热应力等。由于传统设计方法中没有进行详细应力分析,而只有通过简单的计算和对压力容器结构、制造工艺、使用运行等的强制规定以及使用经验的总结来保护设备,这种片面的做法可能给设备带来不安全因素。目前,在工程领域涌现出越来越多先进的压力容器计算分析方法,数值分析法就是其中之一,数值分析法广泛全面地引进新规范,这就要求工程设计人员对数值分析有个正确的认识。本文从压力容器的应力分类、分析评定准则及有限元分析方法进行论述,重点讲述了压力容器的分析评定过程,为从事压力容器的工程设计人员及制造厂家提供参考学习应用。 1 应力的分类 对压力容器进行应力分类是应力分析设计法中的一部分,也是决定其分析正确与否的关键部分,因此在对压力容器进行应力分类时,必须要掌握应力分类的一些基础知识。应力可以从不同角度进行分类:就其范围而言,可以分为总体应力和局部应力;按照沿壁厚的分布情况可以分为均匀分布(薄膜应力)、线性分布(弯曲应力)和非线性分布的应力;按其性质可以分为一次应力、二次应力和峰值应力,这些应力往往又相互交叉包含。 1.1 薄膜应力检测 一次应力是由外加机械荷载作用而产生的正应力或剪应力。它必须满足外载与内力的平衡关系,且无自限性,当其大大超过材料的屈服极限时,将引起容器的过渡变形而破坏。一次应力又分为一次总体薄膜应力Pm、一次弯曲应力Pb和一次局部薄膜应力P1。一次总体薄膜应力是指影响范围遍及整个结构的一次薄膜应力,由于内压在圆筒形或球形壳体中产生的薄膜应力就属于此应力;一次弯曲应力是平衡压力或其他机械荷载所需的沿厚度线性分布的应力;一次局部薄膜应力是由于内压或机械荷载在容器局部范围内所引起的薄膜应力,其应力水平大于一次总体薄膜应力,但影响范围只限于局部区域。 1.2 部件约束引起的热应力 二次应力是由于部件的自身约束或相邻部件的约束而产生的正应力或剪应力,它必须满足变形协调的要求,其主要特征是具有自限性,热应力就属于此类应力。 1.3 应力集中与结构 峰值应力是由荷载或几何形状的突变载接管根部、小的圆角半径或小孔边缘等处所引起的应力集中,其基本特征是在此应力作用下,结构不产生任何形式的变形,它仅是疲劳破坏或脆性断裂的可能根源。 压力容器典型零部件中的应力分类见JB 4732—1995《钢制压力容器——分析设计标准》中表4-1。 2 应力应力场评定线弹性分析在做 应力分析和应力分类的强度评定通常采用第三强度理论,进行评定时,首先在需要评定处(危险截面处)选取穿过容器壁厚的评定线,即定义路径,然后将线弹性分析得到评定线上各种应力分解为薄膜应力、弯曲应力和峰值应力,再求取应力强度(最大主应力和最小主应力的代数差值),按照不同的原则进行评定,如表1所示。表1中,Sm为设计应力强度,Sa为由疲劳曲线得到的许用应力强度幅。 3 应力的分类及评定 随着有限元分析软件的迅速发展,特别是其强大功能的完善,有限元分析越来越成为实施弹性应力分析的主要手段,如何对有限元计算得到的应力结果进行应力分类及评定是目前迫待解决的关键问题,也是国内外压力容器界讨论的热烈话题,在这方面,等效线性化处理方法是一个重要的进展,已为世界各国普遍采用,进入了通用有限元软件。等效线性化处理方法的基本思想来自材料力学和板壳理论中薄膜应力和弯曲应力沿截面均匀分布和线性分布及非线性分布的峰值应力的现象。 3.1 薄膜应力和等效弯曲应力的确定 由于任意主横截面上的任何部位的几何特征都可由该部位的中面曲率半径R和厚度S来描述,那么沿厚度上的各点就可通过中面法向坐标t来确定,如图1所示。 设沿截面厚度方向的实际分布应力为σ,平均应力(薄膜应力)和等效弯曲应力分别为σm和σb=a+bt,这里a和b为待定系数,连同σm一起共3个未知数。根据静力等效、静弯矩等效以及弯曲应力σb沿截面的合力为零的原则得: 式中,h=S/2。由式(1)可得σm的表达式;将σb=a+bt代入式(2)、式(3)可以求出a和b,并得到: 式中,k为中面曲率,k=1/R。 3.2 anasas仿真模型 有限元分析并不能直接得到应力沿路径分布的公式表达式,因此进行线性化时,应首先通过拟合路径上各点的应力得到应力分布曲线,然后积分得到I1、I2、I3的值,计算出各项应力。也可以通过分段数值积分得到,ANSYS软件就是通过分段数值积分各项应力值。在线性化时,ANSYS软件自动将所定义的路径平均分割为48份,如图2所示。用σi(i=1,2,…,49)表示各分段点位置σx、σy、σz、τxy、τyz、τzx6个单项应力中的一个,则沿路径的薄膜应力表达