3. Modélisation A#
3.1. Caractéristiques de la modélisation#
La modélisation 3D est mise en œuvre. On teste les matériaux isotrope transverse et orthotrope (avec éventuellement prise en compte de déformations d’origines thermiques)
Remarques:
L’isotropie transverse n’est pas testée pour les contraintes planes car ce cas correspond à l’isotropie.
Pour le cas axisymétrique le champ de contraintes dépend du point de calcul.
Ce point est choisi au point d’intégration du triangle (i.e. c’est le centre de gravité du triangle).
On rappelle que l’orthotropie dans un repère quelconque n’est pas disponible pour la modélisation en Fourier car il y a alors couplage de toutes les composantes du tenseur de contraintes:
La mise en œuvre actuelle permet de n’utiliser que les composantes symétriques à partir desquelles on peut retrouver les composantes antisymétriques mais pour que ce soit possible, il ne faut pas que les glissements induisent des contraintes de traction.
3.2. Caractéristiques du maillage#
On a un élément tétraèdre à 4 nœuds \(\mathrm{ABCD}\) .
3.3. Valeurs testées#
Identification |
Référence |
Cas de l’isotropie transverse 3D |
|
nom du résultat: \(\mathit{Mest1}\) |
|
champ de déplacement |
|
\(\mathit{dy}(c)\) |
21 |
champ EPSI_ELGA |
|
\(\mathit{EPXY}\) |
3 |
\(\mathit{EPXZ}\) |
4 |
\(\mathit{EPYZ}\) |
6 |
champ SIEF_ELGA |
|
\(\mathit{SIXX}\) |
43310.760 |
\(\mathit{SIYY}\) |
72798.710 |
\(\mathit{SIZZ}\) |
62459.356 |
\(\mathit{SIXY}\) |
39567.891 |
\(\mathit{SIXZ}\) |
31078.597 |
\(\mathit{SIYZ}\) |
84049.301 |
champ SIGM_ELNO |
|
\(\mathit{SIXX}\) |
43310.760 |
champ emel–elga Ep |
1.19123 E6 |
Champ emel–elno–elga Ep |
1.19123 E6 |
Cas de l’orthotropie 3D |
|
nom du résultat: \(\mathit{Mest2}\) |
|
champ de déplacement |
|
\(\mathit{dy}(c)\) |
21 |
champ EPSI_ELGA |
|
\(\mathit{EPXY}\) |
3 |
\(\mathit{EPXZ}\) |
4 |
\(\mathit{EPYZ}\) |
6 |
champ SIEF_ELGA |
|
\(\mathit{SIXX}\) |
601.8754 |
\(\mathit{SIYY}\) |
80053.665 |
\(\mathit{SIZZ}\) |
78596.607 |
\(\mathit{SIXY}\) |
83948.263 |
\(\mathit{SIXZ}\) |
17339.093 |
\(\mathit{SIYZ}\) |
126571.71 |
champ enel–elga \(\mathit{Ep}\) |
1.55286.106 |
champ enel–elno–elga \(\mathit{Ep}\) |
1.55286.106 |
Cas de l’orthotropie avec prise en compte des déformations thermiques (commande STAT_NON_LINE ) |
|
nom du résultat: \(\mathit{Mest3}\) |
|
champ de déplacement |
|
\(\mathit{dy}(c)\) |
21 |
Identification |
Référence |
champ EPSI_ELGA |
|
\(\mathit{EPXY}\) |
3 |
\(\mathit{EPXZ}\) |
4 |
\(\mathit{EPYZ}\) |
6 |
champ SIEF_ELGA |
|
\(\mathit{SIXX}\) |
1226.2014 |
\(\mathit{SIYY}\) |
78597.064 |
\(\mathit{SIZZ}\) |
76585.792 |
\(\mathit{SIXY}\) |
83710.907 |
\(\mathit{SIXZ}\) |
17255.703 |
\(\mathit{SIYZ}\) |
126657.367 |
Cas de l’orthotropie avec prise en compte des déformations thermiques (commande MECA_STATIQUE ) |
|
nom du résultat: \(\mathit{Mest4}\) |
|
champ de déplacement |
|
\(\mathit{dy}(c)\) |
21 |
champ EPSI_ELGA |
|
\(\mathit{EPXY}\) |
3 |
\(\mathit{EPXZ}\) |
4 |
\(\mathit{EPYZ}\) |
6 |
champ SIEF_ELGA |
|
\(\mathit{SIXX}\) |
1226.2014 |
\(\mathit{SIYY}\) |
78597.064 |
\(\mathit{SIZZ}\) |
76585.792 |
\(\mathit{SIXY}\) |
83710.907 |
\(\mathit{SIXZ}\) |
17255.703 |
\(\mathit{SIYZ}\) |
126657.367 |