From b5050b4f21e9ac9ca5de29af7aba986360c8fd5d Mon Sep 17 00:00:00 2001
From: cbsilver <christian.silbermann@ifgt.tu-freiberg.de>
Date: Tue, 13 Jun 2023 14:17:53 +0200
Subject: [PATCH] Bugfix algorithmic tangent, only relevant for theta<1
---
.../MFront/ModCamClay_semiExplParaInitNLnu_inc.mfront | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/MaterialLib/SolidModels/MFront/ModCamClay_semiExplParaInitNLnu_inc.mfront b/MaterialLib/SolidModels/MFront/ModCamClay_semiExplParaInitNLnu_inc.mfront
index 21943f1cd8b..828d2ca5197 100644
--- a/MaterialLib/SolidModels/MFront/ModCamClay_semiExplParaInitNLnu_inc.mfront
+++ b/MaterialLib/SolidModels/MFront/ModCamClay_semiExplParaInitNLnu_inc.mfront
@@ -200,7 +200,7 @@ v.setEntryName("VolumeRatio"); // Total volume per solid volume = 1 + pore numb
dfrpc_ddlp = dfrpc_ddeplV * ntr;
dfrpc_ddeel = dfrpc_ddeplV * dlp * (id2 | dn_ddeel);
- dfrpc_ddrpc = 1 + (deplV * the + dfrpc_ddeplV * dlp * trace(dn_ddrpc)) * theta;
+ dfrpc_ddrpc = 1 + deplV * the * theta + dfrpc_ddeplV * dlp * trace(dn_ddrpc);
}
// compute resulting stresses
--
GitLab