小直徑單立柱上的波浪荷載研究

韓非非，張 瑋，陳 愷

（河海大學(xué)港口海岸與近海工程學(xué)院，南京210098）

由于海洋環(huán)境的特殊性，海洋建筑物會受到波浪、海流、海冰等動力作用，承受著與陸地環(huán)境不同的外部荷載。波浪荷載計算是海洋建筑物設(shè)計中不可或缺的一個環(huán)節(jié)，已發(fā)展成為多種計算理論［1］。目前，研究海洋建筑物上的波浪荷載時，通常采用2種方法［2］：一種是“設(shè)計波”法，用一個給定波浪周期和波高的波浪，代表一定環(huán)境下的最大波，此法只能研究規(guī)則波；另一種是“設(shè)計譜”法，研究不規(guī)則波浪的長期作用規(guī)律，其中波譜是給定波高和周期相應(yīng)的波譜。

無論采用“設(shè)計波”法還是“設(shè)計譜”法，數(shù)值模擬已成為一種趨勢，在波浪荷載研究中被廣泛應(yīng)用。然而在數(shù)值模擬的基礎(chǔ)上，最好能有一定的實測數(shù)據(jù)或可作為驗證資料的數(shù)據(jù)，以此來驗證數(shù)值模擬的準(zhǔn)確性，這是目前波浪荷載研究中的一個空缺。然而用實測數(shù)據(jù)和實驗數(shù)據(jù)作為驗證資料，需要投入大量的人力物力，且人為因素會造成許多誤差，有很大的局限性。用理論計算值作為驗證資料則方便、準(zhǔn)確。

本文以微幅波波浪理論為例，采用“設(shè)計波”法計算小直徑立柱上的波浪荷載，利用計算小直徑立柱上波浪荷載的常用公式——Morison公式，得到波浪荷載的理論計算值［3］，以此作為驗證資料，建立有限元模型，將模型計算值與理論計算值進行對比，進而研究計算參數(shù)計算時長（time）和時間步長（substeps）的最佳值。

1 單立柱上波浪荷載的理論計算方法

1.1 Morison方程

目前，計算作用于小尺度立柱上的波浪荷載時，最常用的公式是Morison方程［2］，雖然只是計算波浪荷載的半經(jīng)驗公式，但已在工程中得到廣泛應(yīng)用。

根據(jù)Morison方程，作用在圓柱上的作用力可以分為兩部分：一部分為由于水質(zhì)點運動而產(chǎn)生沖擊樁的拖曳力，由流體質(zhì)點水平速度引起；另一部分為波浪水質(zhì)點運動引起的對樁柱的慣性力，由流體質(zhì)點水平加速度項引起［3］。公式的物理意義在于將波浪作用于樁柱結(jié)構(gòu)上的力分解為速度分力和慣性分力，再按力矢合成原理，將速度分力和慣性分力疊加，其合力是波浪對樁柱結(jié)構(gòu)的作用力［4］。

小尺度構(gòu)件垂直于其軸線方向單位長度上的波浪力F為

式中：FD為拖曳力；FI為慣性力。其分別可表示為

式中：D為結(jié)構(gòu)物的特征尺度；A為結(jié)構(gòu)物橫截面面積；ρ為流體的密度；ux為垂直于構(gòu)件軸線水質(zhì)點速度分量。當(dāng)海流與波浪聯(lián)合作用時，為波浪水質(zhì)點的速度矢量與海流速度矢量之和在垂直于構(gòu)件方向上的分量；?ux/?t為垂直于構(gòu)件軸線水質(zhì)點加速度分量。當(dāng)海流與波浪聯(lián)合作用時，為波浪水質(zhì)點的速度矢量與海流加速度矢量之和在垂直于構(gòu)件方向上的分量；CD為阻力系數(shù)；CM為慣性系數(shù)（或質(zhì)量系數(shù)）；對圓形構(gòu)件，可取CD＝0.6～1.2，CM＝1.3～2.0。

任意相位時作用于整個圓柱體上的總水平波浪力為

1.2 微幅波表示的波浪理論

式（2）和式（3）中的ux根據(jù)選取的波浪理論具有不同的表達形式，常用的波浪理論有線性波理論和非線性波浪理論。線性波理論主要是Airy波浪理論，非線性波理論主要包括Stocks波浪理論、橢圓余弦波理論和流函數(shù)理論［5］。

本文以微幅波理論為例，得到微幅波理論表示的單個直立柱體上的水平總波浪力和水平總波浪力力矩，如下式所示

1.3 計算實例

采用微幅波波浪理論表示的Morison公式，計算一個直立樁柱上的波浪荷載，基本參數(shù)見表1。計算時代入不同相位角對應(yīng)的K1、K2、K3和K4值，得到理論計算結(jié)果（表2）。

表1 計算基本數(shù)據(jù)Tab.1 Basic data for calculation

2 用ANSYS計算單立柱上的波浪荷載

2.1 單立柱上的總波浪力

本文按照上述參數(shù)建立有限元模型（圖1）。采用微幅波波浪理論，入射角為0°，輸入水流參數(shù)后，計算得到單立柱上的總波浪力與總波浪力力矩。

圖1 有限元模型Fig.1 Finite element model

表2 單立柱上波浪荷載的理論計算結(jié)果Tab.2 Theoretical calculation results of wave load on single pile

一般數(shù)值模擬的誤差來源于3個方面：模型誤差、舍入誤差和截斷誤差。

數(shù)值計算時通常會對實際問題進行簡化。由簡化后的數(shù)學(xué)模型與實際問題的差異產(chǎn)生的誤差即為模型誤差。計算時用精度有限的有理數(shù)代替精度無限的無理數(shù)而產(chǎn)生的誤差即為舍入誤差，此誤差由計算精度引起。截斷誤差是由不含舍入誤差的積分步長產(chǎn)生，減小積分步長可以減小截斷誤差，但是會增加計算時間，且高階算法將會增加數(shù)值的不穩(wěn)定性。因此在數(shù)值模擬時，在不計模型誤差的基礎(chǔ)上，時間步長的選取存在一個恰當(dāng)值［6］。

表3 不同計算時長的相對誤差范圍Tab.3 Relative error range of different calculation times %

本文通過與理論計算值的比較，發(fā)現(xiàn)對小直徑立柱的波浪作用進行數(shù)值模擬時，計算時長和計算步長對計算結(jié)果的影響較大，并存在一個可使計算結(jié)果趨于穩(wěn)定的恰當(dāng)值。計算結(jié)果與理論值的誤差范圍見表3和表4。從表3和表4中可以看出，當(dāng)計算時長大于15倍波浪周期，時間步長約為波浪周期的1/20時，計算結(jié)果與理論值的相對誤差較小。

由計算結(jié)果（表5）可以看出，用本文的計算模式計算單立柱上的波浪荷載時，其結(jié)果十分接近理論計算值，相對誤差基本維持在3%以內(nèi)。

表4 不同時間步長的相對誤差范圍Tab.4 Relative error range of different calculation time steps %

2.2 不同水深時單立柱上的波浪荷載分布

表5 模型計算值與理論計算值的比較Tab.5 Comparison of model calculation and theoretical calculation value

計算深水區(qū)單立柱上的波浪荷載時，發(fā)現(xiàn)立柱剛度較小時，計算結(jié)果會不穩(wěn)定，嚴重偏離理論值（圖2）。若增大立柱剛度，可使計算結(jié)果與理論值符合良好。

圖3和圖4分別給出100 m水深和40 m水深時對應(yīng)的波浪力分布，觀察總波浪力最大值對應(yīng)的單立柱上波浪荷載的分布，可以看出立柱上波浪荷載隨水深呈現(xiàn)遞減的趨勢，且集中在靜水面附近。

不同的是，深水時波浪力沿水深急劇減小，減小到一定值后基本保持不變，并在水底附近達到最小值；淺水時波浪荷載沿水深逐漸減小，變化較為均勻，在水底處達到最小值，且此時底部管樁所受的波浪力遠大于深水底部管樁所受的波浪力。由波浪理論可知，此變化規(guī)律符合質(zhì)點水平速度隨水深的分布規(guī)律［7-8］。

以上分析表明，在采用ANSYS模擬單立柱上的波浪作用時，其計算結(jié)果合理可信，可將此波浪力計算方法推廣到工程應(yīng)用中。

圖2 小剛度立柱上的波浪荷載Fig.2 Wave load on small stiffness pile

圖3 100 m水深單立柱上波浪荷載沿水深的分布Fig.3 Distribution of wave load on single pile on 100 m depth

圖4 40 m水深單立柱上波浪荷載沿水深的分布Fig.4 Distribution of wave load on single pile on 40 m depth

3 結(jié)論

（1）目前計算作用在小直徑立柱上的波浪力時，常用Morison公式，本文用微幅波波浪理論表示的Morison公式計算得到隨相位角變化的單立柱上的總波浪力。

（2）本文用ANSYS模擬小直徑單立柱上的波浪作用，在不考慮模型誤差的基礎(chǔ)上，發(fā)現(xiàn)計算總時長和時間步長對計算結(jié)果影響較大，并存在一個可使計算結(jié)果趨于穩(wěn)定的恰當(dāng)值，即當(dāng)計算時長大于15倍波浪周期，時間步長約為波浪周期的1/20時，計算結(jié)果較為準(zhǔn)確。

（3）在計算深水區(qū)單立柱上的波浪荷載時，發(fā)現(xiàn)立柱剛度較小時，會使計算結(jié)果不穩(wěn)定，嚴重偏離理論值。若增大立柱剛度進行計算，可使其與理論值符合良好。

（4）由本文推薦計算模式得到的計算結(jié)果與理論計算值十分接近，波浪力分布沿水深的分布符合波浪理論，可將此計算波浪力的方法推廣到工程實踐中。

本文計算的波浪條件比較理想，計算模型較為簡單，而在工程實踐中，計算模型的簡化以及波浪理論的選取等都會對計算結(jié)果產(chǎn)生影響。因此對于較復(fù)雜的實際問題，還有待更深入的研究。

［1］趙麗娟.構(gòu)筑物與波浪相互作用的研究［J］.水運工程，2003（12）：21-24.ZHAO L J.An Approach to the Interaction between Structures and Wave［J］.Port&Waterway Engineering，2003（12）：21-24.

［2］張玲.小直徑樁柱波浪荷載的計算分析及程序?qū)崿F(xiàn)［D］.青島：中國海洋大學(xué)，2004.

［3］竺艷蓉.海洋工程波浪力學(xué)［M］.天津：天津大學(xué)出版社，1991.

［4］徐偉國.樁柱結(jié)構(gòu)波浪力計算新探［R］.寧波：第三航務(wù)工程局第四工程處，1985.

［5］徐進，石兆東，張康.基于 Morison 方程的海底管道疲勞壽命分析［J］.承德石油高等?？茖W(xué)校學(xué)報，2009，11（4）：21-23.XU J，SHI Z D，ZHANG K.Morison Equation Based Fatigue Life Analysis of Submarine Pipeline［J］.Journal of Chengde Petroleum College，2009，11（4）：21-23.

［6］李慶揚，王能超，易大義.數(shù)值分析［M］.北京：清華大學(xué)出版社，2001.

［7］張寧川，俞聿修.小直徑樁柱上波浪力峰值的短期分布［J］.海洋工程，1993，11（3）：20-27.ZHANG N C，YU Y X.The short-term distribution of the peak value of wave forces on a pile［J］.The Ocean Engineering，1993，11（3）：20-27.

［8］邱大洪.樁群上的最大總波浪力［J］.海洋學(xué)報，1981，3（1）：157-164.QIU D H.The Maximum Total Wave Forces on Pile Groups［J］.Acta Oceanologic Sinica，1981，3（1）：157-164.