隨著電子技術的發(fā)展，電子、電氣設備越來越廣泛地應用于航空飛行器上，惡劣的電磁環(huán)境往往使電子或電氣設備不能正常工作，導致航空通信、控制系統性能的降低。因此，電磁兼容性就成為工程設計中的一個重要課題。只有在整個系統的初期設計時，對系統的電磁兼容性進行預測，對不滿足電磁兼容要求的部分進行優(yōu)化處理，才能較終實現整個系統的高可靠性。

　　近年來，計算電磁學發(fā)展迅速，理論日趨成熟和完善，并廣泛應用在電磁兼容領域。利用計算電磁學理論解決工程中電磁兼容問題，不僅提高了計算的準確性，而且降低了成本，為在研發(fā)初期掌握產品可能出現的電磁兼容問題提供了有效的手段和途徑。因此，仿真計算對于那些利用實驗和測試方法難以解決的電磁兼容問題是一種很好的解決方法。

　　本文介紹了幾種常用的數值算法，并對各種算法進行了分析。針對航空用電設備常見的電磁兼容問題，利用計算電磁學理論進行分析和數值仿真，實現了對用電設備電磁兼容性有效的預測，為系統的初期電磁兼容性設計提供了理論依據。

　　1 電磁場數值方法的介紹和比較

　　計算電磁學是現代電磁場理論、現代數學方法和現代計算機技術相結合而產生的一門新興交叉學科。其主要任務是通過大型計算解決各相關領域中提出的各種極復雜的電磁場問題。

　　當前電磁學中使用較多的數值方法主要有兩類，一類是以電磁場問題的微分方程為基礎的數值方法，如有限元法(FEM)、時域有限差分法(FDTD)等;另一類是以電磁場問題的積分方程為基礎的數值方法，如矩量法(MOM)、多層快速多極子法(MLFMM)等。

　　1.1 有限元法

　　有限元方法是近似求解數理邊值問題的一種數值技術，較早于20世紀40年代提出，在六七十年代被引進到電磁場問題的求解中。該方法的原理是用許多子域來代表整個連續(xù)區(qū)域，在子域中未知函數用帶有未知系數的簡單插值函數來表示，利用里茲變分法或伽略金方法得到一組代數方程，較后通過求解這組方程得到原邊值問題的近似解。原邊值問題可表示為：

　　有限元法可以方便地分析具有復雜幾何結構和非均勻介質材料的電磁問題，因此，這種方法在各種復雜的靜態(tài)場問題、導波問題、電磁輻射和散射問題中得到了廣泛的應用。

1.2 矩量法

　　自從20世紀60年代Harrington提出矩量法基本概念以來，它在理論上日臻完善，并廣泛用于工程之中，特別是在電磁兼容領域，矩量法更顯示出 其獨特的優(yōu)越性。它的思想是將待求的積分或微分問題轉化為一個矩陣方程問題，借助于計算機，求得其數解。很多電磁場問題都歸結為這樣一個算子方程：

　　式中：L為算子;g為已知激勵函數;f為未知響應函數。展開未知函數f為有限個線性無關的已知簡單函數fn之和：

　　式中：an是展開系數;f1，f2，…，fn為展開函數或基函數。將式(8)代入式(7)，再應用算子L的線性，可以得到：

　　選一組線性無關的函數ωm(m=1，2，…，N)，分別與式(9)兩邊作內積。

　　因為m=1，2，…，N，所以得到N個方程，解出f。

　　矩量法就是這樣一種將算子方程轉化為矩陣方程的一種離散方法。

　　矩量法更適合于求解具有表面電流分布的各種幾何體，如計算天線遠、近場輻射場強、方向圖等。它的算法簡單，不需要設置邊界條件，而且對于適當的尺寸，求解速度較快。

　　1.3 時域有限差分法

　　K.S.Yee于1966年提出求解電磁問題的時域有限差分法，其原理非常簡單，即直接將時域Maxwell方程組的兩個旋度方程中關于空間變量和時間變量的偏導數用差商近似，從而轉換為離散網絡節(jié)點上的時域有限差分方程。

　　為了建立差分方程，首先要將求解空間離散化。通常是以一定形式的網格來劃分求解空間，Yee提出了如圖1所示的差分網格單元，其特點是在同一網格中，E和H的各分量在空間取值點交叉放置，使每個坐標面上的的四周由分量環(huán)繞，同時每個場四周由場環(huán)繞。這樣，配置符合Maxwell方程的基本要求，也符合電磁波空間的傳播規(guī)律，使電磁波的時域特性被直接反映出來，直接給出非常豐富的電磁場問題的時域信息。

　　時域有限差分法在天線輻射特性計算、微波電路分析、散射體雷達散射截面等方面有廣泛的應用，對于計算孔縫對屏蔽效能的影響具有優(yōu)越性。

　　1.4 幾種方法的比較

　　矩量法既要面對繁難的積分方程，又要注意基函數的恰當選取;既要耐心處理奇異點，又要巧妙構思快速求解技術。相對而言，實施有限元要容易些，只需要注意基函數選取及稀疏矩陣存儲方式即可。至于時域有限差分就更容易了。因此一般說來矩量法實施較難，有限元次之，時域有限差分法較易。

　　就通用性而論，有限元與時域有限差分相近，都很通用，矩量法則稍差。矩量法通用性的不足從某種程度上說換來了高精度、高效率。雖然原則上說，三種方法精度相當，然而實際計算表明，矩量法精度較高，有限元次之，時域有限差分較差。其原因是矩量法沒有數值色散誤差，其他兩種都有。時域有限差分不僅有數值色散誤差，且模擬復雜幾何形狀的誤差一般也要大于其他兩種數值方法。

　　2 計算電磁學在電磁兼容領域的應用

　　航空用電設備電磁兼容問題類主要有天線間的耦合干擾、場線耦合，機箱的屏蔽性能等。若在用電設備設計初期就能考慮電磁兼容性問題，就能比較容易地滿足飛行器上電磁兼容性的要求，同時也節(jié)約了產品測試成本。因此，利用數值方法對電磁兼容問題進行仿真計算，并對計算結果進行分析，就可以有效地對設備的電磁兼容性進行預測和評估，為設備電磁兼容設計提供可靠的理論依據，并具有重要的實際意義。

　　2.1 天線耦合問題

　　飛行器上的天線形式多樣，工作時要求能夠與其他天線兼容，因此天線的布局設計尤為重要。若要減小天線間的耦合，就要選擇使兩天線間耦合系數盡可能小地點放置天線，但又要考慮天線本身的輻射特性較佳等諸多因素。因此實際天線布局設計是一個綜合性的調整過程，預先計算天線間的耦合系數，對于系統發(fā)揮較大功效并電磁兼容性良好來說，是非常重要的。

　　當天線端口匹配時，天線間的耦合度可以表示為：

　　式中：Pr為接收天線接收到的功率;Pt為發(fā)射天線的輸入功率。

　　算例分析：兩個喇叭天線相對放置，兩天線間放置一塊金屬板，可以有效去除直線上的直接耦合。采用有限元法計算兩個喇叭天線的耦合系數。圖2為喇叭天線的仿真模型。圖3為計算得到天線間的耦合系數。

2.2 開縫箱體屏蔽效能計算

　　一般情況下，航空電子設備都是用金屬箱體來屏蔽外界電磁干擾的。機箱上可能在蓋板、通風散熱孔、電源信號線處存在孔縫隙，電磁能量可通過屏蔽機箱上這些孔縫直接進入電子設備，孔縫耦合作用嚴重影響了機箱的屏蔽性能，降低設備或系統的可靠性。因此需要對機箱的屏蔽效能進行數值仿真，使敏感器件避開場的峰值區(qū)域，提高電子設備的抗干擾能力。

　　屏蔽體的好壞用屏蔽效能來描述，屏蔽效能表現了屏蔽體對電磁波的衰減程度，定義為屏蔽前某點的場強與屏蔽后該點的場強之比。用公式表示為：

　　式中：E0為屏蔽前某點的電場強度;Es為屏蔽后某點的電場強度。

　　算例分析：某一用電設備機箱，其正面中心處有一縫隙，并面對入射波方向。在小孔面積相同的情況下，考慮了圓形孔、正方形孔和兩個尺寸不同的矩形孔