From 3d5a85e4afa99959352af8d95f1b20fdce06bed8 Mon Sep 17 00:00:00 2001 From: Luv Sharma Date: Mon, 6 Jul 2015 08:47:40 +0000 Subject: [PATCH] corrected implementation of deltaState for isoBrittle damage --- code/source_damage_isoBrittle.f90 | 10 ++++------ 1 file changed, 4 insertions(+), 6 deletions(-) diff --git a/code/source_damage_isoBrittle.f90 b/code/source_damage_isoBrittle.f90 index 1e918aa8f..25795d893 100644 --- a/code/source_damage_isoBrittle.f90 +++ b/code/source_damage_isoBrittle.f90 @@ -277,7 +277,6 @@ subroutine source_damage_isoBrittle_deltaState(C, Fe, ipc, ip, el) instance = source_damage_isoBrittle_instance(phase) sourceOffset = source_damage_isoBrittle_offset(phase) - stiffness = C do mech = 1_pInt, phase_NstiffnessDegradations(phase) select case(phase_stiffnessDegradation(mech,phase)) case (STIFFNESS_DEGRADATION_porosity_ID) @@ -290,15 +289,14 @@ subroutine source_damage_isoBrittle_deltaState(C, Fe, ipc, ip, el) strainenergy = 2.0_pReal*sum(strain*math_mul66x6(stiffness,strain))/ & source_damage_isoBrittle_critStrainEnergy(instance) - if (strainenergy > sourceState(phase)%p(sourceOffset)%subState0(1,constituent)) then + if (strainenergy > sourceState(phase)%p(sourceOffset)%state(1,constituent)) then sourceState(phase)%p(sourceOffset)%deltaState(1,constituent) = & strainenergy - sourceState(phase)%p(sourceOffset)%state(1,constituent) - else + else sourceState(phase)%p(sourceOffset)%deltaState(1,constituent) = & - sourceState(phase)%p(sourceOffset)%subState0(1,constituent) - & - sourceState(phase)%p(sourceOffset)%state(1,constituent) + 0.0_pReal endif - + end subroutine source_damage_isoBrittle_deltaState !--------------------------------------------------------------------------------------------------