HydroMechanicsProcessData.h

Go to the documentation of this file.

 
#pragma once
 
#include <Eigen/Core>
#include <memory>
#include <utility>
#include <variant>
 
#include "MaterialLib/MPL/MaterialSpatialDistributionMap.h"
#include "MaterialLib/MPL/VariableType.h"
#include "ParameterLib/Parameter.h"
#include "ProcessLib/Common/HydroMechanics/InitialStress.h"
 
namespace MaterialLib
{
namespace Solids
{
template <int DisplacementDim>
struct MechanicsBase;
}
}  // namespace MaterialLib
namespace ProcessLib
{
namespace HydroMechanics
{
 
struct Monolithic
{
};
 
struct Staggered
{
    double const fixed_stress_stabilization_parameter;
 
    bool const fixed_stress_over_time_step;
};
 
using CouplingScheme = std::variant<Monolithic, Staggered>;
 
template <int DisplacementDim>
struct HydroMechanicsProcessData
{
    MeshLib::PropertyVector<int> const* const material_ids = nullptr;
 
    MaterialPropertyLib::MaterialSpatialDistributionMap media_map;
 
    std::map<int, std::unique_ptr<
                      MaterialLib::Solids::MechanicsBase<DisplacementDim>>>
        solid_materials;
 
    InitialStress const initial_stress;
 
    Eigen::Matrix<double, DisplacementDim, 1> const specific_body_force;
 
    bool const mass_lumping;
 
    CouplingScheme coupling_scheme;
 
    int const hydraulic_process_id;
 
    int const mechanics_related_process_id;
 
    const bool use_taylor_hood_elements;
 
    MaterialPropertyLib::Variable phase_pressure;
 
    MeshLib::PropertyVector<double>* pressure_interpolated = nullptr;
    std::array<MeshLib::PropertyVector<double>*, 3> principal_stress_vector = {
        nullptr, nullptr, nullptr};
    MeshLib::PropertyVector<double>* principal_stress_values = nullptr;
 
    MeshLib::PropertyVector<double>* permeability = nullptr;
 
    bool isMonolithicSchemeUsed() const
    {
        return std::holds_alternative<Monolithic>(coupling_scheme);
    }
 
    EIGEN_MAKE_ALIGNED_OPERATOR_NEW;
};
 
}  // namespace HydroMechanics
}  // namespace ProcessLib