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

bugfix: change of behavior 

negative values for the resolved stress do not make sense.
The paper does not take this into account (eq (14), Cereceda et. al
2016). According to my understanding, only the non-thermal contributions should be
substracted, so abs(tau_pos)/abs(tau_neg) would not be sufficient.
This commit is contained in:
Martin Diehl 2021-10-30 23:03:27 +02:00
parent e20b705f54
commit f5fe0b9dca
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/phase_mechanical_plastic_dislotungsten.f90
View File

@ -491,7 +491,7 @@ pure subroutine kinetics(Mp,T,ph,en, &
t_n = prm%b_sl*exp(BoltzmannRatio*(1.0_pReal-StressRatio_p) ** prm%q) & t_n = prm%b_sl*exp(BoltzmannRatio*(1.0_pReal-StressRatio_p) ** prm%q) &
/ (prm%omega*effectiveLength) / (prm%omega*effectiveLength)
t_k = effectiveLength * prm%B /(2.0_pReal*prm%b_sl*tau_pos) t_k = effectiveLength * prm%B /(2.0_pReal*prm%b_sl*tau_eff) ! corrected eq. (14)
dot_gamma_pos = dot_gamma_0 * sign(prm%h/(t_n + t_k),tau_pos) * 0.5_pReal dot_gamma_pos = dot_gamma_0 * sign(prm%h/(t_n + t_k),tau_pos) * 0.5_pReal
else where significantPositiveTau else where significantPositiveTau
@ -521,7 +521,7 @@ pure subroutine kinetics(Mp,T,ph,en, &
t_n = prm%b_sl*exp(BoltzmannRatio*(1.0_pReal-StressRatio_p) ** prm%q) & t_n = prm%b_sl*exp(BoltzmannRatio*(1.0_pReal-StressRatio_p) ** prm%q) &
/ (prm%omega*effectiveLength) / (prm%omega*effectiveLength)
t_k = effectiveLength * prm%B /(2.0_pReal*prm%b_sl*tau_pos) t_k = effectiveLength * prm%B /(2.0_pReal*prm%b_sl*tau_eff) ! corrected eq. (14)
dot_gamma_neg = dot_gamma_0 * sign(prm%h/(t_n + t_k),tau_neg) * 0.5_pReal dot_gamma_neg = dot_gamma_0 * sign(prm%h/(t_n + t_k),tau_neg) * 0.5_pReal
else where significantNegativeTau else where significantNegativeTau