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

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