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

should be always set, even if no twinning is enabled 

fixes NaN in dotState
This commit is contained in:
Martin Diehl 2018-09-12 12:08:03 +02:00
parent 253f318e55
commit c9b5b3fb7b
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/plastic_phenopowerlaw.f90
View File

@ -220,10 +220,6 @@ subroutine plastic_phenopowerlaw_init
prm%gdot0_twin = config_phase(p)%getFloat('gdot0_twin')
prm%n_twin = config_phase(p)%getFloat('n_twin')
prm%spr = config_phase(p)%getFloat('s_pr')
prm%twinB = config_phase(p)%getFloat('twin_b')
prm%twinC = config_phase(p)%getFloat('twin_c')
prm%twinD = config_phase(p)%getFloat('twin_d')
prm%twinE = config_phase(p)%getFloat('twin_e')
prm%h0_TwinTwin = config_phase(p)%getFloat('h0_twintwin')
endif
@ -233,6 +229,11 @@ subroutine plastic_phenopowerlaw_init
prm%h0_TwinSlip = config_phase(p)%getFloat('h0_twinslip')
endif
prm%twinB = config_phase(p)%getFloat('twin_b',defaultVal=1.0_pReal)
prm%twinC = config_phase(p)%getFloat('twin_c',defaultVal=0.0_pReal)
prm%twinD = config_phase(p)%getFloat('twin_d',defaultVal=0.0_pReal)
prm%twinE = config_phase(p)%getFloat('twin_e',defaultVal=0.0_pReal)
prm%aTolResistance = config_phase(p)%getFloat('atol_resistance',defaultVal=1.0_pReal)
prm%aTolShear = config_phase(p)%getFloat('atol_shear',defaultVal=1.0e-6_pReal)
prm%aTolTwinfrac = config_phase(p)%getFloat('atol_twinfrac',defaultVal=1.0e-6_pReal)