管道超結(jié)構(gòu)軸向帶隙特性分析及實驗研究

李孟昶，郭少杰,張紅艷

(1.長安大學(xué)理學(xué)院，西安 710064；2.西安交通大學(xué)機(jī)械結(jié)構(gòu)強(qiáng)度與振動國家重點(diǎn)實驗室，西安 710049)

0 引 言

近年來，超材料逐漸引起研究人員的關(guān)注[1]，并被應(yīng)用到多個領(lǐng)域，如聲學(xué)超材料、光學(xué)超材料、電磁超材料等[2-5]。聲學(xué)超材料典型的特征是具有亞波長的尺寸和異常的動態(tài)等效參數(shù)[6]，當(dāng)彈性波在這類結(jié)構(gòu)中傳播時會受到結(jié)構(gòu)的調(diào)制，從而使彈性波在某些頻率范圍內(nèi)無法傳播。

國內(nèi)外學(xué)者最初的研究主要集中于布拉格(Bragg)散射型帶隙[3]。2000年武漢大學(xué)劉正猷教授提出了局域共振型帶隙結(jié)構(gòu)[4]。通常帶隙的形成是結(jié)構(gòu)周期性和單散射體共同作用的結(jié)果，在布拉格散射型帶隙中，各元胞之間相互作用起主導(dǎo)作用，而在局域共振帶隙中，單個散射體的共振特性起主導(dǎo)作用[7-8]。許多學(xué)者對一維聲學(xué)超材料的彈性波帶隙特性展開研究。溫激鴻等[9]通過將一維聲子晶體中的元胞簡化為集中質(zhì)量，引入了一種基于集中質(zhì)量的彈性波集中質(zhì)量計算方法。Xu等[10]用集中質(zhì)量法研究了一維聲子晶體的階頻帶隙特性，分析了集中質(zhì)量數(shù)、元胞尺寸、彈性模量等因素對帶隙的影響。郁殿龍等[11]研究了一維細(xì)長桿的彎曲振動特性,使用平面波展開法計算了桿的扭轉(zhuǎn)帶隙分布情況。鄭玲等[12]研究了鋁/橡膠聲子晶體對橫波的機(jī)械濾波特性,采用有限元法和傳遞矩陣法計算了能帶結(jié)構(gòu)，并且與實驗進(jìn)行了對比。Yan[13]研究了不同材料對一維聲子晶體帶隙的影響，認(rèn)為高密度的材料更容易得到較寬的彈性波帶隙。Gao等[14-16]深入研究了沿軸向排列的單槽、雙槽及多層超材料板，對其低頻帶隙機(jī)理進(jìn)行了分析。Xuan等[17]提出了一種可調(diào)節(jié)一維聲子晶體亞型結(jié)構(gòu)，實現(xiàn)彎曲波相位在0～2π移動。He等[18]研究了彈性支撐條件下一維聲子晶體的帶隙特性，并分析了不同工況對帶隙的影響。

已有研究成果多針對單一類型帶隙結(jié)構(gòu)，如布拉格散射型帶隙或局域共振型帶隙進(jìn)行研究。本文結(jié)合兩種類型帶隙機(jī)理特性提出一種管道超結(jié)構(gòu)構(gòu)型，其元胞由管道、變徑凸臺和局域共振振子構(gòu)成。使用傳遞矩陣法和有限元法計算得到了結(jié)構(gòu)軸向帶隙分布圖，建立該管道超結(jié)構(gòu)的振動實驗平臺，測得了該結(jié)構(gòu)的傳輸特性曲線，并將實驗值和計算結(jié)果進(jìn)行了對比驗證，進(jìn)一步分析了幾何尺寸、材料參數(shù)等因素對結(jié)構(gòu)帶隙特性的影響。

1 管道超結(jié)構(gòu)構(gòu)型及帶隙特性分析

1.1 管道超結(jié)構(gòu)構(gòu)型

結(jié)合布拉格散射型和局域共振型兩種帶隙機(jī)理特性，本文提出了復(fù)合帶隙管道超結(jié)構(gòu)元胞構(gòu)型，其結(jié)構(gòu)示意圖如圖1所示。該元胞由直管段、變徑凸臺和局域共振振子共同構(gòu)成。直徑較大的凸臺嵌套在直管上，凸臺上安裝局域共振振子。局域共振振子由內(nèi)部橡膠環(huán)和外部不銹鋼組成，其元胞側(cè)視平面圖及參數(shù)如圖2所示。

圖1 管道超結(jié)構(gòu)元胞構(gòu)型示意圖Fig.1 Schematic diagram of the cell configuration of the pipeline metastructure

圖2 管道超結(jié)構(gòu)元胞側(cè)視平面圖Fig.2 Side view of the pipeline metastructure

1.2 傳遞矩陣法

圖3 管道超結(jié)構(gòu)簡化示意圖Fig.3 Simplified pipeline metastructure diagram

為研究所構(gòu)建管道超結(jié)構(gòu)帶隙特性，本文分別采用傳遞矩陣法與有限元法分析了其能帶結(jié)構(gòu)。傳遞矩陣法[19]從結(jié)構(gòu)的連續(xù)狀態(tài)參數(shù)基本方程出發(fā)，推導(dǎo)得到管道的傳遞矩陣，并根據(jù)相鄰介質(zhì)位移、應(yīng)力等連續(xù)性獲得邊界條件，最終給出結(jié)構(gòu)的解析解。本文所提出的管道超結(jié)構(gòu)元胞簡化示意圖如圖3所示。A段為直管段，變徑凸臺段以等效質(zhì)量及等效剛度的形式計入B段，該部分材料記為等效材料B，局域共振振子用彈簧質(zhì)量系統(tǒng)表示。

為研究管道軸向振動特性，設(shè)A段彈性模量、密度及截面積分別為EA、ρA和S；B段彈性模量、密度分別為EB、ρB；a1和a2分別表示A段和B段的長度，元胞尺寸a=a1+a2。

管道中縱波的一維波動方程為:

(1)

式中：u(x,t)為x處的位移。其通解可寫成如下形式：

u(x,t)=U(x)eiωt

(2)

由式(1)可得，第n個元胞中A段振幅為：

(3)

同理可得，第n個元胞中B段振幅為：

(4)

局域共振振子n的振動方程為：

Vn(x)=Bneiωt

(5)

式中：Vn為振子的軸向位移，根據(jù)振子處基體位移、振子位移和彈簧回復(fù)力之間的關(guān)系，可知：

(6)

(7)

式中：ω為彈性波的角頻率；k為局域共振振子的體積模量；E和μ為振子橡膠部分的彈性模量及泊松比；m為局域共振振子的質(zhì)量。

在第n個元胞材料A和材料B連接處，由位移和應(yīng)力連續(xù)性條件可以得到：

(8)

(9)

將式(8)、(9)寫成矩陣形式：

KAψnA=HAψnA

(10)

(11)

(12)

將式(11)、(12)寫成矩陣形式：

KBψnA=HBψ(n-1)B

(13)

ψnB=Tψ(n-1)B

(14)

由于該結(jié)構(gòu)在縱向具有周期性，利用Bloch定理[20]可以得到：

ψnB=eiqaψ(n-1)B

(15)

式中：q為Bloch波矢。將式(14)帶入到式(15)中可以轉(zhuǎn)化為求解標(biāo)準(zhǔn)矩陣的特征值問題：

|T-eiqaI|=0

(16)

求解矩陣T的特征值，即可得到波矢q與頻率ω之間的色散關(guān)系。對于每一個給定的ω，波矢q由實部和虛部兩部分組成，其中實部為該結(jié)構(gòu)的色散曲線，虛部則表示彈性波在帶隙內(nèi)衰減能力。

1.3 帶隙特性分析

為研究該結(jié)構(gòu)帶隙分布規(guī)律，本文分別采用傳遞矩陣法和有限元法展開計算。結(jié)構(gòu)直管段為鋁合金，凸臺為不銹鋼，振子內(nèi)環(huán)為聚氨酯橡膠，外環(huán)為不銹鋼，振子寬度為a3，各部分材料參數(shù)和幾何參數(shù)如表1和表2所示。用COMSOL軟件建立元胞結(jié)構(gòu)的有限元模型，如圖4所示。該結(jié)構(gòu)劃分為3 160個域單元，兩端施加Floquet周期性邊界條件，同時限制其徑向位移，使其僅產(chǎn)生軸向運(yùn)動。應(yīng)用特征頻率求解器對布里淵區(qū)邊界進(jìn)行波矢k掃描計算，可得到該元胞的軸向振動能帶結(jié)構(gòu)圖。

表1 材料參數(shù)Table 1 Material parameters

表2 結(jié)構(gòu)幾何參數(shù)Table 2 Structural geometric parameters

圖4 管道超結(jié)構(gòu)有限元模型圖Fig.4 Finite element model of the pipeline metastructure

圖5分別給出了傳遞矩陣法和有限元法計算得到的結(jié)構(gòu)軸向振動帶隙分布圖，可以看出在0～2500 Hz頻率范圍內(nèi)共有兩階明顯的軸向振動帶隙,其位置分布基本一致。傳遞矩陣法所得實波矢部分為該結(jié)構(gòu)的能帶結(jié)構(gòu)圖。其中，第一階帶隙為191～371 Hz，第二階帶隙為1 181～2 141 Hz，帶寬分別為180 Hz和960 Hz；圖5(b)的虛波矢部分反映了結(jié)構(gòu)對彈性波的衰減能力，從圖中可以看出，第一階帶隙較第二階帶隙對彈性波的衰減能力更強(qiáng)。有限元法所得帶隙頻率范圍為185～381 Hz和1 157～2 080.5 Hz，對應(yīng)帶寬分別為196 Hz和923 Hz。兩種方法所得兩階帶隙位置相近，但均存在一定的差異，這是由于推導(dǎo)傳遞矩陣時，對變徑凸臺及局域共振振子進(jìn)行等效簡化從而產(chǎn)生了誤差。兩者結(jié)果表明，該管道超結(jié)構(gòu)元胞在0～1 000 Hz頻段和1 000～2 200 Hz頻段均有帶隙分布。

圖5 管道超結(jié)構(gòu)元胞能帶圖Fig.5 Band diagram of the pipeline metastructure cell

為進(jìn)一步對比分析其帶隙產(chǎn)生機(jī)理，元胞結(jié)構(gòu)可看作由I型基胞(見圖6(a))和Ⅱ型基胞(見圖6(b))構(gòu)成，圖7給出了二種基胞和元胞結(jié)構(gòu)的帶隙分布圖?？梢钥闯觯艿莱Y(jié)構(gòu)第一階帶隙位置與基胞I型的第一階帶隙位置相近。綜合圖5和圖7分析可知，該帶隙對彈性波衰減能力較強(qiáng)，且頻率較低，為局域共振型帶隙。由于振子的共振現(xiàn)象，該頻率范圍內(nèi)的彈性波被局域化，從而導(dǎo)致彈性波無法正常傳播，該帶隙位置主要隨振子內(nèi)環(huán)彈性模量和外環(huán)質(zhì)量的改變而變化。管道超結(jié)構(gòu)第二階帶隙位置與基胞I型的第二階帶隙和基胞Ⅱ型的第一階帶隙頻率范圍相近。同時圖5中該帶隙虛波矢部分呈現(xiàn)對稱特性，應(yīng)為布拉格散射型帶隙[21]。周期出現(xiàn)的變徑凸臺及局域共振振子相對于直管段均發(fā)生尺寸變化，而這種尺寸改變會對彈性波進(jìn)行反射，導(dǎo)致入射波與反射波相互疊加，致使部分頻率的波無法傳播，從而產(chǎn)生布拉格散射帶隙。

圖6 基胞示意圖Fig.6 Schematic diagram of basic cells

1.4 有限周期結(jié)構(gòu)傳輸特性

考慮到實際工程結(jié)構(gòu)大多為有限周期結(jié)構(gòu)，這將導(dǎo)致帶隙頻率范圍內(nèi)某些彈性波的能量衰減不徹底而仍然能透過有限結(jié)構(gòu)。為了更好地研究該管道超結(jié)構(gòu)的實際應(yīng)用效果，需進(jìn)一步計算其傳輸特性，常用傳輸系數(shù)TL來描述一個結(jié)構(gòu)對振動傳播的隔離能力，如下式:

TL=20×lg(u0/u1)

(17)

式中:u0為傳遞到輸出端的位移;u1為輸入位移。

應(yīng)用COMSOL軟件建立包含4個元胞結(jié)構(gòu)的有限元模型，在其一端施加不同頻率的軸向位移激勵，在另一端獲得軸向位移響應(yīng)。圖8給出了管道超結(jié)構(gòu)傳輸特性曲線并與其帶隙分布圖進(jìn)行了對比。從圖中可以看出，傳輸特性曲線共有兩處明顯的凹陷，其頻率范圍與前述計算得到的前兩階帶隙基本一致。

圖7 基胞帶隙與管道超結(jié)構(gòu)帶隙對比Fig.7 Comparison of the band gap between the basic cell and the pipeline metastructure

圖8 管道超結(jié)構(gòu)傳輸特性曲線Fig.8 Transmission curve of the pipeline metastructure

2 實 驗

為了對比驗證所提出的管道超結(jié)構(gòu)帶隙分布和傳輸特性，本文建立了該結(jié)構(gòu)的傳輸特性測試實驗平臺。

2.1 實驗方案

實驗中，管道超結(jié)構(gòu)由4個元胞構(gòu)成，變徑凸臺及振子周期排列嵌套在直管上，用彈性繩將管道超結(jié)構(gòu)懸掛在固定支架上，使用力錘在管道一端沿水平方向敲擊，誘發(fā)該管道的軸向振動。同時在管道兩端分別布設(shè)YD9820電渦流位移傳感器，應(yīng)用B&K3160-A-042數(shù)據(jù)采集儀測量輸入與輸出位移，其實驗裝置及結(jié)構(gòu)局部圖如圖9所示。

圖9 管道超結(jié)構(gòu)傳輸特性實驗裝置圖Fig.9 Experimental devices for transmission characteristics of pipeline metastructure

2.2 實驗結(jié)果分析

實驗中，分別對1～1 000 Hz和1 000～2 200 Hz兩個頻率范圍進(jìn)行振動位移響應(yīng)測試。經(jīng)過傅里葉變換，將所得到的結(jié)果從時域轉(zhuǎn)換到頻域，并與有限元仿真結(jié)果進(jìn)行對比。圖10分別為管道超結(jié)構(gòu)在兩個頻段的傳輸特性曲線，可以看出，該結(jié)構(gòu)出現(xiàn)了兩處振動明顯衰減的帶隙，其中第一階帶隙出現(xiàn)在160～220 Hz，與數(shù)值模擬結(jié)果位置接近；第二階帶隙出現(xiàn)在1 900～2 200 Hz，與數(shù)值模擬結(jié)果1 157～2 080.5 Hz存在一定差異。實驗結(jié)果中1 000～1 800 Hz范圍內(nèi)波的衰減較弱，帶隙并不明顯，其主要原因是第二階帶隙對波的衰減特性較第一階弱，其頻率較高，振幅較小，因此實驗過程中受到干擾因素較多。其次，實驗中管道結(jié)構(gòu)由彈性繩懸空吊起，其實際邊界條件與數(shù)值理論模型存在一定差異。同時周期排列的振子質(zhì)量較大，極易引起管道的彎曲變形，對結(jié)構(gòu)軸向振動產(chǎn)生影響，難以實現(xiàn)嚴(yán)格意義上的軸向振動。綜上原因?qū)е聦嶒灲Y(jié)果與數(shù)值結(jié)果在第二階帶隙范圍處出現(xiàn)一定差異。

圖10 結(jié)構(gòu)傳輸特性曲線Fig.10 Transmission curve of the structure

3 帶隙影響因素討論

通常影響布拉格散射型帶隙的主要因素為晶格結(jié)構(gòu)、尺寸等；而局域共振帶隙則主要強(qiáng)調(diào)單個散射體的設(shè)計。因此，本節(jié)將討論元胞尺寸、凸臺尺寸、橡膠及鋼環(huán)厚度等參數(shù)對帶隙分布的影響規(guī)律。

圖11分別給出了元胞尺寸a和凸臺尺寸a2對結(jié)構(gòu)帶隙的影響分布圖。從圖中可以看出，第一階帶隙以局域共振為主型帶隙，元胞尺寸變化對其影響較??；隨著凸臺長度的增加，其上振子與其他直管段的相對質(zhì)量比有所降低，帶隙寬度減小。第二階帶隙以布拉格散射為主型帶隙，凸臺長度及其與直管段的相對比值對帶隙起始頻率和帶隙寬度影響較大，隨著凸臺長度增加，凸臺內(nèi)部傳播的彈性波波長增大，帶隙頻率先降低；當(dāng)凸臺長度增加至超過直管段時，直管段長度減小將導(dǎo)致在直管段內(nèi)部傳播的彈性波長減小，帶隙起始頻率出現(xiàn)升高。因而帶隙起始頻率和寬度呈拋物線變化，當(dāng)凸臺長度增加到約80 cm時，帶隙寬度達(dá)到峰值。

圖11 元胞和凸臺尺寸對帶隙的影響Fig.11 Effect of the cell and boss size on band gaps

分別改變局域共振振子中橡膠和鋼環(huán)厚度，帶隙分布如圖12所示?？梢钥闯觯S著橡膠厚度的增加，元胞兩階帶隙位置均逐漸降低，帶隙寬度隨橡膠厚度增加而有所減小。鋼環(huán)厚度增加降低了第一階帶隙位置，對第二階帶隙無顯著影響。這是由于橡膠和鋼環(huán)的尺寸增加，其等效剛度和等效質(zhì)量也隨之增加，固有頻率降低，導(dǎo)致第一階帶隙頻率降低。

綜上所述，元胞第一階帶隙主要為局域共振型帶隙，凸臺和振子尺寸對其影響較大，元胞尺寸對其影響較小。第二階帶隙主要為布拉格散射型帶隙，該帶隙分布隨元胞和凸臺長度改變而改變。合理設(shè)計結(jié)構(gòu)各部分幾何尺寸，可滿足工程中特定頻段同時抑振的需求。

圖12 橡膠和鋼環(huán)厚度對帶隙的影響Fig.12 Effect of thickness of rubber and ring on band gaps

4 結(jié) 論

本文設(shè)計了一種同時包含局域共振型和布拉格散射型帶隙的新型復(fù)合帶隙管道超結(jié)構(gòu)構(gòu)型。分別采用傳遞矩陣法和有限元法計算了該結(jié)構(gòu)的帶隙分布及傳輸特性曲線；搭建了包含4個元胞的管道超結(jié)構(gòu)的實驗平臺，采用力錘敲擊法進(jìn)行振動測試，并與數(shù)值計算結(jié)果進(jìn)行了對比驗證，兩者取得了較好的一致性，得到如下結(jié)論：

(1)本文提出的新型管道超結(jié)構(gòu)構(gòu)型，其能帶分布結(jié)合了局域共振型帶隙和布拉格散射型帶隙特點(diǎn)，且第二階帶隙頻率范圍較寬。

(2)該元胞結(jié)構(gòu)的第一階帶隙主要為局域共振型帶隙，凸臺和振子的幾何尺寸對其影響較大，元胞尺寸對其影響較小。其帶隙寬度隨振子寬度的增加而增加，頻率隨鋼環(huán)厚度增加而降低。

(3)該元胞結(jié)構(gòu)的第二階帶隙主要為布拉格散射型帶隙，元胞長度、凸臺長度和橡膠厚度對其影響較大，隨著元胞尺寸的增加，其帶隙中心位置逐漸降低。合理設(shè)計結(jié)構(gòu)各部分幾何尺寸，可滿足工程中特定頻段同時抑振的需求。