ASPECT
drucker_prager_power.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2019 - 2024 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_power_h
22 #define _aspect_material_model_rheology_drucker_prager_power_h
23 
24 #include <aspect/global.h>
26 
30 
31 namespace aspect
32 {
33  namespace MaterialModel
34  {
35  using namespace dealii;
36 
37  namespace Rheology
38  {
39 
40  template <int dim>
42  {
56  public:
58 
62  static
63  void
64  declare_parameters (ParameterHandler &prm);
65 
76  void
77  parse_parameters (ParameterHandler &prm,
78  const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition = nullptr);
79 
88  compute_drucker_prager_parameters (const unsigned int composition,
89  const std::vector<double> &phase_function_values = std::vector<double>(),
90  const std::vector<unsigned int> &n_phase_transitions_per_composition = std::vector<unsigned int>()) const;
91 
95  double
96  compute_yield_stress (const double cohesion,
97  const double angle_internal_friction,
98  const double pressure,
99  const double max_yield_stress) const;
100 
104  double
105  compute_viscosity (const double cohesion,
106  const double angle_internal_friction,
107  const double pressure,
108  const double effective_strain_rate,
109  const double max_yield_stress) const;
110 
115  std::pair<double, double>
116  compute_strain_rate_and_derivative (const double stress,
117  const double pressure,
118  const DruckerPragerParameters &p) const;
119 
124  std::pair<double, double>
125  compute_log_strain_rate_and_derivative (const double log_stress,
126  const double pressure,
127  const DruckerPragerParameters &p) const;
128 
129  private:
134  std::vector<double> angles_internal_friction;
135 
139  std::vector<double> cohesions;
140 
145 
151 
157 
162 
163  };
164  }
165  }
166 }
167 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:59
Definition: compat.h:42