From b33d7e05857bfd7fd8d4a195761436f06e6c63df Mon Sep 17 00:00:00 2001 From: Fengbo Han Date: Thu, 27 Jul 2017 16:28:33 +0200 Subject: [PATCH] implemented yield stop criteria --- src/DAMASK_spectral.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/DAMASK_spectral.f90 b/src/DAMASK_spectral.f90 index 29d9b672a..877970e74 100755 --- a/src/DAMASK_spectral.f90 +++ b/src/DAMASK_spectral.f90 @@ -699,21 +699,21 @@ program DAMASK_spectral call utilities_calcPlasticity(yieldStressNew, plasticStrainNew, eqStressNew, eqTotalStrainNew, & eqPlasticStrainNew, plasticWorkNew) - if(stopFlag .eq. 'totalStrain') then + if(stopFlag == 'totalStrain') then if(eqTotalStrainNew > yieldStopValue) then yieldStress = yieldStressOld * (eqTotalStrainNew - yieldStopValue)/(eqTotalStrainNew - eqTotalStrainOld) & ! linear interpolation of stress values + yieldStressNew * (yieldStopValue - eqTotalStrainOld)/(eqTotalStrainNew - eqTotalStrainOld) plasticStrainRate = (plasticStrainNew - plasticStrainOld)/(time - time0) ! calculate plastic strain rate yieldStopSatisfied = .True. endif - elseif(stopFlag .eq. 'plasticStrain') then + elseif(stopFlag == 'plasticStrain') then if(eqPlasticStrainNew > yieldStopValue) then yieldStress = yieldStressOld * (eqPlasticStrainNew - yieldStopValue)/(eqPlasticStrainNew - eqPlasticStrainOld) & + yieldStressNew * (yieldStopValue - eqPlasticStrainOld)/(eqPlasticStrainNew - eqPlasticStrainOld) plasticStrainRate = (plasticStrainNew - plasticStrainOld)/(time - time0) yieldStopSatisfied = .True. endif - elseif(stopFlag .eq. 'plasticWork') then + elseif(stopFlag == 'plasticWork') then if(plasticWorkNew > yieldStopValue) then yieldStress = yieldStressOld * (plasticWorkNew - yieldStopValue)/(plasticWorkNew - plasticWorkOld) & + yieldStressNew * (yieldStopValue - plasticWorkOld)/(plasticWorkNew - plasticWorkOld)