v6.04.179 SSNV179 - Cube sous fluage via la loi LEMA_SEUIL#
Résumé:
Ce test a pour but de valider la loi LEMA_SEUIL dérivée de la loi de LEMAITRE classique. En particulier, nous nous attarderons sur l’activation du seuil propre à cette loi. On réalise donc un essai de fluage sur une géométrie simple à savoir ici un cube.
Une modélisation 3D avec des éléments HEXA8 est disponible actuellement.
Solution de référence#
Méthode de calcul utilisée pour la solution de référence#
Le but de la solution de référence est de calculer analytiquement la valeur du seuil à partir duquel le fluage apparaît.
Quelques résultats de non régression sur les déplacements au dernier pas de temps sont rajoutés pour vérifier la rigidité globale du système.
Pour le calcul analytique du seuil on a:
Tant que la structure reste élastique et du fait des conditions aux limites, le tenseur des contraintes s’écrit:
Pour le premier pas de temps compris entre \(0\) et \(106s\)
, t en seconde et
en \(\mathrm{MPa}\) . Les autres composantes du tenseur sont nulles.
Pour les autres pas de temps,
Or on a:
Soit comme dans ce cas
, on obtient de façon immédiate en résolvant \(D=1\) la valeur du temps à partir duquel se déclare le fluage: \({t}_{1}=4.08181818{10}^{6}s\)
Ainsi pour un temps égal à \({t}_{1}\) , les déformations visqueuses sont nulles et \(D\) vaut 1.
Résultats de référence#
Variable interne \(\mathrm{V1}\) et \(\mathrm{V2}\) au point \(A\) , \(B\) , \(C\) et \(E\) ainsi que le déplacement du point \(B\) au dernier pas de temps
Références bibliographiques#
DE BONNIERES: Intégration des relations viscoélastiques dans STAT_NON_LINE [R5.03.08] février 2001
Modélisation A#
Caractéristiques de la modélisation#
Eléments 3D ( HEXA8 )
Un seul élément a été modélisé pour représenter le cube
C
D
z
y
A
B
x
H
G
E
F
Suivant l’axe \(Z\) : |
1 couche d’éléments |
Epaisseur totale : |
1 |
Conditions limites: |
|
Nœud \(A\) Nœud \(E\) Nœud \(D\) Nœud \(H\) |
\(\mathrm{DX}=\mathrm{DY}=\mathrm{DZ}=0.\) \(\mathrm{DX}=\mathrm{DY}=0.\) \(\mathrm{DX}=0.\) \(\mathrm{DX}=0.\) |
pression sur la face \(\mathrm{BCFG}\) |
\(F=220.\mathrm{MPa}\) |
Caractéristiques du maillage#
Nombre de nœuds: 8
Nombre de mailles et types: 1 HEXA8 et 1 QUAD4 (faces peau externe).
Grandeurs testées et résultats#
Localisation |
Grandeur |
Référence |
Aster |
% différence |
\(A\) |
\(\mathrm{V1}\) |
0.00000000 |
0.00000000 |
0% |
\(\mathrm{V2}\) |
1.00000000 |
9.9999999949239 10-01 |
-5.08 10-8% |
|
\(B\) |
\(\mathrm{V1}\) |
0.00000000 |
0.00000000 |
0% |
\(\mathrm{V2}\) |
1.00000000 |
9.9999999949239 10-01 |
-5.08 10-8% |
|
\(\mathrm{DX}\) |
/ |
7.8380536694423 10-2 |
/ |
|
\(\mathrm{DY}\) |
/ |
1.2984914071992 10-16 |
/ |
|
\(\mathrm{DZ}\) |
/ |
1.2984914071992 10-16 |
/ |
|
\(C\) |
\(\mathrm{V1}\) |
0.00000000 |
0.00000000 |
0% |
\(\mathrm{V2}\) |
1.00000000 |
9.9999999949239 10-01 |
-5.08 10-8% |
|
\(E\) |
\(\mathrm{V1}\) |
0.00000000 |
0.00000000 |
0% |
\(\mathrm{V2}\) |
1.00000000 |
9.9999999949239 10-01 |
-5.08 10-8% |
Synthèse des résultats#
L’objectif de ce cas test est complètement rempli puisque l’activation du fluage via la détection du seuil est très précise (de l’ordre de \({10}^{-8}\text{\%}\) d’erreur). Bien sûr, puisqu’il s’agit de phénomène visqueux, la discrétisation en temps joue un rôle très important notamment au alentour de l’activation du seuil. Dans ce cas test un soin tout particulier a été pris d’entourer le temps où le seuil est atteint de façon précise pour obtenir des résultats satisfaisants.