淺談PCB電磁場(chǎng)求解方法及仿真軟件

　　商業(yè)化的射頻EDA軟件于上世紀90年代大量的涌現，EDA是計算電磁學(xué)和數學(xué)分析研究成果計算機化的產(chǎn)物，其集計算電磁學(xué)、數學(xué)分析、虛擬實(shí)驗方 法為一體，通過(guò)仿真的方法可以預期實(shí)驗的結果，得到直接直觀(guān)的數據。“興森科技-安捷倫聯(lián)合實(shí)驗室”經(jīng)常會(huì )接到客戶(hù)咨詢(xún)，如何選擇PCB電磁場(chǎng)仿真軟件的 問(wèn)題。那么，在眾多電磁場(chǎng)EDA軟件中，我們如何“透過(guò)現象看本質(zhì)”，知道每種軟件的優(yōu)缺點(diǎn)呢?需要了解此問(wèn)題，首先得從最最基本的求解器維度說(shuō)起。

　　本文旨在工程描述一些電磁場(chǎng)求解器基本概念和市場(chǎng)主流PCB仿真EDA軟件，更為深入的學(xué)習可以參考計算電磁學(xué)相關(guān)資料。

　　電路算法

　　談 到電磁場(chǎng)的算法，不要把場(chǎng)的算法和路的方法搞混，當然也有場(chǎng)路結合的方法。電路算法主要針對線(xiàn)性無(wú)源集總元件和非線(xiàn)性有源器件組成的網(wǎng)絡(luò )，采用頻域 SPICE和純瞬態(tài)電路方程方法進(jìn)行仿真。這類(lèi)仿真的特性是無(wú)需三維實(shí)體模型、線(xiàn)性和非線(xiàn)性器件時(shí)域或頻域模型(SPICE和IBIS等)、仿真速度快、 電壓電流的時(shí)域信號和頻譜為初級求解量。電路仿真簡(jiǎn)稱(chēng)路仿真，主要用于端口間特性的仿真，就是說(shuō)當端口內的電磁場(chǎng)對網(wǎng)絡(luò )外其他部分沒(méi)有影響或者影響可以忽 略時(shí)，則可以采用路仿真;采用路仿真的必要條件是電路的物理尺寸遠小于波長(cháng)。換言之，當電路板的尺寸可以和電路上最高頻率所對應的波長(cháng)相比擬時(shí)，則必須使 用電磁場(chǎng)理論對該電路板進(jìn)行分析。舉例說(shuō)明，一塊PCB尺寸為10*10cm，工作的最高頻率是3GHz，3GHz對應的真空波長(cháng)是10cm，此時(shí)PCB 的尺寸也是10cm，則我們必須使用電磁場(chǎng)理論對此板進(jìn)行分析，否則誤差將很大，而無(wú)法接受。一般工程上，PCB的尺寸是工作波長(cháng)的1/10時(shí)，就需要采 用電磁場(chǎng)理論來(lái)分析了。對于上面的那塊板子，當板上有300MHz的信號時(shí)，就需要場(chǎng)理論來(lái)析了。

　　電磁場(chǎng)求解器分類(lèi)

　　電子產(chǎn)品設計中，對于不同的結構和要求，可能會(huì )用到不同的電磁場(chǎng)求解器。電磁場(chǎng)求解器(Field Solver)以維度來(lái)分：2D、2.5D、3D;逼近類(lèi)型來(lái)分：靜態(tài)、準靜態(tài)、TEM波和全波。

　　1、準靜電磁算法

　　它 需要三維結構模型。所謂“準靜”就是指系統一定支持靜電場(chǎng)和穩恒電流存在，表現為靜電場(chǎng)和靜磁場(chǎng)的場(chǎng)型，更精確地講，磁通變化率或位移電流很小，故在麥克 斯韋方程組中分別可以忽略B和D對時(shí)間的偏導項，對應的麥克斯韋方程分別被稱(chēng)之為準靜電和準靜磁。由此推導出的算法就被稱(chēng)之為準靜電算法和準靜磁算法。這 類(lèi)算法主要用于工頻或低頻電力系統或電機設備中的EMC仿真。如：變流器母線(xiàn)與機柜間分布參數的提取便可采用準靜電磁算法完成。對于高壓絕緣裝置顯然可采 用準靜電近似，而大電流設備，如變流器、電機、變壓器等，采用準靜磁算法是較可取的。

　　2、全波電磁算法

　　簡(jiǎn) 單地講就是求解麥克斯韋方程完整形式的算法。全波算法又分時(shí)域和頻域算法。有限差分法(FD)、有限積分法(FI)、傳輸線(xiàn)矩陣法(TLM)、有限元法 (FEM)、邊界元法(BEM)、矩量法(MoM)和多層快速多極子法(MLFMM)均屬于全波算法。所有的全波算法均需要對仿真區域進(jìn)行體網(wǎng)格或面網(wǎng)格 分割。前三種方法(FD、FI和TLM法)主要是時(shí)域顯式算法，且稀疏矩陣，仿真時(shí)間與內存均正比于網(wǎng)格數一次方;后四種方法(FEM、BEM、MoM和 MLFMM)均為頻域隱式算法。FEM也為稀疏矩陣，仿真時(shí)間和內存正比于網(wǎng)格數的平方;而B(niǎo)EM和MoM由于是密集矩陣，所以時(shí)間與內存正比是網(wǎng)格數的 三次方。FD、FI、TLM和FEM適用于任意結構任意介質(zhì)，BEM和MoM適用于任意結構但須均勻非旋介質(zhì)分布，而MLFMM則主要適用于金屬凸結構， 盡管MLFMM具有超線(xiàn)性的網(wǎng)格收斂性，即大家熟知的NlogN計算量。