Télécharger g_theta.notice

Retour à la liste thématique

Afficher cette notice en

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

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