胡 坤，張 平，黃海峰，孫 奎

(海軍潛艇學院，山東 青島 266199)

1 引言

潛艇水下懸停操縱是潛艇停車后，通過均衡手段使?jié)撏腋∮谒履骋簧疃鹊牟倏v方法。水下懸停是一種以戰(zhàn)術為目的的操縱控制方法，它與潛坐海底、潛坐液體海底和水下錨泊并稱為潛艇水下停泊方式之一[1]。懸停的特點是：潛艇可在其工作深度內的任意深度上實施停泊，并能從停泊狀態(tài)立即轉入戰(zhàn)術機動。潛艇水下懸?？梢韵菪龢肼暫退畡恿υ肼?，降低機械振動噪聲，提高聲納作用距離，節(jié)省能源，對提高潛艇的戰(zhàn)斗力和具有重大的現(xiàn)實意義和戰(zhàn)術價值。應東對潛艇水下懸停運動進行了建模分析，研究了不同因素對懸停操縱的影響[2]；俞科云等分析了懸停運動的特征，并對潛艇水下懸停運動的控制進行了仿真研究[3]；郝英澤等分析了海水密度變化對潛艇懸停垂向運動的影響，得出了穩(wěn)定深度與海水密度變化率的量的關系[4]；熊瑛等針對潛艇水下懸?？刂圃O計了一種模糊控制器，通過仿真表明模糊控制器具有超調量小、調整時間段及魯棒性好的特點[5]。本文通過分析影響潛艇懸停的密度、溫度、鹽度、海流等相關海洋環(huán)境因素，以模型潛艇為研究對象，在潛艇懸停運動數(shù)學模型的基礎上進行仿真計算，研究懸停海域不同密度梯度類型和不同密度梯度值對潛艇懸停操縱的影響，在仿真結果的基礎上，分析了潛艇懸停操縱的海洋環(huán)境因素和懸停穩(wěn)定性之間的邏輯關系。

2 影響潛艇懸停的海洋環(huán)境因素

2.1 海水密度變化對潛艇浮力的影響

海水密度變化是引起潛艇在懸停過程中浮力改變的最根本原因之一。潛艇在海水中懸停時，其受力是平衡的，即潛艇的重力等于艇體所排海水產生的浮力，海水密度的改變將導致潛艇浮力的變化，進而產生浮力差，破壞潛艇懸停的平衡受力狀態(tài)，最終導致潛艇在浮力差的作用下產生垂直面內的加速度和速度，使?jié)撏x初始懸停深度。

假設潛艇在密度為ρ1的海域均衡好以后，過渡到密度為ρ2的海域，引起的浮力變化量可用下式確定：

ΔB1=(ρ2-ρ1)V

(1)

式中，V為潛艇水下容積排水量；ρ1為原海區(qū)的海水密度；ρ2為新海區(qū)的海水密度。

海水密度的變化主要是由海水溫度、鹽度和靜壓力變化引起的，在水深小于75m的大陸沿海，自水面至40m深度區(qū)間，海水的溫、鹽度對海水密度的影響十分突出，在7～10月份，還經常形成密度躍變層。在臺灣東、南部海域，在50～300m的深度區(qū)間，靜壓力和鹽度的變化對海水密度的影響起主要作用。

海水密度在垂直面的變化梯度通常有三種典型的形式：第一種是海水密度隨深度增加而增大，即密度呈正梯度變化，如圖1(a)所示；第二種是海水密度不隨深度變化，即密度均勻層，如圖1(b)所示；第三種是海水密度隨深度的增加而減小，即密度呈負梯度變化，如圖1(c)所示。

圖1 海水密度變化的三種典型梯度類型

2.2 海水溫度變化對潛艇浮力的影響

溫度是引起海水密度變化的主要因素，根據(jù)熱力學原理可知，溫度增高，海水密度減小，溫度降低，海水密度增大。海水溫度除了受季節(jié)和海洋曖流影響外，在水深30～50m以上的海洋表面還受風浪和日照的影響。

實驗表明，海水溫度和海水密度之間并不是簡單的線型關系，當海水溫度范圍不同時，其變化規(guī)律也不盡相同。

當t=20℃～10℃范圍時，海水密度增加率?=0.014%；當t=10℃～4℃范圍時，海水密度增加率?=0.003%。

由此引起的浮力變化量為

ΔB2=?j(t1-t2)ρV

(2)

式中，ρ為原海域海水密度；V為潛艇水下容積排水量；t1為原海域水溫；t2為新海域水溫。

海水的溫度變化對海水密度的影響十分突出，海水溫度的變化引起海水密度的改變，海水密度的改變又會引起浮力變化，導致浮力差的產生，從而影響潛艇懸停穩(wěn)定。

海水溫度的變化引起潛艇排水容積的改變由實驗資料可知，水溫降低1℃，艇殼收縮率為：αt=-4×10-5，由此引起的浮力變化量可表示成

ΔB3=αt(t1-t2)ρV

(3)

2.3 海水鹽度變化對潛艇浮力的影響

海水鹽度對潛艇浮力的影響是通過海水密度表現(xiàn)出來的。海水鹽度是溫度和海水靜壓力的函數(shù)，當海水鹽度改變時，引起海水密度相應的變化，進而引起浮力變化量ΔB4。根據(jù)經驗，在特定的溫度和壓力條件下，當海水鹽度增加1‰時，海水密度增加約1.34‰。

2.4 海水壓力變化對潛艇浮力的影響

潛深增加，水壓力也增大，也會引起海水密度的增加，并使艇殼容積壓縮，導致潛艇浮力改變。

海水壓力變化首先會導致海水密度的改變，雖然通常認為海水是不可壓縮的，但實際上每增加一個大氣壓(對應的深度增加約10米)，海水密度將增加αρ1=5×10-6倍。當潛深由H1米增加到H2米，引起的浮力變化量為ΔB5，并表示成

ΔB5=αρ1(H2-H1)ρV

(4)

式中，H1為原深度；H2為新深度。

海水壓力的變化更是引起潛艇排水容積的變化的重要原因，由實驗資料得知，海水壓力每增加一個大氣壓引起潛艇耐壓殼體容積壓縮率為αρ2=2.0～2.5×10-4。當潛深由H1米增加到H2米，由此引起的浮力變化量可以表示成

ΔB6=αρ2(H2-H1)ρV

(5)

2.5 海流對潛艇懸停的影響

海水受潮汐、海洋地藐、地球自轉、表層風生流、海水與空氣和海底間摩擦阻力的影響，流場在不同深度是存在一定差異的。如潛艇在有起伏較大或有暗礁分布的海區(qū)近海底懸停時，懸停處的流場會因海底不規(guī)則的地藐阻擋、折射而發(fā)生紊亂，紊亂的流場，使艇體產生的動水力和力矩呈不規(guī)則的變化，這種不規(guī)則變化的動水力和力矩不能用均衡的方法所消除，潛艇在不穩(wěn)定的流場中懸停，其深度難以保持穩(wěn)定。

圖2中a為定海流，懸停時，潛艇伴隨定海流一起向下流方向漂移，且對海底產生相對運動。位于潛艇上、下方的流b、c，因風生流和海水對空氣及海底的粘滯阻力作用，大小不盡相同，其產生的動水力，被在懸停初期均衡時所消除，潛艇仍能處于受力平衡狀態(tài)，保持懸停深度不變。當潛艇隨海流漂移至凸起海底上部時，根據(jù)流體運動特性，流c大于流b，從而在艇體上、下部產生壓力差，導致潛艇因原受力平衡狀態(tài)被破壞而使?jié)撏彝Ｉ疃劝l(fā)生改變。因此，為使?jié)撏彝Ｉ疃确€(wěn)定，應避免海面風生流和海底地藐變化產生非均勻流場對潛艇的影響。從保持穩(wěn)定懸停的角度，懸停海區(qū)應選擇在海洋寬闊，海底地藐平坦的海區(qū)進行。而且懸停深度要與海底和海面保持恰當?shù)木嚯x，以減弱流場變化對懸停深度的影響。

圖2 潛艇在不規(guī)則流中的懸停

3 懸停操縱運動控制數(shù)學模型

3.1 懸停運動基本數(shù)學模型

潛艇水下懸停時，先降到經航工況，經準確均衡后再停車。根據(jù)潛艇垂直面運動規(guī)律，水下懸停類似于潛艇在垂直面的慣性潛浮運動，因此，由潛艇垂直面操縱運動非線性方程式[6][7]，經過合理簡化后，潛艇水下懸停運動的基本數(shù)學模型為

(6)

3.2 海洋環(huán)境干擾力模型

為方便研究，對干擾力模型進行了一定的簡化是必要的。由于海水的溫度、鹽度對潛艇懸停的影響是通過海水密度的變化反映出來的，同時海水壓力對潛艇懸停的影響是通過艇體壓縮反映出來的，本文將干擾力模型統(tǒng)一為海水密度變化模型和艇體壓縮模型。

1)艇體壓縮模型

艇體壓縮量隨深度增加而變化設為0.029t/m，艇體壓縮產生的力按下式計算

P3=0.029(ζ-ζ0)

(7)

式中：ζ為當前深度；ζ0為初始深度。

2)海水密度變化模型

海水密度變化產生的力按下計算

P4=-V(ρ-ρ0)

(8)

式中：V為潛艇水下全排水量；ρ為潛艇初始深度海水密度；ρ0為潛艇當前深度海水密度。

均勻層情況下取海水密度為一個定值，即認為海水密度不變化，ρ=ρ0=1.025t/m3。

正梯度情況密度變化公式為

(9)

負梯度情況，密度變化公式為

ρ(ζ)=1.025-0.00001(ζ-ζ0)

(10)

(4)式、(5)式中：ρ(ζ)為當前深度下的海水密度；ζ為當前深度；ζ0為初始深度。

3.3 懸停水艙注排水控制模型

潛艇懸停時所用到的控制執(zhí)行機構主要是懸停專用水艙。懸停水艙實際水量容積Q為考慮相應控制規(guī)律、水艙狀態(tài)初始條件Q0和最大可能的注排水量Qmax后的水艙狀態(tài)方程的積分結果。

懸停專用水艙狀態(tài)方程形式如下

(11)

懸停水艙水量Q1滿足如下條件

(12)

(13)

3.4 懸停水艙閥門開啟模型

在初始狀態(tài)水艙閥門是全關的，閥門開啟系數(shù)Ko=0，在接收到注排水控制信號后閥門在響應時間內按線性過程由全關變化到全開。

注排水時，水艙閥門是全開的，閥門開啟系數(shù)Ko=1，在接收到停止注排水或調水控制信號后閥門在響應時間內按線性過程由全開變化到全關。

3.5 模型有效性分析

為了驗證懸停操縱運動控制數(shù)學模型的有效性，對負梯度條件下的潛艇懸停進行了仿真，仿真條件盡可能與模型潛艇操縱試驗時的海洋環(huán)境接近。具體的仿真環(huán)境為：海水密度梯度類型為弱負梯度，梯度值-0.000015，懸停穩(wěn)定深度40m，潛艇初始不均衡量0.5t，懸停水艙注排水速率10L/s。

表1 數(shù)學模型有效性比對

分析表1數(shù)據(jù)可知，在相似海洋環(huán)境條件下，懸停仿真值與試驗值非常接近，深度偏差、縱傾偏差和單位時間內的注排水量的誤差均在5%以內，特別是單位時間的注排水次數(shù)幾乎相等?？梢钥闯?，本文采用的潛艇懸停運動數(shù)學模型能夠較真實的反映潛艇的實際懸停操縱過程，從而驗證了潛艇懸停運動數(shù)學模型的有效性。

4 數(shù)值仿真分析

前面分析可知，海洋環(huán)境對潛艇懸停的影響是非常復雜的。由于海水的溫度、鹽度、壓力對潛艇懸停的影響是通過海水密度的變化最終反映出來的，因此分析海洋環(huán)境因素對懸停的影響可以落腳于懸停海域的密度梯度類型和不同梯度類型下的梯度值對懸停操縱的影響。

基于懸停運動基本數(shù)學模型、干擾力模型和懸停水艙注排水控制模型，以模型潛艇為研究對象，采用C#語言編寫了潛艇水下懸停操縱運動仿真軟件[9][10]，通過懸停仿真進一步分析密度梯度類型和不同梯度類型下的梯度值對潛艇懸停操縱的影響。

4.1 密度梯度類型對懸停的影響

密度梯度值決定海水密度的變化方式，海水密度變化是影響潛艇懸停穩(wěn)定的最主要因素之一，因此明確懸停海區(qū)密度梯度的情況對潛艇懸停系統(tǒng)的使用至關重要。圖3～圖5分別為海水密度梯度類型為負梯度(梯度值-0.00001)、均勻層和負梯度(梯度值0.00001)時的懸停仿真圖。仿真圖中的四條曲線分別是深度、主動控制力，、累計注排水量和累計注排水次數(shù)與時間的關系。仿真條件設置：懸停水艙注排水速率10L/s，目標懸停深度100m，懸停水艙閥門開啟時間6s。

圖3 負梯度條件下的懸停仿真

圖4 均勻層條件下的懸停仿真

圖5 正梯度條件下的懸停仿真

表2為不同密度梯度類型對潛艇懸停影響的仿真結果。需要特別指出的是，表1及后續(xù)表格中的平均注排水量表示單位時間內的平均注排水量，平均注排水時間表示單位時間內的平均注排水時間，平均注排水次數(shù)表示單位時間內的平均注排水次數(shù)，挽回深度指的是當潛艇存在初始不均衡量時，通過懸停水艙的一次注排水，首次將潛艇深度挽回時，這時的潛艇深度與初始時刻時的潛艇深度之間的差值。

從仿真數(shù)據(jù)可以看出，海洋環(huán)境對潛艇懸停的影響是非常顯著的。海水密度類型為正梯度時，潛艇的懸停要相對于海水密度類型為負梯度時容易的多。正梯度條件下的單位時間內平均注排水次數(shù)和平均注排水量也要遠小于負梯度條件。因此，在選擇懸停深度和海區(qū)方面，從便于操縱控制的角度出發(fā)，應該盡量選擇在海水密度呈正梯度變化的海域進行懸停。

若潛艇懸停海區(qū)的密度梯度類型為正梯度，且密度梯度值大于臨界梯度值時，懸停操縱最為簡單，因為這時潛艇處于穩(wěn)定平衡狀態(tài)，基本上只需經過一到兩次注排水消除潛艇的浮力差后(見表1數(shù)據(jù))，僅僅依靠潛艇所處的正梯度海況條件，就能將潛艇的深度控制住。若潛艇懸停海區(qū)的密度梯度類型為負梯度，且密度梯度值較大時，可能導致潛艇懸停操縱失敗。

4.2 負梯度條件下密度梯度值對懸停的影響

從前面的仿真分析可以看出，在負梯度條件下，潛艇進行懸停操縱時最為困難的。本節(jié)通過仿真分析負梯度條件下不同密度梯度值對懸停操縱的影響。

圖6～圖8分別為海水密度梯度值為-4×10-6、-2×10-5、-8×10-5時的懸停仿真圖。仿真條件設置：海水密度梯度類型為負梯度，懸停穩(wěn)定深度100m，注排水速率10L/s，懸停水艙閥門延遲時間6s，初始不均衡量0.5t。

圖7 密度梯度值為-2×10-5時懸停仿真

圖8 密度梯度值為-8×10-5時懸停仿真

表3為負梯度條件下不同密度梯度值對潛艇懸停的影響仿真數(shù)據(jù)。

表3 負梯度條件下密度梯度值對潛艇懸停影響

通過仿真數(shù)據(jù)可以看出：

1)負梯度條件下，潛艇懸停受海水密度梯度值的變化影響較大。挽回深度、單位時間內平均注排水量、平均注排水時間和平均注排水次數(shù)隨密度梯度值絕對值的增大迅速增大。

2)負梯度條件下，當梯度值大到一定程度時，對控制的要求越高，如果不能夠及時地注排水，即便是一點小小的誤差，也可能造成潛艇深度無法控制從而導致懸停操縱的失敗。

4.3 正梯度條件下密度梯度值對懸停的影響

正梯度條件下，潛艇進行懸停操縱時相對較為容易。本節(jié)通過仿真分析正梯度條件下不同密度梯度值對懸停操縱的影響。圖9～11分別為海水密度梯度值為10-6、9.6×10-6(臨界梯度)、2×10-5時的懸停仿真圖。仿真條件設置：海水密度梯度類型為正梯度，目標懸停深度100m，注排水速率10L/s，初始不均衡量0.5t。

圖9 密度梯度值為10-6時的懸停仿真

圖10 密度梯度值為臨界值的懸停仿真

圖11 密度梯度值為2×10-5時的懸停仿真

表4為海水密度梯度為正梯度條件下密度梯度值對潛艇懸停影響的仿真結果。

表4 正梯度條件下密度梯度值對潛艇懸停影響

通過仿真數(shù)據(jù)可以看出：

1)正梯度條件下，潛艇的懸停操縱要比均勻層和負梯度容易的多，挽回深度、單位時間內平均注排水量、平均注排水時間和平均注排水次數(shù)均較小。只要操縱得當，潛艇可以穩(wěn)定的懸停在目標深度。

2)正梯度條件下，會出現(xiàn)臨界情況，即潛艇的艇體造成的壓縮量所形成的負浮力，和潛艇所處的正梯度所造成的浮力相互抵消，這時潛艇處于隨遇狀態(tài)。經過計算和仿真驗證，臨界梯度值等于9.6×10-6。當梯度值小于臨界梯度值時，相當于潛艇仍然處于不穩(wěn)定平衡狀態(tài)，因此這時的單位時間內平均注排水次數(shù)仍較大；當梯度值超過臨界梯度值以后，潛艇則處于穩(wěn)定平衡狀態(tài)，這時一般只需進行兩次注排水就已足夠將潛艇穩(wěn)定在目標深度上。

5 結論

本文通過對不同海洋環(huán)境條件下潛艇懸停操縱進行了仿真，分析了密度梯度類型和密度值對潛艇懸停的影響。研究表明，海洋環(huán)境對潛艇能否成功進行懸停操縱至關重要。不同的密度梯度類型對潛艇懸停影響非常顯著，海水密度類型為正梯度時，潛艇的懸停要相對于海水密度類型為負梯度時容易的多。正梯度條件下的單位時間內平均注排水次數(shù)和平均注排水量要遠小于負梯度條件下的值。因此，在選擇懸停深度和海區(qū)方面，從便于操縱控制的角度出發(fā)，應該盡量選擇在海水密度呈正梯度變化的海域進行懸停。

需要特別指出的是，懸停系統(tǒng)并不是在任意海洋環(huán)境下都能夠將潛艇成功懸停住，當潛艇處于海水密度為負梯度條件下，且密度梯度值較大時，潛艇深度的微小變化都會產生一個較大的力使得潛艇偏離原來的深度，這時單單利用懸停系統(tǒng)已經很難對潛艇進行懸停控制，必須同時啟動全艇的均衡系統(tǒng)或者使用動力系統(tǒng)和操舵系統(tǒng)控制，否則潛艇會面臨失控的危險。