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

fixed bug in elasticity matrix calculation

This commit is contained in:
Pratheek Shanthraj 2013-01-22 15:48:47 +00:00
parent 652a8366fe
commit 50e874f887
6 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
code/constitutive_dislotwin.f90
View File

@ -543,7 +543,7 @@ do i = 1_pInt,maxNinstance
!* Elasticity matrix and shear modulus according to material.config
constitutive_dislotwin_Cslip_66(1:6,1:6,i) = lattice_symmetrizeC66(constitutive_dislotwin_structureName(i),&
constitutive_dislotwin_Cslip_66)
constitutive_dislotwin_Cslip_66(:,:,i))
constitutive_dislotwin_Gmod(i) = &
0.2_pReal*(constitutive_dislotwin_Cslip_66(1,1,i)-constitutive_dislotwin_Cslip_66(1,2,i)) &
+0.3_pReal*constitutive_dislotwin_Cslip_66(4,4,i)

2
code/constitutive_j2.f90
View File

@ -305,7 +305,7 @@ subroutine constitutive_j2_init(myFile)
constitutive_j2_sizeState(i) = 1_pInt
constitutive_j2_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_j2_structureName(i),&
constitutive_j2_Cslip_66)
constitutive_j2_Cslip_66(:,:,i))
constitutive_j2_Cslip_66(1:6,1:6,i) = &
math_Mandel3333to66(math_Voigt66to3333(constitutive_j2_Cslip_66(1:6,1:6,i)))

2
code/constitutive_none.f90
View File

@ -182,7 +182,7 @@ subroutine constitutive_none_init(myFile)
constitutive_none_sizeState(i) = 1_pInt
constitutive_none_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_none_structureName(i),&
constitutive_none_Cslip_66)
constitutive_none_Cslip_66(:,:,i))
constitutive_none_Cslip_66(:,:,i) = &
math_Mandel3333to66(math_Voigt66to3333(constitutive_none_Cslip_66(:,:,i)))

2
code/constitutive_nonlocal.f90
View File

@ -855,7 +855,7 @@ do i = 1,maxNinstance
!*** elasticity matrix and shear modulus according to material.config
constitutive_nonlocal_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_nonlocal_structureName(i),&
constitutive_nonlocal_Cslip_66)
constitutive_nonlocal_Cslip_66(:,:,i))
constitutive_nonlocal_Gmod(i) = 0.2_pReal * ( constitutive_nonlocal_Cslip_66(1,1,i) - constitutive_nonlocal_Cslip_66(1,2,i) &
+ 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
code/constitutive_phenopowerlaw.f90
View File

@ -451,7 +451,7 @@ subroutine constitutive_phenopowerlaw_init(myFile)
myStructure = constitutive_phenopowerlaw_structure(i)
constitutive_phenopowerlaw_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_phenopowerlaw_structureName(i),&
constitutive_phenopowerlaw_Cslip_66)
constitutive_phenopowerlaw_Cslip_66(:,:,i))
! assign elasticity tensor
constitutive_phenopowerlaw_Cslip_66(:,:,i) = &
math_Mandel3333to66(math_Voigt66to3333(constitutive_phenopowerlaw_Cslip_66(:,:,i)))

2
code/constitutive_titanmod.f90
View File

@ -868,7 +868,7 @@ do i = 1_pInt,maxNinstance
!* Elasticity matrix and shear modulus according to material.config
constitutive_titanmod_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_titanmod_structureName(i),&
constitutive_titanmod_Cslip_66)
constitutive_titanmod_Cslip_66(:,:,i))
constitutive_titanmod_Gmod(i) = &
0.2_pReal*(constitutive_titanmod_Cslip_66(1,1,i)-constitutive_titanmod_Cslip_66(1,2,i))&
+ 0.3_pReal*constitutive_titanmod_Cslip_66(4,4,i)