一種平衡銜鐵受話器數(shù)值仿真分析方法

蔣 鈺，陸 曉，溫周斌

(1.中國科學(xué)院聲學(xué)研究所東海研究站，上海 201815；2.中國科學(xué)院聲學(xué)研究所，北京 100049；3.中國科學(xué)院大學(xué)，北京 100190；4.浙江中科電聲研發(fā)中心，浙江 嘉善 314115)

0 引言

平衡銜鐵受話器(Balanced Armature Receiver,BAR)，具有尺寸小、電聲轉(zhuǎn)換效率高和靈敏度高等特點(diǎn)[1-2]，已被廣泛用于助聽設(shè)備、入耳式耳機(jī)和各種軍用頭戴式耳機(jī)等產(chǎn)品中[3]。近年來，隨著消費(fèi)者對聲品質(zhì)要求提高，平衡銜鐵受話器的應(yīng)用更為廣泛。

數(shù)值仿真分析是一種通過計(jì)算機(jī)進(jìn)行數(shù)值計(jì)算求解一組滿足一定邊界條件的微分方程，并得到其近似解的方法[4]，已廣泛應(yīng)用于動圈式揚(yáng)聲器單體及系統(tǒng)的主要特性的仿真分析[5-8]。近10年來，國內(nèi)外學(xué)者圍繞平衡銜鐵受話器的非線性、聲壓級和失真等主要特性進(jìn)行了深入研究。Jensen等基于集總參數(shù)模型深入研究了平衡銜鐵受話器的非線性和失真等特性[9]?；诩倕?shù)模型的仿真分析方法概念清晰、方法簡便，但所得到的仿真結(jié)果在中高頻段精度不夠高[10]。Jensen還采用數(shù)值仿真分析方法研究了不同平衡銜鐵磁導(dǎo)率對應(yīng)的平衡銜鐵上的磁通密度分布[11]。Hu等創(chuàng)建了平衡銜鐵受話器的振動系統(tǒng)與聲場的數(shù)值仿真分析模型，驗(yàn)證了其所提出的多模態(tài)模型的有效性[10]。Ziolkowski等將平衡銜鐵受話器磁路簡化為2D模型，通過耦合電磁場與結(jié)構(gòu)，得到了平衡銜鐵的位移[12]。Xu等提出了平衡銜鐵受話器電磁場的瞬態(tài)分析和基于模態(tài)疊加法的平衡銜鐵受話器結(jié)構(gòu)與聲場的頻域分析相結(jié)合的方法，利用前者得到的非線性參數(shù)和后者得到的力與聲壓級的傳遞函數(shù)，再通過數(shù)值迭代最終得到平衡銜鐵受話器輸出的聲壓級[13-14]。

本文首先介紹平衡銜鐵受話器及其集總參數(shù)模型，分析指出平衡銜鐵運(yùn)動會對揚(yáng)聲器特性帶來重大影響；然后提出一種基于COMSOL Multiphysics數(shù)值仿真分析軟件的平衡銜鐵受話器聲壓級和阻抗的仿真分析方法。該方法首先對平衡銜鐵受話器的磁場特性進(jìn)行穩(wěn)態(tài)分析，計(jì)算得到可隨平衡銜鐵運(yùn)動變化的驅(qū)動力和反電動勢；然后對電磁場、結(jié)構(gòu)和聲場耦合模型進(jìn)行頻域分析，最終計(jì)算得到平衡銜鐵受話器的聲壓級和阻抗。

1 平衡銜鐵受話器及其集總參數(shù)模型

圖1是一款平衡銜鐵受話器的半剖幾何模型[1]。由圖1可知，平衡銜鐵受話器由線圈(Coil)、傳動桿(Drive pin)、振膜(Membrane)、永磁體(Magnets)、磁軛(Magnet house)、平衡銜鐵(Armature)、前蓋(Cover)、外殼(Case)和出聲孔(Sound outlet)等組成。平衡銜鐵一端被固定在外殼上，另一端懸空。平衡銜鐵通過驅(qū)動桿與振膜連接。振膜被粘于上蓋，其中三邊有折環(huán)，另外一邊是金屬。平衡銜鐵和振膜的振動都是類似于懸臂梁的振動。振膜和上蓋之間的空氣域定義為前腔，振膜另一面的空氣域?yàn)楹笄弧?/p>

圖1 平衡銜鐵受話器半剖幾何模型[1]Fig.1 Section model of BAR[1]

當(dāng)平衡銜鐵受話器音圈中無電流通過時(shí)，平衡銜鐵在磁場中處于受力平衡狀態(tài)，因而穩(wěn)定在兩個(gè)永磁體中間。當(dāng)線圈中有電流通過時(shí)，在線圈周圍除了穩(wěn)定的磁場之外，還有因電流變化所感應(yīng)出的擾動磁場。平衡銜鐵會被擾動磁場磁化，它的磁化方向也會隨著流經(jīng)線圈中電流方向的改變而改變。因此，平衡銜鐵所受到的力將不斷變化，這個(gè)力會使得平衡銜鐵上下往復(fù)運(yùn)動。平衡銜鐵往復(fù)運(yùn)動將通過傳動桿帶動振膜上下振動，而振膜的上下振動又會使得前腔的空氣壓縮或膨脹，從而產(chǎn)生聲波，聲波再從出聲孔傳遞出去。

建立平衡銜鐵受話器的集總參數(shù)模型，可分析得到諸多重要參數(shù)，進(jìn)而還可分析得到這些重要參數(shù)對揚(yáng)聲器特性的影響。圖2(a)給出了一個(gè)平衡銜鐵受話器的簡化的磁回路。將磁阻、磁場力和磁通量分別等效為電阻Rg、電壓源FM和電流Φ，可建立該磁回路的集總參數(shù)模型，圖2(b)是圖2(a)的等效電路圖[10]。利用該等效電路圖，可以得到平衡銜鐵中的磁通量和上下磁隙中的磁通量，進(jìn)而計(jì)算得到線圈上所受到的反電動勢和平衡銜鐵所受的力。

圖2 平衡銜鐵受話器簡化的磁回路及其等效電路[9]Fig.2 Simplified magnetic circuit and equivalent circuit of BAR[9]

反電動勢Uback可表示為[11]

平衡銜鐵所受的力，該力F可表示為[11]

式中，kΦ(x)為非線性補(bǔ)償剛度系數(shù)，Tme(x)和Tme，d(x，i)分別為非線性傳遞系數(shù)和失真?zhèn)鬟f系數(shù)。由式(2)可知，當(dāng)揚(yáng)聲器工作時(shí)，平衡銜鐵所受的力不僅與流經(jīng)線圈的電流有關(guān)，還與平衡銜鐵的位移相關(guān)。與位移相關(guān)的力可稱為動生磁場力。由于平衡銜鐵是導(dǎo)磁材料，其運(yùn)動會導(dǎo)致?lián)P聲器磁場特性變化，磁場特性的變化會導(dǎo)致平衡銜鐵受到動生磁場力的作用。

平衡銜鐵所受的力隨位移變化的特性會對揚(yáng)聲器性能產(chǎn)生重要影響。平衡銜鐵的有效剛度是其機(jī)械剛度與非線性補(bǔ)償剛度系數(shù)的差值，非線性補(bǔ)償剛度系數(shù)增加，會減小平衡銜鐵的有效剛度，從而減小揚(yáng)聲器的共振頻率；該系數(shù)的減小則會使振膜振幅增大，聲壓也會增大[11]。而該系數(shù)的增加會增加揚(yáng)聲器的不穩(wěn)定性，這是因?yàn)橛行偠葴p小會使得平衡銜鐵可能無法回到平衡位置。

2 數(shù)值仿真分析方法

由前文平衡銜鐵受話器集總參數(shù)模型分析可知，平衡銜鐵的運(yùn)動會導(dǎo)致?lián)P聲器的磁場特性發(fā)生變化，在線圈上產(chǎn)生一個(gè)額外的反電動勢，并使得平衡銜鐵受到一個(gè)額外的磁場力作用。采用數(shù)值仿真分析方法仿真計(jì)算揚(yáng)聲器的聲場特性時(shí)一般選用頻域分析，直接的頻域分析則無法考慮平衡銜鐵運(yùn)動對磁場特性的影響。

為此，本文提出一種基于COMSOL Multiphysics數(shù)值仿真分析軟件的平衡銜鐵受話器數(shù)值仿真分析方法，分為兩個(gè)步驟：

(1)對平衡銜鐵受話器的磁場特性進(jìn)行穩(wěn)態(tài)分析，并利用動網(wǎng)格技術(shù)，計(jì)算得到可隨平衡銜鐵運(yùn)動變化的作用在線圈上的反電動勢(速度電動勢)和施加在平衡銜鐵上的動生磁場力。

(2)創(chuàng)建平衡銜鐵受話器電磁場、結(jié)構(gòu)和聲場的耦合分析模型，并利用前文穩(wěn)態(tài)求解得到的可隨平衡銜鐵速度變化的速度電動勢和可隨平衡銜鐵位移變化的動生磁場力，再進(jìn)行頻率分析，最終計(jì)算得到揚(yáng)聲器的聲壓級和阻抗。

無論是穩(wěn)態(tài)分析，還是頻域分析，數(shù)值仿真分析主要包括建立幾何模型、定義材料參數(shù)、設(shè)置邊界條件、網(wǎng)格劃分、計(jì)算和結(jié)果后處理等步驟。

由于影響仿真分析結(jié)果準(zhǔn)確性的主要因素包括幾何模型、材料參數(shù)和仿真分析方法等三個(gè)部分，因此，下文以一款平衡銜鐵受話器為例，首先簡要介紹它的幾何模型和材料參數(shù)，然后再詳細(xì)介紹它的聲壓級和阻抗的數(shù)值仿真分析方法。

2.1 幾何模型和材料參數(shù)

圖3給出了一款平衡銜鐵受話器的幾何模型。平衡銜鐵受話器具有幾何對稱性，因此只需建立其1/2幾何模型即可。為了提高計(jì)算效率，模型中將帶有花紋的金屬球頂簡化為平面，磁隙間距為0.1mm。

圖3 一款平衡銜鐵受話器的幾何模型(1/2模型)Fig.3 Half geometry model of BAR

表1和表2分別給出了該揚(yáng)聲器各個(gè)部件的磁性材料參數(shù)和力學(xué)材料參數(shù)。由于在模型中簡化了振膜的花紋，所以需要適當(dāng)提高振膜的楊氏模量，使得簡化后的振膜剛度可以等效包含花紋的振膜剛度。圖4給出了平衡銜鐵受話器磁軛(軟磁材料)的BH曲線。其中，B為磁通密度，H為磁場強(qiáng)度，線圈匝數(shù)為350匝。

表1 平衡銜鐵受話器部件的磁性材料參數(shù)Table 1 Magnet material properties of BAR

表2 平衡銜鐵受話器部件的力學(xué)材料參數(shù)Table 2 Mechanical material properties of BAR

圖4 磁軛(軟磁材料)的BH曲線Fig.4 BH curve of soft iron

假設(shè)所有材料是各向同性的。對于磁性材料，需要給出磁通密度B和磁場強(qiáng)度H之間的本構(gòu)關(guān)系，用剩磁Br和回復(fù)磁導(dǎo)率μrec定義永磁體本構(gòu)關(guān)系[13-14]：

式中：μ0為真空磁導(dǎo)率。磁軛的磁通密度和磁場強(qiáng)度之間的關(guān)系是非線性的，需用完整的BH曲線來定義磁軛的本構(gòu)關(guān)系。平衡銜鐵的磁通密度和磁場強(qiáng)度之間的關(guān)系也是非線性的，但本文將其簡化，僅用相對磁導(dǎo)率μr(常數(shù)值)來定義平衡銜鐵的本構(gòu)關(guān)系[13-14]：

對于平衡銜鐵、傳動桿和振膜等振動系統(tǒng)部件，均采用線彈性材料模型定義它們的應(yīng)力應(yīng)變關(guān)系。

2.2 平衡銜鐵運(yùn)動對揚(yáng)聲器特性的影響

為了仿真分析平衡銜鐵運(yùn)動對磁場特性的影響，平衡銜鐵受話器的幾何模型只需包含磁路和平衡銜鐵，忽略振膜和傳動桿等其他部件；假設(shè)金屬外殼有一定的磁屏蔽作用，因而忽略外殼以外的空氣域，僅建立外殼內(nèi)的空氣域。材料參數(shù)也只需定義磁路系統(tǒng)各個(gè)部件的磁性材料參數(shù)和平衡銜鐵的力學(xué)材料參數(shù)。在平衡銜鐵與外殼相接面上施加固定約束條件，在對稱面上，通過定義磁矢勢的面內(nèi)分量為0和位移的法向分量為0，來分別設(shè)定磁場和結(jié)構(gòu)的對稱條件。在劃分網(wǎng)格時(shí)，需對平衡銜鐵適當(dāng)加密網(wǎng)格。

在位于磁隙區(qū)域的平衡銜鐵表面上添加力載荷(如圖5所示)，平衡銜鐵受力后會發(fā)生形變，可根據(jù)力的平衡方程計(jì)算得到平衡銜鐵的位移大小。平衡銜鐵作為導(dǎo)磁材料，其形變會改變磁場分布。為了模擬平衡銜鐵運(yùn)動(形變)對磁場分布的影響，模型中需要使用動網(wǎng)格技術(shù)，即網(wǎng)格將隨著平衡銜鐵的運(yùn)動(形變)而變化。基于麥克斯韋方程組計(jì)算得到磁通密度分布。在平衡銜鐵表面上加載不同大小的力，便可得到平衡銜鐵運(yùn)動到不同位置時(shí)的磁通密度分布。

圖5 平衡銜鐵的加載面圖Fig.5 The loading surface diagram of the balanced armature

圖6(a)和圖6(b)分別給出了當(dāng)平衡銜鐵處于初始位置時(shí)及其從初始位置向上移動至磁隙間距的1/10位置時(shí)線圈所包圍區(qū)域的磁通密度分布圖。由圖6可以看出，平衡銜鐵位置(位移)的變化，會導(dǎo)致該區(qū)域磁通密度發(fā)生明顯改變。磁通密度改變，會影響作用在線圈上感應(yīng)電動勢和平衡銜鐵上所受到的力。

圖6 平衡銜鐵處于不同位置時(shí)的磁通密度分布圖Fig.6 Diagrams of the magnetic flux density distribution when the balanced armature is at different positions

2.2.1 作用在線圈上的速度電動勢

由式(1)可知，作用在線圈上的速度電動勢Uv為

根據(jù)電磁感應(yīng)定律，速度電動勢Uv是磁通量Φ(x)對時(shí)間t的導(dǎo)數(shù)：

在平衡銜鐵的位移比較小的情況下，假設(shè)線圈中的磁通量Φ(x)與平衡銜鐵位移x呈線性關(guān)系；根據(jù)式(5)和式(6)，可以將Φ(x)與x的關(guān)系表示為

式中，系數(shù)Tme和a均為常數(shù)。當(dāng)位移x較小時(shí)，近似認(rèn)為非線性傳遞系數(shù)Tme(x)是一個(gè)不隨x改變的常數(shù)Tme。

從前述穩(wěn)態(tài)求解得到的仿真結(jié)果中選取幾個(gè)平衡銜鐵的位置，并根據(jù)平衡銜鐵處于這些位置時(shí)的磁場分布計(jì)算線圈所包圍面積的磁通量，再通過式(7)進(jìn)行函數(shù)擬合，可以得到系數(shù)Tme和a。

將式(7)代入式(6)可得：

2.2.2 平衡銜鐵上所受到的動生磁場力

平衡銜鐵表面的電磁力密度矢量分量f為[14]：

式中，n是表面單位法向矢量，H2=H2x+H2y+H2z。

圖7給出了平衡銜鐵在平衡位置時(shí)其所受的表面電磁力密度矢量分量的分布圖，由圖7可知，該作用力主要集中在磁隙區(qū)域?qū)?yīng)的表面上。

圖7 平衡銜鐵上表面應(yīng)力分布Fig.7 Surface stress distribution of balanced armature

當(dāng)平衡銜鐵運(yùn)動時(shí)，它與上下永磁體之間的距離會發(fā)生變化，距離越小，磁通密度越大，靠近永磁體的平衡銜鐵表面應(yīng)力會增大，遠(yuǎn)離永磁體的平衡銜鐵表面應(yīng)力會減小。

由式(9)可計(jì)算得到平衡銜鐵表面應(yīng)力(表面上每個(gè)節(jié)點(diǎn)的應(yīng)力)隨著位移變化的變化量，再對該表面應(yīng)力的變化量求積分，即可得到(整個(gè))平衡銜鐵上所受到的力的變化量，也即動生磁場力FΦ：

式中：Δf是應(yīng)力的變化量，S表示平衡銜鐵表面。

由式(2)可知，平衡銜鐵上受到的動生磁場力與位移相關(guān)，當(dāng)平衡銜鐵的位移非常小時(shí)，假設(shè)動生磁場力FΦ與位移x呈線性關(guān)系，即：

式中：系數(shù)kΦ為常數(shù)。當(dāng)位移x較小時(shí)，近似認(rèn)為非線性補(bǔ)償剛度系數(shù)kΦ(x)是一個(gè)不隨x改變的常數(shù)kΦ。

同理，從前述穩(wěn)態(tài)求解得到的仿真結(jié)果選取幾個(gè)平衡銜鐵的位置，并根據(jù)式(10)計(jì)算對應(yīng)的動生磁場力，再通過式(11)進(jìn)行函數(shù)擬合，計(jì)算得到系數(shù)kΦ。

2.3 磁路、振動系統(tǒng)和聲場的耦合分析

創(chuàng)建包含磁路和振動系統(tǒng)所有部件的平衡銜鐵受話器的幾何模型，因?yàn)樗哂袔缀螌ΨQ性，所以只建立其1/2幾何模型即可。由于本案例的平衡銜鐵受話器是作為受話器使用，因此，在仿真模型中，其聲輸出方式為711耦合器(壓力場環(huán)境)。為提高計(jì)算效率，可采用特定的“阻抗邊界”條件替代“711耦合器聲腔(模型)”的簡化方法，因而在仿真模型中并不需要建立耦合器聲腔模型。圖8(a)是包含耦合器聲腔的全模型，圖8(b)則是簡化后的幾何模型。對于圖8(b)，只需要包含平衡銜鐵受話器和小部分聲腔，這樣可以大幅降低計(jì)算量。采用表1和表2中的材料參數(shù)定義所創(chuàng)建模型中揚(yáng)聲器各個(gè)部件的材料參數(shù)。

圖8 平衡銜鐵受話器幾何模型(含711耦合器)Fig.8 Geometry model of BAR(including 711 coupler)

在平衡銜鐵與外殼相接面上施加固定約束。在對稱面上，通過定義磁矢勢的面內(nèi)分量為0、位移的法向分量為0和加速度的法向分量為0，分別設(shè)定磁場、結(jié)構(gòu)和聲場的對稱條件?？諝庥蛑谐伺c振動系統(tǒng)耦合面以外的邊界面，全部假設(shè)為全反射邊界。在小部分聲腔與711耦合器的連接面上設(shè)定聲阻抗，可等效711耦合器空氣域?qū)吔缑嫔下晧寒a(chǎn)生的影響，這個(gè)聲阻抗Zin可以預(yù)先利用圖8(a)全模型仿真分析得到，計(jì)算公式為

式中：Pin為輸入端面處的聲壓有效值，Zin是與輸入端面處空氣體積流速Q(mào)in相關(guān)的聲阻抗。合理劃分網(wǎng)格，以保證計(jì)算頻率范圍內(nèi)一個(gè)波長至少包含4～5個(gè)網(wǎng)格。

平衡銜鐵受話器的磁路、振動系統(tǒng)和聲場是相互影響的，在模型中需要正確定義它們的耦合作用。

對于磁場和結(jié)構(gòu)這兩個(gè)物理場，實(shí)際加載在線圈上的電壓U是由外部加載電壓U0和反電動勢Uback兩部分組成，而反電動勢包含由于電流變化帶來的反電動勢Ui和由于平衡銜鐵運(yùn)動帶來的反電動勢(速度電動勢)Uv，即：

平衡銜鐵所受到的力F是由線圈電流產(chǎn)生的感生磁場力Fi和因其自身運(yùn)動而產(chǎn)生的動生磁場力FΦ兩部分組成的，即：

對于結(jié)構(gòu)和聲場這兩個(gè)物理場，振膜表面的加速度即為聲源。同時(shí)，作用在振膜表面上的空氣壓強(qiáng)會使得振動系統(tǒng)受到額外的表面應(yīng)力。

在線圈兩端加載諧波擾動的電壓U0，電壓有效值為0.1 V。根據(jù)電磁感應(yīng)定律，變化的電場會感應(yīng)出磁場，磁路系統(tǒng)中除了原本的靜態(tài)磁場以外，會再疊加一個(gè)諧波擾動的磁場。根據(jù)麥克斯韋方程組求解得到電場分布和擾動的磁場分布。根據(jù)式(12)可計(jì)算得到平衡銜鐵所受到的感生磁場力Fi。平衡銜鐵受力產(chǎn)生振動，根據(jù)結(jié)構(gòu)運(yùn)動方程[15]求解得到振動幅值。根據(jù)振動系統(tǒng)與聲場的邊界面上振幅(或加速度)，可計(jì)算得到聲場分布。一般空氣域的聲場分析建立在理想流體媒質(zhì)的波動方程[15]之上。而對于平衡銜鐵受話器，在它的內(nèi)部存在一些狹窄區(qū)域，不能忽略聲波在狹窄區(qū)域中的粘滯特性。

如前文所述，為提高計(jì)算效率，幾何模型中并沒有包含711耦合器空氣域，因此無法直接獲取麥克風(fēng)面處的聲壓級。利用已建立的幾何模型聲腔邊界面(也即711耦合器輸入端面)上的空氣體積流速Q(mào)in和711耦合器空氣域的轉(zhuǎn)移阻抗Ztrans，可計(jì)算得到麥克風(fēng)表面位置的聲壓pmic：

圖9給出了該平衡銜鐵受話器在頻率為1 kHz信號激勵(lì)下的聲壓級分布。由圖9可知，在振膜與上蓋之間的腔體內(nèi)的聲壓最高，聲波通過出聲孔向外(711耦合器)傳遞。

圖9 平衡銜鐵受話器中的聲壓級分布圖(@1 kHz)Fig.9 Sound pressure level(SPL)of BAR(@1 kHz)

3 仿真分析結(jié)果

通過本文提出的數(shù)值仿真分析方法可以較快地仿真計(jì)算得到平衡銜鐵受話器的聲壓級和阻抗曲線。與此同時(shí)，作者還采用B&K Pulse系統(tǒng)測量了所研究的平衡銜鐵受話器的聲壓級和阻抗曲線，為保證測量與仿真分析條件的一致，測量是在IEC711耦合器(愛宏AWA6162耳模擬器)壓力場條件下進(jìn)行的，符合IEC60318-4標(biāo)準(zhǔn)[16]。被測揚(yáng)聲器通過夾具安裝在耳模擬器上，揚(yáng)聲器和夾具之間使用少量橡皮泥以防止泄露，圖10(a)是被測平衡銜鐵受話器，外尺寸為6.3 mm×4.3 mm×3.0 mm。圖10(b)是被測揚(yáng)聲器通過夾具安裝在耳模擬器上的照片。圖10(c)給出了測量示意圖，揚(yáng)聲器出聲孔平面應(yīng)與耳模擬器參考平面齊平。

圖10 平衡銜鐵受話器的測量Fig.10 Measurement of BAR

圖11(a)和圖11(b)分別給出了平衡銜鐵受話器的聲壓級和阻抗的仿真分析結(jié)果和測量結(jié)果。由圖11可知，在4 kHz以下，仿真分析結(jié)果和測量結(jié)果基本吻合；在高頻段，兩者趨勢一致，數(shù)值存在一定差異。

圖11 仿真結(jié)果和測量結(jié)果Fig.11 Comparison between measurement and simulation

由圖11(a)可知，在低頻段阻抗值基本與直流阻抗一致；而在高頻段阻抗值隨頻率增大而增大。這是因?yàn)橛呻娏髯兓a(chǎn)生的感應(yīng)電動勢增大的緣故。感應(yīng)電動勢越大，電流越小，阻抗就越大。2 kHz頻率附近的峰以及4、7 kHz頻率附近的峰都與速度電動勢有關(guān)，由于振動系統(tǒng)在這些頻率附近發(fā)生共振(或諧振)，平衡銜鐵的振幅相對較大，平衡銜鐵運(yùn)動引起的位置改變所帶來的速度電動勢也增大。同樣，速度電動勢越大，電流越小，阻抗也就越大。

由圖11(b)可知，在低頻段聲壓級基本為一穩(wěn)定值，當(dāng)頻率大于1 kHz，聲壓級曲線出現(xiàn)多個(gè)諧振峰。低頻聲壓級與振動系統(tǒng)的有效剛度密切相關(guān)，而動生剛度是有效剛度的重要組成部分。動生剛度越大，有效剛度越小，振幅則越大，低頻聲壓級也就越高。2 kHz和7 kHz附近的諧振峰也與振動系統(tǒng)的共振(或諧振)相關(guān)，這些峰所對應(yīng)的頻率高低與動生剛度有關(guān)。動生剛度越大，有效剛度則越小，這些峰對應(yīng)的頻率也會越低。而14 kHz頻率附近的共振峰則與此無關(guān)，它是由于711耦合腔的共振所引起的。

本文所提方法較為準(zhǔn)確地仿真了可隨平衡銜鐵速度變化的速度電動勢和可隨平衡銜鐵位移變化的動生磁場力。分析得到了阻抗共振峰和諧振峰，最終較為準(zhǔn)確地仿真得到了低頻聲壓級和聲壓級頻響曲線上的各個(gè)峰。

4 結(jié)論

本文基于平衡銜鐵受話器的集總參數(shù)模型分析，指出平衡銜鐵運(yùn)動會產(chǎn)生額外的反電動勢和磁場力，提出了一種分步進(jìn)行的平衡銜鐵受話器的數(shù)值仿真方法。首先對平衡銜鐵受話器的磁場特性進(jìn)行穩(wěn)態(tài)分析，計(jì)算得到可隨平衡銜鐵運(yùn)動而變化的磁場力和反電動勢；然后，建立電磁場、結(jié)構(gòu)和聲場耦合分析模型，并利用前面的穩(wěn)態(tài)求解方法得到的可隨平衡銜鐵運(yùn)動變化的速度電動勢和動生磁場力，再進(jìn)行頻域分析，最終計(jì)算得到平衡銜鐵受話器的聲壓級和阻抗曲線。仿真結(jié)果和測量結(jié)果基本一致。

平衡銜鐵受話器的數(shù)值仿真分析方法還可以在以下幾個(gè)方面進(jìn)一步優(yōu)化和完善：(1)深入研究平衡銜鐵材料的磁特性，以建立更準(zhǔn)確的非線性模型；(2)細(xì)致分析小體積聲腔中的空氣壓強(qiáng)對振動系統(tǒng)的影響；(3)非線性和失真特性機(jī)理和高效仿真分析方法研究。