洪 明 雷 川 崔洪宇

（大連理工大學 船舶工程學院 大連116024）

環(huán)境激勵下船舶結(jié)構(gòu)模態(tài)參數(shù)識別方法綜述

洪 明 雷 川 崔洪宇

（大連理工大學 船舶工程學院 大連116024）

基于模態(tài)分析理論，對環(huán)境激勵下的模態(tài)參數(shù)識別方法進行歸納總結(jié)，討論了不同模態(tài)參數(shù)識別方法的基本思想和優(yōu)缺點，并著重介紹船舶航行狀態(tài)下的模態(tài)參數(shù)識別。文章對國內(nèi)模態(tài)參數(shù)識別的研究熱點問題進行了闡述，為相關(guān)科研人員研究環(huán)境激勵下模態(tài)參數(shù)識別提供了有益參考。

環(huán)境激勵；模態(tài)參數(shù)識別；船舶結(jié)構(gòu)；運行模態(tài)分析

引 言

船舶航行過程中，結(jié)構(gòu)會受到來自螺旋槳、主機及風浪流等激勵作用，不可避免地會產(chǎn)生結(jié)構(gòu)振動。船舶振動不僅會影響船上人員生活的舒適性，還會對船上設(shè)備及船舶結(jié)構(gòu)自身強度帶來危害，因此船舶振動問題逐漸引起人們的重視。

一般的振動問題由激勵（輸入）、結(jié)構(gòu)（系統(tǒng)）和響應(yīng)（輸出）三部分組成，根據(jù)研究目的不同，通常將振動問題分為圖1所示三種基本類型。

圖1 振動問題分類

振動結(jié)構(gòu)模型可以用三種形式來描述，分別是物理參數(shù)模型、模態(tài)參數(shù)模型和非參數(shù)模型［1］。解決振動問題時最為關(guān)注的是模態(tài)參數(shù)模型，它包括結(jié)構(gòu)的模態(tài)頻率、振型和阻尼比等參數(shù)，而模態(tài)分析理論就是以振動理論為基礎(chǔ)，以模態(tài)參數(shù)為目標的分析方法。根據(jù)模態(tài)分析方法的不同，模態(tài)分析理論可以分為理論模態(tài)分析、實驗模態(tài)分析（EMA）和運行模態(tài)分析（OMA）［2］，由于運行模態(tài)分析的應(yīng)用特點及實施理論難度，近年來基于環(huán)境激勵下的模態(tài)分析方法越來越受到青睞，以下本文主要討論運行模態(tài)。

1 模態(tài)參數(shù)識別方法

傳統(tǒng)的模態(tài)參數(shù)識別方法，通過在結(jié)構(gòu)上人工施加激勵，已知激勵和響應(yīng)，采用信號處理技術(shù)計算系統(tǒng)的頻響函數(shù)或脈沖響應(yīng)函數(shù)，再運用參數(shù)識別方法得到系統(tǒng)的模態(tài)參數(shù)。傳統(tǒng)的模態(tài)分析方法按識別域的不同可分為頻域法和時域法，頻域法是以頻響函數(shù)為基礎(chǔ)的參數(shù)識別方法，時域法則是以脈沖響應(yīng)函數(shù)或自由振動響應(yīng)為基礎(chǔ)的參數(shù)識別方法；按處理各階耦合模態(tài)所用方法的不同，可分為單自由度法和多自由度法，單自由度法基本思想是將具有多階模態(tài)的多自由度系統(tǒng)視作多個單自由度系統(tǒng)進行識別，適用于阻尼小且相鄰模態(tài)分離較遠的系統(tǒng)，多自由度法則是對頻響函數(shù)曲線上的各共振峰同時擬合，在感興趣的頻帶內(nèi)考慮各階模態(tài)的耦合與相互影響；按照使用激勵和響應(yīng)的數(shù)目不同，又可以分為單輸入單輸出（SISO）、單輸入多輸出（SIMO）及多輸入多輸出（MIMO）識別方法，其中SISO屬于局部識別方法，SIMO和MIMO屬于整體識別方法［3］。

以船舶結(jié)構(gòu)為代表的大型結(jié)構(gòu)物傳統(tǒng)實驗模態(tài)分析方法需要同時已知激勵與響應(yīng)，才能通過構(gòu)造頻響函數(shù)對系統(tǒng)參數(shù)進行識別，然而實驗過程中很難對船體結(jié)構(gòu)施加有效激勵。激勵過小，能量不足以激起船舶結(jié)構(gòu)的振動；激勵過大，容易導致局部結(jié)構(gòu)發(fā)生塑性變形等強度問題。由此，基于環(huán)境激勵的運行模態(tài)分析方法顯現(xiàn)出優(yōu)勢，僅需要船舶在航行狀態(tài)測量的振動響應(yīng)信號，便可識別系統(tǒng)的模態(tài)參數(shù)?；诃h(huán)境激勵的運行模態(tài)分析技術(shù)，一方面可識別船體的振動模態(tài)，判斷其是否與船上激勵源發(fā)生共振，還可與設(shè)計預報結(jié)果進行對比，積累預報經(jīng)驗；另一方面可對振動量級過高區(qū)域進行分析，得到該區(qū)域有害振動的主要頻率成分，對比激勵源頻率可以確定有害振動來源，并可以根據(jù)分析得到的模態(tài)振型，對振幅較大位置采取減振措施［4］。

2 環(huán)境激勵下模態(tài)參數(shù)識別方法概述

基于環(huán)境激勵模態(tài)參數(shù)識別方法，按照信號識別域可分為頻域法、時域法以及聯(lián)合時頻域法［5］。頻域法大多利用經(jīng)典譜估計，對于輸入信號已知的實驗狀態(tài)進行頻域模態(tài)分析，根據(jù)頻響函數(shù)在系統(tǒng)固有頻率峰值特征進行模態(tài)識別；基于環(huán)境激勵下的模態(tài)分析，由于無法直接求得頻響函數(shù)，所以一般使用功率譜密度函數(shù)近似地代替頻響函數(shù)，基于功率譜密度函數(shù)進行模態(tài)識別的方法主要有峰值拾取法和頻域分解法。

時域法屬于現(xiàn)代譜分析的范疇，其算法思想大致有兩種：第一種直接利用響應(yīng)信號求系統(tǒng)模態(tài)參數(shù)，如基于數(shù)據(jù)的隨機子空間法以及ARMA時序分析法；第二種先需要對信號進行處理，配合自然激勵技術(shù)或隨機減量技術(shù)，得到中間時間序列如相關(guān)函數(shù)、近似的自由響應(yīng)或脈沖響應(yīng)信號等，然后利用經(jīng)典時域模態(tài)識別方法進行參數(shù)識別，如ITD法、最小二乘復指數(shù)法、多參考點復指數(shù)法、基于協(xié)方差的隨機子空間法以及特征系統(tǒng)實現(xiàn)法等［6］。

無論頻域法還是時域法都是基于傅里葉變換進行模態(tài)識別的，但由于傅里葉變換自身不能反映信號的瞬時性即信號在時頻域內(nèi)的變化規(guī)律，因此能夠在時頻域內(nèi)反映信號能量分布情況的時頻分析開始成為學者研究的熱點，其中小波變換和Hilbert-Huang變換是最為經(jīng)典的兩種時頻分析方法，它們都能夠根據(jù)信號的時間-尺度特征自適應(yīng)地對系統(tǒng)進行模態(tài)識別，彌補了頻域法和時域法在處理非穩(wěn)態(tài)環(huán)境信號和非線性系統(tǒng)模態(tài)識別的不足，擴展了模態(tài)參數(shù)識別的應(yīng)用領(lǐng)域［7］。圖2所示為常見的基于環(huán)境激勵的模態(tài)參數(shù)識別方法以及模態(tài)識別過程中應(yīng)用到的關(guān)鍵技術(shù)，本文將給出這些算法的數(shù)學模型和算法思想，并重點闡述這些算法的適用性和優(yōu)缺點。

圖2 環(huán)境激勵下模態(tài)參數(shù)識別方法分類

2.1 頻域法

2.1.1 峰值拾取法

峰值拾取法［8］的基本思想是系統(tǒng)的頻響函數(shù)在其固有頻率處會出現(xiàn)峰值點，由于環(huán)境激勵前提為白噪聲信號，白噪聲信號的功率譜密度函數(shù)在一定頻率范圍內(nèi)近似均勻分布，因此可以利用輸出響應(yīng)的功率譜密度函數(shù)近似代替頻響函數(shù)。另外該方法假定隨機響應(yīng)的功率譜密度函數(shù)峰值僅由一個模態(tài)確定，這樣系統(tǒng)的固有頻率可以由功率譜密度函數(shù)的峰值得到，利用識別結(jié)果撓度分布就可近似表示模態(tài)振型。峰值拾取法的優(yōu)點是算法簡單且識別效率高；缺點是環(huán)境激勵需滿足白噪聲特征，且不能識別密集模態(tài)和阻尼比，由于對功率譜密度函數(shù)峰值選取存在主觀性，所以其識別精度受到影響。

2.1.2 頻域分解法

頻域分解法是峰值法的改進算法，主要用來解決峰值法難以處理密集模態(tài)的問題。頻域分解法對功率譜密度函數(shù)矩陣進行奇異值分解（SVD），將系統(tǒng)響應(yīng)的功率譜密度函數(shù)分解為對應(yīng)多階模態(tài)的一系列單自由度系統(tǒng)功率譜密度函數(shù)，然后利用峰值法識別模態(tài)參數(shù)。頻域分解法是峰值法的一種延伸，它不僅繼承了峰值法算法簡單、識別效率高的優(yōu)點，還能識別密集結(jié)構(gòu)的模態(tài)參數(shù)，而且由于采用了SVD技術(shù)，所以對噪聲也有一定的抗干擾能力，識別精度較峰值法有所提高。

張毅剛等［9］針對峰值法進行空間網(wǎng)格結(jié)構(gòu)模態(tài)參數(shù)識別時會遇到遺漏模態(tài)及重復頻率無法篩選的問題，對峰值法進行了改進，提出在獲得所有測點平均正則化功率譜基礎(chǔ)上，結(jié)合結(jié)構(gòu)的理論振型特點計算輔助正則化功率譜，取兩者并集為模態(tài)識別結(jié)果，通過典型結(jié)構(gòu)監(jiān)測結(jié)果驗證了改進后的功率譜峰值法的有效性，圖3為實驗測得的平均正則化及輔助正則化功率譜密度函數(shù)曲線［9］。

圖3 平均正則化及輔助正則化功率譜密度函數(shù)曲線

夏祥麟［10］對簡支梁分別進行了峰值拾取法和頻域分解法研究，發(fā)現(xiàn)頻域分解法在保留了峰值法識別快速直觀的基礎(chǔ)上，摒除了峰值由單一模態(tài)貢獻的假定條件，具有一定的抗噪性能，但是對于非白噪聲信號及大阻尼系統(tǒng)，頻域分解法的識別效果并不準確；Brincker等［11］在頻域分解法的基礎(chǔ)上對算法進行改進，其基本思想是將分解后的單自由度功率譜密度函數(shù)進行逆傅里葉變換，在時域內(nèi)計算得到相關(guān)函數(shù)后，利用指數(shù)衰減法計算頻率和阻尼比，改進后的頻域分解法可直接在時域內(nèi)完成頻率和阻尼比的計算，使得對頻率的求解不再依賴FFT分辨率的影響。

2.2 時域法

2.2.1 基于原始響應(yīng)信號的時域法

（1）ARMA時序分析法

1969年Akaike［12］首次使用自回歸滑動平均模型進行白噪聲環(huán)境激勵下的模態(tài)參數(shù)識別。ARMA時序分析法的基本思想是利用差分方程和Z變換，分別建立強迫振動方程傳遞函數(shù)與ARMA模型間等價關(guān)系，直接使用隨機激勵和響應(yīng)信號，并采用ARMA模型結(jié)合最小二乘法來識別系統(tǒng)模態(tài)參數(shù)。ARMA時序分析法屬于SISO參數(shù)識別法，用一個測點就可以識別出系統(tǒng)各階模態(tài)參數(shù)。ARMA參數(shù)模型包括AR自回歸模型、MA滑動平均模型和ARMA自回歸滑動平均模型。

對確定性系統(tǒng)，系統(tǒng)輸入和輸出關(guān)系為：

該式稱為系統(tǒng)的自回歸滑動平均模型。式中：ai和bi分別稱為自回歸系數(shù)和滑動平均系數(shù)，如果bi= 0，則上式變?yōu)榧碅R自回歸模型；如果ai= 0，則即MA滑動平均模型，所以AR模型和MA模型是ARMA模型的兩種特殊情形。

ARMA時序分析法不受觀測樣本長度限制，不存在能量泄露的問題，且識別精度較高；缺點是在建立時序模型時系統(tǒng)的階次難以確定。郭永剛等［13］利用單位脈沖響應(yīng)函數(shù)與ARMA模型Green函數(shù)等價的特點，通過脈沖響應(yīng)函數(shù)來估計系統(tǒng)響應(yīng)的自相關(guān)系數(shù)；然后建立推廣的Yule-Walker方程以求得ARMA模型自回歸系數(shù)并進行參數(shù)識別，在其算例研究中發(fā)現(xiàn)為保證識別結(jié)果的精度需建立較高階數(shù)的ARMA模型，但模型階數(shù)過高又會導致識別計算量劇增；Smail等［14］針對基于ARMA模型的模態(tài)識別方法系統(tǒng)階次難以確定的問題，提出利用輸出信號相關(guān)矩陣特征值個數(shù)來確定系統(tǒng)階次，取得了較好的識別結(jié)果。

（2）數(shù)據(jù)驅(qū)動的隨機子空間法

數(shù)據(jù)驅(qū)動的隨機子空間法是由Overschee等于1991年提出的［15］，算法直接將原始響應(yīng)數(shù)據(jù)作為輸入部分形成Hankel矩陣，利用QR分解和SVD技術(shù)獲得擴展的可觀測矩陣，并利用卡爾曼濾波狀態(tài)序列結(jié)合最小二乘法識別出系統(tǒng)的模態(tài)參數(shù)。

數(shù)據(jù)驅(qū)動的隨機子空間法相對于協(xié)方差驅(qū)動的隨機子空間法最大區(qū)別有兩點：一是不需要進行相關(guān)函數(shù)的計算，提高了算法的識別效率；二是采用QR分解技術(shù)，研究表明QR分解能夠很大程度上提高隨機子空間法的識別精度。辛峻峰等［16］在理論上探討了協(xié)方差驅(qū)動和數(shù)據(jù)驅(qū)動兩種隨機子空間法的不同，并通過相應(yīng)的數(shù)值模擬驗證基于QR分解的數(shù)據(jù)驅(qū)動隨機子空間法無論計算精度或?qū)^弱勢模態(tài)的識別能力均明顯優(yōu)于協(xié)方差驅(qū)動隨機子空間法，辛峻峰等［17］還對隨機子空間法深入分析噪聲與數(shù)據(jù)驅(qū)動的隨機子空間法Hankel矩陣維數(shù)之間的關(guān)系，并提出一種評估隨機子空間法矩陣維數(shù)選擇優(yōu)劣的方法；常軍等［18］就隨機子空間法中存在的虛假模態(tài)問題進行分析，表明產(chǎn)生虛假模態(tài)的原因主要有兩方面：一方面是由于隨機子空間法的算法本身而導致；另一方面是由于輸入信號不滿足白噪聲假定或者輸出信號受到環(huán)境干擾而引起；章國穩(wěn)等［19］針對數(shù)據(jù)驅(qū)動的隨機子空間法計算效率低下的問題，提出一種基于特征值分解的隨機子空間法，通過模型研究證明該方法在保持計算精度的基礎(chǔ)上提高了隨機子空間法的計算效率。數(shù)據(jù)驅(qū)動的隨機子空間法適用于線性結(jié)構(gòu)平穩(wěn)激勵下的模態(tài)參數(shù)識別，對輸出噪聲也有一定的抗干擾能力，但仍不可避免地會遇到虛假模態(tài)的問題；此外，狀態(tài)空間方程中系統(tǒng)階次的確定是限制所有基于狀態(tài)空間方程方法運算速度和精度的最主要因素，基于狀態(tài)空間方程的隨機子空間法也存在同樣的問題。

2.2.2 基于自由響應(yīng)或脈沖響應(yīng)信號的時域法

（1）隨機減量技術(shù)（RDT）

對于船舶等大型結(jié)構(gòu)物來說，工程實際中容易得到的是隨機響應(yīng)信號，但是很多時域模態(tài)識別方法是以自由振動響應(yīng)為數(shù)學模型的，所以需要從隨機響應(yīng)信號中提取出某種自由振動響應(yīng)信號。RDT法是利用樣本平均的方法，通過設(shè)置一定的觸發(fā)條件去除響應(yīng)中的隨機成分，獲得初始激勵下的自由振動響應(yīng)，然后通過相應(yīng)的時域法識別出系統(tǒng)的模態(tài)參數(shù)。Cole［20］首先提出RDT法，并將其成功應(yīng)用于空間飛行器模型結(jié)構(gòu)的模態(tài)識別；Ibrahim［21］將RDT法同ITD法相結(jié)合用于結(jié)構(gòu)運行模態(tài)參數(shù)的識別，并對RDT法中采用的觸發(fā)條件進行詳細闡述，同時在其論文中強調(diào)RDT法僅適用于白噪聲激勵。由于RDT法建立在線性系統(tǒng)疊加原理的基礎(chǔ)上，所以不能用來分析非線性系統(tǒng)；另外RDT法利用了白噪聲激勵下系統(tǒng)產(chǎn)生的平穩(wěn)響應(yīng)方差為零這一特性，所以RDT法僅能用于處理平穩(wěn)隨機信號。

（2）自然激勵技術(shù)（NExT）

NExT法的基本思想是白噪聲環(huán)境激勵下結(jié)構(gòu)兩點間響應(yīng)的互相關(guān)函數(shù)和脈沖響應(yīng)函數(shù)有近似的表達式，求得兩點間響應(yīng)的互相關(guān)函數(shù)后，將其作為輸入信號進行模態(tài)參數(shù)識別。美國SADIA國家實驗室的James等［22］在1994年提出NExT法，并且將該方法運用于汽輪機葉片在工作狀態(tài)下的固有頻率和模態(tài)阻尼的測試。

對自由度為n的線性系統(tǒng)，當系統(tǒng)k點受脈沖激勵時，i點的脈沖響應(yīng)寫成如式（2）所示：

從上式可以發(fā)現(xiàn)，i和j兩點之間的互相關(guān)函數(shù)可看成一系列衰減簡諧函數(shù)的組合，同兩點間脈沖響應(yīng)函數(shù)具有相同表達形式，故可將NExT法與基于脈沖響應(yīng)時域方法結(jié)合起來進行模態(tài)參數(shù)識別。

羅奎［23］分別用NExT結(jié)合ITD法以及NExT結(jié)合ERA法對環(huán)境激勵下三跨混凝土連續(xù)橋進行振動試驗模態(tài)參數(shù)識別，發(fā)現(xiàn)NExT法由于對響應(yīng)信號進行互相關(guān)函數(shù)的計算，其對噪聲的抗干擾能力較強；紀曉東等［24］綜合自然激勵技術(shù)和特征系統(tǒng)實現(xiàn)算法，進行模擬環(huán)境激勵下結(jié)構(gòu)的時域模態(tài)參數(shù)識別，通過改變結(jié)構(gòu)質(zhì)量和特征靈敏度分析，得到質(zhì)量歸一化振型，通過模型實驗驗證了NExTERA法能應(yīng)用于結(jié)構(gòu)在線模態(tài)參數(shù)識別。

（3）ITD法

ITD法［25］是Ibrahim提出的一種SIMO參數(shù)識別方法，基本思想是直接使用各測點的自由響應(yīng)信號，通過三次不同延時采樣，構(gòu)造自由響應(yīng)采樣數(shù)據(jù)的增廣矩陣，然后建立特征方程，求解出特征對后再估算各階模態(tài)參數(shù)，ITD法的最大特點是同時使用全部測點的自由響應(yīng)數(shù)據(jù)。Ibrahim在ITD法的基礎(chǔ)上提出省時的STD法，STD法相對ITD法的計算量大為降低，節(jié)省內(nèi)存和運算時間，在識別精度上也有所提高。

楊佑發(fā)等［26］結(jié)合隨機子空間法對ITD法進行算法改進，提出將隨機子空間法中計算得到Toeplitz矩陣作為輸入傳遞給ITD法，避免使用RDT法和NExT法進行前處理引起的誤差，實驗表明該改進后ITD算法相對于隨機子空間法識別精度沒有降低，而且縮短了計算時間，相對于傳統(tǒng)的ITD法，識別精度明顯有所提高；周凱［27］采用ITD法對環(huán)境載荷下一個近海導管架平臺結(jié)構(gòu)模型進行模態(tài)參數(shù)識別，研究發(fā)現(xiàn)增多ITD法中虛擬測點的個數(shù)，增大特征矩陣的階數(shù)，不僅能增加識別出的模態(tài)數(shù)量，而且有助于降低其對噪聲的敏感度，虛擬測點個數(shù)對系統(tǒng)模態(tài)識別結(jié)果影響如圖4所示［27］。

圖4 虛擬測點個數(shù)對系統(tǒng)模態(tài)識別結(jié)果影響

（4）最小二乘復指數(shù)法（LSCE）和多參考點復指數(shù)法（PRCE）

LSCE法［28］又稱Prony多項式法，與ITD法相比，LSCE法在識別模態(tài)頻率和阻尼比時只用一個測點的脈沖響應(yīng)數(shù)據(jù)，因而LSCE法屬于局部識別法。LSCE法的基本思想是以Z變換因子表示脈沖響應(yīng)，Z變換因子中包含待識別的復頻率，構(gòu)造Prony多項式，使其零點等于Z變換因子的值，將求解Z變換因子就轉(zhuǎn)化為求解Prony多項式系數(shù)：

為求解Prony多項式系數(shù)，構(gòu)造脈沖響應(yīng)數(shù)據(jù)序列的自回歸模型，自回歸系數(shù)即Prony多項式系數(shù)，通過在不同起始點采樣，得到關(guān)于自回歸系數(shù)的線性方程組，用最小二乘法計算得到自回歸系數(shù)的解，可得Prony多項式的根。再由脈沖響應(yīng)數(shù)據(jù)序列構(gòu)造該測點各階脈沖響應(yīng)幅值（留數(shù)）的線性方程組，用最小二乘法求解，對各測點作識別，便可得到各階模態(tài)矢量，脈沖響應(yīng)序列的自回歸模型如下：

式中：為以l為起始采樣點的含2n個采樣值的時序樣本，為自回歸系數(shù)列陣，即Prony多項式的系數(shù)。

為克服LSCE法僅能應(yīng)用于單輸出狀態(tài)的不足，在LSCE法的基礎(chǔ)上建立MIMO整體識別的PRCE法［29］，其數(shù)學模型為基于MIMO的脈沖響應(yīng)函數(shù)矩陣。

（5）協(xié)方差驅(qū)動的隨機子空間法

其基本思想是由輸出響應(yīng)數(shù)據(jù)構(gòu)建Hankel矩陣，計算協(xié)方差序列形成Toeplitz矩陣，Toeplitz矩陣主要作用是保持信號原有信息的情況下縮減數(shù)據(jù)量，然后對其作SVD得到系統(tǒng)矩陣，進而識別結(jié)構(gòu)的模態(tài)參數(shù)。其同樣存在系統(tǒng)階次難以確定的問題，無論是根據(jù)非零奇異值點個數(shù)或是利用奇異值的跳躍都難以在實際工程中確定系統(tǒng)的階次，目前一般采用隨機子空間結(jié)合穩(wěn)定圖的方法進行模態(tài)參數(shù)的識別［30］，穩(wěn)定圖是目前普遍認可的相對較好的確定系統(tǒng)階次的方法［31］。

（6）協(xié)方差驅(qū)動的隨機子空間法特征系統(tǒng)實現(xiàn)法（ERA）

特征系統(tǒng)實現(xiàn)算法是美國航天航空局的Langley研究部門于1984年提出的一種時域模態(tài)分析方法［32］。ERA法以由MIMO得到的脈沖響應(yīng)函數(shù)為基本模型，利用結(jié)構(gòu)的脈沖響應(yīng)信號函數(shù)矩陣來構(gòu)建廣義Hankel矩陣，對其進行SVD，從非零奇異值中提取得到系統(tǒng)的特征量，通過對系統(tǒng)矩陣進行特征值分解識別出系統(tǒng)的模態(tài)參數(shù)。因其在計算中利用了最小實現(xiàn)原理，大幅度減小矩陣維數(shù)和計算量，在實際工程中應(yīng)用十分廣泛。

ERA法通常結(jié)合NExT法用來識別環(huán)境激勵（白噪聲激勵）下的系統(tǒng)模態(tài)參數(shù)，NExT-ERA法程序流程圖如下頁圖5所示。許俊臣和洪明［33］利用NExT-ERA法對一鋼質(zhì)船體梁模型進行了模態(tài)參數(shù)識別，如圖6 -圖8［33］所示。

圖5 NExT-ERA法程序流程圖

圖6 船模結(jié)構(gòu)示意圖

圖7 水池船模結(jié)構(gòu)實驗圖

圖8 穩(wěn)態(tài)白噪聲激勵下的系統(tǒng)穩(wěn)定圖

研究中對存在單頻激勵時對應(yīng)的簡諧頻率的影響，信噪比的影響以及在不同裝載情況下的船模模態(tài)參數(shù)進行大量的實驗對比和系統(tǒng)的分析研究。實驗結(jié)果顯示只要簡諧激勵的頻率與固有頻率沒有離得太近，其頻率和信噪比對模態(tài)分析的影響不大，當兩者頻率較為接近，可能會發(fā)生頻率混疊的現(xiàn)象；局部結(jié)構(gòu)上附加質(zhì)量的存在對局部模態(tài)振型和識別結(jié)果會造成一定影響；萬嶺和洪明等人［34］采用NExT-ERA法對簡支梁模型和船體模型分別進行模態(tài)參數(shù)識別，驗證了NExT-ERA法進行運行模態(tài)識別能夠得到良好效果。文章對ERA法中系統(tǒng)階次的確定進行深入討論，指出選擇較低階次容易丟失系統(tǒng)的模態(tài)信息，而選擇較高階次又會引入噪聲模態(tài)；楊和振、李華軍等人［35］采用NExT-ERA法對海洋環(huán)境載荷下海洋平臺結(jié)構(gòu)模擬模型進行了參數(shù)識別，并將識別結(jié)果同峰值法進行比較，比較結(jié)果顯示NExT-ERA法更加適用于環(huán)境激勵下大型結(jié)構(gòu)物的模態(tài)參數(shù)識別。

2.3 聯(lián)合時頻域法

2.3.1 小波變換法

小波變換是一種基于信號時間-尺度的分析方法［36］，具有多分辨率分析的特點，即在低頻部分具有較高的頻率分辨率和較低的時間分辨率，在高頻部分具有較高的時間分辨率和較低的頻率分辨率，在時域和頻域都具有表征信號局部特征的能力，很適合探測正常信號中夾帶的瞬態(tài)反常現(xiàn)象，因此被廣泛地應(yīng)用于對非平穩(wěn)信號的分析和處理中［37］。

對于一個自由度為n的線性系統(tǒng)，其自由振動響應(yīng)可以表示為：

對上式進行小波變換可以得到：

式中：a和b分別稱為小波變換的尺度參數(shù)與位移參數(shù)。

當尺度參數(shù)為ai時，在小波變換過程中，只有與ai相關(guān)的模態(tài)起主要貢獻，其余的模態(tài)可以忽略不計，則有：

由上式可以分別得到其幅值和相位值：

通過上式便可以識別出系統(tǒng)的模態(tài)參數(shù)。小波分析自身也存在一定的局限性，對具體信號依據(jù)何種原則來選擇小波基，目前在理論上和實際應(yīng)用中尚沒有統(tǒng)一的方法，而且小波變換還可能由于其有限長度的特點出現(xiàn)能量泄露等問題，另外，小波變換的實質(zhì)是一種線性變換，不能用于處理非線性問題。

朱宏平等［38］針對土木工程結(jié)構(gòu)前幾階自振頻率處于低頻區(qū)域以及環(huán)境激勵下結(jié)構(gòu)響應(yīng)信號信噪比很低的特點，著重論述采用小波方法抑制原始測量信號中的高頻噪聲部分，從而突出結(jié)構(gòu)低頻特性的基本原理，并用其對一實際高層結(jié)構(gòu)進行模態(tài)參數(shù)識別，識別結(jié)果顯示基于小波變換的模態(tài)分析方法較傳統(tǒng)模態(tài)分析方法識別精度更高更有效。邵艷秋［39］采用RDT法獲得環(huán)境激勵下結(jié)構(gòu)的自由響應(yīng)，通過改進的Morlet小波時頻分析，分解多自由度系統(tǒng)為若干個單自由度系統(tǒng)的組合，根據(jù)小波系數(shù)和系統(tǒng)固有頻率和阻尼比間關(guān)系識別出環(huán)境激勵下結(jié)構(gòu)模態(tài)參數(shù)，并應(yīng)用GARTEUR模型驗證小波變換方法的有效性；李夕兵等［40］分別應(yīng)用小波變換和HHT變換對實測爆破震動信號進行時頻分析，發(fā)現(xiàn)HHT較小波變換避免了選擇小波基的困難，在分析非平穩(wěn)信號時更具適應(yīng)性，而且小波譜由于受小波基有限長度及測不準原理的限制，容易造成頻譜擴散和能量泄露的問題。

2.3.2 Hilbert-Huang變換法

HHT作為一種針對非線性系統(tǒng)非平穩(wěn)信號的自適應(yīng)時頻處理方法，由Huang等人于1998年首次提出［41］。HHT主要由經(jīng)驗模態(tài)分解（EMD）和希爾伯特變換（HT）組成。經(jīng)驗模態(tài)分解是處理時域信號的方法，可將初始信號分解成一系列本征模態(tài)函數(shù)（IMF）和殘余項的疊加，實現(xiàn)對原信號的強制平穩(wěn)化處理，其中分解得到的IMF分量能夠表征信號在某一特征尺度上的振動模態(tài)；然后根據(jù)Hilbert變換，可以完成瞬時模態(tài)參數(shù)的求解（主要指瞬時頻率），其算法流程如圖9所示。

圖9 HHT算法流程圖

相對于傳統(tǒng)的時域模態(tài)參數(shù)識別方法，HHT將參數(shù)識別擴展到時頻域范圍，打破了傳統(tǒng)方法基于傅里葉變換只能在頻域或時域描述信號的缺陷，提出了瞬時頻率的概念，能夠用來解決非穩(wěn)態(tài)環(huán)境下的模態(tài)參數(shù)問題；相對于小波變換方法，Hilbert譜在時頻域的分辨率都要高于小波譜，具有更強的局部特性，而且能夠反映系統(tǒng)的非線性特性。但是對基于HHT的模態(tài)識別方法的研究也存在很多不足，其理論推導并不完善，對特定信號容易出現(xiàn)端點效應(yīng)，尤其對于模態(tài)頻率較為接近的結(jié)構(gòu)容易發(fā)生模態(tài)混淆的現(xiàn)象。

姚熊亮和張阿漫［42］基于EMD對水下爆炸結(jié)構(gòu)沖擊信號進行分析，通過研究每一個本征模態(tài)函數(shù)，從本質(zhì)上分析了水下爆炸沖擊加速度信號的組成成分及特點；陳雋和徐幼麟［43］結(jié)合青馬橋的實測動力響應(yīng)，研究了HHT方法在結(jié)構(gòu)模態(tài)參數(shù)識別中的應(yīng)用，并將其分別用于平穩(wěn)和非平穩(wěn)的實測記錄以識別結(jié)構(gòu)的模態(tài)信息，驗證了HHT相對傳統(tǒng)模態(tài)識別方法在處理非平穩(wěn)性數(shù)據(jù)方面具有明顯優(yōu)勢；湯寶平和何啟源等人［44］為克服噪聲信號對EMD分解的影響，提出基于小波去噪和HHT相結(jié)合的模態(tài)參數(shù)識別方法，該方法先利用小波進行信號去噪，以減少EMD分解過程的計算量和分解層數(shù)，然后利用HHT進行模態(tài)參數(shù)的識別，簡支梁系統(tǒng)實驗證明該算法對阻尼的識別精度明顯提高；韓建平和李達文［45］針對傳統(tǒng)模態(tài)識別方法對非線性和非穩(wěn)態(tài)信號處理能力差，且對阻尼比識別精度較低的現(xiàn)象，采用基于NExT和HHT相結(jié)合運行模態(tài)參數(shù)識別方法，通過對一個12層鋼筋混凝土框架模型進行模態(tài)分析，驗證了這種算法對阻尼比的識別效果較傳統(tǒng)基于傅里葉變換的半功率帶寬法有所改進，但對識別精確性仍難以確認；劉曉冰［46］采用基于NExT和HHT相結(jié)合運行模態(tài)參數(shù)識別方法，對船模結(jié)構(gòu)進行了模態(tài)識別研究，驗證了該方法有效避免NExTERA法中的虛假模態(tài)和系統(tǒng)定階問題，適用于船舶結(jié)構(gòu)參數(shù)識別及狀態(tài)監(jiān)測。

3 環(huán)境激勵下模態(tài)參數(shù)識別研究熱點

3.1 復雜密頻結(jié)構(gòu)的模態(tài)參數(shù)識別

傳統(tǒng)的模態(tài)參數(shù)識別方法，普遍要求被測結(jié)構(gòu)的各階模態(tài)頻率較分散，阻尼比較小，即模態(tài)混疊不嚴重，但是對于工程中的大型復雜結(jié)構(gòu)，特別是具有質(zhì)量或剛度突變的結(jié)構(gòu)，普遍含有密集模態(tài)，難以通過激振使結(jié)構(gòu)按某一階固有振型振動。產(chǎn)生密集模態(tài)的原因有兩個：一是各階模態(tài)的固有頻率過于靠近；二是模態(tài)阻尼比較大。如果兩個模態(tài)的固有頻率很接近，且各自的阻尼比也較大，則在頻譜上會出現(xiàn)如圖10［49］所示的模態(tài)混疊現(xiàn)象，無法用普通的傅里葉帶通濾波來實現(xiàn)對多模態(tài)的分離，進而影響多模態(tài)參數(shù)識別的精度。

圖10 模態(tài)混疊示意圖

傳統(tǒng)的基于傅里葉變換的模態(tài)識別方法不能直接用來識別復雜密頻結(jié)構(gòu)，因此時頻分析方法被更為廣泛地應(yīng)用于密集模態(tài)結(jié)構(gòu)的模態(tài)識別。陳雋等［47］采用HHT法對密頻結(jié)構(gòu)的阻尼進行識別，對兩個自由度的密頻和非密頻結(jié)構(gòu)分別采用HHT法和半功率帶寬法進行參數(shù)識別，結(jié)果表明HHT法對密頻結(jié)構(gòu)模態(tài)識別的識別精度明顯優(yōu)于半功率帶寬法，但是HHT法對密頻結(jié)構(gòu)的識別精度隨頻率密集度的提高而明顯降低，為解決密頻結(jié)構(gòu)的模態(tài)混疊問題，作者選擇采用在EMD分解過程中引用一種稱為間歇檢測準則的方法來消除混合現(xiàn)象；劉俊斐和李華軍［48］結(jié)合海洋平臺物理模型的動力響應(yīng)數(shù)據(jù)，采用加入偽信號技術(shù)的HHT方法對模態(tài)密集結(jié)構(gòu)進行模態(tài)分析，并對加入偽信號的幅值和頻率選擇進行改進，進而使用峰值法和HHT法結(jié)合偽信號技術(shù)對力錘激勵下的海洋平臺物理模型進行模態(tài)參數(shù)識別，得到結(jié)構(gòu)的自振頻率和阻尼比；黃應(yīng)來和董大偉等人［49］提出了一種針對密集模態(tài)結(jié)構(gòu)模態(tài)參數(shù)識別的方法，其基本思想是先用EMD分解化多模態(tài)參數(shù)識別問題為單模態(tài)參數(shù)識別問題，考慮到適當減小各階模態(tài)的阻尼比可以相應(yīng)降低各階模態(tài)間的混疊程度，采用加逆衰減指數(shù)窗與帶通濾波相結(jié)合的方法，將各個模態(tài)近似提取出來，再使用先進的時頻分析方法識別系統(tǒng)模態(tài)參數(shù)。圖11［49］和圖12［49］分別為密集模態(tài)結(jié)構(gòu)加速度響應(yīng)信號加逆衰減指數(shù)窗前后的時域圖和頻譜圖。

圖11 加逆衰減指數(shù)窗前后時域圖

圖12 加逆衰減指數(shù)窗前后時域圖

可以看出，通過加入適當?shù)哪嫠p指數(shù)窗，可以有效降低模態(tài)密集度，利用帶通濾波器就可以將多模態(tài)信號中的各階模態(tài)分離出來，然后利用Hilbert變換等方法識別系統(tǒng)的模態(tài)參數(shù)；曹軍宏和莊飚等人［50］提出采用基于時延優(yōu)化的二階盲辨識算法對模態(tài)密集結(jié)構(gòu)進行參數(shù)識別，其基本思想是利用量子遺傳算法對二階盲辨識算法中時延的選擇進行優(yōu)化，并對觀測信號的最優(yōu)時延二階協(xié)方差矩陣作對角化處理，同時得到模態(tài)振型和單自由度信號，然后基于FFT法提取出對應(yīng)各階模態(tài)的模態(tài)頻率和阻尼比；孫鵬和丁幼亮等人［51］從小波變換公式推導的角度對結(jié)構(gòu)模態(tài)混疊的機理進行研究分析，明確密集模態(tài)的定義并給出結(jié)構(gòu)模態(tài)密集度指標，即相鄰模態(tài)的頻率比，研究表明當相鄰模態(tài)的頻率比一定時，隨著阻尼比以及相鄰模態(tài)阻尼比差值的增加，模態(tài)混疊現(xiàn)象會愈加嚴重。

3.2 非穩(wěn)態(tài)環(huán)境激勵下的模態(tài)參數(shù)識別

頻域法和時域法中假設(shè)環(huán)境激勵需要近似表示為平穩(wěn)白噪聲激勵，這是因為頻域法和時域法都建立在傅里葉變換的基礎(chǔ)上，而傅里葉變換的全局特性限制其不能夠用來分析非平穩(wěn)信號。對于船舶等大型水上結(jié)構(gòu)，通常將環(huán)境激勵特征歸納為平穩(wěn)的強線譜與弱隨機譜的疊加，這也是對船舶環(huán)境激勵的一種簡化或者近似，真實的環(huán)境激勵應(yīng)當多為非穩(wěn)態(tài)的信號。常軍等［52］對隨機子空間法產(chǎn)生虛假模態(tài)和模態(tài)遺漏現(xiàn)象的原因進行深入分析，發(fā)現(xiàn)產(chǎn)生虛假模態(tài)的主要原因是由于輸入信號不滿足白噪聲假設(shè)，并由此提出分段平穩(wěn)隨機信號參數(shù)識別方法；李愛群等［53］在潤揚大橋的健康監(jiān)測中發(fā)現(xiàn)在風荷載及車輛荷載下大橋的動力響應(yīng)在一定的時間尺度內(nèi)是非平穩(wěn)的，不滿足傳統(tǒng)模態(tài)參數(shù)識別方法的白噪聲假設(shè)，因此以HHT和小波變換為代表的聯(lián)合時頻域法成為分析非穩(wěn)態(tài)環(huán)境激勵下模態(tài)參數(shù)識別的主要手段。

李中付和華宏星［54］認為任意隨機激勵信號都可以表示成白噪聲和非白噪聲信號的疊加，其一般形式寫作：

由此，我們可推導出線性系統(tǒng)響應(yīng)之間的相關(guān)函數(shù)由兩部分組成，一部分與脈沖響應(yīng)具有相同的數(shù)學形式，另一部分為其他形式；然后基于EMD分解，把非穩(wěn)態(tài)環(huán)境激勵下多自由度線性系統(tǒng)的模態(tài)參數(shù)辨識問題轉(zhuǎn)化為類似于已知各個單自由度系統(tǒng)的脈沖響應(yīng)進行參數(shù)辨識問題。這種算法主要有兩個優(yōu)點：一是無論白噪聲激勵、穩(wěn)態(tài)隨機激勵還是非穩(wěn)態(tài)隨機激勵，都可以僅根據(jù)結(jié)構(gòu)的響應(yīng)識別出線性結(jié)構(gòu)的模態(tài)參數(shù)；二是能夠有效地識別出環(huán)境激勵中的周期成分。

Chiang和Lin［55］假設(shè)外部激勵是由平穩(wěn)白噪聲信號和一個調(diào)幅信號乘積組成新的非平穩(wěn)隨機信號，他證明在這種假設(shè)條件下測得的響應(yīng)歷程的均方根函數(shù)中包含了激勵信號中的調(diào)幅成分，并可以通過曲線擬合將調(diào)幅信號提取出來，這樣就可以將非平穩(wěn)信號轉(zhuǎn)變?yōu)槠椒€(wěn)信號，然后計算不同測點間響應(yīng)信號的互相關(guān)函數(shù)，得到的信號便是一個可以進行模態(tài)識別的自由衰減信號。作者使用ITD法對其進行模態(tài)識別，并用一個六自由度的鏈式結(jié)構(gòu)進行數(shù)值實驗對該方法進行驗證。向律楷［56］在文章中將環(huán)境激勵由白噪聲擴展到MA（q）階模型，假設(shè)非平穩(wěn)激勵由d-1階多項式趨勢項與MA（q）階模型之和構(gòu)成，并使用NExT-ERA法識別出系統(tǒng)的模態(tài)參數(shù)。杜秀麗和汪鳳泉［57］將環(huán)境激勵視作均勻調(diào)制的萊維隨機激勵，并基于連續(xù)時間AR模型提出時域模態(tài)分解方法，數(shù)值實驗證明該算法能夠精確識別滿足這種激勵模型下的系統(tǒng)模態(tài)參數(shù)。

本文基于EMD分解，對非穩(wěn)態(tài)環(huán)境激勵下的三自由度系統(tǒng)進行模態(tài)參數(shù)識別，具體參數(shù)如式（13）-式（15）所示：

式中：m為模型質(zhì)量，kg；k為剛度，N/m；c為阻尼，N·s / m。

分別對2號質(zhì)點施加四種不同形式的激勵：平穩(wěn)白噪聲信號、平穩(wěn)白噪聲加簡諧信號、平穩(wěn)白噪聲加線性調(diào)頻信號（頻率從0到30 Hz隨時間線性增加）以及平穩(wěn)白噪聲和調(diào)幅信號（e-0.1t）的乘積，采樣頻率為500 Hz，采樣時間為30 s，用Newmark法求得三個質(zhì)點的加速度響應(yīng)，分別對任意兩點做互相關(guān)，然后對得到的互相關(guān)函數(shù)進行EMD分解，并用Hilbert變換得到該三自由度結(jié)構(gòu)在四種不同激勵下的模態(tài)參數(shù)，識別結(jié)果如下頁表1所示。由識別結(jié)果所示，說明該算法不僅能夠識別穩(wěn)態(tài)環(huán)境激勵下的系統(tǒng)模態(tài)參數(shù)，且可正確識別出各種不同形式非穩(wěn)態(tài)環(huán)境激勵下的系統(tǒng)模態(tài)參數(shù)。

3.3 時變系統(tǒng)的模態(tài)參數(shù)識別

目前針對船舶結(jié)構(gòu)的模態(tài)參數(shù)識別主要是將其視作線性時不變系統(tǒng)來進行識別，但視作線性時變系統(tǒng)更合理，這時傳統(tǒng)的運行模態(tài)分析方法便不再適用。然而目前對線性時變結(jié)構(gòu)的模態(tài)識別研究，也多是將線性時不變結(jié)構(gòu)參數(shù)識別理論基于“凍結(jié)時間”思想的直接推廣，即假設(shè)每一時刻的時間被瞬間“凍結(jié)”，這樣結(jié)構(gòu)就變?yōu)闀r不變結(jié)構(gòu)，然后按時不變結(jié)構(gòu)參數(shù)識別理論識別出每一時刻的模態(tài)參數(shù)，再用曲線擬合得到隨時間變化的模態(tài)參數(shù)［58］。李雷［59］提出一種基于“短時時不變假設(shè)”結(jié)合工況模態(tài)分析的識別方法，對飛行器舵片模型進行模態(tài)辨識，分析得到舵片模型在受熱情況下動力學參數(shù)的時變規(guī)律。但是這種做法忽略了結(jié)構(gòu)參數(shù)隨時間變化給結(jié)構(gòu)動力特性帶來的影響，造成某些情況下的分析結(jié)果不符合實際，因此這種模態(tài)參數(shù)一般被稱為“偽模態(tài)參數(shù)”［7］。

表1 四種不同工況下模態(tài)參數(shù)識別結(jié)果

為了準確描述時變系統(tǒng)模態(tài)參數(shù)隨時間變化的特性，以HHT和小波變換為代表的聯(lián)合時頻域法被更多地應(yīng)用到時變系統(tǒng)的模態(tài)參數(shù)識別過程中。劉建軍［60］提出基于HHT由線性時變SDOF結(jié)構(gòu)的自由振動響應(yīng)識別結(jié)果模態(tài)參數(shù)的新方法，并將其擴展到了MDOF結(jié)構(gòu)及強迫振動條件下結(jié)構(gòu)的模態(tài)識別，并用可調(diào)附加質(zhì)量的懸臂梁實驗驗證HHT用于識別時變系統(tǒng)模態(tài)參數(shù)的有效性；宋斌華［61］在劉建軍工作基礎(chǔ)上研究HHT和LMD［62］兩種新的時頻分析方法在時變結(jié)構(gòu)模態(tài)參數(shù)識別中的應(yīng)用；沈林［63］基于Daubechies小波對線性時變系統(tǒng)進行物理參數(shù)識別，并對線性時變系統(tǒng)在參數(shù)連續(xù)平滑變化、周期變化和突變的工況下進行仿真算例研究，驗證了Daubechies小波能夠正確有效地識別時變系統(tǒng)的物理參數(shù)。

我們基于HHT變換對一剛度和阻尼隨時間連續(xù)變化的兩自由度線性時變系統(tǒng)進行數(shù)值模型研究，采樣頻率為500 Hz，采樣時間為3 s，該數(shù)值模型的具體參數(shù)分別取為：

式中：m為模型質(zhì)量，kg；c為阻尼，N·s / m；k為剛度，N/m 。

將算例中的質(zhì)點m1在t = 0時刻拉離平衡位置2mm，另一個質(zhì)點保持不動，然后從靜止自由釋放，計算質(zhì)點m1的加速度響應(yīng)，對其作EMD分解，并對包含系統(tǒng)單階模態(tài)的瞬態(tài)響應(yīng)分量作Hilbert變換，識別出系統(tǒng)的模態(tài)瞬時頻率和阻尼比，識別結(jié)果如圖13、圖14所示。

圖13 模態(tài)瞬時頻率識別結(jié)果

圖14 模態(tài)阻尼比識別結(jié)果

在前面兩自由度時變系統(tǒng)數(shù)值實驗的基礎(chǔ)上，本文對一移動質(zhì)量-懸臂梁耦合時變系統(tǒng)進行實驗研究，實驗模型如圖15所示。

圖15 移動質(zhì)量-懸臂梁耦合系統(tǒng)實驗模型

將基于HHT變換的識別結(jié)果同傳統(tǒng)的特征系統(tǒng)實現(xiàn)法識別結(jié)果及有限元分析結(jié)果進行研究對比，驗證了HHT變換不但可以識別出線性時不變結(jié)構(gòu)的模態(tài)信息，而且能夠準確地識別得到時變結(jié)構(gòu)的模態(tài)參數(shù)，識別結(jié)果如圖16、圖17所示。

圖16 一階模態(tài)頻率識別結(jié)果

圖17 二階模態(tài)頻率識別結(jié)果

3.4 非線性系統(tǒng)的模態(tài)參數(shù)識別

非線性模態(tài)參數(shù)的概念最早由Rosenberg提出，他將模態(tài)分析思想擴展至非線性系統(tǒng)領(lǐng)域，Rosenberg給出描述非線性系統(tǒng)模態(tài)的頻響函數(shù)模型，并基于系統(tǒng)的響應(yīng)數(shù)據(jù)通過曲線擬合的方法識別出系統(tǒng)的非線性模態(tài)參數(shù)［64］，其典型特征有兩個：一是非線性系統(tǒng)激勵與響應(yīng)之間的關(guān)系無法事先確定，必須選擇正確的動力學模型對非線性系統(tǒng)進行描述；二是非線性系統(tǒng)的固有頻率及阻尼比等模態(tài)參數(shù)都是時間的函數(shù)，因此能夠用來刻畫系統(tǒng)模態(tài)參數(shù)隨時間變化規(guī)律的聯(lián)合時頻域法成為分析非線性系統(tǒng)模態(tài)參數(shù)的最主要手段。張也弛［65］綜合HHT和CxA方法，對強非線性系統(tǒng)進行了模態(tài)參數(shù)識別，同時指出這種算法能夠準確對非線性系統(tǒng)進行模態(tài)識別的前提是動力學模型能夠真實描述原物理過程；劉衛(wèi)華等［66］建立了含非線性參數(shù)單自由度振動系統(tǒng)的力學模型，如圖18所示，系統(tǒng)運動微分方程可寫成：

圖18 非線性單自由度系統(tǒng)力學模型

式中：α為弱非線性系數(shù)。

研究將已知非線性系統(tǒng)產(chǎn)生的混沌響應(yīng)作為系統(tǒng)激勵，并假定該響應(yīng)由若干個不穩(wěn)定周期軌道組成，從混沌響應(yīng)的狀態(tài)空間中提取出近似周期軌道，并采用諧波平衡法識別出非線性系統(tǒng)的模態(tài)參數(shù)；針對目前非線性系統(tǒng)模態(tài)參數(shù)識別的研究主要停留在穩(wěn)態(tài)環(huán)境激勵下進行，王鳳利和趙德有［67］提出非穩(wěn)態(tài)環(huán)境激勵下基于局域波法的非線性系統(tǒng)模態(tài)參數(shù)識別方法，算法思想是通過局域波分解得到等效線性系統(tǒng)，并通過Hilbert變換識別出轉(zhuǎn)子系統(tǒng)的模態(tài)參數(shù)。盡管已有很多學者開始致力于研究非線性系統(tǒng)的模態(tài)參數(shù)識別方法，但對非線性系統(tǒng)的模態(tài)分析還遠沒有達到線性模態(tài)理論研究那樣深入和透徹。

4 結(jié) 論

本文對環(huán)境激勵下的模態(tài)參數(shù)識別方法進行了歸納分類，重點分析這些經(jīng)典模態(tài)識別算法的理論思想和適用條件。主要對目前國內(nèi)研究熱點問題進行了綜述，并從數(shù)值算例和實驗的角度對這些問題進行了分析，為研究和應(yīng)用提供參考。通過上面對國內(nèi)相關(guān)文獻進行分類總結(jié)，本文關(guān)于環(huán)境激勵下模態(tài)參數(shù)識別的研究現(xiàn)狀有如下結(jié)論：

（1）對于環(huán)境激勵，目前普遍需要將激勵信號視作白噪聲或適當寬松地視為平穩(wěn)信號，即使在部分分析非穩(wěn)態(tài)環(huán)境激勵的文獻中也是將非穩(wěn)態(tài)環(huán)境激勵近似看成一種固定的表達形式，但是，在實際工程應(yīng)用中存在的激勵信號是任意形式的非平穩(wěn)信號，如何正確地給出環(huán)境激勵的近似表達形式尚待研究。

（2）如何從實驗模態(tài)中有效地甄別和提取噪聲模態(tài)，尤其對于運行狀態(tài)下的船舶結(jié)構(gòu)，環(huán)境激勵包括近似白噪聲激勵、有害噪聲以及船上螺旋槳主機引起的單頻簡諧激勵等，如何考慮這些簡諧激勵和噪聲模態(tài)給船舶結(jié)構(gòu)模態(tài)分析帶來的影響？也是目前可以研究的一個新方向。

（3）不同于實驗模態(tài)分析，環(huán)境激勵下的運行模態(tài)分析很可能會出現(xiàn)激勵不夠充分的現(xiàn)象，從導致部分階次的模態(tài)難以識別，因此在實驗過程中如何對結(jié)構(gòu)物正確施加激勵？也是工程實際中需要考慮的問題。

（4）目前對環(huán)境激勵下的結(jié)構(gòu)模態(tài)參數(shù)識別尚沒有一種統(tǒng)一完備的方法，各種方法都在一定程度存在缺陷，頻域法和時域法不能夠用來分析非平穩(wěn)信號，且存在由噪聲引起的虛假模態(tài)問題，小波分析不能解決非線性問題且存在分辨率無法保證和能量泄露等問題，而HHT法中針對瞬時頻率的數(shù)學定義在信號處理領(lǐng)域尚未達成共識，IMF分量之間的正交性也沒有完整的數(shù)學證明，且存在虛假模態(tài)和端點效應(yīng)等問題。

（5）模態(tài)識別過程中對系統(tǒng)階次確定是十分關(guān)鍵和棘手的問題，盡管已經(jīng)有學者提出采用模態(tài)置信度及穩(wěn)定圖法來確定系統(tǒng)階次，但噪聲模態(tài)的影響仍不能做到完全消除，識別精度還有待改善。

（6）目前對系統(tǒng)固有頻率已能準確識別，但是對結(jié)構(gòu)阻尼的識別效果往往不好，仍是亟待解決的問題。

［1］ 曹樹謙，張文德，蕭龍翔.振動結(jié)構(gòu)模態(tài)分析［M］.天津：天津大學出版社，2001.

［2］ 萬嶺.環(huán)境激勵下的船舶結(jié)構(gòu)模態(tài)參數(shù)識別［D］.大連：大連理工大學，2010.

［3］ 姜大正.環(huán)境激勵下船舶結(jié)構(gòu)模態(tài)分析實驗與理論研究［D］.大連：大連理工大學，2009.

［4］ 許俊臣.船舶加筋結(jié)構(gòu)的運行模態(tài)分析及動力減振優(yōu)化［D］.大連：大連理工大學，2013.

［5］ 續(xù)秀忠，華宏星，陳兆能.基于環(huán)境激勵的模態(tài)參數(shù)辨識方法綜述［J］.振動與沖擊，2002（3）：1-5.

［6］ 劉宇飛，辛克貴，樊健生，等.環(huán)境激勵下結(jié)構(gòu)模態(tài)參數(shù)識別方法綜述［J］.工程力學，2014（4）：1-5.

［7］ 王學敏.基于Hilbert-Huang變換的橋梁監(jiān)測信號分析與處理和時變模態(tài)參數(shù)識別［D］.長沙：中南大學，2008.

［8］ Roeck G De，Peeters B，Ren W X. Benchmark Study on System Identification through Ambient Vibration Measurements［C］. Proceedings of the 18th International Modal Analysis Conference（IMAC），2000：1106-1112.

［9］ 張毅剛，劉才瑋，吳金志，等.適用空間網(wǎng)格結(jié)構(gòu)模態(tài)識別的改進功率譜峰值法［J］.振動與沖擊，2013（9）：10-15.

［10］ 夏祥麟.環(huán)境激勵模態(tài)分析方法的比較［D］.長沙：中南大學，2013.

［11］ Brincker R，Zhang Ling-mi，Andersen P. Modal Identification from Ambient Responses Using Frequency Domain Decomposition［C］. Proceedings of the 18th International Modal Analysis Conference（IMAC），2000：625-630.

［12］ Akaike H. Power Spectrum Estimation Through Autoregressive Model Fitting［J］. Annals of the Institute of Statistical Mathematics，1969，21：407.

［13］ 郭永剛，許亮華，水小平.基于脈沖響應(yīng)數(shù)據(jù)的ARMA法建模以及模態(tài)參數(shù)識別［J］.地震工程與工程振動，2006（5）：167-171.

［14］ Smail M，Thomas M，Lakis A A.Assessment of Optimal ARMA Model Orders for Modal Analysis［J］.Mechanical Systems and Signal Processing，1999（5）：803-819.

［15］ Van Overschee P，De Moor B. Subspace Algorithms for the Stochastic Identification Problem［C］. Proceedings of the 30th IEEE Conference on Decision and Control，1991，1321-1326.

［16］ 辛峻峰，王樹青，劉福順.數(shù)據(jù)驅(qū)動與協(xié)方差驅(qū)動隨機子空間法差異化分析［J］.振動與沖擊，2013（9）：1-5.

［17］ 辛峻峰，盛進路，張永波.數(shù)據(jù)驅(qū)動隨機子空間法矩陣維數(shù)選擇與噪聲問題研究［J］.振動與沖擊，2013（16）：152-157.

［18］ 常軍，張啟偉，孫利民.隨機子空間產(chǎn)生虛假模態(tài)及模態(tài)遺漏的原因分析［J］.工程力學，2007（11）：57-62.

［19］ 章國穩(wěn)，湯寶平，孟利波.基于特征值分解的隨機子空間算法研究［J］.振動與沖擊，2012（7）：74-78.

［20］ Cole H A. On-Line Failure Detection and Damping Measurement of Aerospace Structures by the Random Decrement Signatures［J］.1973，NASACR-2205.

［21］ Ibrahim S R. Efficient Random Decrement Computation for Identification of Ambient Responses［C］.Proceedings of the 19th IMAC，F(xiàn)lorida，USA，2001：1-6.

［22］ James G H，Came T G，Edmunds R S. STARS Missile-Modal Analysis of First-Flight Data Using the Natural Excitation Technique［C］. Proceedings of the 12th International Modal Analysis Conference（IMAC），1994.

［23］ 羅奎.基于自然激勵技術(shù)的結(jié)構(gòu)模態(tài)參數(shù)識別應(yīng)用研究［D］.武漢：武漢理工大學，2010.

［24］ 紀曉東，錢稼茹，徐龍河.模擬環(huán)境激勵下結(jié)構(gòu)模態(tài)參數(shù)識別試驗研究［J］.清華大學學報，2006（6）：769-772.

［25］ Ibrahim S R，Mikulcik E. A Method for the Direct Identification of Vibration Parameters from the Free Response［J］.The Shock and Vibraton Bulletin，1977（4）：183-198.

［26］ 楊佑發(fā)，李帥，李海龍.環(huán)境激勵下結(jié)構(gòu)模態(tài)參數(shù)識別的改進ITD法［J］.振動與沖擊，2014（1）：194-199.

［27］ 周凱.虛擬測點對ITD法識別結(jié)果的影響［J］.中國水運，2009（2）：184-186.

［28］ Brown D L，Allemang R J，Zimmerman R，et al. Parameter Estimation Techniques for Modal Analysis［J］. SAE Paper，No.790221，SAE Trans.，1979（1）：828-846.

［29］ Vold H，Rocklin G T. The Numerical Implementation of a Multi-input Modal Estimation Method for Minicomputers［C］. Orlando，F(xiàn)L，USA：Proceedings of the 1st International Modal Analysis Conference（IMAC），1982：542-548.

［30］ 常軍，張啟偉，孫利民.穩(wěn)定圖方法在隨機子空間識別模態(tài)參數(shù)中的應(yīng)用［J］.工程力學，2007（2）：39-44.

［31］ Ibrahim S R. Computation of Normal Modes from Identified Complex Modes［J］.AIAAJOURNAL，1982（3）：446-451.

［32］ Juang J N，Pappa R S.An Eigensystem Realization Algorithm（ERA）for Modal Parameter Identification and Model Reduction［J］.NASA/JPL，1984（5）：620-627.

［33］ Xu Junchen，Hong Ming，Liu Xiaobing. Operational Modal Analysis of a Ship Model in the Presence of Harmonic Excitation［J］. J. Marine Sci. Appl.，2013（1）：38-44.

［34］ 萬嶺，洪明，許俊臣.基于NExT/ERA方法緩解激勵下的船體結(jié)構(gòu)模型總體模態(tài)參數(shù)識別［J］.船舶力學，2013（7）：774-784.

［35］ 楊和振，李華軍，黃維平.海洋平臺結(jié)構(gòu)環(huán)境激勵的實驗模態(tài)分析［J］.振動與沖擊，2005（2）：129-132.

［36］ Grossmann A，Morlet J. Decomposition of Hardy Function into Square Integrable Wavelet of Constant Shape［J］. SIAM J. of Math，Anal，1984（4）：723-736.

［37］ 譚冬梅，姚三，瞿偉廉.振動模態(tài)的參數(shù)識別綜述［J］.華中科技大學學報，2002（3）：73-78.

［38］ 朱宏平，翁順.運用小波分析方法進行結(jié)構(gòu)模態(tài)參數(shù)識別［J］.振動與沖擊，2007（4）：1-5.

［39］ 邵艷秋.環(huán)境激勵下系統(tǒng)模態(tài)參數(shù)識別方法研究［D］.南京：南京航空航天大學，2008.

［40］ 李夕兵，張義平，劉志祥，等.爆破震動信號的小波分析與HHT變換［J］.爆炸與沖擊，2005（6）：528-534.

［41］ Huang N E，Zheng S，Steven R L. The empirical mode decomposition and the Hilbert spectrum for nonlinear random-stationary time series analysis［C］.Proc R Soc London Ser A，1998，454：903-95.

［42］ 姚熊亮，張阿漫.經(jīng)驗模態(tài)分解方法在結(jié)構(gòu)沖擊信號分析中的應(yīng)用［J］.中國艦船研究，2006（4）：11-15.

［43］ 陳雋，徐幼麟. HHT方法在結(jié)構(gòu)模態(tài)參數(shù)識別中的應(yīng)用［J］.振動工程學報，2003（3）：383-388.

［44］ 湯寶平，何啟源，蔣恒恒，等.利用小波去噪和HHT的模態(tài)參數(shù)識別［J］.振動、測試與診斷，2009（2）：197-200.

［45］ 韓建平，李達文.基于Hilbert-Huang變換和自然激勵技術(shù)的模態(tài)參數(shù)識別［J］.工程力學，2010（8）：54-59.

［46］ 劉曉冰.基于HHT方法船舶結(jié)構(gòu)運行模態(tài)參數(shù)識別［D］.大連：大連理工大學，2014.

［47］ 陳雋，徐幼麟，李杰. Hilbert-Huang變換在密頻結(jié)構(gòu)阻尼識別中的應(yīng)用［J］.地震工程與工程振動，2003（4）：34-42.

［48］ 劉俊斐，李華軍.HHT方法在模態(tài)密集結(jié)構(gòu)模態(tài)參數(shù)識別中的應(yīng)用研究［J］.中國海洋大學學報（自然科學版），2007（S2）：181-188.

［49］ 黃應(yīng)來，董大偉，閆兵.密集模態(tài)分離及其參數(shù)識別方法研究［J］.機械強度，2009（1）：8-13.

［50］ 曹軍宏，莊飚，韋灼彬.時延優(yōu)化二階盲辨識應(yīng)用于密集模態(tài)結(jié)構(gòu)參數(shù)識別［C］.第22屆全國結(jié)構(gòu)工程學術(shù)會議論文集（第1冊），2013.

［51］ 孫鵬，丁幼亮，張勁泉，等.基于Morlet小波變換的結(jié)構(gòu)密集模態(tài)參數(shù)識別［J］.東南大學學報，2012（2）：339-345.

［52］ 常軍，孫利民，張啟偉.一種分段平穩(wěn)隨機信號的參數(shù)識別方法［J］.振動與沖擊，2007（6）：17-20.

［53］ 李愛群，丁幼亮，費慶國，等.潤揚大橋斜拉橋模態(tài)頻率識別的環(huán)境變異性［J］.東南大學學報，2007（2）：245-250.

［54］ 李中付，華宏星，宋漢文，等.非穩(wěn)態(tài)環(huán)境激勵下線性結(jié)構(gòu)的模態(tài)參數(shù)辨識［J］.振動工程學報，2002（2）：139-143.

［55］ Chiang D Y，Lin C S. Identification of Modal Parameters from Nonstationary Ambient Vibration Data Using Correlation Technique［J］. AIAA Journal，2008（11）：2752-2759.

［56］ 向律楷.非平穩(wěn)環(huán)境激勵下大橋結(jié)構(gòu)的模態(tài)參數(shù)識別［D］.北京：清華大學，2009.

［57］ 杜秀麗，汪鳳泉.非平穩(wěn)環(huán)境激勵下模態(tài)參數(shù)識別的連續(xù)時間AR方法［J］.中國科學，2009（10）：1736-1742.

［58］ 于開平，鄒經(jīng)湘，龐世偉.結(jié)構(gòu)系統(tǒng)模態(tài)參數(shù)識別方法研究進展［J］.世界科技研究與發(fā)展，2005（6）：22-30.

［59］ 李雷.工況模態(tài)分析在時變航天器結(jié)構(gòu)中的應(yīng)用［D］.上海：復旦大學，2011.

［60］ 劉建軍. Hilbert-Huang變換及其在線性時變結(jié)構(gòu)模態(tài)參數(shù)識別中的應(yīng)用［D］.長沙：中南大學，2007.

［61］ 宋斌華.基于Hilbert-Huang變換和局部均值分解的時變結(jié)構(gòu)模態(tài)參數(shù)識別［D］.長沙：中南大學，2009.

［62］ 程軍圣，張亢，楊宇，等.局部均值分解與經(jīng)驗模式分解的對比研究［J］.振動與沖擊，2009（5）：13-16.

［63］ 沈林.基于小波方法的線性時變系統(tǒng)參數(shù)識別［D］.南京：南京航空航天大學，2006.

［64］ Rosenberg R M. Normal Mode in Nonlinear Dual Mode Systems［J］. J Appl Mech，1960，27：263-268.

［65］ 張也弛.基于Hilbert-Huang變換的非線性系統(tǒng)參數(shù)辨識［D］.哈爾濱：哈爾濱工業(yè)大學，2008.

［66］ 劉衛(wèi)華，丁旺才，田海勇.混沌激勵下振動系統(tǒng)的非線性參數(shù)識別［J］.振動與沖擊，2009（5）：80-83.

［67］ 王鳳利，趙德有.基于局域波法的非線性系統(tǒng)模態(tài)參數(shù)識別研究［J］.大連理工大學學報，2011（1）：56-60.

Review of modal identifi cation of ship structure based on ambient excitation

HONG Ming LEI Chuan CUI Hong-yu
(School of Naval Architecture, Dalian University of Technology, Dalian 116024, China)

Based on the theory of modal analysis, this paper summarizes the modal identifi cation methods under ambient excitation, and discusses fundamental concepts and the advantages and disadvantages of these methods. It especially introduces the modal parameter identifi cation of ship structure under the ship sailing condition and hot issues about its domestic researches, which provide benefi cial references for the relevant researchers.

ambient excitation; modal parameters identifi cation; ship structure; operational modal analysis

U661.4

A

1001-9855（2014）05-0001-16

國家自然基金資助“基于運行模態(tài)識別的艦船加筋結(jié)構(gòu)動力優(yōu)化方法及實驗研究（51079027）”。

2014-08-18

洪 明（1959-），男，教授，研究方向：結(jié)構(gòu)失效、振動和噪聲機理、預報及控制研究。

雷 川（1990-），男，碩士，研究方向：結(jié)構(gòu)分析及模態(tài)參數(shù)識別。

崔洪宇（1977-），男，講師，研究方向：結(jié)構(gòu)振動控制及故障診斷。