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

Renamed variables to denote fcc to bcc transformation

This commit is contained in:
Su Leen Wong 2015-06-22 08:33:48 +00:00
parent 5114e0b43e
commit ad17c40d92
2 changed files with 21 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
code/lattice.f90
View File

@ -209,7 +209,7 @@ module lattice
],pInt),[lattice_fcc_Ntwin,lattice_fcc_Ntwin],order=[2,1]) !< Twin--twin interaction types for fcc ],pInt),[lattice_fcc_Ntwin,lattice_fcc_Ntwin],order=[2,1]) !< Twin--twin interaction types for fcc
real(pReal), dimension(4,LATTICE_fcc_Ntrans), parameter, private :: & real(pReal), dimension(4,LATTICE_fcc_Ntrans), parameter, private :: &
LATTICE_fcc_systemTrans = reshape([& LATTICE_fccTobcc_systemTrans = reshape([&
0.0, 1.0, 0.0, 10.26, & ! Pitsch OR (Ma & Hartmaier 2014, Table 3) 0.0, 1.0, 0.0, 10.26, & ! Pitsch OR (Ma & Hartmaier 2014, Table 3)
0.0, 1.0, 0.0, -10.26, & 0.0, 1.0, 0.0, -10.26, &
0.0, 0.0, 1.0, 10.26, & 0.0, 0.0, 1.0, 10.26, &
@ -225,7 +225,7 @@ module lattice
],[ 4_pInt,LATTICE_fcc_Ntrans]) ],[ 4_pInt,LATTICE_fcc_Ntrans])
integer(pInt), dimension(9,LATTICE_fcc_Ntrans), parameter, private :: & integer(pInt), dimension(9,LATTICE_fcc_Ntrans), parameter, private :: &
LATTICE_fcc_bainVariant = reshape(int( [& LATTICE_fccTobcc_bainVariant = reshape(int( [&
1, 0, 0, 0, 1, 0, 0, 0, 1, & ! Pitsch OR (Ma & Hartmaier 2014, Table 3) 1, 0, 0, 0, 1, 0, 0, 0, 1, & ! Pitsch OR (Ma & Hartmaier 2014, Table 3)
1, 0, 0, 0, 1, 0, 0, 0, 1, & 1, 0, 0, 0, 1, 0, 0, 0, 1, &
1, 0, 0, 0, 1, 0, 0, 0, 1, & 1, 0, 0, 0, 1, 0, 0, 0, 1, &
@ -241,7 +241,7 @@ module lattice
],pInt),[ 9_pInt, LATTICE_fcc_Ntrans]) ],pInt),[ 9_pInt, LATTICE_fcc_Ntrans])
real(pReal), dimension(4,LATTICE_fcc_Ntrans), parameter, private :: & real(pReal), dimension(4,LATTICE_fcc_Ntrans), parameter, private :: &
LATTICE_fcc_bainRot = reshape([& LATTICE_fccTobcc_bainRot = reshape([&
1.0, 0.0, 0.0, 45.0, & ! Rotate fcc austensite to bain variant 1.0, 0.0, 0.0, 45.0, & ! Rotate fcc austensite to bain variant
1.0, 0.0, 0.0, 45.0, & 1.0, 0.0, 0.0, 45.0, &
1.0, 0.0, 0.0, 45.0, & 1.0, 0.0, 0.0, 45.0, &
@ -257,7 +257,7 @@ module lattice
],[ 4_pInt,LATTICE_fcc_Ntrans]) ],[ 4_pInt,LATTICE_fcc_Ntrans])
real(pReal), dimension(LATTICE_fcc_Ntrans,LATTICE_fcc_Ntrans), parameter, private :: & ! Matrix for projection of shear from slip system to fault-band (twin) systems real(pReal), dimension(LATTICE_fcc_Ntrans,LATTICE_fcc_Ntrans), parameter, private :: & ! Matrix for projection of shear from slip system to fault-band (twin) systems
LATTICE_fcc_projectionTrans = reshape(real([& ! For ns = nt = nr LATTICE_fccTobcc_projectionTrans = reshape(real([& ! For ns = nt = nr
0, 1,-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, & 0, 1,-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, &
-1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, & -1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, &
1,-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, & 1,-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, &
@ -273,13 +273,13 @@ module lattice
],pReal),[LATTICE_fcc_Ntrans,LATTICE_fcc_Ntrans],order=[2,1]) ],pReal),[LATTICE_fcc_Ntrans,LATTICE_fcc_Ntrans],order=[2,1])
real(pReal), parameter, private :: & real(pReal), parameter, private :: &
LATTICE_fcc_projectionTransFactor = sqrt(3.0_pReal/4.0_pReal) LATTICE_fccTobcc_projectionTransFactor = sqrt(3.0_pReal/4.0_pReal)
real(pReal), parameter, public :: & real(pReal), parameter, public :: &
LATTICE_fcc_shearCritTrans = 0.0224 LATTICE_fccTobcc_shearCritTrans = 0.0224
integer(pInt), dimension(2_pInt,LATTICE_fcc_Ntrans), parameter, public :: & integer(pInt), dimension(2_pInt,LATTICE_fcc_Ntrans), parameter, public :: &
LATTICE_fcc_transNucleationTwinPair = reshape(int( [& LATTICE_fccTobcc_transNucleationTwinPair = reshape(int( [&
4, 7, & 4, 7, &
1, 10, & 1, 10, &
1, 4, & 1, 4, &
@ -1440,13 +1440,13 @@ subroutine lattice_initializeStructure(myPhase,CoverA,CoverA_trans,a_fcc,a_bcc)
select case(trans_lattice_structure(myPhase)) select case(trans_lattice_structure(myPhase))
case (LATTICE_bcc_ID) ! fcc to bcc transformation case (LATTICE_bcc_ID) ! fcc to bcc transformation
do i = 1_pInt,myNtrans do i = 1_pInt,myNtrans
Rtr(1:3,1:3,i) = math_axisAngleToR(lattice_fcc_systemTrans(1:3,i), & ! Pitsch rotation Rtr(1:3,1:3,i) = math_axisAngleToR(lattice_fccTobcc_systemTrans(1:3,i), & ! Pitsch rotation
lattice_fcc_systemTrans(4,i)*INRAD) lattice_fccTobcc_systemTrans(4,i)*INRAD)
Btr(1:3,1:3,i) = math_axisAngleToR(lattice_fcc_bainRot(1:3,i), & ! Rotation of fcc to Bain coordinate system Btr(1:3,1:3,i) = math_axisAngleToR(lattice_fccTobcc_bainRot(1:3,i), & ! Rotation of fcc to Bain coordinate system
lattice_fcc_bainRot(4,i)*INRAD) lattice_fccTobcc_bainRot(4,i)*INRAD)
xtr(1:3,i) = real(LATTICE_fcc_bainVariant(1:3,i),pReal) xtr(1:3,i) = real(LATTICE_fccTobcc_bainVariant(1:3,i),pReal)
ytr(1:3,i) = real(LATTICE_fcc_bainVariant(4:6,i),pReal) ytr(1:3,i) = real(LATTICE_fccTobcc_bainVariant(4:6,i),pReal)
ztr(1:3,i) = real(LATTICE_fcc_bainVariant(7:9,i),pReal) ztr(1:3,i) = real(LATTICE_fccTobcc_bainVariant(7:9,i),pReal)
Utr(1:3,1:3,i) = 0.0_pReal ! Bain deformation Utr(1:3,1:3,i) = 0.0_pReal ! Bain deformation
if ((a_fcc > 0.0_pReal) .and. (a_bcc > 0.0_pReal)) then if ((a_fcc > 0.0_pReal) .and. (a_bcc > 0.0_pReal)) then
Utr(1:3,1:3,i) = (a_bcc/a_fcc)*math_tensorproduct(xtr(1:3,i), xtr(1:3,i)) + & Utr(1:3,1:3,i) = (a_bcc/a_fcc)*math_tensorproduct(xtr(1:3,i), xtr(1:3,i)) + &
@ -1470,7 +1470,8 @@ subroutine lattice_initializeStructure(myPhase,CoverA,CoverA_trans,a_fcc,a_bcc)
lattice_interactionSlipTwin(1:myNslip,1:myNtwin,myPhase) = lattice_fcc_interactionSlipTwin lattice_interactionSlipTwin(1:myNslip,1:myNtwin,myPhase) = lattice_fcc_interactionSlipTwin
lattice_interactionTwinSlip(1:myNtwin,1:myNslip,myPhase) = lattice_fcc_interactionTwinSlip lattice_interactionTwinSlip(1:myNtwin,1:myNslip,myPhase) = lattice_fcc_interactionTwinSlip
lattice_interactionTwinTwin(1:myNtwin,1:myNtwin,myPhase) = lattice_fcc_interactionTwinTwin lattice_interactionTwinTwin(1:myNtwin,1:myNtwin,myPhase) = lattice_fcc_interactionTwinTwin
lattice_projectionTrans(1:myNtrans,1:myNtrans,myPhase) = LATTICE_fcc_projectionTrans*LATTICE_fcc_projectionTransFactor lattice_projectionTrans(1:myNtrans,1:myNtrans,myPhase) = LATTICE_fccTobcc_projectionTrans*&
LATTICE_fccTobcc_projectionTransFactor
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! bcc ! bcc

10
code/plastic_dislotwin.f90
View File

@ -1687,8 +1687,8 @@ subroutine plastic_dislotwin_dotState(Tstar_v,Temperature,ipc,ip,el)
lattice_mu, & lattice_mu, &
lattice_structure, & lattice_structure, &
lattice_fcc_twinNucleationSlipPair, & lattice_fcc_twinNucleationSlipPair, &
lattice_fcc_transNucleationTwinPair, & lattice_fccTobcc_transNucleationTwinPair, &
lattice_fcc_shearCritTrans, & lattice_fccTobcc_shearCritTrans, &
LATTICE_fcc_ID LATTICE_fcc_ID
implicit none implicit none
@ -1900,8 +1900,8 @@ subroutine plastic_dislotwin_dotState(Tstar_v,Temperature,ipc,ip,el)
!* Probability rate of fault band intersection for strain-induced martensite nucleation !* Probability rate of fault band intersection for strain-induced martensite nucleation
select case(lattice_structure(ph)) select case(lattice_structure(ph))
case (LATTICE_fcc_ID) case (LATTICE_fcc_ID)
a = lattice_fcc_transNucleationTwinPair(1,j) a = lattice_fccTobcc_transNucleationTwinPair(1,j)
b = lattice_fcc_transNucleationTwinPair(2,j) b = lattice_fccTobcc_transNucleationTwinPair(2,j)
sa = sign(1_pInt, a) sa = sign(1_pInt, a)
sb = sign(1_pInt, b) sb = sign(1_pInt, b)
ssa = int(sign(1.0_pReal, shearrate_trans(a)),pInt) ssa = int(sign(1.0_pReal, shearrate_trans(a)),pInt)
@ -1912,7 +1912,7 @@ subroutine plastic_dislotwin_dotState(Tstar_v,Temperature,ipc,ip,el)
else else
probrate_trans(j) = 0.0_pReal probrate_trans(j) = 0.0_pReal
endif endif
probrate_trans(j) = probrate_trans(j)/lattice_fcc_shearCritTrans**2 probrate_trans(j) = probrate_trans(j)/lattice_fccTobcc_shearCritTrans**2
!* Dotstate for strain-induced martensite volume fraction !* Dotstate for strain-induced martensite volume fraction
plasticState(ph)%dotState(3_pInt*ns+2_pInt*nt+nr+j, of) = plastic_dislotwin_Cnuc(instance)*& plasticState(ph)%dotState(3_pInt*ns+2_pInt*nt+nr+j, of) = plastic_dislotwin_Cnuc(instance)*&