From 35f9e91e735c231fc1ba80168bc9ac0570a5c7db Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Fri, 6 Mar 2015 13:12:32 +0000 Subject: [PATCH] some corrections to the damage models --- code/damage_anisoBrittle.f90 | 4 +++- code/damage_phaseField.f90 | 2 +- 2 files changed, 4 insertions(+), 2 deletions(-) diff --git a/code/damage_anisoBrittle.f90 b/code/damage_anisoBrittle.f90 index 3d728928a..ef94bb6f0 100644 --- a/code/damage_anisoBrittle.f90 +++ b/code/damage_anisoBrittle.f90 @@ -390,7 +390,9 @@ subroutine damage_anisoBrittle_microstructure(Tstar_v, subdt, ipc, ip, el) localDamage = max(residualStiffness,1.0_pReal/damageState(phase)%state(2,constituent)) damageState(phase)%state(1,constituent) = & - localDamage + localDamage + & + (damageState(phase)%subState0(1,constituent) - localDamage)* & + exp(-subdt/lattice_DamageMobility(phase)) end subroutine damage_anisoBrittle_microstructure diff --git a/code/damage_phaseField.f90 b/code/damage_phaseField.f90 index a7ea882e9..47114b2fb 100644 --- a/code/damage_phaseField.f90 +++ b/code/damage_phaseField.f90 @@ -331,7 +331,7 @@ subroutine damage_phaseField_microstructure(C, Fe, Cv, subdt, ipc, ip, el) (Cv*damage_phaseField_specificVacancyFormationEnergy(instance) + & sum(abs(stress*strain)))/damage_phaseField_surfaceEnergy(instance) damageState(phase)%state(2,constituent) = & - (1.0_pReal - Cv)*(1.0_pReal - Cv)/drivingForce + max(residualStiffness,(1.0_pReal - Cv)*(1.0_pReal - Cv)/drivingForce) damageState(phase)%state(1,constituent) = & damageState(phase)%state(2,constituent) + &