elasticity.h

Go to the documentation of this file.

// ------------------------------------------------------
//
// SPDX-License-Identifier: GPL-2.0-or-later
//
// SPDX-FileCopyrightText: Copyright (C) 2019-2026 by the ASPECT authors.
//
// This file is part of ASPECT.
//
// Detailed license information governing the source code and
// contributions can be found in the folder LICENSES and in
// LICENSE and CONTRIBUTING.md at the top level directory.
//
// ------------------------------------------------------

#ifndef _aspect_material_model_rheology_elasticity_h
#define _aspect_material_model_rheology_elasticity_h

#include <aspect/global.h>
#include <aspect/material_model/interface.h>
#include <aspect/simulator_access.h>

#include <deal.II/matrix_free/fe_point_evaluation.h>

namespace aspect
{
  namespace MaterialModel
  {
    template <int dim>
    class ElasticAdditionalOutputs : public NamedAdditionalMaterialOutputs<dim>
    {
      public:
        explicit ElasticAdditionalOutputs(const unsigned int n_points);

        std::vector<double> get_nth_output(const unsigned int idx) const override;

        std::vector<double> elastic_shear_moduli;

        std::vector<double> elastic_viscosity;

        std::vector<SymmetricTensor<2,dim>> deviatoric_stress;
    };



    namespace Rheology
    {
      template <int dim>
      class Elasticity : public ::aspect::SimulatorAccess<dim>
      {
        public:
          static
          void
          declare_parameters (ParameterHandler &prm);

          void
          parse_parameters (ParameterHandler &prm);

          void
          create_elastic_additional_outputs (MaterialModel::MaterialModelOutputs<dim> &out) const;

          void
          fill_elastic_outputs (const MaterialModel::MaterialModelInputs<dim> &in,
                                const std::vector<double> &average_elastic_shear_moduli,
                                MaterialModel::MaterialModelOutputs<dim> &out) const;

          void
          fill_elastic_additional_outputs (const MaterialModel::MaterialModelInputs<dim> &in,
                                           const std::vector<double> &average_elastic_shear_moduli,
                                           MaterialModel::MaterialModelOutputs<dim> &out) const;

          void
          fill_reaction_outputs (const MaterialModel::MaterialModelInputs<dim> &in,
                                 const std::vector<double> &average_elastic_shear_moduli,
                                 MaterialModel::MaterialModelOutputs<dim> &out) const;

          void
          fill_reaction_rates (const MaterialModel::MaterialModelInputs<dim> &in,
                               const std::vector<double> &average_elastic_shear_moduli,
                               MaterialModel::MaterialModelOutputs<dim> &out) const;

          const std::vector<double> &
          get_elastic_shear_moduli () const;

          double
          calculate_elastic_viscosity (const double shear_modulus) const;

          double
          calculate_viscoelastic_viscosity (const double viscosity,
                                            const double shear_modulus) const;

          SymmetricTensor<2,dim>
          calculate_viscoelastic_strain_rate (const SymmetricTensor<2,dim> &strain_rate,
                                              const SymmetricTensor<2, dim> &stress_0_advected,
                                              const SymmetricTensor<2, dim> &stress_old,
                                              const double viscosity_pre_yield,
                                              const double shear_modulus) const;

          bool
          require_unrotated_viscoelastic_stress() const;

          double
          elastic_timestep () const;

          double
          calculate_timestep_ratio() const;

        private:
          std::vector<SymmetricTensor<2, dim>>
          retrieve_stress_previous_timestep (const MaterialModel::MaterialModelInputs<dim> &in,
                                             const std::vector<Point<dim>> &quadrature_positions) const;

          double elastic_damper_viscosity;

          std::vector<double> elastic_shear_moduli;

          bool use_fixed_elastic_time_step;

          double fixed_elastic_time_step;

          double stabilization_time_scale_factor;

          mutable std::unique_ptr<FEPointEvaluation<dim, dim>> evaluator;
          static constexpr unsigned int n_independent_components = SymmetricTensor<2, dim>::n_independent_components;
          mutable std::unique_ptr<FEPointEvaluation<n_independent_components, dim>> evaluator_composition;

      };
    }
  }
}
#endif