基于ANSYS的管道屈曲临界载荷分析

康习锋, 张宏

康习锋, 张宏. 基于ANSYS的管道屈曲临界载荷分析[J]. 油气储运, 2017, 36(3): 262-266. DOI: 10.6047/j.issn.1000-8241.2017.03.004
引用本文: 康习锋, 张宏. 基于ANSYS的管道屈曲临界载荷分析[J]. 油气储运, 2017, 36(3): 262-266. DOI: 10.6047/j.issn.1000-8241.2017.03.004
KANG Xifeng, ZHANG Hong. Analysis on critical buckling load of pipelines based on ANSYS[J]. Oil & Gas Storage and Transportation, 2017, 36(3): 262-266. DOI: 10.6047/j.issn.1000-8241.2017.03.004
Citation: KANG Xifeng, ZHANG Hong. Analysis on critical buckling load of pipelines based on ANSYS[J]. Oil & Gas Storage and Transportation, 2017, 36(3): 262-266. DOI: 10.6047/j.issn.1000-8241.2017.03.004

基于ANSYS的管道屈曲临界载荷分析

详细信息
    作者简介:

    康习锋, 男, 1977年生, 讲师, 2008年硕士毕业于武警学院作战指挥学专业, 现主要从事管道安全工程、消防装备等技术的研究和教学工作。地址: 河北省廊坊市西昌路220号中国人民武装警察部队学院消防指挥系技术装备教研室, 065000。电话: 13833663121, Email: kangxifeng@outlook.com

    通讯作者:

    张宏, 男, 1963年生, 教授, 2003年博士毕业于中国石油大学(北京)机械设计及理论专业, 现主要从事油气管道的安全评价方面研究工作。地址: 北京市昌平区府学路18号中国石油大学(北京)机械与储运工程学院, 102249。电话: 010-89733274, Email: zhanghong630215@139.com

  • 中图分类号: TE88

Analysis on critical buckling load of pipelines based on ANSYS

  • 摘要: 为了研究管道初始横向位移对其临界屈曲载荷的影响,采用数值分析方法,进行屈曲分析。首先利用特征值法进行屈曲求解,获得屈曲模态,并以此作为初始位移输入的基本参数;其次结合管道的几何参数,应用ANSYS有限元软件建立管道模型,并求解非线性大变形模型,分析初始横向位移对管道屈曲临界载荷的影响。结果表明:初始横向位移与管道的屈曲临界载荷是负相关的,初始位移越大,临界屈曲载荷越小,且几乎呈线性关系。研究结果可为今后的管道设计和工程建设提供技术支持。
    Abstract: In this paper, buckling analysis was performed by numerical analysis method to analyze the influence of initial transversal displacement of pipelines on the critical buckling load. Firstly, buckling solution was obtained by the eigenvalue method, and buckling modal was established and taken as the basic input parameter of initial displacement. Secondly, based on the geometric parameters of pipelines, a pipeline model was set up by ANSYS finite element software. And finally, nonlinear large-deformation model was solved and the influence of initial transversal displacement on the critical buckling load of pipelines was analyzed. It is shown that the initial transverse displacement is negatively correlated with the critical buckling load of pipelines and their relationship is nearly linear. The larger the initial transverse displacement is, the lower the critical buckling load is. The research results can provide the technical support for pipeline design and engineering construction in the future.
  • 油气管道一般由钢材制成,具有热膨胀性,投入运行后,受内部油品传递热量和压力的共同作用,会在轴向发生膨胀。膨胀过程中,一旦端部存在约束,或受周边土壤摩擦力作用,管壁会在轴向产生压缩应力集中,当其达到一定值时,将导致管道发生屈曲[1-2]。在管道建设过程中,复杂的地形特征、施工技术及管理水平等均有可能造成初始位移的产生[3]

    随着冻土区管道的不断建设,管道屈曲现象也时常出现,对屈曲机理及量化评估模型的研究亟待解决[4-6]。由于埋地管道涉及管道与土壤的相互作用,而土壤成分多样,且随季节变化含水量不同,土壤的力学特性也随之改变。这些困难导致管道屈曲机理的研究短时间内难以取得突破,采用解析方法求得精确理论解的难度也较大。

    通过简化,早期有学者找到了屈曲模型的解析解[5-10]。Hobbs在铁轨屈曲分析的基础上,进行了改进,将梁的屈曲模型应用在管道分析[11-13],采用了小变形的假设,建立了二阶微分方程,求解得出轴向载荷P与屈曲段波长L成平方反比关系。在实际工程应用中,由于现场环境复杂、现有测量方法的限制,屈曲波长难以确定[14-16]。Allen通过小尺寸实验进行了模拟分析,在确定边界条件求解方程时,引入了初始位移,得出载荷的平方与初始缺陷成反比关系[17-18]。Heedo等通过研究初始缺陷形状特征与临界载荷的关系发现,缺陷形式不仅对临界载荷产生影响,还可能影响管道的屈曲方式[19-21]。在此采用有限元数值计算方法,分析初始位移对管道临界载荷的影响。

    管道建模采用ANSYS中的beam188梁单元,为三维两节点单元,适于对细长杆件、柱体进行线性或大应变的非线性分析。材料的弹性模量200 GPa,泊松比为0.3,共有400个单元。

    对管道受力情况进行分析(图 1),其中p为管道受到的轴向载荷;V0为管道初始缺陷,即横向最大位移;L为屈曲段管道的长度;g为线密度;v为管道沿y轴的横向位移;k为土壤对管道向下的约束反力系数。土壤产生的约束反力通常是与位移相关且呈非线性,在综合考虑平衡简化与模型精度的前提下,对土壤反力进行分段线性化处理(图 2)。V0点土壤约束力达到最大值,随后线性下降至C点,土壤约束力为0。V0取管道埋深的中点,C为管道埋深。

    图  1  含缺陷管道受力示意图
    图  2  土壤约束反力模型

    管道模型长200 m,直径为700 mm(表 1图 3)。节点1的转动和平移自由度均设定为0,末节点401的x方向无约束,其余节点具有xy轴的平移及z轴的转动自由度,载荷作用于末节点401,为负x轴方向。根据特征值法求解得到的管道位移沿轴向分布情况曲线(图 4),在中间节点201处横向位移达到最大值,为1.2 mm。在后续的非线性大变形稳定性分析中,将该位移曲线作为输入缺陷参数的基础,乘以一系列比例因子而获得最终的缺陷参数。

    表  1  管道-土壤系统参数
    下载: 导出CSV 
    | 显示表格
    图  3  管道几何参数及节点位置
    图  4  特征值法分析管道节点位移的分布曲线

    特征值屈曲分析得到管道临界载荷为0.181×1012 N,利用特征值分析得到的屈曲形态,作为管道模型的非线性分析的输入参数,模拟管道的初始缺陷。由于特征值分析中得到的位移值较小,通过比例因子进行放大。比例因子FI根据其值大小分为3组,每组选7个值,得到最终初始缺陷(表 2)。

    表  2  通过比例因子放大得到节点201的几何缺陷
    下载: 导出CSV 
    | 显示表格

    按照特征值方法分析的结果,中间节点201的位移值取得最大(图 5)。在大变形分析中,此节点所受的弯矩最大,位移也最大,因此,将该点的位移作为临界载荷判据。

    图  5  3组比例因子对应的中间节点201的位移-载荷曲线

    分析上述3组结果可以发现,比例因子越小,位移变化越剧烈。当比例因子为10时,位移-载荷曲线有明显拐点,出现在载荷为17.5×1010 N附近,表明该载荷处位移迅速增大,该点即为管道的屈曲载荷;当比例因子为20时,拐点也可以大致确定,但是其位置左移,在载荷为17×1010 N附近,说明此时屈曲载荷较比例因子为10的工况下降了约0.5×1010 N。随着比例因子的增加,初始缺陷也逐渐放大,位移曲线趋于平缓,无明显的拐点。

    该曲线单调增加,且形状下凹,表明该曲线各阶导数也单调增加,因此,不存在极值点。但是对该曲线应用拉格朗日中值定理,发现曲线中值点位置与位移拐点位置几乎重合,因此,将曲线的拉格朗日中值点确定为临界载荷(表 3图 6)。

    表  3  不同比例因子及其相应的临界载荷
    下载: 导出CSV 
    | 显示表格
    图  6  不同比例因子对应的临界载荷曲线

    经Matlab线性拟合临界载荷曲线,得到最终的屈曲载荷计算公式:

    式中:F为临界屈曲载荷,1010 N;X为特征值分析得到的欧拉屈曲载荷,取18.1×1010 N。

    通过建立有限元模型,对受轴向载荷的埋地管道进行非线性大变形分析,结果表明,初始缺陷与管道的屈曲载荷有较高的关联度。随着初始缺陷的不断增大,屈曲载荷不断降低。经过对临界载荷曲线的线性拟合,得到管道的屈曲临界载荷公式,管道的屈曲载荷与初始缺陷具有负线性相关性。

  • 图  1   含缺陷管道受力示意图

    图  2   土壤约束反力模型

    图  3   管道几何参数及节点位置

    图  4   特征值法分析管道节点位移的分布曲线

    图  5   3组比例因子对应的中间节点201的位移-载荷曲线

    图  6   不同比例因子对应的临界载荷曲线

    表  1   管道-土壤系统参数

    下载: 导出CSV

    表  2   通过比例因子放大得到节点201的几何缺陷

    下载: 导出CSV

    表  3   不同比例因子及其相应的临界载荷

    下载: 导出CSV
  • [1] 张宏. 基于应变的管道强度设计方法的适用性[J]. 油气储运, 2012, 31(12): 952-954. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201212026.htm

    ZHANG H. The applicability of strain-based pipeline strength design method[J]. Oil & Gas Storage and Transportation, 2012, 31(12): 952-954. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201212026.htm

    [2] 何瑞霞, 金会军, 吕兰芝, 等. 格尔木拉萨成品油管道沿线冻土工程和环境问题及其防治对策[J]. 冰川冻土, 2010, 32(1): 18-27. https://www.cnki.com.cn/Article/CJFDTOTAL-BCDT201001004.htm

    HE R X, JIN H J, LYU L Z, et al. Permafrost and environmental problems along the Golmud to Lhasa oil product and their mitigation[J]. Journal of Glaciology and Geocryology, 2010, 32(1): 18-27. https://www.cnki.com.cn/Article/CJFDTOTAL-BCDT201001004.htm

    [3] 刘啸奔, 张宏, 李勐, 等. 断层作用下埋地管道应变分析方法研究进展[J]. 油气储运, 2016, 35(8): 799-807. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201608001.htm

    LIU X B, ZHANG H, LI M, et al. Research progress of strain analysis methods for buried pipelines under faulting process[J]. Oil & Gas Storage and Transportation, 2016, 35(8): 799-807. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201608001.htm

    [4] 蔡克, 杨锋平, 罗恒金. 天然气管道的屈曲变形模拟[J]. 油气储运, 2013, 32(4): 402-405. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201304017.htm

    CAI K, YANG F P, LUO H J. Simulation on buckling deformation of natural gas pipeline[J]. Oil & Gas Storage and Transportation, 2013, 32(4): 402-405. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201304017.htm

    [5] 刘啸奔, 张宏, 夏梦莹. 冻土区管道上浮屈曲临界载荷试验研究[J]. 石油化工设备, 2014, 43(4): 11-14. https://www.cnki.com.cn/Article/CJFDTOTAL-SYSB201404003.htm

    LIU X B, ZHANG H, XIA M Y. Research on experimental method of upheaval buckling of pipeline in permafrost zone[J]. Petro-Chemical Equipment, 2014, 43(4): 11-14. https://www.cnki.com.cn/Article/CJFDTOTAL-SYSB201404003.htm

    [6] 姬贺炯, 白长青, 韩省亮. 输流管道动力有限元建模及实验研究[J]. 应用力学学报, 2013, 30(3): 422-427. https://www.cnki.com.cn/Article/CJFDTOTAL-YYLX201303022.htm

    JI H J, BAI C Q, HAN S L. Dynamic finite element modeling and experimental research of the fluid-filled pipeline[J]. Chinese Journal of Applied Mechanics, 2013, 30(3): 422-427. https://www.cnki.com.cn/Article/CJFDTOTAL-YYLX201303022.htm

    [7]

    WANG L Z, SHI R W, YUAN F, et al. Global buckling of pipelines in the vertical plane with a soft seabed[J]. Applied Ocean Research, 2011, 33(3): 130-136.

    [8] 王丽. 冻土区输油管道轴向稳定性分析[D]. 北京: 中国石油大学(北京), 2009: 70-78.

    WANG L. The axial stability analysis of the oil pipeline in permafrost areas[D]. Beijing: China University of Petroleum (Beijing), 2009: 70-78.

    [9] 余建星, 李智博, 杜尊峰, 等. 深海管道非线性屈曲理论计算方法[J]. 海洋工程, 2013, 31(1): 54-60. https://www.cnki.com.cn/Article/CJFDTOTAL-HYGC201301011.htm

    YU J X, LI Z B, DU Z F, et al. Theoretical calculation method of the nonlinear buckling of deep sea pipes[J]. The Ocean Engineering, 2013, 31(1): 54-60. https://www.cnki.com.cn/Article/CJFDTOTAL-HYGC201301011.htm

    [10] 李璞, 陶燕丽, 周建. 基于应变设计管道局部屈曲应变极限值的计算[J]. 天然气工业, 2013, 33(7): 101-107. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQG201307022.htm

    LI P, TAO Y L, ZHOU J. A study of the ultimate compressive strain of local buckling in strain-based design of pipelines[J]. Natural Gas Industry, 2013, 33(7): 101-107. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQG201307022.htm

    [11] 龚顺风, 陈源, 金伟良, 等. 高静水压力作用下深海油气管道的局部屈曲[J]. 浙江大学学报(工学版), 2012, 46(1): 14-19. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC201201004.htm

    GONG S F, CHEN Y, JIN W L, et al. Local buckling of deep water oil-gas pipeline under high hydrostatic pressure[J]. Journal of Zhejiang University, 2012, 46(1): 14-19. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC201201004.htm

    [12] 金浏, 王苏, 杜修力. 场地沉陷作用下埋地管道屈曲反应分析[J]. 世界地震工程, 2011, 27(2): 142-147. https://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201102022.htm

    JIN L, WANG S, DU X L. Buckling response analysis of buried pipelines subjected to the site soil settlement[J]. World Earthquake Engineering, 2011, 27(2): 142-147. https://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201102022.htm

    [13] 金浏, 李鸿晶. 逆冲断层作用下埋地管道屈曲分析[J]. 工程力学, 2011, 28(12): 98-104. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201112016.htm

    JIN L, LI H J. Buckling analysis of buried pipeline subject to reverse fault crossings[J]. Engineering Mechanics, 2011, 28(12): 98-104. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201112016.htm

    [14] 林镇诗, 余建星, 王永更, 等. 深水海底管道屈曲试验数据采集技术研究[J]. 石油矿场机械, 2011, 40(12): 62-66. https://www.cnki.com.cn/Article/CJFDTOTAL-SKJX201112017.htm

    LIN Z S, YU J X, WANG Y G, et al. Data acquisition technology of deep water submarine pipeline buckling experiments[J]. Oil Field Equipment, 2011, 40(12): 62-66. https://www.cnki.com.cn/Article/CJFDTOTAL-SKJX201112017.htm

    [15] 车小玉, 段梦兰, 曾霞光, 等. 双层管道整体屈曲实验研究及数值模拟[J]. 应用数学和力学, 2014, 35(2): 188-201. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSX201402007.htm

    CHE X Y, DUAN M L, ZENG X G, et al. Experimental study and numerical simulation of global buckling of pipe-in-pipe systems[J]. Applied Mathematics and Mechanics, 2014, 35(2): 188-201. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSX201402007.htm

    [16]

    ZENG X G, DUAN M L, CHE X. Critical upheaval buckling forces of pipelines[J]. Applied Ocean Research, 2014, 45(2): 33-39.

    [17] 刘羽霄, 葛涛, 李昕, 等. 初始几何缺陷在海底管道整体屈曲数值建模中的引入方法[J]. 水资源与水工程学报, 2011, 22(3): 32-35. https://www.cnki.com.cn/Article/CJFDTOTAL-XBSZ201103009.htm

    LIU Y X, GE T, LI X, et al. Introduction method of initial geometric imperfection in the numeral model of subsea pipeline global buckling[J]. Journal of Water Resources & Water Engineering, 2011, 22(3): 32-35. https://www.cnki.com.cn/Article/CJFDTOTAL-XBSZ201103009.htm

    [18] 刘羽霄, 葛涛, 李昕, 等. 管-土相互作用对海底管道横向屈曲的影响[J]. 石油机械, 2011, 39(8): 7-10. https://www.cnki.com.cn/Article/CJFDTOTAL-SYJI201108004.htm

    LIU Y X, GE T, LI X, et al. The effect of pipe-soil interaction on the lateral buckling of seafloor pipeline[J]. China Petroleum Machinery, 2011, 39(8): 7-10. https://www.cnki.com.cn/Article/CJFDTOTAL-SYJI201108004.htm

    [19] 张萌, 余建星, 孙震洲, 等. Patran与Abaqus技术在海底管道屈曲分析中的应用研究[J]. 海洋工程, 2016, 34(3): 55-62. https://www.cnki.com.cn/Article/CJFDTOTAL-HYGC201603007.htm

    ZHANG M, YU J X, SUN Z Z, et al. Application research of Patran and Abaqus in offshore pipeline buckling analysis[J]. The Ocean Engineering, 2016, 34(3): 55-62. https://www.cnki.com.cn/Article/CJFDTOTAL-HYGC201603007.htm

    [20] 闫宏生, 王铿, 李怀亮, 等. 海底管道临时管跨允许最大跨长的计算[J]. 油气储运, 2015, 34(9): 1015-1018. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201509020.htm

    YAN H S, WANG K, LI H L, et al. Calculation of the allowable maximum span of temporarily-suspended submarine pipeline section[J]. Oil & Gas Storage and Transpotation, 2015, 34(9): 1015-1018. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201509020.htm

    [21] YE H, JIN Z J, LIANG Y. Review of research on local buckling of subsea pipeline[J]. Ship Engineering, 2014, 36(3): 1-5. https://www.cnki.com.cn/Article/CJFDTOTAL-DGGC202303001.htm
图(6)  /  表(3)
计量
  • 文章访问数: 
  • HTML全文浏览量:  0
  • PDF下载量: 
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-08-12
  • 修回日期:  2017-02-05
  • 网络出版日期:  2023-08-20
  • 发布日期:  2017-01-03
  • 刊出日期:  2017-03-24

目录

/

返回文章
返回