From: Daniel Arndt Date: Tue, 13 Mar 2018 23:49:35 +0000 (+0100) Subject: Fix Quasi_static_Finite_strain_Quasi_incompressible_ViscoElasticity X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=218e076dfce6e3debefe0632ddd947d19c51ab1c;p=code-gallery.git Fix Quasi_static_Finite_strain_Quasi_incompressible_ViscoElasticity --- diff --git a/Quasi_static_Finite_strain_Quasi_incompressible_ViscoElasticity/viscoelastic_strip_with_hole.cc b/Quasi_static_Finite_strain_Quasi_incompressible_ViscoElasticity/viscoelastic_strip_with_hole.cc index 73772de..b69a651 100644 --- a/Quasi_static_Finite_strain_Quasi_incompressible_ViscoElasticity/viscoelastic_strip_with_hole.cc +++ b/Quasi_static_Finite_strain_Quasi_incompressible_ViscoElasticity/viscoelastic_strip_with_hole.cc @@ -49,7 +49,6 @@ #include #include #include -#include #include #include #include @@ -295,7 +294,7 @@ namespace ViscoElasStripHole prm.enter_subsection("Linear solver"); { prm.declare_entry("Solver type", "cg", - Patterns::Selection(SolverSelector::get_solver_names()), + Patterns::Selection(SolverSelector::get_solver_names()), "Type of solver used to solve the linear system"); prm.declare_entry("Residual", "1e-6", Patterns::Double(0.0), @@ -512,8 +511,8 @@ namespace ViscoElasStripHole } void update_internal_equilibrium(const Tensor<2, dim> &F, - const double &p_tilde, - const double &J_tilde) + const double &/*p_tilde*/, + const double &/*J_tilde*/) { const double det_F = determinant(F); const SymmetricTensor<2,dim> C_bar = std::pow(det_F, -2.0 / dim) * Physics::Elasticity::Kinematics::C(F); @@ -584,7 +583,7 @@ namespace ViscoElasStripHole + Physics::Elasticity::StandardTensors::dev_P * c_bar * Physics::Elasticity::StandardTensors::dev_P; } - SymmetricTensor<4, dim> get_c_bar(const Tensor<2,dim> &F) const + SymmetricTensor<4, dim> get_c_bar(const Tensor<2,dim> &/*F*/) const { // Elastic Neo-Hookean + Linder2011 eq 56 return -2.0*mu_v*((time.get_delta_t()/tau_v)/(1.0 + time.get_delta_t()/tau_v))*Physics::Elasticity::StandardTensors::S; @@ -643,7 +642,7 @@ namespace ViscoElasStripHole material->update_end_timestep(); } private: - std_cxx11::shared_ptr< Material_Compressible_Three_Field_Linear_Viscoelastic > material; + std::shared_ptr< Material_Compressible_Three_Field_Linear_Viscoelastic > material; }; template class Solid @@ -896,7 +895,7 @@ namespace ViscoElasStripHole for (unsigned int d=0; d