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

correction of hardening coefficients also enabled for bcc 

for fcc: now applied to all hardening coefficients, not only hirth, lomer, and glissile
This commit is contained in:
Christoph Kords 2013-08-14 11:44:30 +00:00
parent 659c498aea
commit e909b825cc
1 changed files with 5 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
code/constitutive_nonlocal.f90
View File

@ -1306,12 +1306,10 @@ integer(pInt) neighboring_el, & ! element number o
neighboring_ns, & ! total number of active slip systems at neighboring material point neighboring_ns, & ! total number of active slip systems at neighboring material point
c, & ! index of dilsocation character (edge, screw) c, & ! index of dilsocation character (edge, screw)
s, & ! slip system index 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-) t, & ! index of dilsocation type (e+, e-, s+, s-, used e+, used e-, used s+, used s-)
dir, & dir, &
n, & n, &
nRealNeighbors, & ! number of really existing neighbors nRealNeighbors ! number of really existing neighbors
interactionCoefficient
integer(pInt), dimension(2) :: neighbor integer(pInt), dimension(2) :: neighbor
real(pReal) detFe, & real(pReal) detFe, &
detFp, & detFp, &
@ -1385,26 +1383,19 @@ forall (s = 1_pInt:ns) &
!*** calculate the threshold shear stress for dislocation slip !*** 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 = 0.0_pReal
myInteractionMatrix(1:ns,1:ns) = interactionMatrixSlipSlip(1:ns,1:ns,instance) 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 do s = 1_pInt,ns
myRhoForest = max(rhoForest(s),significantRho(instance)) myRhoForest = max(rhoForest(s),significantRho(instance))
correction = ( 1.0_pReal - linetensionEffect(instance) & correction = ( 1.0_pReal - linetensionEffect(instance) &
+ linetensionEffect(instance) & + linetensionEffect(instance) &
* log(0.35_pReal * burgers(s,instance) * sqrt(myRhoForest)) & * log(0.35_pReal * burgers(s,instance) * sqrt(myRhoForest)) &
/ log(0.35_pReal * burgers(s,instance) * 1e6_pReal)) ** 2.0_pReal / log(0.35_pReal * burgers(s,instance) * 1e6_pReal)) ** 2.0_pReal
do s2 = 1_pInt,ns myInteractionMatrix(s,1:ns) = correction * myInteractionMatrix(s,1: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
enddo enddo
endif endif
forall (s = 1_pInt:ns) & forall (s = 1_pInt:ns) &