Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
1/16
Organisme(s) : EDF-R&D/MMC
Manuel de Validation
Fascicule V6.04 : Statique non-linéaire des structures volumiques
Document : V6.04.163
SSNV163 - Calcul de fluage propre
avec le modèle UMLV
Résumé :
Ce test permet de valider le modèle de fluage propre UMLV. Les résultats de ce test sont comparés avec la
solution analytique pour trois modélisations : 3D, axisymétrique et contraintes planes.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
2/16
1
Problème de référence
1.1 Géométrie
z
e
h
y
x
l
hauteur :
h = 1,00 [m]
largeur :
l = 1,00 [m]
épaisseur :
e = 1,00 [m]
1.2
Propriétés du matériau
E = 31 GPa
= ,
0 2
Ici on renseigne aussi la courbe sorption-désorption qui relie la teneur en eau C à l'hygrométrie h .
Dans ce cas on a supposé que les valeurs numériques de C et de h sont les mêmes.
Paramètres spécifiques au fluage propre :
k sr = ,
2 0E + 5 [MPa]
partie sphérique : rigidité apparente associée au squelette formé
par des blocs d'hydrates à l'échelle mésoscopique
k si = 0
,
5 E + 4 [MPa]
partie sphérique : rigidité apparente associée intrinsèquement
aux hydrates à l'échelle microscopique
k dr = ,
5 0E + 4 [MPa]
partie déviatorique : rigidité associée à la capacité de l'eau
adsorbée à transmettre des charges (load bearing water)
s
r = ,
4 0E +10 [MPa.s]
partie sphérique :viscosité apparente associée au mécanisme
de diffusion au sein de la porosité capillaire
s
i = 0
,
1 E +11 [MPa.s]
partie sphérique : viscosité apparente associée au mécanisme
de diffusion interlamellaire
d
r = ,
1 0E +10 [MPa.s]
partie déviatorique :viscosité associée à l'eau adsorbée par les
feuillets d'hydrates
d
i = 0
,
1 E +11 [MPa.s]
partie déviatorique : viscosité de l'eau libre.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
3/16
1.3
Conditions aux limites et chargements
Dans cet essai, on crée un champ de séchage homogène invariant dans la structure, l'humidité vaut
100% (condition d'une éprouvette scellée). Le chargement mécanique correspond à une compression
unidirectionnelle suivant le direction verticale (z en 3D ou y en 2D) ; son intensité est de 1 [MPa]. La
charge est appliquée en 1s et est maintenue constante pendant 100 jours.
1.4 Conditions
initiales
Le début du calcul est supposé l'instant 1. A cet instant il n'y a ni champ de séchage, ni contrainte
mécanique.
A l'instant 0, on applique un champ de séchage correspondant à 100 % d'hygrométrie.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
4/16
2
Solution de référence
2.1
Méthode de calcul
Cette section présente la résolution analytique complète du problème d'un corps d'épreuve soumis à
un champ de contraintes homogènes et unidirectionnelles appliquées instantanément à l'instant initial
et maintenues constantes par la suite (cas d'un essai de fluage en compression simple) :
= e
0 z e z
éq 2.1-1
Dont la décomposition en partie sphérique et déviatorique s'écrit :
= 1
2
1
01 + e
0
z ez -
0 (ex ex + e y e y )
éq
2.1-2
3 3
2
1
1
3
4
4
4
4
4
4
4
3 2
4
4
4
4
4
4
4
3
partie
dé
partie
viatorique
sphérique
En opérant une décomposition sphérique/déviatorique identique à celle des contraintes, la
déformation axiale s'écrit sous la forme :
= fs
zz
(
fd
0
)3+ (20 )3
éq
2.1-3
Il faut donc résoudre successivement la réponse à un échelon de contrainte sphérique et à un
échelon de contraintes déviatoriques.
2.2
Résolution des équations constitutives du fluage sphérique [bib2]
Le processus de déformation sphérique du fluage est gouverné par le système d'équations couplées
suivant (équations [éq 2.2-1] et [éq 2.2-2], cf. [R7.01.06]) :
fs
& = 1 h - k -
éq
2.2-1
s [
s
s
fs
r
r ]
fs
i
&
r
où s
kr désigne la rigidité apparente associée au squelette formé par des blocs d'hydrates à l'échelle
mésoscopique ;
et s
r la viscosité apparente associée au mécanisme de diffusion au sein de la porosité capillaire.
fs
1
i
& =
k
k
k
h
k
éq
2.2-1
s [ sr
fs - ( sr + si) fsi]-[ s - s fs
r
r ] +
i
où
s
ki désigne la rigidité apparente associée intrinsèquement aux hydrates à l'échelle
microscopique ;
et s
i la viscosité apparente associée au mécanisme de diffusion interfoliaire.
+
+
1
Dans [éq 2.2-2], les crochets
désignent l'opérateur de Mac Cauley : x
= (x + x )
2
La résolution du système d'équations couplées précédent nécessite de distinguer deux cas selon le
signe de la quantité comprise entre les crochets de Mac Cauley. Dans la suite, on présente la
résolution analytique de la réponse à un échelon de contrainte
s
. L'humidité relative est supposée
invariante ; le milieu est saturé en eau.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
5/16
2.2.1 Cas du fluage à court terme
A l'instant initial, t = 0, on applique une contrainte sphérique s
positive. Les déformations de fluage
réversibles et irréversibles sont égales à zéro (conditions initiales). L'équation du système [éq 2.2-2]
s'écrit donc :
fs
&
i (
+
+
t = )
1
0 =
k
k
éq
2.2.1-1
s [2 sr 0 - s
i -
s ]
1
0
=
[- s
s
] =0
i
i
La vitesse de déformation de fluage irréversible est donc égale à zéro. On en déduit que la
déformation de fluage irréversible est aussi égale à zéro. La vitesse de déformation irréversible reste
égale à zéro jusqu'à l'instant t = t0 , défini par la relation [éq 2.2.1-2] :
s
s
fs
2 kr r (t0)
s
fs
- = 0 r (t0) =
éq
2.2.1-2
s
2 kr
Jusqu'à l'instant t = t0 , la déformation de fluage réversible est définie par la relation suivante :
s
s
1
t
&r =
k
t
1 exp
éq
2.2.1-3
s [ s - sr sr ]
sr ( ) =
-
-
s
k
s
r
r
r
s
s
r
r =
est le temps caractéristique associé à la déformation de fluage réversible. L'instant t est
s
0
kr
donc défini par la relation [éq 2.2.1-4] :
fs
t
0
r (t0 )
s
s
=
=
1- exp -
t0 = ln(2) s
s
r 0.69 r éq
2.2.1-4
s
s
s
2 k
r
kr
r
Les déformations de fluage réversibles et irréversibles sont donc déterminées par :
s
fs (
t)
t
=
r
s
1- exp - s
k
éq
2.2.1-5
r
r
fs
(t)
= 0
i
Lors du calcul des déformations de fluage pour t > t0 , les nouvelles conditions initiales sont donc :
s
fs (
t0)
=
r
s
2 kr éq
2.2.1-6
fs
(t0)
= 0
i
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
6/16
2.2.2 Cas du fluage à long terme
En exprimant les vitesses de déformations de fluage réversibles et irréversibles en fonction des
déformations de fluage, on obtient alors la relation :
s
s
k
k
s
k
fs
r
r
fs
i
fs
1
2
& = -
- 4
r + 2
i
+
+
s
r
s
s
s
s
s
r
i
i
r i
éq 2.2.2-1
s
k
k
fs
r
fs
i
fs
1
&
2
i
=
r + -
i + -
s
s
s
s
i
i
i
Afin de simplifier les calculs, on définit les variables intermédiaires suivantes :
kr
1
ki
1
kr
urr =
:
=
, uii =
:
=
et uri =
:
éq
2.2.2-2
r r
i i
i
Le système d'équations [éq 2.2.2-1] peut se mettre alors sous la forme matricielle suivante :
-
- 4
2
fs
& fs
fs
u
u
u
1
u
+ 2u
&
r
rr
ri
ii
r
s
rr
ri
= =
+
éq
2.2.2-3
& fs
fs
s
2 u
- u
k
- u
i
1
4
4
4
ri2
4
4
4
3
ii
i
3
2
1
1r 4
4 2 4
4ri 3
A
s
B
C'est-à-dire :
fs
& = A fs
s
+ B éq
2.2.2-4
Supposons que la matrice A soit diagonalisable (cette propriété sera vérifiée par la suite) :
1
-
A = P D P où D désigne la matrice diagonale des valeurs propres de la matrice A , P
la matrice des vecteurs propres de la matrice A et
1
-
P la matrice inverse de la matrice P . En
effectuant terme à terme le produit par la quantité
1
-
P , [éq 2.2.2-4] peut se mettre sous la
forme :
fs,*
&
= D
fs,*
s
+ B*
fs,*
avec
= P-1
fs
et B* = P-1 B éq
2.2.2-5
Soit 1
et 2 les valeurs propres de la matrice A . On définit les quantités :
*
*
fs
b
,*
=
: 1
*
et B =
: 1
*
2
*2
b
[éq 2.2.2-5] s'écrit alors :
*
s
1& (t) =
1 *
1 (t ) +
*
b
1
éq
2.2.2-6
*
s
&2 (t) = 2 *
2 (t)+ *
b2
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
7/16
Système dont la solution s'écrit :
s
*
*
b
1 (t) = -
1 + 1
µ exp( 1
t)
1
1
0
éq
2.2.2-7
s
*
b
2 (t )
*
0
= -
2 + µ2
2
exp(2 t)
2
On peut alors revenir à l'espace initial, par le biais de la matrice de passage ; les déformations de
fluage réversibles et irréversibles sont des combinaisons linéaires de *
1 et *
2 . Les valeurs propres
de la matrice A , 1
et 2 sont obtenues en résolvant :
det(A - i )
1 = 0
- urr - 4 uri -
2
i
u
éq 2.2.2-8
ii
= 0
2
i + (urr + 4uri + uii ) i + urr uii = 0
2 uri
- uii - i
En remarquant que urr , uri et uii sont strictement positifs, le discriminant est donc toujours
strictement positif. Les valeurs propres sont donc réelles et distinctes, la matrice A est donc
diagonalisable. Par ailleurs, aucune des deux valeurs propres n'est égale à zéro
( = urr uii 0
1
2
). Les deux valeurs propres sont définies par :
- (urr + 4uri + uii )-
1
=
2
éq
2.2.2-9
- (urr + 4uri + uii )+
2 =
2
On peut montrer que les deux valeurs propres sont effectivement négatives. Montrons que la
deuxième valeur propre est négative. La déformation de fluage sphérique est donc asymptotique,
hypothèse émise dans le modèle de fluage propre sphérique [bib1]. Déterminons maintenant une
base des vecteurs propres (X 1, X 2 ) associés aux valeurs propres 1
et 2. Elle se détermine en
résolvant l'équation (A - i )
1 X i = 0 .
Une base particulière de vecteurs propres s'écrit :
x1
1
1
+ uii
2 uri
X 1 =
et X
2 =
avec x
1 =
et x2 =
éq
2.2.2-10
1
x
2
2 uri
2 + uii
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
8/16
Après avoir vérifié que P peut effectivement être inversée, on en déduit la solution dans l'espace
physique :
b*
b*
fs
1
2
r (t) =
s
x
-
1
+
+ x1 1
µ exp( 1
t)+ µ2 exp(2 t)
1
2
b*
b*
fs
1
2
i (t) =
s
x
-
+ 2
+ 1
µ exp( 1
t)+ x2 µ2 exp(2 t)
1
2
avec
éq 2.2.2-11
1
b*1 =
[x2 (urr + 2uri )+ uri ]
x1 x2 -1
1
b*2 =
[- (urr + 2uri )- x1 uri ]
x1 x2 -1
Enfin, 1
µ et µ2 sont définis par les relations :
1
1
1
µ1 = -
x2 exp 2 t0 - exp 2 t0
( 1x x2 - )1exp ([ 1
+ 2 )t0 ]
(
)
(
)
2 k
k
r
i
éq
2.2.2-12
1
1
1
µ2 = -
-
exp 1
t0 + 1
x exp 1
t0
( 1x x2 - )1exp ([ 1
+ 2 )t0 ]
(
)
(
)
2 k
k
r
i
2.3
Résolution des équations constitutives du fluage déviatorique
Les contraintes déviatoriques comportent une partie réversible et une partie irréversible
(cf. [R7.01.06]) :
fd
fd
fd
=
+
éq
2.3-1
{
{
r
{
i
contribution
déformation
contribution
eau
ue
déviatoriq
eau
absobée
totale
libre
La jème composante principale de la déformation déviatorique totale est régie par les équations
[éq 2.3-2] et [éq 2.3-3] :
d d , j
d d , j
d j
r r& + k
,
r r
= h
éq
2.3-2
où d
kr désigne la rigidité associée à la capacité de l'eau adsorbée à transmettre des charges (load
bearing water) ;
et d
r la viscosité associée à l'eau adsorbée par les feuillets d'hydrates.
d d , j
d j
i i& = h
,
éq 2.3-3
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
9/16
où
d
i désigne la viscosité de l'eau libre. Le système d'équations [éq 2.3-2] et [éq 2.3-3] est plus
simple à résoudre que celui régissant le comportement sphérique du fait qu'il est découplé. On
suppose toujours que l'humidité reste égale à 1 durant tout le chargement. L'équation [éq 2.3-2]
correspond au modèle visco-élastique de Kelvin dont le réponse à un échelon de contrainte est de
type exponentiel. Quant à l'équation [éq 2.3-3], la réponse en déformation est linéaire avec la temps.
La déformation de fluage totale s'écrit donc comme la somme de la contribution d'une chaîne de
Kelvin et de la contribution d'un amortisseur et série :
k d
r
-
t
d , j
t
1
d
(t) =
+
1
d ,
- e
j
r
H (t)
éq
2.3-4
d
k d
i
r
2.4
Récapitulatif de la solution analytique
Pour un chargement uniaxial les solutions analytiques des deux composantes de déformation sont
connues. La contribution de la partie déviatorique s'écrit :
d
2
t
1
k t
fd
(t)
r
= 0
+
1- exp -
éq
2.4-1
3
d
d
k
d
i
r
r
Quant à la contribution de la partie sphérique, la solution est définie sur deux intervalles :
s
k t
s
0
r
1- exp -
t r
s
s
ln 2
s
fs
(t) 3k
r
r
k
=
r
éq 2.4-2
s
0 1
1
r
+
+ µ1(1+ 1
x )exp( 1
t)+ µ2 (1+ x2 )exp(2t) t >
ln 2
3 s
s
s
k
k
r
i
kr
La déformation axiale est une fonction linéaire des deux contributions précédentes :
= fs
zz
(
fd
0
)3+ (20 )3
éq
2.4-3
2.5
Grandeurs et résultats de référence
L'essai est homogène. On teste la déformation en un noeud quelconque.
2.6
Incertitudes sur la solution
Résultat analytique exact.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
10/16
2.7 Références
bibliographiques
[1]
BENBOUDJEMA, F. : Modélisation des déformations différées du béton sous sollicitations
biaxiales. Application aux bâtiments réacteurs de centrales nucléaires, Mémoire de D.E.A.
Matériaux Avancés Ingénierie des Structures et des Enveloppes, 38 p. (+ annexes) (1999).
[2]
BENBOUDJEMA, F., MEFTAH, F., HEINFLING, G., LE PAPE, Y. : Étude numérique et
analytique de la partie sphérique du modèle de fluage propre UMLV pour le béton, note
technique HT-25/02/040/A, 56 p (2002).
[3]
LE PAPE, Y. : Relation de comportement UMLV pour le fluage propre du béton,
Documentation de Référence de Code_Aster [R7.01.06], 16 p (2002).
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
11/16
3 Modélisation
A
3.1
Caractéristiques de la modélisation
Modélisation 3D
z
NO7
NO6
NO8
NO5
NO2
NO3
y
NO4
NO1
x
3.2
Caractéristiques du maillage
Nombre de noeuds : 8
Nombre de mailles : 1 de type HEXA 8
6 de type QUAD 4
On définit les mailles suivantes :
S_ARR
NO3 NO7 NO8 NO4
S_AVT
NO1 NO2 NO6 NO5
S_DRT
NO1 NO5 NO8 NO4
S_GCH
NO3 NO2 NO6 NO7
S_INF
NO1 NO2 NO3 NO4
S_SUP
NO5 NO6 NO7 NO8
Les conditions aux limites en déplacement imposées sont :
Sur les noeuds NO1, NO2, NO3 et NO4 : DZ = 0
Sur les noeuds NO3, NO7, NO8 et NO4 : DY = 0
Sur les noeuds NO2, NO6, NO7 et NO8 : DX = 0
Le chargement est constitué du même champ de séchage et de la même force nodale 1/4 appliquée
sur les quatre noeuds de S_SUP.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
12/16
3.3 Fonctionnalités
testées
Commandes
Mot clé
DEFI_MATERIAU ELAS_FO
FONC_DESORP
UMLV_FP
K_RS
K_IS
K_RD
V_RS
V_IS
V_RD
V_ID
CREA_CHAMP `AFFE'
NOM_CMP='TEMP'
TYPE_CHAM='NOEU_TEMP_R'
AFFE_CHAR_MECA DDL_IMPO
LIAISON_UNIF
FORCE_NODALE
SECH_CALCULEE
STAT_NON_LINE COMP_INCR
RELATION=`UMLV_FP'
3.4
Grandeurs testées et résultats
La composante au noeud NO6 a été testée.
xx
Instant Référence
Aster %
différence
0. 0. 0.
-
1.0000E+00 3.225814D-05
3.225810D-05
1.37E-04
9.7041E+04 3.867143D-05
3.867140D-05
8.95E-05
1.8389E+06 6.088552D-05
6.088554D-05
3.25E-05
8.6400E+06 1.100478D-04
1.100473D-04
7.27E-06
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
13/16
4 Modélisation
B
4.1
Caractéristiques de la modélisation
Modélisation 2D axisymétrique.
y
N3
N4
N1
N2
x
4.2
Caractéristiques du maillage
Nombre de noeuds : 4
Nombre de mailles : 1 de type QUAD 4
4 de type SEG2
On définit les mailles suivantes :
L_INF NO1
NO2
L_DRT NO2
NO4
L_SUP NO4
NO3
L_GCH NO3
NO1
Les conditions aux limites en déplacement imposées sont :
Sur L_GCH : DY = 0
Sur L_INF : DX = 0
Le chargement est constitué du même champ de séchage et de la même force nodale 1/2 appliquée
sur les deux noeuds de L_SUP.
4.3 Fonctionnalités
testées
Commandes
Mot clé
DEFI_MATERIAU ELAS_FO
FONC_DESORP
UMLV_FP
K_RS
K_IS
K_RD
V_RS
V_IS
V_RD
V_ID
CREA_CHAMP `AFFE'
NOM_CMP='TEMP'
TYPE_CHAM='NOEU_TEMP_R'
AFFE_CHAR_MECA DDL_IMPO
LIAISON_UNIF
FORCE_NODALE
SECH_CALCULEE
STAT_NON_LINE COMP_INCR
RELATION=`UMLV_FP'
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
14/16
4.4
Grandeurs testées et résultats
La composante yy au noeud NO3 a été testée
Instant Référence
Aster %
différence
0. 0. 0.
-
1.0000E+00 3.225814D-05
3.225810D-05
1.37E-04
9.7041E+04 3.867143D-05
3.867140D-05
8.95E-05
1.8389E+06 6.088552D-05
6.088554D-05
3.25E-05
8.6400E+06 1.100478D-04
1.100473D-04
7.27E-06
5 Modélisation
C
5.1
Caractéristiques de la modélisation
Modélisation en Contraintes Planes.
y
N3
N4
N1
N2
x
5.2
Caractéristiques du maillage
Nombre de noeuds : 4
Nombre de mailles : 1 de type QUAD 4
4 de type SEG2
On définit les mailles suivantes :
L_INF NO1
NO2
L_DRT NO2
NO4
L_SUP NO4
NO3
L_GCH NO3
NO1
Les conditions aux limites en déplacement imposées sont :
Sur L_GCH : DY = 0
Sur L_INF : DX = 0
Le chargement est constitué du même champ de séchage et de la même force nodale 1/2 appliquée
sur les deux noeuds de L_SUP.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
15/16
5.3 Fonctionnalités
testées
Commandes
Mot clé
DEFI_MATERIAU ELAS_FO
FONC_DESORP
UMLV_FP
K_RS
K_IS
K_RD
V_RS
V_IS
V_RD
V_ID
CREA_CHAMP `AFFE'
NOM_CMP='TEMP'
TYPE_CHAM='NOEU_TEMP_R'
AFFE_CHAR_MECA DDL_IMPO
LIAISON_UNIF
FORCE_NODALE
SECH_CALCULEE
STAT_NON_LINE COMP_INCR
RELATION=`UMLV_FP'
ALGO_C_PLAN='DEBORST'
5.4
Grandeurs testées et résultats
La composante yy au noeud NO3 a été testée
Instant Référence
Aster %
différence
0. 0. 0.
-
1.0000E+00 3.225814D-05
3.225810D-05
1.40E-04
9.7041E+04 3.867143D-05
3.867140D-05
9.225E-05
1.8389E+06 6.088552D-05
6.088554D-05
3.08E-05
8.6400E+06 1.100478D-04
1.100478D-04
8.22E-06
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Code_Aster ®
Version
6.4
Titre :
SSNV163 - Calcul de fluage propre avec le modèle UMLV
Date :
01/10/03
Auteur(s) :
Y. LE PAPE Clé
:
V6.04.163-A Page :
16/16
6
Synthèse des résultats
Les valeurs obtenues avec le Code_Aster sont en accord avec les valeurs de la solution analytique de
référence. Ce même test a été tourné avec Castem au Laboratoire de Mécanique à L'Université de
Marne la Vallée, les mêmes résultats ont été obtenus.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-26/03/023/A
Document Outline