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

small corrections on changes related to non-schmid systems

This commit is contained in:
Martin Diehl 2013-01-22 10:04:15 +00:00
parent e691ff1ea1
commit 2b319f02e6
2 changed files with 5 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
code/constitutive_nonlocal.f90
View File

@ -860,7 +860,7 @@ do i = 1,maxNinstance
+ 3.0_pReal*constitutive_nonlocal_Cslip_66(4,4,i) ) ! (C11iso-C12iso)/2 with C11iso=(3*C11+2*C12+4*C44)/5 and C12iso=(C11+4*C12-2*C44)/5
constitutive_nonlocal_nu(i) = ( constitutive_nonlocal_Cslip_66(1,1,i) + 4.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) &
- 2.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) ) &
/ ( 4.0_pReal*constitutive_nonlocal_Cslip_66(1,1,i) + 6.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) &
/ ( 4.0_pReal*constitutive_nonlocal_Cslip_66(1,1,i) + 6.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) &
+ 2.0_pReal*constitutive_nonlocal_Cslip_66(4,4,i) )
constitutive_nonlocal_Cslip_66(1:6,1:6,i) = math_Mandel3333to66(math_Voigt66to3333(constitutive_nonlocal_Cslip_66(1:6,1:6,i)))
constitutive_nonlocal_Cslip_3333(1:3,1:3,1:3,1:3,i) = math_Voigt66to3333(constitutive_nonlocal_Cslip_66(1:6,1:6,i))

12
code/lattice.f90
View File

@ -752,19 +752,15 @@ pure function lattice_symmetrizeC66(structName,C66)
select case(structName(1:3))
case ('iso')
forall(k=1_pInt:3_pInt)
forall(j=1_pInt:3_pInt)
lattice_symmetrizeC66(k,j) = C66(1,2)
end forall
forall(j=1_pInt:3_pInt) lattice_symmetrizeC66(k,j) = C66(1,2)
lattice_symmetrizeC66(k,k) = C66(1,1)
lattice_symmetrizeC66(k+3,k+3) = 0.5_pReal*(C66(1,1)-C66(1,2))
end forall
case ('fcc','bcc')
forall(k=1_pInt:3_pInt)
forall(j=1_pInt:3_pInt)
lattice_symmetrizeC66(k,j) = C66(1,2)
lattice_symmetrizeC66(k,k) = C66(1,1)
lattice_symmetrizeC66(k+3_pInt,k+3_pInt) = C66(4,4)
end forall
forall(j=1_pInt:3_pInt) lattice_symmetrizeC66(k,j) = C66(1,2)
lattice_symmetrizeC66(k,k) = C66(1,1)
lattice_symmetrizeC66(k+3_pInt,k+3_pInt) = C66(4,4)
end forall
case ('hex')
lattice_symmetrizeC66(1,1) = C66(1,1)