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

Detailed Description

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

Definition at line 22 of file ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/SolidMechanics.h.

#include <SolidMechanics.h>

Public Member Functions
	SolidMechanicsModel (SolidConstitutiveRelation< DisplacementDim > const &solid_material)
void	eval (const SpaceTimeData &x_t, TemperatureData< DisplacementDim > const &T_data, CapillaryPressureData< DisplacementDim > const &p_cap_data, StrainData< DisplacementDim > const &eps_data, PrevState< StrainData< DisplacementDim > > const &eps_prev_data, MaterialStateData< DisplacementDim > &mat_state, PrevState< TotalStressData< DisplacementDim > > const &total_stress_data_prev, TotalStressData< DisplacementDim > &total_stress_data, EquivalentPlasticStrainData &equiv_plast_strain_data, SolidMechanicsDataStateless< DisplacementDim > &current_stateless, PrevState< SaturationData > const &S_L_prev_data, SaturationData &S_L_data, SaturationDataDeriv &dS_L_data) const

Static Public Member Functions
static SolidMechanicsModel	create (SolidConstitutiveRelation< DisplacementDim > const &solid_material)

Private Attributes
SolidConstitutiveRelation< DisplacementDim > const &	solid_material_
MSM::OGSMFrontTangentOperatorBlocksView< DisplacementDim, MSM::ForcesGradsCombinations< boost::mp11::mp_list< MSM::Strain, MSM::LiquidPressure >, boost::mp11::mp_list< MSM::Stress, MSM::Saturation >, boost::mp11::mp_list< MSM::Temperature > >::type >	tangent_operator_blocks_view_

Constructor & Destructor Documentation

SolidMechanicsModel()

template<int DisplacementDim>

ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::SolidMechanicsModel< DisplacementDim >::SolidMechanicsModel ( SolidConstitutiveRelation< DisplacementDim > const & solid_material )

inlineexplicit

Definition at line 24 of file ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/SolidMechanics.h.

        : solid_material_(solid_material),
          tangent_operator_blocks_view_{
              solid_material.createTangentOperatorBlocksView()}
    {
    }

References solid_material_, and tangent_operator_blocks_view_.

Referenced by create().

Member Function Documentation

create()

template<int DisplacementDim>

SolidMechanicsModel ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::SolidMechanicsModel< DisplacementDim >::create ( SolidConstitutiveRelation< DisplacementDim > const & solid_material )

inlinestatic

Definition at line 46 of file ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/SolidMechanics.h.

    {
        return SolidMechanicsModel{solid_material};
    }

References SolidMechanicsModel().

eval()

template<int DisplacementDim>

void ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::SolidMechanicsModel< DisplacementDim >::eval	(	const SpaceTimeData &	x_t,
		TemperatureData< DisplacementDim > const &	T_data,
		CapillaryPressureData< DisplacementDim > const &	p_cap_data,
		StrainData< DisplacementDim > const &	eps_data,
		PrevState< StrainData< DisplacementDim > > const &	eps_prev_data,
		MaterialStateData< DisplacementDim > &	mat_state,
		PrevState< TotalStressData< DisplacementDim > > const &	total_stress_data_prev,
		TotalStressData< DisplacementDim > &	total_stress_data,
		EquivalentPlasticStrainData &	equiv_plast_strain_data,
		SolidMechanicsDataStateless< DisplacementDim > &	current_stateless,
		PrevState< SaturationData > const &	S_L_prev_data,
		SaturationData &	S_L_data,
		SaturationDataDeriv &	dS_L_data ) const

Definition at line 11 of file ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/SolidMechanics.cpp.

{
    namespace MPL = MaterialPropertyLib;
 
    double const T_prev = T_data.T_prev;
    auto const& eps_total = eps_data.eps;
    auto const& eps_total_prev = eps_prev_data->eps;
    auto const& sigma_total_prev = total_stress_data_prev->sigma_total;
 
    // current state
    MPL::VariableArray variables;
    {
        // gradients
        // TODO currently we always pass strain via mechanical_strain
        variables.mechanical_strain = eps_total;
        variables.liquid_phase_pressure = -p_cap_data.p_cap;
 
        // external state variables
        variables.temperature = T_data.T;
    }
 
    // previous state
    MPL::VariableArray variables_prev;
    {
        // thermodynamic forces
        variables_prev.stress = sigma_total_prev;
        variables_prev.liquid_saturation = S_L_prev_data->S_L;
 
        // gradients
        // TODO currently we always pass strain via mechanical_strain
        variables_prev.mechanical_strain = eps_total_prev;
        variables_prev.liquid_phase_pressure = -p_cap_data.p_cap_prev;
 
        // external state variables
        variables_prev.temperature = T_prev;
    }
 
    auto solution = solid_material_.integrateStress(
        variables_prev, variables, x_t.t, x_t.x, x_t.dt,
        *mat_state.material_state_variables);
 
    if (!solution)
    {
        OGS_FATAL("Computation of local constitutive relation failed.");
    }
 
    auto& tdyn_forces_data = std::get<0>(*solution);
 
    auto const view = solid_material_.createThermodynamicForcesView();
 
    total_stress_data.sigma_total = view.block(MSM::stress, tdyn_forces_data);
    S_L_data.S_L = view.block(MSM::saturation, tdyn_forces_data);
    mat_state.material_state_variables = std::move(std::get<1>(*solution));
 
    auto const& tangent_operator_data = std::get<2>(*solution);
 
    current_stateless.stiffness_tensor = tangent_operator_blocks_view_.block(
        MSM::stress, MSM::strain, tangent_operator_data);
 
    dS_L_data.dS_L_dp_cap = -tangent_operator_blocks_view_.block(
        MSM::saturation, MSM::liquid_pressure, tangent_operator_data);
 
    current_stateless.J_uT_BT_K_N = tangent_operator_blocks_view_.block(
        MSM::stress, MSM::temperature, tangent_operator_data);
 
    current_stateless.J_up_BT_K_N = tangent_operator_blocks_view_.block(
        MSM::stress, MSM::liquid_pressure, tangent_operator_data);
 
    equiv_plast_strain_data.equivalent_plastic_strain =
        mat_state.material_state_variables->getEquivalentPlasticStrain();
}

References ProcessLib::ThermoRichardsMechanics::SaturationDataDeriv::dS_L_dp_cap, ProcessLib::ConstitutiveRelations::SpaceTimeData::dt, ProcessLib::ConstitutiveRelations::StrainData< DisplacementDim >::eps, ProcessLib::ThermoRichardsMechanics::EquivalentPlasticStrainData::equivalent_plastic_strain, ProcessLib::ThermoRichardsMechanics::SolidMechanicsDataStateless< DisplacementDim >::J_up_BT_K_N, ProcessLib::ThermoRichardsMechanics::SolidMechanicsDataStateless< DisplacementDim >::J_uT_BT_K_N, MaterialPropertyLib::VariableArray::liquid_phase_pressure, MaterialPropertyLib::VariableArray::liquid_saturation, ProcessLib::ThermoRichardsMechanics::MaterialStateData< DisplacementDim >::material_state_variables, MaterialPropertyLib::VariableArray::mechanical_strain, OGS_FATAL, ProcessLib::ThermoRichardsMechanics::CapillaryPressureData< DisplacementDim >::p_cap, ProcessLib::ThermoRichardsMechanics::CapillaryPressureData< DisplacementDim >::p_cap_prev, ProcessLib::ThermoRichardsMechanics::SaturationData::S_L, ProcessLib::ThermoRichardsMechanics::TotalStressData< DisplacementDim >::sigma_total, solid_material_, ProcessLib::ThermoRichardsMechanics::SolidMechanicsDataStateless< DisplacementDim >::stiffness_tensor, MaterialPropertyLib::VariableArray::stress, ProcessLib::ThermoRichardsMechanics::TemperatureData< DisplacementDim >::T, ProcessLib::ConstitutiveRelations::SpaceTimeData::t, ProcessLib::ThermoRichardsMechanics::TemperatureData< DisplacementDim >::T_prev, tangent_operator_blocks_view_, MaterialPropertyLib::VariableArray::temperature, and ProcessLib::ConstitutiveRelations::SpaceTimeData::x.

Member Data Documentation

solid_material_

template<int DisplacementDim>

SolidConstitutiveRelation<DisplacementDim> const& ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::SolidMechanicsModel< DisplacementDim >::solid_material_

private

Definition at line 53 of file ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/SolidMechanics.h.

Referenced by SolidMechanicsModel(), and eval().

tangent_operator_blocks_view_

template<int DisplacementDim>

MSM::OGSMFrontTangentOperatorBlocksView< DisplacementDim, MSM::ForcesGradsCombinations< boost::mp11::mp_list<MSM::Strain, MSM::LiquidPressure>, boost::mp11::mp_list<MSM::Stress, MSM::Saturation>, boost::mp11::mp_list<MSM::Temperature> >::type> ProcessLib::ThermoRichardsMechanics::ConstitutiveStressSaturation_StrainPressureTemperature::SolidMechanicsModel< DisplacementDim >::tangent_operator_blocks_view_

private

Definition at line 61 of file ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/SolidMechanics.h.

Referenced by SolidMechanicsModel(), and eval().

---

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

• ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/SolidMechanics.h
• ThermoRichardsMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/SolidMechanics.cpp