海洋樁基水動力系數(shù)的確定方法研究

(山東科技大學(xué) 山東 青島 266590)

1 引言

樁基結(jié)構(gòu)是海洋風(fēng)電基礎(chǔ)中適用于近海海域的結(jié)構(gòu)形式[1]。在樁基穩(wěn)定性計算中，基樁的波浪力計算至關(guān)重要。目前計算小直徑樁的波浪力一般采用莫里森方程，其中的關(guān)鍵是合理確定慣性力系數(shù)CM和拖曳力系數(shù)CD[2]。目前確定CM和CD的主要有三種，一是利用瞬時值計算數(shù)CM和CD，二是利用傅里葉分析法計算CM和CD，三是利用最小二乘法計算CM和CD。本文就針對這三種計算方法并結(jié)合《海港水文規(guī)范》[3]以及美國API中所規(guī)定的CM和CD，分析這幾種方法的優(yōu)劣。

2 慣性力系數(shù)CM和拖曳力系數(shù)CD的確定

2.1 瞬時值計算CM和CD。根據(jù)線性波理論，當(dāng)cosωt=0時，速度u=0，加速度a有極值，利用該瞬時波浪力資料可計算CM值；當(dāng)sinωt=0時，加速度a=0，速度u有極值，利用該瞬時波浪力資料可計算CD的值。

這種方法雖然計算簡單，但必須保證所取瞬時資料的速度及加速度兩個量中，一個值為極值，同時另一個為0.由于波浪的非線性效應(yīng)，這一條件不一定能夠成立。此外，還應(yīng)考慮觀測資料的某些不規(guī)律性。為了減小這種不規(guī)律性對結(jié)果所造成的誤差，有時對此瞬時附近的一個小時段取多個值，分析后進(jìn)行平均。

2.2 傅里葉分析法計算CM和CD。在穩(wěn)定流中，試驗證明CM和CD為雷諾數(shù)Re的函數(shù)，而在波浪運(yùn)動中雷諾數(shù)是呈周期性變化的。美國的Keulegan和Carpenter發(fā)現(xiàn)系數(shù)CM和CD與拖曳力和慣性力的比值有關(guān)，并定義KC數(shù)為

式中Umax為波浪水質(zhì)點(diǎn)運(yùn)動的最大速度；T為波的周期。

目前人們大多認(rèn)為在波動水流中系數(shù)CM和CD是與Re以及KC兩個參數(shù)有關(guān)。

Keulegan和Carpenter提出由于在波浪水流中流速及加速度兩值均隨時間而周期性地變動，則CM和CD也應(yīng)在一個周期中統(tǒng)一地加以考慮，即可以在一個周期中用傅里葉分析法來取平均值，該值稱為傅里葉分析均值。他們認(rèn)為波浪力應(yīng)為一奇函數(shù)，即半周期內(nèi)異號，則：

+ B1cosθ + B3cos3θ + B5cos5θ + …)

改寫成與莫里森方程相對應(yīng)的表達(dá)式可得到：

-CD|cosθ|cosθ + 2(B3cos3θ + B5cos5θ + …)

如果僅取和莫里森方程相同的項可得：

應(yīng)特別指出的是，在此認(rèn)為CM和CD與相位角θ無關(guān)，在雷諾數(shù)Re和KC數(shù)已知時為一常數(shù)，且當(dāng)n>3時，An和Bn各項均為0。

將(1)式的等號兩側(cè)各乘以cosθ并對θ在(0,2π)的范圍內(nèi)積分可得

將(1)式的等號兩側(cè)各乘以sinθ并對θ在(0,2π)的范圍內(nèi)積分可得

由上述兩式就可求得按傅里葉分析法所得的CM和CD的均值。另外也可根據(jù)最小二乘法原理，采用該種方法對CM和CD進(jìn)行計算，此時所選用的CM和CD值應(yīng)該符合以下條件：E2對CM和CD微分為0。

2.3 各國規(guī)范對CM和CD的建議值。美國API規(guī)定為：CD為0.6～1.0，CM為1.5～2.1；《海港水文規(guī)范》規(guī)定為：CD為1.2，CM為2.0；挪威DNV規(guī)定為： CD為0.5～1.0，CM為2.0；美國DTI規(guī)定為采用可靠的試驗結(jié)果。

從上述數(shù)據(jù)可以看出，各國對CM和CD的定義不盡相同，美國API針對桿件粗糙和光滑，定義了不同的CM和CD值，國內(nèi)并沒有針對圓柱光滑度的規(guī)定，取值相對來說比較保守。最小二乘法(又稱最小平方法)是一種數(shù)學(xué)優(yōu)化技術(shù)。它通過最小化誤差的平方和尋找數(shù)據(jù)的最佳函數(shù)匹配。利用最小二乘法可以簡便地求得未知的數(shù)據(jù)，并使得這些求得的數(shù)據(jù)與實際數(shù)據(jù)之間誤差的平方和為最小，最小二乘法也可用于參數(shù)的計算。

3 波浪要素的計算

3.1 水質(zhì)點(diǎn)的速度和加速度的計算方法.

(1)海浪譜。隨機(jī)波浪與頻率的選取有著較為密切的關(guān)系，頻率和波高、波浪頻率等波浪特征值有密切的聯(lián)系。實用海上波譜材料為根據(jù)大量實測數(shù)據(jù)統(tǒng)計繪制的。擬合相應(yīng)的隨機(jī)波浪荷載是為了反映相類似的真實環(huán)境中海浪的情況。為了反映工程中海上風(fēng)機(jī)所受的實際波浪力而在此采用我國的海浪譜——會戰(zhàn)譜[6]。其形式如下：

上式中：ω為圓頻率；HS為隨機(jī)波浪有效波高。

(2)波面方程。將波浪視為平穩(wěn)隨機(jī)過程，假定其是由多個頻率不同、振幅不同、波數(shù)不同、初始隨機(jī)相位不同的線性波組成的，文中認(rèn)為隨機(jī)波浪由若干余弦波擬合形成。任意組成波的形式如下式：

ηi(t)=Aicos(kix-ωit+εi)

式中：η(t)為波面方程；Ai為任一組成波的振幅；ki為任一組成波的波數(shù)；ωi為任一組成波的圓頻率；εi為任一組成波的v初始隨機(jī)相位，其值平均分布在[0,2π]范圍內(nèi)。

根據(jù)Airy波理論，采用諧波疊加法，則隨機(jī)波浪的波面時程方程如下所示：

任意深度d處的速度時程方程和加速度時程方程如下所示：

3.2 波浪力的計算方法。在海洋工程中，計算小直徑樁波浪力時，經(jīng)常采用Morison方程進(jìn)行計算。Morison方程是1950年由美國加州大學(xué)伯克利分校的Morison等提出的，是一個半理論半經(jīng)驗公式。該方法認(rèn)為結(jié)構(gòu)上的波浪力可以分為兩部分，一部分是波浪質(zhì)點(diǎn)速度所引起的速度力，另一部分是由波浪質(zhì)點(diǎn)的加速度所引起的慣性力。作用于單位長度樁體上的波浪力表達(dá)式如下：

式中各符號含義為：CD為拖曳力系數(shù)；CM為慣性力系數(shù)；ρ為海水密度；D為單樁直徑。

從上式中我們不難發(fā)現(xiàn)拖曳力是非線性的，所以需要通過Borgman L.E方程進(jìn)行線性化處理。σu是線性化的關(guān)鍵，其計算式中積分上下限的取值應(yīng)參照實際工程。

在實際計算中，可以采用等分頻率法將截斷頻率分段形成若干個不同大小的線性波浪，采用上述參數(shù)，通過現(xiàn)行疊加法得到相應(yīng)的時程方程。

3.3 工程概況

(1)基本概況。上海東海大橋100MW海上風(fēng)電示范項目工程位于上海市東海大橋東部海域，總裝機(jī)容重102MW，安裝34臺華銳風(fēng)電科技有限公司生產(chǎn)的單擊容量3MW的SL3000離岸型風(fēng)電機(jī)組。樁柱直徑D=1.7m，水深20.6m。

(2)水文特征。潮位特征(國家85高程)為：平均海平面為0.23米，平均高潮面為1.86米，平均低潮面為-1.34米，設(shè)計高潮面為2.55米，設(shè)計低潮面為-2.09米，極端高潮面為3.68米，極端低潮面為-2.93米。

設(shè)計波浪要素為：重現(xiàn)期為50年，平均波高為2.83米，平均周期長為7.76秒，波長為74.1米，H1%為5.81米，，H4%為5.06米，H5%為4.92米，H13%為4.24米。

4 計算結(jié)果分析

從會戰(zhàn)譜曲線圖中我們可以看出，譜值在圓頻率達(dá)到2.5以后的變化量幾乎很小，所以取0.2<ω<2.5范圍內(nèi)進(jìn)行計算。與此同時，我們還可以采用3.1中的計算方法，求得波面方程時程曲線。

由瞬時值法、傅里葉分析法、最小二乘法和《海港水文規(guī)范》、美國API確定的CM和CD如表4所示。

表4 CM和CD取值

CDCM規(guī)范法1.22最小二乘法0.581.42傅里葉分析法0.621.56瞬時值法0.811.12

據(jù)上述幾種方法確定的CM和CD，代入Morison方程計算得到的單位高度上的波浪力時程曲線，如圖5所示。

從圖中我們可以看出，瞬時值法和《海港水文規(guī)范》確定的CM和CD值偏差較大，而通過傅里葉分析法和最小二乘法計算得到的CM和CD值比較接近且更為合理。

5 結(jié)論

由瞬時值法計算得到的CM和CD值，由于其計算取值的條件過于苛刻，而實際情況并不十分符合計算原理，所以得到的計算結(jié)果誤差較大。而《海港水文規(guī)范》中所規(guī)定的CM和CD值的取值相對比較保守，根據(jù)國內(nèi)規(guī)范計算得出的拖曳力和慣性力較大。傅里葉分析法和最小二乘法計算得到的CM和CD值比較合理，在實際工程中可以根據(jù)已有的波浪數(shù)據(jù)選擇適合的方法進(jìn)行計算。