223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

more sensible sorting of sections; additional polish

This commit is contained in:
Philip Eisenlohr 2023-09-05 13:14:20 -04:00
parent 8399a11703
commit 2434c322ca
1 changed files with 34 additions and 34 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

68
examples/config/numerics.yaml
View File

@ -1,28 +1,10 @@
# Default values of all available numerical parameters # Default values of all available numerical parameters
# Please note that keys are case sensitive # Please note that keys are case sensitive
homogenization:
mechanical:
RGC:
eps_abs_P: 1.0e+4 # absolute tolerance of RGC residuum (in Pa)
eps_rel_P: 1.0e-3 # relative ...
eps_abs_max: 1.0e+10 # absolute upper-limit of RGC residuum (in Pa)
eps_rel_max: 1.0e+2 # relative ...
Delta_a: 1.0e-7 # perturbation for computing penalty tangent
relevant_mismatch: 1.0e-5 # minimum threshold of mismatch
viscosity_exponent: 1.0e+0 # power (sensitivity rate) of numerical viscosity in RGC scheme
viscosity_modulus: 0.0e+0 # stress modulus of RGC numerical viscosity (zero = without numerical viscosity)
# suggestion: larger than the aTol_RGC but still far below the expected flow stress of material
dot_a_ref: 1.0e-3 # reference rate of relaxation (about the same magnitude as straining rate, possibly a bit higher)
dot_a_max: 1.0e+0 # threshold of maximum relaxation vector increment (if exceed this then cutback)
Delta_V_max: 1.0e-5 # maximum allowable relative volume discrepancy
Delta_V_modulus: 1.0e+12
Delta_V_exponent: 5.0
solver: solver:
grid: grid:
N_staggered_iter_max: 10 # max number of field-level staggered iterations N_staggered_iter_max: 10 # max number of field-level staggered iterations
N_cutback_max: 3 # maximum cut-back level (0: 1, 1: 0.5, 2: 0.25, etc) N_cutback_max: 3 # maximum cutback level (0: 1, 1: 0.5, 2: 0.25, etc)
damage: damage:
N_iter_max: 100 # maximum iteration number N_iter_max: 100 # maximum iteration number
@ -41,7 +23,7 @@ solver:
eps_rel_div(P): 5.0e-4 # relative tolerance for fulfillment of stress equilibrium eps_rel_div(P): 5.0e-4 # relative tolerance for fulfillment of stress equilibrium
eps_abs_P: 1.0e3 # absolute tolerance for fulfillment of stress BC eps_abs_P: 1.0e3 # absolute tolerance for fulfillment of stress BC
eps_rel_P: 1.0e-3 # relative tolerance for fulfillment of stress BC eps_rel_P: 1.0e-3 # relative tolerance for fulfillment of stress BC
update_gamma: false # update Gamma-operator with current dPdF (not possible if FFT:memory_efficient=true) update_gamma: false # update Gamma-operator with current dPdF (not possible if FFT: memory_efficient == true)
FFT: FFT:
memory_efficient: true # precalculate Gamma-operator (81 doubles per point) memory_efficient: true # precalculate Gamma-operator (81 doubles per point)
@ -50,21 +32,28 @@ solver:
FFTW_plan_mode: FFTW_MEASURE # planning-rigor flags, see manual at https://www.fftw.org/fftw3_doc/Planner-Flags.html FFTW_plan_mode: FFTW_MEASURE # planning-rigor flags, see manual at https://www.fftw.org/fftw3_doc/Planner-Flags.html
FFTW_timelimit: -1.0 # time limit of plan creation for FFTW, see manual on www.fftw.org. (-1.0: disable time limit) FFTW_timelimit: -1.0 # time limit of plan creation for FFTW, see manual on www.fftw.org. (-1.0: disable time limit)
PETSc_options: -snes_type ngmres -snes_ngmres_anderson # PETSc solver options PETSc_options: -snes_type ngmres -snes_ngmres_anderson # PETSc solver options
alpha: 1.0 # polarization scheme parameter 0.0 < alpha < 2.0 (alpha = 1.0 ==> AL scheme, alpha = 2.0 ==> accelerated scheme) alpha: 1.0 # polarization scheme parameter 0.0 < alpha < 2.0 (1.0: AL scheme, 2.0: accelerated scheme)
beta: 1.0 # polarization scheme parameter 0.0 < beta < 2.0 ( beta = 1.0 ==> AL scheme, beta = 2.0 ==> accelerated scheme) beta: 1.0 # polarization scheme parameter 0.0 < beta < 2.0 (1.0: AL scheme, 2.0: accelerated scheme)
eps_abs_curl(F): 1.0e-10 # absolute tolerance for fulfillment of strain compatibility eps_abs_curl(F): 1.0e-10 # absolute tolerance for fulfillment of strain compatibility
eps_rel_curl(F): 5.0e-4 # relative tolerance for fulfillment of strain compatibility eps_rel_curl(F): 5.0e-4 # relative tolerance for fulfillment of strain compatibility
mesh: homogenization:
maxCutBack: 3 # maximum cut back level (0: 1, 1: 0.5, 2: 0.25, etc) mechanical:
maxStaggeredIter: 10 # max number of field level staggered iterations RGC:
structorder: 2 # order of displacement shape functions eps_abs_P: 1.0e+4 # absolute tolerance of RGC residuum (in Pa)
bbarstabilisation: false eps_rel_P: 1.0e-3 # relative ...
integrationorder: 2 # order of quadrature rule required eps_abs_max: 1.0e+10 # absolute upper-limit of RGC residuum (in Pa)
itmax: 250 # maximum iteration number eps_rel_max: 1.0e+2 # relative ...
itmin: 2 # minimum iteration number Delta_a: 1.0e-7 # perturbation for computing penalty tangent
eps_struct_atol: 1.0e-10 # absolute tolerance for mechanical equilibrium relevant_mismatch: 1.0e-5 # minimum threshold of mismatch
eps_struct_rtol: 1.0e-4 # relative tolerance for mechanical equilibrium viscosity_exponent: 1.0e+0 # power (sensitivity rate) of numerical viscosity in RGC scheme
viscosity_modulus: 0.0e+0 # stress modulus of RGC numerical viscosity (0: without numerical viscosity)
# suggestion: larger than the aTol_RGC but still far below the expected flow stress of material
dot_a_ref: 1.0e-3 # reference rate of relaxation (about the same magnitude as straining rate, possibly a bit higher)
dot_a_max: 1.0e+0 # threshold of maximum relaxation vector increment (if exceed this then cutback)
Delta_V_max: 1.0e-5 # maximum allowable relative volume discrepancy
Delta_V_modulus: 1.0e+12
Delta_V_exponent: 5.0
phase: phase:
mechanical: mechanical:
@ -75,7 +64,7 @@ phase:
N_iter_state_max: 10 # state loop limit N_iter_state_max: 10 # state loop limit
plastic: plastic:
r_linesearch_Lp: 0.5 # factor to decrease the step due to non-convergence in Lp calculation r_linesearch_Lp: 0.5 # factor to decrease the step if Lp calculation fails to converge
eps_rel_Lp: 1.0e-6 # relative tolerance in Lp residuum eps_rel_Lp: 1.0e-6 # relative tolerance in Lp residuum
eps_abs_Lp: 1.0e-8 # absolute tolerance in Lp residuum eps_abs_Lp: 1.0e-8 # absolute tolerance in Lp residuum
N_iter_Lp_max: 40 # stress loop limit for Lp N_iter_Lp_max: 40 # stress loop limit for Lp
@ -83,7 +72,7 @@ phase:
integrator_state: FPI # integration method (FPI = Fixed Point Iteration, Euler = Euler, AdaptiveEuler = Adaptive Euler, RK4 = classical 4th order Runge-Kutta, RKCK45 = 5th order Runge-Kutta Cash-Karp) integrator_state: FPI # integration method (FPI = Fixed Point Iteration, Euler = Euler, AdaptiveEuler = Adaptive Euler, RK4 = classical 4th order Runge-Kutta, RKCK45 = 5th order Runge-Kutta Cash-Karp)
eigen: eigen:
r_linesearch_Li: 0.5 # factor to decrease the step due to non-convergence in Li calculation r_linesearch_Li: 0.5 # factor to decrease the step if Li calculation fails to converge
eps_rel_Li: 1.0e-6 # relative tolerance in Li residuum eps_rel_Li: 1.0e-6 # relative tolerance in Li residuum
eps_abs_Li: 1.0e-8 # absolute tolerance in Li residuum eps_abs_Li: 1.0e-8 # absolute tolerance in Li residuum
N_iter_Li_max: 40 # stress loop limit for Li N_iter_Li_max: 40 # stress loop limit for Li
@ -95,3 +84,14 @@ commercialFEM:
generic: generic:
random_seed: 0 # fixed seeding for pseudo-random number generator (0: use random seed) random_seed: 0 # fixed seeding for pseudo-random number generator (0: use random seed)
phi_min: 1.0e-6 # non-zero residual damage phi_min: 1.0e-6 # non-zero residual damage
mesh:
maxCutBack: 3 # maximum cutback level (0: 1, 1: 0.5, 2: 0.25, etc.)
maxStaggeredIter: 10 # max number of field-level staggered iterations
structorder: 2 # order of displacement shape functions
bbarstabilisation: false
integrationorder: 2 # order of quadrature rule required
itmax: 250 # maximum iteration number
itmin: 2 # minimum iteration number
eps_struct_atol: 1.0e-10 # absolute tolerance for mechanical equilibrium
eps_struct_rtol: 1.0e-4 # relative tolerance for mechanical equilibrium