From a2037db9c6b0295cf752c199161d4d67460d21a4 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Fri, 12 Feb 2010 12:54:00 +0000 Subject: [PATCH] corrected sign of screw part of dislocation stress --- code/constitutive_nonlocal.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index a3268284a..c1b5b27e0 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -915,11 +915,11 @@ do n = 1,FE_NipNeighbors(mesh_element(2,el)) sigma(1,1) = constitutive_nonlocal_nu(myInstance) * (sigma(2,2) + sigma(3,3)) - sigma(1,2) = gb * neighboring_Nscrew(s) * z / (x**2.0_pReal + z**2.0_pReal) + sigma(1,2) = - gb * neighboring_Nscrew(s) * z / (x**2.0_pReal + z**2.0_pReal) sigma(2,3) = gb * ( neighboring_Nedge(s) / (1.0_pReal-constitutive_nonlocal_nu(myInstance)) & * y * (y**2.0_pReal - z**2.0_pReal) / (y**2.0_pReal + z**2.0_pReal)**2.0_pReal & - - neighboring_Nscrew(s) * x / (x**2.0_pReal + z**2.0_pReal) ) + + neighboring_Nscrew(s) * x / (x**2.0_pReal + z**2.0_pReal) ) sigma(2,1) = sigma(1,2) sigma(3,2) = sigma(2,3)