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