From 16d6dfcab6993b26d99c17e65fb9eb9f8bd39cc8 Mon Sep 17 00:00:00 2001 From: abbatux <31911482+abbatux@users.noreply.github.com> Date: Fri, 8 Dec 2017 11:40:44 +1100 Subject: [PATCH] Bugfix for effective plastic strain rate calc --- src/USER-SMD/pair_smd_tlsph.cpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/USER-SMD/pair_smd_tlsph.cpp b/src/USER-SMD/pair_smd_tlsph.cpp index d6ddc9c72d..8dac846017 100644 --- a/src/USER-SMD/pair_smd_tlsph.cpp +++ b/src/USER-SMD/pair_smd_tlsph.cpp @@ -775,9 +775,9 @@ void PairTlsph::AssembleStress() { // compute a characteristic time over which to average the plastic strain double tav = 1000 * radius[i] / (Lookup[SIGNAL_VELOCITY][itype]); - eff_plastic_strain_rate[i] -= eff_plastic_strain_rate[i] / tav; - eff_plastic_strain_rate[i] += (plastic_strain_increment / dt) / tav; - eff_plastic_strain_rate[i] = MAX(0.0, eff_plastic_strain_rate[i]); + eff_plastic_strain_rate[i] -= eff_plastic_strain_rate[i] * dt / tav; + eff_plastic_strain_rate[i] += plastic_strain_increment / tav; + eff_plastic_strain_rate[i] = MAX(0.0, eff_plastic_strain_rate[i]); /* * assemble total stress from pressure and deviatoric stress