From c38e5b84a6c245a2980491cd61f79cc662cb4895 Mon Sep 17 00:00:00 2001
From: labourasse <labourassee@gmail.com>
Date: Tue, 5 Nov 2024 10:35:59 +0100
Subject: [PATCH] add platicity criterion
---
src/scheme/HyperplasticSolver.cpp | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/src/scheme/HyperplasticSolver.cpp b/src/scheme/HyperplasticSolver.cpp
index a0d47d8a2..e56c12c39 100644
--- a/src/scheme/HyperplasticSolver.cpp
+++ b/src/scheme/HyperplasticSolver.cpp
@@ -190,8 +190,12 @@ class HyperplasticSolverHandler::HyperplasticSolver final : public HyperplasticS
{
const R3x3 Id = identity;
const R3x3 Sbar = S - 1. / 3 * trace(S) * Id;
- const double normS = 0;
- return normS;
+ double chi = 0;
+ const double normS = frobeniusNorm(Sbar);
+ if ((normS - std::sqrt(2. / 3) * yield) > 0.) {
+ chi = 1;
+ }
+ return chi;
}
NodeValuePerCell<const Rdxd>
--
GitLab