absolute tolerance was just fine -- residuum is in Lp not in stress. Hence, unit is strain and not Pa...

This commit is contained in:
Philip Eisenlohr 2009-10-19 12:53:56 +00:00
parent ba27947142
commit a94e551c5f
2 changed files with 4 additions and 4 deletions

View File

@ -10,9 +10,9 @@ nCryst 20 # crystallite loop limit (only for debuggi
nState 50 # state loop limit
nStress 200 # stress loop limit
subStepMin 1.0e-3 # minimum (relative) size of sub-step allowed during cutback in crystallite
rTol_crystalliteState 1.0e-6 # relative tolerance in crystallite state loop
rTol_crystalliteStress 1.0e-6 # relative tolerance in crystallite stress loop
aTol_crystalliteStress 1.0e+0 # absolute tolerance in crystallite stress loop
rTol_crystalliteState 1.0e-6 # relative tolerance in crystallite state loop (abs tol provided by constitutive law)
rTol_crystalliteStress 1.0e-6 # relative tolerance in crystallite stress loop (Lp residuum)
aTol_crystalliteStress 1.0e-8 # absolute tolerance in crystallite stress loop (Lp residuum!)
aTol_RGC 1.0e+4 # absolute tolerance of RGC residuum (in Pa)
rTol_RGC 1.0e-3 # relative ...

View File

@ -84,7 +84,7 @@ subroutine numerics_init()
rTol_crystalliteState = 1.0e-6_pReal
rTol_crystalliteTemperature = 1.0e-6_pReal
rTol_crystalliteStress = 1.0e-6_pReal
aTol_crystalliteStress = 1.0e+0_pReal
aTol_crystalliteStress = 1.0e-8_pReal ! residuum is in Lp (hence strain on the order of 1e-8 here)
!* RGC parameters: added <<<updated 31.07.2009>>>
absTol_RGC = 1.0e+3