基于ANSYS-AQWA的墊擋船碼頭靠泊特性分析

袁培銀,李渝鋒,趙 宇,張 哲

(1. 重慶交通大學 航運與船舶工程學院,重慶 400074;2. 重慶交通大學 建筑與城市規(guī)劃學院,重慶 400074)

0 引 言

隨著經(jīng)濟迅猛發(fā)展,三峽庫區(qū)通航船舶數(shù)量逐步增加,碼頭附近靠泊作業(yè)的船舶噸位陡增。往來船只眾多,?？看芭c碼頭發(fā)生碰撞事故時有發(fā)生,碰撞事故不僅會增加了碼頭維修費用,還會對船員和碼頭工作人員安全產(chǎn)生隱患。三峽大壩建成后,三峽庫區(qū)水位上升,水面變寬,水域面積增大;另外,水位抬高導致水深和過流面積的大幅度增加,水流明顯變緩[1]。隨著航道優(yōu)化,三峽庫區(qū)的船舶也逐步向標準化、大型化的方向發(fā)展。因此,開展三峽庫區(qū)系泊船舶水動力性能分析,對于工程建設、防災減災具有十分重要意義。

學界對碼頭系泊船舶運動特性進行了相關的研究。在考慮影響船舶碼頭系泊系統(tǒng)的外界因素時,不能只考慮某一單個因素,影響船舶撞擊的環(huán)境因素應考慮風、流、波浪和水面高度等[2]。高峰等[3]對碼頭系泊過程中的環(huán)境條件、船體運動及護舷碰撞力進行了模擬計算,發(fā)現(xiàn)運動響應幅值、纜繩張力和碰撞力會隨著頻率變化而改變;胡毅等[4]利用多體水動力計算軟件對船碼頭系泊時的運動特性進行了研究,分析了LNG船碼頭系泊時的整體運動響應及系纜繩所受張力變化規(guī)律;周豐等[5]設計了船舶碼頭系泊系統(tǒng),優(yōu)化了系泊系統(tǒng)方式,使得船舶碼頭系泊系統(tǒng)的特征量計算更加地高效;R.GRANT等[6]通過數(shù)值計算,研究了碼頭系泊船舶六自由度運動響應與波浪載荷的關系;TAN Huiming等[7]利用模型試驗方法,研究了在風浪流作用下的泊位長度對LNG船舶系泊影響,試驗表明:波浪力是影響船舶運動的最大因素,泊位長度對船舶運動影響不明顯,但較短的泊位長度有助于降低纜索的張力;P.ROSA-SANTOS等[8]從減少系泊船舶運動、改善泊位操作和安全條件方面出發(fā),分析了護舷類型對系泊纜張力的影響;I.TOUZON等[9]利用運動學方程,模擬了支撐結(jié)構(gòu)與系泊系統(tǒng)之間的相互作用;ZHU Feng等[10]為了解中長周期波浪對船舶動力響應的影響,基于勢流理論,對3種典型停泊集裝箱船在波浪相互作用下進行了水動力分析;劉明維等[11]基于材料力學的基本理論,推導得到了船舶系纜力與浮式系船柱柱身應變之間的定量關系,并結(jié)合系纜力作用下浮式系船柱有限元模型,計算得出了浮式系船柱所受系纜力,用來評估過閘船舶的系纜安全;張婧等[12]以單點系泊FPSO為研究對象,從6自由度運動和系纜張力響應出發(fā),比較了張緊式和懸鏈線式系泊系統(tǒng)對于平臺水動力性能影響,研究結(jié)果表明選取合適的內(nèi)轉(zhuǎn)塔位置和系泊纜與海底夾角布置方案,可有效地降低船體運動響應及系泊纜張力。

筆者基于三峽庫區(qū)水域的實測資料,考慮三峽庫區(qū)淺水深影響,探討了船舶的系泊方式,分析了船舶水動力性能,評估了風、浪、流聯(lián)合作用下船舶的系泊安全。

1 計算理論

1.1 三維勢流理論

假設流體為無黏性、無旋、不可壓縮,則可引入速度勢φ(x,y,z,t)來描述流域內(nèi)的運動[13]。當海洋結(jié)構(gòu)物以自由面為基準時,速度勢滿足Laplace方程,如式(1):

?2φ(x,y,z,t)=0

(1)

Laplace方程和描述物體運動的速度勢需要進行線性化處理。假定波浪和結(jié)構(gòu)物的運動都較小,而流場中的速度勢由入射波速度勢、繞射勢和輻射勢疊加而成,如式(2)。

φ(x,y,z,t)=φI(x,y,z,t)+φD(x,y,z,t)+

φR(x,y,z,t)

(2)

式中:φI為入射波速度勢,表明流場中速度分布的情況;φD為繞射勢,表明結(jié)構(gòu)物對流場內(nèi)的速度產(chǎn)生的影響;φR為輻射勢,表明結(jié)構(gòu)物6個自由度的運動及振蕩對流場的影響。

1.2 系纜力計算方法

系泊纜在外力作用下發(fā)生形變,所產(chǎn)生的力為非線性[14],其計算如式(3):

(3)

式中:FR為纜繩拉力,kN;d為纜繩直徑,m;Kc為纜繩彈性常數(shù),對尼龍繩而言,Kc=1.56×104MPa,對鋼絲纜,Kc=2.75×105MPa;n為材料特性相關的指數(shù),尼龍繩n=1.5;ΔS/S為纜繩相對伸長比。

如果帶纜樁的坐標為(X1,Y1,Z1),船舶在初始船位時導纜孔的坐標為(X2,Y2,Z2),則纜繩應變[14]計算如式(4):

(4)

式中:Δ為導纜孔到系泊絞車的距離;Lo為船舶在初始位置時的纜繩原長。

2 模型建立

2.1 模型參數(shù)

建立模型時,預先對坐標系進行定義:船長為X軸方向,X正方向為船尾指向船艏;船寬為Y軸方向,Y正方向為右舷指向左舷;Z軸向上為正,水線面處Z=0,坐標原點為船尾水線面左右對稱點。環(huán)境載荷作用的方向定義是與X軸正向的夾角,逆時針旋轉(zhuǎn)為正。船舶六自由度運動響應中,沿X軸方向運動為縱蕩(surge);沿Y軸運動為橫蕩(sway);沿Z軸運動為垂蕩(heave);繞X軸轉(zhuǎn)動為橫搖(roll);繞Y軸轉(zhuǎn)動為縱搖(pitch);繞Z軸轉(zhuǎn)動為艏搖(yaw)。表1為墊擋船與碼頭的主尺度參數(shù)。

表1 模型主尺度參數(shù)

2.2 模型驗證

在進行水動力分析之前,為保證文中所建立的水動力模型數(shù)值模擬結(jié)果準確性,筆者對碼頭-船舶模型劃分不同密度的網(wǎng)格,探究在網(wǎng)格密度變化情況下,船舶波浪響應變化趨勢是否滿足要求。方案1為網(wǎng)格允許公差0.1 m,最小網(wǎng)格尺寸為1 m的船舶有限元模型;方案2為網(wǎng)格允許公差0.15 m,最小網(wǎng)格尺寸為1.2 m的船舶有限元模型;方案3為網(wǎng)格允許公差0.2 m,最小網(wǎng)格尺寸為1.5 m的船舶有限元模型。在這3種不同網(wǎng)格劃分下,選取相同的環(huán)境參數(shù),分別對其進行水動力性能數(shù)值分析。

圖1為波浪入射方向為45°時,兩種網(wǎng)格劃分下的縱蕩波浪幅值響應算子。

圖1 網(wǎng)格密度驗證

由圖1可知:這3種網(wǎng)格劃分下的縱蕩波浪幅值響應算子變化趨勢以及具體數(shù)值幾乎一致,其峰值相差僅為1%,故可用原來的水動力模型進行相關的水動力性能分析。

筆者最終采用ANSYS軟件建立船舶-碼頭的濕表面模型。網(wǎng)格允許公差0.1 m,最小網(wǎng)格尺寸為1 m;主船體劃分907網(wǎng)格,碼頭劃分1 766網(wǎng)格,共計2 673網(wǎng)格。碼頭-船舶模型如圖2。

圖2 碼頭和船舶模型示意

3 幅值響應計算分析

基于勢流理論,筆者分析了浪向角對幅值響應算子(response amplitude operator, RAO)影響。規(guī)定以船尾來向為0°浪向角,船首來向為180°浪向角,45°為間隔,共設置5個波浪方向,運用AQWA-LINE模塊進行計算,計算時長為180 s,以5～30 s為典型時間段進行分析。

圖3為浪向角對船舶各自由度幅值響應算子影響的時間歷程。

圖3 浪向角對幅值響應算子的影響

由圖3可知:浪向角對縱蕩幅值響應算子的影響較大,呈現(xiàn)出逐漸上升的趨勢;其中:0°和180°方向的浪向角對縱蕩幅值響應算子的影響最明顯,90°的浪向角對橫蕩、垂蕩和橫搖響應較大,對縱蕩、縱搖和艏搖影響較小。各方向的浪向角對船舶橫搖幅值響應算子的影響趨勢相似,呈現(xiàn)出先增大,后減小。

4 系泊系統(tǒng)設計及水動力性能分析

4.1 環(huán)境參數(shù)

筆者以三峽庫區(qū)典型環(huán)境條件作為模擬工況依據(jù),具體參數(shù)設置如表2。

表2 環(huán)境參數(shù)

4.2 系泊系統(tǒng)設計

將頻域計算的RAO結(jié)果導入AQWA-DRIFT模塊,采用橫向和斜向兩個方向的系泊纜索限制船舶在橫向和縱向的移動,計算得到2-7-2(下文簡稱11系泊)和4-4-4(下文簡稱12系泊)兩種系泊方式下的時間歷程曲線和系泊纜拉力變化規(guī)律。按照3 h回歸周期得到船舶各自由度響應的最大統(tǒng)計值,每組工況模擬3次,將統(tǒng)計值取平均,以60～160 s為典型時間段進行分析,系泊方式如圖4、圖5。

圖4 11系泊方式示意

圖5 12系泊方式示意

4.3 水動力性能分析

筆者對船舶系泊系統(tǒng)進行設計,開展了3 h不規(guī)則波作用下的數(shù)值計算。圖6為兩種系泊方式下船舶六個自由度運動幅值的時間歷程曲線。

由圖6可知:由于外環(huán)境的不規(guī)則性,船舶六自由度呈現(xiàn)出不規(guī)則的變化趨勢,通過兩種系泊方式對比發(fā)現(xiàn),11系泊方式下的垂蕩平均距離為0.505 m,且在運動幅值在0.3～0.6 m內(nèi)變化,振蕩效果明顯;12系泊方式下的垂蕩平均距離為0.772 m,在60～95 s內(nèi)運動頻率較大,振蕩范圍在0.65～0.95 m。11系泊方式下的縱蕩平均距離為0.563 m;12系泊方式下縱蕩平均距離為0.579 m,運動頻率較小,趨于穩(wěn)定。12系泊方式下橫蕩距離隨時間的增長而增加,11系泊方式下橫蕩平均距離為0.047 m且浮動范圍較小。研究表明:11系泊方式下船舶的橫蕩和垂蕩振幅較小,艏搖角度較小,可起到較好的控制效果。

圖6 船舶六自由度運動時歷曲線

在設置系泊纜參數(shù)時,依據(jù)尼龍纜繩材料屬性,纜繩長度為5 m,直徑為50 mm,破斷載荷1 400 kN。表3為兩種系泊方式下各系纜繩最大拉力和平均拉力統(tǒng)計。由表3可知:各自系泊方式下船舶系纜繩最大受力均發(fā)生在船艏和船艉處。11系泊方式下最大拉力出現(xiàn)在Line 2纜繩,最大值為249.3 kN;12系泊方式下最大拉力出現(xiàn)在Line 2纜繩,最大值為244.7 kN。根據(jù)石油公司國際海事論壇(OCIMF)“Mooring Equipment Guidelines(2008)”規(guī)定:“對于材料為尼龍的系泊纜繩,其所受拉力不應大于其最小破斷力的45%”,即該系纜繩拉力不宜超過630 kN,11、12兩種系泊方式下系纜繩的拉力低于規(guī)范中破斷載荷要求,滿足規(guī)范要求。

表3 各系纜繩最大拉力和平均拉力

5 結(jié) 論

筆者以三峽庫區(qū)船舶為研究對象,基于頻域勢流理論,采用ANSYS-AQWA軟件,對船舶及其系泊系統(tǒng)的水動力性能進行分析,通過對船舶在不同浪向角和不同系泊方式下的運動特性進行分析,主要結(jié)論如下:

1)0°和180°方向的浪向角對縱蕩幅值響應算子的影響最明顯,90°的浪向角對橫蕩、垂蕩和橫搖響應較大,對縱蕩、縱搖和艏搖影響較小;

2)各個方向的浪向角對船舶橫搖幅值響應算子的影響趨勢相似,呈現(xiàn)出先增大,后減小;

3)11系泊方式對船舶的橫蕩、垂蕩、艏搖三個自由度的運動起到較好的控制效果;

4)在系泊纜拉力計算中,2-7-2系泊方式下Line 2的拉力大于其他位置的纜繩,最大拉力達到249 kN,平均拉力達到225 kN;4-4-4系泊方式下纜繩拉力最大出現(xiàn)在Line 2,最大值為245 kN,平均拉力達到215 kN,滿足規(guī)范的相應要求。