250 likes | 391 Views
Résultats partiels du benchmark GDR-MOMAS. René Chambon, Robert Charlier, Frédéric Collin, XiangLing Li, Abderrahim Lambarki (ENIM) Liège – Grenoble – SCK/CEN. Hypothèses générales (1/3). Conservation de la quantité de mouvement en quasistatique. Coefficient de Biot constant.
E N D
Résultats partiels du benchmark GDR-MOMAS René Chambon, Robert Charlier, Frédéric Collin, XiangLing Li, Abderrahim Lambarki (ENIM) Liège – Grenoble – SCK/CEN
Hypothèses générales (1/3) • Conservation de la quantité de mouvement en quasistatique Coefficient de Biot constant Petites déformations Grandes déformations Paris 23/11/2004 - GDRMOMAS2
Hypothèses générales (2/3) • Conservation de la masse d’eau Ancienne formulation Conservée ! Paris 23/11/2004 - GDRMOMAS3
Hypothèses générales (3/3) • Equations de bilan sous forme faible • Résolution : Newton-Raphson Paris 23/11/2004 - GDRMOMAS4
Cavité cylindrique (1D) (1/9) • Cavité cylindrique non-soutenue • Etat de contrainte isotrope Rayon intérieur 3 m • Dimensions du problème : Rayon extérieur 20 m Choix : Calcul unidimensionnel 72 éléments – 353 noeuds Paris 23/11/2004 - GDRMOMAS5
Cavité cylindrique (1D) (2/9) • Conditions initiales • Conditions limites et chargement Rext : Rint : T = 1.5 Ms (17 jours) Ttotal = 300 Ms (9.5 ans) Paris 23/11/2004 - GDRMOMAS6
Cavité cylindrique (1D) (3/9) • Modèle adoucissant avec perméabilité variable • Problème de convergence • Modifications proposées : Loi d’adoucissement Relation permabilité-porosité Condition de suintement Paris 23/11/2004 - GDRMOMAS7
Cavité cylindrique (1D) (4/9) • Elément de suintement pEnv Paris 23/11/2004 - GDRMOMAS8
Cavité cylindrique (1D) (5/9) • Elément de suintement Paris 23/11/2004 - GDRMOMAS9
Cavité cylindrique (1D) (6/9) • Elément de suintement Paris 23/11/2004 - GDRMOMAS10
Cavité cylindrique (1D) (7/9) • Elément de suintement Paris 23/11/2004 - GDRMOMAS11
Cavité cylindrique (1D) (8/9) • Comparaisons Paris 23/11/2004 - GDRMOMAS12
Cavité cylindrique (1D) (9/9) • Comparaisons Paris 23/11/2004 - GDRMOMAS13
Cavité cylindrique (2D) (1/12) • Cavité cylindrique non-soutenue • Etat de contrainte anisotrope Rayon intérieur 3 m • Dimensions du problème : Demi-longueur 60 m Choix : Calcul sur un quart de cavité (Symétrie supposée) 894 éléments – 2647 nœuds – 7941 ddl Maillage M1 Paris 23/11/2004 - GDRMOMAS14
Cavité cylindrique (2D) (2/12) 246 él. T = 1.5 Ms (17 jours) Ttotal = 300 Ms (9.5 ans) Paris 23/11/2004 - GDRMOMAS15
Cavité cylindrique (2D) (3/12) • Calcul couplé – Modèle plasticité parfaite Paris 23/11/2004 - GDRMOMAS16
Cavité cylindrique (2D) (4/12) • Calcul couplé – Modèle plasticité parfaite Paris 23/11/2004 - GDRMOMAS17
Cavité cylindrique (2D) (5/12) • Calcul couplé – Modèle plasticité parfaite – Perm. Var. Paris 23/11/2004 - GDRMOMAS18
Cavité cylindrique (2D) (6/12) • Calcul couplé – Modèle plasticité adoucissante Paris 23/11/2004 - GDRMOMAS19
Cavité cylindrique (2D) (7/12) • Calcul couplé – Modèle plasticité adoucissante Critère de bifurcation Indicateur de plasticité Paris 23/11/2004 - GDRMOMAS20
Cavité cylindrique (2D) (8/12) • Calcul couplé – Modèle plasticité adoucissante Incrément de déf. Déviatorique Déf. Déviatorique Paris 23/11/2004 - GDRMOMAS21
Cavité cylindrique (2D) (9/12) • Calcul couplé – Modèle plasticité adoucissante – Perm. Var. – suint. Paris 23/11/2004 - GDRMOMAS22
Cavité cylindrique (2D) (10/12) • Calcul couplé – Modèle plasticité adoucissante – Perm. Var. – suint. Critère de bifurcation Indicateur de plasticité Paris 23/11/2004 - GDRMOMAS23
Cavité cylindrique (2D) (11/12) • Calcul couplé – Modèle plasticité adoucissante – Perm. Var. – suint. Incrément de déf. Déviatorique Déf. Déviatorique Paris 23/11/2004 - GDRMOMAS24
Cavité cylindrique (2D) (12/12) Comparaisons : Paris 23/11/2004 - GDRMOMAS25