From ec377a6e8ee9d710568f67c88852af4314b320c0 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Wed, 21 Aug 2013 08:39:43 +0000 Subject: [PATCH] removed obsolete "dead zone" scaling --- code/config/material.config | 1 - code/constitutive_nonlocal.f90 | 20 ++++---------------- 2 files changed, 4 insertions(+), 17 deletions(-) diff --git a/code/config/material.config b/code/config/material.config index 54c2813fb..4eb32ebf8 100644 --- a/code/config/material.config +++ b/code/config/material.config @@ -325,7 +325,6 @@ significantRho 1e6 # minimum dislocation density c #significantN 0.1 # minimum dislocation number per ip considered relevant aTol_density 1e4 # absolute tolerance for dislocation density in m/m**3 aTol_shear 1e-20 # absolute tolerance for plasgtic shear -deadZone 0 # switch for the modification of the shearrate in presence of dead dislocations randomMultiplication 0 # switch for probabilistic extension of multiplication rate diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index 4cd006d5e..a7e51df87 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -194,7 +194,6 @@ nonSchmidCoeff logical, dimension(:), allocatable, private :: & shortRangeStressCorrection, & !< flag indicating the use of the short range stress correction by a excess density gradient term -deadZoneScaling, & probabilisticMultiplication public :: & @@ -351,7 +350,6 @@ allocate(rhoSglRandomBinning(maxNinstance)) allocate(surfaceTransmissivity(maxNinstance)) allocate(grainboundaryTransmissivity(maxNinstance)) allocate(shortRangeStressCorrection(maxNinstance)) -allocate(deadZoneScaling(maxNinstance)) allocate(probabilisticMultiplication(maxNinstance)) allocate(CFLfactor(maxNinstance)) allocate(fEdgeMultiplication(maxNinstance)) @@ -388,7 +386,6 @@ fEdgeMultiplication = 0.0_pReal linetensionEffect = 0.0_pReal edgeJogFactor = 1.0_pReal shortRangeStressCorrection = .false. -deadZoneScaling = .false. probabilisticMultiplication = .false. allocate(rhoSglEdgePos0(lattice_maxNslipFamily,maxNinstance)) @@ -581,8 +578,6 @@ do while (trim(line) /= '#EOF#') do f = 1_pInt, lattice_maxNnonSchmid nonSchmidCoeff(f,i) = IO_floatValue(line,positions,1_pInt+f) enddo - case('deadzonescaling','deadzone','deadscaling') - deadZoneScaling(i) = IO_floatValue(line,positions,2_pInt) > 0.0_pReal case('probabilisticmultiplication','randomsources','randommultiplication','discretesources') probabilisticMultiplication(i) = IO_floatValue(line,positions,2_pInt) > 0.0_pReal case default @@ -1771,8 +1766,7 @@ real(pReal), dimension(totalNslip(phase_plasticityInstance(material_phase(g,ip,e real(pReal), dimension(totalNslip(phase_plasticityInstance(material_phase(g,ip,el)))) :: & gdotTotal, & !< shear rate tauBack, & !< back stress from dislocation gradients on same slip system - tauThreshold, & !< threshold shear stress - deadZoneSize + tauThreshold !< threshold shear stress !*** initialize local variables @@ -1851,15 +1845,9 @@ forall (s = 1_pInt:ns, t = 5_pInt:8_pInt, rhoSgl(s,t) * v(s,t-4_pInt) < 0.0_pRea !*** Calculation of gdot and its tangent -deadZoneSize = 0.0_pReal -if (deadZoneScaling(myInstance)) then - forall(s = 1_pInt:ns, sum(abs(rhoSgl(s,1:8))) > 0.0_pReal) & - deadZoneSize(s) = maxval(abs(rhoSgl(s,5:8)) / (rhoSgl(s,1:4) + abs(rhoSgl(s,5:8)))) -endif -gdotTotal = sum(rhoSgl(1:ns,1:4) * v, 2) * burgers(1:ns,myInstance) * (1.0_pReal - deadZoneSize) -do t = 1_pInt,4_pInt - dgdot_dtau(:,t) = rhoSgl(1:ns,t) * dv_dtau(1:ns,t) * burgers(1:ns,myInstance) * (1.0_pReal - deadZoneSize) -enddo +gdotTotal = sum(rhoSgl(1:ns,1:4) * v, 2) * burgers(1:ns,myInstance) +forall (t = 1_pInt:4_pInt) & + dgdot_dtau(1:ns,t) = rhoSgl(1:ns,t) * dv_dtau(1:ns,t) * burgers(1:ns,myInstance) !*** Calculation of Lp and its tangent