aboutsummaryrefslogtreecommitdiffstats
path: root/buch/papers/kreismembran
diff options
context:
space:
mode:
Diffstat (limited to '')
-rw-r--r--buch/papers/kreismembran/main.tex26
-rw-r--r--buch/papers/kreismembran/references.bib51
-rw-r--r--buch/papers/kreismembran/teil0.tex88
-rw-r--r--buch/papers/kreismembran/teil1.tex164
-rw-r--r--buch/papers/kreismembran/teil2.tex146
-rw-r--r--buch/papers/kreismembran/teil3.tex108
-rw-r--r--buch/papers/kreismembran/teil4.tex16
7 files changed, 446 insertions, 153 deletions
diff --git a/buch/papers/kreismembran/main.tex b/buch/papers/kreismembran/main.tex
index 67b436c..f6000a1 100644
--- a/buch/papers/kreismembran/main.tex
+++ b/buch/papers/kreismembran/main.tex
@@ -3,34 +3,16 @@
%
% (c) 2020 Hochschule Rapperswil
%
-\chapter{Thema\label{chapter:kreismembran}}
-\lhead{Thema}
+\chapter{Schwingungen einer kreisförmigen Membran\label{chapter:kreismembran}}
+\lhead{Schwingungen einer kreisförmigen Membran}
\begin{refsection}
-\chapterauthor{Hans Muster}
-
-Ein paar Hinweise für die korrekte Formatierung des Textes
-\begin{itemize}
-\item
-Absätze werden gebildet, indem man eine Leerzeile einfügt.
-Die Verwendung von \verb+\\+ ist nur in Tabellen und Arrays gestattet.
-\item
-Die explizite Platzierung von Bildern ist nicht erlaubt, entsprechende
-Optionen werden gelöscht.
-Verwenden Sie Labels und Verweise, um auf Bilder hinzuweisen.
-\item
-Beginnen Sie jeden Satz auf einer neuen Zeile.
-Damit ermöglichen Sie dem Versionsverwaltungssysteme, Änderungen
-in verschiedenen Sätzen von verschiedenen Autoren ohne Konflikt
-anzuwenden.
-\item
-Bilden Sie auch für Formeln kurze Zeilen, einerseits der besseren
-Übersicht wegen, aber auch um GIT die Arbeit zu erleichtern.
-\end{itemize}
+\chapterauthor{Andrea Mozzini Vellen und Tim Tönz}
\input{papers/kreismembran/teil0.tex}
\input{papers/kreismembran/teil1.tex}
\input{papers/kreismembran/teil2.tex}
\input{papers/kreismembran/teil3.tex}
+\input{papers/kreismembran/teil4.tex}
\printbibliography[heading=subbibliography]
\end{refsection}
diff --git a/buch/papers/kreismembran/references.bib b/buch/papers/kreismembran/references.bib
index 0b6a683..acf8f90 100644
--- a/buch/papers/kreismembran/references.bib
+++ b/buch/papers/kreismembran/references.bib
@@ -4,6 +4,25 @@
% (c) 2020 Autor, Hochschule Rapperswil
%
+@online{kreismembran:Duden:Membran,
+ title = {Duden:Membran},
+ url = {https://www.duden.de/rechtschreibung/Membran},
+ date = {2022-07-20},
+ year = {2022},
+ month = {7},
+ day = {20}
+}
+
+@online{kreismembran:wellengleichung_herleitung,
+ title = {Derivation of the 2D Wave Equation},
+ author = {Dr. Christopher Lum},
+ url = {https://www.youtube.com/watch?v=KAS7JBztw8E&t=0s},
+ date = {2022-07-20},
+ year = {2022},
+ month = {7},
+ day = {20}
+}
+
@online{kreismembran:bibtex,
title = {BibTeX},
url = {https://de.wikipedia.org/wiki/BibTeX},
@@ -24,7 +43,7 @@
}
@article{kreismembran:mendezmueller,
- author = { Tabea Méndez and Andreas Müller },
+ author = { Tabea Méndez and Andreas Müller },
title = { Noncommutative harmonic analysis and image registration },
journal = { Appl. Comput. Harmon. Anal.},
year = 2019,
@@ -33,3 +52,33 @@
url = {https://doi.org/10.1016/j.acha.2017.11.004}
}
+@book{lokenath_debnath_integral_2015,
+ edition = {Third Edition},
+ title = {Integral Tansforms and Their Applications},
+ publisher = {{CRC} Press},
+ author = {{Lokenath Debnath} and Dambaru Bhatta},
+ date = {2015},
+}
+
+@thesis{nishanth_p_vibrations_2018,
+ title = {Vibrations of a Circular Membrane - Some Undergraduadte Exercises},
+ type = {phdthesis},
+ author = {{Nishanth P.} and {Udayanandan K. M.}},
+ date = {2018},
+}
+
+@thesis{prof_dr_horst_knorrer_kreisformige_2013,
+ title = {Kreisförmige Membranen},
+ institution = {{ETHZ}},
+ type = {phdthesis},
+ author = {{Prof. Dr. Horst Knörrer}},
+ date = {2013},
+}
+
+@thesis{kreismembran:membrane_vs_thin_plate,
+ title = {Modeling and Control of SPIDER Satellite Components},
+ institution = {{faculty of the Virginia Polytechnic Institute and State University}},
+ type = {Dissertation},
+ author = {{Eric John Ruggiero Doctor of Philosophy In Mechanical Engineering}},
+ date = {2005},
+} \ No newline at end of file
diff --git a/buch/papers/kreismembran/teil0.tex b/buch/papers/kreismembran/teil0.tex
index e4b1711..bb8188d 100644
--- a/buch/papers/kreismembran/teil0.tex
+++ b/buch/papers/kreismembran/teil0.tex
@@ -3,20 +3,80 @@
%
% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
%
-\section{Teil 0\label{kreismembran:section:teil0}}
-\rhead{Teil 0}
-Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam
-nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam
-erat, sed diam voluptua \cite{kreismembran:bibtex}.
-At vero eos et accusam et justo duo dolores et ea rebum.
-Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum
-dolor sit amet.
+\section{Einleitung\label{kreismembran:section:teil0}}
+\rhead{Membran}
+Eine Membran oder selten ein Schwingblatt ist laut Duden \cite{kreismembran:Duden:Membran} ein "dünnes Blättchen aus Metall, Papier o. Ä., das durch seine Schwingungsfähigkeit geeignet ist, Schallwellen zu übertragen ...".
+Ein dünnes Blättchen aus Metall zeig jedoch nicht die selben dynamischen Eigenschaften auf, wie ein gespanntes Stück Papier.
+Beschreibt man das dynamische Verhalten, muss zwischen einer dünnen Platte und einer Membrane unterschieden werden \cite{kreismembran:membrane_vs_thin_plate}.
+Eine dünne Platte zum Beispiel aus Metall, wirkt selbst entgegen ihrer Deformation sobald sie gekrümmt wird.
+Eine Membran auf der anderen Seite besteht aus einem Material welches sich ohne Kraftaufwand verbiegen lässt wie zum Beispiel Papier.
+Bevor Papier als schwingende Membran betrachtet werden kann wird jedoch noch eine Spannung $ T $ benötigt welche das Material daran hindert aus der Ruhelage gebracht zu werden.
-Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam
-nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam
-erat, sed diam voluptua.
-At vero eos et accusam et justo duo dolores et ea rebum. Stet clita
-kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit
-amet.
+Ein geläufiges Beispiel einer Kreismembran ist eine runde Trommel.
+Sie besteht herkömmlicher weise aus einem Leder (Fell), welches auf einen offenen Zylinder (Zargen) aufgespannt wird.
+Das Leder alleine erzeugt nach einem Aufschlag keine hörbaren Schwingungen.
+Sobald das Fell jedoch über den Zargen gespannt wird, kann das Fell auf verschiedensten weisen weiter schwingen, was für den Klang der Trommel verantwortlich ist.
+Wie genau diese Schwingungen untersucht werden können wird in der Folgenden Arbeit Diskutiert.
+\paragraph{Annahmen}
+Um die Wellengleichung herzuleiten \cite{kreismembran:wellengleichung_herleitung}, muss ein Modell einer Membran definiert werden.
+Das untersuchte Modell einer Membrane Erfüllt folgende Eigenschaften:
+\begin{enumerate}[i]
+ \item Die Membran ist homogen.
+ Dies bedeutet, dass die Membran über die ganze Fläche die selbe Dichte $ \rho $ und Elastizität hat.
+ Durch die konstante Elastizität ist die ganze Membran unter gleichmässiger Spannung $ T $.
+ \item Die Membran ist perfekt flexibel.
+ Daraus folgt, dass die Membran ohne Kraftaufwand verbogen werden kann.
+ Die Membran ist dadurch nicht allein stehend schwing-fähig, hierzu muss sie gespannt werden mit einer Kraft $ T $.
+ \item Die Membran kann sich nur in Richtung ihrer Normalen in kleinem Ausmass Auslenken.
+ Auslenkungen in der ebene der Membran sind nicht möglich.
+ \item Die Membran erfährt keine Art von Dämpfung.
+ Die Membran wird also nicht durch ihr umliegendes Medium abgebremst noch erfährt sie Wärmeverluste durch Deformation.
+ Die resultierende Schwingung wird daher nicht gedämpft sein.
+
+\end{enumerate}
+
+\subsection{Wellengleichung} Um die Wellengleichung einer Membran herzuleiten wird vorerst eine schwingende Saite betrachtet.
+Es lohnt sich das Verhalten einer Saite zu beschreiben da eine Saite das selbe Verhalten wie eine Membran aufweist mit dem Unterschied einer fehlenden Dimension.
+Die Verbindung zwischen Membran und Saite ist intuitiv ersichtlich, stellt man sich einen Querschnitt einer Trommel vor.
+%Wie analog zur Membran kann eine Saite erst unter Spannung schwingen.
+
+Abbildung \ref{TODO} ist ein infinitesimales Stück einer Saite mit Länge $ dx $ skizziert.
+Wie für die Membran ist die Annahme iii gültig, keine Bewegung in die Richtung $ \hat{x} $.
+Um dies zu erfüllen muss der Punkt $ P_1 $ gleich stark in Richtung $ -\hat{x} $ gezogen werden wie der Punkt $ P_2 $ in Richtung $ \hat{x} $ gezogen wird. Ist $ T_1 $ die Kraft welche mit Winkel $ \alpha $ auf Punkt $ P_1 $ wirkt sowie $ T_2 $ und $ \beta$ das analoge für Punkt $ P_2 $ ist, so können die Kräfte
+\begin{equation}\label{kreismembran:eq:no_translation}
+ T_1 \cos \alpha = T_2 \cos \beta = T
+\end{equation}
+gleichgesetzt werden.
+Das dynamische verhalten der senkrechten Auslenkung $ u(x,t) $ muss das newtonsche Gesetz
+\begin{equation*}
+ \sum F = m \cdot a
+\end{equation*}
+befolgen. Die senkrecht wirkenden Kräfte werden mit $ T_1 $ und $ T_2 $ ausgedrückt, die Masse als Funktion der Dichte $ \rho $ und die Beschleunigung in Form der zweiten Ableitung als
+\begin{equation*}
+ T_2 \sin \beta - T_1 \sin \alpha = \rho dx \frac{\partial^2 u}{\partial t^2} .
+\end{equation*}
+Die Gleichung wird durch $ T $ dividiert, wobei $ T $ nach \ref{kreismembran:eq:no_translation} geschickt gewählt wird. Somit kann
+\begin{equation*}
+ \frac{T_2 \sin \beta}{T_2 \cos \beta} - \frac{T_1 \sin \alpha}{T_1 \cos \alpha} = \frac{\rho dx}{T} \frac{\partial^2 u}{\partial t^2}
+\end{equation*}
+vereinfacht als
+\begin{equation*}
+ \tan \beta - \tan \alpha = \frac{\rho dx}{T} \frac{\partial^2 u}{\partial t^2}
+\end{equation*}
+geschrieben werden.
+Der $ \tan \alpha $ entspricht der örtlichen Ableitung von $ u(x,t) $ an der Stelle $ x_0 $ und analog der $ \tan \beta $ für die Stelle $ x_0 + dx $.
+Die Gleichung wird dadurch zu
+\begin{equation*}
+ \frac{\partial u}{\partial x} \big\vert_{x_0 + dx} - \frac{\partial u}{\partial x} \big\vert_{x_0} = \frac{\rho dx}{T} \frac{\partial^2 u}{\partial t^2}.
+\end{equation*}
+Durch die Division mit $ dx $ entsteht
+\begin{equation*}
+ \frac{1}{dx} \bigg[\frac{\partial u}{\partial x} \big\vert_{x_0 + dx} - \frac{\partial u}{\partial x} \big\vert_{x_0}\bigg] = \frac{\rho}{T}\frac{\partial^2 u}{\partial t^2}.
+\end{equation*}
+Auf der Linken Seite der Gleichung wird die Differenz der Steigungen durch die Intervall-Länge geteilt, in anderen Worten die zweite Ableitung von $ u(x,t) $ nach $ x $ berechnet. Der Term $ \frac{\rho}{T} $ wird mit $ c^2 $ ersetzt, da der Bruch für eine gegebene Membran eine positive Konstante sein muss. Somit resultiert die, in der Literatur gebräuchliche Form
+\begin{equation}
+ \frac{1}{c^2}\frac{\partial^2u}{\partial t^2} = \Delta u.
+\end{equation}
+In dieser Form ist die Gleichung auch gültig für eine Membran. Für den Fall einer Membran muss lediglich die Ableitung in zwei Dimensionen gerechnet werden. \ No newline at end of file
diff --git a/buch/papers/kreismembran/teil1.tex b/buch/papers/kreismembran/teil1.tex
index b715075..39ca598 100644
--- a/buch/papers/kreismembran/teil1.tex
+++ b/buch/papers/kreismembran/teil1.tex
@@ -2,54 +2,120 @@
% teil1.tex -- Beispiel-File für das Paper
%
% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
-%
-\section{Teil 1
-\label{kreismembran:section:teil1}}
-\rhead{Problemstellung}
-Sed ut perspiciatis unde omnis iste natus error sit voluptatem
-accusantium doloremque laudantium, totam rem aperiam, eaque ipsa
-quae ab illo inventore veritatis et quasi architecto beatae vitae
-dicta sunt explicabo.
-Nemo enim ipsam voluptatem quia voluptas sit aspernatur aut odit
-aut fugit, sed quia consequuntur magni dolores eos qui ratione
-voluptatem sequi nesciunt
-\begin{equation}
-\int_a^b x^2\, dx
-=
-\left[ \frac13 x^3 \right]_a^b
-=
-\frac{b^3-a^3}3.
-\label{kreismembran:equation1}
-\end{equation}
-Neque porro quisquam est, qui dolorem ipsum quia dolor sit amet,
-consectetur, adipisci velit, sed quia non numquam eius modi tempora
-incidunt ut labore et dolore magnam aliquam quaerat voluptatem.
-
-Ut enim ad minima veniam, quis nostrum exercitationem ullam corporis
-suscipit laboriosam, nisi ut aliquid ex ea commodi consequatur?
-Quis autem vel eum iure reprehenderit qui in ea voluptate velit
-esse quam nihil molestiae consequatur, vel illum qui dolorem eum
-fugiat quo voluptas nulla pariatur?
-
-\subsection{De finibus bonorum et malorum
-\label{kreismembran:subsection:finibus}}
-At vero eos et accusamus et iusto odio dignissimos ducimus qui
-blanditiis praesentium voluptatum deleniti atque corrupti quos
-dolores et quas molestias excepturi sint occaecati cupiditate non
-provident, similique sunt in culpa qui officia deserunt mollitia
-animi, id est laborum et dolorum fuga \eqref{000tempmlate:equation1}.
-
-Et harum quidem rerum facilis est et expedita distinctio
-\ref{kreismembran:section:loesung}.
-Nam libero tempore, cum soluta nobis est eligendi optio cumque nihil
-impedit quo minus id quod maxime placeat facere possimus, omnis
-voluptas assumenda est, omnis dolor repellendus
-\ref{kreismembran:section:folgerung}.
-Temporibus autem quibusdam et aut officiis debitis aut rerum
-necessitatibus saepe eveniet ut et voluptates repudiandae sint et
-molestiae non recusandae.
-Itaque earum rerum hic tenetur a sapiente delectus, ut aut reiciendis
-voluptatibus maiores alias consequatur aut perferendis doloribus
-asperiores repellat.
+\section{Lösungsmethode 1: Separationsmethode 
+ \label{kreismembran:section:teil1}}
+\rhead{Lösungsmethode 1: Separationsmethode}
+An diesem Punkt bleibt also nur noch die Lösung der partiellen Differentialgleichung. In diesem Kapitel wird sie mit Hilfe der Separationsmethode gelöst.
+
+\subsection{Aufgabestellung\label{sub:aufgabestellung}}
+Wie im vorherigen Abschnitt gezeigt, lautet die partielle Differentialgleichung, die die Schwingungen einer Membran beschreibt:
+\begin{equation*}
+ \frac{1}{c^2}\frac{\partial^2u}{\partial t^2} = \Delta u.
+\end{equation*}
+Da es sich um eine Kreisscheibe handelt, werden Polarkoordinaten verwendet, so dass sich der Laplaceoperator
+\begin{equation*}
+ \Delta
+ =
+ \frac{\partial^2}{\partial r^2}
+ +
+ \frac1r
+ \frac{\partial}{\partial r}
+ +
+ \frac{1}{r 2}
+ \frac{\partial^2}{\partial\varphi^2}
+ \label{buch:pde:kreis:laplace}
+\end{equation*}
+ergibt.
+
+Es wird eine runde elastische Membran berücksichtigt, die das Gebiet $\Omega$ abdeckt und am Rand $\Gamma$ befestigt ist.
+Es wird daher davon ausgegangen, dass die Membran aus einem homogenen Material von vernachlässigbarer Dicke gefertigt ist.
+Die Membran kann verformt werden, aber innere elastische Kräfte wirken den Verformungen entgegen. Es wirken keine äusseren Kräfte. Es handelt sich somit von einer kreisförmligen eingespannten homogenen schwingenden Membran.
+
+Daher ist die Membranabweichung im Punkt $(r,\varphi)$ $\in$ $\overline{\rm \Omega}$ zum Zeitpunkt $t$:
+\begin{align*}
+ u: \overline{\rm \Omega} \times \mathbb{R}_{\geq 0} &\longrightarrow \mathbb{R}\\
+ (r,\varphi,t) &\longmapsto u(r,\varphi,t)
+\end{align*}
+Da die Membran am Rand befestigt ist, kann es keine Schwingungen geben, so dass die \textit{Dirichlet-Randbedingung} \cite{prof_dr_horst_knorrer_kreisformige_2013}
+\begin{equation*}
+ u\big|_{\Gamma} = 0 \quad \text{für} \quad 0 \leq \varphi \leq 2\pi,\quad t \geq 0
+\end{equation*}
+gilt.
+
+Um eine eindeutige Lösung bestimmen zu können, werden die folgenden Anfangsbedingungen festgelegt:
+\begin{align*}
+ u(r,\varphi, 0) &= f(r,\varphi)\\
+ u_t(r,\varphi, 0) &= g(r,\varphi).
+\end{align*}
+
+\subsection{Lösung\label{sub:lösung1}}
+\subsubsection{Ansatz der Separation der Variablen\label{subsub:ansatz_separation}}
+Daher muss an dieser Stelle von einer Separation der Variablen ausgegangen werden:
+\begin{equation*}
+ u(r,\varphi, t) = F(r)G(\varphi)T(t)
+\end{equation*}
+Dank der Randbedingungen kann also gefordert werden, dass $F(R)=0$ ist, und natürlich, dass $G(\varphi)$ $2\pi$ periodisch ist. Eingesetzt in der Differenzialgleichung ergibt sich:
+\begin{equation*}
+ \frac{1}{c^2}\frac{T''(t)}{T(t)}=\frac{F''(r)}{F(r)}+\frac{1}{r}\frac{F'(r)}{F(r)}+\frac{1}{r^2}\frac{G''(\varphi)}{G(\varphi)}.
+\end{equation*}
+Da die linke Seite nur von $t$ und die rechte Seite nur von $r$ und $\varphi$ abhängt, müssen sie gleich einer reellen Zahl sein. Aus physikalischen Gründen suchen wir nach Lösungen, die weder exponentiell in der Zeit wachsen noch exponentiell abklingen. Dies bedeutet, dass die Konstante negativ sein muss, also schreibt man $k=-k^2$. Daraus ergeben sich die folgenden zwei Gleichungen:
+\begin{align*}
+ T''(t) + c^2\kappa^2T(t) &= 0\\
+ r^2\frac{F''(r)}{F(r)} + r \frac{F'(r)}{F(r)} +\kappa^2 r^2 &= - \frac{G''(\varphi)}{G(\varphi)}.
+\end{align*}
+In der zweiten Gleichung hängt die linke Seite nur von $r$ ab, während die rechte Seite nur von $\varphi$ abhängt. Sie müssen also wiederum gleich einer reellen Zahl $\nu$ sein. Also das:
+\begin{align*}
+ r^2F''(r) + rF'(r) + (\kappa^2 r^2 - \nu)F(r) &= 0 \\
+ G''(\varphi) &= \nu G(\varphi).
+\end{align*}
+
+\subsubsection{Lösung für $G(\varphi)$\label{subsub:lösung_G}}
+Da für die Zweite Gelichung Lösungen von Schwingungen erwartet werden, für die $G''(\varphi)=-\omega^2 G(\varphi)$ gilt, schreibt die gemeinsame Konstante als $-\nu^2$, was die Formeln später vereinfacht. Also:
+\begin{equation*}
+ G(\varphi) = C_n \cos(\varphi) + D_n \sin(\varphi)
+ \label{eq:cos_sin_überlagerung}
+\end{equation*}
+
+\subsubsection{Lösung für $F(r)$\label{subsub:lösung_F}}
+Die Gleichung für $F$ hat die Gestalt
+\begin{align}
+ r^2F''(r) + rF'(r) + (\kappa^2 r^2 - n^2)F(r) = 0
+ \label{eq:2nd_degree_PDE}
+\end{align}
+Wir bereits in Kapitel \ref{buch:differntialgleichungen:section:bessel} gezeigt, sind die Besselfunktionen
+\begin{equation*}
+ J_{\nu}(x) = r^\nu \displaystyle\sum_{m=0}^{\infty} \frac{(-1)^m x^{2m}}{2^{2m+\nu}m! \Gamma (\nu + m+1)}
+\end{equation*}
+Lösungen der Besselschen Differenzialgleichung
+\begin{equation*}
+ x^2 y'' + xy' + (x^2 - \nu^2)y = 0
+\end{equation*}
+Die Funktionen $F(r) = J_n(\kappa r)$ lösen also die Differentialgleichung \eqref{eq:2nd_degree_PDE}. Die
+Randbedingung $F(R)=0$ impliziert, dass $\kappa R$ eine Nullstelle der Besselfunktion
+$J_n$ sein muss. Man kann zeigen, dass die Besselfunktionen $J_n, n \geq 0$, alle unendlich
+viele Nullstellen
+\begin{equation*}
+ \alpha_{1n} < \alpha_{2n} < ...
+\end{equation*}
+haben, und dass $\underset{\substack{m\to\infty}}{\text{lim}} \alpha_{mn}=\infty$. Somit ergibt sich, dass $\kappa = \frac{\alpha_{mn}}{R}$ für ein $m\geq 1$, und dass
+\begin{equation*}
+ F(r) = J_n (\kappa_{mn}r) \quad \text{mit} \quad \kappa_{mn}=\frac{\alpha_{mn}}{R}
+\end{equation*}
+
+\subsubsection{Lösung für $T(t)$\label{subsub:lösung_T}}
+Die Differenzialgleichung $T''(t) + c^2\kappa^2T(t) = 0$, wird auf ähnliche Weise gelöst wie $G(\varphi)$.
+
+\subsubsection{Zusammenfassung der Lösungen\label{subsub:zusammenfassung_lösungen}}
+Durch Überlagerung aller Ergebnisse erhält man die Lösung
+\begin{align}
+ u(r, \varphi, t) = \displaystyle\sum_{m=1}^{\infty}\displaystyle\sum_{n=0}^{\infty} J_n (k_{mn}r)[a_{mn}\cos(n\varphi) + b_{mn}\sin(n\varphi)](n\varphi)[c_{mn}\cos(c \kappa_{mn} t)+d_{mn}\sin(c \kappa_{mn} t)]
+ \label{eq:lösung_endliche_generelle}
+\end{align}
+
+Dabei sind $m$ und $n$ ganze Zahlen, wobei $m$ für die Anzahl der Knotenkreise und $n$
+für die Anzahl der Knotenlinien steht. Es gibt bestimmte Bereiche auf der Membran, in denen es keine Bewegung oder Vibration gibt. Wenn der nicht schwingende Bereich ein Kreis ist, nennt man ihn einen Knotenkreis, und wenn er eine Linie ist, nennt man ihn ebenfalls eine Knotenlinie. $Jn(\kappa_{mn}r)$ ist die Besselfunktion $n$-ter Ordnung, wobei $\kappa mn$ die Wellenzahl und $r$ der Radius ist. $a_{mn}$ und $b_{mn}$ sind die zu bestimmenden Konstanten.
+
+
+An diesem Punkt stellte sich die Frage, ob es möglich wäre, die partielle Differentialgleichung mit einer anderen Methode als der der Trennung der Variablen zu lösen. Nach einer kurzen Recherche wurde festgestellt, dass die beste Methode die Transformationsmethode ist, genauer gesagt die Anwendung der Hankel-Transformation. Im nächsten Kapitel wird daher diese Integraltransformation vorgestellt und entwickelt, und es wird erläutert, warum sie für diese Art von Problem geeignet ist.
diff --git a/buch/papers/kreismembran/teil2.tex b/buch/papers/kreismembran/teil2.tex
index 7ed217f..6efda49 100644
--- a/buch/papers/kreismembran/teil2.tex
+++ b/buch/papers/kreismembran/teil2.tex
@@ -1,40 +1,116 @@
%
-% teil2.tex -- Beispiel-File für teil2
-%
% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
%
-\section{Teil 2
-\label{kreismembran:section:teil2}}
-\rhead{Teil 2}
-Sed ut perspiciatis unde omnis iste natus error sit voluptatem
-accusantium doloremque laudantium, totam rem aperiam, eaque ipsa
-quae ab illo inventore veritatis et quasi architecto beatae vitae
-dicta sunt explicabo. Nemo enim ipsam voluptatem quia voluptas sit
-aspernatur aut odit aut fugit, sed quia consequuntur magni dolores
-eos qui ratione voluptatem sequi nesciunt. Neque porro quisquam
-est, qui dolorem ipsum quia dolor sit amet, consectetur, adipisci
-velit, sed quia non numquam eius modi tempora incidunt ut labore
-et dolore magnam aliquam quaerat voluptatem. Ut enim ad minima
-veniam, quis nostrum exercitationem ullam corporis suscipit laboriosam,
-nisi ut aliquid ex ea commodi consequatur? Quis autem vel eum iure
-reprehenderit qui in ea voluptate velit esse quam nihil molestiae
-consequatur, vel illum qui dolorem eum fugiat quo voluptas nulla
-pariatur?
-
-\subsection{De finibus bonorum et malorum
-\label{kreismembran:subsection:bonorum}}
-At vero eos et accusamus et iusto odio dignissimos ducimus qui
-blanditiis praesentium voluptatum deleniti atque corrupti quos
-dolores et quas molestias excepturi sint occaecati cupiditate non
-provident, similique sunt in culpa qui officia deserunt mollitia
-animi, id est laborum et dolorum fuga. Et harum quidem rerum facilis
-est et expedita distinctio. Nam libero tempore, cum soluta nobis
-est eligendi optio cumque nihil impedit quo minus id quod maxime
-placeat facere possimus, omnis voluptas assumenda est, omnis dolor
-repellendus. Temporibus autem quibusdam et aut officiis debitis aut
-rerum necessitatibus saepe eveniet ut et voluptates repudiandae
-sint et molestiae non recusandae. Itaque earum rerum hic tenetur a
-sapiente delectus, ut aut reiciendis voluptatibus maiores alias
-consequatur aut perferendis doloribus asperiores repellat.
+\section{Die Hankel Transformation \label{kreismembran:section:teil2}}
+\rhead{Die Hankel Transformation}
+
+Hermann Hankel (1839--1873) war ein deutscher Mathematiker, der für seinen Beitrag zur mathematischen Analysis und insbesondere für die nach ihm benannte Transformation bekannt ist.
+Diese Transformation tritt bei der Untersuchung von Funktionen auf, die nur von der Entfernung des Ursprungs abhängen.
+Er studierte auch Funktionen, jetzt Hankel- oder Bessel- Funktionen genannt, der dritten Art.
+Die Hankel-Transformation, die die Bessel-Funktion enthält, taucht natürlich bei achsensymmetrischen Problemen auf, die in zylindrischen Polarkoordinaten formuliert sind.
+In diesem Abschnitt werden die Theorie der Transformation und einige Eigenschaften der Grundoperationen erläutert.
+
+\subsubsection{Hankel-Transformation \label{subsub:hankel_tansformation}}
+Wir führen die Definition der Hankel-Transformation \cite{lokenath_debnath_integral_2015} aus der zweidimensionalen Fourier-Transformation und ihrer Umkehrung ein, die durch:
+\begin{align}
+ \mathscr{F}\{f(x,y)\} & = F(k,l)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-i( \bm{\kappa}\cdot \mathbf{r})}f(x,y) \; dx dy,\label{equation:fourier_transform}\\
+ \mathscr{F}^{-1}\{F(x,y)\} & = f(x,y)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{i(\bm{\kappa}\cdot \mathbf{r}))}F(k,l) \; dx dy \label{equation:inv_fourier_transform}
+\end{align}
+wo $\mathbf{r}=(x,y)$ und $\bm{\kappa}=(k,l)$. Polarkoordinaten sind für diese Art von Problemen am besten geeignet, mit $(x,y)=r(\cos\theta,\sin\theta)$ und $(k,l)=\kappa(\cos\phi,\sin\phi)$ findet man $\bm{\kappa}\cdot\mathbf{r}=\kappa r(\cos(\theta-\phi))$ und danach:
+\begin{align}
+ F(k,\phi)=\frac{1}{2\pi}\int_{0}^{\infty}r \; dr \int_{0}^{2\pi}e^{-ikr\cos(\theta-\phi)}f(r,\theta) \; d\phi.
+ \label{equation:F_ohne_variable_wechsel}
+\end{align}
+Dann wird angenommen dass, $f(r,\theta)=e^{in\theta}f(r)$, was keine strenge Einschränkung ist, und es wird eine Änderung der Variabeln vorgenommen $\theta-\phi=\alpha-\frac{\pi}{2}$, um \eqref{equation:F_ohne_variable_wechsel} zu reduzieren:
+\begin{align}
+ F(k,\phi)=\frac{1}{2\pi}\int_{0}^{\infty}rf(r) \; dr \int_{\phi_{0}}^{2\pi+\phi_{0}}e^{in(\phi-\frac{\pi}{2})+i(n\alpha-kr\sin\alpha)} \; d\alpha,
+ \label{equation:F_ohne_bessel}
+\end{align}
+wo $\phi_{0}=(\frac{\pi}{2}-\phi)$.
+
+Unter Verwendung der Integraldarstellung der Besselfunktion vom Ordnung $n$ \eqref{buch:fourier:eqn:bessel-integraldarstellung}
+\begin{equation*}
+ J_n(\kappa r)=\frac{1}{2\pi}\int_{\phi_{0}}^{2\pi + \phi_{0}}e^{i(n\alpha-\kappa r \sin \alpha)} \; d\alpha
+ \label{equation:bessel_n_ordnung}
+\end{equation*}
+\eqref{equation:F_ohne_bessel} wird sie zu:
+\begin{align}
+ F(k,\phi)&=e^{in(\phi-\frac{\pi}{2})}\int_{0}^{\infty}rJ_n(\kappa r) f(r) \; dr \nonumber \\
+ &=e^{in(\phi-\frac{\pi}{2})}\tilde{f}_n(\kappa),
+ \label{equation:F_mit_bessel_step_2}
+\end{align}
+wo $\tilde{f}_n(\kappa)$ ist die \textit{Hankel-Transformation} von $f(r)$ und ist formell definiert durch:
+\begin{align}
+ \mathscr{H}_n\{f(r)\}=\tilde{f}_n(\kappa)=\int_{0}^{\infty}rJ_n(\kappa r) f(r) \; dr.
+ \label{equation:hankel}
+\end{align}
+
+\subsubsection{Inverse Hankel-Transformation \label{subsub:inverse_hankel_tansformation}}
+Ähnlich verhält es sich mit der inversen Fourier Transformation in Form von polaren Koordinaten unter der Annahme $f(r,\theta)=e^{in\theta}f(r)$ mit \eqref{equation:F_mit_bessel_step_2}, wird die inverse Fourier Transformation \eqref{equation:inv_fourier_transform}:
+
+\begin{align*}
+ e^{in\theta}f(r)&=\frac{1}{2\pi}\int_{0}^{\infty}\kappa \; d\kappa \int_{0}^{2\pi}e^{i\kappa r \cos (\theta - \phi)}F(\kappa,\phi) \; d\phi \\
+ &= \frac{1}{2\pi}\int_{0}^{\infty}\kappa \tilde{f}_n(\kappa) \; d\kappa \int_{0}^{2\pi}e^{in(\phi - \frac{\pi}{2})- i\kappa r \cos (\theta - \phi)} \; d\phi,
+\end{align*}
+was durch den Wechsel der Variablen $\theta-\phi=-(\alpha+\frac{\pi}{2})$ und $\theta_0=-(\theta+\frac{\pi}{2})$,
+
+\begin{align*}
+ &= \frac{1}{2\pi}\int_{0}^{\infty}\kappa \tilde{f}_n(\kappa) \; d\kappa \int_{\theta_0}^{2\pi+\theta_0}e^{in(\theta + \alpha - i\kappa r \sin\alpha)} \; d\alpha \\
+ &= e^{in\theta}\int_{0}^{\infty}\kappa J_n(\kappa r) \tilde{f}_n(\kappa) \; d\kappa,
+\end{align*}
+
+von \eqref{equation:bessel_n_ordnung} also ist, die inverse \textit{Hankel-Transformation} so definiert:
+\begin{align}
+ \mathscr{H}^{-1}_n\{\tilde{f}_n(\kappa)\}=f(r)=\int_{0}^{\infty}\kappa J_n(\kappa r) \tilde{f}_n(\kappa) \; d\kappa.
+ \label{equation:inv_hankel}
+\end{align}
+
+Anstelle von $\tilde{f}_n(\kappa)$, wird häufig für die Hankel-Transformation verwendet, indem die Ordnung angegeben wird.
+\eqref{equation:hankel} und \eqref{equation:inv_hankel} Integralen existieren für eine grosse Klasse von Funktionen, die normalerweise in physikalischen Anwendungen benötigt werden.
+Alternativ kann auch die berühmte Hankel-Transformationsformel verwendet werden,
+
+\begin{align*}
+ f(r) = \int_{0}^{\infty}\kappa J_n(\kappa r) \; d\kappa \int_{0}^{\infty} p J_n(\kappa p)f(p) \; dp,
+ \label{equation:hankel_integral_formula}
+\end{align*}
+um die Hankel-Transformation \eqref{equation:hankel} und ihre Inverse \eqref{equation:inv_hankel} zu definieren.
+Insbesondere die Hankel-Transformation der nullten Ordnung ($n=0$) und der ersten Ordnung ($n=1$) sind häufig nützlich, um Lösungen für Probleme mit der Laplace Gleichung in einer achsensymmetrischen zylindrischen Geometrie zu finden.
+
+\subsection{Operative Eigenschaften der Hankel-Transformation\label{sub:op_properties_hankel}}
+In diesem Kapitel werden die operativen Eigenschaften der Hankel-Transformation aufgeführt. Der Beweis für ihre Gültigkeit wird jedoch nicht analysiert.
+
+\begin{satz}{Skalierung:}
+ Wenn $\mathscr{H}_n\{f(r)\}=\tilde{f}_n(\kappa)$, dann:
+
+ \begin{equation*}
+ \mathscr{H}_n\{f(ar)\}=\frac{1}{a^{2}}\tilde{f}_n \left(\frac{\kappa}{a}\right), \quad a>0.
+ \end{equation*}
+\end{satz}
+
+\begin{satz}{Persevalsche Relation (Skalarprodukt bleibt erhalten):}
+Wenn $\tilde{f}(\kappa)=\mathscr{H}_n\{f(r)\}$ und $\tilde{g}(\kappa)=\mathscr{H}_n\{g(r)\}$, dann:
+
+\begin{equation*}
+ \int_{0}^{\infty}rf(r)g(r) \; dr = \int_{0}^{\infty}\kappa\tilde{f}(\kappa)\tilde{g}(\kappa) \; d\kappa.
+\end{equation*}
+\end{satz}
+
+\begin{satz}{Hankel-Transformationen von Ableitungen:}
+Wenn $\tilde{f}_n(\kappa)=\mathscr{H}_n\{f(r)\}$, dann:
+
+\begin{align*}
+ &\mathscr{H}_n\{f'(r)\}=\frac{\kappa}{2n}\left[(n-1)\tilde{f}_{n+1}(\kappa)-(n+1)\tilde{f}_{n-1}(\kappa)\right], \quad n\geq1, \\
+ &\mathscr{H}_1\{f'(r)\}=-\kappa \tilde{f}_0(\kappa),
+\end{align*}
+bereitgestellt dass $[rf(r)]$ verschwindet als $r\to0$ und $r\to\infty$.
+\end{satz}
+
+\begin{satz}
+Wenn $\mathscr{H}_n\{f(r)\}=\tilde{f}_n(\kappa)$, dann:
+\begin{equation*}
+ \mathscr{H}_n \left\{ \left( \nabla^2 - \frac{n^2}{r^2} f(r)\right)\right\}= \mathscr{H}_n\left\{\frac{1}{r}\frac{d}{dr}\left(r\frac{df}{dr}\right) - \frac{n^2}{r^2}f(r)\right\}=-\kappa^2\tilde{f}_{n}(\kappa),
+\end{equation*}
+bereitgestellt dass $rf'(r)$ und $rf(r)$ verschwinden für $r\to0$ und $r\to\infty$.
+\end{satz}
diff --git a/buch/papers/kreismembran/teil3.tex b/buch/papers/kreismembran/teil3.tex
index 73dee0f..7d5648a 100644
--- a/buch/papers/kreismembran/teil3.tex
+++ b/buch/papers/kreismembran/teil3.tex
@@ -3,38 +3,82 @@
%
% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
%
-\section{Teil 3
+\section{Lösungsmethode 2: Transformationsmethode
\label{kreismembran:section:teil3}}
-\rhead{Teil 3}
-Sed ut perspiciatis unde omnis iste natus error sit voluptatem
-accusantium doloremque laudantium, totam rem aperiam, eaque ipsa
-quae ab illo inventore veritatis et quasi architecto beatae vitae
-dicta sunt explicabo. Nemo enim ipsam voluptatem quia voluptas sit
-aspernatur aut odit aut fugit, sed quia consequuntur magni dolores
-eos qui ratione voluptatem sequi nesciunt. Neque porro quisquam
-est, qui dolorem ipsum quia dolor sit amet, consectetur, adipisci
-velit, sed quia non numquam eius modi tempora incidunt ut labore
-et dolore magnam aliquam quaerat voluptatem. Ut enim ad minima
-veniam, quis nostrum exercitationem ullam corporis suscipit laboriosam,
-nisi ut aliquid ex ea commodi consequatur? Quis autem vel eum iure
-reprehenderit qui in ea voluptate velit esse quam nihil molestiae
-consequatur, vel illum qui dolorem eum fugiat quo voluptas nulla
-pariatur?
-
-\subsection{De finibus bonorum et malorum
-\label{kreismembran:subsection:malorum}}
-At vero eos et accusamus et iusto odio dignissimos ducimus qui
-blanditiis praesentium voluptatum deleniti atque corrupti quos
-dolores et quas molestias excepturi sint occaecati cupiditate non
-provident, similique sunt in culpa qui officia deserunt mollitia
-animi, id est laborum et dolorum fuga. Et harum quidem rerum facilis
-est et expedita distinctio. Nam libero tempore, cum soluta nobis
-est eligendi optio cumque nihil impedit quo minus id quod maxime
-placeat facere possimus, omnis voluptas assumenda est, omnis dolor
-repellendus. Temporibus autem quibusdam et aut officiis debitis aut
-rerum necessitatibus saepe eveniet ut et voluptates repudiandae
-sint et molestiae non recusandae. Itaque earum rerum hic tenetur a
-sapiente delectus, ut aut reiciendis voluptatibus maiores alias
-consequatur aut perferendis doloribus asperiores repellat.
+\rhead{Lösungsmethode 2: Transformationsmethode}
+Die Hankel-Transformation wird dann zur Lösung der Differentialgleichung verwendet. Es müssen jedoch einige Änderungen an dem Problem vorgenommen werden, damit es mit den Annahmen übereinstimmt, die für die Verwendung der Hankel-Transformation erforderlich sind. Das heisst, dass die Funktion $u$ nur von der Entfernung zum Ausgangspunkt abhängt.
+
+\subsubsection{Transformation und Reduktion auf eine algebraische Gleichung\label{subsub:transf_reduktion}}
+Führt man also das Konzept einer unendlichen und achsensymmetrischen Membran ein:
+\begin{equation*}
+ \frac{\partial^2u}{\partial t^2}
+ =
+ c^2 \left(\frac{\partial^2 u}{\partial r^2}
+ +
+ \frac{1}{r}
+ \frac{\partial u}{\partial r} \right), \quad 0<r<\infty, \quad t>0
+ \label{eq:PDE_inf_membane}
+\end{equation*}
+
+\begin{align}
+ u(r,0)=f(r), \quad u_t(r,0) = g(r), \quad \text{für} \quad 0<r<\infty
+ \label{eq:PDE_inf_membane_RB}
+\end{align}
+
+Mit Anwendung der Hankel-Transformation nullter Ordnung in Abhängigkeit von $r$ auf die Gleichungen \eqref{eq:PDE_inf_membane} und \eqref{eq:PDE_inf_membane_RB}:
+
+\begin{align}
+ \tilde{u}(\kappa,t)=\int_{0}^{\infty}r J_0(\kappa r)u(r,t) \; dr,
+\end{align}
+bekommt man:
+
+\begin{equation*}
+ \frac{d^2 \tilde{u}}{dt^2} + c^2\kappa^2\tilde{u}=0,
+\end{equation*}
+
+\begin{equation*}
+ \tilde{u}(\kappa,0)=\tilde{f}(\kappa), \quad
+ \tilde{u}_t(\kappa,0)=\tilde{g}(\kappa).
+\end{equation*}
+Die allgemeine Lösung für diese Transformation lautet, wie in Gleighung \eqref{eq:cos_sin_überlagerung} gesehen, wie folgt
+
+\begin{equation*}
+ \tilde{u}(\kappa,t)=\tilde{f}(\kappa)\cos(c\kappa t) + \frac{1}{c\kappa}\tilde{g}(\kappa)\sin(c\kappa t).
+\end{equation*}
+Wendet man an nun die inverse Hankel-Transformation an, so erhält man die formale Lösung
+
+\begin{align}
+ u(r,t)=\int_{0}^{\infty}\kappa\tilde{f}(\kappa)\cos(c\kappa t) J_0(\kappa r) \; d\kappa +\frac{1}{c}\int_{0}^{\infty}\tilde{g}(\kappa)\sin(c\kappa t)J_0(\kappa r) \; d\kappa.
+ \label{eq:formale_lösung}
+\end{align}
+
+\subsubsection{Erfüllung der Anfangsbedingungen\label{subsub:erfüllung_AB}}
+Es wird in Folgenden davon ausgegangen, dass sich die Membran verformt und zum Zeitpunkt $t=0$ freigegeben wird
+
+\begin{equation*}
+ u(r,0)=f(r)=Aa(r^2 + a^2)^{-\frac{1}{2}}, \quad u_t(r,0)=g(r)=0
+\end{equation*}
+so dass $\tilde{g}(\kappa)\equiv 0$ und
+\begin{equation*}
+ \tilde{f}(\kappa)=Aa\int_{0}^{\infty}r(a^2 + r^2)^{-\frac{1}{2}} J_0 (\kappa r) \; dr=\frac{Aa}{\kappa}e^{-a\kappa}
+\end{equation*}
+Die formale Lösung \eqref{eq:formale_lösung} lautet also
+\begin{align*}
+ u(r,t)&=Aa\int_{0}^{\infty}e^{-a\kappa} J_0(\kappa r)\cos(c\kappa t) \; dk=AaRe\int_{0}^{\infty}e^{-\kappa(a+ict)} J_0(\kappa r) \; dk\\
+ &=AaRe\left\{r^2+\left(a+ict\right)^2\right\}^{-\frac{1}{2}}
+\end{align*}
+
+Nimmt man jedoch die allgemeine Lösung mit Summationen,
+
+\begin{align}
+ u(r, t) = \displaystyle\sum_{m=1}^{\infty} J_0 (k_{m}r)[a_{m}\cos(c \kappa_{m} t)+b_{m}\sin(c \kappa_{m} t)]
+ \label{eq:lösung_unendliche_generelle}
+\end{align}
+kann man die Lösungsmethoden 1 und 2 vergleichen.
+
+\subsection{Vergleich der Analytischen Lösungen
+\label{kreismembran:vergleich}}
+Bei der Analyse der Gleichungen \eqref{eq:lösung_endliche_generelle} und \eqref{eq:lösung_unendliche_generelle} fällt sofort auf, dass die Gleichung \eqref{eq:lösung_unendliche_generelle} nicht mehr von $m$ und $n$ abhängt, sondern nur noch von $n$ \cite{nishanth_p_vibrations_2018}. Das macht Sinn, denn $n$ beschreibt die Anzahl der Knotenlinien, und in einer unendlichen Membran gibt es keine. Tatsächlich werden $a_{m0}$, $b_{m0}$ und $\kappa_{m0}$ in $a_m$, $b_m$ bzw. $\kappa_m$ umbenannt. Die beiden Termen $\cos(n\varphi)$ und $\sin(n\varphi)$ verschwinden ebenfalls, da für $n=0$ der $\cos(n\varphi)$ gleich 1 und der $\sin(n \varphi)$ gleich 0 ist.
+Die Funktion hängt also nicht mehr von der Besselfunktionen $n$-ter Ordnung ab, sondern nur von der $0$-ter Ordnung.
diff --git a/buch/papers/kreismembran/teil4.tex b/buch/papers/kreismembran/teil4.tex
new file mode 100644
index 0000000..c124354
--- /dev/null
+++ b/buch/papers/kreismembran/teil4.tex
@@ -0,0 +1,16 @@
+%
+% einleitung.tex -- Beispiel-File für die Einleitung
+%
+% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
+%
+\section{Lösungsmethode 3: Simulation
+ \label{kreismembran:section:teil4}}
+\paragraph{TODO Einleitung}
+
+Um numerisch das Verhalten einer Membran zu ermitteln, muss eine numerische Darstellung definiert werden.
+Die Membran wird hier in Form der Matrix $ A $ digitalisiert.
+Jedes Element $ A_{ij} $ steh für die Auslenkung der Membran $ u(x,y,t) $ an der Stelle $ \{x,y\}=\{i,j\} $.
+Die zeitliche Dimension wird in Form des Array $ X[] $ aus $ v \times A $ Matrizen dargestellt.
+Das Element auf Zeile $ i $, Spalte $ j $ der $ w $-ten Matrix von $ X[] $ also $ X[w]_{ij} $ entspricht der Auslenkung $ u(i,j,w) $.
+
+\paragraph{title} \ No newline at end of file