ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveSetting< DisplacementDim > Struct Template Reference

Detailed Description

template<int DisplacementDim>
struct ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveSetting< DisplacementDim >

Definition at line 21 of file ConstitutiveSetting.h.

#include <ConstitutiveSetting.h>

Public Member Functions
void	init (ConstitutiveModels< DisplacementDim > &, double const, double const, ParameterLib::SpatialPosition const &, MediaData const &, TemperatureData< DisplacementDim > const &, StatefulData< DisplacementDim > &, StatefulDataPrev< DisplacementDim > &) const
void	eval (ConstitutiveModels< DisplacementDim > &models, double const t, double const dt, ParameterLib::SpatialPosition const &x_position, MaterialPropertyLib::Medium const &medium, TemperatureData< DisplacementDim > const &T_data, CapillaryPressureData< DisplacementDim > const &p_cap_data, KelvinVector< DisplacementDim > const &eps_arg, StatefulData< DisplacementDim > &state, StatefulDataPrev< DisplacementDim > const &prev_state, MaterialStateData< DisplacementDim > &mat_state, ConstitutiveTempData< DisplacementDim > &tmp, OutputData< DisplacementDim > &out, ConstitutiveData< DisplacementDim > &cd) const
	Evaluate the constitutive setting.

Static Public Member Functions
static KelvinVector< DisplacementDim > const &	statefulStress (StatefulData< DisplacementDim > const &state)
static KelvinVector< DisplacementDim > &	statefulStress (StatefulData< DisplacementDim > &state)
static void	convertInitialStressType (StatefulData< DisplacementDim > &state, StatefulDataPrev< DisplacementDim > &prev_state, KelvinVector< DisplacementDim > const &pore_pressure_part)

Member Function Documentation

convertInitialStressType()

template<int DisplacementDim>

static void ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveSetting< DisplacementDim >::convertInitialStressType ( StatefulData< DisplacementDim > & state, StatefulDataPrev< DisplacementDim > & prev_state, KelvinVector< DisplacementDim > const & pore_pressure_part )

inlinestatic

In case that the input initial data for state.s_mech_data.sigma_total are effective stress values, state.s_mech_data.sigma_total is reset to total stress.

Definition at line 59 of file ConstitutiveSetting.h.

    {
        auto& sigma_total =
            std::get<TotalStressData<DisplacementDim>>(state).sigma_total;
        sigma_total.noalias() -= pore_pressure_part;
 
        (std::get<PrevState<TotalStressData<DisplacementDim>>>(prev_state)
             ->sigma_total)
            .noalias() = sigma_total;
    }

eval()

template<int DisplacementDim>

void ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveSetting< DisplacementDim >::eval	(	ConstitutiveModels< DisplacementDim > &	models,
		double const	t,
		double const	dt,
		ParameterLib::SpatialPosition const &	x_position,
		MaterialPropertyLib::Medium const &	medium,
		TemperatureData< DisplacementDim > const &	T_data,
		CapillaryPressureData< DisplacementDim > const &	p_cap_data,
		KelvinVector< DisplacementDim > const &	eps_arg,
		StatefulData< DisplacementDim > &	state,
		StatefulDataPrev< DisplacementDim > const &	prev_state,
		MaterialStateData< DisplacementDim > &	mat_state,
		ConstitutiveTempData< DisplacementDim > &	tmp,
		OutputData< DisplacementDim > &	out,
		ConstitutiveData< DisplacementDim > &	cd ) const

Evaluate the constitutive setting.

Definition at line 20 of file ConstitutiveSetting.cpp.

{
    namespace G = ProcessLib::Graph;
 
    auto const aux_data = std::tuple{SpaceTimeData{x_position, t, dt},
                                     MediaData{medium}, T_data, p_cap_data};
 
    auto const mat_state_tuple = std::tie(mat_state);
 
    // TODO will eps lag one iteration behind? (since it's not updated after
    // solving the global equation system)
    std::get<StrainData<DisplacementDim>>(state).eps.noalias() = eps_arg;
 
    G::eval(models.biot_model, aux_data, tmp);
 
    G::eval(models.s_mech_model, aux_data, tmp, state, prev_state,
            mat_state_tuple, cd);
 
    G::eval(models.solid_compressibility_model, aux_data, tmp, cd);
 
    G::eval(models.bishops_model, aux_data, state, tmp);
    // TODO why not ordinary state tracking?
    G::eval(models.bishops_prev_model, aux_data, prev_state, tmp);
    G::eval(models.poro_model, aux_data, tmp, state, prev_state);
 
    // TODO move to local assembler?
    {
        auto const& biot_data = std::get<BiotData>(tmp);
        auto const& poro_data = std::get<PorosityData>(state);
 
        if (biot_data() < poro_data.phi)
        {
            OGS_FATAL(
                "ThermoRichardsMechanics: Biot-coefficient {} is smaller than "
                "porosity {} in element/integration point {}/{}.",
                biot_data(), poro_data.phi, *x_position.getElementID(),
                *x_position.getIntegrationPoint());
        }
    }
 
    G::eval(models.rho_L_model, aux_data, out);
    G::eval(models.rho_S_model, aux_data, state, out);
    G::eval(models.grav_model, state, out, tmp, cd);
    G::eval(models.mu_L_model, aux_data, out);
    G::eval(models.transport_poro_model, aux_data, tmp, state, prev_state);
    G::eval(models.perm_model, aux_data, state, out, cd, tmp);
    G::eval(models.th_osmosis_model, aux_data, out, cd);
    G::eval(models.darcy_model, aux_data, out, tmp, cd);
    G::eval(models.heat_storage_and_flux_model, aux_data, out, state, tmp, cd);
    G::eval(models.vapor_diffusion_model, aux_data, out, state, tmp, cd);
 
    // TODO Not needed for solid mechanics (solid thermal expansion is computed
    // by the solid material model), but for fluid expansion. This duplication
    // should be avoided in the future.
    G::eval(models.s_therm_exp_model, aux_data, tmp);
 
    G::eval(models.f_therm_exp_model, aux_data, tmp, state, out);
    G::eval(models.storage_model, aux_data, tmp, state, out, prev_state, cd);
    G::eval(models.eq_p_model, aux_data, state, tmp, out, cd);
    G::eval(models.eq_T_model, cd);
}

References ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::biot_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::bishops_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::bishops_prev_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::darcy_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::eq_p_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::eq_T_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::f_therm_exp_model, ParameterLib::SpatialPosition::getElementID(), ParameterLib::SpatialPosition::getIntegrationPoint(), ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::grav_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::heat_storage_and_flux_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::mu_L_model, OGS_FATAL, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::perm_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::poro_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::rho_L_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::rho_S_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::s_mech_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::s_therm_exp_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::solid_compressibility_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::storage_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::th_osmosis_model, ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::transport_poro_model, and ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveModels< DisplacementDim >::vapor_diffusion_model.

init()

template<int DisplacementDim>

void ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveSetting< DisplacementDim >::init ( ConstitutiveModels< DisplacementDim > & , double const , double const , ParameterLib::SpatialPosition const & , MediaData const & , TemperatureData< DisplacementDim > const & , StatefulData< DisplacementDim > & , StatefulDataPrev< DisplacementDim > &  ) const

inline

Definition at line 23 of file ConstitutiveSetting.h.

    {
    }

statefulStress() [1/2]

template<int DisplacementDim>

static KelvinVector< DisplacementDim > & ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveSetting< DisplacementDim >::statefulStress ( StatefulData< DisplacementDim > & state )

inlinestatic

Definition at line 50 of file ConstitutiveSetting.h.

    {
        return std::get<TotalStressData<DisplacementDim>>(state).sigma_total;
    }

statefulStress() [2/2]

template<int DisplacementDim>

static KelvinVector< DisplacementDim > const & ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::ConstitutiveSetting< DisplacementDim >::statefulStress ( StatefulData< DisplacementDim > const & state )

inlinestatic

Definition at line 45 of file ConstitutiveSetting.h.

    {
        return std::get<TotalStressData<DisplacementDim>>(state).sigma_total;
    }

---

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

• ProcessLib/ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/ConstitutiveSetting.h
• ProcessLib/ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/ConstitutiveSetting.cpp