Afficher cette notice en
1 : $$$$ PASAPAS NOTICE MB234859 24/08/22 21:15:02 11987 2 : DATE 24/08/22 3 : 4 : Procedure PASAPAS Voir aussi : PECHE 5 : _________________ EXPLORER 6 : CHARTHER 7 : PAS_SAUV 8 : 9 : TAB1. ACCELERATIONS MODELE 10 : AMORTISSEMENT MOVA 11 : AUGMENTATION_AUTOMATIQUE MTOL 12 : AUTOCRIT NB_BOTH 13 : AUTORESU NITERINTER_MAX 14 : AUTOMATIQUE NITER_KTANGENT 15 : AUTOPAS NMAXSUBSTEPS 16 : BCSTH NPAS_TRACKING 17 : BLOCAGES_DIFFUSIONS NRMAX 18 : BLOCAGES_MECANIQUES ORDRE 19 : BLOCAGES_THERMIQUES PARAMETRE_DE_PILOTAGE 20 : CAPACITE_CONSTANTE PAS_AJUSTE 21 : CARACTERISTIQUES PAS_MAX 22 : CELSIUS PILOTAGE_INDIRECT 23 : CHARGEMENT 24 : CONCENTRATIONS PRECISINTER 25 : CONDUCTIVITE_CONSTANTE PRECISION 26 : PRECSOUSITERATION 27 : CONN PREDICTEUR 28 : CONSOLIDATION PROCEDURE_CHARMECA 29 : CONTRAINTES PROCEDURE_CHARTHER 30 : CONVERGENCE_FORCEE PROCEDURE_PARATHER 31 : CONVERGENCE_MEC_THE PROCEDURE_PERSO1 32 : CONVERGENCE_MONOTONE PROCEDURE_PERSO2 33 : CRITERE_COHERENCE PROCEDURE_REEV_MEC 34 : CTE_STEFAN_BOLTZMANN PROCEDURE_REEV_THE 35 : CTOL PROCEDURE_THERMIQUE 36 : DEFORMATIONS_INELASTIQUES PROCESSEURS 37 : DELTAITER PROJECTION 38 : DEPLACEMENTS PROPORTIONS_PHASE 39 : DEPLACEMENTS_PILOTES REACTIONS 40 : DYNAMIQUE REACTIONS_DIFFUSIONS 41 : ECONOMIQUE REACTIONS_THERMIQUES 42 : FEFP_FORMULATION REAC_GRANDS 43 : FORCES_PILOTEES REEQUILIBRAGE 44 : FREA1 RELAXATION_DUPONT 45 : FTOL RELAXATION_NONCONV 46 : GRANDS_DEPLACEMENTS RELAXATION_THETA 47 : HYPOTHESE_DEFORMATIONS REPRISE 48 : INITIALISATION RENORMALISATION 49 : K_SIGMA RIGIDITE_AUGMENTEE 50 : K_TANGENT RIGIDITE_CONSTANTE 51 : K_TANGENT_ITER0 SOUS_INCREMENT 52 : K_TANGENT_PERT SOUS_RELAXATION 53 : K_TANGENT_SYME STABILITE 54 : K_TANG_PERT_C1 SUBSTEPPING 55 : LAGRANGIEN TEMPERATURES 56 : LBC TEMPS 57 : LINESEARCH TEMPS_ADAPTATION_MODELE 58 : MAN TEMPS_CALCULES 59 : MASSE_CONSTANTE TEMPS_SAUVEGARDES 60 : MAXDEFOR TEMPS_SAUVES 61 : MAXSOUSPAS TRACKING 62 : MAXISOUSPAS TYP_TRAC 63 : MAXITERATION TTOL 64 : MAXSOUSITERATION 65 : MES_SAUVEGARDES UPDATE_LAGRANGIAN 66 : VARIABLES_INTERNES 67 : VITESSES 68 : ZONE_DE_PILOTAGE 69 : 70 : 71 : 72 : 73 : Objet : 74 : _______ 75 : 76 : MECANIQUE : 77 : 78 : Cette procedure permet d'effectuer un calcul non lineaire incremental 79 : La non linearite peut provenir, soit du materiau (plasticite), soit 80 : des grands deplacements soit des deux a la fois. 81 : Les resultats sont calcules a des valeurs du parametre d'evolution 82 : (pseudo temps ou temps reel) definies par l'utilisateur. 83 : Sous l'option MODE FREQ, la procedure resout l'équation dynamique 84 : sur la base modale étendue selon l'approche dite spectrale ou 85 : fréquentielle. Implicitement un instant est interprété comme une 86 : fréquence pour l'objet CHARGEMENT sans changement de terminologie. 87 : La procédure propose un balayage de fréquences par défaut lorsque la 88 : liste des TEMPS_CALCULES n'est pas précisée. 89 : 90 : 91 : THERMIQUE : 92 : 93 : Cette procedure permet d'effectuer un calcul lineaire et 94 : non-lineaire en tenant compte de la conduction, de la convection 95 : et du rayonnement. 96 : 97 : DIFFUSION : 98 : 99 : Cette procedure permet de resoudre un proble lineaire ou non-lineaire 100 : de diffusion. 101 : 102 : Il est possible d'effectuer un calcul couplant MECANIQUE, THERMIQUE 103 : et DIFFUSION. THERMIQUE et DIFFUSION sont resolues simulatnement. 104 : 105 : Commentaire : 106 : _____________ 107 : 108 : En entree, TAB1 sert a definir les options et les parametres du calcul. 109 : Les indices de l'objet TAB1 sont des mots (a ecrire en toutes lettres, 110 : et en majuscules s'ils sont mis entre cotes) dont voici la liste : 111 : 112 : BLOCAGES_DIFFUSIONS : blocages de diffusion (type RIGIDITE ou 113 : CHARGEMENT de nom BLOD). 114 : 115 : BLOCAGES_MECANIQUES : blocages mecaniques (type RIGIDITE ou 116 : CHARGEMENT de nom BLOM). 117 : 118 : BLOCAGES_THERMIQUES : blocages thermiques (type RIGIDITE ou 119 : CHARGEMENT de nom BLOT). 120 : 121 : CARACTERISTIQUES : Champ de caracteristiques materielles et 122 : eventuellement geometriques si necessaire 123 : (type MCHAML, sous-type CARACTERISTIQUES, 124 : ou CHARGEMENT de nom MATE) 125 : Ses composantes peuvent etre de type : 126 : 1) FLOTTANT si la composante est 127 : est constante sur toute la 128 : structure; 129 : 2) MCHAML si la composante depend 130 : uniquement des points de la 131 : structure; 132 : 3) EVOLUTION si la composante 133 : varie en fonction d'un seul 134 : parametre. 135 : 4) NUAGE si la composante est 136 : decrite par une courbe de type 137 : EVOLUTION dependant d'un seul 138 : parametre. 139 : 140 : CHARGEMENT : definition du chargement en fonction du parametre 141 : d'evolution (type CHARGEME) (cree par l'operateur CHAR). 142 : Les chargements elementaires sont obligatoirement 143 : nommes: 144 : 145 : - la temperature T 146 : - les concentrations imposees CIMP 147 : - les deplacements imposes DIMP 148 : - les temperatures imposees TIMP 149 : - Les autres chargements (meca) MECA 150 : - Flux (en consolidation) FLUX 151 : - Les flux de diffusion (diffusion) QCO 152 : - Les flux de chaleur (thermique) Q 153 : - Les temperatures ext (convection) TECO 154 : - Les temperatures ext (rayonnement)TERA 155 : - Une deformation imposee DEFI 156 : - Des blocages mecaniques BLOM 157 : - Des blocages thermiques BLOT 158 : - Des blocages diffusions BLOD 159 : - Le modele MODE 160 : - Ses carateristiques materielles MATE 161 : - Des parametres externes de nom... MOT1 162 : 163 : 164 : "Les autres chargements" concernent uniquement la mecanique 165 : et representent un champ a ajouter au second membre 166 : (type = CHPOINT) 167 : 168 : "Les flux de diffusion" ont pour nom de composante QCO par 169 : defaut. Toutefois, l'utilisateur a la possibilite de le 170 : surcharger (voir notice operateur MODE). Dans ce cas, 171 : il convient de mettre le nom de composante choisi. 172 : 173 : "Les variables exterieures" sont les variables dont la donnee 174 : est indispensable pour instancier le champ de materiau 175 : (type conseille = MCHAML). Il peut s'agir en particulier des 176 : parametres externes du modele, s'il en existe. 177 : 178 : Dans le cas de chargement devant etre evalues a la fin du pas 179 : (par exemple les pressions suiveuses, le rayonnement) 180 : il faut utiliser la possibilite offerte par PROCEDURE_CHARMECA 181 : ou PROCEDURE_CHARTHER. 182 : 183 : Dans le cas de "variables exterieures" devant etre mises a jour 184 : au cours des iterations thermiques, afin de prendre en compte 185 : leur evolution sur le champ de materiau, il faut utiliser la 186 : possibilite offerte par PROCEDURE_PARATHER. 187 : 188 : 189 : PROCEDURE_CHARMECA : Logique. VRAI s'il faut evaluer une partie 190 : du chargement mecanique au cours des iterations 191 : d'equilibrage. Ce calcul se fera a travers une 192 : procedure CHARMECA dont l'appel est : 193 : TAB2 = CHARMECA TAB1 TIME; 194 : TAB1 est la table passee a PASAPAS et TIME est 195 : l'instant pour lequel on veut la charge 196 : TAB2.'ADDI_SECOND'. Il faut evaluer le chargement 197 : sur la configuration courante. Voir exemples 198 : gdep2 et gdep3. Ils montrent comment modeliser 199 : une pression suiveuse. 200 : 201 : PROCEDURE_CHARTHER : Logique. VRAI s'il faut evaluer une partie du 202 : chargement thermique au cours des iterations pour 203 : le calcul d'un pas de temps.Ce calcul se fera a 204 : travers une procedure CHARTHER dont l'appel est : 205 : TAB2 = CHARTHER TAB TIME; 206 : TAB est la table passee a PASAPAS et TIME est 207 : l'instant pour lequel on veut la charge 208 : TAB2.'ADDI_SECOND'. Il faut evaluer le chargement 209 : sur la configuration courante. 210 : 211 : PROCEDURE_PARATHER : Logique a VRAI s'il faut mettre a jour des 212 : "variables exterieures" a chaque iteration d'un 213 : calcul thermique. Cette mise a jour se fera via 214 : la procedure PARATHER dont l'appel est "PARATHER 215 : TAB1 INST ;", TAB1 etant la table passee a 216 : PASAPAS et INST l'instant pour lequel on evalue 217 : les "variables exterieures" necessaires a la 218 : thermique. (Un exemple de procedure PARATHER est 219 : disponible dans le cas-test 220 : "exemple_parather.dgibi".) 221 : 222 : 223 : TEMPS_ADAPTATION_MODELE : LISTREEL, instants ou le maillage du modele 224 : est adapte (voir ADAPTE). 225 : 226 : TEMPS_CALCULES : definition des valeurs du parametre d'evolution 227 : (ou du temps) pour lesquelles on effectue le calcul 228 : (type LISTREEL) (cree par l'operateur PROG). 229 : Dans le cas ou cette donnee est absente, le code 230 : passe en ajustement automatique de pas et considere 231 : comme liste des temps a calculer la liste des temps 232 : sauves donnee sous l'indice TEMPS_SAUVES. 233 : (voir PAS_AJUSTE) 234 : 235 : TEMPS_SAUVES : definition des valeurs du 236 : (ou du temps) pour lesquelles on met les resultats dans 237 : les tables de resultats 238 : (type LISTREEL) (cree par l'operateur PROG) (facultatif) 239 : 240 : TEMPS_SAUVEGARDES : definition des valeurs du parametre d'evolution 241 : pour lesquelles PASAPAS appelle l'operateur SAUV 242 : avec, en argument, la table passee a PASAPAS. La 243 : liste fournie doit etre une sous liste des l'indice 244 : TEMPS_SAUVES. Il faut avoir defini le fichier de 245 : sauvegarde (par OPTION SAUV '...';) avant l'appel 246 : a PASAPAS. 247 : REMARQUE IMPORTANTE : la sauvegarde n'est effectuee 248 : que s'il s'est ecoule plus de 5 minutes depuis la 249 : precedente sauvegarde afin de limiter les acces disque. 250 : MES_SAUVEGARDES : Table permettant de demander de garder certains 251 : resultats en plus des ceux sauves automatiquement. 252 : 253 : Les indices permis sont : 254 : 255 : -'DEFIN'=VRAI pour les deformations inelastiques 256 : -'DEFTO'=VRAI pour les deformations. 257 : Attention, pas necessairement egal aux 258 : deformations totales (voir 'DEFORMATIONS'). 259 : -'DEFAP'=VRAI pour les deformations initiales de 260 : l'increment d'apport de matiere 261 : (necessite chargement TAPP) 262 : 263 : Lorsque la table MES_SAUVEGARDES est definie, les 264 : indices non renseignes sont pris egal a FAUX. 265 : 266 : On peut de plus utiliser cette table pour faire 267 : transiter vers la procedure PAS_SAUV des informations 268 : supplementaires. Cette procedure, qui realise la mise 269 : des resultats dans les tables de resultats, peut etre 270 : surchargee par l'utilisateur (voir la notice de la 271 : procedure PAS_SAUV) 272 : 273 : ECONOMIQUE : Logique VRAI si on veut que pasapas recupere la place 274 : memoire, prise par les resultats, au fur et a mesure 275 : des sauvegardes. 276 : 277 : MODELE : Modele (objet MMODEL ou CHARGEMENT de nom MODE) englobant 278 : toute la structure, excepte celle correspondant a la 279 : RIGIDITE_CONSTANTE. 280 : 281 : indice type objet pointe commentaires 282 : 283 : HYPOTHESE_DEFORMATIONS MOT Indique l'hypothese faite sur le calcul du tenseur des deformations. 284 : Les valeurs possibles sont : 285 : LINEAIRE (par defaut) : deformation infinitesimale 286 : QUADRATIQUE : deformation de Green-Lagrange 287 : TRUESDELL : deformation de Truesdell 288 : JAUMANN : deformation de Jaumann 289 : UTILISATEUR : deformation définie par l'utilisateur 290 : GRANDS_DEPLACEMENTS LOGIQUE VRAI pour reactualiser la geometrie a chaque 291 : iteration dans le cas des grands deplacements. 292 : En standard, cette option entraine : 293 : -- K_SIGMA = VRAI 294 : peut etre impose a FAUX en l'ecrivant 295 : -- HYPOTHESE_DEFORMATIONS = QUADRATIQUE 296 : peut etre modifie en l'ecrivant 297 : -- La sauvegarde des deformations totales. 298 : DYNAMIQUE LOGIQUE VRAI si calcul dynamique 299 : TRACKING LOGIQUE VRAI si le trcaking est active 300 : (methode EFEM). 301 : TYP_TRAC MOT FIXED_HC : tracking par resolution 302 : d un probleme de conduction 303 : VARIABLE_HC : tracking par resolution 304 : d un probleme de convection/diffusion 305 : isotrope 306 : VARIABLE_HCD_SU : tracking par resolution 307 : d un probleme de convection/diffusion 308 : anisotrope 309 : MASSE_CONSTANTE RIGIDITE matrice de masse constante. 310 : Par exemple masses lumpees. 311 : La masse coherente est calculee 312 : automatiquement a partir du modele 313 : et du champ de materiau 314 : AMORTISSEMENT RIGIDITE matrice d'amortissement, 315 : facultative en calcul dynamique 316 : RIGIDITE_AUGMENTEE RIGIDITE matrice d'augmentation de la raideur. 317 : N'intervenant pas dans le calcul 318 : du residu, elle ne change pas la 319 : solution mais peut influer sur la 320 : convergence. 321 : La masse ou l'amortissement sont 322 : classiquement utilisees. 323 : AUGMENTATION_AUTOMATIQUE LOGIQUE Si la rigidite augmentee est fournie, calcul 324 : automatique du coefficient multiplicateur. 325 : Sinon calcul automatique de la matrice 326 : d'augmentation. 327 : RIGIDITE_CONSTANTE RIGIDITE matrice de rigidite constante ajoutee 328 : a celle de la structure. Par exemple, 329 : une rigidite decrivant des appuis 330 : elastiques (voir operateur APPUI). 331 : Le champ de forces nodales associe 332 : est calcule en multipliant le champ 333 : de deplacement par cette matrice. 334 : CAPACITE_CONSTANTE RIGIDITE matrice de capacite ajoutee a celle 335 : du modele thermique. 336 : Les flux thermiques nodaux associes 337 : sont calcules en multipliant le champ 338 : de temperature par cette matrice. 339 : Pour les problemes de DIFFUSION, une 340 : matrice de capacite de diffusion peut 341 : aussi etre ajoutee par cette entree. 342 : CONDUCTIVITE_CONSTANTE RIGIDITE matrice de conductivite ajoutee a 343 : celle du modele thermique. 344 : Les flux thermiques nodaux associes 345 : sont calcules en multipliant le champ 346 : de temperature par cette matrice. 347 : Pour les problemes de DIFFUSION, une 348 : matrice de diffusivite peut aussi etre 349 : ajoutee par cette entree. 350 : PILOTAGE_INDIRECT LOGIQUE VRAI si on desire activer le pilotage 351 : indirect du chargement mecanique. Il 352 : est possible de specifier differents 353 : equations de pilotage en modifiant la 354 : procedure PILOINDI. 355 : Voir l'exemple : pilotage_indirect_1.dgibi 356 : DEPLACEMENTS_PILOTES CHPOINT Pour le PILOTAGE_INDIRECT, il s'agit du 357 : CHPOINT, issu de DEPI, indiquant la 358 : direction du chargement en deplacements 359 : imposes (peut etre unitaire) 360 : FORCES_PILOTEES CHPOINT Pour le PILOTAGE_INDIRECT, il s'agit du 361 : CHPOINT de forces nodales indiquant la 362 : direction du chargement en forces 363 : imposees (peut etre unitaire) 364 : PARAMETRE_DE_PILOTAGE EVOLUTION Evolution en fonction du temps du parametre 365 : de pilotage a etre respecte via PILOINDI. 366 : ZONE_DE_PILOTAGE MAILLAGE Pour le PILOTAGE_INDIRECT, il s'agit de 367 : preciser la region ou la deformation 368 : maximale doit etre cherchee (valable pour 369 : le critere de pilotage standard de PILOINDI). 370 : Cet indice est facultatif. S'il n'est pas 371 : renseigne, l'option par default est activee 372 : (recherche sur l'ensemble du maillage). 373 : AUTOMATIQUE LOGIQUE VRAI si on desire piloter le 374 : chargement en fonction d'un 375 : critere sur le champ de 376 : deplacement. Par defaut le pas 377 : de charge est choisi pour que 378 : l'increment de deformation 379 : correspondant soit proche de la 380 : valeur TAB1.'AUTOCRIT'que l'on 381 : nomme le critere de pilotage. 382 : Il est possible de specifier un 383 : autre type de pilotage en modi 384 : fiant la procedure AUTOPILO . 385 : Dans ce cas il faut encore 386 : fournir le listreel LREEL1 mais 387 : le programme limite l'increment 388 : du temps de facon a respecter 389 : le critere sur epsilon total. 390 : AUTOCRIT FLOTTANT critere de pilotage (par defaut 391 : 0.001) 392 : AUTORECA ENTIER nombre d'iteration entre chaque recalcul du 393 : coefficient de pilotage (par defaut 1) 394 : AUTOPAS ENTIER nombre maximum de sous pas 395 : autorises pour atteindre le 396 : temps demande. 397 : AUTORESU ENTIER valeur precisant tous les combien de 398 : pas il faut garder les resultats 399 : dans le cas de calcul AUTOMATIQUE. 400 : BCSTH RIGIDITE conditions aux limites pour le 401 : tracking (methode EFEM) 402 : PAS_AJUSTE LOGIQUE permet d'imposer un choix de 403 : strategie de calcul de pas de temps. 404 : VRAI on veut faire des pas de temps 405 : ajustes en fonction de la difficulte 406 : de convergence du pas precedent. 407 : CONVERGENCE_FORCEE LOGIQUE VRAI si on force la convergence en 408 : effectuant des sous-pas non converges 409 : (le defaut) 410 : MAXDEFOR FLOTTANT increment de deformation maximum en 411 : convergence forcee (1e-3 par defaut) 412 : CONSOLIDATION LOGIQUE VRAI si calcul de consolidation 413 : DELTAITER ENTIER nombre de pas d'ecart pour le test 414 : de non convergence 415 : K_SIGMA LOGIQUE VRAI si on souhaite ajouter 416 : la matrice KSIGMA a l'operateur 417 : d'iterations. 418 : K_TANGENT LOGIQUE iterations avec la matrice tangente 419 : coherente (appel a KTAN) (FAUX par 420 : defaut) 421 : K_TANGENT_ITER0 MOT mot indiquant la matrice a utiliser 422 : pour demarrer les iterations 423 : mecaniques d'un pas 'MAT_ELASTIQUE' 424 : rigidite elastique 'MAT_TANGENTE' 425 : rigidite tangente en debut de pas 426 : (avec DT = 0.). Par defaut, on 427 : utilise la matrice calculee au pas 428 : precedent. 429 : K_TANGENT_PERT LOGIQUE iterations avec la matrice tangente 430 : calculee par perturbation (appel a 431 : KTAN PERT) (FAUX par defaut) 432 : K_TANGENT_SYME LOGIQUE La matrice tangente calculee est 433 : symetrisee (FAUX par defaut) 434 : K_TANG_PERT_C1 FLOTTANT matrice tangente par perturbation 435 : coefficient multiplicatif applique a 436 : l'increment de deformation pour 437 : determiner la perturbation 438 : (1.D-3 par defaut) 439 : K_TANG_PERT_C2 FLOTTANT matrice tangente par perturbation 440 : valeur minimale de l'increment de 441 : deformation perturbe (C1/100. par 442 : defaut) 443 : NITER_KTANGENT ENTIER la matrice tangente est recalculee 444 : toutes les NITER_KTANGENT iterations 445 : (1 par defaut) 446 : SUBSTEPPING LOGIQUE Substepping pour l'integration locale 447 : avec la matrice tangente consistante. 448 : Modeles: J2, RH_COULOMB, MRS_LADE. 449 : NMAXSUBSTEPS ENTIER nombre maximum de sous-pas pour 450 : SUBSTEPPING. 451 : NITERINTER_MAX ENTIER nombre max d'iterations au niveau 452 : local (avec SUBSTEPPING). Modeles: 453 : J2, RH_COULOMB, MRS_LADE. 454 : PRECISINTER FLOTTANT precision pour le probleme local 455 : (integration des lois constitutives) 456 : 1D-8 par defaut 457 : FEFP_FORMULATION LOGIQUE retour exponentiel avec line-search 458 : au niveau local et matrice tangente 459 : consistante. Modeles: VMT_FEFP, 460 : RHMC_FEFP, POWDER_FEFP, 461 : POWDERCAP_FEFP. 462 : UPDATE_LAGRANGIAN LOGIQUE VRAI pour le traitement de FEFP par 463 : Lagrangien Augmente (par defaut) 464 : PRECISION FLOTTANT Valeur de la precision du calcul, 465 : quelle que soit la physique (par 466 : default 1e-4). 467 : FTOL FLOTTANT l'equilibre des efforts sera 468 : verifie a FTOL pres 469 : MTOL FLOTTANT l'equilibre des moments sera 470 : verifie a MTOL pres 471 : TTOL FLOTTANT l'equilibre en THERMIQUE sera verifie a TTOL pres (Valeur absolue) 472 : sur les variables primales 473 : CTOL FLOTTANT/CHPOINT l'equilibre en DIFFUSION sera verifie a CTOL pres (Valeur absolue) 474 : sur les variables primales 475 : PROCESSEURS MOT 'MONO_PROCESSEUR' pour imposer une 476 : execution avec un seul processeur 477 : 'COMPORTEMENT' pour imposer que seul 478 : le comportement est envoye sur 479 : plusieurs processeurs 'AUTOMATIQUE' 480 : pour imposer, l'utilisation maximale 481 : des processeurs. 482 : PREDICTEUR HPP on utilise un calcul HPP pour 483 : initialiser le calcul en grands 484 : deplacements. 485 : LINESEARCH LOGIQUE VRAI pour accéleration grands 486 : déplacements 487 : STABILITE LOGIQUE VRAI pour test stabilité 488 : ACCELERATION ENTIER on utilise l'acceleration de 489 : convergence tous les 490 : TAB1.'ACCELERATION' pas 491 : (2 par defaut). 492 : MAN LOGIQUE VRAI pour demander l'utilisation de 493 : la methode asymptotique numerique 494 : comme predicteur du deplacements en 495 : cas d'etudes en grands deplacements. 496 : Pour que cela soit vraiment interes- 497 : -sant il faut un comportement 498 : elastique lineaire des materiaux. 499 : LBC MAILLAGE support geometrique des conditions 500 : aux limites pour l'algorithme de 501 : tracking. 502 : ORDRE ENTIER prermet de preciser l'ordre du 503 : developpement de la MAN. Par defaut 504 : il vaut 12. 505 : MAXITERATION ENTIER nombre maximum d'iterations 506 : valeur prise par defaut: 49 507 : MAXISOUSPAS ENTIER nombre maximum de sous- 508 : increments en viscoplasticite ou 509 : en fluage. 510 : MAXSOUSITERATION ENTIER nombre maximum de sousiteration pour 511 : la resolution des equations d Helmholtz 512 : non lineaires utilisables en non local 513 : PRECSOUSITERATION FLOTTANT precision relative pour la convergence 514 : de l equation d HelmHoltz utilisable 515 : en non local 516 : MAXSOUSPAS ENTIER nombre maximum de sous-pas 517 : en convergence forcee 518 : valeur prise par defaut: 200 519 : NPAS_TRACKING FLOTTANT Nombre de pas a considerer pour 520 : realiser le tracking de la fissure 521 : (method EFEM) 522 : NRMAX FLOTTANT Nombre maximum d elements susceptibles 523 : de voir une fissure s initier 524 : (methode EFEM). 525 : CONN MCHAML Champ de connectivites construit par 526 : l'operateur CONN, pour tous les cas 527 : incluant des symetries.a taille de 528 : la plus grande heterogeneite). 529 : MOVA MOT Precise le nom de la variable interne 530 : sur laquelle on teste le nombre de 531 : points plastifies (EPSE par defaut) 532 : CONTRAINTES.0 MCHAML contraintes au debut du pas 533 : (0. par defaut) 534 : NB : cf. remarques. 535 : DEPLACEMENTS.0 CHPOINT deplacements au debut du pas 536 : (0. par defaut) 537 : NB : cf. remarques. 538 : VARIABLES_INTERNES.0 MCHAML variables internes au debut du pas 539 : (0. par defaut) 540 : NB : cf. remarques. 541 : DEFORMATIONS_INELASTIQUES.0 deformations inelastiques au debut 542 : MCHAML du pas (0. par defaut) 543 : NB : cf. remarques. 544 : VITESSES.0 CHPOINT vitesses au debut du pas en 545 : dynamique (0. par defaut) 546 : ACCELERATIONS.0 CHPOINT accelerations au debut du pas en 547 : dynamique (0. par defaut) 548 : TEMPS.0 FLOTTANT temps au debut du pas 549 : (0. par defaut) 550 : FREA1 CHPOINT forces de reactions au debut du pas 551 : (0. par defaut). Ce champ n'est pris 552 : en compte que lorsque un des indices 553 : DYNAMIQUE ou CONSOLIDATION de la 554 : table TAB1 est VRAI. 555 : SOUS_INCREMENT ENTIER Nombre de sous increment, pour l'integration 556 : de la loi de comportement, par pas de 557 : calcul (1 par defaut). 558 : En dynamique, si la contrainte varie 559 : fortement pendant le pas il peut etre 560 : important de suivre au mieux les 561 : contraintes afin d'adapter 562 : l'algorithme dynamique. 563 : CONCENTRATIONS.0 CHPOINT concentrations au debut du pas. 564 : (0. par defaut pour le 1er pas). 565 : TEMPERATURES.0 CHPOINT temperature au debut du pas. 566 : (0. par defaut pour le 1er pas). 567 : PROCEDURE_THERMIQUE MOT Nom de la procedure thermique a 568 : utiliser : 569 : NONLINEAIRE : nonlineaire a un pas de temps 570 : (theta-methode) 571 : LINEAIRE : lineaire (theta-methode) 572 : DUPONT : nonlineaire a deux pas de 573 : temps (methode DUPONT2) 574 : PAS_MAX FLOTTANT pas maximum autorise en thermique. 575 : 576 : CTE_STEFAN_BOLTZMANN FLOTTANT La valeur de la constante de Stefan- 577 : Boltzmann est egale par defaut a 578 : 5.673E-8 (la valeur en unites SI). Si le 579 : calcul est fait dans un autre systeme 580 : d'unites il convient de mettre ici la 581 : bonne valeur de cette constante 582 : CELSIUS LOGIQUE VRAI si le calcul de fait en degres 583 : Celsius. Lorsque cette variable est VRAI 584 : on rajoute 273 aux champs de temperature 585 : avant d'en deduire des quantites liees 586 : au rayonnement. 587 : 588 : RELAXATION_DUPONT FLOTTANT Valeur du coefficient de relaxation 589 : pour DUPONT2 (0.25 pas defaut) 590 : RELAXATION_NONCONV FLOTTANT Valeur du coefficient de relaxation 591 : pour la non convergence (1 par defaut) 592 : qui determine le sou-pas non converge 593 : RELAXATION_THETA FLOTTANT Valeur du coefficient de relaxation 594 : pour theta-methode (1.0 par defaut) 595 : SOUS_RELAXATION FLOTTANT Valeur du coefficient de sous- 596 : relaxation (1.0 par defaut) 597 : 598 : PROCEDURE_PERSO1 LOGIQUE VRAI si on doit appeler une procedure 599 : utlisateur apres le calcul mecanique. 600 : Nom de cette procedure : PERSO1. 601 : ATTENTION ! Cette procedure est appe- 602 : lee apres chaque increment, mais non 603 : au cours des iterations. 604 : Voir aussi les remarques. 605 : PROCEDURE_PERSO2 LOGIQUE VRAI si on doit appeler une procedure 606 : utlisateur apres le calcul thermique. 607 : Nom de cette procedure : PERSO2. 608 : ATTENTION ! Cette procedure est appe- 609 : lee apres chaque increment, mais non 610 : au cours des iterations. 611 : Voir aussi les remarques. 612 : 613 : PROJECTION LOGIQUE VRAI si on effectue un calcul couple 614 : et si la mecanique et la thermique ne 615 : s'appuient pas sur le meme maillage. 616 : 617 : PROCEDURE_REEV_MEC LOGIQUE VRAI si on utilise une procedure de 618 : reactualisation apres le calcul 619 : mecanique (par ex pour reevaluer les 620 : donnees du calcul thermique si elles 621 : dependent de la configuration 622 : mecanique) 623 : PROCEDURE_REEV_THE LOGIQUE VRAI si on utilise une procedure de 624 : reactualisation apres le calcul 625 : thermique (par ex pour reevaluer des 626 : donnees du calcul mecanique si elles 627 : dependent de la configuration 628 : thermique). 629 : 630 : CONVERGENCE_MEC_THE LOGIQUE VRAI si on souhaite iterer le schema 631 : thermique mecanique en cas de depen- 632 : -dance mutuellle des deux problemes. 633 : FAUX par defaut. 634 : Voir aussi les remarques. 635 : CRITERE_COHERENCE FLOTTANT critere de convergence si on itere 636 : l'alternance thermique-mecanique en 637 : cas de dependance mutuelle. Le test 638 : est effectue sur les resulats de la 639 : thermique et il porte sur l'erreur 640 : relative maximale entre les deux 641 : dernieres solutions. 642 : Vaut 'PRECISION' par defaut. 643 : NB_BOTH ENTIER Nombre d'iterations maximale de la 644 : boucle de coherence 645 : mecanique/thermique. 646 : Par defaut : 10 647 : 648 : REAC_GRANDS FLOTTANT increment de deformation pour lequel 649 : on reactualise la matrice de raideur 650 : en grands_deplacements 651 : (valeur par defaut : 10.E-2) 652 : CONVERGENCE_MONOTONE LOGIQUE VRAI si on souhaite imposer la 653 : decroissance monotone du residu au 654 : cours des iterations en recalculant 655 : la matrice de raideur et en limitant 656 : l'itere. 657 : (FAUX par defaut) 658 : INITIALISATION LOGIQUE FAUX si on ne souhaite pas que le 659 : premier residu soit initialise avec 660 : les resultats du pas precedent. 661 : (VRAI par defaut) 662 : RENORMALISATION LOGIQUE VRAI si on souhaite limiter l'itere 663 : a une valeur inferieure a MAXDEFOR. 664 : (FAUX par defaut) 665 : LAGRANGIEN MOT MI_PAS (ou REACTUALISE) (par defaut) pour appliquer le comportement 666 : sur la geometrie au milieu du pas. Meilleur ecoulement plastique. 667 : FIN_PAS pour appliquer le comportement sur la geometrie en fin 668 : du pas. Meilleures contraintes. 669 : TOTAL pour appliquer le comportement sur la geometrie initiale. 670 : Meilleur retour a zero des contraintes si on a un 671 : chargement cyclique. 672 : 673 : REPRISE FLOTTANT instant depuis lequel on souhaite reprendre le calcul 674 : 675 : REEQUILIBRAGE LOGIQUE Ce mot-cle ne sert que pour la poursuite/reprise d'un calcul. 676 : Il ne s'applique qu'a la mecanique. 677 : VRAI : si le dernier instant calcule fait partie de la liste 678 : TEMPS_CALCULES alors un pas de reequilibrage est fait. Le calcul 679 : reprend sur un pas de temps de taille nulle. 680 : 681 : En sortie, TAB1 permet de retrouver les resultats. Ceux 682 : ci sont mis dans des tables dont les indices sont des 683 : entiers (0 1 2 ... N) correspondants aux numero de 684 : sauvegarde des resultats (0 : temps initial). L'utilisateur peut 685 : intervenir de deux façons sur le contenu des resultats. Soit il utilise 686 : la table referencee dans l'indice 'MES_SAUVEGARDES' soit il surcharge 687 : directement la procedure PAS_SAUV. 688 : 689 : 690 : indice type objet commentaires 691 : pointe 692 : 693 : ERREUR LOGIQUE logique de valeur VRAI en cas 694 : d'erreur au cours de l'execution 695 : de la procedure. 696 : CONV LOGIQUE VRAI si calcul converge avant ou 697 : pour le nombre maximum d'itera- 698 : tions, FAUX sinon. 699 : TEMPS TABLE cette table contient les temps 700 : (automatique) de calcul (correspond a 701 : TEMPS_SAUVES) 702 : CONCENTRATIONS TABLE cette table contient les 703 : (automatique) concentrations. 704 : TEMPERATURES TABLE cette table contient les 705 : (automatique) temperatures. 706 : DEPLACEMENTS TABLE cette table contient les 707 : (automatique) deplacements. 708 : CONTRAINTES TABLE cette table contient les 709 : (automatique) contraintes. 710 : VARIABLES_INTERNES TABLE cette table contient les 711 : (automatique) variables internes en 712 : nonlineaire materiau. 713 : DEFORMATIONS_INELASTIQUES TABLE cette table contient les 714 : deformations inelastiques 715 : en nonlineaire materiau. 716 : DEFORMATIONS_APPORT TABLE cette table contient les 717 : deformations initiales de l'increment 718 : d'apport de matiere (voir DEFAP) 719 : DEFORMATIONS TABLE cette table contient les 720 : deformations (DEF) passees a COMP. 721 : DEF = DEFT - DEFTH - DEFDE - DEFIMP 722 : DEFT = deformations totale calculees 723 : selon HYPOTHESE_DEFORMATIONS 724 : DEFTH = deformations thermiques 725 : DEFDE = deformations associees a la 726 : variation de rigidite (E(T)) 727 : DEFIMP = deformations imposees 728 : MODELES TABLE Avec un modele de type CHARGEMENT, 729 : contient le modele utilise au PAS indicé. 730 : TEMPERATURES_APPORT temperature interpolee lors de l'apport 731 : (automatique si MODELE de matiere. Vaut TEMPERATURES en dehors 732 : de type CHARGEMENT) des pas d'apport de matière. 733 : VITESSES TABLE cette table contient les 734 : (automatique si dynamique) vitesses. 735 : ACCELERATIONS TABLE cette table contient les 736 : (automatique si dynamique) accelerations. 737 : REACTIONS TABLE cette table contient les 738 : (automatique) forces de reactions. 739 : REACTIONS_DIFFUSIONS TABLE cette table contient les quantites 740 : (automatique) nodales d'especes chimiques (CHPOINT) 741 : resultant des BLOCAGES_DIFFUSIONS. 742 : REACTIONS_THERMIQUES TABLE cette table contient les quantites 743 : (automatique) de chaleur nodales (CHPOINT) 744 : resultant des BLOCAGES_THERMIQUES. 745 : PROPORTIONS_PHASE TABLE cette table contient les propor- 746 : (automatique) -tions de phase si changement de 747 : phase thermique 748 : 749 : Exemple : pour lister le CHPOINT de deplacements calcules pour 750 : la valeur du parametre d'evolution 2.5, il faudra 751 : coder : 752 : LIST ( PECHE TAB1 DEPLACEMENTS 2.5 ) ; 753 : ou si on connait son numero de sauvegarde i. 754 : LIST (TAB1.DEPLACEMENTS.i); 755 : 756 : ------------------------------ REMARQUES ------------------------------- 757 : 758 : Poursuite / reprise de calcul 759 : _____________________________ 760 : 761 : 762 : l'issue d'un calcul effectue avec PASAPAS il est possible de le 763 : poursuivre ou de le reprendre en utilisant la meme table TAB1. Pour cela, 764 : il faut definir les LISTREEL des instants a calculer et a sauvegarder et 765 : mettre a jour les indices 'TEMPS_CALCULES' et 'TEMPS_SAUVES' de la table TAB1. 766 : 767 : Poursuite de calcul (sans le mot-cle 'REPRISE') 768 : +++++++++++++++++++++++++++++++++++++++++++++++ 769 : Le calcul se poursuit en repartant du dernier instant calcule, note t(i). 770 : Par defaut, ou si le mot-cle REEQUILIBRAGE est mis a FAUX, le premier 771 : pas se fera entre les instants t(i) et t(i+1). En revanche, si le 772 : mot-cle REEQUILIBRAGE est mis a VRAI, alors un reequilibrage sera effectue 773 : a l'instant t(i), puis le calcul se fera entre les instants t(i) et t(i+1). 774 : 775 : Reprise de calcul (avec le mot-cle 'REPRISE') 776 : +++++++++++++++++++++++++++++++++++++++++++++++ 777 : Le calcul reprend a partir de l'instant t(i) = TAB1.'REPRISE'. Cet instant 778 : doit faire partie de la liste des instants sauves (ou sauvegardes selon que 779 : la reprise se fasse ou non dans la meme session Cast3M). Tous les resultats 780 : posterieurs a cet instant sont supprimes de la table TAB1. 781 : Par defaut, ou si le mot-cle REEQUILIBRAGE est mis a FAUX, le premier 782 : pas se fera entre les instants t(i) et t(i+1). En revanche, si le 783 : mot-cle REEQUILIBRAGE est mis a VRAI, alors un reequilibrage sera effectue 784 : a l'instant t(i), puis le calcul se fera entre les instants t(i) et t(i+1). 785 : 786 : 787 : La continuation du calcul peut etre demandee soit dans la meme session 788 : de Cast3M (sans quitter le programme), soit dans une nouvelle session. 789 : Dans ce dernier cas, la table TAB1 doit etre prealablement sauvee avec 790 : la directive SAUV, et restituee dans la nouvelle session avec la 791 : directive REST. 792 : 793 : Etat mecanique initial 794 : ______________________ 795 : 796 : Lorsqu'on fournit un état mécanique initial (indices CONTRAINTES.0, 797 : DEPLACEMENTS.0 et, si besoin, DEFORMATIONS_INELASTIQUES.0, etc.), 798 : celui-ci ne doit pas necessairement etre équilibre (si ce n'est pas le 799 : cas, PASAPAS tiendra compte du déséquilibre initial) mais il doit etre 800 : cohérent avec la loi de comportement (par ex., contraintes linéairement 801 : dépendantes du gradient du déplacement) et avec les hypothèses de calcul 802 : (par ex., en grands déplacement, contraintes correctement transportées 803 : dans la configuration déformée). Voir aussi l'exemple "gdep4.dgibi". 804 : 805 : Procedures PERSO 806 : ________________ 807 : 808 : Une valeur de type logique peut etre positionnee dans la table (avec 809 : l'indice 'ARRET') pour arreter le calcul. Ceci permet a l'interieur 810 : des procedures PERSO1 et PERSO2 de faire sortir de la procedure 811 : PASAPAS. 812 : 813 : Il est possible dans les procedures PERSO1 et PERSO2 de changer la 814 : liste des pas de temps a calculer. Pour cela il faut donner un 815 : objet de type LISTREEL que l'on transmet a PASAPAS via la table TAB1 816 : a l'indice 'A_CALCULER'. 817 : 818 : Convergence Thermique-Mecanique 819 : _______________________________ 820 : 821 : La boucle BO_BOTH a ete mise en place dans PASAPAS pour pouvoir tenir 822 : compte des effets mecaniques sur les donnees de thermique (par exemple, 823 : dans le cas d'un transfert de chaleur dans un jeu qui se ferme). 824 : Deux nouvelles procedures ont ete ajoutees pour modifier la table de 825 : PASAPAS au cours des iterations de BO_BOTH : REEV_THE & REEV_MEC. 826 : 827 : | REEV_THE TAB1 0 828 : | REEV_MEC TAB1 0 829 : | 830 : | BOUCLE Sur les pas de temps 831 : | BOUCLE JUSQU'A NB_BOTH 832 : | | 833 : | | PREPARATION CALCUL THERMIQUE 834 : | | PROCEDURE DE CALCUL THERMIQUE : Resultats stocke dans la table 835 : | | intermediaire CHTER 836 : | | 837 : | | TAB1.'CHTER' = CHTER 838 : | | REEV_THE TAB1 1 839 : | | 840 : | | PREPARATION CALCUL MECANIQUE 841 : | | PROCEDURE DE CALCUL MECANIQUE : Resultat stocke dans la table 842 : | | intermediaire TT 843 : | | 844 : | | TAB1.'TT' = TT 845 : | | REEV_MEC TAB1 1 846 : | | 847 : | | TEST DE CONVERGENCE THERMIQUE-MECANIQUE. 848 : | | Le critere est base sur l'ecart relatif maximal entre les deux 849 : | | derniers resultats issus des procedures de resolution de pro- 850 : | | -blemes transitoires (TRANSNON, TRANSLIN ou DUPONT2) 851 : | | La precision est fixee par l'indice CRITERE_COHERENCE de la 852 : | | table d'entree de PASAPAS (1.E-2 par defaut). 853 : | | 854 : | Reactualisation de TAB1 par les contenus de CHTER et TT 855 : | 856 : 857 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales