From a867fe8f5f6d2e5ce7fa62bad0b181d3a2d2fe47 Mon Sep 17 00:00:00 2001 From: jaspervolmer Date: Tue, 19 Dec 2017 11:42:03 +0100 Subject: [PATCH] Small change in NR derivative --- small-strain/laminate/visco-plasticity.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/small-strain/laminate/visco-plasticity.py b/small-strain/laminate/visco-plasticity.py index f576db6..1e9192d 100644 --- a/small-strain/laminate/visco-plasticity.py +++ b/small-strain/laminate/visco-plasticity.py @@ -130,7 +130,7 @@ def viscoplastic(eps,eps_t,epse_t,ep_t,dt): res = -gamma0*dt*(sigeq_s/sig0)**(1./n) while np.linalg.norm(np.abs(res))/mu>1.e-6: - dres = 1.+3.*mu*gamma0*dt/n*((sigeq_s-3.*mu*dgamma)/sig0)**(1./n-1.) + dres = 1.+3.*mu*gamma0*dt/(n*sig0)*((sigeq_s-3.*mu*dgamma)/sig0)**(1./n-1.) dgamma -= res/dres dgamma[Z] = 0. res = dgamma-gamma0*dt*((sigeq_s-3.*mu*dgamma)/sig0)**(1./n);