李艷麗, 黃維平, 常 爽 , 柳振海,3
(1.上海船舶運輸科學(xué)研究所 a.航運技術(shù)交通行業(yè)重點實驗室;b.航運技術(shù)與安全國家重點實驗室,上海 200135;2.中國海洋大學(xué) 工程學(xué)院,山東 青島 266100;3.浙江省電力設(shè)計院,杭州 310012)
深水頂張緊式立管在深水油氣開發(fā)領(lǐng)域有著廣泛應(yīng)用,其結(jié)構(gòu)在形式上表現(xiàn)為內(nèi)層油管分別被內(nèi)套管和外套管包裹,套管用于提供隔離保護作用,套管之間有扶正器作為支撐,以避免發(fā)生碰撞事故損傷立管。目前國內(nèi)外學(xué)者已對多層立管進行大量研究。在雙套管立管方面:羅坤洪等[1]以兩端鉸支的多層耦合立管為研究對象,建立了多層立管耦合模型,運用MATLAB軟件對深水頂張式立管的全耦合動力響應(yīng)和參數(shù)敏感性進行了分析。在單套管立管方面:康莊等[2]運用Orcaflex軟件研究了立管頂部預(yù)張力對立管結(jié)構(gòu)強度的影響;張崎等[3]運用Abaqus軟件對單套管立管的非線性時域運動進行了研究;MAN等[4]和YUE等[5]提出了一個數(shù)學(xué)模型,研究了重力、壓力和熱膨脹等因素在立管縱向引起的應(yīng)變;ZHOU等[6]對單套管管中管模型的非線性應(yīng)力進行了研究;趙炳星[7]運用Abaqus立管分析軟件提供的間隙模型模擬了管中管模型,分別結(jié)合頂張力系數(shù)、內(nèi)外管間隙、海流和波浪等參數(shù)對立管進行了參數(shù)敏感性分析。
目前,在已有的深水頂張式立管設(shè)計研究中,多采用等效管(按照截面剛度等效)進行模擬分析,沒有考慮立管間復(fù)雜的耦合作用;同時,對雙套管立管的研究大多以兩端鉸支的多層耦合立管為研究對象,沒有考慮上部平臺和管間間隙對立管耦合運動的影響。在實際設(shè)計中,立管采用雙套管結(jié)構(gòu),且頂端與平臺鏈接,底端鉸支,因此有必要對該立管進行建模分析。
本文首先根據(jù)假設(shè)條件建立數(shù)學(xué)模型,其次給出考慮立管管間間隙的各層立管動力響應(yīng)運動方程的迭代求解過程,然后分析不同扶正器布置間距下,各層立管相同節(jié)點處的耦合動力響應(yīng),繪制各層立管相同節(jié)點處的相對位移與立管管間間隙的比值的時程圖,最后分析各層立管的彎曲應(yīng)力,繪制各層立管相同節(jié)點處的彎曲應(yīng)力比值時程圖。
建立立管模型(見圖1),其中:立管的頂端與平臺連接,底端固定于海底;僅考慮波浪載荷的作用,波浪方向為x方向;坐標(biāo)零點位于立管軸線底端中線處。
圖1 立管系統(tǒng)左側(cè)剖面模型
平臺運動對各層立管的動力響應(yīng)有很大影響,不可忽略。在Orcaflex軟件的輔助下,按照彎曲剛度等效的原則模擬計算等效管和平臺在環(huán)境載荷下的運動響應(yīng),有
EIeq=EIo+EIi+EIy
(1)
式(1)中:E為楊氏模量;Io、Ii和Iy分別為外套管、內(nèi)套管和油管的慣性矩;Ieq為等效管慣性矩。
平臺沿x軸的運動響應(yīng)與各層立管頂端節(jié)點的響應(yīng)一致,可將各層立管頂端節(jié)點的響應(yīng)變換為作用于各層立管頂端的頂張力和橫向作用力。
對于長細(xì)形頂張式立管,可忽略其剪切變形。因此,該類型立管的橫向彎曲振動問題可采用Euler-Bernoulli梁的復(fù)雜彎曲理論求解,各層立管均為均質(zhì)等截面管。沿x軸建立各層立管的有限元模型,采用newmark-β法的增量型式進行動力響應(yīng)分析,可得到基本形式的3層立管的運動控制方程為
(2)
管柱間隙的存在會對耦合方程產(chǎn)生影響,根據(jù)圖1所示模型,相鄰立管之間的相對位移為
(3)
(σt+Δt)2=(xi)t+Δt-(xy)t+h2
(4)
式(3)和式(4)中:(σt+Δt)1為外套管與內(nèi)套管在扶正器位置處的相對位移;(σt+Δt)2為內(nèi)套管與油管在扶正器位置處的相對位移;(xo)t+Δt為外套管在t+Δt時刻的位移;(xi)t+Δt為內(nèi)套管在t+Δt時刻的位移;h1為外套管與內(nèi)套管之間扶正器的高度;h2為內(nèi)套管與油管之間扶正器的高度;Di為內(nèi)套管的外徑;di為內(nèi)套管的內(nèi)徑。
在t時刻,外套管和內(nèi)套管在扶正器處的接觸條件為
(5)
式(5)中:上、下2個不等式分別為外套管和內(nèi)套管在軸線左側(cè)接觸和右側(cè)接觸時的判斷條件。
在t時刻,內(nèi)套管和油管在扶正器處的接觸條件為
(6)
式(6)中:Dy為油管的外徑;上、下2個不等式分別為內(nèi)套管和油管在軸線左側(cè)接觸和右側(cè)接觸時的判斷條件。
當(dāng)外套管與內(nèi)套管在扶正器s處的相對位移在t+Δt時刻滿足式(5),或內(nèi)套管與油管在扶正器s處的相對位移在t+Δt時刻滿足式(6)時,接觸節(jié)點s處的位移需通過迭代求得。
由Newmark-β法公式,有
(7)
(8)
(9)
式(7)~式(9)中:γ和β為與精度和穩(wěn)定性有關(guān)的參數(shù)。將式(2)變換為式(7)的形式,可得
(10)
圖2為立管系統(tǒng)左側(cè)剖面接觸示意。當(dāng)外套管與內(nèi)套管在t+Δt時刻的相對位移滿足式(5)時,從圖2a中可看出,外套管與內(nèi)套管在t+Δt時刻的相對位移滿足
圖2 立管系統(tǒng)左側(cè)剖面接觸示意
(11)
(12)
式(11)為在t時刻滿足式(5)時外套管與內(nèi)套管的位移關(guān)系;式(12)為在t時刻不滿足式(5)時外套管與內(nèi)套管的位移關(guān)系。
同理,當(dāng)內(nèi)套管與油管在t+Δt時刻的相對位移滿足式(6)時,從圖2b中可看出,外套管與內(nèi)套管在t+Δt時刻的相對位移滿足
(13)
(14)
式(13)為在t時刻滿足式(6)時內(nèi)套管與油管的位移關(guān)系;式(14)為在t時刻不滿足式(6)時內(nèi)套管與油管的位移關(guān)系。
迭代終止條件為
(15)
式(15)中:ε為迭代終止條件最大值。當(dāng)滿足式(15)時,停止迭代。迭代結(jié)束之后,即可計算各層立管在接觸點s處的速度和加速度增量,計算式為
(16)
(17)
計算立管其他節(jié)點i的位移增量、速度增量和加速度增量,有
(18)
(19)
(20)
(21)
計算模型如圖1所示。張力腿平臺參數(shù):立柱間距為86.25 m;工作水深為1 500 m;排水量為5.45×107kg;平臺縱搖轉(zhuǎn)動慣量為8.237×108kN·m2;平臺艏搖轉(zhuǎn)動慣量為9.807×108kN·m2;張力腿總預(yù)張力為1.373×105kN(頂張力按各層立管的軸向剛度比率分配)。管中管結(jié)構(gòu)主要參數(shù)見表1。采用JONSWAP譜模擬隨機海浪,特征波高Hs=8.3 m,峰值周期Tp=9.7 s。計算采用本文開發(fā)的Tube2D軟件,將立管模型劃分為1 500個單元,時間步長為0.1 s。
表1 管中管結(jié)構(gòu)主要參數(shù)
由于本文的主要研究對象為二維平面的頂張式立管,因此僅給出波浪入射角為0°時,張力腿平臺的縱蕩響應(yīng)時程圖、縱搖響應(yīng)時程圖、縱蕩速度時程圖和縱搖角速度時程圖,見圖3~圖6。
圖3 平臺縱蕩響應(yīng)時程圖
圖4 平臺縱搖響應(yīng)時程圖
圖5 平臺縱蕩速度時程圖
圖6 平臺縱搖角速度時程圖
圖7~圖10為管間扶正器的布置間距為25 m時,立管20 m和1 485 m處的響應(yīng)時程圖和相鄰立管的無量綱化相對位移曲線(其中Δx和p分別為相鄰2層立管的相對位移和間隙),頂張力系數(shù)為1.2。圖11~圖14為扶正器的布置間距為40 m時,立管20 m和1 485 m處的響應(yīng)時程圖和相鄰立管的無量綱化相對位移曲線,頂張力系數(shù)為1.2。
圖7 耦合立管20 m處響應(yīng)時程圖
圖8 耦合立管20 m處無量綱化相對位移曲線
圖9 耦合立管1 485 m處響應(yīng)時程圖
圖10 耦合立管1 485 m處無量綱化相對位移曲線
圖11 耦合立管20 m處響應(yīng)時程圖
圖12 耦合立管20 m處無量綱化相對位移曲線
圖13 耦合立管1 485 m處響應(yīng)時程圖
圖14 耦合立管1 485 m處無量綱化相對位移曲線
由圖7和圖11可知,扶正器的布置間距對立管頂端20 m處響應(yīng)幅值的影響不大。由圖8和圖12可知,當(dāng)扶正器的布置間距為40 m時,內(nèi)套管與油管20 m處的相對位移曲線幅值大于0且小于2,說明此時內(nèi)套管與油管管體發(fā)生碰撞。這與文獻[1]中的立管兩端鉸支分析有所不同,在文獻[1]中,當(dāng)扶正器間距為40 m時,外管與內(nèi)管在20 m處發(fā)生碰撞。造成這種差別的原因是立管上部運動受平臺運動的影響較大。
由圖9和圖13可知,扶正器的布置間距對立管1 485 m處響應(yīng)幅值的影響較小。由圖10和圖14可知,當(dāng)扶正器的布置間距增大時,立管間的相對位移有所減小,這是因為扶正器數(shù)量減少,其傳遞的載荷相對減小,平臺運動對立管底端的影響減小。
圖15為扶正器布置間距為25 m、頂張力系數(shù)為1.2時,立管85 m處彎曲應(yīng)力時程圖。圖16~圖18分別為各層立管彎曲應(yīng)力與總彎曲應(yīng)力的比值時程圖,分別繪制了外套管、內(nèi)套管和油管的彎曲應(yīng)力占三者總彎曲應(yīng)力的比例隨時間的變化情況,用以研究各層立管所受彎曲應(yīng)力是否按立管彎曲剛度之比分配,并分析其原因。
圖15 立管85 m處彎曲應(yīng)力時程圖
圖16 外套管彎曲應(yīng)力與總彎曲應(yīng)力的比值時程圖
圖17 內(nèi)套管彎曲應(yīng)力與總彎曲應(yīng)力的比值時程圖
圖18 油管彎曲應(yīng)力與總彎曲應(yīng)力的比值時程圖
從圖15中可看出,各層立管85 m節(jié)點處所受彎曲應(yīng)力具有隨機性,外套管所受彎曲應(yīng)力最大,內(nèi)套管次之,油管最小,這是因為外套管直接承受環(huán)境載荷,起到保護內(nèi)管的作用,而內(nèi)套管和油管僅受扶正器處相互作用力的作用,扶正器產(chǎn)生的相互作用力屬于被動力。
從圖16中可看出,外套管所受彎曲應(yīng)力與總彎曲應(yīng)力的比值在0.506附近波動;從圖17中可看出,內(nèi)套管所受彎曲應(yīng)力與總彎曲應(yīng)力的比值在0.332附近波動;從圖18中可看出,油管所受彎曲應(yīng)力占總彎曲應(yīng)力的比值在0.175附近波動;根據(jù)彎曲剛度等效分配原理分別求出內(nèi)套管、外套管和油管的彎曲剛度占比分別為0.772 7、0.207 5和0.019 8。
上述分析結(jié)果顯示,將彎曲應(yīng)力按彎曲剛度之比分配所得結(jié)果與實際計算的結(jié)果不符,理論值會高估外管承受的彎曲應(yīng)力,會低估內(nèi)套管和油管承受的彎曲應(yīng)力。因此,在實際設(shè)計立管過程中,應(yīng)根據(jù)耦合立管模型的計算結(jié)果對各層立管進行強度校核。
本文在文獻[1]的基礎(chǔ)上,對深水頂張式立管分析模型進行了改進。一方面,考慮了上部平臺運動對立管運動響應(yīng)的影響;另一方面,對立管管間間隙在計算方程中的處理方式進行了討論。通過對深水頂張式立管的運動響應(yīng)和彎曲應(yīng)力進行分析,主要得出以下結(jié)論:
1)在分析頂張式立管時,不能忽視扶正器和管間間隙對立管運動響應(yīng)的影響,扶正器的幾何參數(shù)、布置間距和管間間隙都會對立管耦合運動產(chǎn)生重要影響。與等效立管模型相比,頂張式立管耦合模型能更好地模擬實際立管模型。在對頂張式立管進行動力學(xué)研究時,應(yīng)充分考慮上部平臺的影響。
2)通過對耦合立管的彎曲應(yīng)力進行分析發(fā)現(xiàn),將彎曲應(yīng)力按彎曲剛度之比分配所得結(jié)果與實際計算的結(jié)果不符,理論值會高估外套管承受的彎曲應(yīng)力,會低估內(nèi)套管與油管承受的彎曲應(yīng)力。因此,在實際設(shè)計立管過程中,應(yīng)根據(jù)耦合立管模型的計算結(jié)果對各層立管進行強度校核,對今后的立管設(shè)計有一定的參考意義。