This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
lagamex:auto [2019/07/10 13:51] helene |
lagamex:auto [2024/01/25 15:59] (current) arthur |
||
---|---|---|---|
Line 13: | Line 13: | ||
|:::|= 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 [[appendices:a19|appendix 19]])| | |:::|= 6|GMRES method coupling with incomplete LU preconditionner, \\ Morse storage (see [[appendices:a19|appendix 19]])| | ||
+ | |:::|= 8|Parallel solver SOLVE_DSS| | ||
|:::|= 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 [[appendices:a22|appendix 22]]| | |:::|= ±10|Iterative solver with mixed constraint preconditioner (GMRES or BiCGstab methods) → see [[appendices:a22|appendix 22]]| | ||
Line 30: | Line 31: | ||
|:::|= 2|Superposition of spectral bands| | |:::|= 2|Superposition of spectral bands| | ||
|IARCL<sup>35</sup>|= 0|No effect| | |IARCL<sup>35</sup>|= 0|No effect| | ||
- | |:::|= 1|Method with spherical step| | + | |:::|> 0|Method with spherical step (see [[lagamex:autosph|spherical steps]])| |
|ICCOR<sup>40</sup>|= 0|Nothing| | |ICCOR<sup>40</sup>|= 0|Nothing| | ||
|IOPT<sup>45</sup>|= 0|Normal analyze| | |IOPT<sup>45</sup>|= 0|Normal analyze| | ||
Line 75: | Line 76: | ||
|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| | ||
+ | |:::|= 2 ALSAV defined cyclically (see [[lagamex:auto#Last lines|Last lines]])| | ||
=== (1) Signification of JSTEP === | === (1) Signification of JSTEP === | ||
Line 129: | Line 131: | ||
|Cranck-Nicolson |θ = 1/2| β = 1/4 | | |Cranck-Nicolson |θ = 1/2| β = 1/4 | | ||
|Implicit| θ = 1| β = 0 | | |Implicit| θ = 1| β = 0 | | ||
+ | |||
+ | The integration scheme is made of two parts : \\ | ||
+ | - the first one (θ) is about the referencial time for the energy balance \\ | ||
+ | - the second one (β) is about the referencial time for the material parameters\\ | ||
+ | |||
+ | which gives : | ||
+ | \[ σ^{t+1} = [(1-θ-β) ALAMX + β ALBMX] ε^{t} + [β ALAMX + (θ-β) ALBMX] ε^{t+1} \] | ||
+ | |||
+ | Where ALAMX represent material properties at time t and ALBMX at t+1 | ||
+ | |||
+ | |||
=== For seepage pollutant flow transport with fixed mesh and Eulerian Lagrangian method (ICRIT=3): === | === For seepage pollutant flow transport with fixed mesh and Eulerian Lagrangian method (ICRIT=3): === | ||
STRAT(4)<sup>40</sup> to STRAT(7)<sup>70</sup> | STRAT(4)<sup>40</sup> to STRAT(7)<sup>70</sup> | ||
Line 195: | Line 208: | ||
---- | ---- | ||
- | ===== Last lines (2G10.0, I5, G10.0) ===== | + | ===== Last lines ===== |
- | Repeat max 15000 times. :!: Over 15000 times this can lead to problems in the execution such as wrong value of DMAXMU. | + | Repeat max 15000 times. :!: Over 15000 times this can lead to problems in the execution such as wrong value of DMAXMU. |
+ | ==== IF ILSAV = 0 (2G10.0, I5, G10.0) ==== | ||
|ALSAV<sup>10</sup>|Multiplier of the imposed forces or displacements for which one wishes a printing according to IOPT and a saving on NRESU (oto file)| | |ALSAV<sup>10</sup>|Multiplier of the imposed forces or displacements for which one wishes a printing according to IOPT and a saving on NRESU (oto file)| | ||
|DELTAT_OTO<sup>20</sup>|If ≠ 0, new value of the time step after a saving on NRESU| | |DELTAT_OTO<sup>20</sup>|If ≠ 0, new value of the time step after a saving on NRESU| | ||
Line 204: | Line 218: | ||
After the last multiplier (when the time is bigger than the largest value of ALSAV), it saves all the steps. \\ | After the last multiplier (when the time is bigger than the largest value of ALSAV), it saves all the steps. \\ | ||
:!: do NOT put a zero in the list (otherwise, the rest of the list will be skipped and every step will be saved) | :!: do NOT put a zero in the list (otherwise, the rest of the list will be skipped and every step will be saved) | ||
+ | ==== IF ILSAV = 1 (G15.0, G10.0, I5, G10.0) ==== | ||
+ | |ALSAV<sup>10</sup>|Multiplier of the imposed forces or displacements for which one wishes a printing according to IOPT and a saving on NRESU (oto file)| | ||
+ | |DELTAT_OTO<sup>20</sup>|If ≠ 0, new value of the time step after a saving on NRESU| | ||
+ | |INIT_OTO<sup>25</sup>|If = 1, reset the nodal speed to 0 until good convergence| | ||
+ | |DMAXMU_OTO<sup>35</sup>|If ≠ 0, new value of the DMAXMU (max. value of the multiplier, see 6th line, col. 40) after a saving on NRESU \\ If = 0, DMAXMU is restored to its initial value| | ||
+ | A blank line indicates the end of the data. \\ | ||
+ | After the last multiplier (when the time is bigger than the largest value of ALSAV), it saves all the steps. \\ | ||
+ | :!: do NOT put a zero in the list (otherwise, the rest of the list will be skipped and every step will be saved) | ||
+ | ==== If ILSAV = 2 ==== | ||
+ | ILSAV = 2 allows to define the printing times cyclically. This is particularly useful for cyclic loadings in combination with the use cyclic definition in the .LOA or .DEP file (see [[appendices:a2|Appendix 2: Non radial loading paths (1)]]). | ||
+ | ^1st line (G10.0, I5)^^ | ||
+ | |Period| Period of the cyclic definition for printing| | ||
+ | |NTIME| Number of printings per period| | ||
+ | ^Line 2 to NTIME + 1 (G10.0)^^ | ||
+ | |t<sub>i</sub>|Time at which one wishes a printing. \\ A printing will be made for every time t<sub>i</sub>+K*Period as long as that time<ALAMBF and K*NTIME<50000| | ||