From e909b825cc67b0e3b9d582e581860baeb18d6599 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Wed, 14 Aug 2013 11:44:30 +0000 Subject: [PATCH] correction of hardening coefficients also enabled for bcc for fcc: now applied to all hardening coefficients, not only hirth, lomer, and glissile --- code/constitutive_nonlocal.f90 | 19 +++++-------------- 1 file changed, 5 insertions(+), 14 deletions(-) diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index 32e3e3275..05b95d844 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -1306,12 +1306,10 @@ integer(pInt) neighboring_el, & ! element number o neighboring_ns, & ! total number of active slip systems at neighboring material point c, & ! index of dilsocation character (edge, screw) s, & ! slip system index - s2, & ! slip system index t, & ! index of dilsocation type (e+, e-, s+, s-, used e+, used e-, used s+, used s-) dir, & n, & - nRealNeighbors, & ! number of really existing neighbors - interactionCoefficient + nRealNeighbors ! number of really existing neighbors integer(pInt), dimension(2) :: neighbor real(pReal) detFe, & detFp, & @@ -1385,26 +1383,19 @@ forall (s = 1_pInt:ns) & !*** calculate the threshold shear stress for dislocation slip +!*** coefficients are corrected for the line tension effect +!*** (see Kubin,Devincre,Hoc; 2008; Modeling dislocation storage rates and mean free paths in face-centered cubic crystals) myInteractionMatrix = 0.0_pReal myInteractionMatrix(1:ns,1:ns) = interactionMatrixSlipSlip(1:ns,1:ns,instance) -if (latticeStruct == 1_pInt) then ! in case of fcc: coefficients are corrected for the line tension effect (see Kubin,Devincre,Hoc; 2008; Modeling dislocation storage rates and mean free paths in face-centered cubic crystals) +if (latticeStruct < 3_pInt) then ! only fcc and bcc do s = 1_pInt,ns myRhoForest = max(rhoForest(s),significantRho(instance)) correction = ( 1.0_pReal - linetensionEffect(instance) & + linetensionEffect(instance) & * log(0.35_pReal * burgers(s,instance) * sqrt(myRhoForest)) & / log(0.35_pReal * burgers(s,instance) * 1e6_pReal)) ** 2.0_pReal - do s2 = 1_pInt,ns - interactionCoefficient = & - lattice_interactionSlipSlip(slipSystemLattice(s,instance), & - slipSystemLattice(s2,instance), & - latticeStruct) - select case(interactionCoefficient) - case(4_pInt,5_pInt,6_pInt) ! only correct junction forming interactions (4,5,6) - myInteractionMatrix(s,s2) = correction * myInteractionMatrix(s,s2) - endselect - enddo + myInteractionMatrix(s,1:ns) = correction * myInteractionMatrix(s,1:ns) enddo endif forall (s = 1_pInt:ns) &