From 5ffe78ff1f821d759e63dad1c76e6d7ea70c2e66 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Thu, 28 Mar 2019 09:47:03 +0100 Subject: [PATCH] bugfix: strain rate was missing dot_gamma is per definition positive, abs only confuses the reader --- src/plastic_isotropic.f90 | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/src/plastic_isotropic.f90 b/src/plastic_isotropic.f90 index be1204d75..a98f8613e 100644 --- a/src/plastic_isotropic.f90 +++ b/src/plastic_isotropic.f90 @@ -403,7 +403,7 @@ subroutine plastic_isotropic_dotState(Mp,instance,of) dot_gamma = prm%dot_gamma_0 * (sqrt(1.5_pReal) * norm_Mp /(prm%M*stt%xi(of))) **prm%n - if (abs(dot_gamma) > 1e-12_pReal) then + if (dot_gamma > 1e-12_pReal) then if (dEq0(prm%c_1)) then xi_inf_star = prm%xi_inf else @@ -412,9 +412,10 @@ subroutine plastic_isotropic_dotState(Mp,instance,of) )**(1.0_pReal / prm%c_3) & / prm%c_4 * (dot_gamma / prm%dot_gamma_0)**(1.0_pReal / prm%n) endif - dot%xi(of) = ( prm%h0 + prm%h_ln * log(dot_gamma) ) & - * abs( 1.0_pReal - stt%xi(of)/xi_inf_star )**prm%a & - * sign(1.0_pReal, 1.0_pReal - stt%xi(of)/xi_inf_star) + dot%xi(of) = dot_gamma & + * ( prm%h0 + prm%h_ln * log(dot_gamma) ) & + * abs( 1.0_pReal - stt%xi(of)/xi_inf_star )**prm%a & + * sign(1.0_pReal, 1.0_pReal - stt%xi(of)/xi_inf_star) else dot%xi(of) = 0.0_pReal endif