ASPECT
drucker_prager.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2019 - 2023 by the authors of the ASPECT code.
3 
4  This file is part of ASPECT.
5 
6  ASPECT is free software; you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation; either version 2, or (at your option)
9  any later version.
10 
11  ASPECT is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with ASPECT; see the file LICENSE. If not see
18  <http://www.gnu.org/licenses/>.
19 */
20 
21 #ifndef _aspect_material_model_rheology_drucker_prager_h
22 #define _aspect_material_model_rheology_drucker_prager_h
23 
24 #include <aspect/global.h>
26 
29 namespace aspect
30 {
31  namespace MaterialModel
32  {
33  using namespace dealii;
34 
35  namespace Rheology
36  {
42  {
47 
51  double cohesion;
52 
57  };
58 
59  template <int dim>
61  {
62  public:
63  DruckerPrager();
67  static
68  void
69  declare_parameters (ParameterHandler &prm);
70 
81  void
82  parse_parameters (ParameterHandler &prm,
83  const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition = nullptr);
84 
93  compute_drucker_prager_parameters (const unsigned int composition,
94  const std::vector<double> &phase_function_values = std::vector<double>(),
95  const std::vector<unsigned int> &n_phase_transitions_per_composition = std::vector<unsigned int>()) const;
96 
100  double
101  compute_yield_stress (const double cohesion,
102  const double angle_internal_friction,
103  const double pressure,
104  const double max_yield_stress) const;
105 
112  double
113  compute_viscosity (const double cohesion,
114  const double angle_internal_friction,
115  const double pressure,
116  const double effective_strain_rate,
117  const double max_yield_stress,
118  const double non_yielding_viscosity = std::numeric_limits<double>::infinity()) const;
119 
124  std::pair<double, double>
125  compute_strain_rate_and_derivative (const double stress,
126  const double pressure,
127  const DruckerPragerParameters p) const;
128 
132  double
133  compute_derivative (const double angle_internal_friction,
134  const double effective_strain_rate) const;
135 
136  private:
137 
138  std::vector<double> angles_internal_friction;
139  std::vector<double> cohesions;
141 
147 
152  };
153  }
154  }
155 }
156 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:42