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