21 #ifndef _aspect_material_model_rheology_drucker_prager_power_h 22 #define _aspect_material_model_rheology_drucker_prager_power_h 33 namespace MaterialModel
77 parse_parameters (ParameterHandler &prm,
78 const std::unique_ptr<std::vector<unsigned int>> &expected_n_phases_per_composition =
nullptr);
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;
96 compute_yield_stress (
const double cohesion,
97 const double angle_internal_friction,
99 const double max_yield_stress)
const;
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;
115 std::pair<double, double>
116 compute_strain_rate_and_derivative (
const double stress,
117 const double pressure,
124 std::pair<double, double>
125 compute_log_strain_rate_and_derivative (
const double log_stress,
126 const double pressure,
double drucker_prager_edot_ref
void declare_parameters(ParameterHandler &prm)
double drucker_prager_stress_exponent
std::vector< double > cohesions
double drucker_prager_log_edot_ref
std::vector< double > angles_internal_friction