COMSOLMultiphysics微型揚(yáng)聲器分析設(shè)計項目解決方案-郭枝權(quán)

..oftWare軟件世界..COMSOLMultiphysics微型揚(yáng)聲器分析設(shè)計解決方案□中仿科技〔CnTech公司郭枝權(quán)安琳..一、背景分析在中國乃至世界范圍內(nèi),消費(fèi)數(shù)碼產(chǎn)業(yè)在過去十年間發(fā)展迅猛。20XX前后,MP3出現(xiàn)。MP3引起的風(fēng)潮刺激了市場,帶來了需求,也催生了更多諸如MP4、MP5這樣的產(chǎn)品。今天,智能手機(jī)、平板電腦、筆記本電腦都把影音播放功能作為賣點(diǎn)之一,價格幾百乃至幾千的耳塞／耳機(jī)在公交車或者地鐵上隨處可見。就連手機(jī)的外放和配套耳機(jī),也開始注重音效。另一方面,汽車消費(fèi)在中國的強(qiáng)勁增長,更是大大推動了小型揚(yáng)聲系統(tǒng)產(chǎn)業(yè)的發(fā)展和擴(kuò)大。當(dāng)前世界正處于刺激消費(fèi)、拉動增長的主旋律之中。中國市場如果排除投資性增長的因素,國民消費(fèi)幾乎全部集中于房產(chǎn)、教育、醫(yī)療、汽車、數(shù)碼電子這幾大領(lǐng)域。微型揚(yáng)聲器產(chǎn)業(yè)面對的市場異常龐大而且充滿機(jī)會。1.揚(yáng)聲器分析技術(shù)背景揚(yáng)聲器設(shè)計離不開聲學(xué)分析,從技術(shù)上說,這是一個波動問題的求解過程,在我們選擇分析方法的時候,頻率〔或者是波長是一個很重要的參考量。如果分析的結(jié)構(gòu)尺寸遠(yuǎn)遠(yuǎn)大于波長,這個時候波動性的體現(xiàn)并不明顯,工程計算一般采用路徑追跡的方法,也才有人把這稱為幾何聲學(xué)的方法。即認(rèn)為聲波在均勻介質(zhì)中沿直線傳播,并按照已知的立體角發(fā)散；遇有介質(zhì)突變的界面時,聲波按照可預(yù)測的角度發(fā)生反射,同時存在可描述的功率衰減。這種追跡的方法,業(yè)內(nèi)通常喜歡使用Matlab一類的工具,用幾行代碼即可實現(xiàn)。如果分析的結(jié)構(gòu)尺寸遠(yuǎn)遠(yuǎn)小于波長,這個時候的波動性也不明顯。彈性波動現(xiàn)象在這個時候弱化為彈性振動,問題變成了結(jié)構(gòu)的動態(tài)力學(xué)分析過程。只有當(dāng)分析的結(jié)構(gòu)尺寸與波長可比擬時,波動效應(yīng)〔衍射、散射體現(xiàn)得最為突出?？諝庵械穆曀贋?40m/s,一般的樂聲和人聲的頻率范圍約為200~4000Hz,換算成波長即為85mm~1.7m。這正是我們常見的微型、小型揚(yáng)聲器,乃至家用級音箱的結(jié)構(gòu)尺寸范圍。這時,需要求解波動方程來進(jìn)行分析。除聲學(xué)之外,揚(yáng)聲器中還涉及電磁線圈的頻率分析、功耗分析和結(jié)構(gòu)失效分析等。2SOLMultiphysics應(yīng)用背景微型揚(yáng)聲器設(shè)計工作涉及三個基礎(chǔ)領(lǐng)域：電磁設(shè)計、結(jié)構(gòu)設(shè)計和聲學(xué)設(shè)計。業(yè)內(nèi)絕大多數(shù)的廠商,其研發(fā)中心也基本照此劃分為三個團(tuán)隊。顯然,這是一個多物理耦合的工程應(yīng)用。另外,工業(yè)級的應(yīng)用對計算分析在很多細(xì)節(jié)上有著很高的要求。比如耳機(jī)結(jié)構(gòu)中存在動圈結(jié)構(gòu)和音膜,二者尺寸跨度很大。同時音膜的結(jié)構(gòu)往往非常復(fù)雜,各部分厚度可能不同、材質(zhì)可能不同,同時存在多層膜復(fù)合的情況,這些都需要仿真工具有專門的解決方案。COMSOLMultiphysics就是這樣一個專門的解決方案,而且就目前而言,是世界上唯一能夠完美解決微型揚(yáng)聲器設(shè)計中所有問題的工業(yè)級解決方案。在微型揚(yáng)聲器設(shè)計領(lǐng)域,COMSOLMultiphysics是行業(yè)標(biāo)準(zhǔn)。二、微型揚(yáng)聲器設(shè)計分析要點(diǎn)從工業(yè)應(yīng)用的角度,用戶的典型需求如下。〔1結(jié)構(gòu)微型化：小型電子設(shè)備〔如手機(jī)中使用的微型揚(yáng)聲器的尺寸相對來說較小,在毫米〔局部微米級的結(jié)構(gòu)上是否能夠?qū)崿F(xiàn)多種物理現(xiàn)象的模擬？能否保證精度？〔2對于電磁線圈的模擬與分析：在揚(yáng)聲器的線圈中通入電流后,電場會產(chǎn)生磁場,同原有永磁體的磁場相互影響,同時線圈存在感應(yīng)電動勢。這一系列的模擬能否提供明確的模擬方法和解決方案？能否得到符合實際的驅(qū)動力？能否給出線圈的阻抗、電感的頻率曲線？〔3在揚(yáng)聲器中所使用的音膜為復(fù)合材料,是由兩層或多層不同材料的薄膜構(gòu)成,每層膜的不同局部,采用不同的厚度和材質(zhì),如何在厚度很小的情況下模擬出符合實際情況的音膜整體性能？..32idnovo..欄目主持：侯琳投稿信箱：houlidnovo〔4音膜在實際生產(chǎn)中使用時,在其表面會有些通過模擬薄膜產(chǎn)生聲場在無限大開放環(huán)境中的傳播。通過進(jìn)行頻蝕刻或成型制造出來的紋路,通過紋路的位置和深淺來調(diào)率掃描計算,得到結(jié)構(gòu)內(nèi)部空氣響應(yīng)和振動模態(tài)。根據(jù)這一節(jié)音膜的響應(yīng)頻率,在模擬中是否能實現(xiàn)對存在紋路的音思路,設(shè)計人員可以調(diào)整任何參數(shù),如改變驅(qū)動電信號、選膜模擬？擇材料和設(shè)計結(jié)構(gòu)等,直到達(dá)到最佳的結(jié)果?！?聲學(xué)模擬：揚(yáng)聲器聲學(xué)特性的模擬,如何實現(xiàn)？下面通過一系列案例說明這一過程。由于我們需要保無限遠(yuǎn)的區(qū)域如何實現(xiàn)？護(hù)相關(guān)廠商的技術(shù)秘密,因此本方案中的模型均為簡化過的〔6CAD模型導(dǎo)入：由于微型揚(yáng)聲器結(jié)構(gòu)的較為復(fù)模型,并有意把模型分離成各小部分分別演示,同時隱去廠雜,尤其是薄膜結(jié)構(gòu)蝕刻或成型后的結(jié)構(gòu),如果使用不是很商的信息和重要的技術(shù)參數(shù),如結(jié)構(gòu)細(xì)節(jié)尺寸、各部分材料熟悉的軟件進(jìn)行建模是比較困難的,同時對于已有CAD模型特性等,旨在說明COMSOLMultiphysics如何進(jìn)行分析。文件,是否能夠直接導(dǎo)入CAE軟件中進(jìn)行計算？對于通過讀1.揚(yáng)聲器驅(qū)動分析取CAD文件導(dǎo)入到CAE軟件中的模型能否在軟件中進(jìn)行修改？圖1展示了一個典型的動圈揚(yáng)聲器的基本結(jié)構(gòu)。永磁如何修改？修改后的模型是否會影響模擬的結(jié)果？鐵產(chǎn)生的磁場被極片、頂板和底板所收集,然后在氣隙處〔7經(jīng)驗方程或特殊材料屬性的輸入：對于一些研究集中釋放。氣隙中的音圈被一個面彈簧懸掛。當(dāng)音圈通有人員自己的一些經(jīng)驗方程能否直接輸入到軟件中,實現(xiàn)求變化的電流〔音頻信號時,磁場產(chǎn)生的羅侖茲力和面彈解？對于一些特殊的或新研發(fā)使用的材料,其屬性是軟件材簧共同作用會導(dǎo)致音圈上下振動,從而帶動音膜振動,聲料庫所無法提供的,在模擬時這些數(shù)據(jù)是否能夠直接輸入？波在空氣中傳播出去,其聲場空間分布受到喇叭輪廓的影〔8結(jié)構(gòu)變形失效研究：結(jié)構(gòu)工作時會有振動〔變響,如圖2所示。形的發(fā)生,需要研究在長時間工作后是否會發(fā)生變形失效問題,導(dǎo)致響應(yīng)頻率的改變,這種研究軟件中是否能夠通過模擬來實現(xiàn)？〔9散熱問題的模擬：在某些稍大型揚(yáng)聲器設(shè)計中,隨著功率的增加,其發(fā)熱量顯著增大,如何模擬熱量對于結(jié)構(gòu)的影響〔如結(jié)構(gòu)熱變形,受熱材料屬性發(fā)生變化,及結(jié)構(gòu)散熱的仿真？..三、COMSOLMultiphyscis微型揚(yáng)聲器設(shè)計由于微型揚(yáng)聲器模擬仿真中所涉及問題的復(fù)雜性,其多物理場耦合的解決是常規(guī)單場分析軟件所無法做到的。在這種情況下,多物理場耦合耦合分析軟件COMSOLMultiphysics能夠發(fā)揮很大的作用。對于多物理場耦合系統(tǒng)的建模和仿真分析,COMSOLMultiphysics提供了一個科學(xué)的軟件分析環(huán)境?？蛇x擇的模塊功能包括電磁學(xué)、聲學(xué)、熱傳遞、微機(jī)電、化學(xué)工程反應(yīng)實驗室以及結(jié)構(gòu)力學(xué)等專業(yè)學(xué)科的專門工具。用戶可以通過求解微型結(jié)構(gòu)的力-電-磁-聲-熱耦合問題,來精確模擬微型揚(yáng)聲器工作中的各種物理現(xiàn)象,進(jìn)一步預(yù)測分析揚(yáng)聲器的工作性能及可靠性。COMSOLMultiphysics集前處理器、求解器和后處理器于一體,在同一個圖形化操作界面中可以完成幾何建模、網(wǎng)格剖分、方程和邊界參數(shù)設(shè)定、求解以及后處理。用戶可以在圖形化界面中構(gòu)建自己的幾何模型,或者直接導(dǎo)入SolidWorks、Inventor和Pro/ENGINEER等CAD軟件構(gòu)建的幾何模型。COMSOL還內(nèi)置了專門的材料庫,支持用戶自定義材料參數(shù)。在設(shè)計中,能夠方便的定義線圈的載荷〔如電流密度等,創(chuàng)建含有永磁體的動圈結(jié)構(gòu),以及很薄的多層復(fù)合音膜。此外,完美匹配層〔PML、預(yù)定義的聲-殼耦合和應(yīng)用阻尼本構(gòu)模型等,使得COMSOLMultiphysics可以自如地圖1動圈揚(yáng)聲器結(jié)構(gòu)示意圖2二維軸對稱建模,加入空氣域和外圍PML考慮無限大聲場模擬仿真分為兩步進(jìn)行?！?電磁分析。使用AC/DC模塊中的磁場〔MagneticFields接口,求解音圈附近的磁場方程,如圖3所示。分析可以獲得驅(qū)動力BL；并且通過計算音圈在一定頻率范圍內(nèi)的時諧場,得到閉合線圈的阻抗Zb?！?聲學(xué)-結(jié)構(gòu)耦合分析。在聲學(xué)-結(jié)構(gòu)耦合分析中使用第一步計算得到的驅(qū)動力BL和阻抗Zb,計算聲壓分布,及同一頻率上的驅(qū)動器總電阻。..CAD/CAM與制造業(yè)信息化·20XX第10期33..oftWare軟件世界..靠近表面。這導(dǎo)致電感及阻抗的電阻部分隨著頻率而變化。圖6顯示的是閉合線圈的電感隨頻率變化的曲線。..圖3音圈附近的磁場分布〔磁場強(qiáng)度的模極片和頂板中的鐵芯是非線性磁性材料,其BH曲線是一組離散的數(shù)據(jù),COMSOLMultiphyscis可以直接導(dǎo)入這些數(shù)據(jù)來定義符合工程實際的材料。圖4清楚地顯示了工程中實際的頂板和極片的磁導(dǎo)率分布。圖4受磁場影響下的頂板和極片的磁導(dǎo)率分布軟件同時可以給出電場的計算結(jié)果,如閉合線圈的阻抗,以及頂板和極片上的感應(yīng)渦旋電流,如圖5所示。圖5極片和頂板分別在52Hz和90Hz時的感應(yīng)電流在較高頻率處,趨膚效應(yīng)〔skineffect使得電流更圖6閉合線圈的電感時頻率的函數(shù)通過聲學(xué)-結(jié)構(gòu)耦合分析,圖7顯示的是3500Hz時聲壓級分布。最小值位于擋板上方45°附近。在較低頻率,聲壓級是平均分布的,但是峰值是在軸向上。圖73500Hz時的聲壓級分布通過進(jìn)一步的掃頻計算,我們可以獲得揚(yáng)聲器的靈敏度,即不同頻率情況下,軸向某點(diǎn)處聲壓級的變化,如圖8所示。首選的工作范圍是在響應(yīng)較平均的區(qū)域,也就是說,通常在100~1500Hz的區(qū)域內(nèi)。圖8揚(yáng)聲器靈敏度,測量距離1米處的軸向聲壓級〔dB得到總電阻,定義為Z=V0/I,如圖9所示。出現(xiàn)在40Hz的峰值是由于恰逢機(jī)械共振；在這個頻率上,阻抗的電阻部分是..34idnovo..欄目主持：侯琳投稿信箱：houlidnovo從電感式到電容式的轉(zhuǎn)換標(biāo)志。在大多數(shù)工作頻率范圍內(nèi)阻抗主要是電阻。在100Hz~1kHz,電阻變化較小。在頻率高于1kHz時,阻抗持續(xù)增加,這里音圈的電感開始扮演重要作用。..圖9揚(yáng)聲器電阻抗是頻率的函數(shù)2.聲-殼耦合分析上例演示了動圈揚(yáng)聲器整個分析的流程。在實際的工業(yè)應(yīng)用中,常常揚(yáng)聲器并不是軸對稱結(jié)構(gòu),因此就必需進(jìn)行三維建模分析。在三維結(jié)構(gòu)中,音膜相比動圈結(jié)構(gòu)而言非常薄,所以整個模型是一個跨尺度分析問題。為了將計算量控制在可接受的范圍內(nèi),同時保證計算的精度,業(yè)內(nèi)通常使用殼結(jié)構(gòu)來描述音膜。所以在實際工業(yè)應(yīng)用中,業(yè)內(nèi)公認(rèn)使用聲-殼耦合來模擬薄膜發(fā)聲,而不是上例所使用的聲-固耦合。也就是把上例中的仿真第二步用本例所描述的方法進(jìn)行分析。下例展示了這項功能。圖10展示了揚(yáng)聲器中的錐體音膜及防塵蓋的模型。其材料設(shè)定為鈦合金,厚度為0.5mm,在模擬仿真中我們使用殼體來描述。在模型外部建立空氣域,計算在空氣域中的聲壓分布。圖10使用殼體來分析音膜振動和發(fā)生假設(shè)電磁分析已計算完畢,電磁驅(qū)動力的大小和頻率已知。音膜結(jié)構(gòu)由于內(nèi)部音圈的受力而振動。我們略去掃頻的過程,僅選擇1kHz和6kHz兩個頻率進(jìn)行頻域分析作為演示,如圖11~圖13所示。圖111kHz和6kHz的錐體瞬時位移〔以z方向的位移來描述圖121kHz和6kHz時空氣域中的聲壓場切面圖,及同零相位相比的位移變形情況圖131kHz和6kHz時空氣域中的聲壓級切面圖,及同零相位相比的位移變形情況..CAD/CAM與制造業(yè)信息化·20XX第10期35..oftWare軟件世界..在圖12和13中可以看到,在6kHz時,在Z軸30°左右存在兩個低值區(qū)域。在同樣的模型上施加同樣的初始位移,可以看到在高頻時整體的聲壓級更高。3.復(fù)合薄膜性能分析有了殼分析仍然不能完美滿足工業(yè)界的需求。在業(yè)內(nèi),音膜是整個揚(yáng)聲器中最重要的部件之一。頂級揚(yáng)聲器設(shè)計,音膜總是最令工程師絞盡腦汁的地方。多層膜復(fù)合,每層膜中各個部分使用不同的厚度和材質(zhì),并且在特定的位置加入各種花紋和凹槽,改善音膜的頻響。在業(yè)內(nèi),各廠商多有自己的CAD環(huán)境進(jìn)行這些音膜的設(shè)計。本例演示了一個Pro/ENGINEER設(shè)計的復(fù)雜雙層復(fù)合薄膜,導(dǎo)入到COMSOLMultiphysics中進(jìn)行仿真分析,如圖14、圖15所示。圖14雙層復(fù)合音膜〔中間藍(lán)色部分為音圈圖15每層膜各個部分往往采用不同的厚度設(shè)計和材質(zhì)設(shè)計COMSOLMultiphysics為這種復(fù)雜膜分析提供了專門的解決方案。軟件中可以設(shè)定任意多層膜復(fù)合,同時分別指定每層膜中各局部的細(xì)部幾何和材料屬性。我們演示了典型的頻域掃略分析,設(shè)定一定的頻率范圍〔一般為200Hz~10kHz,并且在低頻部分設(shè)定更多的采樣點(diǎn)。特別說明：我們使用一臺4核2.66G的臺式計算機(jī),共掃描了從200Hz~10kHz范圍內(nèi)的共137個頻率,運(yùn)算時間為13小時。最終可以給出復(fù)合音膜揚(yáng)聲器的靈敏度曲線〔類似圖8。由于工業(yè)保密原因,本曲線在本方案中略去。4.結(jié)構(gòu)失效分析在工業(yè)界微型揚(yáng)聲器的設(shè)計制造中,一些結(jié)構(gòu)〔如揚(yáng)聲器固定彈片在長時間工作中可能會發(fā)生結(jié)構(gòu)屈服失效的現(xiàn)象。這里我們建立一個如圖16所示模型。圖16彈片屈服分析模型由三部分組成,在接觸面處我們建立接觸對,通過設(shè)定法向接觸系數(shù)來連接模型,如圖17所示。圖17彈片一段固定在結(jié)構(gòu)體中,一端懸空圖18顯示了屈服的情況,圖中繪制的是mises應(yīng)力,紅色部分表示應(yīng)力超限,材料會屈服。同時結(jié)果顯示出,僅在表皮應(yīng)力超限,而不是整體屈服。圖18彈片屈服應(yīng)力分布四、小結(jié)COMSOLMultiphysics是一款大型的高級數(shù)值仿真軟件,廣泛應(yīng)用于各個領(lǐng)域的科學(xué)研究及工程計算。模