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

Read in elastic constants for transformed phase

This commit is contained in:
Su Leen Wong 2015-07-07 13:23:56 +00:00
parent 4435244681
commit 66564c1f75
1 changed files with 27 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
code/lattice.f90
View File

@ -1012,12 +1012,14 @@ module lattice
],pReal),[ 3_pInt + 3_pInt,LATTICE_ortho_Ncleavage]) ],pReal),[ 3_pInt + 3_pInt,LATTICE_ortho_Ncleavage])
real(pReal), dimension(:,:,:), allocatable, public, protected :: & real(pReal), dimension(:,:,:), allocatable, public, protected :: &
lattice_C66 lattice_C66, lattice_trans_C66
real(pReal), dimension(:,:,:,:,:), allocatable, public, protected :: & real(pReal), dimension(:,:,:,:,:), allocatable, public, protected :: &
lattice_C3333 lattice_C3333, lattice_trans_C3333
real(pReal), dimension(:), allocatable, public, protected :: & real(pReal), dimension(:), allocatable, public, protected :: &
lattice_mu, & lattice_mu, &
lattice_nu lattice_nu, &
lattice_trans_mu, &
lattice_trans_nu
real(pReal), dimension(:,:,:), allocatable, public, protected :: & real(pReal), dimension(:,:,:), allocatable, public, protected :: &
lattice_thermalConductivity33, & lattice_thermalConductivity33, &
lattice_thermalExpansion33, & lattice_thermalExpansion33, &
@ -1322,6 +1324,8 @@ subroutine lattice_init
allocate(trans_lattice_structure(Nphases),source = LATTICE_undefined_ID) allocate(trans_lattice_structure(Nphases),source = LATTICE_undefined_ID)
allocate(lattice_C66(6,6,Nphases), source=0.0_pReal) allocate(lattice_C66(6,6,Nphases), source=0.0_pReal)
allocate(lattice_C3333(3,3,3,3,Nphases), source=0.0_pReal) allocate(lattice_C3333(3,3,3,3,Nphases), source=0.0_pReal)
allocate(lattice_trans_C66(6,6,Nphases), source=0.0_pReal)
allocate(lattice_trans_C3333(3,3,3,3,Nphases), source=0.0_pReal)
allocate(lattice_thermalConductivity33 (3,3,Nphases), source=0.0_pReal) allocate(lattice_thermalConductivity33 (3,3,Nphases), source=0.0_pReal)
allocate(lattice_thermalExpansion33 (3,3,Nphases), source=0.0_pReal) allocate(lattice_thermalExpansion33 (3,3,Nphases), source=0.0_pReal)
allocate(lattice_damageDiffusion33 (3,3,Nphases), source=0.0_pReal) allocate(lattice_damageDiffusion33 (3,3,Nphases), source=0.0_pReal)
@ -1345,6 +1349,8 @@ subroutine lattice_init
allocate(lattice_mu(Nphases), source=0.0_pReal) allocate(lattice_mu(Nphases), source=0.0_pReal)
allocate(lattice_nu(Nphases), source=0.0_pReal) allocate(lattice_nu(Nphases), source=0.0_pReal)
allocate(lattice_trans_mu(Nphases), source=0.0_pReal)
allocate(lattice_trans_nu(Nphases), source=0.0_pReal)
allocate(lattice_NnonSchmid(Nphases), source=0_pInt) allocate(lattice_NnonSchmid(Nphases), source=0_pInt)
allocate(lattice_Sslip(3,3,1+2*lattice_maxNnonSchmid,lattice_maxNslip,Nphases),source=0.0_pReal) allocate(lattice_Sslip(3,3,1+2*lattice_maxNnonSchmid,lattice_maxNslip,Nphases),source=0.0_pReal)
@ -1447,6 +1453,24 @@ subroutine lattice_init
lattice_C66(5,5,section) = IO_floatValue(line,positions,2_pInt) lattice_C66(5,5,section) = IO_floatValue(line,positions,2_pInt)
case ('c66') case ('c66')
lattice_C66(6,6,section) = IO_floatValue(line,positions,2_pInt) lattice_C66(6,6,section) = IO_floatValue(line,positions,2_pInt)
case ('c11_trans')
lattice_trans_C66(1,1,section) = IO_floatValue(line,positions,2_pInt)
case ('c12_trans')
lattice_trans_C66(1,2,section) = IO_floatValue(line,positions,2_pInt)
case ('c13_trans')
lattice_trans_C66(1,3,section) = IO_floatValue(line,positions,2_pInt)
case ('c22_trans')
lattice_trans_C66(2,2,section) = IO_floatValue(line,positions,2_pInt)
case ('c23_trans')
lattice_trans_C66(2,3,section) = IO_floatValue(line,positions,2_pInt)
case ('c33_trans')
lattice_trans_C66(3,3,section) = IO_floatValue(line,positions,2_pInt)
case ('c44_trans')
lattice_trans_C66(4,4,section) = IO_floatValue(line,positions,2_pInt)
case ('c55_trans')
lattice_trans_C66(5,5,section) = IO_floatValue(line,positions,2_pInt)
case ('c66_trans')
lattice_trans_C66(6,6,section) = IO_floatValue(line,positions,2_pInt)
case ('covera_ratio','c/a_ratio','c/a') case ('covera_ratio','c/a_ratio','c/a')
CoverA(section) = IO_floatValue(line,positions,2_pInt) CoverA(section) = IO_floatValue(line,positions,2_pInt)
case ('c/a_trans','c/a_martensite','c/a_mart') case ('c/a_trans','c/a_martensite','c/a_mart')