ConstitutiveSetting.h

Go to the documentation of this file.

 
#pragma once
 
#include "ConstitutiveData.h"
#include "ConstitutiveModels.h"
 
namespace ProcessLib::ThermoRichardsMechanics
{
namespace ConstitutiveStress_StrainTemperature
{
template <int DisplacementDim>
struct ConstitutiveSetting
{
    void init(ConstitutiveModels<DisplacementDim>& models, double const t,
              double const dt, ParameterLib::SpatialPosition const& x_position,
              MediaData const& media_data,
              TemperatureData<DisplacementDim> const& T_data,
              StatefulData<DisplacementDim>& state,
              StatefulDataPrev<DisplacementDim>& prev_state) 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;
 
    static KelvinVector<DisplacementDim> const& statefulStress(
        StatefulData<DisplacementDim> const& state)
    {
        return std::get<EffectiveStressData<DisplacementDim>>(state).sigma_eff;
    }
    static KelvinVector<DisplacementDim>& statefulStress(
        StatefulData<DisplacementDim>& state)
    {
        return std::get<EffectiveStressData<DisplacementDim>>(state).sigma_eff;
    }
 
    static void convertInitialStressType(
        StatefulData<DisplacementDim>& state,
        StatefulDataPrev<DisplacementDim>& prev_state,
        KelvinVector<DisplacementDim> const& pore_pressure_part)
    {
        auto& sigma_eff =
            std::get<EffectiveStressData<DisplacementDim>>(state).sigma_eff;
        sigma_eff.noalias() += pore_pressure_part;
 
        (std::get<PrevState<EffectiveStressData<DisplacementDim>>>(prev_state)
             ->sigma_eff)
            .noalias() = sigma_eff;
    }
};
 
extern template struct ConstitutiveSetting<2>;
extern template struct ConstitutiveSetting<3>;
}  // namespace ConstitutiveStress_StrainTemperature
}  // namespace ProcessLib::ThermoRichardsMechanics