ProcessLib::ThermoRichardsMechanics::ConstitutiveStress_StrainTemperature::SwellingModel< DisplacementDim > Struct Template Reference

Detailed Description

template<int DisplacementDim>
struct ProcessLib::ThermoRichardsMechanics::ConstitutiveStress_StrainTemperature::SwellingModel< DisplacementDim >

Definition at line 44 of file Swelling.h.

#include <Swelling.h>

Public Member Functions

void

eval (SpaceTimeData const &x_t, MediaData const &media_data, ElasticTangentStiffnessData< DisplacementDim > const &C_el_data, StrainData< DisplacementDim > const &eps_data, PrevState< StrainData< DisplacementDim > > const &eps_prev_data, SaturationData const &S_L_data, SaturationDataDeriv const &dS_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::ThermoRichardsMechanics::ConstitutiveStress_StrainTemperature::SwellingModel< DisplacementDim >::eval	(	SpaceTimeData const &	x_t,
		MediaData const &	media_data,
		ElasticTangentStiffnessData< DisplacementDim > const &	C_el_data,
		StrainData< DisplacementDim > const &	eps_data,
		PrevState< StrainData< DisplacementDim > > const &	eps_prev_data,
		SaturationData const &	S_L_data,
		SaturationDataDeriv const &	dS_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 22 of file Swelling.cpp.

{
    namespace MPL = MaterialPropertyLib;
 
    static constexpr int kelvin_vector_size =
        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
    using Invariants = MathLib::KelvinVector::Invariants<kelvin_vector_size>;
 
    auto const& solid_phase = media_data.solid;
 
    if (!solid_phase.hasProperty(MPL::PropertyType::swelling_stress_rate))
    {
        out.eps_m.setZero();
        out.J_up_BT_K_N.setZero();
        return;
    }
 
    MPL::VariableArray variables;
    variables.liquid_saturation = S_L_data.S_L;
 
    MPL::VariableArray variables_prev;
    variables_prev.liquid_saturation = S_L_prev_data->S_L;
 
    auto const& identity2 = MathLib::KelvinVector::Invariants<
        MathLib::KelvinVector::kelvin_vector_dimensions(
            DisplacementDim)>::identity2;
 
    // Swelling and possibly volumetric strain rate update.
 
    // If there is swelling, compute it. Update volumetric strain rate,
    // s.t. it corresponds to the mechanical part only.
    state.sigma_sw = prev_state->sigma_sw;
 
    auto const sigma_sw_dot =
        MathLib::KelvinVector::tensorToKelvin<DisplacementDim>(
            MPL::formEigenTensor<3>(
                solid_phase[MPL::PropertyType::swelling_stress_rate].value(
                    variables, variables_prev, x_t.x, x_t.t, x_t.dt)));
    state.sigma_sw += sigma_sw_dot * x_t.dt;
 
    auto const C_el_inv = C_el_data.C_el.inverse().eval();
 
    // !!! Misusing volumetric strain for mechanical volumetric
    // strain just to update the transport porosity !!!
    variables.volumetric_strain =
        Invariants::trace(eps_data.eps) +
        identity2.transpose() * C_el_inv * state.sigma_sw;
    variables_prev.volumetric_strain =
        Invariants::trace(eps_prev_data->eps) +
        identity2.transpose() * C_el_inv * prev_state->sigma_sw;
 
    out.eps_m.noalias() = C_el_inv * (state.sigma_sw - prev_state->sigma_sw);
 
    using DimMatrix = Eigen::Matrix<double, 3, 3>;
    auto const dsigma_sw_dS_L =
        MathLib::KelvinVector::tensorToKelvin<DisplacementDim>(
            solid_phase.property(MPL::PropertyType::swelling_stress_rate)
                .template dValue<DimMatrix>(variables, variables_prev,
                                            MPL::Variable::liquid_saturation,
                                            x_t.x, x_t.t, x_t.dt));
 
    out.J_up_BT_K_N = -dsigma_sw_dS_L * dS_L_data.dS_L_dp_cap;
}

References ProcessLib::ThermoRichardsMechanics::ConstitutiveStress_StrainTemperature::ElasticTangentStiffnessData< DisplacementDim >::C_el, ProcessLib::ThermoRichardsMechanics::SaturationDataDeriv::dS_L_dp_cap, ProcessLib::ConstitutiveRelations::SpaceTimeData::dt, ProcessLib::ConstitutiveRelations::StrainData< DisplacementDim >::eps, ProcessLib::ThermoRichardsMechanics::ConstitutiveStress_StrainTemperature::SwellingDataStateless< DisplacementDim >::eps_m, MaterialPropertyLib::formEigenTensor(), ProcessLib::ThermoRichardsMechanics::ConstitutiveStress_StrainTemperature::SwellingDataStateless< DisplacementDim >::J_up_BT_K_N, MathLib::KelvinVector::kelvin_vector_dimensions(), MaterialPropertyLib::liquid_saturation, MaterialPropertyLib::VariableArray::liquid_saturation, ProcessLib::ThermoRichardsMechanics::SaturationData::S_L, ProcessLib::ThermoRichardsMechanics::ConstitutiveStress_StrainTemperature::SwellingDataStateful< DisplacementDim >::sigma_sw, ProcessLib::ThermoRichardsMechanics::MediaData::solid, MaterialPropertyLib::swelling_stress_rate, ProcessLib::ConstitutiveRelations::SpaceTimeData::t, MathLib::KelvinVector::tensorToKelvin(), MaterialPropertyLib::VariableArray::volumetric_strain, and ProcessLib::ConstitutiveRelations::SpaceTimeData::x.

---

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

• ProcessLib/ThermoRichardsMechanics/ConstitutiveStress_StrainTemperature/Swelling.h
• ProcessLib/ThermoRichardsMechanics/ConstitutiveStress_StrainTemperature/Swelling.cpp