1 : $$$$ PASAPAS NOTICE SP204843 24/02/28 21:15:01 11852 2 : DATE 24/02/28 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 PRECDECHARGE 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 : AUTOPAS ENTIER nombre maximum de sous pas 393 : autorises pour atteindre le 394 : temps demande. 395 : BCSTH RIGIDITE conditions aux limites pour le 396 : tracking (methode EFEM) 397 : AUTORESU ENTIER valeur precisant tous les combien de 398 : pas il faut garder les resultats 399 : dans le cas de calcul AUTOMATIQUE. 400 : PAS_AJUSTE LOGIQUE permet d'imposer un choix de 401 : strategie de calcul de pas de temps. 402 : VRAI on veut faire des pas de temps 403 : ajustes en fonction de la difficulte 404 : de convergence du pas precedent. 405 : CONVERGENCE_FORCEE LOGIQUE VRAI si on force la convergence en 406 : effectuant des sous-pas non converges 407 : (le defaut) 408 : MAXDEFOR FLOTTANT increment de deformation maximum en 409 : convergence forcee (1e-3 par defaut) 410 : CONSOLIDATION LOGIQUE VRAI si calcul de consolidation 411 : DELTAITER ENTIER nombre de pas d'ecart pour le test 412 : de non convergence 413 : K_SIGMA LOGIQUE VRAI si on souhaite ajouter 414 : la matrice KSIGMA a l'operateur 415 : d'iterations. 416 : K_TANGENT LOGIQUE iterations avec la matrice tangente 417 : coherente (appel a KTAN) (FAUX par 418 : defaut) 419 : K_TANGENT_ITER0 MOT mot indiquant la matrice a utiliser 420 : pour demarrer les iterations 421 : mecaniques d'un pas 'MAT_ELASTIQUE' 422 : rigidite elastique 'MAT_TANGENTE' 423 : rigidite tangente en debut de pas 424 : (avec DT = 0.). Par defaut, on 425 : utilise la matrice calculee au pas 426 : precedent. 427 : K_TANGENT_PERT LOGIQUE iterations avec la matrice tangente 428 : calculee par perturbation (appel a 429 : KTAN PERT) (FAUX par defaut) 430 : K_TANGENT_SYME LOGIQUE La matrice tangente calculee est 431 : symetrisee (FAUX par defaut) 432 : K_TANG_PERT_C1 FLOTTANT matrice tangente par perturbation 433 : coefficient multiplicatif applique a 434 : l'increment de deformation pour 435 : determiner la perturbation 436 : (1.D-3 par defaut) 437 : K_TANG_PERT_C2 FLOTTANT matrice tangente par perturbation 438 : valeur minimale de l'increment de 439 : deformation perturbe (C1/100. par 440 : defaut) 441 : NITER_KTANGENT ENTIER la matrice tangente est recalculee 442 : toutes les NITER_KTANGENT iterations 443 : (1 par defaut) 444 : SUBSTEPPING LOGIQUE Substepping pour l'integration locale 445 : avec la matrice tangente consistante. 446 : Modeles: J2, RH_COULOMB, MRS_LADE. 447 : NMAXSUBSTEPS ENTIER nombre maximum de sous-pas pour 448 : SUBSTEPPING. 449 : NITERINTER_MAX ENTIER nombre max d'iterations au niveau 450 : local (avec SUBSTEPPING). Modeles: 451 : J2, RH_COULOMB, MRS_LADE. 452 : PRECISINTER FLOTTANT precision pour le probleme local 453 : (integration des lois constitutives) 454 : 1D-8 par defaut 455 : FEFP_FORMULATION LOGIQUE retour exponentiel avec line-search 456 : au niveau local et matrice tangente 457 : consistante. Modeles: VMT_FEFP, 458 : RHMC_FEFP, POWDER_FEFP, 459 : POWDERCAP_FEFP. 460 : UPDATE_LAGRANGIAN LOGIQUE VRAI pour le traitement de FEFP par 461 : Lagrangien Augmente (par defaut) 462 : PRECISION FLOTTANT Valeur de la precision du calcul, 463 : quelle que soit la physique (par 464 : default 1e-4). 465 : PRECDECHARGE FLOTTANT valeur de la precision en decharge 466 : sous pilotage. defaut: 1.e-2 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 En dynamique et si la contrainte peut 556 : varier fortement pendant le pas il 557 : peut etre important de suivre au mieu 558 : les contraintes afin d'adapter 559 : l'algoritme dynamique. L'entier 560 : fourni est le nombre de 561 : sous_increment par pas de calcul. 562 : CONCENTRATIONS.0 CHPOINT concentrations au debut du pas. 563 : (0. par defaut pour le 1er pas). 564 : TEMPERATURES.0 CHPOINT temperature au debut du pas. 565 : (0. par defaut pour le 1er pas). 566 : PROCEDURE_THERMIQUE MOT Nom de la procedure thermique a 567 : utiliser : 568 : NONLINEAIRE : nonlineaire a un pas de temps 569 : (theta-methode) 570 : LINEAIRE : lineaire (theta-methode) 571 : DUPONT : nonlineaire a deux pas de 572 : temps (methode DUPONT2) 573 : PAS_MAX FLOTTANT pas maximum autorise en thermique. 574 : 575 : CTE_STEFAN_BOLTZMANN FLOTTANT La valeur de la constante de Stefan- 576 : Boltzmann est egale par defaut a 577 : 5.673E-8 (la valeur en unites SI). Si le 578 : calcul est fait dans un autre systeme 579 : d'unites il convient de mettre ici la 580 : bonne valeur de cette constante 581 : CELSIUS LOGIQUE VRAI si le calcul de fait en degres 582 : Celsius. Lorsque cette variable est VRAI 583 : on rajoute 273 aux champs de temperature 584 : avant d'en deduire des quantites liees 585 : au rayonnement. 586 : 587 : RELAXATION_DUPONT FLOTTANT Valeur du coefficient de relaxation 588 : pour DUPONT2 (0.25 pas defaut) 589 : RELAXATION_NONCONV FLOTTANT Valeur du coefficient de relaxation 590 : pour la non convergence (1 par defaut) 591 : qui determine le sou-pas non converge 592 : RELAXATION_THETA FLOTTANT Valeur du coefficient de relaxation 593 : pour theta-methode (1.0 par defaut) 594 : SOUS_RELAXATION FLOTTANT Valeur du coefficient de sous- 595 : relaxation (1.0 par defaut) 596 : 597 : PROCEDURE_PERSO1 LOGIQUE VRAI si on doit appeler une procedure 598 : utlisateur apres le calcul mecanique. 599 : Nom de cette procedure : PERSO1. 600 : ATTENTION ! Cette procedure est appe- 601 : lee apres chaque increment, mais non 602 : au cours des iterations. 603 : Voir aussi les remarques. 604 : PROCEDURE_PERSO2 LOGIQUE VRAI si on doit appeler une procedure 605 : utlisateur apres le calcul thermique. 606 : Nom de cette procedure : PERSO2. 607 : ATTENTION ! Cette procedure est appe- 608 : lee apres chaque increment, mais non 609 : au cours des iterations. 610 : Voir aussi les remarques. 611 : 612 : PROJECTION LOGIQUE VRAI si on effectue un calcul couple 613 : et si la mecanique et la thermique ne 614 : s'appuient pas sur le meme maillage. 615 : 616 : PROCEDURE_REEV_MEC LOGIQUE VRAI si on utilise une procedure de 617 : reactualisation apres le calcul 618 : mecanique (par ex pour reevaluer les 619 : donnees du calcul thermique si elles 620 : dependent de la configuration 621 : mecanique) 622 : PROCEDURE_REEV_THE LOGIQUE VRAI si on utilise une procedure de 623 : reactualisation apres le calcul 624 : thermique (par ex pour reevaluer des 625 : donnees du calcul mecanique si elles 626 : dependent de la configuration 627 : thermique). 628 : 629 : CONVERGENCE_MEC_THE LOGIQUE VRAI si on souhaite iterer le schema 630 : thermique mecanique en cas de depen- 631 : -dance mutuellle des deux problemes. 632 : FAUX par defaut. 633 : Voir aussi les remarques. 634 : CRITERE_COHERENCE FLOTTANT critere de convergence si on itere 635 : l'alternance thermique-mecanique en 636 : cas de dependance mutuelle. Le test 637 : est effectue sur les resulats de la 638 : thermique et il porte sur l'erreur 639 : relative maximale entre les deux 640 : dernieres solutions. 641 : Vaut 'PRECISION' par defaut. 642 : NB_BOTH ENTIER Nombre d'iterations maximale de la 643 : boucle de coherence 644 : mecanique/thermique. 645 : Par defaut : 10 646 : 647 : REAC_GRANDS FLOTTANT increment de deformation pour lequel 648 : on reactualise la matrice de raideur 649 : en grands_deplacements 650 : (valeur par defaut : 10.E-2) 651 : CONVERGENCE_MONOTONE LOGIQUE VRAI si on souhaite imposer la 652 : decroissance monotone du residu au 653 : cours des iterations en recalculant 654 : la matrice de raideur et en limitant 655 : l'itere. 656 : (FAUX par defaut) 657 : INITIALISATION LOGIQUE FAUX si on ne souhaite pas que le 658 : premier residu soit initialise avec 659 : les resultats du pas precedent. 660 : (VRAI par defaut) 661 : RENORMALISATION LOGIQUE VRAI si on souhaite limiter l'itere 662 : a une valeur inferieure a MAXDEFOR. 663 : (FAUX par defaut) 664 : LAGRANGIEN MOT MI_PAS (ou REACTUALISE) (par defaut) pour appliquer le comportement 665 : sur la geometrie au milieu du pas. Meilleur ecoulement plastique. 666 : FIN_PAS pour appliquer le comportement sur la geometrie en fin 667 : du pas. Meilleures contraintes. 668 : TOTAL pour appliquer le comportement sur la geometrie initiale. 669 : Meilleur retour a zero des contraintes si on a un 670 : chargement cyclique. 671 : 672 : REPRISE FLOTTANT instant depuis lequel on souhaite reprendre le calcul 673 : 674 : REEQUILIBRAGE LOGIQUE Ce mot-cle ne sert que pour la poursuite/reprise d'un calcul. 675 : Il ne s'applique qu'a la mecanique. 676 : VRAI : si le dernier instant calcule fait partie de la liste 677 : TEMPS_CALCULES alors un pas de reequilibrage est fait. Le calcul 678 : reprend sur un pas de temps de taille nulle. 679 : 680 : En sortie, TAB1 permet de retrouver les resultats. Ceux 681 : ci sont mis dans des tables dont les indices sont des 682 : entiers (0 1 2 ... N) correspondants aux numero de 683 : sauvegarde des resultats (0 : temps initial). L'utilisateur peut 684 : intervenir de deux façons sur le contenu des resultats. Soit il utilise 685 : la table referencee dans l'indice 'MES_SAUVEGARDES' soit il surcharge 686 : directement la procedure PAS_SAUV. 687 : 688 : 689 : indice type objet commentaires 690 : pointe 691 : 692 : ERREUR LOGIQUE logique de valeur VRAI en cas 693 : d'erreur au cours de l'execution 694 : de la procedure. 695 : CONV LOGIQUE VRAI si calcul converge avant ou 696 : pour le nombre maximum d'itera- 697 : tions, FAUX sinon. 698 : TEMPS TABLE cette table contient les temps 699 : (automatique) de calcul (correspond a 700 : TEMPS_SAUVES) 701 : CONCENTRATIONS TABLE cette table contient les 702 : (automatique) concentrations. 703 : TEMPERATURES TABLE cette table contient les 704 : (automatique) temperatures. 705 : DEPLACEMENTS TABLE cette table contient les 706 : (automatique) deplacements. 707 : CONTRAINTES TABLE cette table contient les 708 : (automatique) contraintes. 709 : VARIABLES_INTERNES TABLE cette table contient les 710 : (automatique) variables internes en 711 : nonlineaire materiau. 712 : DEFORMATIONS_INELASTIQUES TABLE cette table contient les 713 : deformations inelastiques 714 : en nonlineaire materiau. 715 : DEFORMATIONS_APPORT TABLE cette table contient les 716 : deformations initiales de l'increment 717 : d'apport de matiere (voir DEFAP) 718 : DEFORMATIONS TABLE cette table contient les 719 : deformations (DEF) passees a COMP. 720 : DEF = DEFT - DEFTH - DEFDE - DEFIMP 721 : DEFT = deformations totale calculees 722 : selon HYPOTHESE_DEFORMATIONS 723 : DEFTH = deformations thermiques 724 : DEFDE = deformations associees a la 725 : variation de rigidite (E(T)) 726 : DEFIMP = deformations imposees 727 : MODELES TABLE Avec un modele de type CHARGEMENT, 728 : contient le modele utilise au PAS indicé. 729 : TEMPERATURES_APPORT temperature interpolee lors de l'apport 730 : (automatique si MODELE de matiere. Vaut TEMPERATURES en dehors 731 : de type CHARGEMENT) des pas d'apport de matière. 732 : VITESSES TABLE cette table contient les 733 : (automatique si dynamique) vitesses. 734 : ACCELERATIONS TABLE cette table contient les 735 : (automatique si dynamique) accelerations. 736 : REACTIONS TABLE cette table contient les 737 : (automatique) forces de reactions. 738 : REACTIONS_DIFFUSIONS TABLE cette table contient les quantites 739 : (automatique) nodales d'especes chimiques (CHPOINT) 740 : resultant des BLOCAGES_DIFFUSIONS. 741 : REACTIONS_THERMIQUES TABLE cette table contient les quantites 742 : (automatique) de chaleur nodales (CHPOINT) 743 : resultant des BLOCAGES_THERMIQUES. 744 : PROPORTIONS_PHASE TABLE cette table contient les propor- 745 : (automatique) -tions de phase si changement de 746 : phase thermique 747 : 748 : Exemple : pour lister le CHPOINT de deplacements calcules pour 749 : la valeur du parametre d'evolution 2.5, il faudra 750 : coder : 751 : LIST ( PECHE TAB1 DEPLACEMENTS 2.5 ) ; 752 : ou si on connait son numero de sauvegarde i. 753 : LIST (TAB1.DEPLACEMENTS.i); 754 : 755 : ------------------------------ REMARQUES ------------------------------- 756 : 757 : Poursuite / reprise de calcul 758 : _____________________________ 759 : 760 : 761 : l'issue d'un calcul effectue avec PASAPAS il est possible de le 762 : poursuivre ou de le reprendre en utilisant la meme table TAB1. Pour cela, 763 : il faut definir les LISTREEL des instants a calculer et a sauvegarder et 764 : mettre a jour les indices 'TEMPS_CALCULES' et 'TEMPS_SAUVES' de la table TAB1. 765 : 766 : Poursuite de calcul (sans le mot-cle 'REPRISE') 767 : +++++++++++++++++++++++++++++++++++++++++++++++ 768 : Le calcul se poursuit en repartant du dernier instant calcule, note t(i). 769 : Par defaut, ou si le mot-cle REEQUILIBRAGE est mis a FAUX, le premier 770 : pas se fera entre les instants t(i) et t(i+1). En revanche, si le 771 : mot-cle REEQUILIBRAGE est mis a VRAI, alors un reequilibrage sera effectue 772 : a l'instant t(i), puis le calcul se fera entre les instants t(i) et t(i+1). 773 : 774 : Reprise de calcul (avec le mot-cle 'REPRISE') 775 : +++++++++++++++++++++++++++++++++++++++++++++++ 776 : Le calcul reprend a partir de l'instant t(i) = TAB1.'REPRISE'. Cet instant 777 : doit faire partie de la liste des instants sauves (ou sauvegardes selon que 778 : la reprise se fasse ou non dans la meme session Cast3M). Tous les resultats 779 : posterieurs a cet instant sont supprimes de la table TAB1. 780 : Par defaut, ou si le mot-cle REEQUILIBRAGE est mis a FAUX, le premier 781 : pas se fera entre les instants t(i) et t(i+1). En revanche, si le 782 : mot-cle REEQUILIBRAGE est mis a VRAI, alors un reequilibrage sera effectue 783 : a l'instant t(i), puis le calcul se fera entre les instants t(i) et t(i+1). 784 : 785 : 786 : La continuation du calcul peut etre demandee soit dans la meme session 787 : de Cast3M (sans quitter le programme), soit dans une nouvelle session. 788 : Dans ce dernier cas, la table TAB1 doit etre prealablement sauvee avec 789 : la directive SAUV, et restituee dans la nouvelle session avec la 790 : directive REST. 791 : 792 : Etat mecanique initial 793 : ______________________ 794 : 795 : Lorsqu'on fournit un état mécanique initial (indices CONTRAINTES.0, 796 : DEPLACEMENTS.0 et, si besoin, DEFORMATIONS_INELASTIQUES.0, etc.), 797 : celui-ci ne doit pas necessairement etre équilibre (si ce n'est pas le 798 : cas, PASAPAS tiendra compte du déséquilibre initial) mais il doit etre 799 : cohérent avec la loi de comportement (par ex., contraintes linéairement 800 : dépendantes du gradient du déplacement) et avec les hypothèses de calcul 801 : (par ex., en grands déplacement, contraintes correctement transportées 802 : dans la configuration déformée). Voir aussi l'exemple "gdep4.dgibi". 803 : 804 : Procedures PERSO 805 : ________________ 806 : 807 : Une valeur de type logique peut etre positionnee dans la table (avec 808 : l'indice 'ARRET') pour arreter le calcul. Ceci permet a l'interieur 809 : des procedures PERSO1 et PERSO2 de faire sortir de la procedure 810 : PASAPAS. 811 : 812 : Il est possible dans les procedures PERSO1 et PERSO2 de changer la 813 : liste des pas de temps a calculer. Pour cela il faut donner un 814 : objet de type LISTREEL que l'on transmet a PASAPAS via la table TAB1 815 : a l'indice 'A_CALCULER'. 816 : 817 : Convergence Thermique-Mecanique 818 : _______________________________ 819 : 820 : La boucle BO_BOTH a ete mise en place dans PASAPAS pour pouvoir tenir 821 : compte des effets mecaniques sur les donnees de thermique (par exemple, 822 : dans le cas d'un transfert de chaleur dans un jeu qui se ferme). 823 : Deux nouvelles procedures ont ete ajoutees pour modifier la table de 824 : PASAPAS au cours des iterations de BO_BOTH : REEV_THE & REEV_MEC. 825 : 826 : | REEV_THE TAB1 0 827 : | REEV_MEC TAB1 0 828 : | 829 : | BOUCLE Sur les pas de temps 830 : | BOUCLE JUSQU'A NB_BOTH 831 : | | 832 : | | PREPARATION CALCUL THERMIQUE 833 : | | PROCEDURE DE CALCUL THERMIQUE : Resultats stocke dans la table 834 : | | intermediaire CHTER 835 : | | 836 : | | TAB1.'CHTER' = CHTER 837 : | | REEV_THE TAB1 1 838 : | | 839 : | | PREPARATION CALCUL MECANIQUE 840 : | | PROCEDURE DE CALCUL MECANIQUE : Resultat stocke dans la table 841 : | | intermediaire TT 842 : | | 843 : | | TAB1.'TT' = TT 844 : | | REEV_MEC TAB1 1 845 : | | 846 : | | TEST DE CONVERGENCE THERMIQUE-MECANIQUE. 847 : | | Le critere est base sur l'ecart relatif maximal entre les deux 848 : | | derniers resultats issus des procedures de resolution de pro- 849 : | | -blemes transitoires (TRANSNON, TRANSLIN ou DUPONT2) 850 : | | La precision est fixee par l'indice CRITERE_COHERENCE de la 851 : | | table d'entree de PASAPAS (1.E-2 par defaut). 852 : | | 853 : | Reactualisation de TAB1 par les contenus de CHTER et TT 854 : | 855 : 856 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales