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

some corrections to the damage models

This commit is contained in:
Pratheek Shanthraj 2015-03-06 13:12:32 +00:00
parent bbb5ff6ae9
commit 35f9e91e73
2 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
code/damage_anisoBrittle.f90
View File

@ -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)) localDamage = max(residualStiffness,1.0_pReal/damageState(phase)%state(2,constituent))
damageState(phase)%state(1,constituent) = & damageState(phase)%state(1,constituent) = &
localDamage localDamage + &
(damageState(phase)%subState0(1,constituent) - localDamage)* &
exp(-subdt/lattice_DamageMobility(phase))
end subroutine damage_anisoBrittle_microstructure end subroutine damage_anisoBrittle_microstructure

2
code/damage_phaseField.f90
View File

@ -331,7 +331,7 @@ subroutine damage_phaseField_microstructure(C, Fe, Cv, subdt, ipc, ip, el)
(Cv*damage_phaseField_specificVacancyFormationEnergy(instance) + & (Cv*damage_phaseField_specificVacancyFormationEnergy(instance) + &
sum(abs(stress*strain)))/damage_phaseField_surfaceEnergy(instance) sum(abs(stress*strain)))/damage_phaseField_surfaceEnergy(instance)
damageState(phase)%state(2,constituent) = & 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(1,constituent) = &
damageState(phase)%state(2,constituent) + & damageState(phase)%state(2,constituent) + &