含集中參數彈性梁振動特性解析與實驗識別*

含集中參數彈性梁振動特性解析與實驗識別*

王 壯1，2，洪 明2，許俊臣2，崔洪宇2

（1.中國艦船研究設計中心 武漢，430064）（2.大連理工大學船舶工程學院 大連，116024）

基于Laplace變換推導出含集中質量與集中剛度彈性梁的振型函數和典型邊界條件對應的頻率方程，針對有一個集中質量與一個集中剛度的懸臂梁，求解出其固有特性，并利用基于自然激勵技術（natural excitation technique，簡稱NEx T）的特征系統實現算法（eigensystem realization algorithm，簡稱ERA），即NEx T-ERA法對相應結構系統進行了模態識別。通過對比解析結果和實驗結果，分別討論了集中質量與集中剛度大小變化和位置變化時對梁振動特性的影響，得出了集中質量和集中剛度在懸臂梁上位置和大小變化時，懸臂梁固有頻率的相應變化規律，為工程中具有集中質量和集中剛度等直彈性梁的振動分析方法和集中參數布置設計提供了參考。

彈性梁；集中質量；集中剛度；自然激勵技術的特征系統實現算法（NEx T-ERA）；模態識別

引 言

工程中一些結構可以簡化成具有集中質量和集中剛度的等直彈性梁，設計中布置集中質量和集中剛度，追求優良的振動特性是非常有意義的。

具有集中質量與集中剛度的彈性梁的固有振動特性引起了國內外學者的研究興趣［1-2］。Rossit等［3］用解析法推導出了懸臂梁自由端含有一個彈性支撐質量情況下的頻率方程。Banerjee［4］采用動剛度法分析了相同問題。Chang［5］考慮了集中質量轉動慣性的影響，分析了含有一個集中質量簡支梁的振動特性。Low［6］采用特征分析法和Rayleigh法對比分析了僅含有集中質量的歐拉梁的固有頻率。對同樣的模型，Low［7］用Laplace變換的方法推導出了頻率方程，并編寫了計算程序。夏季等［8］采用Laplace變換推導出同時含有任意多個集中質量和集中剛度的彈性梁的振動特性的解析表達式，但表達式較為復雜，只適用于計算含有少量集中參數的情況。彭獻等［9］在Low的基礎上采用符號運算，推導出了便于計算的頻率方程解析表達式，可計算含有集中質量較多情況下的頻率方程，但未考慮含有集中剛度的情況。Wu等［10］采用解析和數值分析相結合的方法求出了含有任意集中質量和集中剛度的歐拉梁的固有頻率和振型。Li等［11］用Ritz法求解了梁兩自由端含有集中質量和彈性支撐的柔性梁的固有特性。以上研究只針對固有特性的求解，沒有探討集中質量和集中剛度大小或者位置布放對梁振動特性的影響，而這個研究認識對指導工程中相應結構減振是非常有必要的。

筆者基于Laplace變換的方法，針對典型邊界條件，推導了含有任意集中質量與集中剛度的彈性梁彎曲振動的固有頻率及振型的解析表達式。采用NEx T-ERA法［12］對含有集中參數的梁結構進行了模態識別實驗，理論解析結果和實驗結果吻合較好。利用驗證的理論算法系統地討論了集中質量和集中剛度的大小以及它們在梁上位置的改變對懸臂梁固有特性的影響，為工程中具有集中質量和集中剛度的等直梁提供了振動分析方法。

1 振型函數及頻率方程推導

對于均勻等直彈性梁，其結構質量與剛度分布是均勻連續的。如果梁上含有集中質量與集中剛度，那么其系統質量與剛度在某一點會產生集中突變，此時可利用δ函數來描述集中質量與集中剛度對均勻梁的振動微分方程的影響，由于方程中δ函數的存在，求解將變得復雜。

1.1 采用δ函數描述集中質量與剛度的梁固有振型函數推導

圖1所示為含有r個集中質量和s個線彈簧的

等直彈性梁，其自由振動微分方程為

圖1 含集中質量與集中剛度彈性梁示意圖Fig.1 Elastic beam with lumped masses and springs

其中：前兩項為均勻直梁的彈性與慣性項；后兩項為集中質量和集中剛度的影響項。

設解為w（xˉ，t）=Y（x）sinωt，采用無量綱形

對式（2）作拉普拉斯逆變換，得到振型函數表達式為

其中：

其中：H（x）為單位階躍函數。

為推導簡便，作如下形式的符號代換，令式（3）中

則滿足以下導數關系

1.2 集中質量和彈簧在梁中位置時梁的頻率方程

為得到頻率方程，考慮梁上含有一個集中質量和一個集中剛度并位于梁中的簡單情況，即式（6）中r=1，s=1，xM1=0.5和xK1=0.5，下面分別推導懸臂梁和簡支梁振動的頻率方程。

懸臂梁的邊界條件為Y（0）=Y＇（0）=0，Y″（1）=Y?（1）=0，將xM1和xK1分別代入振型函數式（6）的對應階導函數式，令x=1，整理可得

其中：

根據邊界條件，式（7）和式（8）可以寫成矩陣形式

由此得到相應的頻率方程

同上推導，即可得出其他簡單邊界條件的頻率方程，這里不再贅述。

1.3 集中質量和彈簧在任意位置時梁的頻率方程

在任意位置處含有一個質量和一個集中剛度的懸臂梁如圖2所示。

同理可得簡支梁的頻率方程為

圖2 集中質量與集中剛度位于梁任意位置示意圖Fig.2 Beam with lumped masses and springs in random position

根據式（6）及懸臂梁的邊界條件得振型函數為

其中：

當ξ1＞ξ2時，得到頻率方程

其中：

當ξ1＜ξ2，則有

2 模態參數識別方法

采用基于自然激勵技術的特征系統實現算法（NEx T-ERA）對梁模態進行識別。

2.1 自然激勵技術

自然激勵技術針對平穩白噪聲激勵下的線性系統，結構上兩點響應信號的互相關或自相關函數與脈沖響應函數具有相近的表達式［13］。因此可應用互相關或自相關函數代替脈沖響應函數數據。

假設一個N自由度的線性系統，若在系統的第l階上作用一個穩態的純白噪聲激勵信號fl（t），則在結構上第n階和第p階的響應信號xnl（t）和xpl（t）的互相關函數［14］為

其中：ψnr為第r階模態振型的第n階成分；λr為第r階特征值；apr為一常量。

根據模態理論，若在系統的第p階自由度上作用一個脈沖激勵力，那么在系統的第n階自由度處的脈沖響應為

其中：bpr為第r階模態的參與因子。

比較式（15）和式（16），發現二者形式相同，僅差兩個常數項apr和bpr。應用NEx T-ERA方法識別系統模態時，可用互相關函數代替脈沖響應函數。

2.2 特征系統實現算法

特征系統實現算法利用結構的脈沖響應函數矩陣構建廣義Hankel矩陣，對其進行奇異值分解得到系統的最小實現，以此識別系統的模態參數［15］。

對N自由度的線性系統，其離散時間的系統狀態空間方程在tk+1時刻表示為

其中：k為離散時間步；x（k）為狀態向量；u（k）= u（ kΔt）為輸入向量；y（k）=y（ kΔt）為輸出向量；A，B和C分別為系統矩陣、輸入矩陣和輸出矩陣。

假設系統在S個自由度上有脈沖輸入，在M個自由度上有傳感器輸出，則ERA算法的數學模型為

設已經得到脈沖響應矩陣h（k），構造廣義Hankel矩陣

其中：P=［C CA…CAα-1］T為αM×2N階可觀矩陣；Q=［B AB…Aβ-1B］為2N×βS階可控矩陣；α和β分別為可觀指數和可控指數。

當k=1時，對H（0）進行奇異值分解，得

將式（18）代入式（19）中，整理可得

其中：Σ=［Σn0］為2N×2N階矩陣。0 0同理，當k=2時，有

設兩個輔助矩陣

其中：IM和0M，IS和0S分別為M階和S階單位陣和零矩陣。

得到系統的最小實現為

對于離散時間系統，設系統矩陣A的特征值矩陣為Z，特征向量矩陣為ψ，對A做特征值分解，得

對于連續時間系統，Ac的特征向量與A的特征向量相同，特征值滿足

其中：zi為特征值矩陣Z的對角元素；λi為Ac的特征值。

根據模態理論，系統的模態頻率、阻尼比和振型為

3 數值算例及模態實驗

3.1 集中質量和集中剛度大小對固有頻率的影響

表1 梁各項參數Tab.1 Material properties of the beam

表2 3種集中質量和彈簧參數Tab.2 Properties of the masses and springs

表3 α=0.075 68和β取不同值時的前4階固有頻率ωiTab.3 The first four natural frequenciesωiwhenα=0.075 68，βhas different values Hz

表4 α=0.123 39和β取不同值時的前4階固有頻率ωiTab.4 The first four natural frequenciesωiwhenα=0.123 39，βhas different values Hz

表5 α=0.185 09和β取不同值時的前4階固有頻率ωiTab.5 The first four natural frequenciesωiwhenα=0.185 09，βhas different values Hz

實驗時彈簧兩端設計有固定裝置，可以安裝在懸臂梁任意位置處，質量塊可用磁座吸附在梁上不同位置。在接近懸臂梁固定端處施加白噪聲激勵，梁上布置有傳感器可測得各點的加速度響應值。實驗儀器及模型示意圖如圖3，4所示。

圖3 實驗模型示意圖Fig.3 The experimental model

圖4 實驗現場Fig.4 Testing ground

分析表3～5可以看出，隨著β值的增大，各階固有頻率ωi逐漸增大，即在集中質量和集中剛度位置不變的前提下，若保持集中質量大小不變，當彈簧剛度增加時，各階固有頻率增大。對比分析表3～5，隨著α值的增大，各階固有頻率ωi逐漸減小，即當彈簧剛度保持不變而集中質量增大時，各階固有頻率呈減小趨勢。通常對于一個單自由度系統，其固有頻率與系統剛度成正比，與系統質量成反比。可見，對含有一個集中質量和一個集中剛度的懸臂梁結構，其固有頻率變化仍符合這一規律。

將實驗結果與理論解對比，很好地驗證了這一規律。根據式（12）計算可得懸臂梁在r=1，s=1，xM1=1和xK1=1/2時振型為

在α=0.123 39和β=21.761 2時，梁前4階模態振型如圖5所示。實驗識別得到的模態振型如圖6所示。對比解析結果與實驗結果，二者十分吻合。

圖5 解析計算得懸臂梁前4階模態振型Fig.5 The analytical first four mode shapes

圖6 實驗識別得懸臂梁前4階模態振型Fig.6 The first four mode shapes obtained from experiment

3.2 集中質量和集中剛度位置對固有頻率的影響

為了探討集中剛度在梁上位置變化時對梁固有頻率的影響，在懸臂梁自由端固定一集中質量，令式（13），（14）中ξ1=1，集中質量系數α1=0.123 39，集中剛度系數β1=21.761 2，當彈簧在梁上的位置為ξ2=0.1，0.2，…，0.9時，分別用本解析方法和實驗識別得到了梁的前4階固有頻率。為便于分析，以集中剛度在梁上的相對位置ξ2為橫坐標，以解析和實驗得到的每一階固有頻率ωi（i=1，2，3，4）為縱坐標，得到ωi的變化曲線如圖7所示。其中，圖左側為解析解縱坐標，右側為實驗值縱坐標。

同理，探討集中質量在梁上位置變化時對梁固有頻率的影響。將彈簧固定在梁中，當ξ2=0.5時，令集中質量在梁上移動，取ξ1=0.1，0.2，…，0.9時，通過本解析法和實驗得到的各階固有頻率ωi（i=1，2，3，4）隨集中質量在梁上相對位置ξ1的變化規律如圖8所示。

圖7 前四階固有頻率隨集中剛度位置的變化曲線Fig.7 The variation of the first four natural frequencies with the spring in different positions

圖8 前四階固有頻率隨集中質量位置的變化曲線Fig.8 The variation of the first four natural frequencies with the mass in different positions

圖9 移動彈簧對梁模態剛度的影響Fig.9 The modal stiffness when the spring is in different positions

對集中質量的位置而言，它對固有頻率的影響和集中剛度相反。如圖10所示，對于一階模態，當集中質量從固定端向自由端移動時，使梁的振動慣性越來越大，導致梁的一階模態質量M1越來越大，因而其一階固有頻率ω1逐漸減小，而二階固有頻率ω2先減小后增大。

圖10 移動質量對梁模態質量的影響Fig.10 The modal mass when the mass is in different positions

對比實驗結果和解析解，實驗得到的各階固有頻率的變化趨勢和解析結果完全吻合，進一步驗證了本研究結論的正確性。

4 結束語

工程中一些結構可以簡化成具有集中質量和集中剛度的等直彈性梁，針對這一模型，通過解析方法得出了其在不同邊界條件下的固有頻率方程和振型函數。通過數值算例求得了含有不同大小集中質量和集中剛度懸臂梁的固有特性解析解，并用NEx TERA法對實際懸臂梁結構的模態進行識別。

將解析結果和實驗結果進行了對比，實驗結果和解析結構互相吻合，驗證了本解析方法的正確性和模態識別方法的有效性。討論了集中質量和集中剛度大小變化以及它們在梁上的位置變化時對梁固有頻率的影響，為工程中具有集中質量和集中剛度的等直彈性梁的振動分析方法和集中參數布置設計提供了參考。

［1］ Gurgoze M.On the eigenfrequencies of a cantilever beam with attached tip mass and a spring-mass system ［J］.Journal of Sound and Vibration，1996，190：149-162.

［2］Lin H Y，Tsai Y C.Free vibration analysis of a uniform multi-span beam carrying multiple spring-mass systems［J］.Journal of Sound and Vibration，2007，302：442-456.

［3］ Rossit C A，Laura P A A.Free vibrations of a cantilever beam with a spring-mass system attached to the free end［J］.Ocean Engineering，2001，28：933-939.

［4］ Banerjee J R.Free vibration of beams carrying springmass systems-A dynamic stiffness approach［J］.Computers and Structures，2012，21-26：104-105.

［5］ Chang C H.Free vibration of a simply supported beam carrying A rigid mass at the middle［J］.Journal of Sound and Vibration，2000，237：733-744.

［6］ Low K H.Frequencies of beams carrying multiple masses：rayleigh estimation versus eigen-analysis solutions［J］.Journal of Sound and Vibration，2003，268：843-853.

［7］ Low K H.On the methods to derive frequency equations of beams carrying multiple masses［J］.International Journal of Mechanical Sciences，2001，43：871-881.

［8］ 夏季，朱目成，馬德毅，帶集中質量和彈性支承梁的橫向固有振動分析［J］.西南工學院學報，1999，14：1-4.Xia Ji，Zhu Mucheng，Ma Deyi，Analysis of lateral natural vibration of beams with lumped masses and elastic supports［J］.Journal of Southwest Institute of Technology，1999，14：1-4.（in Chinese）

［9］彭獻，彭凡，帶集中質量的梁彎曲振動特性新的解析表達式［J］.湖南大學學報，2002，29：44-48.Peng Xian，Peng Fan，New analytical expressions of lateral vibration characteristics of a beam with lumped masses［J］.Journal of Hunan University，2002，29：44-48.（in Chinese）

［10］Wu J S，Chou H M.Free vibration analysis of a cantilever beam carrying any number of elastically mounted point masses with the analytical and numerical combined method［J］.Journal of Sound and Vibration，1998，213：317-332.

［11］Li Jun，Liu Chusheng，Li Feng，et al，Modal analysis of a flexible beam attaching multiple absorbers［J］.Applied Mechanics and Materials，2012，226-228：154-157.

［12］孫熙平，王元戰，趙炳皓，環境激勵下高樁碼頭物理模型模態實驗［J］.振動、測試與診斷，2013，33（2）：263-268.Sun Xiping，Wang Yuanzhan，Zhao Binghao，Modal experiment of physical model for high-piled wharf under ambient excitation［J］.Journal of Vibration，Measurement＆amp；Diagnosis，2013，33（2）：263-268.（in Chinese）

［13］蔣建國，李勤.大型船閘人字閘門工作模態試驗分析［J］.振動、測試與診斷，2008，28（4）：390-394.Jiang Jianguo，Li Qin.Operation modal analysis on miter gates of large lock［J］.Journal of Vibration，Measurement＆amp；Diagnosis，2008，28（4）：390-394.（in Chinese）

［14］Xu Junchen，Hong Ming，Cui Hongyu.The Contrast experimental study on operational modal analysis of ship structural model［J］.Applied Mechanics and Materials，2012，226-228：241-246.

［15］Mohanty P，Rixen D J.Modified ERA method for operational modal analysis in the presence of harmonic excitations［J］.Mechanical Systems and Signal Processing，2006，20：114-130.

O326；O329；TH113

10.16450/j.cnki.issn.1004-6801.2015.04.022

王壯，男，1988年9月生，工程師。主要研究方向為船舶振動噪聲測試與控制。曾發表《含集中質量與彈性支撐壓電層合梁振動控制》（《聲學技術》2015年第34卷第2期）等論文。E-mail：waitorz@qq.com

*國家自然科學基金資助項目（51109034）

2013-12-03；

2014-01-25