劉成龍,許衛(wèi)鍇,呂樹辰,祁武超

(沈陽航空航天大學(xué)，遼寧省飛行器復(fù)合材料結(jié)構(gòu)分析與仿真重點實驗室，沈陽 110136)

0 引 言

聲子晶體(phononic crystal, PnCs)[1-3]是一種具有周期結(jié)構(gòu)的人工材料，由于其所具有的特殊色散特性，能對帶隙頻率范圍內(nèi)的聲波/彈性波產(chǎn)生抑制或調(diào)控作用，從而在波導(dǎo)、頻率濾波器、聲學(xué)隱身和振動控制裝置等實際裝置的設(shè)計中具有潛在的應(yīng)用前景[4-8]。聲子晶體帶隙產(chǎn)生方式有Bragg散射機(jī)理[9]和局域共振機(jī)理[10]兩種，無論哪種機(jī)理基本都是通過合理的構(gòu)建周期性微結(jié)構(gòu)以實現(xiàn)其帶隙特性。相比于自然材料，人工超材料能夠通過合理的微結(jié)構(gòu)設(shè)計以滿足不同的設(shè)計需求[11-12]，因此成為各個領(lǐng)域研究的熱點。

當(dāng)前科技的發(fā)展對材料的設(shè)計提出了更嚴(yán)苛的要求。對于帶隙材料來說，主要體現(xiàn)在帶隙的可調(diào)諧性和多功能性。一方面，當(dāng)聲子晶體的微結(jié)構(gòu)確定后，其帶隙頻率也隨之確定，無法適應(yīng)振動環(huán)境的改變。目前的研究主要集中于通過外加電、磁、熱等物理場使材料的參數(shù)或形狀發(fā)生改變以實現(xiàn)帶隙性能的調(diào)控[13]。另一方面，材料的工作環(huán)境往往存在光、電、磁、熱、力等多場耦合的情況，因此對材料的不同屬性提出了新的要求。例如，聲光子晶體(Phoxonic crystal, PxC)可以同時實現(xiàn)聲帶隙和光帶隙[14]。

在航空航天等領(lǐng)域中，結(jié)構(gòu)與材料在服役過程中往往需要應(yīng)對多種物理場的耦合作用，如微電機(jī)系統(tǒng)、高超聲速熱防護(hù)系統(tǒng)，結(jié)構(gòu)與材料既需要應(yīng)對溫度的變化，又需要在一定程度上應(yīng)對振動和噪聲的挑戰(zhàn)。因此，在設(shè)計具有特定熱膨脹功能的結(jié)構(gòu)時，考慮其是否具有帶隙特性逐漸成為一個重要的科學(xué)問題。目前對特定熱膨脹系數(shù)材料的實現(xiàn)已經(jīng)得到了廣泛的關(guān)注[15-17]?？紤]到實際工程承載需求，Lehman和Lakes等采用三角形超材料單胞(Lehman-Lakes單胞)實現(xiàn)高剛度的設(shè)計[15,17-18]。然而，Lehman-Lakes單胞的連接方式為鉸接，不易制備，Zhang等[19-20]將其改進(jìn)為固定連接(JTCLM單胞)并通過合理的材料/尺寸設(shè)計提高了結(jié)構(gòu)的剛度特性。王炳達(dá)等[21]進(jìn)一步揭示了JTCLM單胞的許用溫變設(shè)計。

迄今為止，研究人員在熱膨脹材料的設(shè)計領(lǐng)域取得了豐碩的成果，但仍缺乏對熱膨脹材料的帶隙分析研究?？梢钥闯觯徽撌荓ehman-Lakes單胞還是JTCLM單胞，其拓?fù)錁?gòu)型都與三角形蜂窩材料類似，即有可能會存在帶隙特性。另一方面，單胞參數(shù)的改變將有可能會改變帶隙的頻率范圍。因此，對考慮熱膨脹系數(shù)設(shè)計的JTCLM單胞進(jìn)行帶隙研究，并討論相關(guān)尺寸對膨脹系數(shù)的變化和帶隙的影響成為必須要考慮的關(guān)鍵問題。

本文通過有限元方法分析了JTCLM單胞的帶隙特性，并討論了單胞的可設(shè)計幾何參數(shù)對能帶的影響。研究表明，肋條的彎曲是JTCLM單胞產(chǎn)生帶隙的重要因素，通過調(diào)整單胞的幾何參數(shù)，可以實現(xiàn)點陣材料的帶隙特征調(diào)控。

1 模型與方法

1.1 JTCLM單胞拓?fù)錁?gòu)型

本文選取文獻(xiàn)[19]中的JTCLM單胞作為研究對象，如圖1所示。其中材料1是因瓦，一種低熱膨脹系數(shù)的合金材料；材料2是鋼，一種熱膨脹系數(shù)較高的材料。其參數(shù)如表1所示。其中L=50 mm，L0=30 mm；θ和r分別為彎曲角度和彎曲半徑，代表了肋條雙材料部分的彎曲程度。

表1 復(fù)合材料中組分材料的參數(shù)

1.2 能帶結(jié)構(gòu)計算

彈性波在線彈性、各向同性的無源介質(zhì)中的控制方程為[22]

ρ(r)-1{▽[μ(r)▽×u(r)]-▽[(λ(r)+2μ(r))(▽·u(r))]}=ω2u(r)

(1)

式中：ω為角頻率；r(x,y,z)為位移矢量；u(r)為位移向量；▽=(?/?x,?/?y,?/?z)為微分算子；λ(r)和μ(r)分別為介質(zhì)的拉梅常數(shù)；ρ(r)為介質(zhì)密度。這三個材料參數(shù)和周期結(jié)構(gòu)具有相同的周期性。

根據(jù)布洛赫定理，有

(2)

式中：k=(kx,ky)為第一布里淵區(qū)內(nèi)的波矢；G為正格矢任一倒格矢。由于平移周期性，一個晶胞的空間內(nèi)場分布與整個晶格結(jié)構(gòu)的周期性是相同的。在晶體學(xué)中，倒格點陣中從某一格點出發(fā),做與相鄰格點的垂直平分面,所圍成的空間就是第一布里淵區(qū)，將第一布里淵區(qū)對稱后得到的區(qū)域稱為不可約布里淵區(qū)。對于由JTCLM單胞組成的點陣超材料，其Wigner-Seitz原胞可取為圖2(a)所示六邊形，其不可約布里淵區(qū)如圖2(b)陰影部分所示。對于帶隙分析，只需要使波矢k歷遍不可約布里淵區(qū)的邊界，即可得到該聲子晶體的帶隙。

2 結(jié)果與討論

對于圖1所示的JTCLM單胞，其等效熱膨脹系數(shù)可由q，m，n，α1、α2，θ六個參數(shù)確定。一旦材料被選定后，則熱膨脹系數(shù)α1、α2，模量比n也將確定，此時單胞的帶隙特性主要受其形狀影響。因此，本文將重點討論幾何參數(shù)對JTCLM單胞帶隙的影響。由于JTCLM單胞組成的熱膨脹點陣材料只在x-y平面內(nèi)體現(xiàn)零膨脹特性，因此本文只考慮二維帶隙即x-y平面內(nèi)帶隙。

2.1 彎曲曲率的影響

JTCLM單胞的彎曲部分是圓弧。當(dāng)彎曲曲率不同時，其曲率半徑和圓心位置也將改變，不易對雙材料彎曲部分進(jìn)行直觀表述和操作。因此可沿肋條方向定義起點M，終點N。點Q是曲線過兩端點的切線的交點，將其稱為控制點。利用點Q縱坐標(biāo)的變化可以更加直觀的表述θ的變化，如圖3所示。由幾何關(guān)系可知

(3)

可以看出，y的值越大，彎曲角度θ也越大。

圖4顯示了點Q縱坐標(biāo)y=12 mm，m=1，q=0.6時的前16階帶隙圖。可以看出，在5.35～7.74 kHz、10.42～14.39 kHz和17.93～20.96 kHz出現(xiàn)了3條較寬的帶隙。為了進(jìn)一步描述不同的彎曲曲率對帶隙的影響，通過計算不同點縱坐標(biāo)的帶隙進(jìn)行規(guī)律對比。對比圖如圖5所示。

由圖5可以看出，y=0 mm時，與常規(guī)的三角點陣結(jié)構(gòu)一樣，并未出現(xiàn)完全帶隙[23]。此時肋條部分可以等效為一直桿。當(dāng)肋條開始變彎后，意味著結(jié)構(gòu)的對稱性降低，更易于產(chǎn)生帶隙[24]，因此在1.8 kHz、3.5 kHz、7.2 kHz、13.6 kHz、21 kHz附近開始出現(xiàn)帶隙，但帶寬較窄。隨著點Q縱坐標(biāo)y的增加，6～10 kHz頻段內(nèi)的帶隙寬度先變寬，在5 mm處達(dá)到最大寬度，隨后帶隙寬度變窄，帶隙頻率范圍整體有下降的趨勢；當(dāng)y增大到7 mm后，10～14 kHz頻段和21 kHz頻率附近的帶隙寬度明顯變寬，且?guī)额l率范圍整體有下降的趨勢，尤其是帶隙的下限下降明顯。3.5 kHz以下的帶隙頻率范圍沒有明顯變化。這是由于q保持不變時，彎曲程度的改變相當(dāng)于振子(彎曲部分)等效質(zhì)量的改變，y值越大，振子的等效質(zhì)量就越大，從而出現(xiàn)禁帶頻域整體下移的趨勢。

2.2 厚度比的影響

圖6顯示了y=12 mm且q=0.6時帶隙結(jié)構(gòu)隨厚度比的變化規(guī)律。由于m表示的是兩種材料的厚度比，為了更好的描述材料在整體厚度的占比，用無量綱參數(shù)m/(m+1)表示厚度比值的變化。這里分別計算了厚度比從0.2到0.9的帶隙情況?？梢钥闯觯S著厚度比的增大，帶隙頻段出現(xiàn)了整體上移的趨勢。這是由于m/(m+1)代表的是材料1的厚度占比，其值越大表示結(jié)構(gòu)的剛度越大而振子(材料2)質(zhì)量越小，因此出現(xiàn)禁帶頻域單邊上升的趨勢。

2.3 長度比的影響

圖7顯示了固定y=12 mm且m=1時帶隙結(jié)構(gòu)隨q值的變化情況。嚴(yán)格地講，q值的變化所引起的單胞構(gòu)型差異化相較前兩節(jié)的曲率或厚度要大得多。當(dāng)L0改變時，雙層材料部分的彎曲程度也會受到影響，因此，更復(fù)雜的單胞構(gòu)型差別使得它們的帶隙特征沒有出現(xiàn)明顯規(guī)律，如圖7所示。根據(jù)Zhang等[19-20]的研究，考慮結(jié)構(gòu)的剛度設(shè)計時q不應(yīng)過小，因此在進(jìn)行材料設(shè)計時q的取值范圍為0.6～0.9。

3 結(jié) 論

JTCLM單胞點陣復(fù)合材料在滿足特定熱膨脹系數(shù)設(shè)計的同時，也具有聲子晶體的帶隙特性。本文對該單胞的帶隙特性進(jìn)行了研究，討論了不同幾何參數(shù)對帶隙特征的影響，并得到以下結(jié)論。

(1) JTCLM單胞點陣復(fù)合材料屬于彎曲主導(dǎo)型的熱膨脹材料，當(dāng)肋條存在彎曲曲率時會產(chǎn)生帶隙，且?guī)短匦允苄螤钣绊戄^大。

(2)針對不同曲率、不同厚度比和不同長度比的單胞帶隙分析，給出帶隙隨幾何參數(shù)改變的分布規(guī)律圖，為實現(xiàn)熱膨脹/帶隙雙目標(biāo)的超材料設(shè)計提供了理論依據(jù)。

(3)通過合理的選材和形狀設(shè)計，有望實現(xiàn)特定膨脹性質(zhì)和帶隙設(shè)計的雙目標(biāo)共贏，使材料具有更好的可調(diào)諧性和多功能性。