3. Opérandes#

3.1. Opérande MATER#

♦ MATER= / mat , [mater]

/ l_mat, [l_mater]

Ce mot-clé permet de renseigner le nom du matériau (mater) défini par DEFI_MATERIAU [U4.43.01], où sont fournis les paramètres nécessaires au comportement choisi.

Dans le cas des polycristaux, on peut avoir à donner plusieurs matériaux (cf. ssnv194c).

3.2. Mot-cléCOMPORTEMENT#

La syntaxe de ce mot clé est décrite dans le document [U4.51.11].

Certaines lois de comportement utilisent des variables d’état de type «intrinsèque», c’est-à-dire des variables d’état (de commande) qui ne sont pas définies par AFFE_VARC. Il y a deux cas:

  • Les variables d’état correspondant à la taille d’un élément. Il s’agit des comportements BETON_DOUBLE_DP, ENDO_PORO_BETON, FLUA_ENDO_PORO et RGI_BETON. Ces comportements ne sont pas utilisables dans SIMU_POINT_MAT

  • La variable d’état correspondant au gradient de vélocité utilisée dans MONOCRISTAL. Par défaut, on met cette variable à zéro

Remarque: dans le mode SUPPORT=”POINT”, on suppose que la les coordonnées du point de Gauss sont \((\mathrm{0,0,0})\) .

3.3. MotsclésINCREMENT / ARCHIVAGE/ NEWTON#

La syntaxe de ces mots-clés est décrite dans le document [U4.51.03].

Le mot-clé INCREMENT définit les intervalles de temps pris dans la méthode incrémentale.

Le mot-clé ARCHIVAGE définit les instants où sont stockés les résultats dans la table tabres. Dans le cas SUPPORT=”POINT”, ces instants ne peuvent être définis que par le mot-clé LIST_INST avec la précision relative PRECISION.

Le mots-cléNEWTON, facultatif, permet de modifier les valeurs par défaut des paramètres de convergence de la méthode de Newton.

3.4. 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.

3.4.1. Opérande RESI_GLOB_RELA/RESI_GLOB_MAXI#

◊    |RESI_GLOB_RELA =    resrel,                    [R]

L’algorithme continue les itérations globales tant que:

\(\underset{i=1,\dots ,\mathrm{nbddl}}{\max}\mid {F}_{i}^{n}\mid >\text{resrel}.\max\mid L\mid\)

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 un 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 =    resmax,                    [R]

L’algorithme continue les itérations globales tant que:

\(\underset{i=1,\dots ,\mathrm{nbddl}}{\max}\mid {F}_{i}^{n}\mid >\text{resmax}\)

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.

3.4.2. Opérande ITER_GLOB_MAXI#

◊    ITER_GLOB_MAXI =    /10                        [DEFAUT]

/maglob

Nombre d’itérations maximum effectué pour résoudre le problème global à chaque instant (10 par défaut).

3.5. MotcléRECH_LINEAIRE#

La syntaxe de ces mots-clés est décrite dans le document [U4.51.03].

Le mot-clé RECH_LINEAIRE permet, dans le cas SUPPORT=”ELEMENT”, d’activer la recherche linéaire pour aider à la convergence de l’algorithme de Newton. Cette fonctionnalité n’est pas disponible pour SUPPORT=”POINT”, car elle ne semble pas nécessaire.

3.6. MotcléMODELISATION#

Le mot-clé MODELISATION permet, dans le cas SUPPORT=”ELEMENT”, d’effectuer le calcul sur un élément 3D ou sur un élément 2D, en contraintes planes ou en déformations planes. Il n’est pas disponible dans le cas SUPPORT=”POINT”, car il suffit d’imposer une valeur nulle aux composantes correspondantes aux contraintes planes ou aux déformations planes pour obtenir le même résultat.

Ce mot-clé permet de définir la dimension du problème traité: 3D (par défaut) ou bien 2D: déformation plane ou contrainte plane. Dans le cas 2D, les composantes des tenseurs fournis sous les mots-clés SIGM_IMPOSE, EPSI_IMPOSE, SIGM_INIT, EPSI_INIT sont au nombre de 4: XX, YY, ZZ, XY.

3.7. Opérande ANGLE#

Ce mot-clé permet de spécifier un angle (en degrés) pour effectuer une rotation d’ensemble autour de \(Z\) appliquée à la fois au chargement, au maillage, et au dépouillement. Ceci permet surtout de vérifier la fiabilité de l’intégration du comportement, comme dans les tests COMP001, COMP002.

Par défaut, la rotation est identiquement nulle.

Dans le cas de matériaux possédant une orientation intrinsèque (orthotropie, comportements cristallins), il convient d’utiliser également le mot-clé MASSIF, avec une première valeur d’angle identique à celle fournie sous ANGLE.

3.8. Mot clé MASSIF#

3.8.1. Opérandes ANGL_EULER / ANGL_REP#

Ces mot-clés permettent de définir une orientation intrinsèque au matériau (orthotropie, comportements cristallins), et permettent de faire appel dans la macro-commande au mot-clé MASSIF de AFFE_CARA_ELEM[U4.42.01] .

Par défaut, l’orientation est nulle, et on ne fait pas appel à AFFE_CARA_ELEM.

3.9. Mots clés SIGM_INIT/EPSI_INIT/VARI_INIT#

Ces mots clés permettent de définir un état initial par la donnée:

  1. des composantesdes contraintes initiales (toutes les composantes ne sont pas nécessaires, par défaut on prend la valeur 0),

  2. des composantes des déformations initiales (si le mot clé EPSI_INIT est présent, il faut fournir toutes les composantes des déformations initiales: 4 en 2D, et 6 en 3D)

  3. l’ensemble des variables internes initiales pour le comportement utilisé.

Cette fonctionnalité est illustrée dans le test SSNV160E.

3.10. Mots clés SIGM_IMPOSE/EPSI_IMPOSE#

3.10.1. OpérandesSIXX, SIYY, SIZZ, SIXY, SIXZ, SIYZ#

Ces mot-clés permettent de définir des composantes du tenseur de contraintes imposées au point matériel, par l’intermédiaire de fonctions du temps. Ces fonctions peuvent être définies à l’aide de DEFI_FONCTION [U4.31.02] ou à l’aide de FORMULE [U4.31.05].

Par défaut, les composantes non affectées sont identiquement nulles.

3.10.2. OpérandesEPXX, EPYY, EPZZ, EPXY, EPXZ, EPYZ#

Ces mot-clés permettent de définir des composantes du tenseur de déformation imposées au point matériel, par l’intermédiaire de fonctions du temps. Ces fonctions peuvent être définies à l’aide de DEFI_FONCTION[U4.31.02] ou à l’aide de FORMULE [U4.31.05] .

Par défaut, les composantes non affectées sont laissées sans valeur (pas de déformation imposée).

Il convient de noter que, dans le cas du modèle de déformation PETIT_REAC, il n’est pas possible d’imposer exactement la déformation à l’aide d’EPSI_IMPOSE. En effet, à cause de la nature incrémentale de ce modèle, les déformations obtenues par ce modèle à la fin du calcul SIMU_POINT_MAT seront différentes de ce qui aura été imposé, sauf en petites déformations. Pour les grandes déformations il donc est préférable d’utiliser le modèle GDEF_LOG qui ne souffre pas de cet inconvénient.

3.11. Mots clé GRAD_IMPOSE#

3.11.1. Opérandes F11, F12, F13, F21, F22, F23,F31, F32, F33#

Ces mot-clés permettent de définir toutes les composantes du tenseur gradient de transformation imposé, en grandes déformations (DEFORMATION=”SIMO_MIEHE”) cf. test ssnd113).

3.12. Mots clésMATR_C1 / MATR_C2 / VECT_IMPO#

Ces mot-clés permettent, dans le cas SUPPORT=”POINT”, de définir directement les coefficients des matrices \(\mathit{C1}\) , \(\mathit{C2}\) et du vecteur \(g\) décrits au § 4.2 : cela permet ainsi de définir des conditions linéaires sur les inconnues (contraintes et déformations du point matériel) plus générales que les composantes imposées par les mots clés SIGM_IMPOSE/EPSI_IMPOSE. Tous les termes des matrices \(\mathit{C1}\) et \(\mathit{C2}\) non précisés sont nuls. Pour un exemple d’utilisation, voir le test WTNV134B [V7.31.134].

3.13. Opérande AFFE_VARC#

Ce mot-clé permet de spécifier une variable de commande (cf . [U4.43.03]) dont le nom est défini sous le mot clé NOM_VARC ; la fonction définissant l’évolution temporelle de cette variable de commande est fournie via le mot clé VALE_FONC. La valeur éventuelle de référence vref est donnée par VALE_REF.

Dans le cas SUPPORT=”ELEMENT”, toutes les variables de commande sont autorisées. De plus, pour M_ZIRC (resp. M_ACIER), il faut fournir les évolutions des 4 (resp. 7) phases métallurgiques en fonction du temps.

Dans le cas SUPPORT=”POINT”,seules les variables de commande “TEMP”, “SECH” et “IRRA” sont autorisées.

3.14. MotcléNB_VARI_TABLE#

Le mot-clé NB_VARI_TABLE permet, dans le cas SUPPORT=”POINT”, de limiter le nombre de variables internes écrites dans la table. En effet pour les milieux polycristallins, celui-ci peut atteindre plusieurs milliers. On limite alors le nombre de colonnes de la table à nvar. Par contre les calculs sont bien sûr effectués avec la totalité des variables internes : celles-ci ne sont tronquées que dans la table en résultat.

3.15. MotcléFORMAT_TABLE#

◊ FORMAT_TABLE = /‘CMP_COLONNE’ [DEFAUT]

/‘CMP_LIGNE’

Le mot-clé FORMAT_TABLE permet, dans le cas SUPPORT=”POINT”, de définir le mode de stockage des grandeurs dans la table résultat (le test SSNV194C illustre ces deux formats) . Dans le cas où le nombre de variables internes dépasse le maximum de colonnes autorisé pour une table (9999, cf. D4,02,05), le format bascule automatiquement en : FORMAT_TABLE = /‘CMP_LIGNE’.

FORMAT_TABLE = /‘CMP_LIGNE’ :

… V845 V846 V847 V848 NB_ITER

… 1.16186E-17 1.32359E-17 1.11751E-17 1.00000E+00 1.00000E+00

… 1.29473E-16 1.47341E-16 1.24474E-16 1.00000E+00 1.00000E+00

… 8.90739E-16 1.00875E-15 8.55093E-16 1.00000E+00 1.00000E+00

… 4.40109E-15 4.92817E-15 4.21938E-15 1.00000E+00 1.00000E+00

… 1.70332E-14 1.87022E-14 1.63484E-14 1.00000E+00 1.00000E+00

… 5.44940E-14 5.80870E-14 5.25904E-14 1.00000E+00 1.00000E+00

FORMAT_TABLE = /‘CMP_COLONNE’ :

… INST GRANDEUR CMP VALEUR

4.97867E-03 VARI V845 -1.43828E+01

4.97867E-03 VARI V846 -2.63548E+01

4.97867E-03 VARI V847 2.80907E+01

4.97867E-03 VARI V848 1.00000E+00

5.00000E-03 EPSI EPXX -2.20535E-03

5.00000E-03 EPSI EPYY -1.96506E-03

5.00000E-03 EPSI EPZZ 5.00000E-03

5.00000E-03 EPSI EPXY -1.98892E-04

5.00000E-03 EPSI EPXZ -2.11427E-04

5.00000E-03 EPSI EPYZ -3.00870E-04

5.00000E-03 V7.32.119SIGM SIXX 1.67146E-04

5.00000E-03 SIGM SIYY 2.78713E-05

5.00000E-03 SIGM SIZZ 3.01140E+02

5.00000E-03 SIGM SIXY 9.15194E-05

5.00000E-03 SIGM SIXZ 1.62000E-04

5.00000E-03 SIGM SIYZ 6.86376E-05

5.00000E-03 SIEQ VMIS 3.01140E+02

5.00000E-03 SIEQ TRACE 3.01140E+02

5.00000E-03 VARI V1 -1.58316E-03

5.00000E-03 VARI V2 -1.34287E-03

5.00000E-03 VARI V3 2.92604E-03

5.00000E-03 VARI V4 -2.81276E-04

3.16. MotcléOPER_TANGENT#

◊ OPER_TANGENT = /‘NON’ [DEFAUT]

/‘OUI’

Le mot-clé OPER_TANGENT permet, dans le cas SUPPORT=”POINT”, d’ajouter à la table résultat les 36 valeurs de l’opérateur tangent issu du comportement.

3.17. OpérandeINFO#

◊    INFO =    inf

Permet d’effectuer dans le fichier message diverses impressions intermédiaires.