Télécharger g_theta.notice

Retour à la liste thématique

Afficher cette notice en

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

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