矢量标量混合位有限元结合修正的IC预条件技术在电磁仿真中地研究.pdfVIP

矢量标量混合位有限元结合修正的IC预条件技术 在电磁仿真中的研究 芮平亮朱剑平学伟丁跃华陈如山 (南京理工大学电子工程与光电技术学院) 摘要：基于电场强度的切向矢量有限元方法(￡一TVFIY．{)被广泛用于电磁问题的有限元分析中。然 而．随着操作频率的降低或者网格剖分的细化，由此产生的有限元刚阵的性态会逐渐变差，从而不 利于对有限元刚阵的迭代求解。本文中采用矢量标量混合位有限元(AV—FEll)方法从本质上弥补 E—TVFEM构造上的缺陷．同时应用一种修正的Ic预条件技术来进一步改善剐阵的性态，提高迭代 法的收敛速度。数值结果表明．在对一些典型的电磁结构进行分析时，显示r新方法的优越性。 关键词：切向矢量有限元；矢量标量混合位有限元；蟹芷蛩Ic预条件技术 一L 1 引言 有限元工具对于建模分析任意形状的复杂几何体和任意特性的材料媒质具有强大的优势， 因此在电磁领域中得到了广泛的应用。其中最常用的是基于电场强度构造的切向矢量有限元 (琶-TVFEM)．简称为E场有限元。它克服了基于节点有限元的几个严重的缺陷：(1)避免了 非物理的或伪解的缺陷：(2)方便了材料表面和导体表面边界条件的强加：(3)克服了处理导 体和介质边缘及角的困难性，避免了与这些结构有关的场的奇异性。虽然暮场有限元有限元 在模型分析方面取得了极大的成功，但是在分析导行波的传播上还存在这一些问题。如文献【l】 中所指出，君场有限元在这种情况下的迭代法求解的收敛速度很慢。此外，这种构造也不适 用于低频的情况。舌场有限元的这两种缺陷归根揭底都是由其构造上的一种缺陷所产生的， 即没有考虑一种不依赖频率的规范条件的强加。幸运的是．有很多的矢量标量位的方法 (万一FEM)，简称刀有限元，来解决这个问题。【2】中提出了一种树型的规范条件，但是由 此产生的刚阵的条件数在很大的程度上取决于具体的树型的选择。我们这里考虑的是引入一种 拉格朗日乘子【3】，来实现非标准的规范条件的强加，从而达到优化有限元刚阵性态的目的． -同时，为了进一步加速迭代算法的收敛速度，前人提出过很多的预条件方法．其中最著名的 是不完全Lu分解预条件技术．将之应用于有限元产生的对称阵时，则称之为Ic预条件技术．由 于LU分解本身的不稳定性。给IC预条件技术的应用造成了很大的困难．这里我们采用了一种修 正补偿的办法，来弥补不完全乔里斯基分解预条件技术的缺陷，这种修正过的Ic预条件我们称 之为MIC．文中将万有限元构造与MIC预条件技术相结合，应用到共扼梯度法(CG)迭代求解电 磁问题．数值结构表明，对于一些典型的结构，我们的方法的效果显著优越于罾场有限元预条 件的迭代算法． 2耖有限元构造： 考虑如下的边值问题： inQ(i) V×[／J，rV×西一《啡】雷=6 34 onF￡。其中露E是强加的表面电流，FH}[1 及边界条件耳×i=露￡onFⅣ和露×h=6 【1，分别指的是理想的磁壁和电壁。引入磁矢量位A和电位V使得 D v×j=一Jr／0Lur]万，fly．[勺】j=一y 自由空间中的特征阻抗。系数∥，y的取值决定了所采用的规范条件的形式，从而对应不同的 j秒有限元。这里我们取∥=ko，y=1．对应于非标准的规范条件的强加。将上式带入由(1) 中定义的边值问题，我们可以得到如下新的边值阀题： V×【肼】…1V×i_{卸】(簖j+盘Ⅳy)=6【 (2) V·[erl(flj+r2vv)--ol 满足边界丢‘‘肋寸‘?j皿五=露￡。nrⅣ和j篙}。nl采脯限元分析批，定义 蕞(【印】矿)=0 的边值问题，我们可以得到代数方程： ㈤