strain_dependent.h

Go to the documentation of this file.

/*
  Copyright (C) 2019 - 2021 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_strain_dependent_h
#define _aspect_material_model_rheology_strain_dependent_h

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

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

namespace aspect
{
  namespace MaterialModel
  {
    using namespace dealii;

    namespace Rheology
    {
      enum WeakeningMechanism
      {
        none,
        finite_strain_tensor,
        total_strain,
        plastic_weakening_with_plastic_strain_only,
        plastic_weakening_with_total_strain_only,
        plastic_weakening_with_plastic_strain_and_viscous_weakening_with_viscous_strain,
        viscous_weakening_with_viscous_strain_only
      };

      enum HealingMechanism
      {
        no_healing,
        temperature_dependent
      };

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

          void
          parse_parameters (ParameterHandler &prm);

          std::array<double, 3>
          compute_strain_weakening_factors(const unsigned int j,
                                           const std::vector<double> &composition) const;

          double
          calculate_strain_healing (const MaterialModel::MaterialModelInputs<dim> &in,
                                    const unsigned int j) const;

          std::pair<double, double>
          calculate_plastic_weakening (const double strain_ii,
                                       const unsigned int j) const;

          double
          calculate_viscous_weakening (const double strain_ii,
                                       const unsigned int j) const;

          void compute_finite_strain_reaction_terms (const MaterialModel::MaterialModelInputs<dim> &in,
                                                     MaterialModel::MaterialModelOutputs<dim> &out) const;
          void fill_reaction_outputs (const MaterialModel::MaterialModelInputs<dim> &in,
                                      const int i,
                                      const double min_strain_rate,
                                      const bool plastic_yielding,
                                      MaterialModel::MaterialModelOutputs<dim> &out) const;

          ComponentMask get_strain_composition_mask() const;

          WeakeningMechanism
          get_weakening_mechanism () const;

          HealingMechanism
          get_healing_mechanism () const;

        private:

          WeakeningMechanism weakening_mechanism;

          HealingMechanism healing_mechanism;

          std::vector<double> start_plastic_strain_weakening_intervals;

          std::vector<double> end_plastic_strain_weakening_intervals;

          std::vector<double> cohesion_strain_weakening_factors;

          std::vector<double> friction_strain_weakening_factors;

          std::vector<double> start_viscous_strain_weakening_intervals;

          std::vector<double> end_viscous_strain_weakening_intervals;

          std::vector<double> viscous_strain_weakening_factors;

          double strain_healing_temperature_dependent_recovery_rate;

          double strain_healing_temperature_dependent_prefactor;

          mutable std::unique_ptr<FEPointEvaluation<dim, dim>> evaluator;
          mutable std::vector<std::unique_ptr<FEPointEvaluation<1, dim>>> composition_evaluators;
      };
    }
  }
}
#endif