模態(tài)跳躍現(xiàn)象對振動分析的影響研究

張　淼,　于　瀾,　鞠　偉

(1.長春工程學(xué)院 理學(xué)院,吉林 長春　130012; 2.中國第一汽車股份有限公司 技術(shù)中心,吉林 長春　130011)

?

模態(tài)跳躍現(xiàn)象對振動分析的影響研究

張淼1,于瀾1,鞠偉2

(1.長春工程學(xué)院 理學(xué)院,吉林 長春130012; 2.中國第一汽車股份有限公司 技術(shù)中心,吉林 長春130011)

摘要:文章揭示了模態(tài)的跳躍現(xiàn)象對振動分析的影響,并給出了一種新的結(jié)構(gòu)優(yōu)化設(shè)計的方法。首先通過特征分析獲得了頻率密集程度不同的3種密頻狀態(tài)系統(tǒng),并給出了接近密頻與接近重頻系統(tǒng)的概念,同時指出密、重頻系統(tǒng)存在模態(tài)跳躍現(xiàn)象;其次對設(shè)計參數(shù)取不同值時所對應(yīng)的多個密、重頻系統(tǒng),利用頻響函數(shù)矩陣計算結(jié)構(gòu)在簡諧激勵下的動力響應(yīng);最后通過響應(yīng)曲線的對比分析,得出有關(guān)模態(tài)跳躍現(xiàn)象對振動分析影響的結(jié)論。

關(guān)鍵詞:振動分析;密集頻率;模態(tài)跳躍;動力響應(yīng);結(jié)構(gòu)優(yōu)化

在實際工程中,很多結(jié)構(gòu)由于自身原因或外界因素的影響,當(dāng)結(jié)構(gòu)參數(shù)發(fā)生微小的變化時,就可能引起結(jié)構(gòu)的部分模態(tài)信息出現(xiàn)急劇變化,即產(chǎn)生模態(tài)跳躍現(xiàn)象,這種現(xiàn)象最早是美國國家航天局(NASA)的專家在彈性薄板構(gòu)件的瞬態(tài)屈曲試驗過程中發(fā)現(xiàn)的,隨后一些學(xué)者做了后續(xù)的研究[1-2]。近年來,這種現(xiàn)象在力學(xué)、土木工程和航空航天等領(lǐng)域逐漸得到關(guān)注,而研究對象也由柔性板等簡單構(gòu)件逐漸向更復(fù)雜結(jié)構(gòu)轉(zhuǎn)變[3-4],同時在密、重頻系統(tǒng)的研究過程中[5-6],很多研究人員發(fā)現(xiàn)頻率密集極可能會引起模態(tài)跳躍現(xiàn)象的發(fā)生[7-8],但這種現(xiàn)象會對結(jié)構(gòu)產(chǎn)生何種影響，目前相關(guān)研究較少。

為了使用模態(tài)分析法來求解結(jié)構(gòu)的精確響應(yīng),必須獲得其固有頻率及模態(tài)信息,這樣解耦的模態(tài)方程才可以很容易地進(jìn)行積分。模態(tài)分析法非常適用于結(jié)構(gòu)動力學(xué)問題的求解,例如在地震載荷的作用下,只需要利用結(jié)構(gòu)的幾個低階模態(tài),就可以近似得到具有較好精度的結(jié)構(gòu)動力響應(yīng)結(jié)果,況且由于高階模態(tài)通常對系統(tǒng)的實際響應(yīng)影響有限,同時有限元法得到的高階模態(tài)又與實際相差很大,因此,從這個意義上講,求解高階模態(tài)的意義不大,但低階模態(tài)是必需的。從上述分析可以看出模態(tài)是計算響應(yīng)的基礎(chǔ),模態(tài)跳躍現(xiàn)象會對結(jié)構(gòu)的動力響應(yīng)產(chǎn)生影響。

隨著靜力條件下的結(jié)構(gòu)優(yōu)化設(shè)計理論的日益成熟及工程設(shè)計要求的逐步提高,結(jié)構(gòu)動力學(xué)優(yōu)化設(shè)計受到越來越多的關(guān)注,結(jié)構(gòu)動力學(xué)優(yōu)化包括結(jié)構(gòu)動力學(xué)特性優(yōu)化和動態(tài)響應(yīng)優(yōu)化2類,相比而言,由于結(jié)構(gòu)動力響應(yīng)的控制屬于多目標(biāo)、多約束優(yōu)化的數(shù)學(xué)模型,因而更為復(fù)雜、難度更大[9-10]。本文結(jié)合一個二自由度阻尼系統(tǒng)的算例,分析模態(tài)的跳躍現(xiàn)象對振動響應(yīng)的影響,并提出一種新的結(jié)構(gòu)優(yōu)化設(shè)計的方法。即首先分析固有頻率隨設(shè)計參數(shù)變化而產(chǎn)生變化的規(guī)律,然后研究當(dāng)固有頻率為密集或重復(fù)狀態(tài)時,相應(yīng)系統(tǒng)的模態(tài)隨之發(fā)生跳躍變化的規(guī)律,最后利用穩(wěn)態(tài)響應(yīng)來對比和分析這些系統(tǒng)的振動特性,從而揭示模態(tài)發(fā)生跳躍時系統(tǒng)響應(yīng)的變化規(guī)律。這樣就可以實現(xiàn)通過討論結(jié)構(gòu)動力響應(yīng)與設(shè)計參數(shù)之間的關(guān)系,來確定在某種動力響應(yīng)要求下的設(shè)計參數(shù)的最優(yōu)值,從而獲得最優(yōu)的設(shè)計方案。

1模態(tài)跳躍現(xiàn)象及其原因

N自由度的線性離散振動系統(tǒng)的自由振動方程為:

(1)

其中,M、C、K分別為系統(tǒng)的質(zhì)量、阻尼和剛度矩陣。其相應(yīng)的強(qiáng)迫振動方程為:

(2)

(3)

設(shè)每個實模態(tài)的規(guī)范化常數(shù)為ai,即

(4)

(K-λiM)φi=0

(5)

實際上特征方程(5)式是關(guān)于矩陣M和K的廣義特征問題。λi稱為固有頻率,與之對應(yīng)的規(guī)范化后的無阻尼模態(tài)φi又可稱為振型或?qū)嵞B(tài)。

當(dāng)結(jié)構(gòu)參數(shù)發(fā)生變化時,研究系統(tǒng)的固有頻率和模態(tài)隨之發(fā)生變化的規(guī)律,既可以使用特征靈敏度作為參考指標(biāo),也可以通過直接研究特征曲線來完成,數(shù)據(jù)處理技術(shù)在這2種方法的應(yīng)用過程中的有效性及可靠性參見文獻(xiàn)[11-12]。本文將使用該方法來實現(xiàn)結(jié)構(gòu)的固有頻率分析及實模態(tài)分析。

1.1　固有頻率分析

考慮一個二自由度的阻尼系統(tǒng),如圖1所示。

圖1　二自由度阻尼振動系統(tǒng)

其中,質(zhì)量、剛度和阻尼矩陣分別為:

(5)式為其無阻尼特征方程,λi(i=1,2)為此系統(tǒng)的特征值,φi(i=1,2)為此系統(tǒng)的特征向量。解得其特征值為:

(6)

其中,p=m1m2;h=-(k1+k2)m2-(k3+k2)m1;c=k1k2+k1k3+k3k2。

令m1=m2=1,k3=1,而取設(shè)計參數(shù)為k1,其變化區(qū)間取為0.6～1.4,分別在k2=0.05,k2=0.005和k2=0.000 5 這3種情況下,利用(6)式得到λ1和λ2關(guān)于設(shè)計參數(shù)k1的關(guān)系圖,并獲得了3個頻率密集程度不同的密頻狀態(tài)系統(tǒng),如圖2所示。

由圖2可知,對于k2取不同值的3種情況下,系統(tǒng)均在k1=1處發(fā)生2個特征值λ1與λ2接近的現(xiàn)象,此時相應(yīng)的系統(tǒng)稱為密頻系統(tǒng),而其余情況下的對應(yīng)系統(tǒng)視為接近密頻系統(tǒng)(或單頻系統(tǒng))。由圖2c可知,第3種密頻狀態(tài)中,當(dāng)k2=0.000 5時,系統(tǒng)在k1=1處出現(xiàn)了重頻現(xiàn)象,即2個特征值相等,此時相應(yīng)的系統(tǒng)稱為重頻系統(tǒng),其余情況下所對應(yīng)的系統(tǒng)視為接近重頻系統(tǒng) (或單頻系統(tǒng))。

圖2　特征值關(guān)于設(shè)計參數(shù)k1的關(guān)系曲線

由圖2可知,k2取值越小,在k1=1處所生成的系統(tǒng)密頻程度越高,甚至出現(xiàn)了重頻現(xiàn)象。

1.2　實模態(tài)分析

相應(yīng)地,按(4)式規(guī)范化后的特征向量為:

i=1,2

(7)

與固有頻率分析的過程類似,令m1=m2=1,k3=1,而取設(shè)計參數(shù)為k1,其變化區(qū)間取為0.6～1.4,分別研究在k2=0.05,k2=0.005和k2=0.000 5時由(7)式得到的φ1和φ2與設(shè)計參數(shù)k1的關(guān)系后發(fā)現(xiàn),當(dāng)k2從0.05向0.000 5變化時,φ1和φ2均在k1=1處兩側(cè)發(fā)生了急轉(zhuǎn)和跳躍,且其變化程度越來越劇烈,其中最為劇烈的情形如圖3所示。這說明φ1與φ2的不穩(wěn)定性均發(fā)生在密頻點k1=1處兩側(cè),系統(tǒng)的密頻程度越高,系統(tǒng)模態(tài)的跳躍性越強(qiáng),換言之,密頻的發(fā)生與模態(tài)跳躍現(xiàn)象呈現(xiàn)出高度的相關(guān)性。

圖3　k2=0.000 5時,φ1、φ2與k1的關(guān)系

為了進(jìn)一步說明這種現(xiàn)象,繪制出φ1和φ2的第1維分量關(guān)于設(shè)計參數(shù)k1的一階導(dǎo)數(shù)圖(或稱為一階靈敏度圖),如圖4所示。事實上,它們的第2維分量關(guān)于設(shè)計參數(shù)k1的一階導(dǎo)數(shù)也大體相同,只是符號相反。

由圖4可知,φ1與φ2在k1∈[0.6,0.9]∪[1.1,1.4]時的變化率幾乎為0,說明此時模態(tài)隨設(shè)計參數(shù)的變化不大,較為穩(wěn)定,但在k1∈[0.9,1.1]時,特別是在k1=1附近,φ1與φ2變化率發(fā)生了較大的波動,并且k2的值越小,變化率的變化幅度越大,這顯示出模態(tài)在這個范圍內(nèi)發(fā)生了一定程度的跳躍,說明密、重頻對靈敏度分析的影響較大,模態(tài)的跳躍現(xiàn)象確實是存在的,而用靈敏度來展示模態(tài)的跳躍性較為直觀有效。

圖4　k2取不同值時,φ11與φ21關(guān)于設(shè)計參數(shù)k1的一階靈敏度

2模態(tài)跳躍對動響應(yīng)分析的影響

根據(jù)上述分析可知,由于設(shè)計參數(shù)變化所產(chǎn)生的密頻系統(tǒng),其密頻程度越高,模態(tài)的跳躍性就越強(qiáng),下面研究和比較模態(tài)的跳躍性對這些系統(tǒng)的動力響應(yīng)分析的影響。

2.1　第1種密頻狀態(tài)下的系統(tǒng)響應(yīng)

對于設(shè)計參數(shù)k2=0.05,k1分別取0.6、0.8、1.2、1.4時所形成的幾種接近密頻的系統(tǒng),同樣使用上述方法便可獲得與圖5相同的簡諧激勵下的穩(wěn)態(tài)響應(yīng)曲線。但是為了便于比較,只是給出時間在0～2.5 s內(nèi)的響應(yīng)曲線,如圖6所示,其他部分的曲線與圖5a的情況類似都是呈周期性變化的,下文與此情況類似的部分不再贅述。

由圖4可知,當(dāng)k2=0.05時,只有設(shè)計參數(shù)k1在密頻點處時,相應(yīng)的密頻系統(tǒng)的實模態(tài)靈敏度最大,即模態(tài)的跳躍性最強(qiáng),但由圖5、圖6對比可知,密頻系統(tǒng)的響應(yīng)幅度較其他接近密頻系統(tǒng)的響應(yīng)而言卻最小,這說明密頻系統(tǒng)相對穩(wěn)定。

圖5　k1=1，k2=0.05時的密頻系統(tǒng)響應(yīng)曲線

圖6　k2=0.05,k1 取不同值時接近密頻的系統(tǒng)響應(yīng)曲線

2.2　第2種密頻狀態(tài)下的系統(tǒng)響應(yīng)

圖7　k1=1,k2=0.005時的密頻系統(tǒng)響應(yīng)曲線

圖8　k2=0.005,k1 取不同值時接近密頻的系統(tǒng)響應(yīng)曲線

由圖4可知,當(dāng)k2=0.005時,只有設(shè)計參數(shù)k1在密頻點處時,相應(yīng)的密頻系統(tǒng)的實模態(tài)靈敏度最大,即模態(tài)的跳躍性最強(qiáng),但通過圖7、圖8的對比,并結(jié)合固有頻率圖2b可知,密頻系統(tǒng)的響應(yīng)幅度最小,這也說明密頻系統(tǒng)相對穩(wěn)定。

2.3　第3種密頻狀態(tài)下的系統(tǒng)響應(yīng)

由于圖2c中呈現(xiàn)出了重頻狀態(tài),其響應(yīng)問題最值得關(guān)注。采用與2.1及2.2中同樣的處理方法得到重頻狀態(tài)及接近重頻狀態(tài)的結(jié)構(gòu)動力響應(yīng)曲線,如圖9、圖10所示。

圖9　k1=1,k2=0.000 5時的重頻系統(tǒng)響應(yīng)曲線

圖10　k2=0.000 5,k1 取不同值時接近重頻的系統(tǒng)響應(yīng)曲線

由實模態(tài)分析圖4可知,當(dāng)k2=0.000 5,只有當(dāng)設(shè)計參數(shù)k1在重頻點處,相應(yīng)的重頻系統(tǒng)的模態(tài)跳躍性最強(qiáng),再通過圖9和圖10對比,并結(jié)合固有頻率圖2c分析可知,此時的重頻系統(tǒng)的響應(yīng)幅度卻最小,說明重頻系統(tǒng)最穩(wěn)定。

圖6、圖8、圖10分別對應(yīng)于圖2a、圖2b、 圖2c中設(shè)計參數(shù)分別取不同值時多個接近密頻的系統(tǒng)響應(yīng)曲線,通過橫向比較可知,振動狀態(tài)區(qū)別不是很大。而圖5、圖7、圖9所對應(yīng)的系統(tǒng)分別對應(yīng)于圖2a、圖2b、 圖2c中設(shè)計參數(shù)k1取最密頻點時,而k2分別取為0.05、0.005及0.000 5所形成的3個系統(tǒng),其密頻程度隨k2減少而增強(qiáng),再由圖4可知,其模態(tài)的跳躍程度也隨之增強(qiáng),因此同時比較這3個系統(tǒng)響應(yīng)顯得尤其必要。從圖5、圖7、圖9中難以發(fā)現(xiàn)細(xì)微差別,因此采用響應(yīng)數(shù)據(jù)來進(jìn)行說明,見表1、表2所列。

表1　k1=1時系統(tǒng)的第1自由度的響應(yīng)值x1(t)

表2　k1=1時系統(tǒng)的第2自由度的響應(yīng)值x2(t)

由表1、表2可知,在各時間點處,頻率越密集的系統(tǒng),其2個自由度的響應(yīng)幅度越小,而重頻系統(tǒng)的響應(yīng)幅度最小,這說明重頻系統(tǒng)更為穩(wěn)定。

從以上研究中可以得出結(jié)構(gòu)優(yōu)化設(shè)計中的一個重要的結(jié)果,即一個設(shè)計參數(shù)k2取所需的值0.05、0.005、0.000 5時,另一個設(shè)計參數(shù)k1只有取1時,系統(tǒng)的響應(yīng)幅度最小,這就是以動力響應(yīng)為優(yōu)化目標(biāo)的最優(yōu)設(shè)計方案。

將上述的研究過程推廣到一般情況,便得到了一種新的結(jié)構(gòu)優(yōu)化設(shè)計方法。即在結(jié)構(gòu)設(shè)計的初始階段,對于某個重要的設(shè)計參數(shù),可以在與其相關(guān)度較大的某個設(shè)計參數(shù)取某些特殊值的前提下,在其可行域內(nèi)計算結(jié)構(gòu)的固有頻率及實模態(tài)的靈敏度,進(jìn)而分析它們隨設(shè)計參數(shù)變化而變化的特點。然后通過計算相應(yīng)結(jié)構(gòu)的響應(yīng),分析其穩(wěn)定性。最后確定滿足某種動力響應(yīng)要求的該設(shè)計參數(shù)的最優(yōu)值，從而獲得最優(yōu)的設(shè)計方案。

3結(jié)論

本文針對由設(shè)計參數(shù)變化所產(chǎn)生的密頻及重頻、接近密頻及接近重頻等多個系統(tǒng),施加相同的簡諧激勵后,對其產(chǎn)生的穩(wěn)態(tài)響應(yīng)進(jìn)行分析,可得到如下結(jié)論:密頻或重頻現(xiàn)象的發(fā)生，對靈敏度分

析所帶來的影響是巨大的，但對結(jié)構(gòu)動力響應(yīng)分析的影響并不大；密頻程度越高時，模態(tài)的跳躍性越強(qiáng)，而系統(tǒng)響應(yīng)的幅度卻越?。恢仡l系統(tǒng)是所有振動系統(tǒng)中最穩(wěn)定的。

同時,本文給出了一種基于動力響應(yīng)的結(jié)構(gòu)設(shè)計的新的優(yōu)化方案。

[參考文獻(xiàn)]

[1]Eduard R,Charle C R,Francis A B.On the solution of mode jumping phenomena in thin-walled shell structures [J].Computer Methods in Applied Mechanics and Engineering,1996,136:59-92.

[2]Everall P R,Hunt G W.Mode jumping in the buckling of structs and plates:a comparatives study [J].International Journal of Non-Linear Mechanics,2000,35:1067-1079.

[3]于巖磊,高維成,劉偉,等.密集模態(tài)結(jié)構(gòu)模態(tài)躍遷分析的簡化攝動法[J].工程力學(xué),2012,29(3):33-39.

[4]蔣友寶,馮健,孟少平.結(jié)構(gòu)損傷識別中模態(tài)躍遷的研究[J].工程力學(xué),2006,23(6):35-40,76.

[5]張淼,于瀾,鞠偉.基于松馳技術(shù)的重頻密頻結(jié)構(gòu)模態(tài)靈敏度分析[J].合肥工業(yè)大學(xué)學(xué)報:自然科學(xué)版,2012,35(12):1605-1609.

[6]張淼,于瀾,鞠偉.重頻系統(tǒng)的頻率靈敏度分析算法研究[J].華南師范大學(xué)學(xué)報:自然科學(xué)版,2014,46(3):39-43.

[7]劉瀟翔,胡軍.包含密集模態(tài)的空間結(jié)構(gòu)的模糊主動振動控制[J].空間控制技術(shù)與應(yīng)用,2010,36(4):18-24.

[8]劉利軍,樊江玲,張志誼,等.密頻系統(tǒng)模態(tài)參數(shù)辯識及其振動控制的研究進(jìn)展[J].振動與沖擊,2007,26(4):109-115.

[9]譚平,卜國雄,劉紅軍,等.帶TMD結(jié)構(gòu)的隨機(jī)地震響應(yīng)分析的新方法[J].北京理工大學(xué)學(xué)報,2010,30(4):390-394.

[10]陳臻林.大型結(jié)構(gòu)動力響應(yīng)的狀態(tài)方程的Krylov精細(xì)時程積分法[J].力學(xué)與實踐,2010,32(2):76-81.

[11]張淼.數(shù)值技術(shù)在結(jié)構(gòu)特征值分析中的應(yīng)用研究[J].長春理工大學(xué)學(xué)報:自然科學(xué)版,2014,37(2):124-127.

[12]張淼.討論密頻結(jié)構(gòu)實模態(tài)參數(shù)的特征曲線法[J].長春師范學(xué)院學(xué)報:自然科學(xué)版,2014,33(2):1-4.

[13]張淼,于瀾,鞠偉.基于頻響函數(shù)矩陣計算阻尼系統(tǒng)動力響應(yīng)的新方法[J].振動與沖擊,2014,33(4):161-166.

(責(zé)任編輯張镅)

Research on the effect of modal jumping phenomenon on vibration analysis

ZHANG Miao1,YU Lan1,JU Wei2

(1.School of Science, Changchun Institute of Technology, Changchun 130012, China; 2.R&D Center, China FAW Group Corporation, Changchun 130011, China)

Abstract:The effect of modal jumping phenomenon on vibration analysis was discussed, and a new method of structural optimization design was proposed. Firstly, three kinds of closed-frequency systems with different closed extents were obtained by eigenproblem analysis, and the concepts of quasi-closed-frequency system and quasi-multiple-frequency system were defined. It was stated that the modal jumping phenomena exist exactly in closed-frequency and multiple-frequency systems. Secondly, in terms of several closed-frequency and multiple-frequency systems corresponding to different values of design parameters, the dynamic responses of them under simple harmonic excitation were computed based on frequency response matrix. Finally, some conclusions about the effect of modal jumping phenomenon on vibration analysis were drawn through contrasting different response curves obtained.

Key words:vibration analysis; closed frequency; modal jumping; dynamic response; structural optimization

doi:10.3969/j.issn.1003-5060.2015.11.024 10.3969/j.issn.1003-5060.2015.11.018

作者簡介：曾偉(1980-),女,四川德陽人,合肥工業(yè)大學(xué)副教授. 汪蓮(1959-),安徽宣城人,合肥工業(yè)大學(xué)副教授,碩士生導(dǎo)師.

基金項目：國家自然科學(xué)青年基金資助項目(71201044)；安徽省級優(yōu)秀青年人才基金重點資助項目(2013SQRW003ZD)和中央高校基本科研業(yè)務(wù)費(fèi)專項資金資助項目(J2014HGXJ0069；2015HGZX0021) 安徽省交通科技進(jìn)步計劃資助項目(2014-35)

收稿日期：2014-11-13；修回日期：2015-01-05 2014-11-13；修回日期：2014-12-24

中圖分類號:O321

文獻(xiàn)標(biāo)識碼：A

文章編號：1003-5060(2015)11-1524-07