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

divide by 3.0 instead of 3 as it is _pReal

This commit is contained in:
Franz Roters 2009-07-27 09:23:39 +00:00
parent 53cff7458c
commit 1cd0b74b08
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
code/constitutive_j2.f90
View File

@ -482,7 +482,7 @@ pure function constitutive_j2_postResults(Tstar_v, Temperature, dt, state, g, ip
matID = phase_constitutionInstance(material_phase(g,ip,el))
! calculate deviatoric part of 2nd Piola-Kirchhoff stress and its norm
Tstar_dev_v(1:3) = Tstar_v(1:3) - sum(Tstar_v(1:3))/3_pReal
Tstar_dev_v(1:3) = Tstar_v(1:3) - sum(Tstar_v(1:3))/3.0_pReal
Tstar_dev_v(4:6) = Tstar_v(4:6)
norm_Tstar_dev = dsqrt(math_mul6x6(Tstar_dev_v,Tstar_dev_v))