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