張亞明

（上海振華重工（集團）股份有限公司，上海 200125）

0 引 言

油田守護船是集拖航、起拋錨、供應、油田守護、救助和服務(wù)等多種功能于一體的遠洋船，一般工作于無遮蔽海域，工作環(huán)境受風、浪、流的影響較大。上海振華重工（集團）股份有限公司（以下簡稱“振華重工”）承接有某型6500HP（1hp=735.49875W）油田守護船的設(shè)計、建造合同，由于其2014年建造交付的某型15000HP三用工作拖船的減搖水艙效果不佳，在試航時遭遇7級大風，橫搖角度過大，導致船舶出現(xiàn)甲板物資掉落海中、甲板上浪嚴重等現(xiàn)象，因此在6500HP油田守護船設(shè)計前期，擬選用可自動調(diào)節(jié)液面減搖水艙系統(tǒng)優(yōu)化減搖效果，將船舶的減搖角度和減搖周期控制在最佳范圍內(nèi)，提升船員的舒適度和船舶的作業(yè)效率，保障船舶作業(yè)時人員和設(shè)備的安全[1]。本文針對6500HP油田守護船可自動調(diào)節(jié)液面減搖水艙系統(tǒng)的設(shè)計方案，采用水池模型試驗和計算預報的方法評估減搖水艙的性能；根據(jù)實船試航數(shù)據(jù)，真實客觀地評估減搖水艙的性能和水艙設(shè)計方案的可行性。

1 矩形減搖水艙的理論支撐

1.1 船舶與減搖水艙的力學模型

將在航船舶的船體簡化為自由漂浮的剛體，賦予其6個自由度，分別為橫蕩、縱蕩、艏搖、橫搖、縱搖和垂蕩；減搖水艙就是基于“雙共振”的思路設(shè)計的，即減搖水艙內(nèi)水流的振蕩周期與船舶的固有橫搖周期相同（諧搖現(xiàn)象），且波浪對船舶的橫搖力矩（擾動力矩）和減搖水艙內(nèi)水流對船舶的減搖力矩兩者的相位相反，互相抵消，從而取得較好的減搖效果。圖1為雙重共振時船舶橫搖角、波浪擾動力矩和水艙內(nèi)水流對船舶的減搖力矩三者之間的相位角關(guān)系。

當水艙內(nèi)的流體與船體復合運動時，假定船體繞過重心的縱向軸線做單純的橫搖運動，且運動呈線性，可用無因次線性微分方程描述該線性運動。

1) 對于水艙內(nèi)的壓載水，有

圖1 雙重共振時船舶橫搖角、波浪擾動力矩和水艙內(nèi)水流對船舶的減搖力矩三者之間的相位關(guān)系

2) 對于船殼，有

式(1)和式(2)中：φ為船的橫搖角度；2φν為船的橫搖衰減系數(shù)；φω為橫搖自然頻率；μ為水艙自由液面穩(wěn)性修正因子；ω為來浪圓頻率；eθ為來浪有效波傾角；δ為水艙中流體相對于船體的有效傾角；tω為水艙內(nèi)液面自然頻率；2tν為水艙內(nèi)液體搖蕩衰減系數(shù)；H為水艙內(nèi)流體平均液面與橫搖縱軸線間的距離。式（1）中的˙為水艙中流體的平移加速度，由船體橫搖加速度引起。

定義φa/θe為橫搖放大因數(shù)，εφξ為船體橫搖滯后于波浪的相位角，εφδ為水艙中流體滯后于船體運動的相位角。由式(1)和式(2)可求得橫搖放大因數(shù)和滯后相位角為

式(3)和式(4)中的P和Q分別滿足

在獲取船體和水艙的主尺度要素之后，可由式(3)～式(7)求出船體的橫搖響應幅值算子。

船舶的橫搖自然頻率φω既可由船模試驗確定，又可由式(8)確定，即

式(8)中：D為船舶排水量；h為船舶初穩(wěn)性高；為船體自身慣性矩與附加慣性矩之和，即總慣性矩。

古德雷奇法[2]認為ωφ和ωt的計算式可表示為

1.2 貝胥法研究矩形減搖水艙

若本文中減搖水艙內(nèi)的流體對船體的作用力矩由船模試驗測定，則水艙內(nèi)流體的運動方程式無須考慮；在貝胥法中，當僅考慮船體繞通過重心的縱向中心軸橫搖運動時，可在船體的橫搖運動方程式中體現(xiàn)水艙中流體力矩項，即

式(11)中：I1為船體質(zhì)量慣性矩（含附加質(zhì)量）；2Nφ為船體橫搖阻尼系數(shù)；Ks為船體橫搖復原系數(shù)；φ為船舶橫搖角；ω為來浪圓頻率；Mω為波浪擾動力矩；Mt為水艙力矩。

式(13)中：Kt為水艙內(nèi)流體對船體復原力矩的貢獻系數(shù)；2Nφδ為水艙內(nèi)流體對船體的橫搖阻尼貢獻系數(shù)。將式(12)和(13)代入式(11)中，定義kt=Kt/Ks，消去式中的I1項，整理后得到

式(14)中：2νφ為船體橫搖衰減系數(shù)；2νφδ為水艙流體橫搖衰減系數(shù)。由式(14)可得減搖水艙裝載后的橫搖幅值放大因數(shù)為

由式(15)可知，在求得 2 （νφ+νφδ）和kt之后，即可得到船體橫搖幅值與來浪有效波傾角的比值，即船體橫搖幅值放大因數(shù)。船體自身的橫搖衰減周期和減搖水艙內(nèi)流體的橫搖衰減周期均可通過船模試驗測得。

2 減搖水艙設(shè)計

2.1 水艙布置和內(nèi)部設(shè)置

由式(13)可知：水艙內(nèi)流體對船體復原力矩的貢獻系數(shù)Kt越大，抵消橫搖力矩就越大，因此一般將減搖水艙設(shè)置在遠離橫搖中心軸的位置；同時，需考慮初穩(wěn)性高的因素，即船體重心位置越高，就越低；綜合考慮，一般將矩形減搖水艙布置在主甲板下，且為獲得較大的φ˙，使減搖水艙沿船寬方向的尺度盡可能地大。

本文中6500HP油田守護船主尺度和相應船模尺度參數(shù)見表1。

表1 6500HP油田守護船主尺度和相應船模尺度參數(shù)

減搖水艙的基本結(jié)構(gòu)和安裝位置見圖2。該減搖水艙為矩形，沿船長方向并列布置2個水艙；水艙的尺寸為3.0m×3.0m×14.8m，在主甲板下貫穿船寬方向。減搖水艙模型采用透明的有機玻璃制作，以便清晰地觀察艙內(nèi)液體的運動狀況；船底兩側(cè)舭部帶有舭龍骨。

圖2 減搖水艙的基本結(jié)構(gòu)和安裝位置

模型試驗的內(nèi)容包括減搖水艙模型的岸上試驗、船模的靜水衰減試驗、減搖水艙開啟和閉合狀態(tài)的橫浪規(guī)則波試驗及減搖性能的計算預報。試驗模型的縮尺比為1:20，船模采用木質(zhì)材料制造，有機玻璃減搖水艙模型嵌入式安裝在船模內(nèi)部。模型試驗在上海交通大學的船模拖曳水池中完成。

2.2 減搖水艙周期和阻尼系數(shù)

設(shè)計方案通過減搖水艙岸上試驗（搖擺臺試驗）獲取減搖水艙固有周期和阻尼系數(shù)；只對一個水艙加水，另一個保持空艙。在試驗過程中對艙內(nèi)自由液面施加初始擾動，通過安裝在艙內(nèi)的水位計測量并記錄自由液面的運動衰減曲線，可獲得減搖水艙的固有周期和阻尼系數(shù)，推導如下。

假定水箱橫搖衰減曲線見圖3。在圖3中，假定阻尼線性變化，從t1到t2=t1+Tφ′/ 2半個周期的時間間隔內(nèi)，橫搖幅值的絕對值變化[9]為

圖3 水箱橫搖衰減曲線

搖擺臺試驗是將減搖水艙固定安裝在減搖臺上，通過電機控制完成規(guī)定周期和規(guī)定幅度的簡諧搖擺運動，測量得到減搖水艙對轉(zhuǎn)動軸的橫搖力矩和相位。將該數(shù)據(jù)與船舶原有橫搖RAO曲線[3]相結(jié)合，可預報減搖水艙裝船之后的減搖RAO曲線。試驗時在設(shè)定安裝高度的搖擺臺（轉(zhuǎn)動軸即為船舶重心）上固定住減搖水艙，在不同的設(shè)定周期下用電機驅(qū)動搖擺臺做經(jīng)過水平角度為原點的小角度簡諧振蕩，檢測并記錄3組時歷數(shù)據(jù)，分別為搖擺臺橫搖角度、水艙內(nèi)兩側(cè)水位高和搖擺臺產(chǎn)生的橫搖力矩。

對上述數(shù)據(jù)進行后處理，并將結(jié)果換算到實船尺度。試驗設(shè)置3組水艙水位和3組激勵橫搖角，試驗狀態(tài)見表2。

表2 試驗水位和激勵橫搖角記錄

每組測試設(shè)置10個轉(zhuǎn)速，對應實船周期為5.7～10.86s；本文在6個狀態(tài)的試驗結(jié)果中選取狀態(tài)1情況，具體見表3和圖4。

表3 狀態(tài)1試驗結(jié)果

圖4 狀態(tài)1減搖力矩和相位差隨橫搖周期變化情況

根據(jù)不同艙內(nèi)水位下的橫搖衰減曲線，由式（17）得到的實尺度自由衰減周期和無因次阻尼結(jié)果見表4。

表4 減搖水艙自由衰減測試數(shù)據(jù)

2.3 船模橫搖固有周期和橫搖黏性系數(shù)

船模通過壓載調(diào)整和傾斜試驗，保證吃水、浮態(tài)和橫穩(wěn)性高度符合表1的要求。減搖水艙內(nèi)無自由液面，完成橫搖靜水衰減試驗，得到船模在無減搖狀態(tài)下的固有周期和阻尼數(shù)據(jù)。測試得到的船模橫搖衰減曲線見圖5，得到實尺度橫搖固有周期為8.13s，實尺度橫搖黏性系數(shù)υ=0.022，相關(guān)定義[3]為

3 減搖水艙系統(tǒng)效果評估

圖5 船模橫搖衰減曲線

3.1 規(guī)則橫搖波預報

規(guī)則橫搖波在上海交通大學船模拖曳水池中實現(xiàn)。該拖曳水池全長110m，寬6m，深3m。水池一端的造波裝置可生成規(guī)則長峰波和非規(guī)則長峰波，船模與水池中線垂直放置，來波為橫向規(guī)則波。船模艏艉兩端分別用柔性繩索與水池側(cè)壁連接，限制船模艏搖，但具有橫搖和橫蕩的自由度。

船模橫搖信號采用安裝在模型上的傾角儀測量；波浪信號采用安裝在離船2m位置處的浪高儀測量。所有測量器件都在試驗前標定，零值對應于靜水狀態(tài)輸出值。試驗測量時，所有信號都以50Hz頻率采樣記錄。

模型尺度數(shù)據(jù)可根據(jù)相似律換算成實船數(shù)據(jù)，試驗中各物理量的換算關(guān)系見表5，所有試驗數(shù)據(jù)都以實船數(shù)據(jù)的形式表達。

表5 實船數(shù)據(jù)換算關(guān)系

3.2 橫搖RAO分析

采用頻域分析方法進行規(guī)則橫搖波試驗，在系列規(guī)則波中試驗結(jié)果的一般表達為橫搖響應幅值算子RAO[3]，定義為橫搖幅值與波浪幅值的比值。橫搖RAO值隨波浪周期或波頻變化的曲線稱為RAO曲線。橫搖 RAO曲線既能用來預報不同波浪譜中的橫搖統(tǒng)計值，又能直觀反映減搖水艙的減搖性能。

橫搖規(guī)則波試驗狀態(tài)見表6，在同一浪向下分別完成減搖水艙關(guān)閉和開啟的2組規(guī)則波試驗，得到橫浪狀態(tài)下的橫搖RAO曲線對比見圖6。2種狀態(tài)下的橫搖規(guī)則波試驗結(jié)果分別見表7和表8。

圖6 橫浪狀態(tài)下的橫搖RAO曲線對比

表6 橫搖規(guī)則波試驗狀態(tài)

表7 狀態(tài)1的橫搖規(guī)則波試驗結(jié)果

表8 狀態(tài)2的橫搖規(guī)則波試驗結(jié)果

3.3 非規(guī)則波中橫搖響應預報

根據(jù)船舶響應的線性假設(shè)，可采用頻譜計算方法預報運動響應的有義雙幅值。當波浪譜和運動響應算子RAO已知的時，運動響應的頻譜的計算式[5]可表示為

若運動響應頻譜曲線與橫軸之間封閉區(qū)域的面積為m0，則運動響應的有義雙幅值Re的計算式可表示為

根據(jù)表7和表8中的數(shù)據(jù)，可計算得到在不同海浪中減搖關(guān)閉和減搖開啟2種狀態(tài)對應的橫搖有義值對比見表9。

表9 橫搖有義雙幅值對比

3.4 試航數(shù)據(jù)評估非規(guī)則波橫搖響應

試航海域為舟山北部海域，浪高1.4～2.0m（4級，中浪），譜峰周期7.2s；傾角通過懸掛于船舶駕駛室的傾斜儀測量，抗橫搖水艙系統(tǒng)控制面板位于船舶駕駛室前駕控臺，通過放置在水艙內(nèi)的液位計和陀螺儀實時監(jiān)測船體的狀態(tài)和水艙液位狀態(tài)；經(jīng)數(shù)據(jù)整理，選取0～4kn航速數(shù)據(jù)和12kn航速數(shù)據(jù)，得到試航減搖效果見表10。

表10 試航減搖效果

由表10可知：在低航速（0～4kn）下，減搖水艙系統(tǒng)實際減搖效果為50%，與非規(guī)則波橫搖響應預報結(jié)果53%非常接近；在高航速（12kn）下，減搖水艙系統(tǒng)實際減搖效果47%，與理論預報結(jié)果53%相差較小，且親身體驗水艙工作時的橫搖角度在可接受的范圍內(nèi)。由此可知：6500HP項目的減搖水艙系統(tǒng)設(shè)計是成功的。

4 結(jié) 語

根據(jù)模型試驗結(jié)果和橫搖預報結(jié)果，可得到以下結(jié)論：

1) 從橫搖RAO曲線對比中可看出，減搖水艙開啟可有效改善船舶橫搖性能。

2) 從橫搖預報結(jié)果來看，在4級浪和5級浪狀態(tài)下，減搖水艙開啟后的減搖效果可達到50%以上；從實際試航結(jié)果來看，減搖效果非常接近50%。

3) 從減搖水艙自由衰減結(jié)果中可看出，減搖水艙自由衰減周期對艙內(nèi)水位高度相當敏感。因此，在船舶裝載狀態(tài)變化的情況下，必須及時調(diào)整水位，以保證減搖水艙周期與船舶橫搖固有周期匹配。

【 參 考 文 獻 】

[1] 李積德. 船舶耐波性[M]. 哈爾濱：哈爾濱工程大學出版社，1992.

[2] 陶堯森. 船舶耐波性[M]. 上海：上海交通大學 1985.

[3] 肖龍飛，彭濤，楊建民，等. 被動式減搖水艙減搖效果的評定[J]. 上海：上海交通大學學報，2001 (10): 1552-1556.

[4] 肖麗娜. 平面被動槽型減搖水艙的研究[D]. 大連：大連理工大學，2005.

[5] 謝廣東. 船舶減搖水艙艙液振蕩周期和阻尼特性的研究[D]. 哈爾濱：哈爾濱工程大學，2011.

[6] 洪超，陳瑩霞. 船舶減搖技術(shù)現(xiàn)狀及發(fā)展趨勢[J]. 船舶工程，2012, 34 (S2): 236-245.

[7] 張松濤. 減搖水艙優(yōu)化設(shè)計及試驗研究[D]. 哈爾濱：哈爾濱工程大學，2006.

[8] 吳文成，鄒康. 深水供應船耐波性及減搖效果研究[J]. 船舶與海洋工程，2016, 32 (1): 36-40.

[9] 盛振邦，劉應中. 船舶原理[M]. 上海：上海交通大學 2003.