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