diff --git a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
index c3b73149f8fdbafb87a618be78ccf700cb9ae26d..daa526c5dcce444798c33b06d0d5531fd2c1be4d 100644
--- a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
+++ b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
@@ -88,6 +88,8 @@ struct IntegrationPointData final
         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>;
 
diff --git a/ProcessLib/RichardsMechanics/IntegrationPointData.h b/ProcessLib/RichardsMechanics/IntegrationPointData.h
index 86c64530ab6ce0a1857896b849312300cb066b98..72d1fa8f6251b67752e027a29ff1d9564460612b 100644
--- a/ProcessLib/RichardsMechanics/IntegrationPointData.h
+++ b/ProcessLib/RichardsMechanics/IntegrationPointData.h
@@ -105,6 +105,8 @@ struct IntegrationPointData final
         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>;
 
diff --git a/ProcessLib/TH2M/IntegrationPointData.h b/ProcessLib/TH2M/IntegrationPointData.h
index 06cdbbed3990b628679805308fec1291970b8213..5733815007a60368a4689b1d144f92a01e8d2ecb 100644
--- a/ProcessLib/TH2M/IntegrationPointData.h
+++ b/ProcessLib/TH2M/IntegrationPointData.h
@@ -215,6 +215,8 @@ struct IntegrationPointData final
         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>;
 
diff --git a/ProcessLib/ThermoHydroMechanics/IntegrationPointData.h b/ProcessLib/ThermoHydroMechanics/IntegrationPointData.h
index c93c697e42fa9b441970bdcc2f0c6315df098556..89306893bdb877084f5623f6d13339f12eba450a 100644
--- a/ProcessLib/ThermoHydroMechanics/IntegrationPointData.h
+++ b/ProcessLib/ThermoHydroMechanics/IntegrationPointData.h
@@ -84,6 +84,8 @@ struct IntegrationPointData final
         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>;
 
diff --git a/ProcessLib/ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/ElasticTangentStiffnessModel.cpp b/ProcessLib/ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/ElasticTangentStiffnessModel.cpp
index 8e0dbdc015fc0dfaf2b6a2df552110f89f4999a4..846595e7e7925228fc4eb0860aae3cffa104119f 100644
--- a/ProcessLib/ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/ElasticTangentStiffnessModel.cpp
+++ b/ProcessLib/ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/ElasticTangentStiffnessModel.cpp
@@ -24,6 +24,7 @@ void ElasticTangentStiffnessModel<DisplacementDim>::eval(
     namespace MPL = MaterialPropertyLib;
 
     auto const null_state = solid_material_.createMaterialStateVariables();
+    solid_material_.initializeInternalStateVariables(x_t.t, x_t.x, *null_state);
 
     using KV = KelvinVector<DisplacementDim>;
 
diff --git a/ProcessLib/ThermoRichardsMechanics/ConstitutiveStress_StrainTemperature/ElasticTangentStiffnessModel.cpp b/ProcessLib/ThermoRichardsMechanics/ConstitutiveStress_StrainTemperature/ElasticTangentStiffnessModel.cpp
index c4613ee651965eacae9c9240c5f345f426b49067..9b29a24479e3fd247e4a047fc93294462e193206 100644
--- a/ProcessLib/ThermoRichardsMechanics/ConstitutiveStress_StrainTemperature/ElasticTangentStiffnessModel.cpp
+++ b/ProcessLib/ThermoRichardsMechanics/ConstitutiveStress_StrainTemperature/ElasticTangentStiffnessModel.cpp
@@ -27,6 +27,7 @@ void ElasticTangentStiffnessModel<DisplacementDim>::eval(
     MPL::VariableArray variable_array_prev;
 
     auto const null_state = solid_material_.createMaterialStateVariables();
+    solid_material_.initializeInternalStateVariables(x_t.t, x_t.x, *null_state);
 
     using KV = KelvinVector<DisplacementDim>;
 
diff --git a/Tests/ProcessLib/RichardsMechanics/MicroporosityComputation.cpp b/Tests/ProcessLib/RichardsMechanics/MicroporosityComputation.cpp
index a43a2e27bc7f4f5c58cae82a036e75298162cb48..5a383de0ac63eca8ee4dc2d2760bb644438b6631 100644
--- a/Tests/ProcessLib/RichardsMechanics/MicroporosityComputation.cpp
+++ b/Tests/ProcessLib/RichardsMechanics/MicroporosityComputation.cpp
@@ -37,6 +37,7 @@ computeElasticTangentStiffness(
     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>;