From 39d45f00411a27092271e210c1a5d2365cfacdf3 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Fri, 3 Feb 2023 16:23:52 +0100
Subject: [PATCH] [PL/TH2M] Move C_el comp. before bulk modulus
This is required by MFront models for evaluation of the bulk modulus.
---
ProcessLib/TH2M/TH2MFEM-impl.h | 11 +++++------
1 file changed, 5 insertions(+), 6 deletions(-)
diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index ec0601f539b..2e6b36fb275 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -146,6 +146,8 @@ std::vector<ConstitutiveVariables<DisplacementDim>> TH2MLocalAssembler<
_element, Nu);
double const T = NT.dot(temperature);
+ double const T_dot = NT.dot(temperature_dot);
+ double const T_prev = T - T_dot * dt;
double const pGR = Np.dot(gas_pressure);
double const pCap = Np.dot(capillary_pressure);
double const pLR = pGR - pCap;
@@ -161,7 +163,9 @@ std::vector<ConstitutiveVariables<DisplacementDim>> TH2MLocalAssembler<
vars.liquid_phase_pressure = pLR;
// medium properties
- auto const K_S = ip_data.solid_material.getBulkModulus(t, pos);
+ auto const C_el =
+ ip_data.computeElasticTangentStiffness(t, pos, dt, T_prev, T);
+ auto const K_S = ip_data.solid_material.getBulkModulus(t, pos, &C_el);
ip_data.alpha_B = medium.property(MPL::PropertyType::biot_coefficient)
.template value<double>(vars, pos, t, dt);
@@ -273,9 +277,6 @@ std::vector<ConstitutiveVariables<DisplacementDim>> TH2MLocalAssembler<
// isotropic solid phase volumetric thermal expansion coefficient
ip_data.beta_T_SR = Invariants::trace(ip_data.alpha_T_SR);
- double const T_dot = NT.dot(temperature_dot);
- double const T_prev = T - T_dot * dt;
-
MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const
dthermal_strain = ip_data.alpha_T_SR * T_dot * dt;
@@ -287,8 +288,6 @@ std::vector<ConstitutiveVariables<DisplacementDim>> TH2MLocalAssembler<
if (solid_phase.hasProperty(MPL::PropertyType::swelling_stress_rate))
{
- auto const C_el =
- ip_data.computeElasticTangentStiffness(t, pos, dt, T_prev, T);
eps_m.noalias() +=
C_el.inverse() * (ip_data.sigma_sw - ip_data.sigma_sw_prev);
}
--
GitLab