From 6cd6172c0cc42e7fa5c29f042059078af8260f57 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 1 Feb 2013 15:44:50 +0000 Subject: [PATCH] fixed bug in dotState causing strange hardening for certain parameters --- code/constitutive_phenopowerlaw.f90 | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/code/constitutive_phenopowerlaw.f90 b/code/constitutive_phenopowerlaw.f90 index 7cbbdb816..5cb478fd5 100644 --- a/code/constitutive_phenopowerlaw.f90 +++ b/code/constitutive_phenopowerlaw.f90 @@ -798,9 +798,11 @@ function constitutive_phenopowerlaw_dotState(Tstar_v,Temperature,state,ipc,ip,el j = j+1_pInt left_SlipSlip(j) = 1.0_pReal ! no system-dependent left part left_SlipTwin(j) = 1.0_pReal ! no system-dependent left part - right_SlipSlip(j) = (1.0_pReal-state(ipc,ip,el)%p(j) / & + right_SlipSlip(j) = abs(1.0_pReal-state(ipc,ip,el)%p(j) / & (constitutive_phenopowerlaw_tausat_slip(f,matID)+ssat_offset)) & - **constitutive_phenopowerlaw_a_slip(matID) + **constitutive_phenopowerlaw_a_slip(matID)& + *sign(1.0_pReal,1.0_pReal-state(ipc,ip,el)%p(j) / & + (constitutive_phenopowerlaw_tausat_slip(f,matID)+ssat_offset)) right_TwinSlip(j) = 1.0_pReal ! no system-dependent part !--------------------------------------------------------------------------------------------------