From 3017bcdf6e7029b3816706022883129efebe3562 Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Tue, 4 Sep 2018 11:31:47 +0200
Subject: [PATCH] [BGRa] A bug fixing: return try stress if the norm of its
deviator is zero
---
MaterialLib/SolidModels/CreepBGRa.cpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/MaterialLib/SolidModels/CreepBGRa.cpp b/MaterialLib/SolidModels/CreepBGRa.cpp
index fbfe7b29ee8..78da877423e 100644
--- a/MaterialLib/SolidModels/CreepBGRa.cpp
+++ b/MaterialLib/SolidModels/CreepBGRa.cpp
@@ -57,7 +57,7 @@ CreepBGRa<DisplacementDim>::integrateStress(
// In case |s_{try}| is zero and _n < 3 (rare case).
if (norm_s_try < std::numeric_limits<double>::epsilon() * C(0, 0))
{
- return {std::make_tuple(sigma_prev, createMaterialStateVariables(), C)};
+ return {std::make_tuple(sigma_try, createMaterialStateVariables(), C)};
}
ResidualVectorType solution = sigma_try;
--
GitLab