elasticity.h

Go to the documentation of this file.

/*
  Copyright (C) 2019 - 2022 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
  {
    using namespace dealii;

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



    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_outputs (MaterialModel::MaterialModelOutputs<dim> &out) const;

          void
          fill_elastic_force_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;

          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;

          double
          calculate_viscoelastic_strain_rate (const SymmetricTensor<2,dim> &strain_rate,
                                              const SymmetricTensor<2,dim> &stored_stress,
                                              const double shear_modulus) const;

          double
          elastic_timestep () const;

        private:
          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;
      };
    }
  }
}
#endif