Sidebar
lagamex:nauto
Table of Contents
Execution data for non automatic strategy
1st line (14I5) - Starting and saving files
| Variable | Values | Description |
|---|---|---|
| NTIN5 | =2,3, or 16 | File from which data are recovered at the beginning of the present execution |
| NTOUT10 | =2 or 3 | File on which the results of the present execution will be saved (NTIN and NTOUT can be the same file) |
| KNSYM15 | =1 | Matrice d'itération non symétrique stockée en mémoire centrale (NSYSOL) - Méthode peu performante |
| = 2 | Matrice symétrique, hauteur de colonne active stockée en mémoire centrale (COLSOL) - Méthode optimale pour les problèmes symétriques | |
| = 3 | Morse storage, METIS renumbering, Direct Solver LU symbolic and real factorization | |
| = 4 | Skyline storage, Direct solver LU | |
| = 5 | Idem as method 3, with parallel real factorization (CAESAR library) | |
| = 6 | GMRES method coupling with incomplete LU preconditionner, Morse storage (see appendix 19) |
|
| = 9 | Morse storage, PARDISO renumbering, Direct Solver LU symbolic and real factorization | |
| = ±10 | Iterative solver with mixed constraint preconditioner (GMRES or BiCGstab methods) → see appendix 22 | |
| IPRES20 | = 0 | if FMULT = DMULT = 0 : the strategy is based on time If FMULT = 0 and DMULT ≠ 0 : the strategy is based on displacements If FMULT ≠ 0 et DMULT = 0 : the strategy is based on forces If FMULT ≠ 0 and DMULT ≠ 0 : impossible |
| ≠ 0 | The strategy is based on time | |
| = 1, 3, 13, 18, 38 | read imposed DOF on file 31 ({namdat}.dep) (see appendix 2) Rem.: If you add 100 (e.g. 101 or 103), idem but with a special strategy adapted to large DEP files (computation time always smaller but does not work with the periodic loading or other special cases) |
|
| = 2, 3, 28, 38 | read imposed forces on file 32 ({namdat}.loa) (see appendix 2) | |
| = 8, 18, 28, 38 | read force multiplier FMULT and displacement multiplier DMULT on file 33 ({namdat}.lic) (see appendix 3) | |
| = 10, 11 or 13 | read imposed relations between the D.O.F. of generalized plane strain state on file 36 (see appendix 13). | |
| = 11, 12 or 13 | read geometry of cylinders in generalized plane strain state on file 35 (see appendix 13) | |
| = 19 | read macroscopic strain (or vector L) for periodic limit boundary conditions. (see appendix 16) → file *.DEM (n°30) | |
| IDENT25 | = 0 | No call of PRISUM and PRISIG and OCASFO routines |
| ≠ 0 | Call of routines PRISUM and PRISIG (see Note) + Lagamine inverse if IOPT≠0 | |
| = 9999 | Computes and saves the reactions following the 3 DOF (X, Y, T, in X, Y, Z according to the analysis type) for the first THREE foundations, according to time (format: 10(1PG10.3), compatible with DESCOU) | |
| = 9998 | Computes and saves the sum of the reactions DOF by DOF | |
| IDSPE30 | = 0 | No effect |
| = 2 | Superposition of spectral bands | |
| IARCL35 | = 0 | No effect |
| = 1 | Method with spherical step (Not available for non automatic strategy) | |
| ICCOR40 | = 0 | Nothing |
| IOPT45 | = 0 | Normal analyze |
| = 1 or 2 | Optimization of parameters by inverse method | |
| IADREM50 | = 0 | Nothing |
| = 1 | Adaptative remeshing | |
| NOWAR55 | =0 | Normal printing in ex.out |
| =1 | Do not print warning ( WARNING - ELEMB TERME DIAGONAL NO XXX NUL OU NEGATIF) in ex.out. This option is especially useful when using switch to reduce the size of ex.out |
Note
PRISUM: to print and to treat nodal values (coord, velocities, reactions, …) according to the user's personal choice
PRISIG: to print and to treat values at integration points (stresses, state variables, …) according to the user's personal choice
2nd line (18A4) - Title
3rd line (14I5) - Strategy control data
| ISTEP5 | N° of the first step to be computed in the present execution |
| JSTEP10 | Number of the last step to be computed in the present execution |
| ISAME15 | Step similitude index |
| = 0 if all the steps are identical, only one step defined | |
| ≠ 0 each step must be defined independently | |
| ISAVE20 | = -1 only the last computed step is saved on NTOUT (the previous ones are deleted) |
| = 0 No result saved or saving defined step by step (ISAUV) | |
| = n > 0 the results are saved every n step on NTOUT | |
| = 0 no result save. The configuration at the beginning of the step is saved in the CONECOLD vector (NUMNP, NSPAC, NCONE) | |
| = n saving on the file NTOUT once every n steps. Only the last saved step is kept (the previous ones are deleted). The purpose is to speed up the simulation by reducing the access to the disk. In that case, PRISIG and PRISUM are only called when a saving on the file NTOUT occurs |
|
| IOLD25 | N° of the set of results from which the present execution will start (in case of NTIN=OTO) |
| ISTR(1)30 | The tangent matrix is recomputed at the following iterations: ISTR1, ISTR2, ISTR2+ISTR3, ISTR2+2*ISTR3, …, ISTR2+K*ISTR3 |
| ISTR(2)35 | |
| ISTR(3)40 | |
| MAXIT45 | Maximum number of iterations per step. Default value : 5 |
| NSWIT50 | = 0 No switch |
| = 1 read switch data on NTSWI file (see appendix 4) | |
| NEXPT55 | = 0 in dynamic analysis, implicit scheme |
| < 0 in dynamic analysis, explicit scheme | |
| = n > 0 in dynamic analysis, mixed scheme; steps n, 2n, 3n are implicit, the others are explicit | |
| NPRIT60 | = 0 no particular printing of nodal values (on file .IPN), element integration points values (on file .IPE) or reactions values (on file .IPR) |
| = 1 reading of the file .PRI, for printing concerning nodal values (on file .IPN), element integration points values (on file .IPE) or reactions values (on file .IPR) - see appendix 9 | |
| ILSAV65 | = 0 nothing |
| = 1 change of ALSAV format to G15.0 |
4th line (14I5)
| ISTRA(3)5 | elepls → unity = 1 integration with initial VGRAD (default value) = 2 integration with actualized VGRAD = 3 integration with average value of VGRAD on the step (optimal value) Be carreful: = -1 if law Dafalias; -3 if law Cloe not used for PLA2EA and ELLI2EA because defined as parameter in the law |
| Law ARB → tens Mechanical or coupled ARB law integrated according to: = 0 return mapping method = 1 Li Kaiping method = 2 old method → hundreds Coupled ARB law = 1 parameters computed at the average temperature of the step = 0 parameters computed at the average temperature of the sub-interval |
|
| Law EVP-NH (Norton-Hoff 2D or 3D) → unity = 0 Resolution method with the analytical matrix = 1 Resolution method by perturbations → tens = 0 Resolution method using VGRADmean = 1 Resolution method using VGRADinitial |
|
| JET2D → hundreds = 0 Jaumann corotational method = 1 GREEN-NAGHDI corot. method (total corot. description) = 2 GREEN-NAGHDI corot. method (updated corot. description) |
|
| ISTRA(4)10 | Only for contact elements → unity: Penalty method = 0 Tangent stiffness matrix computed analytically = 1 Tangent stiffness matrix computed by perturbation = 2 In case of contact between two deformable solids, the skyline and the size of the global stiffness matrix are updated according to the evolution of the effective contact surface and the analytical method is used = 3 same as 2, but the numerical matrix is used → tens: Augmented lagrangian method for the CNTCP = 0 augmentation is realized at the beginning of the step |
| ISTRA(5)15 | For plane shell elements only AND mechanical analysis = 0 no effect = 1 elastic unloading (laws ILYSH and HILSH)  In case of thermo-mechanical analysis, ISTRA(5) is used for automatic step reduction (ISTRA(5)=1 leads to step reduction) |
| ISTRA(6)20 | NTNOR: Choice of a norm for convergence on nodal forces |
| = 0 average of the squares: \[\sqrt{\frac{\displaystyle\sum_{DIM}\frac{\displaystyle\sum_{DIR} \frac{FHE^2}{NEQUA}}{\displaystyle\sum_{DIR} \frac{REACT^2}{NREAC}}}{NDIM}}\] | |
| = 1 average of the absolute values: \[\frac{\displaystyle\sum_{DIM}\frac{\displaystyle\sum_{DIR} \frac{|FHE|}{NEQUA}}{\displaystyle\sum_{DIR} \frac{|REACT|}{NREAC}}}{NDIM}\] | |
| = 2 average of the maximum values: \[\frac{\displaystyle\sum_{DIM}\frac{\displaystyle\max_{DIR} \frac{|FHE|}{NEQUA}}{\displaystyle\max_{DIR} \frac{|REACT|}{NREAC}}}{NDIM}\] |
5th line (7G10.0)
For mechanical analysis
| STRAT(1)10 | ≈Δε | Perturbation for the computation of compliance matrix or stiffness matrix by the perturbation technique (use 10-5 to 10-8) |
For thermal analysis
| STRAT(2)20 | θ | Time integration parameter ∈ [0,1]. Integration scheme is stable if θ ≥ 0.5 |
| STRAT(3)30 | β | Second time integration parameter |
Note:
| Galerkin | θ = 2/3 | β = 1/6 or 2/9 |
| Zienkiewicz | θ = 1/2 | β= 1/6 |
| Cranck-Nicolson | θ = 1/2 | β = 1/4 |
| Implicit | θ = 1 | β = 0 |
For seepage pollutant flow transport with fixed mesh and Eulerian Lagrangian method (ICRIT=3):
STRAT(4)40 to STRAT(7)70
For dynamic analysis
Newmark parameters β and γ (see appendix 10)
If β<0: $\beta=(1+\alpha)^2/4$ and $\gamma=0.5+\alpha k$ (better compromise between stability, numerical damping and frequency distortion 0 ≤ α ≤ 1; 0 ≤ k ≤ 0.5)
| STRAT(2)20 | β or -α | |
| STRAT(3)30 | γ or k | |
| STRAT(4)40 | Damping coefficent α | C = αM + βK For dynamic implicit (or mixed) simulations, used in routines DYJT3D for 8-nodes brick elements and DYCQJ4 for COQJ4 |
| STRAT(4)40 | Damping coefficent β | |
For any type of analysis
| STRAT(6)60 | REDUF - out-of-balance forces reduction parameters = numerical damping, use to achieve convergence in difficult cases with cyclic oscillations; 0,2 means that 80 % of the out-of-balance forces are taken into account. REDUF can decrease in function of the iteration number. If STRAT(6) < 1, REDUF remains constant. Else the integer part gives the final number of iteration before the REDUF becomes null and the decimal part gives the initial value of the STRAT(6) Example : If STRAT(6) =6,6 ⇒ STRAT(6)(ITER=1)=0.6 and STRAT(6)(ITER=6)=0.0 and REDUF = 1 - STRAT(6)  |
| STRAT(7)70 | DMAX - maximum penetration distance in case of contact (used to avoid “false contact”) Default value: 1D+30 |
| DMIN_FOUND80 | In 3D contact (CFI3D element or contact element associated with shell element COQJ4), this parameter is used as the minimum size of the foundations limits (i.e. if the element is not in the limits, the routine is not called). Default value: maximum size of foundation limit divided by 100 |
6th line - Automatic strategy and errors (7G10.0)
STRAT(8)10 to STRAT(11)40: Reserved for automatic strategy.
| STRAT(12)50 | If ICRIT = 2 or 5: % max new “elements” to be remeshed from the last remeshing Else ETASM %max admissible error on stresses |
| STRAT(13)60 | ETAGM % max admissible error on geometric remeshing criteria |
7th line - Convergence parameters (7G10.0)
In the computation of the norm for convergence, the reference efforts (REACT in norm on nodal forces) or reference displacement (DNOR in norm on displacements) can be defined using STRAT(15) to STRAT(20).
If these values are equal to zero, the convergence norms are relative, it is generally the case.
If these values are negative, the convergence norms are generally relatives except if the reaction norm (dimension per dimension) is smaller than the value introduced hereafter. In this case, the norm becomes absolute for the considered dimension. It is then a limit value, a minimum of the reaction norms.
For a detailed explanation, see appendix 15.
| STRAT(15)10 | Mechanical force |
| STRAT(16)20 | Nodal flux/Fluid flux in soil mechanics |
| STRAT(17)30 | Moments/Thermal flux in soil mechanics |
| STRAT(18)40 | Displacement |
| STRAT(19)50 | Variation of temperature/Fluid pressure variations in soil mechanics |
| STRAT(20)60 | Rotation/Temperature variations in soil mechanics |
| STRAT(21)70 | In case of remeshing, minimum interpenetration distance |
8th line (and 9th if required) - Printing control (15I5/14I5)
Only the data structure is described hereafter. The proposed possibilities are detailed in appendix 1. The 9th line is only necessary if IOPT(I) equals to 2, 3 or 5.
| IOPT(I) I=1,15 | Printing option |
| LISTE(14) | Definition of the list of nodes or elements selected for printing, if required |
9th & 10th (or 10th & 11th) lines - Definition of increments (5G10.0, 6I5 / 14I5)
These two lines must be repeated for each defined increment.
| DELTAT10 | Time increment |
| FMULT20 | Increment of the imposed or implicit force multiplier |
| DMULT30 | Increment of the imposed displacements multiplier |
| PRECU40 | Precision for convergence on displacements increments |
| PRECF50 | Precision for out-of-balance forces |
| INITV55 | Initialization of speeds at the beginning of the step = 0 speeds of the preceding step = 1 all the speeds are equal to zero at the beginning of the step until the time increment increases For dynamic analysis, see appendix 10 |
| ISTR(I)60, 65, 70 I=1,3 | Iterations where the stiffness matrix is computed. If the values are equal to zero, the initial definitions are used again |
| IPRECT75 | = 0 Convergence is obtained as far as one of the criteria (force or displacement) is achieved = 1 Convergence is obtained when both criteria are achieved = 2 Convergence is obtained when the displacement criterion is achieved = 3 Convergence is obtained when the force criterion is achieved = -1 in linear analysis, convergence is imposed after a resolution |
| ISAUV80 | = 0 no saving on the NTOUT file for this time step with ISAVE = 0 = 1 saving on the NTOUT file for this time step with ISAVE = 0 |
| IOPT(I) I =1,15 | If IOPT = 0, the initial definitions for IOPT(I) , I= 1 to 10 are used again and the last given value for IOPT(I), I=11 to 15 are kept. According to the given value, read the lists (14I5) if necessary |
lagamex/nauto.txt · Last modified: 2020/08/25 15:46 (external edit)
