玻尔百科电磁学中的有限元方法 是一种通过将连续物理区域划分为简单有限单元网格来近似解决复杂问题的数值方法。该方法利用简单函数表示每个单元内部的物理行为，并将其汇编成由刚度矩阵描述的全局方程组，以求解电磁场分布。作为一种通用的多物理场模拟工具，它支持处理非线性材料属性、边界条件应用以及电与热传导等耦合现象。
电磁学，由麦克斯韦方程组所统治的优美领域，在理论上为我们描绘了电场、磁场与物质相互作用的完整图景。然而，当我们将目光从理想化的教科书问题转向现实世界——从不规则形状的天线，到材料复杂的电子芯片——这些优雅的方程往往变得难以甚至不可能求得解析解。为了跨越理论与应用之间的鸿沟，有限元方法（Finite Element Method, FEM）应运而生，它以“分而治之”的巧妙思想，成为了现代工程与科学研究中不可或缺的数值模拟工具。
尽管FEM功能强大，但许多使用者仅将其视为一个“黑箱”。本文旨在揭开这个黑箱，带领读者深入理解其内在的物理直觉和数学构造。我们将通过两个核心章节来探索FEM的世界。在第一章中，我们将剖析其原理与机制，从根本的离散化思想，到单元内部的场近似，再到全局方程的组装过程。在第二章中，我们将展示其广泛的应用与跨学科连接，看它如何解决从电子元器件设计到天体物理学的各类实际问题。通过本文，你将不仅学会使用FEM，更能深刻领会其精髓。
准备好了吗？现在就让我们卷起袖子，像打开一台精密钟表的后盖一样，一探有限元方法内部精巧绝伦的“原理与机制”。
现在，让我们卷起袖子，像打开一台精密钟表的后盖一样，一探其内部精巧绝伦的“原理与机制”。我们不仅要看它如何工作，更要理解它为何如此工作——这其中蕴含着深刻的物理直觉和数学之美。
想象一下，你想要精确计算一个形状不规则的湖泊的总体积。你不可能用一个简单的公式来完成任务，因为湖的边界和深度变化万千。一个聪明的做法是什么？你可以将湖面划分为许多小的、规则的方格。在每个方格内，你假设湖底是平坦的，然后测量这个方格中心的深度，用“面积 × 深度”来估算这一个小方块的体积。最后，将所有小方块的体积加起来，你就得到了整个湖泊体积的一个很好的近似值。方格划分得越细，结果就越精确。
这，就是有限元方法的核心思想：离散化 (Discretization)。面对一个连续复杂的电磁场问题——比如计算一个奇形怪状的电容器内部的电场分布——我们无法直接得到一个解析解。于是，我们采取一种“分而治之”的策略：将整个求解空间（无论是二维平面还是三维空间）切割成成千上万个小的、简单的、互不重叠的几何单元，比如三角形（二维）或四面体（三维）。我们将这些小单元称为“有限元”。
一旦我们将复杂的世界分解为简单的单元，下一步就是思考在每一个单元内部，物理定律是如何表现的。有限元法的“魔法”在于，它假设在每一个微小的单元内部，物理量的变化是极其简单的。例如，在求解静电场问题时，我们关心的是电势 的分布。在一个小小的三角形单元内，最简单的假设莫过于电势 是一个线性函数，就像一块倾斜的平玻璃板。
这个线性假设 带来了巨大的好处。我们知道，电场 是电势的负梯度，即 。对一个线性函数求梯度，结果是什么？是一个常数向量！这意味着，根据我们的模型，在每一个三角形单元内部，电场是均匀的、恒定不变的。整个复杂的、处处变化的电场，被我们近似成了一系列“马赛克”，每一块“马赛克”内部的电场都简单到不能再简单——均匀分布。
当然，我们知道真实世界的电场并非如此。但是，就像用足够多的平面小镜子可以拼凑出一个光滑的曲面反射镜一样，只要我们的三角形足够小，这种分片常数（piecewise-constant）的电场近似就能以极高的精度逼近真实的电场分布。例如，如果我们通过计算得知了一个三角形单元三个顶点（我们称之为“节点”）的电势分别为 ，我们就能唯一地确定这个倾斜的“电势平面”，从而计算出这个单元内部那个恒定的电场 。
好了，我们知道了在单个单元里发生了什么。但这些单元并非孤岛，它们通过共享的节点和边界连接在一起。物理世界的连续性要求在相邻单元的边界上，电势必须是连续的。这个看似简单的要求，正是将成千上万个孤立的单元方程“缝合”成一个宏大的全局方程组的关键。
这个过程被称为组装 (Assembly)。想象一下，每个三角形单元都对应一个小的 矩阵，我们称之为“单元刚度矩阵”。这个矩阵描述了该单元内部三个节点之间的相互“影响”——给一个节点施加一定的电势，会对其他节点的电荷产生怎样的影响。
现在，我们要构建一个巨大的“全局刚度矩阵”，它描述了整个系统中所有节点之间的相互作用。这个全局矩阵的某个元素 代表了节点 和节点 之间的相互影响。这个影响是如何计算的呢？非常直观：它等于所有同时包含节点 和节点 的单元贡献的总和。这就像在一个社交网络中，两个人之间的关系强度，是他们共同参与的所有小圈子（家庭、公司、兴趣小组）中关系强度的叠加。通过这个“对号入座”并累加的过程，我们将所有局部信息汇集成一个描述整个系统行为的宏伟蓝图。
最终，我们得到一个形式极其简洁的线性方程组：
这里， 是一个列向量，包含了所有待求解的节点电势；而 是源向量，代表了系统中的电荷源或施加的边界条件。解出这个方程组，我们就得到了整个系统在所有节点上的电势分布，这也就意味着我们近似地得到了整个区域的电势场。
求出所有节点的电势 只是第一步。有限元法的真正魅力在于，我们可以从这些数字中提取出丰富的物理信息。例如，我们可以计算存储在每个小单元内的静电能量。其表达式异常优美：，这里 是该单元的节点电势向量。整个系统的总能量就是所有单元能量的总和。 这再次体现了“分而治之”的思想，复杂的整体能量被分解为简单的局部能量之和。
然而，我们必须时刻保持清醒：FEM 的解是一个近似解。一个至关重要的问题是：我们的解有多好？我们如何改进它？
一个朴素的想法是，不断加密网格，让三角形越来越小。通常，当单元的特征尺寸 减小时，计算结果与真实解之间的误差 会以幂函数的形式下降，即 。这里的指数 被称为“收敛阶”，它量化了我们每付出一份努力（加密网格），能够换来多少回报（精度提升）。通过在不同网格密度下进行计算，我们可以验证我们的仿真是否在正确地“收敛”于真实解。
更进一步，我们能否让计算机变得更“聪明”，只在“最需要”的地方加密网格呢？答案是肯定的，这就是自适应网格加密 (Adaptive Mesh Refinement) 的思想。在物理上，误差大的地方通常是场变化剧烈的地方，比如导体尖端。在我们的近似解中，这些地方会表现出一种“不和谐”：物理量在跨越单元边界时会产生不应有的“跳变”。例如，电位移矢量 的法向分量在没有表面电荷的介质交界面上本应是连续的，但在我们的分片近似解中，它在单元边界上可能会出现一个不为零的“跳变”。这个跳变的大小，就成了一个绝佳的后验误差指示器 (a posteriori error indicator)。我们可以指示计算机：去寻找那些 法向分量跳变最大的边界，然后在这些边界附近自动加密网格。 这样，计算资源就被精准地投入到最关键的区域，好比一位经验丰富的艺术家，只在画作的关键细节上精雕细琢。导体尖端的“避雷针效应”——电场急剧增强——就可以通过这种智能的网格加密被精确地捕捉到。
有限元方法的强大不仅在于其数学核心，还在于其处理复杂现实世界场景的灵活性。
模拟“无限”：许多电磁问题，如天线辐射，都发生在开放空间中。我们无法用有限的单元去填充整个宇宙。怎么办？我们可以在离研究对象足够远的地方画一个“人工边界”，然后在这个边界上施加一个特殊的边界条件，用来模拟电磁波“一去不复返”的效果。这些被称为“吸收边界条件”或“辐射边界条件”，它们像是计算区域的“消音墙”，防止电磁波反射回来污染计算结果。选择不同的边界条件，比如简单的强制电势为零（Dirichlet 条件）或更复杂的 Robin 条件，会影响到我们对无限远处的模拟精度。
利用对称性：如果一个问题具有几何和激励上的对称性，比如一个方形的亥姆霍兹线圈，我们完全没必要计算整个三维空间。我们可以只计算其中的八分之一（一个卦限），然后在对称面上施加正确的对称边界条件。比如，在一个平面上，如果磁场线总是与之平行穿过（法向分量为零），我们就施加 Neumann 型边界条件；如果磁场线总是垂直于该平面（切向分量为零），则施加 Dirichlet 型边界条件。正确利用对称性，能将计算量减少数倍甚至数十倍，这是每个优秀数值工程师的必备技巧。
处理悬浮导体：在一个电路中，我们可能会遇到一块孤立的导体，它不与任何电源相连，但我们知道它上面带有总共为 的净电荷。它的电势 是多少呢？这是一个“悬浮电势”问题。在有限元框架中处理它异常简单：我们只需将控制这个导体节点的那一行方程的右端项（源项）直接设为 即可。这体现了 FEM 框架的巨大灵活性，它能自然地将关于电荷的约束直接整合到求解电势的方程组中。
当我们从静力场问题迈向更复杂的电磁波和谐振问题时，比如求解一个三维谐振腔的本征模式，一个诡异的“幽灵”出现了。如果我们天真地将处理电势（标量）的方法直接套用在电场强度（矢量）上，即独立地对电场的三个分量 使用我们之前讨论的“节点”基函数进行离散化，计算结果中会混入大量物理上根本不存在的“伪解”(spurious modes)。这些伪解像噪声一样污染了真实的谐振频谱，让结果变得毫无意义。
这个问题的根源非常深刻。它出在一个被我们忽略的矢量分析恒等式上：任意标量场 的梯度的旋度恒为零，即 。我们那种朴素的、基于节点的矢量场离散化方法，破坏了这一神圣的数学结构！它构建的离散空间中，存在一些“离散旋度”不为零的“离散梯度”场，或者说，存在一些“离散旋度”为零但又不是任何“离散梯度”场的“幽灵场”。这些不满足 这一基本法则的“幽灵”正是伪解的来源。
为了驱逐这些幽灵，数学家和物理学家们发明了一种更精妙的有限元——Nedelec 边元 (Edge Elements)。它的思想堪称革命性：不再将自由度（待求未知数）放在几何体的“节点”上，而是放在“边”上！对于矢量场 ，我们求解的不再是它在顶点上的分量，而是它沿着每个单元的棱边的切向分量的积分。
这种看似奇怪的做法，却有着神奇的效果。它天生就能保证电场 的切向分量在单元与单元之间是连续的——这恰好是麦克斯韦方程组要求的物理边界条件！更重要的是，基于“边”的函数空间完美地重建了 这一离散化的矢量恒等式。它从根本上铲除了“伪解”幽灵滋生的土壤，使得计算结果干净、可靠。
从简单的三角形近似，到组装宏大的矩阵方程，再到与“伪解”幽灵的巧妙搏斗，有限元方法向我们展示了一幅从简单直觉出发，通过严谨的数学构造，最终能够精确描绘复杂物理世界的壮丽画卷。它不仅是一个计算工具，更是一种思想，一种将无限复杂性分解为有限简单性的智慧结晶。
在前一章中，我们已经拆解了有限元方法（FEM）的内部机制——就像钟表匠审视齿轮和弹簧一样。我们看到，通过将复杂问题分解为微小的、可管理的单元，我们能够将描述物理世界的偏微分方程转化为计算机能够求解的代数方程组。这是一个强大而优美的思想。但是，一个工具的真正价值在于它能用来做什么。一个钟表匠的技艺高超与否，最终要看他能否造出一块精准报时的钟表。
现在，我们将开启一段新的旅程，去探索有限元方法在广阔的电磁学世界中所扮演的角色。我们将不再仅仅关注于“如何计算”，而是转向“能够解决什么”。你会发现，这个在上一章看起来有些抽象的数学工具，实际上是我们理解和改造物理世界的一把钥匙。它连接着我们日常生活的方方面面，从你口袋里的手机到医院里的精密仪器，从地球深处的磁场起源到探索宇宙的射电望远镜。让我们一同踏上这段发现之旅，看看有限元方法如何将抽象的理论锻造成工程师手中实在的工具，并点亮物理学前沿的探索之路。
现代科技世界是由无数精心设计的电子元器件构成的。这些元器件性能的优劣，直接决定了从超级计算机到智能手机等一切设备的成败。有限元方法在这里扮演着一个基础而关键的角色——它是一位虚拟的“工匠”，在产品被制造出来之前，就在计算机中精确地预测、分析和优化其电磁特性。
想象一下你电脑主板上那些密如蛛网的微带线（microstrip line），它们是信息高速奔跑的“赛道”。信号的传输速度和质量，严重依赖于这些传输线的电容和电感。如果电容计算不准，信号可能会延迟或失真，导致整个系统崩溃。特别是当传输线经过不同材料区域时，其电容分布会变得非常复杂，难以用简单的公式计算。有限元方法通过将传输线横截面划分为一个个微小的三角形单元，可以精确地计算出每个单元对总电容的贡献，即使介质分布再不规则也不在话下。同样，对于现代电路中用于控制信号完整性的屏蔽电缆，其单位长度的自感是多少？有限元方法可以通过计算电流流过导体时周围空间存储的总磁场能量，反推出这一关键参数。
在微观世界里，有限元方法同样大显身手。微机电系统（MEMS）是一种将机械结构、传感器和电子器件集成在微米尺度芯片上的技术。其中一种常见的驱动器——梳齿驱动器（comb-drive actuator）——就是利用静电力来产生微小位移。当电压施加在两组相互交错的“梳齿”电极上时，静电力会拉动其中一组活动电极，而这个力的大小恰恰取决于系统电容随位移的变化率。有限元方法能够精确模拟梳齿间复杂的电场分布，并计算出电容如何随活动电极的移动而改变，从而为这类微型驱动器的设计提供了根本性的指导。
除了让设备“动起来”，电磁屏蔽——即保护敏感设备免受外界电磁场干扰——同样至关重要。例如，医院里的核磁共振成像（MRI）设备对外界磁场极其敏感，任何微小的干扰都可能导致图像模糊，影响诊断。工程师们会使用高磁导率的材料（例如“坡莫合金”）制成磁屏蔽室或屏蔽罩。有限元方法能够模拟一个屏蔽罩在外部磁场作用下的行为，精确地计算出磁场线是如何被“引导”进屏蔽材料内部，从而使得罩内区域的磁场变得极其微弱，达到保护敏感设备的目的。
当电磁现象不再是静止的场，而是以波的形式在空间中传播时，问题变得更加动态和复杂。天线、雷达、光纤——这些与波和辐射息息相关的技术，正是有限元方法施展才华的另一片广阔天地。
一个核心问题是“开放区域”问题。天线向四面八方辐射能量，或者雷达波从远方射来，被目标散射。我们的计算区域不可能无限大，必须在某个地方“截断”。如何让电磁波在抵达计算边界时，能像在真实无垠空间中那样自由离去，而不是被生硬地反射回来，干扰计算结果？这是计算电磁学中的一个圣杯。早期的解决方案是所谓的“吸收边界条件”（Absorbing Boundary Conditions, ABCs），它在边界上模拟一个向外的行波，虽然不完美，但在很多情况下足够有效。例如，在计算一个物体的雷达散射截面（Radar Cross-Section, RCS）——通俗讲就是它在雷达上看起来有多“亮”——时，我们就可以结合有限元和吸收边界条件，来模拟入射的雷达波是如何被飞机或导弹等目标散射的。
更进一步，物理学家和工程师们发明了一种更为精妙的“计算吸波材料”——完美匹配层（Perfectly Matched Layer, PML）。PML是一种人为设计的、具有特殊复数介电常数和磁导率的虚拟材料层，它可以让来自任何角度、任何频率的电磁波在进入该层后迅速衰减，且几乎不产生任何反射。它就像一个电磁波的“黑洞”，为我们在有限的计算空间里模拟无限宇宙提供了一个近乎完美的解决方案。借助PML，我们可以用有限元方法极其精确地模拟介质天线的辐射方向图，优化其性能，让信号能更集中地发射到我们想要的方向。
除了向外辐射，我们还关心如何将波约束在特定路径内传输，就像用水管输水一样。光纤和微波波导就是电磁波的“管道”。对于一个给定横截面形状的波导，只有特定模式（modes）的电磁波才能在其中稳定传播，而其他模式则会迅速衰减。找到这些传播模式及其对应的传播常数，是设计光纤通信和微波电路的基础。当波导的截面形状复杂（例如，光子晶体光纤）时，解析求解变得不可能。此时，有限元方法再次展现威力，它将亥姆霍兹方程转化为一个矩阵特征值问题。这个特征值问题的美妙之处在于，它的解——特征值和特征向量——直接对应着波导中允许传播的模式及其场分布。求解这个矩阵，就如同找到了光波在这条“高速公路”上行驶的所有“合法车道”。
世界的奇妙之处在于各种物理现象并非孤立存在，而是相互耦合、相互影响。有限元方法的真正威力也体现在它能够构建一个统一的框架，模拟这种“多物理场”的耦合现象，并跨越学科的边界，解决其他领域的问题。
一个绝佳的例子是电热耦合。考虑一个我们日常生活中最常见的元件——保险丝。它的工作原理很简单：当电流过大时，它会熔断以保护电路。这个看似简单的过程背后，却是一个典型的电热耦合问题。电流流过具有电阻的保险丝，根据焦耳定律产生热量。热量使保险丝的温度升高。温度的升高又可能改变其电阻率，进而影响电流和产热。这是一个相互作用的循环。有限元方法可以轻松地应对这种挑战：一个电学求解器计算出电流分布和焦耳热源的分布，然后这个热源分布被无缝地传递给一个热学求解器，作为其输入。热学求解器再计算出新的温度分布。如果材料属性依赖于温度，这个新的温度分布又可以反馈给电学求解器，如此迭代，直至系统达到稳态或发生瞬态变化（如熔断）。这种耦合分析能力对于设计大功率电子器件、电动机乃至核聚变反应堆都至关重要。
另一个有趣的耦合现象是时变磁场在导体中感应出的涡旋电流，即涡流（eddy currents）。当你在电磁炉上烧水时，正是炉盘内线圈产生的高频磁场在金属锅底感应出强大的涡流，涡流的焦耳热效应将锅加热。在大型游乐场的跳楼机或过山车上，一种安全可靠的刹车方式就是让强大的永磁体扫过一块固定的铜板，在铜板中感应出的涡流会产生一个与运动方向相反的洛伦兹力，从而实现平稳减速。有限元方法通过求解包含电导率和频率的复数形式亥姆霍兹方程，能够精确地模拟涡流的分布和大小，从而帮助我们设计高效的感应加热设备或磁力刹车系统。
有限元方法的触角甚至延伸到了生命科学领域。在生物医学工程中，有一种称为“生物电阻抗分析”（Bioelectrical Impedance Analysis）的技术，通过向人体（例如手臂）施加微弱的安全交流电，并测量其表面的电压分布，来推断身体的组成成分（如脂肪和肌肉的比例）或检测病变（如肿瘤）。因为不同生物组织的电导率不同，电流在体内的分布路径也会不同。我们可以将人体组织近似为具有特定电学属性的介质，然后利用有限元方法建立模型，通过模拟电流在复杂组织结构中的流动，帮助医生从表面测量数据中解读出深层的生理信息。在这里，麦克斯韦方程组成为了探索生命奥秘的工具。
到目前为止，我们看到的有限元方法主要扮演着一个“分析师”的角色：给定一个设计，它能告诉我们其性能如何。但在科学和工程的前沿，它的角色正在发生深刻的变革——它不仅能分析，还能“创造”和“发现”。
现实世界远非线性。例如，用于制造电机和变压器铁芯的磁性材料，其磁导率并不是一个常数，而是会随着磁场的增强而减小，这种现象称为“磁饱和”。这种非线性行为使得传统的线性分析方法完全失效。然而，有限元方法框架具有极好的扩展性。我们可以将它与牛顿-拉夫逊等迭代算法相结合，来求解这些非线性方程。每一次迭代，都根据当前的场强更新材料的磁导率，然后重新求解，直到找到一个自洽的解。这使得精确设计高效能的电机和防止变压器因饱和而过热成为可能。
更具革命性的是，我们可以将有限元分析“反过来”用。传统的做法是工程师先构思一个设计，然后用有限元分析其性能。而在一个称为“拓扑优化”的领域，我们只告诉计算机目标是什么——比如，设计一个天线，使其在某个方向上的增益最大，而在其他方向上尽可能小——以及可用的材料和设计空间。然后，计算机会在一个巨大的设计可能性空间中进行搜索，自动地“生长”出最优的材料分布形态。这种方法常常会产生一些人类工程师意想不到、但性能卓越的仿生或异形结构。例如，我们可以用这种思想来自动设计出最佳的八木天线（Yagi-Uda antenna）的引向器和反射器布局。这标志着从“计算机辅助分析”到“计算机驱动设计”的飞跃。
有限元方法也让我们能够探索一些曾经只存在于理论物理学家想象中的奇异现象。在20世纪60年代，苏联物理学家Victor Veselago预言，如果一种材料同时具有负的介电常数（）和负的磁导率（），那么光在这种材料中将会发生许多奇特的行为，比如负折射。一个由这种材料制成的平板透镜，原则上可以实现“完美成像”，即突破传统光学中的衍射极限，将一个点光源发出的所有光（包括会衰逝的近场波）重新聚焦成一个完美的点。在几十年后，随着“超材料”（metamaterials）的出现，这种想法不再是天方夜谭。而有限元方法，由于其对材料属性的普遍适应性，可以毫不费力地模拟电磁波在这些奇异材料中的传播，亲眼“看”到完美透镜是如何工作的。
最令人惊叹的是，有限元方法有时甚至能揭示出物理学最深层次的统一性。考虑一个带电物体以接近光速的速度（即相对论速度）匀速运动。根据狭义相对论，运动物体在运动方向上会发生长度收缩。那么，它产生的电场和磁场会是怎样的呢？这个问题可以通过求解一组耦合的势方程来描述。当我们用有限元方法来求解这些方程时，一个惊人的结果自然而然地出现了：在推导单元刚度矩阵的过程中，代表空间各向异性的矩阵 中，沿着运动方向的分量自然地包含了因子 ，其中 。这正是洛伦兹变换的核心因子！这意味着，我们仅仅是忠实地遵循麦克斯韦方程组和有限元方法的数学规则，就不可避免地与狭义相对论的深刻内涵迎头相遇。数值方法本身，以一种无可辩驳的方式，向我们展示了时空结构的奥秘。
最后，让我们将目光投向更宏大的尺度——我们的地球和太阳。它们为何拥有强大的磁场？这个问题的答案被认为是“发电机理论”（dynamo theory）：导电流体（地核中的液态铁或太阳中的等离子体）的复杂运动，能够像发电机一样自我维持并放大磁场。这是一个极其复杂的“运动诱导”问题，涉及到流体力学和电磁学的深度耦合。有限元方法，连同其他数值技术，是科学家们用来模拟这一过程、试图解开行星和恒星磁场起源之谜的有力武器之一。他们构建整个球体的模型，设定可能的流体运动模式，然后观察一个初始的微弱磁场是否能够被放大和维持，从而与天文观测数据进行比对。
从设计一个微小的芯片，到探索宇宙尺度的奥秘，有限元方法如同一座桥梁，连接着纯粹的数学理论与纷繁复杂的现实世界。它不仅是工程师手中精密的绘图笔，更是物理学家探索未知的望远镜。通过它，我们不仅能够回答“是什么”和“为什么”，甚至开始有能力去回答“什么才是最好的”和“还有什么可能性”。这正是科学与工程相结合时，最激动人心的篇章。
在我们能够进行计算之前，我们必须正确地构建问题。这个练习将重点放在关键的第一步：将一个静磁装置的物理描述转化为适用于每个区域的正确偏微分方程。这确保了我们的数值模拟是基于坚实的物理定律，并为有限元分析奠定了基础。
问题: 一位工程师的任务是使用有限元法 (FEM) 对一个长电磁执行器中的磁场分布进行建模。该执行器在 平面内的横截面由三个不同的区域组成：一个中心铁芯 ()，铁芯内一个非同心的圆柱形气隙 ()，以及围绕铁芯的铜线圈区域 ()。铁芯由一种具有恒定相对磁导率 的软磁材料制成，因此其磁导率为 。气隙和铜线圈区域的磁导率假定为真空磁导率 。线圈承载着稳定的均匀直流电 (DC)，在区域 内产生轴向电流密度 ，而在其他地方电流密度为零。
由于系统的长度很长，该问题可以被视为二维问题。磁场由电流 产生，需要求解的合适的势是磁矢量势的 分量, 。整个执行器在一个大的计算边界框内进行建模，在该边界框的外边界 上，假定磁矢量势为零，即 。
为了设置仿真，工程师必须为三个区域（、、）中的每一个区域定义 的控制偏微分方程。下列哪组方程正确描述了需要通过 FEM 软件求解的物理问题？
A. 问题在所有区域均由拉普拉斯方程描述：
B. 问题由一个泊松方程描述，其中源项 存在于所有区域：
C. 问题由一组方程描述，其中源项 仅在线圈区域非零：
D. 问题由一个不正确的公式描述，该公式不当地使用了磁导率而非磁阻率：
E. 问题由单一形式的泊松方程 描述，并应用于所有区域：
有限元法的核心在于将一个大型的方程组分解为许多简单的、局部的“单元”矩阵，然后再将它们组装起来。这个练习提供了一个亲手计算准静态问题中一个典型单元矩阵项的机会。通过这个过程，你将理解材料属性、几何形状和工作频率是如何被编码进模型的数值核心中。
问题: 在采用有限元法（FEM）对导体中的涡流进行数值建模时，一项关键的计算任务是评估对导体横截面进行网格划分的微小离散三角形单元上的耦合积分。这些积分构成了单元特征矩阵的元。
考虑一个这样的线性三角形单元，它由三个顶点处的节点定义，坐标分别为：节点1位于 ，节点2位于 ，节点3位于 。该单元内的一个物理量使用线性基函数（也称为形函数）、和来近似。基函数在节点处的值为1，在另外两个节点处的值为0。
你的任务是计算此单元上基函数和之间的复值耦合积分。该积分定义为： 其中是磁导率，是电导率，是角频率，，积分区域为三角形单元的面积。
使用以下非磁性铜导体的参数来计算的数值：
你的答案应以包含实部和虚部的行矩阵形式给出，格式为。数值必须以SI单位表示，并四舍五入到四位有效数字。最终答案中不要包含单位。
对电势场或磁场的求解通常只是第一步，真正的工程价值来自于从这些场中推导出如力或转矩等物理量。这个练习演示了如何对静电模拟的结果进行“后处理”来计算作用在介电元件上的力。这是设计执行器、传感器和其他机电设备时至关重要的一步。
问题: 在静电学中，不仅非均匀电场可以对电介质材料施加力，当材料本身具有空间变化的介电常数时，即使在均匀场中也会受到力的作用。这种力的物理来源由开尔文力密度描述，，其中 是电场， 是位置 处的介电常数。
考虑一个使用有限元法（FEM）建模的二维静电学问题。问题域被离散化为线性三角形单元。我们关注其中一个单元，其厚度为均匀的 。该三角形的顶点（节点）标记为 1、2 和 3，分别位于坐标 、 和 处。三角形面的面积记为 。通过先前的有限元全局求解得到的这些节点上的电势为 和 。
在该特定单元内，介电常数梯度被假定为一个恒定矢量，由 给出。使用线性单元的一个关键特征是，近似电场 在整个单元内是恒定的。作用在该单元上的合力矢量 是通过对力密度在单元体积上积分得到的，并且可以证明它与介电常数梯度矢量 成正比，即 。
您的任务是推导标量比例因子 的符号表达式。请用节点电势 、节点坐标 、单元面积 及其厚度 来表示您的答案。