Afficher cette notice en
1 : $$$$ PASAPAS NOTICE PV090527 26/08/21 21:15:05 12621 2 : DATE 26/08/21 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 : Peut etre necessaire en post flambage. 325 : RIGIDITE_CONSTANTE RIGIDITE matrice de rigidite constante ajoutee 326 : a celle de la structure. Par exemple, 327 : une rigidite decrivant des appuis 328 : elastiques (voir operateur APPUI). 329 : Le champ de forces nodales associe 330 : est calcule en multipliant le champ 331 : de deplacement par cette matrice. 332 : CAPACITE_CONSTANTE RIGIDITE matrice de capacite ajoutee a celle 333 : du modele thermique. 334 : Les flux thermiques nodaux associes 335 : sont calcules en multipliant le champ 336 : de temperature par cette matrice. 337 : Pour les problemes de DIFFUSION, une 338 : matrice de capacite de diffusion peut 339 : aussi etre ajoutee par cette entree. 340 : CONDUCTIVITE_CONSTANTE RIGIDITE matrice de conductivite ajoutee a 341 : celle du modele thermique. 342 : Les flux thermiques nodaux associes 343 : sont calcules en multipliant le champ 344 : de temperature par cette matrice. 345 : Pour les problemes de DIFFUSION, une 346 : matrice de diffusivite peut aussi etre 347 : ajoutee par cette entree. 348 : PILOTAGE_INDIRECT LOGIQUE VRAI si on desire activer le pilotage 349 : indirect du chargement mecanique. Il 350 : est possible de specifier differents 351 : equations de pilotage en modifiant la 352 : procedure PILOINDI. 353 : Voir l'exemple : pilotage_indirect_1.dgibi 354 : DEPLACEMENTS_PILOTES CHPOINT Pour le PILOTAGE_INDIRECT, il s'agit du 355 : CHPOINT, issu de DEPI, indiquant la 356 : direction du chargement en deplacements 357 : imposes (peut etre unitaire) 358 : FORCES_PILOTEES CHPOINT Pour le PILOTAGE_INDIRECT, il s'agit du 359 : CHPOINT de forces nodales indiquant la 360 : direction du chargement en forces 361 : imposees (peut etre unitaire) 362 : PARAMETRE_DE_PILOTAGE EVOLUTION Evolution en fonction du temps du parametre 363 : de pilotage a etre respecte via PILOINDI. 364 : ZONE_DE_PILOTAGE MAILLAGE Pour le PILOTAGE_INDIRECT, il s'agit de 365 : preciser la region ou la deformation 366 : maximale doit etre cherchee (valable pour 367 : le critere de pilotage standard de PILOINDI). 368 : Cet indice est facultatif. S'il n'est pas 369 : renseigne, l'option par default est activee 370 : (recherche sur l'ensemble du maillage). 371 : AUTOMATIQUE LOGIQUE VRAI si on desire piloter le 372 : chargement en fonction d'un 373 : critere sur le champ de 374 : deplacement. Par defaut le pas 375 : de charge est choisi pour que 376 : l'increment de deformation 377 : correspondant soit proche de la 378 : valeur TAB1.'AUTOCRIT'que l'on 379 : nomme le critere de pilotage. 380 : Il est possible de specifier un 381 : autre type de pilotage en modi 382 : fiant la procedure AUTOPILO . 383 : Dans ce cas il faut encore 384 : fournir le listreel LREEL1 mais 385 : le programme limite l'increment 386 : du temps de facon a respecter 387 : le critere sur epsilon total. 388 : AUTOCRIT FLOTTANT critere de pilotage (par defaut 389 : 0.001) 390 : AUTORECA ENTIER nombre d'iteration entre chaque recalcul du 391 : coefficient de pilotage (par defaut 1) 392 : AUTOPAS ENTIER nombre maximum de sous pas 393 : autorises pour atteindre le 394 : temps demande. 395 : AUTORESU ENTIER valeur precisant tous les combien de 396 : pas il faut garder les resultats 397 : dans le cas de calcul AUTOMATIQUE. 398 : BCSTH RIGIDITE conditions aux limites pour le 399 : tracking (methode EFEM) 400 : PAS_AJUSTE LOGIQUE permet d'imposer un choix de 401 : strategie de calcul de pas de temps. 402 : VRAI on veut faire des pas de temps 403 : ajustes en fonction de la difficulte 404 : de convergence du pas precedent. 405 : CONVERGENCE_FORCEE LOGIQUE VRAI si on force la convergence en 406 : effectuant des sous-pas non converges 407 : (le defaut) 408 : MAXDEFOR FLOTTANT increment de deformation maximum en 409 : convergence forcee (1e-3 par defaut) 410 : CONSOLIDATION LOGIQUE VRAI si calcul de consolidation 411 : DELTAITER ENTIER nombre de pas d'ecart pour le test 412 : de non convergence 413 : K_SIGMA LOGIQUE VRAI si on souhaite ajouter 414 : la matrice KSIGMA a l'operateur 415 : d'iterations. 416 : K_TANGENT LOGIQUE iterations avec la matrice tangente 417 : coherente (appel a KTAN) (FAUX par 418 : defaut) 419 : K_TANGENT_ITER0 MOT mot indiquant la matrice a utiliser 420 : pour demarrer les iterations 421 : mecaniques d'un pas 'MAT_ELASTIQUE' 422 : rigidite elastique 'MAT_TANGENTE' 423 : rigidite tangente en debut de pas 424 : (avec DT = 0.). Par defaut, on 425 : utilise la matrice calculee au pas 426 : precedent. 427 : K_TANGENT_PERT LOGIQUE iterations avec la matrice tangente 428 : calculee par perturbation (appel a 429 : KTAN PERT) (FAUX par defaut) 430 : K_TANGENT_SYME LOGIQUE La matrice tangente calculee est 431 : symetrisee (FAUX par defaut) 432 : K_TANG_PERT_C1 FLOTTANT matrice tangente par perturbation 433 : coefficient multiplicatif applique a 434 : l'increment de deformation pour 435 : determiner la perturbation 436 : (1.D-3 par defaut) 437 : K_TANG_PERT_C2 FLOTTANT matrice tangente par perturbation 438 : valeur minimale de l'increment de 439 : deformation perturbe (C1/100. par 440 : defaut) 441 : NITER_KTANGENT ENTIER la matrice tangente est recalculee 442 : toutes les NITER_KTANGENT iterations 443 : (1 par defaut) 444 : SUBSTEPPING LOGIQUE Substepping pour l'integration locale 445 : avec la matrice tangente consistante. 446 : Modeles: J2, RH_COULOMB, MRS_LADE. 447 : NMAXSUBSTEPS ENTIER nombre maximum de sous-pas pour 448 : SUBSTEPPING. 449 : NITERINTER_MAX ENTIER nombre max d'iterations au niveau 450 : local (avec SUBSTEPPING). Modeles: 451 : J2, RH_COULOMB, MRS_LADE. 452 : PRECISINTER FLOTTANT precision pour le probleme local 453 : (integration des lois constitutives) 454 : 1D-8 par defaut 455 : FEFP_FORMULATION LOGIQUE retour exponentiel avec line-search 456 : au niveau local et matrice tangente 457 : consistante. Modeles: VMT_FEFP, 458 : RHMC_FEFP, POWDER_FEFP, 459 : POWDERCAP_FEFP. 460 : UPDATE_LAGRANGIAN LOGIQUE VRAI pour le traitement de FEFP par 461 : Lagrangien Augmente (par defaut) 462 : PRECDEFO FLOTTANT Valeur de la precision du calcul, pour le critere de Cauchy 463 : sur les deformations (par default 1e-4). 464 : PRECISION FLOTTANT Valeur de la precision du calcul, 465 : quelle que soit la physique (par 466 : default 1e-4). 467 : FTOL FLOTTANT l'equilibre des efforts sera 468 : verifie a FTOL pres 469 : MTOL FLOTTANT l'equilibre des moments sera 470 : verifie a MTOL pres 471 : TTOL FLOTTANT l'equilibre en THERMIQUE sera verifie a TTOL pres (Valeur absolue) 472 : sur les variables primales 473 : CTOL FLOTTANT/CHPOINT l'equilibre en DIFFUSION sera verifie a CTOL pres (Valeur absolue) 474 : sur les variables primales 475 : PROCESSEURS MOT 'MONO_PROCESSEUR' pour imposer une 476 : execution avec un seul processeur 477 : 'COMPORTEMENT' pour imposer que seul 478 : le comportement est envoye sur 479 : plusieurs processeurs 'AUTOMATIQUE' 480 : pour imposer, l'utilisation maximale 481 : des processeurs. 482 : PREDICTEUR HPP on utilise un calcul HPP pour 483 : initialiser le calcul en grands 484 : deplacements. 485 : LINESEARCH LOGIQUE VRAI pour accéleration grands 486 : déplacements 487 : STABILITE LOGIQUE VRAI pour test stabilité 488 : ACCELERATION ENTIER on utilise l'acceleration de 489 : convergence tous les 490 : TAB1.'ACCELERATION' pas 491 : (2 par defaut). 492 : MAN LOGIQUE VRAI pour demander l'utilisation de 493 : la methode asymptotique numerique 494 : comme predicteur du deplacements en 495 : cas d'etudes en grands deplacements. 496 : Pour que cela soit vraiment interes- 497 : -sant il faut un comportement 498 : elastique lineaire des materiaux. 499 : LBC MAILLAGE support geometrique des conditions 500 : aux limites pour l'algorithme de 501 : tracking. 502 : ORDRE ENTIER prermet de preciser l'ordre du 503 : developpement de la MAN. Par defaut 504 : il vaut 12. 505 : MAXITERATION ENTIER nombre maximum d'iterations 506 : valeur prise par defaut: 49 507 : MAXISOUSPAS ENTIER nombre maximum de sous- 508 : increments en viscoplasticite ou 509 : en fluage. 510 : MAXSOUSITERATION ENTIER nombre maximum de sousiteration pour 511 : la resolution des equations d Helmholtz 512 : non lineaires utilisables en non local 513 : PRECSOUSITERATION FLOTTANT precision relative pour la convergence 514 : de l equation d HelmHoltz utilisable 515 : en non local 516 : MAXSOUSPAS ENTIER nombre maximum de sous-pas 517 : en convergence forcee 518 : valeur prise par defaut: 200 519 : NPAS_TRACKING FLOTTANT Nombre de pas a considerer pour 520 : realiser le tracking de la fissure 521 : (method EFEM) 522 : NRMAX FLOTTANT Nombre maximum d elements susceptibles 523 : de voir une fissure s initier 524 : (methode EFEM). 525 : CONN MCHAML Champ de connectivites construit par 526 : l'operateur CONN, pour tous les cas 527 : incluant des symetries.a taille de 528 : la plus grande heterogeneite). 529 : MOVA MOT Precise le nom de la variable interne 530 : sur laquelle on teste le nombre de 531 : points plastifies (EPSE par defaut) 532 : CONTRAINTES.0 MCHAML contraintes au debut du pas 533 : (0. par defaut) 534 : NB : cf. remarques. 535 : DEPLACEMENTS.0 CHPOINT deplacements au debut du pas 536 : (0. par defaut) 537 : NB : cf. remarques. 538 : VARIABLES_INTERNES.0 MCHAML variables internes au debut du pas 539 : (0. par defaut) 540 : NB : cf. remarques. 541 : DEFORMATIONS_INELASTIQUES.0 deformations inelastiques au debut 542 : MCHAML du pas (0. par defaut) 543 : NB : cf. remarques. 544 : VITESSES.0 CHPOINT vitesses au debut du pas en 545 : dynamique (0. par defaut) 546 : ACCELERATIONS.0 CHPOINT accelerations au debut du pas en 547 : dynamique (0. par defaut) 548 : TEMPS.0 FLOTTANT temps au debut du pas 549 : (0. par defaut) 550 : FREA1 CHPOINT forces de reactions au debut du pas 551 : (0. par defaut). Ce champ n'est pris 552 : en compte que lorsque un des indices 553 : DYNAMIQUE ou CONSOLIDATION de la 554 : table TAB1 est VRAI. 555 : SOUS_INCREMENT ENTIER Nombre de sous increment, pour l'integration 556 : de la loi de comportement, par pas de 557 : calcul (1 par defaut). 558 : En dynamique, si la contrainte varie 559 : fortement pendant le pas il peut etre 560 : important de suivre au mieux les 561 : contraintes afin d'adapter 562 : l'algorithme dynamique. 563 : DEPS_MAX FLOTTANT Increment maximum de deformation a partir 564 : duquel on sous decoupe le pas pour le 565 : calcul des deformations. 566 : CONCENTRATIONS.0 CHPOINT concentrations au debut du pas. 567 : (0. par defaut pour le 1er pas). 568 : TEMPERATURES.0 CHPOINT temperature au debut du pas. 569 : (0. par defaut pour le 1er pas). 570 : PROCEDURE_THERMIQUE MOT Nom de la procedure thermique a 571 : utiliser : 572 : NONLINEAIRE : nonlineaire a un pas de temps 573 : (theta-methode) 574 : LINEAIRE : lineaire (theta-methode) 575 : DUPONT : nonlineaire a deux pas de 576 : temps (methode DUPONT2) 577 : PAS_MAX FLOTTANT pas maximum autorise en thermique. 578 : 579 : CTE_STEFAN_BOLTZMANN FLOTTANT La valeur de la constante de Stefan- 580 : Boltzmann est egale par defaut a 581 : 5.673E-8 (la valeur en unites SI). Si le 582 : calcul est fait dans un autre systeme 583 : d'unites il convient de mettre ici la 584 : bonne valeur de cette constante 585 : CELSIUS LOGIQUE VRAI si le calcul de fait en degres 586 : Celsius. Lorsque cette variable est VRAI 587 : on rajoute 273 aux champs de temperature 588 : avant d'en deduire des quantites liees 589 : au rayonnement. 590 : 591 : RELAXATION_DUPONT FLOTTANT Valeur du coefficient de relaxation 592 : pour DUPONT2 (0.25 pas defaut) 593 : RELAXATION_NONCONV FLOTTANT Valeur du coefficient de relaxation 594 : pour la non convergence (1 par defaut) 595 : qui determine le sou-pas non converge 596 : RELAXATION_THETA FLOTTANT Valeur du coefficient de relaxation 597 : pour theta-methode (1.0 par defaut) 598 : SOUS_RELAXATION FLOTTANT Valeur du coefficient de sous- 599 : relaxation (1.0 par defaut) 600 : 601 : PROCEDURE_PERSO1 LOGIQUE VRAI si on doit appeler une procedure 602 : utlisateur apres le calcul mecanique. 603 : Nom de cette procedure : PERSO1. 604 : ATTENTION ! Cette procedure est appe- 605 : lee apres chaque increment, mais non 606 : au cours des iterations. 607 : Voir aussi les remarques. 608 : PROCEDURE_PERSO2 LOGIQUE VRAI si on doit appeler une procedure 609 : utlisateur apres le calcul thermique. 610 : Nom de cette procedure : PERSO2. 611 : ATTENTION ! Cette procedure est appe- 612 : lee apres chaque increment, mais non 613 : au cours des iterations. 614 : Voir aussi les remarques. 615 : 616 : PROJECTION LOGIQUE VRAI si on effectue un calcul couple 617 : et si la mecanique et la thermique ne 618 : s'appuient pas sur le meme maillage. 619 : 620 : PROCEDURE_REEV_MEC LOGIQUE VRAI si on utilise une procedure de 621 : reactualisation apres le calcul 622 : mecanique (par ex pour reevaluer les 623 : donnees du calcul thermique si elles 624 : dependent de la configuration 625 : mecanique) 626 : PROCEDURE_REEV_THE LOGIQUE VRAI si on utilise une procedure de 627 : reactualisation apres le calcul 628 : thermique (par ex pour reevaluer des 629 : donnees du calcul mecanique si elles 630 : dependent de la configuration 631 : thermique). 632 : 633 : CONVERGENCE_MEC_THE LOGIQUE VRAI si on souhaite iterer le schema 634 : thermique mecanique en cas de depen- 635 : -dance mutuellle des deux problemes. 636 : FAUX par defaut. 637 : Voir aussi les remarques. 638 : CRITERE_COHERENCE FLOTTANT critere de convergence si on itere 639 : l'alternance thermique-mecanique en 640 : cas de dependance mutuelle. Le test 641 : est effectue sur les resulats de la 642 : thermique et il porte sur l'erreur 643 : relative maximale entre les deux 644 : dernieres solutions. 645 : Vaut 'PRECISION' par defaut. 646 : NB_BOTH ENTIER Nombre d'iterations maximale de la 647 : boucle de coherence 648 : mecanique/thermique. 649 : Par defaut : 10 650 : 651 : REAC_GRANDS FLOTTANT increment de deformation pour lequel 652 : on reactualise la matrice de raideur 653 : en grands_deplacements 654 : (valeur par defaut : 10.E-2) 655 : CONVERGENCE_MONOTONE LOGIQUE VRAI si on souhaite imposer la 656 : decroissance monotone du residu au 657 : cours des iterations en recalculant 658 : la matrice de raideur et en limitant 659 : l'itere. 660 : (FAUX par defaut) 661 : INITIALISATION LOGIQUE FAUX si on ne souhaite pas que le 662 : premier residu soit initialise avec 663 : les resultats du pas precedent. 664 : (VRAI par defaut) 665 : RENORMALISATION LOGIQUE VRAI si on souhaite limiter l'itere 666 : a une valeur inferieure a MAXDEFOR. 667 : (FAUX par defaut) 668 : LAGRANGIEN MOT MI_PAS (ou REACTUALISE) (par defaut) pour appliquer le comportement 669 : sur la geometrie au milieu du pas. Meilleur ecoulement plastique. 670 : FIN_PAS pour appliquer le comportement sur la geometrie en fin 671 : du pas. Meilleures contraintes. 672 : TOTAL pour appliquer le comportement sur la geometrie initiale. 673 : Meilleur retour a zero des contraintes si on a un 674 : chargement cyclique. 675 : 676 : REPRISE FLOTTANT instant depuis lequel on souhaite reprendre le calcul 677 : 678 : REEQUILIBRAGE LOGIQUE Ce mot-cle ne sert que pour la poursuite/reprise d'un calcul. 679 : Il ne s'applique qu'a la mecanique. 680 : VRAI : si le dernier instant calcule fait partie de la liste 681 : TEMPS_CALCULES alors un pas de reequilibrage est fait. Le calcul 682 : reprend sur un pas de temps de taille nulle. 683 : 684 : En sortie, TAB1 permet de retrouver les resultats. Ceux 685 : ci sont mis dans des tables dont les indices sont des 686 : entiers (0 1 2 ... N) correspondants aux numero de 687 : sauvegarde des resultats (0 : temps initial). L'utilisateur peut 688 : intervenir de deux façons sur le contenu des resultats. Soit il utilise 689 : la table referencee dans l'indice 'MES_SAUVEGARDES' soit il surcharge 690 : directement la procedure PAS_SAUV. 691 : 692 : 693 : indice type objet commentaires 694 : pointe 695 : 696 : ERREUR LOGIQUE logique de valeur VRAI en cas 697 : d'erreur au cours de l'execution 698 : de la procedure. 699 : CONV LOGIQUE VRAI si calcul converge avant ou 700 : pour le nombre maximum d'itera- 701 : tions, FAUX sinon. 702 : TEMPS TABLE cette table contient les temps 703 : (automatique) de calcul (correspond a 704 : TEMPS_SAUVES) 705 : CONCENTRATIONS TABLE cette table contient les 706 : (automatique) concentrations. 707 : TEMPERATURES TABLE cette table contient les 708 : (automatique) temperatures. 709 : DEPLACEMENTS TABLE cette table contient les 710 : (automatique) deplacements. 711 : CONTRAINTES TABLE cette table contient les 712 : (automatique) contraintes. 713 : VARIABLES_INTERNES TABLE cette table contient les 714 : (automatique) variables internes en 715 : nonlineaire materiau. 716 : DEFORMATIONS_INELASTIQUES TABLE cette table contient les 717 : deformations inelastiques 718 : en nonlineaire materiau. 719 : DEFORMATIONS_APPORT TABLE cette table contient les 720 : deformations initiales de l'increment 721 : d'apport de matiere (voir DEFAP) 722 : DEFORMATIONS TABLE cette table contient les 723 : deformations (DEF) passees a COMP. 724 : DEF = DEFT - DEFTH - DEFDE - DEFIMP + DEFAP 725 : DEFT = deformations totale calculees 726 : selon HYPOTHESE_DEFORMATIONS 727 : DEFTH = deformations thermiques 728 : DEFDE = deformations associees a la 729 : variation de rigidite (E(T)) 730 : DEFIMP = deformations imposees 731 : DEFAP = deformations apport de matiere 732 : MODELES TABLE Avec un modele de type CHARGEMENT, 733 : contient le modele utilise au PAS indicé. 734 : TEMPERATURES_APPORT temperature initiale dans la partie ajou- 735 : (automatique si MODELE -tee du maillage lors d'un pas d'apport 736 : de type CHARGEMENT) de matiere (evolution de modele). Vaut 737 : TREF dans partie initiale du modele. 738 : VITESSES TABLE cette table contient les 739 : (automatique si dynamique) vitesses. 740 : ACCELERATIONS TABLE cette table contient les 741 : (automatique si dynamique) accelerations. 742 : REACTIONS TABLE cette table contient les 743 : (automatique) forces de reactions. 744 : REACTIONS_DIFFUSIONS TABLE cette table contient les quantites 745 : (automatique) nodales d'especes chimiques (CHPOINT) 746 : resultant des BLOCAGES_DIFFUSIONS. 747 : REACTIONS_THERMIQUES TABLE cette table contient les quantites 748 : (automatique) de chaleur nodales (CHPOINT) 749 : resultant des BLOCAGES_THERMIQUES. 750 : PROPORTIONS_PHASE TABLE cette table contient les propor- 751 : (automatique) -tions de phase si changement de 752 : phase thermique 753 : 754 : Exemple : pour lister le CHPOINT de deplacements calcules pour 755 : la valeur du parametre d'evolution 2.5, il faudra 756 : coder : 757 : LIST ( PECHE TAB1 DEPLACEMENTS 2.5 ) ; 758 : ou si on connait son numero de sauvegarde i. 759 : LIST (TAB1.DEPLACEMENTS.i); 760 : 761 : ------------------------------ REMARQUES ------------------------------- 762 : 763 : Poursuite / reprise de calcul 764 : _____________________________ 765 : 766 : 767 : l'issue d'un calcul effectue avec PASAPAS il est possible de le 768 : poursuivre ou de le reprendre en utilisant la meme table TAB1. Pour cela, 769 : il faut definir les LISTREEL des instants a calculer et a sauvegarder et 770 : mettre a jour les indices 'TEMPS_CALCULES' et 'TEMPS_SAUVES' de la table TAB1. 771 : 772 : Poursuite de calcul (sans le mot-cle 'REPRISE') 773 : +++++++++++++++++++++++++++++++++++++++++++++++ 774 : Le calcul se poursuit en repartant du dernier instant calcule, note t(i). 775 : Par defaut, ou si le mot-cle REEQUILIBRAGE est mis a FAUX, le premier 776 : pas se fera entre les instants t(i) et t(i+1). En revanche, si le 777 : mot-cle REEQUILIBRAGE est mis a VRAI, alors un reequilibrage sera effectue 778 : a l'instant t(i), puis le calcul se fera entre les instants t(i) et t(i+1). 779 : 780 : Reprise de calcul (avec le mot-cle 'REPRISE') 781 : +++++++++++++++++++++++++++++++++++++++++++++++ 782 : Le calcul reprend a partir de l'instant t(i) = TAB1.'REPRISE'. Cet instant 783 : doit faire partie de la liste des instants sauves (ou sauvegardes selon que 784 : la reprise se fasse ou non dans la meme session Cast3M). Tous les resultats 785 : posterieurs a cet instant sont supprimes de la table TAB1. 786 : Par defaut, ou si le mot-cle REEQUILIBRAGE est mis a FAUX, le premier 787 : pas se fera entre les instants t(i) et t(i+1). En revanche, si le 788 : mot-cle REEQUILIBRAGE est mis a VRAI, alors un reequilibrage sera effectue 789 : a l'instant t(i), puis le calcul se fera entre les instants t(i) et t(i+1). 790 : 791 : 792 : La continuation du calcul peut etre demandee soit dans la meme session 793 : de Cast3M (sans quitter le programme), soit dans une nouvelle session. 794 : Dans ce dernier cas, la table TAB1 doit etre prealablement sauvee avec 795 : la directive SAUV, et restituee dans la nouvelle session avec la 796 : directive REST. 797 : 798 : Etat mecanique initial 799 : ______________________ 800 : 801 : Lorsqu'on fournit un état mécanique initial (indices CONTRAINTES.0, 802 : DEPLACEMENTS.0 et, si besoin, DEFORMATIONS_INELASTIQUES.0, etc.), 803 : celui-ci ne doit pas necessairement etre équilibre (si ce n'est pas le 804 : cas, PASAPAS tiendra compte du déséquilibre initial) mais il doit etre 805 : cohérent avec la loi de comportement (par ex., contraintes linéairement 806 : dépendantes du gradient du déplacement) et avec les hypothèses de calcul 807 : (par ex., en grands déplacement, contraintes correctement transportées 808 : dans la configuration déformée). Voir aussi l'exemple "gdep4.dgibi". 809 : 810 : Procedures PERSO 811 : ________________ 812 : 813 : Une valeur de type logique peut etre positionnee dans la table (avec 814 : l'indice 'ARRET') pour arreter le calcul. Ceci permet a l'interieur 815 : des procedures PERSO1 et PERSO2 de faire sortir de la procedure 816 : PASAPAS. 817 : 818 : Il est possible dans les procedures PERSO1 et PERSO2 de changer la 819 : liste des pas de temps a calculer. Pour cela il faut donner un 820 : objet de type LISTREEL que l'on transmet a PASAPAS via la table TAB1 821 : a l'indice 'A_CALCULER'. 822 : 823 : Convergence Thermique-Mecanique 824 : _______________________________ 825 : 826 : La boucle BO_BOTH a ete mise en place dans PASAPAS pour pouvoir tenir 827 : compte des effets mecaniques sur les donnees de thermique (par exemple, 828 : dans le cas d'un transfert de chaleur dans un jeu qui se ferme). 829 : Deux nouvelles procedures ont ete ajoutees pour modifier la table de 830 : PASAPAS au cours des iterations de BO_BOTH : REEV_THE & REEV_MEC. 831 : 832 : | REEV_THE TAB1 0 833 : | REEV_MEC TAB1 0 834 : | 835 : | BOUCLE Sur les pas de temps 836 : | BOUCLE JUSQU'A NB_BOTH 837 : | | 838 : | | PREPARATION CALCUL THERMIQUE 839 : | | PROCEDURE DE CALCUL THERMIQUE : Resultats stocke dans la table 840 : | | intermediaire CHTER 841 : | | 842 : | | TAB1.'CHTER' = CHTER 843 : | | REEV_THE TAB1 1 844 : | | 845 : | | PREPARATION CALCUL MECANIQUE 846 : | | PROCEDURE DE CALCUL MECANIQUE : Resultat stocke dans la table 847 : | | intermediaire TT 848 : | | 849 : | | TAB1.'TT' = TT 850 : | | REEV_MEC TAB1 1 851 : | | 852 : | | TEST DE CONVERGENCE THERMIQUE-MECANIQUE. 853 : | | Le critere est base sur l'ecart relatif maximal entre les deux 854 : | | derniers resultats issus des procedures de resolution de pro- 855 : | | -blemes transitoires (TRANSNON, TRANSLIN ou DUPONT2) 856 : | | La precision est fixee par l'indice CRITERE_COHERENCE de la 857 : | | table d'entree de PASAPAS (1.E-2 par defaut). 858 : | | 859 : | Reactualisation de TAB1 par les contenus de CHTER et TT 860 : | 861 : 862 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales