王勝勇，王冬姣*，季勇志

（1.華南理工大學土木與交通學院，廣東廣州510640；2.中交廣州航道局有限公司，廣東廣州510221）

絞吸船水流力的數(shù)值模擬研究

王勝勇1，王冬姣1*，季勇志2

（1.華南理工大學土木與交通學院，廣東廣州510640；2.中交廣州航道局有限公司，廣東廣州510221）

利用CFD方法對7025型絞吸船遭遇不同來流攻角時的周圍流場進行了數(shù)值模擬，得出作用在絞吸船上的水流力系數(shù)隨來流攻角和水深的變化關系。結果表明：水深對橋架縱向水流力影響較大而對橫向水流力和縱搖力矩系數(shù)的影響較小。結果可為絞吸船施工作業(yè)提供指導。

絞吸船；水流力系數(shù)；數(shù)值模擬

0 引言

絞吸式挖泥船是一種結構復雜、技術含量高的水上施工設備[1]。挖泥船施工作業(yè)時，變化的水流力會影響船舶的操控甚至可能會危及船舶的安全[2]。在作業(yè)過程中，船體所受各種載荷非常復雜，如絞刀切削力及風、浪、流對船體的作用力等。目前大多采用OCIMF[3]提供的VLCC水流力系數(shù)的試驗曲線計算作用在絞吸船上的水流載荷，沒有考慮首尾凹槽對流載荷的影響[1，4-5]。由于絞吸船上有從甲板貫穿到船底的凹槽，與一般流線型運輸船有截然不同的阻力性能[6]，因此最好能夠根據(jù)絞吸船的實際線型確定挖泥作業(yè)時船體所受到的流載荷。分析作用在船體上的流載荷，對于保證絞吸船施工安全，提高作業(yè)能力至關重要[7]。在挖泥工況，受流載荷作用的有船體、定位樁和絞刀橋架。作用在橋架上的流載荷將通過船體傳遞到定位樁上，因此本文主要討論絞吸船船體和鉸刀橋架所受到的流載荷。應用Fluent軟件中VOF模型分析來流攻角和水深吃水比對絞吸船水流力系數(shù)的影響。7025型絞吸船允許施工的設計工況條件為最大風速15 m/s、最大流速2 m/s、最大波高1 m，波浪周期5 s，絞吸船船尾必須與波浪或水流的方向一致，挖寬約為80 m，擺角約66°[5]，因此絞吸船遭遇尾流作用時其來流攻角范圍約為-33°≤θc≤33°（挖深5 m，鋼樁位于船尾）。利用數(shù)值模擬得到的水流力系數(shù)可作為7025型絞吸船挖泥工況環(huán)境載荷中水流力系數(shù)的輸入依據(jù)，對施工設計具有指導意義。

1 控制方程

假定船舶靜止不動，水流勻速流過船體。將ox軸取在船體縱中剖面，指向船首，平行于水線面；oy軸與縱中剖面垂直，指向右舷；oz軸垂直向上為正。同時假定流體是不可壓縮的，任何流動都要滿足質量守恒方程：

和動量方程

式中：ρ為流體密度；p為流體壓強；μ為流體動力黏性系數(shù)；ρg為重力；F 為外力項；V為流體速度。

自由液面采用VOF法模擬，體積分數(shù)αq指單元內(nèi)第q相流體所占體積與該單元總體積之比。本文采用氣液兩相模擬，αq滿足如下方程：

可定義空氣為基本相，水為第二相，基本相的體積分數(shù)方程不需要求解，而從式（4）得到：

定義αq=0.5為自由液面，流場中每個單元的流體密度ρ由式（5）確定。

每個單元內(nèi)的動力黏性系數(shù)μ，與密度的定義一致。

設作用在船體上的流載荷在x方向的分量為Fx，y方向的分量為Fy，相對于重心的艏搖力矩為Mxy，則水流力系數(shù)可表示為[3]：

式中：L為船長；T為吃水；ρw為水的密度；V為水流速度。

2 計算模型驗證

取長4 m，寬0.843 m，型深0.256 m，吃水0.243 m的某沉管實驗模型，拖曳水池寬度為6 m，實驗段水深0.322 m[8]。在數(shù)值模擬時取一半流場進行分析，取計算域x=-6~10 m，y=0~3 m，z= -0.322~0.2 m。邊界條件設置為：x=-6 m處為壓力入口，x=10 m處為壓力出口，y=0及z= 0.2 m處為對稱邊界，在沉管表面、池壁y=3 m處及池底z=-0.322 m處均為壁面邊界。全流場采用六面體結構化網(wǎng)格，一半流場生成的網(wǎng)格總數(shù)為38.28萬，選取RNG k-ε模型，在浮運縱拖工況利用數(shù)值模擬得到的沉管管段模型靜水阻力曲線與模型試驗結果的比較見圖1，由圖可知，兩者吻合良好。

圖1 沉管管段模型縱拖阻力曲線Fig.1Longitudinal towing resistance of an immersed tube element

3 數(shù)值模擬

取7025型絞吸船為研究對象進行建模計算，其主尺度為：船長66.5 m，型寬15.4 m，型深4 m，設計滿載吃水2.68 m，艏開槽（長×寬）24.5 m× 5 m，艉開槽（長×寬）12.6 m×4 m，船底艏部結構向上傾斜，重心距船尾30.95 m，其工作水深為5~25 m，來流攻角定義為來流速度與x軸正向的夾角。

流體計算域x軸方向從艏、艉向上、下游各延伸2倍船長，y軸方向從舷側左右方向各延長2倍船長，空氣域上界面距離水面為5 m。使用ANSYS中Design Model模塊建立模型并使用ICEM網(wǎng)格劃分軟件對計算域進行全結構化網(wǎng)格劃分，并對船體附近區(qū)域網(wǎng)格加密。采用速度入口（Velocity Inlet），壓力出口（Pressure Outlet）邊界條件，船體和海底采用壁面（Wall）邊界條件，空氣流體域上表面采用對稱（Symmetry）邊界條件。

在挖泥工況，受流載荷作用的除船體外、還有定位樁和絞刀橋架，作用在橋架上的水流力也將通過船體傳遞到定位樁上。絞吸船船體和橋架的幾何模型見圖2。

圖2 船體和橋架幾何模型Fig.2Geometric model of barge and ladder

7025型絞吸船的橋架外形為方柱狀，橋架面板上有許多開孔，為簡化計算，建模時沒有考慮橋架面板上的開孔，橋架面板上的開孔會使作用在橋架上的水流載荷有所減小，因此參照文獻[9]，取折減系數(shù)0.75來考慮橋架面板上開孔對橋架流載荷的影響。

4 結果分析

4.1 絞吸船船體水流力系數(shù)

取流速2 m/s，計算了絞吸船船體在0°~180°來流攻角范圍，3種水深d=5 m（d/T=1.866）、8.04 m（d/T=3）、16.08 m（d/T=6）對應的流載荷，水流力系數(shù)曲線如圖3所示，其中艏搖力矩系數(shù)相對于重心。

圖3 艏艉開槽船體上的水流力系數(shù)Fig.3Current forces coefficients on barge with notches at fore and aft

由圖3可知，由于絞吸船首部船底有一斜升角，使絞吸船遭遇頂流作用時的縱向流載荷小于尾流工況。在0°~180°來流攻角范圍內(nèi)，絞吸船船體的縱、橫向水流力系數(shù)、艏搖力矩系數(shù)隨著水深吃水比的增加均呈明顯下降趨勢。這表明，水深對流載荷影響很大，淺水會導致絞吸船所受流載荷增加。

水深吃水比d/T=3時非流線型絞吸船船體和VLCC隨來流攻角變化的水流力系數(shù)對比如圖4所示，其中艏搖力矩相對于船舯。

圖4 艏艉開槽船體和VLCC流載荷比較Fig.4Comparison of current loads on barge with notches at fore and aft and on VLCC

由圖4可知，絞吸船船體和VLCC縱向水流力系數(shù)差距較大，這是由于7025型絞吸船采用非流線型設計，水流阻力大，而VLCC為流線型船，水流阻力小。絞吸船和VLCC橫向水流力系數(shù)和艏搖水流力系數(shù)變化規(guī)律一致。

4.2 橋架對水流力的影響

在挖泥工況，整個橋架幾乎都淹沒在水中，橋架對流載荷有較大影響。圖2所示幾何模型在5 m和16.08 m工作水深包含橋架在內(nèi)的水流力系數(shù)如圖5和圖6所示。

圖5 絞吸船流載荷系數(shù)（d=5 m）Fig.5Coefficients of current forces and moment on cutter suction dredger（d=5 m）

圖6 絞吸船流載荷系數(shù)（d/T=6，d=16.08 m）Fig.6Coefficients of current forces and moment on cutter suction dredger（d/T=6，d=16.08 m）

由圖5和圖6可知，隨著作業(yè)水深的增加，對橋架橫向水流力系數(shù)和艏搖力矩系數(shù)的影響較小，但橋架在橫剖面上的投影面積隨之增加，導致作用在橋架上的縱向水流力隨之增加。作業(yè)水深較大時作用在橋架的縱向水流力將大于船體上的水流力。因此，在施工作業(yè)中，要特別注意橋架產(chǎn)生的水流力隨水深的變化。

5 結語

挖泥工況需同時考慮風、浪、流環(huán)境載荷及定位樁、絞刀架、橫移鋼絲繩、絞刀等對絞吸船運動的影響。為了得到較為精確的結果，有必要利用絞吸船的實際線型確定作用在船體上的環(huán)境載荷。本文采用Fluent軟件中RNG k-ε模型對7025型絞吸船對不同來流攻角、不同水深吃水比時船體周圍的流場進行了數(shù)值模擬計算，分析了橋架對流載荷的影響。主要結論如下：

1）絞吸船船體部分所受水流力系數(shù)大小與來流攻角及作業(yè)水深有關。橫向流載荷遠大于縱向，淺水會導致橫向流載荷大幅增加。

2）傳遞到定位樁上的流載荷除來自船體之外、還有絞刀橋架。橋架在橫剖面上的投影面積隨水深而增加，導致作用在橋架上的縱向流載荷隨作業(yè)水深增加，甚至超過船體上的縱向水流力；作用在橋架上的橫向水流力和艏搖力矩隨作業(yè)水深變化的影響較小。

[1]張文宜.絞吸式挖泥船主要載荷分析和計算軟件開發(fā)[D].上海：上海交通大學，2009. ZHANG Wen-yi.Main loads analysis for cutter suction dredger and compilation of calculation software[D].Shanghai:Shanghai JiaoTong University,2009.

[2]柯軍.絞吸式挖泥船三錨定位系統(tǒng)的設計[D].武漢：武漢理工大學.2007. KE Jun.The design of"The Christmas Tree"system on a cutter suction dredger[D].Wuhan：Wuhan Universityof Technology,2007.

[3]OCIMF.Prediction of wind and current loads on VLCCs[Z].1994.

[4]陳新權，楊啟，顧敏童，等.絞吸挖泥船環(huán)境載荷及其對作業(yè)能力的影響[J].中國港灣建設，2012（6）：1-3. CHEN Xin-quan,YANG Qi,GU Ming-tong,et al.Environment load and its effects on capability of cutter suction dredger[J].China Harbour Engineering.2012(6):1-3.

[5]季勇志，吳秀艷，吳小明.風浪流對絞吸船定位樁的沖擊研究[J].水運工程，2013（2）：31-36. JI Yong-zhi,WU Xiu-yan,WU Xiao-ming.On spud impacted by wind,wave and current[J].Port and Waterway Engineering,2013(2): 31-36.

[6]石珣，陳新權，譚家華.基于計算流體力學方法的水線面大開口工程船阻力性能分析[J].上海交通大學學報，2012，46（8）：1 178-1 183. SHI Xun,CHEN Xin-quan,TAN Jia-hua.Study on resistance performance of working ship with notches based on CFD[J].Journal of Shanghai Jiaotong University,2012,46(8):1 178-1 183.

[7]張盼盼，徐立群，倪福生.絞吸挖泥船水流作用力數(shù)值模擬[J].機械制造與自動化.2013.42（3）：139-142. ZHANG Pan-pan,XU Li-qun,NI Fu-sheng.Numerical simulation of current drag force on cutter suction dredger[J].Machine Building &Automation,2013,42(3):139-142.

[8]中國交通建設股份有限公司，中交第四航務工程局有限公司，華南理工大學.港珠澳大橋主體工程初步設計——沉管管段浮運沉放的物理模型試驗研究專題研究報告[R].2009. China Communications Construction Co.,Ltd.,CCCC Fourth Harbor Engineering Co.,Ltd.,South China University of Technology. Preliminary design of the main project of the Hongkong Zhuhai Macao Bridge:A monographic research report on the physical model test of the floating and sinking of immersed tube section[R]. 2009.

[9]季勇志.絞吸船定位樁抗沖擊能力的數(shù)值研究[J].中國港灣建設，2011（5）：16-20. JI Yong-zhi.Numerical research of anti-impact capacity of spud of a cutter suction dredger[J].China Harbour Engineering,2011(5): 16-20.

Numerical simulation of current loads on cutter suction dredge

WANG Sheng-yong1,WANG Dong-jiao1*,JI Yong-zhi2
（1.School of Civil Engineering and Transportation,South China University of Technology,Guangzhou,Guangdong 510640, China;2.CCCC Guangzhou Dredging Co.,Ltd.,Guangzhou,Guangdong 510221,China）

The CFD method was used to simulate the flow around the 7025 type cutter suction dredger at various angles of attack with the incoming flow.The relationship of the current force coefficients on the cutter suction dredger with current angle of attack and the water depth are obtained.It is shown that the water depth has great influence on the longitudinal current force of the ladder,but the influence on the lateral force and current yaw moment is relatively small.The conclusion can provide guidance for dredging operation of cutter suction dredger.

cutter suction dredger;current force coefficient;numerical simulation

U615.35

A

2095-7874（2017）05-0010-05

10.7640/zggwjs201705003

2016-11-19

2017-01-16

廣東省科技計劃項目（2014A020217001）

王勝勇（1991—），男，山西運城人，碩士研究生，船舶與海洋工程專業(yè)。

*通訊作者：王冬姣，E-mail：djwang@scut.edu.cn