![电磁场数值计算及基于FreeFEM的编程实现](https://wfqqreader-1252317822.image.myqcloud.com/cover/797/49070797/b_49070797.jpg)
1.1 电磁场基本理论
1.1.1 麦克斯韦方程组
19世纪中叶,英国物理学家麦克斯韦在总结前人工作的基础上,提出了适用于所有宏观电磁现象的数学模型,称之为麦克斯韦方程组。科学技术发展的实践充分证明,麦克斯韦方程组正确描述了电磁场中各物理量之间的相互关系。麦克斯韦方程组有微分形式与积分形式,本节根据电磁场有限元数值分析格式推导的需要,仅给出其微分形式。
(1)Maxwell方程组的微分形式如下:
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/16_01.jpg?sign=1738869895-uvp4N0dhdJXA1ZW5voXQaE975vCycLqN-0-3382d08f81c6880ab37fe8eb9afe45c8)
式中,为电场强度(V/m);
为电通量密度(C/m2);
为磁场强度(A/m);
为磁通量密度(T;Wb/m2);
为传导电流密度(A/m2);ρ为电荷密度(C/m3)。另一个基本方程是电流连续性方程,即
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/16_07.jpg?sign=1738869895-xIjnWesezHhQd6tiT5efDacF8puDpWND-0-5ec7aedf5cf2cf8091e9c8e83190d1ae)
(2)静电场、恒定电场与静磁场 当电磁场量不随时间变化时,所有与时间导数相关的项为0,此时,式(1-1)、式(1-2)与式(1-5)可写作
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/16_08.jpg?sign=1738869895-xnU2gmK1n1pV7TbgGxwdM1cqnu47iw92-0-294c15f909bfe46b1a01c46409449303)
而电介质中的高斯定理式(1-3)、磁场高斯定理式(1-4)保持不变。此时,电场和磁场之间不存在耦合关系。将电磁场解耦,即得到静电场、恒定电场及静磁场的控制方程。其中,式(1-3)和式(1-6)描述了静电场,式(1-6)与式(1-8)描述了恒定电场,式(1-4)与式(1-7)描述了静磁场。
(3)时谐电磁场 当麦克斯韦方程组中的激励源做时谐变化,且所有材料参数均为线性时,所有的电磁场量也在该频率下做时谐变化,因此也叫作时谐电磁场。用复相量法,可将式(1-1)、式(1-2)和式(1-5)写成以下形式:
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/16_09.jpg?sign=1738869895-zGIkPOfnRFwatoICIsv0BCSRX3xFanLN-0-609e5be5e57682edd2adc20aa8840a44)
式中,ω为角频率;j为虚数单位。时谐激励下,电场和磁场同时存在并相互激发。根据激励频率、所分析结构的电尺寸及区域中的材料参数特性,时谐电磁场可进一步分为时谐电场[忽略式(1-9)中的],时谐磁场[忽略式(1-10)中的
]及时谐全波电磁场[比如由式(1-9)和式(1-10)可以导出以
为变量的二阶矢量波动方程[8]]。
(4)本构方程 为了确定所有的电磁场量,还需要列出描述媒质宏观电磁性质的本构关系。电磁场的三个本构关系为
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_04.jpg?sign=1738869895-OlO2AHWfHIPGNqMvKc4TcejuiwyRDrRi-0-489ace3ed046a6149d7e7287eb9e0c5d)
式中,本构参数ε、σ和μ分别表示媒质的介电常数(F/m)、磁导率(H/m)和电导率(S/m)。对于各向异性媒质,这些参数是张量;对于各向同性媒质,它们是标量。对于非均匀媒质,它们是空间坐标的函数;对于均匀媒质,它们不随空间坐标变化,取值为常数。另外,通常使用的电工磁性材料,其磁导率及电导率参数同时也是磁场强度及电场强度的非线性函数。
(5)不同媒质分界面的衔接条件 在电磁场分析中,如果场域内包含有不同的媒质(例如媒质1和媒质2),则一般需要确定电磁场量在不同媒质分界面处的衔接关系(连续性条件)。两种媒质交界面上的连续性条件表达如下:
对于电场有
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_05.jpg?sign=1738869895-K9ix5DVA60WHhsFtt1zL021OoO48DWhc-0-b8568b23b6973a2f1b283c7505c738c0)
对于磁场有
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_06.jpg?sign=1738869895-G5n5qUc1PVXBWIexXhuO0F8amNe96F6x-0-75486110599eebcd490a0c3b1ecae7bc)
式中,为媒质分界面上的单位法向量,由媒质2指向媒质1,如图1-1所示。需要注意的是,以上衔接条件假设分界面处不存在线电流及自由面电荷。若分界面存在线电流密度
和自由面电荷密度ρs,则需对以上方程作以下修改[9]:
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_09.jpg?sign=1738869895-8ZzAITHbvCuUv0q5sUE3rfQMkdRKo2eO-0-b677f060acdd76d87b582b2c8ea49cd1)
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_10.jpg?sign=1738869895-9S7E0ZnupB13f2p0wRJSU4XA48tk4h5W-0-46c7e27cc7d6cecdfe94bd02018b0eb8)
图1-1 两种媒质界面的示意图