From 4980f57fd5b06553afdff323fd7310089cd03a52 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Fri, 20 Jan 2012 10:25:35 +0000 Subject: [PATCH] added total Green-Lagrange strain to crystallite outputs --- code/config/material.config | 1 + code/crystallite.f90 | 7 ++++++- 2 files changed, 7 insertions(+), 1 deletion(-) diff --git a/code/config/material.config b/code/config/material.config index 5baaf2f5e..1b6a1ad3a 100644 --- a/code/config/material.config +++ b/code/config/material.config @@ -77,6 +77,7 @@ crystallite 1 (output) f # deformation gradient tensor; synonyms: "defgrad" (output) fe # elastic deformation gradient tensor (output) fp # plastic deformation gradient tensor +(output) e # total strain as Green-Lagrange tensor (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" diff --git a/code/crystallite.f90 b/code/crystallite.f90 index 1384becf7..2019c7a0d 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -276,7 +276,7 @@ do i = 1,material_Ncrystallite mySize = 4 case('eulerangles') ! Bunge (3-1-3) Euler angles mySize = 3 - case('defgrad','f','fe','fp','lp','ee','p','firstpiola','1stpiola','s','tstar','secondpiola','2ndpiola') + case('defgrad','f','fe','fp','lp','e','ee','p','firstpiola','1stpiola','s','tstar','secondpiola','2ndpiola') mySize = 9 case default mySize = 0 @@ -3210,6 +3210,11 @@ function crystallite_postResults(& case ('defgrad','f') mySize = 9_pInt crystallite_postResults(c+1:c+mySize) = reshape(math_transpose3x3(crystallite_partionedF(1:3,1:3,g,i,e)),(/mySize/)) + case ('e') + mySize = 9_pInt + crystallite_postResults(c+1:c+mySize) = 0.5_pReal * reshape((math_mul33x33( & + math_transpose3x3(crystallite_partionedF(1:3,1:3,g,i,e)), & + crystallite_partionedF(1:3,1:3,g,i,e)) - math_I3),(/mySize/)) case ('fe') mySize = 9_pInt crystallite_postResults(c+1:c+mySize) = reshape(math_transpose3x3(crystallite_Fe(1:3,1:3,g,i,e)),(/mySize/))