劉 洋，宋曉東，李春來

(1. 海軍大連艦艇學(xué)院，遼寧 大連 116000；2. 中國人民解放軍 91421 部隊，山東 青島 266011；3. 海軍工程大學(xué)，湖北 武漢 430033)

潛艇外掛武器模塊后的運動建模與仿真

劉 洋1，宋曉東2，李春來3

(1. 海軍大連艦艇學(xué)院，遼寧 大連 116000；2. 中國人民解放軍 91421 部隊，山東 青島 266011；3. 海軍工程大學(xué)，湖北 武漢 430033)

為了對潛艇外掛武器模塊后的運動性能進行研究，本文用近似估算的方法推算外掛后潛艇的水動力系數(shù)，建立六自由度運動的 Matlab 數(shù)學(xué)模型，并對潛艇外掛武器模塊前后的5種典型運動進行對比仿真分析。仿真結(jié)果表明外掛武器模塊對潛艇的運動性能有一定的影響，外掛后潛艇的操縱性降低、應(yīng)舵性下降，所建立的仿真模型可基本反映潛艇外掛武器模塊后的運動特性，可為下一步潛艇外掛武器模塊論證奠定一定的理論基礎(chǔ)。

潛艇；外掛武器模塊；仿真；水動力系數(shù)估算

0 引 言

潛艇外掛武器模塊（模塊內(nèi)可為 UUV、導(dǎo)彈、魚雷等），對提升我海軍現(xiàn)在及未來水下作戰(zhàn)能力具有重大意義。有關(guān)潛艇外外掛武器的應(yīng)用主要有布雷箱和 UUV[1–3]，在潛艇外部外掛研究方面，國內(nèi)主要對潛器攜帶吊艙的水動力系數(shù)進行計算[4]。而在潛艇空間操縱仿真方面，國外已在潛艇水下操縱性、航行的安全性和攻擊規(guī)避的有效性、陸上操縱仿真訓(xùn)練和綜合顯示集中控制涉及的相關(guān)理論和研究方面達到了一定的研究水平[5–7]，國內(nèi)也取得了很多研究成果[8]。本文將在這些成果的基礎(chǔ)上，針對潛艇外外掛武器模塊后的運動進行建模仿真和初步分析。

1 潛艇運動建模與仿真

1.1 模型的建立與求解

建立如圖1所示的固定坐標(biāo)系 E-ξηζ（簡稱“定系”）和運動坐標(biāo)系 G-xyz（簡稱“動系”）。假設(shè)潛艇重心處相對于地球的速度為 V，V 在動系上的投影為 u（縱向速度）、v（橫向速度）、w（垂向速度）；潛艇以角速度 Ω 轉(zhuǎn)動，Ω 在動系上的投影為p（橫傾角速度）、q（縱傾角速度）、r（偏航角速度）。同時規(guī)定：橫傾角 φ 向右傾為正，首向角 ψ 向右轉(zhuǎn)為正，縱傾角 θ 向尾傾為正。方向舵角 δr右舵為正，首舵角 δb和尾舵角 δs的正負(fù)號按右手定則決定，則有：首舵角 δb為正是上浮舵角，尾舵角 δs為正是下潛舵角。

在所建立的坐標(biāo)系基礎(chǔ)上，依據(jù)葛特勒（Gertler）等提出的“用于潛艇模擬研究的標(biāo)準(zhǔn)運動方程”[9]，省略螺旋槳負(fù)荷的影響，同時將微分項移到方程的左邊，其余項放在方程右邊，從而得到六自由度空間運動方程如下：

縱向力方程

橫向力方程

豎向力方程

橫傾力矩方程

縱傾力矩方程

偏航力矩方程

式中：

上式中各符號意義參見文獻[9]。

聯(lián)立式（1）～ 式（6）求解可以得到加速度矩陣：

式中未賦值 aij的值為 0。

對式（7）矩陣中的加速度進行積分得到速度和角速度，根據(jù)定系和動系之間的轉(zhuǎn)換關(guān)系再次積分可求得潛艇任意時刻的位置和姿態(tài)角。

1.2 模型的驗證

為了檢驗所建立的數(shù)學(xué)模型的準(zhǔn)確性，根據(jù)數(shù)學(xué)模型編制仿真程序，對潛艇的 5 種典型運動進行驗證。

1）水平面定?；剞D(zhuǎn)運動

潛艇直航時，當(dāng)方向舵保持某一固定的舵角，潛艇將作圓周運動，稱為水平面定?；剞D(zhuǎn)運動。對潛艇操 15°，25°，35° 方向舵時的定?；剞D(zhuǎn)運動進行仿真，得到潛艇重心軌跡曲線如圖 2 所示。

2）水平面Z型操舵運動

潛艇在實際的航行中，常常是在保持航向的基礎(chǔ)上執(zhí)行操舵機動，這種中等舵角時的機動性可以通過Z型操舵來研究，其重點是潛艇的應(yīng)舵性。對潛艇在δr/ψ = 10°/10° 時的Z型操舵運動仿真，得到的舵角、首向角變化曲線如圖 3 所示。

3）垂直面定常直線潛浮運動

潛艇在定深航向時，當(dāng)升降舵變?yōu)槟骋还潭ǘ娼?，潛艇將作定常直線潛浮運動。定常潛浮運動能夠反映潛艇的深度機動性。對潛艇操 10° 和 20° 尾升降舵的定常直線潛浮運動仿真，計算所得垂直面重心軌跡曲線如圖 4 所示，縱傾角變化曲線如圖 5 所示。

4）垂直面超越機動運動

垂直面內(nèi)周期性變化潛艇升降舵，潛艇將作超越機動運動，此時可以研究潛艇對升降舵的應(yīng)舵能力。對潛艇在 δs/θ = 20°/10° 時的超越機動運動，計算所得舵角、縱傾角變化曲線如圖 6 所示，潛深變化曲線如圖 7 所示。

5）空間定常螺旋運動

當(dāng)固定潛艇某一方向舵角和升降舵角，或者只固定方向舵角，潛艇將繞 Eζ 軸以等角速度回轉(zhuǎn)，其重心的三維運動軌跡為螺旋線，此時潛艇的運動稱為空間定常螺旋運動?？臻g定常螺旋運動可用來研究潛艇的 空間機動性。對潛艇在操舵條件：δr= 30°，δb= 10° 及δs= 0° 時的空間定常螺旋運動進行仿真，得到的重心軌跡曲線如圖 8 所示。

5 種典型運動的仿真結(jié)果表明：仿真結(jié)果反映出來的潛艇運動規(guī)律基本符合客觀實際，所建立的潛艇六自由度數(shù)學(xué)模型基本正確，可用來預(yù)報外掛武器模塊之后潛艇的運動。

2 潛艇外掛武器模塊后的水動力系數(shù)近似估算

潛艇外掛載武器模塊后的水動力系數(shù)會發(fā)生改變，為此必須要確定外掛的武器模塊的水動力系數(shù)。潛艇外掛武器后的水動力系數(shù)獲取可采用實驗測定、理論計算或近似估算等方法[10]，本文在初步論證階段，采用近似估算方法。近似估算方法[11]是根據(jù)在大量船模試驗結(jié)果基礎(chǔ)上整理得出的半經(jīng)驗的圖譜經(jīng)過估算得到的，能在允許一定誤差的情況下較快的推算出水動力系數(shù)的近似值。這種方法假設(shè)潛艇外掛武器模塊后的水動力系數(shù)等于潛艇和外掛的武器模塊的水動力系數(shù)之和。

2.1 武器模塊的加速度系數(shù)近似估算

潛艇外掛的武器模塊一般可有一個或多個。本文所研究的為潛艇外部外掛一個武器模塊情況（見圖 9）。武器模塊長 L 取 10 m，寬 B 和高 H 取 1.5 m，水下全排水容積 ? 取 8 m3。

流線形的武器模塊可以用 1 個三軸橢球體代替，以主尺度 L，B，H 作為橢球體的三軸 2a，2b，2c。由于橢球體的 λij= 0（i ≠ j），故所需計算的附加質(zhì)量只有 λij= 0（i = j）。

橢球體的附加質(zhì)量系數(shù)曲線是蘭姆（Lamb）根據(jù)勢流理論得到的[11]。當(dāng)已知橢球體的半軸長 a，b，c 時，在圖 10 曲線就可以查得橢球體的附加質(zhì)量。圖中的 K11，K22，K33分別是橢球體附加質(zhì)量 λ11，λ22，λ33的相對值，K55，K66為橢球體附加轉(zhuǎn)動慣量 λ55，λ66的相對值。

由 L/B = 10.0，H/B = 1.0 查橢球體的附加質(zhì)量系數(shù)曲線（見圖 10），得到附加質(zhì)量系數(shù)：K11= 0.02，K22= 0.96，K33= 0.96，K55= 0.87，K66= 0.87，將得到的質(zhì)量系數(shù)通過計算得到武器模塊的加速度系數(shù)：

2.2 武器模塊的速度系數(shù)近似估算

橢球體的阻力系數(shù)取

計算得到武器模塊的速度水動力系數(shù)：

2.3 武器模塊的角速度系數(shù)近似估算

假設(shè)武器模塊的重心 G 在模塊正中，即浮心位置修正系數(shù)為 1，則武器模塊的角速度系數(shù)為：

3 外掛武器模塊后潛艇的運動仿真

根據(jù)估算得到的外掛后潛艇的水動力系數(shù)，結(jié)合建立的潛艇運動數(shù)學(xué)模型，便可以對外掛后的潛艇 5種典型運動進行仿真。

3.1 水平面定?；剞D(zhuǎn)運動對比

對潛艇操 15°，25°，35° 方向舵 3 種情況分別進行運動仿真，這里給出 15° 方向舵時外掛前后的潛艇重心軌跡曲線如圖 13 所示。

通過重心軌跡曲線計算得到外掛武器模塊前后潛艇的操縱性能指標(biāo)如表 1 所示。

表1 外掛前后的水平面操縱性能指標(biāo)對照表Tab. 1 Contrast on performance indexes of maneuverability in horizontal plane before & after carrying weapons module

從表 1 可知：潛艇外掛武器模塊后的回轉(zhuǎn)直徑、潛艇從打舵開始至達到定常狀態(tài)的穩(wěn)定時間、回轉(zhuǎn)周期以及回轉(zhuǎn)角速度等操縱性能指標(biāo)均比操舵前的指標(biāo)大，說明外掛武器模塊后潛艇的水平面操縱性能減低。另外，通過外掛前后的性能指標(biāo)的絕對變化量分析，潛艇在小舵角時外掛前后的操縱性變化較大、在大舵角時的變化較小。

3.2 水平面 Z 型操舵運動對比

仿真得到的潛艇外掛武器模塊前后的首向角對比曲線如圖 14 所示。

計算結(jié)果表明：外掛武器模塊前潛艇的初轉(zhuǎn)期為31.0 s、超越首向角為 11°，外掛后的潛艇的初轉(zhuǎn)期為33.7 s、超越首向角為 12°。外掛前后初轉(zhuǎn)期的絕對變化量為 8.71%、超越首向角的絕對變化量為 4.76%。分析可知外掛武器模塊后潛艇的應(yīng)舵性降低了。

3.3 垂直面定常直線潛浮運動對比

對潛艇（尾升降舵 20°）在外掛武器模塊前后的定常直線潛浮運動進行仿真，得到的重心軌跡對比曲線（見圖 15）和縱傾角對比曲線（見圖 16）。

從圖 15 和圖 16 分析可知：潛艇在外掛前的穩(wěn)定時間為 46.2 s、最大縱傾角為 59°；外掛后的穩(wěn)定時間為 50.1 s、最大縱傾角為 52°。分析可知外掛的武器模塊對潛艇在垂直面上的操縱性能指標(biāo)也有較大的影響，特別是縱傾角變化很大。

3.4 垂直面超越機動運動對比

對潛艇外掛前后在 δs/θ = 20°/10° 時的超越機動運動仿真計算，所得的深度對比曲線和縱傾角對比曲線如圖 17 和圖 18 所示。

計算結(jié)果表明：外掛武器模塊前潛艇的初轉(zhuǎn)期為17.8 s、超越深度為 0.4 m、超越縱傾角為 2.5°，外掛后的初轉(zhuǎn)期為 20.0 s、超越深度為 0.5 m、超越縱傾角為3.2°。計算可知搭載前后初轉(zhuǎn)期的絕對變化量為15.7%、超越深度的絕對變化量為 25.0%、超越縱傾角的絕對變化量為 28.0%。分析各性能指標(biāo)的絕對變化量可知搭載的武器模塊對潛艇垂直面的應(yīng)舵性影響比對水平面的影響更大，潛艇在垂直面的機動能力將損失更多。

3.5 空間定常螺旋運動對比

在 δr= 30°，δb= 10°，δs= 0° 操舵條件下，對潛艇外掛武器模塊前后的空間定常螺旋運動進行仿真，得到的潛艇外掛武器模塊前后的重心運動軌跡對比曲線如圖 19 所示。

仿真得出的潛艇外掛武器模塊前后的回轉(zhuǎn)直徑、回轉(zhuǎn)周期、升距、相對升速等主要特征參數(shù)如表 2所示。

由表 2 可知：潛艇外掛武器模塊前后的空間操縱性能有較大的變化，特別是回轉(zhuǎn)直徑和回轉(zhuǎn)周期，這說明潛艇外掛武器模塊后的空間操縱能力下降了，空間機動能力降低了。

表2 外掛前后空間操縱性能指標(biāo)對照表Tab. 2 Contrast on performance indexes of spatial maneuverability before and after carrying weapons module

4 結(jié) 語

本文結(jié)合 5 種典型運動對外掛武器模塊后的潛艇運動進行了仿真，結(jié)果表明外掛后的潛艇的機動性降低，相關(guān)的操縱性能指標(biāo)下降。仿真結(jié)果符合潛艇本身的運動特性，說明所建立的模型合理。本文的仿真研究為下一步潛艇外掛武器模塊的可行性論證奠定了一定的理論基礎(chǔ)。

[1]范乃忠. 現(xiàn)代潛艇的外掛布雷裝置[J]. 現(xiàn)代艦船, 1995, 6(2): 32–34.

[2]錢東. 美國未來的大型UUV-MANTA[J]. 魚雷技術(shù), 2003(1)41–50. QIAN Dong. The large-scale unmanned undersea vehicles by the US in the future [J]. Torpedl Tochnology, 2003(1) 41–50.

[3]SIRMALIS J E. Pursuing the MANTA vision: recent at-sea technology demonstration results[C]//UDT 2001

[4]楊路春. 潛艇在有外部外掛情況下操縱性水動力導(dǎo)數(shù)的數(shù)值計算方法研究[D]. 哈爾濱：哈爾濱工程大學(xué), 2009.

[5]付延輝, 萬延橙. 潛艇空間運動研究進展[C]//國內(nèi)外潛艇操縱性研究評述, 北京: 國防科工委水動力學(xué)專業(yè)組, 1992: 13–31

[6]夏飛, 陳源. 潛艇操縱性預(yù)報評述[A]. 國內(nèi)外潛艇操縱性研究評述, 北京: 國防科工委水動力學(xué)專業(yè)組, 1992: 1–12

[7]JOHN T H. Development of a six degee of freedom motion Simulation model for use in subma-rine design analysis[R]. ADA078410. USA: DTNSRDC, 1978: 1–20

[8]胡坤, 張洪剛, 徐亦凡. 潛艇水下空間機動仿真與分析[J]. 北京: 計算機仿真, 2006, 23(5): 10–13

[9]MARTLON G, GRANT P H. Standard equations of motion for submarine simulation. AD653861 June 1967, Naval Ship Research and Development Center.

[10]陳厚泰, 徐玉如, 黃錫榮. 關(guān)于深潛艇六自由度運動的流體動力系數(shù)測定和推算[J]. 船工學(xué)報, 1980(1).

[11]李殿璞. 船舶運動與建模(第二版)[M]. 北京: 國防工業(yè)出版社, 2008.

Motion modeling and simulation of submarines carrying external weapons module

LIU Yang1, SONG Xiao-dong2, LI Chun-lai3
(1. Dalian Naval Academy Brigade for Midshipmen, Dalian 116000, China; 2. No. 91421 Unit of PLA, Qingdao 266011, China; 3. Naval University of Engineering, Wuhan 430033, China)

To research the motion performance of submarines carrying external weapons module, this paper uses a way of estimation to calculate the hydrodynamic coefficients of submarines with weapons module, build the Matlab mathematical model of 6D motion, and do comparative simulation analysis of 5 kinds of typical motion of submarines carrying external weapons module before and after. The result of simulation indicates that carrying external weapons module has an effect on a submarine's motion performance to some extent, causing its maneuverability and steering quality to reduce. The simulated model can reflect motion features of submarines with external weapons module, which may provide theoretical foundation for further feasibility analysis.

submarine；external weapons module；simulation；hydrodynamics coefficient estimation

PT391.9

A

1672–7619(2017)05–0053–07

10.3404/j.issn.1672–7619.2017.05.011

2016–12–21

劉洋(1994–)，男，本科生，研究方向為武器發(fā)射與控制。