From c10306567f569f6c36a302f72bf8d8e18fd895e3 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Sun, 21 Feb 2021 21:35:30 +0100
Subject: [PATCH] [PL/RM] Fix error in mechanical strain computation
Only happened in the assemble function, which does not
have a ctest with swelling.
---
ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index 4b587a5e1a3..3e1cc7ac731 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -537,7 +537,7 @@ void RichardsMechanicsLocalAssembler<
//
eps_m.noalias() =
solid_phase.hasProperty(MPL::PropertyType::swelling_stress_rate)
- ? eps - C_el.inverse() * sigma_sw
+ ? eps + C_el.inverse() * sigma_sw
: eps;
variables[static_cast<int>(MPL::Variable::mechanical_strain)]
.emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
--
GitLab