$$$$ PASAPAS NOTICE PASCAL 14/01/27 21:15:00 7920 DATE 14/01/27 Procedure PASAPAS Voir aussi : PECHE _________________ EXPLORER CHARTHER PAS_SAUV TAB1. ACCELERATIONS MES_SAUVEGARDES AMORTISSEMENT MODELE AUTOCRIT MOVA AUTORESU MTOL AUTOMATIQUE NB_BOTH AUTOPAS NITERINTER_MAX BLOCAGES_DIFFUSIONS NITER_KTANGENT BLOCAGES_MECANIQUES NMAXSUBSTEPS BLOCAGES_THERMIQUES NON_LOCAL CONCENTRATIONS ORDRE CONN PAS_AJUSTE CELSIUS PAS_MAX CARACTERISTIQUES PRECDECHARGE CHARGEMENT PRECISINTER CONSOLIDATION PRECISION CONTRAINTES PREDICTEUR CONVERGENCE_FORCEE PROCEDURE_CHAR_MECA CONVERGENCE_MEC_THE PROCEDURE_CHAR_THER CONVERGENCE_MONOTONE PROCEDURE_PARA_THER CRITERE_COHERENCE PROCEDURE_PERSO1 CTE_STEFAN_BOLTZMANN PROCEDURE_PERSO2 DEFORMATIONS_INELASTIQUES PROCEDURE_REAC_MECA DELTAITER PROCEDURE_REAC_THER DEPLACEMENTS PROCEDURE_THERMIQUE DYNAMIQUE PROCESSEURS ECONOMIQUE PROJECTION FEFP_FORMULATION REACTIONS FREA1 REAC_GRANDS FTOL RELAXATION_DUPONT GRANDS_DEPLACEMENTS RELAXATION_NONCONV INITIALISATION RELAXATION_THETA K_SIGMA RENORMALISATION K_TANGENT RIGIDITE_CONSTANTE K_TANGENT_ITER0 SOUS_INCREMENT K_TANGENT_PERT SOUS_RELAXATION K_TANGENT_SYME STABILITE K_TANG_PERT_C1 SUBSTEPPING K_TANG_PERT_C2 TALPHA_REFERENCE LCAR TEMPERATURES LINESEARCH TEMPS LVIAM TEMPS_CALCULES MAN TEMPS_SAUVEGARDES MASSE_CONSTANTE TEMPS_SAUVES MAXDEFOR UPDATE_LAGRANGIAN MAXSOUSPAS VARIABLES_INTERNES MAXISOUSPAS VITESSES MAXITERATION Description : _____________ MECHANIC : This procedure enables the user to carry out an incremental non linear calculation. The non linearity can result either from the material (plasticity) or from the large displacements, or else from both. When the option MODE FREQ is adopted, the procedure solves the dynamic equation according to the spectral method. Thus time is implicitly read as frequency for the CHARGEMENT object with no change in the name of component. The procedure computes following a default frequency list when the TEMPS_CALCULES (computed times) list is not given. THERMIC : This procedure enables the user to carry out a linear and non-linear calculation with conduction, convection and radiation. DIFFUSION This procedure computes the solution of a linear or non-linear problem of diffusion. This procedure enables the user to carry out a coupled THERMAL- MECHANICAL-DIFFUSION calculations. The results are calculated for some values of the evolution parameter (pseudo time or real time) defined by the user. Contents : __________ TAB1 : TABLE type object On input, TAB1 is used to define the options and the calculation parameters. The indices of the TAB1 object are words (written in full words, and in capital letters if they are placed between quotes). Here is their list : index pointed object comments type BLOCAGES_DIFFUSIONS RIGIDITE mortar matrices for diffusion BLOCAGES_MECANIQUES RIGIDITE mortar matrices for mechanics BLOCAGES_THERMIQUES RIGIDITE mortar matrices for thermics CARACTERISTIQUES MCHAML field of material and geometrical (when necessary) characteristics (CARACTERISTIQUES subtype) Its components may be the following types : 1) FLOTTANT if the component is constant through whole structure; 2) MCHAML if the component is independent of any parameters, but depends upon the position of a point; 3) EVOLUTION if the component varies in function of one parameter field; 4) NUAGE if the component is described by a curve of EVOLUTION type depending upon on a parameter field. CHARGEMENT CHARGEME definition of the loading in relation to the evolution parameter. (created by the CHAR operator). The elementary loading are called : - temperature T - imposed concentrations CIMP - imposed displacements DIMP - imposed temperatures TIMP - others loadings in mecanics MECA - flux (consolidation) FLUX - diffusion flux QC0 - thermic flux Q - external tempe convection TECO - external temperat radiation TERA - imposed deformations DEFI - external variables ? The external variables are these which are necessary for the evaluation of the material field (recommended type = MCHAML). They may be the external parameters of the model, if such parameters exist. The name of the component of the "diffusion flux" is named QCO by default. However, users can rename it (see the manual page of the MODE operator). In this case, the chosen name must be used insteadf of QCO. In case of loads that must be evaluated during the equilibrium process (By example following pressure or radiation) You have to take advantage of possibilty offered by PROCEDURE_CHAR_MECA or PROCEDURE_CHAR_THER. If the calculation of some "external variables" are required at each thermal iteration, in order to take into account their evolution on the material field, the use of the PARATHER procedure will be helpful. PROCEDURE_CHAR_MECA : Logical. VRAI(true) If there are loads to be evaluated during equilibrium iteration process. This computation will be done by invoking the CHARMECA procedure. The real call is : TAB2 = CHARMECA TAB TIME; TAB is the input TABLE in PASAPAS and TIME is the time for which loads must be computed. TAB2.'ADDI_SECOND' is the computed loads. Loads must be evaluated on actual configuration. Examples gdep2 and gdep3 show how to introduce following pressure. PROCEDURE_CHAR_THER : Logical. VRAI(true) If there are loads to be evaluated during iteration process to find temperature. This computation will be done by invoking the CHARTHER procedure. The real call is : TAB2 = CHARTHER TAB TIME; TAB is the TABLE object input in PASAPAS and TIME is the time for which loads must be computed. For old syntaxes for radiative flux, the index PROCEDURE_CHARG_THER is set to VRAI (TRUE) by the PASAPAS procedure. The CHARTHER procedure given with the code is an example of use in case of radiative flux using old syntax.(see example rayo-2D-4 and rayo-2D-4-bis) PROCEDURE_PARA_THER Logical. VRAI(true) if some "external variables" need to be evaluated at each thermal process iteration. This update will be released in the PARATHER procedure whose call will be "PARATHER TAB1 INST ;". TAB1 is the input TABLE object for PASAPAS procedure and INST is the time for which the "extranal variables" need to be calculated for the thermial problem. (An example for the PARATHER procedure using is available in the "exemple_parather.dgibi" file.) TEMPS_CALCULES LISTREEL definition of the values of the evolution parameter (or time) for which results are requested. If this data is missing, the code will adopt an adaptative time step and will consider the values given for results recording ( index TEMPS_SAUVES ). see also PAS_AJUSTE. TEMPS_SAUVES LISTREEL definition of the values of the evolution parameter (or time) for which results are recorded. TEMPS_SAUVEGARDES LISTREEL definition of the values of the evolution parameter (or time) for which PASAPAS will invoke the SAUV operator (SAUV TAB1). The values must be part of the values of TEMPS_SAUVES. The file where to keep result must be defined (by OPTION SAUV '.... ';) before invoking PASAPAS. MES_SAUVEGARDES TABLE allow to ask for storing results in addition of those automaticlly stored. Authorized indexes are : -'DEFIN'=VRAI for inelastic strains -'DEFTO'=VRAI for total strains (in corotationnnal axes if 'GRANDES_DEFORMATIONS' of Jaumann's type) -'DEFLO'=VRAI for logaritmic strains if 'GRANDES_DEFORMATIONS' type 'UTILISATEUR' -'ROTAF'=VRAI for rotations of the corotationnal axes in case of 'GRANDES_DEFORMATIONS' of Jaumann's type One can use this table to put informations to be used in the PAS_SAUV procedure. This procedure put results in results' tables and can be overwritten by the user.( see the manual of PAS_SAUV procedure) ECONOMIQUE LOGIQUE ask the code to clean out the memory space with already computed results by the FANTOM operator. MODELE MMODEL model object including the whole stucture excepted the model object corresponding to RIGIDITE_CONSTANTE. GRANDS_DEPLACEMENTS LOGIQUE VRAI (true) if the geometry must be recomputed at each step in the case of large displacements. This option set the two logical values GRANDES_ROTATIONS and K_SIGMA to TRUE. To cancel these options declare the corresponding index to FALSE. DYNAMIQUE LOGIQUE VRAI (true) if dynamic analysis MASSE_CONSTANTE RIGIDITE Constant mass matrix (for example lumped mass matrix). Consistent mass matrix is defined by the given model and material. AMORTISSEMENT RIGIDITE damping matrix optional for dynamic calculation RIGIDITE_CONSTANTE RIGIDITE stiffness of the material being always linear. The associed nodal force is obtained by the displacement multiplied by this constante stiffness matrice. TEMPERATURE_REFERENCE CHPOINT gives the reference field of temperature, i.e. the temperature field which does not create stresses. By default, this temperature chart is considered to equal that of the initial time.If the reference temperature chart (Tref) and that of the initial time (T0) are different, therefore the deformations will be null at the initial time but there will be initial stresses equivalent to : YOUNG*ALPH*( T0 - Tref) AUTOMATIQUE LOGIQUE VRAI (TRUE) if one wants to limit the load step according to a certain criterion on the displacement field. By default the load step is chosen in order that the maximum strain increment be about TAB1.'AUTOCRIT'. It is possible to specify another criterion by changing the content of the procedure AUTOPILO. In this case, the listreel LREEL1 will have to be supplied once again but the program limits the time increment so as to respect the criterion on total epsilon. AUTOCRIT FLOTTANT value used for the criterion above (0.001 by default) AUTOPAS ENTIER maximum number of possible substeps for reaching the requested time. AUTORESU ENTIER value specifying how often the results must be stored in case of AUTOMATIQUE calculation. PAS_AJUSTE LOGIQUE VRAI (true) if an automatic adaptative time step must be computed function of difficulty of converge for the previous time step. CONVERGENCE_FORCEE LOGIQUE VRAI (true) if convergence is forced by relaxation (the default) MAXDEFOR FLOTTANT maximum strain increment in CONVERGENCE_FORCEE (default is 1e-3) DELTAITER ENTIER nomber of steps for the nonconvergence test K_SIGMA LOGIQUE VRAI (true) if the matrix KSIGMA is to be added to the iteration operator ( the stiffness matrix by default ) K_TANGENT LOGIQUE iterations using the consistent stiffness matrix (call to the KTAN operator) (FAUX by default) K_TANGENT_ITER0 MOT word to indicate the tangential stiffness matrix to start the first mechanical iterations at the beginning of a time step 'MAT_ELASTIQUE' elastic stiffness matrix 'MAT_TANGENTE' tangential stiffness matrix at the beginning of a step (with DT = 0.). By default, the last calculated matrix is used. K_TANGENT_PERT LOGIQUE iterations using the tangential stiffness matrix calculated with the disturbing method (call to the KTAN operator with the keyword PERT) (by default FAUX) K_TANGENT_SYME LOGIQUE The used tangential stiffness matrix is symetrized (by default FAUX) K_TANG_C1 FLOTTANT tangential sitffness matrix by disturbance multiplicative factor of the strain increment to calculate the disturbance (by default 1.D-3) K_TANG_C2 FLOTTANT tangential sitffness matrix by disturbance minimum of the disturbed increment strain (by default C1/100.) NITER_KTANGENT ENTIER The k_tangent is recomputed every NITER_KTANGENT iterations (by default 1). SUBSTEPPING LOGIQUE Substepping local integration with consistent tangent matrix. Models: J2, RH_COULOMB, MRS_LADE. NMAXSUBSTEPS ENTIER max. number of substeps for SUBSTEPPING. NITERINTER_MAX ENTIER max. number of iterations at local level (with SUBSTEPPING). Models: J2, RH_COULOMB, MRS_LADE. PRECISINTER FLOTTANT precision of local problems (constitutive laws integration). Default value 1.D-8. FEFP_FORMULATION LOGIQUE exponential return mapping with line-search at local level and cons. tangent matrix. Models: VMT_FEFP, RHMC_FEFP, POWDER_FEFP, POWDERCAP_FEFP. UPDATE_LAGRANGIAN LOGIQUE VRAI (true) for UL treatment of FeFp formulations (the default). PRECISION FLOTTANT Computational precision, for all models (1.E-4 by default). PRECDECHARGE FLOTTANT value of the specification for discharge under AUTOMATIQUE option defaut: 1.e-2 FTOL FLOTTANT the effort equilibrium will be checked to within about FTOL MTOL FLOTTANT the moment equilibrium will be checked to within about MTOL PROCESSEURS MOT 'MONO_PROCESSEUR' to impose a run with only one processor 'COMPORTEMENT' to impose that only the material behavior will use all processors 'AUTOMATIQUE' to ask the code to use as much as possible processors PREDICTEUR HPP Initialise the great displacement schema with a small pertubation analysis. LINESEARCH LOGIQUE VRAI for large displacements acceleration STABILITE LOGIQUE VRAI for stability control ACCELERATION ENTIER We use the convergence acceleration every TAB1.'ACCELERATION' steps (by default 2). MAN LOGICAL VRAI (true) ifin order to use the numerical asymptotic method to compute the prediction of displacements in case of GRANDS_DEPLACEMENTS and GRANDES_ROTATIONS. The interest of this method is only when behaviour of materails is linear elastic. ORDRE ENTIER Allows to define the order of expansion of the MAN. (default value 12) MAXITERATION ENTIER maximum number of iterations default value : 49 MAXISOUSPAS ENTIER maximum number of sub- increments in viscoplasticity or creep default value : 200 MAXSOUSPAS ENTIER maximum number of substeps in forced convergence NON_LOCAL MOT Used in case of non local caclulation and allows to precise the method to be used: 'MOYE' : by default method Then LCAR, CONN and LVIAM indexes will have to be provided. 'SB ' : non local stress-based method Then LCAR, CONN, LVIAM, NLOC_SB_MODL and NLOC_SB_FT indexes will have to be provided LCAR FLOTTANT feature length to be supplied in case of non local calculation (it will be possible to take 3 times the size of the greatest heterogeneity) CONN MCHAML field of connectivities built from the CONN operator, for all the cases including symmetries.a, size of the greatest heterogeneity) LVIAM LISTMOTS names of component to be averaged in non local. NLOC_SB_MODL MMODEL model considered for the 'SB ' regularization method NLOC_SB_FT FLOTTANT value/field corresponding to the limit stress MCHAML in traction used to normalize fields in the 'SB ' regularization method MOVA MOT specifies the name of the internal variable on which the number of plastic points are tested (EPSE by default) CONTRAINTES.0 MCHAML stresses at beginning of step (0. by default) NB : see remarks below. DEPLACEMENTS.0 CHPOINT displacements at beginning of step (0. by default) NB : see remarks below. VARIABLES_INTERNES.0 MCHAML internal variables at beginning of step (0. by default) NB : see remarks below. DEFORMATIONS_INELASTIQUES.0 MCHAML inelastic strains at beginning of step (0. by default) NB : see remarks below. VITESSES.0 CHPOINT speeds at beginning of step (0. by default) ACCELERATIONS.0 CHPOINT accelerations at beginning of step (0. by default) TEMPS.0 FLOTTANT time at beginning of step (0. by default) FREA1 CHPOINT reaction forces at beginning of step (0. by default). This field is taken into account only when one of DYNAMIQUE or CONSOLIDATION indices is equal to VRAI. SOUS_INCREMENT ENTIER For dynamic and non linear material behaviour it maybe important to follow as well as possible srtresses values. The integer value indicates how many substeppings will be performed for each step of calculation. CONCENTRATIONS.0 CHPOINT concentration at the beginning of the step (by default 0. at the 1st step). TEMPERATURES.0 CHPOINT temperature at the beginning of the step (by default 0. at the 1st step). PROCEDURE_THERMIQUE ENTIER Number of the thermic procedure to use NONLINEAIRE : nonlinear , one step of time theta-method LINEAIRE : linear theta-method DUPONT : nonlinear, two steps of time DUPONT2 method PAS_MAX FLOTTANT maximum step authorized in thermics CTE_STEFAN_BOLTZMANN FLOTTANT The value of the Stefan-Boltzmann constant is equal by default to 5.673E-8 (SI units). If another unit system is used, the correct value should be put here. CELSIUS LOGIQUE VRAI (true) if the temperatures are in Celsius degrees. When it's true, 273.15 is added to the temperature fields before the radiation coefficients are computed. RELAXATION_DUPONT FLOTTANT value of the relaxation coefficient for DUPONT2 method (0.25 by default) RELAXATION_NONCONV FLOTTANT value of the relaxation coefficient for forced convergence (1 by default). It defines the unconverged substep RELAXATION_THETA FLOTTANT value of the relaxation coefficient for theta-method (1.0 by default) SOUS_RELAXATION FLOTTANT value of the sub-relaxation coefficient (1.0 by default) PROCEDURE_PERSO1 LOGIQUE VRAI if a personal procedure should be called after each mecanic step; name of this procedure : PERSO1. WARNING ! This procedure will be called after each step, but not during the iterations. See also the remarks. PROCEDURE_PERSO2 LOGIQUE VRAI if a personal procedure should be called after each thermic step; name of this procedure : PERSO2. WARNING ! This procedure will be called after each step, but not during the iterations. See also the remarks. PROJECTION LOGIQUE VRAI if coupled calculation and if mecanic and thermic meshes are different PROCEDURE_REAC_MECA LOGIQUE VRAI if a procedure is used to update data after the mecanichal resolution. (for instance to evaluate the data for the thermic resolution if it depends on the mecanical configuration) PROCEDURE_REAC_THER LOGIQUE VRAI if an update procedure is used after the thermic resolution (for instance to update the mechanical data depending on the temperature -Internal pressure of a container-) CONVERGENCE_MEC_THE LOGIQUE VRAI if one wants to activate the thermic-mecanic loop (see remarks below) in case of mutual dependance of the two problems. FAUX by default. CRITERE_COHERENCE FLOTTANT convergence criterion if the thermic- mechanic loop is activated. The convergence test is based on the relative error of the thermic results. By default, equal to the 'PRECISION'. NB_BOTH ENTIER Maximum number of time the thermic/mecanic loop is executed. Default : 10 REAC_GRANDS ENTIER define at which iteration the stiffness matrix is recomputed on the new geometry in case of "grands_deplacements" (default value : 2) CONVERGENCE_MONOTONE LOGIQUE VRAI if residue is imposed to monotonically decrease during the iterations. This is achieve by recomputing the stiffness and limiting the iterate. (FAUX by default) INITIALISATION LOGIQUE FAUX to prevent the first residue to be initialized with the result of the previous step. (VRAI by default) RENORMALISATION LOGIQUE VRAI to limite the iterate to a linear strain less than MAXDEFOR. (FAUX by default) On output, TAB1 enables the user to find again the results which are stored in tables whose indices are intergers (0 1 2 .... N) corresponding to the number of the results saving. 0 defined the initial time. index pointed object comments type ERREUR LOGIQUE logical result of value VRAI (TRUE) in case of error during the procedure execution. CONV LOGIQUE VRAI (TRUE) if calculation converged before or for the maximum number of iterations, FAUX (FALSE) if not. DEPLACEMENTS TABLE this table contains the computed displacements CONTRAINTES TABLE this table contains the computed stresses VARIABLES_INTERNES TABLE this table contains the internal variables in non linear material DEFORMATIONS_INELASTIQUES TABLE this table contains the inelastic strains in non linear material. DEFORMATIONS TABLE this table contains the strains cumulated in the corotationnal referential in case of large strains with Jaumann's derivative ROTATIONS TABLE this table contains the rotations of the corotationnal referential in case of large strains with Jaumann's derivative REACTIONS TABLE this table contains the computed reaction forces VITESSES TABLE this table contains the computed velocities ACCELERATIONS TABLE this table contains the computed accelerations CONCENTRATIONS TABLE this table contains the computed concentrations TEMPERATURES TABLE this table contains the computed temperatures PROPORTIONS_PHASE TABLE contains proportions of phases in case of thermal phase change. Example : in order to list the CHPOINT of displacements calculated for the value of the evolution parameter 2.5, the following will have to be coded : LIST ( PECHE TAB1 DEPLACEMENTS 2.5 ) ; ------------------------------- REMARKS -------------------------------- Continuation / Restart with PASAPAS ___________________________________ It is possible to continue (restart) a calculus carried out with the PASAPAS procedure by invoking one more time the command : PASAPAS TAB1 ; Before issuing this command, one must redefine the LISTREEL objects pointed on the 'TEMPS_CALCULES' and 'TEMPS_SAUVES' index. If (one of) those new LISTREEL objects contain(s) the values that are less than TAB1.'TEMPS'.NUMPAS, these values will be ignored. The restart of the calculus can be asked in the same session (without exiting the program) or in an other one. In the latter case, the table TAB1 should have been output with the SAUV directive, and reloaded with the REST directive. Initial mechanical state ________________________ When one gives an initial mechanical state (index CONTRAINTES.0, DEPLACEMENTS.0 and, possibly, DEFORMATIONS_INELASTIQUES.0, etc.), it does not have to be equilibrated (if not, PASAPAS will take into account this initial desiquilibrium) but it must be consistant with the material behavior law (for instance, the stress depends linearly on the gradient of the displacements) or with the hypothesis of the calculus, e.g. within the large displacement option, the stress field has to be correctly transported onto the deformed configuration. See also the example "gdep4.dgibi". Procedures PERSO ________________ Within procedures PERSO1 or PERSO2 it is possible to quit the procedure PASAPAS. To do so the table TAB1 at index 'ARRET' (stop) must be set to VRAI (true) ( TAB1.'ARRET' = VRAI;). It is possible to change the list of times step to be computed by seting an object of type LISTREEL giving the new times which must be computed. This object must be parse to the procedure PASAPAS by the table TAB1 at index 'A_CALCULER' (to_compute). Convergence THERMIC-MECHANICS _____________________________ The loop BO_BOTH has been programmed in PASAPAS to take into account the feedback of mechanics on the thermal calculation data (for instance in case of heat transfer in a gap that is progressively closed). Two news procedures have been added to allow user modifications in the PASAPAS table : REEV_THE & REEV_MEC. | REEV_THE TAB1 0 | REEV_MEC TAB1 0 | | LOOP on the time steps | LOOP UNTIL NB_BOTH | | | | THERMAL CALCULATION INITIALIZATION | | THERMAL CALCULATION PROCEDURE : results stored in table CHTER | | | | TAB1.'CHTER' = CHTER | | REEV_THE TAB1 1 | | | | MECHANICAL CALCULATION INITIALIZATION | | MECHANICAL CALCULATION PROCEDURE : results stored in table TT | | | | TAB1.'TT' = TT | | REEV_MEC TAB1 1 | | | | The convergence in BO_BOTH is tested by the procedure PRO_CONV. | | The criterion is based on the relative error beetween the last | | two results calculated in the THERMAL CALCULATION PROCEDURE. | | The precision is set by the 'CRITERE_COHERENCE' input data | | of the PASAPAS table (1.E-2 by default). | | | Update of TAB1 by CHTER et TT | $$$$
© Cast3M 2003 - All rights reserved.
Disclaimer