LargeDeformation/ConstitutiveRelations/Base.h

Go to the documentation of this file.

 
#pragma once
 
#include "BaseLib/StrongType.h"
#include "MaterialLib/MPL/Medium.h"
#include "MaterialLib/MPL/Utils/Tensor.h"
#include "MathLib/KelvinVector.h"
#include "ParameterLib/SpatialPosition.h"
#include "ProcessLib/ConstitutiveRelations/Base.h"
#include "ProcessLib/Reflection/ReflectionData.h"
 
namespace ProcessLib::LargeDeformation
{
 
using namespace ProcessLib::ConstitutiveRelations;
namespace KV = MathLib::KelvinVector;
 
template <int DisplacementDim>
using KelvinVector = KV::KelvinVectorType<DisplacementDim>;
 
template <int DisplacementDim>
using KelvinMatrix = KV::KelvinMatrixType<DisplacementDim>;
 
template <int DisplacementDim>
using GlobalDimVector = Eigen::Vector<double, DisplacementDim>;
 
template <int DisplacementDim>
using GlobalDimMatrix =
    Eigen::Matrix<double, DisplacementDim, DisplacementDim, Eigen::RowMajor>;

template <int D>
constexpr GlobalDimVector<D> DVnan()
{
    return GlobalDimVector<D>::Constant(nan);
}

template <int D>
constexpr GlobalDimMatrix<D> DMnan()
{
    return GlobalDimMatrix<D>::Constant(nan);
}
 
struct MediaData
{
    explicit MediaData(MaterialPropertyLib::Medium const& medium)
        : medium{medium}, solid{medium.phase("Solid")}
    {
    }
 
    MaterialPropertyLib::Medium const& medium;
    MaterialPropertyLib::Phase const& solid;
};
 
using Temperature = BaseLib::StrongType<double, struct TemperatureTag>;
 
template <int DisplacementDim>
struct DeformationGradientData
{
    // TODO Move initialization to the local assembler.
    MaterialPropertyLib::Tensor<DisplacementDim> deformation_gradient =
        MaterialPropertyLib::Tensor<DisplacementDim>::Zero();
    double volume_ratio = 0;
 
    static auto reflect()
    {
        using Self = DeformationGradientData<DisplacementDim>;
 
        return std::tuple{
            ProcessLib::Reflection::makeReflectionData(
                "deformation_gradient", &Self::deformation_gradient),
            ProcessLib::Reflection::makeReflectionData("volume_ratio",
                                                       &Self::volume_ratio)};
    }
};
}  // namespace ProcessLib::LargeDeformation