strain_dependent.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_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
  {
    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 std::vector<double> &composition,
                                           const unsigned int j) const;

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

          std::array<double, 3>
          apply_temperature_dependence_to_strain_weakening_factors(const std::array<double, 3> &weakening_factors,
                                                                   const double temperature,
                                                                   const unsigned int j) 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;

          bool use_temperature_activated_strain_softening;

          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;

          std::vector<double> viscous_strain_weakening_T0;
          std::vector<double> viscous_strain_weakening_T1;
          std::vector<double> viscous_strain_weakening_T2;
          std::vector<double> viscous_strain_weakening_T3;

          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