ASPECT
drucker_prager.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2019 - 2021 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 
86  compute_drucker_prager_parameters (const unsigned int composition,
87  const std::vector<double> &phase_function_values = std::vector<double>(),
88  const std::vector<unsigned int> &n_phases_per_composition = std::vector<unsigned int>()) const;
89 
93  double
94  compute_yield_stress (const double cohesion,
95  const double angle_internal_friction,
96  const double pressure,
97  const double max_yield_stress) const;
98 
102  double
103  compute_viscosity (const double cohesion,
104  const double angle_internal_friction,
105  const double pressure,
106  const double effective_strain_rate,
107  const double max_yield_stress,
108  const double pre_yield_viscosity = std::numeric_limits<double>::infinity()) const;
109 
114  std::pair<double, double>
115  compute_strain_rate_and_derivative (const double stress,
116  const double pressure,
117  const DruckerPragerParameters p) const;
118 
122  double
123  compute_derivative (const double angle_internal_friction,
124  const double effective_strain_rate) const;
125 
126  private:
127 
128  std::vector<double> angles_internal_friction;
129  std::vector<double> cohesions;
131 
137 
142  };
143  }
144  }
145 }
146 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:88