jianmu
/
lammps
forked from lijiext/lammps
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Bugfix for effective plastic strain rate calc

This commit is contained in:
abbatux 2017-12-08 11:40:44 +11:00 committed by GitHub
parent d029cb9002
commit 16d6dfcab6
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/USER-SMD/pair_smd_tlsph.cpp
View File

@ -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