朱軍，劉瑞杰，葛義軍，王智宇，黃昆侖

1海軍工程大學艦船工程系，湖北武漢430033

2海軍裝備研究院，北京100073

規(guī)則波浪中艦船純穩(wěn)性喪失計算研究

朱軍1，劉瑞杰2，葛義軍1，王智宇1，黃昆侖1

1海軍工程大學艦船工程系，湖北武漢430033

2海軍裝備研究院，北京100073

［目的］波浪中純穩(wěn)性喪失是第二代艦船完整穩(wěn)性中穩(wěn)性薄弱性的衡準之一。針對波浪中大傾斜角穩(wěn)性計算發(fā)散問題，［方法］以靜水面坐標系為基準，定義了廣義縱傾角和廣義吃水變量，推導出物理含義清晰的波面方程。在Froude-Krylov假定條件下，結合AutoCAD二維圖形面域計算技術和VBA編程方法，提出了規(guī)則波浪中艦船純穩(wěn)性喪失的計算方法。針對一艘具有階梯式甲板的艦船，計算了規(guī)則波浪中艦船復原力臂曲線，證明大橫傾狀態(tài)具有一致收斂性。［結果］計算得到了規(guī)則波浪中穩(wěn)性變化規(guī)律：縱向波浪中，波面高于甲板或低于船底導致有效水線面消失是穩(wěn)性降低的主要原因；斜浪和正橫浪中，波面相對船體橫剖面傾斜引起的不對稱性，是穩(wěn)性大幅度降低的主要原因。［結論］計算收斂一致性表明，基于廣義縱傾角定義的計算方法可成為評估艦船波浪中純穩(wěn)性喪失的有效手段。

艦船；純穩(wěn)性喪失；規(guī)則波浪

0 引 言

20世紀30年代，Kempf［1］通過不斷調查和分析發(fā)現，當船舯位于波峰附近時，船舶穩(wěn)性相比靜水會有所下降，反之，當船舯位于波谷附近時，穩(wěn)性有所上升。1961年，Paulling等［2-3］研究了波浪中船舶的穩(wěn)性，其后通過在舊金山灣中進行隨浪、尾斜浪和橫浪環(huán)境中的大尺度比船舶試驗，并結合時域模擬分析方法，發(fā)現了船舶在波浪環(huán)境中的傾覆機理。Grim［4］提出了有效波的觀點，將海洋中的隨機波轉化為規(guī)則波。上世紀80年代，Hamamoto等［5-6］對尾斜浪的船舶穩(wěn)性計算方法進行了研究，但是大橫傾角計算時出現發(fā)散情況，隨后他建立了一個新的坐標轉換體系研究船舶在波浪中的操縱性［7］。Fang等［8］在 Froud-Krylov假定的基礎上考慮散射和繞射的影響，計算了船舶在斜浪中的復原力臂。Kreuzer等［9］在考慮船舶橫搖、縱搖及垂蕩的基礎上對船舶傾覆問題進行了研究。

國內孔祥金等［10］為了解決大橫傾角下計算的發(fā)散問題，采用攝動法計算波面下船體橫剖面積和中心；謝微［11］對穩(wěn)性薄弱性的純穩(wěn)性喪失進行了研究；周耀華等［12-13］對純穩(wěn)性喪失薄弱性的一、二層衡準進行了計算分析；胡麗芬等［14］綜述了參數激振橫搖的研究現狀。朱軍等［15］應用AutoCAD的圖形面域技術數值計算了艦船在靜水中的穩(wěn)性與破損穩(wěn)性。

研究發(fā)現，波浪中大橫傾角的穩(wěn)性計算出現發(fā)散的原因，在于繞船體坐標軸定義的縱傾角與縱向傾斜力矩（重力與浮力構成的力矩）方向不一致，橫傾角越大偏差也越大。在橫傾角為90°的極端情況下，縱傾角在水平面內轉動，而縱向傾斜力矩則始終是垂直水平面的。這導致大橫傾角下縱傾力矩平衡方程式的不匹配。

本文將定義繞靜水面軸線轉動的廣義縱傾角，以及垂直靜水面的廣義吃水，保證縱向傾斜力矩平衡方程式的匹配性。

1 船體坐標系下的波面方程

1.1 坐標系定義

靜水面固定坐標系O-xyz如圖1所示，平面Oxy與靜水面重合，Oz坐標軸垂直向上，波面方程為

式中：λ為波長；Hw為波高；χ為浪向角；ξ0為波峰位置。當ξ0=0時，波峰線通過坐標原點，當ξ0=λ/2時，波谷線通過坐標原點。

圖1 靜水面固定坐標系Fig.1 Fixed coordinate system of static surface

船體固定坐標系B-xByBzB如圖2所示，以左右對稱面、龍骨基面和船舯面的3面交點B為坐標系原點，BxB軸指向船艏，ByB軸指向右舷，BzB軸垂直向上。T為靜水船舯處的吃水，得到靜水面固定坐標系與船體坐標系的轉換關系為

圖2 船體坐標系Fig.2 Coordinate system on hull

令靜水面坐標系繞O點轉動?角度（圖3（a）），再令靜水面坐標系繞沿靜水面平移ΔT（圖3（b），相當于船體下沉），最后令靜水面坐標系繞Oy軸旋轉θ角（相當于船體繞靜水線面軸線轉動），則轉動后兩坐標系的轉換關系為

圖3 坐標系旋轉與平移Fig.3 Coordinate systems transformed by rotation and shift

1.2 船體坐標系下波面方程

令式（3）的第3項z=0，則到船體坐標下靜水面方程為

將式（1）代入到式（3）的第3項中，得到船體坐標下的波面方程為

式中，x和y由式（3）的第1項和第2項確定，這是一個隱式方程。

式（5）為波面在船體坐標系下的表達式，當吃水T、吃水增量ΔT、橫傾角?和縱傾角θ一定時，波面方程唯一確定。值得注意的是，T和?與通常的定義是一致的，但吃水增量ΔT是垂直靜水面的變化量，即廣義吃水，θ則是繞船舯剖面與靜水面交線的轉動角度，即廣義縱傾角。

2 復原力臂計算公式

2.1 艦船平衡條件

平衡條件為艦船重量等于瞬時波面下的浮力，且作用在同一條鉛垂線上（垂直靜水線面），平衡方程為

式中：ρ為水的密度；g為重力加速度；為瞬時波面下船體的排水體積；(xBb，yBb，zBb)為船體排水體積中心；W和(xBg，yBg，zBg)為船體的重力和重心。排水體積和浮心坐標按二維切片法計算：

式中：A(xB)為瞬時波線下船體橫剖面面積；yBc，zBc為對應的面積中心。

2.2 瞬時波面下船體復原力臂

當船體橫傾角為?時，在滿足式（7）中的第1式（重力與浮力平衡）和第3式（縱傾力矩平衡）條件下，瞬時波面下船體復原力臂GZ為

上述平衡條件和復原力臂計算式是基于瞬時波面下靜水壓力（Froude-Krylov）的假定，忽略了史密斯的動壓修正項（e-kz）。為了簡化計算，取波面以下的排水重量作為瞬時浮力，將流體壓力近似成p=ρg[z+?zcos(kx)]。該近似相當于在波峰區(qū)域高估了壓力，在波谷區(qū)域低估了壓力，即有效波陡應該比幾何波陡略小，通常的近似方法是修正幾何波陡，即動壓部分取為?ze-kd/2cos(kx)，其中d為吃水。因此，修正將不會影響運動的變化規(guī)律。

3 數值計算及分析

3.1 圖形面域計算技術

本文應用AutoCAD二維圖形面域的計算方法［15］，通過編程語言VBA二次開發(fā)了波浪中艦船穩(wěn)性計算程序。計算步驟如下：

1）繪制船體橫剖面圖，形成船體剖面圖形面域；

2）按式（5）在對應的剖面處繪制波線，形成包含船體剖面的波線圖形面域；

3）提取船體面域和波線圖形面域并集的面域特征值（面積和形心坐標）；

4）按式（7）積分得到排水體積和體積中心坐標，在滿足式（6）第1式和第3式的條件下，由式（8）計算復原力臂。

3.2 波浪中穩(wěn)性數值計算

本文計算艦船的水線長為125 m，艏樓計入穩(wěn)性部分，艉部有階梯型甲板。波浪中穩(wěn)性計算范圍為：浪向角χ分別為 0°，30°，60°和90°；波長船長比λ/L分別為 0.5，1.0和 1.5；波陡Hw/λ為0.01～0.1，波峰相位(2π/λ×ξ0)為0～2π。

3.3 艦船波浪中穩(wěn)性分析

1）縱向波浪中純穩(wěn)性喪失計算分析。

圖4為縱向波浪中（浪向角χ=0°）復原力臂計算曲線圖。波峰相位(2π/λ×ξ0)=0表示波峰位于船舯，(2π/λ×ξ0)=π表示波谷位于船舯。計算結果表明，波高顯著導致穩(wěn)性喪失，復原力臂的最大值和穩(wěn)性范圍均顯著減少；另一個特征是中垂狀態(tài)（波谷位于船舯）也有顯著的穩(wěn)性喪失。

圖4 縱向波浪中復原力臂計算曲線（χ=0°）Fig.4 Calculating curves of righting arm in longitudinal waves（χ=0°）

當λ/L=1.0和Hw/λ=0.06時，船體與波面相交位置如圖5所示。說明波峰位于船舯時，有效水線面（提供穩(wěn)性）位于船體中部，船艏和船艉的底部高出波面導致有效水線面喪失；波谷位于船舯時，只有船舯后和船舯前2部分區(qū)域為有效水線面，船舯部因波谷低于船體而喪失了有效水線面，船艏和船艉上甲板因波浪淹沒也導致有效水線面的喪失。

圖5 船體與波面相交位置圖（λ/L=1.0，Hw/λ=0.06）Fig.5 The interface between hull and waves on crest and trough（λ/L=1.0，Hw/λ=0.06）

2）正橫波浪中純穩(wěn)性喪失計算分析。

圖6為正橫波浪中（浪向角χ=90°）復原力臂計算曲線圖，圖6給出了波峰相位分別為0.5π和0.75π下2種穩(wěn)性喪失較為嚴重的情況。計算結果表明，波高也是導致穩(wěn)性喪失的主要影響因素，穩(wěn)性喪失表現在最大復原力臂值隨波高增加而大為減少，其次是穩(wěn)性范圍也大為減少。由于波面相對船體橫剖面成傾斜狀，類似于船體存在初始橫傾角，復原力臂曲線呈現左右不對稱性，導致最大復原力臂值和穩(wěn)性范圍減少。

3）斜浪中純穩(wěn)性喪失計算分析。

圖7為斜浪中（浪向角χ=30°）復原力臂計算曲線圖。計算結果表明，斜浪狀態(tài)的穩(wěn)性特征介于縱向波浪和正橫浪之間。波長較短時復原力臂左右不對稱性不明顯，這是由于波長較短時，波面相對船體橫剖面的傾斜方向相反，抵消了傾斜的不對稱性。隨著波長的增加，傾斜的不對稱性變得顯著，最大復原力臂值隨波高增加明顯降低。

圖6 正橫波浪中復原力臂計算曲線（χ=90°）Fig.6 Calculating curves of righting arm in beam seas（χ=90°）

圖7 斜浪中復原力臂計算曲線（χ=30°）Fig.7 Calculating curves of righting arm in oblique waves（χ=30°）

4 結 論

本文定義了廣義縱傾角θ和廣義吃水ΔT，推導出了船體坐標系下的波面方程。導出的船體坐標系下的波面方程由5項構成，各物理含義明確，且吃水增量ΔT、橫傾角?和縱傾角θ分別主要影響到浮力、橫向力矩和縱向力矩，從而使得迭代計算具有穩(wěn)定收斂特性。在Froude-Krylov假定條件下，對階梯式甲板艦船波浪中復原力臂計算結果的對比和分析表明：

1）λ/L=1.0是艦船波浪中純穩(wěn)性喪失的主要波浪環(huán)境，波高是導致穩(wěn)性喪失的主要因素。

2）縱向波浪中，波谷位于船舯時（某些波高較大時）穩(wěn)性也會大幅度降低，波面高于甲板或低于船底導致有效水線面消失是其主要原因。

3）正橫浪和斜浪情況下，因波面相對船體橫剖面的傾斜形成的不對稱性是穩(wěn)性大幅度喪失的主要原因，正橫浪情況最為嚴重。

本文提出的規(guī)則波浪中艦船純穩(wěn)性喪失直接計算的方法具有穩(wěn)定收斂特性，且可擴展至船舶搖蕩耦合運動計算方法中。

［1］KEMPF G.Die stabilit?tsbeanspruchung der schiffe durch wellen und schwingungen［J］.Werft Reederei Hafen，1938，19：200-202.

［2］PAULLING J R.The transverse stability of a ship in a longitudinal seaway［J］.Journal of Ship Research，1961，4（4）：37-49.

［3］PAULLING J R，KAESTNER S，SCHAFFRAN S.Ex?perimental studies of capsizing of intact ships in heavy seas［R］.U.S.Coast Guard Technical Report，1972.

［4］GRIM O.Beitrag zu dem problem der Sicherheit des Schiffes im seegang［J］.Schiff und Hafen，1961，6：490-497.

［5］HAMAMOTO M，UMEDA N，SHIGEHIRO R，et al.Transverse stability of a ship in following sea［J］.The Japan Society of Naval Architects and Ocean Engi?neers，1982，185：49-56.

［6］浜本剛実，金潤洙，松田秋彥，等.斜め追波中の船の転覆とその原因の分析［J］.日本造船學會論文集，1992，172：135-145.HAMAMOTO M，KIM Y S，MATSUDA A，et al.An analysis of a ship capsizing in quartering seas［J］.The Society of Naval Architects of Japan，1992，172：135-145（in Japanese）.

［7］浜本剛実，金潤洙.波浪中の操縦運動を記述する新しい座標系とその運動方程式［J］.日本造船學會論文集，1993，173：209-220.HAMAMOTO M，KIM Y S.A new coordinate system and the equations describing manoeuvring motion of a ship in waves［J］.The Society of Naval Architects of Ja?pan，1993，173：209-220（in Japanese）.

［8］FANG M C，LEE C K.On the dynamic stability of a ship advancing in longitudinal waves［J］.International Shipbuilding Progress，1993，40（422）：177-197.

［9］KREUZER E，WENDT M.Ship capsizing analysis us?ing advanced hydrodynamic modelling［J］.Philosophi?cal Transactions of the Royal Society A：Mathemati?cal，Physical and Engineering Sciences，2000，358（1771）：1835-1851.

［10］孔祥金，曹振海.尾斜浪中船舶復原力計算［J］.中國造船，1994（2）：32-41.KONG X J，CAO Z H.Calculation of righting arm for ships in quartering seas［J］.Shipbuilding of China，1994（2）：32-41（in Chinese）.

［11］謝微.船舶純穩(wěn)性喪失薄弱性衡準研究［D］.大連：大連理工大學，2014.XIE W.Research on the draft vulnerability criteria of the pure loss of stability of ship［D］.Dalian：Dalian University of Technology，2014（in Chinese）.

［12］周耀華，張高峰，牛業(yè)興，等.純穩(wěn)性喪失衡準實船分析研究［J］.中國造船，2015，56（增刊 1）：37-47.ZHOU Y H，ZHANG G F，NIU Y X，et al.Study on full-scale ships based on pure loss of stability criteria［J］.Shipbuilding of China，2015，56（Supp 1）：37-47（in Chinese）.

［13］馬坤，甘鄂，謝微.船舶純穩(wěn)性喪失薄弱性衡準樣船計算分析［J］.中國造船，2015，56（增刊 1）：193-200.MA K，GAN E，XIE W.Sample calculations and analysis on vulnerability criteria of pure loss of stabili?ty［J］.Shipbuilding of China，2015，56（Supp 1）：193-200（in Chinese）.

［14］胡麗芬，魯江，顧民.船舶波浪穩(wěn)性研究現狀分析［J］.中國造船，2015，56（增刊1）：211-216.HU L F，LU J，GU M.Status analysis of research on ship's stability in waves［J］.Shipbuilding of China，2015，56（Supp 1）：211-216（in Chinese）.

［15］朱軍，黃昆侖.艦船穩(wěn)性與破損穩(wěn)性的二維圖形面域計算法［C］//第九屆全國海事技術研討會論文集.上海：上海市海事交流協會，2003：99-104.ZHU J，HUANG K L.The calculating method of tow-dimension region for ship stability and dam?aged-stability［C］//Proceedings of the Ninth National Conference on Maritime Technology.Shanghai：ShanghaiMaritimeExchangeAssociation， 2003：99-104（in Chinese）.

Pure loss of stability calculation of naval ship in regular waves

ZHU Jun1，LIU Ruijie2，GE Yijun1，WANG Zhiyu1，HUANG Kunlun1
1 Department of Naval Architecture Engineering，Naval University of Engineering，Wuhan 430033，China
2 Naval Academy of Armament，Beijing 100073，China

The pure loss of stability of a ship is one of the failure modes of the second-generation intact stability criteria.Based on the static surface coordinate system,the generalized pitch angle and draught variable are defined,and a wave equation with clear physical meaning is deduced.Based on the Froude-Krylov hypothesis,combined with AutoCAD 2D surface area computing technology and the VBA programming method,a calculation method of the loss of pure stability of ships in regular waves is proposed.For a naval ship,the righting arm is calculated in regular waves,and the large heel ship state is shown to have an identical convergence.The results of the calculations show that a significant reduction in the meta centric height of the maximum righting arm occurs not just in a wave crest but also in a trough.By analyzing the wave profile under the hull,it can clearly be seen that wave amplitude above the deck or below the bottom of the hull causes the pure loss of stability,and in oblique waves or beam seas,the pure loss of stability is caused the asymmetry of the wave profile on the hull.The coinciding convergence of the calculations shows that the process of the definition of generalized pitch can be employed to assess the pure loss of stability of naval ships in regular waves.

naval ship；pure loss of stability；regular waves

U661.2+2

：ADOI：10.3969/j.issn.1673-3185.2017.03.001

http://kns.cnki.net/kcms/detail/42.1755.TJ.20170512.1303.040.html期刊網址：www.ship-research.com

朱軍，劉瑞杰 ，葛義軍，等.規(guī)則波浪中艦船純穩(wěn)性喪失計算研究［J］.中國艦船研究，2017，12（3）：1-6，15.

ZHU J，LIU R J，GE Y J，et al.Pure loss of stability calculation of naval ship in regular waves［J］.Chinese Journal of Ship Research，2017，12（3）：1-6，15.

2017-02-16< class="emphasis_bold">網絡出版時間

時間：2017-5-12 13:03

朱軍，男，1959年生，博士，教授。研究方向：船舶操縱性。E-mail：zhjun101@sina.com

劉瑞杰（通信作者），男，1987年生，博士，工程師。研究方向：船舶水動力。E-mail：jackrygy@126.com