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

removed obsolete "dead zone" scaling

This commit is contained in:
Christoph Kords 2013-08-21 08:39:43 +00:00
parent f706ba3ff9
commit ec377a6e8e
2 changed files with 4 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
code/config/material.config
View File

@ -325,7 +325,6 @@ significantRho 1e6 # minimum dislocation density c
#significantN 0.1 # minimum dislocation number per ip considered relevant #significantN 0.1 # minimum dislocation number per ip considered relevant
aTol_density 1e4 # absolute tolerance for dislocation density in m/m**3 aTol_density 1e4 # absolute tolerance for dislocation density in m/m**3
aTol_shear 1e-20 # absolute tolerance for plasgtic shear 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 randomMultiplication 0 # switch for probabilistic extension of multiplication rate

20
code/constitutive_nonlocal.f90
View File

@ -194,7 +194,6 @@ nonSchmidCoeff
logical, dimension(:), allocatable, private :: & logical, dimension(:), allocatable, private :: &
shortRangeStressCorrection, & !< flag indicating the use of the short range stress correction by a excess density gradient term shortRangeStressCorrection, & !< flag indicating the use of the short range stress correction by a excess density gradient term
deadZoneScaling, &
probabilisticMultiplication probabilisticMultiplication
public :: & public :: &
@ -351,7 +350,6 @@ allocate(rhoSglRandomBinning(maxNinstance))
allocate(surfaceTransmissivity(maxNinstance)) allocate(surfaceTransmissivity(maxNinstance))
allocate(grainboundaryTransmissivity(maxNinstance)) allocate(grainboundaryTransmissivity(maxNinstance))
allocate(shortRangeStressCorrection(maxNinstance)) allocate(shortRangeStressCorrection(maxNinstance))
allocate(deadZoneScaling(maxNinstance))
allocate(probabilisticMultiplication(maxNinstance)) allocate(probabilisticMultiplication(maxNinstance))
allocate(CFLfactor(maxNinstance)) allocate(CFLfactor(maxNinstance))
allocate(fEdgeMultiplication(maxNinstance)) allocate(fEdgeMultiplication(maxNinstance))
@ -388,7 +386,6 @@ fEdgeMultiplication = 0.0_pReal
linetensionEffect = 0.0_pReal linetensionEffect = 0.0_pReal
edgeJogFactor = 1.0_pReal edgeJogFactor = 1.0_pReal
shortRangeStressCorrection = .false. shortRangeStressCorrection = .false.
deadZoneScaling = .false.
probabilisticMultiplication = .false. probabilisticMultiplication = .false.
allocate(rhoSglEdgePos0(lattice_maxNslipFamily,maxNinstance)) allocate(rhoSglEdgePos0(lattice_maxNslipFamily,maxNinstance))
@ -581,8 +578,6 @@ do while (trim(line) /= '#EOF#')
do f = 1_pInt, lattice_maxNnonSchmid do f = 1_pInt, lattice_maxNnonSchmid
nonSchmidCoeff(f,i) = IO_floatValue(line,positions,1_pInt+f) nonSchmidCoeff(f,i) = IO_floatValue(line,positions,1_pInt+f)
enddo enddo
case('deadzonescaling','deadzone','deadscaling')
deadZoneScaling(i) = IO_floatValue(line,positions,2_pInt) > 0.0_pReal
case('probabilisticmultiplication','randomsources','randommultiplication','discretesources') case('probabilisticmultiplication','randomsources','randommultiplication','discretesources')
probabilisticMultiplication(i) = IO_floatValue(line,positions,2_pInt) > 0.0_pReal probabilisticMultiplication(i) = IO_floatValue(line,positions,2_pInt) > 0.0_pReal
case default 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)))) :: & real(pReal), dimension(totalNslip(phase_plasticityInstance(material_phase(g,ip,el)))) :: &
gdotTotal, & !< shear rate gdotTotal, & !< shear rate
tauBack, & !< back stress from dislocation gradients on same slip system tauBack, & !< back stress from dislocation gradients on same slip system
tauThreshold, & !< threshold shear stress tauThreshold !< threshold shear stress
deadZoneSize
!*** initialize local variables !*** 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 !*** Calculation of gdot and its tangent
deadZoneSize = 0.0_pReal gdotTotal = sum(rhoSgl(1:ns,1:4) * v, 2) * burgers(1:ns,myInstance)
if (deadZoneScaling(myInstance)) then forall (t = 1_pInt:4_pInt) &
forall(s = 1_pInt:ns, sum(abs(rhoSgl(s,1:8))) > 0.0_pReal) & dgdot_dtau(1:ns,t) = rhoSgl(1:ns,t) * dv_dtau(1:ns,t) * burgers(1:ns,myInstance)
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
!*** Calculation of Lp and its tangent !*** Calculation of Lp and its tangent