Télécharger g_theta.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ G_THETA  NOTICE  JB251061  21/12/13    21:15:10     11219          
   2 :                                              DATE     21/12/13
   3 : 
   4 :   Procedure G_THETA                        Voir aussi : G_CALCUL
   5 :     -----------------                                     CH_THETA  
   6 :                                                           CH_THETX  
   7 :                                                           G_CAS  
   8 :                                                           G_AUX  
   9 :     G_THETA SUPTAB ;
  10 : 
  11 :                        SUPTAB.'BLOCAGES_MECANIQUES'    'LEVRE_INFERIEURE'
  12 :                               'CALCUL_CRITERE'         'METH_AUX'
  13 :                               'CARACTERISTIQUES'       'MODELE'
  14 :                               'CHAMP_THETA'            'MODELES_COMPOSITES'
  15 :                               'CHARGEMENTS_MECANIQUES' 'NOEUDS_AVANCES'
  16 :                               'CHPOINT_TRANSFORMATION' 'OBJECTIF'
  17 :                               'CHPO_RESULTATS'         'OPERATEUR'
  18 :                               'COUCHE'                 'PHI'
  19 :                               'CRIT_DECHA_GLOBAL1'     'POINT_CENTRE'
  20 :                               'CRIT_DECHA_GLOBAL2'     'POINT_1'
  21 :                               'CRIT_DECHA_LOCAL1'      'POINT_2'
  22 :                               'CRIT_DECHA_LOCAL2'      'POINT_3'
  23 :                               'DEFORMATIONS_IMPOSEES'  'PRESSION'
  24 :                               'ELEMENT_MULTICOUCHE'    'PSI'
  25 :                               'EPAISSEUR_RESULTATS'    'RESULTATS'
  26 :                               'EVOLUTION_RESULTATS'    'ROTATION_RIGIDIFIANTE'
  27 :                               'FISSURE_2'              'SOLUTION_PASAPAS'
  28 :                               'FRONT_FISSURE'          'SOLUTION_RESO'
  29 :                               'FRONT_FISSURE_2'        'TEMPERATURES'
  30 :                               'LEVRE_SUPERIEURE'
  31 : 
  32 : 
  33 :     Objet :
  34 :     -------
  35 : 
  36 :     Cette procedure a deux objectifs principaux :
  37 : 
  38 :     1 ) calculer les integrales suivantes de la mecanique de la rupture :
  39 :         1.1) l'integrale J (ou G) d'un materiau isotrope, caracteristique
  40 :              en elasto-plastique. Les discontinuites de proprietes ne sont
  41 :              pas encore acceptables dans le cas des elements 3D massifs.
  42 :         1.2) l'integrale J dynamique d'un materiau isotrope,
  43 :              caracteristique en elasto-dynamique. Les discontinuites de
  44 :              proprietes ne sont pas encore acceptables dans le cas des
  45 :              elements 3D massifs.
  46 :         1.3) l'integrale C* d'un materiau isotrope, caracteristique
  47 :              dans le cas de fluage secondaire stationnaire. Le chargement
  48 :              doit etre mecanique et les dicsontinuites de proprietes ne sont
  49 :              pas encore acceptables, ni en 2D ni en 3D.
  50 :         1.4) l'integrale C*(h) d'un materiau isotrope, caracteristique
  51 :              dans le cas de fluage primaire ou tertiaire sous un chargement
  52 :              radial. Le chargement doit etre mecanique et les dicsontinuites
  53 :              de proprietes ne sont pas encore acceptables, ni en 2D ni en 3D.
  54 :         1.5) l'integrale de derivation dJ/da (a : longueur de la fissure)
  55 :              d'un materiau homogene et isotrope, utile pour etudier la
  56 :              stabilite de propagation d'une fissure ou des fissures
  57 :              interagissantes. Ne sont pas encore acceptables les discontinuites
  58 :                          de proprietes en 2D ou en 3D et les elements de coque (mince
  59 :              ou epaisse).
  60 : 
  61 :     2 ) decoupler les modes mixtes d'un solide homogene constitue d'un materiau
  62 :         elastique lineaire isotrope, c'est a dire la separation des facteurs
  63 :         d'intensite des contraintes K1, K2 (et K3 en 3D).
  64 :         Les discontinuites de proprietes en 3D et les elements de coque ne sont
  65 :         pas encore acceptables.
  66 :         Les proprietes materielles doivent etre des constantes.
  67 : 
  68 : 
  69 :     ENTREE :
  70 :     --------
  71 : 
  72 :     En entree, SUPTAB (objet de type TABLE) sert a definir les options
  73 :     et les parametres du calcul. Ses indices sont des objets de type
  74 :     MOT (a ecrire en toutes lettres) dont voici la liste :
  75 : 
  76 :     Arguments obligatoires dans tous les cas
  77 :     ----------------------------------------
  78 : 
  79 :     SUPTAB.'OBJECTIF'
  80 :     = MOT pour preciser le but du calcul, valant :
  81 :       1) 'J'      pour calculer l'integrale J (ou G), caracteristique
  82 :                   en elasto-plastique.
  83 :       2) 'J_DYNA' pour calculer l'integrale J (ou G), caracteristique
  84 :                   en elasto-dynamique.
  85 :       3) 'C*'     pour calculer l'integrale C*, caracteristique
  86 :                   en fluage secondaire stationnaire.
  87 :       4) 'C*H'    pour calculer l'integrale C*(h), caracteristique
  88 :                   en fluage primaire ou tertiaire.
  89 :       5) 'DJ/DA'  pour calculer l'integrale de la derivation dJ/da,
  90 :                   caracteristique pour analyser la stabilite de
  91 :                   propagation d'une fissure ou des fissures
  92 :                   interagissantes.
  93 :       6) 'DECOUPLAGE' pour decouper les modes mixtes, c'est a dire la
  94 :                   separation des facteurs K1, K2 (et K3 et 3D).
  95 : 
  96 :     SUPTAB.'COUCHE'
  97 :     = ENTIER representant le nombre de couches d'elements autour du
  98 :              front de la fissure qui se deplacent pour simuler la
  99 :              propagation de la fissure. Il vaut 0 si seul la pointe de
 100 :              la fissure se deplace, 1 si c'est la premiere couche
 101 :              d'elements entourant la fissure qui se deplace, 2 si c'est
 102 :              l'ensemble des premiere et deuxieme couches d'elements qui
 103 :              se deplace, etc.
 104 :              Il convient veiller a ce que l'ensemble des elements a deplacer
 105 :              n'atteint pas le bord de la structure fissuree.
 106 :              Si COUCHE et CHAMP_THETA sont tous deux donnés, CHAMP_THETA
 107 :              est écrasé (cf.8.).
 108 : 
 109 :     SUPTAB.'FRONT_FISSURE'
 110 :     = type POINT en 2D massif ou 3D coque
 111 :       representant la pointe de la fissure;
 112 :     = type MAILLAGE en 3D massif (elements SEG2 ou SEG3)
 113 :       representant le front de la fissure.
 114 : 
 115 :     Arguments obligatoires avec des elements standards
 116 :     --------------------------------------------------
 117 : 
 118 :     SUPTAB.'LEVRE_SUPERIEURE'
 119 :     = Selon la convention habituelle de definition, cet objet (type
 120 :       MAILLAGE) represente la levre superieure de la fissure.
 121 : 
 122 :     SUPTAB.'LEVRE_INFERIEURE'
 123 :     = Selon la convention habituelle de definition, cet objet (type
 124 :       MAILLAGE) represente la levre inferieure de la fissure.
 125 : 
 126 :     Si une seule levre est modelisee, un des deux mots (LEVRE_SUPERIEURE
 127 :     ou LEVRE_INFERIEURE) sera suffisant pour decrire la fissure entiere.
 128 :     Dans le cas de l'objectif DECOUPLAGE, les deux levres sont obligatoires.
 129 : 
 130 :     Arguments obligatoires avec des elements enrichis (XFEM)
 131 :     --------------------------------------------------------
 132 : 
 133 :     SUPTAB.'PSI' =   1ere level set (CHPOINT) decrivant la fissure dans
 134 :                      le cas ou l'on utilise des elements XFEM .
 135 : 
 136 :     SUPTAB.'PHI' =   2eme level set.
 137 : 
 138 : 
 139 :     Arguments obligatoires avec une solution issue de la procedure PASAPAS
 140 :     ----------------------------------------------------------------------
 141 : 
 142 :     SUPTAB.'SOLUTION_PASAPAS' = TABLE sortant de la procedure PASAPAS.
 143 : 
 144 : 
 145 :     Arguments obligatoires avec une solution issue de l'operateur RESO
 146 :     ------------------------------------------------------------------
 147 : 
 148 :     SUPTAB.'SOLUTION_RESO'
 149 :     = CHPOINT de deplacement issus de RESO.
 150 : 
 151 :     SUPTAB.'CARACTERISTIQUES'
 152 :     = Champ de caractristiques materielles et eventuellement
 153 :       geometriques si necessaire.
 154 : 
 155 :     SUPTAB.'MODELE'
 156 :     = Objet modele (type MMODEL) englobant toute la structure.
 157 : 
 158 :     Arguments optionnels avec une solution issue de l'operateur RESO
 159 :     ----------------------------------------------------------------
 160 : 
 161 :     SUPTAB.'TEMPERATURES'
 162 :     = CHPOINT ou MCHAML de temperature creant une contrainte thermique
 163 :       non nulle si elle existe.
 164 : 
 165 :     SUPTAB.'CHARGEMENTS_MECANIQUES'
 166 :     = CHPOINT representant l'ensemble des forces exterieures
 167 :       (surfaciques, volumiques ou ponctuelles ....) appliquees
 168 :       sur le systeme si elles existent, SAUF la pression sur les
 169 :       levres de la fissure qui doit obligatoirement etre donnee
 170 :       dans l'indice PRESSION et le modele correspondant dans
 171 :       l'indice MODELE.
 172 : 
 173 :     SUPTAB.'BLOCAGES_MECANIQUES'
 174 :     = RIGIDITE representant les blocages mecaniques du probleme, a fournir
 175 :       uniquement pour le calcul de l'integrale de derivation dJ/da.
 176 : 
 177 :     SUPTAB.'PRESSION'
 178 :     = MCHAML de pression, obligatoire si l'on fournit un modele de
 179 :       chargement pression a l'indice MODELE.
 180 : 
 181 :     SUPTAB.'DEFORMATIONS_IMPOSEES'
 182 :     = MCHAML de deformations imposees s'il y en a
 183 : 
 184 :     Autres arguments optionnels
 185 :     ---------------------------
 186 : 
 187 :     1 : Discontinuites de proprietes (2D massif ou 3D coque seulement)
 188 : 
 189 :     SUPTAB.'MODELES_COMPOSITES'
 190 :     = TABLE indicee par des entiers (1 2... M, M = nombre de modeles)
 191 :       pour donner les modeles permettant d'identifier les discontinuites
 192 :       de proprietes materielles ou geometriques. Les proprietes sont donc
 193 :       continues sur chaque modele, et les lieux de discontinuites sont
 194 :       les interfaces entre les modeles.
 195 : 
 196 :     2 : Pour un front de fissure tridimensionnel massif
 197 : 
 198 :     SUPTAB.'NOEUDS_AVANCES'
 199 :     = MAILLAGE de type POI1 pour donner les points du front pour
 200 :       lesquels le calcul sera effectue. Si cet argument est obsent, le
 201 :       calcul sera fait pour tous les noeuds sur le front de la fissure.
 202 : 
 203 :     3 : Calcul des termes croises de la matrice de derivation dJi/daj
 204 :         (i n'est pas egal a j) dans le cas des fisures interagissantes.
 205 : 
 206 :     SUPTAB.'FISSURE_2'
 207 :     = Objet de type MAILLAGE representant une autre fissure (levres
 208 :       superieure + inferieure si toutes les deux levres sont presentes).
 209 : 
 210 :     SUPTAB.'FRONT_FISSURE_2'
 211 :     = POINT ou MAILLAGE reprsentant le front de la fissure 2 telle que
 212 :       decrite ci-dessus.
 213 : 
 214 :     4 : Cas d'une fissure circulaire dans une geometrie plane
 215 : 
 216 :     SUPTAB.'POINT_CENTRE' = centre de la fissure circulaire
 217 : 
 218 :     5 : Cas ou l'extension de la fissure correspond a une simple
 219 :         translation dans un tuyauterie droite (3D). Dans ce cas
 220 :         on effectue dans la procedure CH_THETA une transformation
 221 :         de tuyau en plaque en passant au systeme de coordonnees
 222 :         cylindriques. Il est alors necessaire de fournir :
 223 : 
 224 :     SUPTAB.'POINT_1' = centre du systeme de coordonnees
 225 : 
 226 :     SUPTAB.'POINT_2' = POINT tel que l'axe defini par POINT_1
 227 :                        vers POINT_2 soit l'axe Z poisitif
 228 : 
 229 :     SUPTAB.'POINT_3' = POINT tel que le plan defini par les 3 points
 230 :                        POINT_1 POINT_2 POINT_3 donne l'angle theta nul
 231 : 
 232 :     6 : Cas ou l'extension de la fissure ne correspond
 233 :         pas a une simple translation (3D)
 234 : 
 235 :         6.1 : Fissure dans un tuyauterie droite (3D, Rotation)
 236 : 
 237 :         SUPTAB.'POINT_1' = Objet de type POINT sur l'axe du tuyau
 238 : 
 239 :         SUPTAB.'POINT_2' = Objet de type POINT qui, avec le point
 240 :                            SUPTAB.'POINT_1', defini l'axe du tuyau
 241 : 
 242 :         6.2 : Fissure dans un coude (3D, rotation + transformation)
 243 :               Outre les deux points SUPTAB.'POINT_1' et SUPTAB.'POINT_2'
 244 :               definis en haut on donne encore :
 245 : 
 246 :         SUPTAB.'CHPOINT_TRANSFORMATION'
 247 :         = Objet de type CHPOINT utilise pour transformer un coude en un
 248 :           tuyauterie droite.
 249 :         SUPTAB.'OPERATEUR'
 250 :         = Objet de type MOT valant 'PLUS' ou 'MOIN' pour indiquer
 251 :           l'operateur PLUS ou MOIN a utiliser si l'on veut transformer
 252 :           le coude en un tuyauterie droite.
 253 : 
 254 :     7 : Rotation rigidifiante imposee dans le calcul par PASAPAS
 255 : 
 256 :     SUPTAB.'ROTATION_RIGIDIFIANTE'
 257 :     = table indicee par entiers 0,1,2...donnant les champs de
 258 :       deplacements due a une rotation rigidifiante de la piece autour
 259 :       d'un point. Cette rotation rigidifiante est imposee dans le calcul
 260 :       par PASAPAS en tant d'un calcul en grand deplacement.
 261 : 
 262 :     8 : Cas ou on souhaite donner soi-meme un champ de type Theta
 263 : 
 264 :       SUPTAB.'CHAMP_THETA'
 265 :       = Objet de type CHPOINT caracterisant l'avancee virtuelle de la
 266 :         fissure. Attention : dans le cas d'une demi-eprouvette
 267 :         (condition de symetrie dans le plan de fissure), la norme du
 268 :         champ theta doit logiquement varier de 2. à 0.
 269 :        (au lieu de 1. à 0. dans le cas d'une éprouvette complète).
 270 : 
 271 :     9 : Cas ou on souhaite calculer une integrale dans l epaisseur
 272 :         d une structure en coque (rapport DMT/96-317)
 273 : 
 274 :         On utilise pour cela la technique de multicouche, qui
 275 :         consiste, avant d'appeler la proceduer G_THETA, a :
 276 :         1) Etablir un modele multicouches (cf MODE CONS) sur un ou
 277 :            des element(s) proche(s) de la fissure sachant qu'il faut
 278 :            au moins une couche en peau inferieure, une couche en
 279 :            peau superieure, une couche en ligne moyenne {ces couches
 280 :            doivent avoir une epaisseur inferieure a 1e-4*(epaisseur
 281 :            totale de la coque) et donc 2 couches intermediaires.
 282 :         2) Penser a donner un excentrement et un nom constituant
 283 :            different a ces couches.
 284 :         3) Assembler le modele multicouches avec le modele du reste
 285 :            de la structure.
 286 :         4) Effectuer le calcul des contraintes et des deplacements
 287 :            avec le modele total et le materiau qui en decoule.
 288 :         Le calcul de l'integrale avec la procedure G_THETA sera
 289 :         realise sur UN SEUL element en multicouche et pour toutes les
 290 :         couches dans cet element qui ont une epaisseur inferieure a
 291 :         1e-4*(epaisseur totale de la coque). Un tel element doit
 292 :         etre designe par l'argument suivant :
 293 : 
 294 :         SUPTAB.'ELEMENT_MULTICOUCHE'
 295 :         = Objet MAILLAGE comportant UN SEUL element modelise en
 296 :           multicouche. Il doit etre a l'interieur de la zone THETA,
 297 :           c'est a dire dans la zone definie par le nombre SUPTAB.'COUCHE'.
 298 :           Il ne doit pas etre trop loin, ni trop proche de la pointe de
 299 :           la fissure. Theoriquement, l'integrale a calculer est
 300 :           independante du choix de l'element pres de la fissure, ce qui
 301 :           est numeriquement verifiable en la determinant sur des elements
 302 :           en multicouche differents. NOTA : Cette technique necessite un
 303 :           maillage tres fin dans la zone de la pointe de la fissure.
 304 : 
 305 :     10 : Si on souhaite le calcul de criteres de decharge des contraintes
 306 : 
 307 :        SUPTAB.'CALCUL_CRITERE'
 308 :        = LOGIQUE = VRAI si on veut le calcul,
 309 :                  = FAUX si on ne veut pas le calcul.
 310 :          (Valeur par defaut = VRAI)
 311 : 
 312 :     11 : Calcul des facteurs d'intensité des contraintes (DECOUPLAGE)
 313 : 
 314 :         SUPTAB.'METH_AUX'
 315 :         = MOT pour preciser la methode de calcul des champs auxiliaires :
 316 :           1) 'GANA'      pour utiliser les expressions analytiques du
 317 :                          gradient du deplacement et du champ de contraintes
 318 :                          en deformations planes (champs de Westergaard).
 319 :                          Cette methode n'est pas encore disponible dans
 320 :                          le cas ou la fissure est dans une discontinuite
 321 :                          materielle.
 322 :           2) 'UANA'      pour utiliser l'expression analytique du champ de
 323 :                          deplacement en deformations planes ou contraintes
 324 :                          planes. Le gradient du deplacement et le champ de
 325 :                          contraintes sont ensuite deduits en appelant
 326 :                          GRAD et SIGM respectivement.
 327 :                          Cette methode n'est pas encore disponible en 3D.
 328 :           3) 'MECANIQUE' pour resoudre des problemes elementaires en appelant
 329 :                          RESO.
 330 :                          Cette methode n'est actuellement disponible qu'en 3D
 331 :                          avec des elements standards. De plus, celle-ci est
 332 :                          assez couteuse.
 333 :           Lorsque l'indice n'est pas rempli, on utilise pas defaut la methode
 334 :           G-analytique dans tous les cas ou la fissure n'est pas dans une
 335 :           discontinuite materielle, et la methode U-analytique sinon (en 2D).
 336 : 
 337 : 
 338 :     SORTIE :
 339 :     --------
 340 : 
 341 :     Les resultats du calcul correspondant a un champ THETA specifie
 342 :     par l'objet SUPTAB.'COUCHE' (ou SUPTAB.'CHAMP_THETA' dans le cas
 343 :     ou on souhaite donner soi-meme un champ de type Theta) sont sauves
 344 :     de la maniere suivante :
 345 : 
 346 : 
 347 :     Dans tous les cas de calcul
 348 :     ---------------------------
 349 : 
 350 :     SUPTAB.'CHPO_RESULTATS' = Objet decrivant l'evolution du la grandeur
 351 :                               calculee le long du front de fissure.
 352 :                               Son type est variable selon la solution du
 353 :                               probleme traite :
 354 : 
 355 :                               1) solution provenant de l'operateur RESO
 356 :                                   => CHPOINT
 357 : 
 358 :                               2) solution provenant de la procedure PASAPAS
 359 :                                   => TABLE de CHPOINT indicee par
 360 :                                      .(numero du pas de calcul)
 361 : 
 362 :                               La ou les composantes du ou des CHPOINT sont
 363 :                               déterminées par l'objectif :
 364 :                               'J' ou 'J_DYNA' => J
 365 :                               'DJ/DA'         => DJDA
 366 :                               'C*' ou 'C*H'   => C*
 367 :                               'DECOUPLAGE'    => KI, KII (et KIII en 3D)
 368 : 
 369 :     SUPTAB.'RESULTATS' = Objet contenant la valeur numerique du calcul.
 370 :                          Son type est variable selon qu'on est en 2D ou
 371 :                          3D et selon la solution du probleme traite :
 372 : 
 373 :                          1) valeur de l'integrale de contour dans le cas
 374 :                             d'une solution provenant de l'operateur RESO
 375 :                             2D        => FLOTTANT
 376 :                             3D massif => TABLE indicee par
 377 :                                .(fonction d'interpolation) et
 378 :                                .'GLOBAL' pour une estimation globale
 379 :                             3D coque  => TABLE indicee par
 380 :                                .'SUPERI' en peau superieure
 381 :                                .'INFERI' en peau inferieure
 382 :                                .'MEDIAN' au plan median et
 383 :                                .'GLOBAL' pour une estimation globale
 384 : 
 385 :                          2) valeur de l'integrale de contour a un
 386 :                             certain numero du pas de calcul dans le
 387 :                             cas d'une solution provenant de la
 388 :                             procedure PASAPAS
 389 :                             2D        => TABLE indicee par
 390 :                                .(numero du pas de calcul)
 391 :                             3D massif => TABLE indicees par
 392 :                                .(numero du pas de calcul).(fonction
 393 :                                 d'interpolation) et
 394 :                                .(numero du pas de calcul).'GLOBAL'
 395 :                             3D coque  => TABLE indicees
 396 :                                .(numero du pas de calcul).'SUPERI'
 397 :                                .(numero du pas de calcul).'INFERI'
 398 :                                .(numero du pas de calcul).'MEDIAN' et
 399 :                                .(numero du pas de calcul).'GLOBAL'
 400 : 
 401 :                          3) valeur des F.I.C. (facteurs d'intensite des
 402 :                             contraintes) dans le cas de decouplage des
 403 :                             modes avec une solution provenant de
 404 :                             l'operateur RESO
 405 :                             2D        => TABLE indicee par
 406 :                                .'I'  pour KI et
 407 :                                .'II' pour KII
 408 :                             3D massif => TABLE indicee par
 409 :                                .'I'  .(fonction d'interpolation) et
 410 :                                .'I'  .'GLOBAL' pour KI
 411 :                                .'II' .(fonction d'interpolation) et
 412 :                                .'II' .'GLOBAL' pour KII
 413 :                                .'III'.(fonction d'interpolation) et
 414 :                                .'III'.'GLOBAL' pour KIII
 415 : 
 416 :                          4) valeur des F.I.C. (facteurs d'intensite des
 417 :                             contraintes) a un certain numero du pas de
 418 :                             calcul dans le cas de decouplage des modes
 419 :                             avec une solution provenant de la procedure
 420 :                             PASAPAS
 421 :                             2D        => TABLE indicee par
 422 :                                .'I' .(numero du pas de calcul) pour KI et
 423 :                                .'II'.(numero du pas de calcul) pour KII
 424 :                             3D massif => TABLE indicee par
 425 :                                .'I'  .(numero du pas de calcul).(fonction
 426 :                                  d'interpolation) et
 427 :                                .'I'  .(numero du pas de calcul).'GLOBAL'
 428 :                                  pour KI
 429 :                                .'II' .(numero du pas de calcul).(fonction
 430 :                                  d'interpolation) et
 431 :                                .'II' .(numero du pas de calcul).'GLOBAL'
 432 :                                  pour KII
 433 :                                .'III'.(numero du pas de calcul).(fonction
 434 :                                  d'interpolation
 435 :                                  de fissure) et
 436 :                                .'III'.(numero du pas de calcul).'GLOBAL'
 437 :                                  pour KIII
 438 : 
 439 :                          Les fonctions d'interpolation utilisees pour indicer
 440 :                          les TABLE dans le cas 3D massif sont celles qui sont
 441 :                          utilisees pour reconstituer l'intégrale de contour le
 442 :                          long du front de fissure (pour l'instant il s'agit
 443 :                          necessairement des fonctions de forme des segments du
 444 :                          front).
 445 : 
 446 :     Dans le cas de calcul effectue pas a pas
 447 :     ----------------------------------------
 448 : 
 449 :     SUPTAB.'EVOLUTION_RESULTATS' = Objet contenant l'evolution des
 450 :                                    resultats en fonction du temps.
 451 :                                    Son type est variable selon la
 452 :                                    configuration du probleme traite :
 453 : 
 454 :                                    1) Evolution de l'integrale de contour
 455 :                                       2D ou 3D coque => EVOLUTION
 456 :                                               (dans le cas des coques il y a 4
 457 :                                                SOUS-EVOLUTION :
 458 :                                                SUPE, MEDI, INFE et GLOB)
 459 :                                       3D massif      => TABLE indicee par
 460 :                                          .(points au front de fissure)
 461 :                                          .'GLOBAL' evolution pour une
 462 :                                            estimation globale
 463 : 
 464 :                                    2) Evolution des F.I.C. (facteurs
 465 :                                       d'intensite de contrainte)
 466 :                                       2D        => EVOLUTION a 2 SOUS-EVOLUTION
 467 :                                                    KI et KII
 468 :                                       3D massif => TABLE indicee par
 469 :                                               .(points au front de fissure) et
 470 :                                               .'GLOBAL'
 471 :                                               La table contient des evolutions
 472 :                                               a 3 SOUS-EVOLUTION KI, KII et KIII
 473 : 
 474 : 
 475 :     Dans le cas des elements de coque
 476 :     ---------------------------------
 477 : 
 478 :     SUPTAB.'EPAISSEUR_RESULTATS' = representant l'evolution de la valeur
 479 :                                    des integrales dans l'epaisseur de la
 480 :                                    coque. Son type est variable selon la
 481 :                                    solution du probleme traite :
 482 : 
 483 :                                   1) EVOLUTION dans le cas d'une solution
 484 :                                      provenant de l'operateur RESO
 485 :                                   2) TABLE indicee par
 486 :                                      .(numero du pas de calcul)
 487 :                                      dans le cas d'une solution provenant
 488 :                                      de la procedure PASAPAS
 489 : 
 490 : 
 491 :     Dans le cas de calcul elasto-plastique
 492 :     --------------------------------------
 493 : 
 494 :      En cas de calcul elasto-plastique isotrope ou cinematique,
 495 :      eventuellement thermique :
 496 : 
 497 :     SUPTAB.'CRIT_DECHA_GLOBAL1'
 498 :     = On calcul un critere de decharge des contraintes defini par
 499 :       (si F = courbe de traction ) :   crit = F(EPSeq)/ SIGeq.
 500 :       - crit = 1. si non-decharge,
 501 :       - et crit > 1. si decharge.
 502 :       SUPTAB.'CRIT_DECHA_GLOBAL1' est une table indicee par les numeros
 503 :       du pas de calcul contenant des reels.
 504 :     SUPTAB.'CRIT_DECHA_GLOBAL2'
 505 :     = On calcul un critere de decharge et de changement de direction des
 506 :       contraintes defini par : crit = 2 - SIG:EPSpl/SIGeq_max.EPSE.
 507 :       - crit = 1. si non-decharge et non-changement de direction,
 508 :       - crit > 1. dans le cas contraire.
 509 :       SUPTAB.'CRIT_DECHA_GLOBAL2' est une table indicee par les numeros
 510 :       du pas de calcul contenant des reels.
 511 :     SUPTAB.'CRIT_DECHA_LOCAL1'
 512 :     = On calcul un critere de decharge defini par: crit = (delta SIG)/SIG
 513 :       - crit = 0. si non-decharge,
 514 :       - crit > 0. si decharge.
 515 :       SUPTAB.'CRIT_DECHA_LOCAL1' est une table indicee par les numeros
 516 :       du pas de calcul contenant des CHPOINTs.
 517 :     SUPTAB.'CRIT_DECHA_LOCAL2'
 518 :     = On calcul un critere de changement de direction des contraintes
 519 :       defini par: crit = SIG:(delta SIG)/norme(SIG).norme(delta SIG)
 520 :       - crit = 1. si la contrainte evolue proportionnellement et de facon
 521 :         monotone,
 522 :       - crit = -1. si la contrainte evolue proportionnellement mais dans la
 523 :         direction opposee
 524 :       SUPTAB.'CRIT_DECHA_LOCAL2' est une table indicee par les numeros
 525 :       du pas de calcul contenant des CHPOINTs.
 526 : 
 527 : 
 528 :      Remarque
 529 :      --------
 530 : 
 531 :      La table SUPTAB contient aussi d'autres objets servant aux
 532 :      reprises. C'est cette table qu'il convient de sauver en vue
 533 :      d'une reprise ulterieure du calcul.
 534 : 
 535 : 

© Cast3M 2003 - Tous droits réservés.
Mentions légales