，

(1.江蘇科技大學(xué) 船舶與海洋工程學(xué)院， 江蘇 鎮(zhèn)江 212003； 2. 上海外高橋造船有限公司， 上海 200137 )

0 引 言

導(dǎo)管架平臺(tái)是近些年來(lái)被廣泛應(yīng)用于近海海洋石油開(kāi)采活動(dòng)的一種固定式平臺(tái)，具有適應(yīng)性較強(qiáng)、結(jié)構(gòu)安全可靠、結(jié)構(gòu)簡(jiǎn)單等特點(diǎn)。但是，導(dǎo)管架平臺(tái)長(zhǎng)期處于惡劣的海洋環(huán)境中，不斷承受風(fēng)、浪、流的交替作用或聯(lián)合作用，有時(shí)還承受冰載荷、腐蝕、沖刷、疲勞和低溫的影響，在這些載荷長(zhǎng)期作用下，局部或整體結(jié)構(gòu)的強(qiáng)度或抵抗能力會(huì)下降。導(dǎo)管架結(jié)構(gòu)如果強(qiáng)度不足，導(dǎo)致平臺(tái)倒塌，會(huì)帶來(lái)嚴(yán)重的安全威脅及經(jīng)濟(jì)損失。因此，對(duì)在波浪載荷作用下的導(dǎo)管架平臺(tái)的研究具有重要的意義。

嵇春艷等[1]利用ANSYS分析導(dǎo)管架平臺(tái)在波浪和海流作用下的結(jié)構(gòu)內(nèi)力響應(yīng)，獲得平臺(tái)的頻率特性和動(dòng)力特性。吳家鳴等[2]通過(guò)數(shù)值分析平臺(tái)所有構(gòu)件的最大應(yīng)力分布以及樁基的位移特征，得出主樁的最大應(yīng)力區(qū)域主要發(fā)生在主樁與海底泥面相接的節(jié)點(diǎn)附近，然后對(duì)平臺(tái)結(jié)構(gòu)在特性海況下的安全性進(jìn)行觀測(cè)及研究。黃懷州等[3]研究導(dǎo)管架平臺(tái)結(jié)構(gòu)在波浪載荷作用下的疲勞可靠性，發(fā)現(xiàn)一階模態(tài)對(duì)深水導(dǎo)管架平臺(tái)結(jié)構(gòu)響應(yīng)起主要作用，并求得結(jié)構(gòu)在隨機(jī)應(yīng)力譜下給定疲勞壽命時(shí)的疲勞可靠性指標(biāo)。GAIDAI等[4]和NAESS等[5]基于蒙特卡羅的方法研究隨機(jī)波浪載荷作用下近海固定結(jié)構(gòu)物的極限響應(yīng)問(wèn)題，該方法可以快速、簡(jiǎn)單、有效地預(yù)報(bào)短期及長(zhǎng)期平臺(tái)的極限響應(yīng)歷程，對(duì)于近海平臺(tái)結(jié)構(gòu)的安全性預(yù)測(cè)具有重要作用。

目前對(duì)于導(dǎo)管架平臺(tái)的結(jié)構(gòu)響應(yīng)研究都偏向規(guī)則波，而對(duì)隨機(jī)波的研究則較少。因此，本文選取渤海灣塘沽地區(qū)某導(dǎo)管架平臺(tái)，結(jié)合隨機(jī)波浪理論及Morison方程，研究導(dǎo)管架平臺(tái)在無(wú)流狀態(tài)時(shí)8個(gè)不同浪向的隨機(jī)波浪載荷作用下的結(jié)構(gòu)響應(yīng)，對(duì)于今后平臺(tái)結(jié)構(gòu)設(shè)計(jì)及加強(qiáng)具有重要的意義。

1 波浪頻譜及非線性波浪載荷計(jì)算

1.1 JONSWAP譜

通過(guò)給定的有義波高和峰值周期，利用MATLAB得到JONSWAP譜[6]，即譜密度與圓頻率的關(guān)系為

(1)

1.2 導(dǎo)管架運(yùn)動(dòng)方程

利用有限元法建立導(dǎo)管架平臺(tái)動(dòng)力響應(yīng)模型，該模型假設(shè)水動(dòng)力作用在離散成節(jié)點(diǎn)的平臺(tái)上，作用區(qū)間為甲板到海床之間，該運(yùn)動(dòng)方程為

(2)

1.3 Morison方程

Morison方程是工程上常用的計(jì)算小尺度樁柱上波浪力的計(jì)算方法，具體形式為

(3)

2 導(dǎo)管架平臺(tái)響應(yīng)分析

2.1 導(dǎo)管架平臺(tái)模型概況

本文選取渤海灣塘沽地區(qū)實(shí)際運(yùn)營(yíng)導(dǎo)管架平臺(tái)，該平臺(tái)已經(jīng)運(yùn)營(yíng)數(shù)十年，考慮到運(yùn)營(yíng)安全性和可靠性，有必要對(duì)其進(jìn)行結(jié)構(gòu)分析。

圖1 導(dǎo)管架平臺(tái)有限元模型

利用ANSYS有限元軟件建立多自由度三維導(dǎo)管架平臺(tái)結(jié)構(gòu)模型。該平臺(tái)由4根樁腿和主甲板組成，不考慮波浪對(duì)甲板結(jié)構(gòu)的影響，對(duì)導(dǎo)管架平臺(tái)結(jié)構(gòu)進(jìn)行一定的簡(jiǎn)化，甲板結(jié)構(gòu)利用超單元建模。導(dǎo)管架平臺(tái)如圖1所示。該海域平均水深為17 m，導(dǎo)管架平臺(tái)海底泥線下的樁腿長(zhǎng)度是90 m，平均水位(Mean Water Level,MWL)距底甲板的距離為12 m。本文導(dǎo)管架圓管的許用應(yīng)力為355 MPa。

2.2 不同方向的波浪工況

由于平臺(tái)的截面為矩形，本文考慮8個(gè)不同方向隨機(jī)波浪對(duì)導(dǎo)管架平臺(tái)的作用。查閱相關(guān)資料得到渤海地區(qū)全年不同浪向的隨機(jī)波浪概率分布見(jiàn)表1。

表1 不同浪向的隨機(jī)波浪概率分布

為對(duì)導(dǎo)管架平臺(tái)結(jié)構(gòu)強(qiáng)度進(jìn)行分析，根據(jù)渤海灣地區(qū)波浪的分布情況，分別從每個(gè)浪向中選取1個(gè)經(jīng)典海況，并且該海況中的隨機(jī)波浪為出現(xiàn)次數(shù)最多的工況。 不同工況的波浪參數(shù)具體數(shù)值見(jiàn)表2。

表2 不同工況的波浪參數(shù)

圖2 JONSWAP譜

2.3 數(shù)值仿真及分析

2.3.1 敏感性測(cè)點(diǎn)的選取

根據(jù)渤海灣海域的具體情況，海浪譜選取JONSWAP譜，由表2所示的波浪基本參數(shù)，首先選擇工況8在無(wú)流狀態(tài)下進(jìn)行分析，利用MATLAB軟件對(duì)表2中的8種工況進(jìn)行編程運(yùn)算，輸入有義波高和譜峰周期生成JONSWAP譜，然后分離出多個(gè)規(guī)則波，用來(lái)模擬不規(guī)則的海況。JONSWAP譜如圖2所示。

將MATLAB生成的數(shù)據(jù)導(dǎo)入ANSYS進(jìn)行計(jì)算，得到工況8的數(shù)值仿真結(jié)果。圖3為工況8某一時(shí)刻導(dǎo)管架平臺(tái)整體應(yīng)力和位移變形特征，可以看出：該時(shí)刻導(dǎo)管架平臺(tái)的最大應(yīng)力值達(dá)到54 MPa，最大位移0.49 mm；應(yīng)力的集中區(qū)域主要在樁腿與泥面的交匯處。選取該危險(xiǎn)區(qū)域的2個(gè)敏感性單元以及1個(gè)節(jié)點(diǎn)進(jìn)行詳細(xì)分析，測(cè)點(diǎn)的分布如圖4所示，編號(hào)1～3分別表示Node 2953、Elem 44173、Elem 77150。

圖3 工況8整體結(jié)構(gòu)應(yīng)力和變形云圖

圖4 測(cè)點(diǎn)分布圖

2.3.2 不同浪向的結(jié)構(gòu)響應(yīng)分析

通過(guò)對(duì)工況8的計(jì)算，確定了危險(xiǎn)區(qū)域，為研究該危險(xiǎn)區(qū)域在不同浪向的隨機(jī)波浪作用下的結(jié)構(gòu)響應(yīng)，在該危險(xiǎn)區(qū)域選擇2個(gè)敏感單元進(jìn)行詳細(xì)分析。利用ANSYS進(jìn)行數(shù)值仿真瞬態(tài)分析，得到工況8下測(cè)點(diǎn)2和測(cè)點(diǎn)3的應(yīng)力隨時(shí)間的變化情況，如圖5所示。測(cè)點(diǎn)1的x方向和y方向的位移如圖6所示。

根據(jù)工況8仿真結(jié)果，可以得到測(cè)點(diǎn)2的最大應(yīng)力值為45.029 MPa，最小應(yīng)力值為39.723 MPa，應(yīng)力均值為42.480 MPa。該測(cè)點(diǎn)的應(yīng)力在許用范圍內(nèi)，結(jié)構(gòu)滿足強(qiáng)度要求，但是該組應(yīng)力的標(biāo)準(zhǔn)差卻達(dá)到了0.86 MPa，由此可見(jiàn)該工況對(duì)于結(jié)構(gòu)的應(yīng)力響應(yīng)具有較大的影響，結(jié)構(gòu)的應(yīng)力具有較大的波動(dòng)性，較不穩(wěn)定。同樣的情況，測(cè)點(diǎn)3的應(yīng)力也在許用范圍內(nèi)，但是應(yīng)力值變化幅度較大。在該工況下，測(cè)點(diǎn)1的x方向最大變形為0.33 mm，y方向最大變形為2.52 mm。因此，在工況8即浪向?yàn)楸?、無(wú)流狀態(tài)下，導(dǎo)管架平臺(tái)結(jié)構(gòu)的應(yīng)力變化幅度較大，應(yīng)是重點(diǎn)關(guān)注的海況。導(dǎo)管架平臺(tái)在不同工況下各測(cè)點(diǎn)的應(yīng)力統(tǒng)計(jì)見(jiàn)表3～表5。

圖5 應(yīng)力時(shí)歷圖

圖6 測(cè)點(diǎn)1位移

表3 在不同工況下測(cè)點(diǎn)2的應(yīng)力統(tǒng)計(jì) MPa

表4 在不同工況下測(cè)點(diǎn)3的應(yīng)力統(tǒng)計(jì) MPa

表5 在不同工況下測(cè)點(diǎn)1的位移統(tǒng)計(jì) mm

從表3～表5可以看出：在8種不同工況下，測(cè)點(diǎn)2的應(yīng)力值均達(dá)到40 MPa，應(yīng)力平均值約為42 MPa，所以測(cè)點(diǎn)2的應(yīng)力值均在許用應(yīng)力范圍內(nèi)，結(jié)構(gòu)的強(qiáng)度滿足要求；測(cè)點(diǎn)2在工況1中應(yīng)力最大，但應(yīng)力隨時(shí)間的變化卻十分地平緩，此外工況4、5、7、8的應(yīng)力也相對(duì)較高；測(cè)點(diǎn)3的應(yīng)力均值都在37 MPa左右，也滿足強(qiáng)度要求；在工況7即浪向?yàn)槲鞅狈较驎r(shí)，應(yīng)力達(dá)到最大值40.381 MPa，且應(yīng)力標(biāo)準(zhǔn)差也達(dá)到最大值0.858，雖然測(cè)點(diǎn)3在工況8下，應(yīng)力不是最大的，但是應(yīng)力標(biāo)準(zhǔn)差也達(dá)到0.781，說(shuō)明工況7和工況8對(duì)測(cè)點(diǎn)3的應(yīng)力影響較大。

導(dǎo)管架平臺(tái)在8個(gè)不同工況作用下，位移變化不大，x方向位移最大值為1.44 mm，y方向最大位移為2.52 mm。從圖中可以看出工況4和工況8，測(cè)點(diǎn)3的y方向位移較大，因?yàn)椴ɡ说姆较蛘醚貀軸；同樣地，工況2和工況6下x方向位移較大，說(shuō)明這是由于浪向角的變化，引起x和y方向分量的變化。另外，由于本文選取的海況波高不大，所以結(jié)構(gòu)的位移不明顯。

綜合各種工況，導(dǎo)管架平臺(tái)在給定工況的隨機(jī)波浪載荷作用下結(jié)構(gòu)響應(yīng)比較大，在有義波高為1.51 m，周期為6.9 s(即工況8)時(shí)，敏感性測(cè)點(diǎn)的應(yīng)力變化幅度最大，應(yīng)力值也較大，因此該工況為最危險(xiǎn)的工況。為進(jìn)一步研究該危險(xiǎn)工況對(duì)導(dǎo)管架平臺(tái)的影響，本文還考慮了流速的影響。

2.3.3 波流聯(lián)合作用下結(jié)構(gòu)響應(yīng)分析

本節(jié)考慮流速分別為0 m/s、0.5 m/s和1 m/s 的情況，流速的方向?yàn)闁|南方向，在有義波高為1.51 m，周期為6.9 s(即工況8)時(shí)，針對(duì)危險(xiǎn)區(qū)的敏感性測(cè)點(diǎn)2和3進(jìn)行響應(yīng)分析，并測(cè)得測(cè)點(diǎn)1的位移變化情況，具體數(shù)值如表6～表8所示。

表6 測(cè)點(diǎn)2在不同流速下應(yīng)力統(tǒng)計(jì)

表7 測(cè)點(diǎn)3在不同流速下應(yīng)力統(tǒng)計(jì)

表8 測(cè)點(diǎn)1在不同流速下位移統(tǒng)計(jì)

由表6～表8可以看出：隨著流速的增加，結(jié)構(gòu)的位移量逐漸減小，這是因?yàn)榱魉俚姆较蚺c波浪的方向相反，使得x和y方向的分量逐漸變??；導(dǎo)管架平臺(tái)在波流聯(lián)合作用下，結(jié)構(gòu)的應(yīng)力明顯增大，隨著流速的不斷增加，應(yīng)力也逐漸增加，且應(yīng)力標(biāo)準(zhǔn)差也不斷增大，表明結(jié)構(gòu)的應(yīng)力幅值增大，波流聯(lián)合作用對(duì)導(dǎo)管架平臺(tái)的結(jié)構(gòu)響應(yīng)具有很大的影響。

3 結(jié) 論

(1) 導(dǎo)管架平臺(tái)的危險(xiǎn)區(qū)域主要在樁腿與泥面的交匯處。

(2) 平臺(tái)有義波高為1.51 m，周期為6.9 s，浪向?yàn)?0°時(shí)，結(jié)構(gòu)的應(yīng)力變化幅度較大，屬于危險(xiǎn)工況，值得重點(diǎn)關(guān)注。

(3) 導(dǎo)管架平臺(tái)在波流聯(lián)合作用下，結(jié)構(gòu)的應(yīng)力明顯增大，隨著流速的不斷增加，應(yīng)力也逐漸增加，結(jié)構(gòu)的應(yīng)力變化幅度也增大，因此波流聯(lián)合作用對(duì)導(dǎo)管架平臺(tái)的結(jié)構(gòu)響應(yīng)具有很大的影響。

[1] 嵇春艷, 劉聰. 隨機(jī)波浪載荷作用下深水自升式海洋平臺(tái)動(dòng)力響應(yīng)分析[J]. 造船技術(shù), 2012(05): 27-31.

[2] 吳家鳴, 劉昊宇. 不同浪向的波浪載荷作用下導(dǎo)管架平臺(tái)主樁強(qiáng)度分析[J]. 海洋技術(shù)學(xué)報(bào), 2014, 33(03): 85-93.

[3] 黃懷州, 洪明, 遲少艷. 隨機(jī)波浪載荷作用下導(dǎo)管架平臺(tái)動(dòng)力響應(yīng)及疲勞可靠性分析[J]. 船舶力學(xué), 2006, 10(04): 65-71.

[4] GAIDAI O, NAESS A. Extreme Response Statistics for Drag Dominated Offshore Structures[J]. Probabilistic Engineering Mechanics, 2008(23): 180-187.

[5] NAESS A, GAIDAI O. Monte Carlo Methods for Estimating the Extreme Response of Dynamical Systems[J]. Journal Engineering Mechanics, 2008.

[6] 付昱華. 有效波高、平均周期和峰頻率確定的JONSWA P譜[J]. 中國(guó)海上油氣(工程), 1994, 02(01): 45-47.

[7] WANG Z F, WU K J, ZHOU L M, et al. Wave Characteristics and Extreme Parameters in the Bohai Sea[J]. China Ocean Engineering, 2012, 26(02)：341-350.