SolidCompressibilityModel.h

Go to the documentation of this file.

 
#pragma once
 
#include "ElasticTangentStiffnessData.h"
#include "ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/Biot.h"
#include "ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/SolidCompressibilityData.h"
 
namespace ProcessLib::ThermoRichardsMechanics
{
namespace ConstitutiveStress_StrainTemperature
{
// TODO remove typename SolidMaterial
template <int DisplacementDim, typename SolidMaterial>
struct SolidCompressibilityModel
{
    explicit SolidCompressibilityModel(SolidMaterial const& solid_material)
        : solid_material_(solid_material)
    {
    }
 
    void eval(SpaceTimeData const& x_t,
              BiotData const& biot_data,
              ElasticTangentStiffnessData<DisplacementDim> const& C_el_data,
              SolidCompressibilityData& out) const
    {
        out.beta_SR = (1 - biot_data()) / solid_material_.getBulkModulus(
                                              x_t.t, x_t.x, &C_el_data.C_el);
    }
 
    static SolidCompressibilityModel create(SolidMaterial const& solid_material)
    {
        return SolidCompressibilityModel{solid_material};
    }
 
private:
    SolidMaterial const& solid_material_;
};
}  // namespace ConstitutiveStress_StrainTemperature
}  // namespace ProcessLib::ThermoRichardsMechanics