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

first calculation of Tstar still used hard coded Hooke, now usese constitutive_TandItsTangent(...)

This commit is contained in:
Franz Roters 2012-06-18 07:26:34 +00:00
parent 8033a9d0da
commit 380f40561b
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
code/crystallite.f90
View File

@ -546,7 +546,8 @@ real(pReal), dimension(3,3,3,3) :: dSdFe, &
dFedFdot, & dFedFdot, &
dSdF, & dSdF, &
dSdFdot, & dSdFdot, &
dFp_invdFdot dFp_invdFdot, &
junk2
real(pReal) :: counter real(pReal) :: counter
! --+>> INITIALIZE TO STARTING CONDITION <<+-- ! --+>> INITIALIZE TO STARTING CONDITION <<+--
@ -726,8 +727,7 @@ enddo
#endif #endif
invFp = math_inv33(crystallite_partionedFp0(1:3,1:3,g,i,e)) invFp = math_inv33(crystallite_partionedFp0(1:3,1:3,g,i,e))
Fe_guess = math_mul33x33(crystallite_partionedF(1:3,1:3,g,i,e), invFp) Fe_guess = math_mul33x33(crystallite_partionedF(1:3,1:3,g,i,e), invFp)
Tstar = math_Mandel6to33( math_mul66x6( 0.5_pReal*constitutive_homogenizedC(g,i,e), & call constitutive_TandItsTangent(Tstar, junk2, Fe_guess,g,i,e)
math_Mandel33to6( math_mul33x33(transpose(Fe_guess),Fe_guess) - math_I3 ) ) )
crystallite_P(1:3,1:3,g,i,e) = math_mul33x33(Fe_guess,math_mul33x33(Tstar,transpose(invFp))) crystallite_P(1:3,1:3,g,i,e) = math_mul33x33(Fe_guess,math_mul33x33(Tstar,transpose(invFp)))
endif endif
#ifndef _OPENMP #ifndef _OPENMP