Finite-Integral-Methode

Finite-Integral-Methode

Die Finite-Integral-Methode basiert auf der Finite Integration Theorie (FIT) und ist ein numerisches Simulationsverfahren zur näherungsfreien Lösung der elektromagnetischen Grundgleichungen nach Maxwell. Sie bildet die mathematische Grundlage von Simulationsprogrammen für elektromagnetische Probleme wie z.B. MAFIA und CST MICROWAVE STUDIO®.

Grundlagen

Die Finite-Integral-Methode, zuerst 1976 von Thomas Weiland vorgestellt, löst die nach Maxwell benannten elektromagnetischen Grundgleichungen


\oint_{C(A)}  \vec{E} \cdot d \vec{s} \;=\; - \int_A \frac{\partial}{\partial t}\vec{B} \cdot d \vec{A}

\oint_{C(A)}  \vec{H} \cdot d \vec{s} \;=\; \int_A \left( \frac{\partial}{\partial t}\vec{D} + \vec{J}\right)\cdot d \vec{A}

\oint_{A(V)}  \vec{D} \cdot d \vec{A} \;=\; \int_V \rho  \cdot dV

 \oint_{A(V)}  \vec{B} \cdot d \vec{A} \;=\; 0

näherungsfrei in deren Integralform und mittels Integralapproximation die Materialgleichungen


\vec{D} \;=\; \varepsilon_0 \varepsilon_r \vec{E}

\vec{B} \;=\; \mu_0 \mu_r \vec{H}

\vec{J} \;=\; \kappa \vec{E} + \vec{J_s}

in diskretisierter Form.

Vorgehensweise

Das gesamte Problemgebiet wird in ein erstes (oder primäres) dreidimensionales Netz von einzelnen, möglichst kleinen Gitterzellen (engl.: mesh cells) mit den Materialeigenschaften εrr unterteilt, die jede bzgl. ihrer elektrischen Kantenspannung e_i = E_i \cdot a und ihres magnetischen Fluss durch die Randflächen b_j = B_j \cdot A_j berechnet wird.

Zusätzlich wird ein orthogonal zum ersten Gitternetz angesetztes, zweites (duales) Gitterzellennetz bzgl. der magnetischen Kantenspannung h_i = H_i \cdot a und des elektrischen Fluss durch die Randflächen d_j = D_j \cdot A_j unter Berücksichtigung der Stetigkeitsbedingungen berechnet.

Kurvenintegral über die elektrische Feldstärke in einer Gitterzelle

Durch die quaderförmige Form der Gitterzellen vereinfacht sich das Konturintegral der elektrischen Feldstärke \oint_{C(A)}{\vec{E}\cdot d\vec{s}}=\int_{{{C}_{1}}}{\vec{E}\cdot d\vec{s}}+\int_{{{C}_{2}}}{\vec{E}\cdot d\vec{s}}-\int_{{{C}_{3}}}{\vec{E}\cdot d\vec{s}}-\int_{{{C}_{4}}}{\vec{E}\cdot d\vec{s}} zur Summe \sum_{i=1}^4 e_i der Kantenspannungen einer Quaderwand der Gitterzelle.

Die zeitliche Ableitung des magnetischen Flusses durch die Randfläche der Gitterzelle wird nun dieser Summe gleichgesetzt, so dass sich folgende Gleichung ergibt:


\sum_{i=1}^4 e_i \;=\; e_1 + e_2 - e_3 - e_4 \;=\; - \frac{\partial}{\partial t} b_n

Diese Berechnung muss für alle sechs Randflächen einer Gitterzelle wiederholt werden. In Matrixschreibweise ergibt sich das Gleichungssystem


{{\begin{pmatrix}
\ldots & \ldots & \ldots & \ldots & \ldots & \ldots & \ldots & \ldots \\
\ldots & 1 & \ldots & 1 & \ldots & -1 & \ldots & -1 \\
\ldots & \ldots & \ldots & \ldots & \ldots & \ldots & \ldots &
\ldots
\end{pmatrix}}}
\cdot {{\begin{pmatrix} e_i\\
\ldots \\ e_j \\ \ldots \\ e_k \\ \ldots \\ e_l\end{pmatrix}}}
\;=\; - \frac{\partial}{\partial t} \; {{\begin{pmatrix} \ldots \\b_n
\\ \ldots
\end{pmatrix}}}
\textbf{C} \cdot \vec{e} \;=\; - \frac{\partial}{\partial t} \;\vec{b}

Die beschreibende Matrix \textbf{C} besitzt als Elemente nur die Werte 1, 0, -1.

Analog dazu werden die übrigen Maxwell’schen Gleichungen behandelt. In Matrixschreibweise ergibt sich das Gitter-Maxwell-Gleichungssystem


\textbf{C} \cdot \vec{e} \;=\; - \frac{\partial}{\partial t} \;\vec{b}

\textbf{C}_{\text{Dual}} \cdot \vec{h} \;=\; \frac{\partial}{\partial t} \; \vec{d} + \vec{j}


\textbf{S}_{\text{Dual}} \cdot \vec{d} \;=\; \vec{q}


\textbf{S} \cdot \vec{b} \;=\; \vec{0}

Die Matrix \textbf{C} entspricht dem analytischen Rotations-Operator, die Matrix \textbf{S} entspricht dem analytischen Divergenz-Operator. Der Index Dual weist auf die Berechnung der Kantenspannungen und Flüsse im Dualen Gitter hin.

Die Materialgleichungen werden analog zu den Maxwell’schen Gleichungen diskretisiert.


\vec{d} \;=\; \textbf{M}_{\varepsilon} \;\cdot\; \vec{e}

\vec{b} \;=\; \textbf{M}_{\mu} \;\cdot\; \vec{h}

\vec{j} \;=\; \textbf{M}_{\kappa} \;\cdot\; \vec{e} \;+\; \vec{j}_s,

wobei die Materialgrößen orts-, frequenz- und richtungsabhängig sein können.

Die FIT-Methode ist auf alle elektromagnetische Probleme im Zeit- und Frequenzbereich anwendbar, sowohl in der Elektrostatik, als auch in der Elektrodynamik. Durch den speziellen Zuschnitt der FIT-Methode auf die Maxwell’schen Gleichungen und das daraus entstehende diskrete Analogon sind die Stetigkeitsbedingungen a priori erfüllt und die analytischen Eigenschaften der Vektoroperationen werden beibehalten.

Für elektrodynamische Probleme werden im Frequenzbereich alle zeitlichen Ableitungen durch jω ersetzt. Das Ergebnis einer Simulation im Frequenzbereich liefert die Impulsantwort auf ein monofrequentes Eingangssignal.

Im Zeitbereich ist eine breitbandige Anregung mit freien Signalverläufen gestattet. Die Simulationsrechnung beschreibt in diesem Fall das Frequenzverhalten über einen vorab definierten Frequenzbereich.

Literatur

  • T. Weiland: Eine Methode zur Lösung der Maxwellschen Gleichungen für sechskomponentige Felder auf diskreter Basis, AEÜ, Band 31, Heft 3, pp. 116–120, 1977
  • T. Weiland: A Discretization Method for the Solution of Maxwell’s Equations for Six-Component Fields, Electronics and Communications AEUE, vol. 31, no. 3, pp. 116–120, 1977.

Wikimedia Foundation.

Игры ⚽ Поможем написать курсовую

Schlagen Sie auch in anderen Wörterbüchern nach:

  • Finite-Elemente-Methode — Die Finite Elemente Methode (FEM), auch „Methode der finiten Elemente“ genannt, ist ein numerisches Verfahren zur Lösung von partiellen Differentialgleichungen. Sie ist ein weit verbreitetes modernes Berechnungsverfahren im Ingenieurwesen und ist …   Deutsch Wikipedia

  • Finite-Elemente-Analyse — Die Finite Elemente Methode (FEM) ist ein numerisches Verfahren zur näherungsweisen Lösung, insbesondere elliptischer partieller Differentialgleichungen mit Randbedingungen. Sie ist auch ein weit verbreitetes modernes Berechnungsverfahren im… …   Deutsch Wikipedia

  • Finite-Elemente-Verfahren — Die Finite Elemente Methode (FEM) ist ein numerisches Verfahren zur näherungsweisen Lösung, insbesondere elliptischer partieller Differentialgleichungen mit Randbedingungen. Sie ist auch ein weit verbreitetes modernes Berechnungsverfahren im… …   Deutsch Wikipedia

  • Finite Elemente — Die Finite Elemente Methode (FEM) ist ein numerisches Verfahren zur näherungsweisen Lösung, insbesondere elliptischer partieller Differentialgleichungen mit Randbedingungen. Sie ist auch ein weit verbreitetes modernes Berechnungsverfahren im… …   Deutsch Wikipedia

  • Finite Elementemethode — Die Finite Elemente Methode (FEM) ist ein numerisches Verfahren zur näherungsweisen Lösung, insbesondere elliptischer partieller Differentialgleichungen mit Randbedingungen. Sie ist auch ein weit verbreitetes modernes Berechnungsverfahren im… …   Deutsch Wikipedia

  • Methode der finiten Elemente — Die Finite Elemente Methode (FEM) ist ein numerisches Verfahren zur näherungsweisen Lösung, insbesondere elliptischer partieller Differentialgleichungen mit Randbedingungen. Sie ist auch ein weit verbreitetes modernes Berechnungsverfahren im… …   Deutsch Wikipedia

  • Finite-Volumen-Verfahren — Das Finite Volumen Verfahren ist ein numerisches Verfahren zur Lösung von Erhaltungsgleichungen, also von speziellen, häufig hyperbolischen, partiellen Differentialgleichungen, denen ein Erhaltungssatz zugrunde liegt. Am prominentesten ist der… …   Deutsch Wikipedia

  • Finite Quantenfeldtheorie — Die finite Quantenfeldtheorie (FQFT) ist ein Versuch, mit den klassischen Schwierigkeiten der Quantenfeldtheorie (QFT) fertig zu werden. Eine dieser klassischen Schwierigkeiten ist die UV Katastrophe, die in der klassischen Theorie durch eine… …   Deutsch Wikipedia

  • Galerkin-Methode — Die Galerkin Methode (auch Galerkin Ansatz, nach Boris Galerkin) ist ein numerisches Verfahren zur näherungsweisen Lösung von partiellen Differentialgleichungen. Sie stellt die gebräuchlichste Variante der Methode der gewichteten Residuen dar,… …   Deutsch Wikipedia

  • FEM (Finite Element Method) — Methode der finiten Elemente, i. d. R. auf Verschiebungsgrößenverfahren aufbauende Methode unter Verwendung von Näherungsansätzen für den Verlauf der Verschiebungsgrößen längs des finiten Elements; Unbekannte sind demzufolge Verschiebungsgrößen;… …   Erläuterung wichtiger Begriffe des Bauwesens

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”