蔡計強， 諶志新

(1 上海交通大學(xué)船舶海洋與建筑工程學(xué)院，上海 200240；2 農(nóng)業(yè)部漁業(yè)裝備與工程技術(shù)重點實驗室，中國水產(chǎn)科學(xué)研究院漁業(yè)機械儀器研究所，上海 200092)

?

基于NAPA軟件的漁業(yè)船舶完整穩(wěn)性分析

蔡計強1,2， 諶志新2

(1 上海交通大學(xué)船舶海洋與建筑工程學(xué)院，上海 200240；2 農(nóng)業(yè)部漁業(yè)裝備與工程技術(shù)重點實驗室，中國水產(chǎn)科學(xué)研究院漁業(yè)機械儀器研究所，上海 200092)

為驗證NAPA軟件在漁業(yè)船舶應(yīng)用上的可行性和采用船舶靜力學(xué)原理計算完整穩(wěn)性的精確性，根據(jù)南海海域49.5 m金槍魚延繩釣船的特定作業(yè)工況，分別采用NAPA軟件和表格進行計算。計算結(jié)果表明，兩種計算方法的偏差在2%左右，能夠滿足工程應(yīng)用的精度要求。對氣象衡準(zhǔn)計算，分別采用了漁業(yè)船舶規(guī)則和國際海事組織(IMO)規(guī)則，由于算法不同，其偏差較大。通過比較發(fā)現(xiàn)，漁業(yè)船舶規(guī)則對氣象衡準(zhǔn)的要求高于IMO規(guī)則的要求。在沒有條件使用大型船舶計算軟件，且船舶處于平浮狀態(tài)時，采用船舶靜力學(xué)原理的計算結(jié)果是精確的。采用了船舶靜力學(xué)原理并建立了NAPA計算模型，通過實例計算表明，NAPA軟件具有在漁業(yè)船舶應(yīng)用上的可行性和采用船舶靜力學(xué)原理計算完整穩(wěn)性的精確性，該結(jié)論可為船舶設(shè)計人員選取不同計算方法時提供參考依據(jù)。

完整穩(wěn)性；NAPA；漁業(yè)船舶；氣象衡準(zhǔn)

完整穩(wěn)性是關(guān)系船舶能否安全運營的一個重要指標(biāo)[1-4]。船級社和漁船檢驗機構(gòu)都非常重視對其的檢驗和審查，并頒布了詳細的檢驗衡準(zhǔn)[5-6]。鑒于船舶穩(wěn)性的重要性及近年來海損事故頻發(fā)的現(xiàn)實，國際海事組織(IMO)正在制定完整穩(wěn)性的第二代衡準(zhǔn)，并將其作為2008年國際完整穩(wěn)性規(guī)則(2008 IS CODE)A部分強制性要求的補充和替代方法。國內(nèi)已開展了相關(guān)的研究工作[7-10]。在圖紙審查方面，相較于運輸船舶由船級社審核，我國的漁業(yè)船舶則由漁船檢驗局審核。不同的檢驗機構(gòu)對完整穩(wěn)性衡準(zhǔn)的定義有所不同。運輸船舶一般采用芬蘭的總體性能軟件NAPA進行校核。而漁業(yè)船舶，由于受制于資金不足等問題，長期以來都是自行編制計算表格用于計算。

為提高漁業(yè)船舶的設(shè)計效率，已有多家設(shè)計單位引進了NAPA軟件，并運用于完整穩(wěn)性校核。運輸船舶長期采用NAPA軟件，其計算結(jié)果經(jīng)過設(shè)計單位、船廠和船級社的多次驗證，其建模方法和計算結(jié)果均具有較高的可信度。工程應(yīng)用一般都有一定的可接受誤差范圍，如中國船級社對初穩(wěn)心高GM值的可接受誤差為1%，對復(fù)原力臂GZ的可接受誤差為5%[11]。首次在漁業(yè)船舶上應(yīng)用NAPA軟件，并且由于船型、作業(yè)工況和校核衡準(zhǔn)的差異性，有必要對NAPA軟件進行校核。本文分別采用表格計算和NAPA軟件進行計算。表格計算是指運用船舶原理并采用辛浦生積分法所進行的穩(wěn)性計算(在不產(chǎn)生誤解的前提下，以下簡稱“表格計算”)；而NAPA軟件則采用數(shù)值積分的方法進行計算。

1 漁業(yè)船舶作業(yè)工況配載

1.1 漁業(yè)船舶作業(yè)工況

漁業(yè)船舶有多種作業(yè)工況，不僅要從事捕魚作業(yè)，并且要盡快地將漁獲物運輸上岸以保證魚的品質(zhì)，同時還要具備優(yōu)良的快速性以躲避惡劣海況，實現(xiàn)安全作業(yè)。漁業(yè)船舶規(guī)則要求必須校核滿載出港、捕魚中、滿載返港、滿載到港和空載到港等典型作業(yè)工況?？傮w布置時，漁業(yè)船舶通常將燃油艙布置在船舶的底部，在滿載出港時有利于降低船舶重心。但在滿載到港時，燃料消耗得只剩下10%，且漁獲物有時存放在高處的凍結(jié)間，這將提高船舶的重心。本文計算時將選取較危險的滿載到港作為校核工況。

1.2 空船重量估算與作業(yè)工況配載

本文的考察對象為在南海作業(yè)的49.5 m金槍魚延繩釣船，其垂線間長LPP=43 m，型寬B=8.2 m，型深D=3.55 m，設(shè)計吃水d=3.1 m，排水量Δ=780.78 t?？沾亓康臏?zhǔn)確估算對完整穩(wěn)性計算很重要。對相同船型進行統(tǒng)計回歸后，得到了漁業(yè)船舶重量重心的多項回歸公式[12]：

LW=K×Lpp0.9×B1.2×D0.7

(1)

式中：LW—船舶空船重量，t；K—無因次系數(shù)；Lpp—垂線間長，m；B—型寬，m；D—型深，m。

選取船型相近的母型船后運用該回歸公式并選取相同的K值，可以較為準(zhǔn)確地估算出新船的重量重心。估算得到的空船重量LW=527.5 t，船舶重心縱向坐標(biāo)LCG=19.25 m，船舶重心垂向坐標(biāo)VCG=3.33 m。表格計算時，需要計算空船重量、油水艙和其他配載的重量重心，并匯總垂向力矩和縱向力矩，最后得到該作業(yè)工況下的整船重量重心。用NAPA軟件計算時，在完成線型建模和艙室分艙后，只需要輸入空船重量重心并進行合理配載，就可以完成作業(yè)工況的配載。用到的NAPA語言為：LOAD，load，amount，location，如LOAD，F(xiàn)W，*1，R1.01，其含義為在艙室代號為R1.01的No.1 淡水艙內(nèi)加載滿艙的淡水。用以上2種方法分別進行配載，其計算結(jié)果是一致的，排水量766.6 t，重心距艉柱20.504 m，重心高度3.071m。

2 自由液面修正

2.1 自由液面對初穩(wěn)心高的修正

根據(jù)漁業(yè)船舶規(guī)則[6]的要求，凡存在自由液面的液體艙，均應(yīng)計算裝載50%艙容液體的自由液面影響，如果艙的形狀特殊，存在更為不利的自由液面影響，則應(yīng)按更不利的情況計算自由液面影響，而IS CODE[5]要求按照形成最大慣性矩的液面進行計算。為了船舶安全、提高穩(wěn)性裕度，船級社和漁船檢驗局一般都按形成最大慣性矩的液面進行初穩(wěn)心高修正。本文的兩種計算方法均考慮最大自由液面影響。計算時對表1中的液艙進行自由液面修正。表格計算時，自由液面修正力矩為26.66 t·m，自由液面修正值為0.034 8 m；應(yīng)用NAPA軟件時，自由液面修正力矩為26.68 t·m，自由液面修正值為0.034 8 m?？梢钥闯?，兩種算法的計算結(jié)果是一致的[13-14]，按照船舶靜力學(xué)原理，其初穩(wěn)心高修正公式如下：

δGM=∑ρnIn/Δ

(2)

式中：δGM—初穩(wěn)心高修正值，m；ρn—某一個液艙內(nèi)的液體密度，t/m3；In—某一個液艙自由液面慣性矩，m4；n—全船應(yīng)計及自由液面影響的液艙的序號。

GM=KM-KG-δGM

(3)

式中：KM—橫穩(wěn)心高，m；KG—作業(yè)工況下的重心高度，m。

通過公式(3)，可以算出該作業(yè)工況下的初穩(wěn)心高GM。

2.2 自由液面對穩(wěn)性曲線的修正

考慮自由液面對穩(wěn)性曲線的修正時，應(yīng)根據(jù)漁業(yè)船舶規(guī)則[6]的要求進行核算，計算某一個液艙自由液面修正力矩，公式如下：

(4)式中：Mφ—不同橫傾角的自由液面修正力矩，t·m；V—艙柜的最大容積，m3；b—艙柜的最大寬度，m；ρ—艙柜中液體的密度，t/m3；k—無因次系數(shù)(按照b/h進行計算)；δ—艙柜的方形系數(shù)；h—艙柜的最大高度，m；。

(5)

(6)

(7)

采用公式(3)～(7)，算出不同橫傾角的自由液面修正力矩Mφ，并用下面公式計算復(fù)原力臂減小值δGZ：

(8)

表1為按照上面的計算方法并用表格計算得到的自由液面修正力矩。表中所列艙室在船上的位置按圖1所示。

表1 自由液面修正力矩

R1.01P&S—No.1淡水艙(左&右)；R1.02P&C&S—No.2 淡水艙(左&中&右)；R2.01P&C&S—No.1 燃油艙(左&中&右)；R2.02P&S—No.2燃油艙(左&右)；R2.03—日用油柜；R3.01—滑油艙；R3.02—污油艙；R4.01、R4.02、R4.03—No.1、No.2、No.3冷藏魚艙；R5.01—生活污水艙。圖1 NAPA計算模型Fig.1 The NAPA calculation model

用NAPA軟件計算時，按照IS CODE的方法進行修正。NAPA軟件的語法為：FRS，‘(R1.02P R1.02S R1.02C)/ALL IMO 0’，‘(R2.01P R2.01S)/ALL IMO 0’，其含義為R1.02(P&C&S)和R2.01(P&S)按照IMO的規(guī)則進行修正。NAPA軟件沒有詳細列出自由液面修正過程，只能通過改變修正方法來更改。

采用公式(8)進行表格計算得到不同橫傾角度的復(fù)原力臂減小值δGZ，并與NAPA軟件計算得到的修正值進行比較(表2)，兩種計算方法的最大差異比值為3.17%。復(fù)原力臂減小值δGZ對靜穩(wěn)性曲線的準(zhǔn)確繪制具有至關(guān)重要的作用，差異越小的δGZ表明靜穩(wěn)性曲線的差異也越小。雖然NAPA軟件沒有給出自由液面的具體計算公式，但計算結(jié)果表明，船舶原理對表格計算或NAPA軟件均具有普適性。由于δGZ的數(shù)值為毫米級別，其計算結(jié)果的差異更多來自截斷誤差。

表2 不同橫傾角下的復(fù)原力臂減小值δGZ

3 氣象衡準(zhǔn)校核

3.1 漁業(yè)船舶規(guī)則對氣象衡準(zhǔn)的規(guī)定

氣象衡準(zhǔn)數(shù)K值體現(xiàn)了船舶在風(fēng)浪聯(lián)合作用下的穩(wěn)性水平，一般情況下K值越大穩(wěn)性越好[15-17]。采用漁業(yè)船舶規(guī)則[6]進行校核時，穩(wěn)性衡準(zhǔn)數(shù)K應(yīng)符合下式要求：

(9)

式中：lq—最小傾復(fù)力臂，m；lf—風(fēng)壓傾側(cè)力臂，m。

風(fēng)壓傾側(cè)力臂lf按下式計算：

(10)

式中：Av—船舶受風(fēng)面積，m2；Z—計算風(fēng)力作用力臂，m；Δ—所核算裝載工況下船舶的排水量，t；P—單位計算風(fēng)壓，Pa。該工況在遠海航區(qū)時，其風(fēng)力作用力臂為2.689 m，受風(fēng)面積為231.54 m2，對應(yīng)的單位計算風(fēng)壓P為1 062.3 Pa。表格計算時，由圖1可得最小傾覆力臂，該值大小與船舶橫搖角和進水角有關(guān)系。如圖1所示且經(jīng)計算，最小傾覆力臂lq=0.171 3 m，風(fēng)壓傾側(cè)力臂lf=0.088，K=1.948。

3.2 國際海事組織對氣象衡準(zhǔn)的規(guī)定

我國漁業(yè)船舶規(guī)則與國際海事組織(IMO)規(guī)則在氣象衡準(zhǔn)計算上所用的方法是不一樣的。前者采用力臂比，后者采用面積比，但其所要求的K值都要大于1。從表3可以直觀地看出兩種衡準(zhǔn)的區(qū)別[2]。

表3 氣象衡準(zhǔn)對照表

用NAPA進行氣象衡準(zhǔn)計算時，其定義語法如下：

CRIT, IMOWEATHER, ′IMO weather criterion′

TYPE, ARATIO1

REQ, 1

RANG, EQ-ROLL, MIN(50, FAUN)

MOM, IMOWIND

OK

該條語法的意思是：無限航區(qū)時采用IMO定義的風(fēng)壓，即504 Pa，其橫傾角范圍為從平衡點向左到橫搖角，向右到50度或進水點的較小值。經(jīng)計算，b=0.150 7 m·rad，a=0.048 8 m·rad，最終計算結(jié)果為面積比K=b/a=3.088。

4 完整穩(wěn)性衡準(zhǔn)

4.1 穩(wěn)性衡準(zhǔn)計算方法介紹

衡量漁業(yè)船舶穩(wěn)性的優(yōu)劣性，主要從以下幾條衡準(zhǔn)來判斷：不同橫傾角時GZ曲線下的面積；初穩(wěn)性高、橫傾30度時的復(fù)原力臂；最大復(fù)原力臂對應(yīng)角和氣象衡準(zhǔn)數(shù)。表格計算時，將靜穩(wěn)性臂和動穩(wěn)性臂畫在CAD圖上(圖2)，再從圖中量取所需的面積值、復(fù)原力臂和角度等。

用NAPA計算，則是采用編制NAPA BASIC語言實現(xiàn)穩(wěn)性衡準(zhǔn)的定義。

圖2 表格計算完整穩(wěn)性Fig.2 The table calculation of intact stability

4.2 完整穩(wěn)性差異比較

為了直觀地反映完整穩(wěn)性衡準(zhǔn)的差異性，對國內(nèi)規(guī)則和國際規(guī)則對完整穩(wěn)性的不同要求進行了匯總，并表列了采用漁業(yè)船舶規(guī)則的表格計算和采用2008國際完整穩(wěn)性規(guī)則(2008 IS CODE)

的NAPA計算的驗算結(jié)果，列出了這兩種計算方法的差異比值(表4)。

由表4可以看出，除了氣象衡準(zhǔn)差異較大外，其余衡準(zhǔn)數(shù)的偏差都在2%以內(nèi)。表格計算時，K=1.948，而NAPA計算時，K=3.088。其計算結(jié)果的差異更多地來源于計算公式的不一致，如表3所示的橫搖角計算公式。表格計算時橫搖角φa=17.974，而NAPA計算時φa=18.714。由表4也可看出，采用漁業(yè)船舶規(guī)則所得的K值要比采用IMO規(guī)則所得的數(shù)值小，即漁業(yè)船舶規(guī)則更加嚴格[2]，與通常認為的國際法規(guī)比國內(nèi)法規(guī)嚴格是有矛盾的。

漁業(yè)船舶規(guī)則[6]上寫明：如果船舶符合IMO的穩(wěn)性衡準(zhǔn)要求，驗船部門可予以同意。表4也對完整穩(wěn)性進行了對比，其不同點在于對初穩(wěn)心高GM的要求不同，對常規(guī)運輸船GM最小值為0.15 m，對漁業(yè)船舶則為0.35 m。通過本案例分析可以發(fā)現(xiàn)，采用IMO規(guī)則進行校核時，特別是氣象衡準(zhǔn)，考慮到船舶安全性，應(yīng)該留有一定的余量，這可為漁業(yè)船舶規(guī)則的修正提供參考依據(jù)。

表4 完整穩(wěn)性差異匯總表

5 結(jié)論

分別采用表格計算和NAPA軟件進行計算，可以起到相互驗證的作用。在驗證過程中，設(shè)計人員能更好地掌握軟件的應(yīng)用方法和其中參數(shù)的具體含義，避免由于參數(shù)設(shè)定不正確而造成計算結(jié)果的錯誤。在作業(yè)工況較少的情況下，采用表格計算是沒問題的，但在作業(yè)工況較多的情況下，會比較繁瑣，且容易出錯。如果第二代完整穩(wěn)性衡準(zhǔn)生效，將加大表格計算的難度，也將影響計算結(jié)果的精確性。

NAPA軟件目前已經(jīng)得到了廣泛的應(yīng)用，并且有部分船級社直接將其作為校核軟件，如日本船級社NK、意大利船級社RINA和英國LR船級社等。該軟件具有優(yōu)良的拓撲關(guān)系，這將提高多工況的計算效率，也能快速地比較多個設(shè)計方案的優(yōu)劣性，可縮短計算時間。采用NAPA軟件計算具有快速性和準(zhǔn)確性的特點，該軟件應(yīng)用于漁業(yè)船舶，能夠提高設(shè)計效率、減少出錯率，有利于漁業(yè)船舶工程快速發(fā)展。

□

[1] 李曉峰,周瑋. 淺析DR67對耙吸挖泥船完整穩(wěn)性的要求 [J]. 船舶, 2007 (2)：11-15.

[2] 劉金剛. 海洋救助船完整穩(wěn)性特點及校核方法 [J].船舶設(shè)計通訊, 2012 (2)：1-5.

[3] 李培勇,裘泳銘,顧敏童等. 多體船型完整穩(wěn)性計算 [J].上海交通大學(xué)學(xué)報, 2002,36 (11)：1560-1563.

[4] 石珣, 張高峰, 周耀華. 新型拖船完整穩(wěn)性衡準(zhǔn)應(yīng)對研究 [J].中國造船, 2015,56 (1)：230-239.

[5] 中華人民共和國海事局. 國際航行海船法定檢驗技術(shù)規(guī)則[M].北京：人民交通出版社, 2007：159-176.

[6] 漁業(yè)船舶檢驗局. 漁業(yè)船舶法定檢驗規(guī)則[M].北京：人民交通出版社, 2000：231-241.

[7] 顧民,魯江,王志榮. 國內(nèi)第二代完整穩(wěn)性衡準(zhǔn)研究主要成果和重點發(fā)展方向 [J].中國造船, 2015,56 (1)：1-8.

[8] 馬坤,劉飛,張如凱. 癱船穩(wěn)性第二層薄弱性衡準(zhǔn)技術(shù)研究 [J].中國造船, 2015,56 (1)：97-105.

[9] 付麗麗,鄭雙燕,柳衛(wèi)東,等. 第二代完整穩(wěn)性參數(shù)橫搖薄弱性衡準(zhǔn)發(fā)展綜述 [J].中國造船, 2015,56 (1)：217-223.

[10]李曉,朱永峨,周耀華.第二代完整穩(wěn)性衡準(zhǔn)的驗證與改進研究 [J].船舶與海洋工程, 2012 (4)：50-53.

[11]中國船級社.鋼質(zhì)海船入級與建造規(guī)范 [M].北京：人民交通出版社,2009：481-482.

[12]梁建生.玻璃鋼漁船總體設(shè)計中的重量控制[J].漁業(yè)現(xiàn)代化, 2000,37(1)：7-9.

[13]中國船舶工業(yè)集團公司.船舶設(shè)計手冊[M].北京：國防工業(yè)出版社, 2013：83-115.

[14]盛振邦.船舶靜力學(xué)[M].北京：國防工業(yè)出版社, 1992：73-106.

[15]齊紹江.關(guān)于IMO完整穩(wěn)性要求中氣象衡準(zhǔn)探討 [J].世界海運, 2008,31 (2)：30-32.

[16]張文新,包輝陽. NAPA軟件在鋪管船完整穩(wěn)性計算的應(yīng)用[J]. 船舶, 2008 (2)：62-63.

[17]馮振玉.國際國內(nèi)氣象衡準(zhǔn)的抗風(fēng)等級[J].中國船檢, 2015 (8)：27-30.

Analysis on the intact stability of the fishing vessel >based on NAPA software

CAI Jiqiang1,2， CHEN Zhixin1

(1 School of Naval Architecture, Ocean & Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China ; 2 Key Laboratory of Fishery Equipment and Engineering, Ministry of Agriculture, Fishery Machinery and Instrument Research Institute, Chinese Academy of Fishery Sciences, Shanghai 200092, China)

For the purpose of verifying the feasibility of NAPA software in fishing vessel and the accuracy of the stability calculation with the ship’s hydrostatical theory, according to the special loading condition of the 49.5 m tuna longliner ship, the NAPA software and the table calculation were used respectively. The calculation results showed that the deviation of the two methods was about 2%, which could satisfy the accuracy requirements of the engineering application. About the weather criterion calculation, the Fishing Vessel Statutory Inspection Rules and the International Maritime Organization (IMO) Codes were used respectively. Because of the method difference, the result deviation was a bit large. The comparison results indicated that the requirements of Fishing Vessel Statutory Inspection Rules to weather criterion was stricter than that of the IMO Codes. In the absence of using the large ship calculation software and the ship was untrimmed, the results using the ship’s hydrostatical theory were reliable. As the large ship calculation software is used frequently by the designer, so the calculation principle of the NAPA software is necessary for analysis. The ship’s hydrostatical theory was adopted and the NAPA calculation model was established. Through the real example calculation, the results showed the feasibility of NAPA software in fishing vessel and the accuracy of the stability calculation with the ship’s hydrostatical theory. In addition, the results could provide reference foundation when the ship designer selects different calculation methods.

intact stability; NAPA; fishing vessel; weather criterion

2015-09-03

2016-01-10

公益性行業(yè)(農(nóng)業(yè))科研專項“漁業(yè)節(jié)能關(guān)鍵技術(shù)研究與重大裝備開發(fā)(201003024)”；上海市科技興農(nóng)項目 “漁船船型標(biāo)準(zhǔn)化及節(jié)能技術(shù)集成示范與推廣( 滬農(nóng)科推字2012第2-3號)”

蔡計強(1985—)，男，碩士研究生，研究方向：船舶總體設(shè)計。E-mail：cjq-ssy@163.com

諶志新(1969—)，男，研究員，研究方向：海洋漁業(yè)工程裝備。E-mail：chenzhixin@fmiri.ac.cn

10.3969/j.issn.1007-9580.2016.01.008

U674.4

A

1007-9580(2016)01-041-06