ASPECT
Public Member Functions | Static Public Member Functions | Private Attributes | List of all members
aspect::MaterialModel::Rheology::DruckerPrager< dim > Class Template Reference
Inheritance diagram for aspect::MaterialModel::Rheology::DruckerPrager< dim >:
Inheritance graph
[legend]

Public Member Functions

 DruckerPrager ()
 
void parse_parameters (ParameterHandler &prm, const std::unique_ptr< std::vector< unsigned int >> &expected_n_phases_per_composition=nullptr)
 
const DruckerPragerParameters compute_drucker_prager_parameters (const unsigned int composition, const std::vector< double > &phase_function_values=std::vector< double >(), const std::vector< unsigned int > &n_phases_per_composition=std::vector< unsigned int >()) const
 
double compute_yield_stress (const double cohesion, const double angle_internal_friction, const double pressure, const double max_yield_stress) const
 
double compute_viscosity (const double cohesion, const double angle_internal_friction, const double pressure, const double effective_strain_rate, const double max_yield_stress, const double pre_yield_viscosity=std::numeric_limits< double >::infinity()) const
 
std::pair< double, double > compute_strain_rate_and_derivative (const double stress, const double pressure, const DruckerPragerParameters p) const
 
double compute_derivative (const double angle_internal_friction, const double effective_strain_rate) const
 
- Public Member Functions inherited from aspect::SimulatorAccess< dim >
 SimulatorAccess ()
 
 SimulatorAccess (const Simulator< dim > &simulator_object)
 
virtual ~SimulatorAccess ()=default
 
virtual void initialize_simulator (const Simulator< dim > &simulator_object)
 
template<typename PostprocessorType >
PostprocessorType * find_postprocessor () const
 
const Introspection< dim > & introspection () const
 
const Simulator< dim > & get_simulator () const
 
const Parameters< dim > & get_parameters () const
 
SimulatorSignals< dim > & get_signals () const
 
MPI_Comm get_mpi_communicator () const
 
TimerOutput & get_computing_timer () const
 
const ConditionalOStream & get_pcout () const
 
double get_time () const
 
double get_timestep () const
 
double get_old_timestep () const
 
unsigned int get_timestep_number () const
 
unsigned int get_nonlinear_iteration () const
 
const parallel::distributed::Triangulation< dim > & get_triangulation () const
 
double get_volume () const
 
const Mapping< dim > & get_mapping () const
 
std::string get_output_directory () const
 
bool include_adiabatic_heating () const
 
bool include_latent_heat () const
 
bool include_melt_transport () const
 
int get_stokes_velocity_degree () const
 
double get_adiabatic_surface_temperature () const
 
double get_surface_pressure () const
 
bool convert_output_to_years () const
 
unsigned int get_pre_refinement_step () const
 
unsigned int n_compositional_fields () const
 
void get_refinement_criteria (Vector< float > &estimated_error_per_cell) const
 
void get_artificial_viscosity (Vector< float > &viscosity_per_cell, const bool skip_interior_cells=false) const
 
void get_artificial_viscosity_composition (Vector< float > &viscosity_per_cell, const unsigned int compositional_variable) const
 
const LinearAlgebra::BlockVectorget_current_linearization_point () const
 
const LinearAlgebra::BlockVectorget_solution () const
 
const LinearAlgebra::BlockVectorget_old_solution () const
 
const LinearAlgebra::BlockVectorget_old_old_solution () const
 
const LinearAlgebra::BlockVectorget_reaction_vector () const
 
const LinearAlgebra::BlockVectorget_mesh_velocity () const
 
const DoFHandler< dim > & get_dof_handler () const
 
const FiniteElement< dim > & get_fe () const
 
const LinearAlgebra::BlockSparseMatrixget_system_matrix () const
 
const LinearAlgebra::BlockSparseMatrixget_system_preconditioner_matrix () const
 
const MaterialModel::Interface< dim > & get_material_model () const
 
const GravityModel::Interface< dim > & get_gravity_model () const
 
const InitialTopographyModel::Interface< dim > & get_initial_topography_model () const
 
const GeometryModel::Interface< dim > & get_geometry_model () const
 
const AdiabaticConditions::Interface< dim > & get_adiabatic_conditions () const
 
bool has_boundary_temperature () const
 
DEAL_II_DEPRECATED const BoundaryTemperature::Interface< dim > & get_boundary_temperature () const
 
const BoundaryTemperature::Manager< dim > & get_boundary_temperature_manager () const
 
const BoundaryHeatFlux::Interface< dim > & get_boundary_heat_flux () const
 
bool has_boundary_composition () const
 
DEAL_II_DEPRECATED const BoundaryComposition::Interface< dim > & get_boundary_composition () const
 
const BoundaryComposition::Manager< dim > & get_boundary_composition_manager () const
 
const std::map< types::boundary_id, std::unique_ptr< BoundaryTraction::Interface< dim > > > & get_boundary_traction () const
 
DEAL_II_DEPRECATED const InitialTemperature::Interface< dim > & get_initial_temperature () const
 
std::shared_ptr< const InitialTemperature::Manager< dim > > get_initial_temperature_manager_pointer () const
 
const InitialTemperature::Manager< dim > & get_initial_temperature_manager () const
 
DEAL_II_DEPRECATED const InitialComposition::Interface< dim > & get_initial_composition () const
 
std::shared_ptr< const InitialComposition::Manager< dim > > get_initial_composition_manager_pointer () const
 
const InitialComposition::Manager< dim > & get_initial_composition_manager () const
 
const std::set< types::boundary_id > & get_fixed_temperature_boundary_indicators () const
 
const std::set< types::boundary_id > & get_fixed_heat_flux_boundary_indicators () const
 
const std::set< types::boundary_id > & get_fixed_composition_boundary_indicators () const
 
const std::set< types::boundary_id > & get_mesh_deformation_boundary_indicators () const
 
const BoundaryVelocity::Manager< dim > & get_boundary_velocity_manager () const
 
const HeatingModel::Manager< dim > & get_heating_model_manager () const
 
const MeshRefinement::Manager< dim > & get_mesh_refinement_manager () const
 
const MeltHandler< dim > & get_melt_handler () const
 
const VolumeOfFluidHandler< dim > & get_volume_of_fluid_handler () const
 
const NewtonHandler< dim > & get_newton_handler () const
 
const MeshDeformation::MeshDeformationHandler< dim > & get_mesh_deformation_handler () const
 
const LateralAveraging< dim > & get_lateral_averaging () const
 
const AffineConstraints< double > & get_current_constraints () const
 
bool simulator_is_past_initialization () const
 
double get_pressure_scaling () const
 
bool pressure_rhs_needs_compatibility_modification () const
 
bool model_has_prescribed_stokes_solution () const
 
TableHandler & get_statistics_object () const
 
template<typename PostprocessorType >
DEAL_II_DEPRECATED PostprocessorType * find_postprocessor () const
 
const Postprocess::Manager< dim > & get_postprocess_manager () const
 
const Particle::World< dim > & get_particle_world () const
 
Particle::World< dim > & get_particle_world ()
 
bool is_stokes_matrix_free ()
 
const StokesMatrixFreeHandler< dim > & get_stokes_matrix_free () const
 
RotationProperties< dim > compute_net_angular_momentum (const bool use_constant_density, const LinearAlgebra::BlockVector &solution, const bool limit_to_top_faces=false) const
 

Static Public Member Functions

static void declare_parameters (ParameterHandler &prm)
 
- Static Public Member Functions inherited from aspect::SimulatorAccess< dim >
static void get_composition_values_at_q_point (const std::vector< std::vector< double >> &composition_values, const unsigned int q, std::vector< double > &composition_values_at_q_point)
 

Private Attributes

std::vector< double > angles_internal_friction
 
std::vector< double > cohesions
 
double max_yield_stress
 
bool use_plastic_damper
 
double damper_viscosity
 

Detailed Description

template<int dim>
class aspect::MaterialModel::Rheology::DruckerPrager< dim >

Definition at line 60 of file drucker_prager.h.

Constructor & Destructor Documentation

§ DruckerPrager()

Member Function Documentation

§ declare_parameters()

template<int dim>
static void aspect::MaterialModel::Rheology::DruckerPrager< dim >::declare_parameters ( ParameterHandler &  prm)
static

Declare the parameters this function takes through input files.

§ parse_parameters()

template<int dim>
void aspect::MaterialModel::Rheology::DruckerPrager< dim >::parse_parameters ( ParameterHandler &  prm,
const std::unique_ptr< std::vector< unsigned int >> &  expected_n_phases_per_composition = nullptr 
)

Read the parameters this class declares from the parameter file. If expected_n_phases_per_composition points to a vector of unsigned integers, this is considered the number of phases for each compositional field (plus possibly a background field) and this number will be checked against the parsed parameters.

Parameters
[in]prmThe ParameterHandler to read from.
expected_n_phases_per_compositionOptional list of number of phases.

§ compute_drucker_prager_parameters()

template<int dim>
const DruckerPragerParameters aspect::MaterialModel::Rheology::DruckerPrager< dim >::compute_drucker_prager_parameters ( const unsigned int  composition,
const std::vector< double > &  phase_function_values = std::vector< double >(),
const std::vector< unsigned int > &  n_phases_per_composition = std::vector< unsigned int >() 
) const

§ compute_yield_stress()

template<int dim>
double aspect::MaterialModel::Rheology::DruckerPrager< dim >::compute_yield_stress ( const double  cohesion,
const double  angle_internal_friction,
const double  pressure,
const double  max_yield_stress 
) const

Compute the plastic yield stress based on the Drucker Prager yield criterion.

§ compute_viscosity()

template<int dim>
double aspect::MaterialModel::Rheology::DruckerPrager< dim >::compute_viscosity ( const double  cohesion,
const double  angle_internal_friction,
const double  pressure,
const double  effective_strain_rate,
const double  max_yield_stress,
const double  pre_yield_viscosity = std::numeric_limits< double >::infinity() 
) const

Compute the plastic viscosity with the yield stress and effective strain rate.

§ compute_strain_rate_and_derivative()

template<int dim>
std::pair<double, double> aspect::MaterialModel::Rheology::DruckerPrager< dim >::compute_strain_rate_and_derivative ( const double  stress,
const double  pressure,
const DruckerPragerParameters  p 
) const

Compute the strain rate and first stress derivative as a function of stress based on the damped Drucker-Prager flow law.

§ compute_derivative()

template<int dim>
double aspect::MaterialModel::Rheology::DruckerPrager< dim >::compute_derivative ( const double  angle_internal_friction,
const double  effective_strain_rate 
) const

Compute the derivative of the plastic viscosity with respect to pressure.

Member Data Documentation

§ angles_internal_friction

template<int dim>
std::vector<double> aspect::MaterialModel::Rheology::DruckerPrager< dim >::angles_internal_friction
private

Definition at line 128 of file drucker_prager.h.

§ cohesions

template<int dim>
std::vector<double> aspect::MaterialModel::Rheology::DruckerPrager< dim >::cohesions
private

Definition at line 129 of file drucker_prager.h.

§ max_yield_stress

template<int dim>
double aspect::MaterialModel::Rheology::DruckerPrager< dim >::max_yield_stress
private

Definition at line 130 of file drucker_prager.h.

§ use_plastic_damper

template<int dim>
bool aspect::MaterialModel::Rheology::DruckerPrager< dim >::use_plastic_damper
private

Whether to add a plastic damper in the computation of the drucker prager plastic viscosity.

Definition at line 136 of file drucker_prager.h.

§ damper_viscosity

template<int dim>
double aspect::MaterialModel::Rheology::DruckerPrager< dim >::damper_viscosity
private

Viscosity of a damper used to stabilize plasticity

Definition at line 141 of file drucker_prager.h.


The documentation for this class was generated from the following file: