大庆石油学院学报第35卷第3期2011年6月JOU RNAL OF DAQING PE TRO LE UM INS TIT UT E V o l.35No.3Jun.2011
刘羽霄1，葛 涛2，李 昕3，周 晶3
(1.山东工商学院管理科学与工程学院，山东烟台 264005; 2.山东工商学院工商管理学院，山东烟台 264005;
3.大连理工大学海岸与近海工程国家重点实验室，辽宁大连 116023)
摘 要：为了研究初始几何缺陷对海底管道前、后屈曲的影响，基于ANS YS软件建立平坦海床上裸铺管道的非线性
有限元模型，分析建模过程中初始几何缺陷的引入方法，得到缺陷幅值及缺陷长度对管道前屈曲及后屈曲的影响规律.
结果表明：临界温差随着初始几何缺陷幅值的增加而降低，随着缺陷段长度的增加呈现先降低后增加的趋势；缺陷幅值
及缺陷段长度对管道后屈曲的变形、弯矩及轴向应变影响不大.这对海洋管道工程设计提供指导.
关 键 词：海底管道;ANS YS软件；横向屈曲；引入方法；缺陷幅值；缺陷长度
中图分类号:TE973.11 文献标识码:A 文章编号:1000-1891(2011)03-0076-05
0 引言
随着深水和边际海洋油气田的开发，越来越多高温/高压海底管道应用于工程实际.海底管道在内压和温度等工作荷载作用下，易发生整体屈曲.对于埋设在海床以下的管道，易发生垂向隆起屈曲；对于直接铺设在海床上的管道，易发生水平横向屈曲.
H obbs R E给出理想直管道隆起及横向屈曲临界荷载的计算公式，其研究基于小坡角、线弹性假设，且没有考虑初始几何缺陷的影响[1].Taylo r N等对海底管道的隆起屈曲进行理论和实验研究，针对3种不同的缺陷，分别给出管道的临界屈曲荷载，并通过实验结果对理论研究进行验证[2].James G A C采用简化模型研究隆起热屈曲，考虑初始几何缺陷及残余应力影响，给出管道隆起屈曲段的临界荷载的计算公式[3].刑静忠等研究部分悬跨管道的悬跨段和埋地段的临界屈曲荷载和变形[4].刘润等基于编制的软件探讨地基土摩阻力、管道上覆土厚度对埋地管道安全温度的影响及膨胀弯的作用[5].贾旭等研究海底管道在静水压力下的局部屈曲问题[6].苗同臣等分析初始弯曲变形对A型石油井架稳定性的影响[7].
整体屈曲理论公式仅用于判断初始几何缺陷对管道临界屈曲荷载的影响，无法用于判断初始几何缺陷对管道后屈曲的变形、应变及弯矩的影响.对于初始几何缺陷对管道后屈曲的影响分析需要采用有限元方法，如基于有限元分析软件ANSYS、ABAQUS等建立海底管道整体屈曲分析的有限元模型.初始几何缺陷是海底管道整体屈曲分析中的关键因素，基于有限元分析软件建立的海底管道整体屈曲分析的有限元模型，必须引入几何缺陷才能进行整体屈曲分析.关于海底管道横向屈曲问题的相关研究较少，基于ANSYS软件建立平坦海床上裸铺管道的非线性有限元模型，分析建模过程中初始几何缺陷的引入方法，以及初始几何缺陷幅值及缺陷长度对管道前屈曲及后屈曲的影响，以对海洋管道工程设计提供指导.
1 裸铺管道非线性有限元模型
1.1 管道模型
为了研究管道材料非线性的影响，选择PIPE20单元模拟管道.整体屈曲分析主要研究屈曲段的变形，因此管道屈曲段单元长度取1倍管道外直径，屈曲段之外单元长度逐渐增加.
管道钢的材料应力应变关系采用Ramberg-Osg ood本构关系[8-9]，可以表示为
收稿日期:2010-12-15；审稿人：罗 敏；编辑：任志平
基金项目：国家自然科学基金重点项目(50439010)；教育部科学技术研究重大项目(305003)
作者简介：刘羽霄(1978-)，女，博士，讲师，主要从事海底管线破坏机理方面的研究.
ε(σ)=σE
1+37σσR n +1,
十二五科技成就展>中华研修大学(1)
n =lg
ε2-σ2
E ε1-σ1
E
lg σ2σ1,
(2)
σR =
σ173ε1E σ1
-11
法国国庆日
n +1
,(3)
图1 X65号钢应力-应变关系
式中:n 为材料的硬化系数;σR 为Ramberg -Osgo od 应力；下标1和2分别代表屈服强度和极限抗拉强度对
应的应力-应变点;E 为弹性模量.X65号钢的应力-应变关系见图1.1.2 土体模型
由于土体的高度非线性及管道易变形，管-土之间的相互作用复杂.所建立的非线性有限元模型中，土体采用弹簧单元模拟.对于平坦海床上裸铺管道的横向屈曲，因管道的屈曲主要在水平面内发生，因此只需要建立二维平面模型.为了模拟土体的非线性摩擦特性，选用COM BIN39号弹簧单元.
海床与管道之间的非线性弹塑性的轴向/横向相互作用力采用弹簧单元模拟，土弹簧单元的力-位移关系见图2.其中:F ly 为屈服力,Δly 为屈服位移.海床与管道之间相互作用的轴向/横向摩擦力采用库仑定律计算
[10]
:
F =μW ,
(4)
式中:F 为土弹簧对应的屈服力;μ为管道与海床之间的摩擦因数;W 为单位长度管道的有效重力.1.3 管道与海床的相互作用模型
基于ANS YS
软件，建立平坦海床上裸铺管道的非线性有限元模型[9-12]
，管-土相互作用模型见图3.主要研究管道横向的屈曲，模型中约束管道竖向的运动.根据实际模拟情况，管道两端采用完全约束或自由边界条件，主要考虑压力及温度荷载.引入初始几何缺陷后的管道有限元剖分见图4.
图2 土弹簧的力-位移关系图3 管-土相互作用模型
第3期 刘羽霄等：初始几何缺陷对海底管道横向屈曲的影响
2 初始几何缺陷
2.1 形状
数值分析中采用管道初始几何缺陷的形状[3,11]见图5.
图4 管道有限元剖分图5 管道初始几何缺陷形状
y 0=w 012+12cos 2πx L 0
,
(5)
式中:y 0为缺陷段相对于x 轴的横向变形;w 0为初始几何缺陷幅值;L 0为初始几何缺陷长度.
屈曲幅值与屈曲段长度之间的关系[2]为
w m /L 4
=2.407×10-3
μl W /EI =w 0/L 4
0,
(6)即
w 0=μl W L 4
0415EI
,
(7)
式中:μl 为管道与海床间的横向摩擦因数;w m 为屈曲幅值,w m =y
x =0.2.2 引入方法
2.2.1 建模引入
国宝沉浮录
图6 含初始几何缺陷的管道模型
此种方法在建立管道有限元模型时引入初始缺陷，即管道初始轴向力为零[11]
(见图6).在计算过程中，忽略管道形成初始几何缺陷时在管道中产生的残余应力的影响.此种方法适用于在铺设管道时人为设置缺陷，且管道设置缺陷后其轴向力为零的情况.2.2.2 加载引入
此种方法在加载过程中引入初始几何缺陷，考虑管道形成初始缺陷时引起的残余应力影响.它适用于图7 管道非线性有限元分析流程管道在铺设时为直管道，在铺设过程中或在管道运行过程中形成随机几何缺陷的情况.加载及引入缺陷的
流程[9]
:
(1)施加管道淹没重力、静水压力和有效残余力.(2)通过移动管道节点和弹簧节点施加初始缺陷.(3)释放缺陷段管道节点的约束.(4)释放管末端.在步骤(3)之后，步骤(4)之前，约束轴向弹簧.
(5)施加内压荷载.(6)施加温度荷载.2.3 横向屈曲数值分析
结合初始几何缺陷的引入方法，给出管道非线性有限元分析流程(见图7).2.4 结果分析2.4.1 缺陷幅值
缺陷是引发管道整体屈曲的关键因素，这里所取缺陷长度为100m ，缺陷幅值分别为0.8,1.0,1.2，大 庆 石 油 学 院 学 报 第35卷 2011年
1.5m .数值分析时管道相关参数见表1.
表1 管道相关参数
管道直径/mm
管道壁厚/mm
热膨胀因数管道淹没重力/(N ·m
-
1)
屈服应力/M Pa
极限强度/M Pa
3001411.7×10-6900448550横向摩擦因数
轴向摩擦因数
管道长度/m 运行压力/M Pa
运行温度/℃
环境温度/℃
0.75
0.5
2800
20
95
图8 不同初始缺陷幅值时的温差-最大横向位
移曲线
用户环境
不同初始几何缺陷幅值时的温差-最大横向位移曲线见图8.由图8可见，初始几何缺陷幅值对管道前屈曲临界温差有明显影响，临界温差随着初始几何缺陷幅值的增加而降低.
达到目标温度时不同初始缺陷幅值对管道后屈曲响应的影响曲线见图9.由图9可见，初始几何缺陷幅值对管道后屈曲的变形、弯矩、轴向应变影响不大.2.4.2
缺陷长度
缺陷幅值为
0.8m ，初始几何缺陷长度为52,64,80,100,110m 时临界温差变化曲线见10.由图10可见，临界温差随着缺陷长度的增加呈现先降低后增加的趋势.
图9 不同初始几何缺陷幅值对后屈曲响应的影响曲线
图10 不同初始几何缺陷长度时的临界温差变化曲线
达到目标温度时初始几何缺陷长度对管道后屈曲响应的影响曲线见图11.由图11可见，在缺陷幅值一定时，改变缺陷长度对管道后屈曲的变形、弯矩及轴向应变影响不大.
第3期 刘羽霄等：初始几何缺陷对海底管道横向屈曲的影响
大 庆 石 油 学 院 学 报 第35卷 2011年
图11 不同初始缺长度对后屈曲响应的影响
3 结论
迪布韦克
(1)基于ANSYS软件建立平坦海床上裸铺管道的非线性有限元模型，分析建模过程中初始几何缺陷的引入方法，以及初始几何缺陷幅值及缺陷长度对管道前屈曲及后屈曲响应的影响.
(2)改变初始几何缺陷幅值对管道前屈曲临界温差有明显影响；临界温差随着初始几何缺陷幅值的增加而降低，对管道后屈曲的变形、弯矩、轴向应变影响不大.
(3)临界温差随着缺陷段长度的增加呈现先降低后增加的趋势；在缺陷幅值一定时，改变缺陷段长度对管道后屈曲的变形、弯矩及轴向应变影响不大.
参考文献:
[1] Hobb s R E.In-service buckling of heated pipeline[J].J.Transport Eng,1984,110(2):175-189.
[2] Taylor N,T ran V.Experimental and theoretical studies in subsea pipeline b uckling[J].M arine Structu res,1996,9(2):211-257.
[3] J ames G A C.A sim plified model of u pheaval thermal buck ling of sub sea pipelines[J].T hin-Walled S tructures,1997,29:59-78.
[4] 刑静忠，柳春图，徐永君.埋设悬跨海底管道的屈曲分析[J].工程力学,2006,23(2):173-162.
[5] 刘润，闫澍旺，孙国民.温度应力下海底管线屈曲分析方法的改进[J].天津大学学报,2005,38(2):124-128.
[6] 贾旭，侯静，田英辉.静水压力作用下海底单层保温管道压溃屈曲分析[J].中国海上油气,2006,15(5):342-343.
[7] 苗同臣，张艳艳，赖永星，等.初始弯曲变形对A型石油井架稳定性的影响[J].大庆石油学院学报,2008,32(3):59-60.
[8] Ramberg W,Osgood W R.Description of s tress-s train curves by th ree parameters[C].NAC A Technical Note,1943.
[9] S ren H,Bai Y.Use of finite elemen t analysis for local bu cklin g design of pipelines[C].Lisbon,Proceedings of the17th International
C onferen ce on Offshore M ech anics and Arctic Engineering,1998:1-9.
[10] Lars C.Dis placement control in lateral buckling of“sh ort”pipelines[C].S eoul,Proceeding s of th e15th International Offshore and
Polar Engineering Conference,2005:84-92.
[11] Nader Yoosef-Ghodsi.Analy sis of bu ried pipelines w ith therm al applications[D].Canada:University of Alb erta,2002.
[12] Enrico T,Luigino V,Erik L.S naking of sub marine pipelines res tin g on flat sea bottom u s ing finite elemen t method[C].Brest,Pro-
ceedings of the9th International Offs hore an d Polar E ngineering Conference,1999,2:34-44.
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议