OGS
ProcessLib::TH2M::ConstitutiveRelations::SwellingModel< DisplacementDim > Struct Template Reference

Detailed Description

template<int DisplacementDim>
struct ProcessLib::TH2M::ConstitutiveRelations::SwellingModel< DisplacementDim >

Definition at line 35 of file TH2M/ConstitutiveRelations/Swelling.h.

#include <Swelling.h>

Public Member Functions

void eval (SpaceTimeData const &x_t, MediaData const &media_data, ElasticTangentStiffnessData< DisplacementDim > const &C_el_data, SaturationData const &S_L_data, PrevState< SaturationData > const &S_L_prev_data, PrevState< SwellingDataStateful< DisplacementDim > > const &prev_state, SwellingDataStateful< DisplacementDim > &state, SwellingDataStateless< DisplacementDim > &out) const

Member Function Documentation

◆ eval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::SwellingModel< DisplacementDim >::eval ( SpaceTimeData const & x_t,
MediaData const & media_data,
ElasticTangentStiffnessData< DisplacementDim > const & C_el_data,
SaturationData const & S_L_data,
PrevState< SaturationData > const & S_L_prev_data,
PrevState< SwellingDataStateful< DisplacementDim > > const & prev_state,
SwellingDataStateful< DisplacementDim > & state,
SwellingDataStateless< DisplacementDim > & out ) const

Definition at line 15 of file TH2M/ConstitutiveRelations/Swelling.cpp.

23{
24 namespace MPL = MaterialPropertyLib;
25
26 if (media_data.swelling_stress_rate_solid == nullptr)
27 {
28 out.eps_m.setZero();
29 state.sigma_sw.setZero();
30 return;
31 }
32
33 MPL::VariableArray variables;
34 variables.liquid_saturation = S_L_data.S_L;
35
36 MPL::VariableArray variables_prev;
37 variables_prev.liquid_saturation = S_L_prev_data->S_L;
38
39 // Swelling and possibly volumetric strain rate update.
40
41 // If there is swelling, compute it. Update volumetric strain rate,
42 // s.t. it corresponds to the mechanical part only.
43 state.sigma_sw = prev_state->sigma_sw;
44
45 auto const sigma_sw_dot =
46 MathLib::KelvinVector::tensorToKelvin<DisplacementDim>(
47 MPL::formEigenTensor<3>(
48 media_data.swelling_stress_rate_solid->value(
49 variables, variables_prev, x_t.x, x_t.t, x_t.dt)));
50 state.sigma_sw += sigma_sw_dot * x_t.dt;
51
52 auto const C_el_inv = C_el_data.stiffness_tensor.inverse().eval();
53
54 out.eps_m.noalias() = C_el_inv * (state.sigma_sw - prev_state->sigma_sw);
55}
MathLib::KelvinVector::tensorToKelvin
KelvinVectorType< DisplacementDim > tensorToKelvin(Eigen::Matrix< double, 3, 3 > const &m)
ProcessLib::TH2M::ConstitutiveRelations::SwellingModel::eval
void eval(SpaceTimeData const &x_t, MediaData const &media_data, ElasticTangentStiffnessData< DisplacementDim > const &C_el_data, SaturationData const &S_L_data, PrevState< SaturationData > const &S_L_prev_data, PrevState< SwellingDataStateful< DisplacementDim > > const &prev_state, SwellingDataStateful< DisplacementDim > &state, SwellingDataStateless< DisplacementDim > &out) const
Definition TH2M/ConstitutiveRelations/Swelling.cpp:15

References ProcessLib::ConstitutiveRelations::SpaceTimeData::dt, ProcessLib::TH2M::ConstitutiveRelations::SwellingDataStateless< DisplacementDim >::eps_m, MaterialPropertyLib::formEigenTensor(), MaterialPropertyLib::VariableArray::liquid_saturation, ProcessLib::TH2M::ConstitutiveRelations::SaturationData::S_L, ProcessLib::TH2M::ConstitutiveRelations::SwellingDataStateful< DisplacementDim >::sigma_sw, ProcessLib::TH2M::ConstitutiveRelations::ElasticTangentStiffnessData< DisplacementDim >::stiffness_tensor, ProcessLib::TH2M::ConstitutiveRelations::MediaData::swelling_stress_rate_solid, ProcessLib::ConstitutiveRelations::SpaceTimeData::t, MathLib::KelvinVector::tensorToKelvin(), MaterialPropertyLib::Property::value(), and ProcessLib::ConstitutiveRelations::SpaceTimeData::x.

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