Homec4science
Diffusion µSpectre (18_OSX_fixes) Browse

muSpectre/src/materials/material_hyper_elastic1.hh
1f24dca2c71c18_OSX_fixes

Last Change
Raw File

material_hyper_elastic1.hh
View Options

/**
* file material_hyper_elastic1.hh
*
* @author Till Junge <till.junge@epfl.ch>
*
* @date 13 Nov 2017
*
* @brief Implementation for hyperelastic reference material like in de Geus
* 2017. This follows the simplest and likely not most efficient
* implementation (with exception of the Python law)
*
* @section LICENCE
*
* Copyright (C) 2017 Till Junge
*
* µSpectre is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License as
* published by the Free Software Foundation, either version 3, or (at
* your option) any later version.
*
* µSpectre is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with GNU Emacs; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* Boston, MA 02111-1307, USA.
*/
#include "materials/material_muSpectre_base.hh"
#ifndef MATERIAL_HYPER_ELASTIC1_H
#define MATERIAL_HYPER_ELASTIC1_H
namespace muSpectre {
template<Dim_t DimS, Dim_t DimM>
class MaterialHyperElastic1;
template <>
template <Dim_t DimS, Dim_t DimM>
struct MaterialMuSpectre_traits<MaterialHyperElastic1<DimS, DimM>>:
public MaterialMuSpectre_traits<void> {
using Parent = MaterialMuSpectre_traits<void>;
using InternalVariables = typename Parent::DefaultInternalVariables;
};
//! DimS spatial dimension (dimension of problem
//! DimM material_dimension (dimension of constitutive law)
template<Dim_t DimS, Dim_t DimM>
class MaterialHyperElastic1:
public MaterialMuSpectre<MaterialHyperElastic1<DimS, DimM>, DimS, DimM>
{
public:
using Parent = MaterialMuSpectre<MaterialHyperElastic1, DimS, DimM>;
using NeedTangent = typename Parent::NeedTangent;
using GFieldCollection_t = typename Parent::GFieldCollection_t;
// declare what type of strain measure your law takes as input
constexpr static auto strain_measure{StrainMeasure::GreenLagrange};
// declare what type of stress measure your law yields as output
constexpr static auto stress_measure{StressMeasure::PK2};
// declare whether the derivative of stress with respect to strain is uniform
constexpr static bool uniform_stiffness = true;
// declare the type in which you wish to receive your strain measure
using StrainMap_t = MatrixFieldMap<GFieldCollection_t, Real, DimM, DimM, true>;
using StressMap_t = MatrixFieldMap<GFieldCollection_t, Real, DimM, DimM>;
using TangentMap_t = T4MatrixFieldMap<GFieldCollection_t, Real, DimM>;
using Strain_t = typename StrainMap_t::const_reference;
using Stress_t = typename StressMap_t::reference;
using Tangent_t = typename TangentMap_t::reference;
using Stiffness_t = Eigen::TensorFixedSize
<Real, Eigen::Sizes<DimM, DimM, DimM, DimM>>;
using InternalVariables = typename Parent::DefaultInternalVariables;
//! Default constructor
MaterialHyperElastic1() = delete;
//! Copy constructor
MaterialHyperElastic1(const MaterialHyperElastic1 &other) = delete;
//! Construct by name, Young's modulus and Poisson's ratio
MaterialHyperElastic1(std::string name, Real young, Real poisson);
//! Move constructor
MaterialHyperElastic1(MaterialHyperElastic1 &&other) noexcept = delete;
//! Destructor
virtual ~MaterialHyperElastic1() noexcept = default;
//! Copy assignment operator
MaterialHyperElastic1& operator=(const MaterialHyperElastic1 &other) = delete;
//! Move assignment operator
MaterialHyperElastic1& operator=(MaterialHyperElastic1 &&other) noexcept = delete;
template <class s_t>
inline decltype(auto) evaluate_stress(s_t && E);
template <class s_t>
inline decltype(auto) evaluate_stress_tangent(s_t && E);
const Tangent_t & get_stiffness() const;
protected:
const Real young, poisson, lambda, mu;
const Stiffness_t C;
private:
};
/* ---------------------------------------------------------------------- */
template <Dim_t DimS, Dim_t DimM>
template <class s_t>
decltype(auto)
MaterialHyperElastic1<DimS, DimM>::evaluate_stress(s_t && E) {
return E.trace()*lambda * Strain_t::Identity() + 2*mu*E;
}
/* ---------------------------------------------------------------------- */
template <Dim_t DimS, Dim_t DimM>
template <class s_t>
decltype(auto)
MaterialHyperElastic1<DimS, DimM>::evaluate_stress_tangent(s_t && E) {
return std::make_tuple(std::move(this->evaluate_stress(std::move(E))),
std::move(Tangent_t(const_cast<double*>(this->C.data()))));
}
} // muSpectre
#endif /* MATERIAL_HYPER_ELASTIC1_H */
c4science · Help