1 引言

在電磁領域里多層介質結構有著廣泛的應用，例如具有夾芯層的雷達罩壁結構、多層吸波屏等。對于這類三維復雜結構介質體的電磁仿真分析，VIE－MoM具有很高的計算精度，它與高頻方法的結合可以高效地分析結構復雜的電大尺寸的雷達罩的電磁性能。矩量法在計算時需要對積分方程中未知函數(shù)的定義域進行離散，即需要對所分析的目標形體（即未知函數(shù)定義域）的幾何模型進行適當?shù)木W(wǎng)格離散。所建立的幾何模型是否較好地符合實際目標形體，離散后的模型在幾何上是否具有對模型的良好擬合效果，將影響到最終矩量法計算結果的精確性。因此，實現(xiàn)有效的三維形體建模與網(wǎng)格離散是成功實施分析的關鍵問題之一。

現(xiàn)有的商用建模軟件有較強的適應性，對于比較規(guī)則的形體，都能快速地建立起精確的幾何模型，并且能夠實現(xiàn)對幾何模型完成典型離散單元的快速網(wǎng)格剖分，獲得滿足電磁特性分析要求的網(wǎng)格參數(shù)文件。但是，對于雷達罩中常見的多層介質的復雜結構，直接應用這些商用軟件實體建模功能得到的實體，在四面體網(wǎng)格離散之后很難解決多層結構之間網(wǎng)格匹配等問題，進而會影響到矩量法分析結果的準確性。因此，本文將在商用軟件較強的復雜曲面建模功能的基礎上，針對多夾層結構的天線罩，拓展其在多層結構之間網(wǎng)格匹配方面的建模功能，從而為多層結構介質體電磁性能的分析提供高質量的網(wǎng)格數(shù)據(jù)文件。

2 實體建模和四面體網(wǎng)格離散

現(xiàn)有CAD軟件在建模功能和網(wǎng)格離散功能方面均很強大，可以使用它們來實現(xiàn)對復雜形體結構的幾何建模和快速的自動網(wǎng)格剖分。

以MSC.Patran為例，所分析的三維介質區(qū)域用實體模型（solid）來實現(xiàn)，對實體可以進行以以四面體為單元的自動剖分。在實體建模時，為便于VIE－MoM的分析與計算，將介質參數(shù)不同的區(qū)域用不同的實體來表示，這樣可以在網(wǎng)格離散后對各個區(qū)域的離散單元分配其相應的介電參數(shù)。

四面體能夠較好地離散擬合任意復雜形狀實體，本文采用四面體單元(Tet)對實體模型進行離散，全局尺度設置為工作波長的1/10左右。自動網(wǎng)格剖分軟件會根據(jù)形體邊界的曲率和尺寸因素對剖分尺度進行自適應調(diào)整，以達到較好的形體擬合效果。但是CAD軟件靈活性比較差，尤其是有多個實體組合的情況下，網(wǎng)格剖分后往往在各個實體的交界面處會出現(xiàn)不匹配情況，即不同實體在分界面上形成不同的三角網(wǎng)格分布。采用低頻方法分析多層介質體電磁問題時時，無論是矢量有限元法還是基于體積分方程的矩量法，網(wǎng)格的不連續(xù)性會造成計算結果的不準確性，因此必須要保證在兩個不同實體交界面上的網(wǎng)格匹配性。

3 四面體網(wǎng)格匹配性的檢查

對實體進行四面體網(wǎng)格剖分后，四面體網(wǎng)格分布在實體的內(nèi)部，其在分界面處的網(wǎng)格情況是不易被觀察到的，如圖1所示。

圖1 實體模型及其四面體網(wǎng)格剖分效果圖

為能夠有效檢查不同實體的四面體網(wǎng)格在其交界面處的匹配情況，本文通過對網(wǎng)格數(shù)據(jù)文件進行處理，將不同實體分界面處的網(wǎng)格匹配情況提取出來，使其可視化，便于檢查和修改。圖2給出的是圖1所示的模型在不同網(wǎng)格剖分情況下的交界面上網(wǎng)格匹配情況。

（a）完全匹配 （b）部分匹配

圖2 網(wǎng)格剖分后提取到分界面處的信息

從圖2可以看到，如果在交界面是匹配的，則得到的交界面一個由三角面元組成的完整方形平面。反之，當交界面不完整時，表示在其缺損的區(qū)域出現(xiàn)了網(wǎng)格不匹配的情況。造成這種不匹配的情況主要有兩個原因造成的：一種是在交界面兩側的實體在交界面分別形成了不一致的結點體系，如圖3(a)所示：實體1在分界面上產(chǎn)生的結點分別為1,2,3,7,8,9,10；實體2在分界面上產(chǎn)生的結點為4,5,6,7,8,9,10。另一種情況是雖然結點體系是一致，但是其組成三角網(wǎng)格的組合是不同的，如圖3(b)所示，在實體1中的四面體（1356）與實體2中的四面體（2456）在交界面上就形成了兩個交叉的三角面元，從而造成了兩側實體在分界面上由結點3、4、5和6圍成的區(qū)域上內(nèi)是不匹配的。前一種情況當兩個實體結構不相似時，或者兩邊的剖分尺度不同時，比較容易產(chǎn)生，而后一種情況是比較難以被預先估計的。

a）結點體系不一致的情況 b）結點組合不一致的情況

圖3 分界面處網(wǎng)格不匹配的成因