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

corrections to crystallite_push33ToRef function

This commit is contained in:
Pratheek Shanthraj 2015-08-03 11:07:19 +00:00
parent 1ec3639cf1
commit a01d7d8d66
1 changed files with 10 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
code/crystallite.f90
View File

@ -3400,19 +3400,25 @@ end function crystallite_stateJump
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
function crystallite_push33ToRef(g,i,e, tensor33) function crystallite_push33ToRef(g,i,e, tensor33)
use math, only: & use math, only: &
math_inv33 math_mul33x33, &
math_inv33, &
math_transpose33, &
math_EulerToR
use material, only: &
material_EulerAngles
implicit none implicit none
real(pReal), dimension(3,3) :: crystallite_push33ToRef real(pReal), dimension(3,3) :: crystallite_push33ToRef
real(pReal), dimension(3,3), intent(in) :: tensor33 real(pReal), dimension(3,3), intent(in) :: tensor33
real(pReal), dimension(3,3) :: invFp real(pReal), dimension(3,3) :: T
integer(pInt), intent(in):: & integer(pInt), intent(in):: &
e, & ! element index e, & ! element index
i, & ! integration point index i, & ! integration point index
g ! grain index g ! grain index
invFp = math_inv33(crystallite_Fp(1:3,1:3,g,i,e)) T = math_mul33x33(math_EulerToR(material_EulerAngles(1:3,g,i,e)), &
crystallite_push33ToRef = matmul(invFp,matmul(tensor33,transpose(invFp))) math_transpose33(math_inv33(crystallite_subF(1:3,1:3,g,i,e))))
crystallite_push33ToRef = math_mul33x33(math_transpose33(T),math_mul33x33(tensor33,T))
end function crystallite_push33ToRef end function crystallite_push33ToRef