电磁场数值计算之1-西安交通大学电气工程学院.doc

第一章 电磁场基本概念11 Maxwell 方程组（一） maxwell方程微分形式 积分形式全电流定律 （1-1） 电磁感应定律 （1-2）高斯定律 （1-3）磁通连续性原理 (1-4)电流连续性方程 （1-5）说明：1、四个方程的物理意义，电生磁，磁生电，预言电磁波；积分形式（环量与旋度，通量与散度之间的关系）、复数形式（可作为稳态场计算）；梯度、散度、旋度的概念（描述“点”上电磁场的性质）。2、方程（1-1）、（1-2）、（1-5）是一组独立方程，其它两个方程可以由此推出。但独立方程有6个变量（），因此，方程数少于未知量，是非定解方式，必须加本构方程才为定解形式，对于简单媒质，本构方程为 (1-6)3、材料性质材料是均匀的 ， ，材料是非均匀： ，材料是各向异性：材料参数用张量形式表示 ， 材料为非线性：材料参数是未知函数的函数 ， （1-7） 4、直接求解矢量偏微分方程不易：一般矢量方程要转化为标量方程才能求解，另外，在边界上不易写出场量边界条件，因此，常化为位函数的定解问题（位函数容易确定边界条件），通过位函数与场量的关系 （1-8）得到场量。12 偏微分方程的基本概念1.2.1 偏微分方程的基本概念微分方程分为常微分方程和偏微分方程（又分为描述不同物理现象的椭圆型方程、双曲型方程、抛物型方程及其线性和非线性方程），电磁场问题多为偏微分方程问题。1、 常微分方程未知函数是一元函数（即一个变量的函数）的微分方程（组）。如R、L、C串联电路是两阶常系数非齐次微分方程， （1-9）对于一个n阶场微分方程，通常可将其分解为有n个任意常数的通解形式，根据初始条件解出常数。2、 偏微分方程未知函数是多元函数的微分方程，如 。又分为线性和非线性偏微分方程，除了极有限的问题可以用分离变量法求解外，多数问题难以用解析表达式表示。（1） 线性偏微分方程设 ，（如：，如：），则 （1-10）中，如果a,b,c,d,e,r,s与p无关，只是x,y的函数，则称式(1-10)为线性微分方程。（2） 非线性微分方程a,b,c, d,e,r,s,f中只要有一项不满足上述条件，或未知函数及其偏导数是非线性的微分方程，则都称为非线性微分方程。如恒定磁场中的定解问题 如：在电磁场中，若，或媒质不均匀时，均为线性方程。若，或，则为非线性方程。1.2.2 偏微分方程的分类宏观电磁场都是二阶微分方程，下面以二阶电磁场偏微分方程为例，看偏微分方程的不同类型所反映的物理现象。以二元函数为例，y可以是时间变量t，那么偏微分方程的普遍形式为 最高阶项称为主部，主部决定着公式所代表的物理特性： 椭圆型方程，如 ， 双曲型方程，如 ， 抛物型方程，如 ，1、椭圆型方程如泊松方程、拉普拉斯方程 （ 与椭圆方程 形象对比）特点：所有二阶偏导数的系数同符号，描述的物理现象：描述平衡、定常的稳定状态，因此方程与时间无关，定解条件中只有边界条件，没有初始条件。如重力场、静电场、恒定电场、恒定磁场、稳定温度分布过程。2、双曲型方程如波动方程 无损耗，无激励源（ 与双曲型方程 形象对比）特点：对时间的偏导数系数与对空间偏导数的系数相差一负号。描述波的传播过程，它具有对时间可逆的性质（用（-t）代入方程后，方程不变）如：弦振动、膜振动、声波、电磁波。3、抛物型方程如，热传导方程 扩散率或导温系数涡流方程 ， ， （ 与双曲型方程 形象对比）特点：对时间变量的二阶导数为零。描述各种场的扩散过程，它具有对时间不可逆的性质。1.2.3 定解问题1、初值问题只有初始条件，没有边界条件的定解问题。如电路中的过渡过程问题、无界空间电磁波传播问题等。2、边值问题只有边界条件，没有初始条件的定解问题。如静电场、恒定电场、恒定磁场等问题。3、 混合问题既有边界条件，又有初始条件的定解问题，又称定解问题。如电气设备中的瞬态电磁场问题等。4、 解的稳定性问题如果定解条件的微小变化只引起方程的解在整个定义域中的微小变化，称其解是稳定的。反之称为不稳定解。（第1次课）13 电磁场中的定解问题定解问题 = 泛定方程+定解条件（初始条件+边界条件）下面先介绍各种场的泛定方程，然后介绍各类边界条件。1.3.1 静态、稳态电磁场中的泛定方程1、静电场方程静电场的基本方程 泊松方程 三维方程 若是均匀、各向同性介质，上式为 椭圆型方程静电场方程是椭圆型方程，只有边值问题。2、 稳态电流场问题稳态（直流）电流场满足的基本方程： 说明在导电媒质中，电流不会自成闭合回路（从电源正极出发到电源负极终止），电位满足拉普拉斯方程 椭圆型方程若是均匀、线性、各向同性介质，上式为 产生该电流场的源往往需要借助边界条件引入。3、 稳态磁场稳态（直流）电流产生的磁场满足的基本方程 （1） 标量磁位的泊松方程当求解区域内，那么，必定存在一个标量函数，使得根据，上式为拉普拉斯方程 椭圆型方程 上述方程只能用于的单连通域（见雷银照教材），因此应用的局限性较小。当磁场区域内存在铁磁质时，展开后为非线性方程为：若为线性，则为拉普拉斯方程：若已知磁化强度M，那么代入，得到此时，媒质的磁导率为。（2） 矢量磁位的泊松方程根据，有双旋度方程取库伦规范，及矢量恒等式 ，得 矢量泊松方程 若为线性、均匀媒质 若存在铁磁质，可将其作用等效为磁化电流的作用，它与磁化强度的关系为磁矢位A的方程可以写为真空中的泊松方程1.3.2 交变电磁场中的泛定方程时变场中，（下面分段没有绝对的分界线）缓慢变化 ( f 10 KHz )快速变化准静态场准静态场 电磁波MQS：EQS： MQS场求解时，磁场可以用稳态磁场的方法求解，然后用上述公式求电场；EQS场求解时，电场可以用静电场的方法求解，然后用上述公式求磁场。1、 扩散方程（抛物型方程）忽略位移电流，MQS场的方程为由此得到的扩散方程为（对第一式再取旋度） 非线性介质 ， 线性介质 ， 若为正弦交变场，扩散方程为 ， 涡流损耗是引起导体发热的主要原因。2、 波动方程（双曲型方程）一般不考虑非线性问题，因为如果在铁磁材料中传播电磁波，高频下的涡流损耗及磁滞损耗很大，电磁波很快衰减，能量不可能传递很远。因此，场量的波动方程取洛伦兹规范，则位函数满足的波动方程1.3.3 定解条件1、初始条件（柯西问题）在瞬态电磁场中，初始条件是整个系统初始状态的表达式扩散方程初始条件： 波动方程初始条件： 如：初始的速度、电流、电压等。2、边界条件（1）第一类边界条件（Dirichlet 狄利赫里条件） 强加边界条件例1 铁磁体的磁场和电容器的电场（二维） 图1-1第一类边界条件 （a）磁场问题；（b）静电问题在距离磁体足够远的地方，设磁力线平行于边界，因此可以假设。在距离电容器足够远的地方，设等位线平行于边界，可以假设。关键问题是第一类边界条件取得多远，才能保证计算精度。例2 电机的磁场图1-2（a）、（b）：需要考虑定子外的漏磁，因此，第一类边界条件取在大于定子外径20%之处，磁力线于边界平行，可以设A=0。图1-2（c）、（d）：如果定子深度饱和，漏磁很小，可以忽略，可将第一类边界条件取在定子外径，减少计算量。图1-2（e）、（f）：如果要分析远场，第一类边界条件可以取在大于定子外径56倍之处，如图（e）所示。或者用开于边界条件，如Kelvin-transformation边界（后面介绍），边界可以小一些，如图（f）所示。 (a) (b) (c) (d) (e) (f)图1-2 电机的磁场计算（第一类边界条件） （第二次课）（2）第二、三类边界条件（Neumann 聂以曼条件） 自然边界条件 自然边界条件（先复习平行平面场、轴对称场及两种场的等A线方程，球对称场）已知边界上的法向导数，代表着几种物理意义：(a) 已知边界上的激励源：面电荷分布或面电流分布 已知面电荷分布 已知面电流分布 以平行平面场为例：如：面电流在边界XOZ平面上，沿x方向，计算场域在区域，区域场为零（如电流方向相反的一对平行平板电流）（）那么 则 第二类非齐次边界条件如果面电流在YOZ平面，电流方向不变，沿y方向，计算场域在区域那么 则 第二类非齐次边界条件(b) 已知边界为场分布的对称线（面）， 第二类齐次边界条件，自然边界条件（提问1：磁场对称问题。若平行平面场中，对称线上磁力线与其平行，是哪一类边界条件？（对A而言，由于磁力线为等A线，所以是第一类边界条件，可以设为参考位。对而言，是第二类边界条件，在Ansoft中称为偶对称）。垂直呢？（如由于磁力线垂直于对称线，即其切线分量等于零，所以，所以是第二类齐次边界条件，如例3。对而言，是第一类边界条件，。在Ansoft中称为奇对称）。轴对称场中磁场沿轴对称？（如螺线管，磁力线垂直于对称线，所以。如磁力线平行于对称线，则，是第一类边界条件）。提问2：电场对称问题。若对称线上电力线与其平行，是哪一类边界条件（齐次第二类边界条件）？垂直呢（第一类，可以设为参考位）？如例4例3 如图1-3所示，对称线上磁力线与对称边界垂直，即对称线上只有法向分量根据，且二维场中，则有 由于，所以，在对称边界上 图1-3 齐次第二类边界条件例4 如图1-4（）所示，对称线上电力线与之平行，即只有切向分量，等位线与之垂直，根据 由于，所以，在对称边界上。还可进一步简化，再利用另一对称边界，如图1-4()所示，电力线与之垂直。在100V电压之间的0V等位线可以作为第一类边界条件。这样只需计算四分之一区域，计算量大大减少。 (a) (b)图1-4对称线边界条件例5 计算E型、U型电磁铁的磁场分布（近似为二维场，）两种模型无限远处的磁感应强度为零，取计算场域足够大时，可以认为模型的截断边界上磁感应强度为零，所以可以取A=0。对于E型电磁铁，对称轴为y轴，磁力线与y轴重合，是等磁位线。且电磁铁左右两部分中的电流方向相反，两边A值相等，符号相反，故y轴上的A=0，是第一类边界条件。图1-5 E型、U型电磁铁对于U型电磁铁，对称轴为y轴，但磁力线与y轴垂直，即y轴与等A线垂直（平行平面场中B线就是等A线），因此在y轴上，在有限元法中，这类边界节点不需要处理，按内部节点对待（Ansys中是默认值）。如果这两种电磁铁的结构还具有上下对称的特点，那么，B线与x轴垂直，在x轴上满足第二类齐次边界条件。例6： U型电磁铁如图1-7所示（三维），只需计算四分之一区间，即， 区间（第一象限），三维边界条件分为（三个分量都应有表达式） 1、无限远边界条件 无限远处，所以取三个截断边界面上：2、边界条件（也可以是的边界）U型磁铁对平面呈对称，即在边界面yoz上，只有磁场垂直于分界面，即只有分量，故Ax=0，且，根据，展开 （1）所以 第一式按强加边界处理，后两式是自然边界条件，且为齐次，不需要专门处理。 图1-7 U型电磁铁边界条件 2、边界条件（也可以是的边界）磁场强度平行于边界面xoz（），根据A与B的正交关系，所以在该平面上。但因为，无法从A的旋度中写出Ay分量的边界条件（只有是已知的，而。反之，如果），所以，考虑A的散度，稳态场取库伦规范，即由于在边界面xoz上，Ax、Az分量处处为零，因此其偏导数也为零，也是齐次第二类边界条件，最终边界条件为第一、三式是第一类边界条件，第二式是齐次第二类边界条件，不需专门处理。例7 用三维静磁磁矢位有限元法对CJ20-25交流接触器E形磁系统进行计算，如图1-8(a)所示。由于几何结构前后、左右对称，因此只需取四分之一空间计算，如图1-8(b)所示。同上分析，无限远处，所以 （a） （b） 图1-8 CJ20-25交流接触器E形磁系统 B平行于对称面yoz平面，因此该平面上，即，从A的旋度中无法得到齐次第二类边界条件，从A的散度中得到Ax的表述如图示。B又平行于对称面xoz平面，因此该平面上，即，从A的旋度中无法得到齐次第二类边界条件，从A的散度中得到Ay的表述如图示。例8 圆盘形电磁铁的磁场计算。采用圆柱坐标系，设矢量磁位A和电流密度J垂直于zor平面，只有方向的分量，即，。微分方程为令，则轴对称场的微分方程可以写