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

following paper

This commit is contained in:
Martin Diehl 2022-04-01 12:40:00 +02:00
parent 4f646f03ee
commit fb633798af
5 changed files with 6 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
PRIVATE

@ -1 +1 @@
Subproject commit 459326e9840c843ade72b04cf28e50889c9779f1 Subproject commit dd7f4efd64ebadac74cab84f713a738c71d8b737

2
examples/config/phase/damage/anisobrittle_cubic.yaml
View File

@ -9,5 +9,5 @@ s_crit: [0.006666]
dot_o: 1.e-3 dot_o: 1.e-3
q: 20 q: 20
D_11: 1.0 l_c: 1.0
mu: 0.001 mu: 0.001

4
examples/config/phase/damage/isobrittle_generic.yaml
View File

@ -2,6 +2,6 @@ type: isobrittle
output: [f_phi] output: [f_phi]
W_crit: 1400000.0 G: 1400000.0
D_11: 1.0 l_c: 1.0
mu: 0.001 mu: 0.001

4
src/phase_damage.f90
View File

@ -105,9 +105,7 @@ module subroutine damage_init
damage_active = .true. damage_active = .true.
source => sources%get(1) source => sources%get(1)
param(ph)%mu = source%get_asFloat('mu') param(ph)%mu = source%get_asFloat('mu')
param(ph)%D(1,1) = source%get_asFloat('D_11') param(ph)%D = math_I3 * source%get_asFloat('l_c')**2
if (any(phase_lattice(ph) == ['hP','tI'])) param(ph)%D(3,3) = source%get_asFloat('D_33')
param(ph)%D = lattice_symmetrize_33(param(ph)%D,phase_lattice(ph))
end if end if
end do end do

2
src/phase_damage_isobrittle.f90
View File

@ -63,7 +63,7 @@ module function isobrittle_init() result(mySources)
associate(prm => param(ph), dlt => deltaState(ph), stt => state(ph)) associate(prm => param(ph), dlt => deltaState(ph), stt => state(ph))
src => sources%get(1) src => sources%get(1)
prm%W_crit = src%get_asFloat('W_crit') prm%W_crit = src%get_asFloat('G')/src%get_asFloat('l_c')
#if defined (__GFORTRAN__) #if defined (__GFORTRAN__)
prm%output = output_as1dString(src) prm%output = output_as1dString(src)