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

modify fomular in pheno+ to remove dampping

This commit is contained in:
zhangc43 2016-04-25 15:18:46 -04:00
parent bc4f04a1c5
commit b168317a17
1 changed files with 4 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/plastic_phenoplus.f90
View File

@ -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