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

updated tests

This commit is contained in:
Martin Diehl 2021-05-24 17:39:09 +02:00
parent f44edb31fc
commit f525999f52
2 changed files with 10 additions and 24 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
PRIVATE

@ -1 +1 @@
Subproject commit 14e754308e81545a231220f9d993fb4571729ced Subproject commit 51fd422131bbef9ea37a8e63429d907473b778cd

32
src/lattice.f90
View File

@ -462,10 +462,8 @@ subroutine lattice_init
lattice_C66(4,4,ph) = elasticity%get_asFloat('C_44') lattice_C66(4,4,ph) = elasticity%get_asFloat('C_44')
lattice_C66(1,3,ph) = elasticity%get_asFloat('C_13',defaultVal=0.0_pReal) lattice_C66(1,3,ph) = elasticity%get_asFloat('C_13',defaultVal=0.0_pReal)
lattice_C66(2,2,ph) = elasticity%get_asFloat('C_22',defaultVal=0.0_pReal)
lattice_C66(2,3,ph) = elasticity%get_asFloat('C_23',defaultVal=0.0_pReal) lattice_C66(2,3,ph) = elasticity%get_asFloat('C_23',defaultVal=0.0_pReal)
lattice_C66(3,3,ph) = elasticity%get_asFloat('C_33',defaultVal=0.0_pReal) lattice_C66(3,3,ph) = elasticity%get_asFloat('C_33',defaultVal=0.0_pReal)
lattice_C66(5,5,ph) = elasticity%get_asFloat('C_55',defaultVal=0.0_pReal)
lattice_C66(6,6,ph) = elasticity%get_asFloat('C_66',defaultVal=0.0_pReal) lattice_C66(6,6,ph) = elasticity%get_asFloat('C_66',defaultVal=0.0_pReal)
select case(phase%get_asString('lattice')) select case(phase%get_asString('lattice'))
@ -1715,34 +1713,22 @@ function applyLatticeSymmetryC66(C66,structure) result(C66_sym)
select case(structure) select case(structure)
case ('cF','cI') case ('cF','cI')
C66_sym(1,1) = C66(1,1) C66_sym(1,1) = C66(1,1); C66_sym(2,2) = C66(1,1); C66_sym(3,3) = C66(1,1)
C66_sym(2,2) = C66(1,1) C66_sym(1,2) = C66(1,2); C66_sym(1,3) = C66(1,2); C66_sym(2,3) = C66(1,2)
C66_sym(3,3) = C66(1,1) C66_sym(4,4) = C66(4,4); C66_sym(5,5) = C66(4,4); C66_sym(6,6) = C66(4,4) ! isotropic C_44 = (C_11-C_12)/2
C66_sym(1,2) = C66(1,2)
C66_sym(1,3) = C66(1,2)
C66_sym(2,3) = C66(1,2)
C66_sym(4,4) = C66(4,4) ! isotropic C_44 = (C_11-C_12)/2
C66_sym(5,5) = C66(4,4)
C66_sym(6,6) = C66(4,4)
case ('hP') case ('hP')
C66_sym(1,1) = C66(1,1) C66_sym(1,1) = C66(1,1); C66_sym(2,2) = C66(1,1)
C66_sym(2,2) = C66(1,1)
C66_sym(3,3) = C66(3,3) C66_sym(3,3) = C66(3,3)
C66_sym(1,2) = C66(1,2) C66_sym(1,2) = C66(1,2)
C66_sym(1,3) = C66(1,3) C66_sym(1,3) = C66(1,3); C66_sym(2,3) = C66(1,3)
C66_sym(2,3) = C66(1,3) C66_sym(4,4) = C66(4,4); C66_sym(5,5) = C66(4,4)
C66_sym(4,4) = C66(4,4)
C66_sym(5,5) = C66(4,4)
C66_sym(6,6) = 0.5_pReal*(C66(1,1)-C66(1,2)) C66_sym(6,6) = 0.5_pReal*(C66(1,1)-C66(1,2))
case ('tI') case ('tI')
C66_sym(1,1) = C66(1,1) C66_sym(1,1) = C66(1,1); C66_sym(2,2) = C66(1,1)
C66_sym(2,2) = C66(1,1)
C66_sym(3,3) = C66(3,3) C66_sym(3,3) = C66(3,3)
C66_sym(1,2) = C66(1,2) C66_sym(1,2) = C66(1,2)
C66_sym(1,3) = C66(1,3) C66_sym(1,3) = C66(1,3); C66_sym(2,3) = C66(1,3)
C66_sym(2,3) = C66(1,3) C66_sym(4,4) = C66(4,4); C66_sym(5,5) = C66(4,4)
C66_sym(4,4) = C66(4,4)
C66_sym(5,5) = C66(4,4)
C66_sym(6,6) = C66(6,6) C66_sym(6,6) = C66(6,6)
case default case default
call IO_error(137,ext_msg='applyLatticeSymmetryC66: '//trim(structure)) call IO_error(137,ext_msg='applyLatticeSymmetryC66: '//trim(structure))