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

polishing, fixed typo

This commit is contained in:
Martin Diehl 2019-02-13 07:34:16 +01:00
parent aa8d218ce7
commit b3e705e628
1 changed files with 5 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/source_damage_isoBrittle.f90
View File

@ -235,10 +235,9 @@ subroutine source_damage_isoBrittle_deltaState(C, Fe, ipc, ip, el)
phase_NstiffnessDegradations, & phase_NstiffnessDegradations, &
phase_stiffnessDegradation phase_stiffnessDegradation
use math, only : & use math, only : &
math_sym33to6, &
math_mul33x33, & math_mul33x33, &
math_mul66x6, & math_mul66x6, &
math_Mandel33to6, &
math_transpose33, &
math_I3 math_I3
implicit none implicit none
@ -254,7 +253,6 @@ subroutine source_damage_isoBrittle_deltaState(C, Fe, ipc, ip, el)
phase, constituent, instance, sourceOffset, mech phase, constituent, instance, sourceOffset, mech
real(pReal) :: & real(pReal) :: &
strain(6), & strain(6), &
stiffness(6,6), &
strainenergy strainenergy
phase = phaseAt(ipc,ip,el) !< phase ID at ipc,ip,el phase = phaseAt(ipc,ip,el) !< phase ID at ipc,ip,el
@ -263,11 +261,11 @@ subroutine source_damage_isoBrittle_deltaState(C, Fe, ipc, ip, el)
instance = source_damage_isoBrittle_instance(phase) !< instance of damage_isoBrittle source instance = source_damage_isoBrittle_instance(phase) !< instance of damage_isoBrittle source
sourceOffset = source_damage_isoBrittle_offset(phase) sourceOffset = source_damage_isoBrittle_offset(phase)
stiffness = C
strain = 0.5_pReal*math_Mandel33to6(math_mul33x33(math_transpose33(Fe),Fe)-math_I3)
strainenergy = 2.0_pReal*sum(strain*math_mul66x6(stiffness,strain))/ & strain = 0.5_pReal*math_sym33to6(math_mul33x33(transpose(Fe),Fe)-math_I3)
param(instances)%critStrainEnergy
strainenergy = 2.0_pReal*sum(strain*math_mul66x6(C,strain))/ &
param(instance)%critStrainEnergy
if (strainenergy > sourceState(phase)%p(sourceOffset)%subState0(1,constituent)) then if (strainenergy > sourceState(phase)%p(sourceOffset)%subState0(1,constituent)) then
sourceState(phase)%p(sourceOffset)%deltaState(1,constituent) = & sourceState(phase)%p(sourceOffset)%deltaState(1,constituent) = &
strainenergy - sourceState(phase)%p(sourceOffset)%state(1,constituent) strainenergy - sourceState(phase)%p(sourceOffset)%state(1,constituent)