From 1ae252c1ce8cdc16f2af773f4697136f1ce11f24 Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Thu, 30 Sep 2010 17:31:19 +0000 Subject: [PATCH] altered hexagonal slip systems (needs comparison to former set, which is potentially wrong) --- code/lattice.f90 | 48 ++++++++++++++++++++++++------------------------ 1 file changed, 24 insertions(+), 24 deletions(-) diff --git a/code/lattice.f90 b/code/lattice.f90 index 165d94954..5060b4b26 100644 --- a/code/lattice.f90 +++ b/code/lattice.f90 @@ -429,37 +429,37 @@ integer(pInt), allocatable, dimension(:,:,:) :: lattice_interactionSlipSlip, & integer(pInt), parameter :: lattice_hex_Ntwin = 24 ! sum(lattice_hex_NtwinSystem) integer(pInt) :: lattice_hex_Nstructure = 0_pInt - !* sorted by YJ.Ro and Philip + !* sorted by A. Alankar & P. Eisenlohr real(pReal), dimension(4+4,lattice_hex_Nslip), parameter :: lattice_hex_systemSlip = & reshape((/& ! Basal systems <1120>{0001} (independent of c/a-ratio, Bravais notation (4 coordinate base)) 2, -1, -1, 0, 0, 0, 0, 1, & - 1, 1, -2, 0, 0, 0, 0, 1, & -1, 2, -1, 0, 0, 0, 0, 1, & + -1, -1, 2, 0, 0, 0, 0, 1, & ! 1st type prismatic systems <1120>{1010} (independent of c/a-ratio) - 2, -1, -1, 0, 0, -1, 1, 0, & - 1, 1, 2, 0, 1, -1, 0, 0, & - -1, 2, -1, 0, 1, 0, -1, 0, & -! 1st type 1st order pyramidal systems <1120>{1011} - 2, -1, -1, 0, 0, -1, 1, 1, & - 1, 1, -2, 0, 1, -1, 0, 1, & - -1, 2, -1, 0, 1, 0, -1, 1, & - -2, 1, 1, 0, 0, 1, -1, 1, & - -1, -1, 2, 0, -1, 1, 0, 1, & - 1, -2, 1, 0, -1, 0, 1, 1, & + 2, -1, -1, 0, 0, 1, -1, 0, & + -1, 2, -1, 0, -1, 0, 1, 0, & + -1, -1, 2, 0, 1, -1, 0, 0, & +! 1st type 1st order pyramidal systems <1120>{1011} -- plane normals depend on the c/a-ratio + 2, -1, -1, 0, 0, 1, -1, 1, & + 1, 1, -2, 0, -1, 1, 0, 1, & + -1, 2, -1, 0, -1, 0, 1, 1, & + -2, 1, 1, 0, 0, -1, 1, 1, & + -1, -1, 2, 0, 1, -1, 0, 1, & + 1, -2, 1, 0, 1, 0, -1, 1, & ! pyramidal system: c+a slip <2113>{1011} -- plane normals depend on the c/a-ratio - 2, -1, -1, -3, 1, -1, 0, 1, & - 1, 1, -2, -3, 1, 0, -1, 1, & - -1, 2, -1, -3, 1, -1, 0, 1, & - -2, 1, 1, -3, 1, 0, -1, 1, & - -1, -1, 2, -3, 0, 1, -1, 1, & - 1, -2, 1, -3, 0, -1, 1, 1, & - -2, 1, 1, -3, -1, 0, 1, 1, & - -1, -1, 2, -3, 0, -1, 1, 1, & - 1, -2, 1, -3, 0, 1, -1, 1, & - 2, -1, -1, -3, -1, 1, 0, 1, & - 1, 1, -2, -3, -1, 0, 1, 1, & - -1, 2, -1, -3, -1, 1, 0, 1 & + -1, 2, -1, 3, 0, 1, -1, 1, & + 1, 1, -2, 3, 0, 1, -1, 1, & + -2, 1, 1, 3, -1, 1, 0, 1, & + -1, 2, -1, 3, -1, 1, 0, 1, & + -1, -1, 2, 3, -1, 0, 1, 1, & + -2, 1, 1, 3, -1, 0, 1, 1, & + 1, -2, 1, 3, 0, -1, 1, 1, & + -1, -1, 2, 3, 0, -1, 1, 1, & + 2, -1, -1, 3, 1, -1, 0, 1, & + 1, -2, 1, 3, 1, -1, 0, 1, & + 1, 1, -2, 3, 1, 0, -1, 1, & + 2, -1, -1, 3, 1, 0, -1, 1, & /),(/4+4,lattice_hex_Nslip/)) real(pReal), dimension(4+4,lattice_hex_Ntwin), parameter :: lattice_hex_systemTwin = &