Homec4science
Diffusion µSpectre (feat/repair_intersection_octree) Browse

muSpectre/src/materials/material_anisotropic.hh
80eb1f984de2feat/repair_intersectio…

Last Change
Raw File

material_anisotropic.hh
View Options

/**
* @file material_anisotropic.hh
*
* @author Ali Falsafi<ali.falsafi@epfl.ch>
*
* @date 9 Jul 2018
*
* @brief Base class for materials (constitutive models)
*
* Copyright © 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.
*/
#ifndef MATERIAL_ANISOTROPIC_H
#define MATERIAL_ANISOTROPIC_H
#include "materials/material_base.hh"
#include "materials/material_muSpectre_base.hh"
#include "common/common.hh"
#include "common/T4_map_proxy.hh"
namespace muSpectre {
template <Dim_t DimS, Dim_t DimM>
class MaterialAnisotropic;
//traits for anisotropic material
template <Dim_t DimS, Dim_t DimM>
struct MaterialMuSpectre_traits<MaterialAnisotropic<DimS, DimM>> {
using Parent = MaterialMuSpectre_traits<void>;//!< base for elasticity
//! global field collection
using GFieldCollection_t = typename
MaterialBase<DimS, DimM>::GFieldCollection_t;
//! expected map type for strain fields
using StrainMap_t = MatrixFieldMap<GFieldCollection_t, Real, DimM, DimM, true>;
//! expected map type for stress fields
using StressMap_t = MatrixFieldMap<GFieldCollection_t, Real, DimM, DimM>;
//! expected map type for tangent stiffness fields
using TangentMap_t = T4MatrixFieldMap<GFieldCollection_t, Real, DimM>;
//! 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};
//! anisotropicity without internal variables
using InternalVariables = std::tuple<>;
};
/**
* Material implementation for anisotropic constitutive law
*/
template <Dim_t DimS, Dim_t DimM=DimS>
class MaterialAnisotropic:
public MaterialMuSpectre< MaterialAnisotropic<DimS, DimM>, DimS, DimM>
{
//! base class
using Parent = MaterialMuSpectre<MaterialAnisotropic, DimS, DimM>;
using Stiffness_t =T4Mat<Real, DimM>;
//! traits of this material
using traits = MaterialMuSpectre_traits<MaterialAnisotropic>;
//! this law does not have any internal variables
using InternalVariables = typename traits::InternalVariables;
public:
using parent = MaterialMuSpectre
<MaterialAnisotropic<DimS, DimM>, DimS, DimM>;
//! global field collection
using GFieldCollection_t = typename
MaterialBase<DimS, DimM>::GFieldCollection_t;
//! expected map type for tangent stiffness fields
using Tangent_t = T4MatrixFieldMap<GFieldCollection_t, Real, DimM>;
//! Default constructor
MaterialAnisotropic() = delete;
// constructor
// a std::vector is utilized as the input of the constructor to
//enable us to check its length so to prevent user mistake
MaterialAnisotropic(std::string name, std::vector<Real> input_c);
//! Copy constructor
MaterialAnisotropic(const MaterialAnisotropic &other) = delete;
//! Move constructor
MaterialAnisotropic(MaterialAnisotropic &&other) = delete;
//! Destructor
virtual ~MaterialAnisotropic() = default;
template <class s_t>
inline decltype(auto) evaluate_stress(s_t && E);
/**
* evaluates both second Piola-Kirchhoff stress and stiffness given
* the Green-Lagrange strain (or Cauchy stress and stiffness if
* called with a small strain tensor) and the local stiffness tensor.
*/
template <class s_t>
inline decltype(auto)
evaluate_stress_tangent(s_t && E);
/**
* return the empty internals tuple
*/
InternalVariables & get_internals() {
return this->internal_variables;};
// takes the elements of the C and makes it:
virtual auto c_maker(std::vector<Real> input) -> Stiffness_t;
protected:
Stiffness_t C; //!< stiffness tensor
//! empty tuple
InternalVariables internal_variables{};
};
/* ---------------------------------------------------------------------- */
template <Dim_t DimS, Dim_t DimM>
template <class s_t>
decltype(auto)
MaterialAnisotropic<DimS, DimM>::
evaluate_stress(s_t && E) {
return Matrices::tensmult(this->C, E);
}
/* ---------------------------------------------------------------------- */
template <Dim_t DimS, Dim_t DimM>
template <class s_t>
decltype(auto)
MaterialAnisotropic<DimS, DimM>::
evaluate_stress_tangent(s_t && E) {
return std::make_tuple
(evaluate_stress(E), this->C);
}
}//muSpectre
#endif
c4science · Help