摘要： 超長梁結構的振動頻率跨度大，為實現(xiàn)超長梁結構的寬頻振動控制，在梁上周期布置含有阻尼的振幅增強型動力吸振器。振幅放大裝置通過人為地放大受控點處的振幅，從而提高吸振器的工作能力。為能考慮阻尼的影響，建立了復能帶分析模型，并基于人工彈簧模型和能量法提出了一種復能帶計算方法。利用該方法分析了振幅放大裝置形式、吸振器阻尼和放大系數(shù)對復能帶的影響。并研究了非接地式振幅放大裝置連接點相對位置對其減振性能的影響。結果表明，振幅放大裝置形式、吸振器阻尼和放大系數(shù)對復能帶的影響很大；合適的相對位置能夠大幅度提高非接地式振幅放大裝置的工作能力。

關鍵詞： 超材料梁； 動力吸振器； 復能帶； 振幅放大機制； 人工彈簧

中圖分類號： O327； O328 " "文獻標志碼： A " "文章編號： 1004-4523（2025）03-0499-08

DOI：10.16385/j.cnki.issn.1004-4523.2025.03.006

Complex band structure characterization of elastic waves in amplitude enhanced metamaterial beams

GUO Wenjie， HONG Xian， LUO Wenjun， YAN Jianwei， NEI Biao

（State Key Laboratory of Performance Monitoring and Protecting of Rail Transit Infrastructure， East China Jiaotong University， Nanchang 330013， China）

Abstract： Ultra-long beam structure has large vibration frequency span. To control the wide frequency vibration of ultra-long beam structures， amplitude-enhanced dynamic vibration absorbers with damping are arranged periodically on the beam. The amplitude magnification device artificially magnifies the amplitude at the controlled point， thereby increasing the operating ability of the absorber. To be able to consider the effect of damping， a complex band structure analysis model is established， and a new complex band structure calculation method is proposed based on the artificial spring model and the energy method. This method is used to analyze the effect of amplification device type， absorber damping and magnification factor on the complex band structure in detail. The effect of the relative position of the connection points of the ungrounded magnification device on its vibration damping performance is studied. The results show that the magnitude amplification device type， absorber damping and magnification factor have great influence on the complex band structure. The suitably relative position can significantly improve the working ability of the ungrounded magnification device.

Keywords： metamaterial beam； dynamic vibration absorber； complex band structure； amplitude magnification； artificial spring

超長梁結構在工程中具有廣泛的應用，如橋梁、鐵路軌道、管道等。這些超長梁結構受到風荷載、車荷載、地震波等［1?2］動態(tài)荷載會產(chǎn)生振動，有效抑制結構振動帶來的危害至關重要。然而這些振動的頻率通常并不固定，且振動頻率跨度很大，因此超長梁結構的寬頻減振具有重要意義。

動力吸振器（dynamic vibration absorber， DVA）是工程上常用的減振裝置之一，DVA通過與主體結構振動頻率發(fā)生共振轉移和耗散結構中的能量，從而降低結構的振動［3］。振幅放大機構能夠通過放大連接處的位移來提高DVA的吸振效果，杠桿是一種簡易且有效的放大機構。LEE等［4］設計了一種杠桿式DVA，并證明了它的有效性。李春祥等［5］研究了杠桿多重DVA的動力特性，結果表明該裝置能夠提供較好的有效性和魯棒性。楊曉彤等［6］提出了一種含振幅放大機構和慣容的接地剛度動力吸振器，并依據(jù)H∞優(yōu)化準則對吸振器進行了參數(shù)優(yōu)化，結果表明其能夠大幅降低主系統(tǒng)的振幅，拓寬減振頻帶。從上述研究中可知，振幅增強型DVA能實現(xiàn)寬頻減振的目的。

周期結構具有特殊的濾波特性，能夠讓某些頻率波在結構中傳播，而其他頻率波則會產(chǎn)生衰減，該特性也被稱為帶隙特性。LIU等［7］提出的局域共振的概念，拓寬了周期結構運用場景。將吸振裝置作為局域共振器，周期布設于結構中，以抑制波在結構中的傳播。這一想法也應用在梁結構的減振中［8?9］。XIAO等［10］研究了周期布設DVA的梁的傳輸特性，并通過推導和物理模型解釋了周期性系統(tǒng)帶隙的形成機制。這種周期性梁結構也被稱為超材料梁。張垚等［11］設計了基于磁流變彈性體的超材料梁，可以在不改變結構的情況下實現(xiàn)對超材料梁帶隙的調節(jié)。朱學治等［12］將轉動振子周期布置于基體梁上形成超材料梁，并分析了轉動振子的轉動慣量和轉動剛度對帶隙的影響。本文也將利用周期結構的帶隙特性，通過周期布設振幅增強型DVA形成超材料梁，以實現(xiàn)對大跨度梁的寬頻振動調控。

周期結構的帶隙計算方法根據(jù)計算模式的不同可以分為實能帶和復能帶。實能帶是通過給定波數(shù)得到頻率，現(xiàn)今很多方法都是基于該模式，如有限元法［11］、平面波展開法［13］、多重散射法［14］和集中質量法［15］等。而復能帶與實能帶相反，其是通過給定頻率求解波數(shù)。相較于實能帶，復能帶既能得到帶隙的頻率又能獲得衰減特性，有更大的分析優(yōu)勢。但復能帶求解方法更加復雜，因此計算方法也相對有限。如傳遞矩陣法［12］、擴展平面波展開法［16］、擴展微分正交法［17］和擴展有限元法［18］等。傳遞矩陣法是一種天然的復能帶分析方法，但適用性較差，一般僅用于一維周期結構的計算；擴展平面波展開法計算流程直觀，但收斂性差；擴展微分正交法也僅在一維結構中使用；擴展有限元法的適用性很好，但計算結果精度與網(wǎng)格劃分精度相關。上述方法在其使用范圍內都有很好的應用，也有各自的缺陷。本文將提出一種復能帶的分析方法，該方法基于能量法，其計算模式與有限元相似，但無需劃分網(wǎng)格。

在之前的工作中［19］提出了基于人工彈簧模型（artificial spring model，ASM）和能量法的實能帶計算方法。該方法利用人工彈簧模擬周期邊界條件，使得能量法的形函數(shù)無需滿足復雜的Bloch理論，并且將波數(shù)僅存在于周期邊界條件的彈簧勢能中，掃描波數(shù)時無需計算全部矩陣，大大提高了能量法計算帶隙的效率。本文進一步改進ASM，將其擴展到復能帶的計算中。

1 理論計算

1.1 裝置模型介紹

動力吸振器的吸振效果與吸振器的質量和剛度相關［20］，放大受控點的振幅，相當于增加吸振器的質量和剛度，從而提高吸振器的工作能力。文獻［21］分析了振幅放大裝置的工作機制。杠桿是一種常見、簡易和有效的振幅放大（amplitude magnification， AM）裝置。本文研究利用杠桿連接吸振器形成的振幅放大吸振器（AMDVA）。根據(jù)杠桿支撐的方式不同，可以分為接地式和非接地式。

1.1.1 接地式AMDVA模型

如圖1所示，接地式AMDVA的支撐點和地面（固定點）相連，接地式AMDVA能夠獲得很好的減振效果。其中，質量塊提供AMDVA的質量mDVA，橡膠層提供AMDVA的剛度kDVA和阻尼cDVA，振幅放大系數(shù)表示為α=l2/l1，其中，l1為支撐點到受控點的桿長，l2為吸振器與桿的連接處到支撐點的桿長。假設受控點處基體梁的垂向位移為w1，此時DVA處的位移為αw1。

1.1.2 非接地式AMDVA模型

如圖2所示，非接地式AMDVA的受控點和支撐點均在連接的基體梁上，因此非接地式AMDVA在裝配上更加靈活。設支撐點處的垂向位移為w2，那么傳導到DVA上的振幅則為α（w1-w2）。

1.2 基于ASM的復能帶計算

提出的基于ASM和能量法的實能帶計算方法［19］，能夠簡單高效地得到結構的實能帶曲線。在工程中，往往需要得到結構的衰減特性，而這是實能帶所不能做到的。因此，建立復能帶分析模型，能更全面地分析周期結構的振動特性。本文將基于ASM和能量法提出一種復能帶計算方法。該方法的基本計算思路是：利用ASM模擬周期邊界條件，波數(shù)k僅存在人工彈簧的彈性勢能中。通過處理ASM的剛度矩陣，將波數(shù)k和剛度矩陣解耦，進而給定頻率ω求解波數(shù)k。

如圖3所示，在梁下周期布設AMDVA，間隔為l，梁的寬度和高度分別為b和h。梁的垂向振動基于歐拉梁理論。根據(jù)能量法，梁的垂向振動可以由與時間相關的未知系數(shù)a（t）和與位置相關的形函數(shù)f（x）組成：

2 數(shù)值分析

2.1 模型參數(shù)

在本文中，梁的材料為鋼，密度為7850 kg/m3，彈性模量為210 GPa，泊松比為0.3，梁的高度和寬度均為0.03 m。AMDVA的布設間距為0.6 m，質量塊和一個周期單元梁的質量比為0.03，吸振頻率fDVA=500 Hz，剛度取值采用單自由度的固有頻率計算公式得到：k_DVA=m_DVA （2πf_DVA ）^2。AMDVA的受控點位于一個周期單元的中間，受控點到支撐點的桿長為l1=0.05 m，支撐點到DVA連接處的桿長為l2=0.15 m。在本節(jié)中，采用接地式AMDVA作為研究對象，暫不考慮阻尼的影響。后續(xù)分析中，除非另外說明，否則參數(shù)保持不變。

2.2 收斂性分析

在使用基于ASM的能量法時，復能帶的計算精度與模擬位移的形函數(shù)個數(shù)和虛擬彈簧的剛度相關。本節(jié)將利用控制變量法，研究二者取值的收斂性。

2.2.1 形函數(shù)個數(shù)

本文方法是利用N個形函數(shù)來模擬位移，N的取值也決定了復能帶的精度。將位移和轉角的人工彈簧剛度均取為1014，其單位分別為N/m和N·m/rad。以頻率為2000、 4000、 6500 Hz時，波數(shù)的實數(shù)部分取值為參考。圖4為形函數(shù)個數(shù)的收斂曲線。從圖4中可以看到，當形函數(shù)個數(shù)取為14時，6500 Hz以內的復能帶就可以收斂。

2.2.2 虛擬彈簧剛度

正如1.2節(jié)中所述，使用ASM模擬周期邊界時，虛擬彈簧的剛度取值為無窮大，因此在計算過程中，需要用一個“大”值來代替。ASM計算結果的準確性與剛度密切相關，因此需要分析虛擬彈簧剛度取值的收斂性。取形函數(shù)個數(shù)N=15，同樣以頻率為2000、 4000、 6500 Hz時，波數(shù)的實數(shù)部分取值為參考。圖5為虛擬彈簧剛度取值的收斂曲線。從圖5中可以看出，當虛擬彈簧剛度取值為1013 N/m時，6500 Hz以內的復能帶結構就可以收斂。綜上，在后續(xù)的分析中，將形函數(shù)個數(shù)為15，虛擬彈簧剛度取值為1013。

2.3 準確性驗證

COMSOL能夠提供周期邊界條件，得到周期結構的帶隙特性。本節(jié)將通過與有限元計算的結果進行對比，研究本文方法的準確性。利用COMSOL的固體力學模塊建立如圖3所示的有限元模型，并利用特征頻率研究得到結構的復能帶。圖6為COMSOL的有限元模型，圖中的藍色部分為基體梁，下方的紅色部分為剛性桿，桿的末尾黑色部分為吸振器，黑點的位置為彈簧阻尼器的連接位置。如圖7所示，梁和桿之間采用彈簧阻尼器連接，彈簧阻尼器的垂向剛度設置為1015 N/m，以模擬剛性連接；支撐點處采用彈性基礎模擬，彈性基礎各向彈簧均為1015 N/m，以模擬剛性固定約束；質量塊和桿之間采用彈性薄層連接，剛度取值與吸振器剛度取值一致。

可以看出FEM的實能帶結果和ASM計算得到的復能帶實部結果在4000 Hz以內吻合地較好，而對于4000 Hz以上，因為梁模型選用的不同會產(chǎn)生誤差，其中最大的誤差為4.04%，可以認為該段誤差在可接受范圍內，由此可以證明ASM的正確性。從圖8（a）中可以看出，設有AMDVA的周期性梁的前三段帶隙分別為：155～195 Hz、481～615 Hz、782～907 Hz。根據(jù)Bragg散射機理，當波長滿足條件：2l=nr（n=1， 2， 3，…），其中，r為波長，l為單元長度。Bragg散射頻率的估算公式為：f_Bragg=√（EI/ρA） "（n^2 π）/（2l^2 ），將2.1節(jié)中的參數(shù)代入此公式中，當n=1時，fBragg=195.4 Hz；當n=2時，fBragg=781.7 Hz。由此可知155～195 Hz和782～907 Hz為Bragg帶隙。圖8（b）為結構的衰減特性曲線，可以看出由AMDVA共振產(chǎn)生的帶隙具有最大的衰減。而由Bragg散射產(chǎn)生的帶隙的衰減相對較小。

3 復能帶特性研究

復能帶既能得到周期結構的振動傳播規(guī)律，還能獲得衰減特性，相較于實能帶具有更大的分析優(yōu)勢。本文將分析分別布設了接地式AMDVA和非接地式MDVA的周期性梁的復能帶特性。

3.1 裝置形式的復能帶

在1.1節(jié)中介紹了接地式和非接地式的AMDVA，并分析了兩者的特點，兩種裝置的不同結構形式的減振效果也需要進行定量分析。在實際工程中往往會加入阻尼來提升DVA的減振效果，因此在本節(jié)中，設置阻尼系數(shù)cDVA=100 N/（m/s）。

設置的振幅放大吸振器的吸振頻率為500 Hz，因此在分析復能帶時，分析頻段范圍為0～1500 Hz。由圖9可知，AMDVA和傳統(tǒng)DVA在吸振頻率500 Hz均會產(chǎn)生帶隙。接地式AMDVA具有遠超非接地式AMDVA的吸振頻率與吸振效果，并且非接地式AMDVA比傳統(tǒng)DVA有更好的工作能力。從圖9（a）中可以看出，當DVA中含有阻尼后，復能帶的實部很難正確地反映出振動的傳播規(guī)律。因此唯有復能帶的虛部才能準確地分析出含阻尼DVA的振動傳播與衰減規(guī)律。

3.2 阻尼對復能帶的影響

阻尼的耗散作用對工程結構的減振具有重要意義，而這是實能帶所不能準確描述的。在上一節(jié)的分析中也看出復能帶實部并不能準確地描述振動的傳播與衰減，本節(jié)將利用復能帶虛部研究阻尼的影響。

如圖10（a）所示，無阻尼的非接地式AMDVA僅在500 Hz左右小范圍內有強烈衰減，而添加阻尼的非接地式AMDVA能夠擴大減振頻率，但也降低了衰減程度。在190～196 Hz有一段很小的衰減區(qū)域，這段帶隙區(qū)域幾乎不受阻尼的影響，這段帶隙是由Bragg散射產(chǎn)生，而吸振器對Bragg散射的影響很小，也就導致阻尼對該段頻率的影響不大。圖10（b）中的154～196 Hz處的帶隙不變的原因與非接地式AMDVA一致。而且因為接地式AMDVA的良好放大效果，梁結構產(chǎn)生了寬頻帶隙且有較好的衰減。綜上所述，阻尼的滯后效應拓寬了局域共振帶隙，并降低了衰減程度；原來的通帶頻率會產(chǎn)生衰減，且阻尼越大衰減程度越大。

3.3 放大系數(shù)對復能帶的影響

振幅放大裝置的放大系數(shù)很大程度上決定了AMDVA減振效果。放大系數(shù)越大，即對杠桿的材料和空間要求更高，同時也具有傾覆的風險。對于放大系數(shù)的分析是十分必要的。圖11為放大系數(shù)對復能帶虛部的影響，其中云圖的顏色表示復能帶振動的衰減程度，顏色越紅表示衰減越大，顏色越藍表示衰減越小，云圖顏色條上的值表示復能帶的虛部取值。

從圖11中可以看出，隨著放大系數(shù)的增大，這兩類AMDVA的減振頻率和減振能力均會得到增強，且接地式AMDVA的工作能力具有明顯的增強。對于非接地式AMDVA，其最大的衰減頻率始終維持在500 Hz。而接地式AMDVA的最大衰減頻率會隨著放大系數(shù)的增加而增加。其原因在于：放大系數(shù)會增大DVA的等效質量，從而導致DVA的等效質量相較于主體結構變得不可忽視，單自由度的計算方法不再準確。

4 受控點和支撐點之間的相對位置

對非接地式AMDVA的影響

非接地式AMDVA的布設靈活，能夠很好地適應工程實際，但在上述研究中也發(fā)現(xiàn)，非接地式AMDVA的減振效果相對有限，如何在不改變DVA參數(shù)的前提下提升非接地式AMDVA減振效果是本節(jié)的研究重點。受控點和支撐點的相對位置會影響非接地式AMDVA的工作能力。因為結構具有周期性，因此本節(jié)將受控點置于一個周期單元的起點，通過移動支撐點的位置來研究受控點和支撐點之間的相對位置對非接地式AMDVA的影響，并設受控點和支撐點的相對位置為Δx。將Δx從0增大到0.6，結果如圖12所示，云圖的顏色條上的值表示復能帶的虛部取值。

從圖12中可以看出：受控點和支撐點的相對位置會顯著影響非接地式AMDVA的減振效果，也就意味著調節(jié)兩點的相對位置能夠提高非接地式AMDVA的工作能力；當兩連接點的相對位置處于0.2～0.35時，非接地式AMDVA具有較寬的減振頻率以及很好的衰減率。

提取Δx=0.6和Δx=0.33時非接地式AMDVA的復能帶虛部，如圖13所示?？傮w上，Δx=0.33時非接地式AMDVA的減振效果比Δx=0.6時的更加優(yōu)異。而且，Δx=0.33時，非接地式AMDVA的減振效果甚至超過了接地式AMDVA。為研究產(chǎn)生該現(xiàn)象的原因，繪出振動頻率為500 Hz時的基體梁的振型圖，如圖14所示。當Δx=0.6時，受控點和支撐點之間的位移差很小，意味著放大效果不明顯，也就導致減振效果提升不夠；當Δx=0.33時，兩點之間的位移差更大，也就有著更好的減振效果。為更好地使用非接地式AMDVA，在確定吸振頻率后，可以根據(jù)基體結構的振型布設受控點和支撐點，以獲得更好減振效果。

5 結 "論

為降低梁結構的振動，本文將兩種振幅放大型吸振器分別周期布設于梁上，并從復能帶的角度研究兩種裝置的減振效果。本文基于ASM和能量法，提出了一種復能帶的計算方法，并詳細說明了該方法的計算流程。通過控制變量法，確定了本方法虛擬彈簧剛度取值和形函數(shù)個數(shù)的收斂性，并利用FEM驗證了提出方法的正確性。之后研究了兩種裝置形式、材料阻尼和振幅放大系數(shù)對復能帶的影響。得到了以下結論：

（1）接地式AMDVA和非接地式AMDVA比傳統(tǒng)DVA有著更好的減振寬度和衰減性能。

（2）復能帶能夠準確地反映阻尼對振動傳遞的影響，阻尼能夠拓寬局域共振帶隙，并降低其衰減程度；而原本的通帶頻率在阻尼作用下會產(chǎn)生衰減，且阻尼越大衰減程度越大。

（3）放大系數(shù)的增加會提高AMDVA的工作能力，且對接地式AMDVA的工作能力提高得更加明顯。

（4）非接地式AMDVA的減振效果與受控點和支撐點的相對位置相關，通過調整相對位置，能夠大幅度提高其工作能力。

參考文獻：

［1］ 王少欽， 夏禾， 郭薇薇， 等. 考慮橋梁幾何非線性的風-車-橋耦合振動分析［J］. 工程力學， 2013， 30（4）： 122-128.

WANG Shaoqin， XIA He， GUO Weiwei， et al. Coupling vibration analysis of wind-train-bridge system considering geometric nonlinearity of bridge［J］. Engineering Mechanics， 2013， 30（4）： 122-128.

［2］ 王亞偉， 鄭凱鋒， 熊籽躒， 等. 地震與風聯(lián)合作用下大跨橋梁車-橋耦合振動分析［J］. 中國公路學報， 2021， 34（2）： 298-308.

WANG Yawei， ZHENG Kaifeng， XIONG Ziluo， et al. Coupled vibration analysis of vehicle-bridge for long-span bridge under wind and earthquake action［J］. China Journal of Highway and Transport， 2021， 34（2）： 298-308.

［3］ SUN C， NAGARAJAIAH S. Study on semi-active tuned mass damper with variable damping and stiffness under seismic excitations［J］. Structural Control and Health Monitoring， 2014， 21（6）： 890-906.

［4］ LEE E T， EUN H C. Lever-type tuned mass damper for alleviating dynamic responses［J］. Advances in Civil Engineering， 2019（1）： 5824972.

［5］ 李春祥， 王肇民. 杠桿式多重調諧質量阻尼器的動力特性［J］. 固體力學學報， 2003，24（4）： 451-455.

LI Chunxiang， WANG Zhaomin. Dynamic characteristics of lever-type multiple tuned mass dampers［J］. Chinese Journal of Solid Mechanics， 2003，24（4）： 451-455.

［6］ 楊曉彤， 申永軍， 王俊鋒. 一種含放大機構、慣容和接地剛度的動力吸振器的參數(shù)優(yōu)化［J］. 振動與沖擊， 2022， 41（21）： 308-315.

YANG Xiaotong， SHEN Yongjun， WANG Junfeng. Parametric optimization of a DVA with amplifying mechanism， inerter and grounded stiffness［J］. Journal of Vibration and Shock， 2022， 41（21）： 308-315.

［7］ LIU Z Y， ZHANG X X， MAO Y W， et al. Locally resonant sonic materials［J］. Science， 2000， 289（5485）： 1734-1736.

［8］ LIU Y Z， YU D L， LI L， et al. Design guidelines for flexural wave attenuation of slender beams with local resonators［J］. Physics Letters A， 2007， 362（5-6）： 344-347.

［9］ YU D L， LIU Y Z， WANG G， et al. Flexural vibration band gaps in Timoshenko beams with locally resonant structures［J］. Journal of Applied Physics， 2006， 100（12）： 124901.

［10］ XIAO Y， WEN J H， YU D L， et al. Flexural wave propagation in beams with periodically attached vibration absorbers： band-gap behavior and band formation mechanisms［J］. Journal of Sound and Vibration， 2013， 332（4）： 867-893.

［11］ 張垚， 李金強， 宋智廣， 等. 利用磁流變彈性體對超材料梁進行帶隙調控［J］. 哈爾濱工程大學學報， 2022， 43（9）： 1271-1276.

ZHANG Yao， LI Jinqiang， SONG Zhiguang， et al. Bandgap control of metamaterial beam using magnetorheological elastomers［J］. Journal of Harbin Engineering University， 2022， 43（9）： 1271-1276.

［12］ 朱學治， 陳照波， 焦映厚， 等. 含有周期分布轉動振子的聲子晶體梁的彎曲振動帶隙研究［J］. 振動與沖擊， 2017， 36（21）： 85-91.

ZHU Xuezhi， CHEN Zhaobo， JIAO Yinghou， et al. Flexural vibration band gaps in a phononic crystal beam containing rotational resonators［J］. Journal of Vibration and Shock， 2017， 36（21）： 85-91.

［13］ DAL POGGETTO V F， SERPA A L. Elastic wave band gaps in a three-dimensional periodic metamaterial using the plane wave expansion method［J］. International Journal of Mechanical Sciences， 2020， 184： 105841.

［14］ 房學謙， 胡超， 杜善義. 半無限功能梯度材料結構中圓孔對彈性波的多重散射［J］. 工程力學， 2007，24（12）： 14-18.

FAGN Xueqian， HU Chao， DU Shanyi. Multiple scattering of elastic waves in semi-infinite functionally graded material with a circular cavity［J］. Engineering Mechanics， 2007，24（12）： 14-18.

［15］ 溫激鴻， 王剛， 劉耀宗， 等. 基于集中質量法的一維聲子晶體彈性波帶隙計算［J］. 物理學報， 2004，53（10）： 3384-3388.

WEN Jihong， WANG Gang， LIU Yaozong， et al. Lumped-mass method on calculation of elastic band gaps of one-dimensional phononic crystals［J］. Acta Physica Sinica， 2004，53（10）： 3384-3388.

［16］ LAUDE V， ACHAOUI Y， BENCHABANE S， et al. Evanescent Bloch waves and the complex band structure of phononic crystals［J］. Physical Review B， 2009， 80（9）： 092301.

［17］ CHENG Z B， SHI Z F， MO Y L. Complex dispersion relations and evanescent waves in periodic beams via the extended differential quadrature method［J］. Composite Structures， 2018， 187： 122-136.

［18］ VERES I A， BERER T. Complexity of band structures： Semi-analytical finite element analysis of one-dimensional surface phononic crystals［J］. Physical Review B， 2012， 86（10）： 104304.

［19］ 馮青松， 楊舟， 郭文杰， 等. 基于人工彈簧模型的周期結構帶隙計算方法研究［J］. 力學學報， 2021， 53（6）： 1684-1697.

FENG Qingsong， YANG Zhou， GUO Wenjie， et al. Research on band gap calculation method of periodic structure based on artificial spring model［J］. Chinese Journal of Theoretical and Applied Mechanics， 2021， 53（6）： 1684-1697.

［20］ ?ZYAR O， YILMAZ ?. A self-tuning adaptive-passive lever-type vibration isolation system［J］. Journal of Sound and Vibration， 2021， 505： 116159.

［21］ FENG Q S， YANG Z， GUO W J， et al. Enhanced suppression of vibrational wave transmission in structures with periodically attached absorbers exploiting amplitude magnification mechanisms［J］. European Journal of Mechanics-A/Solids， 2022， 95： 104666.

第一作者： 郭文杰（1991―），男，博士，副教授。

E-mail： 739633869@qq.com

通信作者： 羅文?。?979―），女，博士，教授。

E-mail： lwj06051979@163.com