From 651e9eecce4ff945720c1ee835a58af1232ed946 Mon Sep 17 00:00:00 2001 From: ling <l.wu@ulg.ac.be> Date: Tue, 30 May 2023 10:11:20 +0200 Subject: [PATCH] bug --- NonLinearSolver/materialLaw/mlawHyperelastic.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/NonLinearSolver/materialLaw/mlawHyperelastic.cpp b/NonLinearSolver/materialLaw/mlawHyperelastic.cpp index 5ea592ab2..94bb565dc 100644 --- a/NonLinearSolver/materialLaw/mlawHyperelastic.cpp +++ b/NonLinearSolver/materialLaw/mlawHyperelastic.cpp @@ -142,8 +142,8 @@ double mlawHyperViscoElastic::deformationEnergy(const IPHyperViscoElastic& q) co double Av, dAv,B_d; static STensor3 dB_d; - evaluatePhiPCorrection(trEe, devEpr, A_v, dAprdEepr, B_d, dB_vddevEe); - Psy = getUpdatedBulkModulus(&q)*A_v*(0.5*trEe*trEe)+getUpdatedShearModulus(&q)*B_d*STensorOperation::doubledot(devEe,devEe); //this is not correct we should do an integral + evaluatePhiPCorrection(trEe, devEe, Av, dAv, B_d, dB_d); + Psy = getUpdatedBulkModulus(&q)*Av*(0.5*trEe*trEe)+getUpdatedShearModulus(&q)*B_d*STensorOperation::doubledot(devEe,devEe); //this is not correct we should do an integral for (int i=0; i<_N; i++){ -- GitLab