223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

added some output variables in crystallite: F, Fe, Fp, Ee, P, S

This commit is contained in:
Christoph Kords 2010-05-18 07:57:13 +00:00
parent 15328e51cf
commit d96780967b
2 changed files with 37 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

36
code/crystallite.f90
View File

@ -226,7 +226,7 @@ subroutine crystallite_init(Temperature)
mySize = 3 mySize = 3
case('grainrotation') ! Deviation from initial grain orientation in axis-angle form (angle in degrees) case('grainrotation') ! Deviation from initial grain orientation in axis-angle form (angle in degrees)
mySize = 4 mySize = 4
case('defgrad') case('defgrad','f','fe','fp','ee','p','firstpiola','1stpiola','s','tstar','secondpiola','2ndpiola')
mySize = 9 mySize = 9
case default case default
mySize = 0 mySize = 0
@ -1671,7 +1671,10 @@ function crystallite_postResults(&
pReal pReal
use math, only: math_QuaternionToEuler, & use math, only: math_QuaternionToEuler, &
math_QuaternionToAxisAngle, & math_QuaternionToAxisAngle, &
inDeg math_mul33x33, &
math_I3, &
inDeg, &
math_Mandel6to33
use mesh, only: mesh_element use mesh, only: mesh_element
use material, only: microstructure_crystallite, & use material, only: microstructure_crystallite, &
crystallite_Noutput, & crystallite_Noutput, &
@ -1724,10 +1727,31 @@ function crystallite_postResults(&
math_QuaternionToAxisAngle(crystallite_rotation(1:4,g,i,e)) ! grain rotation away from initial orientation as axis-angle math_QuaternionToAxisAngle(crystallite_rotation(1:4,g,i,e)) ! grain rotation away from initial orientation as axis-angle
crystallite_postResults(c+4) = inDeg * crystallite_postResults(c+4) ! angle in degree crystallite_postResults(c+4) = inDeg * crystallite_postResults(c+4) ! angle in degree
c = c + 4_pInt c = c + 4_pInt
case ('defgrad') case ('defgrad','f')
forall (k=0:2,l=0:2) crystallite_postResults(c+1+k*3+l) = & mySize = 9_pInt
crystallite_partionedF(k+1,l+1,g,i,e) crystallite_postResults(c+1:c+1+mySize) = reshape(crystallite_partionedF(:,:,g,i,e),(/mySize/))
c = c+9_pInt c = c + mySize
case ('fe')
mySize = 9_pInt
crystallite_postResults(c+1:c+1+mySize) = reshape(crystallite_Fe(:,:,g,i,e),(/mySize/))
c = c + mySize
case ('ee')
Ee = 0.5_pReal * (math_mul33x33(transpose(crystallite_Fe(:,:,g,i,e)), crystallite_Fe(:,:,g,i,e)) - math_I3)
mySize = 9_pInt
crystallite_postResults(c+1:c+1+mySize) = reshape(Ee(:,:),(/mySize/))
c = c + mySize
case ('fp')
mySize = 9_pInt
crystallite_postResults(c+1:c+1+mySize) = reshape(crystallite_Fp(:,:,g,i,e),(/mySize/))
c = c + mySize
case ('p','firstpiola','1stpiola')
mySize = 9_pInt
crystallite_postResults(c+1:c+1+mySize) = reshape(crystallite_P(:,:,g,i,e),(/mySize/))
c = c + mySize
case ('s','tstar','secondpiola','2ndpiola')
mySize = 9_pInt
crystallite_postResults(c+1:c+1+mySize) = reshape(math_Mandel6to33(crystallite_Tstar_v(:,g,i,e)),(/mySize/))
c = c + mySize
end select end select
enddo enddo

9
code/material.config
View File

@ -66,8 +66,13 @@ crystallite 1
(output) volume (output) volume
(output) orientation (output) orientation
(output) eulerangles (output) eulerangles
(output) grainrotation # deviation from initial orientation as axis (1-3) and angle in degree (4) (output) grainrotation # deviation from initial orientation as axis (1-3) and angle in degree (4)
(output) defgrad (output) f # deformation gradient tensor; synonyms: "defgrad"
(output) fe # elastic deformation gradient tensor
(output) fp # plastic deformation gradient 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"
#-------------------# #-------------------#
<phase> <phase>