From c4e54793fca6fe1f80c7926752f27217e7e063d5 Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Fri, 7 Nov 2014 12:10:10 +0000 Subject: [PATCH] missing factor of 2 in damage driving force --- code/damage_anisoDuctile.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/code/damage_anisoDuctile.f90 b/code/damage_anisoDuctile.f90 index 980fd8266..c442d0b77 100644 --- a/code/damage_anisoDuctile.f90 +++ b/code/damage_anisoDuctile.f90 @@ -244,7 +244,7 @@ subroutine damage_anisoDuctile_init(fileUnit) if (any(numerics_integrator == 5_pInt)) & allocate(damageState(phase)%RKCK45dotState (6,sizeDotState,NofMyPhase),source=0.0_pReal) - call damage_anisoDuctile_stateInit(phase,instance) + call damage_anisoDuctile_stateInit(phase) call damage_anisoDuctile_aTolState(phase,instance) endif @@ -254,14 +254,14 @@ end subroutine damage_anisoDuctile_init !-------------------------------------------------------------------------------------------------- !> @brief sets the relevant state values for a given instance of this damage !-------------------------------------------------------------------------------------------------- -subroutine damage_anisoDuctile_stateInit(phase,instance) +subroutine damage_anisoDuctile_stateInit(phase) use material, only: & damageState use math, only: & math_I3 implicit none - integer(pInt), intent(in) :: phase, instance !< number specifying the phase of the damage + integer(pInt), intent(in) :: phase !< number specifying the phase of the damage real(pReal), dimension(damageState(phase)%sizeState) :: tempState @@ -332,7 +332,7 @@ subroutine damage_anisoDuctile_dotState(nSlip, accumulatedSlip, ipc, ip, el) damageState(phase)%dotState(index+1,constituent) = & -damage_anisoDuctile_sdot_0(instance)* & max(0.0_pReal, & - damageState(phase)%state(index+1,constituent)* & + 2.0_pReal*damageState(phase)%state(index+1,constituent)* & accumulatedSlip(index)/damage_anisoDuctile_critAccShear(f,instance) - & nonlocalFactor)**damage_anisoDuctile_N(instance) index = index + 1_pInt