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