From adbd9c100c68915e3b6be462dc33cef36b571af2 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Thu, 5 Sep 2019 14:00:13 +0200
Subject: [PATCH] [MatL] MCAS; Use incremental form.
MFront's computations are absolute, and incorporation
of initial stresses requires computation of initial
strain.
---
MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront | 5 +++++
1 file changed, 5 insertions(+)
diff --git a/MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront b/MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront
index f9f7bbf7e87..97d5cfa1b82 100644
--- a/MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront
+++ b/MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront
@@ -82,6 +82,11 @@ a.setEntryName("TensionCutOffParameter");
cos_3_lodeT = cos(3. * lodeT);
sin_3_lodeT = sin(3. * lodeT);
tan_3_lodeT = tan(3. * lodeT);
+
+ // Compute initial elastic strain
+ const auto S = invert(D);
+ eel = S*sig;
+
// elastic prediction
const auto sig_el = computeElasticPrediction();
const auto s_el = deviator(sig_el);
--
GitLab