A modified slip-twin interaction matrices for FCC and BCC are introduced.

The shear twin value for BCC is corrected.
An introduction of a third parameter at the interactionSlipTwin in the material.config for dislotwin materials is essential.
This commit is contained in:
Nader Zaafarani 2012-11-14 11:13:16 +00:00
parent b988caa6ba
commit b0464ad176
1 changed files with 80 additions and 97 deletions

View File

@ -178,19 +178,23 @@ module lattice
integer(pInt), dimension(lattice_fcc_Nslip,lattice_fcc_Ntwin), target, private :: & integer(pInt), dimension(lattice_fcc_Nslip,lattice_fcc_Ntwin), target, private :: &
lattice_fcc_interactionSlipTwin = transpose(reshape(int( [& lattice_fcc_interactionSlipTwin = transpose(reshape(int( [&
1,1,1,2,2,1,1,2,2,2,1,2, & ! ---> twin 1,1,1,3,3,3,2,2,2,3,3,3, & ! ---> twin
1,1,1,2,2,1,1,2,2,2,1,2, & ! | 1,1,1,3,3,3,3,3,3,2,2,2, & ! |
1,1,1,2,2,1,1,2,2,2,1,2, & ! | 1,1,1,2,2,2,3,3,3,3,3,3, & ! |
2,2,1,1,1,1,2,1,2,1,2,2, & ! v slip 3,3,3,1,1,1,3,3,3,2,2,2, & ! v slip
2,2,1,1,1,1,2,1,2,1,2,2, & 3,3,3,1,1,1,2,2,2,3,3,3, &
2,2,1,1,1,1,2,1,2,1,2,2, & 2,2,2,1,1,1,3,3,3,3,3,3, &
1,2,2,2,1,2,1,1,1,2,2,1, & 2,2,2,3,3,3,1,1,1,3,3,3, &
1,2,2,2,1,2,1,1,1,2,2,1, & 3,3,3,2,2,2,1,1,1,3,3,3, &
1,2,2,2,1,2,1,1,1,2,2,1, & 3,3,3,3,3,3,1,1,1,2,2,2, &
2,1,2,1,2,2,2,2,1,1,1,1, & 3,3,3,2,2,2,3,3,3,1,1,1, &
2,1,2,1,2,2,2,2,1,1,1,1, & 2,2,2,3,3,3,3,3,3,1,1,1, &
2,1,2,1,2,2,2,2,1,1,1,1 & 3,3,3,3,3,3,2,2,2,1,1,1 &
],pInt),[lattice_fcc_Ntwin,lattice_fcc_Nslip])) ],pInt),[lattice_fcc_Ntwin,lattice_fcc_Nslip]))
!< Interaction types
!< 1 --- coplanar interaction
!< 2 --- colinenar interaction
!< 3 --- hardened interaction
integer(pInt), dimension(lattice_fcc_Ntwin,lattice_fcc_Nslip), target, private :: & integer(pInt), dimension(lattice_fcc_Ntwin,lattice_fcc_Nslip), target, private :: &
lattice_fcc_interactionTwinSlip = 0_pInt lattice_fcc_interactionTwinSlip = 0_pInt
@ -259,30 +263,30 @@ module lattice
-1, 1, 1, 1,-1, 2, & -1, 1, 1, 1,-1, 2, &
1, 1, 1, 1, 1,-2, & 1, 1, 1, 1, 1,-2, &
! Slip system <111>{123} meaningful sorting ? ! Slip system <111>{123} meaningful sorting ?
1, 1,-1, 1, 2, 3, & ! 1, 1,-1, 1, 2, 3, &
1,-1, 1, -1, 2, 3, & ! 1,-1, 1, -1, 2, 3, &
-1, 1, 1, 1,-2, 3, & !-1, 1, 1, 1,-2, 3, &
1, 1, 1, 1, 2,-3, & ! 1, 1, 1, 1, 2,-3, &
1,-1, 1, 1, 3, 2, & ! 1,-1, 1, 1, 3, 2, &
1, 1,-1, -1, 3, 2, & ! 1, 1,-1, -1, 3, 2, &
1, 1, 1, 1,-3, 2, & ! 1, 1, 1, 1,-3, 2, &
-1, 1, 1, 1, 3,-2, & !-1, 1, 1, 1, 3,-2, &
1, 1,-1, 2, 1, 3, & ! 1, 1,-1, 2, 1, 3, &
1,-1, 1, -2, 1, 3, & ! 1,-1, 1, -2, 1, 3, &
-1, 1, 1, 2,-1, 3, & !-1, 1, 1, 2,-1, 3, &
1, 1, 1, 2, 1,-3, & ! 1, 1, 1, 2, 1,-3, &
1,-1, 1, 2, 3, 1, & ! 1,-1, 1, 2, 3, 1, &
1, 1,-1, -2, 3, 1, & ! 1, 1,-1, -2, 3, 1, &
1, 1, 1, 2,-3, 1, & ! 1, 1, 1, 2,-3, 1, &
-1, 1, 1, 2, 3,-1, & !-1, 1, 1, 2, 3,-1, &
-1, 1, 1, 3, 1, 2, & !-1, 1, 1, 3, 1, 2, &
1, 1, 1, -3, 1, 2, & ! 1, 1, 1, -3, 1, 2, &
1, 1,-1, 3,-1, 2, & ! 1, 1,-1, 3,-1, 2, &
1,-1, 1, 3, 1,-2, & ! 1,-1, 1, 3, 1,-2, &
-1, 1, 1, 3, 2, 1, & !-1, 1, 1, 3, 2, 1, &
1, 1, 1, -3, 2, 1, & ! 1, 1, 1, -3, 2, 1, &
1, 1,-1, 3,-2, 1, & ! 1, 1,-1, 3,-2, 1, &
1,-1, 1, 3, 2,-1 & ! 1,-1, 1, 3, 2,-1 &
],pReal),[ 3_pInt + 3_pInt ,lattice_bcc_Nslip]) ],pReal),[ 3_pInt + 3_pInt ,lattice_bcc_Nslip])
! twin system <111>{112} ! twin system <111>{112}
@ -306,18 +310,18 @@ module lattice
real(pReal), dimension(lattice_bcc_Ntwin), parameter, private :: & real(pReal), dimension(lattice_bcc_Ntwin), parameter, private :: &
lattice_bcc_shearTwin = reshape([& lattice_bcc_shearTwin = reshape([&
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal, & 0.7071067812_pReal, &
0.123_pReal & 0.7071067812_pReal &
],[lattice_bcc_Ntwin]) ! Twin system {111}<112> just a dummy ],[lattice_bcc_Ntwin]) ! Twin system {111}<112> just a dummy
!> slip--slip interactions for BCC structures (2) !> slip--slip interactions for BCC structures (2)
@ -376,56 +380,35 @@ module lattice
!> slip--twin interactions for BCC structures (2) MISSING: not implemented yet !> slip--twin interactions for BCC structures (2) MISSING: not implemented yet
integer(pInt), dimension(lattice_bcc_Nslip,lattice_bcc_Ntwin), target, private :: & integer(pInt), dimension(lattice_bcc_Nslip,lattice_bcc_Ntwin), target, private :: &
lattice_bcc_interactionSlipTwin = transpose(reshape(int( [& lattice_bcc_interactionSlipTwin = transpose(reshape(int( [&
0,0,0,0,0,0,0,0,0,0,0,0, & ! ---> twin 3,3,3,2,2,3,3,3,3,2,3,3, & ! ---> twin
0,0,0,0,0,0,0,0,0,0,0,0, & ! | 3,3,2,3,3,2,3,3,2,3,3,3, & ! |
0,0,0,0,0,0,0,0,0,0,0,0, & ! | 3,2,3,3,3,3,2,3,3,3,3,2, & ! |
0,0,0,0,0,0,0,0,0,0,0,0, & ! v slip 2,3,3,3,3,3,3,2,3,3,2,3, & ! v slip
0,0,0,0,0,0,0,0,0,0,0,0, & 2,3,3,3,3,3,3,2,3,3,2,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,2,3,3,2,3,3,2,3,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,2,3,3,3,3,2,3,3,3,3,2, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,3,2,2,3,3,3,3,2,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 2,3,3,3,3,3,3,2,3,3,2,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,3,2,2,3,3,3,3,2,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,2,3,3,3,3,2,3,3,3,3,2, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,2,3,3,2,3,3,2,3,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 1,3,3,3,3,3,3,2,3,3,2,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,1,3,3,3,3,2,3,3,3,3,2, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,1,3,3,2,3,3,2,3,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,3,1,2,3,3,3,3,2,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,3,2,1,3,3,3,3,2,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,2,3,3,1,3,3,2,3,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,2,3,3,3,3,1,3,3,3,3,2, &
0,0,0,0,0,0,0,0,0,0,0,0, & 2,3,3,3,3,3,3,1,3,3,2,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,2,3,3,2,3,3,1,3,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,3,3,2,2,3,3,3,3,1,3,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 2,3,3,3,3,3,3,2,3,3,1,3, &
0,0,0,0,0,0,0,0,0,0,0,0, & 3,2,3,3,3,3,2,3,3,3,3,1 &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0, &
0,0,0,0,0,0,0,0,0,0,0,0 &
],pInt),[lattice_bcc_Ntwin,lattice_bcc_Nslip])) ],pInt),[lattice_bcc_Ntwin,lattice_bcc_Nslip]))
!< Interaction types
!< 1 --- coplanar interaction
!< 2 --- colinenar interaction
!< 3 --- hardened interaction
!>twin--slip interactions for BCC structures (2) MISSING: not implemented yet !>twin--slip interactions for BCC structures (2) MISSING: not implemented yet
integer(pInt), dimension(lattice_bcc_Ntwin,lattice_bcc_Nslip), target, private :: & integer(pInt), dimension(lattice_bcc_Ntwin,lattice_bcc_Nslip), target, private :: &