048_用有限元方法计算电准静态场

2006年用户年会论文

用有限元方法计算电准静态场

石磊，金大志，唐平瀛

（中国工程物理研究院电子工程研究所，四川绵阳，621900）

[ 摘要 ] 本文介绍了电准静态场EQS(Electric Quasi-Static)的控制方程，用通用的大型有限元软件ANSYS计算电准静态场的方法。通过典型的时谐场和暂态场例题，显示理论分析和数值计算

结果完全符合，可以用到实际的工程问题。

[ 关键词]电准静态，有限元 ,ANSYS

Computation of Quasi-Static Electric Fields with Fea

SHI Lei, JIN Da-zhi,TANG Ping-ying

(Institute of Electronic Engineering, China Academy of Engineering Physics,

Mianyang 621900, China)

[ Abstract ] The controlled equation of Electric Quasi-Static is described. An application of software ANSYS with FEA to the computation of quasi-static electric fields for systems with perfectly

conducting and lossy components is presented. This method shows good agreement

between numerical and analytical calculations for benchmark model, and has been

effectively applied to partical problems.

[ Keyword ] EQS,FEA,ANSYS

1前言

电场强度的计算是高压绝缘设计中必不可少和极为重要的一部分，只考虑介质电导σ的恒流场计算和只考虑绝缘介质介电常数ε的静电场计算早已为人们熟悉[1]。然而实际的绝缘介质总是有一定的电导，特别是人为的通过一层薄的半导电层来调整绝缘子表面的电

位分布[2]，这时电导σ和介电常数ε共同起作用，单纯用恒流场或静电场计算已无能为力，必须采用准静态电场的计算方法。

2控制方程和计算方法

2.1电准静态场控制方程

2006年用户年会论文

电磁场是一种特殊的物质形态，Maxwell 方程组完全描述了电磁场的宏观性质。电磁场

理论有两个分支：一是从无线电频率到光频的电磁波传播问题，其中特征时间τ比电磁波通过系统的特征长度L 所需时间要短得多即τ<>L/c ,也就是准静态电磁场问题[5]。实际上很多电工（包括电力工程）学科中的电磁场问题，大多是准静态的。准静态场绝不是静态场，它是与时间有关的，并且不论那一类准静态场，都是既含有电场又含有磁场，只是电场和磁场之间没有耦合。

当磁感应强度随时间的变化可以忽略时，即

0≈∂∂t B

r ,称为电准静态，简称EQS ，这

时Maxwell 方程组退化为 ⎪⎪⎪⎩⎪⎪⎪⎨⎧=⋅∇=⋅∇≈×∇∂∂+=×∇ρD B E t D J H r r r r r r 00

对

t D J H ∂∂+=×∇r r r 求散度，可以有连续性方程:0(=∂∂+⋅∇t D J r r 由于0≈×∇E r ，电场E r 是无旋的，仍可以用标量位来描述V E −∇=r ，电场与磁场已经是非耦合的，当仅关心电场时，考虑到物性方程：E D E J r r r r εσ==,

可得同时考虑了材料的电导σ和介电常数ε的电准静态场的控制方程为： 0()(=∂∂∇⋅∇−∇⋅∇−t V V εσ

忽略电位随时间的变化它就成了恒流场的控制方程：0)(=∇⋅∇−V σ 在谐性电场分析中，控制方程为可写为：0

)()(=∇⋅∇+

∇⋅∇−V j V σωε

2.2 ANSYS 计算方法

ANSYS 软件是一个大型通用有限元分析软件，在我国有着广泛的应用，它也可以用来

分析电磁领域多方面的问题，自8.1版本后它增加了准静态电场的计算功能。

对于准静态时谐电场的计算，ANSYS 一方面扩充了原有的基于电荷的静电场单元

PLANE121,SOLID122, SOLID123的能力，使它们也能进行准静态时谐电场的计算；另一方面新增了基于电流的准静态电场单元，即PLANE230，SOLID231，SOLID232。对于准静态暂态电场的计算，ANSYS 只能用新增的准静态电场单元。

2006年用户年会论文

在准静态电场的计算中，要求输入的材料特性是介电常数PERX 和电阻率RSVX ，要

注意的是电阻率是电导率的倒数，而电工书上经常给出的是材料的电导率。

3 计算比较

3.1 时谐场

在均匀的时谐场中，放入一个介质球体，它的表面覆盖有一层薄的导电层，如图1所

示，

在这种情形下，如图1所示最大电场强度的解析解是[3]

))((2)2)(2()(6)2(303132313321333133max

εεεεεεεεεεεεεε−−−++−++=r r d d E E

式中： )3,2

,1(/)(0=+=+=x j r

r d d rx

x x r εωεσε

导电层中的功率损耗为：

a r E P 32036πσ=Δ

图１ 表面覆盖有一层薄的导电层的介质球体在均匀时谐场中

式中：

3313213323

2312322321

)

)((2)2)(2()11()

(2)1()

2(r r r d b d b d b a εεεεεεεεεεεεεε−−−++=−−−−+= 取表面厚度d=1mm, 对电导在10-4～10-8S/m 范围内的Emax 解析解的计算结果与有限元

的计算结果见图2，导电层中的功率损耗解析解的计算结果与有限元的计算结果见图3，电导在10-6时的等位线实部和虚部见图4