Ansys Maxwell2——二维电磁场理论和有限元基础

对于偏微分方程,使其解成为唯一的辅助条件可分为两种:一种是表达场的边界所处的物理情况、称为边界条件;另一种是确定场的初始状态,称为初始条件。边界条件和初始条件合称为定解条件。

2.1 二维电磁场基本理论

电磁场的经典描述是麦克斯韦方程组,电机电磁场分析—般采用位函数表示,位函数

比场量本身更容易建立边界条件。位函数包括磁矢位 $A$ 和磁标位$\varphi $。

2.1.1 麦克斯韦方程

麦克斯韦方程组实际上由似个定律组成,它们分别为:安培环路定律、法拉第电磁感应定律、高斯电通定律和高斯磁通定律(亦称磁通连续定律)。

\[\left\{ \begin{array}{l}
\nabla  \times \vec H = \vec J + \frac{{\partial \vec D}}{{\partial t}}\\
\nabla  \times \vec E =  - \frac{{\partial \vec B}}{{\partial t}}\\
\nabla  \cdot \vec D = \rho \\
\nabla  \cdot \vec B = 0
\end{array} \right.\]

对于线性媒质:

\[\left\{ \begin{array}{l}
\vec D = \varepsilon \vec E\\
\vec B = \mu \vec H\\
\vec J = \sigma \vec E
\end{array} \right.\]

2.1.2 位函数及其微分方程

在无旋场(即旋度为零的场)中可以采用标量位函数,而在有旋场中,则必须用矢量位函数。

因此可以引入标量电位或标量磁位。它们与场强的关系是:

\[\left\{ \begin{array}{l}
\vec E =  - \nabla \varphi \\
\vec H =  - \nabla {\varphi _m}
\end{array} \right.\]

Ansoft 中常用的求解方程有:

二维、三维静电场求解器所满足的泊松方程

\[\nabla \varphi  =  - \frac{\rho }{\varepsilon }\]

二维稳恒电场求解器所满足的拉普拉斯方程

\[\nabla \varphi  = 0\]

二维交变电场求解器所满足的复数拉普拉斯方程

二维静磁场求解器所满足的非齐次标量波动方程

二维涡流场求解器所满足的波动方程

二维轴向磁场求解器所满足的齐次波动方程

三维静磁场和涡流场求解器所满足的齐次波动方程

2.1.3 电磁场中的边界条件

电磁场求解过程中有各种各样的边界条件,具体包括以下几类:

狄利克莱边界条件

\[{\left. \varphi  \right|_\Gamma } = g(\Gamma )\]

诺依曼边界条件

\[\left. {\frac{{\partial \varphi }}{{\partial n}}} \right|\Gamma  + {\left. {f(\Gamma )\varphi } \right|_\Gamma } = h(\Gamma )\]

自然边界条件

对称边界条件

周期边界条件

气球边界条件

阻抗边界条件

2.2 二维有限元理论初步

2.2.1 二维有限元法

将有限元法的过程简要地归纳成如下几个步骤:

Step1 列出与偏微分方程边值问题等价的条件变分问题。

Step 2 将区域作三角形单元剖分,并在单元中,构造出线性插值函数。

Step 3 将能量泛函的极值问题转化为能量函数的极值问题,建立线性代数方程组。

Step 4 求解线性代数方程组。

2.2.2 电磁场求解后处理

通过有限元法求解出的节点的电势或者磁势值是远远不够的,在实际的问题当中我们还需要得到许多其它物理量,如磁感应强度、电位移通量、电磁场能量、电磁力和力矩、电感和电容等。这些量是比较容易通过求得的电势和磁势量导出的,这种导出过程称之为有限元解的后处理。

电场储能:

\[\begin{array}{*{20}{l}}
{{W_e} = \frac{1}{2}\int_\Omega  {\vec D \cdot \vec E{\rm{d}}\Omega }  = \frac{\varepsilon }{2}\int_\Omega  {{{\left| {\nabla \varphi } \right|}^2}{\rm{d}}\Omega } }\\
{\;\;\;\; = \frac{\varepsilon }{2}\sum\limits_{e = 1}^n {\int_{{\Omega ^e}} {[{{(\frac{{\partial \varphi }}{{\partial x}})}^2} + {{(\frac{{\partial \varphi }}{{\partial y}})}^2}]} } {\rm{d}}\Omega }
\end{array}\]

磁场储能

\[\begin{array}{*{20}{l}}
{{W_m} = \frac{1}{2}\int_\Omega  {\vec B \cdot \vec H{\rm{d}}\Omega }  = \frac{1}{{2\mu }}\int_\Omega  {{{\left| {\vec B} \right|}^2}{\rm{d}}\Omega } }\\
{\;\;\;\; = \frac{1}{{2\mu }}\sum\limits_{e = 1}^n {\int_{{\Omega ^e}} {[{{(\frac{{\partial A}}{{\partial x}})}^2} + {{(\frac{{\partial A}}{{\partial y}})}^2}]} } {\rm{d}}\Omega }
\end{array}\]

电感

\[L = \frac{{2{W_m}}}{{{I^2}}}\]

电容

\[C = \frac{{2{W_e}}}{{{V^2}}}\]

电磁力\[F = \frac{{\partial {W_m}^\prime }}{{\partial S}}\]

力矩\[{T_m} = P\frac{{\partial {W_m}^\prime }}{{\partial \theta }}\]

posted @ 2020-09-18 09:47  叮叮当当sunny  阅读(3058)  评论(0编辑  收藏  举报