局部壁厚减薄CNG储气瓶应变分布与爆破压力研究.pdf

石油管材与仪器28PETROLEUM TUBULAR GOODS&INSTRUMENTS2023年12 月试验与研究局部壁厚减薄CNG储气瓶应变分布与爆破压力研究*李丽锋12,夏若琛”，廖凯”，许多林”，柳（1.中国石油集团工程材料研究院有限公司，石油管材及装备材料服役行为与结构安全国家重点实验室2.中国石油大学（华东）机电工程学院山东青岛2 6 6 58 0；3.中国石油西南油气田分公司蜀南气矿4.中国石油天然气股份有限公司玉门油田分公司机械厂甘肃酒泉7 352 0 0）摘要：针对钢质压缩天然气(CNG)储气瓶因局部壁厚减薄而发生强度失效的问题，采用有限元方法结合材料真实本构关系和弧长法分析了CNG储气瓶的爆破压力和减薄区域的应力应变分布，并与现有经验公式的计算结果对比分析。结果表明，随内压升高，减薄区域中心外壁面处首先发生塑性变形，而后塑性区不断向周围扩展；临近爆破压力时，减薄区塑性应变迅速增长，且最大塑性应变由外壁面向壁厚中心处移动；随着减薄深度增加,CNG储气瓶爆破压力下降幅度逐渐增加；对于局部减薄容器爆破压力的计算，经验公式相比有限元法的结果均偏保守，CB/T19624标准与有限元法的预测结果较为接近。关键词：CNG储气瓶；局部减薄；有限元分析；应变分析；爆破压力中图法分类号：TE972D0I:10.19459/ki.61-1500/te.2023.06.006(1.State Key Laboratory for Performance and Structure Safety of Petroleum Tubular Goods and Equipment Materials,CNPC Tubular Goods Research Institute,Xian,Shaanxi 710077,China;2.College of Mechanicaland Electrical Engineering,China University of Petroleum(East China),Qingdao,Shandong 266580,China;3.Shu nan Gas-mine Field,PetroChina Southwest Oil&Gasfield Company,Luzhou,Sichuan 646000,China;4.PetroChina Yumen Oilfield Company Machine Plant,Jiuquan,Gansu 735200,China)Abstract:In view of the strength failure of compressed natural gas(CNG)cylinder due to local wall thinning,the stress and strain distribu-tion in the thinned area and the blasting pressure of CNG cylinder are analyzed by using finite element method combined with the real consti-tutive relation of materials and arc length method,and the results are compared with those of existing empirical formulas.The results showthat with the increase of internal pressure,the plastic deformation occurs at the center of the thinning zone,and then the plastic zone ex-pands to the surrounding area.When the bursting pressure is about to be reached,the plastic strain in the thinning zone increases rapidlyand the maximum plastic strain moves from the outer wall to the center of the thickness.With the increasing of thinning depth,the decrea-sing range of bursting pressure of CNG storage cylinder increases gradually.For the calculation of blursting pressure in locally reduced ves-sels,the empirical formula is more conservative than that of the finite element method,and the prediction results of GB/T 19624 standardand the finite element method are close to each other.Key words:CNG storage cylinder;local thinning;finite element analysis;strain analysis;bursting pressure印磊”，杜洋，王维鑫4陕西西安710077;四川泸州646000;文献标识码：A文章编号：2 0 9 6-0 0 7 7(2 0 2 3)0 6-0 0 2 8-0 7Strain Distribution and Bursting Pressure of CNGCylinder with Local Wall ThinningLI Lifeng,XIA Ruochen,LIAO Kai,XU Duolin,LIU Lei,DU Yang,WANG Weixin初投稿收稿日期：2 0 2 1-0 7-18；修改稿收稿日期：2 0 2 3-0 7-2 5*基金项目：国家重点研发计划“典型危险化学品储存设施安全预警与防护一体化关键技术研究与应用示范”（编号：2 0 16 YFC0801200）；国家自然科学基金“基于率型损伤演化与流固耦合的内爆炸下管道裂纹动态扩展机理与预测研究”（编号：518 0 5544）第一作者简介：李丽锋，男，19 8 3年生，高级工程师，2 0 0 8 年毕业于北京航空航天大学材料学专业，现从事油气管道、储气库和储罐安全评估方向的研究工作。E-mail:2023年第9 卷第6 期0引言钢质储气瓶作为最为常见的天然气储存容器，其内的天然气通常以高压的形式储存。由于 CNG 储气瓶长期处于高压状态且天然气具有强烈的易燃易爆特性，通常认为CNG储气瓶具有发生泄漏燃爆事故的安全隐患 。钢质储气瓶的失效模式主要包括强度失效和泄漏失效两类,其中强度失效主要由于制造缺陷、材料性能不合格、疲劳、应力腐蚀开裂、管体腐蚀、机械损伤等因素导致承压能力下降致失效，同时包括充压过载失效。钢质储气瓶一旦发生强度失效，轻则造成财产损失、环境污染，重则导致人员伤亡。因此，开展储气瓶爆破压力预测研究对于其安全设计和安全评定尤为重要。当前容器爆破压力的研究主要有3种方法：理论研究、实验研究、有限元分析。随着容器结构及类型的复杂多样化,理论分析已越来越困难。实验研究需要进行大量爆破实验，成本较为高昂且由于焊接质量或材料自身缺陷的存在,结果存在一定随机性。近年来,有限元分析以高效、经济、较为准确等特点,逐步成为重要的压力容器爆破压力分析方法2-6 。刘萌等7-9 利用ABAQUS软件对复合材料罐缠绕层进行变厚度变角度的各向应力分布模拟，预测罐的爆破压强，优化设计得到的缠绕线型参数可有效提高复合材料罐的结构力学性能。D S.Kushan 等10-11 针对带接管压力容器，应用ANSYS软件分析接管对容器爆破压力的影响，提出计算容器爆破压力的拟合公式。目前多数研究集中在完整压力容器的性能评估与应力分析方面，在含局部缺陷容器的爆破压力分析方面12-16 ,Vladimir等用ANSYS模拟腐蚀缺陷对圆柱形容器爆破压力的影响,并对容器进行水压试验,试验结果与模拟结果吻合较好。在实际生产和运输过程中,压力容器可能会因机容积/m32外形尺寸（长度外径）/mm10 780559材料30CrMoCNG储气瓶在实际使用过程中，会因机械损伤或人为因素造成外部划伤，造成局部壁厚减薄。本研究中CNG储气瓶局部壁厚减薄设定：沿平行于储气瓶轴线的平面切除管道外表面部分材料，对切除形成的直角部分李丽锋等：局部壁厚减薄CNG储气瓶应变分布与爆破压力研究样,对容器的影响也难以评估。Nidhi Dwivedi 等10 曾总结现有压力容器的爆破压力计算公式，通过与有限元分析结果对比发现,没有一种理论公式可以适用于所有容器的爆破压力预测，应用有限元法估算容器的爆破压力是一条合理的途径。本文针对钢质CNG储气库因局部壁厚减薄而发生强度失效的问题，提出基于有限元计算局部减薄压力容器爆破压力的方法，建立局部壁厚减薄CNG储气瓶模型，将有限元模型分析结果与爆破试验结果进行对比，验证模型的准确性,基于此模型分析不同局部壁厚减薄深度下钢质CNG储气瓶爆破压力及应力应变分布的规律,并与现有经验公式的计算结果对比分析，为钢质CNG储气瓶的安全设计和安全评定提供支持。1研究对象与有限元模型1.1研究对象研究对象是一个容积为2 m的钢质CNG储气瓶,符合TSG212016标准，储存介质为压缩天然气，主要成分是甲烷（CH4），结构如图1所示，由前后端塞、排污管及瓶体3部分组成，设计参数见表1,CNG储气瓶材料为30CrMo,其力学性能参数见表2。10.78036523.846STLLLLLLL图1CNG储气瓶结构图表1CNG储气瓶设计参数设计压力/MPa工作压力/MPa27.525充装介质体积/m3充装质量/kg435.8291表2 30 CrMo力学性能数据表弹性模量/(Nm-2)泊松比2.11 10110.279倒圆角。具体尺寸参数为，减薄的轴向长度L等于CNG储气瓶半径,即2 7 9.5mm；壁厚减薄深度d分别取值为4.812、16、2 0 m m;倒圆半径与减薄深度一致。建立的局部壁厚减薄CNG储气瓶模型如图2 所示（本文中局部29械撞击摩擦和人为等因素产生外部划伤，造成容器局部壁厚减薄,影响容器的安全性及使用寿命。划伤形式多设计壁厚/mm23.84水压试验压力/MPa34.4屈服强度/MPa712SR279.5排污管设计温度/60气密封压力/MPa27.5抗拉强度/MPa845端塞30.减薄深度均为壁厚减薄的深度）。图2 局部壁厚减薄CNG储气瓶模型1.2爆破压力数值分析模型采用有限元分析软件ANSYS对局部减薄CNG储气瓶进行爆破压力分析,减薄CNG储气瓶从受压到爆破过程中,局部材料经历了弹性屈服强化断裂4个阶段的过程,本研究中采用多线性等向强化模型（MISO)将材料完整应力一应变曲线输人程序。根据ASMEVIII-22017（美国机械工程师协会锅炉和压力容器规范），可由材料的弹性模量、屈服强度以及抗拉强度拟合出材料的真实应力应变曲线17 。拟合公式如式(1)所示。(1)8t+1+2E式中：8,为真实应变；o,为真实应力，MPa;E,为弹性模%(1.0+2tanhH）;8 为微应变区的真实塑性应变；8 2 为宏观应变区的真实塑性应变；H为曲线拟合参数（详见ASMEVIII-2 2017-P202)。将30 CrMo的相关参数代人计算后即可得到材料的的真应力应变曲线，如图3所示。1000r800600400200图330 CrMo真实应力应变曲线CNG储气瓶的网格划分采用2 0 节点六面体高阶单元Solid186,Solid186可以具有任意的空间各向异性，支持塑性、超弹性、蠕变、应力钢化、大变形和大应变能力8 。局部壁厚减薄CNG储气瓶属于不规则模型,将其切分为3部分：储气瓶主体，减薄区，过渡区。其中储气瓶主体和减薄区采用六面体网格划分，过渡区采用四面体网格。因需要重点关注减薄区在加压过程中的应力应变分布及变化情况，对减薄区网格进行加密，具体的划石油管材与仪器分结果如图4所示。CNG储气瓶内部处于高压状态时,端塞或管线会对CNG储气瓶端部施加拉应力，在有限元分析中也应施加该等效拉应力,其大小可由式(2)计算：P=P.Sr-S.S2式中:P,为内壁面施加内压，MPa;S,为瓶口外径面积，mm;S,为瓶口内径面积,mm。图4局部壁厚减薄CNG储气瓶网格划分结果由于在加压至容器爆破的过程中,容器材料大部分区域进人塑性阶段,容器中产生明显的塑性变形,因此在进行有限元分析时需要考虑几何非线性，为准确捕捉到CNG储气瓶的塑性塌压力,CNG储气瓶爆破压力计算方法采用弧长法计算。弧长法是一种非线性求解的迭代控制方法，可以在荷载和位移增量均不确定的情况下，生成变化的增量值，在迭代求解过程中自动调节增量步长，跟踪各种复杂的非线性路径全过程，它是目前结构非线性分析中数值计算最为稳定、计算效率最高且最为可靠的迭代控制方法19-2 0 对局部壁厚减薄CNG储气瓶加载求解完毕后，选取减薄区域中心等效应变最大处节点（根据有限元分析结果该节点位于壁厚中心处），定义该节点等效塑形应变为自变量,内压载荷为因变量，绘制图形。如图5所示，当曲线开始出现下降段时,说明此时CNG储气瓶该位置发生了塑性垮塌,曲线最高点对应压力即为 CNG储气瓶爆破压力。0.050.1真实应变2023年12 月(2)0.15100r80F6040H20F%图5CNG储气瓶减薄区中心点等效塑性应变-压力曲线爆破压力0.10.2等效塑性应变0.32023年第9 卷第6 期李丽锋等：局部壁厚减薄CNG储气瓶应变分布与爆破压力研究31MN2模型验证在将上述有限元方法及参数设置应用于局部减薄CNG储气瓶爆破压力分析之前，首先需验证其有效性。由于缺乏该钢质CNG储气瓶的爆破压力实验数据，本文以浙江某单位生产的高压氧气瓶（几何特征与CNG储气瓶一致)为对象来进行模型与方法的验证。该高压氧气瓶几何参数如图6 所示，其材料为37Mn,工作压力为15MPa，容积为7 L。根据材料力学性能试验结果,该氧气瓶材料屈服强度,为535MPa,抗拉强度,为7 40 MPa,断后延伸率为22.3%，弹性模量为2 10 MPa,泊松比为0.3。材料真实本构关系由式（1)拟合得到,并采用多线性等向强化模型（MISO)输人ANSYS软件。600367SR5884.5图6 高压氧气瓶结构图对该氧气瓶建立1/4模型，并采用六面体单元进行网格划分,其中厚度方向采用4层网格，模型共划分为43.580个单元，如图7 所示。图7 高压氧气瓶网格划分结果使用福贝尔公式等理论方法估算出氧气瓶的爆破压力后，在内壁面施加压力载荷6 0 MPa，根据式（2）,在氧气瓶端部施加等效拉应力38.46 MPa，在模型对称面施加对称边界条件。为保证受力平衡不发生刚性位移，选取氧气瓶底部一节点施加轴向位移约束。求解方式选择弧长法，计算结束后绘制氧气瓶临近爆破时的等效应变云图，如图8 所示。选取氧气瓶等效塑性应变最大节点，绘制此节点等效塑性应变内压曲线，如图9 所示。当内压增至40 MPa左右时,塑性应变迅速增长,内压曲线对应的最大值即为氧气瓶的爆破压力，大小为51.9 MPa。MX0.0474450.1418290.253E-030.094637图：临近爆破时氧气瓶的等效应变云图60r50F4030H2010F70101025LSR800.2362130.1890210.283405爆破压力0.20.4.等效塑形应变图9 氧气瓶内压与塑性应变关系图该氧气瓶的水压爆破试验结果显示其爆破压力为51.55MPa,本文所建立有限元模型结果的相对偏差为0.68%，说明本文建立的压力容器爆破压力数值分析模型与方法具有较高的准确度。3结果分析与讨论3.1局部减薄区域塑性应变分析对不同局部减薄深度的CNG储气瓶内表面施加内压载荷，并在CNG储气瓶端部表面施加等效拉应力，同时约束其刚性位移并在模型对称面施加对称边界条件，对模型进行求解。为分析加压过程中局部减薄区域应变分布演化特征与规律,选取减薄深度为8 mm的CNG储气瓶为对象，分别绘制减薄区域沿轴向（路径1）、环向（路径2）和径向（路径3）的Mises等效塑性应变分布，路径选取如图10所示，等效应变分布分别如图11至图13所示。内压载荷为51.1MPa至爆破压力6 7.8 MPa。12图10 路径选取位置0.3305970.3777890.60.80.4249811.0320.250.200.150.100.05图11局部减薄区域轴向等效塑性应变0.25环向缺陷边缘位置0.200.150.100.050.0050图12局部减薄区域环向等效塑性应变0.25斤-P=67.8 MPa-P=67.5 MPa0.20P=66.7MPaP=65.5MPa0.15-P=64.1 MPa-P=62.4 MPa-P=60.4MPa0.10-P=57.7 MPaP=54.5MPa0.05-P=51.1 MPa0.005缺陷处径向距离/mm图13局部减薄区域径向等效塑性应变由轴向、环向路径等效塑性应变分布图（图11、图12)可知，当内压达到51.1MPa左右时,减薄区域中心处已发生塑性变形，随着内压逐渐增大，减薄区域塑性应变迅速增长,且塑性区范围由中心不断向外扩展。临近爆破压力时，减薄中心处塑性应变对内压变化极为敏感，此时该处等效塑性应变远大于远离减薄处筒体塑性应变。由拉美公式2 1 可知,对于仅受内压作用的完整圆筒压力容器,圆筒内壁面到外壁面的径向应力与环向应力逐渐减小,轴向应力保持不变，即完整圆筒从内壁面到外壁面等效应力逐渐减小，等效应变也逐渐减小。绘制CNG储气瓶局部减薄区域对称位置等效塑性应变沿径向的分布,局部减薄区域对称位置属于完整圆筒,如图14所示,可知CNG储气瓶未减薄区域内壁面到外壁面等效塑性应变逐渐减小，与一般完整圆筒的等效石油管材乌仪器应变分布规律相同。轴向缺陷-P=67.8 MPa边缘位置-P=67.5 MPa-P=66.7MPa-P=65.5 MPa-P=64.1 MPa-P=62.4 MPa-P=60.4 MPa-P=57.7MPa-P=54.5 MPa-P=51.1 MPa50100150200250300350400距缺陷中心轴向距离/mm-P=67.8MPa-P=67.5MPaP=66.7 MPa-P=65.5 MPa-P=64.1 MPa.-P=62.4MPa-P=60.4MPa-P=57.7 MPa-P=54.5MPa-P=51.1MPa2040距缺陷中心环向距离/mm10152023年12 月7104-P=67.8MPa6F.-P=67.5 MPaP=66.7 MPa5P=65.5MPa-P=62.4MPa-P=51.1 MPa25图14局部减薄区域对称位置径向等效塑性应变由径向路径等效塑性应变分布图（图13）可知，局部减薄处的径向等效塑性应变从内壁面到外壁面逐渐增大,外壁面首先开始发生塑性变形,且临近 CNG储气瓶爆破压力时，减薄中心的塑性应变最大值由外壁面逐渐向罐壁中心处移动。这与完整圆筒的等效应变分布规608010径向路径/mm1001201402015律存在明显不同,原因在于局部减薄导致了该区域几何特征的改变,进而导致了应力、应变分布的变化。3.2壁厚减薄对爆破压力影响规律分析使用有限元方法计算出不同局部壁厚减薄深度下CNG储气瓶的爆破压力，见表3。绘制出CNG储气瓶爆破压力随局部减薄深度的变化情况，如图15所示。表3不同减薄深度下CNG储气瓶爆破压力表局部减薄施加内压/MPa深度/mm048121620图15不同局部减薄深度下CNG储气瓶爆破压力由表3与图15可知，随着局部壁厚减薄深度的增大,CNG储气瓶爆破压力逐渐减小；当减薄深度为12 mm20最大求解时间/：1100.736900.834850.797800.697650.639300.8151008060402025爆破压力/MPa80.9675.0667.7555.7641.5424.451015壁厚减薄深度/mm1520252023年第9 卷第6 期（约1/2 壁厚）时，爆破压力下降幅度为（8 0.9 655.76)/80.96=31.1%；壁厚减薄深度越大，爆破压力的下降速率越快。该CNG储气瓶的设计压力为2 7.5MPa，当CNG储气瓶局部减薄深度达到2 0 mm时，已不满足设计要求。为分析理论方法在计算局部壁厚减薄容器爆破压CNG储气瓶局部减薄深度/mm0福贝尔公式P,=.2-)inK李生昌公式爆破压力P,=m(0.876gb+0.155g,)lnK/MPa史文逊公式2TGB/T 196242019GB/T19624一2 0 19 中计算容器爆破压力的公式如式(6)所示:g.+lnk (0.95-Po)P,=V30.013a)0.85+)cmin(3,a)bc1/3(6)6式中:P,为爆破压力，MPa;K为外径与内径之比;n为薄壁减薄系数，当K2时,值取1T为容器壁厚,mm;R为容器内径,mm;a,b,c为缺陷的相对轴向长度、相对环向宽度、相对深度。以有限元计算结果为基准，得出4种理论方法计算局部减薄CNG储气瓶爆破压力的偏差，如图16 所示。10%0-10%-20%-30%-40%-50%图16 4种理论方法计算偏差对比由表4和图16 可知，理论公式计算完整CNG储气瓶爆破压力偏差较小，在7%范围内；但对于局部减薄CNG储气瓶爆破压力的计算存在较大误差,且随着减薄深度的增大，偏差显著增大。以有限元计算结果为基准，3种经验公式和GB/T19624国家标准在不同减薄深度下计算结果的标准差分别为6.54、4.58、6.38、2.56，李丽锋等：局部壁厚减薄CNG储气瓶应变分布与爆破压力研究容器安全评定4种方法2-2 6 计算不同减薄深度下CNG储气瓶的爆破压力,运算时以内壁面到减薄区的最短距离为壁厚,内径保持不变，中心轴线到减薄区的最短距离为外径，计算结果见表4。表44种理论公式计算结果484.8471.1275.8363.5782.8968.9980.1469.37GB/T19624国家标准偏差最小，原因是GB/T19624在一定程度上考虑了CNG储气瓶局部减薄的相对尺寸。4种理论计算结果普遍比实际值偏小,主要原因是公式中采取的CNG储气瓶壁厚是局部减薄区壁厚的最小值,而局部减薄区的轴向中间位置壁厚最薄,两侧逐渐加厚,公式计算的CNG储气瓶壁厚比实际壁厚小,因此得到的CNG储气瓶爆破压力值也会普遍偏小。4结论本文建立了局部减薄钢质CNG储气瓶爆破压力有限元分析模型，结合高压氧气瓶爆破压力实验结果验证了模型的有效性,采用该模型分析了局部减薄深度对CNG储气瓶爆破压力的影响规律,以及随内压升高局部减薄区域应变分布的演化规律，可为储气瓶安全设计和安全评定提供技术支持。主要结论如下：1)随着局部壁厚减薄深度的增大,钢质CNG储气瓶福贝尔公式爆破压力逐渐减小，且减薄深度越大，爆破压力下降速李牛富公式使支逐公式CB/T19624十048121620储罐局部减薄深度/mm33力方面的准确度及局限性，采用福贝尔公式、李生昌公式、史文逊公式和CB/T196242019在用含缺陷压力857.2051.1355.0859.75率越快。对于2 5MPa设计压力，直径559 mm的钢质CNG储气瓶减薄深度为1/2 壁厚时，爆破压力下降31%24左右。2)随着内压升高,局部减薄区域中心外壁面处首先发生塑性变形,而后塑性区不断向周围扩展；临近爆破压力时，减薄区塑性应变迅速增长，且最大等效塑性应变由外壁面向壁厚中心位置移动,这与完整圆筒应力应变分布规律存在明显不同。3）福贝尔等经验公式法主要适用于完整圆筒爆破压力的计算，对于局部减薄压力容器存在适用局限性。1243.0838.5141.1748.281628.7425.6927.2635.322014.1912.6813.3519.2334.GB/T19624标准考虑了减薄区域的具体尺寸，相比其他经验公式,与有限元预测的爆破压力较为接近。对于局部减薄压力容器，经验公式及GB/T19624标准计算的爆破压力值均偏保守。参考文献1李丽锋,罗金恒,张超，等.CNG储气瓶物理爆炸裂纹动态扩展特性J.油气储运,2 0 2 2,41（10）：118 1-118 8.2马刚,许林滔,顾荣见，等.简单压力容器储气罐爆破压力公式研究J.中国化工装备,2 0 17,19（5）：31-34.3何廷梅.呼吸器氧气瓶爆破数值模拟J.煤矿安全，2 0 15，46(4):112-114.4】杨超，惠虎，黄淞.超高压爆破片设计爆破压力的理论计算与数值模拟的对比研究J.压力容器,2 0 2 0,37（3）：2 1-26,40.5杨蓉遵，李剑容，严伟华.有限元分析方法在压力容器检验中的应用J.化工装备技术,2 0 0 7,2 8（1：1417.6彦辉，李晓帅,王昱，等.基于显式有限元分析的X80钢压力容器爆破压力研究J.河北工业大学学报,2 0 17,46(2):42-48.7刘萌,祖磊,李书欣,等.复合材料球形气瓶非测地线缠绕线型设计和强度分析J.玻璃钢/复合材料,2 0 18（2）：8-14.8祖磊，王继辉,李书欣，等.金属内衬碳纤维全缠绕气瓶的有限元分析J.玻璃钢/复合材料,2 0 14（11）：30-35.9 PRAMOD R,KRISHNADASAN C K,SIVA SHANMUGAMN.Design and finite element analysis of metal-elastomerlined composite over wrapped spherical pressure vessel J.Composite Structures,2019,224:111028.10 KUSHAN D S,SANYAL S,BHOWMICK S.Parametric studyof interaction effect between closely-spaced nozzles in athin cylindrical pressure vessel J.International Journal ofPressure Vessels and Piping,2018,165:34-42.11姜雅洲.压力容器爆破压力数值模拟与试验研究D.杭州：浙江工业大学，2 0 15.12卢沛,张吉,李璟，等.储气罐爆破压力数值模拟分析J.化学工程与装备,2 0 18（9）：2 4-2 5,2 0.13陈凤.带接管压力容器极限压力及爆破压力的研究D.南京：南京工业大学,2 0 0 5.石油管材仪器14 CHMELKO V,BERTA I.Analytical solution of the pipe burstpressure using bilinear stress-strain model and influence ofcorrosion defects on it J.Procedia Structural Integrity,2019,18:600-607.15】徐君臣，吴云龙，沈，等.带几何缺陷拱顶罐球壳的稳定性分析J.化工装备技术,2 0 17,38（6)：2 8-32.16 DWIVEDI N,KUMAR V,SHRIVASTAVA A,et al.Burstpressure assessment of pressure vessel using finite element a-nalysis:a review J Journal of Pressure Vessel Technology,2013,135(4):044502.17 ASME.Section VIII-rules for construction of pressure ves-sels-division 2-alternative rules:BPVC-VIII-2-2017 S.NewYork:American Society of MechanicalEngineers,2017.18梁赛，杨冰，吴亚运，等.有限元方法中实体单元选择策略研究J.机械制造与自动化,2 0 19,48（2）：7 9-8 3.19 CHEN D,CAI Y H,LUO Y,et al.Measurement and analysison the thermal and mechanical transfer characteristics ofGTA arc affected by arc length J.Measurement,2020,153:107431.20 KHOE DINH D,LEE D H,SONG Y H,et al.Arc length con-trol for efficiency enhancement of energy usage in plasma dryreforming process J.Journal of C02 Utilization,2018,28:274-282.21】黄泽,郑津洋，刘鹏飞，等.内压圆筒厚度计算公式分析讨论J.压力容器,2 0 12,2 9(8):18-2 1,59.22郑传祥，文棋.低碳钢压力容器爆破试验及爆破压力公式研究J.压力容器,2 0 0 2,19(9):9-12.23刘小宁.压力容器爆破压力的计算J.化学世界,19 9 1(4):176-181.24刘小宁,刘岑，吴元祥，等.超高压圆筒形容器爆破压力计算公式的比较J.机械强度,2 0 15,37（2）：37 3-37 6.25王飞,邢璐，赵哲明，等.基于有限元法的压力容器爆破压力预测J.化工装备技术,2 0 16,37（4):42-44.26国家市场监督管理总局，中国国家标准化管理委员会.在用含缺陷压力容器安全评定：GB/T19624-2019S.北京：中国标准出版社，2 0 19.（编辑：马小芳）2023 年12 月