From 3d4b2d2f527398816dd673b509f979ca38cb70fa Mon Sep 17 00:00:00 2001 From: Luv Sharma Date: Fri, 21 Nov 2014 13:12:54 +0000 Subject: [PATCH] tolerance for damage can be any positive value. Set to 0.001 if found negative. --- code/damage_anisoBrittle.f90 | 2 +- code/damage_anisoDuctile.f90 | 2 +- code/damage_isoBrittle.f90 | 2 +- code/damage_isoDuctile.f90 | 2 +- 4 files changed, 4 insertions(+), 4 deletions(-) diff --git a/code/damage_anisoBrittle.f90 b/code/damage_anisoBrittle.f90 index 5fd350a5e..d1efa231f 100644 --- a/code/damage_anisoBrittle.f90 +++ b/code/damage_anisoBrittle.f90 @@ -224,7 +224,7 @@ subroutine damage_anisoBrittle_init(fileUnit) min(lattice_NcleavageSystem(1:lattice_maxNcleavageFamily,phase),& ! limit active cleavage systems per family to min of available and requested damage_anisoBrittle_Ncleavage(1:lattice_maxNcleavageFamily,instance)) damage_anisoBrittle_totalNcleavage(instance) = sum(damage_anisoBrittle_Ncleavage(:,instance)) ! how many cleavage systems altogether - if (damage_anisoBrittle_aTol_damage(instance) >= 1.0e-3_pReal) & + if (damage_anisoBrittle_aTol_damage(instance) < 0.0_pReal) & damage_anisoBrittle_aTol_damage(instance) = 1.0e-3_pReal ! default absolute tolerance 1e-3 if (damage_anisoBrittle_aTol_disp(instance) >= 1.0e-3_pReal) & damage_anisoBrittle_aTol_disp(instance) = 1.0e-3_pReal ! default absolute tolerance 1e-3 diff --git a/code/damage_anisoDuctile.f90 b/code/damage_anisoDuctile.f90 index 5ba3f87a9..f518afa16 100644 --- a/code/damage_anisoDuctile.f90 +++ b/code/damage_anisoDuctile.f90 @@ -203,7 +203,7 @@ subroutine damage_anisoDuctile_init(fileUnit) min(lattice_NslipSystem(1:lattice_maxNslipFamily,phase),& ! limit active cleavage systems per family to min of available and requested damage_anisoDuctile_Nslip(1:lattice_maxNslipFamily,instance)) damage_anisoDuctile_totalNslip(instance) = sum(damage_anisoDuctile_Nslip(:,instance)) - if (damage_anisoDuctile_aTol_damage(instance) >= 1.0e-3_pReal) & + if (damage_anisoDuctile_aTol_damage(instance) < 0.0_pReal) & damage_anisoDuctile_aTol_damage(instance) = 1.0e-3_pReal ! default absolute tolerance 1e-3 if (any(damage_anisoDuctile_critAccShear(:,instance) < 0.0_pReal)) & call IO_error(211_pInt,el=instance,ext_msg='critical_accshear ('//LOCAL_DAMAGE_anisoDuctile_LABEL//')') diff --git a/code/damage_isoBrittle.f90 b/code/damage_isoBrittle.f90 index e18ed8ca0..fa62f4503 100644 --- a/code/damage_isoBrittle.f90 +++ b/code/damage_isoBrittle.f90 @@ -179,7 +179,7 @@ subroutine damage_isoBrittle_init(fileUnit) NofMyPhase=count(material_phase==phase) instance = phase_damageInstance(phase) ! sanity checks - if (damage_isoBrittle_aTol(instance) >= 1.0e-3_pReal) & + if (damage_isoBrittle_aTol(instance) < 0.0_pReal) & damage_isoBrittle_aTol(instance) = 1.0e-3_pReal ! default absolute tolerance 1e-3 if (damage_isoBrittle_critStrainEnergy(instance) <= 0.0_pReal) & call IO_error(211_pInt,el=instance,ext_msg='critical_strain_energy ('//LOCAL_DAMAGE_isoBrittle_LABEL//')') diff --git a/code/damage_isoDuctile.f90 b/code/damage_isoDuctile.f90 index e5ce8155f..00130d0ba 100644 --- a/code/damage_isoDuctile.f90 +++ b/code/damage_isoDuctile.f90 @@ -177,7 +177,7 @@ subroutine damage_isoDuctile_init(fileUnit) NofMyPhase=count(material_phase==phase) instance = phase_damageInstance(phase) ! sanity checks - if (damage_isoDuctile_aTol(instance) >= 1.0e-3_pReal) & + if (damage_isoDuctile_aTol(instance) < 0.0_pReal) & damage_isoDuctile_aTol(instance) = 1.0e-3_pReal ! default absolute tolerance 1e-3 if (damage_isoDuctile_critpStrain(instance) <= 0.0_pReal) & call IO_error(211_pInt,el=instance,ext_msg='critical_plastic_strain ('//LOCAL_DAMAGE_isoDuctile_LABEL//')')