From 71c6179ee5ed6db5610857ce8568ea89a7b1cbd9 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Tue, 8 Jun 2010 09:34:23 +0000 Subject: [PATCH] dislocation stress had wrong sign in some components --- code/constitutive_nonlocal.f90 | 18 +++++++++--------- 1 file changed, 9 insertions(+), 9 deletions(-) diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index 837fa6c4c..ec7dc8c4d 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -984,21 +984,21 @@ do n = 1,FE_NipNeighbors(mesh_element(2,el)) do i = 1,2 r2_2 = lambda(i)**2.0_pReal + r1_2 - sigma(1,1,i) = z * lambda(i) / dsqrt(r2_2) * ( - 2.0_pReal * constitutive_nonlocal_nu(myInstance) / r1_2 & + sigma(1,1,i) = - z * lambda(i) / dsqrt(r2_2) * ( 2.0_pReal * constitutive_nonlocal_nu(myInstance) / r1_2 & * ( 1.0_pReal + a**2.0_pReal/r1_2 + 0.5_pReal*a**2.0_pReal/r2_2 ) & - + 1.0_pReal / r2_2 ) + - 1.0_pReal / r2_2 ) - sigma(2,2,i) = z * lambda(i) / ( r1_2 * dsqrt(r2_2) ) & + sigma(2,2,i) = - z * lambda(i) / ( r1_2 * dsqrt(r2_2) ) & * ( 1.0_pReal + 2.0_pReal*(y**2.0_pReal+a**2.0_pReal)/r1_2 + (y**2.0_pReal+a**2.0_pReal)/r2_2 ) - sigma(3,3,i) = z * lambda(i) / ( r1_2 * dsqrt(r2_2) ) & + sigma(3,3,i) = + z * lambda(i) / ( r1_2 * dsqrt(r2_2) ) & * ( 1.0_pReal - 2.0_pReal*(z**2.0_pReal+a**2.0_pReal)/r1_2 - (z**2.0_pReal+a**2.0_pReal)/r2_2 ) - sigma(1,2,i) = y * z / ( r2_2 * dsqrt(r2_2) ) + sigma(1,2,i) = + y * z / ( r2_2 * dsqrt(r2_2) ) - sigma(2,3,i) = y * lambda(i) / ( r1_2 * dsqrt(r2_2) ) * ( 1.0_pReal - 2.0_pReal*z**2.0_pReal/r1_2 - z**2.0_pReal/r2_2 ) + sigma(2,3,i) = + y * lambda(i) / ( r1_2 * dsqrt(r2_2) ) * ( 1.0_pReal - 2.0_pReal*z**2.0_pReal/r1_2 - z**2.0_pReal/r2_2 ) - sigma(1,3,i) = 1.0_pReal / dsqrt(r2_2) * ( constitutive_nonlocal_nu(myInstance) - z**2.0_pReal/r2_2 & + sigma(1,3,i) = + 1.0_pReal / dsqrt(r2_2) * ( constitutive_nonlocal_nu(myInstance) - z**2.0_pReal/r2_2 & - 0.5_pReal*(1.0_pReal-constitutive_nonlocal_nu(myInstance))*a**2.0_pReal/r2_2 ) enddo @@ -1015,10 +1015,10 @@ do n = 1,FE_NipNeighbors(mesh_element(2,el)) do i = 1,2 r2_2 = lambda(i)**2.0_pReal + r1_2 - sigma(1,2,i) = z * (1.0_pReal - constitutive_nonlocal_nu(myInstance)) * lambda(i) / ( r1_2 * dsqrt(r2_2) ) & + sigma(1,2,i) = - z * (1.0_pReal - constitutive_nonlocal_nu(myInstance)) * lambda(i) / ( r1_2 * dsqrt(r2_2) ) & * ( 1.0_pReal + a**2.0_pReal/r1_2 + 0.5_pReal*a**2.0_pReal/r2_2 ) - sigma(2,3,i) = - x * (1.0_pReal - constitutive_nonlocal_nu(myInstance)) * lambda(i) / ( r1_2 * dsqrt(r2_2) ) & + sigma(2,3,i) = + x * (1.0_pReal - constitutive_nonlocal_nu(myInstance)) * lambda(i) / ( r1_2 * dsqrt(r2_2) ) & * ( 1.0_pReal + a**2.0_pReal/r1_2 + 0.5_pReal*a**2.0_pReal/r2_2 ) enddo