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

enhanced comment on dislocation kinetics

This commit is contained in:
Christoph Kords 2011-02-08 12:48:50 +00:00
parent 6731659aa4
commit 17a6c0e59b
1 changed files with 12 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
code/constitutive_nonlocal.f90
View File

@ -1088,8 +1088,16 @@ if (Temperature > 0.0_pReal) then
tau(s) = math_mul6x6(Tstar_v + Tdislocation_v, & tau(s) = math_mul6x6(Tstar_v + Tdislocation_v, &
lattice_Sslip_v(:,constitutive_nonlocal_slipSystemLattice(s,myInstance),myStructure)) lattice_Sslip_v(:,constitutive_nonlocal_slipSystemLattice(s,myInstance),myStructure))
!*** only if the resolved shear stress exceeds the threshold stress, dislocations are able to cut the dislocation forest !*** Only if the resolved shear stress exceeds the threshold stress, dislocations are able to cut the dislocation forest.
!*** the forest can't be overcome by thermal activation !*** In contrast to small atomic obstacles the forest can't be overcome by thermal activation.
!***
!*** mean travel distance
!*** The mean dislocation velocity is calculated as: --------------------------
!*** dwell time + travel time
!***
!*** with : mean travel distance = inverse of the root of forest density
!*** dwell time = inverse of attack frequency times probability of success
!*** travel time = mean travel distance over velocity of sound
tauRel = (abs(tau(s)) - tauThreshold(s)) / constitutive_nonlocal_tauObs(myInstance) tauRel = (abs(tau(s)) - tauThreshold(s)) / constitutive_nonlocal_tauObs(myInstance)
if (tauRel > 0.0_pReal .and. tauRel < 1.0_pReal) then if (tauRel > 0.0_pReal .and. tauRel < 1.0_pReal) then
@ -1107,8 +1115,8 @@ if (Temperature > 0.0_pReal) then
/ ((1.0_pReal - tauRel) * (abs(tau(s)) - tauThreshold(s))) / ((1.0_pReal - tauRel) * (abs(tau(s)) - tauThreshold(s)))
endif endif
!*** if resolved stress exceeds even threshold plus obstacle stress, then allow for maximum travel velocity !*** If resolved stress exceeds threshold plus obstacle stress, the probability for thermal activation is 1.
!*** the tangent is zero, since no dependency of tau !*** The tangent is zero, since no dependency of tau.
elseif (tauRel >= 1.0_pReal) then elseif (tauRel >= 1.0_pReal) then
constitutive_nonlocal_v(s,:,g,ip,el) = sign(constitutive_nonlocal_vs(myInstance), tau(s)) & constitutive_nonlocal_v(s,:,g,ip,el) = sign(constitutive_nonlocal_vs(myInstance), tau(s)) &