多點激勵功率譜再現(xiàn)振動試驗控制研究

第一作者王海東男，高級工程師，1975年6月生

多點激勵功率譜再現(xiàn)振動試驗控制研究

王海東1，欒強利2，陳章位2，賀惠農3(1.上海航天精密機械研究所，上海200233；2.浙江大學流體動力與機電系統(tǒng)國家重點實驗室，杭州310027；3.杭州億恒科技有限公司，杭州310015)

摘要：通過對多點激勵功率譜再現(xiàn)振動試驗控制算法研究.設計基于偏相干分析理論的振動試驗系統(tǒng)頻響矩陣辨識策略，針對系統(tǒng)頻響矩陣存在奇異點及系統(tǒng)頻響矩陣為長方矩陣情形，設計基于求解頻響矩陣廣義逆和矩陣最小范數最小二乘解的Moore-Penrose逆系統(tǒng)解耦算法。針對傳統(tǒng)差分修正驅動譜控制算法中存在系統(tǒng)功率譜自譜為負數或零值問題，通過引進比例均方根反饋修正算法，設計改進的功率譜均衡控制策略，有效避免功率譜均衡過程中自譜產生負值或零值問題。多點激勵功率譜再現(xiàn)振動試驗表明，改進的功率譜均衡控制策略對多點激勵系統(tǒng)具有可靠、高精度的控制效果。

關鍵詞：偏相干分析；系統(tǒng)解耦；比例均方根反饋修正；功率譜均衡

收稿日期：2014-01-02修改稿收到日期：2014-04-29

中圖分類號:O324; TB535文獻標志碼：A

Control designs for multi-exciter power spectral density replication vibration tests

WANGHai-dong1,LUANQiang-li2,CHENZhang-wei2,HEHui-nong3(1. Shanghai Spaceflight Precision Machinery Institute，Shanghai 200233，China;2.The State Key Lab of Fluid Power Transmission and Control, Zhejiang University, Hangzhou 310027, China;3.Hangzhou ECON Science and Technology Co., LTD, Hangzhou 310015, China)

Abstract:Through studying the control algorithm for multi-exciter PSD (power spectral density) replication vibration tests, an identification strategy for the frequency response matrix of a vibration test system based on partial coherence analysis was designed. For situations of singular points existing in a system frequency response matrix and the matrix being rectangular, different system decoupling algorithms were desgined based on the solutions to the matrix generalized inverse and its Moore-Penrose inverse with the minimum norm and least square methods. Aiming at the situation that the system APSD (auto-power spectral density) might be negative or zero in the traditional difference correction of the drive spectral control algorithm, an improved PSD equalization control strategy was designed by introducing a proportional RMS feedback correction algorithm. The algorithm could effectively avoid the APSD being negative and zero in the process of PSD equalization. The multi-exciter PSD replication vibration tests showed that the improved PSD equalization control strategy can achieve a reliable, high-precision control effect on a multi-exciter vibration system.

Key words:partial coherence analysis; system decoupling; proportional RMS feedback correction; PSD equalization

振動試驗通過模擬產品工作過程中可能承受的各種振動環(huán)境，可有效評估產品設計可靠性。由于對產品工作可靠性要求越來越苛刻，尤其航空航天領域，各類元器件、組合件及整機均需進行振動試驗，可及時暴露系統(tǒng)中各環(huán)節(jié)可能存在的缺陷及不足。以便快速整改及優(yōu)化，從而保證滿足產品工作可靠性[1-3]要求。功率譜再現(xiàn)試驗廣泛用于各類產品的振動試驗，通過在振動臺上復現(xiàn)規(guī)定參考譜，檢驗產品及關鍵零部件可靠性。

振動試驗據激勵點不同分為單、多點激勵。其中多點激勵能提供更大推力，充分考慮產品不同部位振動響應，尤其適用于航空航天類產品復雜化、大型化、重型化結構特點。多點激勵振動試驗[4-7]對多點激勵振動控制相關理論發(fā)展、完善具有開創(chuàng)性貢獻。Underwood等[8-10]則在多點激勵正弦試驗自適應控制、頻響函數矩陣奇異點控制及多輸入多輸出(Multiple-Input Multiple-Output)振動控制方面進行深入研究。國內研究已取得諸多理論及實踐突破[11-15]。本文在已有成果基礎上，利用偏相干分析理論進行多點激勵振動試驗系統(tǒng)辨識，并在功率譜均衡過程中通過采用改進比例均方根修正算法實現(xiàn)驅動信號更新；通過雙點同向激勵振動試驗對控制算法進行驗證。試驗表明，控制算法對振動試驗系統(tǒng)控制效果較好。

1傳遞函數辨識

傳遞函數是整個振動試驗系統(tǒng)的特性描述，因此其參數辨識的好壞對整個試驗系統(tǒng)控制效果有重要影響。傳遞函數偏相干分析理論通過將信號中與其它信號相干部分去掉，計算剩余信號對輸出信號影響[16-18]。多點激勵問題可表示為多個簡單多輸入單輸出(MISO)問題。簡單多輸入單輸出系統(tǒng)見圖1。

圖1　原始多輸入單輸出系統(tǒng) Fig.1 The original multiple-input single-output (MISO) system

該系統(tǒng)表達式為

(1)

式中：Xi為第i{i=1,2…q}個輸入信號頻譜；Y為輸出信號頻譜；Hiy為輸入Xi到輸出Y的頻響函數；N為系統(tǒng)中存在的噪聲干擾。

將條件輸入信號代替圖1中原始輸入信號，所得條件多輸入單輸出系統(tǒng)見圖2，其為原始多輸入單輸出系統(tǒng)的等效模型。

圖2　條件多輸入單輸出系統(tǒng) Fig.2 The conditional multiple-input single-output (MISO) system

該系統(tǒng)表達式為

(2)

式中：Xi(i-1)!為原始輸入信號Xi去除X1…Xi-1線性影響后殘差，且兩兩互不相關；Liy為等效模型頻響函數，不同于原始頻響Hiy。

定義條件系統(tǒng)(圖2)產生最小系統(tǒng)均方誤差的傳遞函數為

(3)

式中：Siy(i-1)!為Xi(i-1)!與Y的互譜；Sii(i-1)!為Xi(i-1)!的自譜。

將Y用Xj代替，得

(4)

由式(1)、(3)、(4)，并經推導得

(5)

式(5)為原始多輸入單輸出系統(tǒng)與條件多輸入單輸出系統(tǒng)之關系。原始系統(tǒng)頻響可由條件輸入系統(tǒng)頻響求得，Lij為對角線為1的上三角矩陣。將多輸入單輸出系統(tǒng)頻響估計方法用于多輸入多輸出系統(tǒng)，可實現(xiàn)對MIMO系統(tǒng)頻響估計，M輸入N輸出系統(tǒng)的頻響估計流程見圖3。

圖3　多輸入多輸出系統(tǒng)頻響估計流程 Fig.3 The flow chart for the estimation of system frequency response matrix of a multiple-input multiple-output (MIMO) system

2系統(tǒng)解耦

多點激勵功率譜再現(xiàn)過程中需求得系統(tǒng)頻響矩陣的逆，即需對系統(tǒng)進行解耦。當系統(tǒng)頻響矩陣為方陣(激勵點數等于響應點數)且非奇異時，系統(tǒng)解耦補償矩陣唯一確定，即系統(tǒng)頻響矩陣的逆。

系統(tǒng)頻響矩陣為方陣時，振動試驗過程中系統(tǒng)結構共振點處頻響矩陣易發(fā)生奇異，導致較大控制誤差，因此可采用奇異值分解方法求解奇異矩陣的廣義逆。當系統(tǒng)頻響矩陣為長方矩陣(激勵點數不等于響應點數)時，需求解長方矩陣具有最小范數最小二乘解的Moore-Penrose逆。對分塊長方矩陣L=(A,B)，則長方矩陣最小范數最小二乘解的Moore-Penrose逆為

L*=[A(I,A*B)]*=(I∶A*B)*A*=

(6)

3控制方法

3.1驅動信號生成

多點激勵振動試驗系統(tǒng)中驅動信號與控制信號頻譜可表示為

D(f)=[H(f)]-1C(f)

(7)

式中：D(f)，C(f)分別為驅動信號、控制信號頻譜；H(f)為系統(tǒng)頻響矩陣。

多點激勵功率譜再現(xiàn)振動試驗中[GRR]，[GDD]，[GCC]分別表示參考信號、驅動信號及控制信號功率譜密度，多點激勵功率譜再現(xiàn)目的為使系統(tǒng)控制信號功率譜[GCC]等于系統(tǒng)參考信號功率譜[GRR]。由式(7)得

[GDD]=[Z][GRR][Z]H

(8)

式中：[Z]=[H]-1為系統(tǒng)解耦補償矩陣；[GRR]為參考譜矩陣，可用n個獨立白噪聲向量{W}及1個下三角矩陣[L]表示，同時對驅動譜矩陣[GDD]作Cholesky分解，得

[GDD]=[D][D]H

(9)

DDH=[Z][L]{W}{W}H[L]H[Z]H

(10)

由此獲得系統(tǒng)驅動信號頻譜，即

D=[Z][L]{W}

(11)

3.2迭代修正算法

為實現(xiàn)多點激勵振動試驗系統(tǒng)高精度功率譜再現(xiàn)，需對試驗系統(tǒng)采用閉環(huán)控制策略。由于傳統(tǒng)的差分修正驅動譜控制算法據控制點位置不同賦予不同優(yōu)先權，導致系統(tǒng)參考譜中位置靠前的自譜控制效果較好，而位置靠后的自譜控制效果較差，甚至導致自譜負數。為此，本文采用比例均方根反饋修正算法[19]，即

[L]i+1=[L]i[Δ]i

(12)

(13)

(14)

式中：i為迭代次數。

當參考譜與響應譜存在誤差時，采用比例均方根修正算法對下三角矩陣[L]進行修正，可避免功率譜對角元素出現(xiàn)負值或零值；但比例均方根修正算法主要針對功率譜對角元素修正，并不單獨考慮非對角元素，互譜控制在修正自譜中自行完成，將導致互譜控制無法實現(xiàn)準確、高精度控制。因此，采用振動系統(tǒng)頻響矩陣實時修正算法，實現(xiàn)功率譜矩陣的自譜、互譜控制。改進的功率譜再現(xiàn)控制流程見圖4。為避免由于控制譜相干系數過高導致功率譜矩陣奇異，當目標控制譜相干系數小于0.95時，可采用頻響矩陣修正方法；當相干系數大于0.95時采用修正參考譜矩陣下三角矩陣[L]方式實現(xiàn)驅動信號修正。本文綜合考慮兩種控制修正方法，設計出改進功率譜均衡控制策略，以逆補各自不足，實現(xiàn)多點激勵功率譜再現(xiàn)振動控制的高精度。

圖4　改進的功率譜再現(xiàn)控制框圖 Fig.4 The control block diagram of the improved PSD replication

4振動試驗

為進一步驗證控制算法對多點激勵功率譜再現(xiàn)振動試驗的控制效果，建立振動試驗系統(tǒng)見圖5。系統(tǒng)中選基于控制算法的ECON多點激勵振動控制器，并采用雙點激勵、響應控制。振動臺及加速度傳感器參數見表1。

圖5　雙輸入雙輸出振動試驗系統(tǒng) Fig.5 A double-input double-output vibration test system

功率譜再現(xiàn)振動試驗中，信號頻率范圍20～2 000 Hz，自譜相同，互譜相位差為0°，相干系數設為1。系統(tǒng)的解耦補償(阻抗)曲線見圖6。由圖6看出，系統(tǒng)在160 Hz、602 Hz、1 245 Hz附近有奇異點。系統(tǒng)自譜、互譜控制曲線分別見圖7、圖8。由圖7、圖8可知，在系統(tǒng)具有明顯反共振點情況下(圖6)，控制算法控制效果較好，自功率譜幅值控制在±1.5 dB范圍，互功率譜相干函數及相位波動控制在有限波動范圍，保證系統(tǒng)自功率譜的高精度再現(xiàn)。

表1　振動試驗系統(tǒng)參數

圖6　振動試驗系統(tǒng)阻抗曲線 Fig.6 The impedance curves of the vibration test system

圖7　振動試驗系統(tǒng)自譜控制曲線 Fig.7 The APSD control curves of the vibration test system

圖8　振動試驗系統(tǒng)互譜控制曲線 Fig.8 The CPSD control curves of the vibration test system

5結論

(1)通過對多點激勵功率譜再現(xiàn)振動試驗相關理論、算法深入研究，設計出基于偏相干理論的系統(tǒng)頻響矩陣估計方法；針對系統(tǒng)頻響矩陣中奇異現(xiàn)象及系統(tǒng)長方矩陣存在，設計出相關的解耦控制算法；針對傳統(tǒng)頻譜均衡修正算法中存在可能導致系統(tǒng)功率譜自譜為負數或零值問題，設計出改進的功率譜均衡控制策略。

(2)通過對雙點激勵、控制振動試驗系統(tǒng)進行隨機振動試驗表明，振動系統(tǒng)存在奇異點時，通過控制策略可實現(xiàn)振動系統(tǒng)可靠、準確的功率譜再現(xiàn)。

參考文獻

[1]王述成. 振動試驗實時控制系統(tǒng)的研究[D]. 杭州：浙江大學，2006.

[2]齊曉軍. 航天器振動試驗控制技術研究[D]. 長沙：國防科學技術大學，2011.

[3]王珂晟,雷勇軍,朱曉瑩. 系統(tǒng)級產品振動試驗的探討與研究[J]. 振動與沖擊, 2004, 23(4): 112-115.

WANG Ke-sheng, LEI Yong-jun, ZHU Xiao-ying. Investigation on vibration testing for system_level equipments[J]. Journal of Vibration and Shock,2004, 23(4):112-115.

[4]Smallwood D O. Multiple shaker random control with cross coupling [A]. Proceedings of the IES[C].1978: 341-347.

[5]Smallwood D O. Random vibration testing of a single test item with a multiple input control system[A]. Proceedings of the IES[C].1982: 42-49.

[6]Smallwood D O, Woodall T D, Buksa E J. Minimum driver requirements for a multiple input multiple output linear system [A]. Proceedings of the IES[C].1986: 295-301.

[7]Smallwood D O. Multiple shaker random vibration control-an update [A]. Proceedings of the IEST[C].1999:212-221.

[8]Underwood M A, Keller T. Recent system developments for multi-actuator vibration control[J]. Sound and Vibration, 2001, 35(10):16-23.

[9]Underwood M A. Adaptive control method for multiexciter sine tests[P]. United States Patient:5299459, 1994-04-05.

[10]Underwood M A,Keller T. Apparatus and method for adaptive closed loop control of shock testing system[P]. United States Patent, 5517426, 1996-05-14.

[11]陳家焱. 多點激勵振動試驗系統(tǒng)的控制策略研究[D]. 杭州：浙江大學,2011.

[12]賀旭東,陳懷海. 多點隨機振動控制中的互譜矩陣研究[J].南京航空航天大學學報,2004, 36(6): 744-747.

HE Xu-dong, CHEN Huai-hai. Cross-spectra in multiple shaker random vibration test [J]. Journal of Nanjing University of Aeronautics & Astronautics,2004,36(6):744-747.

[13]賀旭東,陳懷海. 一種多點隨機振動試驗控制的新方法研究[J]. 振動工程學報,2004, 17(1): 49-52.

HE Xu-dong, CHEN Huai-hai. A new method for the control of multi-shakers in random vibration tests[J]. Journal of Vibration Engineering, 2004, 17(1): 49-52.

[14]葉建華,李傳日. 多點隨機振動試驗控制技術[J]. 系統(tǒng)工程與電子技術,2008, 30(1):124-127.

YE Jian-hua, LI Chuan-ri. Systems engineering and electronics[J]. Journal of Systems Engineering and Electronics, 2008, 30(1):124-127.

[15]陳家焱,王海東,周建川,等. 多點激勵振動試驗控制技術進展[J]. 振動與沖擊,2011, 30(3): 69-73.

CHEN Jia-yan, WANG Hai-dong, ZHOU Jian-chuan, et al. Progress in multi-exciter vibration testing control technology [J]. Journal of Vibration and Shock, 2011, 30(3): 69-73.

[16]馬忠成. 偏相干方法分析及其工程應用[J]. 應用聲學, 1994, 13(3): 25-30.

MA Zhong-cheng. Analysis and engineering application of the partial coherence method[J]. Applied Acoustics,1994,13(3):25-30.

[17]趙海瀾,汪鴻振. 偏相干分析識別噪聲源的計算[J].噪聲與振動控制,2005(5):31-33.

ZHAO Hai-lan, WANG Hong-zhen. Calculation of identifying noise source based on partial coherence analysis [J]. Noise and Vibration Control,2005(5): 31-33.

[18]貝達特J S，皮爾索A G. ,著.凌福根，譯.相關分析和譜分析的工程應用[M]. 北京：國防工業(yè)出版社，1983.

[19]賀旭東,陳懷海,申凡,等. 雙振動臺隨機振動綜合控制研究[J]. 振動工程學報, 2006, 19(2): 145-149.

HE Xu-dong, CHEN Huai-hai, SHEN Fan, et al. Study on dual-shaker random vibration test control[J]. Journal of Vibration Engineering, 2006, 19(2): 145-149.