From 0185eda888879f9875a5366ba30f3d11fc9d8143 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sat, 7 Jan 2023 13:03:12 +0100 Subject: [PATCH] correct calculation of strain was introduced when changing to Voigt notation in fa8218124ace56b392ec1b639379106e7738fce7 --- src/phase_damage_isobrittle.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/phase_damage_isobrittle.f90 b/src/phase_damage_isobrittle.f90 index ff91265b9..0e00974e3 100644 --- a/src/phase_damage_isobrittle.f90 +++ b/src/phase_damage_isobrittle.f90 @@ -93,7 +93,7 @@ end function isobrittle_init !-------------------------------------------------------------------------------------------------- -!> @brief calculates derived quantities from state +!> @brief !-------------------------------------------------------------------------------------------------- module subroutine isobrittle_deltaState(C, Fe, ph,en) @@ -109,11 +109,11 @@ module subroutine isobrittle_deltaState(C, Fe, ph,en) r_W - epsilon = math_33toVoigt6_strain(matmul(transpose(Fe),Fe)-math_I3) + epsilon = math_33toVoigt6_strain(0.5_pReal*(matmul(transpose(Fe),Fe)-math_I3)) associate(prm => param(ph), stt => state(ph), dlt => deltaState(ph)) - r_W = (0.5_pReal*dot_product(epsilon,matmul(C,epsilon)))/prm%W_crit + r_W = (2.0_pReal*dot_product(epsilon,matmul(C,epsilon)))/prm%W_crit dlt%r_W(en) = merge(r_W - stt%r_W(en), 0.0_pReal, r_W > stt%r_W(en)) end associate @@ -122,7 +122,7 @@ end subroutine isobrittle_deltaState !-------------------------------------------------------------------------------------------------- -!> @brief writes results to HDF5 output file +!> @brief Write results to HDF5 output file. !-------------------------------------------------------------------------------------------------- module subroutine isobrittle_results(phase,group)