Télécharger g_theta.notice

Retour à la liste thématique

Afficher cette notice en

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

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