From b168317a17cd7148da46eaecd8e2f2cd69f6cbdd Mon Sep 17 00:00:00 2001 From: zhangc43 Date: Mon, 25 Apr 2016 15:18:46 -0400 Subject: [PATCH] modify fomular in pheno+ to remove dampping --- src/plastic_phenoplus.f90 | 10 ++++------ 1 file changed, 4 insertions(+), 6 deletions(-) diff --git a/src/plastic_phenoplus.f90 b/src/plastic_phenoplus.f90 index b12d476ab..1646f5f41 100644 --- a/src/plastic_phenoplus.f90 +++ b/src/plastic_phenoplus.f90 @@ -979,13 +979,11 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar ne_mprimes(n) = maxval(m_primes) ENDIF ENDDO LOOPMYNEIGHBORS - !*******check if one of the neighbor already can provide a kick for this slip system - IF ( maxval(ne_mprimes) > mprime_cut ) THEN - plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal + 0.5_pReal*maxval(ne_mprimes) - ELSE - plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal - 0.1_pReal*maxval(ne_mprimes) - ENDIF + + plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal & + + 0.15_pReal*(1.0_pReal + ERF(10.0_pReal*(maxval(ne_mprimes)-0.8_pReal))) ENDDO LOOPMYSLIP + ENDIF ENDIF