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