王文濤，查建平，許寧

中國直升機設(shè)計研究所，江西 景德鎮(zhèn) 333001

隨著社會的不斷發(fā)展，軍用/民用直升機的需求日益增加。使用環(huán)境及任務(wù)需求的多樣化，使得直升機日益復(fù)雜化，機載任務(wù)設(shè)備更加多元化，大尺寸的外載外掛式設(shè)備（如桅桿式雷達、外掛吊艙等）也逐漸增多。目前，國內(nèi)在裝備的研制中更多地關(guān)注裝備的整體性能、可靠性等，而對外掛設(shè)備外形導(dǎo)致的整機氣動特性關(guān)注不足，可能產(chǎn)生額外的氣動力作用到機體上；同時，由于分析和驗證手段的局限性，當(dāng)直升機上出現(xiàn)低頻的振動問題時很難進行分析和處理。NH90、EH101、AH-64D“長弓阿帕奇”“科曼奇”等直升機在試飛初期均出現(xiàn)過“Tail shake/Buffet”等類似的低頻振動現(xiàn)象[1，2]，國內(nèi)在直升機載荷測試、旋翼頂部加裝大尺寸設(shè)備等項目的試飛中也多次出現(xiàn)過低頻振動問題。

早期低頻振動在對振動測試數(shù)據(jù)的常規(guī)分析中無法發(fā)現(xiàn)，僅靠飛行員的感受獲知存在低頻振動問題。飛行員的描述往往存在較大的主觀性，其感受會受到個體敏感性、試飛科目、直升機狀態(tài)等因素的影響產(chǎn)生較大的差異，同時飛行員在試飛時處于多振源的復(fù)雜振動和噪聲環(huán)境中，因此僅靠飛行員描述對低頻振動問題進行定位和分析無法給出解決問題的有效措施。

飛行振動測試是分析直升機振動的有效手段，直升機的振動特性[3，4]決定了振動測試數(shù)據(jù)中包含了主旋翼1Ω、kΩ（k為槳葉片數(shù)）等周期激勵成分以及復(fù)雜的隨機成分，數(shù)據(jù)分析發(fā)現(xiàn)這些成分對應(yīng)低頻振動的幅值相對較小，且沒有明顯的規(guī)律性，目前國內(nèi)在直升機試飛中還沒有有效的手段分析該問題。以傅里葉變換為代表的譜分析方法以周期信號為分析對象，而對非平穩(wěn)信號將失去原有的意義，而希爾伯特變換（HT）是處理非平穩(wěn)信號的有效方法，可用于對調(diào)制信號的解調(diào)[5]，其衍生方法被廣泛應(yīng)用到小突變信號檢測[6]、旋轉(zhuǎn)機械信號分析[7]、軸承、齒輪等旋轉(zhuǎn)部件的故障診斷[8，9]等。

本文通過對直升機中的低頻振動現(xiàn)象及機理分析，基于希爾伯特變換，形成了直升機低頻振動現(xiàn)象的分析方法，建立了飛行員描述與振動數(shù)據(jù)之間的對應(yīng)關(guān)系，為低頻振動問題的分析和處理提供了有效的客觀評價手段。本文建立的方法已成功用于多個直升機型號的低頻振動問題分析，取得了很好的效果。

1 現(xiàn)象及機理

直升機的振源一般包含主旋翼、尾槳、發(fā)動機及傳動系統(tǒng)等旋轉(zhuǎn)部件引起的周期激勵和氣動環(huán)境中的隨機激勵。對飛行員及乘員而言，主旋翼轉(zhuǎn)速頻率1Ω、一階通過頻率kΩ（k為槳葉片數(shù)）及其倍頻NkΩ為其所感受到的主要振動源。

直升機的低頻振動不同于上述常規(guī)的振動源，表現(xiàn)為“抖動”“晃動”“篩動”等現(xiàn)象，類似于坐船或騎馬等，方向多為橫側(cè)向或垂向，有能量逐步積累然后釋放的過程，無明顯的規(guī)律性，時間間隔為1～5s。低頻振動一般出現(xiàn)在某個或某些飛行狀態(tài)或飛行姿態(tài)，與一個或多個因素相關(guān)聯(lián)。因低頻振動是無規(guī)律性的，“低頻”僅指單位時間內(nèi)發(fā)生的次數(shù)較少。

參考文獻[9]指出直升機存在一種特殊的振動問題，其來源于主槳轂等區(qū)域的脫落渦對尾部結(jié)構(gòu)的氣動影響，如圖1所示。影響因素有主槳轂、發(fā)動機進氣及排氣、整流罩等位置的氣動外形以及旋翼到機身的距離等，表現(xiàn)為機體側(cè)向一階彎曲模態(tài)下的響應(yīng)。而機體側(cè)向一階彎曲模態(tài)頻率一般大于5Hz，這時需要與時間間隔1～5s建立關(guān)聯(lián)，為分析低頻振動問題的機理，引入經(jīng)典的“拍現(xiàn)象”[10]。兩個頻率相近、幅值相同的簡諧振動用式（1）來描述：

其中：a、ω和φ分別為簡諧振動的幅值、頻率和相位，其合成振動可以寫成調(diào)制波形的形式：

其中：

即合成振動是頻率為(ω2+ω1)/2的變幅振動，其幅值為以頻率(ω2-ω1)/2變化，是典型的幅值調(diào)制過程，實際呈現(xiàn)出的是幅值調(diào)制函數(shù)(t)的特性，即幅值包絡(luò)線，如圖2所示。

考慮任意非周期函數(shù)a(t)對頻率為ω0的簡諧振動調(diào)制的情況，即：

可描述典型的低頻振動現(xiàn)象，其中被調(diào)制信號sin(ω0+φ)的頻率ω0為機體側(cè)向和垂向一階彎曲固有頻率。

因此，直升機中的低頻振動為非周期調(diào)制函數(shù)對周期振動的調(diào)制，周期振動的頻率為機體低階固有頻率，飛行員所感受的為非周期調(diào)制函數(shù)的變化規(guī)律，即周期函數(shù)的幅值包絡(luò)。低頻指單位時間內(nèi)超出某一閾值的次數(shù)相對較少，時間間隔相對較長，如圖3中出現(xiàn)幅值超出0.1g的時間間隔為2～3s。

圖1 脫落渦對尾部干擾示意Fig.1 Turbulent wake effects on tail

圖2 兩個同幅簡諧振動合成的拍Fig.2 Beats with two simple harmonic vibrations

圖3 典型低頻振動信號Fig.3 Typical low frequency vibration signal

2 分析方法

2.1 希爾伯特變換

設(shè)實函數(shù) x(t)，t(-∞，∞)，其希爾伯特變換[11，12]為函數(shù)x(t)與 1/πt的卷積，即：

由傅里葉變換的性質(zhì)可知，函數(shù)h(t)的傅里葉變換為函數(shù)x(t)的傅里葉變換X( f )與1/πt的傅里葉變換F（1/πt）的乘積，即：

其中

引入解析函數(shù)：

即：

式中：A(t)為希爾伯特變換的包絡(luò)

φ(t)為瞬時相位瞬時頻率定義為

考慮式（4）所示的調(diào)制波形函數(shù)，根據(jù)希爾伯特變換的性質(zhì)有：

解析函數(shù)：

其包絡(luò)為：

因此，通過希爾伯特變換，可以得到低頻振動的幅值包絡(luò)。

2.2 參數(shù)ü幅值圖

直升機低頻振動現(xiàn)象由飛行員提出，一般為直升機在某個特定的飛行狀態(tài)或飛行姿態(tài)下出現(xiàn)，如某個速度下平飛、某個姿態(tài)爬升等，其感受為主觀的，隨飛行員的狀態(tài)、關(guān)注點不同而有所差異。

為相對客觀地反映低頻振動與飛行狀態(tài)、直升機姿態(tài)等參數(shù)的關(guān)系，建立參數(shù)p與振動幅值A(chǔ)之間的關(guān)系，形成參數(shù)ü幅值（püA）圖，其中參數(shù)p可以為飛行狀態(tài)參數(shù)如速度、壓力高度、升降率等，可以為俯仰角、滾轉(zhuǎn)角、側(cè)滑角等直升機姿態(tài)參數(shù)，亦可為總距、周期變距等操縱量，同時根據(jù)需要可以為幾個量的組合；振動幅值A(chǔ)為利用希爾伯特變換得到的幅值包絡(luò)。通過püA圖能夠直觀地看出振動幅值A(chǔ)隨參數(shù)p的變化關(guān)系。

3 實例驗證

為證明本文方法的適用性和有效性，應(yīng)用本文方法對某型機試飛中的橫側(cè)向抖動問題進行分析。某型機調(diào)整試飛過程中，機組人員反映在速度170～220km/h平飛狀態(tài)下機體存在較為明顯的抖動現(xiàn)象，方向為橫側(cè)向，在1～2s出現(xiàn)一次，并存在類似能量累積然后釋放的過程。利用本文方法分析，確定該低頻振動現(xiàn)象為槳轂整流罩后部無規(guī)律的脫落渦引起了機體側(cè)向一階彎曲模態(tài)下的響應(yīng)，是無規(guī)律激勵對周期振動的調(diào)制，典型時間段低頻振動數(shù)據(jù)如圖4所示，利用希爾伯特變換獲得振動幅值的包絡(luò)如圖5所示。對槳轂整流罩進行優(yōu)化設(shè)計后，該低頻振動現(xiàn)象得到明顯改善。

利用本文方法對整個架次數(shù)據(jù)進行處理，以飛行速度作為參數(shù)p，得到速度與振動幅值（püA）的關(guān)系圖，如圖6所示。圖中顯示在速度170～220km/h狀態(tài)，振動峰值的分布明顯大于其他速度狀態(tài)，與飛行員反映的規(guī)律一致。更改后的püA的關(guān)系如圖7所示，圖中顯示在速度170～220 km/h狀態(tài)，振動幅值明顯小于初始狀態(tài)，與飛行員所述“低頻振動現(xiàn)象明顯改善”相吻合。

由此可以證明，基于希爾伯特變換得到振動幅值的包絡(luò)，然后與參數(shù)p形成的püA圖，能夠直觀地呈現(xiàn)飛行員所反映的低頻振動現(xiàn)象。püA圖建立了振動測試數(shù)據(jù)與飛行員主觀感受的對應(yīng)關(guān)系，能夠客觀地描述直升機的低頻振動現(xiàn)象，可以作為分析、定位和解決低頻振動問題的有效手段。

圖4 低頻振動信號Fig.4 Low frequency vibration signal

圖5 低頻振動的幅值包絡(luò)Fig.5 Envelope of low frequency vibration signal

圖6 初始狀態(tài)püA關(guān)系圖Fig.6 p-A relationship in initial status

圖7 更改狀態(tài)püA關(guān)系圖Fig.7 p-A relationship in improved status

4 結(jié)論

直升機中的“抖動”“晃動”“篩動”等低頻振動現(xiàn)象為直升機低階模態(tài)下響應(yīng)的幅值包絡(luò)。通過分析直升機低頻振動的現(xiàn)象和機理，基于希爾伯特變換，引入püA圖，形成了直升機低頻振動現(xiàn)象的分析方法，相對客觀地反映出機組人員的主觀感受，為分析、定位和解決直升機的低頻現(xiàn)象提供有效的手段。通過對某直升機橫側(cè)向抖動問題的分析，證明本文方法在分析直升機低頻振動現(xiàn)象的適用性和有效性。

[1] Pieter D W， Michel T. Tail shake vibration in fl ight： objective comparison of aerodynamic configurations in a subjective environment[C]// AHS 52nd Annual Forum， Montreal， Quebec，1999： 2306-2316.

[2] Hassan A A， Thompson T， Duque E P， et al. Resolution of tail buffet phenomena for AH-64D Longbow Apache： a case study on the application of a modern CFD analysis tool for improved helicopter aerodynamic design[C]// AHS 53rd Annual Forum，Virginia Beach， Virginia， 1997： 1028-1048.

[3] 航空航天工業(yè)部科學(xué)技術(shù)研究院.直升機動力學(xué)手冊[M].北京：航空工業(yè)出版社，1991.Science and Technology Institute in Ministry of Aerospace Industry. Handbook in helicopter dynamics[M]. Beijing：Aviation Industry Press， 1991.（in Chinese）

[4] 張曉谷.直升機動力學(xué)設(shè)計[M].北京：航空工業(yè)出版社，1995.ZHANG Xiaogu. Helicopter dynamics design [M]. Beijing：Aviation Industry Press， 1995.（in Chinese）

[5] 胡異丁，任偉新，楊棟，等.基于希爾伯特變換的非平穩(wěn)調(diào)幅信號解調(diào) [J].振動與沖擊，2013，32（10）：181-183.HU Yiding， REN Weixin， YANG Dong， et al. Demodulation of nonstationary amplitude modulated signal based on Hilbert transform [J]. Journal of Vibration and Shock， 2013， 32（10）：181-183.（in Chinese）

[6] 王幫峰，林劍祥，蘆吉云.基于EEMD-HT的飛行數(shù)據(jù)小突變信號檢測 [J].振動、測試與診斷，2013，33（3）：388-392.WANG Bangfeng， LIN Jianxiang， LU Jiyun. Detection of small mutation in fl ight data based on EEMD-HT [J]. Journal of Vibration， Measurement & Diagnosis， 2013， 33（3）：388-392.（in Chinese）

[7] 胡勁松，楊世錫，吳昭同，等.基于EMD和HT的旋轉(zhuǎn)機械振動信號時頻分析[J].振動、測試與診斷，2004，24（2）：106-110.HU Jingsong， YANG Shixi， WU Zhaotong， et al. Timefrequency analysis of vibration signal in rotating machinery[J].Journal of Vibration， Measurement & Diagnosis， 2004， 24（2）：106-110.（in Chinese）

[8] 楊宇，于德介，程軍圣.基于Hilbert邊際譜的滾動軸承故障診斷方法 [J].振動與沖擊，2005，24（1）：70-72.YANG Yu， YU Dejie， CHENG Junsheng. Roller bearing fault diagnosis based on Hilbert marginal spectrum [J]. Journal of Vibration and Shock， 2005， 24（1）：70-72.（in Chinese）

[9] 夏均忠，劉遠宏，李樹珉，等.應(yīng)用Hilbert變換和ZFFT提取變速器齒輪故障特征[J].振動與沖擊，2013，32（6）：63-66.XIA Junzhong， LIU Yuanhong， LI Shumin， et al. Gearbox fault detection using Hilbert transform and ZFFT [J]. Journal of Vibration and Shock， 2013， 32（6）： 63-66.（in Chinese）

[10] Bramwell A R， George D， David B. Bramwell’s helicopter dynamics[M]. Oxford： Butterworth-Heinemann， 2001.

[11] 胡海巖，孫久厚，陳懷海.機械振動與沖擊[M].北京：航空工業(yè)出版社，1998.HU Haiyan， SUN Jiuhou， CHEN Huaihai. Mechanical vibration and shock[M]. Beijing： Aviation Industry Press， 1998.（in Chinese）

[12] Hahn S L. Hilbert transform in signal processing[M].Boston：Artech House， 1996.