表征人體對通過空氣傳播的聲音信號的影響并不是一項簡單的任務。為了對受人體影響的聲學裝置的性能進行準確評估，我們必須直面這項挑戰。原位電聲測試是一種標準的評估方法，然而軀干、頭部、耳廓和耳道形狀不規則，并且因人而異，常常給聲學工程師帶來極大的困擾。

為了確保采集的聲學數據盡可能準確，聲音與振動測量設備制造商 Brüel  & Kj?r開發出了頭部和軀干模擬器（Head and Torso Simulator，簡稱 HATS）。這一模型可以準確地再現普通成人頭部和軀干的幾何結構與聲學性能（圖 1）。此裝置配有耳口模擬器，是對耳機、免提通話裝置、助聽器、頭盔和聽力保護器進行電聲測試的理想選擇。耳模擬器通常是一個與前置放大器相連的半英寸耳機，而口模擬器是一臺高順性揚聲器，它通過聲壓分布來模擬普通成人講話時的情況。

圖 1. 使用 HATS（頭和軀干模擬器）測試集成式頭盔系統。

多年來，丹麥瑞聲達（GN Resound）公司一直使用 HATS 作為測量裝置，用于研究耳朵幾何結構對耳膜感知聲音的影響，并測定由耳朵和頭部幾何結構而產生的聲散射的方向性特征。近期，他們使用 COMSOL Multiphysics? 軟件構建 HATS 的計算模型。高級工程師 René Christensen 利用軟件中的邊界元法（boundary element method，簡稱 BEM）功能完成了這項任務。模型完成后，許多同事都發現各自的研究均可以采用此仿真工具替代人體模型。

為何選擇邊界元法？

軀干、頭部、耳廓和耳道對空氣傳播聲音信號以及頭部周圍聲場的綜合影響，被稱之為頭部相關傳輸函數。 “在遠離頭部的空間中指定一點，你可以描述該點發出的聲音到耳膜的傳輸過程。”Christensen 解釋道。這種確定性模擬方法可以對麥克風擺放在不同位置時產生的聲音（示例興趣點見圖 2）進行可視化研究，由此協助工程師制定工程決策，并探索頭部和軀體對總聲壓的影響（圖3）。以米為量級的尺度下進行聲學模擬的挑戰在于，常用的有限元法（finite element method，簡稱 FEM）的計算量會很大并且占用大量內存。

圖 2. 上和中：“近耳”點，在這些位置計算了口對點的近耳傳輸函數，這是一種更具體的關于頭部的傳輸函數。下：3200 Hz 頻率下的聲壓分布。紅色表示高正壓，藍色表示高負壓。

圖 3. 3200 Hz 頻率、1 m 半徑的總聲壓極坐標圖。

對于在筆記本電腦上執行大部分模擬工作的 Christensen 而言，計算量與內存需求是一個嚴重的桎梏，但是 COMSOL Multiphysics? 軟件讓他能夠隨意選擇適合的方法。針對這一案例，他可以充分利用軟件中的聲學邊界元法，實現更為高效的模擬。雖然邊界元法在每個自由度上的計算量超過了有限元法，但邊界元法在大體積范圍內實現相同的精度時，所需的自由度要比有限元法少很多。利用邊界元法，可以在域內的任一點上提取聲壓值，而僅需對表面進行網格劃分和計算。顯而易見，這項功能具備實質性的計算優勢。有限元法需要對整個體進行網格劃分，更加適用于近場分析。通過使用邊界元法，Christensen得以減少模擬的計算量，讓筆記本電腦的計算資源足以應付此類仿真運算?！斑吔缭ê唵味嗔?。”Christensen 評價道，“如果你擁有完整的幾何結構，那就只需要創建表面網格，這種工作方法簡直是一種享受。通常人們會使用體網格，他們不得不離散空氣域，還得應用輻射條件來限制反射，但這些特征都已經包含在邊界元法當中?！?/p>

講話與傾聽

HATS 的優點在于，它可以模擬包含多個聲源的環境。舉例來說，用戶在來訪時的說話聲，就像在忙碌的辦公室中戴著耳機工作的職員說話時的情況一樣。而助聽器中由于配置有麥克風和揚聲器，復雜的結構可能會引起聲反饋?！爸犉髦醒b有兩個麥克風，而且它們均可感知聲音和振動，這讓問題更加復雜。”Christensen 解釋道?！爸犉鳟a生的聲反饋現象有時會造成令人難以忍受的噪音?！?/p>

因為邊界元法需要極少的物理場設置，所以建模相對簡單。在上述案例中，研究人員對口腔采用加速度邊界條件以模擬揚聲器。假設 HATS 自身表現為剛性，并且邊界元法（BEM）可自動添加任何輻射條件，確保聲波不會反射?；谏鲜鰲l件，計算頭部和軀干周圍的聲壓場，運行仿真后可以通過后處理對多個麥克風位置進行研究?；谥犉鞯膸缀谓Y構，機械設計師可以指出所有可行的位置，并且根據仿真結果推測最佳位置。“借助邊界元法，我們確實不需要考慮過多物理場設置。”Christensen 說道，“因為我們感興趣的是口與耳的關系，所以在口的位置施加了一個加速度條件，并在口的前方，或耳朵內部和表面上提取聲壓，僅此而已。我們需要關心的唯一物理場就是聲場。”

未來的應用

Christensen 的 HATS 模型對于他和同事來說具有極高的價值，因為該模型一旦完成求解，就可以通過各種后處理工具提取海量的有用信息。他們可以對不同的配置（比如有助于助聽器緊密貼合的拱形結構，以及連接內耳和外耳的端口）進行分析研究，針對各種位置推導出新的傳輸函數，也可以與其他軟件搭配使用。

Christensen 還表示，他將繼續使用 COMSOL Multiphysics? 改進他的聲學模型，不僅僅是 HATS，他還將使用該軟件研究各種振動聲學問題。“我喜歡這款軟件，因為它具有易于使用且直觀的界面，還支持根據自己的需求添加或修改方程，我經常使用這項功能?！彼u價道，“這款軟件對于像我這樣試圖通過非現有的數學方法來描述物理現象的人來說非常有用。我經常自問：‘你可以給出描述此問題的方程嗎？’如果答案是肯定的，那就能夠在 COMSOL 中建立并求解這個問題?！?br>

本文內容來自多物理場仿真:《IEEE Spectrum》特刊 2018