From e8b5a6dfe40113c65015cf3c2b9b892149e3863a Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Tue, 10 Mar 2020 14:02:09 +0100 Subject: [PATCH] bugfix: full stress tensor should be used --- src/kinematics_slipplane_opening.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/kinematics_slipplane_opening.f90 b/src/kinematics_slipplane_opening.f90 index 3f77ecbbb..ab9035918 100644 --- a/src/kinematics_slipplane_opening.f90 +++ b/src/kinematics_slipplane_opening.f90 @@ -141,9 +141,9 @@ subroutine kinematics_slipplane_opening_LiAndItsTangent(Ld, dLd_dTstar, S, ipc, dLd_dTstar = 0.0_pReal do i = 1, prm%totalNslip - traction_d = math_mul33xx33(S,prm%P_d(3,3,i)) - traction_t = math_mul33xx33(S,prm%P_t(3,3,i)) - traction_n = math_mul33xx33(S,prm%P_n(3,3,i)) + traction_d = math_tensordot(S,prm%P_d(1:3,1:3,i)) + traction_t = math_tensordot(S,prm%P_t(1:3,1:3,i)) + traction_n = math_tensordot(S,prm%P_n(1:3,1:3,i)) traction_crit = prm%critLoad(i)* damage(homog)%p(damageOffset) ! degrading critical load carrying capacity by damage