吉洪蕾，趙輝, 陳仁良，吳文韜

1. 重慶大學(xué) 航空航天學(xué)院，重慶 400030 2. 南京航空航天大學(xué) 航空宇航學(xué)院 直升機(jī)旋翼動(dòng)力學(xué)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室，南京 210016 3.中國(guó)船舶工業(yè)系統(tǒng)工程研究院， 北京 100094

直升機(jī)艦面起降的風(fēng)限圖(Wind-Over-Deck,WOD)是指某一特定機(jī)/艦組合在不同海況、晝夜環(huán)境和不同進(jìn)場(chǎng)方向等條件下，直升機(jī)在其載艦上安全起降的最大相對(duì)風(fēng)速、風(fēng)向限制[1]。風(fēng)限圖是保障直升機(jī)艦面起降安全的重要依據(jù)，在執(zhí)行艦面起降任務(wù)之前駕駛員觀(guān)察艦船相對(duì)風(fēng)向、風(fēng)速是否位于風(fēng)限圖內(nèi)判定該任務(wù)能否完成。如果實(shí)際使用的風(fēng)限圖大于真實(shí)的風(fēng)向風(fēng)速限制，會(huì)危及直升機(jī)艦面起降安全；反之，則會(huì)降低機(jī)/艦動(dòng)態(tài)配合的任務(wù)效能。目前確定風(fēng)限圖的主要手段仍然是直升機(jī)與載艦的動(dòng)態(tài)配合試飛，盡管該方法能夠考慮各方面因素影響、試飛結(jié)果真實(shí)，但是耗資多、周期長(zhǎng)、風(fēng)險(xiǎn)大、受環(huán)境條件限制[2]，迫使各軍事技術(shù)強(qiáng)國(guó)尋找替代方法。采用理論方法計(jì)算風(fēng)限圖，不僅能夠克服機(jī)/艦動(dòng)態(tài)配合試飛的諸多缺點(diǎn)，還能用于直升機(jī)艦面起降安全的影響因素分析，獲得廣泛關(guān)注。

孫文勝[3]采用二維定常CFD方法計(jì)算艦船尾流，較早地開(kāi)展了直升機(jī)艦面起降風(fēng)限圖的理論計(jì)算研究。艦面起降的危險(xiǎn)狀況常常發(fā)生于直升機(jī)懸停在艦面起降點(diǎn)上方跟進(jìn)艦船運(yùn)動(dòng)時(shí)，根據(jù)這一特點(diǎn)，建立了耦合艦船尾流的飛行動(dòng)力學(xué)模型，提出了直升機(jī)艦面起降的最大姿態(tài)、操縱和功率約束等風(fēng)限圖計(jì)算判據(jù)，通過(guò)求解直升機(jī)在起降點(diǎn)上方某一危險(xiǎn)高度艦船尾流中的飛行動(dòng)力學(xué)平衡方程獲得直升機(jī)狀態(tài)和駕駛員操縱量，與風(fēng)限圖計(jì)算判據(jù)對(duì)比得到最大的風(fēng)速風(fēng)向限制，繪制風(fēng)限圖，形成了直升機(jī)穩(wěn)態(tài)配平的風(fēng)限圖計(jì)算方法。孫傳偉等[4]進(jìn)一步探索了艦面上層建筑“陡壁效應(yīng)”引起的艦船尾流三維空間變化對(duì)直升機(jī)平衡的影響，指出飛行甲板上方由于機(jī)庫(kù)導(dǎo)致的垂向氣流及其分布對(duì)艦載直升機(jī)懸停操縱特性影響較大，而側(cè)向分量影響較小。趙唯義和王占勇[5]采用縮比艦船模型的風(fēng)洞試驗(yàn)數(shù)據(jù)分析了艦船尾流下洗、側(cè)沖以及渦流對(duì)直升機(jī)操穩(wěn)性能的作用機(jī)理。李軍亮等[6]研究了側(cè)風(fēng)對(duì)艦載直升機(jī)懸停狀態(tài)和駕駛員操縱的作用，發(fā)現(xiàn)小風(fēng)速時(shí)側(cè)風(fēng)對(duì)直升機(jī)平衡狀態(tài)的作用較小，然而大風(fēng)速時(shí)容易導(dǎo)致尾槳操縱量超過(guò)安全操縱極限。基于大量的直升機(jī)岸基懸停試驗(yàn)數(shù)據(jù)和駕駛員感受描述，Hoencamp和Pavel[7]總結(jié)了直升機(jī)艦面起降風(fēng)限圖的理論計(jì)算判據(jù)。徐廣等[8-9]求解穩(wěn)態(tài)艦船尾流中的直升機(jī)平衡方程，研究了UH-60A直升機(jī)在國(guó)外某航母甲板空氣流場(chǎng)中的平衡與著艦策略，指出駕駛員通過(guò)調(diào)整直升機(jī)航向姿態(tài)保持迎風(fēng)方向進(jìn)場(chǎng)能夠減少氣流的干擾，提升著艦的安全性。

上述風(fēng)限圖的穩(wěn)態(tài)配平計(jì)算方法固然能夠捕捉到起降點(diǎn)上方艦船尾流沿空間非均勻變化的作用，然而直升機(jī)艦面起降是一個(gè)動(dòng)態(tài)過(guò)程，當(dāng)直升機(jī)快速飛過(guò)隨時(shí)間和空間變化的艦船尾流場(chǎng)時(shí)，可能導(dǎo)致飛行狀態(tài)和需要的駕駛員操縱發(fā)生劇烈變化，駕駛員對(duì)快速變化的直升機(jī)狀態(tài)的應(yīng)對(duì)能力也有重要影響，現(xiàn)有方法由于無(wú)法考慮到這些因素的作用，導(dǎo)致風(fēng)限圖計(jì)算的準(zhǔn)確度受到限制。雖然也有研究人員開(kāi)展了艦船非定常尾流的數(shù)值模擬[10-14]和直升機(jī)艦面起降的動(dòng)態(tài)仿真研究[15-22]，但由于集成艦船非定常尾流的直升機(jī)飛行動(dòng)力學(xué)模型計(jì)算效率低，尚未用于風(fēng)限圖的理論計(jì)算。

基于耦合艦船非定常尾流的直升機(jī)飛行動(dòng)力學(xué)高效率計(jì)算模型，本文首先建立適于多軸協(xié)同操控的駕駛員模型，控制直升機(jī)飛行動(dòng)力學(xué)模型完成艦面起降動(dòng)態(tài)仿真，捕捉駕駛員感知、決策延誤和操控頻率限制的作用；然后發(fā)展多種艦面起降軌跡的數(shù)學(xué)描述及生成方法，以考慮不同起降方式導(dǎo)致艦船尾流時(shí)空變化的影響；最后將飛行動(dòng)力學(xué)模型、駕駛員模型和艦面起降軌跡生成算法集成，形成全面計(jì)及艦船非定常尾流環(huán)境、直升機(jī)運(yùn)動(dòng)和駕駛員操控特點(diǎn)綜合影響的直升機(jī)艦面起降動(dòng)態(tài)仿真方法，用于風(fēng)限圖計(jì)算，提升風(fēng)限圖計(jì)算的準(zhǔn)確度。

1 耦合艦船非定常尾流的飛行動(dòng)力學(xué)模型

耦合艦船尾流的飛行動(dòng)力學(xué)模型是直升機(jī)艦面起降動(dòng)態(tài)仿真和風(fēng)限圖計(jì)算的基礎(chǔ)，難點(diǎn)在于隨時(shí)間和空間變化的艦船尾流與直升機(jī)飛行動(dòng)力學(xué)模型的動(dòng)態(tài)耦合建模。采用CFD方法計(jì)算艦面非定常流場(chǎng)，獲得離散數(shù)據(jù)樣本，建立艦船非定常尾流與飛行動(dòng)力學(xué)模型的高效率動(dòng)態(tài)耦合方法，形成耦合艦船非定常尾流的直升機(jī)飛行動(dòng)力學(xué)模型。詳細(xì)的艦船尾流模擬、耦合艦船尾流的飛行動(dòng)力學(xué)建模及驗(yàn)證見(jiàn)文獻(xiàn)[20-21]。以UH-60直升機(jī)為例，最終得到的直升機(jī)艦面起降飛行動(dòng)力學(xué)模型共21個(gè)運(yùn)動(dòng)自由度，具體包括機(jī)體的6個(gè)剛體運(yùn)動(dòng)自由度、旋翼4片槳葉的揮舞和擺振自由度、旋翼槳葉動(dòng)態(tài)扭轉(zhuǎn)的1個(gè)自由度、旋翼動(dòng)態(tài)入流的3個(gè)自由度、尾槳?jiǎng)討B(tài)入流的1個(gè)自由度、機(jī)身的1個(gè)側(cè)洗和下洗自由度。簡(jiǎn)化后的直升機(jī)飛行動(dòng)力學(xué)方程組可以表示為

(1)

2 直升機(jī)艦面起降動(dòng)態(tài)仿真方法

2.1 適用于多軸協(xié)同操控的駕駛員模型

實(shí)際的直升機(jī)艦面起降是一種復(fù)雜的多軸協(xié)同操控任務(wù)，然而艦載直升機(jī)往往具有良好的角速率穩(wěn)定性和姿態(tài)操控性，即使直升機(jī)各軸響應(yīng)間存在交叉耦合，駕駛員也可通過(guò)對(duì)各軸響應(yīng)的補(bǔ)償實(shí)現(xiàn)解耦操控。因此，假設(shè)在艦面起降過(guò)程中駕駛員對(duì)直升機(jī)各軸響應(yīng)的補(bǔ)償操控是獨(dú)立的，各軸操控的駕駛員模型可以單獨(dú)設(shè)計(jì)?；谏鲜黾僭O(shè)，設(shè)計(jì)了適用于多軸協(xié)同操控的駕駛員模型，其結(jié)構(gòu)如圖1所示，x,y,z為直升機(jī)縱向、橫向和垂向的位置坐標(biāo)，vz為直升機(jī)垂向運(yùn)動(dòng)速度，YpΦ、YpΘ、YpΨ分別為滾轉(zhuǎn)、俯仰和偏航姿態(tài)駕駛員控制模型，Ypw為垂向速率駕駛員控制模型，Ypx、Ypy、Ypz分別為縱向、橫向和垂向位置駕駛員控制模型，Φcom、Θcom分別為滾轉(zhuǎn)和俯仰姿態(tài)駕駛員指令，vzcom為垂向速度駕駛員指令，xcom、ycom、zcom為要求駕駛員跟蹤的直升機(jī)縱向、橫向和垂向位置坐標(biāo)，Ψcom為要求駕駛員跟蹤的直升機(jī)偏航角。

駕駛員的縱、橫向操控行為類(lèi)似，模型分別由內(nèi)環(huán)和外環(huán)兩個(gè)回路構(gòu)成。其中，內(nèi)環(huán)為姿態(tài)回路，主要功能是對(duì)直升機(jī)艦面起降過(guò)程的姿態(tài)角偏差進(jìn)行補(bǔ)償。本文采用精確駕駛員模型[23-24]模擬這一過(guò)程，俯仰、滾轉(zhuǎn)姿態(tài)的駕駛員控制模型YpΘ(s)和YpΦ(s)為

(2)

圖1 多軸協(xié)同操控的駕駛員模型Fig.1 Pilot model for multi-axis coordinated controls

式中：KpΘ和KpΦ分別為俯仰、滾轉(zhuǎn)姿態(tài)駕駛員控制模型增益，用來(lái)模擬駕駛員對(duì)控制變量誤差的應(yīng)對(duì)能力。通過(guò)調(diào)整駕駛員控制模型增益KpΘ和KpΦ獲得期望的駕駛員與直升機(jī)組合系統(tǒng)穿越頻率ωcΘ和ωcΦ，使組合系統(tǒng)的幅值YpΘ/ΦYc(s)s=jω在穿越頻率處為1，即

YpΘ/ΦYcjωcΘ/Φ=1

Yc(s)為直升機(jī)飛行動(dòng)力學(xué)模型；TLs+1/TDs+1為駕駛員模型的超前滯后環(huán)節(jié)，TL為超前時(shí)間常數(shù)，反映了駕駛員預(yù)測(cè)操控輸入的能力，TD為滯后時(shí)間常數(shù)，描述了駕駛員產(chǎn)生需要操控輸入的難易程度。根據(jù)駕駛員模型理論[23]，TL和TD的調(diào)整應(yīng)使YpΘ/ΦYc(s)在穿越頻率附近具有以下形式:

式(2)中余下的兩項(xiàng)描述了駕駛員生理特征的限制，其中e-τs是一個(gè)延遲時(shí)間為τ的純時(shí)間延遲環(huán)節(jié)，用來(lái)模擬駕駛員決策過(guò)程導(dǎo)致的時(shí)間延誤，1/(TNs+1)描述了駕駛員操控的肌肉動(dòng)力學(xué)特性。

外環(huán)為軌跡跟蹤控制，主要功能是跟蹤預(yù)定的飛行軌跡完成艦面起降任務(wù)。根據(jù)Hess等[25-26]的研究，縱、橫位置駕駛員控制模型YpX(s)和YpY(s)的選取應(yīng)使直升機(jī)與駕駛員組合系統(tǒng)的縱、橫向位置響應(yīng)在穿越頻率ωcX和ωcY附近具有與K/s類(lèi)似的特性，本研究中YpX(s)和YpY(s)的形式為

(3)

式中：駕駛員控制模型增益系數(shù)KPX和KPY的選取應(yīng)使外環(huán)位置回路穿越頻率為內(nèi)環(huán)姿態(tài)回路穿越頻率的1/3[25]，即

(4)

縱、橫向內(nèi)環(huán)姿態(tài)駕駛員控制模型的穿越頻率由直升機(jī)艦面起降任務(wù)的難易程度決定。預(yù)先對(duì)直升機(jī)艦面起降任務(wù)進(jìn)行動(dòng)態(tài)仿真，獲得起降過(guò)程中內(nèi)環(huán)姿態(tài)響應(yīng)的最大滾轉(zhuǎn)和俯仰角變化ΔΦmax和ΔΘmax，以及最大滾轉(zhuǎn)和俯仰角速率變化Δpmax和Δqmax，俯仰、滾轉(zhuǎn)操控的穿越頻率為[27]

(5)

對(duì)于垂向操控，以垂向速度回路為內(nèi)環(huán)，駕駛員操控模型同樣采用精確駕駛員模型理論建立；垂向位置回路為外環(huán)，駕駛員模型僅由一個(gè)常數(shù)增益項(xiàng)構(gòu)成。航向操控僅由偏航姿態(tài)駕駛員控制模型完成。垂向和航向的駕駛員操控模型設(shè)計(jì)過(guò)程與縱、橫向?qū)?yīng)部分類(lèi)似，這里不再贅述。

2.2 艦面起降軌跡的數(shù)學(xué)描述和生成

艦面起降空間狹小，起降甲板周?chē)嬖谥T多障礙物，且艦船尾部氣流環(huán)境復(fù)雜，根據(jù)起降場(chǎng)地和風(fēng)況的不同應(yīng)當(dāng)因地制宜發(fā)展多種艦面起降方式，增加艦面起降的安全性。目前，國(guó)內(nèi)外常用的直升機(jī)艦面起降方式包括從后向前起降、迎風(fēng)起降、橫向起降、從前向后起降、艦艉起降和斜向起降6種[28]。上述6種起降方式的軌跡都可由穩(wěn)定下滑、橫向平移、航向調(diào)整和起飛爬升這4個(gè)基本飛行軌跡組合而成。

1) 穩(wěn)定下滑

在穩(wěn)定下滑過(guò)程，駕駛員實(shí)施必要的操控使直升機(jī)從遠(yuǎn)方到達(dá)艦船附近一點(diǎn)跟進(jìn)艦船運(yùn)動(dòng)，直升機(jī)初始前飛速度降低至與艦船相同的速度，同時(shí)降低飛行高度，在下滑過(guò)程中還要保持相對(duì)恒定的下滑角。此任務(wù)可分為3個(gè)階段：初始階段，直升機(jī)在一定高度上作水平均速飛行，穩(wěn)定下滑開(kāi)始后，在垂向逐步過(guò)渡到以一定的加速度下降，在水平方向減速逐步達(dá)到設(shè)定的水平速度；穩(wěn)定下滑階段，直升機(jī)在垂向保持恒定的加速度下滑，在水平方向保持恒定的加速度減速；終止階段，在水平和垂直方向的加速度逐漸降為零，垂向速度也逐漸降為零，直升機(jī)下降到一定高度跟進(jìn)艦船運(yùn)動(dòng)，與艦船相對(duì)位置保持不變。

假設(shè)穩(wěn)定下滑初始階段的直升機(jī)水平速度為Vg0，垂向位置為zg0，艦船航行速度為Vs，直升機(jī)和艦船的水平距離為Δxg0，設(shè)定的下滑角為γd，終止階段的直升機(jī)垂向位置為zgd，3個(gè)階段的開(kāi)始時(shí)刻為t0、t1和t2，結(jié)束時(shí)刻為t1、t2和t3，每個(gè)階段所用時(shí)間為Δt1、 Δt2和Δt3，3個(gè)階段所用總時(shí)間為Δtf。沿慣性坐標(biāo)系x軸的水平加速度的數(shù)學(xué)描述為

ax(t)=

(6)

沿慣性坐標(biāo)系z(mì)軸的垂向加速度數(shù)學(xué)描述為

(7)

假設(shè)Δt1=Δt3=ctΔtf,ct為穩(wěn)定下滑初始階段時(shí)間Δt1和終止階段時(shí)間Δt3所占穩(wěn)定下滑總時(shí)間Δtf的比例系數(shù)，式(7)和式(8)均為參數(shù)ax2、az1、az2、ct和Δtf的函數(shù)，通過(guò)穩(wěn)定下滑階段的下滑角約束、終止階段的水平和垂直方向速度和位置約束可求解上述參數(shù)，對(duì)加速度時(shí)間歷程積分得到直升機(jī)穩(wěn)定下滑的速度和位移時(shí)間歷程。

2) 橫向平移

橫向平移開(kāi)始前直升機(jī)位于艦船甲板起降點(diǎn)側(cè)方與艦船相對(duì)位置保持不變。橫向平移加速階段，直升機(jī)逐漸加速，直至達(dá)到設(shè)定的側(cè)向平移速度；然后是均速階段，保持設(shè)定的側(cè)向平移速度，直至接近艦面起降點(diǎn)上方；最后是減速階段，逐步將側(cè)向平移速度降為零，使直升機(jī)最終位于起降點(diǎn)上方跟進(jìn)艦船運(yùn)動(dòng)。

假設(shè)直升機(jī)橫向平移加速、勻速和減速階段的開(kāi)始時(shí)刻為t0、t1和t2，結(jié)束時(shí)刻為t1、t2和t3，3個(gè)階段所用時(shí)間為Δt1、 Δt2和Δt3。橫向平移加速度的數(shù)學(xué)描述為

ay(t)=

(8)

假設(shè)3個(gè)階段所用時(shí)間相同，即Δt1=Δt2=Δt3=Δtf/3，式(8)僅為ay2和Δtf的函數(shù)，給定最大橫向平移加速度ay2和橫向移動(dòng)距離Δy，可以解得上述參數(shù)，對(duì)加速度時(shí)間歷程積分可得直升機(jī)橫向平移的速度和位移時(shí)間歷程。

3) 航向調(diào)整

當(dāng)直升機(jī)位于艦面起降點(diǎn)上方時(shí)往往需要調(diào)整航向姿態(tài)，以適應(yīng)直升機(jī)起飛或者降落時(shí)的風(fēng)向限制。航向調(diào)整的偏航角加速度的數(shù)學(xué)描述為

α(t)=

(9)

假設(shè)Δt1=Δt3=Δt2/2=Δtf/4，給定航向調(diào)整的偏航角變化ΔΨ和最大偏航角加速度αΨ，可以解得Δtf，代入式(9)，積分可得偏航角速率和角度時(shí)間歷程。

4) 起飛爬升

典型的直升機(jī)艦面起飛爬升包括3個(gè)階段：初始階段，駕駛員操控直升機(jī)在垂向加速達(dá)到要求的爬升率，水平方向加速達(dá)到要求的水平加速度；穩(wěn)定爬升階段，直升機(jī)保持恒定的爬升率和水平加速度飛行；終止階段，將爬升率和水平加速度降為零，直升機(jī)在指定高度穩(wěn)定平飛。

假設(shè)直升機(jī)起飛爬升3個(gè)階段的開(kāi)始時(shí)刻為t0、t1和t2，結(jié)束時(shí)刻為t1、t2和t3，所用時(shí)間為Δt1、 Δt2和Δt3。沿慣性坐標(biāo)系x軸的加速度的數(shù)學(xué)描述為

ax(t)=

(10)

沿慣性坐標(biāo)系z(mì)軸的加速度的數(shù)學(xué)描述為

az(t)=

(11)

假設(shè)Δt1=Δt3=ctΔtf，式(10)和式(11)均為參數(shù)ax2、az2、ct和Δtf的函數(shù)。給定允許的最大爬升加速度az2，通過(guò)起飛爬升終止階段的水平速度、垂向速度和高度約束可以解得上述參數(shù)，對(duì)加速度時(shí)間歷程積分可得直升機(jī)起飛爬升的速度和位移時(shí)間歷程。

2.3 直升機(jī)艦面起降的數(shù)值仿真方法

將駕駛員模型和艦面起降軌跡生成算法與耦合艦船非定常尾流的飛行動(dòng)力學(xué)模型集成，構(gòu)成圖2的直升機(jī)艦面起降數(shù)值仿真流程。其中，耦合艦船非定常尾流的直升機(jī)飛行動(dòng)力學(xué)模型能夠反映艦船非定常尾流和直升機(jī)運(yùn)動(dòng)的作用，駕駛員模型能夠模擬駕駛員感知、決策延誤和操控頻率限制的作用，而采用多種艦面起降方式仿真可以研究飛行軌跡導(dǎo)致艦船尾流時(shí)空變化對(duì)直升機(jī)艦面起降安全的影響。本文采用一種變步長(zhǎng)、變階數(shù)的Adams方法[29]對(duì)式(1)的直升機(jī)飛行動(dòng)力學(xué)方程組積分模擬直升機(jī)的操縱響應(yīng)，并采用雙線(xiàn)性變換方法將駕駛員模型與飛行動(dòng)力學(xué)模型耦合求解，實(shí)現(xiàn)直升機(jī)艦面起降動(dòng)態(tài)仿真，輸出直升機(jī)運(yùn)動(dòng)狀態(tài)和駕駛員操縱量，用于風(fēng)限圖計(jì)算。

圖2 直升機(jī)艦面起降數(shù)值仿真流程Fig.2 Numerical simulation procedure of helicopter shipboard flight

3 直升機(jī)艦面起降風(fēng)限圖計(jì)算方法

3.1 風(fēng)限圖計(jì)算判據(jù)

直升機(jī)艦面起降安全的影響因素眾多，涉及艦船尾流環(huán)境、直升機(jī)運(yùn)動(dòng)、駕駛員操控和選擇的起降方式三方面。風(fēng)限圖是針對(duì)特定的機(jī)/艦組合制定的，不同型號(hào)直升機(jī)的運(yùn)動(dòng)特性有差異，而不同型號(hào)艦船由于幾何外形和尺寸的差異，艦船尾流的時(shí)空分布也不同。因此，不同機(jī)/艦組合的動(dòng)態(tài)配合特性不同，風(fēng)限圖也不同。本文以UH-60直升機(jī)和LHA兩棲登陸艦組合為例計(jì)算風(fēng)限圖，其中，UH-60直升機(jī)的飛行動(dòng)力學(xué)建模見(jiàn)文獻(xiàn)[20-21]，LHA兩棲登陸艦?zāi)Ｐ腿鐖D3所示，OSHiSHjSHKSH為艦船坐標(biāo)系，具體的流場(chǎng)計(jì)算及驗(yàn)證見(jiàn)文獻(xiàn)[20-22]。針對(duì)UH-60直升機(jī)的特點(diǎn)和相關(guān)規(guī)定[30]，選定如下約束條件作為該機(jī)/艦組合風(fēng)限圖的計(jì)算判據(jù)。

圖3 艦船模型及起降點(diǎn)位置Fig.3 Ship model and locations of its landing spots

判據(jù)1在任何方向的風(fēng)速超過(guò)22.5 m/s、尾槳正側(cè)風(fēng)大于17.5 m/s時(shí)禁止直升機(jī)艦面起降。

判據(jù)2直升機(jī)總距桿、縱/橫向操縱桿和尾槳腳蹬操縱余量均不小于10%。

判據(jù)3在艦船附近飛行時(shí)，直升機(jī)滾轉(zhuǎn)角不大于8°，俯仰角后倒不超過(guò)7°、前傾不超過(guò)4°。

判據(jù)4發(fā)動(dòng)機(jī)可用功率余量不小于10%。

3.2 風(fēng)限圖計(jì)算方法

根據(jù)3.1節(jié)建立的風(fēng)限圖計(jì)算判據(jù)，首先采用穩(wěn)態(tài)配平方法計(jì)算風(fēng)限圖，然后采用直升機(jī)艦面起飛和降落的動(dòng)態(tài)仿真進(jìn)一步縮小風(fēng)限圖范圍，提升風(fēng)限圖計(jì)算的準(zhǔn)確度。以穩(wěn)態(tài)配平方法的風(fēng)限圖計(jì)算結(jié)果為基礎(chǔ)，基于動(dòng)態(tài)仿真的風(fēng)限圖計(jì)算流程如圖4所示，具體過(guò)程如下：

步驟1從風(fēng)向角θw=-90°、風(fēng)速Uinf=Umθw開(kāi)始，執(zhí)行直升機(jī)艦面起飛和降落的動(dòng)態(tài)仿真，將駕駛員操縱量和直升機(jī)運(yùn)動(dòng)狀態(tài)與風(fēng)限圖計(jì)算判據(jù)對(duì)比，判斷是否超出。其中，θw為風(fēng)向角，左舷為負(fù)，右舷為正；Umθw為穩(wěn)態(tài)配平確定的風(fēng)向角為θw的直升機(jī)懸停風(fēng)速限制。如果是，將風(fēng)速Uinf減小2.5 m/s，重新執(zhí)行直升機(jī)艦面起降動(dòng)態(tài)仿真；如果否，進(jìn)入下一步。

圖4 基于直升機(jī)艦面起降動(dòng)態(tài)仿真的風(fēng)限圖計(jì)算流程Fig.4 WOD envelope calculation procedure based on dynamic simulation of helicopter shipboard operations

步驟2在風(fēng)向角θw的最大風(fēng)速限制為Umθw=Uinf。將風(fēng)向角θw增加15°，Uinf重新設(shè)為Umθw，進(jìn)入下一步。

步驟3判斷風(fēng)向角θw是否大于90°。如果是，終止計(jì)算；如果否，進(jìn)入步驟1繼續(xù)計(jì)算。

4 計(jì)算結(jié)果與分析

首先采用斜向方式執(zhí)行直升機(jī)艦面起飛和降落的動(dòng)態(tài)仿真，并確定風(fēng)限圖。斜向起飛和降落的軌跡由穩(wěn)定下滑、航向調(diào)整和起飛爬升3個(gè)基本飛行軌跡組合而成，如圖5所示。圖6為直升機(jī)從艦面起降點(diǎn)上方斜向起飛的地速時(shí)間歷程，其中，地速坐標(biāo)系的x和z軸與圖3(a)中艦船坐標(biāo)系的x和z軸方向相反，兩個(gè)坐標(biāo)系的y軸指向相同。

圖7～圖9為UH-60直升機(jī)在不同風(fēng)況艦船尾流中斜向起飛的駕駛員操縱桿量、直升機(jī)姿態(tài)角和需用功率Preq時(shí)間歷程。圖7中，δlat為橫向操縱桿量，向右壓桿為正；δlon為縱向操縱桿量，向后拉桿為正；δcol為總距操縱桿量，總距增大桿量為正；δped為尾槳腳蹬操縱量，右踩為正。直升機(jī)起飛的起點(diǎn)位于艦面起降點(diǎn)8上方5 m，總重為8 500 kg。動(dòng)態(tài)仿真的初始狀態(tài)由直升機(jī)在艦船尾流中的配平計(jì)算獲得。從圖中可以看出，艦船尾流對(duì)直升機(jī)艦面起飛操縱量、姿態(tài)和需用功率的影響非常大，風(fēng)況變化使直升機(jī)在艦面起降過(guò)程中經(jīng)歷不同的艦船尾流速度時(shí)空變化，導(dǎo)致直升機(jī)飛行狀態(tài)和需要的駕駛員操縱補(bǔ)償變化劇烈。雖然在圖中所有風(fēng)況艦船尾流環(huán)境中直升機(jī)保持懸停的飛行狀態(tài)和駕駛員操縱量都位于風(fēng)限圖判據(jù)內(nèi)，但當(dāng)風(fēng)向角為左舷和右舷45°風(fēng)況時(shí)，直升機(jī)艦面起降的腳蹬操縱和滾轉(zhuǎn)角分別超出了風(fēng)限圖計(jì)算判據(jù)，表明在此風(fēng)況下直升機(jī)雖然能夠安全地在起降點(diǎn)8上方懸停，但是以斜向方式起飛過(guò)程中受到艦船尾流時(shí)空間變化的強(qiáng)烈干擾，無(wú)法安全地完成起飛任務(wù)。

圖5 斜向起降方式示意圖Fig.5 Diagram of oblique procedure

圖6 斜向起飛地速時(shí)間歷程Fig.6 Time history of ground velocities of oblique takeoff

圖7 斜向起飛駕駛員操縱時(shí)間歷程Fig.7 Time history of pilot controls of oblique takeoff

除了艦面起飛過(guò)程，還分析了直升機(jī)以斜向方式在艦面降落的動(dòng)態(tài)響應(yīng)，得到相似的結(jié)論，即：雖然在某些風(fēng)況下直升機(jī)能夠安全地在艦面起降點(diǎn)上方保持懸停，但在降落過(guò)程中受到艦船尾流時(shí)空變化的強(qiáng)烈干擾，導(dǎo)致無(wú)法安全地完成降落任務(wù)。與起飛過(guò)程相比，在降落過(guò)程中當(dāng)直升機(jī)逐漸接近艦船時(shí)，受到艦船尾流非定常時(shí)空變化的干擾更加強(qiáng)烈，導(dǎo)致直升機(jī)飛行狀態(tài)和需要的駕駛員操縱變化更劇烈，對(duì)直升機(jī)艦面起降的飛行安全和風(fēng)限圖計(jì)算準(zhǔn)確度的影響也更大。因此，風(fēng)限圖的理論計(jì)算應(yīng)同時(shí)考慮直升機(jī)艦面起飛和降落過(guò)程中艦船尾流時(shí)空變化的干擾。

圖8 斜向起飛姿態(tài)角時(shí)間歷程Fig.8 Time history of attitudes of oblique takeoff

圖9 斜向起飛需用功率時(shí)間歷程Fig.9 Time history of required power of oblique takeoff

圖10為采用不同方法計(jì)算UH-60直升機(jī)在LHA艦起降點(diǎn)5、6和8上方艦面起降的風(fēng)限圖，圖中極徑坐標(biāo)為艦船前方來(lái)流風(fēng)速Uinf，極角坐標(biāo)為艦船前方來(lái)流風(fēng)向角θw。穩(wěn)態(tài)配平方法取自文獻(xiàn)[2]，并采用斜向起降方式執(zhí)行直升機(jī)艦面起飛和降落的動(dòng)態(tài)仿真，計(jì)入起飛和降落過(guò)程中艦船尾流的干擾計(jì)算風(fēng)限圖。艦面起降過(guò)程中影響直升機(jī)飛行安全的因素也標(biāo)于圖中，“Col.”表示總距操縱限制，“Ped.”表示腳蹬操縱限制，“Roll”表示滾轉(zhuǎn)角限制。兩種方法均采用本文提出的判據(jù)進(jìn)行風(fēng)限圖計(jì)算。起降點(diǎn)5、6和8的位置見(jiàn)圖3(b)。其中，起降點(diǎn)5位于艦面島型建筑左前方，起降點(diǎn)6位于島型建筑左側(cè)中間位置，起降點(diǎn)8位于島型建筑左后方。

從圖10中可以看出，總距、腳蹬操縱量和滾轉(zhuǎn)角是限制UH-60直升機(jī)在LHA艦安全起降的重要因素。由于艦面上層建筑的作用，不同起降點(diǎn)上方艦船尾流的空間分布存在差異，導(dǎo)致在起降點(diǎn)上方懸停和動(dòng)態(tài)起降的直升機(jī)飛行狀態(tài)和駕駛員操縱量不同，直升機(jī)在3個(gè)起降點(diǎn)的風(fēng)限圖也不同。穩(wěn)態(tài)配平方法和動(dòng)態(tài)仿真方法都能捕捉不同起降點(diǎn)上方艦船尾流的空間分布特征對(duì)風(fēng)限圖理論計(jì)算的影響。

圖10 不同方法計(jì)算的直升機(jī)風(fēng)限圖Fig.10 Predicted WOD envelopes of helicopter by different methods

與穩(wěn)態(tài)配平方法計(jì)算的風(fēng)限圖相比，采用艦面起降動(dòng)態(tài)仿真方法計(jì)算的風(fēng)限圖范圍更小。這是因?yàn)榉€(wěn)態(tài)配平方法只能考慮起降點(diǎn)上方艦船尾流的空間非均勻變化對(duì)直升機(jī)懸停的作用，而直升機(jī)艦面起降動(dòng)態(tài)仿真方法不僅能考慮起降點(diǎn)上方艦船尾流的作用，還能捕捉直升機(jī)艦面起飛和降落機(jī)動(dòng)過(guò)程中經(jīng)歷的艦船尾流時(shí)空變化對(duì)直升機(jī)飛行狀態(tài)和駕駛員操控的影響。從上文分析可知，某些風(fēng)況下直升機(jī)在艦面起降過(guò)程中受到艦船尾流的干擾遠(yuǎn)大于在起降點(diǎn)上方懸停時(shí)受到的作用。因此，直升機(jī)艦面起降動(dòng)態(tài)仿真方法顯著提升了現(xiàn)有風(fēng)限圖計(jì)算方法的準(zhǔn)確度。

為分析不同起降方式對(duì)直升機(jī)艦面起降風(fēng)向、風(fēng)速的限制，除了斜向方式，本文還采用從后向前方式進(jìn)行直升機(jī)艦面起降動(dòng)態(tài)仿真，得到直升機(jī)在起降點(diǎn)8上方艦面起降的風(fēng)限圖，如圖11所示。圖中，“Power”代表直升機(jī)艦面起降過(guò)程的需用功率限制，作為對(duì)比，斜向起降的風(fēng)限圖也繪于圖中。與斜向起降的風(fēng)限圖相比，對(duì)于UH-60直升機(jī)和LHA艦船組合，從總體上看，采用從后向前起降方式能夠擴(kuò)大直升機(jī)艦面起降的風(fēng)速限制。這是因?yàn)閺暮笙蚯捌鸾祵⑽ｋU(xiǎn)的直升機(jī)降落終點(diǎn)和起飛起點(diǎn)移至艦船側(cè)方，有效避開(kāi)了起降點(diǎn)上方劇烈的艦船尾流時(shí)空變化對(duì)直升機(jī)飛行狀態(tài)和駕駛員操控的影響。然而，當(dāng)艦面來(lái)流的風(fēng)向角為右舷30°時(shí)，采用從后向前方式執(zhí)行艦面起降任務(wù)容易導(dǎo)致直升機(jī)發(fā)動(dòng)機(jī)的可用功率不足，這是因?yàn)楫?dāng)風(fēng)向角為右舷30°風(fēng)況時(shí)，直升機(jī)從后向前方式艦面起降的起點(diǎn)和終點(diǎn)恰巧位于島型建筑的下洗氣流內(nèi)。因此，采用不同起降方式完成艦面起降使直升機(jī)經(jīng)歷不同的艦船尾流時(shí)空變化，各有利弊。對(duì)于UH-60直升機(jī)和LHA艦組合，與斜向方式相比，采用從后向前方式執(zhí)行直升機(jī)艦面起降的風(fēng)速限制范圍更大，能夠提升機(jī)/艦動(dòng)態(tài)配合的任務(wù)效能。采用從后向前方式進(jìn)行直升機(jī)艦面起降動(dòng)態(tài)仿真時(shí)，直升機(jī)總重為8 500 kg，從后向前起飛和降落的軌跡由穩(wěn)定下滑、橫向平移和爬升飛行3個(gè)基本飛行軌跡組合而成，如圖12所示。

圖11 直升機(jī)以不同起降方式在起降點(diǎn)8的風(fēng)限圖Fig.11 Predicted WOD envelopes of helicopter in landing spot 8 with different procedures

圖12 從后向前起降方式示意圖Fig.12 Diagram of fore/aft procedure

5 結(jié) 論

本文基于耦合艦船非定常尾流的飛行動(dòng)力學(xué)模型，發(fā)展了適于多軸協(xié)同操控的駕駛員模型，建立了艦面起降軌跡的數(shù)學(xué)描述與生成方法，最終形成了直升機(jī)艦面起降的動(dòng)態(tài)仿真方法。在此基礎(chǔ)上，建立了基于直升機(jī)艦面起降動(dòng)態(tài)仿真的風(fēng)限圖計(jì)算方法。

1) 常用直升機(jī)艦面起降方式的軌跡可由穩(wěn)定下滑、橫向平移、航向調(diào)整和起飛爬升這4種基本飛行軌跡組合而成。本文的艦面起降軌跡數(shù)學(xué)描述與生成方法能夠給出光滑的飛行軌跡。建立的多軸協(xié)同操控駕駛員模型能夠控制飛行動(dòng)力學(xué)模型跟蹤生成的飛行軌跡實(shí)現(xiàn)直升機(jī)艦面起降動(dòng)態(tài)仿真。

2) 某些風(fēng)況下直升機(jī)在艦面起降過(guò)程中受到艦船尾流的干擾遠(yuǎn)大于在起降點(diǎn)上方保持平衡時(shí)受到的作用。基于直升機(jī)艦面起降動(dòng)態(tài)仿真的風(fēng)限圖計(jì)算方法不僅能考慮起降點(diǎn)上方艦船尾流沿空間非線(xiàn)性變化對(duì)直升機(jī)平衡狀態(tài)的作用，還能捕捉到不同起降方式導(dǎo)致艦船尾流時(shí)空變化對(duì)直升機(jī)飛行狀態(tài)和駕駛員操控的干擾，顯著提升了風(fēng)限圖計(jì)算的準(zhǔn)確度。

3) 采用不同起降方式在不同艦面起降點(diǎn)執(zhí)行起飛和降落任務(wù)使直升機(jī)經(jīng)歷不同的艦船尾流時(shí)空變化，直升機(jī)艦面起降的風(fēng)向、風(fēng)速限制也不同。對(duì)于UH-60直升機(jī)和LHA艦組合，與斜向起降方式相比，從后向前方式使直升機(jī)避開(kāi)在起降點(diǎn)上方時(shí)空變化劇烈的艦船尾流中執(zhí)行起飛和降落任務(wù)，風(fēng)速限制范圍更大，能夠提升直升機(jī)艦面起降的任務(wù)效能。