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

lower bound for damage is residual stiffness

This commit is contained in:
Shaokang Zhang 2015-06-22 14:27:32 +00:00
parent 12bc4f01ef
commit a40d247f7f
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
code/spectral_damage.f90
View File

@ -257,6 +257,8 @@ end function spectral_damage_solution
!> @brief forms the spectral damage residual vector !> @brief forms the spectral damage residual vector
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine spectral_damage_formResidual(in,x_scal,f_scal,dummy,ierr) subroutine spectral_damage_formResidual(in,x_scal,f_scal,dummy,ierr)
use numerics, only: &
residualStiffness
use mesh, only: & use mesh, only: &
gridLocal gridLocal
use math, only: & use math, only: &
@ -325,7 +327,7 @@ subroutine spectral_damage_formResidual(in,x_scal,f_scal,dummy,ierr)
call utilities_fourierGreenConvolution(D_ref, mobility_ref, params%timeinc) call utilities_fourierGreenConvolution(D_ref, mobility_ref, params%timeinc)
call utilities_FFTscalarBackward() call utilities_FFTscalarBackward()
where(scalarField_realMPI > 1.0_pReal) scalarField_realMPI = 1.0_pReal where(scalarField_realMPI > 1.0_pReal) scalarField_realMPI = 1.0_pReal
where(scalarField_realMPI < 0.0_pReal) scalarField_realMPI = 0.0_pReal where(scalarField_realMPI < residualStiffness) scalarField_realMPI = residualStiffness
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! constructing residual ! constructing residual