gsMaterialMatrixUtils.h

Go to the documentation of this file.

#pragma once

#include <gsKLShell/src/gsMaterialMatrixBase.h>
#include <gsIO/gsOptionList.h>

namespace gismo
{


enum class MatIntegration : short_t
{
    NotIntegrated   = 0,
    Integrated      = 1,
    Constant        = 2,
    Linear          = 3,
};


enum class Material : short_t
{
    SvK = 0,
    NH  = 1,
    NH_ext  = 2,
    MR  = 3,
    OG  = 4
};


enum class Implementation : short_t
{
    Composite   = 0,
    Analytical  = 1,
    Generalized = 2,
    Spectral    = 3
};


enum class MaterialOutput : short_t
{
    Generic = 0,
    Density = 1,
    VectorN = 2,
    VectorM = 3,
    CauchyVectorN = 4,
    CauchyVectorM = 5,
    MatrixA  = 6,
    MatrixB  = 7,
    MatrixC  = 8,
    MatrixD  = 9,
    PStress  = 10,
    PStressDir= 11,
    PStressN = 12,
    PStressM = 13,
    PCauchyStressN = 14,
    PCauchyStressM = 15,
    PStrainN = 16,
    PStrainM = 17,
    Stretch = 18,       // ONLY ON MID-PLANE
    StretchDir = 19,    // ONLY ON MID-PLANE
    TensionField = 20,  // Tension field indicator (1: slack, 0: wrinkled, -1: taut)
    Theta = 21,  // Tension field indicator (1: slack, 0: wrinkled, -1: taut)
    Gamma = 22,  // Tension field indicator (1: slack, 0: wrinkled, -1: taut)
    Strain  = 23,
    StrainN = 24,
    StrainM = 25,
    Stress  = 26,
    StressN = 27,
    StressM = 28,
    CauchyStress  = 29,
    CauchyStressN = 30,
    CauchyStressM = 31,
    Spec2CovTransform = 101,  // Transformation matrix from spectral to covariant basis
    Spec2ConTransform = 102,  // Transformation matrix from spectral to contravariant basis
    Cov2CartTransform = 103,  // Transformation matrix from covariant basis to cartesian basis
    Con2CartTransform = 104,  // Transformation matrix from contravariant basis to cartesian basis
    StretchTransform  = 105,  // Transformation matrix from principal stretch axes to contravariant axes, such that E_con = E_p * T
    PStressTransform  = 106,  // Transformation matrix from principal stress axes to covariant axes, such that S_cov = S_p * T
    Thickness = 1000,
    Parameters = 1001,
    Deformation = 1002,
    // FINISH VON MISES WITH THE CAUCHY STRESS!!!
};


template<enum Material material, enum Implementation implementation>
struct encodeMat_id {
  static const constexpr short_t id = 10*(short_t)implementation + (short_t)material;
};

template<short_t id>
struct decodeMat_id {
  static const constexpr enum Material          material = (enum Material)(id%10);
  static const constexpr enum Implementation    implementation = (enum Implementation)((id/10)%10);
};

} // namespace