elasticity.h

Go to the documentation of this file.

/*
  Copyright (C) 2019 - 2024 by the authors of the ASPECT code.

  This file is part of ASPECT.

  ASPECT 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 2, or (at your option)
  any later version.

  ASPECT 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 ASPECT; see the file LICENSE.  If not see
  <http://www.gnu.org/licenses/>.
*/

#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;

          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