From b8a3e5d958d147e1ced7dba67fd58580bdda975b Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Fri, 3 Mar 2023 17:21:29 +0100
Subject: [PATCH] [PL/SD] Calculate eps prev instead of storing it
---
ProcessLib/SmallDeformation/SmallDeformationFEM.h | 7 ++-----
1 file changed, 2 insertions(+), 5 deletions(-)
diff --git a/ProcessLib/SmallDeformation/SmallDeformationFEM.h b/ProcessLib/SmallDeformation/SmallDeformationFEM.h
index 79f5d2942c9..82380c40d77 100644
--- a/ProcessLib/SmallDeformation/SmallDeformationFEM.h
+++ b/ProcessLib/SmallDeformation/SmallDeformationFEM.h
@@ -54,7 +54,7 @@ struct IntegrationPointData final
}
typename BMatricesType::KelvinVectorType sigma, sigma_prev;
- typename BMatricesType::KelvinVectorType eps, eps_prev;
+ typename BMatricesType::KelvinVectorType eps;
double free_energy_density = 0;
MaterialLib::Solids::MechanicsBase<DisplacementDim> const& solid_material;
@@ -68,7 +68,6 @@ struct IntegrationPointData final
void pushBackState()
{
- eps_prev = eps;
sigma_prev = sigma;
material_state_variables->pushBackState();
}
@@ -164,7 +163,6 @@ public:
// Previous time step values are not initialized and are set later.
ip_data.sigma_prev.resize(kelvin_vector_size);
- ip_data.eps_prev.resize(kelvin_vector_size);
_secondary_data.N[ip] = shape_matrices[ip].N;
}
@@ -245,7 +243,6 @@ public:
MPL::VariableArray variables_prev;
MPL::VariableArray variables;
- auto const& eps_prev = ip_data.eps_prev;
auto const& sigma_prev = ip_data.sigma_prev;
auto& eps = ip_data.eps;
@@ -270,7 +267,7 @@ public:
sigma_prev);
variables_prev.mechanical_strain
.emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
- eps_prev);
+ B * (u - u_dot * dt));
double const T_ref =
_process_data.reference_temperature
--
GitLab