u4.82.07 Opérateur POST_KCP correction β#

Syntaxe#

Détail de la syntaxe
table = POST_KCP(
    ◆ MODELE = modele,
    ◆ MATER_REV = mater,
    ◆ EPAIS_REV = float,
    ◆ MATER_MDB = mater,
    ◆ EPAIS_MDB = float,
    ◆ UNITE_LONGUEUR = / "M" (par défaut),
                       / "MM",
    ◆ FISSURE = _F(
         ◇ FORM_FISS = / "SEMI_ELLIPTIQUE",
                       / "BANDE",
         # Si: equal_to("FORM_FISS", 'SEMI_ELLIPTIQUE')
             ◆ LONGUEUR = float,
             ◆ CORRECTION = / "BETA_3D" (par défaut),
                           / "BETA_2D",
         ◆ PROFONDEUR = float,
         ◆ ORIENTATION = / "CIRC",
                         / "LONGI",
      ),
    ◆ K1D = _F(
         ◆ TABL_MECA= table,
         ◆ TABL_THER = table,
      ),
)


◆ : obligatoire
◇ : optionnel
⟐ : présent par défaut
& : ensemble
/ : un seul parmi
| : plusieurs choix possibles

Opérandes#

L’opérateur POST_KCP peut-être utilisé dans le cadre de modélisations axisymétriques ou tridimensionnelles.

Sur la figure ci-dessous nous présentons une schématisation d’une tranche de cuve avec la présence d’un défaut semi-elliptique.

../../../../_images/Cadre310.gif

Défaut semi-elliptique situé sous le revêtement

Opérande MODELE#

Concept de type modele.

Opérande MATER_REV#

Nom du concept de type materiau définissant le matériau constituant le revêtement. Nécessaire pour récupérer :

  • Le module d’Young, pour le calcul des coefficients d’influence (défaut semi-elliptique),

  • Les limites élastiques, pour la correction plastique des facteurs d’intensité de contraintes.

Opérande EPAIS_REV#

Épaisseur du revêtement. Nécessaire pour :

  • Le calcul des facteurs d’intensité de contrainte dans le cas de la méthode des coefficients d’influence (défaut semi-elliptique)

  • la correction plastique.

Opérande MATER_MDB#

Nom du concept de type materiau définissant le matériau constituant le métal de base. Nécessaire pour récupérer le module d’Young, pour le calcul des coefficients d’influence (défaut semi-elliptique).

Opérande EPAIS_MDB#

Épaisseur du métal de base, nécessaire pour la correction des facteurs d’intensité de contraintes par la méthode des coefficients d’influence (défaut semi-elliptique).

Opérande UNITE_LONGUEUR#

Permet d’utiliser le millimètre comme unité de longueur.

Si l’utilisateur spécifie UNITE_LONGUEUR = “MM” dans la commande POST_KCP, cela signifie que l’utilisateur souhaite manipuler des unités de longueur en millimètres, et la commande POST_KCP effectue la transformation des unités.

Mot clé FISSURE#

Mot-clé facteur pour la caractérisation géométrique du défaut. Il ne peut être utilisé qu’une seule fois.

Opérande FORM_FISS#

Forme du défaut:

“BANDE” pour un défaut bande a/c = 0,

“SEMI_ELLIPTIQUE” pour un défaut semi-elliptique,

Opérande CORRECTION#

CORRECTION:

“beta_2d” pour un défaut bande,

“beta_3d” pour un défaut semi-elliptique,

Opérande LONGUEUR#

Deuxième dimension du défaut (axiale ou orthoradiale suivant l’orientation du défaut) en SI.

Opérande PROFONDEUR#

Dimension radiale du défaut en SI.

Opérande ORIENTATION#

Caractérisation de l’orientation du défaut:

“CIRC” pour un défaut circonférentiel,

“LONGI” pour un défaut longitudinal.

Mot clé K1D#

Mot-clé facteur pour la caractérisation du transitoire thermomécanique. La répétition de ce mot-clé est possible.

Opérandes TABL_MECA#

La table TABL_MECA est obligatoire, elle fournit le transitoire des contraintes subies par la cuve au cours de l’histoire du chargement, respectivement côté revêtement (de la pointe A du défaut : interface entre le revêtement et le métal de base) et côté métal de base (à la pointe B du défaut).

Les paramètres communs nécessaires à ces tables sont :

  • “INST”[R],

  • “ABSC_CURV”[R]

  • “COOR_X”[R],

  • “COOR_Y”[R]

  • “COOR_Z”[R] (en 3D)

  • “SIXX”[R] (en 3D avec “LONGI”),

  • “SIYY”[R] (en 3D avec “LONGI” ou en 2D avec “CIRC”),

  • “SIXY”[R] (en 3D avec “LONGI”),

  • “SIZZ”[R] (en 2D avec “LONGI” ou en 3D avec “CIRC”).

Le lecteur est invité à consulter la documentation des commandes POST_RELEVE_T [U4.81.21] et MACR_LIGN_COUPE [U4.81,13] pour davantage d’informations sur le sens des paramètres.

La macro-commande MACR_LIGN_COUPE doit être utilisée avec les mots clés TYPE=”SEGMENT” et NB_POINTS=5. La méthode des coefficients d’influence nécessite l’évaluation d’un polynôme de degré 5 représentant l’évolution de la contrainte normale au segment représentant la fissure [R7.02.10]. La macro-commande MACR_LIGN_COUPE permet à l’utilisateur de définir les lieux où la contrainte est relevée et ceci indépendamment du maillage utilisé lors de l’étude thermo-mécanique.

Opérande TABL_THER#

Table fournissant le transitoire thermique dans la cuve au cours de l’histoire du chargement de la pointe A à la pointe B.

Les paramètres nécessaires de cette table sont:

  • “INST”[R],

  • “ABSC_CURV”[R],

  • “TEMP”[R].

Le lecteur est invité à consulter la documentation des commandes POST_RELEVE_T [U4.81.21] et MACR_LIGN_COUPE [U4.81.13] pour davantage d’informations sur le sens des paramètres.

Table produite#

Les paramètres de la table produite sont décrits dans le tableau suivant :

Forme du défaut : FORM_FISS = “SEMI_ELLIPTIQUE”

PARAMETRE

TYPE

DESCRIPTION

INST

R

instant

TEMP_A

R

température à la pointe de fissure côté revêtement (pointe A)

TEMP_B

R

température à la pointe de fissure côté métal de base (pointe B)

K1_A

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté revêtement (pointe A)

K1_B

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté métal de base (pointe B)

K1_C

R

facteur d’intensité des contraintes à la pointe de fissure côté revêtement (pointe C)

K2_A

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté revêtement (pointe A)

K2_B

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté métal de base (pointe B)

K2_C

R

facteur d’intensité des contraintes à la pointe de fissure côté revêtement (pointe C)

K3_A

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté revêtement (pointe A)

K3_B

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté métal de base (pointe B)

K3_C

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté revêtement (pointe C)

Keq_A

R

facteur d’intensité des contraintes élastiques équivalent à la pointe de fissure côté revêtement (pointe A)

Keq_B

R

facteur d’intensité des contraintes élastiques équivalent à la pointe de fissure côté métal de base (pointe B)

Keq_C

R

facteur d’intensité des contraintes élastiques équivalent à la pointe de fissure côté revêtement (pointe C)

KCP_A

R

facteur d’intensité des contraintes avec correction plastique à la pointe de fissure côté revêtement (pointe A)

KCP_B

R

facteur d’intensité des contraintes avec correction plastique à la pointe de fissure côté métal de base (pointe B)

KCP_C

R

facteur d’intensité des contraintes avec correction plastique à la pointe de fissure côté revêtement (pointe C)

Forme du défaut : FORM_FISS = “BANDE”

PARAMETRE

TYPE

DESCRIPTION

INST

R

instant

TEMP_A

R

température à la pointe de fissure côté revêtement (pointe A)

TEMP_B

R

température à la pointe de fissure côté métal de base (pointe B)

K1_A

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté revêtement (pointe A)

K1_B

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté métal de base (pointe B)

K2_A

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté revêtement (pointe A)

K2_B

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté métal de base (pointe B)

K3_A

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté revêtement (pointe A)

K3_B

R

facteur d’intensité des contraintes élastiques à la pointe de fissure côté métal de base (pointe B)

Keq_A

R

facteur d’intensité des contraintes élastiques équivalent à la pointe de fissure côté revêtement (pointe A)

Keq_B

R

facteur d’intensité des contraintes élastiques équivalent à la pointe de fissure côté métal de base (pointe B)

KCP_A

R

facteur d’intensité des contraintes avec correction plastique à la pointe de fissure côté revêtement (pointe A)

KCP_B

R

facteur d’intensité des contraintes avec correction plastique à la pointe de fissure côté métal de base (pointe B)

Exemples#

Des exemples d’utilisation de la commande POST_KCP sont fournis dans le cas test EPICU03.

Avant l’utilisation de la commande POST_KCP, il est nécessaire de relever les contraintes et les températures le long du défaut.

Les tables TABL_MECA et TABL_THER sont obtenues directement en utilisant la macro-commande MACR_LIGN_COUPE sur le segment support du défaut allant de la pointe A du défaut côté revêtement jusqu’à la pointe B côté métal de base.

S1G = MACR_LIGN_COUPE(RESULTAT=SIG,
                      VIS_A_VIS=_F(GROUP_MA_1='M'),
                      NOM_CHAM='SIGM ELNO',
                      LIGN_COUPE=_F(TYPE='SEGMENT',
                                    NB_POINTS=5,
                                    COOR_ORIG=(Pointe_A, 0., 0.),
                                    COOR_EXTR=(Pointe_B, 0., 0.)))

TEMP_G = MACR_LIGN_COUPE(RESULTAT=TEMP,
                         NOM_CHAM='TEMP',
                         LIST_INST=INST_TH,
                         LIGN_COUPE=_F(TYPE='SEGMENT',
                                       NB_POINTS=5,
                                       COOR_ORIG=(Pointe_A, 0., 0.),
                                       COOR_EXTR=(Pointe_B, 0., 0.)))

Après avoir effectué ces différents relevés, le calcul du facteur d’intensité de contrainte peut se faire effectivement à l’aide de la commande POST_KCP.

KBETA = POST_KCP(MODELE=MO,
                    MATER_REV=MAME_RE2,
                    MATER_MDB=MAME_MET,
                    EPAIS_REV=EPREV,
                    EPAIS_MDB=EPMDB,
                    UNITE_LONGUEUR='M',
                    FISSURE=_F(FORM_FISS='SEMI_ELLIPTIQUE',
                               PROFONDEUR=0.005,
                               LONGUEUR=0.025,
                               ORIENTATION='LONGI'),
                    K1D=_F(TABL_MECA=S1G,
                           TABL_THER=TEMP_G),)