From a011b0a4a26e743710d97d3789791ce1453da839 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Thu, 16 Jan 2014 09:53:45 +0000 Subject: [PATCH] corrected tangent of plastic velocity gradient dLp_dT for slip (twin volume fraction was missing) --- code/constitutive_dislotwin.f90 | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/code/constitutive_dislotwin.f90 b/code/constitutive_dislotwin.f90 index f63a6657b..b461fb767 100644 --- a/code/constitutive_dislotwin.f90 +++ b/code/constitutive_dislotwin.f90 @@ -1224,8 +1224,9 @@ subroutine constitutive_dislotwin_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v,Temperat constitutive_dislotwin_v0PerSlipSystem(j,matID) !* Shear rates due to slip - gdot_slip(j) = DotGamma0*exp(-BoltzmannRatio*(1-StressRatio_p)**constitutive_dislotwin_q(matID))*& - sign(1.0_pReal,tau_slip(j)) + gdot_slip(j) = (1.0_pReal - sumf) * DotGamma0 & + * exp(-BoltzmannRatio*(1-StressRatio_p) ** constitutive_dislotwin_q(matID)) & + * sign(1.0_pReal,tau_slip(j)) !* Derivatives of shear rates dgdot_dtauslip(j) = & @@ -1234,7 +1235,7 @@ subroutine constitutive_dislotwin_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v,Temperat StressRatio_pminus1*(1-StressRatio_p)**(constitutive_dislotwin_q(matID)-1.0_pReal) !* Plastic velocity gradient for dislocation glide - Lp = Lp + (1.0_pReal - sumf)*gdot_slip(j)*lattice_Sslip(:,:,1,index_myFamily+i,structID) + Lp = Lp + gdot_slip(j)*lattice_Sslip(:,:,1,index_myFamily+i,structID) !* Calculation of the tangent of Lp forall (k=1_pInt:3_pInt,l=1_pInt:3_pInt,m=1_pInt:3_pInt,n=1_pInt:3_pInt) &