范倪圣，陳穩(wěn)杰，周振綱

(嘉興學(xué)院建筑工程學(xué)院，浙江 嘉興 314001)

隨機子空間方法[1-2]只需要根據(jù)測點響應(yīng)數(shù)據(jù)即可識別結(jié)構(gòu)工作模態(tài)參數(shù)，但該方法進(jìn)行模態(tài)識別時可能存在模態(tài)遺漏和出現(xiàn)虛假模態(tài)問題[3-5]，對于此，國內(nèi)外學(xué)者通過采用穩(wěn)定圖方法或模糊聚類方法等對其進(jìn)行了改進(jìn)，在一定程度上避免了出現(xiàn)虛假模態(tài)問題，較好地解決了隨機子空間方法系統(tǒng)階數(shù)n取值問題。然而，該方法的識別精度實際還會受到方法中相關(guān)參數(shù)(如Hankel矩陣行數(shù)i、測點數(shù)L，識別用到的測點響應(yīng)數(shù)據(jù)量s)取值的影響。國內(nèi)個別學(xué)者就Hankel矩陣行數(shù)i對隨機子空間方法模態(tài)參數(shù)識別精度問題做了初步研究[6-7]。但同時考慮參數(shù)i和s不同取值組合對隨機子空間方法模態(tài)參數(shù)識別精度的影響規(guī)律研究還尚未見相關(guān)報道。因此，為了將隨機子空間方法更好地應(yīng)用于環(huán)境激勵下實際工程結(jié)構(gòu)的模態(tài)參數(shù)識別中，有必要對參數(shù)s和i不同取值組合對隨機子空間方法模態(tài)參數(shù)識別精度的影響規(guī)律問題進(jìn)行研究。

1 隨機子空間方法模態(tài)參數(shù)識別方法

采用基于協(xié)方差的隨機子空間方法進(jìn)行結(jié)構(gòu)模態(tài)參數(shù)識別時，首先根據(jù)各測點的加速度響應(yīng)測試數(shù)據(jù)構(gòu)造如式(1)所示的Hankel矩陣。

(1)

其中,H為Hankel矩陣；Hp為Hankel矩陣過去分塊矩陣；Hf為Hankel矩陣未來分塊矩陣；yk(k=0,…,2i+j-2)為根據(jù)測點響應(yīng)數(shù)據(jù)構(gòu)造的列向量，上述構(gòu)造的列向量的行數(shù)為測點數(shù)L；i和j分別為Hankel矩陣分塊矩陣的行數(shù)和列數(shù)，Hankel矩陣的分塊矩陣的行數(shù)和列數(shù)應(yīng)滿足關(guān)系式s=2i+j-1，s為計算時選取的測點響應(yīng)數(shù)據(jù)數(shù)量。

根據(jù)生成的Hankel矩陣過去分塊矩陣和未來分塊矩陣按式(2)計算得到Toeplitz矩陣T：

(2)

對T矩陣按式(3)進(jìn)行奇異值分解：

(3)

其中,U,S,V均為T矩陣經(jīng)奇異值分解得到的矩陣，U和V為正交矩陣，S為對角矩陣；S1為S矩陣的分塊矩陣，其為對角矩陣，對角元素T為矩陣異值分解得到的正奇異值，S1∈Rn×n；U1和U2為U矩陣的分塊矩陣，U1∈RLi×n，U2∈RLi×(Li-n)；V1和V2為V矩陣的分塊矩陣，V1∈RLi×n，V2∈RLi×(Li-n)。

系統(tǒng)的觀測矩陣Γ及其逆矩陣Γ-1可表示為：

(4)

(5)

其中，?！蔙Li×n；Γ-1∈Rn×Li。

根據(jù)式(6)，式(7)計算得到矩陣P1和P2。

P1=Γ-1(1∶n,1∶L(i-1))

(6)

P2=Γ(L+1∶Li,1∶n)

(7)

其中，P1∈Rn×(Li-L+1)；P2∈R(Li-L)×n；Γ-1(1∶n,1∶L(i-1))為取Γ-1矩陣第1行至第n行和第1列至第L(i-1)列范圍內(nèi)元素的構(gòu)成的矩陣；Γ(L+1∶Li,1∶n)為取Γ矩陣第L+1行至第Li行和第1列至第n列范圍內(nèi)元素的構(gòu)成的矩陣。

系統(tǒng)的狀態(tài)矩陣A和輸出矩陣Q可表示為：

A=P1P2

(8)

Q=Γ(1∶L,1∶n)

(9)

其中，A∈Rn×n;Q∈RL×n;Γ(1∶L,1∶n)為取Γ矩陣第1行至第L行和第1列至第n列范圍內(nèi)元素的構(gòu)成的矩陣。

對計算得到的狀態(tài)矩陣A按式(10)進(jìn)行特征值分解：

A=EDE-1

(10)

其中，E為特征向量構(gòu)成的矩陣，E∈Rn×n；D為特征值構(gòu)成的對角矩陣，D∈Rn×n。矩陣E和矩陣D的元素都為復(fù)數(shù)。

提取矩陣D的對角線元素構(gòu)成列向量Z，即：

Z={z1z2…zn}T

(11)

其中，zk為Z向量的第k行元素值，其大小等于D矩陣第k行的對角元素值；Z∈Rn×1。

根據(jù)式(12)計算得到識別模態(tài)頻率列向量Ft。

Ft=abs(ln(Z)·fs)/2π

(12)

其中，fs為測點響應(yīng)測量的采樣頻率，Hz；π為圓周率；ln(·)函數(shù)為對復(fù)向量元素進(jìn)行自然對數(shù)運算；abs(·)函數(shù)為對復(fù)向量元素進(jìn)行取模運算；Ft∈Rn×1。

根據(jù)輸出矩陣Q和狀態(tài)矩陣A特征值分解得到的特征向量構(gòu)成的矩陣E按式(13)計算復(fù)矩陣G。

G=QE

(13)

按式(14)計算得到識別模態(tài)振型向量構(gòu)成的振型矩陣Mt。

Mt=abs(G)*sign(180·(|angle(G)|-π/2)/π)

(14)

其中，angle(·)函數(shù)為對復(fù)矩陣元素進(jìn)行取相位運算；abs(·)函數(shù)為對復(fù)矩陣元素進(jìn)行取模運算；sign(·)函數(shù)為符號函數(shù)，當(dāng)因變量分別為正值、零和負(fù)值時，函數(shù)值分別取1，0和-1。

2 隨機子空間方法模態(tài)參數(shù)識別精度評價方法

本文主要就參數(shù)i和s取值對基于協(xié)方差的隨機子空間方法的模態(tài)頻率和模態(tài)振型的識別精度影響規(guī)律進(jìn)行研究。以下對本文模態(tài)頻率和模態(tài)振型識別精度計算方法進(jìn)行簡介。

2.1 模態(tài)頻率識別精度評價方法

前述隨機子空間方法識別得到的結(jié)構(gòu)的各階模態(tài)頻率構(gòu)成的列向量Ft可表示為：

(15)

其中，ftk為根據(jù)試驗測量數(shù)據(jù)由隨機子空間方法識別得到的第k階模態(tài)的模態(tài)頻率，k=1,2,…,m，m為識別模態(tài)數(shù)。

與隨機子空間方法識別得到的各階模態(tài)對應(yīng)的結(jié)構(gòu)真實的模態(tài)頻率構(gòu)成的列向量Fa可表示為：

(16)

其中，F(xiàn)a∈Rn×1；fak為與隨機子空間方法識別得到的第k階模態(tài)對應(yīng)的結(jié)構(gòu)真實的模態(tài)頻率，k=1,2,…,m，m為識別模態(tài)數(shù)。

根據(jù)式(17)可計算得到隨機子空間方法識別得到的各階模態(tài)的模態(tài)頻率識別精度。

(17)

其中，eFk為與隨機子空間方法識別得到的第k階模態(tài)的模態(tài)頻率識別精度。

參數(shù)i和s取值對隨機子空間方法不同階模態(tài)的模態(tài)頻率識別精度的影響有一定差異，故本文采用按式(18)計算得到的模態(tài)頻率平均識別精度指標(biāo)EM來對隨機子空間方法模態(tài)頻率識別精度進(jìn)行評價。EM值越大則隨機子空間法的模態(tài)頻率整體識別精度越高。

(18)

其中，m為識別模態(tài)數(shù)。

2.2 模態(tài)振型識別精度評價方法

前述隨機子空間方法識別得到的結(jié)構(gòu)的各階模態(tài)振型構(gòu)成的振型矩陣Mt可表示為：

(19)

其中，φtk為根據(jù)試驗測量數(shù)據(jù)由隨機子空間方法識別得到的第k階模態(tài)的模態(tài)振型列向量，k=1,2,…,m。

與隨機子空間方法識別得到的各階模態(tài)對應(yīng)的結(jié)構(gòu)真實的模態(tài)振型向量構(gòu)成的振型矩陣Ma可表示為：

(20)

其中，φak為與隨機子空間方法識別得到的第k階模態(tài)對應(yīng)的結(jié)構(gòu)真實的模態(tài)振型列向量，k=1,2,…,m。

模態(tài)振型識別的精度可采用識別模態(tài)振型與真實模態(tài)振型的模態(tài)保證準(zhǔn)則MAC值來評價。反映隨機子空間方法各階模態(tài)的模態(tài)振型識別精度的MAC值可按式(21)計算得到。

(21)

其中，MACk為隨機子空間方法識別得到的第k階模態(tài)模態(tài)振型的MAC值。

參數(shù)i和s取值對隨機子空間方法不同階模態(tài)的模態(tài)振型識別精度的影響有一定差異，故本文采用按式(22)計算得到的模態(tài)振型平均模態(tài)模態(tài)保證準(zhǔn)則MACM值對隨機子空間方法模態(tài)振型的識別精度進(jìn)行評價。MACM值越大則隨機子空間法的模態(tài)振型整體識別精度越高。

(22)

3 五質(zhì)點的彈簧-質(zhì)量-阻尼系統(tǒng)模態(tài)識別算例

五質(zhì)點的彈簧-質(zhì)量-阻尼系統(tǒng)算例模型如圖1所示。算例模型的質(zhì)量、剛度和阻尼參數(shù)如表1所示。

表1 算例模型參數(shù)

在該算例模型的節(jié)點1和節(jié)點2處施加豎向白噪聲激振力f1(t)和f2(t)來模擬算例模型系統(tǒng)受到的環(huán)境激勵，通過計算獲得各節(jié)點豎向加速度響應(yīng)數(shù)據(jù)作為本文隨機子空間方法模態(tài)識別的測點響應(yīng)數(shù)據(jù)。

根據(jù)各節(jié)點的動力平衡關(guān)系可建立算例模型受上述白噪聲激勵時的動力方程[8]，如式(23)所示。

(23)

(24)

(25)

(26)

(27)

(28)

(29)

(30)

在節(jié)點1和節(jié)點2處均輸入相同的高斯白噪聲激勵力，高斯白噪聲激振力時程及功率譜曲線如圖2所示，激振力時長為10 s，步長為0.001 s。采用Newmark-β法用Matlab軟件編程求解算例模型動力響應(yīng)時程。限于篇幅，此處僅給出節(jié)點1加速度響應(yīng)時程及功率譜曲線，如圖3所示。

算例模型的模態(tài)頻率和歸一化模態(tài)振型的真實值(理論值)可通過對式(23)忽略阻尼項和無激振力時得到的無阻尼自由振動方程進(jìn)行特征值和特征向量分析計算得到。算例模型的模態(tài)頻率和模態(tài)振型真實值分別如表2，圖4～圖8所示。

表2 i和s不同取值時隨機子空間法模態(tài)頻率和模態(tài)振型識別平均精度初步比較

在測點數(shù)L=5和系統(tǒng)階次n=12取值一定時，對i=38，s=6 542和i=50，s=2 838兩種參數(shù)取值情況根據(jù)前述基于協(xié)方差的隨機子空間模態(tài)參數(shù)識別方法對算例模型進(jìn)行了模態(tài)參數(shù)識別計算，識別得到的各階模態(tài)頻率和模態(tài)振型如表2，圖4～圖8所示。同時，給出了兩種參數(shù)取值情況時模態(tài)頻率平均識別精度指標(biāo)EM值和模態(tài)振型平均識別精度指標(biāo)MACM值，如表2所示。

由表2，圖4～圖8可知：i=38，s=6 542時隨機子空間方法識別得到的算例模型各階模態(tài)頻率和模態(tài)振型識別值與其真實值非常接近，EM和MACM都接近于1；而i=50，s=2 838時模態(tài)頻率和模態(tài)振型的平均識別精度要低于i=38，s=6 542取值情況。由此可見，參數(shù)i和s取值對隨機子空間方法模態(tài)參數(shù)識別精度有一定影響，故有必要對參數(shù)i和s取值對隨機子空間方法模態(tài)參數(shù)識別精度的影響規(guī)律進(jìn)行研究。

4 i和s取值對模態(tài)參數(shù)識別精度影響規(guī)律研究

為了進(jìn)一步研究參數(shù)i和s取值對隨機子空間方法模態(tài)參數(shù)識別精度的影響規(guī)律，對i和s不同取值組合時采用隨機子空間方法對前述算例模型的模態(tài)頻率和模態(tài)振型進(jìn)行了識別計算，并計算得到了各取值組合時模態(tài)頻率平均識別精度指標(biāo)EM值和模態(tài)振型平均識別精度指標(biāo)MACM值。參數(shù)i和s不同取值組合構(gòu)造時參數(shù)i取值從30到158，增量為2，參數(shù)s取值從2 000到8 998，增量為2，共構(gòu)造了227 500個i和s不同取值組合。識別計算中測點數(shù)L=5和系統(tǒng)階次n=12保持不變。

圖9給出了i和s不同取值組合時采用隨機子空間法識別得到的算例模型的模態(tài)頻率平均識別精度指標(biāo)EM值變化情況。

圖10給出了i和s不同取值組合時采用隨機子空間法識別得到的算例模型的模態(tài)振型平均識別精度指標(biāo)MACM值變化情況。

由圖9可知：參數(shù)i和s不同取值組合時，隨機子空間法模態(tài)頻率的平均識別精度指標(biāo)EM值均較大，EM值在0.971～0.998范圍內(nèi)變化，表明隨機子空間法識別算例模型的模態(tài)頻率具有很高的精度。參數(shù)i和s不同取值對隨機子空間法模態(tài)頻率的平均識別精度指標(biāo)EM有一定影響，主要表現(xiàn)為：

1)當(dāng)參數(shù)i在42～110取值范圍內(nèi)時，參數(shù)s取值對模態(tài)頻率平均識別精度指標(biāo)EM的影響相對較明顯。

2)參數(shù)i取值對模態(tài)頻率平均識別精度指標(biāo)值EM的影響相對明顯，參數(shù)i取值在36～42及110～158取值范圍內(nèi)時，模態(tài)頻率平均識別精度更高。

由圖10可知：參數(shù)i和s不同取值組合時，隨機子空間法模態(tài)振型的平均識別精度指標(biāo)MACM值變化較大，MACM值在0.283～0.999范圍內(nèi)變化，表明參數(shù)i和s取值對隨機子空間法識別算例模型的模態(tài)振型的精度影響較大。參數(shù)i和s不同取值對隨機子空間法模態(tài)振型的平均識別精度指標(biāo)MACM的影響主要表現(xiàn)為：

1)當(dāng)參數(shù)i在30～96和130～136取值范圍內(nèi)時，參數(shù)s取值對模態(tài)振型平均識別精度指標(biāo)MACM的影響較明顯。

2)參數(shù)i取值對模態(tài)振型平均識別精度指標(biāo)MACM的影響明顯，參數(shù)i取值在22～24,32～36,60～62,72～78,82～84,102～120,136～138,142～148及156～158取值范圍內(nèi)時，模態(tài)振型平均識別精度較高，MACM值基本在0.7及以上。

綜上可知：

1)參數(shù)i和s取值對隨機子空間法模態(tài)頻率的平均識別精度指標(biāo)EM有一定影響，但影響總體相對較小。

2)參數(shù)i和s取值對隨機子空間法模態(tài)振型的平均識別精度指標(biāo)MACM影響較大。

3)參數(shù)i和s在若干取值范圍內(nèi)取值組合時隨機子空間法的模型頻率及振型均具有較高的識別精度，如i取36～38,s取2 500～7 000時，模態(tài)頻率的平均識別精度可達(dá)0.99以上，模態(tài)振型的平均識別精度可達(dá)0.8以上。

5 結(jié)論及展望

本文以一五質(zhì)點的彈簧-質(zhì)量-阻尼系統(tǒng)的模態(tài)參數(shù)識別算例為例，對參數(shù)i和s取值對基于協(xié)方差的隨機子空間方法的模態(tài)頻率及模態(tài)振型的識別精度影響進(jìn)行了研究。主要結(jié)論如下：

1)參數(shù)i和s不同取值組合時隨機子空間法模態(tài)頻率的平均識別精度指標(biāo)EM值在0.971～0.998范圍內(nèi)變化，變化幅度較小，表明參數(shù)i和s取值對隨機子空間法模態(tài)頻率的平均識別精度影響相對較小，但仍有一定影響，同時也表明隨機子空間法識別模態(tài)頻率具有很高的精度。

2)參數(shù)i和s不同取值組合時隨機子空間法模態(tài)振型的平均識別精度指標(biāo)MACM在0.283～0.999范圍內(nèi)變化，變化幅度較大，表明參數(shù)i和s取值對隨機子空間法模態(tài)振型的平均識別精度影響較大。

3)參數(shù)i和s在若干取值范圍內(nèi)取值組合時隨機子空間法的模型頻率及振型均具有較高的識別精度。

本文對參數(shù)i和s不同取值組合時隨機子空間法模態(tài)頻率和模型振型的平均識別精度的影響進(jìn)行了初步研究。為了更好地將隨機子空間方法應(yīng)用于環(huán)境激勵下實際工程結(jié)構(gòu)的模態(tài)參數(shù)識別中，后續(xù)將結(jié)合相關(guān)數(shù)學(xué)方法對參數(shù)i和s合理取值方法進(jìn)行進(jìn)一步的研究。