![电磁场数值计算及基于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=1739291079-gjWvGy9MElS2XU2fnSHw8LItIbstGylh-0-e4b371f1cb1e06c841e2d44061ac56d9)
式中,为电场强度(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=1739291079-uuVpAfco4V9jAJgkdQDXmVYZMnp0xcoN-0-1986fcd6d56bc334a68399dd2bcd66ae)
(2)静电场、恒定电场与静磁场 当电磁场量不随时间变化时,所有与时间导数相关的项为0,此时,式(1-1)、式(1-2)与式(1-5)可写作
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/16_08.jpg?sign=1739291079-3fe43AMp7Z9ur0RPztjQlQ5bMVVJNFzk-0-baeacc71c4e0947ba8716a87ae9d2b73)
而电介质中的高斯定理式(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=1739291079-i1nB64mSv8Ds3CLMEjXKwIXwoBKdvNal-0-252f2bef047613ea4079192490c7fc21)
式中,ω为角频率;j为虚数单位。时谐激励下,电场和磁场同时存在并相互激发。根据激励频率、所分析结构的电尺寸及区域中的材料参数特性,时谐电磁场可进一步分为时谐电场[忽略式(1-9)中的],时谐磁场[忽略式(1-10)中的
]及时谐全波电磁场[比如由式(1-9)和式(1-10)可以导出以
为变量的二阶矢量波动方程[8]]。
(4)本构方程 为了确定所有的电磁场量,还需要列出描述媒质宏观电磁性质的本构关系。电磁场的三个本构关系为
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_04.jpg?sign=1739291079-nhHQ9Ma2XFqQ7rsBd6k5S3078fCP5Xvp-0-6e760eab4de5eab5fae81226b00f8d63)
式中,本构参数ε、σ和μ分别表示媒质的介电常数(F/m)、磁导率(H/m)和电导率(S/m)。对于各向异性媒质,这些参数是张量;对于各向同性媒质,它们是标量。对于非均匀媒质,它们是空间坐标的函数;对于均匀媒质,它们不随空间坐标变化,取值为常数。另外,通常使用的电工磁性材料,其磁导率及电导率参数同时也是磁场强度及电场强度的非线性函数。
(5)不同媒质分界面的衔接条件 在电磁场分析中,如果场域内包含有不同的媒质(例如媒质1和媒质2),则一般需要确定电磁场量在不同媒质分界面处的衔接关系(连续性条件)。两种媒质交界面上的连续性条件表达如下:
对于电场有
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_05.jpg?sign=1739291079-kv7vDcl64dGGUaGK5M0nGfhm7wS52EC5-0-e0878c8b8b62c91a4c8c5de244b5a767)
对于磁场有
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_06.jpg?sign=1739291079-isrEXCne4PeuciAY8qJJSlc4PZIjmobD-0-4c62d4d0cc137c9186ad77cb37264eb8)
式中,为媒质分界面上的单位法向量,由媒质2指向媒质1,如图1-1所示。需要注意的是,以上衔接条件假设分界面处不存在线电流及自由面电荷。若分界面存在线电流密度
和自由面电荷密度ρs,则需对以上方程作以下修改[9]:
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_09.jpg?sign=1739291079-r6d58WHOIT33WWCMyiKHO9zBV4zublPg-0-deb13a9e5ab41d9e00f49865b75056ef)
![](https://epubservercos.yuewen.com/F7EA69/28831410007864506/epubprivate/OEBPS/Images/17_10.jpg?sign=1739291079-SFwQJnK0SqQh08uMjjlMQ2ZXbTzbYwhM-0-31855ee7e3ff53ff0707fb1c133e4d04)
图1-1 两种媒质界面的示意图