韓雷 季宏麗 裘進浩

摘要：針對具有深度亞波長厚度的聲學超構(gòu)吸聲板，研究了其聲學特性和物理吸聲機制，以及進一步提高其低頻聲學特性的可能性。該吸聲板單元通過將帶有小孔的分隔板內(nèi)置于以阿基米德螺旋形成的共面腔體中形成雙層穿孔吸聲體結(jié)構(gòu)。研究表明，通過調(diào)整分隔板在腔體中的位置，可以在較寬的頻帶范圍內(nèi)改變吸聲板的共振頻率，而不改變其近似完美吸聲的聲學性能。對該吸聲板建立了完整的用于描述其聲學特性的理論模型，并通過數(shù)值模擬的方法揭示了其物理吸聲機制。同時，研究并分析了吸聲板結(jié)構(gòu)參數(shù)對其聲學特性的影響，為進一步改善其低頻聲學性能提供了基礎(chǔ)。為了驗證研究結(jié)論的有效性，采用3D打印技術(shù)對所設(shè)計的樣品進行打印，測量結(jié)果與理論計算和數(shù)值模擬得到的結(jié)果具有很好的一致性。與未加入分隔板的吸聲結(jié)構(gòu)相比，該吸聲板可以在更低的頻率下獲得近全吸聲的效果。

關(guān)鍵詞：可調(diào)近全吸聲;聲學超構(gòu)吸聲板;亞波長厚度;分隔板

中圖分類號：0422.4;TB535

文獻標志碼：A

文章編號：10044523（ 2022）01-022008

DOI： 10.1638 5/j .cnki.issn.10044523.2022.01.024

引 言

傳統(tǒng)的吸聲材料如多孔材料和纖維材料等，通常對于高頻噪聲具有良好的吸聲性能，而對于低頻噪聲吸聲效率較低[1-5]，需要增加吸聲材料的厚度才能達到更好的降噪效果6-8]。同樣，對于基于聲共振原理的傳統(tǒng)吸聲裝置，如Helmholtz吸聲器[9-13]、微穿孔板吸聲器[14-20]以及兩者組合而成的復合吸聲器[21-23]，需要使其背腔厚度達到與工作波長相當?shù)某叽缌考壊拍軐崿F(xiàn)高效的吸聲。因此，在空間及重量受到嚴格限制的飛行器結(jié)構(gòu)中，無法應用上述手段進行有效的降噪。聲學超材料的出現(xiàn)為新型噪聲控制裝置的設(shè)計帶來了新的前景，也使利用輕薄、緊湊的結(jié)構(gòu)實現(xiàn)對飛行器噪聲的控制成為可能。

近年來，基于馬大猷提出的微穿孔板聲學理論[14.24-25]，一些學者提出了采用卷曲空間構(gòu)成具有近全吸聲性能的深度亞波長吸聲板[26-28]。Cai等[27]和Li等[28]介紹了由共面螺旋管和穿孔板組成的超薄吸聲板，將傳統(tǒng)的穿孔吸聲體的垂直背腔采用共面螺旋管替代，極大地縮減了吸聲體的總體厚度，并且在特定的頻率下可以實現(xiàn)完美吸聲。Chen等[29]設(shè)計了一種由兩段具有不同截面積的同軸螺旋管串聯(lián)形成的超薄吸聲板.可以實現(xiàn)兩個吸聲峰的疊加，從而改善吸聲帶寬。一些可調(diào)[30]和多階聲學吸聲結(jié)構(gòu)[31-32]在此基礎(chǔ)上也陸續(xù)出現(xiàn)。因此，聲學超材料吸聲結(jié)構(gòu)對于飛行器噪聲控制具有較高的應用價值。

與文獻[29]中利用類似結(jié)構(gòu)實現(xiàn)多階吸聲峰從而構(gòu)建連續(xù)的寬頻帶吸收結(jié)構(gòu)不同，本文主要對單個結(jié)構(gòu)在引入分隔板前后的聲學特性變化進行研究，旨在保持結(jié)構(gòu)參數(shù)不變的情況下僅通過改變分隔板在腔體中的位置，實現(xiàn)在較寬的頻帶范圍內(nèi)調(diào)節(jié)吸聲峰頻率而不改變其吸聲性能。相對于未加入分隔板的吸聲結(jié)構(gòu)，可實現(xiàn)該吸聲板在更低頻率下的完美吸聲。

1 結(jié)構(gòu)設(shè)計及理論分析

本文詳盡分析了一種由穿孔前面板和內(nèi)部帶有穿孔分隔板的共面螺旋腔體組成的超薄吸聲板的聲學特性。其中，螺旋腔體采用阿基米德曲線形成，穿孔分隔板垂直放置在螺旋腔體中并將其分隔為兩部分。吸聲板由周期性排列的螺旋腔體和穿孔前面板單元形成（周期為l1），周期單元的內(nèi)部結(jié)構(gòu)如圖1所示。用于形成螺旋腔體的阿基米德曲線的初始和末端半徑分別為1 mm和22 mm。螺旋腔體的總有效長度leff約為95.3 mm，分隔板置于從阿基米德螺旋起點沿著螺旋方向旋轉(zhuǎn)θ=90°的位置。前面板及分隔板的小孔直徑均為di=3 mm（i=1，2），厚度為ti=1 mm。由前面板和框架組成的吸聲板的總厚度H=11 mm，其中分隔板兩側(cè)螺旋腔體的寬度wi和深度hi分別為8.5 mm和8 mm。腔體壁厚b=2 mm。

螺旋腔體在理論分析中可以簡化為由兩段同軸串聯(lián)的直通管組成[29-31]。根據(jù)電一聲類比法，吸聲單元的簡化模型及對應的等效電路如圖2所示。其中，Ri和yi分別為穿孔板等效的聲阻和聲抗，ZDi為分隔板兩側(cè)腔體等效的聲順。z0= pc/Si為聲波入射至截面積為S1的前面板一側(cè)的空氣特性阻抗，p和c分別為空氣的密度和聲速。P為人射聲波的聲壓幅值。根據(jù)等效電路圖，可以得到吸聲單元的總的聲阻抗可以表示為：

可見在加入穿孔分隔板后，形成了兩個共振結(jié)構(gòu)，即以穿孔前面板與腔1以及穿孔分隔板與腔2組成的共振結(jié)構(gòu)。假設(shè)前者共振頻率較高。當頻率升高時，余切項cot（k1l1）的值逐漸變小，式（7）的第三項逐漸可以忽略不計，該結(jié)構(gòu)的聲學特性和穿孔前面板與腔1組成的共振結(jié)構(gòu)相近，高頻阻抗特性及其吸收特性受到的影響不大。而在頻率較低時，根據(jù)馬大猷雙層微穿孔板結(jié)構(gòu)中的結(jié)論[14]，可知該結(jié)構(gòu)的共振頻率等于或低于由下式?jīng)Q定的值：

比起無穿孔分隔板的結(jié)構(gòu)，其共振頻率更低。因此通過引入穿孔分隔板改善結(jié)構(gòu)的低頻聲學性能具有一定的可行性。

2 聲學特性及物理吸聲機制

為了展示所設(shè)計結(jié)構(gòu)的聲學特性，并揭示該結(jié)構(gòu)的物理吸聲機制，采用COMSOL有限元仿真軟件對吸聲板進行數(shù)值仿真與分析。通過數(shù)值仿真以及理論模型計算得到的圖1所示吸聲板（圖3中左側(cè)插圖）的吸聲系數(shù)曲線如圖3所示。其中，理論計算得到的吸聲峰幅值為1.0，峰值頻率為648 Hz，數(shù)值模擬獲得的吸聲峰幅值為0.99，峰值頻率同為648Hz。該吸聲板的總體厚度為H=11 mm，僅為共振頻率對應波長的1/48。此外，為了與未加入分隔板的吸聲結(jié)構(gòu)進行對比，在結(jié)構(gòu)參數(shù)完全相同的情況下，對腔體中未加入穿孔分隔板的吸聲板（圖3中右側(cè)插圖）的吸聲系數(shù)曲線也進行了計算。無分隔板的吸聲結(jié)構(gòu)其共振頻率為722 Hz，明顯高于本文所設(shè)計的吸聲板的共振頻率。由此可見，在結(jié)構(gòu)參數(shù)不變的情況下，穿孔分隔板的引入能夠在不影響其吸聲性能的情況下顯著降低其工作頻率，從而改善吸聲板的低頻聲學性能。gzslib202204041831

為了研究吸聲板的物理吸聲機制，在共振頻率648 Hz下通過數(shù)值仿真得到的吸聲板內(nèi)部空氣域質(zhì)點速度以及聲壓場的分布如圖4所示。在單位幅值的聲波垂直入射的情況下，分隔板兩側(cè)的腔體存在較大的聲壓差，對腔內(nèi)的氣體形成一定的推拉效應。吸聲板前面板以及分隔板附近的空氣質(zhì)點速度的方向均向小孔處匯集，且小孔處的平均質(zhì)點速度比其他地方的質(zhì)點速度高出數(shù)十倍。因而小孔中快速移動的空氣所產(chǎn)生的較大的摩擦力極大地耗散了入射聲能，形成完美吸聲。

為了進一步理解在共振頻率下吸聲板的物理吸聲機制，吸聲板表面歸一化聲阻抗Zi/z0的實部和虛部曲線如圖5所示。為了獲得完美吸聲效果，吸聲器的聲學阻抗需要與空氣阻抗獲得完美匹配，即相對于空氣的歸一化阻抗的虛部Im（Zi/Z0）趨于0的同時，其實部Re（Zt/zo）必須趨于1。由圖5所示，吸聲板的聲阻抗虛部在共振頻率648 Hz處穿過零點，此時聲阻抗的實部趨于1（約為1.053），因而所設(shè)計的吸聲板的聲學阻抗?jié)M足阻抗匹配條件，能夠獲得完美吸聲性能。

3 結(jié)構(gòu)參數(shù)對聲學特性的影響

根據(jù)以上對吸聲系數(shù)的理論推導可知，吸聲板的聲學特性受到許多結(jié)構(gòu)參數(shù)的影響，比如孔的厚度和直徑、腔體的有效長度以及其截面積等。圖6所示為在頻率50～1500 Hz范圍內(nèi)，當其他結(jié)構(gòu)參數(shù)不變時，吸聲板的吸聲峰幅值a-peak以及共振頻率fo因上述其中一個結(jié)構(gòu)參數(shù)變化所受到的影響。可以發(fā)現(xiàn)，影響吸聲峰幅值的主要因素為孔的直徑，其次為腔體2的截面積。吸聲峰幅值對于其他結(jié)構(gòu)參數(shù)的變化并不敏感，所受到的影響基本可以忽略不計。此外，只有當孔的直徑變大時，吸聲板的共振頻率會向高頻方向偏移，而其他結(jié)構(gòu)參數(shù)變大均會使得共振頻率向低頻偏移。相比之下，腔體2的結(jié)構(gòu)參數(shù)對共振頻率的影響比腔體1的影響更大。前面板和分隔板上的孔對吸聲系數(shù)的影響較為相近，同時值得注意的是，腔體1和腔體2的長度變化均可以在不改變其他結(jié)構(gòu)參數(shù)以及保證吸聲板近似完美吸聲的情況下，在很大的頻率范圍內(nèi)調(diào)節(jié)吸聲板的共振頻率。

通過以上對吸聲板結(jié)構(gòu)參數(shù)的分析，腔體1和腔體2長度的變化主要影響吸聲板的共振頻率，而對其吸聲峰幅值的影響則可以基本忽略。同時，腔體2的結(jié)構(gòu)參數(shù)相比于腔體1對共振頻率產(chǎn)生的影響更大。因此，考慮在腔體總長度l1+l2=95.3 mm不變的情況下，研究分隔板在腔體中的不同位置對聲學特性產(chǎn)生的影響。圖7（a）所示為在保持吸聲板結(jié)構(gòu)參數(shù)不變的情況下，分隔板在腔體中的位置變化對吸聲峰幅值以及共振頻率所產(chǎn)生的影響，對應的吸聲系數(shù)曲線如圖7（b）所示。從圖中所示曲線可以觀察到，改變分隔板在腔體中的位置可以對共振頻率在較大的頻帶范圍內(nèi)進行調(diào)節(jié)，同時吸聲板的吸聲性能基本不變。

4 低頻聲學性能的改善

當吸聲結(jié)構(gòu)產(chǎn)生共振時，對應的聲阻抗虛部，即聲抗為0，此時吸聲板相對于空氣的歸一化聲阻抗的實部（聲阻）越接近于1，其產(chǎn)生的吸聲系數(shù)越接近完美吸聲。因此，為了探究圖7所示結(jié)果的根本原因，分析了隔板在腔體中不同位置時吸聲板的歸一化聲學阻抗。如圖8所示，當改變隔板在腔體中的位置時，吸聲板聲阻抗的虛部基本保持在0值不變，說明隔板的位置并不影響結(jié)構(gòu)產(chǎn)生共振，即吸聲系數(shù)始終為極大值。同時，虛部為0值時對應的頻率（即結(jié)構(gòu)共振頻率）如圖中點劃線所示，表明改變隔板的位置所導致的共振頻率的變化，使得其共振頻率可以在較寬的頻帶范圍內(nèi)調(diào)節(jié)，為改善結(jié)構(gòu)的低頻聲學特性提供可能。此外，當隔板在腔體中不同位置時，歸一化聲阻始終在1附近，保證了吸聲板始終具有近似完美吸聲的聲學性能。

相應地，根據(jù)以上結(jié)構(gòu)參數(shù)對吸聲板聲學特性的影響，可以在不改變吸聲性能的情況下進一步降低吸聲板的共振頻率。根據(jù)圖7和8中得到的結(jié)論，在不改變其他結(jié)構(gòu)參數(shù)的情況下將分隔板旋轉(zhuǎn)角度設(shè)置為θ=0°，此時其共振頻率下降為f0=628 Hz，并同時具有完美吸聲性能，如圖9（a）中曲線A所示。在此基礎(chǔ)上，進一步增加分隔板的厚度以及腔體1和2的深度至t2=2 mm，h1=h29mm，其吸聲系數(shù)曲線如圖9（a）中曲線B所示，此時共振頻率降至f0=572 Hz，且其吸聲峰幅值保持不變。與圖3所示吸聲系數(shù)曲線對應的吸聲板相比，此時的吸聲板總體厚度僅增加了1 mm，而共振頻率相對降低了76 Hz。對于未加入分隔板的吸聲板，共振頻率降低了150 Hz。圖9（b）所示為曲線A和B對應結(jié)構(gòu)的歸一化聲阻抗的虛部和實部曲線。二者的相對聲阻抗虛部分別在頻率628 Hz和572 Hz處穿過零點，滿足阻抗匹配的第一個條件，此時相對聲阻抗的實部分別為1.092和1.191，均接近阻抗匹配的第二個條件，能夠得到近似完美的吸聲效果。因此，穿孔分隔板的引入可以很好地改善吸聲板的低頻聲學性能，降低結(jié)構(gòu)的工作頻率。

5 實驗驗證

為了驗證上述結(jié)論的有效性，采用雙傳聲器傳遞函數(shù)法[36]，利用阻抗管系統(tǒng)對吸聲板的吸聲系數(shù)進行實驗測量。聲望公司生產(chǎn)的阻抗管直徑為100mm，為了盡量減小由于邊界條件對吸聲系數(shù)產(chǎn)生的影響，在3D打印的單個實驗樣品中包含有4個完全相同的吸聲單元。本文設(shè)計了2種不同的實驗樣品，樣品1的結(jié)構(gòu)參數(shù)與第1節(jié)中描述的結(jié)構(gòu)相同，樣品2在樣品1的基礎(chǔ)上將分隔板旋轉(zhuǎn)角度調(diào)整為θ=0°，并且將分隔板的厚度和腔體的深度增加至t2=2 mm以及h1=h2=9 mm。利用3D打印技術(shù)制作的樣品1和樣品2如圖10所示，其外部結(jié)構(gòu)均為直徑100 mm的圓盤形，其中前面板上孔的中心間距為35.53 mm。

通過實驗測量、數(shù)值仿真以及理論計算3種方法得到的兩個樣品的吸聲系數(shù)曲線分別如圖11（a）和11（b）所示。其中理論計算和仿真得到的曲線比較一致，而實驗測得的吸聲系數(shù)其峰值和共振頻率與以上二者略有偏差。造成這種情況的原因一方面是由于實驗樣品的制造誤差所產(chǎn)生的影響，尤其是分隔板尺寸的誤差;另一方面是由于實驗條件下樣品的邊界條件與理論和仿真中所設(shè)定的周期性邊界條件不匹配。盡管如此，通過上述3種方法獲得的吸聲系數(shù)仍能夠充分證明所設(shè)計吸聲板的聲學性能，并且結(jié)構(gòu)參數(shù)對吸聲系數(shù)的影響與以上分析完全一致，證實了以上研究的正確性。gzslib202204041831

6 結(jié)論

針對本文提出的內(nèi)置有穿孔分隔板的共面螺旋腔體超薄吸聲結(jié)構(gòu)，在結(jié)構(gòu)參數(shù)不變的情況下，相比未加入分隔板的結(jié)構(gòu)，其吸聲頻率更低。此外，僅通過改變分隔板在腔體中的位置便可在較寬的頻帶范圍內(nèi)調(diào)節(jié)其共振頻率。通過數(shù)值仿真的方法揭示了該吸聲板的物理吸聲機理，并且采用理論模型詳細分析了結(jié)構(gòu)參數(shù)對其聲學特性的影響。基于上述研究，帶孔分隔板的引入可以使得該結(jié)構(gòu)獲得更低的共振頻率，從而改善其低頻聲學性能，并且通過理論分析、數(shù)值仿真以及實驗測量驗證了以上分析的有效性。小巧輕薄的物理結(jié)構(gòu)使得該吸聲板能夠在空間受限的情況下得到良好的應用。此外，該結(jié)構(gòu)可以方便地利用3D打印技術(shù)進行制作，相比于傳統(tǒng)吸聲結(jié)構(gòu)具有更大的實用性。

參考文獻：

[1]Zarek J H B. Sound absorption in flexible porous materi-als [J] . Joumal of Sound and Vibration ， 1978 ， 61 （ 2） ：205234.

[2]Arenas J P. Crocker M J. Recent trends in poroussoundabsorbing materials [J]. Sound and Vibration，2010 ， 44（ 7） ： 12-17.

[3]Yu N， Genevet P， Kats M A， et al. Light propagationwith phase discontinuities ： generalized laws of reflectionand refraction[J] . Science ， 2011 ， 334（ 6054） ： 333337.

[4]Yang J， Lee J S. Kim Y Y. Metaporous layer to over-come the thickness constraint for broadband sound ab-sorption [Jl. Journal of Applied Physics， 2015， 117（17） ： 174-903.

[5]Liu S， Chen W， Zhang Y. Design optimization of po-rous fibrous material for maximizing absorption ofsounds under set frequency bands [J]. Applied Acous-tics ， 2014 . 76 ： 319328.

[6]Norton M P， Nelson F C. Fundamentals of Noise andVibration Analysis for Engineers [M]. 2nd ed. Cam-bridge ： Cambridge University Press ， 2003.

[7]ALRahman L A， Raja R I， Rahman R A. Attenuationof noise by using absorption materials and barriers ： a re-view [ J] . Intemational Journal of Engineering and Tech-nology， 2012 ， 2（ 7） ： 1207-1217.

[8]Delany M E， Bazley E N. Acoustical properties of fi-brous absorbent materials [J]. Applied Acoustics，1970， 3（2） ： 105-116.

[9]Tang S K. On Helmholtz resonators with tapered necks[Jl. Journal of Sound and Vibration， 2005， 279： 10851096.

[10]Tang S K， Ng C H， Lam E Y L. Experimental investi-gation of the sound absorption performance of compart-mented Helmholtz resonators [J]. Applied Acoustics，2012 ， 73 ： 969976.

[11]Jordan V L. The application of Helmholtz resonators tosoundabsorbing structures [Jl. The Journal of theAcoustical Society of America ， 1947 ， 19（ 6） ： 972981.

[12]Kim S ， Kim Y H ， Jang J H. A theoretical model to pre-dict the lowfrequency sound absorption of a Helmholtzresonator array （L） [J]. The Joumal of the AcousticalSociety of America ， 2006 ， 119（4） ： 1933-1936.

[13]Ingard U. On the theory and design of acoustic resona-tors [ J] . The Joumal of the Acoustical Society of Ameri-ca， 1953， 25 （ 6） ： 1037-1061.gzslib202204041831

[14]Maa D Y. Theory and design of microperforated panelsoundabsorbing constructions [J]. Scientia Sinica，1975， 18（1） ： 55-71.

[15]Sakagami K， Nagayama Y， Morimoto M， et al. Pilotstudy on wideband sound absorber obtained by combina-tion of two different microperforated panel （MPP） ab-sorbers [ J] . Acoustical Science and Technology ， 2009 ，30： 154-156.

[16]Li X， Wu Q， Kang L， et al. Design of multiple parallel-arranged perforated panel absorbers for low frequencysound absorption [J] . Materials ， 2019 ， 12（ 13） ： 2099.

[17]Wang C ， Huang L. On the acoustic properties of parallelarrangement of multiple microperforated panel absorb-ers with different cavity depths EJl. The Journal of theAcoustical Society of America ， 2011 ， 130（ 1） ： 208218.

[18]閔鶴群，郭文成，具有并聯(lián)不等深度子背腔序列的微穿孔板吸聲體吸聲特性[J].東南大學學報 ， 2017， 47（1） ： 177-183.

Min Hequn， Guo Wencheng. Absorption characteristicsof microperforated panel sound absorbers with array ofparallel-arranged sub-cavities with different depths E Jl .Journal of Southeast University ， 2017 ， 47 （ 1） ： 177-183.

[19]Zhu X Z， Chen Z B， Jiao Y H， et al. Broadening of thesound absorption bandwidth of the perforated panel us-ing a membranetype resonator [ J] . Journal of Vibrationand Acoustics ， 2018 ， 140 ： 031014.

[20]Lee Y Y. Lee E W M， Ng C F. Sound absorption of afinite flexible microperforated panel backed by an aircavity [J] . Journal of Sound and Vibration . 2005 ， 287 ：227-243.

[21]Gai X L， Xing T， Li X H， et al. Sound absorption ofmicroperforated panel mounted with Helmholtz resona-tors [ J ] . Applied Acoustics ， 2016 ， 114 ： 260265.

[22]Toyoda M， Mu R L， Takahashi D. Relationship be-tween Helmholtzresonance absorption and panel-typeabsorption in finite flexible microperforatedpanel ab-sorbers[ J] . Applied Acoustics ， 2010 ， 71： 315320.

[23]Park S H. Acoustic properties of microperforated panelabsorbers backed by Helmholtz resonators for the im-provement of lowfrequency sound absorptionEJl . Jour-nal of Sound and Vibration ， 2013 ， 332 ： 48954911.

[24]Maa D Y. General theory and design of microperforatedpanel absorbers [J]. Chinese Joumal of Acoustics，1997. 16（3） ： 193202.

[25]Maa D Y. Design of microperforated panel constructions[J]. Chinese Joumal of Acoustics， 1988， 7 （3） ：193200.

[26]Ma F Y ， Chen J Y， Wu J H. Threedimensional acous-tic subdiffraction focusing by coiled metamaterials withstrong absorption [J]. Journal of Materials ChemistryC ， 2019 ， 7 ： 5131-5138.gzslib202204041832

[27]Cai X B， Guo Q Q， Hu G K， et al. Ultrathin lowfre-quency sound absorbing panels based on coplanar spiraltubes or coplanar Helmholtz resonators [J . AppliedPhysics Letters ， 2014 ， 105（ 12） ： 339356.

[28]Li Y， Assouar B M. Acoustic metasurfacebased per-fect absorber with deep subwavelength thickness [J].Applied Physics Letters ， 2016， 108（6） ： 063502.

[29]Chen C， Du Z， Hu G， et al. A lowfrequency sound ab-sorbing material with subwavelength thickness [J] . Ap-plied Physics Letters ， 2017 ， 110 ： 221903.

[30]Wang Y， Zhao H G， Yang H B， et al. A tunable soundabsorbing metamaterial based on coiledup space EJl.Journal of Applied Physics ， 2018 ， 123 ： 185109.

[31]Liu C R， Wu J H， Chen X， et al. A thin lowfrequencybroadband metasurface with multiorder sound absorp-tion[J]. Journal of Physics D Applied Physics， 2019，52（10） ： 105302.

[32]Liu C R， Wu J H， Ma F Y， et al. A thin multiorderHelmholtz metamaterial with perfect broadband acousticabsorption [J]. Applied Physics Express. 2019， 12：084002.

[33]Maa D Y. Potential of microperforated panel absorber[ J] . The Journal of the Acoustical Society of America ，1998， 104（5） ： 2861-2866.

[34]Kirchhoff G. Ueber den Einflufs der Warmeleitung in ei-nem Gase auf die Schallbewegung [J]. Annalen derPhysik ， 1868 ， 134 ： 177-193.

[35]Stinson M R. The propagation of plane sound waves innarrow and wide circular tubes， and generalization touniform tubes of arbitrary crosssectional shape [ J] . TheJournal of the Acoustical Society of America， 1991，89 ： 550558.