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

only Hooke as model 

and this model should take care of the elastic constants, not lattice
This commit is contained in:
Martin Diehl 2021-05-24 17:19:38 +02:00
parent 0d0bc188eb
commit fe281f4592
1 changed files with 21 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

33
src/phase_mechanical_elastic.f90
View File

@ -1,12 +1,11 @@
submodule(phase:mechanical) elastic submodule(phase:mechanical) elastic
enum, bind(c); enumerator :: & type :: tParameters
ELASTICITY_UNDEFINED_ID, & real(pReal), dimension(6,6) :: &
ELASTICITY_HOOKE_ID C66 !< Elastic constants in Voig notation
end enum end type tParameters
integer(kind(ELASTICITY_UNDEFINED_ID)), dimension(:), allocatable :: & type(tParameters), allocatable, dimension(:) :: param
phase_elasticity !< elasticity of each phase
contains contains
@ -23,19 +22,27 @@ module subroutine elastic_init(phases)
mech, & mech, &
elastic elastic
print'(/,a)', ' <<<+- phase:mechanical:elastic init -+>>>' print'(/,a)', ' <<<+- phase:mechanical:elastic init -+>>>'
allocate(phase_elasticity(phases%length), source = ELASTICITY_undefined_ID) allocate(param(phases%length))
do ph = 1, phases%length do ph = 1, phases%length
phase => phases%get(ph) phase => phases%get(ph)
mech => phase%get('mechanical') mech => phase%get('mechanical')
elastic => mech%get('elastic') elastic => mech%get('elastic')
if(IO_lc(elastic%get_asString('type')) == 'hooke') then ! accept small letter h for the moment if (elastic%get_asString('type') /= 'Hooke') call IO_error(200,ext_msg=elastic%get_asString('type'))
phase_elasticity(ph) = ELASTICITY_HOOKE_ID
else associate(prm => param(ph))
call IO_error(200,ext_msg=elastic%get_asString('type')) prm%C66(1,1) = elastic%get_asFloat('C_11')
endif prm%C66(1,2) = elastic%get_asFloat('C_12')
prm%C66(4,4) = elastic%get_asFloat('C_44')
prm%C66(1,3) = elastic%get_asFloat('C_13',defaultVal=0.0_pReal)
prm%C66(2,3) = elastic%get_asFloat('C_23',defaultVal=0.0_pReal)
prm%C66(3,3) = elastic%get_asFloat('C_33',defaultVal=0.0_pReal)
prm%C66(6,6) = elastic%get_asFloat('C_66',defaultVal=0.0_pReal)
end associate
enddo enddo
end subroutine elastic_init end subroutine elastic_init
@ -63,9 +70,9 @@ module subroutine phase_hooke_SandItsTangents(S, dS_dFe, dS_dFi, &
real(pReal), dimension(3,3) :: E real(pReal), dimension(3,3) :: E
real(pReal), dimension(3,3,3,3) :: C real(pReal), dimension(3,3,3,3) :: C
integer :: & integer :: &
d, & !< counter in degradation loop
i, j i, j
C = math_66toSym3333(phase_homogenizedC(ph,en)) C = math_66toSym3333(phase_homogenizedC(ph,en))
C = phase_damage_C(C,ph,en) C = phase_damage_C(C,ph,en)