ASPECT
visco_plastic.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2020 - 2022 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_visco_plastic_h
22 #define _aspect_material_model_rheology_visco_plastic_h
23 
24 #include <aspect/global.h>
37 
38 #include<deal.II/fe/component_mask.h>
39 
40 namespace aspect
41 {
42  namespace MaterialModel
43  {
44  using namespace dealii;
45 
51  template <int dim>
53  {
54  public:
55  PlasticAdditionalOutputs(const unsigned int n_points);
56 
57  std::vector<double> get_nth_output(const unsigned int idx) const override;
58 
64  std::vector<double> cohesions;
65 
71  std::vector<double> friction_angles;
72 
77  std::vector<double> yielding;
78  };
79 
84  {
88  std::vector<double> composition_viscosities;
89 
93  std::vector<bool> composition_yielding;
94 
98  std::vector<double> current_friction_angles;
99 
103  std::vector<double> current_cohesions;
104  };
105 
106  namespace Rheology
107  {
108 
109  template <int dim>
111  {
112  public:
116  ViscoPlastic();
117 
123  calculate_isostrain_viscosities ( const MaterialModel::MaterialModelInputs<dim> &in,
124  const unsigned int i,
125  const std::vector<double> &volume_fractions,
126  const std::vector<double> &phase_function_values = std::vector<double>(),
127  const std::vector<unsigned int> &n_phases_per_composition =
128  std::vector<unsigned int>()) const;
129 
135  void compute_viscosity_derivatives(const unsigned int point_index,
136  const std::vector<double> &volume_fractions,
137  const std::vector<double> &composition_viscosities,
140  const std::vector<double> &phase_function_values = std::vector<double>(),
141  const std::vector<unsigned int> &n_phases_per_composition =
142  std::vector<unsigned int>()) const;
143 
150  ComponentMask get_volumetric_composition_mask() const;
151 
155  static
156  void
157  declare_parameters (ParameterHandler &prm);
158 
166  void
167  parse_parameters (ParameterHandler &prm,
168  const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition = nullptr);
169 
174  void
175  create_plastic_outputs (MaterialModel::MaterialModelOutputs<dim> &out) const;
176 
182  void fill_plastic_outputs(const unsigned int point_index,
183  const std::vector<double> &volume_fractions,
184  const bool plastic_yielding,
187  const IsostrainViscosities &isostrain_viscosities) const;
188 
193 
198 
203 
208 
213 
218 
219 
220  private:
221 
227 
233  std::vector<double> minimum_viscosity;
234  std::vector<double> maximum_viscosity;
235 
241  {
245  composite
246  } viscous_flow_law;
247 
253  {
255  drucker_prager
256  } yield_mechanism;
257 
264 
274 
279  std::vector<double> exponents_stress_limiter;
280 
285 
291  std::unique_ptr<Rheology::FrankKamenetskii<dim>> frank_kamenetskii_rheology;
292 
297 
301  std::unique_ptr<Rheology::PeierlsCreep<dim>> peierls_creep;
302 
309 
310  /*
311  * Object for computing plastic stresses, viscosities, and additional outputs
312  */
314 
315  /*
316  * Input parameters for the drucker prager plasticity.
317  */
319 
320  };
321  }
322  }
323 }
324 #endif
MaterialUtilities::CompositionalAveragingOperation viscosity_averaging
Rheology::FrictionModels< dim > friction_models
std::unique_ptr< Rheology::FrankKamenetskii< dim > > frank_kamenetskii_rheology
Rheology::DruckerPrager< dim > drucker_prager_plasticity
void declare_parameters(ParameterHandler &prm)
Rheology::ConstantViscosityPrefactors< dim > constant_viscosity_prefactors
Rheology::DislocationCreep< dim > dislocation_creep
Rheology::StrainDependent< dim > strain_rheology
Rheology::DiffusionCreep< dim > diffusion_creep
Definition: compat.h:88
std::unique_ptr< Rheology::PeierlsCreep< dim > > peierls_creep
Rheology::DruckerPragerParameters drucker_prager_parameters
Rheology::Elasticity< dim > elastic_rheology