From 7fff4bb048b18c28dc997f65e511a4c8b35802cd Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Mon, 2 Sep 2019 12:12:56 +0200
Subject: [PATCH] [PL] RM; Setting initial stress in FEM.
Was forgotten in the main PR #2635.
---
.../RichardsMechanics/RichardsMechanicsFEM.h | 22 ++++++++++++++++++-
1 file changed, 21 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM.h
index 99e672d7f5f..c90a1b48fee 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM.h
@@ -95,7 +95,27 @@ public:
for (unsigned ip = 0; ip < n_integration_points; ip++)
{
- _ip_data[ip].pushBackState();
+ auto& ip_data = _ip_data[ip];
+
+ /// Set initial stress from parameter.
+ if (_process_data.initial_stress != nullptr)
+ {
+ ParameterLib::SpatialPosition const x_position{
+ boost::none, _element.getID(), ip,
+ MathLib::Point3d(
+ interpolateCoordinates<ShapeFunctionDisplacement,
+ ShapeMatricesTypeDisplacement>(
+ _element, ip_data.N_u))};
+
+ ip_data.sigma_eff =
+ MathLib::KelvinVector::symmetricTensorToKelvinVector<
+ DisplacementDim>((*_process_data.initial_stress)(
+ std::numeric_limits<
+ double>::quiet_NaN() /* time independent */,
+ x_position));
+ }
+
+ ip_data.pushBackState();
}
}
--
GitLab