Loading [MathJax]/jax/output/HTML-CSS/config.js
OGS
IntegrationPointData.h
Go to the documentation of this file.
1
11#pragma once
12
13#include <memory>
14
20
21namespace ProcessLib
22{
23namespace RichardsMechanics
24{
25template <typename BMatricesType, typename ShapeMatrixTypeDisplacement,
26 typename ShapeMatricesTypePressure, int DisplacementDim, int NPoints>
28{
29 typename ShapeMatrixTypeDisplacement::NodalRowVectorType N_u;
30 typename ShapeMatrixTypeDisplacement::GlobalDimNodalMatrixType dNdx_u;
31
32 typename ShapeMatricesTypePressure::NodalRowVectorType N_p;
33 typename ShapeMatricesTypePressure::GlobalDimNodalMatrixType dNdx_p;
34
35 double integration_weight = std::numeric_limits<double>::quiet_NaN();
36
37 MathLib::KelvinVector::
38 KelvinMatrixType<DisplacementDim> static computeElasticTangentStiffness(
39 double const t,
40 ParameterLib::SpatialPosition const& x_position,
41 double const dt,
42 double const temperature,
44 solid_material,
46 MaterialStateVariables const& material_state_variables)
47 {
48 namespace MPL = MaterialPropertyLib;
49
50 MPL::VariableArray variable_array;
51 MPL::VariableArray variable_array_prev;
52
54
55 variable_array.stress.emplace<KV>(KV::Zero());
56 variable_array.mechanical_strain.emplace<KV>(KV::Zero());
57 variable_array.temperature = temperature;
58
59 variable_array_prev.stress.emplace<KV>(KV::Zero());
60 variable_array_prev.mechanical_strain.emplace<KV>(KV::Zero());
61 variable_array_prev.temperature = temperature;
62
63 auto&& solution = solid_material.integrateStress(
64 variable_array_prev, variable_array, t, x_position, dt,
65 material_state_variables);
66
67 if (!solution)
68 {
69 OGS_FATAL("Computation of elastic tangent stiffness failed.");
70 }
71
73 std::move(std::get<2>(*solution));
74
75 return C;
76 }
77
78 static typename BMatricesType::KelvinMatrixType updateConstitutiveRelation(
79 MaterialPropertyLib::VariableArray const& variable_array,
80 double const t,
81 ParameterLib::SpatialPosition const& x_position,
82 double const dt,
83 double const temperature,
85 sigma_eff,
87 DisplacementDim>> const& sigma_eff_prev,
89 DisplacementDim> const&
90 /*eps_m*/,
92 DisplacementDim>> const& eps_m_prev,
94 solid_material,
95 std::unique_ptr<typename MaterialLib::Solids::MechanicsBase<
96 DisplacementDim>::MaterialStateVariables>& material_state_variables)
97 {
98 MaterialPropertyLib::VariableArray variable_array_prev;
99 variable_array_prev.stress = sigma_eff_prev->sigma_eff;
100 variable_array_prev.mechanical_strain
102 eps_m_prev->eps_m);
103 variable_array_prev.temperature = temperature;
104
105 auto&& solution = solid_material.integrateStress(
106 variable_array_prev, variable_array, t, x_position, dt,
107 *material_state_variables);
108
109 if (!solution)
110 {
111 OGS_FATAL("Computation of local constitutive relation failed.");
112 }
113
115 std::tie(sigma_eff.sigma_eff, material_state_variables, C) =
116 std::move(*solution);
117
118 return C;
119 }
120
122};
123
124} // namespace RichardsMechanics
125} // namespace ProcessLib
#define OGS_FATAL(...)
Definition Error.h:26
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), 1, Eigen::ColMajor > KelvinVectorType
Eigen::Matrix< double, kelvin_vector_dimensions(DisplacementDim), kelvin_vector_dimensions(DisplacementDim), Eigen::RowMajor > KelvinMatrixType
virtual std::optional< std::tuple< KelvinVector, std::unique_ptr< MaterialStateVariables >, KelvinMatrix > > integrateStress(MaterialPropertyLib::VariableArray const &variable_array_prev, MaterialPropertyLib::VariableArray const &variable_array, double const t, ParameterLib::SpatialPosition const &x, double const dt, MaterialStateVariables const &material_state_variables) const =0
MathLib::KelvinVector::KelvinVectorType< DisplacementDim > sigma_eff
Represents a previous state of type T.
Definition Base.h:21
ShapeMatrixTypeDisplacement::NodalRowVectorType N_u
ShapeMatricesTypePressure::GlobalDimNodalMatrixType dNdx_p
ShapeMatrixTypeDisplacement::GlobalDimNodalMatrixType dNdx_u
ShapeMatricesTypePressure::NodalRowVectorType N_p
static BMatricesType::KelvinMatrixType updateConstitutiveRelation(MaterialPropertyLib::VariableArray const &variable_array, double const t, ParameterLib::SpatialPosition const &x_position, double const dt, double const temperature, ProcessLib::ConstitutiveRelations::EffectiveStressData< DisplacementDim > &sigma_eff, PrevState< ProcessLib::ConstitutiveRelations::EffectiveStressData< DisplacementDim > > const &sigma_eff_prev, ProcessLib::ConstitutiveRelations::MechanicalStrainData< DisplacementDim > const &, PrevState< ProcessLib::ConstitutiveRelations::MechanicalStrainData< DisplacementDim > > const &eps_m_prev, MaterialLib::Solids::MechanicsBase< DisplacementDim > const &solid_material, std::unique_ptr< typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables > &material_state_variables)
static MathLib::KelvinVector::KelvinMatrixType< DisplacementDim > computeElasticTangentStiffness(double const t, ParameterLib::SpatialPosition const &x_position, double const dt, double const temperature, MaterialLib::Solids::MechanicsBase< DisplacementDim > const &solid_material, typename MaterialLib::Solids::MechanicsBase< DisplacementDim >::MaterialStateVariables const &material_state_variables)