lagamex:nauto
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
lagamex:nauto [2019/06/20 14:35] helene |
lagamex:nauto [2020/08/25 15:46] (current) |
||
|---|---|---|---|
| Line 12: | Line 12: | ||
| |:::|= 4|Skyline storage, Direct solver LU| | |:::|= 4|Skyline storage, Direct solver LU| | ||
| |:::|= 5|Idem as method 3, with parallel real factorization (CAESAR library)| | |:::|= 5|Idem as method 3, with parallel real factorization (CAESAR library)| | ||
| - | |:::|= 6|GMRES method coupling with incomplete LU preconditionner, \\ Morse storage (see appendix 19)| | + | |:::|= 6|GMRES method coupling with incomplete LU preconditionner, \\ Morse storage (see [[appendices:a19|appendix 19]])| |
| |:::|= 9|Morse storage, PARDISO renumbering, Direct Solver LU symbolic and real factorization| | |:::|= 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| | + | |:::|= ±10|Iterative solver with mixed constraint preconditioner (GMRES or BiCGstab methods) → see [[appendices:a22|appendix 22]]| |
| |IPRES<sup>20</sup>| = 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| | |IPRES<sup>20</sup>| = 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| | |:::|≠ 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)| | + | |:::|= 1, 3, 13, 18, 38|read imposed DOF on file 31 ({namdat}.dep) (see [[appendices:a2|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)| | + | |:::|= 2, 3, 28, 38| read imposed forces on file 32 ({namdat}.loa) (see [[appendices:a2|appendix 2]])| |
| - | |:::|= 8, 18, 28, 38| read force multiplier FMULT and displacement multiplier DMULT on file 33 ({namdat}.lic) (see appendix 3)| | + | |:::|= 8, 18, 28, 38| read force multiplier FMULT and displacement multiplier DMULT on file 33 ({namdat}.lic) (see [[appendices:a3|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).| | + | |:::|= 10, 11 or 13 |read imposed relations between the D.O.F. of generalized plane strain state on file 36 (see [[appendices:a13|appendix 13]]).| |
| - | |:::|= 11, 12 or 13|read geometry of cylinders in generalized plane strain state on file 35 (see appendix 13)| | + | |:::|= 11, 12 or 13|read geometry of cylinders in generalized plane strain state on file 35 (see [[appendices:a13|appendix 13]])| |
| - | |:::|= 19|read macroscopic strain (or vector L) for periodic limit boundary conditions. (see appendix 16) → file *.DEM (n°30)| | + | |:::|= 19|read macroscopic strain (or vector L) for periodic limit boundary conditions. (see [[appendices:a16|appendix 16]]) → file *.DEM (n°30)| |
| |IDENT<sup>25</sup>| = 0|No call of PRISUM and PRISIG and OCASFO routines| | |IDENT<sup>25</sup>| = 0|No call of PRISUM and PRISIG and OCASFO routines| | ||
| |:::|≠ 0 | Call of routines PRISUM and PRISIG (see [[lagamex:auto#Note|Note]]) + Lagamine inverse if IOPT≠0| | |:::|≠ 0 | Call of routines PRISUM and PRISIG (see [[lagamex:auto#Note|Note]]) + Lagamine inverse if IOPT≠0| | ||
| Line 68: | Line 68: | ||
| |MAXIT<sup>45</sup>|Maximum number of iterations per step. Default value : 5| | |MAXIT<sup>45</sup>|Maximum number of iterations per step. Default value : 5| | ||
| |NSWIT<sup>50</sup>|= 0 No switch| | |NSWIT<sup>50</sup>|= 0 No switch| | ||
| - | |:::|= 1 read switch data on NTSWI file (see appendix 4bis)| | + | |:::|= 1 read switch data on NTSWI file (see [[appendices:a4|appendix 4]])| |
| |NEXPT<sup>55</sup>|= 0 in dynamic analysis, implicit scheme| | |NEXPT<sup>55</sup>|= 0 in dynamic analysis, implicit scheme| | ||
| |:::|< 0 in dynamic analysis, explicit scheme| | |:::|< 0 in dynamic analysis, explicit scheme| | ||
| |:::|= n > 0 in dynamic analysis, mixed scheme; steps n, 2n, 3n are implicit, the others are explicit| | |:::|= n > 0 in dynamic analysis, mixed scheme; steps n, 2n, 3n are implicit, the others are explicit| | ||
| |NPRIT<sup>60</sup>|= 0 no particular printing of nodal values (on file .IPN), element integration points values (on file .IPE) or reactions values (on file .IPR)| | |NPRIT<sup>60</sup>|= 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| | + | |:::|= 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 [[appendices:a9|appendix 9]]| |
| |ILSAV<sup>65</sup>|= 0 nothing| | |ILSAV<sup>65</sup>|= 0 nothing| | ||
| |:::|= 1 change of ALSAV format to G15.0| | |:::|= 1 change of ALSAV format to G15.0| | ||
| Line 109: | Line 109: | ||
| === For dynamic analysis === | === For dynamic analysis === | ||
| - | Newmark parameters β and γ (see Appendix 10) \\ | + | Newmark parameters β and γ (see [[appendices:a10|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) | 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)<sup>20</sup>|β or -α|| | |STRAT(2)<sup>20</sup>|β or -α|| | ||
| Line 133: | Line 133: | ||
| If these values are equal to zero, the convergence norms are __relative__, it is generally the case. \\ \\ | 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. \\ | 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. | + | For a detailed explanation, see [[appendices:a15|appendix 15]]. |
| |STRAT(15)<sup>10</sup>|Mechanical force| | |STRAT(15)<sup>10</sup>|Mechanical force| | ||
| Line 146: | Line 146: | ||
| ---- | ---- | ||
| ===== 8th line (and 9th if required) - Printing control (15I5/14I5) ===== | ===== 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. | + | Only the data structure is described hereafter. The proposed possibilities are detailed in [[appendices:a1|appendix 1]]. The 9th line is only necessary if IOPT(I) equals to 2, 3 or 5. |
| |IOPT(I) \\ I=1,15| Printing option| | |IOPT(I) \\ I=1,15| Printing option| | ||
| |LISTE(14) |Definition of the list of nodes or elements selected for printing, if required| | |LISTE(14) |Definition of the list of nodes or elements selected for printing, if required| | ||
| Line 160: | Line 160: | ||
| |PRECU<sup>40</sup>|Precision for convergence on displacements increments| | |PRECU<sup>40</sup>|Precision for convergence on displacements increments| | ||
| |PRECF<sup>50</sup>|Precision for out-of-balance forces| | |PRECF<sup>50</sup>|Precision for out-of-balance forces| | ||
| - | |INITV<sup>55</sup>|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| | + | |INITV<sup>55</sup>|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 [[appendices:a10|appendix 10]]| |
| |ISTR(I)<sup>60, 65, 70</sup> \\ I=1,3|Iterations where the stiffness matrix is computed. If the values are equal to zero, the initial definitions are used again| | |ISTR(I)<sup>60, 65, 70</sup> \\ I=1,3|Iterations where the stiffness matrix is computed. If the values are equal to zero, the initial definitions are used again| | ||
| |IPRECT<sup>75</sup>|= 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| | |IPRECT<sup>75</sup>|= 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| | ||
lagamex/nauto.1561034100.txt.gz · Last modified: 2020/08/25 15:34 (external edit)
Page Tools
