孫元璋,陳洲,宋健,李東旭,朱凡述,楊卓懿

(1.招商局金陵船舶(威海)有限公司,山東 威海 264200;2.山東交通學(xué)院 船舶與港口工程學(xué)院,山東 威海 264200)

船舶裝載后,如果作用在船體各橫剖面上的彎矩代數(shù)和不等于零,則船體會(huì)發(fā)生總縱彎曲變形,也稱為拱垂變形,拱垂變形會(huì)對(duì)船舶的排水量產(chǎn)生明顯的影響。部分港口或碼頭,對(duì)船舶的最大吃水有限制,這種情況下,還應(yīng)關(guān)注中垂對(duì)載貨量的影響[1]。船體梁變形還會(huì)通過影響船舶浮力分布而對(duì)其產(chǎn)生一附加的彎矩、剪力。尤為明顯的是,當(dāng)船體梁置于一坦谷波中,波峰/波谷位于船中時(shí),波浪力的作用會(huì)進(jìn)一步放大船體梁的變形,此時(shí)這種“附加力”的作用將會(huì)更明顯[2]。不過從宏觀的角度分析,這種變形對(duì)減少總縱彎矩是有利的。因?yàn)樵陟o水、滿載、中垂的情況下,船舯浮力增大,兩端浮力減少,此時(shí)船中最大中垂彎矩會(huì)有略微減少;更進(jìn)一步,如果此時(shí)不是靜水狀態(tài),而是船舯處于坦谷波的波谷位置,由于坦谷波的“加強(qiáng)”作用,中垂變形越發(fā)增大,較之初始狀態(tài),船舯浮力增加更多,中垂彎矩減少作用更明顯。這也意味著,中垂變形可以降低中垂彎矩的峰值,船體梁的“抗彎”強(qiáng)度能力增強(qiáng)了;中拱的情形亦然。

1 拱垂變形對(duì)排水量的影響

1.1 拱垂作用下靜水力計(jì)算的理論

船舶靜水力參數(shù)的計(jì)算,以排水體積為例。在不考慮船體梁拱垂作用情況下,只要基于平浮水線對(duì)各站線剖面積做積分即可;考慮拱垂的情況下,首先須保證排水量、浮心位置一致(重力與浮力相平衡、重力矩與浮力矩相平衡),然后基于拱垂時(shí)水線對(duì)各站線剖面積進(jìn)行積分。

NAPA軟件在計(jì)算過程中,程序會(huì)進(jìn)行多次迭代,以求得相對(duì)精確解,保證排水量、浮心位置一致;迭代計(jì)算的基礎(chǔ)為拱垂變形后水線面,其既取決于彎矩、剪力分布以及船體梁整體的慣性矩信息,又會(huì)進(jìn)一步影響彎矩、剪力分布。

1.2 靜水力數(shù)據(jù)、拱垂變形量

對(duì)某船舶滿載離港(舯垂)、壓載到港(舯拱)兩種工況進(jìn)行分析[3],基于NAPA計(jì)算得到的拱垂變形量及靜水力相關(guān)數(shù)據(jù)見表1、2。

表2 壓載到港中拱變形

同時(shí),對(duì)拱垂變形數(shù)據(jù)進(jìn)行二階擬合,具體數(shù)據(jù)見圖1～2,其方差非常小。且與業(yè)內(nèi)認(rèn)為船體梁的拱垂變形沿船長(zhǎng)方向近似為拋物線的習(xí)慣看法是一致的。對(duì)于這么高的擬合度分析認(rèn)為,原因如下。

圖1 滿載離港中垂變形

圖2 壓載到港中拱變形

1)滿載離港、壓載到滿2個(gè)工況的配載基本上是均勻過渡的,不存在“隔艙”工況;空船重量沿船長(zhǎng)分布基本上是均勻的,除了艉部的樓子、橫艙壁會(huì)引起空船重量分布的波動(dòng);船殼常規(guī)的型線設(shè)計(jì),浮力沿船長(zhǎng)分布的變化基本上是均勻的。

2)船體梁的慣性矩,沿船長(zhǎng)分布的變化基本上是均勻的。

1.3 船舶運(yùn)行過程中拱垂變形對(duì)排水量的修正

對(duì)于大尺度船而言,船體梁的拱垂變形對(duì)排水量影響的絕對(duì)值是不可忽略的。目前市場(chǎng)上只有較少的裝載儀可對(duì)排水量進(jìn)行拱垂修正,業(yè)內(nèi)對(duì)排水量的修正仍多采用經(jīng)驗(yàn)公式。

1.3.1 經(jīng)驗(yàn)公式法

船舶營(yíng)運(yùn)過程中,船舶的排水量是基于艏、舯、艉水尺觀測(cè)值,修正為船體的艏、舯、艉型吃水,進(jìn)一步進(jìn)行中拱或中垂修正,查靜水力表得到。

平均吃水T的拱垂修正如下。

經(jīng)濟(jì)學(xué)理論顯示，政府監(jiān)管部門和市場(chǎng)主體的利益關(guān)切點(diǎn)不同，它們之間難免存在矛盾和博弈，市場(chǎng)主體進(jìn)入市場(chǎng)的目的是實(shí)現(xiàn)個(gè)體利益的最大化，而政府監(jiān)管部門則力爭(zhēng)通過維護(hù)市場(chǎng)秩序?qū)崿F(xiàn)公共利益的最大化，兩者看似對(duì)立，實(shí)則統(tǒng)一。在實(shí)踐層面，公共利益的過于抽象使得政府為謀取公共利益而制定的政策缺乏說服力，而實(shí)實(shí)在在的社會(huì)總體利益則扮演著公共利益代言人的角色。負(fù)面清單模式可以最大限度地保障市場(chǎng)主體追求個(gè)人利益，個(gè)人利益的最大化也會(huì)反哺社會(huì)，在一定程度上起到了引導(dǎo)社會(huì)總體進(jìn)步和發(fā)展的作用，最終與政府追求的社會(huì)總體利益最大化目標(biāo)相吻合。

T′=(TF+TA+6TM)/8

(1)

式中:TF、TA、TM分別為修正到首垂線、尾垂線、船舯處的吃水。根據(jù)T′可查得船舶排水量。

1.3.2 函數(shù)修正法

由1.2可知,拱垂變形沿船長(zhǎng)方向?yàn)橐粩M合度非常高的二階函數(shù),而恰好可以根據(jù)六面水尺的觀測(cè)值,以3個(gè)點(diǎn)擬合出1條二次曲線?；诖办o力學(xué),給出如下垂拱對(duì)排水量的修正方法。

b·Aw·LCA+c·Aw

(2)

式中:xmin為船長(zhǎng)方向的最尾端;xmax為船長(zhǎng)方向的最首端;y(x)為沿船長(zhǎng)方向的水線面的寬度;h(x)為沿船長(zhǎng)方向的吃水;a、b、c為二階函數(shù)對(duì)應(yīng)的3個(gè)常數(shù);IX為水線面對(duì)y軸的慣性矩;LCA為飄心的縱向坐標(biāo);Aw為水線面面積。

1.3.3 兩種修正方法的誤差

兩種修正方法的計(jì)算數(shù)據(jù)見表3、4。表中,x11、x12、x13分別為假定的艉部、舯部、艏部水尺位置;t1、t2、t3分別為假定的艉部、舯部、艏部水尺讀數(shù),T為型吃水,T′為根據(jù)經(jīng)驗(yàn)公式修正的型吃水;h1、h2、h3分別為艉垂線、舯部、艏垂線處的基于艏、艉吃水邊線的相對(duì)吃水值;a、b、c為二階函數(shù)對(duì)應(yīng)的3個(gè)常數(shù);IX為水線面對(duì)y軸的慣性矩;LCA為飄心的縱向坐標(biāo);AW為水線面面積;Cb為方形系數(shù);TPC為海水中每厘米吃水噸數(shù);ρ為海水密度;disp0為相關(guān)工況的不計(jì)及拱垂變形的排水量;disp1為經(jīng)驗(yàn)公式修正方法計(jì)算得到的排水量;disp20為采用函數(shù)修正方法得到的初始排水量;disp21為采用函數(shù)修正方法得到的修正排水量,disp2為disp20與disp21之和,Δdisp1為disp1與disp0之差,Δdisp2為disp2與disp0之差。

表3 滿載離港拱垂變形對(duì)排水量的修正

表4 壓載到港拱垂變形對(duì)排水量的修正

從數(shù)據(jù)可知,兩種修正方法的誤差相對(duì)值都比較小,最大不過為0.06%。函數(shù)修正方法更精確;經(jīng)驗(yàn)公式修正方法更簡(jiǎn)易。

2 拱垂變形對(duì)總縱強(qiáng)度的影響

2.1 初始工況

船體梁拱垂變形對(duì)總縱彎矩、剪力的影響,量化滿載離港工況下的靜水(工況A)、坦谷波(等效設(shè)計(jì)波,波高H=10.75 m波長(zhǎng)λ=326.6 m)、坦谷波(波斯灣海域波譜,有義波高Hs=1.672 m跨零周期Tz=4.777 s)3種情況下的修正值,且3種浮態(tài)的排水量及浮心位置是一樣的。需要說明的是,用坦谷波模擬的波浪載荷與規(guī)范給出的用于評(píng)估總縱強(qiáng)度的彎矩、剪力是不同的,亦不同于三維勢(shì)流理論預(yù)報(bào)的彎矩、剪力,在此采用坦谷波法僅是為了簡(jiǎn)化計(jì)算、給出相對(duì)定量的數(shù)據(jù)供參考[6-7]。

上述工況C的取選僅為說明在短波峰海況中船體梁的拱垂變形情況,并非反應(yīng)波斯灣25年一遇海況所對(duì)應(yīng)的等效設(shè)計(jì)波。

2.2 3種工況的浮態(tài)特性

吃水沿船長(zhǎng)方向的分布反應(yīng)了浮力的分布,考慮到重力分布不變,因此吃水的分布某種程度上可以反應(yīng)彎矩、剪力的分布特性。

工況A的吃水見表1。工況B的吃水見表5,工況C的吃水見表6。

表5 工況B的中垂變形

表6 工況C的中垂變形

2.3 彎矩、剪力分布

工況A、B、C的彎矩分布匯總見表7。BMREL代表DFL與NO-DFL之差與許用靜水彎矩的比值;剪力分布見表8,SFREL代表DFL與NO-DFL之差與許用靜水剪力的比值,SFREL′SFREL"含義類似。

表7 3種工況的彎矩分布

表8 3種工況的剪力分布

由上述數(shù)據(jù)可知,在靜水工況(工況A)及短波峰(工況C)情況下,單純的舯垂變形引起的彎矩修正會(huì)使彎矩峰值減少3.6%,約對(duì)應(yīng)許靜水彎矩的2.8%;長(zhǎng)波峰(工況B)情況下,舯垂變形引起的彎矩修正會(huì)使彎矩峰值減少3.8%,約對(duì)應(yīng)許用靜水彎矩的11%左右。

需要特別指出的是,工況B在坦谷波的作用下,舯垂變形約為0.7 m,此工況僅作為極端設(shè)計(jì)工況,真實(shí)的船舶營(yíng)運(yùn)過程中幾乎不可能發(fā)生。

3 結(jié)論

1)對(duì)營(yíng)運(yùn)船而言,優(yōu)化船舶裝載方案可以減少拱垂變形,進(jìn)而增加滿載離港時(shí)載貨量。在保證載貨量不變的情況下,可以將中部貨艙的載貨量移至首尾貨艙,減輕中垂影響,優(yōu)化裝載方案。

2)如何在船舶設(shè)計(jì)階段優(yōu)化船體梁的拱垂變形,由于牽涉到分艙布置、船舶的經(jīng)濟(jì)效益等,目前還難以給出有針對(duì)性的設(shè)計(jì)思路。但是在設(shè)計(jì)裝載方案時(shí)力圖減少彎矩的峰值對(duì)解決拱垂變形大的問題是有利的,且有利于優(yōu)化船體梁自身的總縱強(qiáng)度。

3)船體梁的拱垂變形會(huì)“削低”船體梁在靜水、波浪中的彎矩的峰值。如將來有公認(rèn)的方法來量化評(píng)估船體梁的拱垂變形以及對(duì)應(yīng)的彎矩、剪力,對(duì)于船長(zhǎng)250 m以上、且裝載模式相對(duì)單一的大船的總縱強(qiáng)度校核而言非常重要。