一、通用面积分方程构造
为了建立任意多种材料混合结构通用积分方程,先对求解问题作如下约定:
- 模型中的材料按顺序编号,其中0号默认为_PEC_区域,1号为背景空间材料(通常为空气),其它介质材料从2号开始依次排列;
- 每种材料定义为一个区域,任意两种材料的分界面定义为一个边界,如纯金属结构的金属与空气分界面,纯介质结构的介质与空气分界面,介质与介质相交结构的分界面,介质包裹金属结构中的金属介质分界面;
- 金属与任意材料的分界面上只存在等效电流,没有等效磁流;
- 波端口面设置为金属面,并在其上强制放置等效磁流。
1.1、等效原理构造积分方程
如图1.1所示,假设背景空间中存在金属、介质i,介质j,波端口,线结构,且线和金属面相连接。背景空间中存在均匀平面波。
图1.1、任意材料线、面混合结构等效原理
定义金属与背景空间边界为分界面1,介质i与背景空间边界为分界面m,介质j与背景空间边界为分界面n。金属表面放置等效电流密度;;线结构放置等效电流I;线面连接处等效电流密度;波端口面放置等效磁流密度;边界面m外侧放置等效电流、磁流密度,内侧放置等效电流、磁流密度;边界面n外侧放置等效电流、磁流密度。
所有等效源在区域1背景空间中产生的辐射场与外加场总场可表示为:
第i区域内等效源产生的辐射场可表示为:
第j区域内等效源产生的辐射场可表示为:
以下在不同分界面上建立通用积分方程,
1、线结构表面EFIE:
2、金属分界面 1表面上EFIE:
3、连接域上EFIE:
4、介质分界面m表面PMCHW:
其中积分算子简记为:,, 。
5、介质分界面n表面PMCHW:
6、波端口面MFIE:
算法参考"Higher Order MoM Analysis of Traveling-Wave Waveguide Antennas with Matched Waveports," IEEE Transactions onAntennas and Propagation, vol. 63. no. 8.2015
选取合适检验函数在各边界上检验积分方程,将积分方程转化为矩阵方程:
1、线结构表面检验:
2、金属分界面1表面上EFIE:
3、连接域检验:
4、介质分界面m表面PMCHW:
5、介质分界面n表面PMCHW:
6、波端口面MFIE:
1.2、矩阵元素块分类
由矩阵方程可知,通用积分方程矩阵元素主要由线、连接域、金属面、介质面、波端口之间的相互作用构成。可将矩阵元素拆分为以下16个部分分别计算。
A1、线自作用
A2、线、金属面互作用
A3、线、连接域互作用
A4、线、介质面互作用
A5、线、波端口互作用
B1、金属面自作用
B2、金属面、连接域互作用
B3、金属、介质面互作用
B4、金属、波端口面互作用
C1、连接域自作用
C2、连接域、介质面互作用
C3、连接域、波端口面互作用
D1、介质面自作用
D2、不同介质面互作用
D3、介质面、波端口互作用
E1、波端口自作用
以上16个子矩阵块主要由,,,几个算子构成。其中下标i表示材料区域,即当前积分算子作用区域的材料编号;下标m表示分界面编号,表示第m个分界面上的电流或磁流。
1.3、积分方程普适性推广
上述通用积分推导过程建立在各金属、介质结构分离状态下,但该方程依然适用于任意结构。前面推导过程中在介质分界面外侧放置等效电磁流J,内侧放置等效电磁流密度。根据等效原理只需介质面内外侧的电磁流反号,至于哪一侧为正则是相对的。为了便于理解积分方程的通用性,作如图2假设:背景空间1中存在两个介质材料2、3,其中材料2包围材料3。背景空间和材料2分界面记为分界面m,材料2和材料3分界面记为分界面n。分界面m外侧放置等效电磁流、内侧放置等效电磁流,分界面n外侧放置等效电磁流、内侧放置等效电磁流。根据等效原理,,,。
图1.2、包围结构等效原理
首先两个分界面是否产生作用需要根据分界面两侧材料属性进行判断,如图1.2中,分界面m和分界面n材料属性中同时存在2号材料,则这两个分界面产生互作用,且作用算子都是在2号材料区域。在前面推导的积分方程中,阻抗元素公式中对应的源分布函数和检验函数分别是和。至于这两组电磁流的符号完全根据分界面定义时的参考值,确定。也就是说,在前述通用积分方程中,只要保证源函数和检验函数携带自身的符号,方程就可以适用于任意结构。