公路隧道聲環境有限元數值分析

摘要：

公路隧道聲環境復雜，影響行車的舒適性和安全性，在緊急情況下嚴重影響疏散信息的傳遞。采用有限單元法對公路隧道的聲學模態及衰變特性進行數值分析。結果表明，公路隧道的聲學模態在橫向具有對稱性，對稱的區域聲壓幅值相等，相位相同或者相反，模態頻率越高，聲壓幅值最大值越大，隧道長度對模態頻率影響較小。與擴散聲場的線性衰變特性不同，公路隧道內各點的聲壓級隨時間的衰變是非線性的，在聲源停止激勵1.5 s內，衰變速率較快，從停止激勵到聲壓級衰變30 dB的時間需要5.5～7.0 s的時間。

關鍵詞：

公路隧道； 聲場； 模態； 衰變； 數值分析

中圖分類號：

U453.8；TU112.2

文獻標志碼：A

文章編號：1674-4764（2014）05-0061-05

Numerical Analysis of Acoustic Environment with Finite Element Method

Ma Fei， Zhou Xiaoping， Li Zuwei

（College of Civil Engineering， Chongqing University， Chongqing 400045， P. R. China

Abstract：

The acoustic environment of highway tunnel is complex and has influences on the driving safety and comfort. In case of emergency， evacuation information is difficult to transfer. Finite element method is used for calculating highway tunnel acoustic modal and attenuation characteristics. The results show that acoustic mode of highway tunnel has symmetry in the horizontal pressure amplitude of symmetry regional is equal and phase is the same or opposite. The higher modal frequency， the greater sound pressure amplitude is， and the effect of tunnel length on modal frequency is smaller. Decay curves of pressure level at each point in the highway tunnel over time are nonlinear. After the source excitation， in 1.5 s， attenuation rate is faster than that of few seconds later and the attenuation time of sound pressure level 30 dB is 5.5 to 7 s.

Key words：

highway tunnel； acoustic； modal； attenuation； numerical analysis

公路隧道內噪聲源主要有交通噪聲和設備噪聲，其中交通噪聲是主要來源 [1-4]，參照《聲環境質量標準》（GB 3096-2008）的環境評價體系，隧道內噪聲水平嚴重超標。與洞外相比，行駛在隧道高分貝的噪音環境中會導致人產生煩躁、恐懼的心理，進而導致注意力不集中，嚴重影響行車的安全性和舒適性。隧道內一旦發生火災、爆炸、擁堵等應急事件，公共廣播系統發出的疏導信息或者現場指揮信息難以傳達給被困人員，易使場面失控，造成嚴重的人身和財產損失^[5]。研究公路隧道聲場特性對提高行車安全性和舒適性有著長遠意義。

公路隧道的斷面尺寸一般在10 m左右，縱向尺寸一般在100 m以上，有的隧道縱向長度長達18 km，縱向尺寸一般大于橫向尺寸30倍。根據長空間的定義，縱向尺寸大于橫向尺寸6倍以上的建筑物為長空間，因此，公路隧道屬于典型的長空間建筑物。對于一般長空間的聲學特性，采用的研究方法有統計聲學理論、幾何聲線法、有限元法、邊界元法以及經驗公式等[6-15]。已有學者基于幾何聲線法對矩形斷面的長空間進行了深入研究，研究的對象多為走廊或者是地下地鐵站等截面形式簡單，縱向尺寸與斷面尺寸的比例較小的長空間，基于有限單元法的研究則較少。有限單元法可以分析任意幾何斷面的長空間的聲場特性，還可以研究不同吸聲材料對聲場的影響。公路隧道雖與建筑物走廊、地下車站在定義上同屬于長空間，但其長度尺寸與橫斷面尺寸的比值較一般長空間大，且兩端開口。公路隧道的斷面一般為三心圓，具有較強的聲聚焦效應。在目前涉及公路隧道聲環境的研究可以歸納為以下3個方面：一是對公路隧道噪聲水平的調查研究；二是路面與車輛輪胎之間產生噪聲的機理研究；三是通過在隧道內部敷設降噪材料改善隧道聲學環境的研究[16-20]，目前尚缺少對公路隧道聲學基本特性精細化深入研究。

筆者采用有限單元法對公路隧道的聲學特性進行分析，研究公路隧道的基本聲學模態以及聲壓級隨時間的衰變規律。

1 基本假定和聲波方程

聲波是物質波，是媒質中傳播的聲壓p、質點速度v、密度的一種或多種變化，聲音也是這些物理量的變化引起人的聲覺。聲傳播過程必須滿足3個基本的物理定律即平衡方程、質量守恒定律以及描述壓強、溫度與體積等狀態參數關系的物態方程。為了使問題簡化，有如下假設：

1）空氣為理想流體，不存在粘滯性，聲波在傳遞過程中沒有能量損耗。

2）沒有聲擾動時，空氣在宏觀上是靜止的，同時空氣是均勻的，靜態壓強和靜態密度都是常數。

3）聲波傳播時，空氣中稠密和稀疏的過程是絕熱的，沒有能量交換。

聲波的三維波動方程為

2 數值分析

2.1 分析模型

公路隧道分析模型采用的橫斷面尺寸如圖1（a）所示。首先在有限元程序ANSYS建立網格，然后將網格導入聲學程序LMS Virtual.Lab完成聲學的分析。為了盡量降低洞口對洞內聲場的影響，在模態分析中縱向長度取400 m，在時程分析中縱向長度取600 m。采用六面體三維單元，單元最大邊長為0.4 m，在分析中最大分析頻率不大于100 Hz，因此能夠保證一個波長范圍內有6個單元，有限元模型如圖1（b）所示。

公路隧道內混凝土襯砌和瀝青路面都有一定的粗糙度，具有一定的吸聲性能，因此在相應位置施加吸聲邊界條件，吸聲系數取0.04 [22]。隧道兩端洞口的吸聲系數為1，即聲波傳播到此處時全部被外界吸收，沒有聲波反射至洞內。

2.2 聲學模態分析

聲學模態是聲波在空間內傳播時，在不同位置產生不同的聲壓分布，其以空間中各點的聲壓分布為特征。隧道是由襯砌和路面圍成的封閉空間，與其他建筑物一樣具有模態形狀和模態頻率。圖2為通過有限元模態分析得到的幾個代表性的隧道聲學模態。

當頻率小于或等于17.4 Hz時，聲學模態均為垂直于縱軸的平面聲波，橫斷面內的聲壓處處相等，這里定義為橫向模態，沿縱向等間隔出現零聲壓節面，模態頻率越高，縱向出現節面的數目越多。節面兩側相位相反，聲壓幅值逐漸增大。模態頻率17.5 Hz對應的聲學模態以斷面豎直中心線為對稱軸在橫斷面內聲壓左右對稱分布，相位相反，檢修道上方的區域聲壓幅值較大，而隧道豎向中心線附近區域聲壓值較小，這里定義為縱向模態。頻率為25.9 Hz和34.2 Hz所對應的聲學模態也均為縱向聲學模態，可以看出檢修道、隧道路面中心上方區域和拱頂在縱向模態中的聲壓較大。

隨著模態頻率的升高，聲壓同時沿縱向和橫向變化，顯示出橫向聲學模態和縱向聲學模態的疊加特征，節面不規則，在縱向交替出現相位相反的聲壓區域。模態頻率越高，聲壓幅值最大值越大，同時聲學模態越密集，局部模態增多，縱橫模態疊加的特征越明顯，如圖2（e）、（f）所示。

由于隧道斷面在橫向左右對稱，因此各階聲學模態在橫向表現為左右對稱，相位相同或者相反，說明斷面形狀對聲學模態在橫向的聲壓分布有重要影響。

將兩車道隧道有限元模型縱向長度設置為300和200 m，然后計算其聲學模態，并與長度為400 m的隧道模型計算結果比較，三者首次出現縱向模態的頻率均為17.5 Hz，隧道模型長度不同，首次出現相同縱向聲學模態的頻率相同，說明隧道模型長度對模態的影響較小。將公路隧道簡化為圓管，通過管道聲學計算其模態頻率，首次出現縱向模態的頻率為18.3 Hz，與有限元方法計算的數值較接近，說明有限元方法計算隧道聲學模態的可靠性。

2.3 衰變特性分析

采用有限元時程分析法對隧道內的聲壓級隨時間的變化進行研究，即聲壓級為時間的函數。聲源設置在距入口100 m處的橫斷面內，節點坐標為（0，0，100）， 激勵頻率為100 Hz，振幅為5 N/m2。以聲源為起點所在斷面開始，每隔40 m設置一個觀察斷面。在每個斷面內布置6個觀察點，坐標分別為（0，0），（1，0），（2，0）（3，0），（0，1），（0，2）。在隧道內聲場達到穩態后（約20 s）停止聲源激勵，繼續計算過程直至計算結束。然后通過后處理模塊觀察聲壓級隨時間的變化。

公路隧道內聲壓級隨時間的衰變曲線是非線性的，如圖3所示，在聲源停止激勵約1.5 s內所有觀測點聲壓級衰變的速度較快，停止激勵1.5 s后各個觀察點的聲壓級衰變速率逐漸變緩，并且衰變速率趨向一致。

除聲源所在的橫斷面的各個觀測點的衰變時間差異較大外，其它各個橫斷面內的觀測點的衰變時間基本一致。與聲源在同一橫斷面的點的衰變30 dB的時間約為5.5 s，在距聲源200 m橫斷面內所需的時間約6 s，處于兩者之間的橫斷面內的所需時間在5.5～7.0 s之間，與擴散聲場衰變速率處處相等有所差別，在公路隧道內即表現為距聲源越近，聲壓級衰變的速率越快，距聲源越遠，聲壓級衰變的較慢的特點。

從幾何聲線法的原理可以看出，聲源與接收點之間的距離為d，因壁面入射角度θ的不同而存在著數條傳播路徑，如果圖4所示，每一條聲線傳播路徑長度為d/sinθ。假設存在著入射角θ1和θ2的兩條聲線，傳播路徑長度分別為d/sin θ1和d/sin θ2。在聲速相同的情況下，從聲源同時出發，聲線通過路徑1和路徑2到達接收點的時間是不同的，到達時間相差d（1/sinθ1-1/sinθ2）/c0。入射角度不變，聲源與接收點的距離d越大，時間差越長，因此導致衰變時間沿縱向逐漸增大。

根據計算結果可知，公路隧道內的混響時間t60在12 s以上，一般住房的混響時間為0.5 s、大型廳堂的混響時間為1.2 s^[22]。在播放預警信息時，聲波長時間在隧道內回蕩，很難識別疏導信息的內容，應引起足夠的重視。要改善公路隧道內的混響，一是要控制噪聲源的輻射能量，二是在襯砌上敷設降噪材料來加強噪聲能量的吸收。

3 結論

1）公路隧道聲學模態在橫向具有對稱性，高階次的聲學模態表現為橫向聲學模態和縱向聲學模態的疊加。模態頻率越高，聲壓幅值的最大值越大。

2）公路隧道內聲壓級隨時間的衰變曲線是非線性的，大致可以分為兩個階段。擴散聲場各點的混響時間處處相等，而隧道內衰變時間與距聲源的距離有一定的關系，距聲源越近，混響時間越短，混響時間越長，聲源停止激勵后，聲壓級衰變30 dB的時間在5.5～7 s之間。

3）在公路隧道內的聲能衰變時間過長，對聽覺環境造成惡劣干擾，須通過其他措施進行改善。參考文獻：

[1]

王明年，翁漢民，關寶樹.我國3條高速公路隧道污染狀況調查[J].環境科學，1997，18（3），50-53.

Wang M N， Weng H M， Guan B X. Three highway tunnel pollution survey [J]. Environmental Science， 1997， 18（3）： 50-53.

[2]張銳，黃曉明，趙永利，等. 隧道噪聲的調查與分析[J]. 公路交通科技，2006， 23（10）：29-32.

Zhang R， Huang X M， Zhao Y L， et al. Investigation and analysis of tunnel noise [J]. Journal of Highway and Transportation Research and Development， 2006， 23（10）： 29-32.

[3]楊超，王志偉. 公路隧道通風技術現狀及發展趨勢[J].地下空間與工程學報，2011，7（4）：819-824.

Yang C， Wang Z W. The technology status and development trend of highway tunnel ventilation [J]. Chinese Journal of Underground Space and Engineering， 2011， 7（4）： 819-824.

[4]Ma F， Zhang Y X. Experiments and analysis of the acoustics properties of concrete tunnel [J]. Electronic Journal of Geotechnical Engineering， 2013，18： 2031-2040.

[5]趙煒華， 劉浩學， 丁光明. 高速公路特長隧道對駕駛人的影響[J]. 公路，2011，5（1）：108-111.

Zhao W H， Liu H X， Ding G M. Impact of expressway long tunnel on drivers [J]， Highway， 2011， 5（1）： 108-111.

[6]曾向陽.室內聲場計算機模擬發展40年（1968—2008）[J].電聲技術，2008，12（6）：13-23.

Zeng X Y. Forty years of development of room acoustic computer simulation （1968-2008） [J]. Audio Engineering， 2008， 12（6）： 13-23.

[7]Kang J. The unsuitability of the classic room acoustical theory in long enclosures [J]. Architectural Science Review， 1996， 39：89-94.

[8]Allen J B， Berkley D A. Image method for efficiently simulating small-room acoustics [J]. Journal of the Acoustical Society of America， 1979， 65（4）： 943-950.

[9]Borish J. Extension of the image model to arbitrary polyhedral [J]. Journal of the Acoustical Society of America，1984，75：1827-1836.

[10]Kang J. Sound attenuation in long enclosures[J]. Building and Environment， 1996， 31（3）： 245-253.

[11]Kang J. Reverberation in rectangular long enclosures with geometrically reflecting boundaries [J]. Acta Acustica United with Acustica， 1996， 82（3）： 509-516.

[12]Vinokur R Y. Propagation of sound generation by a point source in certain built-up urban environments [J]. Soviet Physical Acoustics， 1979， 25： 377-379.

[13]Redmore T L. A method to predict the transmission of sound through corridors [J]. Applied Acoustics， 1982， 15（2）： 133-146.

[14]Pu H J， Min H Q， Qiu X J， et al. On the sound field of a vault with two open ends [J]. Applied Acoustics， 2010， 71（6）： 556-563.

[15]羅威力，蔡銘，劉濟科.矩形截面隧道交通噪聲計算機仿真[J].振動與沖擊，2012，31（8）：33-37.

Luo W L， Cai M， Liu J K. Computer simulation of traffic noise in rectangular cross-section tunnel [J]. Journal of Vibration and Shock， 2012， 31（8）：33-37.

[16]黃學輝，尚福亮，薛紅亮，等.公路隧道降噪用吸聲材料的研制[J].武漢理工大學學報，2003，25（4）：27-30.

Huang X H， Shang F L， Xue H L， et al. Development of the sound absorption material used for noise reduction in highway tunnel [J]. Journal of Wuhan University of Technology， 2003， 25（4）： 27-30.

[17]師利明，羅德春，鄧順熙.公路隧道內噪聲預測和降噪措施的理論研究[J].中國公路學報，1999，12（5）：101-105.

Shi L M， Luo D C， Deng S X. Theoretical research on noise predication and measures of reducing noise in highway tunnel [J]. China Journal of Highway and Transport， 1999， 12（5）： 101-105.

[18]張金喜，孔靜靜，黃頌昌，等. 不同類型瀝青路面吸聲性能及機理的實驗研究[J]. 北京工業大學學報，2010，36（8）：1084-1089.

Zhang J X， Kong J J， Huang S C， et al. Experimental study on sound absorption of asphalt pavement with different types [J]. Journal of Beijing University of Technology， 2010， 36（8）：1084-1089.

[19]叢卓紅，王莎，鄭南翔. 抗滑降噪型隧道瀝青混合料設計[J]. 中國公路學報，2012，25（4）：37-42.

Cong Z H， Wang S， Zheng N X. Asphalt mixture design with anti-skidding and low noise performance in tunnel pavement [J]. China Journal of Highway and Transport， 2012， 25（4）： 37-42.

[20]詹成根，郝增恒. 隧道路面明色降噪微表處技術研究[J]. 公路工程. 2013， 38（5）：67-71.

Zhan C G， Hao Z H. The tunnel pavement light color noise reduction micro-surfacing technology research [J]， Highway Engineering， 2013， 38（5）：67-71.

[21]商躍進.有限元原理與ANSYS應用指南[M].北京：清華大學出版社，2005.

[22]馬大猷.現代聲學理論基礎[M].北京：科學出版社，2004.