minor change to strain energy calculation

This commit is contained in:
Pratheek Shanthraj 2014-08-14 09:59:16 +00:00
parent c8005247ba
commit ad6467638f
1 changed files with 1 additions and 1 deletions

View File

@ -294,7 +294,7 @@ subroutine damage_local_dotState(Tstar_v, Fe, Lp, ipc, ip, el)
if (pressure < 0.0_pReal) stress = stress - pressure*math_I3
trialDamage = min(1.0_pReal, &
(math_trace33(lattice_surfaceEnergy33(1:3,1:3,phase))/3.0_pReal)/ &
(abs(sum(strain*stress)) + &
(sum(abs(strain*stress)) + &
damageState(phase)%state(1,constituent)))
damageState(phase)%dotState(1,constituent) = &