From 47a748831fad8397cb96a4cc38d5602f3ca41509 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 5 Oct 2018 07:26:13 +0200 Subject: [PATCH 1/2] Cannot write out resolved stress including non schmid Would need to write out tau_pos and tau_neg separately, so limit to resolved stress in the sense of Schmid --- src/plastic_phenopowerlaw.f90 | 10 +--------- 1 file changed, 1 insertion(+), 9 deletions(-) diff --git a/src/plastic_phenopowerlaw.f90 b/src/plastic_phenopowerlaw.f90 index 76e7dd4b4..edba9a71a 100644 --- a/src/plastic_phenopowerlaw.f90 +++ b/src/plastic_phenopowerlaw.f90 @@ -769,8 +769,6 @@ function plastic_phenopowerlaw_postResults(Mp,instance,of) result(postResults) integer(pInt) :: & o,c,i,j - real(pReal) :: & - tau_slip_pos, tau_slip_neg real(pReal), dimension(param(instance)%totalNslip) :: & gdot_slip_pos,gdot_slip_neg @@ -797,13 +795,7 @@ function plastic_phenopowerlaw_postResults(Mp,instance,of) result(postResults) c = c + prm%totalNslip case (resolvedstress_slip_ID) do i = 1_pInt, prm%totalNslip - tau_slip_pos = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i)) - tau_slip_neg = tau_slip_pos - !do j = 1,size(prm%nonSchmidCoeff) - ! tau_slip_pos = tau_slip_pos + math_mul33xx33(S,prm%nonSchmid_pos(1:3,1:3,j,i)) - ! tau_slip_neg = tau_slip_neg + math_mul33xx33(S,prm%nonSchmid_neg(1:3,1:3,j,i)) - !enddo - postResults(c+i) = 0.5_pReal*(tau_slip_pos+tau_slip_neg) + postResults(c+i) = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i)) enddo c = c + prm%totalNslip From 081ca482e96b6e6fc61613cc9f1d16bdd5044006 Mon Sep 17 00:00:00 2001 From: Test User Date: Sat, 6 Oct 2018 02:47:57 +0200 Subject: [PATCH 2/2] [skip ci] updated version information after successful test of v2.0.2-707-g47a74883 --- VERSION | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/VERSION b/VERSION index 79ee48b68..1bfab081d 100644 --- a/VERSION +++ b/VERSION @@ -1 +1 @@ -v2.0.2-638-gd5cf7d21 +v2.0.2-707-g47a74883