From f5fe0b9dca43ff022f9de1154814821a5821b32d Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sat, 30 Oct 2021 23:03:27 +0200 Subject: [PATCH] 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. --- src/phase_mechanical_plastic_dislotungsten.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/phase_mechanical_plastic_dislotungsten.f90 b/src/phase_mechanical_plastic_dislotungsten.f90 index 6a947cc70..9f565df67 100644 --- a/src/phase_mechanical_plastic_dislotungsten.f90 +++ b/src/phase_mechanical_plastic_dislotungsten.f90 @@ -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) & / (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 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) & / (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 else where significantNegativeTau