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

prefactors were not used/not mentioned in the paper

This commit is contained in:
Martin Diehl 2019-03-19 06:21:40 +01:00
parent 7d40d13b4c
commit d9e497bdd8
1 changed files with 9 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
src/plastic_dislotwin.f90
View File

@ -50,7 +50,6 @@ module plastic_dislotwin
qShearBand, & !< q-exponent in shear band velocity qShearBand, & !< q-exponent in shear band velocity
CEdgeDipMinDistance, & !< CEdgeDipMinDistance, & !<
Cmfptwin, & !< Cmfptwin, & !<
Cthresholdtwin, & !<
SolidSolutionStrength, & !<strength due to elements in solid solution SolidSolutionStrength, & !<strength due to elements in solid solution
L0_twin, & !< Length of twin nuclei in Burgers vectors L0_twin, & !< Length of twin nuclei in Burgers vectors
L0_trans, & !< Length of trans nuclei in Burgers vectors L0_trans, & !< Length of trans nuclei in Burgers vectors
@ -67,7 +66,6 @@ module plastic_dislotwin
aTolTransFrac, & !< absolute tolerance for integration of trans volume fraction aTolTransFrac, & !< absolute tolerance for integration of trans volume fraction
deltaG, & !< Free energy difference between austensite and martensite deltaG, & !< Free energy difference between austensite and martensite
Cmfptrans, & !< Cmfptrans, & !<
Cthresholdtrans, & !<
transStackHeight !< Stack height of hex nucleus transStackHeight !< Stack height of hex nucleus
real(pReal), dimension(:), allocatable :: & real(pReal), dimension(:), allocatable :: &
rho_mob_0, & !< initial unipolar dislocation density per slip system rho_mob_0, & !< initial unipolar dislocation density per slip system
@ -141,7 +139,7 @@ module plastic_dislotwin
mfp_slip, & mfp_slip, &
mfp_twin, & mfp_twin, &
mfp_trans, & mfp_trans, &
threshold_stress_slip, & tau_pass, &
threshold_stress_twin, & threshold_stress_twin, &
threshold_stress_trans, & threshold_stress_trans, &
twinVolume, & twinVolume, &
@ -353,7 +351,6 @@ subroutine plastic_dislotwin_init
prm%xc_twin = config%getFloat('xc_twin') prm%xc_twin = config%getFloat('xc_twin')
prm%L0_twin = config%getFloat('l0_twin') prm%L0_twin = config%getFloat('l0_twin')
prm%Cthresholdtwin = config%getFloat('cthresholdtwin', defaultVal=0.0_pReal)
prm%Cmfptwin = config%getFloat('cmfptwin', defaultVal=0.0_pReal) ! ToDo: How to handle that??? prm%Cmfptwin = config%getFloat('cmfptwin', defaultVal=0.0_pReal) ! ToDo: How to handle that???
prm%shear_twin = lattice_characteristicShear_Twin(prm%N_tw,config%getString('lattice_structure'),& prm%shear_twin = lattice_characteristicShear_Twin(prm%N_tw,config%getString('lattice_structure'),&
@ -386,7 +383,6 @@ subroutine plastic_dislotwin_init
prm%burgers_trans = config%getFloats('transburgers') prm%burgers_trans = config%getFloats('transburgers')
prm%burgers_trans = math_expand(prm%burgers_trans,prm%Ntrans) prm%burgers_trans = math_expand(prm%burgers_trans,prm%Ntrans)
prm%Cthresholdtrans = config%getFloat('cthresholdtrans', defaultVal=0.0_pReal) ! ToDo: How to handle that???
prm%transStackHeight = config%getFloat('transstackheight', defaultVal=0.0_pReal) ! ToDo: How to handle that??? prm%transStackHeight = config%getFloat('transstackheight', defaultVal=0.0_pReal) ! ToDo: How to handle that???
prm%Cmfptrans = config%getFloat('cmfptrans', defaultVal=0.0_pReal) ! ToDo: How to handle that??? prm%Cmfptrans = config%getFloat('cmfptrans', defaultVal=0.0_pReal) ! ToDo: How to handle that???
prm%deltaG = config%getFloat('deltag') prm%deltaG = config%getFloat('deltag')
@ -591,7 +587,7 @@ subroutine plastic_dislotwin_init
allocate(dst%invLambdaSlipTwin (prm%totalNslip, NipcMyPhase),source=0.0_pReal) allocate(dst%invLambdaSlipTwin (prm%totalNslip, NipcMyPhase),source=0.0_pReal)
allocate(dst%invLambdaSlipTrans (prm%totalNslip, NipcMyPhase),source=0.0_pReal) allocate(dst%invLambdaSlipTrans (prm%totalNslip, NipcMyPhase),source=0.0_pReal)
allocate(dst%mfp_slip (prm%totalNslip, NipcMyPhase),source=0.0_pReal) allocate(dst%mfp_slip (prm%totalNslip, NipcMyPhase),source=0.0_pReal)
allocate(dst%threshold_stress_slip (prm%totalNslip, NipcMyPhase),source=0.0_pReal) allocate(dst%tau_pass (prm%totalNslip, NipcMyPhase),source=0.0_pReal)
allocate(dst%invLambdaTwin (prm%totalNtwin, NipcMyPhase),source=0.0_pReal) allocate(dst%invLambdaTwin (prm%totalNtwin, NipcMyPhase),source=0.0_pReal)
allocate(dst%mfp_twin (prm%totalNtwin, NipcMyPhase),source=0.0_pReal) allocate(dst%mfp_twin (prm%totalNtwin, NipcMyPhase),source=0.0_pReal)
@ -715,7 +711,7 @@ subroutine plastic_dislotwin_LpAndItsTangent(Lp,dLp_dMp,Mp,Temperature,instance,
0, 1, 1 & 0, 1, 1 &
],pReal),[ 3,6]) ],pReal),[ 3,6])
associate(prm => param(instance), stt => state(instance), dst => microstructure(instance)) associate(prm => param(instance), stt => state(instance))
f_unrotated = 1.0_pReal & f_unrotated = 1.0_pReal &
- sum(stt%twinFraction(1:prm%totalNtwin,of)) & - sum(stt%twinFraction(1:prm%totalNtwin,of)) &
@ -968,17 +964,17 @@ subroutine plastic_dislotwin_dependentState(temperature,instance,of)
dst%mfp_trans(:,of) = prm%Cmfptrans*prm%GrainSize/(1.0_pReal+prm%GrainSize*dst%invLambdaTrans(:,of)) dst%mfp_trans(:,of) = prm%Cmfptrans*prm%GrainSize/(1.0_pReal+prm%GrainSize*dst%invLambdaTrans(:,of))
!* threshold stress for dislocation motion !* threshold stress for dislocation motion
forall (i = 1:prm%totalNslip) dst%threshold_stress_slip(i,of) = & forall (i = 1:prm%totalNslip) dst%tau_pass(i,of) = &
prm%mu*prm%b_sl(i)*& prm%mu*prm%b_sl(i)*&
sqrt(dot_product(stt%rhoEdge(1:prm%totalNslip,of)+stt%rhoEdgeDip(1:prm%totalNslip,of),& sqrt(dot_product(stt%rhoEdge(1:prm%totalNslip,of)+stt%rhoEdgeDip(1:prm%totalNslip,of),&
prm%h_sl_sl(:,i))) prm%h_sl_sl(:,i)))
!* threshold stress for growing twin/martensite !* threshold stress for growing twin/martensite
if(prm%totalNtwin == prm%totalNslip) & if(prm%totalNtwin == prm%totalNslip) &
dst%threshold_stress_twin(:,of) = prm%Cthresholdtwin* & dst%threshold_stress_twin(:,of) = &
(SFE/(3.0_pReal*prm%b_tw)+ 3.0_pReal*prm%b_tw*prm%mu/(prm%L0_twin*prm%b_sl)) ! slip burgers here correct? (SFE/(3.0_pReal*prm%b_tw)+ 3.0_pReal*prm%b_tw*prm%mu/(prm%L0_twin*prm%b_sl)) ! slip burgers here correct?
if(prm%totalNtrans == prm%totalNslip) & if(prm%totalNtrans == prm%totalNslip) &
dst%threshold_stress_trans(:,of) = prm%Cthresholdtrans* & dst%threshold_stress_trans(:,of) = &
(SFE/(3.0_pReal*prm%burgers_trans) + 3.0_pReal*prm%burgers_trans*prm%mu/& (SFE/(3.0_pReal*prm%burgers_trans) + 3.0_pReal*prm%burgers_trans*prm%mu/&
(prm%L0_trans*prm%b_sl) + prm%transStackHeight*prm%deltaG/ (3.0_pReal*prm%burgers_trans) ) (prm%L0_trans*prm%b_sl) + prm%transStackHeight*prm%deltaG/ (3.0_pReal*prm%burgers_trans) )
@ -1052,7 +1048,7 @@ function plastic_dislotwin_postResults(Mp,Temperature,instance,of) result(postRe
enddo enddo
c = c + prm%totalNslip c = c + prm%totalNslip
case (threshold_stress_slip_ID) case (threshold_stress_slip_ID)
postResults(c+1:c+prm%totalNslip) = dst%threshold_stress_slip(1:prm%totalNslip,of) postResults(c+1:c+prm%totalNslip) = dst%tau_pass(1:prm%totalNslip,of)
c = c + prm%totalNslip c = c + prm%totalNslip
case (f_tw_ID) case (f_tw_ID)
@ -1158,7 +1154,7 @@ pure subroutine kinetics_slip(Mp,Temperature,instance,of, &
tau(i) = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i)) tau(i) = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i))
enddo enddo
tau_eff = abs(tau)-dst%threshold_stress_slip(:,of) tau_eff = abs(tau)-dst%tau_pass(:,of)
significantStress: where(tau_eff > tol_math_check) significantStress: where(tau_eff > tol_math_check)
stressRatio = tau_eff/(prm%SolidSolutionStrength+prm%tau_peierls) stressRatio = tau_eff/(prm%SolidSolutionStrength+prm%tau_peierls)