旋翼振動監(jiān)測信號的演化分析與損傷跟蹤方法*

朱旭程

(海軍航空工程學(xué)院飛行器工程系 煙臺，264001)

?

旋翼振動監(jiān)測信號的演化分析與損傷跟蹤方法*

朱旭程

(海軍航空工程學(xué)院飛行器工程系 煙臺，264001)

為獲得一種對損傷敏感而對正常狀態(tài)擾動不敏感的旋翼損傷觀測信號、解決直升機(jī)旋翼在飛行狀態(tài)下的損傷監(jiān)測問題，研究了旋翼在重構(gòu)相空間中的損傷演化特性并提出一種新的損傷跟蹤方法。該方法首先應(yīng)用嵌入技術(shù)將旋翼氣彈有限元模型產(chǎn)生的振動監(jiān)測信號重構(gòu)到維數(shù)更高的相空間中，采用非線性Volterra級數(shù)建立旋翼基準(zhǔn)狀態(tài)預(yù)測模型，以基準(zhǔn)狀態(tài)預(yù)測結(jié)果與實測相軌跡之件的差異作為狀態(tài)預(yù)測殘差，在多個鄰域內(nèi)對其統(tǒng)計平均形成一個損傷觀測特征向量；然后,應(yīng)用奇異譜分解方法從損傷觀測特征向量的時間序列中提取出槳葉損傷演化的維度和趨勢信息，采用雙指數(shù)平滑方法建立損傷演化趨勢預(yù)測模型并估計出槳葉損傷故障剩余壽命。采用槳葉損傷模型和旋翼氣彈模型仿真數(shù)據(jù)驗證了所提方法的可行性和有效性。結(jié)果表明:該方法能充分利用監(jiān)測信號的非線性特性在高維相空間中重建系統(tǒng)動力學(xué)本質(zhì)，以不同的時間尺度來觀測系統(tǒng)的演化特性；具有損傷模式自動識別能力，可用于難以事先確定系統(tǒng)損傷演化模型或維數(shù)信息的場合。

直升機(jī)；旋翼；健康監(jiān)測；故障診斷；相空間

引言

旋翼是直升機(jī)飛行的關(guān)鍵動力部件，在槳葉周期變距和氣動交變載荷作用下槳葉、槳轂、減擺器和變距拉桿等部件會出現(xiàn)磨損和疲勞失效故障。能夠準(zhǔn)確監(jiān)測旋翼的損傷狀況，對于提高直升機(jī)飛行安全性極具現(xiàn)實意義。文獻(xiàn)[1]對旋翼損傷監(jiān)測問題開展了研究，提出了通過槳葉相對位移和機(jī)體振動識別旋翼故障的方案。Ganguli等[2-3]分析了典型損傷對旋翼氣彈響應(yīng)信號的影響，運(yùn)用神經(jīng)網(wǎng)絡(luò)、模糊推理等人工智能技術(shù)建立了一些旋翼故障診斷算法。目前，旋翼在飛行條件下的損傷監(jiān)測技術(shù)還存在瓶頸[4]：在高速旋轉(zhuǎn)的旋翼上加裝各種數(shù)據(jù)采集系統(tǒng)尚存在一些技術(shù)上的限制；此外，由于槳距時變、氣彈耦合和非定常尾流等因素影響，傳統(tǒng)損傷跟蹤方法難以獲得一種對損傷敏感而對正常狀態(tài)擾動不敏感的損傷觀測信號，有效的旋翼損傷跟蹤方法依然缺乏。

筆者應(yīng)用混沌信號嵌入分析理論的相空間重構(gòu)技術(shù)，將旋翼狀態(tài)監(jiān)測數(shù)據(jù)變換到更高維數(shù)的狀態(tài)空間內(nèi)、以不同的時間尺度來研究旋翼系統(tǒng)的損傷演化特性。采用非線性Volterra級數(shù)建立旋翼基準(zhǔn)狀態(tài)預(yù)測模型，用基準(zhǔn)狀態(tài)預(yù)測結(jié)果與實測相軌跡之間的差異作為狀態(tài)預(yù)測殘差，發(fā)展一種新的旋翼損傷跟蹤方法。

1 旋翼在重構(gòu)相空間的損傷演化特性

1.1 旋翼氣彈有限元模型

為研究損傷對槳葉位移、槳轂力等旋翼狀態(tài)監(jiān)測信號的影響，需要建立旋翼氣彈有限元模型。由Hamilton能量變分原理

(1)

其中：Π為系統(tǒng)總能；T為動能；U為應(yīng)變能；W為外力虛功。

氣動載荷采用葉素理論建立，旋翼誘導(dǎo)速度采用White-Blake入流模型[5]計算。

采用有限元方法離散上述能量積分，進(jìn)一步采用分部積分法建立旋翼氣彈有限元方程

(2)

其中：M，C，K，F(xiàn)分別表示慣性矩陣、阻尼矩陣、剛度矩陣和載荷向量，在定常平飛狀態(tài)下它們?yōu)榉轿唤堑闹芷诤瘮?shù)。

旋翼氣彈有限元方程可采用時域有限元方法化為代數(shù)方程求解，或者將其轉(zhuǎn)換為非線性狀態(tài)方程采用數(shù)值積分方法求解。

旋翼的故障模式種類多樣，不同類型的故障對系統(tǒng)參數(shù)會有不同影響。采用復(fù)合材料槳葉剛度疲勞損傷演化模型說明本研究方法的有效性。

(3)

槳葉損傷變量隨使用壽命的擴(kuò)展如圖1所示。

圖1 槳葉剛度損傷變量擴(kuò)展曲線Fig.1 Damage accumulation curve of blade stiffness

目前，直升機(jī)健康和使用監(jiān)測系統(tǒng)(health and usage monitor system,簡稱HUMS)能夠提供的旋翼狀態(tài)監(jiān)測信號主要有二類：a. 槳葉的揮舞、擺動和扭轉(zhuǎn)位移，可通過在機(jī)身上安裝激光多普勒測振儀實時監(jiān)測[6]；b. 槳轂力和力矩，可通過機(jī)身上的應(yīng)變計或振動傳感器結(jié)合飛行參數(shù)得到。下面以槳尖位移、槳轂力和力矩作為旋翼典型狀態(tài)監(jiān)測信號。

將旋翼氣彈有限元模型和槳葉損傷演化模型結(jié)合，分析槳葉損傷對槳尖位移和槳轂載荷的影響。筆者以BO-105直升機(jī)旋翼為參考設(shè)置實例旋翼的基本參數(shù)：槳葉數(shù)為4，旋翼半徑R為4.94 m，旋翼轉(zhuǎn)速Ω為40 rad/s，槳葉密度m為6.46 kg/m，槳葉弦長為0.055R，弦向剛度為0.026 8mΩ2R4，揮向剛度為0.010 8mΩ2R4，扭轉(zhuǎn)剛度為0.006 15mΩ2R4，洛克數(shù)為5.2，槳葉升力線斜率為6.0，型阻系數(shù)為0.002，誘導(dǎo)阻力系數(shù)為0.2，俯仰力矩系數(shù)為0.00。飛機(jī)在前進(jìn)比0.3，旋翼拉力系數(shù)0.004 9的條件下定常平飛。旋翼氣彈模型計算程序采用Matlab軟件建立(如圖2)，槳葉在健康(t/tF=0)和損傷(t/tF=1/3)狀態(tài)下仿真產(chǎn)生的槳尖位移、槳根力和力矩如圖3所示。圖中:w為槳尖揮舞位移；v為槳尖擺動位移；Fz為槳轂力；Mx為槳轂力矩；n為樣本長度。

圖2 旋翼氣彈模型計算流程Fig.2 Flowchart of rotor simulation procedure

圖3 健康和損傷旋翼槳尖位移和槳轂載荷Fig.3 Health and damaged blade tip deflections and Hub loads

1.2 旋翼狀態(tài)監(jiān)測信號的相空間重構(gòu)

旋翼損傷監(jiān)測的一個基本問題是找到一種對損傷敏感而對系統(tǒng)狀態(tài)變化、干擾噪聲不敏感的觀測信號。然而槳葉損傷與氣彈耦合現(xiàn)象會引起氣動載荷沿槳葉結(jié)構(gòu)的分布發(fā)生變化[7]，隨著槳葉損傷程度的擴(kuò)展，在揮舞、擺動、扭轉(zhuǎn)方向的一些固有頻率互相穿越與交叉，使得槳尖位移和槳轂載荷的峰值等特征信號在單調(diào)性關(guān)系上出現(xiàn)突變，因而難以由這些特征信號辨識出損傷的變化趨勢。源自混沌理論的相空間重構(gòu)技術(shù)可以在不同維數(shù)的空間內(nèi)、以不同的時間尺度來研究演化系統(tǒng)的動力學(xué)特性，為非線性振動系統(tǒng)的損傷跟蹤提供了一條新途徑。只要時間尺度和空間維數(shù)選擇合適，可通過相軌跡的移動變化定量跟蹤系統(tǒng)損傷的演化過程[8]。將某監(jiān)測信號的時序數(shù)據(jù){y(i):i=1,2,…,M} 進(jìn)行時間坐標(biāo)延遲，構(gòu)造一狀態(tài)點列 {x(n):n=1,2,…,N}。其中：x(n) = [y1(n),y2(n), …,ym(n)]T；yk(n) =y(n+(k-1)τ)；m為嵌入維數(shù)；τ為延遲時間。該點列描述了系統(tǒng)在m維狀態(tài)空間中的相軌跡，由它構(gòu)成的空間稱為相空間，這個過程稱為相空間重構(gòu)。根據(jù)嵌入理論[9]，如果重構(gòu)系統(tǒng)與原動力學(xué)系統(tǒng)的動態(tài)響應(yīng)可限制在相同維數(shù)的吸引子流形上，重構(gòu)系統(tǒng)與原系統(tǒng)在拓?fù)潢P(guān)系或微分意義上屬于同一等價(微分同胚)，具有相同的動力學(xué)性質(zhì)和幾何性質(zhì)，系統(tǒng)的許多特性如分維、K熵、李指數(shù)等均相同。

圖4 平均互信息法和假近鄰法求重構(gòu)參數(shù)流程Fig.4 Flowchart for time delay and embedding dimension based on average mutual information and FNN percentage

圖5 相軌跡在三維空間中的投影Fig.5 3D views of reconstructed phase space trajectories

為直觀表明旋翼在相空間中的損傷演化過程，將帶槳葉損傷BO-105旋翼的槳尖位移和槳轂載荷數(shù)據(jù)進(jìn)行相空間重構(gòu)。在相空間重構(gòu)時，延遲時間和嵌入維數(shù)是兩個主要參數(shù)，采用平均互信息法和假近鄰法來確定[10]，計算流程如圖4所示。求得延遲時間為30，嵌入維數(shù)為4，信號相軌跡在三維空間中的投影如圖5所示?？梢姡硐到y(tǒng)的穩(wěn)態(tài)相軌跡(吸引子)在重構(gòu)空間中完全展開，通過相軌跡隨槳葉損傷擴(kuò)展而出現(xiàn)的漂移和變形等幾何變化可直觀地了解和跟蹤旋翼損傷演化過程。

2 在重構(gòu)相空間跟蹤旋翼損傷的方法

筆者所提的旋翼損傷跟蹤方法首先在多個鄰域內(nèi)對狀態(tài)預(yù)測殘差統(tǒng)計平均形成一個損傷觀測特征向量，然后應(yīng)用奇異譜分解方法從損傷觀測特征向量的時間序列中提取出槳葉損傷演化的維度和趨勢信息，采用雙指數(shù)平滑(double exponential smoothing，簡稱DES)方法[11]建立損傷演化趨勢預(yù)測模型并估計出槳葉損傷故障剩余壽命。程序主要包括損傷觀測信號生成模塊和損傷狀態(tài)優(yōu)化估計模塊(如圖6)。損傷觀測信號生成模塊事先根據(jù)健康旋翼在相空間內(nèi)的運(yùn)動軌跡建立基準(zhǔn)狀態(tài)預(yù)測模型，運(yùn)行時利用現(xiàn)場采集數(shù)據(jù)和預(yù)測模型產(chǎn)生短時狀態(tài)預(yù)測殘差并統(tǒng)計生產(chǎn)損傷觀測特征向量。損傷分析和剩余壽命預(yù)測模塊采用奇異譜分解方法獲得旋翼槳葉損傷維數(shù)和演化趨勢信息，利用損傷演化趨勢預(yù)測模型估計槳葉剩余壽命。

圖6 旋翼損傷監(jiān)測程序原理圖Fig.6 Rotor health monitoring system framework

2.1 采用Volterra預(yù)測模型生成損傷觀測信號

在旋翼系統(tǒng)狀態(tài)空間中，用損傷旋翼和健康旋翼軌跡之間的短時偏差作為系統(tǒng)的損傷觀測信號。設(shè)旋翼在某任意時刻的狀態(tài)為xn，損傷為φ，經(jīng)過短時間τ后旋翼狀態(tài)變?yōu)閤(xn,τ,φ)。若按健康系統(tǒng)(φ=φR)估計其狀態(tài)，則應(yīng)為x(xn,τ,φR)，兩狀態(tài)之間的殘差為

S(xn,τ)(φ-φR)+O(ε)

(4)

短時狀態(tài)殘差利用系統(tǒng)的動力學(xué)模型去除信號中因運(yùn)行時間和系統(tǒng)狀態(tài)的正常轉(zhuǎn)移而造成的變化，保留了損傷造成的變化，是一種更好的損傷觀測信號。若能忽略高階小量的影響，式(4)可近似為

(5)

(6)

其中:k5(x0;φ)=(k1+2k2+2k3+k4)/6;k1=F(x0;φ)；k2=F(x0+0.5τk1;φ)；k3=F(x0+0.5τk2;φ)；k4=F(x0+τk3;φ)。

應(yīng)用式(6)計算得到短時狀態(tài)殘差后，就可研究靈敏度S(xn,τ)隨系統(tǒng)初始狀態(tài)和損傷大小的變化情況，具體分析高階小量的影響。下面分析實例旋翼槳尖揮舞位移的短時殘差靈敏度隨初始狀態(tài)和槳葉剛度損傷的變化情況。首先，在健康旋翼的狀態(tài)空間中隨機(jī)選擇16個點作為計算初始狀態(tài)，在槳葉壽命周期(0

本文提出了(N,m)維修模型,在充分考慮不同部件的維修優(yōu)先級不同的條件下，對高維修優(yōu)先級部件，采用修復(fù)性維修以及預(yù)防性維修；對低維修優(yōu)先級部件，采用成組維修。不同于傳統(tǒng)的依賴于單一團(tuán)隊的集中式維修模式，本文提出了異構(gòu)多部件系統(tǒng)的競爭式維修模式。以全周期的最低維修費率為目標(biāo)，通過計算不同部件的失效次序統(tǒng)計量，對參數(shù)進(jìn)行優(yōu)化，進(jìn)而得到最優(yōu)維修策略。通過數(shù)值算例驗證了維修策略的有效性，并對參數(shù)敏感度進(jìn)行了分析。

圖7 短時狀態(tài)殘差與壽命和損傷的關(guān)系Fig.7 Curves of short time state error vs. life time and damage variable

(7)

其中：y為傳感器輸出序列；x為相點；xR為預(yù)測器輸出的參考點；ex為狀態(tài)預(yù)測殘差；ek為損傷觀測信號。Pha為相空間重構(gòu)算子；m，τ分別為嵌入維數(shù)和時延參數(shù)；Bi為相空間鄰域(i=1,2，…,M)；Ni為鄰域內(nèi)相點數(shù)目；Vol為Volterra預(yù)測器；p為預(yù)測器階數(shù)，Volterra預(yù)測器的濾波系可根據(jù)健康旋翼的相軌跡采用自適應(yīng)算法確定[12]。

該方法的計算流程如圖8所示。根據(jù)健康旋翼在相空間內(nèi)的運(yùn)動軌跡，采用非線性Volterra級數(shù)建立基準(zhǔn)狀態(tài)短時預(yù)測模型。監(jiān)測時將旋翼的狀態(tài)采集數(shù)據(jù)運(yùn)用延遲嵌入方法重構(gòu)出實測相軌跡，采用Volterra預(yù)測模型計算出參考軌跡，將參考軌跡和實測相軌跡之差統(tǒng)計平均生成損傷觀測信號。由于Volterra預(yù)測模型在不同區(qū)域內(nèi)的精度變化，單個短時狀態(tài)殘差靈敏度具有較大的不確定性。為提高損傷檢測的魯棒性并增加損傷模式的跟蹤數(shù)目，將參考相空間劃分成靈敏度大體相等且互不重疊的幾個區(qū)域Bi(i=1,2，…,M)，狀態(tài)預(yù)測殘差對多種損傷類型的非線性特性可通過其在不同相空間區(qū)域內(nèi)的靈敏度差異來體現(xiàn)，用各個區(qū)域內(nèi)所有殘差的平均值來消除靈敏度矩陣的隨機(jī)誤差，將一條記錄中各區(qū)域內(nèi)殘差均值組合形成一個維數(shù)足夠高的損傷觀測特征向量ek。

圖8 損傷觀測信號的計算流程Fig.8 Flow chart of damage measurement prediction

應(yīng)用該方法對實例旋翼在不同損傷狀態(tài)下的槳葉位移和槳轂載荷進(jìn)行分析。在槳葉壽命周期內(nèi)(0≤t/tF<1)，等時間間隔采集60個損傷程度(0≤φ<1)下的仿真記錄，將參考相空間按照坐標(biāo)分量x1大小劃分成8個等間隔的鄰域，在各鄰域內(nèi)統(tǒng)計每條記錄的短時狀態(tài)殘差，計算所得每個殘差特征向量ek包含32個列，各殘差特征向量隨使用壽命的變化曲線如圖9所示?？梢?，各監(jiān)測信號的殘差特征向量可大致反映槳葉損傷曲線的變化趨勢，但數(shù)據(jù)依然有較大起伏，難以判斷其中含有幾種損傷發(fā)展模態(tài)。

圖9 殘差特征向量隨時間的變化Fig.9 Curves of residual tracking metric(Scaled)

2.2 槳葉損傷信息的提取

理論上可以利用殘差特征向量ek與損傷變量之間的函數(shù)關(guān)系φk=H(ek)來估計槳葉損傷。然而，旋翼的實際損傷情況比較復(fù)雜，損傷產(chǎn)生的原因可能來自疲勞、腐蝕和照射等不確定性因素，損傷變量的實際維數(shù)往往不能事先知道；有些損傷類型不但損傷擴(kuò)展規(guī)律不確定，甚至沒有事先可直接測量的損傷變量，難以建立損傷與觀測之間映射關(guān)系。因此，需要直接應(yīng)用ek序列分析損傷變量的維數(shù)和損傷演化的模態(tài)信息。奇異譜分析方法[13]能夠研究多維時序數(shù)據(jù)在模態(tài)空間中的能量分布情況，可用于分析旋翼的損傷維數(shù)和演化趨勢信息。

(8)

式(8)是一種Rayleigh商的優(yōu)化問題[13]，可化為廣義特征值問題

(9)

其中：

Δt為序列的時間間隔；D為微分算子矩陣，大小為(n-1)n；λ為特征值；q為特征向量。

在Matlab軟件中廣義特征值方程(8)可采用奇異譜分解法求解：(Q,λ) =eig(YTY, (DY)TDY)。矩陣Y在模態(tài)坐標(biāo)系中的投影分量為

φk=Yqk

(10)

其中：qk為第k個廣義模態(tài)向量；φk為第k個模態(tài)投影。

矩陣Y中既包含損傷變量緩慢演化信號，又包含旋翼系統(tǒng)的狀態(tài)變量快速運(yùn)動的信號。如果殘差矩陣中有k個特征值比其它特征值高出一個數(shù)量級，則表明系統(tǒng)存在著k個緩慢的損傷演化模式，這些特征值對應(yīng)的投影可指示損傷演化過程。

采用奇異譜分解算法對BO-105旋翼各監(jiān)測信號的殘差矩陣Y進(jìn)行分析，得到的32個廣義特征值如圖10(a)所示，相應(yīng)的，各殘差矩陣在32個模態(tài)坐標(biāo)系中的投影分量如圖10(b)所示?？梢?，殘差信號的能量按模態(tài)頻率重新分配后，有一個特征值遠(yuǎn)遠(yuǎn)大于其他特征值，說明殘差矩陣中有一種確定的趨勢分量。與之對應(yīng)模態(tài)的投影分量較為清晰地反映槳葉損傷擴(kuò)展的趨勢信息，容易驗證槳葉損傷與其之間有較為顯著的線性關(guān)系。將該模態(tài)投影分量作為一種名義損傷變量，可對旋翼的損傷狀態(tài)進(jìn)行定量跟蹤。

圖10 廣義特征值和模態(tài)坐標(biāo)投影Fig.10 Eigenvalues and modal projections

2.3 槳葉剩余壽命的預(yù)測

為估計槳葉剩余壽命需要建立合理的損傷趨勢預(yù)測模型。槳葉剛度損傷過程是一種非線性隨機(jī)過程，雙指數(shù)平滑方法[11]既能靈活適應(yīng)回歸變量的非線性動態(tài)特性，又能快速計算多變量回歸預(yù)測模型的統(tǒng)計分析參數(shù)，是一種優(yōu)選的自適應(yīng)預(yù)測技術(shù)。DES方法的基本預(yù)測公式為

(11)

當(dāng)給定失效條件φF后，根據(jù)式(10)求得槳葉剩余壽命。采用DES方法對實例旋翼各殘差矩陣的模態(tài)投影分量進(jìn)行分析，得到槳葉剩余壽命估計結(jié)果如圖11所示。在使用時間t≥0.2tF區(qū)間內(nèi)，采用各特征模態(tài)投影分量對槳葉剩余壽命跟蹤的誤差如表1所示。可見，這種方法能較為有效地跟蹤槳葉壽命過程。

圖11 應(yīng)用不同觀測數(shù)據(jù)對剩余壽命估計Fig.11 Prediction of remaining time to failure using different measurements

表1 槳葉剩余壽命的跟蹤誤差

3 結(jié)束語

對結(jié)構(gòu)損傷機(jī)理和材料失效理論的研究通常都是在微觀尺度上進(jìn)行的。這種研究是必要的，但這類研究不能闡明這些損傷過程在整個系統(tǒng)宏觀行為中是如何表現(xiàn)的，也不能提供一種系統(tǒng)級的損傷跟蹤手段。源自混沌信號分析方法的相空間重構(gòu)技術(shù)可以從不同的時空角度觀察系統(tǒng)的損傷行為特性。筆者利用信號重構(gòu)技術(shù)和Volterra預(yù)測模型發(fā)展了一種新的損傷觀測殘差生成方法，對損傷觀測殘差序列分析可提取損傷變量的維度和趨勢信息。所提方法能充分利用被監(jiān)測信號的非線性特性在高維相空間中重建系統(tǒng)動力學(xué)本質(zhì)，以不同的時間尺度來觀測系統(tǒng)的演化特性。該方法具有損傷模式的自動識別能力，可用于損傷演化模型或維數(shù)信息難以事先確定的場合。

[1] Azzam H, Andrew M J.The use of math-dynamic model to aid the development of integrated health and usage monitoring[J].Journal of Aerospace Engineering,1992,206(1):71-96.

[2] Ganguli R,Chopra I,Haas D J.Detection of helicopter rotor system simulated faults using neural networks[J].Journal of American Helicopter Society,1997,42(2):161-171.

[3] Pawar P M, Ganguli R.Fuzzy logic based health monitoring and residual life prediction for composite rotor[J].Journal of Aircraft,2007,44(3): 981-995.

[4] Pawar P M, Ganguli R. Helicopter rotor health monitoring-a review[J]. Journal of Aerospace Engineering,2007,221(1): 631-647.

[5] Johnson W. Comprehensive analytical model for rotorcraft aerodynamics and dynamics volumes ii:components theory [M].Palo Alto:Johnson Aeronautics,2007:280-283.

[6] Miles T J, Lucas M.Torsional and bending vibration measurements on rotors using laser technology[J]. Journal of Sound Vibration, 1999,266(3):1441-1467.

[7] Pawar P M,Ganguli R. On the effect of progressive damage on composite helicopter rotor system behavior[J]. Composite Structures, 2007,78:410-423.

[8] Chelidze D, Liu M.Dynamical systems approach to fatigue damage identification [J]. Journal of Sound Vibration,2005, 281(3):887-904.

[9] Hilborn R C.Chaos and nonlinear dynamics[M].Oxford:Oxford University Press,2006:375-389.

[10]Cao Liangyue.Practical method for determining the minimum embedding dimension of a scalar time series[J]. Physica D:Nonlinear Phenomena,1997,110(1-2):43-50.

[11]Guo Zhihao, Shahdi M, Shaya S. Multi-objective OLSR for proactive routing in manet with delay,energy, and link lifetime predictions[J].Applied Mathematical Modeling,2011,35:1413-1426.

[12]張家樹.混沌信號的非線性自適應(yīng)預(yù)測技術(shù)及其應(yīng)用研究[D].成都：中國電子科技大學(xué), 2001.

[13]錢征文,程禮,李應(yīng)紅.利用奇異值分解的信號降噪方法[J].振動、測試與診斷,2011，31(4):459-463.

Qian Zhengwen,Cheng Li,Li Yinghong. Noise reduction method based on singular value decomposition[J].Journal of Vibration,Measurement & Diagnosis, 2011,31(4):459-463.(in Chinese)

[14]Marquardt D. An algorithm for least-squares estimation of nonlinear parameters [J]. SIAM Journal on Applied Mathematics, 1963, 11: 431-441．

歡迎訂閱《振動、測試與診斷》

《振動、測試與診斷》由工業(yè)和信息化部主管，南京航空航天大學(xué)和全國高校機(jī)械工程測試技術(shù)研究會聯(lián)合主辦，是反映振動、動態(tài)測試及故障診斷學(xué)科領(lǐng)域的科研成果及其應(yīng)用情況的技術(shù)性刊物。主要刊登國內(nèi)外以振動測試與故障診斷為中心的動態(tài)測試?yán)碚?、方法和手段的研究及?yīng)用方面的技術(shù)文獻(xiàn)，包括實驗測試技術(shù)、測試儀器的研制、方法和系統(tǒng)組成、信號分析、數(shù)據(jù)處理、參數(shù)識別與故障診斷以及有關(guān)裝置的設(shè)計、使用、控制、標(biāo)定和校準(zhǔn)等，不拘泥于行業(yè)和測試項目。

本刊為EICompendex數(shù)據(jù)庫收錄期刊和中文核心期刊，雙月刊，每逢雙月末出版，每本定價20元，全年120元。歡迎訂閱和投稿，歡迎在本刊刊登各類廣告和科技信息。

編輯部地址：南京市御道街29號 郵政編碼：210016 電話：(025)84893332

傳真：(025)84893332 E-mail:qchen@nuaa.edu.cn 網(wǎng)址：http:∥zdcs.nuaa.edu.cn

*國家自然科學(xué)基金資助項目(51375490)

2013-01-07；

2013-02-24

10.16450/j.cnki.issn.1004-6801.2015.01.020

V275.1; TH113.1

朱旭程,男，1973年11月，博士、副教授。主要研究方向為航空機(jī)械維修、綜合保障。曾發(fā)表《直升機(jī)旋翼健康監(jiān)測系統(tǒng)地面試驗臺方案設(shè)計與開發(fā)》(《兵工學(xué)報》2008年第29卷第4期)等論文。 E-mail：zhu133@189.cn