馬網(wǎng)扣， 袁飛暉

(1.中國船舶及海洋工程設(shè)計研究院，上海 200011；2.上海外高橋造船有限公司，上海 200137)

流速對自升式平臺動力放大因子的影響研究

馬網(wǎng)扣1， 袁飛暉2

(1.中國船舶及海洋工程設(shè)計研究院，上海 200011；2.上海外高橋造船有限公司，上海 200137)

研究自升式平臺動力放大因子的隨機波動態(tài)響應(yīng)分析法，以某自升式鉆井平臺為例進(jìn)行動力放大因子分析，結(jié)果顯示流速越大，流對平臺動態(tài)和靜態(tài)運動響應(yīng)的抑制越顯著，平臺動力放大因子隨著流速增大而減小。進(jìn)一步進(jìn)行算例平臺的慣性力分析，結(jié)果顯示隨機波動態(tài)響應(yīng)分析法能夠考慮流速的影響，該方法得到的慣性力增大相對平緩。研究表明，流速較大時使用隨機波動態(tài)響應(yīng)分析法確定自升式平臺的動力放大因子和慣性力，能優(yōu)化樁腿設(shè)計，避免無謂的結(jié)構(gòu)加強，有利于增大平臺的可變載荷設(shè)計指標(biāo)。

動力放大因子；隨機波動態(tài)響應(yīng)分析法；流速；慣性力

0 引 言

自升式平臺的作業(yè)特點要求其能在站立狀態(tài)下抵抗一定的風(fēng)暴。風(fēng)暴載荷是計算校核自升式平臺的樁腿強度、傾覆穩(wěn)性和滑移穩(wěn)性的關(guān)鍵參數(shù)。風(fēng)暴載荷組成中波浪激振作用引起的結(jié)構(gòu)慣性力不可忽略，而慣性力大小主要取決于平臺的動力放大因子(DAF)。隨著作業(yè)水深的加大，自升式平臺的慣性力在整體載荷中的比重也越來越大。自升式平臺的動力放大因子可理解為波激作用下平臺的動力響應(yīng)極值與靜態(tài)響應(yīng)極值之比，是反映平臺結(jié)構(gòu)動力放大效應(yīng)的參數(shù)。目前業(yè)界對自升式平臺動力放大因子的研究主要集中在平臺動態(tài)響應(yīng)極值計算、波浪載荷線性化、結(jié)構(gòu)模型簡化、樁基對動力放大因子的影響等幾個方面[1—4]。Zhang等[5]通過對兩個自升式平臺的動態(tài)響應(yīng)極值和動力放大因子算例進(jìn)行比較，討論了4種常用的動力放大因子的動態(tài)響應(yīng)分析法之間的差異，并研究了生成隨機波的隨機種子對動力放大因子計算結(jié)果的敏感性。黃盟宇等[6]通過Stokes五階波方法計算了某自升式平臺的靜態(tài)和動態(tài)基礎(chǔ)剪切力，討論了單自由度法計算動力放大因子的適用性，還比較了樁腿與船體采用齒輪齒條連接和鎖緊機構(gòu)連接間平臺動力放大因子的差異，指出增加平臺剛度能對其動力放大因子起到一定的抑制作用。Saeidtehrani等[7]選擇波斯灣典型海底土質(zhì)，比較了5種不同的樁基條件對自升式平臺動力放大因子的影響。Perry等[8]研究指出當(dāng)使用拖曳/慣性參數(shù)法計算自升式平臺動力放大因子時，應(yīng)將結(jié)果乘以一個與周期比有關(guān)的放大系數(shù)，以免在平臺固有周期與波浪激勵周期接近時低估其動力放大因子。

風(fēng)暴自存工況是自升式平臺的關(guān)鍵控制工況之一，而風(fēng)暴氣象下，高速流往往與狂風(fēng)巨浪同時出現(xiàn)。比如南海東部，熱帶氣旋作用下的海水表層流速最高可達(dá)2.5m/s以上[9]。對自升式平臺來說，當(dāng)使用莫里森公式求解波流聯(lián)合作用下的載荷時，需要將海流速度矢量與波浪水質(zhì)點速度矢量疊加，而波流載荷中的拖曳力部分與速度平方成正比，因而流的影響不可忽略。目前關(guān)于流速，特別是高流速時波流聯(lián)合對自升式平臺動力放大因子影響的研究和文獻(xiàn)資料不多。本文以筆者參與設(shè)計的某400英尺(1英尺≈30.48cm)自升式鉆井平臺作為算例平臺，研究能夠考慮流速影響的動力放大因子隨機動態(tài)響應(yīng)分析法，分析流速對動力放大因子以及慣性力的影響。

1 自升式平臺動力放大因子分析方法

自升式平臺動力放大因子的分析方法主要有單自由度法(SDOF法)和隨機波動態(tài)響應(yīng)分析法兩種[10]。單自由度法將自升式平臺簡化為質(zhì)量-彈簧-阻尼的單自由度系統(tǒng)，假定波浪為規(guī)則波，忽略波浪激振力的隨機性本質(zhì)，且假定波浪激振力作用于質(zhì)心。單自由度法計算公式如下：

(1)

式中：Ω表示周期比，Ω=Tn/Tw，Tn為平臺的自振周期，Tw為波浪周期；ξ為阻尼比，ξ≤0.07。單自由度法形式簡便，目前工程應(yīng)用較多。但單自由度法一個比較明顯的缺點是不能考慮流速的影響。

自升式平臺隨機波動態(tài)響應(yīng)分析方程如下：

(2)

(3)

隨機波動態(tài)響應(yīng)分析法比較復(fù)雜，目前工程應(yīng)用還不是很多。由于波浪拖曳力和自由面的非線性，隨機波動態(tài)響應(yīng)分析得到的自升式平臺時域響應(yīng)過程是非高斯過程，拖曳/慣性參數(shù)法(D/I法)是常用的響應(yīng)極值求解方法之一。

1.1 動力放大因子的拖曳/慣性參數(shù)法

D/I法基于這樣一個假定，即自升式平臺的動態(tài)響應(yīng)極值可由靜態(tài)響應(yīng)極值、慣性響應(yīng)極值以及兩者間的關(guān)聯(lián)系數(shù)組合而得[11]：

Rd_mpm2=(CRsσRs)2+(CRiσRi)2
+2ρR(CRsσRs)(CRiσRi)，

(4)

(5)

式中：Rd_mpm為動態(tài)響應(yīng)最大可能值；CRs、CRi分別為靜態(tài)響應(yīng)和慣性響應(yīng)最大值系數(shù)；σRd、σRs、σRi分別為動態(tài)響應(yīng)、靜態(tài)響應(yīng)和慣性響應(yīng)均方差；ρR為靜態(tài)響應(yīng)和慣性響應(yīng)最大值間的關(guān)聯(lián)系數(shù)。

(6)

(7)

D=8.0σRs(Cm=0)，

M=3.7σRs(Cd=0)，

式中：σRs(Cm=0)為不考慮慣性力的靜態(tài)響應(yīng)均方差；σRs(Cd=0)為不考慮拖曳力的靜態(tài)響應(yīng)均方差。因此，D/I法還需要分別分析不考慮慣性力和不考慮拖曳力的靜態(tài)響應(yīng)。

將式(2)中的質(zhì)量矩陣和阻尼矩陣均設(shè)為零，即為平臺的靜態(tài)響應(yīng)。平臺的慣性響應(yīng)則為動態(tài)響應(yīng)與靜態(tài)響應(yīng)之差。

平臺的動態(tài)響應(yīng)最大可能值確定后，即按下式確定平臺的動態(tài)響應(yīng)最大可能極值Rd_mpme：

Rd_mpme=μRd+Rd_mpm，

(8)

式中：μRd為動態(tài)響應(yīng)均值。

平臺的靜態(tài)響應(yīng)最大可能極值Rs_mpme按下式確定：

Rs_mpme=μRs+CRsσRs，

(9)

式中：μRs為靜態(tài)響應(yīng)均值。

最終按下式確定平臺的動力放大因子：

(10)

式中：FDAF為與周期比相關(guān)的比例放大系數(shù)，

式中： Tp表示波浪譜譜峰周期。

1.2 隨機波生成

自升式平臺隨機波動態(tài)響應(yīng)分析的隨機波可由規(guī)則波隨機線性疊加而得?？刹捎玫确诸l率法劃分頻率區(qū)間，頻率區(qū)間不少于200，每一個頻率區(qū)間對應(yīng)一個規(guī)則波。第i個規(guī)則波的波幅ai為

(11)

(12)

式中：εi為第i個組成波的初相位，是在(0, 2π)上均勻分布的隨機數(shù)。

美國船級社(ABS)推薦使用JONSWAP譜生成隨機波，峰值因子γ取3.3[10]。波浪譜的有效有義波高Hs和峰值周期Tp則由確定性波高極值Hmax和對應(yīng)周期Tass轉(zhuǎn)換而得：

(13)

式中：d表示水深，d≥25m；Hsrp表示風(fēng)暴作用下3h回歸周期對應(yīng)的有義波高，

Hsrp=Hmax/1.75.

需要指出，式(13)引入了Wheeler擴展來修正淺水自由面的非線性效應(yīng)。

2 算例分析

圖1 算例平臺模型Fig.1 Finite element (FE) model of sample unit

圖1所示為算例平臺的有限元分析模型。其中樁腿采用梁單元詳細(xì)建模，保證剛度和水動力系數(shù)與實船一致；主船體、樁靴按剛度等效原則簡化建模；樁腿與主船體的連接使用彈簧單元模擬[12]。模型采用笛卡爾直角坐標(biāo)系，坐標(biāo)原點位于船底板、中縱剖面與尾封板的交點處，縱向以指向船首為正，垂向以向上為正。限于篇幅，本文只考慮0°波流作用方向，如圖2中箭頭所示。圖2中位于船體重心的水平負(fù)剛度彈簧，用來考慮P-Δ效應(yīng)，彈簧剛度為-1 175kN/m。模型邊界條件取為泥線下3m處鉸支。計算波浪激勵力時，忽略結(jié)構(gòu)運動響應(yīng)，阻尼比取0.07。波流載荷使用莫里森公式計算后直接施加到樁腿梁單元上。

圖2 浪向角和計入P-Δ效應(yīng)的水平彈簧Fig.2 Wave heading and horizontal spring for P-Δ effect

表1列出了平臺與波流載荷隨機動力分析有關(guān)的主要參數(shù)。

表1 算例平臺主要參數(shù)Table 1 Main parameters of sample unit

注：De為等效直徑。

本文使用的流速垂向分布如圖3所示，即海底流速為零，流速從海底至瞬時波面處呈線性分布。研究選擇的表面流速Vc=0～2.0m/s，間隔0.5m/s。

圖3 流速垂向分布Fig.3 Vertical distribution of current speed

根據(jù)表1中波高和周期推算得到Hs=10.49m，Tp=15.26s，按γ=3.3的JONSWAP譜生成1h隨機波。選擇頻率ω=0.05～2.5rad/s模擬不規(guī)則波，采用等分頻率法劃分頻率區(qū)間，取Δω=0.025rad/s。算例隨機波的波面升高時歷曲線如圖4所示，對應(yīng)的波面升高平均值為0，均方差為2.62m，峰度為3.05，偏度為-0.029，最大波峰高度8.93m，滿足ABS要求[10]。

圖4 隨機波波面升高時歷Fig.4 Time history of random wave surface elevation

按D/I法的要求，通過式(2)分別進(jìn)行不同流速的動態(tài)響應(yīng)分析、靜態(tài)響應(yīng)分析、不考慮慣性力的靜態(tài)響應(yīng)分析以及不考慮拖曳力的靜態(tài)響應(yīng)分析，時間步長取0.5s。

圖5顯示了Vc=0和2.0m/s的樁基剪力的各個響應(yīng)時歷。從圖5可以看出，隨著流速的增大，平臺的動態(tài)響應(yīng)、靜態(tài)響應(yīng)和慣性響應(yīng)幅值均增大。只有不考慮拖曳力的靜態(tài)響應(yīng)不隨流速變化[圖5(g)、(h)]，這是因為流速對波流載荷的影響體現(xiàn)在拖曳力上。流速只具有單向性，隨著流速的增大，流對動態(tài)和靜態(tài)響應(yīng)的抑制作用越加顯著，平臺動態(tài)和靜態(tài)響應(yīng)均值均逐漸偏向流力作用方向一側(cè)，平臺動態(tài)響應(yīng)極值與靜態(tài)響應(yīng)極值的差異越來越小，這在圖5(b)、 (d)、 (f)中表現(xiàn)得尤為明顯。慣性響應(yīng)時歷則分別由圖5(a)與圖5(c)相減以及圖5(b)與圖5(d)相減而得。計算動力放大因子所需的響

應(yīng)均方差值和平均值已在圖中列出。

圖5 算例平臺樁基剪力響應(yīng)時歷

3 考慮流速影響的動力放大因子及慣性力比較

圖6顯示了按D/I法確定的不同流速下算例平臺樁基剪力(SF)和傾覆力矩(BM)的動力放大因子，以及按SDOF法確定的動力放大因子?？梢钥闯觯S著流速逐漸增大，D/I法確定的樁基剪力和傾覆力矩的動力放大因子均快速減小。而SDOF法確定的動力放大因子不隨流速而變，無法體現(xiàn)流速對平臺動力放大因子的影響。

圖6 不同流速下的動力放大因子比較Fig.6 Comparison of DAF at different current speeds

波流載荷對自升式平臺的動力放大效應(yīng)通過慣性力施加到平臺上。D/I法確定的慣性樁基剪力FI和傾覆力矩MI分別按以下兩式計算：

FI=(Fda_BS-1)Fmax，

(15)

MI=(Fda_BM-1)Mmax，

(16)

式中：Fda_BS和Fda_BM分別為樁基剪力和傾覆力矩的動力放大因子：Fmax和Mmax分別為通過確定性方法計算計入流速后的極端波高而得到的最大波流力和彎矩。

SDOF法確定的慣性樁基剪力FI和傾覆力矩MI分別按以下兩式計算：

FI=(Fda-1)×Fwave amp，

(17)

MI=FI×d，

(18)

式中：d為平臺垂向重心至樁腿約束點的垂向距離；Fwave amp為靜態(tài)波流力幅值，

Fwave amp=0.5(Fmax-Fmin)，

Fmin為通過確定性方法計算計入流速后的極端波高而得到的最小波流力。

按表1中列出的極端波高和對應(yīng)周期，通過Stokes五階波進(jìn)行算例平臺的波流載荷確定性分析。不同流速下的波流載荷最大值和最小值在表2中列出。

表2 不同流速下的波流力、彎矩最大值和最小值Table 2 Maximum and minimum wave forces and overturning moments at each current speed

(a) 樁基剪力

(b) 傾覆力矩圖7 不同流速下的慣性力比較Fig.7 Comparison of inertial force at different current speeds

圖7分別顯示了按D/I法和SDOF法計算得到的慣性樁基剪力和傾覆力矩的差異。從圖7可以看到，按SDOF法計算得到的慣性力顯著大于按D/I法計算得到的慣性力。隨著流速的增大，兩種方法計算得到的慣性力都增大，但SDOF法計算得到的慣性力增大更為劇烈。流速2.0m/s時兩種方法得到的慣性力差值基本達(dá)到不考慮流速時的2倍，這說明隨著流速的增大，SDOF法得到的慣性力越來越偏于保守。進(jìn)一步比較表2中的靜態(tài)波浪力，可以看到隨著流速的增大，F(xiàn)max增大的速度要比Fwave amp快得多。之所以D/I法計算得到的慣性力比SDOF法要小得多，原因是D/I法確定的動力放大因子隨著流速增大而減小，抵消了Fmax快速增大的影響。

4 結(jié) 語

通過隨機波動態(tài)響應(yīng)分析法分析了某400英尺自升式鉆井平臺在不同流速下的動力放大因子和慣性力變化情況，并與單自由度法計算結(jié)果進(jìn)行了簡單比較。研究表明：

(1) 由于流是方向不隨時間而變的定常流，即流載荷具有單向性，故流速越大，它對自升式平臺運動響應(yīng)的抑制越顯著。盡管隨著流速的增大，平臺的動態(tài)、靜態(tài)響應(yīng)極值均增大，但兩者間的差異越來越小，造成平臺的動力放大因子越來越小。這從物理上解釋了自升式平臺動力放大因子隨流速增大而減小的原因。

(2) 在波浪參數(shù)不變的前提下，流速增大，隨機波動態(tài)響應(yīng)分析法和單自由度法得到的自升式平臺慣性力都隨之增大。隨機波動態(tài)響應(yīng)分析法考慮了流速的影響，該方法得到的慣性力增大相對平緩。而單自由度法得到的慣性力增大較為劇烈，越來越偏于保守。兩種方法確定的平臺動力放大因子間的較大差異，是造成兩種方法得到的慣性力間較大差異的主要原因。

因此，為合理評估自升式平臺的動力放大因子和慣性力，盡管過程更為繁瑣，流速較大時還是應(yīng)采用隨機波動態(tài)響應(yīng)分析法。目前，F(xiàn)&G和GustoMSC等世界主流設(shè)計公司都已引入隨機波動態(tài)響應(yīng)分析法進(jìn)行高規(guī)格、大水深自升式平臺的設(shè)計。應(yīng)用隨機波動態(tài)響應(yīng)分析法進(jìn)行自升式平臺的動力放大因子和慣性力分析，能夠優(yōu)化風(fēng)暴環(huán)境條件下的樁腿設(shè)計載荷，避免無謂的樁腿結(jié)構(gòu)加強，有利于控制空船重量，增加平臺的可變載荷設(shè)計指標(biāo)，從而提高平臺在國際市場上的競爭力。這也是本研究的實際工程意義所在。

[1] Jenson J J. Dynamic amplification of jack-up platforms subjected to non-gaussian wave loads [J]. Journal of Research of the National Institute of Standards and Technology, 1998,99(4)： 455.

[2] Lu Y J, Chen Y N, Tan P L, et al. Prediction of most probable extreme values for jack-up dynamic analysis [J]. Marine Structures, 2002,15(1): 15.

[3] Williams M S, Thompson R S G, Houlby G T. A parametric study of the non-linear dynamic behaviour of an offshore jack-up unit [J]. Engineering Structures, 1999,21(5)： 383.

[4] Cassidy M J, Taylor P H, Taylor R E, et al. Evaluation of long-term extreme response statistics of jack-up platforms [J]. Ocean Engineering, 2002,29(13)： 1603.

[5] Zhang X Y, Cheng Z P, Wu J F. Dynamic response of jack-up units—reevaluation of SNAME 5-5A four methods [C]. ISOPE, 2010： I-10-587.

[6] 黃盟宇,唐文勇,汪怡,等.自升式海洋平臺動態(tài)放大因子分析與抑制[J].上海交通大學(xué)學(xué)報,2010,44(6)： 816.

Huang Meng-yu, Tang Wen-yong, Wang Yi, et al. Analysis and suppression on the dynamic amplification factor of self-elevating offshore platform [J]. Journal of Shanghai Jiao Tong University, 2010,44(6)： 816.

[7] Saeidtehrani S, Asgarian B, Daghigh M, et al. Calculation of dynamic amplification factor for jack-up structures [J]. International Geoinformatics Research and Development Journal, 2015,6(3)： 15.

[8] Perry M J, Mobbs B. Factoring to improve conservatism of the drag-inertia method when determining dynamic amplification of jack-up units [C]. OMAE, 2011： 50335.

[9] 陳順楠,喬方利,潘增弟,等.中國南海東部海域氣候特征及風(fēng)浪流極值參數(shù)的研究—LAGFD數(shù)值模式群的運用[J].黃渤海海洋,1998,16(2)： 6.

Chen Shun-nan, Qiao Fang-li, Pan Zeng-di, et al. The study of climatic characters and wind, wave, current extreme parameters in the east of South China Sea： the application of LAGFD numerical models [J]. Journal of Oceanography of Huanghai & Bohai Seas, 1998,16(2)： 6.

[10] American Bureau of Shipping. Guidance notes on dynamic analysis procedure for self-elevating units [S]. 2014.

[11] Maatschappij B V. Dynamic analysis and estimation of extreme response for jack-ups [R]. 1991.

[12] 郭東杰.在役自升式平臺風(fēng)暴自存狀態(tài)強度評估[D].天津： 天津大學(xué),2009.

Guo Dong-jie. Strength evaluation of jack-up in use under storm survival condition [D]. Tianjin： Tianjin University, 2009.

InfluenceofCurrentSpeedonDynamicAmplificationFactorofJack-UpUnits

MA Wang-kou1, YUAN Fei-hui2

(1.MarineDesign&ResearchInstituteofChina,Shanghai200011,China;2.ShanghaiWaigaoqiaoShipbuildingCo.,Ltd.,Shanghai200137,China)

Random wave dynamic response analysis method is studied and used to analyze a sample drilling jack-up unit. The result shows that current suppression on the dynamic and static responses of the unit increases as the current speed increases, but the dynamic amplification factor decreases meanwhile. The result also shows that the inertial force computed by random wave dynamic response analysis method increases relatively slowly since this method can take account of current effect. So random wave dynamic response analysis method should be selected to calculate dynamic amplification factor and inertial force under large current speed for optimizing structure design of legs and increasing variable deck load of jack-up units.

dynamic amplification factor; random wave dynamic response analysis method; current speed; inertial force

2017-06-19

工業(yè)和信息化部高技術(shù)船舶科研項目“自升式鉆井平臺品牌工程(Ⅰ型)”(工信部聯(lián)裝[2014]507號文)

馬網(wǎng)扣(1978—)，男，高級工程師，主要從事海洋工程設(shè)計和研發(fā)工作。

U661.3

A

2095-7297(2017)05-0293-07