v1.01.135 ZZZZ135 - Validation informatique des méthodes SRM dans la macro-commande CALC_STAB_PENTE#
Résumé:
Ce test a pour objectif de valider les fonctionnalités informatiques de la macro-commande CALC_STAB_PENTE.
Dans ce cas test on étudie la stabilité d’un barrage en remblai avec l’infiltration d’eau. Ce test est un test de non-régression permettant de valider l’utilisation de la macro-commande, et donc n’a pas de sens physique ni de solution de référence scientifique.
Il y a 2 modélisations proposées pour ce cas test :
Modélisation A : Calcul du facteur de sécurité de la pente amont avec la loi Mohr-Coulomb. Cette modélisation valide les fonctionnalités suivantes :
L’utilisation des variables de commande dans le cham_mater en entrée.
La vérification de la loi de comportement dans la zone SRM.
L’utilisation de la loi Mohr-Coulomb dans l’algorithme SRM.
Le mot clé GROUP_MA permettant la définition de la zone SRM par les groupes des mailles indépendants.
Le mot clé CHAM_DEFO permettant de visualiser la surface de glissement.
Le mot clé METHODE = “EXPONENTIELLE”
Modélisation B : Calcul du facteur de sécurité du modèle entier avec la loi Drucker-Prager. Cette modélisation valide les fonctionnalités suivantes :
L’utilisation de la loi Drucker-Prager dans l’algorithme SRM.
Le mot-clé METHODE = “LINEAIRE”
Le mot-clé TOUT = “OUI”
Modélisation A#
Caractéristiques de la modélisation#
On calcule le facteur de sécurité de la pente amont du barrage en tenant compte de la présence du champ de pression interstitielle due à l’infiltration d’eau dans le barrage. L’effet mécanique sur le champ hydraulique est négligé. On résout donc le problème hydromécanique par chaînage, et introduit le champ de pression interstitielle dans le calcul mécanique par l’affectation de la variable de commande PTOT dans le cham_mater qui est entré dans la macro-commande CALC_STAB_PENTE.
Le barrage est modélisé en D_PLAN_HS pour le calcul hydraulique et en D_PLAN pour l’analyse de stabilité.
Fig. 480 Définition de la zone SRM#
En raison de la pression hydrostatique, la pente aval est normalement moins stable que la pente amont. Afin de calculer le facteur de sécurité de la pente amont, il faut définir une zone SRM qui exclut la pente en aval. La zone en rouge dans la est définie comme la zone SRM. Cette zone est définie par un groupe de maille non utilisé lors de l’affectation des matériaux, et donc permet de valider l’algorithme de traitement du maillage dans la macro-commande CALC_STAB_PENTE.
Afin de tester la vérification de la loi de comportement dans la zone SRM, on définit la relation de comportement élastique sur la recharge en aval. Cela démontre la possibilité de l’affectation des relations de comportement autres que Mohr-Coulomb et Drucker-Prager en dehors de la zone SRM.
On emploie la loi de variation exponentielle de l’incrément de facteur de sécurité,
avec INCR_INIT = 0,04 et RESI_MAXI = 0,02,
soit 1 raffinement en total avant la convergence de l’algorithme.
La surface de glissement est visualisée via le concept evol_noli produit.
Grandeurs testées et résultats#
On teste le facteur de sécurité et le déplacement maxi total à chaque itération SRM. Les valeurs de références sont montrées dans le Tableau 2.
On teste aussi sur le nœud nommé N1 le composant V1 (déformation plastique cumulée) du champ VARI_NOEU compris dans le concept evol_noli produit. La valeur de référence est montrée dans le Tableau 3.
Tableau 2 : Valeur de référence au nœud N1 du cham_defo
Identification |
Type de référence |
Valeur de référence |
Tolérance |
|---|---|---|---|
V1 |
“NON_REGRESSION” |
0.001287715823488581 |
0,0001% |
Tableau 3 : Valeurs de référence de la table des FS (Modélisation A)
Nb Itération SRM |
Identification |
Type de référence |
Valeur de référence |
Tolérance |
|---|---|---|---|---|
1 |
FS |
“NON_REGRESSION” |
1.84 |
0,0001% |
DISP_TOT_MAXI |
“NON_REGRESSION” |
0.32075657405211994 |
0,0001% |
|
2 |
FS |
“NON_REGRESSION” |
1.86 |
0,0001% |
DISP_TOT_MAXI |
“NON_REGRESSION” |
0.3230683648659598 |
0,0001% |
|
3 |
FS |
“NON_REGRESSION” |
1.87 |
0,0001% |
DISP_TOT_MAXI |
“NON_REGRESSION” |
0.32430147650733704 |
0,0001% |
Modélisation B#
Caractéristiques de la modélisation#
La modélisation B réutilise le champ de pression interstitielle de la modélisation A,
mais applique l’algorithme SRM au modèle entier avec la loi Drucker-Prager.
En outre, on emploie la loi de variation linéaire de l’incrément du facteur de sécurité,
avec INCR_INIT = 0,05, RESI_MAXI = 0,01, et ITER_RAFF_LINE= 1.
L’écart entre l’incrément initial et le résidu maximal étant grand, cette configuration des paramètres de convergence impose que l’algorithme converge à l’issue d’un raffinement de l’incrément.
Grandeurs testées et résultats#
On teste le facteur de sécurité et le déplacement maxi total à chaque itération SRM.
Les valeurs de référence sont montrées dans le Tableau 4.
Tableau 4 : Valeurs de référence de la table des FS (Modélisation B)
Nb Itération SRM |
Identification |
Type de référence |
Valeur de référence |
Tolérance |
|---|---|---|---|---|
1 |
FS |
“NON_REGRESSION” |
1.78 |
0,0001% |
DISP_TOT_MAXI |
“NON_REGRESSION” |
0.3196069268900407 |
0,0001% |
|
2 |
FS |
“NON_REGRESSION” |
1.83 |
0,0001% |
DISP_TOT_MAXI |
“NON_REGRESSION” |
0.4610095633405791 |
0,0001% |
|
3 |
FS |
“NON_REGRESSION” |
1.84 |
0,0001% |
DISP_TOT_MAXI |
“NON_REGRESSION” |
0.6337163549233418 |
0,0001% |