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 }}\]


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM