added check for CFL condition in constitutive_nonlocal_dotState

This commit is contained in:
Christoph Kords 2011-04-06 09:07:36 +00:00
parent e00d073ee3
commit 4ea1fe436b
1 changed files with 12 additions and 0 deletions

View File

@ -1611,6 +1611,17 @@ forall (s = 1:ns, t = 1:4, rhoSgl(s,t+4) * constitutive_nonlocal_v(s,t,g,ip,el)
!****************************************************************************
!*** check CFL condition for flux
if (any(minval(mesh_ipVolume(ip,el) / max(tiny(1.0_pReal),mesh_ipArea(:,ip,el))) &
< constitutive_nonlocal_v(1:ns,1:4,g,ip,el) * timestep)) then
dotState%p = NaN
return
endif
!**************************************************************************** !****************************************************************************
!*** calculate limits for stable dipole height !*** calculate limits for stable dipole height
@ -1627,6 +1638,7 @@ dUpper(1:ns,2) = min( 1.0_pReal / sqrt( sum(abs(rhoSgl),2)+sum(rhoDip,2) ), &
dUpper(1:ns,1) = dUpper(1:ns,2) / ( 1.0_pReal - constitutive_nonlocal_nu(myInstance) ) dUpper(1:ns,1) = dUpper(1:ns,2) / ( 1.0_pReal - constitutive_nonlocal_nu(myInstance) )
!**************************************************************************** !****************************************************************************
!*** dislocation remobilization (bauschinger effect) !*** dislocation remobilization (bauschinger effect)