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

include sanity check

This commit is contained in:
Martin Diehl 2020-03-14 22:53:05 +01:00
parent 554631d554
commit 9ad391ae50
1 changed files with 7 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
src/constitutive_plastic_nonlocal.f90
View File

@ -214,6 +214,8 @@ module subroutine plastic_nonlocal_init
dst => microstructure(phase_plasticityInstance(p)), & dst => microstructure(phase_plasticityInstance(p)), &
config => config_phase(p)) config => config_phase(p))
prm%output = config%getStrings('(output)',defaultVal=emptyStringArray)
prm%atol_rho = config%getFloat('atol_rho', defaultVal=0.0_pReal) prm%atol_rho = config%getFloat('atol_rho', defaultVal=0.0_pReal)
prm%atol_gamma = config%getFloat('atol_shear', defaultVal=0.0_pReal) prm%atol_gamma = config%getFloat('atol_shear', defaultVal=0.0_pReal)
@ -223,7 +225,6 @@ module subroutine plastic_nonlocal_init
prm%mu = lattice_mu(p) prm%mu = lattice_mu(p)
prm%nu = lattice_nu(p) prm%nu = lattice_nu(p)
prm%Nslip = config%getInts('nslip',defaultVal=emptyIntArray) prm%Nslip = config%getInts('nslip',defaultVal=emptyIntArray)
prm%totalNslip = sum(prm%Nslip) prm%totalNslip = sum(prm%Nslip)
slipActive: if (prm%totalNslip > 0) then slipActive: if (prm%totalNslip > 0) then
@ -372,12 +373,10 @@ module subroutine plastic_nonlocal_init
extmsg = trim(extmsg)//' surfaceTransmissivity' extmsg = trim(extmsg)//' surfaceTransmissivity'
if (prm%fEdgeMultiplication < 0.0_pReal .or. prm%fEdgeMultiplication > 1.0_pReal) & if (prm%fEdgeMultiplication < 0.0_pReal .or. prm%fEdgeMultiplication > 1.0_pReal) &
extmsg = trim(extmsg)//' fEdgeMultiplication' extmsg = trim(extmsg)//' fEdgeMultiplication'
endif slipActive endif slipActive
prm%output = config%getStrings('(output)',defaultVal=emptyStringArray)
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! allocate state arrays ! allocate state arrays
NofMyPhase = count(material_phaseAt==p) * discretization_nIP NofMyPhase = count(material_phaseAt==p) * discretization_nIP
@ -483,6 +482,10 @@ module subroutine plastic_nonlocal_init
if (NofMyPhase > 0) call stateInit(p,NofMyPhase) if (NofMyPhase > 0) call stateInit(p,NofMyPhase)
plasticState(p)%state0 = plasticState(p)%state plasticState(p)%state0 = plasticState(p)%state
!--------------------------------------------------------------------------------------------------
! exit if any parameter is out of range
if (extmsg /= '') call IO_error(211,ext_msg=trim(extmsg)//'('//PLASTICITY_NONLOCAL_LABEL//')')
enddo enddo
allocate(compatibility(2,maxval(totalNslip),maxval(totalNslip),nIPneighbors,& allocate(compatibility(2,maxval(totalNslip),maxval(totalNslip),nIPneighbors,&