任 怡 張 昊 張 望 劉海東 李 釗 張滿弓,2)

* (西南科技大學(xué)制造過程測試技術(shù)教育部重點實驗室,四川綿陽 621010)

? (武漢第二船舶設(shè)計研究所,武漢 430064)

** (天津航天瑞萊科技有限公司,天津 300462)

引言

點陣結(jié)構(gòu)一般是由連續(xù)、重復(fù)排列的相互連通的單元組成的三維結(jié)構(gòu),也可以理解為三維空間中由相互連通的支柱和節(jié)點組成的多孔材料結(jié)構(gòu)[1-3].點陣結(jié)構(gòu)的周期和多孔特征使其具有特殊的力學(xué)和聲學(xué)特性,如比強度高、比剛度大、輕質(zhì)和隔振等,在航空航天、艦船等領(lǐng)域具有十分廣泛的應(yīng)用前景[4-8].

對點陣結(jié)構(gòu)的減振和隔聲性能,已經(jīng)有廣泛的理論[9-10]和實驗研究[11-13],點陣結(jié)構(gòu)通過單胞的周期性排列獲得聲子晶體帶隙,在帶隙范圍內(nèi)抑制彈性波的傳播[14-16].An 等[17]提出了一種基于拉伸主導(dǎo)的三維聲學(xué)超材料點陣結(jié)構(gòu),通過振動傳遞試驗驗證了其結(jié)構(gòu)的減振性能.Liu 等[18]將一些輔助結(jié)構(gòu)引入點陣結(jié)構(gòu),使用布洛赫定理和有限元相結(jié)合的方法研究了彈性波在這種改進的二維點陣結(jié)構(gòu)中的傳播,通過對輔助結(jié)構(gòu)的參數(shù)分析,設(shè)計了一種基于改進點陣結(jié)構(gòu)的寬頻帶隔振結(jié)構(gòu).Matlack 等[19]將局域共振結(jié)構(gòu)和點陣結(jié)構(gòu)相結(jié)合,提出了一種彈性超結(jié)構(gòu),可以實現(xiàn)低頻寬帶隙的減振效果.Cao 等[20]通過有限元方法和實驗研究了桿件截面變化的點陣結(jié)構(gòu),討論了形狀參數(shù)對點陣結(jié)構(gòu)的力學(xué)性能和能量吸收能力的影響,并得到了最優(yōu)的桿件形狀參數(shù).文獻[21]在立方體點陣結(jié)構(gòu)上附加了局部諧振器,研究了不同尺寸和質(zhì)量的局部諧振器對帶隙的寬度和位置的影響.Wu 等[22-23]使用譜元法研究了由兩種不同材料組成的晶格結(jié)構(gòu)的振動帶隙特性,參數(shù)分析了結(jié)構(gòu)和材料對帶隙的影響.Syam 等[24]設(shè)計了六種不同的點陣結(jié)構(gòu),并使用有限元方法和實驗研究了這些結(jié)構(gòu)的減振性能.Wang 等[25]使用遺傳算法在最小重量的前提下對夾層板的聲學(xué)和力學(xué)性能進行了優(yōu)化設(shè)計研究,最終得到了一種重量輕、機械效率高的隔音材料.

與目前廣泛研究的基于二維平面和三維長方體的點陣結(jié)構(gòu)不同,點陣圓筒結(jié)構(gòu)是沿軸向和環(huán)向周期性陣列的結(jié)構(gòu).對點陣圓筒結(jié)構(gòu)的研究主要集中于結(jié)構(gòu)的力學(xué)性能[26-28],而缺乏對其低頻寬帶減振的相關(guān)研究.事實上,針對圓筒結(jié)構(gòu)(如管道)等,存在著豐富而又迫切的應(yīng)用場景,如化工管道和艦船管道的低頻減振降噪控制等,合理的骨架結(jié)構(gòu)形式將是能否實現(xiàn)預(yù)期隔振效果的關(guān)鍵.因此,本文以金字塔型單胞為對象,通過環(huán)向和軸向陣列形成點陣圓筒結(jié)構(gòu),針對獨立參數(shù)進行結(jié)構(gòu)設(shè)計,并聯(lián)合ABAQUS 對結(jié)構(gòu)軸向減振性能進行分析.利用無量綱參數(shù)系統(tǒng)分析了金字塔單胞的高和寬,桿件的直徑以及環(huán)向陣列的個數(shù)對金字塔點陣圓筒結(jié)構(gòu)軸向減振性能的影響規(guī)律.在此基礎(chǔ)上,設(shè)計了一種具有低頻寬帶減振性能的金字塔點陣圓筒結(jié)構(gòu),利用3D打印技術(shù)制備了樣件,完成振動實驗測試并與有限元模擬結(jié)果進行對比.

1 金字塔點陣圓筒結(jié)構(gòu)參數(shù)化模型

通過Rhino 和Grasshopper 等軟件建立的金字塔點陣圓筒模型如圖1 所示,主要結(jié)構(gòu)參數(shù)如下: 金字塔單胞軸向和環(huán)向陣列的個數(shù)分別為m和n,結(jié)構(gòu)的內(nèi)接圓和外接圓的半徑分別為rc和Rc,金字塔單胞的四根桿件的交點為底面四邊形的中心,為等腰四棱錐結(jié)構(gòu),桿件為圓柱桿,直徑為d,單胞的長寬高分別為a,b和h,一共有8 個結(jié)構(gòu)參數(shù),它們之間滿足如下關(guān)系

圖1 金字塔點陣圓筒結(jié)構(gòu)示意圖Fig.1 Schematic diagram of the pyramid lattice cylinder structure

所以獨立的結(jié)構(gòu)參數(shù)只有6 個.

Rhino 是目前普遍應(yīng)用的三維NURBS 工業(yè)建模軟件,Grasshopper 是與Rhino 緊密結(jié)合的參數(shù)化設(shè)計插件.Grasshopper 作為一種節(jié)點可視化編程語言,是一種圖示算法編輯器或者可視化腳本編輯器,可高效實現(xiàn)實時數(shù)據(jù)結(jié)構(gòu)調(diào)整并觀察、記錄結(jié)構(gòu)的推敲過程[29].因此,本文使用Rhino 和Grasshopper軟件對結(jié)構(gòu)進行建模,固定選取金字塔點陣圓筒結(jié)構(gòu)的內(nèi)接圓半徑rc=100 mm,通過改變桿件的直徑、金字塔單胞的寬和高、單胞環(huán)向和軸向陣列的個數(shù),實現(xiàn)對金字塔點陣圓筒結(jié)構(gòu)的設(shè)計和模型生成.

2 結(jié)構(gòu)軸向振動有限元模型

本研究關(guān)注金字塔點陣圓筒結(jié)構(gòu)軸向振動的控制情況,為了驗證結(jié)構(gòu)的減振性能,利用有限元仿真分析結(jié)構(gòu)的振動傳輸特性,通過在一定頻率范圍內(nèi)進行諧響應(yīng)分析得到響應(yīng)曲線,從曲線中可以直觀獲得結(jié)構(gòu)對彈性波的衰減能力.傳輸特性函數(shù)(FRF)的定義為

式中,u1和u2為振動傳輸激勵和響應(yīng)的位移值.

利用有限元軟件ABAQUS 進行模擬仿真[30],在金字塔點陣圓筒結(jié)構(gòu)的左端面施加z方向的位移激勵,并在右端面提取位移響應(yīng),有限元模型如圖2 所示,選用梁單元.計算中使用的具體結(jié)構(gòu)參數(shù)如表1所示,使用的材料為樹脂,密度ρ=1180 kg/m3,彈性模量E=2 GPa .圖3 給出了兩個模型(軸向陣列個數(shù)不同)在軸向位移激勵下的振動傳輸特性曲線.從圖中可以看出,在116～466 Hz 范圍內(nèi),兩個模型的軸向振動都有較大衰減,且隨著軸向陣列個數(shù)的增加,振動衰減的強度逐漸增大.

圖2 有限元模型Fig.2 Finite element analysis model

表1 金字塔點陣圓筒模型結(jié)構(gòu)參數(shù)Table 1 Structural parameters of pyramid lattice cylinder model

圖3 振動傳輸特性曲線Fig.3 Vibration transmission characteristic curve

3 金字塔點陣圓筒結(jié)構(gòu)參數(shù)影響分析

量綱分析理論表明,某一物理過程不因其中相關(guān)物理量量綱選擇的任意性而改變,由此可以得到物理量之間的關(guān)系[31].因此,本文將引入歸一化頻率的概念,用于對頻率進行無量綱處理.歸一化頻率定義為

式中,f,rc和分別為結(jié)構(gòu)的頻率、內(nèi)接圓半徑和縱波波速,根據(jù)歸一化頻率的表達式可以得知,歸一化頻率與具體材料選取和結(jié)構(gòu)的內(nèi)接圓半徑無關(guān).

同時為了評估結(jié)構(gòu)振動衰減的能力,使用無量綱參數(shù)—相對帶寬BG%的概念,其定義為

式中,fH和fL分別為帶隙的上、下邊界頻率值,相對帶寬即為帶寬與中心頻率之比[32-34].使用相對帶寬的概念可以降低帶隙效果評價對絕對頻率的依賴,便于在更寬頻率范圍內(nèi)統(tǒng)一評估結(jié)構(gòu)的減振效果,一般認為BG%越大,減振效果越好.

學(xué)科地位的確定，不是學(xué)科建設(shè)的終結(jié)，而是學(xué)科建設(shè)的開始。治安學(xué)雖然學(xué)科地位得以確認，但它還只是正在成長中的新興學(xué)科。審視治安學(xué)學(xué)科的形成與發(fā)展，不難發(fā)現(xiàn)學(xué)科成熟程度低，學(xué)科建設(shè)問題多的現(xiàn)狀。

下面將利用無量綱參數(shù)—歸一化頻率和相對帶寬,系統(tǒng)討論金字塔點陣圓筒結(jié)構(gòu)中單胞環(huán)向陣列的個數(shù)、單胞的寬和高以及桿件的直徑對結(jié)構(gòu)減振效果的影響.

3.1 單胞的寬對金字塔點陣圓筒結(jié)構(gòu)振動傳輸特性的影響

金字塔單胞的寬決定了點陣結(jié)構(gòu)的軸向長度,當(dāng)選取單胞的高為28.4 mm,桿件直徑為1 mm,軸向和環(huán)向陣列個數(shù)分別為15 和20,此時金字塔點陣圓筒結(jié)構(gòu)的外接圓半徑為130 mm.改變單胞的寬,得到不同單胞寬的金字塔點陣圓筒結(jié)構(gòu)振動傳輸特性曲線如圖4 所示,當(dāng)單胞的寬為20.00 mm,33.33 mm和46.67 mm 時,帶隙的起始歸一化頻率分別為:0.010,0.009 和0.008,終止歸一化頻率分別為: 0.047,0.040 和0.032.由此可見,隨著單胞的寬增加,帶隙向低頻移動,帶寬變小,相對帶寬分別為129.82%,126.53%和120.00%.

圖4 不同單胞寬下的振動傳輸特性曲線Fig.4 Vibration transmission characteristic curves of different unit cell widths

3.2 單胞的高對金字塔點陣圓筒結(jié)構(gòu)振動傳輸特性的影響

選取單胞的寬為40 mm,桿件直徑為1 mm,軸向和環(huán)向陣列個數(shù)分別為15 和20 時,改變單胞的高,分別取8.65 mm,18.52 mm 和28.40 mm,此時的金字塔點陣圓筒結(jié)構(gòu)的外接圓半徑分別為110 mm,120 mm 和130 mm.不同單胞高對金字塔點陣圓筒結(jié)構(gòu)振動傳輸特性的影響如圖5 所示,當(dāng)單胞的寬為8.65 mm,18.52 mm 和28.40 mm 時,帶隙的起始歸一化頻率分別為: 0.024,0.013 和0.009,終止歸一化頻率分別為: 0.062,0.045 和0.036.由此可見,隨著單胞高度增加,帶隙向低頻移動,帶寬變小,相對帶寬分別為88.37%,110.34%和120.00%.

圖5 不同單胞高下的振動傳輸特性曲線Fig.5 Vibration transmission characteristic curves of different unit cell heights

綜合考慮金字塔單胞的高和寬對整個結(jié)構(gòu)減振性能的影響,從圖6 可以看出,隨著金字塔單胞的高和寬的增加,帶隙頻率的上限和下限都在降低,帶隙頻率的上限比下限降低的幅度更大,所以結(jié)構(gòu)的帶隙往低頻移動時,整個帶隙的寬度也在變小.從相對帶寬的角度去評價整個結(jié)構(gòu)的低頻寬帶減振性能,如圖6(d)所示,隨著金字塔單胞寬的增加,相對帶寬變小,隨著單胞高度增加,相對帶寬變大.相比于單胞寬的影響,相同的變化幅值下,單胞高對結(jié)構(gòu)的減振性能影響更大,所以在考慮設(shè)計金字塔點陣圓筒結(jié)構(gòu)用于振動控制時,應(yīng)該主要考慮金字塔單胞高對結(jié)構(gòu)減振性能的影響.

圖6 單胞的高和寬對結(jié)構(gòu)減振性能的影響Fig.6 The effect of unit cell height and width on structural vibration reduction performance

圖6 單胞的高和寬對結(jié)構(gòu)減振性能的影響(續(xù))Fig.6 The effect of unit cell height and width on structural vibration reduction performance (continued)

3.3 單胞桿件的直徑對金字塔點陣圓筒結(jié)構(gòu)振動傳輸特性的影響

選取單胞的寬為40 mm,單胞的高為28.4 mm,軸向和環(huán)向陣列個數(shù)分別為15 和20 時,改變金字塔單胞桿件的直徑,分別取0.5 mm,1.0 mm 和1.5 mm,得到圖7 所示的傳遞特性曲線,當(dāng)桿件的直徑為0.5 mm,1.0 mm 和1.5 mm 時,帶隙的起始歸一化頻率分別為: 0.004,0.009 和0.013,終止歸一化頻率分別為: 0.018,0.036 和0.054.由此可見,隨著桿件的直徑增加,帶隙向高頻移動,帶寬變大,相對帶寬分別為127.27%,120.00%和122.39%.桿件的直徑直接影響整個結(jié)構(gòu)的剪切剛度和等效密度,改變整個結(jié)構(gòu)的減振性能,桿件越細,在低頻時的減振效果越好,但是結(jié)構(gòu)的剛度會隨之變小,所以在結(jié)構(gòu)設(shè)計中,應(yīng)該首先滿足結(jié)構(gòu)剛度的需要,然后選擇較小的桿件直徑,以達到低頻減振的效果.

圖7 不同桿件直徑的振動傳輸特性曲線Fig.7 Vibration transmission characteristic curves of different rod diameters

3.4 單胞環(huán)向陣列的個數(shù)對金字塔點陣圓筒結(jié)構(gòu)振動傳輸特性的影響

選取單胞的高和寬分別為28.4 mm 和40 mm,桿件直徑為1 mm,軸向陣列個數(shù)為15,此時金字塔點陣圓筒結(jié)構(gòu)的外接圓半徑為130 mm.改變單胞環(huán)向陣列的個數(shù),得到不同單胞環(huán)向陣列個數(shù)下金字塔點陣圓筒結(jié)構(gòu)振動傳輸特性曲線如圖8 所示,當(dāng)環(huán)向陣列的個數(shù)為10,20,30 和40 時,帶隙的起始歸一化頻率分別為: 0.007,0.009,0.010 和0.010,終止歸一化頻率分別為: 0.027,0.038,0.041 和0.044.由此可見,隨著環(huán)向陣列個數(shù)的增加,帶隙向高頻移動,帶寬變大,相對帶寬分別為121.89%,123.72%,120.93%和124.72%.從中發(fā)現(xiàn),當(dāng)環(huán)向陣列個數(shù)從30 增加到40 時,結(jié)構(gòu)的帶隙范圍變化顯著減小,對結(jié)構(gòu)減振性能的影響減弱.

圖8 不同單胞環(huán)向陣列個數(shù)下的圓筒振動傳輸特性曲線Fig.8 Vibration transmission characteristic curves of different numbers of unit cell circular arrays

3.5 等比例模型參數(shù)分析

本文重點關(guān)注的結(jié)構(gòu)參數(shù)為單胞的高、寬和桿件的直徑.為了進一步探究這些結(jié)構(gòu)參數(shù)之間的關(guān)系,使單胞的高、寬和桿件的直徑控制在一定的比例,如b:h:d=20.0:14.2:0.5,同時改變?nèi)齻€參數(shù)的數(shù)值,如表2 所示,此時固定金字塔單胞軸向、環(huán)向陣列的個數(shù)為15 和20,結(jié)構(gòu)內(nèi)接圓的半徑為 100 mm.得到的振動傳遞曲線和相對帶寬隨著桿件直徑變化的結(jié)果如圖9 所示,當(dāng)單胞的高、寬和桿件的直徑三個參數(shù)的比值為上述定值時,隨著桿件的直徑增加,帶隙向低頻移動,帶寬減小,但是結(jié)構(gòu)的相對帶寬BG%是一樣的,均為122.50%左右.

表2 等比例模型的結(jié)構(gòu)參數(shù)Table 2 Structural parameters of the equal scale model

圖9 等比例模型的帶隙范圍和相對帶寬Fig.9 Band gap and relative bandwidth for the equal scale model

在實際應(yīng)用中,當(dāng)結(jié)構(gòu)內(nèi)接圓的半徑為 100 mm時,為了保證結(jié)構(gòu)整體的協(xié)調(diào)性和力學(xué)性能,單胞的高、寬應(yīng)不大于60 mm,桿件的直徑不大于5 mm,此時單胞的高、寬和桿件的直徑的比值如果相同,帶隙頻率會隨著結(jié)構(gòu)參數(shù)有規(guī)律的變化.定義無量綱結(jié)構(gòu)參數(shù)k為直徑的變化量 Δ 與結(jié)構(gòu)內(nèi)接圓半徑rc的比值,其表達式為k=Δ/rc.如圖9 所示,以桿件直徑為1 mm 的模型為基準,隨著直徑大小改變,無量綱結(jié)構(gòu)參數(shù)k相應(yīng)變化,其帶隙范圍變?yōu)閇2-100k fL0,2-100k fH0],fL0和fH0分別為桿件直徑為1 mm 時的帶隙頻率下限和上限.

點陣結(jié)構(gòu)的固有頻率決定了帶隙寬度和頻率位置,固有頻率f0正比于結(jié)構(gòu)剛度k0與質(zhì)量m0比值的1/2 次方,即[24].結(jié)構(gòu)的剛度越低,參與振動的質(zhì)量越大,那么結(jié)構(gòu)的固有頻率就會越低,帶隙范圍的頻率也就越低[35].圖10 為等比例模型中隨桿件直徑變化的結(jié)果,其反映了結(jié)構(gòu)固有頻率的變化趨勢,隨著等比例模型中桿件直徑的增加,結(jié)構(gòu)固有頻率減小,帶隙向低頻移動,與圖9 中的結(jié)果相符合.

圖10 等比例模型的Fig.10 for the equal scale model

4 實驗驗證

通過本文總結(jié)的參數(shù)規(guī)律,可以根據(jù)實際的減振需要設(shè)計具體的金字塔點陣圓筒結(jié)構(gòu),首先確定一組基準模型結(jié)構(gòu)參數(shù),通過歸一化頻率比較結(jié)構(gòu)的帶隙范圍和需要的帶隙范圍的差距,改變基準模型結(jié)構(gòu)的單胞高、單胞寬和桿件的直徑之間的比例,使結(jié)構(gòu)的相對帶寬和所需要的減振結(jié)構(gòu)的相對帶寬一致,然后再調(diào)整等比例模型的具體尺寸大小,使結(jié)構(gòu)完全滿足所需要的減振效果.

在實驗中,選用增材制造技術(shù)制備樣件,希望結(jié)構(gòu)能夠在500～1500 Hz 內(nèi)能有一個良好的減振性能,以表1 中軸向周期數(shù)為8 的金字塔點陣圓筒結(jié)構(gòu)為基準模型,固定結(jié)構(gòu)內(nèi)接圓半徑rc=100 mm,通過上述的設(shè)計方法,進行參數(shù)調(diào)整,最終確定了滿足減振要求的結(jié)構(gòu)參數(shù)為:rc=100 mm,b=40 mm,h=28.4 mm,d=4 mm,m=8和n=20 .

由于金字塔點陣圓筒結(jié)構(gòu)兩端不便于安裝和測量,所以在結(jié)構(gòu)的兩端加裝一個厚度Hb=5 mm,直徑Db=296 mm 的圓板.3D 打印時使用的材料為光敏樹脂,是一種具有精確和耐久特性的類ABS 的立體光造型樹脂,材料參數(shù)為密度ρ=1180 kg/m3,彈性模量E=2 GPa .蓋板對傳遞特性的影響如圖11所示,發(fā)現(xiàn)蓋板并不影響整體振動衰減的效果,因此使用這種加蓋的結(jié)構(gòu)驗證純金字塔點陣圓筒結(jié)構(gòu)的減振效果是可行的.圖12 為實驗用金字塔圓筒點陣結(jié)構(gòu)的模型圖和實物圖.

圖11 有無蓋板的傳遞特性曲線對比Fig.11 Comparison of transfer characteristic curves with and without cover

圖12 實驗用金字塔點陣圓筒結(jié)構(gòu)Fig.12 Pyramid lattice cylinder structure in experiment

隨機振動的頻帶寬,且有連續(xù)的頻譜,能同時在所有的頻率上對試件進行激勵,相比于簡諧振動僅對某些頻率或連續(xù)掃頻來模擬實際環(huán)境振動的影響,隨機激勵信號的振動實驗測試更能反映結(jié)構(gòu)在一定頻率范圍內(nèi)的真實減振效果,因此在實驗中采用隨機激勵信號.實驗示意圖如圖13(a)所示,樣件的一端通過夾具固定在振動臺上,激勵端的加速度傳感器固定在振動臺上,用于接收響應(yīng)的兩個加速度傳感器固定在樣件的上端面,通過振動臺產(chǎn)生隨機激勵信號,響應(yīng)端的信號被接收后通過數(shù)據(jù)處理,得到實驗振動傳輸曲線如圖13(b)所示,結(jié)果表明,在500～1500 Hz,實驗結(jié)果與有限元模擬的結(jié)果基本一致,結(jié)構(gòu)體現(xiàn)了良好的減振效果,平均的衰減強度達到了50 dB 左右.

圖13 振動測試實驗Fig.13 Vibration test experiment

5 結(jié)論

金字塔點陣圓筒結(jié)構(gòu)具有良好的軸向低頻寬帶減振性能.本文通過參數(shù)化建模,使用有限元方法計算結(jié)構(gòu)的傳遞特性函數(shù),系統(tǒng)研究了單胞的寬、高以及桿件直徑等關(guān)鍵參數(shù)對減振效果的影響,并通過振動測試實驗對結(jié)構(gòu)的減振性能進行了驗證,具體結(jié)論如下.

(1)隨著金字塔單胞寬的增加,帶隙向低頻移動,帶寬和相對帶寬都減小;隨著金字塔單胞高的增加,帶隙向低頻移動,帶寬減小,相對帶寬變大;隨著桿件的直徑和環(huán)向陣列個數(shù)增加,帶隙向高頻移動,帶寬變大.

(2)當(dāng)確定金字塔單胞軸向、環(huán)向陣列的個數(shù)和結(jié)構(gòu)內(nèi)接圓的半徑時,使單胞的高和寬以及桿件的直徑控制在一定的比例,參數(shù)的數(shù)值變大,帶隙向低頻移動,帶寬減小,但結(jié)構(gòu)的相對帶寬BG%不變.同時,在一定的尺寸范圍內(nèi),當(dāng)直徑變化k時,帶隙頻率的上下限為變化前的倍.

(3)通過結(jié)構(gòu)參數(shù)分析得到的規(guī)律,設(shè)計了一種具有低頻寬帶減振性能的金字塔點陣圓筒結(jié)構(gòu),并通過增材制造技術(shù)進行制樣,最后進行振動測試實驗,實驗與有限元結(jié)果基本一致,結(jié)構(gòu)在500～1500 Hz 有明顯減振效果.