von Dirk » So 13. Okt 2013, 14:00
Hallo,
ich möchte die Anisotropie eines Einkristalls mittels des E-Moduls im dreidimensionalen Raum in Abhängigkeit der Kristallrichtungen darstellen.
Dazu soll zunächst die Elastizitätsmatrix von Silizium genutzt werden:
$\boldsymbol{C} =
\begin{bmatrix}
C_{11} & C_{12} & C_{12} &0&0&0 \\
C_{12} & C_{11} & C_{12} &0&0&0 \\
C_{12} & C_{12} & C_{11} &0&0&0 \\
0&0&0&C_{44}&0&0\\
0&0&0&0&C_{44}&0\\
0&0&0&0&0&C_{44}\\
\end{bmatrix}$
Elastischen Konstanten:
$C_{11} = 165,7\,\mathrm{GPa}, \; C_{12} = 63,9\,\mathrm{GPa}, \; C_{44} = 79,6\,\mathrm{GPa}$
Diese Konstanten gelten für ein parallel zu den <100> -Richtungen ausgerichtetes Koordinatensystem. Geht man von einem einachsigen Spannungszustand in einer bestimmten Richtung aus, so kann man den E-Modul für diese Richtung berechnen. Führt man diese Rechnung für beliebige Richtungen durch und trägt die erhaltenen Werte als Abstand um einen Mittelpunkt in ihrer jeweiligen Richtung ab, so entsteht der angehangene Körper.
Zur Umsetzung bin ich auf das Paket 3dplot gestoßen. (
http://www.ctan.org/tex-archive/graphic ... kz-3dplot/) Dazu hangele ich mich an folgendem Minimalbeispiel entlang.
\documentclass{minimal}
\usepackage{verbatim}
\usepackage{tikz}
\usepackage{3dplot}
\begin{document}
\tdplotsetmaincoords{70}{135}
\begin{tikzpicture}[scale=1,line join=bevel,tdplot_main_coords, fill opacity=.7]
\pgfsetlinewidth{.1pt}
\tdplotsphericalsurfaceplot[parametricfill]{72}{36}%
{sqrt(15/2)*sin(\tdplottheta)*cos(\tdplottheta)}{black}{\tdplotphi}%
{\draw[color=black,-latex] (0,0,0) -- (2,0,0) node[anchor=north east]{$x$};}%
{\draw[color=red,-latex] (0,0,0) -- (0,2,0) node[anchor=north west]{$y$};}%
{\draw[color=black,-latex] (0,0,0) -- (0,0,2) node[anchor=south]{$z$};}%
\end{tikzpicture}
\end{document}
Nun ist jedoch die Frage, wie ich meine Elastizitätsmatrix hier implementiere und diese dann auch zum Rotieren um den Koordinatenursprung bringe.
Kann hier jemand einen Tip geben bzw. gibt es alternative Lösungsmöglichkeiten?
Viele Grüße
Dirk
Hallo,
ich möchte die Anisotropie eines Einkristalls mittels des E-Moduls im dreidimensionalen Raum in Abhängigkeit der Kristallrichtungen darstellen.
Dazu soll zunächst die Elastizitätsmatrix von Silizium genutzt werden:
[code]
$\boldsymbol{C} =
\begin{bmatrix}
C_{11} & C_{12} & C_{12} &0&0&0 \\
C_{12} & C_{11} & C_{12} &0&0&0 \\
C_{12} & C_{12} & C_{11} &0&0&0 \\
0&0&0&C_{44}&0&0\\
0&0&0&0&C_{44}&0\\
0&0&0&0&0&C_{44}\\
\end{bmatrix}$
[/code]
Elastischen Konstanten:
[code]
$C_{11} = 165,7\,\mathrm{GPa}, \; C_{12} = 63,9\,\mathrm{GPa}, \; C_{44} = 79,6\,\mathrm{GPa}$
[/code]
Diese Konstanten gelten für ein parallel zu den <100> -Richtungen ausgerichtetes Koordinatensystem. Geht man von einem einachsigen Spannungszustand in einer bestimmten Richtung aus, so kann man den E-Modul für diese Richtung berechnen. Führt man diese Rechnung für beliebige Richtungen durch und trägt die erhaltenen Werte als Abstand um einen Mittelpunkt in ihrer jeweiligen Richtung ab, so entsteht der angehangene Körper.
[img]http://ft.carstenkloehn.de/Joomla/images/Projekte/2005/Knorr_DA05/Si_EMod_Koerper_tp.png[/img]
Zur Umsetzung bin ich auf das Paket 3dplot gestoßen. ([url]http://www.ctan.org/tex-archive/graphics/pgf/contrib/tikz-3dplot/[/url]) Dazu hangele ich mich an folgendem Minimalbeispiel entlang.
[code]
\documentclass{minimal}
\usepackage{verbatim}
\usepackage{tikz}
\usepackage{3dplot}
\begin{document}
\tdplotsetmaincoords{70}{135}
\begin{tikzpicture}[scale=1,line join=bevel,tdplot_main_coords, fill opacity=.7]
\pgfsetlinewidth{.1pt}
\tdplotsphericalsurfaceplot[parametricfill]{72}{36}%
{sqrt(15/2)*sin(\tdplottheta)*cos(\tdplottheta)}{black}{\tdplotphi}%
{\draw[color=black,-latex] (0,0,0) -- (2,0,0) node[anchor=north east]{$x$};}%
{\draw[color=red,-latex] (0,0,0) -- (0,2,0) node[anchor=north west]{$y$};}%
{\draw[color=black,-latex] (0,0,0) -- (0,0,2) node[anchor=south]{$z$};}%
\end{tikzpicture}
\end{document}
[/code]
Nun ist jedoch die Frage, wie ich meine Elastizitätsmatrix hier implementiere und diese dann auch zum Rotieren um den Koordinatenursprung bringe.
Kann hier jemand einen Tip geben bzw. gibt es alternative Lösungsmöglichkeiten?
Viele Grüße
Dirk