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

missing factor of 2 in damage driving force

This commit is contained in:
Pratheek Shanthraj 2014-11-07 12:10:10 +00:00
parent dc0745800f
commit c4e54793fc
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
code/damage_anisoDuctile.f90
View File

@ -244,7 +244,7 @@ subroutine damage_anisoDuctile_init(fileUnit)
if (any(numerics_integrator == 5_pInt)) & if (any(numerics_integrator == 5_pInt)) &
allocate(damageState(phase)%RKCK45dotState (6,sizeDotState,NofMyPhase),source=0.0_pReal) 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) call damage_anisoDuctile_aTolState(phase,instance)
endif endif
@ -254,14 +254,14 @@ end subroutine damage_anisoDuctile_init
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief sets the relevant state values for a given instance of this damage !> @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: & use material, only: &
damageState damageState
use math, only: & use math, only: &
math_I3 math_I3
implicit none 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 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) = & damageState(phase)%dotState(index+1,constituent) = &
-damage_anisoDuctile_sdot_0(instance)* & -damage_anisoDuctile_sdot_0(instance)* &
max(0.0_pReal, & 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) - & accumulatedSlip(index)/damage_anisoDuctile_critAccShear(f,instance) - &
nonlocalFactor)**damage_anisoDuctile_N(instance) nonlocalFactor)**damage_anisoDuctile_N(instance)
index = index + 1_pInt index = index + 1_pInt