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

reworked damage. cleaned up homogenisation. removed VI stuff and using projection instead

This commit is contained in:
Pratheek Shanthraj 2014-07-07 17:50:40 +00:00
parent 826f3d4dc7
commit a68b216773
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
code/damage_gradient.f90
View File

@ -262,7 +262,10 @@ subroutine damage_gradient_microstructure(Tstar_v, Fe, ipc, ip, el)
math_Mandel6to33, & math_Mandel6to33, &
math_mul33x33, & math_mul33x33, &
math_transpose33, & math_transpose33, &
math_trace33, &
math_I3 math_I3
use lattice, only: &
lattice_surfaceEnergy33
implicit none implicit none
integer(pInt), intent(in) :: & integer(pInt), intent(in) :: &
@ -284,8 +287,8 @@ subroutine damage_gradient_microstructure(Tstar_v, Fe, ipc, ip, el)
damage = damageState(phase)%state(3,constituent)*damageState(phase)%state(3,constituent) damage = damageState(phase)%state(3,constituent)*damageState(phase)%state(3,constituent)
damageState(phase)%state(2,constituent) = & damageState(phase)%state(2,constituent) = &
0.5_pReal*sum(math_Mandel6to33(Tstar_v/damage)*(math_mul33x33(math_transpose33(Fe),Fe)-math_I3)) - & (0.5_pReal*sum(math_Mandel6to33(Tstar_v/damage)*(math_mul33x33(math_transpose33(Fe),Fe)-math_I3)) - &
0.5_pReal*damageState(phase)%state(1,constituent) 0.5_pReal*damageState(phase)%state(1,constituent))/math_trace33(lattice_surfaceEnergy33(1:3,1:3,phase))/3.0_pReal
end subroutine damage_gradient_microstructure end subroutine damage_gradient_microstructure