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

only mobile dislocations (t=1,2,3,4) are used

This commit is contained in:
Martin Diehl 2023-12-29 22:30:04 +01:00
parent 3b57143662
commit f697182971
No known key found for this signature in database
GPG Key ID: 1FD50837275A0A9B
1 changed files with 10 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
src/phase_mechanical_plastic_nonlocal.f90
View File

@ -1122,11 +1122,10 @@ function rhoDotFlux(timestep,ph,en)
rho, & rho, &
rho_0, & !< dislocation density at beginning of time step rho_0, & !< dislocation density at beginning of time step
rhoDotFlux !< density evolution by flux rhoDotFlux !< density evolution by flux
real(pREAL), dimension(param(ph)%sum_N_sl,8) :: &
rho_0_sgl_nbr, & !< current single dislocation densities of neighboring ip (positive/negative screw and edge without dipoles)
rho_0_sgl !< single dislocation densities of central ip (positive/negative screw and edge without dipoles)
real(pREAL), dimension(param(ph)%sum_N_sl,4) :: & real(pREAL), dimension(param(ph)%sum_N_sl,4) :: &
v, & !< current dislocation glide velocity rho_0_sgl_mob, & !< mobile dislocation densities of neighboring ip (positive/negative screw and edge)
rho_0_sgl_mob_nbr, & !< mobile dislocation densities of neighboring ip (positive/negative screw and edge)
v, & !< dislocation glide velocity
v_0, & v_0, &
v_0_nbr, & !< dislocation glide velocity of enighboring ip v_0_nbr, & !< dislocation glide velocity of enighboring ip
dot_gamma !< shear rates dot_gamma !< shear rates
@ -1160,7 +1159,7 @@ function rhoDotFlux(timestep,ph,en)
rho = getRho(ph,en) rho = getRho(ph,en)
rho_0 = getRho0(ph,en) rho_0 = getRho0(ph,en)
rho_0_sgl = rho_0(:,sgl) rho_0_sgl_mob = rho_0(:,mob)
v = reshape(stt%v(:,en),[prm%sum_N_sl,4]) !ToDo: MD: I think we should use state0 here v = reshape(stt%v(:,en),[prm%sum_N_sl,4]) !ToDo: MD: I think we should use state0 here
dot_gamma = rho(:,mob) * v * spread(prm%b_sl,2,4) dot_gamma = rho(:,mob) * v * spread(prm%b_sl,2,4)
@ -1229,12 +1228,12 @@ function rhoDotFlux(timestep,ph,en)
forall (s = 1:ns, t = 1:4) forall (s = 1:ns, t = 1:4)
v_0_nbr(s,t) = plasticState(ph_nbr)%state0(iV (s,t,ph_nbr),en_nbr) v_0_nbr(s,t) = plasticState(ph_nbr)%state0(iV (s,t,ph_nbr),en_nbr)
rho_0_sgl_nbr(s,t) = max(plasticState(ph_nbr)%state0(iRhoU(s,t,ph_nbr),en_nbr),0.0_pREAL) rho_0_sgl_mob_nbr(s,t) = max(plasticState(ph_nbr)%state0(iRhoU(s,t,ph_nbr),en_nbr),0.0_pREAL)
endforall endforall
where (rho_0_sgl_nbr * IPvolume0(ip_nbr,el_nbr) ** 0.667_pREAL < prm%rho_min & where (rho_0_sgl_mob_nbr * IPvolume0(ip_nbr,el_nbr) ** 0.667_pREAL < prm%rho_min &
.or. rho_0_sgl_nbr < prm%rho_significant) & .or. rho_0_sgl_mob_nbr < prm%rho_significant) &
rho_0_sgl_nbr = 0.0_pREAL rho_0_sgl_mob_nbr = 0.0_pREAL
normal_neighbor2me_defConf = math_det33(Favg) * matmul(math_inv33(transpose(Favg)), & normal_neighbor2me_defConf = math_det33(Favg) * matmul(math_inv33(transpose(Favg)), &
IPareaNormal0(1:3,n_nbr,ip_nbr,el_nbr)) ! normal of the interface in (average) deformed configuration (pointing neighbor => en) IPareaNormal0(1:3,n_nbr,ip_nbr,el_nbr)) ! normal of the interface in (average) deformed configuration (pointing neighbor => en)
normal_neighbor2me = matmul(transpose(F_e_nbr), normal_neighbor2me_defConf) & normal_neighbor2me = matmul(transpose(F_e_nbr), normal_neighbor2me_defConf) &
@ -1247,7 +1246,7 @@ function rhoDotFlux(timestep,ph,en)
topp = t + mod(t,2) - mod(t+1,2) topp = t + mod(t,2) - mod(t+1,2)
if (v_0_nbr(s,t) * math_inner(m(1:3,s,t), normal_neighbor2me) > 0.0_pREAL & ! flux from my neighbor to en == entering flux for en if (v_0_nbr(s,t) * math_inner(m(1:3,s,t), normal_neighbor2me) > 0.0_pREAL & ! flux from my neighbor to en == entering flux for en
.and. v_0(s,t) * v_0_nbr(s,t) >= 0.0_pREAL ) then ! ... only if no sign change in flux density .and. v_0(s,t) * v_0_nbr(s,t) >= 0.0_pREAL ) then ! ... only if no sign change in flux density
lineLength = rho_0_sgl_nbr(s,t) * v_0_nbr(s,t) & lineLength = rho_0_sgl_mob_nbr(s,t) * v_0_nbr(s,t) &
* math_inner(m(1:3,s,t), normal_neighbor2me) * a ! positive line length that wants to enter through this interface * math_inner(m(1:3,s,t), normal_neighbor2me) * a ! positive line length that wants to enter through this interface
where (dependentState(ph)%compatibility(c,:,s,n,en) > 0.0_pREAL) & where (dependentState(ph)%compatibility(c,:,s,n,en) > 0.0_pREAL) &
rhoDotFlux(:,t) = rhoDotFlux(1:ns,t) & rhoDotFlux(:,t) = rhoDotFlux(1:ns,t) &
@ -1287,7 +1286,7 @@ function rhoDotFlux(timestep,ph,en)
else ! sign change in flux density means sign change in stress which does not allow for dislocations to arive at the neighbor else ! sign change in flux density means sign change in stress which does not allow for dislocations to arive at the neighbor
transmissivity = 0.0_pREAL transmissivity = 0.0_pREAL
end if end if
lineLength = rho_0_sgl(s,t) * v_0(s,t) & lineLength = rho_0_sgl_mob(s,t) * v_0(s,t) &
* math_inner(m(1:3,s,t), normal_me2neighbor) * a ! positive line length of mobiles that wants to leave through this interface * math_inner(m(1:3,s,t), normal_me2neighbor) * a ! positive line length of mobiles that wants to leave through this interface
rhoDotFlux(s,t) = rhoDotFlux(s,t) - lineLength / geom(ph)%v_0(en) ! subtract dislocation flux from current type rhoDotFlux(s,t) = rhoDotFlux(s,t) - lineLength / geom(ph)%v_0(en) ! subtract dislocation flux from current type
rhoDotFlux(s,t+4) = rhoDotFlux(s,t+4) & rhoDotFlux(s,t+4) = rhoDotFlux(s,t+4) &