Theorie - Statische Analyse: Erste vollständige Version eingebaut

Preamble/Main.tex

@@ -367,7 +367,7 @@
 \newcommand*{\ddt}[1]{%
   \overset{\mbox{\large\bfseries .\hspace{-0.1ex}.}}{#1}}
 
-%\usepackage{accents}
+\usepackage{accents} % for tensor 
 \usepackage{ushort}  % for tensor
 %\usepackage{ulem}    % not for math mode, too much space below with uuline; for tensor
 

sections/References.bib

@@ -20,6 +20,17 @@
   isbn = {3-834-80025-2},
   gender={sm},
 }
+@BOOK{becker02,
+  author = {Becker, Wilfried and Gross, Dietmar},
+  title= {Mechanik elastischer Körper und Strukturen},
+  subtitle= {},
+  publisher = {Springer-Verlag},
+  year = {2002},
+  address = {Berlin Heidelberg},
+  edition = {},
+  isbn = {3-540-43511-5},
+  gender={sm},
+}
 @BOOK{tm409,
   author = {Gross, Dietmar and Hauger, Werner and Wriggers, Peter},
   title = {Technische Mechanik},

sections/SymbolsCodes.tex

@@ -9,6 +9,7 @@
 \(L\) & mm & Länge \\
 [0.25cm]
 \(\nu\) &  & Querkontraktionszahl \\
+\(\tensorII{\sigma}\) & MPa & Cauchy Spannungstensor \\
 \end{longtable*}
 
 

sections/Theorie.tex

@@ -24,17 +24,20 @@ Die \ac{FEM} umfasst eine Vielzahl von Methodiken physikalische Fragestellungen
 So sind für strukturmechanische Probleme einer \ac{WEA} beispielsweise die statische Durchbiegung der Rotorblätter aufgrund Eigengewicht und Windlasten 
 als auch Eigenformen und die Anlagenbelastung bei drehenden Rotorblätter von Interesse. 
 In den nachfolgenden Abschnitten wird auf die Grundlagen der statischen und dynamischen Analyse eingegangen.
-Zur weiteren Vertiefung der Themengebiete beziehungsweise bei Interesse zur Lösung von anderen Problemstellungen sei unter anderem auf die Literatur \citep{bathe86} und \cite{klein05} sowie das \acs{ANSYS}"=Programmhandbuch %\cite{ansys}
+Zur weiteren Vertiefung der Themengebiete beziehungsweise bei Interesse zur Lösung von anderen Problemstellungen sei unter anderem auf die Literatur \citep{bathe86}, \cite{klein05} und \cite{becker02} sowie das \acs{ANSYS}"=Programmhandbuch %\cite{ansys}
 verwiesen.
 
 Im Allgemeinen werden mit der \ac{FEM} Differentialgleichungen beziehungsweise Systeme von Differentialgleichungen gelöst.
 Am Beispiel eines Biegebalkens ist es die Biegedifferentialgleichung
 und bei dynamischen Problemen eine Bewegungsdifferentialgleichung.
 Gewonnen werden diese Differentialgleichung durch Gleichgewichtsbetrachtungen eines beliebigen Festkörpers, der entweder statisch bestimmt oder statisch überbestimmt gelagert und mit äußeren Lasten beaufschlagt ist.
-Die \ac{FEM} umfasst zur Lösung der Feldgröße dabei drei grundlegende Schritte mit der die Differentialgleichung als FE"=Gleichungssystem beschrieben wird; die Partitionierung, die Approximation und die Assemblierung.
+Die \ac{FEM} umfasst zur Lösung der Feldgröße dabei drei grundlegende Schritte mit der die Differentialgleichung als FE"=Gleichungssystem beschrieben wird;
+die Partitionierung wobei das Berechnungsgebiet in Elementen mit Knoten eingeteilt wird und das Gebiet als Summe der Elemente anzusehen ist,
+die Approximation bei der die gesuchte Verformung durch Polynome angenähert wird und die Elemente fortan in Abhängigkeit der Knotenverformung beschrieben sind und schließlich
+die Assemblierung welche das globale \glqq diskrete\grqq\ Gleichgewicht (FE"=Gleichungssystem) aufstellt mit der das Gebiet als Summe der beschriebenen Elemente anzusehen ist, wobei kinematische Bedingungen gleicher Verformungen zu berücksichtigen sind.
 
 \subsubsection{Statische Analysen}
-Mit dem Kräftegleichgewicht am infinitesimalen Element ergibt sich für ein beliebigen Körper~\(V\) im allgemeinen dreidimensionalen Fall die klassische elliptische Differentialgleichung mit den Randbedingungen auf der Körperoberfläche~\(\partial V\) -- in differentielle Formulierung -- 
+Mit dem Kräftegleichgewicht am infinitesimalen Element ergibt sich für ein beliebigen Körper~\(V\) im allgemeinen dreidimensionalen Fall die klassische elliptische Differentialgleichung mit den Randbedingungen auf der Körperoberfläche~\(\partial V\) -- in differentielle Tensor"=Formulierung -- 
 \[
   \begin{Array}{rcl}
     \nabla\cdot\tensorII{\sigma} + \tensorI{f} &= \tensorI{0} & \quad\text{in } V, \\
@@ -43,17 +46,124 @@ Mit dem Kräftegleichgewicht am infinitesimalen Element ergibt sich für ein bel
     &&\quad\text{wobei }\partial V_1 \cup \partial V_2 = \partial V.
   \end{Array}
 \]
-Die partiell abgeleiteten inneren Spannungen~\(\tensorII{\sigma}\) stehen mit den  Volumenkräften~\(\tensorI{f}\) im Gleichgewicht. 
-Dabei beschreiben die Randbedingungen zum einen Verschiebungen~\(\tensorI{u}_0\) auf einem Teil der Körperoberfläche~\(\partial V_1\) und zum anderen Belastungen in Form eines Spannungsvektors~\(\tensorI{t}\) auf den restlichen Teil der Körperoberfläche~\(\partial V_2\).
+Die Summe der partiell abgeleiteten inneren Spannungen~\(\sigma_{ji,j}(\tensorI{x})\) stehen mit der ortsabhängigen Volumenkraft~\(f_i(\tensorI{x})\) der jeweiligen Richtung im Gleichgewicht. 
+Dabei beschreiben die Randbedingungen zum einen Verschiebungen~\(\tensorI{u}_0(\tensorI{x})\) auf einem Teil der Körperoberfläche~\(\partial V_1\) und zum anderen Belastungen in Form eines Spannungsvektors~\(\tensorI{t}(\tensorI{x})\) auf den restlichen Teil der Körperoberfläche~\(\partial V_2\).
 Zur näherungsweisen Verarbeitung der Differentialgleichung wird hier das Variationsprinzip -- das sogenannte \emph{Prinzip der virtuellen Verrückung} -- herangezogen, mit der ein Ersatzgleichgewichtsgleichung formuliert wird. 
 
+\paragraph{Globales Gleichgewicht über schwache Formulierung des Randwertproblems}~\\
 Auf Grundlage der Differentialgleichung (starke Formulierung) erfolgt die schwache Formulierung durch Multiplikation einer Testfunktion beziehungsweise virtuellen Verrückung \(\delta\tensorI{u}\), welche die kinematischen beziehungsweise wesentlichen Randbedingungen erfüllen muss (\(\forall \delta\tensorI{u}\in C^1(V)\cap C(\overline{V}),~ \delta\tensorI{u} = \tensorI{u}_0 \text{ auf } \partial V_1\)), mit anschließender Integration über das Berechnungsgebiet \(V\).
 \[
-  \underbrace{\int\limits_{V\vphantom{V_1}}\!\!\!(\tensorII{\sigma} : \nabla\delta\tensorI{u})\dif V }_{\delta W\ti{i}} = \underbrace{\!\!\int\limits_V\!\!\!(\tensorI{f}\cdot\delta\tensorI{u})\dif V + \!\!\int\limits_{\partial V_2}\!\!\!(\tensorI{t}\cdot\delta\tensorI{u})\dif A}_{\delta W\ti{a}}
+  \underbrace{\int\limits_{V\vphantom{V_1}}\!\!\!(\tensorII{\sigma} : \nabla\delta\tensorI{u})\dif V }_{\delta W\ti{i}} = \underbrace{\!\!\int\limits_V\!\!\!(\tensorI{f}\cdot\delta\tensorI{u})\dif V + \!\!\int\limits_{\partial V_2}\!\!\!(\tensorI{t}\cdot\delta\tensorI{u})\dif A 
+%  + \sum\tensorI{F}_i\cdot\delta\tensor{u}_i
+  }_{\delta W\ti{a}}
 \]
-wobei \( \nabla\cdot(\tensorII{\sigma}\cdot\delta\tensorI{u}) = \nabla\cdot\tensorII{\sigma}\cdot\delta\tensorI{u} + \tensorII{\sigma}:\nabla\delta\tensorI{u} \)
-und \( \int_V \nabla\cdot\tensorII{\sigma}\dif V = \int_{\partial V}\tensorII{\sigma}\cdot\tensorI{n}\dif A \)
-(Gauß'sche Integralsatz)
+wobei 
+%\( \nabla\cdot(\tensorII{\sigma}\cdot\delta\tensorI{u}) = \nabla\cdot\tensorII{\sigma}\cdot\delta\tensorI{u} + \tensorII{\sigma}:\nabla\delta\tensorI{u} \)
+%(Produktregel mit symmetrischen Tensor \(\tensorII{\sigma}^\T=\tensorII{\sigma}\))
+%und \( \int_V \nabla\cdot\tensorII{\sigma}\dif V = \int_{\partial V}\tensorII{\sigma}\cdot\tensorI{n}\dif A \)
+%(Gauß'sche Integralsatz)
+\( \int_V \nabla\cdot\tensorII{\sigma}\cdot\delta\tensorI{u}\dif V = \int_A(\tensorII{\sigma}\cdot\tensorI{n})\cdot\delta\tensorI{u}\dif A - \int_V\tensorII{\sigma}:(\nabla\delta\tensorI{u})\dif V \)
+(Green'sche Integralsatz)
+sowie
+\( \nabla\delta\tensorI{u} = \tensorII{\varepsilon}(\delta\tensorI{u}) = \delta\tensorII{\varepsilon} \) im linearen Fall. Die Integration erfolgt im linearen Fall über das unverformte und bei dem nichtlinearen Fall über das verformte Volumen.
+
+\paragraph{Materialgesetz}~\\
+Im linear elastischen Fall besteht zwischen der Spannung und der Dehnung das folgende Materialgesetz
+\[
+  \begin{Array}{rrcll}
+  &\tensorII{\sigma} &=& \tensorIV{C}:\tensorII{\varepsilon}  \\
+  \text{bzw. in Ingeniersnotation}~~~~~ & \tensor{\sigma}_{6\times1} &=& \tensor{C}_{6\times6} \tensor{\varepsilon}_{6\times1} & \quad\text{da } \tensorII{\sigma}^\T = \tensorII{\sigma} ~\wedge~ \tensorII{\varepsilon}^\T = \tensorII{\varepsilon}
+  \end{Array}
+\]
+Im allgemeinem Fall entsprechen die enthaltene Richtungen der Dehnung nicht die enthaltene Richtung der Kraftkomponente. Es kann somit jede Spannungskomponente von jeder Dehnungskomponente abhängen, weshalb entgegen dem vergleichsweise eindimensionalen Fall \(\sigma(x) = E\varepsilon(x)\) keine direkte Proportionalität vorherrscht.
+Damit ist die Elastizitätsbeziehung kein Skalar sondern ein Tensor vierter Stufe.
+
+Für den physikalisch nichtlinearen Fall sind die Spannungen allgemein von dem Verlauf der Verformung abhängig \(\tensorII{\sigma} = \tensorII{\sigma}(\tensorI{u})\).
+Für eine genauere Darstellung wird auf \cite{becker02} verwiesen.
+
+\paragraph{Kinematik}~\\
+Die kinematische Beziehung der örtlich dreidimensionalen Dehnung, für kleine Verschiebungen, 
+%sind die partiellen Verhältnisse der Längenänderung zur ursprünglichen Länge
+ist als Verzerrungstensor gegeben
+\[
+  \begin{Array}{rrll} 
+  &\varepsilon_{ij} &= \displaystyle\frac{1}{2}\left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right) = \frac{1}{2} \left( u_{i,j} + u_{j,i} \right) \\
+  \text{bzw. in Ingeniersnotation}~~~~~ & \tensor{\varepsilon}_{6\times1} &= \tensor{D}_{6\times3}\tensor{u}_{3\times1} \quad \quad\text{da } \tensorII{\varepsilon}^\T = \tensorII{\varepsilon}  \text{ mit } \tensor{D}: \text{Operatormatrix}
+  \end{Array}
+\]
+% \[ \tensor{D} = \begin{bmatrix}
+% \frac{\partial}{\partial x} & 0 & 0 \\
+% 0 & \frac{\partial}{\partial y} & 0 \\
+% 0 & 0 & \frac{\partial}{\partial z} \\
+% \frac{\partial}{\partial y} & \frac{\partial}{\partial x} & 0 \\
+% 0 & \frac{\partial}{\partial z} & \frac{\partial}{\partial y} \\
+% \frac{\partial}{\partial z} & 0 & \frac{\partial}{\partial x}
+% \end{bmatrix} \]
+Für den geometrisch nichtlinearen Fall sind die Dehnungen allgemein von dem Verlauf der Verformung abhängig \(\tensorII{\varepsilon} = \tensorII{\varepsilon}(\tensorI{u})\).
+Für eine genauere Darstellung wird auf \cite{becker02} verwiesen.
+
+\paragraph{Das FE"=Gleichungssystem}~\\
+Das Aufstellen des FE"=Gleichungssystems startet mit der \emph{Vernetzung} beziehungsweise \emph{Partitionierung}
+\[
+  V \approx \bigcup V^{(e)}
+\] 
+Hiermit ist das Gesamtvolumen die Vereinigung der diskreten Einzelvolumen beziehungsweise Elemente.
+
+Die \emph{Approximation} als Überlagerung aller Formfunktionen mit den jeweiligen Knotenverformungen lautet
+\[
+  \tensor{u}\big|_{V^{(e)}} \approx \tensor{u}\ti{fe}\big|_{V^{(e)}} = \tensor{N}\big|_{V^{(e)}} \tensor{\hat{u}} \quad \text{mit } \forall \tensor{u}\ti{fe} \in C^0(V)
+\]
+worin \( u\ti{fe} \) die FE"=Verformung, \(N\) die Formfunktionen und \(\hat{u}\) die Knotenverformung sind. 
+Werden für die Volumenvernetzung acht Knoten Hexaeder-Elemente mit trilinearen Formfunktionen verwendet nimmt die Matrixnotation folgende Darstellung an
+\[
+  \tensor{u}_{3\times1} = \tensor{N}_{\!3\times24}\, \tensor{\hat{u}}_{24\times1}   
+\]
+Mit der \emph{Assemblierung} wird das Gebiet als Summe aller Elementintegrale abgebildet.
+\[
+  \underbrace{\int\limits_{V\vphantom{V_1}}\!\!\!\tensorII{\sigma} : \nabla\delta\tensorI{u} \dif V }_{\delta W\ti{i}} = 
+  \underbrace{\sum\limits_{(e)} \int\limits_{~V^{(e)}}\!\!\!\tensorII{\sigma} : \nabla\delta\tensorI{u} \dif V }_{\delta W\ti{i}\big|_{V^{(e)}}}
+  \quad , \quad 
+  \underbrace{\int\limits_{V\vphantom{V_1}}\!\!\!\tensorI{f}\cdot\delta\tensorI{u}\dif V}_{\delta W\ti{a}} =
+  \underbrace{\sum\limits_{(e)} \int\limits_{~V^{(e)}}\!\!\!\tensorI{f}\cdot\delta\tensorI{u}\dif V}_{\delta W\ti{a}\big|_{V^{(e)}}}
+\]
+Für kleine Verformungen ergibt sich nun, mit den Materialgesetz, der kinematischen Beziehung und der Assemblierung, die virtuelle innere Arbeit zu
+\[
+  \begin{Array}{rll} \displaystyle
+  \sum\limits_{(e)} \int\limits_{~V^{(e)}}\!\!\!\tensorII{\sigma} : \nabla\delta\tensorI{u} \dif V & \displaystyle = 
+  \sum\limits_{(e)} \int\limits_{~V^{(e)}}\!\!\!\tensor{C} \tensor{\varepsilon} \cdot \delta\tensor{\varepsilon} \dif V  = 
+  \sum\limits_{(e)} \delta\tensor{\hat{u}}^\T \!\!\!\!\!\int\limits_{~V^{(e)}}\!\!\! \underbrace{(\tensor{D}\tensor{N})^\T}_{\tensor{B}^\T} \tensor{C} \underbrace{(\tensor{D}\tensor{N})}_{\tensor{B}\vphantom{B^\T}} \dif V \tensor{\hat{u}}  \\
+    & \displaystyle = 
+  \sum\limits_{(e)} \delta\tensor{\hat{u}}^\T \tensor{K}^{(e)} \tensor{\hat{u}} \quad\text{mit } \tensor{K}^{(e)} =  \!\!\int\limits_{~V^{(e)}}\!\!\! \tensor{B}^\T \tensor{C} \tensor{B} \dif V \\
+    & \displaystyle = 
+  \delta\tensor{\hat{u}}^\T \tensor{K} \tensor{\hat{u}} \quad \text{mit } \tensor{K} = \sum\limits_{(e)} \tensor{K}^{(e)}
+  \end{Array}
+\]
+Hierin ist \(\tensor{K}^{(e)}\) die Elementsteifigkeitmatrix und \(\tensor{K}\) die Gesamtsteifigkeitsmatrix. Für die äußere virtuelle Arbeit zu
+\[
+  \begin{Array}{rll} \displaystyle
+  \sum\limits_{(e)} \int\limits_{~V^{(e)}}\!\!\!\tensorI{f}\cdot\delta\tensorI{u}\dif V &\displaystyle =
+  \sum\limits_{(e)} \delta\tensor{\hat{u}}^\T \!\!\!\!\!\int\limits_{~V^{(e)}}\!\!\!\tensor{N}^\T\tensor{f}\dif V \\
+  & \displaystyle =
+  \sum\limits_{(e)} \delta\tensor{\hat{u}}^\T\tensor{\hat{r}}^{(e)} \quad \text{mit } \tensor{\hat{r}}^{(e)} = \!\!\int\limits_{~V^{(e)}}\!\!\!\tensor{N}^\T\tensor{f}\dif V\\
+  & \displaystyle =
+  \delta\tensor{\hat{u}}^\T\tensor{\hat{r}} \quad \text{mit } \tensor{\hat{r}} = \sum\limits_{(e)} \tensor{\hat{r}}^{(e)}
+  \end{Array}
+\]
+Hierin ist \(\tensor{\hat{r}}\) der Knotenlastvektor.
+Mit den beiden virtuellen Arbeitsaussagen lautet das globale Gleichgewicht
+\[
+  \delta\tensor{\hat{u}}^\T \tensor{K} \tensor{\hat{u}} = \delta\tensor{\hat{u}}^\T\tensor{\hat{r}}
+\]
+Abschließend wird die virtuelle Verschiebung ausgeklammert
+\[
+  \delta\tensor{\hat{u}}^\T \tensor{K} \tensor{\hat{u}} = \delta\tensor{\hat{u}}^\T\tensor{\hat{r}} \quad\rightarrow\quad
+  \delta\tensor{\hat{u}}^\T (\tensor{K} \tensor{\hat{u}} - \tensor{\hat{r}}) = 0
+\]
+Aus dem Variationsargument, dass diese Gleichung für beliebige virtuelle Verschiebungen gelten soll, muss der Klammerausdruck zu Null werden.
+\[
+  \tensor{K} \tensor{\hat{u}} = \tensor{\hat{r}}
+\]
+
 
 
 \subsubsection{Dynamische Analysen}