OGS
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits > Class Template Reference

Detailed Description

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
class ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >

Definition at line 31 of file ThermoRichardsMechanicsFEM.h.

#include <ThermoRichardsMechanicsFEM.h>

Inheritance diagram for ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >:
[legend]
Collaboration diagram for ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >:
[legend]

Classes

class  LocalMatrices

Public Types

using ShapeMatricesTypeDisplacement
using ShapeMatricesType
using GlobalDimMatrixType = typename ShapeMatricesType::GlobalDimMatrixType
using GlobalDimVectorType = typename ShapeMatricesType::GlobalDimVectorType
using BMatricesType
using KelvinVectorType = typename BMatricesType::KelvinVectorType
using IpData
using Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>
using SymmetricTensor = Eigen::Matrix<double, KelvinVectorSize, 1>

Public Member Functions

 ThermoRichardsMechanicsLocalAssembler (ThermoRichardsMechanicsLocalAssembler const &)=delete
 ThermoRichardsMechanicsLocalAssembler (ThermoRichardsMechanicsLocalAssembler &&)=delete
 ThermoRichardsMechanicsLocalAssembler (MeshLib::Element const &e, std::size_t const, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > &process_data)
void setInitialConditionsConcrete (Eigen::VectorXd const local_x, double const t, int const process_id) override
void assembleWithJacobian (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, std::vector< double > &local_rhs_data, std::vector< double > &local_Jac_data) override
int getNumberOfVectorElementsForDeformation () const override
void initializeConcrete () override
void computeSecondaryVariableConcrete (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev) override
Eigen::Map< const Eigen::RowVectorXd > getShapeMatrix (const unsigned integration_point) const override
 Provides the shape matrix at the given integration point.
Public Member Functions inherited from ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >
 LocalAssemblerInterface (MeshLib::Element const &e, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > &process_data)
std::size_t setIPDataInitialConditions (std::string_view name, double const *values, int const integration_order)
std::vector< double > getMaterialStateVariableInternalState (std::function< std::span< double >(typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables &)> const &get_values_span, int const &n_components) const
unsigned getNumberOfIntegrationPoints () const
int getMaterialID () const
MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables const & getMaterialStateVariablesAt (unsigned integration_point) const
void postTimestepConcrete (Eigen::VectorXd const &, Eigen::VectorXd const &, double const, double const, int const) override
Public Member Functions inherited from ProcessLib::LocalAssemblerInterface
virtual ~LocalAssemblerInterface ()=default
virtual void setInitialConditions (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, int const process_id)
virtual void initialize (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)
virtual void preAssemble (double const, double const, std::vector< double > const &)
virtual void assemble (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
virtual void assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
virtual void assembleWithJacobianForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)
virtual void computeSecondaryVariable (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id)
virtual void preTimestep (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, double const delta_t)
virtual void postTimestep (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
void postNonLinearSolver (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
virtual Eigen::Vector3d getFlux (MathLib::Point3d const &, double const, std::vector< double > const &) const
virtual Eigen::Vector3d getFlux (MathLib::Point3d const &, double const, std::vector< std::vector< double > > const &) const
 Fits to staggered scheme.
Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual ~ExtrapolatableElement ()=default

Static Public Attributes

static int const KelvinVectorSize
static constexpr auto & N_u_op

Private Member Functions

void assembleWithJacobianSingleIP (double const t, double const dt, ParameterLib::SpatialPosition const &x_position, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, IpData const &ip_data, typename ConstitutiveTraits::ConstitutiveSetting &CS, MaterialPropertyLib::Medium &medium, LocalMatrices &out, typename ConstitutiveTraits::StatefulData &current_state, typename ConstitutiveTraits::StatefulDataPrev const &prev_state, MaterialStateData< DisplacementDim > &mat_state, typename ConstitutiveTraits::OutputData &output_data) const
void addToLocalMatrixData (double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, LocalMatrices const &loc_mat, std::vector< double > &local_rhs_data, std::vector< double > &local_Jac_data) const
void massLumping (LocalMatrices &loc_mat) const
void convertInitialStressType (unsigned const ip, double const t, ParameterLib::SpatialPosition const x_position, MaterialPropertyLib::Medium const &medium, MPL::VariableArray const &variables, double const p_at_ip)

Static Private Member Functions

static constexpr auto localDOF (std::vector< double > const &x)
static auto block_uu (auto &mat)
static auto block_up (auto &mat)
static auto block_uT (auto &mat)
static auto block_pu (auto &mat)
static auto block_pp (auto &mat)
static auto block_pT (auto &mat)
static auto block_Tp (auto &mat)
static auto block_TT (auto &mat)
static auto block_u (auto &vec)
static auto block_p (auto &vec)
static auto block_T (auto &vec)

Private Attributes

std::vector< IpData > ip_data_

Static Private Attributes

static constexpr int temperature_index = 0
static constexpr int temperature_size = ShapeFunction::NPOINTS
static constexpr int pressure_index = temperature_size
static constexpr int pressure_size = ShapeFunction::NPOINTS
static constexpr int displacement_index = 2 * ShapeFunction::NPOINTS
static constexpr int displacement_size

Additional Inherited Members

Static Public Member Functions inherited from ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >
static auto getReflectionDataForOutput ()
Protected Attributes inherited from ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >
ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > & process_data_
std::vector< typename ConstitutiveTraits::StatefulData > current_states_
std::vector< typename ConstitutiveTraits::StatefulDataPrev > prev_states_
std::vector< MaterialStateData< DisplacementDim > > material_states_
NumLib::GenericIntegrationMethod const & integration_method_
MeshLib::Element const & element_
bool const is_axially_symmetric_
ConstitutiveTraits::SolidConstitutiveRelation const & solid_material_
std::vector< typename ConstitutiveTraits::OutputData > output_data_

Member Typedef Documentation

◆ BMatricesType

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::BMatricesType
Initial value:
BMatrixPolicyType<ShapeFunctionDisplacement, DisplacementDim>

Definition at line 53 of file ThermoRichardsMechanicsFEM.h.

◆ GlobalDimMatrixType

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::GlobalDimMatrixType = typename ShapeMatricesType::GlobalDimMatrixType

Definition at line 50 of file ThermoRichardsMechanicsFEM.h.

◆ GlobalDimVectorType

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::GlobalDimVectorType = typename ShapeMatricesType::GlobalDimVectorType

Definition at line 51 of file ThermoRichardsMechanicsFEM.h.

◆ Invariants

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>

Definition at line 63 of file ThermoRichardsMechanicsFEM.h.

◆ IpData

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::IpData
Initial value:
IntegrationPointData<ShapeMatricesTypeDisplacement,
ShapeMatricesType, DisplacementDim,
ShapeFunctionDisplacement::NPOINTS>
ProcessLib::HydroMechanics::HydroMechanicsLocalAssembler::ShapeMatricesTypeDisplacement
ShapeMatrixPolicyType< ShapeFunctionDisplacement, DisplacementDim > ShapeMatricesTypeDisplacement
Definition HydroMechanicsFEM.h:163
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler::ShapeMatricesType
ShapeMatrixPolicyType< ShapeFunction, DisplacementDim > ShapeMatricesType
Definition ThermoRichardsMechanicsFEM.h:47

Definition at line 57 of file ThermoRichardsMechanicsFEM.h.

◆ KelvinVectorType

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::KelvinVectorType = typename BMatricesType::KelvinVectorType

Definition at line 55 of file ThermoRichardsMechanicsFEM.h.

◆ ShapeMatricesType

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::ShapeMatricesType
Initial value:
ShapeMatrixPolicyType<ShapeFunction, DisplacementDim>
ShapeMatrixPolicyType
EigenFixedShapeMatrixPolicy< ShapeFunction, GlobalDim > ShapeMatrixPolicyType
Definition ShapeMatrixPolicy.h:134

Definition at line 47 of file ThermoRichardsMechanicsFEM.h.

◆ ShapeMatricesTypeDisplacement

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::ShapeMatricesTypeDisplacement
Initial value:
ShapeMatrixPolicyType<ShapeFunctionDisplacement, DisplacementDim>

Definition at line 43 of file ThermoRichardsMechanicsFEM.h.

◆ SymmetricTensor

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
using ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::SymmetricTensor = Eigen::Matrix<double, KelvinVectorSize, 1>

Definition at line 65 of file ThermoRichardsMechanicsFEM.h.

Constructor & Destructor Documentation

◆ ThermoRichardsMechanicsLocalAssembler() [1/3]

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::ThermoRichardsMechanicsLocalAssembler ( ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits > const & )
delete

References ThermoRichardsMechanicsLocalAssembler().

Referenced by ThermoRichardsMechanicsLocalAssembler(), ThermoRichardsMechanicsLocalAssembler(), addToLocalMatrixData(), assembleWithJacobianSingleIP(), and massLumping().

◆ ThermoRichardsMechanicsLocalAssembler() [2/3]

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::ThermoRichardsMechanicsLocalAssembler ( ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits > && )
delete

References ThermoRichardsMechanicsLocalAssembler().

◆ ThermoRichardsMechanicsLocalAssembler() [3/3]

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::ThermoRichardsMechanicsLocalAssembler ( MeshLib::Element const & e,
std::size_t const ,
NumLib::GenericIntegrationMethod const & integration_method,
bool const is_axially_symmetric,
ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > & process_data )

Definition at line 32 of file ThermoRichardsMechanicsFEM-impl.h.

40 : LocalAssemblerInterface<DisplacementDim, ConstitutiveTraits>(
41 e, integration_method, is_axially_symmetric, process_data)
42{
43 unsigned const n_integration_points =
44 this->integration_method_.getNumberOfPoints();
45
46 ip_data_.resize(n_integration_points);
47
48 auto const shape_matrices_u =
49 NumLib::initShapeMatrices<ShapeFunctionDisplacement,
50 ShapeMatricesTypeDisplacement,
51 DisplacementDim>(e, is_axially_symmetric,
52 integration_method);
53
54 auto const shape_matrices =
55 NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
56 DisplacementDim>(e, is_axially_symmetric,
57 integration_method);
58
59 for (unsigned ip = 0; ip < n_integration_points; ip++)
60 {
61 auto& ip_data = ip_data_[ip];
62 auto const& sm_u = shape_matrices_u[ip];
63 ip_data_[ip].integration_weight =
64 integration_method.getWeightedPoint(ip).getWeight() *
65 sm_u.integralMeasure * sm_u.detJ;
66
67 ip_data.N_u = sm_u.N;
68 ip_data.dNdx_u = sm_u.dNdx;
69
70 // ip_data.N_p and ip_data.dNdx_p are used for both p and T variables
71 ip_data.N_p = shape_matrices[ip].N;
72 ip_data.dNdx_p = shape_matrices[ip].dNdx;
73 }
74}
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler::ShapeMatricesTypeDisplacement
ShapeMatrixPolicyType< ShapeFunctionDisplacement, DisplacementDim > ShapeMatricesTypeDisplacement
Definition ThermoRichardsMechanicsFEM.h:43
ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface::LocalAssemblerInterface
LocalAssemblerInterface(MeshLib::Element const &e, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > &process_data)
Definition ThermoRichardsMechanics/LocalAssemblerInterface.h:23

References ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::LocalAssemblerInterface(), MathLib::WeightedPoint::getWeight(), NumLib::GenericIntegrationMethod::getWeightedPoint(), NumLib::initShapeMatrices(), ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::integration_method_, and ip_data_.

Member Function Documentation

◆ addToLocalMatrixData()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
void ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::addToLocalMatrixData ( double const dt,
std::vector< double > const & local_x,
std::vector< double > const & local_x_prev,
LocalMatrices const & loc_mat,
std::vector< double > & local_rhs_data,
std::vector< double > & local_Jac_data ) const
private

Definition at line 277 of file ThermoRichardsMechanicsFEM-impl.h.

287{
288 constexpr auto local_matrix_dim =
289 displacement_size + pressure_size + temperature_size;
290 assert(local_x.size() == local_matrix_dim);
291
292 auto local_Jac = MathLib::createZeroedMatrix<
293 typename ShapeMatricesTypeDisplacement::template MatrixType<
294 local_matrix_dim, local_matrix_dim>>(
295 local_Jac_data, local_matrix_dim, local_matrix_dim);
296
297 auto local_rhs =
298 MathLib::createZeroedVector<typename ShapeMatricesTypeDisplacement::
299 template VectorType<local_matrix_dim>>(
300 local_rhs_data, local_matrix_dim);
301
302 local_Jac.noalias() = loc_mat.Jac;
303 local_rhs.noalias() = -loc_mat.res;
304
305 //
306 // -- Jacobian
307 //
308 block_TT(local_Jac).noalias() += loc_mat.M_TT / dt + loc_mat.K_TT;
309 block_Tp(local_Jac).noalias() +=
310 loc_mat.K_Tp + loc_mat.dK_TT_dp + loc_mat.M_Tp / dt;
311
312 block_pT(local_Jac).noalias() += loc_mat.M_pT / dt + loc_mat.K_pT;
313 block_pp(local_Jac).noalias() +=
314 loc_mat.K_pp + loc_mat.storage_p_a_p / dt + loc_mat.storage_p_a_S_Jpp;
315 block_pu(local_Jac).noalias() = loc_mat.M_pu / dt;
316
317 //
318 // -- Residual
319 //
320 auto const [T, p_L, u] = localDOF(local_x);
321 auto const [T_prev, p_L_prev, u_prev] = localDOF(local_x_prev);
322
323 block_T(local_rhs).noalias() -= loc_mat.M_TT * (T - T_prev) / dt +
324 loc_mat.K_TT * T + loc_mat.K_Tp * p_L +
325 loc_mat.M_Tp * (p_L - p_L_prev) / dt;
326 block_p(local_rhs).noalias() -=
327 loc_mat.K_pp * p_L + loc_mat.K_pT * T +
328 (loc_mat.storage_p_a_p + loc_mat.storage_p_a_S) * (p_L - p_L_prev) /
329 dt +
330 loc_mat.M_pu * (u - u_prev) / dt + loc_mat.M_pT * (T - T_prev) / dt;
331}

References ThermoRichardsMechanicsLocalAssembler(), block_p(), block_pp(), block_pT(), block_pu(), block_T(), block_Tp(), block_TT(), MathLib::createZeroedMatrix(), MathLib::createZeroedVector(), displacement_size, localDOF(), pressure_size, and temperature_size.

Referenced by assembleWithJacobian().

◆ assembleWithJacobian()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
void ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::assembleWithJacobian ( double const t,
double const dt,
std::vector< double > const & local_x,
std::vector< double > const & local_x_prev,
std::vector< double > & local_rhs_data,
std::vector< double > & local_Jac_data )
overridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 209 of file ThermoRichardsMechanicsFEM-impl.h.

215{
216 auto& medium =
217 *this->process_data_.media_map.getMedium(this->element_.getID());
218
219 LocalMatrices loc_mat;
220 loc_mat.setZero();
221 LocalMatrices loc_mat_current_ip;
222 loc_mat_current_ip.setZero(); // only to set the right matrix sizes
223
224 typename ConstitutiveTraits::ConstitutiveSetting constitutive_setting;
225
226 for (unsigned ip = 0; ip < this->integration_method_.getNumberOfPoints();
227 ++ip)
228 {
229 ParameterLib::SpatialPosition const x_position{
230 std::nullopt, this->element_.getID(),
231 MathLib::Point3d(
232 NumLib::interpolateCoordinates<ShapeFunctionDisplacement,
233 ShapeMatricesTypeDisplacement>(
234 this->element_, ip_data_[ip].N_u))};
235
236 assembleWithJacobianSingleIP(
237 t, dt, x_position, //
238 local_x, local_x_prev, //
239 ip_data_[ip], constitutive_setting,
240 medium, //
241 loc_mat_current_ip, //
242 this->current_states_[ip], this->prev_states_[ip],
243 this->material_states_[ip], this->output_data_[ip]);
244 loc_mat += loc_mat_current_ip;
245 }
246
247 massLumping(loc_mat);
248
249 addToLocalMatrixData(dt, local_x, local_x_prev, loc_mat, local_rhs_data,
250 local_Jac_data);
251}
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler::addToLocalMatrixData
void addToLocalMatrixData(double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, LocalMatrices const &loc_mat, std::vector< double > &local_rhs_data, std::vector< double > &local_Jac_data) const
Definition ThermoRichardsMechanicsFEM-impl.h:278
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler::assembleWithJacobianSingleIP
void assembleWithJacobianSingleIP(double const t, double const dt, ParameterLib::SpatialPosition const &x_position, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, IpData const &ip_data, typename ConstitutiveTraits::ConstitutiveSetting &CS, MaterialPropertyLib::Medium &medium, LocalMatrices &out, typename ConstitutiveTraits::StatefulData &current_state, typename ConstitutiveTraits::StatefulDataPrev const &prev_state, MaterialStateData< DisplacementDim > &mat_state, typename ConstitutiveTraits::OutputData &output_data) const
Definition ThermoRichardsMechanicsFEM-impl.h:338
ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface::process_data_
ThermoRichardsMechanicsProcessData< DisplacementDim, ConstitutiveTraits > & process_data_
Definition ThermoRichardsMechanics/LocalAssemblerInterface.h:174

References addToLocalMatrixData(), assembleWithJacobianSingleIP(), ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::current_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::element_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::integration_method_, NumLib::interpolateCoordinates(), ip_data_, massLumping(), ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::material_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::output_data_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::prev_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::process_data_, and ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::LocalMatrices::setZero().

◆ assembleWithJacobianSingleIP()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
void ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::assembleWithJacobianSingleIP ( double const t,
double const dt,
ParameterLib::SpatialPosition const & x_position,
std::vector< double > const & local_x,
std::vector< double > const & local_x_prev,
IpData const & ip_data,
typename ConstitutiveTraits::ConstitutiveSetting & CS,
MaterialPropertyLib::Medium & medium,
LocalMatrices & out,
typename ConstitutiveTraits::StatefulData & current_state,
typename ConstitutiveTraits::StatefulDataPrev const & prev_state,
MaterialStateData< DisplacementDim > & mat_state,
typename ConstitutiveTraits::OutputData & output_data ) const
private

Definition at line 337 of file ThermoRichardsMechanicsFEM-impl.h.

355{
356 auto const& N_u = ip_data.N_u;
357 auto const& dNdx_u = ip_data.dNdx_u;
358
359 // N and dNdx are used for both p and T variables
360 auto const& N = ip_data.N_p;
361 auto const& dNdx = ip_data.dNdx_p;
362
363 auto const B =
364 LinearBMatrix::computeBMatrix<DisplacementDim,
365 ShapeFunctionDisplacement::NPOINTS,
366 typename BMatricesType::BMatrixType>(
367 dNdx_u, N_u, (*x_position.getCoordinates())[0],
368 this->is_axially_symmetric_);
369
370 typename ConstitutiveTraits::ConstitutiveData CD;
371
372 auto const [T, p_L, u] = localDOF(local_x);
373 auto const [T_prev, p_L_prev, u_prev] = localDOF(local_x_prev);
374
375 {
376 auto models = ConstitutiveTraits::createConstitutiveModels(
377 this->process_data_, this->solid_material_);
378 typename ConstitutiveTraits::ConstitutiveTempData tmp;
379
380 double const T_ip = N * T;
381 double const T_prev_ip = N * T_prev;
382 GlobalDimVectorType const grad_T_ip = dNdx * T;
383
384 double const p_cap_ip = -N * p_L;
385 double const p_cap_prev_ip = -N * p_L_prev;
386 GlobalDimVectorType const grad_p_cap_ip = -dNdx * p_L;
387
388 KelvinVectorType eps = B * u;
389
390 CS.eval(models, t, dt, x_position, //
391 medium, //
392 {T_ip, T_prev_ip, grad_T_ip}, //
393 {p_cap_ip, p_cap_prev_ip, grad_p_cap_ip}, //
394 eps, current_state, prev_state, mat_state, tmp, output_data,
395 CD);
396 }
397
398 using NodalMatrix = typename ShapeMatricesType::NodalMatrixType;
399 NodalMatrix const NTN = N.transpose() * N;
400 NodalMatrix const dNTdN = dNdx.transpose() * dNdx;
401
402 // TODO is identity2.transpose() * B something like divergence?
403 auto const& identity2 = MathLib::KelvinVector::Invariants<
404 MathLib::KelvinVector::kelvin_vector_dimensions(
405 DisplacementDim)>::identity2;
406 typename ShapeMatricesTypeDisplacement::template MatrixType<
407 displacement_size, pressure_size>
408 BTI2N = B.transpose() * identity2 * N;
409
410 /*
411 * Conventions:
412 *
413 * * use positive signs exclusively, any negative signs must be included in
414 * the coefficients coming from the constitutive setting
415 * * the used coefficients from the constitutive setting are named after the
416 * terms they appear in, e.g. K_TT_X_dNTdN means it is a scalar (X) that
417 * will be multiplied by dNdx.transpose() * dNdx. Placefolders for the
418 * coefficients are:
419 * * X -> scalar
420 * * V -> vector
421 * * K -> Kelvin vector
422 * * the Laplace terms have a special name, e.g., K_TT_Laplace
423 * * there shall be only one contribution to each of the LocalMatrices,
424 * assigned with = assignment; this point might be relaxed in the future
425 * * this method will overwrite the data in the passed LocalMatrices& out
426 * argument, not add to it
427 */
428
429 // residual, order T, p, u
430 block_p(out.res).noalias() =
431 dNdx.transpose() * std::get<EqPData<DisplacementDim>>(CD).rhs_p_dNT_V;
432 block_u(out.res).noalias() =
433 B.transpose() *
434 ProcessLib::Graph::get<TotalStressData<DisplacementDim>>(
435 CD, current_state)
436 .sigma_total -
437 static_cast<int>(this->process_data_.apply_body_force_for_deformation) *
438 N_u_op(N_u).transpose() *
439 std::get<GravityData<DisplacementDim>>(CD).volumetric_body_force;
440
441 // Storage matrices
442 out.storage_p_a_p.noalias() =
443 std::get<EqPData<DisplacementDim>>(CD).storage_p_a_p_X_NTN * NTN;
444 out.storage_p_a_S.noalias() =
445 std::get<TRMStorageData>(CD).storage_p_a_S_X_NTN * NTN;
446 out.storage_p_a_S_Jpp.noalias() =
447 std::get<TRMStorageData>(CD).storage_p_a_S_Jpp_X_NTN * NTN;
448
449 // M matrices, order T, p, u
450 out.M_TT.noalias() =
451 std::get<EqTData<DisplacementDim>>(CD).M_TT_X_NTN * NTN;
452 out.M_Tp.noalias() =
453 std::get<TRMVaporDiffusionData<DisplacementDim>>(CD).M_Tp_X_NTN * NTN;
454
455 out.M_pT.noalias() =
456 std::get<EqPData<DisplacementDim>>(CD).M_pT_X_NTN * NTN;
457 out.M_pu.noalias() =
458 std::get<EqPData<DisplacementDim>>(CD).M_pu_X_BTI2N * BTI2N.transpose();
459
460 // K matrices, order T, p, u
461 out.K_TT.noalias() =
462 dNdx.transpose() *
463 std::get<TRMHeatStorageAndFluxData<DisplacementDim>>(CD)
464 .K_TT_Laplace *
465 dNdx +
466 N.transpose() *
467 (std::get<EqTData<DisplacementDim>>(CD).K_TT_NT_V_dN.transpose() *
468 dNdx) +
469 std::get<TRMVaporDiffusionData<DisplacementDim>>(CD).K_TT_X_dNTdN *
470 dNTdN;
471
472 out.dK_TT_dp.noalias() =
473 N.transpose() *
474 (std::get<TRMHeatStorageAndFluxData<DisplacementDim>>(CD)
475 .K_Tp_NT_V_dN.transpose() *
476 dNdx) +
477 std::get<TRMHeatStorageAndFluxData<DisplacementDim>>(CD).K_Tp_X_NTN *
478 NTN;
479 out.K_Tp.noalias() =
480 dNdx.transpose() *
481 std::get<ThermoOsmosisData<DisplacementDim>>(CD).K_Tp_Laplace *
482 dNdx +
483 std::get<TRMVaporDiffusionData<DisplacementDim>>(CD).K_Tp_X_dNTdN *
484 dNTdN;
485
486 out.K_pp.noalias() =
487 dNdx.transpose() * std::get<EqPData<DisplacementDim>>(CD).K_pp_Laplace *
488 dNdx +
489 std::get<TRMVaporDiffusionData<DisplacementDim>>(CD).K_pp_X_dNTdN *
490 dNTdN;
491 out.K_pT.noalias() =
492 dNdx.transpose() *
493 std::get<ThermoOsmosisData<DisplacementDim>>(CD).K_pT_Laplace * dNdx;
494
495 // direct Jacobian contributions, order T, p, u
496 block_pT(out.Jac).noalias() =
497 std::get<TRMVaporDiffusionData<DisplacementDim>>(CD).J_pT_X_dNTdN *
498 dNTdN;
499 block_pp(out.Jac).noalias() =
500 std::get<TRMStorageData>(CD).J_pp_X_NTN * NTN +
501 std::get<EqPData<DisplacementDim>>(CD).J_pp_X_BTI2NT_u_dot_N *
502 BTI2N.transpose() * (u - u_prev) / dt *
503 N // TODO something with volumetric strain rate?
504 + dNdx.transpose() *
505 std::get<EqPData<DisplacementDim>>(CD).J_pp_dNT_V_N * N;
506
507 block_uT(out.Jac).noalias() =
508 B.transpose() *
509 std::get<SolidMechanicsDataStateless<DisplacementDim>>(CD).J_uT_BT_K_N *
510 N;
511 block_up(out.Jac).noalias() =
512 B.transpose() *
513 std::get<SolidMechanicsDataStateless<DisplacementDim>>(CD)
514 .J_up_BT_K_N *
515 N +
516 N_u_op(N_u).transpose() *
517 std::get<GravityData<DisplacementDim>>(CD).J_up_HT_V_N * N;
518 block_uu(out.Jac).noalias() =
519 B.transpose() *
520 std::get<SolidMechanicsDataStateless<DisplacementDim>>(CD)
521 .stiffness_tensor *
522 B;
523
524 out *= ip_data.integration_weight;
525}
ProcessLib::BMatrixPolicyType< ShapeFunctionDisplacement, DisplacementDim >::BMatrixType
MatrixType< _kelvin_vector_size, _number_of_dof > BMatrixType
Definition BMatrixPolicy.h:48
MathLib::KelvinVector::kelvin_vector_dimensions
constexpr int kelvin_vector_dimensions(int const displacement_dim)
Kelvin vector dimensions for given displacement dimension.
Definition KelvinVector.h:18
EigenFixedShapeMatrixPolicy::NodalMatrixType
MatrixType< ShapeFunction::NPOINTS, ShapeFunction::NPOINTS > NodalMatrixType
Definition ShapeMatrixPolicy.h:67

References ThermoRichardsMechanicsLocalAssembler(), block_p(), block_pp(), block_pT(), block_u(), block_up(), block_uT(), block_uu(), ProcessLib::LinearBMatrix::computeBMatrix(), displacement_size, ProcessLib::Graph::get(), ParameterLib::SpatialPosition::getCoordinates(), ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::is_axially_symmetric_, ProcessLib::ThermoRichardsMechanics::TRMVaporDiffusionData< DisplacementDim >::K_pp_X_dNTdN, ProcessLib::ThermoRichardsMechanics::TRMVaporDiffusionData< DisplacementDim >::K_Tp_X_dNTdN, MathLib::KelvinVector::kelvin_vector_dimensions(), localDOF(), N_u_op, pressure_size, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::process_data_, and ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::solid_material_.

Referenced by assembleWithJacobian().

◆ block_p()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_p ( auto & vec)
inlinestaticprivate

Definition at line 409 of file ThermoRichardsMechanicsFEM.h.

410 {
411 return vec.template segment<pressure_size>(pressure_index);
412 }

References pressure_index.

Referenced by addToLocalMatrixData(), assembleWithJacobianSingleIP(), and computeSecondaryVariableConcrete().

◆ block_pp()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_pp ( auto & mat)
inlinestaticprivate

Definition at line 384 of file ThermoRichardsMechanicsFEM.h.

385 {
386 return mat.template block<pressure_size, pressure_size>(pressure_index,
387 pressure_index);
388 }

References pressure_index.

Referenced by addToLocalMatrixData(), and assembleWithJacobianSingleIP().

◆ block_pT()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_pT ( auto & mat)
inlinestaticprivate

Definition at line 389 of file ThermoRichardsMechanicsFEM.h.

390 {
391 return mat.template block<pressure_size, temperature_size>(
392 pressure_index, temperature_index);
393 }

References pressure_index, and temperature_index.

Referenced by addToLocalMatrixData(), and assembleWithJacobianSingleIP().

◆ block_pu()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_pu ( auto & mat)
inlinestaticprivate

Definition at line 379 of file ThermoRichardsMechanicsFEM.h.

380 {
381 return mat.template block<pressure_size, displacement_size>(
382 pressure_index, displacement_index);
383 }

References displacement_index, and pressure_index.

Referenced by addToLocalMatrixData().

◆ block_T()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_T ( auto & vec)
inlinestaticprivate

Definition at line 413 of file ThermoRichardsMechanicsFEM.h.

414 {
415 return vec.template segment<temperature_size>(temperature_index);
416 }

References temperature_index.

Referenced by addToLocalMatrixData(), and computeSecondaryVariableConcrete().

◆ block_Tp()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_Tp ( auto & mat)
inlinestaticprivate

Definition at line 394 of file ThermoRichardsMechanicsFEM.h.

395 {
396 return mat.template block<temperature_size, pressure_size>(
397 temperature_index, pressure_index);
398 }

References pressure_index, and temperature_index.

Referenced by addToLocalMatrixData().

◆ block_TT()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_TT ( auto & mat)
inlinestaticprivate

Definition at line 399 of file ThermoRichardsMechanicsFEM.h.

400 {
401 return mat.template block<temperature_size, temperature_size>(
402 temperature_index, temperature_index);
403 }

References temperature_index.

Referenced by addToLocalMatrixData().

◆ block_u()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_u ( auto & vec)
inlinestaticprivate

Definition at line 405 of file ThermoRichardsMechanicsFEM.h.

406 {
407 return vec.template segment<displacement_size>(displacement_index);
408 }

References displacement_index.

Referenced by assembleWithJacobianSingleIP(), and computeSecondaryVariableConcrete().

◆ block_up()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_up ( auto & mat)
inlinestaticprivate

Definition at line 369 of file ThermoRichardsMechanicsFEM.h.

370 {
371 return mat.template block<displacement_size, pressure_size>(
372 displacement_index, pressure_index);
373 }

References displacement_index, and pressure_index.

Referenced by assembleWithJacobianSingleIP().

◆ block_uT()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_uT ( auto & mat)
inlinestaticprivate

Definition at line 374 of file ThermoRichardsMechanicsFEM.h.

375 {
376 return mat.template block<displacement_size, temperature_size>(
377 displacement_index, temperature_index);
378 }

References displacement_index, and temperature_index.

Referenced by assembleWithJacobianSingleIP().

◆ block_uu()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::block_uu ( auto & mat)
inlinestaticprivate

Definition at line 364 of file ThermoRichardsMechanicsFEM.h.

365 {
366 return mat.template block<displacement_size, displacement_size>(
367 displacement_index, displacement_index);
368 }

References displacement_index.

Referenced by assembleWithJacobianSingleIP().

◆ computeSecondaryVariableConcrete()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
void ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::computeSecondaryVariableConcrete ( double const t,
double const dt,
Eigen::VectorXd const & local_x,
Eigen::VectorXd const & local_x_prev )
overridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 531 of file ThermoRichardsMechanicsFEM-impl.h.

535{
536 auto const T = block_T(local_x);
537 auto const p_L = block_p(local_x);
538 auto const u = block_u(local_x);
539
540 auto const T_prev = block_T(local_x_prev);
541 auto const p_L_prev = block_p(local_x_prev);
542
543 auto const e_id = this->element_.getID();
544 auto const& process_data = this->process_data_;
545 auto& medium = *process_data.media_map.getMedium(e_id);
546
547 unsigned const n_integration_points =
548 this->integration_method_.getNumberOfPoints();
549
550 typename ConstitutiveTraits::ConstitutiveSetting constitutive_setting;
551
552 auto models = ConstitutiveTraits::createConstitutiveModels(
553 process_data, this->solid_material_);
554 typename ConstitutiveTraits::ConstitutiveTempData tmp;
555 typename ConstitutiveTraits::ConstitutiveData CD;
556
557 for (unsigned ip = 0; ip < n_integration_points; ip++)
558 {
559 auto& current_state = this->current_states_[ip];
560 auto& output_data = this->output_data_[ip];
561
562 auto const& ip_data = ip_data_[ip];
563
564 // N is used for both p and T variables
565 auto const& N = ip_data.N_p;
566 auto const& N_u = ip_data.N_u;
567 auto const& dNdx_u = ip_data.dNdx_u;
568 auto const& dNdx = ip_data.dNdx_p;
569
570 ParameterLib::SpatialPosition const x_position{
571 std::nullopt, this->element_.getID(),
572 MathLib::Point3d(
573 NumLib::interpolateCoordinates<ShapeFunctionDisplacement,
574 ShapeMatricesTypeDisplacement>(
575 this->element_, N_u))};
576
577 auto const x_coord =
578 x_position.getCoordinates().value()[0]; // r for axisymmetry
579 auto const B =
580 LinearBMatrix::computeBMatrix<DisplacementDim,
581 ShapeFunctionDisplacement::NPOINTS,
582 typename BMatricesType::BMatrixType>(
583 dNdx_u, N_u, x_coord, this->is_axially_symmetric_);
584
585 double const T_ip = N * T;
586 double const T_prev_ip = N * T_prev;
587 GlobalDimVectorType const grad_T_ip = dNdx * T;
588
589 double const p_cap_ip = -N * p_L;
590 double const p_cap_prev_ip = -N * p_L_prev;
591 GlobalDimVectorType const grad_p_cap_ip = -dNdx * p_L;
592
593 KelvinVectorType eps = B * u;
594
595 constitutive_setting.eval(models, //
596 t, dt, x_position, //
597 medium, //
598 {T_ip, T_prev_ip, grad_T_ip}, //
599 {p_cap_ip, p_cap_prev_ip, grad_p_cap_ip}, //
600 eps, current_state, this->prev_states_[ip],
601 this->material_states_[ip], tmp, output_data,
602 CD);
603 }
604
605 NumLib::interpolateToHigherOrderNodes<
606 ShapeFunction, typename ShapeFunctionDisplacement::MeshElement,
607 DisplacementDim>(this->element_, this->is_axially_symmetric_, p_L,
608 *process_data.pressure_interpolated);
609 NumLib::interpolateToHigherOrderNodes<
610 ShapeFunction, typename ShapeFunctionDisplacement::MeshElement,
611 DisplacementDim>(this->element_, this->is_axially_symmetric_, T,
612 *process_data.temperature_interpolated);
613}

References block_p(), block_T(), block_u(), ProcessLib::LinearBMatrix::computeBMatrix(), ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::current_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::element_, ParameterLib::SpatialPosition::getCoordinates(), ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::integration_method_, NumLib::interpolateCoordinates(), NumLib::interpolateToHigherOrderNodes(), ip_data_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::is_axially_symmetric_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::material_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::output_data_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::prev_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::process_data_, and ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::solid_material_.

◆ convertInitialStressType()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
void ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::convertInitialStressType ( unsigned const ip,
double const t,
ParameterLib::SpatialPosition const x_position,
MaterialPropertyLib::Medium const & medium,
MPL::VariableArray const & variables,
double const p_at_ip )
private

This function resets the initial stress type according to the input initial stress type, either total or effective. If subtype = StressSaturation_StrainPressureTemperature is bing used in the process setting, the initial effective stress is converted to total stress. Otherwise, the initial total stress is converted to effective stress.

Definition at line 167 of file ThermoRichardsMechanicsFEM-impl.h.

174{
175 bool constexpr is_strain_temperature_constitutive =
176 std::is_same<ConstitutiveStress_StrainTemperature::ConstitutiveTraits<
177 DisplacementDim>,
178 ConstitutiveTraits>::value;
179 if (is_strain_temperature_constitutive &&
180 this->process_data_.initial_stress.type ==
181 InitialStress::Type::Effective)
182 {
183 return;
184 }
185
186 if (!is_strain_temperature_constitutive &&
187 this->process_data_.initial_stress.type == InitialStress::Type::Total)
188 {
189 return;
190 }
191
192 double const alpha_b =
193 medium.property(MPL::PropertyType::biot_coefficient)
194 .template value<double>(variables, x_position, t, 0.0 /*dt*/);
195
196 double const bishop =
197 medium.property(MPL::PropertyType::bishops_effective_stress)
198 .template value<double>(variables, x_position, t, 0.0 /*dt*/);
199
200 ConstitutiveTraits::ConstitutiveSetting::convertInitialStressType(
201 this->current_states_[ip], this->prev_states_[ip],
202 bishop * alpha_b * p_at_ip * Invariants::identity2);
203}
MathLib::KelvinVector::Invariants::identity2
static Eigen::Matrix< double, KelvinVectorSize, 1 > const identity2
Kelvin mapping of 2nd order identity tensor.
Definition KelvinVector.h:101

References MaterialPropertyLib::biot_coefficient, MaterialPropertyLib::bishops_effective_stress, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::current_states_, ProcessLib::InitialStress::Effective, MathLib::KelvinVector::Invariants< KelvinVectorSize >::identity2, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::prev_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::process_data_, MaterialPropertyLib::Medium::property(), and ProcessLib::InitialStress::Total.

Referenced by setInitialConditionsConcrete().

◆ getNumberOfVectorElementsForDeformation()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
int ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::getNumberOfVectorElementsForDeformation ( ) const
inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 244 of file ThermoRichardsMechanicsFEM.h.

245 {
246 return displacement_size;
247 }

References displacement_size.

◆ getShapeMatrix()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
Eigen::Map< const Eigen::RowVectorXd > ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::getShapeMatrix ( const unsigned integration_point) const
inlineoverridevirtual

Provides the shape matrix at the given integration point.

Implements NumLib::ExtrapolatableElement.

Definition at line 345 of file ThermoRichardsMechanicsFEM.h.

347 {
348 auto const& N_u = ip_data_[integration_point].N_u;
349
350 // assumes N is stored contiguously in memory
351 return Eigen::Map<const Eigen::RowVectorXd>(N_u.data(), N_u.size());
352 }

References ip_data_.

◆ initializeConcrete()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
void ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::initializeConcrete ( )
inlineoverridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 272 of file ThermoRichardsMechanicsFEM.h.

273 {
274 unsigned const n_integration_points =
275 this->integration_method_.getNumberOfPoints();
276 auto const time_independent = std::numeric_limits<double>::quiet_NaN();
277 auto const& medium =
278 *this->process_data_.media_map.getMedium(this->element_.getID());
279
280 for (unsigned ip = 0; ip < n_integration_points; ip++)
281 {
282 ParameterLib::SpatialPosition const x_position{
283 std::nullopt, this->element_.getID(),
284 MathLib::Point3d(NumLib::interpolateCoordinates<
285 ShapeFunctionDisplacement,
286 ShapeMatricesTypeDisplacement>(
287 this->element_, this->ip_data_[ip].N_u))};
288 auto& current_state = this->current_states_[ip];
289
290 // Set initial stress from parameter.
291 if (this->process_data_.initial_stress.value)
292 {
293 ConstitutiveTraits::ConstitutiveSetting::statefulStress(
294 current_state) =
295 MathLib::KelvinVector::symmetricTensorToKelvinVector<
296 DisplacementDim>(
297 (*this->process_data_.initial_stress.value)(
298 time_independent, x_position));
299 }
300
301 if (*this->process_data_.initialize_porosity_from_medium_property)
302 {
303 // Initial porosity. Could be read from integration point data
304 // or mesh.
305 std::get<PorosityData>(current_state).phi =
306 medium.property(MPL::porosity)
307 .template initialValue<double>(x_position,
308 time_independent);
309
310 if (medium.hasProperty(MPL::PropertyType::transport_porosity))
311 {
312 std::get<TransportPorosityData>(current_state).phi =
313 medium.property(MPL::transport_porosity)
314 .template initialValue<double>(x_position,
315 time_independent);
316 }
317 else
318 {
319 std::get<TransportPorosityData>(current_state).phi =
320 std::get<PorosityData>(current_state).phi;
321 }
322 }
323
324 double const t = 0; // TODO (naumov) pass t from top
325 this->solid_material_.initializeInternalStateVariables(
326 t, x_position,
327 *this->material_states_[ip].material_state_variables);
328 }
329
330 for (unsigned ip = 0; ip < n_integration_points; ip++)
331 {
332 this->material_states_[ip].pushBackState();
333 }
334
335 for (unsigned ip = 0; ip < n_integration_points; ip++)
336 {
337 this->prev_states_[ip] = this->current_states_[ip];
338 }
339 }

References ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::current_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::element_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::integration_method_, NumLib::interpolateCoordinates(), ip_data_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::material_states_, MaterialPropertyLib::porosity, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::prev_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::process_data_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::solid_material_, MathLib::KelvinVector::symmetricTensorToKelvinVector(), and MaterialPropertyLib::transport_porosity.

◆ localDOF()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
constexpr auto ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::localDOF ( std::vector< double > const & x)
inlinestaticconstexprprivate

Definition at line 357 of file ThermoRichardsMechanicsFEM.h.

358 {
359 return NumLib::localDOF<
360 ShapeFunction, ShapeFunction,
361 NumLib::Vectorial<ShapeFunctionDisplacement, DisplacementDim>>(x);
362 }

References NumLib::localDOF().

Referenced by addToLocalMatrixData(), and assembleWithJacobianSingleIP().

◆ massLumping()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
void ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::massLumping ( LocalMatrices & loc_mat) const
private

Definition at line 257 of file ThermoRichardsMechanicsFEM-impl.h.

261{
262 if (this->process_data_.apply_mass_lumping)
263 {
264 loc_mat.storage_p_a_p =
265 loc_mat.storage_p_a_p.colwise().sum().eval().asDiagonal();
266 loc_mat.storage_p_a_S =
267 loc_mat.storage_p_a_S.colwise().sum().eval().asDiagonal();
268 loc_mat.storage_p_a_S_Jpp =
269 loc_mat.storage_p_a_S_Jpp.colwise().sum().eval().asDiagonal();
270 }
271}

References ThermoRichardsMechanicsLocalAssembler(), and ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::process_data_.

Referenced by assembleWithJacobian().

◆ setInitialConditionsConcrete()

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
void ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::setInitialConditionsConcrete ( Eigen::VectorXd const local_x,
double const t,
int const process_id )
overridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 80 of file ThermoRichardsMechanicsFEM-impl.h.

84{
85 assert(local_x.size() ==
86 temperature_size + pressure_size + displacement_size);
87
88 auto const p_L = local_x.template segment<pressure_size>(pressure_index);
89 auto const T =
90 local_x.template segment<temperature_size>(temperature_index);
91
92 constexpr double dt = std::numeric_limits<double>::quiet_NaN();
93 auto const& medium =
94 *this->process_data_.media_map.getMedium(this->element_.getID());
95
96 MediaData const media_data{medium};
97
98 typename ConstitutiveTraits::ConstitutiveSetting const constitutive_setting;
99 auto models = ConstitutiveTraits::createConstitutiveModels(
100 this->process_data_, this->solid_material_);
101
102 unsigned const n_integration_points =
103 this->integration_method_.getNumberOfPoints();
104 for (unsigned ip = 0; ip < n_integration_points; ip++)
105 {
106 // N is used for both T and p variables.
107 auto const& N = ip_data_[ip].N_p;
108
109 ParameterLib::SpatialPosition const x_position{
110 std::nullopt, this->element_.getID(),
111 MathLib::Point3d(
112 NumLib::interpolateCoordinates<ShapeFunctionDisplacement,
113 ShapeMatricesTypeDisplacement>(
114 this->element_, ip_data_[ip].N_u))};
115
116 double T_ip;
117 NumLib::shapeFunctionInterpolate(T, N, T_ip);
118
119 double p_cap_ip;
120 NumLib::shapeFunctionInterpolate(-p_L, N, p_cap_ip);
121
122 MPL::VariableArray variables;
123
124 if (medium.hasProperty(MPL::PropertyType::saturation))
125 {
126 variables.capillary_pressure = p_cap_ip;
127 variables.liquid_phase_pressure = -p_cap_ip;
128 // setting pG to 1 atm
129 // TODO : rewrite equations s.t. p_L = pG-p_cap
130 variables.gas_phase_pressure = 1.0e5;
131 variables.temperature = T_ip;
132
133 double const S_L =
134 medium.property(MPL::PropertyType::saturation)
135 .template value<double>(variables, x_position, t, dt);
136 std::get<PrevState<SaturationData>>(this->prev_states_[ip])->S_L =
137 S_L;
138
139 variables.liquid_saturation = S_L;
140 }
141
142 // dt = 0 at initialization: there is no time step yet, which yields
143 // the elastic tangent. The function-wide dt is NaN to keep
144 // initialization and integration strictly separated.
145 constitutive_setting.init(models, t, 0.0 /*dt*/, x_position, media_data,
146 {T_ip, 0, {}}, this->current_states_[ip],
147 this->prev_states_[ip]);
148
149 if (this->process_data_.initial_stress.value)
150 {
151 if (!medium.hasProperty(MPL::PropertyType::saturation))
152 {
153 OGS_FATAL(
154 "The medium has no saturation property required to compute "
155 "initial stress.");
156 }
157 convertInitialStressType(ip, t, x_position, medium, variables,
158 -p_cap_ip);
159 }
160 }
161}
OGS_FATAL
#define OGS_FATAL(...)
Definition Error.h:10
ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler::convertInitialStressType
void convertInitialStressType(unsigned const ip, double const t, ParameterLib::SpatialPosition const x_position, MaterialPropertyLib::Medium const &medium, MPL::VariableArray const &variables, double const p_at_ip)
Definition ThermoRichardsMechanicsFEM-impl.h:168
NumLib::detail::shapeFunctionInterpolate
void shapeFunctionInterpolate(const NodalValues &, const ShapeMatrix &)
Definition Interpolation.h:20

References MaterialPropertyLib::VariableArray::capillary_pressure, convertInitialStressType(), ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::current_states_, displacement_size, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::element_, MaterialPropertyLib::VariableArray::gas_phase_pressure, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::integration_method_, NumLib::interpolateCoordinates(), ip_data_, MaterialPropertyLib::VariableArray::liquid_phase_pressure, MaterialPropertyLib::VariableArray::liquid_saturation, OGS_FATAL, pressure_index, pressure_size, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::prev_states_, ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::process_data_, MaterialPropertyLib::saturation, setInitialConditionsConcrete(), NumLib::detail::shapeFunctionInterpolate(), ProcessLib::ThermoRichardsMechanics::LocalAssemblerInterface< DisplacementDim, ConstitutiveTraits >::solid_material_, MaterialPropertyLib::VariableArray::temperature, temperature_index, and temperature_size.

Referenced by setInitialConditionsConcrete().

Member Data Documentation

◆ displacement_index

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
int ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::displacement_index = 2 * ShapeFunction::NPOINTS
staticconstexprprivate

Definition at line 38 of file ThermoRichardsMechanicsFEM.h.

Referenced by block_pu(), block_u(), block_up(), block_uT(), and block_uu().

◆ displacement_size

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
int ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::displacement_size
staticconstexprprivate
Initial value:
=
ShapeFunctionDisplacement::NPOINTS * DisplacementDim

Definition at line 39 of file ThermoRichardsMechanicsFEM.h.

Referenced by addToLocalMatrixData(), assembleWithJacobianSingleIP(), getNumberOfVectorElementsForDeformation(), setInitialConditionsConcrete(), and ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::LocalMatrices::setZero().

◆ ip_data_

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
std::vector<IpData> ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::ip_data_
private

Definition at line 355 of file ThermoRichardsMechanicsFEM.h.

Referenced by ThermoRichardsMechanicsLocalAssembler(), assembleWithJacobian(), computeSecondaryVariableConcrete(), getShapeMatrix(), initializeConcrete(), and setInitialConditionsConcrete().

◆ KelvinVectorSize

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
int const ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::KelvinVectorSize
static
Initial value:
=
MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim)

Definition at line 61 of file ThermoRichardsMechanicsFEM.h.

◆ N_u_op

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
auto& ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::N_u_op
staticconstexpr
Initial value:
= MathLib::eigenBlockMatrixView<
DisplacementDim,
typename ShapeMatricesTypeDisplacement::NodalRowVectorType>
MathLib::eigenBlockMatrixView
constexpr Eigen::CwiseNullaryOp< EigenBlockMatrixViewFunctor< D, M >, typename EigenBlockMatrixViewFunctor< D, M >::Matrix > eigenBlockMatrixView(const Eigen::MatrixBase< M > &matrix)
Definition EigenBlockMatrixView.h:43
EigenFixedShapeMatrixPolicy::NodalRowVectorType
RowVectorType< ShapeFunction::NPOINTS > NodalRowVectorType
Definition ShapeMatrixPolicy.h:71

Definition at line 67 of file ThermoRichardsMechanicsFEM.h.

Referenced by assembleWithJacobianSingleIP().

◆ pressure_index

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
int ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::pressure_index = temperature_size
staticconstexprprivate

Definition at line 36 of file ThermoRichardsMechanicsFEM.h.

Referenced by block_p(), block_pp(), block_pT(), block_pu(), block_Tp(), block_up(), and setInitialConditionsConcrete().

◆ pressure_size

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
int ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::pressure_size = ShapeFunction::NPOINTS
staticconstexprprivate

Definition at line 37 of file ThermoRichardsMechanicsFEM.h.

Referenced by addToLocalMatrixData(), assembleWithJacobianSingleIP(), setInitialConditionsConcrete(), and ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::LocalMatrices::setZero().

◆ temperature_index

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
int ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::temperature_index = 0
staticconstexprprivate

Definition at line 34 of file ThermoRichardsMechanicsFEM.h.

Referenced by block_pT(), block_T(), block_Tp(), block_TT(), block_uT(), and setInitialConditionsConcrete().

◆ temperature_size

template<typename ShapeFunctionDisplacement, typename ShapeFunction, int DisplacementDim, typename ConstitutiveTraits>
int ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::temperature_size = ShapeFunction::NPOINTS
staticconstexprprivate

Definition at line 35 of file ThermoRichardsMechanicsFEM.h.

Referenced by addToLocalMatrixData(), setInitialConditionsConcrete(), and ProcessLib::ThermoRichardsMechanics::ThermoRichardsMechanicsLocalAssembler< ShapeFunctionDisplacement, ShapeFunction, DisplacementDim, ConstitutiveTraits >::LocalMatrices::setZero().

The documentation for this class was generated from the following files: