电磁场有限元分析1VIP

§1电磁场基本概念及理论

1.1电磁场基本理论

1.2运动导体中的电磁场

1.3电磁场能量和力

1.4物质的磁性10/9/202210/9/20222/27/20231M&ISI电荷与电流在空间的每一点产生具有它本身真实性的场。电磁场是由电荷引起的，它对电荷又有作用。在我们的应用中，探讨的是宏观场，我们所探讨的场量是对无穷小的体积和时间间隔内场量的时间和空间的统计平均值。在宏观理论中的无穷小量，是一个抽象的概念，一方面它是小到不影响场中物理量连续变更的充分小，另一方面它与分子、原子、晶胞等微观结构相比又是充分大。麦克斯韦理论所依据的一些试验定律，忽视了物质构造的不匀整性，是把一些有关的物理量进行了统计平均得到的。1.基本概念1.1电磁场基本理论10/9/202210/9/20222/27/20232M&ISI牛顿力学中的长度、质量、时间和力等基本概念都是适用的。电荷的量值与运动无关

|q|=1.6008×10－19C。电荷守恒：

J是电流密度，ρ

是电荷的体积密度。在任何媒质中，电磁扰动以速度：相对于接收者传播。在电磁场中，作用于单位电荷的力为：电荷体密度大小为ρ处单位体积所受的力：2.基本假设1.1电磁场基本理论10/9/202210/9/20222/27/20233M&ISI电极化

电介质中电荷的分布不同于正常状态而发生畸变。可用电偶极子模型来描述,极化的程度可以用极化强度表示：

+-

即每单位体积的电偶极矩,极化后,在介质内部要引起作体分布的束缚电荷,在介质表面出现面分布的束缚电荷。表示电介质极化特性的参数：介电常数（电容率）。电位移矢量和电场强度的关系：当电介质为各向同性且线性时，可简化为

：3.媒质的极化和磁化

1.1电磁场基本理论10/9/202210/9/20222/27/20234M&ISI媒质的磁化可用磁偶极子模型描述。磁化程度可用磁化强度表示：表示媒质磁化特性的参数：磁导率。磁场强度和磁感应强度的关系：

当媒质为各向同性且线性时，可简化为：E和H是由场源产生的电磁场本身的特性，与媒质无关；D和B则是考虑受媒质极化和磁化影响后的电磁场特性。1.1电磁场基本理论10/9/202210/9/20222/27/20235M&ISI以E和B表示的麦克斯韦方程组：第四个方程中多出一项位移电流密度，是因为没有此项，就不能同时满足电荷守恒方程，也即右面三个方程不能同时成立：4.媒质中的电磁场

1.1电磁场基本理论10/9/202210/9/20222/27/20236M&ISI以E、D、B、H表示的麦克斯韦方程组：构成关系：上面七个方程构成了媒质中电磁现象严格的宏观描述的基础。1.1电磁场基本理论10/9/202210/9/20222/27/20237M&ISI赫姆赫兹定理：一个有限区域的矢量场由它的旋度和散度唯一确定矢量磁位A：1）库伦规范：对于静态场：因此矢量位A是唯一确定的。2）洛伦兹规范：

5.用矢量磁位A和标量位φ表示的电磁场方程组

1.1电磁场基本理论10/9/202210/9/20222/27/20238M&ISI对于时变场：1.1电磁场基本理论10/9/202210/9/20222/27/20239M&ISI对于非导电媒质（波动方程）：导电媒质中低频场（扩散方程）：1.1电磁场基本理论10/9/202210/9/20222/27/202310M&ISI用E、D、B、H表示的麦克斯韦方程组对于空间中媒质特性连续点处成立。在穿过不同媒质的分界面时，媒质参数ε、μ、γ要发生突变，因而场矢量会出现相应的不连续。设想用一很薄的过渡层代替媒质1和媒质2的分界面，媒质参数ε、μ、γ很快地但是连续地从媒质1中值变到媒质2中去。在过渡层内，场矢量及其一阶导数是连续、有界的。

6.不同媒质分界面上的边界条件

1.1电磁场基本理论10/9/202210/9/20222/27/202311M&ISI当，柱壁的贡献为无限小；则那么时，1.1电磁场基本理论(1).B满足的边界条件10/9/202210/9/20222/27/202312M&ISIn是分界面的法线单位矢量，由媒质1到媒质2；n0是切向单位矢量，与回路的绕向成右螺旋。1.1电磁场基本理论(2).H满足的边界条件10/9/202210/9/20222/27/202313M&ISI而，因此1.1电磁场基本理论（a）若假定和有界，当，时，式中的面积分都为零，两条侧边的贡献也为零。=?10/9/202210/9/20222/27/202314M&ISI(b)若假定时，原先通过回路的电流总量不变，被挤压到S面上的无限薄层里，则分界面S上，也就是说，必需用面电流密度K来表示：因此，当，时，对应的矢量:1.1电磁场基本理论冲突10/9/202210/9/20222/27/202315M&ISI1.1电磁场基本理论（磁通连续性）B的法向重量连续H的切向重量连续，只有在跨越志向完纯导体或超导体的边界时，K≠0，它才会不连续假如两种媒质都不是志向导体，则：（H的切向重量连续）（B的法向重量连续）因此：10/9/202210/9/20222/27/202316M&ISI1.1电磁场基本理论上式描述了磁场的折射性质：假如媒质1是非磁性的，而媒质2是铁磁性的则，，这意味着，对于随意一个不接近0的角度，在铁磁性媒质中，磁场几乎是与分界面平行；假如媒质2是非磁性的，而媒质1是铁磁性的则，，这意味着，假如磁场起源于铁磁性媒质，则磁通将以近似垂直的角度穿出分界面。10/9/202210/9/20222/27/202317M&ISI1.1电磁场基本理论（4）电场矢量满足的边界条件：（3）用矢量磁位A表示的边界条件：10/9/202210/9/20222/27/202318M&ISI对于低速运动状况下的电磁场，分析计算的方法之一是：选定一个静止的坐标系，运动物体对于该坐标系有相对运动。计算时仍从静止系统的麦克斯韦方程动身，但要计入自由电荷、极化电荷和物质磁化等因素由于运动而引起的附加效应。也就是说，电磁场的场源，除了自由电荷、传导电流和运流电流以外，还将出现若干附加项表示由于运动引起的影响。众所周知，在全部惯性坐标系中，麦克斯韦方程具有相同的形式，在与电磁场有关的问题中，伽利略变换不适用。但是对于低速运动状况下的电磁场，运用伽利略变换仍旧可以得到特别精确的结果。1.2运动导体中的电磁场10/9/202210/9/20222/27/202319M&ISI其中而在两个参照系中的构成方程分别为：对于在电磁场中匀速运动的导体，如图所示，坐标系x1-y1-z固结在转子上，相对于惯性系x-y-z以角速度绕z轴匀速转动，电磁场量在两个系中的伽利略变换为：1.2运动导体中的电磁场10/9/202210/9/20222/27/202320M&ISI从应用角度来说，很多电磁器件就是为了传递能量或进行能量转换而设计的。因此必需关注电磁场中的能量和力。另一方面，从能量的观点考虑问题有时可以很简洁得到场的重要性质。有电介质时的电场能量媒质中的磁场能量1.电磁场能量

1.3电磁场能量和力10/9/202210/9/20222/27/202321M&ISI力学中的应力张量和表面应力应力张量T：2.静态场中的力

1.3电磁场能量和力10/9/202210/9/20222/27/202322M&ISI表面应力（矢量）Tn：力平衡力矩平衡(体积力密度与应力张量间的关系)(对称张量，九个重量，只有六个独立)1.3电磁场能量和力10/9/202210/9/20222/27/202323M&ISI电磁应力张量和表面应力依据场的观点，全部的力都要通过场以连续的方式从一个物体传递到另一个物体。法拉第认为：电场和磁场处于一种特殊的张紧状态。把一个处于静平衡下的系统，用一个闭合面A把它一分为二，则部分2对部分1作用的合力，应按某一方式通过表面A。麦克斯韦证明白：部分1所受到的体积力可以用包围部分1的边界面A上的表面力来表示。力平衡关系：磁场中的表面应力：A211.3电磁场能量和力10/9/202210/9/20222/27/202324M&ISI电磁应力张量：T是正方向为单位面积外法线时的应力张量形式,S是常用的电磁应力张量形式。(两者相差一个负号)表面应力：1.3电磁场能量和力10/9/202210/9/20222/27/202325M&ISI抗磁性：,惰性气体顺磁性：铜、铝、不导磁钢等，铁磁性：铁、硅钢、碳钢等，反铁磁性：其磁化率的大小近乎于强顺磁体，属弱磁性铁氧体磁性：与铁磁性有一些相像的特性，如自发磁化、磁滞、临界温度等。1.4物质的磁性10/9/202210/9/20222/27/202326M&ISI1.4物质的磁性自发磁化：不论是在居里点以上还是以下，铁磁质中都存在有分子场作用，而且相当强。甚至在没有外磁场作用的状况下它使物质内部磁化到饱和。磁畴：在去磁状态下，铁磁质中有很多称为“磁畴”的小区域，每一磁畴都自发磁化到饱和。在不同的磁畴中饱和磁化强度的方向各不相同，宏观看来材料的净磁化强度为0；材料的磁化过程，就是从磁化强度方向各异的多磁畴状态变为与外磁场同方向的单磁畴状态的过程。磁化：将它们放入磁场就会在磁场相同的方向上强力磁化;非线性：磁通密度对磁场强度的变更，除了在小范围内是线性的以外，都是非线性的；饱和、磁滞、剩磁等特性；铁损耗：当经受随时间变更的磁化时要消耗能量。铁磁性的基本特点：微观机理：惠斯在1906年提出了分子场理论。基于对铁磁现象的探讨不断深化，他提出了如下假设：10/9/202210/9/20222/27/202327M&ISI1.4物质的磁性1.磁化曲线的非线性

起始磁导率区：在外磁场的作用下，磁畴壁做可逆移动;不行逆磁化区：这是畴壁做不行逆跳动式移动的过程，磁化强度随着磁场增大急剧增大；旋转磁化区：畴壁位移基本完毕，只有靠磁畴内磁矩的转动来进一步增大磁化强度；接近饱和区：磁畴内磁矩的可逆转动造成，此时增大磁场，磁化强度增加很小。(a)M=0(b)M>0(c)M=Mscosθ (d)M=MsMsMsMsMsHMsHMsH10/9/202210/9/20222/27/202328M&ISI1.4物质的磁性2.磁滞回线

剩磁：铁磁体经过外磁场的磁化达到饱和后，将外磁场移去，其磁化强度不为零，而是停留在Mr处，称之为剩余磁化强度Br；内禀矫顽力要在反方向再加外磁场后才能使磁化强度渐渐回复到零，这时的外磁场称为内禀矫顽力Hc；假如反向磁场再增大，即在与原来相反的方向达到饱和，磁场回到零，这时的磁化状态是处于反方向的剩磁状念(-Mr)再以原来的磁化方向渐渐加大磁场，通过正向矫颓力点再度饱和。10/9/202210/9/20222/27/202329M&ISI1.4物质的磁性3.硬磁材料和软磁材料

工程上常用的磁性材料通常分软磁材料和硬磁材料两类。软磁材料：磁性能的主要特点磁导率高，矫顽力低。对某些软磁材料来说，还要求磁化损耗尽可能低；硬磁材科：磁性能的主要特点是矫顽力高，剩磁感应强度和磁能积也是其最重要的参数；10/9/202210/9/20222/27/202330M&ISI§2ANSYS低频电磁场有限元分析

2.1ANSYS低频电磁场有限元分析概要

2.2分析方法

2.33D载流体的建模

2.4对称性和边界条件

2.52D分析常用单元简介10/9/202210/9/20222/27/202331M&ISI2.1ANSYS低频电磁场有限元分析概要1.主要探讨对象：电机/发电机变压器螺线管作动器2.主要物理量：磁感应强度（磁通密度）磁场强度磁力和力矩阻抗滤波器磁成像系统连接器感抗涡流功耗漏磁10/9/202210/9/20222/27/202332M&ISI3.电磁分析中系统部件按电磁属性的分类1）按电特性分：绝缘体：不存在涡流绞合导体：不存在涡流无涡流效应实心导体包含涡流效应实心导体2）按磁特性分：不导磁体：空气、铜、铝等软磁体：铁、钢硬磁体：永磁材料2.1ANSYS低频电磁场有限元分析概要10/9/202210/9/20222/27/202333M&ISI基本分析方法如下表，依据电流引入的不同方式，依据不同问题选用相应的分析方法。2.2分析方法10/9/202210/9/20222/27/202334M&ISI2.2分析方法10/9/202210/9/20222/27/202335M&ISI2.2分析方法10/9/202210/9/20222/27/202336M&ISI1.电流驱动多股绞线－电流密度已知J：电流密度N：线圈匝数i：单股导线电流A：线圈横截面积通过施加实体载荷给定每个单元的电流密度,吩咐：BFE,,JS2.33D载流体的建模10/9/202210/9/20222/27/202337M&ISI2.SOURC36:电流驱动多股绞线2.33D载流体的建模10/9/202210/9/20222/27/202338M&ISI3.载流实心导体4.电压驱动实心导体5.电路驱动实心导体2.33D载流体的建模10/9/202210/9/20222/27/202339M&ISI

利用物理问题本身的对称性可以简化模型，提高分析效率。轴对称：3D问题可简化为2D奇对称：对称面两边的场分布符号相反偶对称：对称面两边的场分布相同周期对称：周期问题只需分析其中一个周期区域2.4对称性和边界条件10/9/202210/9/20222/27/202340M&ISI边界条件(以2D分析为例)：2.4对称性和边界条件10/9/202210/9/20222/27/202341M&ISI1.PLANE538节点，4自由度：AZ，MMF,CURR,VOLT2.52D分析常用单元简介KEYOPT(1)：单元自由度0--AZ:静态区域,涡流区域1--VOLT,AZ:载流实心导体2--AZ,CURR:电压驱动绞线3--AZ,CURR,EMF:电路耦合绞线4--AZ,CURR,EMF:电路耦合实心导体10/9/202210/9/20222/27/202342M&ISIKEYOPT(2)单元速度:0--单元无相对运动1--经典速度公式(KEYOPT(1)=0，1)KEYOPT(3)单元几何特性:0--平面1--轴对称当包含速度效应时(KEYOPT(2)=1)：单元只能有AZ或AZ、VOLT自由度各向同性电阻率磁雷诺数远大于1时，计算精度会降低KEYOPT(1)≥2或KEYOPT(2)≥1时，系统矩阵是非对称的2.52D分析常用单元简介10/9/202210/9/20222/27/202343M&ISI假设和限制：单元面积必需大于零;在做轴对称分析时，在全局坐标系下的Y轴必需是对称轴（轴向），而且模型必需在左半平面(X>0);当单元只有AZ自由度(KEYOPT