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

last commit did't compile, hope it's working as it should now

This commit is contained in:
Mahesh Balasubramaniam 2013-02-05 20:07:25 +00:00
parent 0346506448
commit 037081831b
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
code/constitutive_j2.f90
View File

@ -404,7 +404,9 @@ pure subroutine constitutive_j2_LpAndItsTangent(Lp, dLp_dTstar_99, Tstar_v, Temp
use prec, only: p_vec use prec, only: p_vec
use math, only: math_mul6x6, & use math, only: math_mul6x6, &
math_Mandel6to33, & math_Mandel6to33, &
math_Plain3333to99 math_Plain3333to99, &
math_deviatoric33, &
math_mul33xx33
use mesh, only: mesh_NcpElems, & use mesh, only: mesh_NcpElems, &
mesh_maxNips mesh_maxNips
use material, only: homogenization_maxNgrains, & use material, only: homogenization_maxNgrains, &
@ -441,7 +443,7 @@ pure subroutine constitutive_j2_LpAndItsTangent(Lp, dLp_dTstar_99, Tstar_v, Temp
! deviatoric part of 2nd Piola-Kirchhoff stress ! deviatoric part of 2nd Piola-Kirchhoff stress
Tstar_dev_33 = math_deviatoric33(math_Mandel6to33(Tstar_v)) Tstar_dev_33 = math_deviatoric33(math_Mandel6to33(Tstar_v))
squarenorm_Tstar_dev = math_mul33x33(Tstar_dev_33,Tstar_dev_33) squarenorm_Tstar_dev = math_mul33xx33(Tstar_dev_33,Tstar_dev_33)
norm_Tstar_dev = sqrt(squarenorm_Tstar_dev) norm_Tstar_dev = sqrt(squarenorm_Tstar_dev)
! Initialization of Lp and dLp_dTstar ! Initialization of Lp and dLp_dTstar