有限元方法求解金屬襯底介質(zhì)片對平面波的反射問題

徐源，孫一月

(1.中國傳媒大學(xué) 理學(xué)院，北京100024；2.中國人民武裝警察部隊工程大學(xué) 理學(xué)院，陜西 西安710086)

有限元方法求解金屬襯底介質(zhì)片對平面波的反射問題

徐源1，孫一月2

(1.中國傳媒大學(xué) 理學(xué)院，北京100024；2.中國人民武裝警察部隊工程大學(xué) 理學(xué)院，陜西 西安710086)

任意平面波可以被分解為只有Z軸向電場的Ez-極化平面波和只有Z軸向磁場的Hz-極化平面波。我們需要研究的是金屬襯底介質(zhì)片對這兩種波的反射能力。根據(jù)標(biāo)量的赫爾姆霍茲方程分別得到相應(yīng)的極化方程，在方程的邊界上采用不同邊界條件。最后我們通過采用有限元的思想來分析問題，通過Matlab語言實現(xiàn)有限元方法的過程來解決這個問題。文章的最后給出了兩個數(shù)值實驗來驗證真解和數(shù)值解之間的誤差關(guān)系保證了算法的正確性和時效性。

極化平面波；有限元；Matlab

1 問題的描述

在圖1中，一束均勻的入射波在非均質(zhì)介質(zhì)片面上發(fā)生了反射。介質(zhì)片的寬度為L，相對電導(dǎo)率為εr，相對的磁導(dǎo)率為μr都視為x的方程。周圍的自由空間介質(zhì)中εr=μr=1，接下來我們需要考慮的是介質(zhì)片的反射能力。

眾所周知任意平面波可以被分解為一個只有Z軸向電場的Ez-極化平面波和一個只有Z軸向磁場的Hz-極化平面波。因此我們可以充分考慮這兩種波的案例。對于Ez-極化平面波例子，其入射波的表達式可以是：

(1)

圖1 金屬介質(zhì)片平面波的反射圖

E0為一常量，指的是入射場的大小，θ代表圖1中的入射角度。顯而易見為了滿足在與x坐標(biāo)軸相垂直交界面處場的連續(xù)性條件，那么總場必須要有一個共同的系數(shù)e-jk0ysinθ。根據(jù)標(biāo)量的赫爾姆霍茲方程使得電場Ez滿足：

(2)

對于邊界條件我們采用的是齊次的狄利克雷邊界條件

Ez|x=0=0

(3)

相似地，對于極化波，入射波可以表示為：

(4)

總的磁場滿足標(biāo)量的赫爾姆霍茲方程：

(5)

對應(yīng)的紐曼邊界條件為：

(6)

2 解的分析

此問題的解決第一步需要把介質(zhì)片分成很多的薄層，定義為M層(圖2)每一層的相對電導(dǎo)率εr和相對磁導(dǎo)率可以近似地看作為常值，用εrm和μrm其中(m=1，2，3，…，M)。因此關(guān)于方程2在第m層的解可以表示為：

Ezm=(Amejkxmx+Bme-jkxmx)e-jk0ysinθ

(7)

圖2

我們強制Ez和Hy在兩層之間的交界面處是連續(xù)的，指出第m層和第(m+1)層的時候，我們有：

(8)

這里的xm+1指的是交界面處的位置，Rm=Bm/Am，并且有：

(9)

由于邊界條件式3，一旦有A1+B1=0成立，那么R1=-1。根據(jù)式8我們又可以計算出R2，R3等等。最后我們就可以計算出我們感興趣的反射系數(shù)RM+1，我們用|RM+1|2來表示介質(zhì)片反射能力的百分比。相類似的過程可以直接用在磁場的Hz極化波上面：

(10)

由邊界條件式6，一旦有A1-B1=0成立，那么R1=1。根據(jù)式8我們又可以計算出R2，R3等等。我們就可以計算出我們感興趣的Hz極化波反射系數(shù)RM+1。

3 有限元的解

首先，讓我們用有限元的方法來解決這個問題。首先確定問題的求解區(qū)域。顯然問題的求解區(qū)域是半無窮的(0≤x≤∞)，遺憾是有限元方法不能直接的用于這種半有界的求解域。于是我們需要適當(dāng)?shù)目s小下求解域，我們可以引進適當(dāng)?shù)娜斯み吔鐥l件。為了滿足解的有效性，人工的邊界和求解域直接應(yīng)該足夠的接近。所以針對此問題最可能的選擇是令x=L，這里的L指的是板的厚度。

接下來要做的是在此邊界上找一個適當(dāng)?shù)倪吔鐥l件。需要指出的是在板的外部，總場可以表示為入射場和反射場的疊加。因此對于Ez極化波我們有：

Ez(x)=E0ejk0xcosθ+RE0e-jk0xcosθxgt;L

(11)

這里的R代表反射系數(shù)以及公共系數(shù)e-jk0xcosθ已經(jīng)表示過了，對Ez進行對x求偏導(dǎo)數(shù)，我們有：

(12)

也可以表示為：

(13)

或者

(14)

此條件適用于在μr=1處介質(zhì)片外表面的邊界上。如果希望將邊界放在介質(zhì)片的內(nèi)表面，那么可以用到式14的邊界條件需要指出的是：

這是針對Ez和Hz連續(xù)邊界條件的結(jié)果，因此我們不難發(fā)現(xiàn)：

(15)

作為邊界條件僅僅在介質(zhì)片的內(nèi)表面是可行的。顯然式13和式14是等價的，它們的選擇并不會影響方程的解。相類似地針對極化波我們同樣可以找到適當(dāng)?shù)倪吔鐥l件強加在Hz如：

(16)

或

(17)

對于式2和式5的不同，我們分別給出了在x=0處的邊界條件和在x=L處推導(dǎo)出的邊界條件。我們就可以繼續(xù)運用有限元方法來求解問題的數(shù)值解。我們首先是把介質(zhì)片或者是區(qū)域(0，L)分成了一定數(shù)量的薄層，再次將M層每一層都看作是一個單元。因此我們就可以通過編寫有限元方法的計算機語言可以高效地解決這個問題。一旦場在x=L處被解決了，那么方程11的反射系數(shù)也可以得到，或者更具體地對于Ez極化波的反射系數(shù)來說：

Hz極化波來說也是類似，就是把Ez用Hz來代替，E0用H0來代替。

4 數(shù)值實驗

為了保證有限元解的有效性，我們令電導(dǎo)率εr=4+(2-0.1j)(1-x/L)2，相對應(yīng)的磁導(dǎo)率μr=2-0.1j。對于介質(zhì)片的厚度采取五倍于自由空間的波長。實驗的結(jié)果已經(jīng)展示在圖3。采用50單元(51個節(jié)點)和100單元(101個節(jié)點)的有限元的數(shù)值解和真解的比較。正如我們所看到的節(jié)點的數(shù)量越多，數(shù)值解就越接近真解。

5 實驗分析

通過Matlab語言實現(xiàn)的有限元算法來求解金屬襯底介質(zhì)板對平面波的反射問題，我們不難發(fā)現(xiàn)，隨著求解單元數(shù)量的增加，數(shù)值解與真實解之間的誤差會越來越小。這達到了我們預(yù)先的估計和目的，用有限元的方法來解決偏微分方程問題的時效性和準確性優(yōu)勢明顯。根據(jù)完整的有限元的求解過程構(gòu)編寫程序語言是本文的最終目的。

(a)電場極化算例 (b)磁場極化算例圖3 金屬襯底介質(zhì)片的反射系數(shù)εr=4+(2-0.1j)(1-x/L)2，μr=2-0.1j，L=5λ

[1]Harris，Edward G.Introduction to modern theoretical physics[M].V I，1975.

[2]Jianming Jin.The finite element method in electromagnetics[M].Canada：John Wiley amp; Sonsinc，2002：83-90.

[3]嚴登俊，黃學(xué)良，胡敏強，勾磊.有限元網(wǎng)格生成技術(shù)分析[J].微特電機，1999(01).

[4]傅為農(nóng)，江建中.關(guān)于電磁場數(shù)值分析的若干認識[J].微特電機，1998(05).

[5]劉春太，楊曉東，申長雨，陳靜波.任意平面區(qū)域的變尺寸有限元網(wǎng)格劃分[J].計算力學(xué)學(xué)報，2000(01).

[6]Gerald Recktenwald.數(shù)值方法和MATLAB實現(xiàn)與應(yīng)用[M].北京：機械工業(yè)出社，2004.

[7]Snyder A W，Love J D.光波導(dǎo)理論[M].北京：人民郵電出版社，1991.

(責(zé)任編輯：宋金寶)

FiniteElementMethodtoSolvePlane-waveReflectionbyaMetal-backedDielectricSlab

XU Yuan1，SUN Yi-yue2

(1.School of Science，Communication University of China，Beijing 100024；2.College of Science，Engineering University of the Chinese People’s Armed Police Force，Xi’an 710086)

It is well known that any plane wave can be decomposed into an Ez-polarized plane wave having only a z-compoment for the electric field and a Hz-polarized plane wave having only a z-compoment for the magnetic field，we need to find the power reflected by the slab for two polarization cases only.With observation，the scalar Helmholtz equation governing the functions，the boundary condition for equations is different.We consider the finite element solution to this problem and write Matlab code to solve the problem. Finally this dissertation gives the running result of this system and analyzes the accuracy of the algorithm and its stabilization.

polarization plane wave；finite element solution；Matlab

O24

A

1673-4793(2017)06-0032-04

2017-07-06

徐源(1990-)，男(漢族)，山東臨沂人，中國傳媒大學(xué)碩士研究生.E-mail：xuyuan@cuc.edu.cn