From d8fecff73f801c34d507185ab8515d6df579f08c Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Mon, 17 Feb 2020 12:38:29 +0100
Subject: [PATCH] [PL/RM] Add porosity to integration points.
And the corresponding push-back.
---
ProcessLib/RichardsMechanics/IntegrationPointData.h | 3 +++
1 file changed, 3 insertions(+)
diff --git a/ProcessLib/RichardsMechanics/IntegrationPointData.h b/ProcessLib/RichardsMechanics/IntegrationPointData.h
index d7558186d89..d8442413c26 100644
--- a/ProcessLib/RichardsMechanics/IntegrationPointData.h
+++ b/ProcessLib/RichardsMechanics/IntegrationPointData.h
@@ -56,6 +56,8 @@ struct IntegrationPointData final
typename ShapeMatricesTypePressure::GlobalDimNodalMatrixType dNdx_p;
double saturation;
+ double porosity = std::numeric_limits<double>::quiet_NaN();
+ double porosity_prev = std::numeric_limits<double>::quiet_NaN();
MaterialLib::Solids::MechanicsBase<DisplacementDim> const& solid_material;
std::unique_ptr<typename MaterialLib::Solids::MechanicsBase<
@@ -67,6 +69,7 @@ struct IntegrationPointData final
{
eps_prev = eps;
sigma_eff_prev = sigma_eff;
+ porosity_prev = porosity;
material_state_variables->pushBackState();
}
--
GitLab