直接用电密J计算导体三维涡流场的FEBE耦合法.pdfVIP

中目计季(E辑) 第32卷第6期 SCIENCEINCHINA(SeriesE) 2002年12月 直接用电密t，计算导体三维涡流场的 FE．BE耦合法 mf n 刘志珍 王衍章 f山东大学电气工程学院．济南250061) (山东大学科技处．济南250061) 贾智平 孙英明 {山东大学计算机科学与技术学院，济南25006 (济南大学计算中心，济南250002 摘要 提出一种计算导体三维涡流场的FE．BE耦台法．在涡流区中直接采用涡流电 密J为求解量的FEM，其公式采用八节点长方体单元进行离散，在非涡流区中采用磁场 强度日为求解量的BEM，其公式采用矩彤单元离散．根据交界面条件，耦合得到一个求 解变量仅为，的代数方程组．该方法的优点就是所求量无需微分，求解变量少，祸合客易， 适用于多连域问题．对两个模型进行了计算，计算结果与测试结果吻台较好． 关键词 导体 涡流电密 磁场强度 三堆涡流场 有限元法 边界元法 耦台 电气工程中大型涡流场求解问题仍是电磁场数值计算的难题之一，研究行之有效的计算 方法成为电磁场数值计算人员追求的目标．目前，在这个领域应用较多且比较完善的计算方法 主要是有限元法(FEM)和边界元法(BEM)，它们是电磁场数值计算的两个非常有用的工 具FEM容易处理封闭区域、非线性、复杂介质问题且系数阵稀疏、正定、对称，但处理开 域问题较困难．BEM易于解决开域、线性问题，尽管其代数方程组较FEM要少，但其系数阵 为不对称满阵．这两种方法对于计算简单的涡流场问题非常有用，但对于电气工程中复杂的涡 流场问题，仅用一种方法往往不是很奏效，所以混合法在涡流场计算中越来越显示出优势 最有优势且应用与研究最多的一种，FEM一般用于涡流区，BEM用于非涡流区，所求解的状态 量也各式各样，如H．似A一铲皿T-．Q-H和A’．口等组合”“．这些方法的优势显而易见，但缺陷 也同样明显，如状态量较多．耗费计算资源：求解涡流分布时需要微分，影响精度：耦合方 法复杂等．这些成为FE．BE耦合法急需解决的问题． 基于FE．BE耦合法存在的问题，本文提出一种计算导体三维涡流场的FE-BE耦合法．在 涡流区直接采用涡流电密．，为求解量的FEM，非涡流医采用磁场强度日为求解量的BEM，根 据交界面条件，耦台得到一个求解变量仅为_，的代数方程组．该方法使用场量求解，无需考虑 唯一性及规范条件问题：涡流区用电流密度_，求解，无需微分，保证了较高的精度；求解变量 较少，耦合容易：适用于多连域问题． 20024)4．09收稿，2002-08-05收修改稿 第6期 刘志珍等直接用电密J计算导体三维涡流场的FE—BE耦合法 频情况下，磁场是准静态的，不考虑位移电流的影响． 图l涡流场模型示意图 (ii)导体各向同性．线性均匀．(jii)温度的影响不计．(iv)各量均按正弦变化， 1．1涡流区皿 在涡流区垃，根据Maxwell微分方程组，有 V×V‘J+jaWoa,J=0， (1) 甲xJ=-jaWo暝,H． (2) 利用加权余量法，取矢权函数砒对(1)式进行加权积分，得 (3) 『nw‘(vxVxj+jm#otYeJ)dO=0- 对(3)式利用散度定理及矢量恒等式进行变换，可得