solidmechanics/hyperelastic/model.hh

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
//
// SPDX-FileCopyrightText: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
// SPDX-License-Identifier: GPL-3.0-or-later
//
#ifndef DUMUX_SOLIDMECHANICS_HYPERELASTIC_MODEL_HH
#define DUMUX_SOLIDMECHANICS_HYPERELASTIC_MODEL_HH
 
#include <type_traits>
 
#include <dumux/common/properties.hh>
#include <dumux/common/properties/model.hh>
#include <dumux/discretization/concepts.hh>
#include <dumux/discretization/gridvariables.hh>
#include <dumux/discretization/cvfe/gridvariablescache_.hh>
#include <dumux/discretization/cvfe/hybrid/gridvariablescache.hh>
#include <dumux/discretization/cvfe/interpolationpointdata.hh>
#include <dumux/discretization/cvfe/variablesadapter.hh>
#include <dumux/discretization/fem/gridvariablescache.hh>
 
#include "spatialparams.hh"
#include "localresidual.hh"
#include "variables.hh"
 
namespace Dumux {
 
template<class PV, class MT>
struct HyperelasticVolumeVariablesTraits
{
    using PrimaryVariables = PV;
    using ModelTraits = MT;
};
 
struct HyperelasticIndices
{
    static constexpr int displacementIdx(int i) { return i; };
    static constexpr int equationIdx(int i) { return i; };
};

template<int dim>
struct HyperelasticModelTraits
{
    using Indices = HyperelasticIndices;
    static constexpr int numEq() { return dim; }
};
 
namespace Properties {
 
// Create new type tags
namespace TTag {
struct Hyperelastic { using InheritsFrom = std::tuple<ModelProperties>; };
} // end namespace TTag
 
template<class TypeTag>
struct ModelTraits<TypeTag, TTag::Hyperelastic>
{ using type = HyperelasticModelTraits<GetPropType<TypeTag, Properties::GridGeometry>::GridView::dimension>; };
 
template<class TypeTag>
struct LocalResidual<TypeTag, TTag::Hyperelastic>
{ using type = HyperelasticLocalResidual<TypeTag>; };
 
template<class TypeTag>
struct VolumeVariables<TypeTag, TTag::Hyperelastic>
{
    using PV = GetPropType<TypeTag, Properties::PrimaryVariables>;
    using MT = GetPropType<TypeTag, Properties::ModelTraits>;
    using Traits = HyperelasticVolumeVariablesTraits<PV, MT>;
    using type = HyperelasticVariables<Traits>;
};

template<class TypeTag>
struct GridVariables<TypeTag, TTag::Hyperelastic>
{
private:
    using GG = GetPropType<TypeTag, Properties::GridGeometry>;
    using ElementDiscretization = typename GG::LocalView;
    using Problem = GetPropType<TypeTag, Properties::Problem>;
    using Variables = Dumux::Detail::CVFE::VariablesAdapter<GetPropType<TypeTag, Properties::VolumeVariables>>;
    using IPDataCache = Dumux::CVFE::LocalBasisInterpolationPointData<GG>;
    static constexpr bool enableCache = getPropValue<TypeTag, Properties::EnableGridVolumeVariablesCache>();
 
    using FECache = Dumux::Experimental::FEGridVariablesCache<
        Dumux::Experimental::FEDefaultGridVariablesCacheTraits<Problem, Variables, IPDataCache>, enableCache>;
    using HybridCache = Dumux::Experimental::CVFE::HybridCVFEGridVariablesCache<
        Dumux::Experimental::CVFE::HybridCVFEDefaultGridVariablesCacheTraits<Problem, Variables, IPDataCache>, enableCache>;
    using CVFECache = Dumux::Experimental::CVFE::CVFEGridVariablesCache<
        Dumux::Experimental::CVFE::CVFEDefaultGridVariablesCacheTraits<Problem, Variables, IPDataCache>, enableCache>;
 
    using Cache = std::conditional_t<Dumux::Experimental::Concepts::FEElementDiscretization<ElementDiscretization>, FECache,
                  std::conditional_t<Dumux::Experimental::Concepts::HybridElementDiscretization<ElementDiscretization>, HybridCache,
                                     CVFECache>>;
public:
    using type = Dumux::Experimental::GridVariables<GG, Cache>;
};
 
// Set the default spatial parameters
template<class TypeTag>
struct SpatialParams<TypeTag, TTag::Hyperelastic>
{
    using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>;
    using Scalar = GetPropType<TypeTag, Properties::Scalar>;
    using type = DefaultHyperelasticSpatialParams<GridGeometry, Scalar>;
};
 
} // end namespace Properties
} // end namespace Dumux
 
#endif