From 2f503f5cdba7652055e531156a5c964864a0a740 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Fri, 25 Feb 2011 08:15:26 +0000 Subject: [PATCH] added "Lp" as crystallite output --- code/crystallite.f90 | 6 +++++- code/material.config | 1 + 2 files changed, 6 insertions(+), 1 deletion(-) diff --git a/code/crystallite.f90 b/code/crystallite.f90 index 230bfcc20..5f1e9e365 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -261,7 +261,7 @@ do i = 1,material_Ncrystallite mySize = 3 case('grainrotation') ! Deviation from initial grain orientation in axis-angle form (angle in degrees) mySize = 4 - case('defgrad','f','fe','fp','ee','p','firstpiola','1stpiola','s','tstar','secondpiola','2ndpiola') + case('defgrad','f','fe','fp','lp','ee','p','firstpiola','1stpiola','s','tstar','secondpiola','2ndpiola') mySize = 9 case default mySize = 0 @@ -2894,6 +2894,10 @@ function crystallite_postResults(& mySize = 9_pInt crystallite_postResults(c+1:c+1+mySize) = reshape(math_transpose3x3(crystallite_Fp(:,:,g,i,e)),(/mySize/)) c = c + mySize + case ('lp') + mySize = 9_pInt + crystallite_postResults(c+1:c+1+mySize) = reshape(math_transpose3x3(crystallite_Lp(:,:,g,i,e)),(/mySize/)) + c = c + mySize case ('p','firstpiola','1stpiola') mySize = 9_pInt crystallite_postResults(c+1:c+1+mySize) = reshape(math_transpose3x3(crystallite_P(:,:,g,i,e)),(/mySize/)) diff --git a/code/material.config b/code/material.config index a20be0e2e..92e316f48 100644 --- a/code/material.config +++ b/code/material.config @@ -73,6 +73,7 @@ crystallite 1 (output) ee # elastic strain as Green-Lagrange tensor (output) p # first Piola-Kichhoff stress tensor; synonyms: "firstpiola", "1stpiola" (output) s # second Piola-Kichhoff stress tensor; synonyms: "tstar", "secondpiola", "2ndpiola" +(output) lp # plastic velocity gradient tensor #-------------------#