From 038adeeb4537a59954917f995af6e7059d3b1a08 Mon Sep 17 00:00:00 2001
From: ludovic noels <l.noels@ulg.ac.be>
Date: Sun, 4 Jun 2023 17:57:55 +0200
Subject: [PATCH] missing basic energy
---
NonLinearSolver/materialLaw/mlawHyperelastic.cpp | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/NonLinearSolver/materialLaw/mlawHyperelastic.cpp b/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
index 94bb565dc..50a6486a3 100644
--- a/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
+++ b/NonLinearSolver/materialLaw/mlawHyperelastic.cpp
@@ -143,9 +143,9 @@ double mlawHyperViscoElastic::deformationEnergy(const IPHyperViscoElastic& q) co
double Av, dAv,B_d;
static STensor3 dB_d;
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
-
-
+ Psy = getUpdatedBulkModulus(&q)*(1.+Av)*(0.5*trEe*trEe)+getUpdatedShearModulus(&q)*(1.+B_d)*STensorOperation::doubledot(devEe,devEe); //this is not correct we should do an integral
+ //
+ //
for (int i=0; i<_N; i++){
static STensor3 devEebranch;
Psy += _Ki[i]*0.5*(q._B[i])*(q._B[i])+_Gi[i]*STensorOperation::doubledot(q._A[i],q._A[i]);
--
GitLab