21 #ifndef _aspect_material_model_rheology_elasticity_h 22 #define _aspect_material_model_rheology_elasticity_h 28 #include <deal.II/matrix_free/fe_point_evaluation.h> 32 namespace MaterialModel
45 std::vector<double>
get_nth_output(
const unsigned int idx)
const override;
89 parse_parameters (ParameterHandler &prm);
121 const std::vector<double> &average_elastic_shear_moduli,
145 const std::vector<double> &average_elastic_shear_moduli,
169 const std::vector<double> &average_elastic_shear_moduli,
187 const std::vector<double> &average_elastic_shear_moduli,
194 const std::vector<double> &
195 get_elastic_shear_moduli ()
const;
202 calculate_elastic_viscosity (
const double shear_modulus)
const;
209 calculate_viscoelastic_viscosity (
const double viscosity,
210 const double shear_modulus)
const;
229 SymmetricTensor<2,dim>
230 calculate_viscoelastic_strain_rate (
const SymmetricTensor<2,dim> &
strain_rate,
231 const SymmetricTensor<2, dim> &stress_0_advected,
232 const SymmetricTensor<2, dim> &stress_old,
233 const double viscosity_pre_yield,
234 const double shear_modulus)
const;
240 elastic_timestep ()
const;
247 calculate_timestep_ratio()
const;
258 std::vector<SymmetricTensor<2, dim>>
260 const std::vector<Point<dim>> &quadrature_positions)
const;
301 mutable std::unique_ptr<FEPointEvaluation<dim, dim>>
evaluator;
302 static constexpr
unsigned int n_independent_components = SymmetricTensor<2, dim>::n_independent_components;
double stabilization_time_scale_factor
std::vector< double > elastic_shear_moduli
std::vector< double > elastic_viscosity
void declare_parameters(ParameterHandler &prm)
double elastic_damper_viscosity
bool use_fixed_elastic_time_step
std::unique_ptr< FEPointEvaluation< n_independent_components, dim > > evaluator_composition
std::unique_ptr< FEPointEvaluation< dim, dim > > evaluator
std::vector< double > elastic_shear_moduli
std::vector< SymmetricTensor< 2, dim > > deviatoric_stress
double fixed_elastic_time_step
std::vector< double > get_nth_output(const unsigned int idx) const override
ElasticAdditionalOutputs(const unsigned int n_points)