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