时域电场积分方程的精确求解.pdfVIP

￡l￡ 婴￡1 2立Q垒Z盛盏 时域电场积分方程的精确求解。 赵延文武胜波聂在平徐建华 (电子科技大学电子工程学院。网jn成都610054) 摘要：用于求解堕堑麴盆直提(TDIE)的时间步进(MOT)算法的后时不稳定性严 重阻碍了TDIE法的广泛应用和发展。利用参数坐标和广义Duffy坐标变换将时域电场积分 方程(TDEFIE)的奇异性积分转换成非奇异性积分．井据时间基函数定义域将积分转换成 能快速精确计算的分区域积分，人幅提高了利_I}j 时稳定性。该方法适用于任意时间基函数并可推广‘到高阶曲面单元和高阶空间基函数情形。 1．引言 近年来，利用时域方法模拟电磁现象越来越受到人们的重视。在分析均匀无限大介质 中的理想导体或均匀介质体的表面时域散射问题时，采用TDIE法比采用时域微分方程 地阻碍了TDIE法的广泛应用和发展。有许多学者对TDIE的时间后期稳定性问题作了广泛 深入的研究：例如，利用混合场积分方程(CFIE)法来克服内谐振所导致的不稳定性IlI： 采用时问平均等滤波方法来克服高频成份的影响12I以及通过避免直流成份在时问迭代过程 的积累来提高时间后期稳定性13l等等。实际上，MOT算法出现后时不稳定性主要是由于在 离散TDIE时采用了不精确的数值方法和不恰当的解析近似所导致的．精确计算TDEFIE的 时域阻抗矩阵元素可以大幅度提高利用MOT算法求解TDEFIE的精度和后时稳定性I．J。 2．堕缝垒堑墼坌复堡及其时间步进算法 假设s为位于均匀介质中的理想导体的表面，则瞬态入射电场￡“(，，f)将在理想导体表 面s上产生感应电流．，(，r)．利用理想导体表面上切向电场为零的边界条件，可得感应电流 1(r，f)所满足的时域电场积分方程为 吼归曙昙』掣册石IVj严警以0es㈩ 其中。下标“tan”表示切向分量，推迟时间f=t一叫c，R=I，一，1为场点，和源点，，间的 距离，s和Ⅳ分别为散射体周围介质的介电常数和磁导率。c为电磁波在其中的传播速度。 利用矩量法求解方程(1)，可以求得电流J(r，，)。为此，将感应电流．，(，，f)_}}j空间基函 数^(，)(=I，2，…，M)和时间基函数巧(f)(，=I，2，…，Ⅳ，)展开为 ⅣⅣ J(r，f)=艺艺C正(，)‘(r) (2) ^=lI=l 其中以为与时空基函数一(，)巧(f)相关的待求加权系数。在本文中，选取时间基函数 Z(f)=T(t一，出)(At为时间步长)为二角形函数、空间基函数正(，)为RWG基函数吼 将展开式(2)代入方程(1)，井在空间上采州伽略金方法、在时间上采_l}J点匹配方法(时 间匹配点为kAt)可以得剑F列矩阵方稃， 田家自然科学幕会资助项目(堆会o， ●-I 执·，I=垆一∑鸭．·，^一r (3) ，．。I 其中 帅∽2刨剖。一舻避地…+ ㈣ 去川。dIP-R／c川】监等型舔，舔 去jnP瑚，J当半拶舔 3．TDEFIE的奇异性积分的精确快速计算 计算时域阻抗矩阵元素的积分(5)为二重面积分。式(s)中对场点，的积分可采用标 准的高斯积分法。若在给定场点，o之后，积分(5)中对源单元的面积分可表示成 t=f矗’7“‘，R)F(r)IR (6) r 其中r表示源点，’所在的平面三角形单元；F表示非奇异性的任意标量或矢量函数，而 7(f．，R)表示时间基函数的时间导数或积分：R=Ire一，1为场点％和源点r’之间的距离。 分当场点乇和源点，，位于同一个三角形单元上时，积分(6)为奇异性积分。为计算奇 异性积分(6)，可以将平面三角形单元，，上的积分，