From cd7b885c959488f18b8c69b0d4aeb8e74f9c25b0 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Thu, 28 Mar 2024 14:21:10 +0100
Subject: [PATCH] [PL/THM] Use current solid state for C_el eval
---
ProcessLib/ThermoHydroMechanics/IntegrationPointData.h | 10 +++-------
1 file changed, 3 insertions(+), 7 deletions(-)
diff --git a/ProcessLib/ThermoHydroMechanics/IntegrationPointData.h b/ProcessLib/ThermoHydroMechanics/IntegrationPointData.h
index 5937aca687e..b85e8da06eb 100644
--- a/ProcessLib/ThermoHydroMechanics/IntegrationPointData.h
+++ b/ProcessLib/ThermoHydroMechanics/IntegrationPointData.h
@@ -98,10 +98,6 @@ struct IntegrationPointData final
MPL::VariableArray variable_array;
MPL::VariableArray variable_array_prev;
- auto const null_state = solid_material.createMaterialStateVariables();
- solid_material.initializeInternalStateVariables(t, x_position,
- *null_state);
-
using KV = MathLib::KelvinVector::KelvinVectorType<DisplacementDim>;
variable_array.stress.emplace<KV>(KV::Zero());
@@ -112,9 +108,9 @@ struct IntegrationPointData final
variable_array_prev.mechanical_strain.emplace<KV>(KV::Zero());
variable_array_prev.temperature = temperature;
- auto&& solution =
- solid_material.integrateStress(variable_array_prev, variable_array,
- t, x_position, dt, *null_state);
+ auto&& solution = solid_material.integrateStress(
+ variable_array_prev, variable_array, t, x_position, dt,
+ *material_state_variables);
if (!solution)
{
--
GitLab