MIMO混合振動試驗中正弦與隨機(jī)信號分離方法

張步云， 陳懷海， 賀旭東

(南京航空航天大學(xué)振動工程研究所， 江蘇 南京 210016)

引 言

振動環(huán)境試驗作為檢驗產(chǎn)品可靠性與耐久性的有效手段，近半個世紀(jì)取得了長足的進(jìn)展。多輸入多輸出振動試驗比單軸振動試驗更準(zhǔn)確地模擬外場振動環(huán)境，能夠克服單軸推力不足、量級不高、試驗耗時較長等缺點(diǎn)，越來越受到研究者們的重視。在正弦、隨機(jī)和沖擊等振動試驗研究中[1～6]，人們已建立了成熟的理論體系，完善了先進(jìn)的振動試驗方法。國內(nèi)外對相應(yīng)的振動試驗都制定了完整的試驗標(biāo)準(zhǔn)與試驗規(guī)范[7，8]，使振動試驗更具實用性與權(quán)威性。但已往的研究大都針對于某一單類振動試驗的研究，而事實上工程機(jī)械或飛行器在運(yùn)作過程中所經(jīng)歷的振動環(huán)境是相當(dāng)復(fù)雜的。如直升機(jī)在飛行中由于漿片的旋轉(zhuǎn)而產(chǎn)生一定頻率的正弦振動，而引擎的工作則產(chǎn)生一個寬帶的隨機(jī)噪聲。若用單一的正弦振動或隨機(jī)振動來模擬機(jī)身所受的振動環(huán)境顯然已經(jīng)不準(zhǔn)確，得到的試驗結(jié)果不能用來作為飛機(jī)疲勞破壞的程度或損傷的判斷標(biāo)準(zhǔn)。正弦加隨機(jī)混合振動試驗可以更好地完成對這種振動環(huán)境的模擬與再現(xiàn)，故其已成為振動試驗研究的新課題[9，10]。

正弦振動試驗中以信號的幅值與相位為控制對象，而隨機(jī)振動試驗則是通過閉環(huán)反饋將控制點(diǎn)的功率譜控制在給定的譜容差范圍內(nèi)。常用的正弦加隨機(jī)控制方法是將正弦信號與隨機(jī)信號分離后分別用已成熟的控制方法加以控制，再將反饋的控制信號疊加輸送到振動臺進(jìn)行激振。正弦與隨機(jī)信號分離的好壞將決定控制的精度，故其成為正弦加隨機(jī)振動控制中最關(guān)鍵的環(huán)節(jié)。目前已有的分離方法有跟蹤濾波法[11]，將混合信號通過一比隨機(jī)信號分析頻率大很多的帶通濾波器分離出正弦信號，此法高頻分離尚可，低頻的控制精度較差。王述成針對低頻控制精度差的問題，提出了基于多分辨譜的控制方法[12]。借助維納-辛欽定理從頻域獲取正弦信號的幅值，但頻域識別存在能量泄露問題，不能保證控制的精度。Engelhardt提出用Vold-Kalman濾波器法[13]，其分離方法比較復(fù)雜，且難以應(yīng)用到振動試驗控制之中。此外，還有計算法、截取法、復(fù)數(shù)二乘法和基于最小二乘法、盲分離等分離算法等[14～17]。這些方法運(yùn)算復(fù)雜、計算量大，增加了控制回路的時間。

相關(guān)積分法可以有效地抑制測量噪聲的影響，能提高測試信號的信噪比，是一種常用的信號處理方法。本文利用相關(guān)積分法的濾波特性，詳細(xì)推導(dǎo)了正弦信號分離公式，有效分離了正弦與隨機(jī)信號，提高了混合振動試驗的控制精度。文中以一懸臂梁為仿真對象，驗證了此方法的有效性。

1 多輸入多輸出混合振動控制

多輸入多輸出正弦加隨機(jī)混合振動控制的整個流程如圖1所示。先將混合參考譜分解為正弦參考值和隨機(jī)功率譜，再將由實測信號分離得到的正弦值和實測隨機(jī)信號的功率譜與參考譜進(jìn)行比較，經(jīng)過控制算法均衡后得到新的正弦信號和隨機(jī)信號，將兩者綜合后輸入激振器進(jìn)行激振。如此循環(huán)，直到控制點(diǎn)的振動滿足參考譜為止。

圖1 正弦加隨機(jī)混合振動控制流程

1.1 隨機(jī)振動試驗控制

頻譜再現(xiàn)式隨機(jī)振動試驗是在較短的時間內(nèi)“同時”激發(fā)給定功率譜譜形的寬帶隨機(jī)振動[17]。設(shè)參考譜為Sr(ω)，根據(jù)振動理論，激勵與響應(yīng)的關(guān)系為

Sy(ω)=H(ω)Sx(ω)H*(ω)

(1)

式中Sy(ω)表示響應(yīng)信號的功率譜密度，Sx(ω)表示激勵信號的功率譜密度，H(ω)表示系統(tǒng)的頻響函數(shù)矩陣，其上標(biāo)*表示矩陣的共軛轉(zhuǎn)置。為方便計，下文省略ω。隨機(jī)振動試驗的控制目標(biāo)就是使Sy=Sr，則要求激勵信號的功率譜滿足

Sx=H-1Sy(H*)-1=ASrA*

(2)

式中A=H-1。對激勵譜Sx和參考譜Sy分別進(jìn)行Cholesky分解得

Sx=XX*，Sr=RR*

(3)

式中X，R均為下三角矩陣，將其帶入式(2)，由Cholesky分解的唯一性可得

X=AR

(4)

將X配以隨機(jī)相位，進(jìn)行逆FFT變換就可得到激勵信號，再進(jìn)行時域隨機(jī)化便可得到純隨機(jī)激勵信號。因系統(tǒng)中的非線性因素以及頻響函數(shù)估計誤差的存在，實際所測頻響函數(shù)G并不是系統(tǒng)真正的頻響，設(shè)它們的誤差可用如下式表示

GA=I+ΔI

(5)

其中I為單位矩陣，ΔI為誤差矩陣。則由式(4)所產(chǎn)生的激勵信號所得到的響應(yīng)譜為

(I+ΔI)(ΔR)R((ΔR)R)*(I+ΔI)*=Sr

(7)

可計算出ΔR對角元素為

(8)

1.2 正弦振動試驗控制

以某幾個固定頻率的正弦信號的疊加為激勵，采集響應(yīng)信號并進(jìn)行幅值識別與參考值進(jìn)行比較，計算出誤差。然后反饋修正激勵信號，使每個頻率點(diǎn)的響應(yīng)都能滿足參考值的要求。一般正弦修正算法為

(9)

式中Ai表示當(dāng)前驅(qū)動信號幅值，Ai-1表示前一次驅(qū)動幅值，E為前一次響應(yīng)幅值與參考值的比，K表示壓縮比，其大小影響響應(yīng)信號收斂的速度。

假設(shè)初始響應(yīng)信號與參考值之比為1.1，取不同的K觀察響應(yīng)收斂程度。分別取K為0.5，1，2，6，其收斂情況如圖2所示。

圖2 壓縮比K對收斂速度的影響

從圖中可以看出收斂速度與K大小成反比，K取值越小，收斂速度越快。

2 混合信號中正弦分離

設(shè)正弦加隨機(jī)混合信號可表示為如下形式

(10)

式中xr表示均值為零的平穩(wěn)隨機(jī)信號；混合信號中共含n個正弦信號，幅值為Ak，頻率為ωk，相位為φk。因隨機(jī)信號與正弦信號是相互獨(dú)立的量，所以對任一頻率的正弦函數(shù)，與隨機(jī)信號xr有如下關(guān)系

(11)

在相關(guān)積分計算過程中，隨機(jī)信號與正弦信號的計算是基于時間域的。它們乘積的相關(guān)積分函數(shù)不嚴(yán)格歸零，但其積分值遠(yuǎn)低于原信號量級，可以忽略。式(10)與待識別頻率正弦函數(shù)sin(ωpt)的相關(guān)積分函數(shù)為

(12)

式中Tp表示sin(ωpt)的周期。同理，式(10)與待識別頻率正弦函數(shù)cos(ωpt)的相關(guān)積分函數(shù)

(13)

由式(12)和(13)可以求得頻率為ωp的正弦信號的幅值與相位為

(14)

按式(12)～(14)分別求出n個正弦信號，然后將其從混合信號中減掉，就可以很容易地得到隨機(jī)信號

(15)

下面以某一直升機(jī)振動環(huán)境為研究對象，分析正弦信號與隨機(jī)信號的分離。根據(jù)文獻(xiàn)[8]，可得該振動環(huán)境如圖3所示。圖中參數(shù)如表1所示。

表1 直升機(jī)振動環(huán)境參數(shù)

圖3 某直升機(jī)振動環(huán)境

由表1和圖3可以看出直升機(jī)的振動環(huán)境是由正弦振動和隨機(jī)振動組合而成的。而且正弦的頻率相對較低，以一個帶寬為500 Hz的低量級隨機(jī)振動作為振動背景。混合振動信號可以寫成

(16)

由式(16)可得時域信號如圖4所示。

圖4 混合時域信號

利用相關(guān)積分法和頻域法同時識別式(16)中4個頻率的正弦信號幅值，以比較兩種方法識別精度的高低。頻域法的識別流程是先對混合信號作FFT變換，然后根據(jù)維納-辛欽定理獲取正弦信號的幅值；再將兩者識別的幅值同時與參考值作比較，計算出識別誤差。經(jīng)過計算，識別結(jié)果如表2所示。

表2 正弦信號幅值識別誤差分析

表2顯示相關(guān)法識別正弦信號結(jié)果精度很高，絕對誤差最高不足0.45%。而頻域法識別結(jié)果誤差很大，最少是相關(guān)法識別誤差的4倍以上。固從表2可知，相關(guān)法識別的精度結(jié)果可進(jìn)一步進(jìn)行振動控制試驗。

3 仿真算例

本文以一懸臂梁作為研究對象，進(jìn)行兩輸入兩輸出混合振動隨機(jī)控制仿真試驗，模型尺寸及示意圖見下圖。用Matlab軟件編制程序進(jìn)行控制仿真試驗，各階阻尼比取0.03。

圖5 試驗?zāi)Ｐ统叽缂笆疽鈭D (單位:mm)

控制點(diǎn)的正弦成分如表1，而隨機(jī)成分的參考譜如圖6所示。

圖6 控制點(diǎn)參考譜譜形

注意到式(2)，在求驅(qū)動信號時須求頻響的逆陣。在結(jié)構(gòu)的反共振峰處，頻響的幅值較小，甚至該頻率處的矩陣是奇異的。若直接求逆陣，可能會發(fā)生病態(tài)，此時須求頻響矩陣的廣義逆H+。所以頻響函數(shù)的準(zhǔn)確估計成為振動控制的一個關(guān)鍵環(huán)節(jié)，本文采用快速正弦掃描法估計頻響函數(shù)，此法能獲得平譜，在整個測試頻段內(nèi)激勵能量相同，信噪比大，能夠提高估計精度。經(jīng)過仿真計算，得到控制結(jié)果如分別如圖7～10所示。

圖7 控制點(diǎn)1隨機(jī)振動控制結(jié)果

圖8 控制點(diǎn)2隨機(jī)振動控制結(jié)果

圖9 控制點(diǎn)1正弦振動控制結(jié)果

圖10 控制點(diǎn)2隨機(jī)振動控制結(jié)果

圖7,8表示混合控制中隨機(jī)振動的控制結(jié)果。圖中虛線表示控制警告線，即參考譜的±3 dB值處；上下兩條實線表示控制停止線，為參考譜的±6 dB值處。經(jīng)過閉環(huán)反饋控制，可以從圖中看出兩點(diǎn)的響應(yīng)譜都被控制在警告線以內(nèi)。控制點(diǎn)1和2參考譜的總均方根值分別為1.77g和2.1g，而圖7，8中兩點(diǎn)響應(yīng)譜的總均方根值為1.78g和2.15g，誤差分別為0.6%和2.4%。綜前，隨機(jī)控制滿足控制要求。

圖8，9表示正弦振動控制結(jié)果。混合振動中含有4個定頻振動，頻率分別為7，18，34和51 Hz，以正弦加速度(g)為控制對象。從圖中可以看出，正弦幅值控制誤差滿足工程要求，即可以實現(xiàn)振動控制。

4 結(jié) 論

本文將相關(guān)積分法運(yùn)用到隨機(jī)加正弦混合振動試驗中，利用其濾波特性將隨機(jī)信號與正弦信號精確分離，從而分別進(jìn)行隨機(jī)振動控制和正弦振動控制，解決了混合振動控制中的關(guān)鍵問題。通過建立兩輸入兩輸出振動控制系統(tǒng)，利用Matlab軟件編制控制程序，對系統(tǒng)進(jìn)行混合控制。仿真結(jié)果表明：隨機(jī)振動信號被控制在參考譜的±3 dB以內(nèi)，正弦信號幅值被控制在參考值的±10%以內(nèi)，實現(xiàn)了良好的控制效果。

參考文獻(xiàn):

[1] Underwood M A. Adaptive control method for multi- exciter sine tests[P]. USA: US5299459 A, 1994,4,5.

[2] Smallwood D O. Multiple shaker random vibration control-an update[R]. Institute of Environmental Sciences and Technology (IEST) Annual Meeting, SAND98 -2044C, USA, US Department of Energy,1999:212—221.

[3] 賀旭東,陳懷海.一種多點(diǎn)隨機(jī)振動試驗控制的新方法研究[J].振動工程學(xué)報,2004,17(1):49—52.He X D, Chen H H. A new method for the control of multi-shakers in random vibration tests[J]. Journal of Vibration Engineering, 2004,17(1):49—52.

[4] 崔旭利,陳懷海,賀旭東,等.全相位Hv頻響函數(shù)估計在MIMO隨機(jī)振動試驗中的應(yīng)用[J].振動工程學(xué)報,2011,24(2):181—185.Cui X L, Chen H H, He X D, et al. All-phase Hvfrequency response function estimation for MIMO random vibration test control[J]. Journal of Vibration Engineering, 2011,24(2):181—185.

[5] Baren J V. What is random vibration testing[J]. Sound and Vibration, 2012,2:9—12.

[6] Underwood M A. Apparatus and method for adaptive closed loop control of shock testing system[P]. USA: US5517426A, 1996,5,14.

[7] Department of Defense of United States of America. MIL-STD-810G-2008[S]. Environmental engineering conside- rations and laboratory tests.

[8] 中國人民解放軍總裝備部.GJB 150.16A-2009.軍用裝備實驗室環(huán)境試驗方法第16部分: 振動試驗[S].Chinese PLA General Armament Department. GJB 150.16A-2009. Laboratory environmental test methods for military materiel, Part 16: Vibration Test[S].

[9] 孫成寬.多輸入多輸出隨機(jī)加正弦振動控制系統(tǒng)的研究[D].南京:南京航空航天大學(xué),2010.Sun C K. Research on MIMO sine on random vibration control system[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2010.

[10] 上海市質(zhì)量監(jiān)督檢驗技術(shù)研究院. GB/T 2423.58-2008,電工電子產(chǎn)品環(huán)境試驗.第2-80部分:試驗方法.試驗Fi:振動.混合模式[S].北京:國家質(zhì)檢總局, 2009.Shanghai Institute of Quality Inspection and Technical Research. GB/T 2423.58-2008, Enviromental testing for electric and electonic products. Part2: Tests methods. Test Fi: Vibration. Mixed Mode[S]. Beijing: CN-GB, 2009.

[11] 袁宏杰,李傳日.正弦加隨機(jī)振動控制技術(shù)的研究[J].航空學(xué)報,2000,21(4):383—384.Yuan H J, Li C R. Study on the sine-on-random vibration control[J]. Acta Aeronautica et Astronautica Sinica, 2000,21(4):383—384.

[12] 王述成,陳章位.基于多分辨譜分析的正弦加隨機(jī)振動試驗控制算法的研究[J].中國機(jī)械工程,2005,16(15):1 335—1 338.Wang S C, Chen Z W. A new algorithm for sine on random vibration test based on multirate spectrum estimation[J]. China Mechanical Engineering, 2005,16(15):1 335—1 338.

[13] Engelhardt C, Baker M, Mouron A, et al. Separation of sine and random components from vibration measurements[J]. Sound and Vibration, 2012,6:6—11.

[14] 龔慶祥,陳應(yīng)汀.頻率分析中正弦和隨機(jī)的分離與判辨[J].洪都科技,1979,4:1—6.Gong Q X, Chen Y T. Separation and analysis of sine and random signals in the frequency analysis[J]. Hongdu Technology, 1979,4:1—6.

[15] 袁宏杰,李傳日. VCS振動控制系統(tǒng)控制算法研究[J].儀器儀表學(xué)報(S), 2002,23(3):846—847.Yuan H J, Li C R. VCS vibration control system control algorithm[J]. Chinese Journal of Scientific Instrument (S), 2002,23(3):846—847.

[16] 江勇,張林柯,何琳,等.兩次去相關(guān)用于振動信號盲分離[J].振動、測試與診斷,2011,31(2):241—245.Jiang Y, Zhang L K, He L, et al. Blind separation of vibration signals using double decorrelation[J]. Journal of Vibration, Measurement & Diagnosis, 2011,31(2):241—245.

[17] 胡志強(qiáng).隨機(jī)振動試驗應(yīng)用技術(shù)[M].北京:中國計量出版社,1998:118.Hu Z Q. Application Technology of Random Vibration Test[M]. Beijing: China Metrology Publishing House, 1998:118.

[18] 賀旭東,陳懷海,申凡,等.雙振動臺隨機(jī)振動綜合控制研究[J].振動工程學(xué)報,2006,19(2):145—149.He X D, Chen H H, Shen F, et al. Study on dual-shaker random vibration test control[J]. Journal of Vibration Engineering, 2006,19(2):145—149.