Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <cstddef>
       7              : #include <memory>
       8              : #include <pup.h>
       9              : 
      10              : #include "DataStructures/Tensor/TypeAliases.hpp"
      11              : #include "Utilities/Gsl.hpp"
      12              : #include "Utilities/Serialization/CharmPupable.hpp"
      13              : #include "Utilities/TMPL.hpp"
      14              : 
      15              : /// \cond
      16              : class DataVector;
      17              : /// \endcond
      18              : 
      19              : namespace Elasticity {
      20              : /*!
      21              :  * \brief Constitutive (stress-strain) relations that characterize the elastic
      22              :  * properties of a material
      23              :  */
      24            1 : namespace ConstitutiveRelations {
      25              : 
      26              : /*!
      27              :  * \brief Base class for constitutive (stress-strain) relations that
      28              :  * characterize the elastic properties of a material
      29              :  *
      30              :  * \details A constitutive relation, in the context of elasticity, relates the
      31              :  * Stress \f$T^{ij}\f$ and Strain \f$S_{ij}=\nabla_{(i}u_{j)}\f$ within an
      32              :  * elastic material (see \ref Elasticity). For small stresses it is approximated
      33              :  * by the linear relation
      34              :  *
      35              :  * \f[
      36              :  * T^{ij} = -Y^{ijkl}S_{kl}
      37              :  * \f]
      38              :  *
      39              :  * (Eq. 11.17 in \cite ThorneBlandford2017) that is referred to as _Hooke's
      40              :  * law_. The constitutive relation in this linear approximation is determined by
      41              :  * the elasticity (or _Young's_) tensor \f$Y^{ijkl}=Y^{(ij)(kl)}=Y^{klij}\f$
      42              :  * that generalizes a simple proportionality to a three-dimensional and
      43              :  * (possibly) anisotropic material.
      44              :  *
      45              :  * \note We assume a Euclidean metric in Cartesian coordinates here (for now).
      46              :  */
      47              : template <size_t Dim>
      48            1 : class ConstitutiveRelation : public PUP::able {
      49              :  public:
      50            0 :   static constexpr size_t volume_dim = Dim;
      51              : 
      52            0 :   ConstitutiveRelation() = default;
      53            0 :   ConstitutiveRelation(const ConstitutiveRelation&) = default;
      54            0 :   ConstitutiveRelation& operator=(const ConstitutiveRelation&) = default;
      55            0 :   ConstitutiveRelation(ConstitutiveRelation&&) = default;
      56            0 :   ConstitutiveRelation& operator=(ConstitutiveRelation&&) = default;
      57            0 :   ~ConstitutiveRelation() override = default;
      58              : 
      59            0 :   WRAPPED_PUPable_abstract(ConstitutiveRelation);  // NOLINT
      60              : 
      61              :   /// Returns a `std::unique_ptr` pointing to a copy of the
      62              :   /// `ConstitutiveRelation`.
      63            1 :   virtual std::unique_ptr<ConstitutiveRelation<Dim>> get_clone() const = 0;
      64              : 
      65              :   /// The constitutive relation that characterizes the elastic properties of a
      66              :   /// material
      67            1 :   virtual void stress(gsl::not_null<tnsr::II<DataVector, Dim>*> stress,
      68              :                       const tnsr::ii<DataVector, Dim>& strain,
      69              :                       const tnsr::I<DataVector, Dim>& x) const = 0;
      70              : };
      71              : 
      72              : }  // namespace ConstitutiveRelations
      73              : }  // namespace Elasticity