u4.90.21 Macro-commande CALC_ESSAI_GEOMECA#
Syntaxe#
Détail de la syntaxe
table = CALC_ESSAI_GEOMECA(
◆ MATER = mater,
◆ | COMPORTEMENT = _F(
◆ RELATION = / "HUJEUX",
/ "DRUCK_PRAGER",
/ "DRUCK_PRAG_N_A",
/ "CAM_CLAY",
/ "CJS",
/ "MOHR_COULOMB",
/ "Iwan",
/ "MohrCoulombAS",
/ "GonfElas",
/ "MCC",
/ "CSSM",
/ "NLH_CSRM",
/ "MFRONT",
# Si: equal_to('RELATION', 'MFRONT')
◆ COMPOR_MFRONT = compor_mgis,
◇ VERI_BORNE = / "ARRET" (par défaut),
/ "SANS",
/ "MESSAGE",
◇ SYME_MATR_TANG = / "OUI" (par défaut),
/ "NON",
# Si: not equal_to('RELATION', 'MFRONT')
◇ DEFORMATION = / "PETIT" (par défaut),
/ "GDEF_LOG",
/ "GREEN_LAGRANGE",
# Si: equal_to('RELATION', 'MFRONT')
◇ RESI_INTE = float,
◇ ITER_INTE_MAXI = int,
# Si: not equal_to('RELATION', 'MFRONT')
◇ RESI_INTE = float (défaut: 1e-06),
◇ ITER_INTE_MAXI = int (défaut: 20),
# Si: is_in('DEFORMATION', ('PETIT','PETIT_REAC','GROT_GDEP'))
◇ ITER_INTE_PAS = int,
◇ ALGO_INTE = / "ANALYTIQUE",
/ "SECANTE",
/ "DEKKER",
/ "NEWTON_1D",
/ "BRENT",
/ "NEWTON",
/ "NEWTON_RELI",
/ "NEWTON_PERT",
/ "RUNGE_KUTTA",
/ "SPECIFIQUE",
/ "SEMI_EXPLICITE",
/ "BASCULE_EXPLICITE",
/ "SANS_OBJET",
◇ TYPE_MATR_TANG = / "PERTURBATION",
/ "VERIFICATION",
# Si: (exists("TYPE_MATR_TANG"))
◇ VALE_PERT_RELA = float (défaut: 1e-05),
◇ PARM_THETA = float (défaut: 1.0),
# Si: not exists("TYPE_MATR_TANG")
◇ RESI_RADI_RELA = float,
),
⟐ CONVERGENCE = _F(
◆ | RESI_GLOB_MAXI = float,
| RESI_GLOB_RELA = float,
◇ ITER_GLOB_MAXI = int (défaut: 10),
),
◆ | ESSAI_TRIA_DR_M_D = _F(
◆ PRES_CONF = float,
◆ EPSI_IMPOSE = float,
◇ NB_INST = int (défaut: 100),
◇ KZERO = float (défaut: 1.0),
◇ TABLE_RESU = co,
◇ GRAPHIQUE = text (défaut: "('P-Q', 'EPS_AXI-Q', 'EPS_AXI-EPS_VOL', 'P-EPS_VOL')"),
◇ PREFIXE_FICHIER = text,
◇ TABLE_REF = table,
◇ NOM_CMP = text,
◇ COULEUR = int,
◇ MARQUEUR = int,
◇ STYLE = int,
),
| ESSAI_TRIA_ND_M_D = _F(
◆ PRES_CONF = float,
◆ EPSI_IMPOSE = float,
◇ BIOT_COEF = float (défaut: 1.0),
◇ KZERO = float (défaut: 1.0),
◇ NB_INST = int (défaut: 100),
◇ TABLE_RESU = co,
◇ GRAPHIQUE = text (défaut: "('P-Q', 'EPS_AXI-Q', 'EPS_AXI-PRE_EAU')"),
◇ PREFIXE_FICHIER = text,
◇ TABLE_REF = table,
◇ NOM_CMP = text,
◇ COULEUR = int,
◇ MARQUEUR = int,
◇ STYLE = int,
),
| ESSAI_CISA_DR_C_D = _F(
◆ PRES_CONF = float,
◆ GAMMA_IMPOSE = float,
◇ GAMMA_ELAS = float (défaut: 1e-07),
◆ NB_CYCLE = int,
◇ KZERO = float (défaut: 1.0),
◇ NB_INST = int (défaut: 25),
◇ TABLE_RESU = co,
◇ GRAPHIQUE = text (défaut: "('GAMMA-SIG_XY', 'GAMMA-G_SUR_GMAX', 'GAMMA-DAMPING', 'G_SUR_GMAX-DAMPING')"),
◇ PREFIXE_FICHIER = text,
◇ TABLE_REF = table,
◇ NOM_CMP = text,
◇ COULEUR_NIV1 = int,
◇ MARQUEUR_NIV1 = int,
◇ STYLE_NIV1 = int,
◇ COULEUR_NIV2 = int,
◇ MARQUEUR_NIV2 = int,
◇ STYLE_NIV2 = int,
◇ TYPE_CHARGE = / "SINUSOIDAL" (par défaut),
/ "TRIANGULAIRE",
),
| ESSAI_TRIA_ND_C_F = _F(
◆ PRES_CONF = float,
◆ SIGM_IMPOSE = float,
◇ BIOT_COEF = float (défaut: 1.0),
◇ KZERO = float (défaut: 1.0),
◆ CRIT_LIQUEFACTION = / "RU_MAX",
/ "EPSI_ABSO_MAX",
/ "EPSI_RELA_MAX",
◆ VALE_CRIT = float,
◇ ARRET_LIQUEFACTION = / "OUI" (par défaut),
/ "NON",
◆ UN_SUR_K = float,
◆ NB_CYCLE = int,
◇ NB_INST = int (défaut: 25),
◇ NB_INST_MONO = int (défaut: 400),
◇ TABLE_RESU = co,
◇ GRAPHIQUE = text (défaut: "('NCYCL-DSIGM', 'P-Q', 'SIG_AXI-PRE_EAU', 'EPS_AXI-PRE_EAU', 'SIG_AXI-RU', 'EPS_AXI-RU', 'EPS_AXI-Q')"),
◇ PREFIXE_FICHIER = text,
◇ TABLE_REF = table,
◇ NOM_CMP = text,
◇ COULEUR_NIV1 = int,
◇ MARQUEUR_NIV1 = int,
◇ STYLE_NIV1 = int,
◇ COULEUR_NIV2 = int,
◇ MARQUEUR_NIV2 = int,
◇ STYLE_NIV2 = int,
◇ TYPE_CHARGE = / "SINUSOIDAL" (par défaut),
/ "TRIANGULAIRE",
),
| ESSAI_TRIA_ND_C_D = _F(
◆ PRES_CONF = float,
◆ EPSI_MINI = float,
◆ EPSI_MAXI = float,
◇ EPSI_ELAS = float (défaut: 1e-07),
◇ NB_INST = int (défaut: 25),
◆ NB_CYCLE = int,
◇ BIOT_COEF = float (défaut: 1.0),
◇ KZERO = float (défaut: 1.0),
◆ UN_SUR_K = float,
◇ RU_MAX = float (défaut: 0.8),
◇ TABLE_RESU = co,
◇ GRAPHIQUE = text (défaut: "('NCYCL-DEPSI', 'DEPSI-RU_MAX', 'DEPSI-E_SUR_EMAX', 'DEPSI-DAMPING', 'P-Q', 'EPS_AXI-EPS_VOL', 'EPS_AXI-Q', 'P-EPS_VOL', 'EPS_AXI-PRE_EAU', 'EPS_AXI-RU', 'P-PRE_EAU')"),
◇ PREFIXE_FICHIER = text,
◇ TABLE_REF = table,
◇ NOM_CMP = text,
◇ COULEUR_NIV1 = int,
◇ MARQUEUR_NIV1 = int,
◇ STYLE_NIV1 = int,
◇ COULEUR_NIV2 = int,
◇ MARQUEUR_NIV2 = int,
◇ STYLE_NIV2 = int,
◇ TYPE_CHARGE = / "SINUSOIDAL" (par défaut),
/ "TRIANGULAIRE",
),
| ESSAI_TRIA_DR_C_D = _F(
◆ PRES_CONF = float,
◆ EPSI_MINI = float,
◆ EPSI_MAXI = float,
◇ EPSI_ELAS = float (défaut: 1e-07),
◇ NB_INST = int (défaut: 25),
◆ NB_CYCLE = int,
◇ KZERO = float (défaut: 1.0),
◇ TABLE_RESU = co,
◇ GRAPHIQUE = text (défaut: "('P-Q', 'EPS_AXI-Q', 'EPS_VOL-Q', 'EPS_AXI-EPS_VOL', 'P-EPS_VOL', 'DEPSI-E_SUR_EMAX', 'DEPSI-DAMPING')"),
◇ PREFIXE_FICHIER = text,
◇ TABLE_REF = table,
◇ NOM_CMP = text,
◇ COULEUR_NIV1 = int,
◇ MARQUEUR_NIV1 = int,
◇ STYLE_NIV1 = int,
◇ COULEUR_NIV2 = int,
◇ MARQUEUR_NIV2 = int,
◇ STYLE_NIV2 = int,
◇ TYPE_CHARGE = / "SINUSOIDAL" (par défaut),
/ "TRIANGULAIRE",
),
| ESSAI_OEDO_DR_C_F = _F(
◆ PRES_CONF = float,
◆ SIGM_IMPOSE = float,
◆ SIGM_DECH = float,
◇ KZERO = float (défaut: 1.0),
◇ NB_INST = int (défaut: 25),
◇ TABLE_RESU = co,
◇ GRAPHIQUE = text (défaut: "('P-EPS_VOL', 'SIG_AXI-EPS_VOL')"),
◇ PREFIXE_FICHIER = text,
◇ TABLE_REF = table,
◇ NOM_CMP = text,
◇ COULEUR = int,
◇ MARQUEUR = int,
◇ STYLE = int,
◇ TYPE_CHARGE = / "SINUSOIDAL" (par défaut),
/ "TRIANGULAIRE",
),
| ESSAI_ISOT_DR_C_F = _F(
◆ PRES_CONF = float,
◆ SIGM_IMPOSE = float,
◆ SIGM_DECH = float,
◇ NB_INST = int (défaut: 25),
◇ TABLE_RESU = co,
◇ GRAPHIQUE = text (défaut: "P-EPS_VOL"),
◇ PREFIXE_FICHIER = text,
◇ TABLE_REF = table,
◇ NOM_CMP = text,
◇ COULEUR = int,
◇ MARQUEUR = int,
◇ STYLE = int,
◇ TYPE_CHARGE = / "SINUSOIDAL" (par défaut),
/ "TRIANGULAIRE",
),
◇ INFO = / 1 (par défaut),
/ 2,
)
◆ : obligatoire
◇ : optionnel
⟐ : présent par défaut
& : ensemble
/ : un seul parmi
| : plusieurs choix possibles
Opérandes#
OpérandeMATER#
♦ MATER = mat [mater]
Ce mot-clé permet de renseigner le nom du matériau défini par DEFI_MATERIAU[U4.43.01], où sont fournis les paramètres nécessaires au comportement choisi.
Mot-cléCOMPORTEMENT#
La syntaxe de ce mot clé est décrite dans le document [U4.51.11].
Dans le cadre de cette macro-commande, l’utilisation de l’opérande RELATIONdu mot-clé COMPORTEMENTest bornée aux lois élastoplastiques de sol suivantes:
“MOHR_COULOMB”
“MohrCoulombAS”
“CAM_CLAY”
“CJS”
“DRUCK_PRAGER”
“DRUCK_PRAG_N_A”
“HUJEUX”
“Iwan”
“MFRONT”
“NLH_CSRM”
“MCC”
“CSSM”
Mot clé CONVERGENCE#
◊CONVERGENCE= _F()
Si aucun des deux opérandes suivants n’est présent, alors tout se passe comme si: RESI_GLOB_RELA = 1.E-6
Opérande RESI_GLOB_RELA/RESI_GLOB_MAXI#
◊| RESI_GLOB_RELA = resi_rela [R]
L’algorithme continue les itérations globales tant que:
\(\underset{i=1,\mathit{nbddl}}{max}|{F}_{i}^{n}|\phantom{\rule{2em}{0ex}}>\phantom{\rule{2em}{0ex}}\text{resi\_rela}\times max|L|\)
où \({F}^{n}\) est le résidu de l’itération \(n\) et \(L\) le vecteur du chargement imposé et des réactions d’appuis (Cf. [R5.03.01] pour plus de détails).
Lorsque le chargement et les réactions d’appui deviennent nuls, c’est-à-dire lorsque \(L\) est nul (par exemple dans le cas d’une décharge totale), on essaie de passer du critère de convergence relatif RESI_GLOB_RELA au critère de convergence absolu RESI_GLOB_MAXI. Cette opération est transparente pour l’utilisateur (message d’alarme émis dans le fichier .mess). Lorsque le vecteur \(L\) redevient différent de zéro, on repasse automatiquement au critère de convergence relatif RESI_GLOB_RELA.
Toutefois, ce mécanisme de basculement ne peut pas fonctionner au premier pas de temps. En effet, pour trouver une valeur de RESI_GLOB_MAXI raisonnable de manière automatique (puisque l’utilisateur ne l’a pas renseigné), on a besoin d’avoir eu au moins un pas convergé sur le mode RESI_GLOB_RELA . Dès lors, si le chargement est nul dès le premier instant, le calcul s’arrête. L’utilisateur doit déjà alors vérifier que le chargement nul est normal du point de vue de la modélisation qu’il réalise, et si tel est le cas, trouver un autre critère de convergence (RESI_GLOB_MAXI par exemple).
Si cet opérande est absent, le test est effectué avec la valeur par défaut, sauf si RESI_GLOB_MAXI est présent.
◊| RESI_GLOB_MAXI = resi_maxi [R]
L’algorithme continue les itérations globales tant que:
\(\underset{i=1,\mathit{nbddl}}{max}|{F}_{i}^{n}|\phantom{\rule{2em}{0ex}}>\phantom{\rule{2em}{0ex}}\text{resi\_maxi}\)
où \({F}^{n}\) est le résidu de l’itération \(n\) (cf. [R5.03.01] pour plus de détails). Si cet opérande est absent, le test n’est pas effectué.
Si RESI_GLOB_RELA et RESI_GLOB_MAXI sont présents tous les deux, les deux tests sont effectués.
Opérande ITER_GLOB_MAXI#
◊ITER_GLOB_MAXI = | 10 [DEFAUT]
| nb_iter
Nombre d’itérations maximum effectué pour résoudre le problème global à chaque instant.
MotcléESSAI_TRIA_DR_M_D#
Ce mot-clé facteur (répétable) permet d’effectuer une série de simulations d’un même essai triaxial drainé monotone à déformation imposée pour lesquelles on fait varier les paramètres de chargement (pression de confinement et déformation axiale imposée), de post-traiter les résultats obtenus et de les écrire sous forme de graphiques (au format xmgrace) et/ou de tables.
Convention de signe des entrées et sorties#
La convention de signe des géomécaniciens s’applique aux paramètres d’entrée en contraintes ou déformations imposées, c’est-à-dire que les valeurs sontpositivesen compression.
Cette convention s’applique de même sur les variables de sorties pré-définies dont la liste complète pour cet essai est la suivante:
INST: instant
EPS_AXI: déformation axiale
EPS_LAT: déformation latérale
EPS_VOL: déformation volumique
SIG_AXI: contrainte effective axiale
SIG_LAT: contrainte effective latérale
P: contrainte effective moyenne
Q: déviateur des contraintes
En revanche, cette convention ne s’applique pas sur les variables non pré-définies demandées en sortie dans NOM_CMP(§ 3.4.6 ).
Opérandes PRES_CONF, EPSI_IMPOSE, NB_INST#
♦ PRES_CONF = l_sigma_conf [l_R]
♦ EPSI_IMPOSE = l_epsi_impo [l_R]
◊NB_INST = | 100 [DEFAUT]
| nbinst [I]
L’opérande PRES_CONF permet de définir la liste des pressions de confinement qui seront maintenuesau coursde chaque essai. De même l’opérande EPSI_IMPOSE permet de définir la liste des valeurs finales du chargement de compression (rampe de déformation axiale imposée).
Pour cet essai, on fait correspondre à chaque pression de confinement une valeur finale pour larampe dedéformation axiale (voir figure ): les listes PRES_CONF et EPSI_IMPOSE doivent donc avoir même cardinal. Ce cardinal correspond au nombre de simulations qui seront exécutées sous ce mot-clé facteur.
Les contraintes et les déformations étant comptées positivement en compression, les valeurs renseignées pour PRES_CONF doivent être strictement positives. Les valeurs renseignées dans EPSI_IMPOSE sont strictement positives ou négatives, une valeur négative indiquant un chargement en extension.
L’opérande NB_INSTpermet de définir la discrétisation temporelle du chargement (voir figure ), avec une valeur par défaut de 100 pas de chargement au cours de la rampe.
.. image:: images/Cadre5.gif
.. _RefSchema_Cadre5.gif:
Figure 3.4.2-a: discrétisation et allure du chargement pour les mots-clés ESSAI_TD et ESSAI_TND
Opérande KZERO#
◊KZERO = | 1 [DEFAUT]
Valeur du coefficient des terres au repos, permet de définir un état de confinement anisotrope: \({\sigma}_{xx}={\sigma}_{yy}={K}_{0}{\sigma}_{zz}={K}_{0}\times \text{PRES\_CONF}\)
Remarque: Lorsque la valeur de KZEROest renseignée différente de 1, la pression de confinement réelle de l’essai n’est plus PRES_CONF, elle devient:
\({P}_{c}=\frac{1+2{K}_{0}}{3}\times \text{PRES\_CONF}\)
Opérande TABLE_RESU#
◊TABLE_RESU= l_tabres [l_CO]
Cet opérande facultatif permet de donner la liste des noms des concepts produits par la macro-commande qui seront alors de type [table]. Chaque table produite contient les résultats bruts et post-traités d’une simulation d’essai: la liste TABLE_RESU doit donc avoir même cardinal que les listes PRES_CONF et EPSI_IMPOSE.
Le titre de chaque table produite est complété par la macro-commande, il comprend:
le nom du mot-clé facteur (ici ESSAI_TRIA_DR_M_D) et son numéro d’occurrence (celui-ci étant répétable);
le couple de valeurs (PRES_CONF, EPSI_IMPOSE)caractérisant lechargement del’essai;
Exemple :
TABRES1=CO(“TRES1”)
TABRES2=CO(“TRES2”)
TABRES3=CO(“TRES3”)
CALC_ESSAI_GEOMECA(
…
ESSAI_TRIA_DR_M_D = _F(
PRES_CONF = (-1.E4, -1.5E4, -2.E4 ),
EPSI_IMPOSE = (-1.E-2, -1.E-2, -1.E-2),
TABLE_RESU = (TABRES1, TABRES2, TABRES3),),
);
TABRES1, TABRES2, et TABRES3sont successivement remplies selon l’ordre des listes PRES_CONF et EPSI_IMPOSE, le tableau ci-dessous précise les résultats d’essais contenus dans chaque table.
EPSI_IMPOSE ® |
-1.0E-2 |
-1.0E-2 |
-1.0E-2 |
PRES_CONF |
|||
-1.0E4 |
TABRES1 |
||
-1.5E4 |
TABRES2 |
||
-2.0E4 |
TABRES3 |
Opérande GRAPHIQUE, PREFIXE_FICHIER#
◊GRAPHIQUE = | ('P-Q', 'EPS_AXI-Q', 'EPS_AXI-EPS_VOL',
'P-EPS_VOL') [DEFAUT]
| l_graphique [l_Kn]
Cet opérande permet de préciser les types des graphiques produits par la macro-commande. Ces graphiques récapitulent les résultats des simulations exécutées sous le mot-clé facteur courant. Une liste par défaut est proposée. Cependant, il est possible de remonter un graphique quelconque en donnant les combinaisons souhaitées de composantes existantes dans l_graphique, de la forme: ABSC-ORDO. Le tiret «-» permet de séparer la composante demandée en abscisse de celle demandée en ordonnée.
Exemple :
l_graphique = (“P-Q”, “EPS_AXI-Q”, “INST-EPS_VOL”, “INST-V23”)
La composante non-prédéfinie V23 (correspondant à la déformation volumique plastique de la loi de Hujeux [R7.01.23]) doit être déclarée dans NOM_CMP ( § 3.4.6 ) pour être prise en compte. Dans le cas contraire, le graphique demandé n’est pas remonté.
Quant aux composantes pré-définies dont la liste complète est donnée dans § 3.4.1 , elles seront affichées éventuellement avec la convention de signe des géomécaniciens, c’est-à-dire avec des valeurs positives en compression.
Par ailleurs, toutes les composantes demandées pour la construction graphique seront affichées dans les tableaux de sorties ( § 3.4.4 ).
Les fichiers contenant ces graphiques sont écrits au format xmgrace dans un même répertoire spécifié par l’utilisateur (type repeen résultat dans astk), et sont nommés de la manière suivante:
“préfixe’_’nom_essai’_’numéro_occurrence’_’absc-ordo”.agr
Un préfixe peut être ajouté au nom du fichier de sortie grâce au mot-clé facultatif:
◊PREFIXE_FICHIER= préfixe [Kn]
Opérande NOM_CMP#
◊NOM_CMP= l_composante [l_Kn]
liste des composantes non-prédéfinies demandées en sortie. Il s'agit de toutes les composantes existantes contenues dans les champs SIGM, EPSI et VARI produits par le calcul. Ces composantes non-prédéfinies seront produites avec la convention de signe par défaut, c'est-à-dire avec des valeurs négatives en compression.
Les composantes non existantes sont ignorées.
Opérande TABLE_REF#
◊ TABLE_REF = l_tabref [l_table]
Cet opérande permet de renseigner des courbes de référence (par exemple, expérimentales) tabulées et stockées sous forme de tables, afin de les superposer aux courbes issues des simulations exécutées sous le mot-clé facteur courant. Ces courbes de référence sont alors incluses dans les fichiers produits par le mot-clé GRAPHIQUE.
Chaque table contenue dans la liste TABLE_REF doit être créée au préalable à l’aide de l’opérateur CREA_TABLE [U4.33.02], et formatée de la manière suivante:
tabref = CREA_TABLE(
LISTE = ( _F(PARA = “TYPE”, LISTE_K = [typgraph,]),
_F(PARA = “LEGENDE”, LISTE_K = [malegend,]),
_F(PARA = “ABSCISSE”, LISTE_R = l_absc),
_F(PARA = “ORDONNEE”, LISTE_R = l_ordo),),);
avec:
typgraph une chaîne de caractères dont la valeur appartient obligatoirement à la liste de valeurs par défaut du mot-clé GRAPHIQUE. Cette valeur permet d’identifier le type de graphique (et donc le fichier) auquel la courbe de référence doit être ajoutée.
malegend une chaîne de caractères qui contient la légende associée à la courbe de référence
l_absc et l_ordo sont des listes python de réels contenant respectivement les abscisses et les ordonnées des points de la courbe de référence. Ces listes doivent donc avoir même cardinal
Opérandes COULEUR, MARQUEUR, STYLE#
◊ COULEUR= l_couleur [l_I]
◊MARQUEUR= l_marqueur [l_I]
◊STYLE = l_style [l_I]
Ces opérandes acceptent une liste d’entiers permettant de définir respectivement la couleur, le type de marqueur et le style des courbes affichées dans les graphiques. L’information codée par chaque entier est donnée dans la documentation sur XMGRACE U2.51.01. La longueur de la liste doit être égale à celle de PRES_CONF.
MotcléESSAI_TRIA_ND_M_D#
Ce mot-clé facteur (répétable) permet d’effectuer une série de simulations d’un même essai triaxial non drainé monotone à déformation imposée (on suppose la saturation totale) pour lesquelles on fait varier les paramètres de chargement (pression de confinement et déformation axiale imposée), de post-traiter les résultats obtenus et de les écrire sous forme de graphiques (au format xmgrace) et/ou de tables.
Convention de signe des entrées et sorties#
La convention de signe des géomécaniciens s’applique aux paramètres d’entrée en contraintes ou déformations imposées, c’est-à-dire que les valeurs sontpositivesen compression.
Cette convention s’applique de même sur les variables de sorties pré-définies dont la liste complète est la suivante:
INST: instant
EPS_AXI: déformation axiale
EPS_LAT: déformation latérale
EPS_VOL: déformation volumique
SIG_AXI: contrainte effective axiale
SIG_LAT: contrainte effective latérale
P: contrainte effective moyenne
Q: déviateur des contraintes
PRE_EAU: pression interstitielle
En revanche, cette convention ne s’applique pas sur les variables non pré-définies demandées en sortie dans NOM_CMP(§ 3.5.7 ).
Opérandes PRES_CONF, EPSI_IMPOSE, NB_INST#
♦ PRES_CONF = l_sigma_conf [l_R]
♦ EPSI_IMPOSE = l_epsi_impo [l_R]
◊NB_INST = | 100 [DEFAUT]
| nbinst [I]
Idem qu’au § 3.4.2 .
Opérande BIOT_COEF#
◊BIOT_COEF = | 1 [DEFAUT]
| biot [R]
Valeur du coefficient de Biot.
Opérande KZERO#
◊KZERO = | 1 [DEFAUT]
Idem qu’au § 3.4.3 .
Opérande TABLE_RESU#
◊TABLE_RESU= l_tabres [l_CO]
Idem qu’au § 3.4.4 .
Opérande GRAPHIQUE, PREFIXE_FICHIER#
◊GRAPHIQUE = |('P-Q', 'EPS_AXI-Q',
'EPS_AXI-PRE_EAU') [DEFAUT]
| l_graphique [l_Kn]
◊PREFIXE_FICHIER = préfixe [Kn]
Idem qu’au § 3.4.5 .
Opérande NOM_CMP#
◊NOM_CMP= l_composante [l_Kn]
Idem qu’au § 3.4.6 .
Opérande TABLE_REF#
◊ TABLE_REF = l_tabref [l_table]
Idem qu’au § 3.4.7 .
Opérandes COULEUR, MARQUEUR, STYLE#
◊ COULEUR= l_couleur [l_I]
◊MARQUEUR= l_marqueur [l_I]
◊STYLE = l_style [l_I]
Idem qu’au § 3.4.8 .
MotcléESSAI_CISA_DR_C_D#
Ce mot-clé facteur (répétable) permet d’effectuer une série de simulations d’un même essai de cisaillement cyclique drainé à déformation de cisaillement imposée pour lesquelles on fait varier les paramètres de chargement (pression de confinement, amplitude de déformation de cisaillement et nombre de cycles), de post-traiter les résultats obtenus et de les écrire sous forme de graphiques (au format xmgrace) et/ou de tables.
Convention de signe des entrées et sorties#
La convention de signe des géomécaniciens s’applique aux paramètres d’entrée en contraintes ou déformations imposées, c’est-à-dire que les valeurs sontpositivesen compression.
Cette convention s’applique de même sur les variables de sortie pré-définies. En revanche, cette convention ne s’applique pas sur les variables non pré-définies demandées en sortie dans NOM_CMP(§ 3.6.7 ).
Dans cet essai, on distingue les variables de sortiepré-définies de niveau 2 qui correspondent aux variables produites par l’essaidont la liste complète est la suivante:
INST: instant
GAMMA: déformation de cisaillement \(\gamma =2{{\varepsilon}}_{xy}\)
SIG_XY: contrainte de cisaillement
Etles variables de sortie pré-définies de niveau 1 qui correspondent à des courbes dont les points représentent le résultat d’un essai (par ex. la courbe :math:`frac{G}{{G}_{max}}-gamma ` ). La liste complète des variables de niveau 1 est la suivante:
G_SUR_GMAX: \(\frac{G}{{G}_{\max}}\)
DAMPING: amortissement hystérétique \(\frac{\Delta W}{\pi W}\)
Opérandes PRES_CONF, GAMMA_IMPOSE, NB_CYCLE, NB_INST, TYPE_CHARGE#
♦ PRES_CONF = l_sigma_conf [l_R]
♦ GAMMA_IMPOSE = l_gamma_impo [l_R]
♦ NB_CYCLE = nbcyc [I]
◊ NB_INST = | 25 [DEFAUT]
|nbinst [I]
◊ TYPE_CHARGE = | 'SINUSOIDAL' [DEFAUT]
|'TRIANGULAIRE' [Kn]
Ces opérandes permettent de définir le chargement de chacune des simulations à exécuter sous le mot-clé facteur courant, ainsi que sa discrétisation. Leur signification est résumée à la figure et détaillée ci-dessous:
PRES_CONF permet de définir la liste des pressions de confinement (strictement positives) qui seront maintenues au cours de chaque essai;
GAMMA_IMPOSE permet de définir la liste des amplitudes (strictement positives) de déformation de cisaillement \(\gamma =2{\epsilon}_{xy}\) du chargement cyclique imposé;
NB_CYCLEcorrespond au nombre de cycles, fixé pour toutes les simulations.
NB_INSTpermet de définir la discrétisation temporelle du chargement, et correspond au nombre de pas de chargement par quart de cycle
TYPE_CHARGE indique le type de chargement souhaité: sinusoïdal ou triangulaire;
Pour chaque pression de confinement PRES_CONF, on effectue autant de simulations qu’il y a d’éléments dans la liste GAMMA_IMPOSE. Contrairement aux essais TRIA_DR_M_Det TRIA_ND_M_D(voir respectivement § 3.4 et § 3.5 ), ces listes ne sont pas en bijection et il y a au total
\(\mathit{card}(\text{PRES\_CONF})\times \mathit{card}(\text{GAMMA\_IMPOSE})\) simulations exécutées.
Figure 3.6.2-a : discrétisation et allure du chargement de type triangulaire pour le mot-clé ESSAI_CISA_DR_C_D pour 3 cycles de chargement
Opérande GAMMA_ELAS#
◊GAMMA_ELAS = | 1.E-7 [DEFAUT]
|gamma_elas [R]
Pour chaque pression de confinement, le module de cisaillement sécant maximal (c’est à dire du matériau sain) est déterminé en simulant une rampe de déformation de cisaillement imposée (en termes de distorsion) dont la valeur finale est GAMMA_ELAS. Cette valeur doit être telle que le matériau reste dans son domaine d’élasticité (linéaire ou non, selon la relation de comportement utilisée). GAMMA_ELAS vaut 1.E-7 par défaut, et toute valeur renseignée par l’utilisateur doit lui être inférieure. Si la valeur renseignée ne permet pas de rester dans le domaine d’élasticité, le code s’arrête en erreur fatale.
Opérande KZERO#
◊KZERO = | 1 [DEFAUT]
Idem qu’au § 3.4.3 .
Opérande TABLE_RESU#
◊TABLE_RESU= l_tabres [l_CO]
Cet opérande facultatif permet de donner la liste des noms des concepts produits par la macro-commande qui seront alors de type [table]. La taille de cette liste doit vérifier:
\(\mathit{card}(\text{TABLE\_RESU})=\mathit{card}(\text{PRES\_CONF})+1\)
En effet, chaque table produite regroupe les résultats bruts de toutes les simulations exécutées pour une même pression de confinement (PRES_CONF), dans laquelle chaque simulation correspond à un paquet de colonnes contiguës dont les titres sont tous indexés par le même entier (indice de la valeur considérée dans la liste GAMMA_IMPOSE). Une table supplémentaire récapitulant les post-traitements réalisés à l’issue de toutes les simulations est également produite. Cette table contient pour chaque pression de confinement (PRES_CONF) les valeurs du module de cisaillement sécant normalisé \(\frac{G}{{G}_{max}}\) et du taux d’amortissement hystérétique \(\frac{\Delta W}{\pi W}\) en vis-à-vis des amplitudes de distorsion imposée (GAMMA_IMPOSE). L’amortissement hystérétique \(\frac{\Delta W}{\pi W}\) pour le dernier cycle simulé est calculé d’après Kohusho [] de la façon suivante (Figure ):
:math:`Delta W={int}_{C}delta {sigma}_{xy}delta gamma ` l’aire de la dernière boucle d’hystérésis;
:math:`W=Delta {sigma}_{xy}Delta gamma ` l’énergie élastique associée
Cette table correspond au nom de concept donné en dernière position dans la liste TABLE_RESU. Des extraits de ces tables sont présentés dans l’exemple ci-dessous.
Figure 3.6.5-a: Définition de l’amortissement hystérétique d’après Kokusho []
Exemple :
TABRES1 = CO(“TRES1”)
TABRES2= CO(“TRES2”)
TABBILA= CO(“TBILA”)
CALC_ESSAI_GEOMECA(
…
ESSAI_CISA_DR_C_D= _F(
PRES_CONF = (1.E5, 2.05E5,),
GAMMA_IMPOSE = (1.E-5, 5.E-5, 1.E-4, 1.E-3),
NB_CYCLE = 3,
TABLE_RESU = (TABRES1, TABRES2,TABBILA),),
…
);
Le tableau ci-dessous précise pour cet exemple les résultats de simulations contenus dans les tables TABRES1et TABRES2, ainsi que l’ordre dans lequel ces tables sont remplies.
GAMMA_IMPOSE |
1.E-5 |
5.E-5 |
1.E-4 |
1.E-3 |
PRES_CONF |
||||
1.E5 |
TABRES1 |
TABRES1 |
TABRES1 |
TABRES1 |
2.05E5 |
TABRES2 |
TABRES2 |
TABRES2 |
TABRES2 |
On présente ci-dessous un extrait de la table TABRES1contenant les résultats bruts des simulations exécutées pour la première valeur de PRES_CONF(TABRES2étant construite de la même manière, pour la seconde valeur de PRES_CONF).
#
#——————————————————————————–
#
#Resultats bruts : ESSAI_CISA_C numero 1 / PRES_CONF = -1.000000E+05
#
GAMMA_IMPOSE_1 INST_1 GAMMA_1 SIG_XY_1 GAMMA_IMPOSE_2 INST_2 GAMMA_2 SIG_XY_2 …
1.00000E-05 0.00000E+00 0.00000E+00 0.00000E+00 5.00000E-05 0.00000E+00 0.00000E+00 0.00000E+00 …
4.00000E-01 -4.00000E-07 -3.79317E+01 - 4.00000E-01 -2.00000E-06 -1.89658E+02 …
8.00000E-01 -8.00000E-07 -7.58633E+01 - 8.00000E-01 -4.00000E-06 -3.79317E+02 …
1.20000E+00 -1.20000E-06 -1.13795E+02 - 1.20000E+00 -6.00000E-06 -5.50536E+02 …
1.60000E+00 -1.60000E-06 -1.51727E+02 - 1.60000E+00 -8.00000E-06 -7.19697E+02 …
2.00000E+00 -2.00000E-06 -1.89658E+02 - 2.00000E+00 -1.00000E-05 -8.88778E+02 …
2.40000E+00 -2.40000E-06 -2.27590E+02 - 2.40000E+00 -1.20000E-05 -1.05778E+03 …
2.80000E+00 -2.80000E-06 -2.65522E+02 - 2.80000E+00 -1.40000E-05 -1.22670E+03 …
3.20000E+00 -3.20000E-06 -3.03453E+02 - 3.20000E+00 -1.60000E-05 -1.39554E+03 …
3.60000E+00 -3.60000E-06 -3.41385E+02 - 3.60000E+00 -1.80000E-05 -1.56429E+03 …
4.00000E+00 -4.00000E-06 -3.79317E+02 - 4.00000E+00 -2.00000E-05 -1.73297E+03 …
4.40000E+00 -4.40000E-06 -4.15150E+02 - 4.40000E+00 -2.20000E-05 -1.90156E+03 …
4.80000E+00 -4.80000E-06 -4.49001E+02 - 4.80000E+00 -2.40000E-05 -2.07007E+03 …
5.20000E+00 -5.20000E-06 -4.82849E+02 - 5.20000E+00 -2.60000E-05 -2.23850E+03 …
5.60000E+00 -5.60000E-06 -5.16693E+02 - 5.60000E+00 -2.80000E-05 -2.40684E+03 …
6.00000E+00 -6.00000E-06 -5.50535E+02 - 6.00000E+00 -3.00000E-05 -2.57510E+03 …
6.40000E+00 -6.40000E-06 -5.84374E+02 - 6.40000E+00 -3.20000E-05 -2.74327E+03 …
6.80000E+00 -6.80000E-06 -6.18209E+02 - 6.80000E+00 -3.40000E-05 -2.91136E+03 …
7.20000E+00 -7.20000E-06 -6.52041E+02 - 7.20000E+00 -3.60000E-05 -3.07937E+03 …
7.60000E+00 -7.60000E-06 -6.85870E+02 - 7.60000E+00 -3.80000E-05 -3.24729E+03 …
8.00000E+00 -8.00000E-06 -7.19695E+02 - 8.00000E+00 -4.00000E-05 -3.41511E+03 …
8.40000E+00 -8.40000E-06 -7.53517E+02 - 8.40000E+00 -4.20000E-05 -3.58284E+03 …
8.80000E+00 -8.80000E-06 -7.87336E+02 - 8.80000E+00 -4.40000E-05 -3.75049E+03 …
9.20000E+00 -9.20000E-06 -8.21152E+02 - 9.20000E+00 -4.60000E-05 -3.91805E+03 …
9.60000E+00 -9.60000E-06 -8.54965E+02 - 9.60000E+00 -4.80000E-05 -4.08436E+03 …
1.00000E+01 -1.00000E-05 -8.88774E+02 - 1.00000E+01 -5.00000E-05 -4.24698E+03 …
1.04000E+01 -9.60000E-06 -8.50842E+02 - 1.04000E+01 -4.80000E-05 -4.05732E+03 …
1.08000E+01 -9.20000E-06 -8.12911E+02 - 1.08000E+01 -4.60000E-05 -3.86767E+03 …
1.12000E+01 -8.80000E-06 -7.74979E+02 - 1.12000E+01 -4.40000E-05 -3.67801E+03 …
1.16000E+01 -8.40000E-06 -7.37047E+02 - 1.16000E+01 -4.20000E-05 -3.48835E+03 …
1.20000E+01 -8.00000E-06 -6.99116E+02 - 1.20000E+01 -4.00000E-05 -3.31514E+03 …
1.24000E+01 -7.60000E-06 -6.61184E+02 - 1.24000E+01 -3.80000E-05 -3.14591E+03 …
1.28000E+01 -7.20000E-06 -6.23252E+02 - 1.28000E+01 -3.60000E-05 -2.97673E+03 …
1.32000E+01 -6.80000E-06 -5.85321E+02 - 1.32000E+01 -3.40000E-05 -2.80759E+03 …
1.36000E+01 -6.40000E-06 -5.47389E+02 - 1.36000E+01 -3.20000E-05 -2.63849E+03 …
1.40000E+01 -6.00000E-06 -5.09458E+02 - 1.40000E+01 -3.00000E-05 -2.46942E+03 …
1.44000E+01 -5.60000E-06 -4.71526E+02 - 1.44000E+01 -2.80000E-05 -2.30040E+03 …
1.48000E+01 -5.20000E-06 -4.33594E+02 - 1.48000E+01 -2.60000E-05 -2.13142E+03 …
1.52000E+01 -4.80000E-06 -3.95663E+02 - 1.52000E+01 -2.40000E-05 -1.96247E+03 …
1.56000E+01 -4.40000E-06 -3.57731E+02 - 1.56000E+01 -2.20000E-05 -1.79357E+03 …
1.60000E+01 -4.00000E-06 -3.19799E+02 - 1.60000E+01 -2.00000E-05 -1.62471E+03 …
1.64000E+01 -3.60000E-06 -2.81868E+02 - 1.64000E+01 -1.80000E-05 -1.45589E+03 …
1.68000E+01 -3.20000E-06 -2.43936E+02 - 1.68000E+01 -1.60000E-05 -1.28711E+03 …
1.72000E+01 -2.80000E-06 -2.06004E+02 - 1.72000E+01 -1.40000E-05 -1.11836E+03 …
1.76000E+01 -2.40000E-06 -1.68073E+02 - 1.76000E+01 -1.20000E-05 -9.49666E+02 …
1.80000E+01 -2.00000E-06 -1.30141E+02 - 1.80000E+01 -1.00000E-05 -7.81009E+02 …
1.84000E+01 -1.60000E-06 -9.23287E+01 - 1.84000E+01 -8.00000E-06 -6.12391E+02 …
… … … … … … … …
Ci dessous, on présente également le contenu de la table supplémentaire TABBILA, récapitulant les post-traitements (\(\frac{G}{{G}_{max}}\) et \(\mathit{DAMPING}\) ) réalisés à l’issue de toutes les simulations. Chaque paquet de colonnes contiguës dont les titres sont indexés par le même entier (indice de la valeur considérée dans la liste PRES_CONF) correspond aux post-traitements réalisés pour une même pression de confinement.
#
#——————————————————————————–
#
#Resultats globaux : ESSAI_CISA_C numero 1
#
PRES_CONF_1 GAMMA_IMPOSE_1 G_SUR_GMAX_1 DAMPING_1 PRES_CONF_2 GAMMA_IMPOSE_2 G_SUR_GMAX_2 DAMPING_2
-1.00000E+05 1.00000E-05 9.37244E-01 1.79484E-02 -2.00000E+05 1.00000E-05 9.65623E-01 1.30248E-02
5.00000E-05 8.95847E-01 7.93301E-03 - 5.00000E-05 9.26152E-01 6.34512E-03
1.00000E-04 7.91534E-01 5.24988E-02 - 1.00000E-04 8.82990E-01 2.26809E-02
1.00000E-03 2.81687E-01 2.16980E-01 - 1.00000E-03 3.69590E-01 1.89750E-01
Opérande GRAPHIQUE, PREFIXE_FICHIER#
◊GRAPHIQUE = |('GAMMA-SIG_XY', 'GAMMA-G_SUR_GMAX',
'GAMMA-DAMPING',
'G_SUR_GMAX-DAMPING') [DEFAUT]
| l_graphique [l_Kn]
◊PREFIXE_FICHIER= préfixe [Kn]
Idem qu’au § 3.4.5 , sauf que contrairement aux graphiques de niveau 2 qui peuvent être quelconques, les graphiques de niveau 1 sont à sélectionner dans la liste:
“GAMMA-G_SUR_GMAX”;
“GAMMA-DAMPING”;
“G_SUR_GMAX-DAMPING”
Opérande NOM_CMP#
◊NOM_CMP= l_composante [l_Kn]
Idem qu’au § 3.4.6 .
Opérande TABLE_REF#
◊ TABLE_REF = l_tabref [l_table]
Idem qu’au § 3.4.7 .
Opérandes COULEUR_NIV1, MARQUEUR_NIV1, STYLE_NIV1, COULEUR_NIV2, MARQUEUR_NIV2, STYLE_NIV2#
◊ COULEUR_NIV1= l_couleur_niv1 [l_I]
◊MARQUEUR_NIV1= l_marqueur_niv1 [l_I]
◊STYLE_NIV1= l_style_niv1 [l_I]
◊ COULEUR_NIV2= l_couleur_niv2 [l_I]
◊MARQUEUR_NIV2= l_marqueur_niv2 [l_I]
◊STYLE_NIV2= l_style_niv2 [l_I]
Ces opérandes acceptent une liste d’entiers permettant de définir respectivement la couleur, le type de marqueur et le style des courbes affichées dans les graphiques. L’information codée par chaque entier est donnée dans la documentation sur XMGRACEU2.51.01. La longueur de la liste de niveau 1 doit être égale à celle de PRES_CONF. La longueur de la liste de niveau 2 doit être égale à celle de GAMMA_IMPOSE.
MotcléESSAI_TRIA_ND_C_F#
Ce mot-clé facteur (répétable) permet d’effectuer une série de simulations d’un même essai triaxial non drainé (on suppose la saturation totale) cyclique à force imposée pour lesquelles on fait varier les paramètres de chargement (pression de confinement, amplitude de contrainte effective axiale imposée, et nombre de cycles), de post-traiter les résultats obtenus et de les écrire sous forme de graphiques (au format xmgrace) et/ou de tables.
Convention de signe des entrées et sorties#
La convention de signe des géomécaniciens s’applique aux paramètres d’entrée en contraintes ou déformations imposées, c’est-à-dire que les valeurs sontpositivesen compression.
Cette convention s’applique de même sur les variables de sortie pré-définies. En revanche, cette convention ne s’applique pas sur les variables non pré-définies demandées en sortie dans NOM_CMP(§ 3.7.9 ).
Dans cet essai, on distingue les variables de sortiepré-définies de niveau 2 qui correspondent aux variables produites par l’essaidont la liste complète est la suivante:
INST: instant
EPS_AXI: déformation axiale
EPS_LAT: déformation latérale
EPS_VOL: déformation volumique
SIG_AXI: contrainte effective axiale
SIG_LAT: contrainte effective latérale
P: contrainte effective moyenne
Q: déviateur des contraintes
PRE_EAU: pression interstitielle
RU: coefficient de pression interstitielle égale à \({r}_{u}=\frac{3}{1+2{K}_{0}}\frac{\Delta {u}_{w}}{\sigma {'}_{v,0}}\)
Etles variables de sortie pré-définies de niveau 1 qui correspondent à des courbes dont les points représentent le résultat d’un essai (par ex. la courbe \(\mathit{CRR}-{N}_{\mathit{cyc}}\) ). La liste complète des variables de niveau 1 est la suivante:
NCYCL: nombre de cycles de chargement à la liquéfaction
DSIGM: contrainte maximaleimposée normalisée par la contrainte effective moyenne initiale \(\mathit{DSIGM}=\mathit{CRR}=\frac{{Q}_{\max}}{P{'}_{0}}\)
Opérandes PRES_CONF, SIGM_IMPOSE, NB_CYCLE, NB_INST, NB_INST_MONO, TYPE_CHARGE#
♦ PRES_CONF = l_sigma_conf [l_R]
♦ SIGM_IMPOSE = l_sigma_impo [l_R]
♦ NB_CYCLE = nbcyc [I]
◊NB_INST = | 25 [DEFAUT]
| nbinst [I]
◊NB_INST_MONO= | 400 [DEFAUT]
| nbinst_mono [I]
◊ TYPE_CHARGE = | 'SINUSOIDAL' [DEFAUT]
|'TRIANGULAIRE' [Kn]
Ces opérandes permettent de définir le chargement de chacune des simulations à exécuter sous le mot-clé facteur courant, ainsi que sa discrétisation. Leur signification est résumée à la figure et détaillée ci-dessous:
PRES_CONF permet de définir la liste des pressions de confinement (strictement positives) qui seront maintenues au cours de chaque essai;
SIGM_IMPOSE permet de définir la liste des amplitudes de contrainte effective axiale du chargement cyclique imposé (avec PRES_CONF la contrainte moyenne). Une valeur strictement positive indique un premier chargement en compression, une valeur strictement négative un premier chargement en extension;
NB_CYCLEcorrespond au nombre de cycles, fixé pour toutes les simulations;
NB_INSTpermet de définir la discrétisation temporelle du chargement, et correspond au nombre de pas de chargement par quart de cycle;
NB_INST_MONO permet de définir la discrétisation temporelle des chargements monotones à déformation imposée après détection d’une instabilité, et correspond au nombre de pas de chargement par quart de cycle. Le nombre d’instants doit être suffisamment grand (>100) pour une détection suffisamment précise de l’atteinte de la consigne en force;
TYPE_CHARGE indique le type de chargement souhaité: sinusoïdal ou triangulaire;
Pour chaque pression de confinement PRES_CONF, on effectue autant de simulations qu’il y a d’éléments dans la liste SIGM_IMPOSE. Contrairement aux essais TRIA_DR_M_Det TRIA_ND_M_D(voir respectivement § 3.4 et § 3.5 ), ces listes ne sont pas en bijection et il y a au total:
\(\mathit{card}(\text{PRES\_CONF})\times \mathit{card}(\text{SIGM\_IMPOSE})\) simulations exécutées.
Figure 3.7.2-a: discrétisation et allure du chargement de type triangulaire pour le mot-clé ESSAI_TRIA_ND_C_D pour 3 cycles de chargement
Remarque:
Pour les sables lâches, le contrôle en contrainte de l’essai pose des difficultés lors de la traversée des deux lignes d’instabilité, représentées en bleue sur la Figure *. En effet, l’imposition d’une consigne de contrainte maximale supérieure à la contrainte maximale admissible sur la ligne d’instabilité conduit soit à une divergence, soit à une solution fausse (saut brutal de contrainte visible sur la Figure* ). En effet, la ligne d’instabilité représente le lieu de tous les maxima admissibles de contraintes d’un essai TRIA_ND_M_D monotone pour différentes valeurs de consolidation initiale (courbe noire).
Le problème ne se pose pas pour un sable dense, car il n’y a pas de maxima de contraintes dans ce cas-là, comme on peut le voir la courbe noire de la Figure .
Pour cet essai, il existe donc une procédure **automatiquede gestion des situations instables. Elle consiste à détecter l’instabilité et à poursuivre l’essai en déformation contrôlée. Les critères de détection sont les suivants:*
non convergence du calcul
\(\frac{\Delta Q}{\Delta P}<0.25\) et \(\frac{\Delta {{\varepsilon}}_{zz}^{+}}{\Delta {{\varepsilon}}_{zz}^{-}}>10\)
On poursuit sur le nombre de cycles restant par enchaînement d’essais TRIA_ND_M_D monotones à déformation contrôlée, à raison de deux essais par cycle ( \(\pm {{\varepsilon}}_{max}\) pour atteindre \(\pm {\sigma}_{max}\) ). La consigne de déformation maximale \({{\varepsilon}}_{max}\) imposée est de 4%, ou de 12% si la consigne précédente était insuffisante. La liste d’instants pour ces essais TRIA_ND_M_D s’échelonne de 0 à 100 secondes par pas temporel de 0.2 secondes (ou de 0.1 secondes si la consigne est \({{\varepsilon}}_{max}=12\text{\%}\) ).
L’enchaînement d’un essai TRIA_ND_M_D à déformation contrôlée à l’autre s’effectue par une reprise du calcul à partir du dernier instant où la consigne en contrainte \(\pm {\sigma}_{max}\) est atteinte.
Sur la Figure , on montre sur un exemple de sable lâche (cas-test comp012c) la solution obtenue avec ou sans la procédure de gestion de l’instabilité.
Figure 3.7.2-b: Résultat des essais TRIA_ND_M_D (noir) et TRIA_ND_C_F (rose) pour un sable dense. L’état de rupture est représentée par la ligne violette, et les lignes d’instabilité sont bleues.
Figure 3.7.2-c : Résultat des essais TRIA_ND_M_D (noir) et TRIA_ND_C_F (rose) pour un sable lâche. L’état de rupture est représentée par la ligne violette, et les lignes d’instabilité sont bleues.
Figure 3.7.2-d : Résultat d’un essai TRIA_ND_C_F pour un sable lâche avec (bleue) ou sans (rouge) la procédure de gestion de l’instabilité.
Opérande CRIT_LIQUEFACTION, VALE_CRIT, ARRET_LIQUEFACTION#
♦CRIT_LIQUEFACTION= (“RU_MAX”, “EPSI_ABSO_MAX”,
“EPSI_RELA_MAX”)
♦VALE_CRIT = l_vale [l_R]
◊ARRET_LIQUEFACTION = | “OUI” [DEFAUT]
|'NON'
Valeur maximale du critère de liquéfaction, à comparer avec : \({r}_{u}=∣\frac{u}{{P}_{0}}∣\)
Dans CRIT_LIQUEFACTION, on choisit une combinaison de critères dans la liste suivante:
RU_MAXcorrespond au coefficient de liquéfaction égal à \({r}_{u}=\frac{3}{1+2{K}_{0}}\frac{\Delta {u}_{w}}{\sigma {'}_{v,0}}\) , nécessairement compris entre \(\left]0,1\right]\) (Figure );
EPSI_ABSO_MAX correspond au seuil en déformation axiale dit «single amplitude», positif en compression et négatif en extension dans l’intervalle \(\left]-5,5\right]\phantom{\rule{2em}{0ex}}(\text{\%})\) (Figure );
EPSI_RELA_MAX correspond au critère en amplitude de déformation axiale sur une période, dit «double amplitude», nécessairement positif et inférieur à 5% (Figure );
La liste des valeurs, de même cardinal, est donnée dans VALE_CRIT. L’ordre des critères doit correspondre à l’ordre des valeurs numériques.
Si plusieurs critères sont choisis, la liquéfaction est stipulée avoir lieu quand tous les critères ont été atteints successivement. Le calcul est arrêté au dernier cycle où le dernier critère est atteint, si c’est possible, si ARRET_LIQUEFACTION = “OUI”. Si ARRET_LIQUEFACTION = “NON”, l’essai est poursuivi jusqu’au nombre de cycles renseigné dans NB_CYCLE.
Figure 3.7.3-a : Critère RU_MAX=0.9
Figure 3.7.3-b : Critères EPSI_ABSO_MAX = 5E-3 et EPSI_RELA_MAX = 1E-4
Opérande BIOT_COEF#
◊BIOT_COEF = | 1 [DEFAUT]
| biot [R]
Valeur du coefficient de Biot.
Opérande KZERO#
◊KZERO = | 1 [DEFAUT]
Idem qu’au § 3.4.3 .
Opérande UN_SUR_K#
♦UN_SUR_K = unsurk [R]
Valeur de l’inverse du module de compressibilité de l’eau.
Opérande TABLE_RESU#
◊TABLE_RESU = l_tabres [l_CO]
Cet opérande facultatif permet de donner la liste des noms des concepts produits par la macro-commande qui seront alors de type [table]. La taille de cette liste doit vérifier:
\(\mathit{card}(\text{TABLE\_RESU})=\mathit{card}(\text{PRES\_CONF})+1\)
En effet, chaque table produite regroupe les résultats bruts de toutes les simulations exécutées pour une même pression de confinement (PRES_CONF), dans laquelle chaque simulation correspond à un paquet de colonnes contiguës dont les titres sont tous indexés par le même entier (indice de la valeur considérée dans la liste SIGM_IMPOSE). Une table supplémentaire récapitulant les post-traitements réalisés à l’issue de toutes les simulations est également produite. Cette table contient pour chaque pression de confinement (PRES_CONF) le nombre de cycles au bout duquel le critère de liquéfaction du sol a été atteint NCYCL, en vis-à-vis des amplitudes de contrainte effective imposée (SIGM_IMPOSE). La liquéfaction est considérée avoir lieu quand tous les critères choisis au paragraphe § 3.7.3 sont atteints successivement.
Cette table correspond au nom de concept donné en dernière position dans la liste TABLE_RESU. Des extraits de ces tables sont présentés dans l’exemple ci-dessous.
Exemple :
TABRES1=CO(“TRES1”)
TABRES2=CO(“TRES2”)
TABRES3=CO(“TRES3”)
TABBILA=CO(“TBILA”)
CALC_ESSAI_GEOMECA(
…
ESSAI_TRIA_ND_C_F = _F(
PRES_CONF = (3.E4, 3.25E4, 3.5E4),
SIGM_IMPOSE = (1.E4, 1.1E4, 1.2E4, 1.3E4, 1.6E4),
UN_SUR_K = 1.E-12,
TABLE_RESU = (TABRES1, TABRES2, TABRES3, TABBILA)),
…
);
Le tableau ci-dessous précise pour cet exemple les résultats de simulations contenus dans les tables TABRES1,TABRES2 et TABRES3, ainsi que l’ordre dans lequel ces tables sont remplies.
SIGM_IMPOSE |
1.E4 |
1.1E4 |
1.2E4 |
1.3E4 |
1.6E4 |
PRES_CONF ¯ |
|||||
3.E4 |
TABRES1 |
TABRES1 |
TABRES1 |
TABRES1 |
TABRES1 |
3.25E4 |
TABRES2 |
TABRES2 |
TABRES2 |
TABRES2 |
TABRES2 |
3.5E4 |
TABRES3 |
TABRES3 |
TABRES3 |
TABRES3 |
TABRES3 |
On présente ci-dessous un extrait de la table TABRES2 contenant les résultats bruts des simulations exécutées pour la deuxième valeur de PRES_CONF.
#
#——————————————————————————–
#
#Resultats bruts : ESSAI_TND_C numero 1 / PRES_CONF = -3.250000E+04
#
SIGM_IMPOSE_1 INST_1 EPS_AXI_1 EPS_LAT_1 … PRE_EAU_1 SIGM_IMPOSE_2 INST_2 …
1.00000E+04 0.00000E+00 0.00000E+00 0.00000E+00 … -0.00000E+00 1.10000E+04 0.00000E+00 …
4.00000E-01 2.37788E-06 -1.18887E-06 … -1.33316E+02 - 4.00000E-01 …
8.00000E-01 5.27009E-06 -2.63491E-06 … -2.66632E+02 - 8.00000E-01 …
1.20000E+00 8.17434E-06 -4.08697E-06 … -3.99948E+02 - 1.20000E+00 …
1.60000E+00 1.10910E-05 -5.54522E-06 … -5.33263E+02 - 1.60000E+00 …
2.00000E+00 1.40203E-05 -7.00984E-06 … -6.66579E+02 - 2.00000E+00 …
2.40000E+00 1.69628E-05 -8.48101E-06 … -7.99895E+02 - 2.40000E+00 …
2.80000E+00 1.99187E-05 -9.95886E-06 … -9.33211E+02 - 2.80000E+00 …
3.20000E+00 2.30002E-05 -1.14996E-05 … -1.06631E+03 - 3.20000E+00 …
3.60000E+00 2.65689E-05 -1.32839E-05 … -1.19796E+03 - 3.60000E+00 …
4.00000E+00 3.06705E-05 -1.53346E-05 … -1.32706E+03 - 4.00000E+00 …
4.40000E+00 3.53282E-05 -1.76634E-05 … -1.45247E+03 - 4.40000E+00 …
4.80000E+00 4.05684E-05 -2.02834E-05 … -1.57305E+03 - 4.80000E+00 …
5.20000E+00 4.64193E-05 -2.32088E-05 … -1.68756E+03 - 5.20000E+00 …
5.60000E+00 5.29170E-05 -2.64576E-05 … -1.79469E+03 - 5.60000E+00 …
6.00000E+00 6.01011E-05 -3.00496E-05 … -1.89303E+03 - 6.00000E+00 …
6.40000E+00 6.80160E-05 -3.40070E-05 … -1.98113E+03 - 6.40000E+00 …
6.80000E+00 7.67173E-05 -3.83576E-05 … -2.05731E+03 - 6.80000E+00 …
7.20000E+00 8.62682E-05 -4.31330E-05 … -2.11975E+03 - 7.20000E+00 …
7.60000E+00 9.67455E-05 -4.83717E-05 … -2.16636E+03 - 7.60000E+00 …
8.00000E+00 1.08244E-04 -5.41211E-05 … -2.19469E+03 - 8.00000E+00 …
8.40000E+00 1.20878E-04 -6.04380E-05 … -2.20193E+03 - 8.40000E+00 …
8.80000E+00 1.34791E-04 -6.73947E-05 … -2.18477E+03 - 8.80000E+00 …
9.20000E+00 1.50170E-04 -7.50842E-05 … -2.13899E+03 - 9.20000E+00 …
9.60000E+00 1.67253E-04 -8.36253E-05 … -2.05939E+03 - 9.60000E+00 …
1.00000E+01 1.86360E-04 -9.31790E-05 … -1.93918E+03 - 1.00000E+01 …
1.04000E+01 1.84117E-04 -9.20575E-05 … -1.80586E+03 - 1.04000E+01 …
1.08000E+01 1.81336E-04 -9.06674E-05 … -1.67255E+03 - 1.08000E+01 …
1.12000E+01 1.78318E-04 -8.91580E-05 … -1.53923E+03 - 1.12000E+01 …
1.16000E+01 1.75296E-04 -8.76474E-05 … -1.40592E+03 - 1.16000E+01 …
1.20000E+01 1.72272E-04 -8.61353E-05 … -1.27260E+03 - 1.20000E+01 …
1.24000E+01 1.69244E-04 -8.46217E-05 … -1.13928E+03 - 1.24000E+01 …
1.28000E+01 1.66214E-04 -8.31066E-05 … -1.00597E+03 - 1.28000E+01 …
1.32000E+01 1.63181E-04 -8.15900E-05 … -8.72652E+02 - 1.32000E+01 …
1.36000E+01 1.60145E-04 -8.00719E-05 … -7.39335E+02 - 1.36000E+01 …
1.40000E+01 1.57105E-04 -7.85523E-05 … -6.06019E+02 - 1.40000E+01 …
1.44000E+01 1.54063E-04 -7.70311E-05 … -4.72703E+02 - 1.44000E+01 …
1.48000E+01 1.51017E-04 -7.55085E-05 … -3.39387E+02 - 1.48000E+01 …
1.52000E+01 1.47969E-04 -7.39842E-05 … -2.06071E+02 - 1.52000E+01 …
1.56000E+01 1.44857E-04 -7.24283E-05 … -7.26273E+01 - 1.56000E+01 …
… … … … … … …
Ci dessous, on présente également le contenu de la table supplémentaire TABBILA, récapitulant les post-traitements (nombre de cycles à la liquéfaction) réalisés à l’issue de toutes les simulations. Chaque paquet de colonnes contiguës dont les titres sont indexés par le même entier (indice de la valeur considérée dans la liste PRES_CONF) correspond aux post-traitements réalisés pour une même pression de confinement.
#
#——————————————————————————–
#
#Resultats globaux : ESSAI_TND_C numero 1
#
PRES_CONF_1 NCYCL_1 SIGM_IMPOSE_1 PRES_CONF_2 NCYCL_2 SIGM_IMPOSE_2 PRES_CONF_3 NCYCL_3 SIGM_IMPOSE_3
-3.00000E+04 1.10000E+01 1.00000E+04 -3.25000E+04 1.50000E+01 1.00000E+04 -3.50000E+04 2.10000E+01 1.00000E+04
7.00000E+00 1.10000E+04 - 1.10000E+01 1.10000E+04 - 1.50000E+01 1.10000E+04
6.00000E+00 1.20000E+04 - 8.00000E+00 1.20000E+04 - 1.10000E+01 1.20000E+04
4.00000E+00 1.30000E+04 - 6.00000E+00 1.30000E+04 - 8.00000E+00 1.30000E+04
3.00000E+00 1.60000E+04 - 3.00000E+00 1.60000E+04 - 4.00000E+00 1.60000E+04
#
Opérande GRAPHIQUE#
◊GRAPHIQUE = |('P-Q', 'SIG_AXI-PRE_EAU',
'SIG_AXI-RU', 'EPS_AXI-PRE_EAU',
'EPS_AXI-Q', 'EPS_AXI-RU',
'NCYCL-DSIGM') [DEFAUT]
| l_graphique [l_Kn]
◊PREFIXE_FICHIER= préfixe [Kn]
Idem qu’au § 3.4.5 , sauf que contrairement aux graphiques de niveau 2 qui peuvent être quelconques, les graphiques de niveau 1 sont à sélectionner dans la liste:
* 'NCYCL-DSIGM'
Opérande NOM_CMP#
◊NOM_CMP= l_composante [l_Kn]
Idem qu’au § 3.4.6 .
Opérande TABLE_REF#
◊ TABLE_REF = l_tabref [l_table]
Idem qu’au § 3.4.7 .
Opérandes COULEUR_NIV1, MARQUEUR_NIV1, STYLE_NIV1, COULEUR_NIV2, MARQUEUR_NIV2, STYLE_NIV2#
◊ COULEUR_NIV1= l_couleur_niv1 [l_I]
◊MARQUEUR_NIV1= l_marqueur_niv1 [l_I]
◊STYLE_NIV1= l_style_niv1 [l_I]
◊ COULEUR_NIV2= l_couleur_niv2 [l_I]
◊MARQUEUR_NIV2= l_marqueur_niv2 [l_I]
◊STYLE_NIV2= l_style_niv2 [l_I]
Idem qu’au § 3.4.8 .
MotcléESSAI_TRIA_DR_C_D#
Ce mot-clé facteur (répétable) permet d’effectuer une série de simulations d’un même essai triaxial drainé cyclique à déformation imposée pour lesquelles on fait varier les paramètres de chargement (pression de confinement, déformation axiale imposée, et nombre de cycles), de post-traiter les résultats obtenus et de les écrire sous forme de graphiques (au format xmgrace) et/ou de tables.
Convention de signe des entrées et sorties#
La convention de signe des géomécaniciens s’applique aux paramètres d’entrée en contraintes ou déformations imposées, c’est-à-dire que les valeurs sontpositivesen compression.
Cette convention s’applique de même sur les variables de sortie pré-définies. En revanche, cette convention ne s’applique pas sur les variables non pré-définies demandées en sortie dans NOM_CMP(§ 3.8.7 ).
Dans cet essai, on distingue les variables de sortiepré-définies de niveau 2 qui correspondent aux variables produites par l’essaidont la liste complète est la suivante:
INST: instant
EPS_AXI: déformation axiale
EPS_LAT: déformation latérale
EPS_VOL: déformation volumique
SIG_AXI: contrainte effective axiale
SIG_LAT: contrainte effective latérale
P: contrainte effective moyenne
Q: déviateur des contraintes
Etles variables de sortie pré-définies de niveau 1 qui correspondent à des courbes dont les points représentent le résultat d’un essai (par ex. la courbe \(\frac{E\text{*}}{E{\text{*}}_{\max}}-\Delta {{\varepsilon}}_{v}\) ). La liste complète des variables de niveau 1 est la suivante:
E_SUR_EMAX: \(\frac{E\text{*}}{E{\text{*}}_{\max}}\) avec \(E\text{*}=\frac{|\Delta ({\sigma}_{v}-{\sigma}_{h})|}{|\Delta {{\varepsilon}}_{v}|}\) le module apparent de l’essai;
DAMPING: amortissement hystérétique \(\frac{\Delta W}{\pi W}\)
Opérandes PRES_CONF, EPSI_MAXI, EPSI_MINI, NB_CYCLE, NB_INST#
♦ PRES_CONF = l_sigma_conf [l_R]
♦ EPSI_MAXI = l_epsi_maxi_impo [l_R]
♦ EPSI_MINI = l_epsi_mini_impo [l_R]
♦ NB_CYCLE = nbcyc [I]
◊ NB_INST = | 25 [DEFAUT]
|nbinst [I]
◊ TYPE_CHARGE = | 'SINUSOIDAL' [DEFAUT]
Ces opérandes permettent de définir le chargement de chacune des simulations à exécuter sous le mot-clé facteur courant, ainsi que sa discrétisation. Leur signification est résumée à la figure et détaillée ci-dessous:
PRES_CONF permet de définir la liste des pressions de confinement (strictement positives) qui seront maintenues au cours de chaque essai;
EPSI_MAXI et EPSI_MINI permettent de définir la liste des déformations axiales maximales et minimales, respectivement, du chargement cyclique imposé. Le cardinal de EPSI_MAXI doit être égal à celui de EPSI_MINI. La condition \(\Delta {{\varepsilon}}_{v}=\mathit{EPSI}\text{\_}\mathit{MAXI}–\mathit{EPSI}\text{\_}\mathit{MINI}>0\) doit être respectée terme à terme;
NB_CYCLEcorrespond au nombre de cycles, fixé pour toutes les simulations.
NB_INSTpermet de définir la discrétisation temporelle du chargement, et correspond au nombre de pas de chargement par quart de cycle
TYPE_CHARGE indique le type de chargement souhaité: sinusoïdal ou triangulaire;
Pour chaque pression de confinement PRES_CONF, on effectue autant de simulations qu’il y a d’éléments dans la liste GAMMA_IMPOSE. Contrairement aux essais TRIA_DR_M_Det TRIA_ND_M_D(voir respectivement § 3.4 et § 3.5 ), ces listes ne sont pas en bijection et il y a au total:
\(\mathit{card}(\text{PRES\_CONF})\times \mathit{card}(\text{EPSI\_MAXI})\) simulations exécutées.
Figure 3.8.2-a : discrétisation et allure du chargement pour le mot-clé ESSAI_TD_A, pour 3 cycles
Opérande EPSI_ELAS#
◊EPSI_ELAS = | 1.E-7 [DEFAUT]
|epsi_elas [R]
Pour chaque pression de confinement, le module de Young cyclique équivalent maximal (c’est à dire du matériau sain) est déterminé en simulant un cycle de chargement alterné contrôlé en déformation axiale imposée jusqu’à la valeur EPSI_ELAS. Cette valeur doit être telle que le matériau reste dans son domaine d’élasticité (linéaire ou non, selon la relation de comportement utilisée). EPSI_ELAS vaut 1.E-7 par défaut, et toute valeur renseignée par l’utilisateur doit lui être inférieure. Si la valeur renseignée ne permet pas de rester dans le domaine d’élasticité, le code s’arrête en erreur fatale.
Opérande KZERO#
◊KZERO = | 1 [DEFAUT]
Idem qu’au § 3.4.3 .
Opérande TABLE_RESU#
◊TABLE_RESU = l_tabres [l_CO]
Cet opérande facultatif permet de donner la liste des noms des concepts produits par la macro-commande qui seront alors de type [table]. La taille de cette liste doit vérifier:
\(\mathit{card}(\text{TABLE\_RESU})=\mathit{card}(\text{PRES\_CONF})+1\)
En effet, chaque table produite regroupe les résultats bruts de toutes les simulations exécutées pour une même pression de confinement (PRES_CONF), dans laquelle chaque simulation correspond à un paquet de colonnes contiguës dont les titres sont tous indexés par le même entier (indice de la valeur considérée dans la liste EPSI_MAXI). Une table supplémentaire récapitulant les post-traitements réalisés à l’issue de toutes les simulations est également produite. Cette table contient pour chaque pression de confinement (PRES_CONF) les valeurs de:
module apparentcyclique équivalent normalisé \(\frac{E\text{*}}{E{\text{*}}_{\max}}\) pour le dernier cycle simulé, en vis-à-vis des amplitudes de déformation imposées (EPSI_MAXI). \(E\text{*}\) et \(E{\text{*}}_{max}\) sont calculés de la manière suivante: \(E\text{*}=\frac{|\Delta ({\sigma}_{v}-{\sigma}_{h})|}{|\Delta {{\varepsilon}}_{v}|}\) avec \(\Delta \sigma ={\sigma}_{max}-{\sigma}_{min}\) et \(\Delta {{\varepsilon}}_{v}={{\varepsilon}}_{v,max}-{{\varepsilon}}_{v,min}\) ;
l’amortissement hystérétique \(\frac{\Delta W}{\pi W}\) pour le dernier cycle simulé, en vis-à-vis des amplitudes de déformation imposées (EPSI_MAXI). \(W\) Est calculé de la façon suivante: \(\Delta W={\int}_{C}\delta ({\sigma}_{v}-{\sigma}_{h}).\delta ({{\varepsilon}}_{v}-{{\varepsilon}}_{h})\) l’aire de la dernière boucle d’hystérésis \(W=\Delta ({\sigma}_{v}-{\sigma}_{h}).\Delta ({{\varepsilon}}_{v}-{{\varepsilon}}_{h})\) l’énergie élastique associée
Cette table correspond au nom de concept donné en dernière position dans la liste TABLE_RESU. Des extraits de ces tables sont présentés dans l’exemple ci-dessous.
Exemple :
TABRES1=CO(“TRES1”)
TABRES2=CO(“TRES2”)
TABBILA=CO(“TBILA”)
CALC_ESSAI_GEOMECA(
…
ESSAI_TRIA_DR_C_D = _F(
PRES_CONF = (3.E4, 5.E4),
EPSI_MAXI = (1.E-4, 5.E-4, 1.E-3, 2.E-3, 5.E-3),
EPSI_MINI= (-1.E-4, -5.E-4, -1.E-3, -2.E-3, -5.E-3),
NB_CYCLE = 3,
TABLE_RESU = (TABRES1, TABRES2, TABBILA),),
…
);
Le tableau ci-dessous précise pour cet exemple les résultats de simulations contenus dans les tables TABRES1,TABRES2, ainsi que l’ordre dans lequel ces tables sont remplies.
EPSI_IMPOSE |
1.E-4 |
5.E-4 |
1.E-3 |
2.E-3 |
5.E-3 |
PRES_CONF ¯ |
|||||
3.E4 |
TABRES1 |
TABRES1 |
TABRES1 |
TABRES1 |
TABRES1 |
5.E4 |
TABRES2 |
TABRES2 |
TABRES2 |
TABRES2 |
TABRES2 |
On présente ci-dessous un extrait de la table TABRES2 contenant les résultats bruts des simulations exécutées pour la deuxième valeur de PRES_CONF.
#
#——————————————————————————–
#
#Resultats bruts : ESSAI_TD_A numero 1 / PRES_CONF = -5.000000E+04
#
EPSI_IMPOSE_1 INST_1 EPS_AXI_1 EPS_LAT_1 … EPSI_IMPOSE_2 INST_2 …
1.00000E-04 0.00000E+00 0.00000E+00 0.00000E+00 … 5.00000E-04 0.00000E+00 …
4.00000E-01 -4.00000E-06 -1.97971E-06 … - 4.00000E-01 …
8.00000E-01 -8.00000E-06 -3.80754E-06 … - 8.00000E-01 …
1.20000E+00 -1.20000E-05 -5.47994E-06 … - 1.20000E+00 …
1.60000E+00 -1.60000E-05 -7.15274E-06 … - 1.60000E+00 …
2.00000E+00 -2.00000E-05 -8.82571E-06 … - 2.00000E+00 …
2.40000E+00 -2.40000E-05 -1.04989E-05 … - 2.40000E+00 …
2.80000E+00 -2.80000E-05 -1.21721E-05 … - 2.80000E+00 …
3.20000E+00 -3.20000E-05 -1.38456E-05 … - 3.20000E+00 …
3.60000E+00 -3.60000E-05 -1.55192E-05 … - 3.60000E+00 …
4.00000E+00 -4.00000E-05 -1.71930E-05 … - 4.00000E+00 …
4.40000E+00 -4.40000E-05 -1.88669E-05 … - 4.40000E+00 …
4.80000E+00 -4.80000E-05 -2.05409E-05 … - 4.80000E+00 …
5.20000E+00 -5.20000E-05 -2.22151E-05 … - 5.20000E+00 …
5.60000E+00 -5.60000E-05 -2.38895E-05 … - 5.60000E+00 …
6.00000E+00 -6.00000E-05 -2.55640E-05 … - 6.00000E+00 …
6.40000E+00 -6.40000E-05 -2.72385E-05 … - 6.40000E+00 …
6.80000E+00 -6.80000E-05 -2.88926E-05 … - 6.80000E+00 …
7.20000E+00 -7.20000E-05 -3.04742E-05 … - 7.20000E+00 …
7.60000E+00 -7.60000E-05 -3.19868E-05 … - 7.60000E+00 …
8.00000E+00 -8.00000E-05 -3.34349E-05 … - 8.00000E+00 …
8.40000E+00 -8.40000E-05 -3.48219E-05 … - 8.40000E+00 …
8.80000E+00 -8.80000E-05 -3.61512E-05 … - 8.80000E+00 …
9.20000E+00 -9.20000E-05 -3.74257E-05 … - 9.20000E+00 …
9.60000E+00 -9.60000E-05 -3.86481E-05 … - 9.60000E+00 …
1.00000E+01 -1.00000E-04 -3.98209E-05 … - 1.00000E+01 …
1.04000E+01 -9.60000E-05 -4.10210E-05 … - 1.04000E+01 …
1.08000E+01 -9.20000E-05 -4.23186E-05 … - 1.08000E+01 …
1.12000E+01 -8.80000E-05 -4.36541E-05 … - 1.12000E+01 …
1.16000E+01 -8.40000E-05 -4.49913E-05 … - 1.16000E+01 …
1.20000E+01 -8.00000E-05 -4.63303E-05 … - 1.20000E+01 …
1.24000E+01 -7.60000E-05 -4.76712E-05 … - 1.24000E+01 …
1.28000E+01 -7.20000E-05 -4.90139E-05 … - 1.28000E+01 …
1.32000E+01 -6.80000E-05 -5.03585E-05 … - 1.32000E+01 …
1.36000E+01 -6.40000E-05 -5.17050E-05 … - 1.36000E+01 …
1.40000E+01 -6.00000E-05 -5.30534E-05 … - 1.40000E+01 …
1.44000E+01 -5.60000E-05 -5.44039E-05 … - 1.44000E+01 …
1.48000E+01 -5.20000E-05 -5.54605E-05 … - 1.48000E+01 …
1.52000E+01 -4.80000E-05 -5.38679E-05 … - 1.52000E+01 …
1.56000E+01 -4.40000E-05 -5.22783E-05 … - 1.56000E+01 …
1.60000E+01 -4.00000E-05 -5.06917E-05 … - 1.60000E+01 …
1.64000E+01 -3.60000E-05 -4.91081E-05 … - 1.64000E+01 …
1.68000E+01 -3.20000E-05 -4.75276E-05 … - 1.68000E+01 …
1.72000E+01 -2.80000E-05 -4.59671E-05 … - 1.72000E+01 …
1.76000E+01 -2.40000E-05 -4.44567E-05 … - 1.76000E+01 …
1.80000E+01 -2.00000E-05 -4.29952E-05 … - 1.80000E+01 …
… … … … … … … …
Ci dessous, on présente également le contenu de la table supplémentaire TABBILA , récapitulant les post-traitements \(\frac{E\text{*}}{E{\text{*}}_{\max}}\) et \(\frac{\Delta W}{\pi W}\) réalisés à l’issue de toutes les simulations. Chaque paquet de colonnes contiguës dont les titres sont indexés par le même entier (indice de la valeur considérée dans la liste PRES_CONF) correspond aux post-traitements réalisés pour une même pression de confinement.
#
#——————————————————————————–
#
#Resultats globaux : ESSAI_TD_A numero 1
#
PRES_CONF_1 EPSI_IMPOSE_1 E_SUR_EMAX_1 PRES_CONF_2 EPSI_IMPOSE_2 E_SUR_EMAX_2
-3.00000E+04 1.00000E-04 4.08330E-01 -5.00000E+04 1.00000E-04 4.30922E-01
5.00000E-04 1.72000E-01 - 5.00000E-04 2.05213E-01
1.00000E-03 1.13241E-01 - 1.00000E-03 1.38724E-01
2.00000E-03 7.38666E-02 - 2.00000E-03 9.14469E-02
5.00000E-03 4.33285E-02 - 5.00000E-03 5.32730E-02
Opérande GRAPHIQUE, PREFIXE_FICHIER#
◊GRAPHIQUE = |('P-Q','EPS_AXI-Q', 'EPS_VOL-Q','EPS_AXI-EPS_VOL', 'P-EPS_VOL',
'DEPSI-E_SUR_EMAX',
'DEPSI-DAMPING') [DEFAUT]
| l_graphique [l_Kn]
◊PREFIXE_FICHIER= préfixe [Kn]
Idem qu’au § 3.4.5 , sauf que contrairement aux graphiques de niveau 2 qui peuvent être quelconques, les graphiques de niveau 1 sont à sélectionner dans la liste:
* 'DEPSI-E_SUR_EMAX'
* 'DEPSI-DAMPING'
Opérande NOM_CMP#
◊NOM_CMP= l_composante [l_Kn]
Idem qu’au § 3.4.6 .
Opérande TABLE_REF#
◊ TABLE_REF = l_tabref [l_table]
Idem qu’au § 3.4.7 .
Opérandes COULEUR_NIV1, MARQUEUR_NIV1, STYLE_NIV1, COULEUR_NIV2, MARQUEUR_NIV2, STYLE_NIV2#
◊ COULEUR_NIV1= l_couleur_niv1 [l_I]
◊MARQUEUR_NIV1= l_marqueur_niv1 [l_I]
◊STYLE_NIV1= l_style_niv1 [l_I]
◊ COULEUR_NIV2= l_couleur_niv2 [l_I]
◊MARQUEUR_NIV2= l_marqueur_niv2 [l_I]
◊STYLE_NIV2= l_style_iv2 [l_I]
Idem qu’au § 3.4.8 .
MotcléESSAI_TRIA_ND_C_D#
Ce mot-clé facteur (répétable) permet d’effectuer une série de simulations d’un même essai triaxial non drainé cyclique à déformation imposée pour lesquelles on fait varier les paramètres de chargement (pression de confinement, déformation axiale imposée, et nombre de cycles), de post-traiter les résultats obtenus et de les écrire sous forme de graphiques (au format xmgrace) et/ou de tables.
Convention de signe des entrées et sorties#
La convention de signe des géomécaniciens s’applique aux paramètres d’entrée en contraintes ou déformations imposées, c’est-à-dire que les valeurs sontpositivesen compression.
Cette convention s’applique de même sur les variables de sortie pré-définies. En revanche, cette convention ne s’applique pas sur les variables non pré-définies demandées en sortie dans NOM_CMP(§ 3.9.10 ).
Dans cet essai, on distingue les variables de sortiepré-définies de niveau 2 qui correspondent aux variables produites par l’essaidont la liste complète est la suivante:
INST: instant
EPS_AXI: déformation axiale
EPS_LAT: déformation latérale
EPS_VOL: déformation volumique
SIG_AXI: contrainte effective axiale
SIG_LAT: contrainte effective latérale
P: contrainte effective moyenne
Q: déviateur des contraintes
PRE_EAU: pression interstitielle
RU: coefficient de pression interstitielle égale à \({r}_{u}=\frac{3}{1+2{K}_{0}}\frac{\Delta {u}_{w}}{\sigma {'}_{v,0}}\)
Etles variables de sortie pré-définies de niveau 1 qui correspondent à des courbes dont les points représentent le résultat d’un essai (par ex. la courbe \(\mathit{CRR}-{N}_{\mathit{cyc}}\) ). La liste complète des variables de niveau 1 est la suivante:
NCYCL: nombre de cycles de chargement à la liquéfaction
DEPSI: amplitude de déformation axiale imposée \(\Delta {{\varepsilon}}_{v}=\text{EPSI\_MAX}-\text{EPSI\_MIN}\)
RU_MAX: maximum de RU
E_SUR_EMAX: \(\frac{E\text{*}}{E{\text{*}}_{\max}}\) avec \(E\text{*}=\frac{|\Delta ({\sigma}_{v}-{\sigma}_{h})|}{|\Delta {{\varepsilon}}_{v}|}\)
DAMPING: amortissement hystérétique \(\frac{\Delta W}{\pi W}\)
Opérandes PRES_CONF, EPSI_MAXI, EPSI_MINI, NB_CYCLE, NB_INST#
♦ PRES_CONF = l_sigma_conf [l_R]
♦ EPSI_MAXI = l_epsi_maxi_impo [l_R]
♦ EPSI_MINI = l_epsi_mini_impo [l_R]
♦ NB_CYCLE = nbcyc [I]
◊ NB_INST = | 25 [DEFAUT]
|nbinst [I]
◊ TYPE_CHARGE = | 'SINUSOIDAL' [DEFAUT]
Idem qu’au § 3.8.2 .
Opérande EPSI_ELAS#
◊EPSI_ELAS = | 1.E-7 [DEFAUT]
|epsi_elas [R]
Idem qu’au § 3.9.3 .
Opérande BIOT_COEF#
◊BIOT_COEF = | 1 [DEFAUT]
| biot [R]
Valeur du coefficient de Biot.
Opérande KZERO#
◊KZERO = | 1 [DEFAUT]
Idem qu’au § 3.4.3 .
Opérande UN_SUR_K#
♦UN_SUR_K = unsurk [R]
Valeur de l’inverse du module de compressibilité de l’eau.
Opérande RU_MAX#
◊RU_MAX = | 0.8 [DEFAUT]
Valeur du critère de liquéfaction en RU.
Opérande TABLE_RESU#
◊ TABLE_RESU = l_tabres, [l_CO]
Cet opérande facultatif permet de donner la liste des noms des concepts produits par la macro-commande qui seront alors de type [table]. La taille de cette liste doit vérifier:
\(\mathit{card}(\text{TABLE\_RESU})=\mathit{card}(\text{PRES\_CONF})+1\)
En effet, chaque table produite regroupe les résultats bruts de toutes les simulations exécutées pour une même pression de confinement (PRES_CONF), dans laquelle chaque simulation correspond à un paquet de colonnes contiguës dont les titres sont tous indexés par le même entier (indice de la valeur considérée dans la liste EPSI_MAXI). Une table supplémentaire récapitulant les post-traitements réalisés à l’issue de toutes les simulations est également produite. Cette table contient pour chaque pression de confinement (PRES_CONF) les valeurs de:
module apparentcyclique équivalent normalisé \(\frac{E\text{*}}{E{\text{*}}_{\max}}\) pour le dernier cycle simulé, en vis-à-vis des amplitudes de déformation imposées (EPSI_MAXI). \(E\text{*}\) et \(E{\text{*}}_{max}\) sont calculés de la manière suivante: \(E\text{*}=\frac{|\Delta ({\sigma}_{v}-{\sigma}_{h})|}{|\Delta {{\varepsilon}}_{v}|}\) avec \(\Delta \sigma ={\sigma}_{max}-{\sigma}_{min}\) et \(\Delta {{\varepsilon}}_{v}={{\varepsilon}}_{v,max}-{{\varepsilon}}_{v,min}\) ;
l’amortissement hystérétique \(\frac{\Delta W}{\pi W}\) pour le dernier cycle simulé, en vis-à-vis des amplitudes de déformation imposées (EPSI_MAXI). \(W\) Est calculé de la façon suivante: \(\Delta W={\int}_{C}\delta ({\sigma}_{v}-{\sigma}_{h}).\delta ({{\varepsilon}}_{v}-{{\varepsilon}}_{h})\) l’aire de la dernière boucle d’hystérésis \(W=\Delta ({\sigma}_{v}-{\sigma}_{h}).\Delta ({{\varepsilon}}_{v}-{{\varepsilon}}_{h})\) l’énergie élastique associée
nombre de cycles à la liquéfactionNCYCL
Cette table correspond au nom de concept donné en dernière position dans la liste TABLE_RESU.
Opérande GRAPHIQUE, PREFIXE_FICHIER#
◊GRAPHIQUE = |('NCYCL-DEPSI', 'DEPSI-RU_MAX',
'DEPSI-E_SUR_EMAX', 'DEPSI-DAMPING',
'P-Q', 'EPS_AXI-EPS_VOL'
'EPS_AXI-Q', 'P-EPS_VOL',
'EPS_AXI-PRE_EAU', 'EPS_AXI-RU',
'P-PRE_EAU') [DEFAUT]
| l_graphique [l_Kn]
◊PREFIXE_FICHIER= préfixe [Kn]
Idem qu’au § 3.4.5 , sauf que contrairement aux graphiques de niveau 2 qui peuvent être quelconques, les graphiques de niveau 1 sont à sélectionner dans la liste:
* 'NCYCL-DEPSI'
* 'DEPSI-RU_MAX'
* 'DEPSI-E_SUR_EMAX'
* 'DEPSI-DAMPING'
Opérande NOM_CMP#
◊NOM_CMP= l_composante [l_Kn]
Idem qu’au § 3.4.6 .
Opérande TABLE_REF#
◊ TABLE_REF = l_tabref [l_table]
Idem qu’au § 3.4.7 .
Opérandes COULEUR_NIV1, MARQUEUR_NIV1, STYLE_NIV1, COULEUR_NIV2, MARQUEUR_NIV2, STYLE_NIV2#
◊ COULEUR_NIV1= l_couleur_niv1 [l_I]
◊MARQUEUR_NIV1= l_marqueur_niv1 [l_I]
◊STYLE_NIV1= l_style_niv1 [l_I]
◊ COULEUR_NIV2= l_couleur_niv2 [l_I]
◊MARQUEUR_NIV2= l_marqueur_niv2 [l_I]
◊STYLE_NIV2= l_style_niv2 [l_I]
Idem qu’au § 3.4.8 .
MotcléESSAI_OEDO_DR_C_F#
Ce mot-clé facteur (répétable) permet d’effectuer une série de simulations d’un même essai oedométrique drainé cyclique à force imposée pour lesquelles on fait varier les paramètres de chargement (pression de consolidation isotrope initiale , amplitude de contrainte effective axiale imposée, et amplitude de contrainte effective axiale en fin de décharge), de post-traiter les résultats obtenus et de les écrire sous forme de graphiques (au format xmgrace) et/ou de tables.
Convention de signe des entrées et sorties#
La convention de signe des géomécaniciens s’applique aux paramètres d’entrée en contraintes ou déformations imposées, c’est-à-dire que les valeurs sontpositivesen compression.
Cette convention s’applique de même sur les variables de sorties pré-définies dont la liste complète pour cet essai est la suivante:
INST: instant
EPS_VOL: déformation volumique
SIG_AXI: contrainte effective axiale
SIG_LAT: contrainte effective latérale
P: contrainte effective moyenne
En revanche, cette convention ne s’applique pas sur les variables non pré-définies demandées en sortie dans NOM_CMP(§ 3.10.6 ).
Opérandes PRES_CONF, SIGM_IMPOSE, SIGM_DECH, NB_CYCLE, NB_INST, TYPE_CHARGE#
♦ PRES_CONF = l_sigma_conf [l_R]
♦ SIGM_IMPOSE = l_sigma_impo [l_R]
♦ SIGM_DECH = l_sigma_décharge [l_R]
♦ NB_CYCLE = nbcyc [I]
◊NB_INST = | 25 [DEFAUT]
| nbinst [I]
◊ TYPE_CHARGE = | 'SINUSOIDAL' [DEFAUT]
|'TRIANGULAIRE' [Kn]
Ces opérandes permettent de définir le chargement de chacune des simulations à exécuter sous le mot-clé facteur courant, ainsi que sa discrétisation. Leur signification est résumée à la figure et détaillée ci-dessous:
PRES_CONF permet de définir la liste des pressions de consolidation verticales initiales (strictement positives). La contrainte moyenne initiale sera donc égale à \(\frac{1+2{K}_{0}}{3}\times \text{PRES\_CONF}\) ;
SIGM_IMPOSE permet de définir la liste des valeurs de contrainte axiale imposées (strictement positives et supérieures à SIGM_DECH);
SIGM_DECH permet de définir la liste des valeurs (strictement positives) de la contrainte axiale de décharge, fixée pour tous les cycles de chargement, à une pression de consolidation initiale donnée. Les listes PRES_CONF et SIGM_DECH doivent avoir le même cardinal;
NB_INSTpermet de définir la discrétisation temporelle du chargement, et correspond au nombre de pas de chargement par moitié de cycle;
TYPE_CHARGE indique le type de chargement souhaité: sinusoïdal ou triangulaire;
Pour chaque pression de consolidation initiale PRES_CONF, et chaque contrainte de décharge SIGM_DECH,on effectue autant de cycles qu’il y a d’éléments dans la liste SIGM_IMPOSE. Contrairement aux essais TRIA_DR_M_D et TRIA_ND_M_D (voir respectivement § 3.4 et § 3.5 ), ces listes ne sont pas en bijection et il y a au total \(\mathit{card}(\text{PRES\_CONF})=\mathit{card}(\text{SIGM\_DECH})\) simulations exécutées, chaque simulation comportant \(\mathit{card}(\text{SIGM\_IMPOSE})\) cycles.
2.NB_INST pas de temps
SIGM_IMPOSE_4
σ’zz
SIGM_IMPOSE_1
SIGM_DECH
t
PRES_CONF
Figure 3.10.2-a: discrétisation et allure du chargement pour le mot-clé ESSAI_OEDO_C, pour 4 cycles
Opérande KZERO#
◊KZERO = | 1 [DEFAUT]
Idem qu’au § 3.4.3 .
Opérande TABLE_RESU#
◊TABLE_RESU = l_tabres [l_CO]
Cet opérande facultatif permet de donner la liste des noms des concepts produits par la macro-commande qui seront alors de type [table]. La taille de cette liste doit vérifier:
\(\mathit{card}(\text{TABLE\_RESU})=\mathit{card}(\text{PRES\_CONF})\)
En effet, chaque table produite regroupe les résultats bruts de la simulation exécutée pour une même pression de consolidation initiale (PRES_CONF) et une même valeur de contrainte en fin de décharge (SIGM_DECH), chaque cycle de cette simulation correspond à une valeur (SIGM_IMPOSE).
Exemple :
TABRES1=CO(“TRES1”)
TABRES2=CO(“TRES2”)
CALC_ESSAI_GEOMECA(
…
ESSAI_OEDO_DR_C_F = _F(
PRES_CONF = (1.E5, 2.E5),
SIGM_DECH = (2.E5, 3.E5),
SIGM_IMPOSE = (3.E5, 4.E5, 5.E5),
…
);
Les trajets de chargement sont illustrés sur la Figure .
Figure 3.10.4-a : Chargements oedométriques
Le tableau ci-dessous précise pour cet exemple les résultats de simulations contenus dans les tables TABRES1 etTABRES2, ainsi que l’ordre dans lequel ces tables sont remplies.
SIGM_IMPOSE |
3.E5 |
4.E5 |
5.E5 |
|
PRES_CONF |
SIGM_DECH |
|||
1.E5 |
2.E5 |
TABRES1 |
TABRES1 |
TABRES1 |
2.E5 |
3.E5 |
TABRES2 |
TABRES2 |
TABRES2 |
On présente ci-dessous un extrait de la table TABRES2 contenant les résultats bruts des simulations exécutées pour le deuxième valeur de PRES_CONF, et de SIGM_DECH
#
#——————————————————————————–
#
# ESSAI_OEDO_C numero 1 / PRES_CONF = -2.000000E+05 / SIGM_DECH = -3.000000E+05
#
SIGM_IMPOSE_2 INST_2 EPS_VOL_2 P_2 SIG_AXI_2 SIG_LAT_2
-3.00000E+05 0.00000E+00 0.00000E+00 -2.00000E+05 -2.00000E+05 -2.00000E+05
-3.00000E+05 4.00000E-01 -7.39723E-05 -2.01077E+05 -2.12000E+05 -1.95616E+05
-3.00000E+05 8.00000E-01 -1.75678E-04 -2.02549E+05 -2.24000E+05 -1.91824E+05
-3.00000E+05 1.20000E+00 -3.18644E-04 -2.04590E+05 -2.36000E+05 -1.88885E+05
-3.00000E+05 1.60000E+00 -4.99908E-04 -2.07143E+05 -2.48000E+05 -1.86714E+05
-3.00000E+05 2.00000E+00 -7.15464E-04 -2.10149E+05 -2.60000E+05 -1.85224E+05
-3.00000E+05 2.40000E+00 -9.61164E-04 -2.13554E+05 -2.72000E+05 -1.84332E+05
-3.00000E+05 2.80000E+00 -1.23276E-03 -2.17306E+05 -2.84000E+05 -1.83959E+05
-3.00000E+05 3.20000E+00 -1.52656E-03 -2.21361E+05 -2.96000E+05 -1.84041E+05
-3.00000E+05 3.60000E+00 -1.83869E-03 -2.25674E+05 -3.08000E+05 -1.84511E+05
-3.00000E+05 4.00000E+00 -2.16608E-03 -2.30211E+05 -3.20000E+05 -1.85317E+05
… … … … … …
… … … … … …
… … … … … …
-3.00000E+05 2.00000E+01 -6.04300E-03 -2.76465E+05 -3.00000E+05 -2.64698E+05
-4.00000E+05 2.04000E+01 -6.11132E-03 -2.80404E+05 -3.16000E+05 -2.62606E+05
-4.00000E+05 2.08000E+01 -6.19529E-03 -2.81924E+05 -3.32000E+05 -2.56886E+05
-4.00000E+05 2.12000E+01 -6.28651E-03 -2.83575E+05 -3.48000E+05 -2.51363E+05
-4.00000E+05 2.16000E+01 -6.39457E-03 -2.85529E+05 -3.64000E+05 -2.46293E+05
-4.00000E+05 2.20000E+01 -6.51836E-03 -2.87762E+05 -3.80000E+05 -2.41643E+05
-4.00000E+05 2.24000E+01 -6.65666E-03 -2.90254E+05 -3.96000E+05 -2.37381E+05
-4.00000E+05 2.28000E+01 -6.80833E-03 -2.92984E+05 -4.12000E+05 -2.33476E+05
-4.00000E+05 2.32000E+01 -6.97209E-03 -2.95933E+05 -4.28000E+05 -2.29900E+05
… … … … … …
… … … … … …
… … … … … …
-4.00000E+05 3.96000E+01 -1.01176E-02 -3.27727E+05 -3.16000E+05 -3.33591E+05
-4.00000E+05 4.00000E+01 -9.76282E-03 -3.20137E+05 -3.00000E+05 -3.30206E+05
-5.00000E+05 4.04000E+01 -9.84598E-03 -3.24672E+05 -3.20000E+05 -3.27007E+05
-5.00000E+05 4.08000E+01 -9.94584E-03 -3.26579E+05 -3.40000E+05 -3.19868E+05
-5.00000E+05 4.12000E+01 -1.00631E-02 -3.28816E+05 -3.60000E+05 -3.13225E+05
-5.00000E+05 4.16000E+01 -1.02034E-02 -3.31485E+05 -3.80000E+05 -3.07228E+05
… … … … … …
Opérande GRAPHIQUE, PREFIXE_FICHIER#
◊GRAPHIQUE = |('P-EPS_VOL', 'SIG_AXI-EPS_VOL') [DEFAUT]
| l_graphique [l_Kn]
◊PREFIXE_FICHIER = préfixe [Kn]
Idem qu’au § 3.4.5 .
Opérande NOM_CMP#
◊NOM_CMP= l_composante [l_Kn]
Idem qu’au § 3.4.6 .
Opérande TABLE_REF#
◊ TABLE_REF = l_tabref [l_table]
Idem qu’au § 3.4.7 .
Opérandes COULEUR, MARQUEUR, STYLE#
◊ COULEUR= l_couleur [l_I]
◊MARQUEUR= l_marqueur [l_I]
◊STYLE = l_style [l_I]
Idem qu’au § 3.4.8 .
MotcléESSAI_ISOT_DR_C#
Ce mot-clé facteur (répétable) permet d’effectuer une série de simulations d’un même essai de compression isotrope drainée cyclique pour lesquelles on fait varier les paramètres de chargement (pression de consolidation isotrope initiale , amplitude de contrainte effective isotrope imposée, et amplitude de contrainte isotrope de décharge), de post-traiter les résultats obtenus et de les écrire sous forme de graphiques (au format xmgrace) et/ou de tables.
Convention de signe des entrées et sorties#
Idem qu’au § 3.10.1 .
Opérandes PRES_CONF, SIGM_IMPOSE, SIGM_DECH, NB_CYCLE, NB_INST, TYPE_CHARGE#
♦ PRES_CONF = l_sigma_conf [l_R]
♦ SIGM_IMPOSE = l_sigma_impo [l_R]
♦ SIGM_DECH = l_sigma_décharge [l_R]
♦ NB_CYCLE = nbcyc [I]
◊NB_INST = | 25 [DEFAUT]
| nbinst [I]
◊ TYPE_CHARGE = | 'SINUSOIDAL' [DEFAUT]
|'TRIANGULAIRE' [Kn]
Idem qu’au § 3.10.2 avec \({K}_{0}=1\) .
Opérande TABLE_RESU#
◊TABLE_RESU = l_tabres [l_CO]
Idem qu’au § 3.10.4 .
Opérande GRAPHIQUE, PREFIXE_FICHIER#
◊GRAPHIQUE = |'P-EPS_VOL' [DEFAUT]
| l_graphique [l_Kn]
◊PREFIXE_FICHIER = préfixe [Kn]
Idem qu’au § 3.4.5 .
Opérande NOM_CMP#
◊NOM_CMP= l_composante [l_Kn]
Idem qu’au § 3.4.6 .
Opérande TABLE_REF#
◊ TABLE_REF = l_tabref [l_table]
Idem qu’au § 3.4.7 .
Opérandes COULEUR, MARQUEUR, STYLE#
◊ COULEUR= l_couleur [l_I]
◊MARQUEUR= l_marqueur [l_I]
◊STYLE = l_style [l_I]
Idem qu’au § 3.4.8 .