張寶華，王元戰(zhàn)，孫建軍

（1.交通部天津水運(yùn)工程科學(xué)研究所水工構(gòu)造物檢測(cè)、診斷與加固技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室，天津300456；2.天津大學(xué)建筑工程學(xué)院，天津 300072）

拋石基床是重力式碼頭或直立式防波堤的基礎(chǔ)。近年來(lái)，作業(yè)船舶的大型化、高速化和操作性能的提高給港口工程建設(shè)和運(yùn)營(yíng)管理提出了新的課題，如大型拖輪和滾裝船螺旋槳產(chǎn)生的水流對(duì)水工建筑物拋石基床的沖刷問(wèn)題。目前，重力式碼頭的拋石基床出現(xiàn)了被大水流掏空的重大安全隱患。

20世紀(jì)70年代以來(lái)，伴隨著引擎功率大幅提升，我國(guó)到港船舶也日趨大型化。特別是近年來(lái)，各港口都在增大船載噸位，大多需要大型拖輪協(xié)助調(diào)頭或靠泊碼頭作業(yè)。因而，各港口的港作拖輪馬力越來(lái)越大，螺旋槳的推力及尾流流速也越來(lái)越大。大型船舶尾流流速超過(guò)了10～100 kg塊石的起動(dòng)流速，導(dǎo)致螺旋槳射流沖刷基床而損壞碼頭。然而，設(shè)計(jì)規(guī)范未涉及港作大拖輪或尾直式滾裝船對(duì)拋石基床沖刷破壞的工況，也沒(méi)有相應(yīng)防范措施。規(guī)范的不足客觀上造成了設(shè)計(jì)的缺陷，因此有必要對(duì)此展開(kāi)研究。

1 有限元模型及其可行性研究

1.1 二維有限元模型的建立

塊石的沖刷穩(wěn)定性問(wèn)題實(shí)質(zhì)上是塊石的穩(wěn)定重量問(wèn)題。目前，國(guó)內(nèi)外相關(guān)的研究工作有很多［1-5］，其中比較著名的有伊茲巴什公式［6］、沙莫夫公式［7］、交通部防波堤施工規(guī)范［8］等。但是絕大多數(shù)研究工作都基于理論推導(dǎo)和物理模型試驗(yàn)展開(kāi)，很少有人通過(guò)建立有限元數(shù)值模型進(jìn)行研究。本文首次通過(guò)建立有限元模型，對(duì)船舶尾流作用下拋石基床塊石穩(wěn)定性問(wèn)題展開(kāi)研究。

根據(jù)等容原則，將形狀不規(guī)則的塊石簡(jiǎn)化成同體積的圓球后，再根據(jù)等容原則將體積為V、直徑為D的圓球，簡(jiǎn)化成為體積為V、長(zhǎng)度為D的圓柱體。這樣使圓柱體的高度方向與斷面的法線方向一致排布，就將一排圓球簡(jiǎn)化成一排圓柱體，簡(jiǎn)化過(guò)程見(jiàn)圖1。在理想狀態(tài)下，由于每個(gè)小圓柱體重量相等，只有受到相同水流力作用時(shí)，才可能發(fā)生位移，也就是說(shuō)當(dāng)水流達(dá)到一定條件后，這一排圓柱體石塊將同時(shí)起動(dòng)。可以將這一排圓柱體石塊看作是一根中間連續(xù)的石柱。這樣就將三維的石球簡(jiǎn)化成了石柱，取其斷面進(jìn)行二維模型分析，有限元模型見(jiàn)圖2。

1.2 有限元模型的可行性研究

為了驗(yàn)證本文所提出二維模型的可行性，首先建立了水平面上3種排布狀態(tài)（圖3）的單層塊石，分別進(jìn)行計(jì)算分析。將所得起動(dòng)流速計(jì)算結(jié)果與伊茲巴什公式、交通部規(guī)范擬合公式做比較，比對(duì)結(jié)果見(jiàn)圖4。

由圖4可以看出，本文有限元計(jì)算結(jié)果與伊茲巴什公式等理論計(jì)算結(jié)果相比，在起動(dòng)流速隨穩(wěn)定重量的變化趨勢(shì)上是一致的，且計(jì)算結(jié)果與以往理論計(jì)算結(jié)果較為吻合。說(shuō)明本文所提出的二維數(shù)值模型是可行的，也驗(yàn)證了通過(guò)有限元建立的數(shù)值模型在研究塊石起動(dòng)流速問(wèn)題上的可行性。

圖1 模型的簡(jiǎn)化過(guò)程圖Fig.1 Simplified model

圖2 拋石基床遭尾流沖刷數(shù)值模型圖Fig.2 Numerical model of rubble bed scoured by ship stern flow

圖3 3種排布狀態(tài)模型圖Fig.3 Three assignments of configuration

圖4 有限元計(jì)算結(jié)果與以往理論成果比較圖Fig.4 Comparison between finite element results and previous theoretical results

2 船舶尾流數(shù)值模擬

目前，國(guó)內(nèi)外在船舶尾流領(lǐng)域的研究已經(jīng)取得了豐碩的成果。本文針對(duì)重力式碼頭明基床形式，分別計(jì)算螺旋槳位于岸壁前3 m、5 m、7.5 m和10 m時(shí)，距離海底1.5 m、2 m、2.5 m共12種情況，運(yùn)用FLOW-3D模擬計(jì)算螺旋槳尾流場(chǎng)分布狀況，其數(shù)值模型及部分計(jì)算結(jié)果見(jiàn)圖5～圖6。

圖5 螺旋槳射流數(shù)值模型圖Fig.5 Numerical model of propeller jet

圖6 槳位于岸壁前5 m、離地1.5 m結(jié)果圖Fig.6 Numerical model of propeller（5 m offshore，1.5 m above the ground）

計(jì)算結(jié)果表明：當(dāng)螺旋槳距離岸壁結(jié)構(gòu)3 m，螺旋槳距離海底1.5 m時(shí)，射流主要沿著拋石基床表面?zhèn)鞑?，遇到岸壁結(jié)構(gòu)后，沿岸壁結(jié)構(gòu)物向上傳播，回流對(duì)拋石基床的沖刷作用不明顯；而當(dāng)螺旋槳距離岸壁結(jié)構(gòu)3 m，螺旋槳距離海底2.0 m和2.5 m時(shí)，螺旋槳射流遇到岸壁結(jié)構(gòu)物后反射，部分發(fā)射回流沿斜向下方向傳播，這部分回流對(duì)基床的沖刷作用顯著。當(dāng)螺旋槳距離岸壁5 m時(shí)，螺旋槳射流遇到岸壁的反射回流依然存在，只是這部分回流主要分布在離岸壁較近的區(qū)域，此時(shí)拋石基床堤頭主要受到螺旋槳射流的直接沖刷作用。當(dāng)螺旋槳距離岸壁結(jié)構(gòu)7.5 m時(shí)，射流場(chǎng)分布與5 m距離時(shí)相比，直接沖刷作用和回流作用都明顯減弱。當(dāng)螺旋槳距離岸壁結(jié)構(gòu)10 m，螺旋槳尾流傳播到基床位置時(shí)，已經(jīng)逐漸趨于平穩(wěn)，基本為水平流向，渦旋作用減弱，此時(shí)螺旋槳離海底的距離對(duì)基床處流場(chǎng)分布的影響已經(jīng)很弱。

對(duì)沿岸壁方向流速分布進(jìn)行分析，結(jié)果表明：當(dāng)螺旋槳距離岸壁結(jié)構(gòu)較近時(shí)，螺旋槳中軸線附近的流速較大；當(dāng)螺旋槳離開(kāi)岸壁一段距離后，流速較大的區(qū)域開(kāi)始向中軸線兩側(cè)偏移，此時(shí)中軸線處的流速相對(duì)小一些；當(dāng)螺旋槳離岸壁較遠(yuǎn)時(shí)，尾流流場(chǎng)已經(jīng)逐漸趨于穩(wěn)定，此時(shí)在拋石基床堤頭處的水流分布變化不大。這樣，當(dāng)給出螺旋槳出流的初速度時(shí)，通過(guò)流速比可以粗略推算出螺旋槳距離岸壁一定距離時(shí)，螺旋槳尾流在拋石基床堤頭處的流速分布狀態(tài)，可用于針對(duì)性的有限元建模計(jì)算，或作為實(shí)際工程安全穩(wěn)定性評(píng)估的參考。

研究發(fā)現(xiàn)，當(dāng)螺旋槳距離岸壁結(jié)構(gòu)物約5 m時(shí)（當(dāng)射流初始速度較大時(shí)，該區(qū)域也會(huì)相應(yīng)增大），且距離海底在1.5 m以上（拋石基床厚度1.0～1.2 m），此時(shí)螺旋槳射流遇到岸壁結(jié)構(gòu)物的反射回流對(duì)拋石基床的沖刷作用顯著，這是拋石基床被大水流掏空的主要原因。當(dāng)螺旋槳距離岸壁結(jié)構(gòu)較遠(yuǎn)時(shí)，螺旋槳射流對(duì)拋石基床的沖刷作用減弱，此時(shí)水流主要沿著拋石基床斜坡面向岸傳播，遇到岸壁結(jié)構(gòu)后產(chǎn)生向上或斜向上的反射作用?？梢?jiàn)，在船舶靠泊和起航的時(shí)候，螺旋槳射流對(duì)基床的沖刷是造成拋石基床沖刷破壞的重要因素，應(yīng)針對(duì)這種工況重點(diǎn)對(duì)拋石基床進(jìn)行設(shè)計(jì)防范。

3 拋石基床遭沖刷穩(wěn)定性研究

圖7 30 kg塊石在2.3 m/s流速條件下起動(dòng)結(jié)果圖Fig.7 Incipient motion for 30 kg stone blocks under 2.3 m/s flow condition

圖8 30 kg塊石在2.5 m/s水平流和2 m/s斜向流共同作用下起動(dòng)結(jié)果圖Fig.8 Incipient motion for 30 kg stone blocks under 2.5 m/s horizontal flow and 2 m/s oblique flow condition

圖9 30 kg塊石在1.3 m/s回流作用下起動(dòng)結(jié)果圖Fig.9 Incipient motion for 30 kg stone blocks under 1.3 m/s backflow condition

船舶尾流流場(chǎng)分布極其復(fù)雜，本文結(jié)合二維模型的斷面結(jié)構(gòu)形式，將復(fù)雜的流場(chǎng)分布簡(jiǎn)化為二維模型中的流速荷載。結(jié)合三維流場(chǎng)中流速在基床面上的分布規(guī)律，忽略徑向流速（即y向流速）的影響，針對(duì)不同螺旋槳位置時(shí)刻的流場(chǎng)分布，選取20 kg、30 kg、40 kg、50 kg 和 75 kg 5 種不同質(zhì)量的塊石建立數(shù)值模型，將流場(chǎng)簡(jiǎn)化為以下3種流速荷載作用情況進(jìn)行研究分析：（1）水平來(lái)流作用下塊石的起動(dòng)（圖7）。這種情況模擬的是螺旋槳距岸壁結(jié)構(gòu)一段距離且在特定深度時(shí)，塊石受到螺旋槳尾流直接作用，發(fā)生沖刷位移的情況；（2）水平流和沿堤頭坡面向上水流共同作用下堤頭塊石的起動(dòng)（圖8）。這種情況模擬的是螺旋槳距岸壁結(jié)構(gòu)一段距離后，螺旋槳射流遇到岸壁結(jié)構(gòu)的回流作用已經(jīng)很弱，不是塊石起動(dòng)的主導(dǎo)作用時(shí)，塊石受到螺旋槳尾流直接作用，發(fā)生沖刷位移的情況；（3）斜向下回流作用下堤頭塊石的起動(dòng)（圖9）。這種情況模擬的是當(dāng)螺旋槳距岸壁結(jié)構(gòu)距離很近，螺旋槳射流遇到岸壁結(jié)構(gòu)的反射回流流速很大，這部分回流在距離岸壁較近的區(qū)域，是斜向下作用在拋石基床表層塊石上的，這時(shí)候堤頭部分的塊石受到斜向下方向水流作用。

將上述拋石基床二維模型的塊石在水平來(lái)流和回流作用下起動(dòng)的有限元計(jì)算結(jié)果，與水平面上單層塊石排布的計(jì)算結(jié)果作比較，對(duì)照結(jié)果見(jiàn)圖10和圖11。

圖10 水平來(lái)流作用下水平面單層塊石與拋石基床堤頭塊石起動(dòng)流速對(duì)照?qǐng)DFig.10 Comparison diagram of stone incipient velocity between horizontal layer and rubble bed under horizontal flow condition

圖11 回流作用下水平面單層塊石與拋石基床堤頭塊石起動(dòng)流速對(duì)照?qǐng)DFig.11 Comparison diagram of stone incipient velocity between horizontal layer and rubble bed under backflow condition

通過(guò)計(jì)算分析拋石基床結(jié)構(gòu)型式下不同水流作用時(shí)塊石的起動(dòng)流速發(fā)現(xiàn)：螺旋槳射流作用下，沿射流方向只有水平流速作用時(shí)，塊石的起動(dòng)流速比水平面上單層塊石的起動(dòng)流速大；沿射流方向，在水平流和沿基床坡面向上的斜向流兩股水流共同作用時(shí)，基床堤頭塊石反而更不易起動(dòng)，塊石的穩(wěn)定性相對(duì)提高。

通過(guò)計(jì)算斜向下的回流作用下基床堤頭塊石的起動(dòng)發(fā)現(xiàn)：回流作用下，塊石的起動(dòng)流速比同種情況水平流作用下小10%以上。也就是說(shuō)，在螺旋槳射流遇岸壁結(jié)構(gòu)的反射回流作用下，塊石更容易發(fā)生失穩(wěn)起動(dòng)，這種情況下的沖刷作用相當(dāng)明顯，從而解釋了目前許多重力式碼頭出現(xiàn)的拋石基床遭船舶尾流沖刷，塊石被水流掏空的現(xiàn)象。針對(duì)回流作用下塊石易發(fā)生失穩(wěn)的情況，建議以后的研究工作可以針對(duì)回流的流場(chǎng)分布特點(diǎn)，以及回流對(duì)塊石沖刷的作用機(jī)理，重點(diǎn)展開(kāi)研究，提出新的加固方案，提高拋石基床的沖刷穩(wěn)定性，為設(shè)計(jì)設(shè)防工作提供有力參考。

4 小結(jié)

采用有限元軟件Flow-3D建立數(shù)值模型，模擬碼頭岸壁附近船舶螺旋槳尾流流場(chǎng)分布情況。結(jié)果表明，當(dāng)螺旋槳距離岸壁結(jié)構(gòu)物5 m左右（當(dāng)射流初始速度較大時(shí)，該區(qū)域也會(huì)相應(yīng)增大），且距離海底在1.5 m以上（拋石基床厚度1.0～1.2 m）時(shí)，螺旋槳射流遇到岸壁結(jié)構(gòu)物的反射回流對(duì)拋石基床的沖刷作用顯著，這是造成拋石基床被大水流掏空的主要因素。

運(yùn)用ADINA軟件，針對(duì)幾種較容易引起塊石起動(dòng)的流速條件進(jìn)行計(jì)算分析，發(fā)現(xiàn)螺旋槳射流作用下，沿著螺旋槳射流方向，當(dāng)水平流速起主導(dǎo)作用時(shí)，塊石的起動(dòng)流速略大于水平面上單層塊石的起動(dòng)流速；當(dāng)塊石在受到水平向和沿堤頭斜坡面向上的兩股水流共同作用時(shí)，塊石的起動(dòng)流速相對(duì)較大，其穩(wěn)定性相對(duì)較高，通過(guò)分析其流速分布發(fā)現(xiàn)，沿堤頭斜坡面向上的水流使水平流的流向發(fā)生了改變，對(duì)堤頭塊石的沖擊作用減弱，堤頭塊石在這兩股水流共同作用下，處于一種更穩(wěn)定的狀態(tài)。通過(guò)計(jì)算塊石在斜向下方向的回流作用下的起動(dòng)情況發(fā)現(xiàn)，同種情況下塊石在回流作用下比水平面上更易起動(dòng)，起動(dòng)流速低10%以上。這一結(jié)論解釋了目前許多重力式碼頭拋石基床，遭到船舶尾流沖刷破壞，出現(xiàn)拋石基床塊石被大水流掏空的工況。

［1］Headquarters US Army Corps of Engineers.Engineer Manual：Hydraulic，Design of Flood Control Channels［M］.Washington D C：Headquarters US Army Corps of Engineers，1994.

［2］Bakker K J，Verheij H J，Groot M B.Design Relationship for Filters in Bed Protection［J］.Journal of Hydraulic Engineering，1994，120（9）：1 082-1 088.

［3］吳蘇舒，張瑋.海岸防護(hù)中護(hù)底塊石穩(wěn)定重量的確定［J］.中國(guó)港灣建設(shè)，2007（6）：13-17.WU S S，ZHANG W.Determination of Stable Weight of Bottom Blocks in Coastal Protection Works［J］.China Harbour Engineering，2007（6）：13-17.

［4］張瑋，沈正潮，鐘春欣.涌潮作用下護(hù)灘塊石的起動(dòng)研究［J］.泥沙研究，2007（1）：24-29.ZHANG W，SHEN Z C，ZHONG C X.Incipient Motion of Rock Blocks for Beach Protection Under the Action of the Tidal Bores［J］.Journal of Sediment Research，2007（1）：24-29.

［5］陳小莉，馬吉明.受渦旋作用的水下塊石的起動(dòng)流速［J］.清華大學(xué)學(xué)報(bào)，2005，45（3）：315-318.CHEN X L，MA J M.Critical Velocity for Initiation Motion of Rocks Under Water Due to Vortex Motion［J］.Journal of Tsinghua University，2005，45（3）：315-318.

［6］陳吉余.中國(guó)圍海工程［M］.北京：中國(guó)水利水電出版社，2000.

［7］王昌杰.河流動(dòng)力學(xué)［M］.北京：人民交通出版社，2001.

［8］JTJ298-98，防波堤設(shè)計(jì)與施工規(guī)范［S］.