ASPECT
visco_plastic.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2011 - 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_visco_plastic_h
22 #define _aspect_material_model_visco_plastic_h
23 
28 
29 #include<deal.II/fe/component_mask.h>
30 
31 namespace aspect
32 {
33  namespace MaterialModel
34  {
35  using namespace dealii;
36 
180  template <int dim>
182  {
183  public:
184 
185  void evaluate(const MaterialModel::MaterialModelInputs<dim> &in,
186  MaterialModel::MaterialModelOutputs<dim> &out) const override;
187 
198  bool is_compressible () const override;
199 
200  static
201  void
202  declare_parameters (ParameterHandler &prm);
203 
204  void
205  parse_parameters (ParameterHandler &prm) override;
206 
207  void
208  create_additional_named_outputs (MaterialModel::MaterialModelOutputs<dim> &out) const override;
209 
210  double get_min_strain_rate() const;
211 
220  bool
221  is_yielding (const double pressure,
222  const double temperature,
223  const std::vector<double> &composition,
224  const SymmetricTensor<2,dim> &strain_rate) const;
225 
231  bool
232  is_yielding (const MaterialModelInputs<dim> &in) const;
233 
234  private:
235 
247  std::unique_ptr<Rheology::ViscoPlastic<dim>> rheology;
248 
249  std::vector<double> thermal_diffusivities;
250 
255 
256  std::vector<double> thermal_conductivities;
257 
262 
267 
268  };
269 
270  }
271 }
272 
273 #endif
MaterialUtilities::PhaseFunction< dim > phase_function
std::unique_ptr< Rheology::ViscoPlastic< dim > > rheology
StructuredDataLookup< dim > DEAL_II_DEPRECATED
EquationOfState::MulticomponentIncompressible< dim > equation_of_state
std::vector< double > thermal_diffusivities
void declare_parameters(ParameterHandler &prm)
std::vector< double > thermal_conductivities
Definition: compat.h:88