diff --git a/examples/config/numerics.yaml b/examples/config/numerics.yaml index 4ff1d3863..e04888055 100644 --- a/examples/config/numerics.yaml +++ b/examples/config/numerics.yaml @@ -1,28 +1,10 @@ # Default values of all available numerical parameters # 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: grid: 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: 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_abs_P: 1.0e3 # absolute 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: 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_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 - alpha: 1.0 # polarization scheme parameter 0.0 < alpha < 2.0 (alpha = 1.0 ==> AL scheme, alpha = 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) + 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 (1.0: AL scheme, 2.0: accelerated scheme) 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 -mesh: - maxCutBack: 3 # maximum cut back 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 +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 (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: mechanical: @@ -75,7 +64,7 @@ phase: N_iter_state_max: 10 # state loop limit 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_abs_Lp: 1.0e-8 # absolute tolerance in Lp residuum 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) 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_abs_Li: 1.0e-8 # absolute tolerance in Li residuum N_iter_Li_max: 40 # stress loop limit for Li @@ -95,3 +84,14 @@ commercialFEM: generic: random_seed: 0 # fixed seeding for pseudo-random number generator (0: use random seed) 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