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

use parameters from param structure

This commit is contained in:
Martin Diehl 2019-02-13 07:24:48 +01:00
parent 61baa66c38
commit 6a0d739d48
2 changed files with 2 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/source_damage_isoBrittle.f90
View File

@ -25,9 +25,6 @@ module source_damage_isoBrittle
integer(pInt), dimension(:), allocatable, target, public :: &
source_damage_isoBrittle_Noutput !< number of outputs per instance of this damage
real(pReal), dimension(:), allocatable, private :: &
source_damage_isoBrittle_critStrainEnergy
enum, bind(c)
enumerator :: undefined_ID, &
damage_drivingforce_ID
@ -141,8 +138,6 @@ subroutine source_damage_isoBrittle_init(fileUnit)
allocate(source_damage_isoBrittle_outputID(maxval(phase_Noutput),Ninstance), source=undefined_ID)
allocate(source_damage_isoBrittle_Noutput(Ninstance), source=0_pInt)
allocate(source_damage_isoBrittle_critStrainEnergy(Ninstance), source=0.0_pReal)
allocate(param(Ninstance))
do p=1, size(config_phase)
@ -196,9 +191,6 @@ subroutine source_damage_isoBrittle_init(fileUnit)
IO_lc(IO_stringValue(line,chunkPos,2_pInt))
end select
case ('isobrittle_criticalstrainenergy')
source_damage_isoBrittle_critStrainEnergy(instance) = IO_floatValue(line,chunkPos,2_pInt)
end select
endif; endif
enddo parsingFile
@ -275,7 +267,7 @@ subroutine source_damage_isoBrittle_deltaState(C, Fe, ipc, ip, el)
strain = 0.5_pReal*math_Mandel33to6(math_mul33x33(math_transpose33(Fe),Fe)-math_I3)
strainenergy = 2.0_pReal*sum(strain*math_mul66x6(stiffness,strain))/ &
source_damage_isoBrittle_critStrainEnergy(instance)
param(instances)%critStrainEnergy
if (strainenergy > sourceState(phase)%p(sourceOffset)%subState0(1,constituent)) then
sourceState(phase)%p(sourceOffset)%deltaState(1,constituent) = &
strainenergy - sourceState(phase)%p(sourceOffset)%state(1,constituent)

9
src/source_damage_isoDuctile.f90
View File

@ -25,8 +25,6 @@ module source_damage_isoDuctile
integer(pInt), dimension(:), allocatable, target, public :: &
source_damage_isoDuctile_Noutput !< number of outputs per instance of this damage
real(pReal), dimension(:), allocatable, private :: &
source_damage_isoDuctile_critPlasticStrain
enum, bind(c)
enumerator :: undefined_ID, &
@ -142,8 +140,6 @@ subroutine source_damage_isoDuctile_init(fileUnit)
allocate(source_damage_isoDuctile_outputID(maxval(phase_Noutput),Ninstance), source=undefined_ID)
allocate(source_damage_isoDuctile_Noutput(Ninstance), source=0_pInt)
allocate(source_damage_isoDuctile_critPlasticStrain(Ninstance), source=0.0_pReal)
allocate(param(Ninstance))
do p=1, size(config_phase)
@ -197,9 +193,6 @@ subroutine source_damage_isoDuctile_init(fileUnit)
IO_lc(IO_stringValue(line,chunkPos,2_pInt))
end select
case ('isoductile_criticalplasticstrain')
source_damage_isoDuctile_critPlasticStrain(instance) = IO_floatValue(line,chunkPos,2_pInt)
end select
endif; endif
enddo parsingFile
@ -262,7 +255,7 @@ subroutine source_damage_isoDuctile_dotState(ipc, ip, el)
sourceState(phase)%p(sourceOffset)%dotState(1,constituent) = &
sum(plasticState(phase)%slipRate(:,constituent))/ &
((damage(homog)%p(damageOffset))**param(instance)%N)/ &
source_damage_isoDuctile_critPlasticStrain(instance)
param(instance)%critPlasticStrain
end subroutine source_damage_isoDuctile_dotState