From 1f77594ca0ae39e02fd659e197a96ddcccb473c2 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Wed, 21 Aug 2013 13:54:36 +0000 Subject: [PATCH] dislocation velocity was not correctly stored to state in last version --- code/constitutive_nonlocal.f90 | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index 723a164f8..a6af4acc8 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -1829,7 +1829,6 @@ call constitutive_nonlocal_kinetics(v(1:ns,1), dv_dtau(1:ns,1), dv_dtauNS(1:ns,1 v(1:ns,2) = v(1:ns,1) dv_dtau(1:ns,2) = dv_dtau(1:ns,1) dv_dtauNS(1:ns,2) = dv_dtauNS(1:ns,1) -state%p(iV(1:ns,2,myInstance)) = v(1:ns,1) !screws if (lattice_NnonSchmid(myStructure) == 0_pInt) then ! no non-Schmid contributions @@ -1837,25 +1836,28 @@ if (lattice_NnonSchmid(myStructure) == 0_pInt) then v(1:ns,t) = v(1:ns,1) dv_dtau(1:ns,t) = dv_dtau(1:ns,1) dv_dtauNS(1:ns,t) = dv_dtauNS(1:ns,1) - state%p(iV(1:ns,t,myInstance)) = v(1:ns,1) endforall else ! take non-Schmid contributions into account do t = 3_pInt,4_pInt call constitutive_nonlocal_kinetics(v(1:ns,t), dv_dtau(1:ns,t), dv_dtauNS(1:ns,t), & tau(1:ns), tauNS(1:ns,t), tauThreshold(1:ns), & 2_pInt , Temperature, g, ip, el) - state%p(iV(1:ns,t,myInstance)) = v(1:ns,t) enddo endif +!*** store velocity in state + +forall (t = 1_pInt:4_pInt) & + state%p(iV(1:ns,t,myInstance)) = v(1:ns,t) + + !*** Bauschinger effect forall (s = 1_pInt:ns, t = 5_pInt:8_pInt, rhoSgl(s,t) * v(s,t-4_pInt) < 0.0_pReal) & rhoSgl(s,t-4_pInt) = rhoSgl(s,t-4_pInt) + abs(rhoSgl(s,t)) - !*** Calculation of Lp and its tangent gdotTotal = sum(rhoSgl(1:ns,1:4) * v, 2) * burgers(1:ns,myInstance)