潘德位,林成新,孫德平,劉志杰,周超玉

(1.大連海事大學(xué) 交通運(yùn)輸裝備與海洋工程學(xué)院,遼寧 大連 116026;2.大連海事大學(xué) 輪機(jī)工程學(xué)院,遼寧 大連 116026)

基于GHS軟件的傾覆船舶扳正計(jì)算與分析

潘德位1,林成新1,孫德平2,劉志杰1,周超玉1

(1.大連海事大學(xué) 交通運(yùn)輸裝備與海洋工程學(xué)院,遼寧 大連 116026;2.大連海事大學(xué) 輪機(jī)工程學(xué)院,遼寧 大連 116026)

在考慮難船浮性和穩(wěn)性的基礎(chǔ)上,研究了傾覆船舶的扳正過程。根據(jù)船舶的受力特點(diǎn),建立了打撈難船的扳正力力學(xué)模型;針對(duì)傳統(tǒng)扳正力編程計(jì)算時(shí)間長(zhǎng)和前期數(shù)據(jù)輸入量大的問題,提出利用GHS軟件的Salvage模塊進(jìn)行分析計(jì)算。并以某傾覆船舶扳正為例,建立三維模型并對(duì)其扳正過程進(jìn)行仿真計(jì)算研究,從計(jì)算結(jié)果可以分析出最大扳正力值及其出現(xiàn)的扳正階段、該船舶的穩(wěn)性狀態(tài)及船舶縱向上的強(qiáng)度分布情況。GHS軟件能夠合理模擬橫傾船舶的扳正過程,適用于制定船舶救助打撈方案。

傾覆船舶;沉船打撈;扳正力;浮力和穩(wěn)性;GHS軟件

船舶沉沒后在水下的姿態(tài)各有不同,有的為正坐即甲板向上,更多的為橫傾,其橫傾角度為0°～180°[1]。打撈時(shí),當(dāng)沉船出現(xiàn)較大橫傾角度時(shí),起浮前應(yīng)對(duì)沉船進(jìn)行扳正[2]。船舶扳正就是施加在難船上的扳正力矩大于阻力矩使船舶向正方向翻轉(zhuǎn)的過程。傾覆于水面的船舶的浮力主要來自于船艙內(nèi)儲(chǔ)備的空氣墊,所以保持船底部的完整或可以快速修補(bǔ)是非常重要的。扳正過程也是扳正力、浮力和船舶重力共同作用的結(jié)果。對(duì)于具有一定儲(chǔ)備浮力和穩(wěn)性的難船,可以通過計(jì)算進(jìn)行調(diào)載和卸載,使船舶橫傾與縱傾恢復(fù)至適于航行的狀態(tài)[3-5]。

本文在研究?jī)A覆船舶扳正方法的基礎(chǔ)上,根據(jù)扳正過程中船舶受力的特點(diǎn),推導(dǎo)出扳正傾覆船舶的力學(xué)公式。針對(duì)常規(guī)方法計(jì)算扳正力復(fù)雜低效的缺點(diǎn),嘗試使用GHS軟件Salvage模塊模擬扳正過程,求解所需最大扳正力和扳正角度,同時(shí)得到扳正各階段沿船舶縱向上的扭矩分布值,對(duì)于合理統(tǒng)籌設(shè)計(jì)扳正施力點(diǎn)具有現(xiàn)實(shí)的指導(dǎo)意義。

1 GHS軟件簡(jiǎn)介

GHS軟件是以應(yīng)用水靜力學(xué)程序?yàn)榛A(chǔ)不斷發(fā)展起來的,它利用水靜力學(xué)理論,分析任意形狀海洋浮體結(jié)構(gòu)在給定吃水和傾角時(shí)所具備的浮性、穩(wěn)性和強(qiáng)度。目前,GHS軟件由4個(gè)模塊組成[6],第1個(gè)模塊為GHS,用于所有類型的船舶及其他浮動(dòng)結(jié)構(gòu)設(shè)計(jì)計(jì)算和評(píng)估的軟件系統(tǒng),它通過使用船只的數(shù)學(xué)/幾何模型計(jì)算其受力狀態(tài)來求解浮性,傾角,穩(wěn)性和強(qiáng)度;第2個(gè)模塊為GHS/Salvage,在打撈船舶前期計(jì)算時(shí),可以不限制縱傾和橫傾,模擬出船舶直立甚至倒立狀態(tài);第3個(gè)模塊為GLM,用于船舶在各種載荷作用下的穩(wěn)定性和強(qiáng)度監(jiān)測(cè),具備快速檢查船舶安全能力和預(yù)期負(fù)載能力;第4個(gè)模塊為BHS,是GHS的衍生程序,具備GHS很多功能,但更易于編輯模型,生成船艙和船舶附體等,特別是可以生成高質(zhì)量的圖表。本文分析計(jì)算時(shí)采用的是GHS/Salvage模塊。

2 數(shù)學(xué)模型

對(duì)于一定傾角的難船進(jìn)行打撈,首先要對(duì)難船進(jìn)行扳正。通常需要扳正的難船呈3種狀態(tài)[7]:①小傾角傾斜,需要施加很小的扳正力,甚至不用施加外力,調(diào)節(jié)艙室載荷即可實(shí)現(xiàn)難船扳正;②大傾角傾斜,傾角接近或超過90°,需要外力輔助扳正;③船舶傾覆,一般出現(xiàn)在相對(duì)船舶較深的水域。

圖1 難船扳正過程的某一階段受力示意Fig.1 Force diagram of righting process at a certain stage

為了便于分析,忽略縱傾的影響,只考慮橫傾的變化。因此,船舶的位置可以用橫傾角和吃水2個(gè)參變量來表示。難船扳正過程的某一階段受力如圖1。

當(dāng)難船自浮水面穩(wěn)定時(shí),根據(jù)重力與浮力平衡可得到:

式中:W為船舶自重;▽為船舶排水量。

假設(shè)在外力F作用時(shí),扳正的過程是分段緩慢進(jìn)行,各階段可近似地看作是穩(wěn)定的,故滿足以下的力平衡式:

式中:Φ為橫傾角;B為船舶型寬。

當(dāng)船舶發(fā)生瞬間傾覆后,貨艙內(nèi)仍充滿著空氣,這一部分空氣會(huì)給船舶提供浮力,稱為空氣墊?？諝鈮|體積可以假設(shè)是等溫變化的,根據(jù)波意耳定律[8]:

式中:h為傾斜之前船內(nèi)自由液面低于船外水面的距離;h0為傾斜之后的船內(nèi)自由液面低于船外水面的距離;V0為初始空氣體積;V為某一橫傾角下的空氣墊體積;P0為大氣壓強(qiáng)。

傾覆船舶空氣墊的存在增加了儲(chǔ)備浮力。但是,隨著扳正過程的進(jìn)行,空氣墊會(huì)逐漸縮小直到消失,同時(shí)難船艙室內(nèi)的自由液面對(duì)扳正過程的阻力也會(huì)下降?？偢×ψ冃?dǎo)致扳正可以利用有效浮力減小,因而需要增加提升沉船力量。否則難船可能會(huì)直接沉底,大大增加了扳正難度,如果有條件可以在難船上施加吊力阻止下沉。如果空氣墊的面積很大但是總的體積卻很小,壓縮空氣的影響可以忽略。對(duì)于大型船舶和空氣墊很大的船舶而言,不能忽略空氣壓縮的影響。

根據(jù)浮力與重力的合力矩平衡可得:

進(jìn)一步根據(jù)式(2)和(6)可得關(guān)于Φ,T的扳正力公式:

通常扳正力需要多次根據(jù)橫傾角和吃水計(jì)算得到。對(duì)于任何一個(gè)給定的橫傾角,可以求出式(7)的解,進(jìn)一步輸入吃水得出扳正力,可以制成圖表以便分析。但是,這樣只能得到某一船只的扳正力,不具備普遍應(yīng)用的條件,打撈計(jì)算起來極為不便,難以適應(yīng)快速打撈的需要。

如果針對(duì)扳正力公式,對(duì)于在任何給定吃水和橫傾狀況下船舶所有靜水力特征的輸出信息編輯程序求解,會(huì)使計(jì)算得到極大的簡(jiǎn)化,從而加快工作速度,但也需要一定的時(shí)間。本文根據(jù)船舶的具體信息,采用GHS軟件模擬其扳正過程,可以迅速求解。

3 模擬計(jì)算實(shí)例

在傳統(tǒng)的打撈工程中,扳正力的計(jì)算采用表格法和經(jīng)驗(yàn)公式或半經(jīng)驗(yàn)公式,所得到的結(jié)果往往與實(shí)際結(jié)果偏差較大,曾出現(xiàn)過因扳正力計(jì)算錯(cuò)誤導(dǎo)致浮吊損毀的情況。所以,傳統(tǒng)計(jì)算的扳正力對(duì)打撈意義有限。扳正現(xiàn)場(chǎng)只能不斷調(diào)整試探,從而安排所需設(shè)備,大大增加打撈成本和時(shí)間,工程效率極低[9]。利用GHS軟件輔助設(shè)計(jì)沉船扳正方案能很好地解決上述問題,合理配置扳正力(PULL命令),使扳正方案更為科學(xué),從而完善打撈方案。

3.1 船體模型建立

本文利用GHS軟件建立某船舶的模型(圖2)。建模完成之后經(jīng)過進(jìn)一步修正得到比較準(zhǔn)確的模型,以便在隨后的計(jì)算中得到最接近實(shí)際狀態(tài)的模擬數(shù)據(jù)。圖2中艙室分別填充紫色、深藍(lán)、淡藍(lán)、褐色、黃色和綠色。假設(shè)該船由于風(fēng)浪作用,以傾覆狀態(tài)漂浮在海面。為了方便扳正難船,卸掉艙室內(nèi)的貨物和燃料。由于褐色艙室頂蓋存在破損,導(dǎo)致該艙進(jìn)水,紅色部分表示進(jìn)水口。難船的傾覆狀態(tài)如圖3,圖中由2個(gè)橫向剖面圖組成,文中a表示船尾方向,f表示船首方向,s表示船舶右側(cè),p表示船舶左側(cè)。

圖2 船舶船體和艙室Fig.2 The hull and cabin of a yacht

圖3 難船初始狀態(tài),左傾180°(上移)Fig.3 The initial state,180 degrees to the left

3.2 靜穩(wěn)性

在圖4中,A,B曲線分別為該船遇難前、后的靜穩(wěn)性曲線。船體遇難后,由于上層艙室進(jìn)水,嚴(yán)重降低了船體的穩(wěn)性,甚至出現(xiàn)負(fù)初穩(wěn)性,船體的最大復(fù)原力臂、極限靜傾角和穩(wěn)矩都不同程度減小,船體的抗傾斜能力必然降低。

3.3 初始狀態(tài)計(jì)算

對(duì)建立的模型進(jìn)行初始狀態(tài)仿真模擬,可以求解出船舶的吃水、縱傾角、橫傾角、重量和重心,從而大大減少打撈前期人工探摸的工作量,提高工作效率。

本次計(jì)算中,船舶總長(zhǎng)L=5.5 m,型寬B=2.18 m,型深D=0.9 m,船體質(zhì)量W=1.1 t,初始狀態(tài)時(shí)橫傾角=180°(p),縱傾角=1.59°(a)。在模型的橫剖面上合理布置施力作用點(diǎn),為簡(jiǎn)化分析,在船體的兩個(gè)橫剖面上合理布置6個(gè)施力作用點(diǎn),確定每個(gè)施力點(diǎn)扳正力的作用方向和角度,見表1。

圖4 船舶靜穩(wěn)性曲線Fig.4 Curves of statical stability of a yacht

表1 初始載荷下某船舶的基本數(shù)據(jù)Tab.1 Basic data of a yacht in initial state

3.4 扳正過程模擬

3.4.1扳正過程 迅速準(zhǔn)確計(jì)算難船所需的最大扳正力能有效減少扳正時(shí)間,對(duì)整個(gè)打撈工程非常必要。通過GHS軟件模擬仿真扳正過程,能很快解出結(jié)果,為打撈初期方案設(shè)計(jì)提供理論參考。打撈工作開始后,適當(dāng)?shù)卦诟髯饔命c(diǎn)施加扳正力,當(dāng)扳正力矩大于阻力矩,難船開始滾動(dòng)。扳正力達(dá)到所需要的最大值并不是一次實(shí)現(xiàn)的,而是通過數(shù)次加載完成,扳正過程也是施加扳正力的過程。為了保證扳正過程的安全、有序,每進(jìn)行一次扳正力的加載,都應(yīng)當(dāng)暫停,通過對(duì)扳正纜和沉船等進(jìn)行觀察,確定無安全隱患后,再繼續(xù)按扳正方案進(jìn)行加載。

由初始狀態(tài)開始對(duì)扳正力適當(dāng)賦值,并逐漸遞增(本例中增力為0.05 t)。難船扳正過程見圖5和6,圖中分別為船體縱向坐標(biāo)為1(f)和1.25(a)的截面。圖6中船舶向右傾斜表示完全扳正后施加的扳正力產(chǎn)生扳正力矩繼續(xù)對(duì)船舶作用的結(jié)果。

圖5 難船扳正過程中任一狀態(tài),左傾167.78°Fig.5 A certain state of righting process,167.78 degrees to the left

圖6 難船完全扳正后狀態(tài),右傾30.93°Fig.6 The state after sitting,30.93 degrees to the right

3.4.2破艙計(jì)算 船舶傾覆后,破損艙室進(jìn)入一定量的海水。隨著船舶的轉(zhuǎn)動(dòng),在船體自重、傾角和艙室內(nèi)剩余氣體的作用下,艙室內(nèi)的進(jìn)水量和空氣占艙室的體積分?jǐn)?shù)會(huì)發(fā)生變化,見表2。

表2 進(jìn)水艙室基本數(shù)據(jù)Tab.2 Basic data of damaged cabin

隨著船體橫傾角度的變化,艙室進(jìn)水口的深度發(fā)生變化,在壓力作用下艙室內(nèi)空氣體積分?jǐn)?shù)也相應(yīng)發(fā)生改變。由于該船自重較小,完整艙室均無載貨,吃水較淺,導(dǎo)致破艙內(nèi)壓力變化較小,故軟件沒有精確顯示出破艙進(jìn)水量小范圍的波動(dòng)。

圖7 扳正過程中扳正力的變化Fig.7 Variation of righting force during uprighting process

3.4.3扳正力 船體扳正前期,扳正力值會(huì)隨著船體橫傾角度的減小而逐漸增加。當(dāng)扳正力值達(dá)到最大后,會(huì)隨著船體轉(zhuǎn)動(dòng)而逐漸降低,直至外力扳正作用消失。實(shí)際打撈工程中,根據(jù)最大扳正力,可以進(jìn)一步預(yù)備相關(guān)扳正設(shè)備和制定詳細(xì)的打撈方案。圖7為所需扳正力值隨橫傾角變化。

圖7中,扳正力變化線由兩組線段組成,其中實(shí)線段表示船體扳正過程中扳正力逐漸增大的過程;虛線表示船體正坐后持續(xù)施加最大扳正力使其傾斜的結(jié)果。虛線段一端較大值表示整個(gè)扳正過程的最大扳正力值為1 t。扳正力的大小與施力點(diǎn)的位置設(shè)計(jì)、扳正力方向和角度相關(guān),有條件的情況可以利用設(shè)備增加扳正力矩,降低扳正力。扳正點(diǎn)的位置應(yīng)根據(jù)實(shí)際情況進(jìn)行調(diào)整,保證扳正過程中船體的安全。

3.4.4船體強(qiáng)度 船舶扳正過程中,由于外力作用,需要慎重考慮船體強(qiáng)度的分布情況。正常航行中,普通船舶極少考慮船體縱向扭轉(zhuǎn)強(qiáng)度,通常在滿足縱向彎曲強(qiáng)度的條件下,扭轉(zhuǎn)強(qiáng)度一般也能滿足。但是,在船舶扳正過程中,則需要適當(dāng)考慮扭轉(zhuǎn)強(qiáng)度。由圖7可知,船體扳正到左傾167.08°后,扳正力達(dá)到最大值。圖8為難船在左傾167.08°時(shí)船體強(qiáng)度分布曲線。

圖8 船體強(qiáng)度分布(船尾←縱向→船頭)Fig.8 Intensity distribution of hull(stern←longitudinal→bow)

船體質(zhì)量分布是難船打撈設(shè)計(jì)的基礎(chǔ)參數(shù)之一,包括船體質(zhì)量和載貨分布。浮力分布與船體質(zhì)量分布、船體型線、船體傾角和吃水密切相關(guān)。剪力分布為船體質(zhì)量分布、浮力分布和扳正力分布的代數(shù)和的積分。而彎矩分布和扭矩分布則為剪力分布的進(jìn)一步積分。船體強(qiáng)度的計(jì)算結(jié)果是整體打撈難船的依據(jù),也是檢驗(yàn)扳正點(diǎn)位置設(shè)計(jì)合理與否的重要依據(jù)。

難船扳正過程中,GHS軟件只能求解出某一階段的扭矩分布值,不能得到整個(gè)扳正過程中的扭矩值。由于GHS軟件能快速計(jì)算船舶強(qiáng)度,特殊工況下,可以減小扳正力的增值,模擬扳正過程中的多個(gè)階段求解扭矩,使各階段之間的差值相對(duì)較小,得到的數(shù)據(jù)盡可能詳細(xì),所以對(duì)于打撈前期的力學(xué)估算依舊具有重要的實(shí)際意義。

4 結(jié) 語

在考慮難船浮性和穩(wěn)性的基礎(chǔ)上,根據(jù)打撈方案,對(duì)難船扳正過程進(jìn)行力學(xué)分析,并推導(dǎo)出扳正過程的扳正力公式。針對(duì)船體打撈計(jì)算復(fù)雜、低效的缺點(diǎn),利用GHS軟件對(duì)扳正過程進(jìn)行仿真模擬,得到扳正該船舶所需的最大扳正力,特別是船舶縱向上的強(qiáng)度分布值,可以檢驗(yàn)初始設(shè)定的扳正施力點(diǎn)合理性。所采用的軟件分析方法,遵循力學(xué)基本原理,對(duì)其他相似情況的船舶扳正工作均適用。

GHS軟件的計(jì)算方式是線性修正計(jì)算,沿船長(zhǎng)方向切分的微體積單元通過線性連接以差分法計(jì)算。因此,該軟件計(jì)算簡(jiǎn)單,所占計(jì)算機(jī)資源較少,運(yùn)行速度更快。GHS軟件能實(shí)現(xiàn)逆向邏輯求解船舶穩(wěn)性狀態(tài),即給定某一船舶狀態(tài),逆向反推回去計(jì)算滿足此狀態(tài)下壓載艙、貨艙以及其他一些艙室的配載情況,這對(duì)迅速響應(yīng)救助打撈工程也非常適用。

參 考 文 獻(xiàn):

[1]張?jiān)雒?打撈工程[M].大連:大連海事大學(xué)出版社,2012:208-211.(ZHANG Zeng-meng.Salvage engineering[M]. Dalian:Dalian Maritime University Press,2012:208-211.(in Chinese))

[2]劉漢明,臧海鵬.沉船打撈技術(shù)設(shè)計(jì)與分析[M].青島:中國海洋大學(xué)出版社,2011:18-21.(LIU Han-ming,ZANG Haipeng.Technical design and analysis on salvage[M].Qingdao:China Ocean University Press,2011:18-21.(in Chinese))

[3]PHILLIP S.Counter-capsizing a battleship[J].Naval Engineers Journal,2005,117(1):101-110.

[4]MORRIS L P.Salvage of the Oklahoma at Pearl Harbor[J].Engineering and Science,1947(11):9-17.

[5]United States Navy.U S Navy salvage manual volume 2:pol offloading[M].American:Direction of Commander,Naval Sea Systems Command,1991:22-24.

[6]Creative Systems.GHS user′s manual[M].American:Creative Systems,Inc.2007.

[7]ZHANG Meng-xin.Discussion on uprighting big-angle tilted sunken ships underwater in salvage operation[J].China Ocean Engineering,1989(2):217-227.

[8]DROBYSHEVSKI Y.A note on uprighting of a ship floating upside-down[J].Ocean Engineering,2004,31(11/12):1447-1467.

[9]林金雄,王仁義,孫召才.“銀河山168”沉船打撈工程[C]∥救撈專業(yè)委員會(huì)2006年學(xué)術(shù)交流會(huì)論文集,2006:83-89. (LIN Jin-xiong,WANG Ren-yi,SUN Zhao-cai.The salvage of the ship Yinheshan l68[C]∥Proceedings of China International Rescue&Salvage Conference 2006(CIRSC),2006:83-89.(in Chinese))

Calculation and analysis of righting capsized vessel based on GHS software

PAN De-wei1,LIN Cheng-xin1,SUN De-ping2,LIU Zhi-jie1,ZHOU Chao-yu1
(1.Transportation Equipment&Ocean Engineering,Dalian Maritime University,Dalian 116026,China; 2.Marine Engineering College,Dalian Maritime University,Dalian 116026,China)

The uprighting process of the capsized vessel considering the buoyancy and stability of the wreck has been studied in this paper.A model of the righting force for the wreck salvage is established according to the mechanical characteristics of the wreck.The analysis and calculation are demonstrated by the salvage module of GHS software in view of the problem of righting force programme,which needs long computing time and large amount of data input.Taking the simulation of righting a capsized ship as an example,a three-dimensional model is developed and the righting process is simulated by using GHS software.The maximum righting force and the occurring phase,angle of vanishing stability and longitudinal distribution of torque values can be obtained from the analysis.And the analysis results can play a crucial role in the uprighting process of the capsized ship.The uprighting process of the inclining ship can be simulated by using GHS software,which is suitable for making salvage schemes.

capsized vessels;wreck salvage;righting force;buoyant force and stability;GHS software

U676.6

A

1009-640X(2014)06-0078-06

2014-04-29

潘德位(1986-),男,遼寧莊河人,博士研究生,主要從事船舶救助打撈的研究。E-mail:wssrwsmt@163.com

林成新(E-mail:lchxin@dlmu.edu.cn)