Variable.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <MGIS/Behaviour/Variable.hxx>
 
#include "BaseLib/cpp23.h"
#include "MaterialLib/MPL/Utils/Tensor.h"
#include "MaterialLib/MPL/VariableType.h"
#include "MathLib/KelvinVector.h"
 
namespace MaterialLib::Solids::MFront
{
template <typename Derived>
struct Variable
{
    template <int DisplacementDim>
    static constexpr std::size_t size()
    {
        return rows<DisplacementDim>() * cols<DisplacementDim>();
    }

    template <int DisplacementDim>
    static constexpr std::size_t rows()
    {
        using T = mgis::behaviour::Variable::Type;
 
        switch (Derived::type)
        {
            case T::SCALAR:
                return 1;
            case T::VECTOR:
                return DisplacementDim;
            case T::STENSOR:
                return MathLib::KelvinVector::kelvin_vector_dimensions(
                    DisplacementDim);
            case T::TENSOR:
                return MaterialPropertyLib::tensorSize(DisplacementDim);
            case T::VECTOR_1D:
            case T::VECTOR_2D:
            case T::VECTOR_3D:
            case T::STENSOR_1D:
            case T::STENSOR_2D:
            case T::STENSOR_3D:
            case T::TENSOR_1D:
            case T::TENSOR_2D:
            case T::TENSOR_3D:
            case T::HIGHER_ORDER_TENSOR:
            case T::ARRAY:
                OGS_FATAL("Unsupported MGIS variable type {}.",
                          BaseLib::to_underlying(Derived::type));
        }
    }

    template <int DisplacementDim>
    static constexpr std::size_t cols()
    {
        using T = mgis::behaviour::Variable::Type;
        switch (Derived::type)
        {
            case T::SCALAR:
                return 1;
            case T::VECTOR:
                return 1;
            case T::STENSOR:
                return 1;
            case T::TENSOR:
                return 1;
            case T::VECTOR_1D:
            case T::VECTOR_2D:
            case T::VECTOR_3D:
            case T::STENSOR_1D:
            case T::STENSOR_2D:
            case T::STENSOR_3D:
            case T::TENSOR_1D:
            case T::TENSOR_2D:
            case T::TENSOR_3D:
            case T::HIGHER_ORDER_TENSOR:
            case T::ARRAY:
                OGS_FATAL("Unsupported MGIS variable type {}.",
                          BaseLib::to_underlying(Derived::type));
        }
    }
};

struct Strain : Variable<Strain>
{
    constexpr static const char* name = "Strain";

    constexpr static mgis::behaviour::Variable::Type type =
        mgis::behaviour::Variable::Type::STENSOR;

    constexpr static auto mpl_var =
        &MaterialPropertyLib::VariableArray::mechanical_strain;
};

static constexpr Strain strain;

struct GreenLagrangeStrain : Variable<GreenLagrangeStrain>
{
    constexpr static const char* name = "GreenLagrangeStrain";

    constexpr static mgis::behaviour::Variable::Type type =
        mgis::behaviour::Variable::Type::STENSOR;

    constexpr static auto mpl_var =
        &MaterialPropertyLib::VariableArray::mechanical_strain;
};

static constexpr GreenLagrangeStrain green_lagrange_strain;

struct DeformationGradient : Variable<DeformationGradient>
{
    constexpr static const char* name = "DeformationGradient";

    constexpr static mgis::behaviour::Variable::Type type =
        mgis::behaviour::Variable::Type::TENSOR;

    constexpr static auto mpl_var =
        &MaterialPropertyLib::VariableArray::deformation_gradient;
};

static constexpr DeformationGradient deformation_gradient;
 
struct LiquidPressure : Variable<LiquidPressure>
{
    constexpr static const char* name = "LiquidPressure";
 
    constexpr static mgis::behaviour::Variable::Type type =
        mgis::behaviour::Variable::Type::SCALAR;
 
    constexpr static auto mpl_var =
        &MaterialPropertyLib::VariableArray::liquid_phase_pressure;
};
 
static constexpr LiquidPressure liquid_pressure;
 
struct Stress : Variable<Stress>
{
    constexpr static const char* name = "Stress";
 
    constexpr static mgis::behaviour::Variable::Type type =
        mgis::behaviour::Variable::Type::STENSOR;
 
    constexpr static auto mpl_var = &MaterialPropertyLib::VariableArray::stress;
};
 
static constexpr Stress stress;
 
struct SecondPiolaKirchhoffStress : Variable<SecondPiolaKirchhoffStress>
{
    constexpr static const char* name = "SecondPiolaKirchhoffStress";
 
    constexpr static mgis::behaviour::Variable::Type type =
        mgis::behaviour::Variable::Type::STENSOR;
 
    constexpr static auto mpl_var = &MaterialPropertyLib::VariableArray::stress;
};
 
static constexpr SecondPiolaKirchhoffStress second_piola_kirchhoff_stress;
 
struct Saturation : Variable<Saturation>
{
    constexpr static const char* name = "Saturation";
 
    constexpr static mgis::behaviour::Variable::Type type =
        mgis::behaviour::Variable::Type::SCALAR;
 
    constexpr static auto mpl_var =
        &MaterialPropertyLib::VariableArray::liquid_saturation;
};
 
static constexpr Saturation saturation;
 
struct Temperature : Variable<Temperature>
{
    constexpr static const char* name = "Temperature";
 
    constexpr static mgis::behaviour::Variable::Type type =
        mgis::behaviour::Variable::Type::SCALAR;
 
    constexpr static auto mpl_var =
        &MaterialPropertyLib::VariableArray::temperature;
};
 
static constexpr Temperature temperature;
}  // namespace MaterialLib::Solids::MFront