淺水低頻長波對海洋核動力平臺動力響應(yīng)的影響分析

2019-09-02 06:52:38梁雙令

兵器裝備工程學(xué)報 2019年8期

梁雙令，程堃

(武漢第二船舶設(shè)計研究所，武漢 430064)

在工程設(shè)計中，波浪譜一般通過理論公式得到，這些理論波浪譜的能量主要集中在0.2～2.0 rad/s頻率范圍內(nèi)，在0～0.2 rad/s的低頻范圍內(nèi)，一般沒有波浪能量成分。但是在實際淺水海域，伴隨著大幅度振蕩的波浪時歷，總會出現(xiàn)較明顯的低頻波浪成分，這即是在淺水中存在的低頻長波現(xiàn)象[1]。淺水低頻長波及其誘導(dǎo)的一階波浪力會導(dǎo)致總的波浪力大幅度增加，從而導(dǎo)致單點系泊船體的運動響應(yīng)以及由此產(chǎn)生的系泊力相比于數(shù)值計算有一定幅度的增大。因此，為保證船體運動響應(yīng)和系泊力預(yù)報的準(zhǔn)確性，在總體設(shè)計流程中必須要考慮淺水低頻長波的影響。

關(guān)于淺水低頻長波的存在及其形成機理，已被許多理論研究和實際觀測所證實[2]。為分析淺水低頻長波對船體縱蕩、垂蕩和縱搖動力響應(yīng)的影響，主要包括波浪力和運動響應(yīng)，許多學(xué)者都采用數(shù)值計算與水池試驗相結(jié)合的方法。肖龍飛等[3]通過頻域數(shù)值計算和模型試驗研究，分析了淺水不規(guī)則波浪譜中低頻成分對低頻縱蕩運動的影響，結(jié)果表明淺水低頻縱蕩運動響應(yīng)主要是一階的，二階成分幾乎可以忽略。李欣等[4]在線性三維勢流理論的基礎(chǔ)上，采用時域計算方法分析了淺水軟剛臂系泊FPSO的觸底情況，并與模型試驗進行了對比。郭彬[5]通過時域數(shù)值計算和模型試驗分析了船體在淺水非共線風(fēng)浪流海況作用下的運動性能與舷側(cè)甲板上浪情況。

本研究針對作業(yè)于渤海淺水海域的海洋核動力平臺示范工程，通過數(shù)值計算和水池試驗兩種方法得到淺水波浪譜密度，通過結(jié)果對比，分析淺水低頻長波對平臺縱蕩、垂蕩和縱搖三個自由度上的波浪力和運動響應(yīng)的影響。

1 理論基礎(chǔ)

1.1 縱蕩運動方程

單點系泊平臺在迎浪不規(guī)則波作用下的縱蕩運動方程可表示為[6]：

(1)

式中：x1為縱蕩運動；a11(μ1)和B11(μ1)分別為固有頻率μ1對應(yīng)的附加質(zhì)量和靜水阻尼系數(shù)；C11為單點系泊系統(tǒng)提供的等效線性縱蕩水平回復(fù)力系數(shù)；F1為縱蕩波浪力，包括一階成分和二階成分。

對于一階波浪力，其譜密度函數(shù)可直接根據(jù)波浪譜密度函數(shù)S(ω)和波浪力幅值響應(yīng)函數(shù)HF1(ω)計算，如式(2)所示，其標(biāo)準(zhǔn)差可由式(3)計算得到。

(2)

(3)

二階平均波浪力決定了平臺的平均位移，可以表示為：

(4)

式中，T(ω,ω)為平均波浪力二次傳遞函數(shù)。

二階波浪慢漂力譜密度函數(shù)如式(5)所示，而其標(biāo)準(zhǔn)差可由式(6)計算得到。

(5)

(6)

式中，T(ω+μ,ω)為波浪慢漂力二次傳遞函數(shù)。一階波浪力和二階波浪慢漂力決定了平臺在平均位移附近的運動幅度。

對于一階縱蕩運動，其譜密度函數(shù)也可由波浪譜密度函數(shù)S(ω)和運動幅值響應(yīng)函數(shù)Hx1(ω)計算，如式(7)所示，其標(biāo)準(zhǔn)差可由式(8)計算得到。

(7)

(8)

以二階低頻縱蕩波浪力作為輸入，以二階低頻縱蕩運動作為輸出，二者的能量譜密度函數(shù)關(guān)系可以通過線性幅值傳遞函數(shù)進行表示，如式(9)所示。

(9)

(10)

由式(10)可知，二階低頻縱蕩運動不僅受二階縱蕩波浪力影響，也與運動阻尼、回復(fù)力系數(shù)有關(guān)。其中回復(fù)力系數(shù)與單點系泊系統(tǒng)的設(shè)計有直接關(guān)系，在本文討論平臺動力響應(yīng)時，假設(shè)其為常量。因此影響縱蕩運動響應(yīng)大小的因素主要有縱蕩波浪力和運動阻尼，而后者又包含靜水阻尼和波浪慢漂阻尼兩部分。

靜水阻尼系數(shù)與黏性有關(guān)，一般根據(jù)靜水衰減試驗結(jié)果計算得到，如式(11)所示[8]。

(11)

(12)

式中：n為縱蕩衰減次數(shù)；xn為第i次縱蕩運動幅值。

縱蕩波浪慢漂阻尼是由波浪引起的，可以通過平均波浪力的二次傳遞函數(shù)計算得到，如式(13)和(14)所示[9]。

(13)

(14)

1.2 垂蕩和縱搖運動方程

對于垂蕩和縱搖運動，由于其固有頻率處在波頻范圍內(nèi)，所以主要表現(xiàn)為一階運動，而且其自身靜水回復(fù)力遠大于系泊系統(tǒng)提供的回復(fù)力，因此一般在頻域計算中不考慮系泊系統(tǒng)的影響。基于線性勢流理論，自由漂浮的平臺在頻率ω的規(guī)則波作用下的一階運動方程如式(15)所示。

(15)

與縱蕩一階波浪力和運動響應(yīng)相同，垂蕩和縱搖一階波浪力的譜密度函數(shù)和標(biāo)準(zhǔn)差可通過式(2)和式(3)計算得到，垂蕩和縱搖運動響應(yīng)的譜密度函數(shù)和標(biāo)準(zhǔn)差可通過式(7)和式(8)計算得到。

2 海洋核動力平臺簡介

海洋核動力平臺是我國核動力水面艦船示范工程，屬浮動式核動力能源保障船式平臺，由軟剛臂單點系泊裝置長期定位于渤海海域，主要解決渤海油氣開采中的電力、淡水等能源需求。

平臺有兩部分組成：船體和軟剛臂單點系泊系統(tǒng)，分別如圖1和圖2所示。平臺作業(yè)場址平均海平面30.6 m，設(shè)計重現(xiàn)期為500年。基于國家海洋環(huán)境預(yù)報中心監(jiān)測數(shù)據(jù)，推算得到重現(xiàn)期500年最低水位相對于平均海平面為-5.1 m，因此平臺的最小作業(yè)水深為25.5 m，屬于淺水海域。

圖1 船體網(wǎng)格劃分示意圖

圖2 軟剛臂系泊系統(tǒng)總體布置

單點系泊系統(tǒng)的縱蕩水平回復(fù)力曲線如圖3所示，圖3中x為縱蕩位移，F(xiàn)x為縱蕩水平回復(fù)力，二者成非線性關(guān)系，尤其是當(dāng)位移比較大時更為明顯。

基層思想政治工作是促和諧、保穩(wěn)定的生命線，是企業(yè)生產(chǎn)經(jīng)營管理實現(xiàn)良性發(fā)展的有力思想保證。基層思想政治工作發(fā)揮著統(tǒng)一思想、提高認(rèn)識，化解矛盾、理順情緒，凝聚力量、振奮精神的重要作用，是協(xié)調(diào)關(guān)系、理順情緒的潤滑油，是宣傳政策、解惑釋疑的穩(wěn)壓器，是統(tǒng)一思想、鼓舞士氣的凝固劑。它雖然不直接產(chǎn)生經(jīng)濟效益，但它猶如水、空氣對于人一樣，是確保基層和諧、穩(wěn)定發(fā)展的生命線。

圖3 單點系泊系統(tǒng)縱蕩水平回復(fù)力曲線

3 數(shù)值計算結(jié)果與比較分析

3.1 不規(guī)則波校驗

本研究選取在上海交通大學(xué)海洋工程水池中模擬的一組不規(guī)則波，Hs=5.0 m，Tp=10.3 s，γ=3，JONSWAP譜。波浪譜密度的計算值與試驗值如圖4所示。

圖4 波浪譜密度計算值與試驗值對比結(jié)果

由圖4可知，波浪譜密度計算值與試驗值吻合良好，但試驗值在0～0.2 rad/s的低頻范圍內(nèi)出現(xiàn)了計算值所沒有的能量成分，這即是在淺水中出現(xiàn)的低頻長波現(xiàn)象。該淺水低頻長波的存在，將對平臺的運動響應(yīng)和單點系泊系統(tǒng)的設(shè)計產(chǎn)生重要的影響，在實際工程設(shè)計中需要給予考慮。

3.2 縱蕩動力響應(yīng)

根據(jù)式(2)、式(4)、式(5)可分別計算得到縱蕩一階波浪力譜、二階平均波浪力和二階波浪慢漂力譜，分別如圖5～圖7所示。為分析淺水低頻長波對縱蕩波浪力的影響，以標(biāo)準(zhǔn)差為判斷標(biāo)準(zhǔn)，在波浪譜計算值和試驗值兩種情況下的縱蕩一階波浪力和波浪漂力標(biāo)準(zhǔn)差無因次比值以及平均波浪力大小的無因次比值如圖8所示。

圖5 縱蕩一階波浪力譜

圖6 縱蕩二階平均波浪力

圖7 縱蕩二階波浪慢漂力

圖8 縱蕩波浪力標(biāo)準(zhǔn)差對比

對比波浪譜計算值和測量值兩種計算結(jié)果：

1) 對比圖5和圖7可知，在0～0.2 rad/s的低頻范圍內(nèi)，一階波浪力為1012量級，二階波浪慢漂力為1011量級，兩者相差十倍至幾十倍。這說明，在淺水海域，由于低頻長波的存在，低頻波浪力的主要貢獻為一階低頻波浪力，這一點與深水海域截然不同。

2) 由圖8可知，淺水低頻長波對縱蕩一階波浪力、二階平均波浪力影響較小，而對二階波浪慢漂力影響較大，考慮淺水低頻長波相比于不考慮淺水低頻長波時二階波浪慢漂力增大一倍左右。

由式(4)計算得到平均波浪力大小為186.8 kN，重現(xiàn)期500年對應(yīng)迎風(fēng)迎流條件下的風(fēng)載荷和流載荷分別為340.4 kN和33.8 kN，因此對應(yīng)平均系泊力為561.0 kN。由圖3可知，在該平均系泊力下，平臺的平均位移為-1.70 m，此時系泊系統(tǒng)的縱蕩水平回復(fù)力系數(shù)C11(位移處取斜率)為3.48×105N/m，則縱蕩固有頻率μ1可由式(16)計算得到。

(16)

式中，a11(∞)為頻率為∞時的縱蕩附加質(zhì)量，此時μ1為0.104 rad/s。

根據(jù)式(11)和式(13)分別計算得到B11和Bwdd為1 001.1 kN·s/m和39.8 kN·s/m。因此可由式(8)和式(10)計算得到平臺的一階縱蕩運動、二階低頻縱蕩運動和總的縱蕩運動的標(biāo)準(zhǔn)差，如圖9所示。