李宗威，紀(jì)剛，周其斗

海軍工程大學(xué) 艦船與海洋學(xué)院，湖北 武漢 430033

0 引 言

多自由度系統(tǒng)是指在任意時候都需要2 個或者更多的廣義坐標(biāo)才能完全確定其位置的系統(tǒng)。工程結(jié)構(gòu)與多跨梁具有相似的特征：每個單元都表現(xiàn)為多自由度振動，屬于多自由度特征的系統(tǒng)。開展多自由度系統(tǒng)的局域化研究對于工程振動控制設(shè)計具有特別重要的意義。此外工程結(jié)構(gòu)都是有限尺寸的結(jié)構(gòu)，需要將無限結(jié)構(gòu)的波傳遞響應(yīng)推廣至有限尺寸結(jié)構(gòu)的振動傳遞響應(yīng)。通過有界媒質(zhì)中的強迫振動分析[1]，表明有界媒質(zhì)的振動響應(yīng)可視為由擾動源發(fā)出的左、右傳播波及其在邊界多次反射波疊加的結(jié)果。有界媒質(zhì)中，遠離擾動源處的響應(yīng)特征被歸結(jié)為無界媒質(zhì)中自由行進波的特征。

周期系統(tǒng)在不同的頻段具有不同的振動特性。振動傳遞若在某些頻段內(nèi)不隨傳播距離的增加而衰減，則該頻段稱為通頻帶；而在另一些頻段，振動隨著傳播距離的增加而呈指數(shù)衰減，此時，振動能量被局限于振源附近，即振動被“局域化”，該頻段稱為止頻帶。與周期系統(tǒng)不同，非周期系統(tǒng)的振動在所有頻段都表現(xiàn)為局域化特征，即在所有頻帶內(nèi)，振動傳遞隨傳播距離的增加而衰減，因而不存在通頻帶。1958 年，固態(tài)物理學(xué)家Anderson[2]在研究不規(guī)則條件下電子傳播規(guī)律時，首次發(fā)現(xiàn)了局域化現(xiàn)象。Hodges[3]通過對耦合振子鏈模型的研究，首次證實了在結(jié)構(gòu)動力學(xué)領(lǐng)域中也存在局域化效應(yīng)，發(fā)現(xiàn)無序結(jié)構(gòu)振動能量局限于激勵源附近，振動幅值隨傳播距離的增加呈指數(shù)衰減。Pierre 等[4-7]對有限彈簧振子耦合擺模型進行分析，從模態(tài)角度分析解釋了無序結(jié)構(gòu)的局域化效應(yīng)，說明振動局域化效應(yīng)在模態(tài)上體現(xiàn)為局部模態(tài)，并給出了振動強局域化發(fā)生的條件。國內(nèi)學(xué)者劉文璽和譚路等[8-9]通過具體的實例研究了圓柱殼結(jié)構(gòu)不等間距布置對圓柱殼振動特性與聲學(xué)性能的影響，發(fā)現(xiàn)結(jié)構(gòu)不等間距布置具有一定的減振作用。

從理論分析方法上看，上述學(xué)者多采用波動分析方法，因為這些模型結(jié)構(gòu)大多是一維波傳播模型，動力學(xué)表達簡單，容易通過理論解析定量給出局域化因子。但是，若直接采取理論分析方法對具有多自由度特征的實際工程系統(tǒng)進行振動局域化的研究，尚有一定的困難：一是波動控制方程在表達上比較復(fù)雜，因為解析方法給出的傳遞遞歸關(guān)系需要引入假定，在理論推導(dǎo)上存在一定的困難；二是較難表示有限尺寸結(jié)構(gòu)的振動場；三是對于平板、圓柱殼等具有二維振動特征的結(jié)構(gòu)無法直接運用針對一維模型得出的理論。

為解決多自由度系統(tǒng)理論分析存在的難題，建立一種等效的方法，即將多自由度系統(tǒng)的振動等效為彈簧質(zhì)量鏈系統(tǒng)的振動，然后利用彈簧質(zhì)量鏈系統(tǒng)的局域化公式預(yù)報多自由度系統(tǒng)局域化因子，進而開展局域化規(guī)律的研究。該方法可以有效規(guī)避對多自由度系統(tǒng)進行復(fù)雜解析的困難，實現(xiàn)對多自由度系統(tǒng)開展局域化定量分析。

本文將以梁模型為研究對象開展上述等效方法的研究。首先，介紹等效方法的原理； 然后，利用等效方法對一維梁系統(tǒng)進行局域化分析，預(yù)報局域化因子； 最后，通過比較等效方法與統(tǒng)計擾動法的計算結(jié)果，以及與具體實例的有限元模型衰減特征對比，來驗證等效方法的有效性。

1 單通帶近似法

對于多自由度周期系統(tǒng)，波數(shù)kw與頻率ω 的色散曲線可用圖1（a）所示的曲線表示。與單自由度系統(tǒng)具有唯一的通頻帶不同，多自由度系統(tǒng)具有多個通頻帶。Photiadis 和Houston[10]指出對于具有多個通頻帶的系統(tǒng)，可以通過單通帶近似（single band approximation，SBA）的等效方法進行簡化分析。通過采用彈簧質(zhì)量鏈系統(tǒng)來代替實際系統(tǒng)，將結(jié)構(gòu)在各階通頻帶的振動視為固有頻率不同的振子系統(tǒng)來進行研究，即分析單一通頻帶。因此，只要給出了相應(yīng)的振動等效參數(shù)，再根據(jù)彈簧質(zhì)量鏈系統(tǒng)的局域化理論結(jié)果，就可以預(yù)報實際模型的局域化因子。

單通帶近似法[10]的本質(zhì)在于將多自由度系統(tǒng)的色散曲線中的某一根色散曲線視為等效耦合彈簧質(zhì)量鏈系統(tǒng)的色散曲線。耦合彈簧質(zhì)量鏈系統(tǒng)的色散曲線可以表示為圖1（b）所示曲線，可根據(jù)曲線的特征給出該通頻帶的帶寬和該通頻帶的中心頻率。若圖中通頻帶的上、下限頻率為 ω+和ω- ，則可以得出等效的耦合度 V 和，其中 ω0為中心頻率。

針對如圖2 所示的彈簧質(zhì)量鏈系統(tǒng)進行分析，其中彈簧振子質(zhì)量為 m，振子間耦合剛度為 k，振子對地剛度為 kj。

圖 1 波數(shù)與頻率的色散曲線Fig. 1 Dispersion curve of wave number and frequency

圖 2 彈簧質(zhì)量鏈系統(tǒng)Fig. 2 Mass-spring chain system

以第 j 個彈簧質(zhì)量振子為例，可以得到其運動關(guān)系式為

2 等間距梁模型的振動參數(shù)識別與等效

等間距梁模型如圖3 所示。梁跨距為l0，跨間扭矩彈簧剛度為 C，彈性模量為E，梁的線密度為m0，梁的截面面積慣性矩為I。梁的有限元模型結(jié)構(gòu)尺寸和材料屬性如表1 所示。

圖 3 等間距梁模型Fig. 3 Model of beam with equal spacing

表 1 有限元梁模型結(jié)構(gòu)相關(guān)參數(shù)Table 1 Structural parameters of finite element beam model

當(dāng)K=1 時，等間距梁的色散曲線如圖4 所示。圖中，實線表示梁的通頻帶。由圖可以看出，與彈簧質(zhì)量鏈系統(tǒng)具有唯一的通頻帶不同，梁系統(tǒng)具有多個通頻帶。

圖 4 耦合梁模型的通頻帶及色散關(guān)系（K=1）Fig. 4 The pass band and dispersion relations of coupled beam model（K=1）

通過單通帶近似法將梁的某一根色散曲線視作等效耦合振子鏈的色散曲線，如圖5 所示。

圖 5 等間距梁的色散曲線Fig. 5 Dispersion curves of beam with equal sapcing

針對某一具體通頻帶，通過上、下限頻率ω+和 ω-，可以計算得到相應(yīng)的振動等效參數(shù)。

3 不等間距梁的局域化因子預(yù)報

4 等效方法有效性的驗證

通過對比統(tǒng)計擾動方法與等效方法計算的局域化因子結(jié)果以及通過不等間距梁有限元模型實例來驗證等效方法預(yù)報結(jié)果的有效性 。

4.1 等效方法預(yù)報結(jié)果與統(tǒng)計擾動方法結(jié)果比較

文獻[7]中，由統(tǒng)計擾動方法分析得到的局域化因子為

以K=1， σ2=0.005 77 ，的 計 算 模 型 為 算例，計算出各通頻帶的帶寬和中心頻率，以此作為輸入，利用式（17）計算等效局域化因子。

利用統(tǒng)計擾動法和單通帶近似法計算出的局域化因子結(jié)果對比如圖6 所示。其中：散點為等效方法計算得到的局域化因子，實線為統(tǒng)計擾動方法計算的局域化因子值。從圖中可以看出，通過單通帶近似的方法將梁系統(tǒng)等效為彈簧質(zhì)量鏈系統(tǒng)預(yù)報得到的局域化因子同統(tǒng)計擾動方法給出的結(jié)果具有相同的量級，而且規(guī)律基本符合：隨著頻率的增加，局域化因子具有增加的趨勢。這說明了運用等效方法預(yù)報多自由度系統(tǒng)的局域化因子具有可行性。

圖 6 統(tǒng)計擾動法和單通帶近似法計算的局域化因子結(jié)果對比Fig. 6 Local factors comparison calculated by statistical perturbation method and single pass band approximation method

4.2 不等間距梁有限元模型驗證

圖 7 梁的有限元模型Fig. 7 Finite element model of the beam

圖 8 振動幅值同理論局域化因子衰減線比較Fig. 8 Comparison of the vibration amplitude by numerical and theoretical results

5 結(jié) 論

本文以梁模型為研究對象，給出了預(yù)報不等間距梁局域化因子的等效方法，驗證了預(yù)報方法的有效性，得到如下主要結(jié)論：

1）對不等間距梁的色散曲線進行分析，運用單通帶近似方法將單一的色散曲線等效為耦合彈簧質(zhì)量鏈系統(tǒng)的色散曲線，采用參數(shù)等效手段將梁的振動參數(shù)等效為彈簧質(zhì)量鏈系統(tǒng)的振動參數(shù)，基于波動理論給出不等間距梁的無序度參數(shù)，進而利用彈簧質(zhì)量鏈模型的局域化因子公式預(yù)報出不等間距梁的局域化因子。該方法有效規(guī)避了對多自由度系統(tǒng)進行復(fù)雜解析的困難，實現(xiàn)了對多自由度系統(tǒng)開展無序局域化定量分析的工程化使用。

2）通過對比統(tǒng)計擾動方法與單通帶近似等效方法的局域化因子結(jié)果驗證了等效方法的正確性。

3）通過對比不等間距梁FEM 模型證實了采用不等間距布置實現(xiàn)振動局域化的可行性，驗證了采用單通帶近似等效方法對多自由度系統(tǒng)進行局域化分析的有效性。