龔曉怡,莫建新,李存興
(中交第二航務(wù)工程勘察設(shè)計院有限公司,湖北 武漢 430060)
隨著“一帶一路”倡議的拓展,我國工程企業(yè)越來越多地參與到國際項目的實施中,其中碼頭的設(shè)計是其中一項。對沿海的港口碼頭而言,受當(dāng)?shù)亟?jīng)濟(jì)發(fā)展水平等條件制約,不可能每個項目都建設(shè)防波堤進(jìn)行掩護(hù)。對于沒有掩護(hù)的高樁碼頭結(jié)構(gòu)而言,外海波浪是高樁結(jié)構(gòu)設(shè)計中需要考慮的最重要因素之一,其中波浪對結(jié)構(gòu)的浮托力,不僅影響碼頭面高程的確定,也影響到結(jié)構(gòu)的受力、樁基的選型和布置。如果一味提高碼頭面使之不受波浪力,會導(dǎo)致工程造價顯著增加,并且可能影響到碼頭裝卸使用的方便性;如果結(jié)構(gòu)高程過低,則結(jié)構(gòu)底部將承受一定的波浪力,波浪荷載甚至為樁拔力的控制荷載。波浪浮托力計算是否合理,將直接影響結(jié)構(gòu)的安全性和造價的合理性。
我國水運(yùn)行業(yè)規(guī)范中,有一套基于物理模型試驗的高樁碼頭波浪浮托力計算經(jīng)驗公式,滿足規(guī)范約定的條件時,可用于工程計算。但是多數(shù)海外工程的設(shè)計,合同約定的參考規(guī)范沒有中國規(guī)范,業(yè)主的咨詢工程師也不接受中國規(guī)范,需要選擇合適的規(guī)范和手冊用于計算。同時項目受到條件限制不能開展物理模型試驗時,為了保證結(jié)構(gòu)的安全性,往往需要采用多種規(guī)范計算并與中國規(guī)范手冊進(jìn)行對比。在計算高樁碼頭上部結(jié)構(gòu)浮托力的計算時,不同計算方法得出的結(jié)果不同,因此,選擇合適的規(guī)范顯得尤為重要。
本文對我國水運(yùn)行業(yè)規(guī)范、設(shè)計手冊有關(guān)波浪力對上部結(jié)構(gòu)浮托力的計算方法進(jìn)行了研究,并結(jié)合常用國外規(guī)范有關(guān)計算方法進(jìn)行對比分析,對波浪浮托力計算公式的適用性和合理性進(jìn)行了探討。
根據(jù)《海港設(shè)計手冊》中冊[1],浮托力p采用壓制波理論,可按式(1)汁算:
p=βρg(η-hi)
(1)
式中:η為波峰面在靜水面以上的高度,計算方法見手冊;hi為結(jié)構(gòu)底部的高程;β為壓力反力系數(shù),取1.5~2.0;ρ為水的密度。
PIANCRECOMWG149(part IV)(2017)[2](《國際航運(yùn)協(xié)會碼頭設(shè)計指南(第4部分)》)中,對于碼頭的波浪浮托力也采用壓制波理論,建議β=2.0。
1.1節(jié)和1.2節(jié)均為壓制波理論,本文對比計算時,β統(tǒng)一取值1.5。
根據(jù)JTS 145—2015《港口與航道水文規(guī)范》[3]第10.5節(jié),波浪浮托力計算公式為:
(2)
其中各變量用以下公式計算:
(3)
(4)
且
η1%H1%≤0.7
(5)
(6)
(7)
PF=KFP1%
(8)
(9)
式中:P1%為累積頻率為1% 的波高作用下面板底部縱向單位長度上的最大總浮托力(kNm);ρ為水的密度( tm3);B為沿波浪傳播方向的面板寬度(m);x1%為面板底部的波浪作用寬度(m),即均布壓強(qiáng)分布寬度,當(dāng)x1%>B時取為B;H1%為累積頻率為1%的入射波波高(m);K1為 面板寬度影響系數(shù);Δh為 面板底部在靜水面以上的高度(m);η1%為H1%波高對應(yīng)的波峰在靜水面以上的高度(m);Ls為 有效波波長(m);B為沿波浪傳播方向的面板寬度(m);d為碼頭前沿水深(m);g為重力加速度(ms2);Ts為有效波周期(s);PF為累積頻率為F的波高作用下面板底部縱向單位長度上的最大浮托力(kNm);KF為高樁碼頭面板底部不同累積頻率下的波浪總浮托力換算系數(shù);為累積頻率下的波高作用下面板底部均布壓強(qiáng)。
根據(jù)《海港碼頭設(shè)計手冊》[4],波浪作用壓力可分為靜水壓力和動水壓力,分別按照式(9)和式(10)計算,其中靜水壓力按照波峰面高出結(jié)構(gòu)底面的靜水壓強(qiáng)來考慮。
靜水壓力:
p=ρg(η-hi)
(10)
動水壓力:
(11)
式中:ρ為水的密度;u為水質(zhì)點速度豎向分量。
美國碼頭設(shè)計有關(guān)波浪力主要根據(jù)CoastalEngineeringManual(CEM,《美國陸軍海岸工程手冊》) Part VI-5[5],波浪浮托力可以分解為靜水壓力和動水壓力。
靜水壓力:
(12)
動水壓力:
(13)
式中:FL為向上的托力;ps為結(jié)構(gòu)表面上產(chǎn)生的靜水壓力(不考慮由于波面高度變化引起的靜水壓力差);dA為結(jié)構(gòu)體的表面積微分(其方向向外垂直于表面);n為向上的法向單位矢量;CL為經(jīng)驗上托力系數(shù),與KC數(shù)和雷諾數(shù)有關(guān),KC數(shù)與直徑、波浪周期、水質(zhì)點最大水平速度有關(guān);An為投影面積;ρw為水的密度;g為重力加速度;u為水質(zhì)點速度豎向分量。
1.5和1.4節(jié)計算方法類似,對比計算時兩者統(tǒng)一為靜水壓力和動水壓力方法,按照美標(biāo)CEM計算,系數(shù)CL取特定值1.7。
德國碼頭設(shè)計規(guī)范[6](EAU 2012)沒有給出波浪浮托力的計算,僅僅給出波浪沖擊力的計算,由靜水壓力和動水壓力組成,如圖1所示,計算見式(14)~(16)。
圖1 波浪浮托力(沖擊)定義
波浪沖擊力壓強(qiáng):
p=p0+ps
(14)
靜水壓強(qiáng)分量:
p0=C0ρwg(η-RC)
(15)
動水壓強(qiáng)分量:
Ps=0.5CSρww2
(16)
式中:表面水質(zhì)點向上的速度分量w=csinβ;ρw為水的密度;波速c=LT,為波長與周期的比值;β為接觸點處波浪與碼頭面的夾角;靜水圧力系數(shù)C0取1.0,動水沖擊系數(shù)
(17)
單位寬度的總波浪力:
(18)
其中:l′=min(l,lcr,b),l為接觸區(qū)域的長度,lcr為初始接觸點波峰與碼頭的水平距離,b為平臺寬度。
本節(jié)計算方法為波浪沖擊力的計算,在碼頭面高程較低、波長較長時,得到的波浪沖擊力巨大,因此,原則上應(yīng)避免結(jié)構(gòu)上浪。本節(jié)不參與對比。
根據(jù)日本港工技術(shù)規(guī)范(OCDI)[7],對于行進(jìn)波,波浪上托力計算見式(19)~(22)。
(19)
F*=bwb1p2
(20)
p2=(ηmax-cl)ρg
(21)
(22)
a、b是與結(jié)構(gòu)構(gòu)件有關(guān)的系數(shù)。無掩護(hù)時a=0.82、b=0.61;有掩護(hù)時,對板a=0.71、b=0.71,對梁a=0.82、b=0.61。其余詳細(xì)內(nèi)容見OCDI相關(guān)規(guī)定。
BS 6349-1-2:2016(A1+2017)[8]是英國標(biāo)準(zhǔn)港口工程結(jié)構(gòu)設(shè)計的總則,包含了荷載計算等基礎(chǔ)設(shè)計,但是有關(guān)波浪浮托力的計算沒有直接給出,而是推薦參照文獻(xiàn)[9],計算圖示為圖2,波浪力計算公式為(23)~(24)。
頂部壓力:
p1=[ηmax-(bh+cl)]ρg
(23)
底部壓力:
p2=(ηmax-cl)ρg
(24)
單寬總浮托力為:
Fv=p2bwbl
(25)
平均浮托力公式同式(19),該計算方法與OCDI相近,對比時統(tǒng)一按照OCDI公式。
圖2 波浪浮托力計算圖示
ISO 21650—2007《海岸工程結(jié)構(gòu)波浪與水流作用》[10]有關(guān)波浪浮托力的計算沒有直接給出,而是同樣推薦參照文獻(xiàn)[9]。
印尼某水泥廠配套碼頭項目的重件滾裝碼頭,波浪條件見表1。
對于重件滾裝碼頭工作平臺,受工藝等因素制約,碼頭面高程定為3.5 m,墩臺底部高程定為2.00 m,前沿海底高程為-9.0 m。采用不同規(guī)范手冊進(jìn)行計算,波浪浮托力計算結(jié)果見圖3a)。
越南某海港項目的油碼頭工作平臺,波浪條件見表1。
對于油碼頭工作平臺受工藝裝卸條件等因素制約,碼頭面高程定為6.8 m,墩臺底部高程定為4.80 m,前沿海底高程為-4.0 m。采用不同規(guī)范手冊進(jìn)行計算,波浪浮托力計算結(jié)果見圖3b)。
表1 50 a一遇設(shè)計波要素
注:δhHs為結(jié)構(gòu)底部高出靜水面的高度與有效波高的比值。圖3 波浪浮托力計算結(jié)果對比
對于兩個算例,計算結(jié)果分析如下:1)結(jié)構(gòu)底部高出水位0.4~0.5倍有效波高時,幾種計算方法結(jié)果相對較為接近。2)采用JTS 145—2015標(biāo)準(zhǔn)計算,結(jié)構(gòu)高出水面約0.4~ 0.45倍波高時,波浪浮托力出現(xiàn)峰值;在結(jié)構(gòu)底部高出水位0.2~0.7倍波高區(qū)間時,波浪浮托力相對較大,區(qū)間內(nèi)變化較?。凰贿M(jìn)一步升高時,波浪浮托力隨之減少。3)采用海港設(shè)計手冊壓制波理論,或靜水+動水壓力公式,則隨著面板底部超高δh的減小,波浪浮托力迅速增加,計算結(jié)果明顯大于其他計算方法;而底部超高δh較大時,該方法計算結(jié)果較小。 4)采用OCDI標(biāo)準(zhǔn)進(jìn)行計算,其結(jié)果介于上述規(guī)范之間,水位升高時,波浪浮托力緩緩增加,結(jié)果比JTS 145—2015大,而水位較低時,計算結(jié)果又較JTS 145—2015標(biāo)準(zhǔn)的計算結(jié)果要小。
由于該項目合同中,允許的規(guī)范列表沒有中國規(guī)范,而OCDI規(guī)范計算的波浪浮托力最大值與JTS 145—2015規(guī)范接近,因此設(shè)計計算時采用了OCDI標(biāo)準(zhǔn),既能滿足設(shè)計質(zhì)量控制,又能滿足合同要求。
為了進(jìn)一步分析不同規(guī)范手冊計算之間的差異,本文選取幾個工程常用的波浪情況,進(jìn)行計算對比分析。有效波高為3 m,水深分別取10、25 m,平均周期分別按照5、6、7、9、12 s進(jìn)行計算,有效周期取平均周期的1.15倍,相應(yīng)的有效波長根據(jù)有效周期和對應(yīng)水深進(jìn)行計算,1%波高取H1%=1.5Hs,最大波高取Hmax=1.8Hs。
不同水深、不同波周期、不同規(guī)范(手冊)波浪浮托力的計算對比見表2。
根據(jù)表2分析可知:1)其余條件相同,水深越深,波浪浮托力越小;水深越淺,浮托力越大。2)其余條件相同,波周期越短,波浪浮托力越??;波周期越長,浮托力越大。3)總體而言,在結(jié)構(gòu)底部高出水面0.3~0.5倍波高時,幾種計算方法的結(jié)果相對比較接近,但波浪周期很長或很短時,幾種方法沒有明顯的接近值。4)其他條件相同,在結(jié)構(gòu)底部高程較低(水位較高)時,幾種計算方法結(jié)果大致為:海港設(shè)計手冊波壓制理論> 靜水+動水總壓力方法> OCDI標(biāo)準(zhǔn)(無掩護(hù))> OCDI標(biāo)準(zhǔn)(有掩護(hù))> 海港水文規(guī)范JTS 145—2015。5)其他條件相同,在結(jié)構(gòu)底部高程較高(水位較低)時,除海港設(shè)計手冊波壓制理論計算結(jié)果最小,其余幾種方案沒有明顯規(guī)律。6)對于較短周期波浪或陡波,海港水位規(guī)范JTS 145—2015計算結(jié)果相對較小,OCDI等計算結(jié)果相對較大。7)對于較長周期波浪或坦波,在結(jié)構(gòu)底面較高時,海港水位規(guī)范JTS 145—2015計算結(jié)果相對較大,其他方法計算結(jié)果較小。
表2 工程常用波浪要素下不同規(guī)范(手冊)的波浪力浮托力計算對比
1)壓制波理論方法,在結(jié)構(gòu)底面超高較大時,浮托力計算結(jié)果最??;但結(jié)構(gòu)底面超高較小(水位附近)時,計算結(jié)果最大。波壓力反力經(jīng)驗參數(shù)確定困難,不宜作為計算波浪浮托力的首選。
2)靜水壓力+動水壓力方法,對于短周期波浮托力結(jié)果偏大,長周期波偏小。與其他方法結(jié)果相差較大,也不宜做計算首選。
3)JTS 145 —2015公式作為現(xiàn)行行業(yè)規(guī)范,在其適用范圍內(nèi),可以作為設(shè)計首選;但是超出其使用范圍時,應(yīng)慎重。
4)OCDI規(guī)范公式,與其他國際規(guī)范計算方法基本一致,在國際上有一定影響力,對于國外工程,如果合同不能采用中國規(guī)范,則OCDI是良好的替代規(guī)范。