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

- added sanity check for temperature tolerance 

- added temperature tolerance in numerics.config
This commit is contained in:
Christoph Kords 2009-07-02 16:37:48 +00:00
parent 1c35103524
commit fe2f3e4d36
3 changed files with 56 additions and 51 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
code/IO.f90
View File

@ -852,7 +852,9 @@ endfunction
case (274) case (274)
msg = 'Non-positive relative maximum value (upper bound) for GIA residual' msg = 'Non-positive relative maximum value (upper bound) for GIA residual'
case (275) case (275)
msg = 'Limit for GIA iteration too small' msg = 'Limit for GIA iteration too small'
case (276)
msg = 'Non-positive relative tolerance for temperature'
case (300) case (300)
msg = 'This material can only be used with elements with three direct stress components' msg = 'This material can only be used with elements with three direct stress components'
case (500) case (500)

33
code/numerics.config
View File

@ -1,20 +1,21 @@
### numerical parameters ### ### numerical parameters ###
relevantStrain 1.0e-7 # strain increment considered significant relevantStrain 1.0e-7 # strain increment considered significant
iJacoStiffness 1 # frequency of stiffness update iJacoStiffness 1 # frequency of stiffness update
iJacoLpresiduum 1 # frequency of Jacobian update of residuum in Lp iJacoLpresiduum 1 # frequency of Jacobian update of residuum in Lp
pert_Fg 1.0e-6 # strain perturbation for FEM Jacobi pert_Fg 1.0e-6 # strain perturbation for FEM Jacobi
nHomog 10 # homogenization loop limit nHomog 10 # homogenization loop limit
nCryst 20 # crystallite loop limit (only for debugging info, real loop limit is "subStepMin") nCryst 20 # crystallite loop limit (only for debugging info, real loop limit is "subStepMin")
nState 10 # state loop limit nState 10 # state loop limit
nStress 40 # stress loop limit nStress 40 # stress loop limit
subStepMin 1.0e-3 # minimum (relative) size of sub-step allowed during cutback in crystallite 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_crystalliteState 1.0e-6 # relative tolerance in crystallite state loop
rTol_crystalliteStress 1.0e-6 # relative tolerance in crystallite stress loop rTol_crystalliteTemperature 1.0e-6 # relative tolerance in crystallite temperature loop
aTol_crystalliteStress 1.0e-8 # absolute tolerance in crystallite stress loop rTol_crystalliteStress 1.0e-6 # relative tolerance in crystallite stress loop
aTol_crystalliteStress 1.0e-8 # absolute tolerance in crystallite stress loop
resToler 1.0e-4 # relative tolerance of residual in GIA iteration resToler 1.0e-4 # relative tolerance of residual in GIA iteration
resAbsol 1.0e+2 # absolute tolerance of residual in GIA iteration (corresponds to ~1 Pa) resAbsol 1.0e+2 # absolute tolerance of residual in GIA iteration (corresponds to ~1 Pa)
resBound 1.0e+1 # relative maximum value (upper bound) for GIA residual resBound 1.0e+1 # relative maximum value (upper bound) for GIA residual
NRiterMax 24 # maximum number of GIA iteration NRiterMax 24 # maximum number of GIA iteration

70
code/numerics.f90
View File

@ -68,6 +68,7 @@ subroutine numerics_init()
! nStress ! nStress
! subStepMin ! subStepMin
! rTol_crystalliteState ! rTol_crystalliteState
! rTol_crystalliteTemperature
! rTol_crystalliteStress ! rTol_crystalliteStress
! aTol_crystalliteStress ! aTol_crystalliteStress
! resToler ! resToler
@ -80,23 +81,23 @@ subroutine numerics_init()
write(6,*) write(6,*)
! initialize all parameters with standard values ! initialize all parameters with standard values
relevantStrain = 1.0e-7_pReal relevantStrain = 1.0e-7_pReal
iJacoStiffness = 1_pInt iJacoStiffness = 1_pInt
iJacoLpresiduum = 1_pInt iJacoLpresiduum = 1_pInt
pert_Fg = 1.0e-6_pReal pert_Fg = 1.0e-6_pReal
nHomog = 10_pInt nHomog = 10_pInt
nCryst = 20_pInt nCryst = 20_pInt
nState = 10_pInt nState = 10_pInt
nStress = 40_pInt nStress = 40_pInt
subStepMin = 1.0e-3_pReal subStepMin = 1.0e-3_pReal
rTol_crystalliteState = 1.0e-6_pReal rTol_crystalliteState = 1.0e-6_pReal
rTol_crystalliteTemperature = 1.0e-6_pReal rTol_crystalliteTemperature = 1.0e-6_pReal
rTol_crystalliteStress = 1.0e-6_pReal rTol_crystalliteStress = 1.0e-6_pReal
aTol_crystalliteStress = 1.0e-8_pReal aTol_crystalliteStress = 1.0e-8_pReal
resToler = 1.0e-4_pReal resToler = 1.0e-4_pReal
resAbsol = 1.0e+2_pReal resAbsol = 1.0e+2_pReal
resBound = 1.0e+1_pReal resBound = 1.0e+1_pReal
NRiterMax = 24_pInt NRiterMax = 24_pInt
! try to open the config file ! try to open the config file
if(IO_open_file(fileunit,numerics_configFile)) then if(IO_open_file(fileunit,numerics_configFile)) then
@ -132,7 +133,7 @@ subroutine numerics_init()
subStepMin = IO_floatValue(line,positions,2) subStepMin = IO_floatValue(line,positions,2)
case ('rtol_crystallitestate') case ('rtol_crystallitestate')
rTol_crystalliteState = IO_floatValue(line,positions,2) rTol_crystalliteState = IO_floatValue(line,positions,2)
case ('rtol_crystalliteTemperature') case ('rtol_crystallitetemperature')
rTol_crystalliteTemperature = IO_floatValue(line,positions,2) rTol_crystalliteTemperature = IO_floatValue(line,positions,2)
case ('rtol_crystallitestress') case ('rtol_crystallitestress')
rTol_crystalliteStress = IO_floatValue(line,positions,2) rTol_crystalliteStress = IO_floatValue(line,positions,2)
@ -178,23 +179,24 @@ subroutine numerics_init()
write(6,'(a24,x,i8)') 'NRiterMax: ',NRiterMax write(6,'(a24,x,i8)') 'NRiterMax: ',NRiterMax
write(6,*) write(6,*)
! sanity check (Temperature check missing!!!!!!!) ! sanity check
if (relevantStrain <= 0.0_pReal) call IO_error(260) if (relevantStrain <= 0.0_pReal) call IO_error(260)
if (iJacoStiffness < 1_pInt) call IO_error(261) if (iJacoStiffness < 1_pInt) call IO_error(261)
if (iJacoLpresiduum < 1_pInt) call IO_error(262) if (iJacoLpresiduum < 1_pInt) call IO_error(262)
if (pert_Fg <= 0.0_pReal) call IO_error(263) if (pert_Fg <= 0.0_pReal) call IO_error(263)
if (nHomog < 1_pInt) call IO_error(264) if (nHomog < 1_pInt) call IO_error(264)
if (nCryst < 1_pInt) call IO_error(265) if (nCryst < 1_pInt) call IO_error(265)
if (nState < 1_pInt) call IO_error(266) if (nState < 1_pInt) call IO_error(266)
if (nStress < 1_pInt) call IO_error(267) if (nStress < 1_pInt) call IO_error(267)
if (subStepMin <= 0.0_pReal) call IO_error(268) if (subStepMin <= 0.0_pReal) call IO_error(268)
if (rTol_crystalliteState <= 0.0_pReal) call IO_error(269) if (rTol_crystalliteState <= 0.0_pReal) call IO_error(269)
if (rTol_crystalliteStress <= 0.0_pReal) call IO_error(270) if (rTol_crystalliteTemperature <= 0.0_pReal) call IO_error(276)
if (aTol_crystalliteStress <= 0.0_pReal) call IO_error(271) if (rTol_crystalliteStress <= 0.0_pReal) call IO_error(270)
if (resToler <= 0.0_pReal) call IO_error(272) if (aTol_crystalliteStress <= 0.0_pReal) call IO_error(271)
if (resAbsol <= 0.0_pReal) call IO_error(273) if (resToler <= 0.0_pReal) call IO_error(272)
if (resBound <= 0.0_pReal) call IO_error(274) if (resAbsol <= 0.0_pReal) call IO_error(273)
if (NRiterMax < 1_pInt) call IO_error(275) if (resBound <= 0.0_pReal) call IO_error(274)
if (NRiterMax < 1_pInt) call IO_error(275)
endsubroutine endsubroutine