From c9c286b10173d1b94cba1cc867f951af4b703a68 Mon Sep 17 00:00:00 2001
From: Van Dung NGUYEN <vdg.nguyen@gmail.com>
Date: Tue, 10 Jan 2023 22:35:40 +0100
Subject: [PATCH] fix bug when elastic stiffnes is not filled for DG
formulation
---
dG3D/src/dG3DMaterialLaw.cpp | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/dG3D/src/dG3DMaterialLaw.cpp b/dG3D/src/dG3DMaterialLaw.cpp
index fcb797284..c568c2a25 100644
--- a/dG3D/src/dG3DMaterialLaw.cpp
+++ b/dG3D/src/dG3DMaterialLaw.cpp
@@ -2965,7 +2965,10 @@ void HyperViscoElasticDG3DMaterialLaw::stress(IPVariable* ipv, const IPVariable*
HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw::HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw(const int num, const double rho, const double E,const double nu, const double tol ,
const bool matrixbyPerturbation, const double pert):
- dG3DMaterialLaw(num,rho), _viscoLaw(num,E,nu,rho,tol,matrixbyPerturbation,pert){};
+ dG3DMaterialLaw(num,rho), _viscoLaw(num,E,nu,rho,tol,matrixbyPerturbation,pert)
+{
+ fillElasticStiffness(E, nu, elasticStiffness);
+};
HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw::HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw(const HyperViscoElastoPlasticPowerYieldDG3DMaterialLaw& src): dG3DMaterialLaw(src),
_viscoLaw(src._viscoLaw){};
--
GitLab