Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | Static Private Attributes | List of all members
oomph::AxisymFoepplvonKarmanEquations Class Referenceabstract

#include <axisym_displ_based_fvk_elements.h>

+ Inheritance diagram for oomph::AxisymFoepplvonKarmanEquations:

Public Types

typedef void(* AxisymFoepplvonKarmanPressureFctPt) (const double &r, double &f)
 
typedef void(* AxisymFoepplvonKarmanPressureFctPt) (const double &r, double &f)
 
- Public Types inherited from oomph::FiniteElement
typedef void(* SteadyExactSolutionFctPt) (const Vector< double > &, Vector< double > &)
 
typedef void(* UnsteadyExactSolutionFctPt) (const double &, const Vector< double > &, Vector< double > &)
 

Public Member Functions

 AxisymFoepplvonKarmanEquations ()
 
 AxisymFoepplvonKarmanEquations (const AxisymFoepplvonKarmanEquations &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const AxisymFoepplvonKarmanEquations &)=delete
 Broken assignment operator. More...
 
const double & nu () const
 Poisson's ratio. More...
 
double *& nu_pt ()
 Pointer to Poisson's ratio. More...
 
const double & eta () const
 FvK parameter. More...
 
double *& eta_pt ()
 Pointer to FvK parameter. More...
 
virtual unsigned nodal_index_fvk (const unsigned &i=0) const
 
void output (std::ostream &outfile)
 Output with default number of plot points. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 
void output (FILE *file_pt)
 C_style output with default number of plot points. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 
void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output exact soln: r,w_exact at n_plot plot points. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Get error against and norm of exact solution. More...
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Dummy, time dependent error checker. More...
 
AxisymFoepplvonKarmanPressureFctPt & pressure_fct_pt ()
 Access function: Pointer to pressure function. More...
 
AxisymFoepplvonKarmanPressureFctPt pressure_fct_pt () const
 Access function: Pointer to pressure function. Const version. More...
 
virtual void get_pressure_fvk (const unsigned &ipt, const double &r, double &pressure) const
 
void get_gradient_of_deflection (const Vector< double > &s, Vector< double > &gradient) const
 Get gradient of deflection: gradient[i] = dw/dr_i *‍/. More...
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Fill in the residuals with this element's contribution. More...
 
double interpolated_w_fvk (const Vector< double > &s) const
 Return FE representation of transverse displacement. More...
 
double interpolated_u_fvk (const Vector< double > &s) const
 Return FE representation of radial displacement. More...
 
bool interpolated_stress (const Vector< double > &s, double &sigma_r_r, double &sigma_phi_phi) const
 
unsigned self_test ()
 Self-test: Return 0 for OK. More...
 
void use_linear_bending_model ()
 
 AxisymFoepplvonKarmanEquations ()
 
 AxisymFoepplvonKarmanEquations (const AxisymFoepplvonKarmanEquations &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const AxisymFoepplvonKarmanEquations &)=delete
 Broken assignment operator. More...
 
const double & eta () const
 FvK parameter. More...
 
double *& eta_pt ()
 Pointer to FvK parameter. More...
 
virtual unsigned nodal_index_fvk (const unsigned &i=0) const
 
void output (std::ostream &outfile)
 Output with default number of plot points. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 
void output (FILE *file_pt)
 C_style output with default number of plot points. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 
void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output exact soln: r,w_exact at n_plot plot points. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Get error against and norm of exact solution. More...
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Dummy, time dependent error checker. More...
 
AxisymFoepplvonKarmanPressureFctPt & pressure_fct_pt ()
 Access function: Pointer to pressure function. More...
 
AxisymFoepplvonKarmanPressureFctPt pressure_fct_pt () const
 Access function: Pointer to pressure function. Const version. More...
 
AxisymFoepplvonKarmanPressureFctPt & airy_forcing_fct_pt ()
 Access function: Pointer to Airy forcing function. More...
 
AxisymFoepplvonKarmanPressureFctPt airy_forcing_fct_pt () const
 Access function: Pointer to Airy forcing function. Const version. More...
 
virtual void get_pressure_fvk (const unsigned &ipt, const double &r, double &pressure) const
 
virtual void get_airy_forcing_fvk (const unsigned &ipt, const double &r, double &airy_forcing) const
 
void get_gradient_of_deflection (const Vector< double > &s, Vector< double > &gradient) const
 Get gradient of deflection: gradient[i] = dw/dr_i *‍/. More...
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Fill in the residuals with this element's contribution. More...
 
double interpolated_w_fvk (const Vector< double > &s) const
 
bool interpolated_stress (const Vector< double > &s, double &sigma_r_r, double &sigma_phi_phi)
 
unsigned self_test ()
 Self-test: Return 0 for OK. More...
 
void use_linear_bending_model ()
 
- Public Member Functions inherited from oomph::FiniteElement
void set_dimension (const unsigned &dim)
 
void set_nodal_dimension (const unsigned &nodal_dim)
 
void set_nnodal_position_type (const unsigned &nposition_type)
 Set the number of types required to interpolate the coordinate. More...
 
void set_n_node (const unsigned &n)
 
int nodal_local_eqn (const unsigned &n, const unsigned &i) const
 
double dJ_eulerian_at_knot (const unsigned &ipt, Shape &psi, DenseMatrix< double > &djacobian_dX) const
 
 FiniteElement ()
 Constructor. More...
 
virtual ~FiniteElement ()
 
 FiniteElement (const FiniteElement &)=delete
 Broken copy constructor. More...
 
virtual bool local_coord_is_valid (const Vector< double > &s)
 Broken assignment operator. More...
 
virtual void move_local_coord_back_into_element (Vector< double > &s) const
 
void get_centre_of_gravity_and_max_radius_in_terms_of_zeta (Vector< double > &cog, double &max_radius) const
 
virtual void local_coordinate_of_node (const unsigned &j, Vector< double > &s) const
 
virtual void local_fraction_of_node (const unsigned &j, Vector< double > &s_fraction)
 
virtual double local_one_d_fraction_of_node (const unsigned &n1d, const unsigned &i)
 
virtual void set_macro_elem_pt (MacroElement *macro_elem_pt)
 
MacroElement * macro_elem_pt ()
 Access function to pointer to macro element. More...
 
void get_x (const Vector< double > &s, Vector< double > &x) const
 
void get_x (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void get_x_from_macro_element (const Vector< double > &s, Vector< double > &x) const
 
virtual void get_x_from_macro_element (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void set_integration_scheme (Integral *const &integral_pt)
 Set the spatial integration scheme. More...
 
Integral *const & integral_pt () const
 Return the pointer to the integration scheme (const version) More...
 
virtual void shape (const Vector< double > &s, Shape &psi) const =0
 
virtual void shape_at_knot (const unsigned &ipt, Shape &psi) const
 
virtual void dshape_local (const Vector< double > &s, Shape &psi, DShape &dpsids) const
 
virtual void dshape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids) const
 
virtual void d2shape_local (const Vector< double > &s, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
virtual void d2shape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
virtual double J_eulerian (const Vector< double > &s) const
 
virtual double J_eulerian_at_knot (const unsigned &ipt) const
 
void check_J_eulerian_at_knots (bool &passed) const
 
void check_jacobian (const double &jacobian) const
 
double dshape_eulerian (const Vector< double > &s, Shape &psi, DShape &dpsidx) const
 
virtual double dshape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsidx) const
 
virtual double dshape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsi, DenseMatrix< double > &djacobian_dX, RankFourTensor< double > &d_dpsidx_dX) const
 
double d2shape_eulerian (const Vector< double > &s, Shape &psi, DShape &dpsidx, DShape &d2psidx) const
 
virtual double d2shape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsidx, DShape &d2psidx) const
 
virtual void assign_nodal_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void describe_local_dofs (std::ostream &out, const std::string &current_string) const
 
virtual void describe_nodal_local_dofs (std::ostream &out, const std::string &current_string) const
 
virtual void assign_all_generic_local_eqn_numbers (const bool &store_local_dof_pt)
 
Node *& node_pt (const unsigned &n)
 Return a pointer to the local node n. More...
 
Node *const & node_pt (const unsigned &n) const
 Return a pointer to the local node n (const version) More...
 
unsigned nnode () const
 Return the number of nodes. More...
 
virtual unsigned nnode_1d () const
 
double raw_nodal_position (const unsigned &n, const unsigned &i) const
 
double raw_nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position (const unsigned &n, const unsigned &i) const
 
double nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double dnodal_position_dt (const unsigned &n, const unsigned &i) const
 Return the i-th component of nodal velocity: dx/dt at local node n. More...
 
double dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
virtual void disable_ALE ()
 
virtual void enable_ALE ()
 
virtual unsigned required_nvalue (const unsigned &n) const
 
unsigned nnodal_position_type () const
 
bool has_hanging_nodes () const
 
unsigned nodal_dimension () const
 Return the required Eulerian dimension of the nodes in this element. More...
 
virtual unsigned nvertex_node () const
 
virtual Node * vertex_node_pt (const unsigned &j) const
 
virtual Node * construct_node (const unsigned &n)
 
virtual Node * construct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
virtual Node * construct_boundary_node (const unsigned &n)
 
virtual Node * construct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
int get_node_number (Node *const &node_pt) const
 
virtual Node * get_node_at_local_coordinate (const Vector< double > &s) const
 
double raw_nodal_value (const unsigned &n, const unsigned &i) const
 
double raw_nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
unsigned dim () const
 
virtual ElementGeometry::ElementGeometry element_geometry () const
 Return the geometry type of the element (either Q or T usually). More...
 
virtual double interpolated_x (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated coordinate x[i] at local coordinate s. More...
 
virtual double interpolated_x (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
virtual void interpolated_x (const Vector< double > &s, Vector< double > &x) const
 Return FE interpolated position x[] at local coordinate s as Vector. More...
 
virtual void interpolated_x (const unsigned &t, const Vector< double > &s, Vector< double > &x) const
 
virtual double interpolated_dxdt (const Vector< double > &s, const unsigned &i, const unsigned &t)
 
virtual void interpolated_dxdt (const Vector< double > &s, const unsigned &t, Vector< double > &dxdt)
 
unsigned ngeom_data () const
 
Data * geom_data_pt (const unsigned &j)
 
void position (const Vector< double > &zeta, Vector< double > &r) const
 
void position (const unsigned &t, const Vector< double > &zeta, Vector< double > &r) const
 
void dposition_dt (const Vector< double > &zeta, const unsigned &t, Vector< double > &drdt)
 
virtual double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
void interpolated_zeta (const Vector< double > &s, Vector< double > &zeta) const
 
void locate_zeta (const Vector< double > &zeta, GeomObject *&geom_object_pt, Vector< double > &s, const bool &use_coordinate_as_initial_guess=false)
 
virtual void node_update ()
 
virtual void identify_field_data_for_interactions (std::set< std::pair< Data *, unsigned >> &paired_field_data)
 
virtual void identify_geometric_data (std::set< Data * > &geometric_data_pt)
 
virtual double s_min () const
 Min value of local coordinate. More...
 
virtual double s_max () const
 Max. value of local coordinate. More...
 
double size () const
 
virtual double compute_physical_size () const
 
virtual void point_output_data (const Vector< double > &s, Vector< double > &data)
 
void point_output (std::ostream &outfile, const Vector< double > &s)
 
virtual unsigned nplot_points_paraview (const unsigned &nplot) const
 
virtual unsigned nsub_elements_paraview (const unsigned &nplot) const
 
void output_paraview (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_output_offset_information (std::ofstream &file_out, const unsigned &nplot, unsigned &counter) const
 
virtual void write_paraview_type (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_offsets (std::ofstream &file_out, const unsigned &nplot, unsigned &offset_sum) const
 
virtual unsigned nscalar_paraview () const
 
virtual void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
virtual std::string scalar_name_paraview (const unsigned &i) const
 
virtual void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, const SolutionFunctorBase &exact_soln) const
 Output a time-dependent exact solution over the element. More...
 
virtual void get_s_plot (const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
 
virtual std::string tecplot_zone_string (const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
virtual unsigned nplot_points (const unsigned &nplot) const
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_abs_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error)
 
void integrate_fct (FiniteElement::SteadyExactSolutionFctPt integrand_fct_pt, Vector< double > &integral)
 Integrate Vector-valued function over element. More...
 
void integrate_fct (FiniteElement::UnsteadyExactSolutionFctPt integrand_fct_pt, const double &time, Vector< double > &integral)
 Integrate Vector-valued time-dep function over element. More...
 
virtual void build_face_element (const int &face_index, FaceElement *face_element_pt)
 
virtual unsigned get_bulk_node_number (const int &face_index, const unsigned &i) const
 
virtual int face_outer_unit_normal_sign (const int &face_index) const
 Get the sign of the outer unit normal on the face given by face_index. More...
 
virtual unsigned nnode_on_face () const
 
void face_node_number_error_check (const unsigned &i) const
 Range check for face node numbers. More...
 
virtual CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt (const int &face_index) const
 
virtual BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt (const int &face_index) const
 
- Public Member Functions inherited from oomph::GeneralisedElement
 GeneralisedElement ()
 Constructor: Initialise all pointers and all values to zero. More...
 
virtual ~GeneralisedElement ()
 Virtual destructor to clean up any memory allocated by the object. More...
 
 GeneralisedElement (const GeneralisedElement &)=delete
 Broken copy constructor. More...
 
void operator= (const GeneralisedElement &)=delete
 Broken assignment operator. More...
 
Data *& internal_data_pt (const unsigned &i)
 Return a pointer to i-th internal data object. More...
 
Data *const & internal_data_pt (const unsigned &i) const
 Return a pointer to i-th internal data object (const version) More...
 
Data *& external_data_pt (const unsigned &i)
 Return a pointer to i-th external data object. More...
 
Data *const & external_data_pt (const unsigned &i) const
 Return a pointer to i-th external data object (const version) More...
 
unsigned long eqn_number (const unsigned &ieqn_local) const
 
int local_eqn_number (const unsigned long &ieqn_global) const
 
unsigned add_external_data (Data *const &data_pt, const bool &fd=true)
 
bool external_data_fd (const unsigned &i) const
 
void exclude_external_data_fd (const unsigned &i)
 
void include_external_data_fd (const unsigned &i)
 
void flush_external_data ()
 Flush all external data. More...
 
void flush_external_data (Data *const &data_pt)
 Flush the object addressed by data_pt from the external data array. More...
 
unsigned ninternal_data () const
 Return the number of internal data objects. More...
 
unsigned nexternal_data () const
 Return the number of external data objects. More...
 
unsigned ndof () const
 Return the number of equations/dofs in the element. More...
 
void dof_vector (const unsigned &t, Vector< double > &dof)
 Return the vector of dof values at time level t. More...
 
void dof_pt_vector (Vector< double * > &dof_pt)
 Return the vector of pointers to dof values. More...
 
void set_internal_data_time_stepper (const unsigned &i, TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data)
 
void assign_internal_eqn_numbers (unsigned long &global_number, Vector< double * > &Dof_pt)
 
void describe_dofs (std::ostream &out, const std::string &current_string) const
 
void add_internal_value_pt_to_map (std::map< unsigned, double * > &map_of_value_pt)
 
virtual void assign_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void complete_setup_of_dependencies ()
 
virtual void get_residuals (Vector< double > &residuals)
 
virtual void get_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
virtual void get_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void get_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void get_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void get_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void get_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void get_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void get_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
virtual void compute_norm (Vector< double > &norm)
 
virtual void compute_norm (double &norm)
 
virtual unsigned ndof_types () const
 
virtual void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- Public Member Functions inherited from oomph::GeomObject
 GeomObject ()
 Default constructor. More...
 
 GeomObject (const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim, TimeStepper *time_stepper_pt)
 
 GeomObject (const GeomObject &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const GeomObject &)=delete
 Broken assignment operator. More...
 
virtual ~GeomObject ()
 (Empty) destructor More...
 
unsigned nlagrangian () const
 Access function to # of Lagrangian coordinates. More...
 
unsigned ndim () const
 Access function to # of Eulerian coordinates. More...
 
void set_nlagrangian_and_ndim (const unsigned &n_lagrangian, const unsigned &n_dim)
 Set # of Lagrangian and Eulerian coordinates. More...
 
TimeStepper *& time_stepper_pt ()
 
TimeStepper * time_stepper_pt () const
 
virtual void position (const double &t, const Vector< double > &zeta, Vector< double > &r) const
 
virtual void dposition (const Vector< double > &zeta, DenseMatrix< double > &drdzeta) const
 
virtual void d2position (const Vector< double > &zeta, RankThreeTensor< double > &ddrdzeta) const
 
virtual void d2position (const Vector< double > &zeta, Vector< double > &r, DenseMatrix< double > &drdzeta, RankThreeTensor< double > &ddrdzeta) const
 

Protected Member Functions

virtual double dshape_and_dtest_eulerian_axisym_fvk (const Vector< double > &s, Shape &psi, DShape &dpsidr, Shape &test, DShape &dtestdr) const =0
 
virtual double dshape_and_dtest_eulerian_at_knot_axisym_fvk (const unsigned &ipt, Shape &psi, DShape &dpsidr, Shape &test, DShape &dtestdr) const =0
 
virtual double dshape_and_dtest_eulerian_axisym_fvk (const Vector< double > &s, Shape &psi, DShape &dpsidr, Shape &test, DShape &dtestdr) const =0
 
virtual double dshape_and_dtest_eulerian_at_knot_axisym_fvk (const unsigned &ipt, Shape &psi, DShape &dpsidr, Shape &test, DShape &dtestdr) const =0
 
- Protected Member Functions inherited from oomph::FiniteElement
virtual void assemble_local_to_eulerian_jacobian (const DShape &dpsids, DenseMatrix< double > &jacobian) const
 
virtual void assemble_local_to_eulerian_jacobian2 (const DShape &d2psids, DenseMatrix< double > &jacobian2) const
 
virtual void assemble_eulerian_base_vectors (const DShape &dpsids, DenseMatrix< double > &interpolated_G) const
 
template<unsigned DIM>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double invert_jacobian_mapping (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping_diagonal (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual void dJ_eulerian_dnodal_coordinates (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<unsigned DIM>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
virtual void d_dshape_eulerian_dnodal_coordinates (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<unsigned DIM>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
virtual void transform_derivatives (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
void transform_derivatives_diagonal (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
virtual void transform_second_derivatives (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
virtual void fill_in_jacobian_from_nodal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_nodal_by_fd (DenseMatrix< double > &jacobian)
 
virtual void update_before_nodal_fd ()
 
virtual void reset_after_nodal_fd ()
 
virtual void update_in_nodal_fd (const unsigned &i)
 
virtual void reset_in_nodal_fd (const unsigned &i)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Zero-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 One-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Two-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
- Protected Member Functions inherited from oomph::GeneralisedElement
unsigned add_internal_data (Data *const &data_pt, const bool &fd=true)
 
bool internal_data_fd (const unsigned &i) const
 
void exclude_internal_data_fd (const unsigned &i)
 
void include_internal_data_fd (const unsigned &i)
 
void clear_global_eqn_numbers ()
 
void add_global_eqn_numbers (std::deque< unsigned long > const &global_eqn_numbers, std::deque< double * > const &global_dof_pt)
 
virtual void assign_internal_and_external_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void assign_additional_local_eqn_numbers ()
 
int internal_local_eqn (const unsigned &i, const unsigned &j) const
 
int external_local_eqn (const unsigned &i, const unsigned &j)
 
void fill_in_jacobian_from_internal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_internal_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
virtual void update_before_internal_fd ()
 
virtual void reset_after_internal_fd ()
 
virtual void update_in_internal_fd (const unsigned &i)
 
virtual void reset_in_internal_fd (const unsigned &i)
 
virtual void update_before_external_fd ()
 
virtual void reset_after_external_fd ()
 
virtual void update_in_external_fd (const unsigned &i)
 
virtual void reset_in_external_fd (const unsigned &i)
 
virtual void fill_in_contribution_to_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void fill_in_contribution_to_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void fill_in_contribution_to_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void fill_in_contribution_to_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 

Protected Attributes

double * Eta_pt
 Pointer to FvK parameter. More...
 
AxisymFoepplvonKarmanPressureFctPt Pressure_fct_pt
 Pointer to pressure function: More...
 
double * Nu_pt
 Pointer to Poisson's ratio. More...
 
bool Linear_bending_model
 
AxisymFoepplvonKarmanPressureFctPt Airy_forcing_fct_pt
 Pointer to Airy forcing function. More...
 
- Protected Attributes inherited from oomph::FiniteElement
MacroElement * Macro_elem_pt
 Pointer to the element's macro element (NULL by default) More...
 
- Protected Attributes inherited from oomph::GeomObject
unsigned NLagrangian
 Number of Lagrangian (intrinsic) coordinates. More...
 
unsigned Ndim
 Number of Eulerian coordinates. More...
 
TimeStepper * Geom_object_time_stepper_pt
 

Static Private Attributes

static double Default_Physical_Constant_Value = 0.0
 Default value for physical constants. More...
 

Additional Inherited Members

- Static Public Attributes inherited from oomph::FiniteElement
static double Tolerance_for_singular_jacobian = 1.0e-16
 Tolerance below which the jacobian is considered singular. More...
 
static bool Accept_negative_jacobian = false
 
static bool Suppress_output_while_checking_for_inverted_elements
 
- Static Public Attributes inherited from oomph::GeneralisedElement
static bool Suppress_warning_about_repeated_internal_data
 
static bool Suppress_warning_about_repeated_external_data = true
 
static double Default_fd_jacobian_step = 1.0e-8
 
- Static Protected Attributes inherited from oomph::FiniteElement
static const unsigned Default_Initial_Nvalue = 0
 Default value for the number of values at a node. More...
 
static const double Node_location_tolerance = 1.0e-14
 
static const unsigned N2deriv [] = {0, 1, 3, 6}
 
- Static Protected Attributes inherited from oomph::GeneralisedElement
static DenseMatrix< double > Dummy_matrix
 
static std::deque< double * > Dof_pt_deque
 

Detailed Description

A class for all isoparametric elements that solve the axisYm Foeppl von Karman equations in a displacement based formulation.

This contains the generic maths. Shape functions, geometric mapping etc. must get implemented in derived class.

A class for all isoparametric elements that solve the axisum Foeppl von Karman equations.

This contains the generic maths. Shape functions, geometric mapping etc. must get implemented in derived class.

Member Typedef Documentation

◆ AxisymFoepplvonKarmanPressureFctPt [1/2]

typedef void(* oomph::AxisymFoepplvonKarmanEquations::AxisymFoepplvonKarmanPressureFctPt) (const double &r, double &f)

Function pointer to pressure function fct(r,f(r)) – r is a Vector!

◆ AxisymFoepplvonKarmanPressureFctPt [2/2]

typedef void(* oomph::AxisymFoepplvonKarmanEquations::AxisymFoepplvonKarmanPressureFctPt) (const double &r, double &f)

Function pointer to pressure function fct(r,f(r)) – r is a Vector!

Constructor & Destructor Documentation

◆ AxisymFoepplvonKarmanEquations() [1/4]

oomph::AxisymFoepplvonKarmanEquations::AxisymFoepplvonKarmanEquations ( )
inline

Constructor (must initialise the Pressure_fct_pt). Also set physical parameters to their default values.

60 : Pressure_fct_pt(0), Nu_pt(0) {}
oomph::AxisymFoepplvonKarmanEquations::Nu_pt
double * Nu_pt
Pointer to Poisson's ratio.
Definition: axisym_displ_based_fvk_elements.h:378
oomph::AxisymFoepplvonKarmanEquations::Pressure_fct_pt
AxisymFoepplvonKarmanPressureFctPt Pressure_fct_pt
Pointer to pressure function:
Definition: axisym_displ_based_fvk_elements.h:375

◆ AxisymFoepplvonKarmanEquations() [2/4]

oomph::AxisymFoepplvonKarmanEquations::AxisymFoepplvonKarmanEquations ( const AxisymFoepplvonKarmanEquations &  dummy)
delete

Broken copy constructor.

◆ AxisymFoepplvonKarmanEquations() [3/4]

oomph::AxisymFoepplvonKarmanEquations::AxisymFoepplvonKarmanEquations ( )
inline

Constructor (must initialise the Pressure_fct_pt and Airy_forcing_fct_pt to null). Also set physical parameters to their default values.

62  : Pressure_fct_pt(0), Airy_forcing_fct_pt(0)
63  {
64  // Set all the physical constants to the default value (zero)
65  Eta_pt = &Default_Physical_Constant_Value;
66  Linear_bending_model = false;
67  }
oomph::AxisymFoepplvonKarmanEquations::Default_Physical_Constant_Value
static double Default_Physical_Constant_Value
Default value for physical constants.
Definition: axisym_displ_based_fvk_elements.h:386
oomph::AxisymFoepplvonKarmanEquations::Airy_forcing_fct_pt
AxisymFoepplvonKarmanPressureFctPt Airy_forcing_fct_pt
Pointer to Airy forcing function.
Definition: axisym_fvk_elements.h:365
oomph::AxisymFoepplvonKarmanEquations::Eta_pt
double * Eta_pt
Pointer to FvK parameter.
Definition: axisym_displ_based_fvk_elements.h:372
oomph::AxisymFoepplvonKarmanEquations::Linear_bending_model
bool Linear_bending_model
Definition: axisym_displ_based_fvk_elements.h:382

References Default_Physical_Constant_Value, Eta_pt, and Linear_bending_model.

◆ AxisymFoepplvonKarmanEquations() [4/4]

oomph::AxisymFoepplvonKarmanEquations::AxisymFoepplvonKarmanEquations ( const AxisymFoepplvonKarmanEquations &  dummy)
delete

Broken copy constructor.

Member Function Documentation

◆ airy_forcing_fct_pt() [1/2]

AxisymFoepplvonKarmanPressureFctPt& oomph::AxisymFoepplvonKarmanEquations::airy_forcing_fct_pt ( )
inline

Access function: Pointer to Airy forcing function.

182  {
183  return Airy_forcing_fct_pt;
184  }

References Airy_forcing_fct_pt.

◆ airy_forcing_fct_pt() [2/2]

AxisymFoepplvonKarmanPressureFctPt oomph::AxisymFoepplvonKarmanEquations::airy_forcing_fct_pt ( ) const
inline

Access function: Pointer to Airy forcing function. Const version.

188  {
189  return Airy_forcing_fct_pt;
190  }

References Airy_forcing_fct_pt.

◆ compute_error() [1/4]

void oomph::AxisymFoepplvonKarmanEquations::compute_error ( std::ostream &  outfile,
FiniteElement::SteadyExactSolutionFctPt  exact_soln_pt,
double &  error,
double &  norm 
)
virtual

Get error against and norm of exact solution.

Validate against exact solution

Solution is provided via function pointer. Plot error at a given number of plot points.

Reimplemented from oomph::FiniteElement.

438  {
439  // Initialise
440  error = 0.0;
441  norm = 0.0;
442 
443  // Vector of local coordinates
444  Vector<double> s(1);
445 
446  // Vector for coordintes
447  Vector<double> r(1);
448 
449  // Find out how many nodes there are in the element
450  unsigned n_node = nnode();
451 
452  Shape psi(n_node);
453 
454  // Set the value of n_intpt
455  unsigned n_intpt = integral_pt()->nweight();
456 
457  // Tecplot
458  outfile << "ZONE" << std::endl;
459 
460  // Exact solution Vector (here a scalar)
461  // Vector<double> exact_soln(1);
462  Vector<double> exact_soln(1);
463 
464  // Loop over the integration points
465  for (unsigned ipt = 0; ipt < n_intpt; ipt++)
466  {
467  // Assign values of s
468  s[0] = integral_pt()->knot(ipt, 0);
469 
470  // Get the integral weight
471  double w = integral_pt()->weight(ipt);
472 
473  // Get jacobian of mapping
474  double J = J_eulerian(s);
475 
476  // Premultiply the weights and the Jacobian
477  double W = w * J;
478 
479  // Get r position as Vector
480  interpolated_x(s, r);
481 
482  // Get FE function value
483  double w_fe = interpolated_w_fvk(s);
484 
485  // Get exact solution at this point
486  (*exact_soln_pt)(r, exact_soln);
487 
488  // Output r,error
489  outfile << r[0] << " ";
490  outfile << exact_soln[0] << " " << exact_soln[0] - w_fe << std::endl;
491 
492  // Add to error and norm
493  norm += exact_soln[0] * exact_soln[0] * W;
494  error += (exact_soln[0] - w_fe) * (exact_soln[0] - w_fe) * W;
495  }
496 
497  {
498  // Initialise
499  error = 0.0;
500  norm = 0.0;
501 
502  // Vector of local coordinates
503  Vector<double> s(1);
504 
505  // Vector for coordintes
506  Vector<double> r(1);
507 
508  // Find out how many nodes there are in the element
509  unsigned n_node = nnode();
510 
511  Shape psi(n_node);
512 
513  // Set the value of n_intpt
514  unsigned n_intpt = integral_pt()->nweight();
515 
516  // Tecplot
517  outfile << "ZONE" << std::endl;
518 
519  // Exact solution Vector (here a scalar)
520  // Vector<double> exact_soln(1);
521  Vector<double> exact_soln(1);
522 
523  // Loop over the integration points
524  for (unsigned ipt = 0; ipt < n_intpt; ipt++)
525  {
526  // Assign values of s
527  s[0] = integral_pt()->knot(ipt, 0);
528 
529  // Get the integral weight
530  double w = integral_pt()->weight(ipt);
531 
532  // Get jacobian of mapping
533  double J = J_eulerian(s);
534 
535  // Premultiply the weights and the Jacobian
536  double W = w * J;
537 
538  // Get r position as Vector
539  interpolated_x(s, r);
540 
541  // Get FE function value
542  double w_fe = interpolated_w_fvk(s);
543 
544  // Get exact solution at this point
545  (*exact_soln_pt)(r, exact_soln);
546 
547  // Output r error
548  outfile << r[0] << " ";
549  outfile << exact_soln[0] << " " << exact_soln[0] - w_fe << std::endl;
550 
551  // Add to error and norm
552  norm += exact_soln[0] * exact_soln[0] * W;
553  error += (exact_soln[0] - w_fe) * (exact_soln[0] - w_fe) * W;
554  }
555  }
556  }
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
w
RowVector3d w
Definition: Matrix_resize_int.cpp:3
oomph::AxisymFoepplvonKarmanEquations::interpolated_w_fvk
double interpolated_w_fvk(const Vector< double > &s) const
Return FE representation of transverse displacement.
Definition: axisym_displ_based_fvk_elements.h:249
oomph::FiniteElement::interpolated_x
virtual double interpolated_x(const Vector< double > &s, const unsigned &i) const
Return FE interpolated coordinate x[i] at local coordinate s.
Definition: elements.cc:3962
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
oomph::FiniteElement::integral_pt
Integral *const & integral_pt() const
Return the pointer to the integration scheme (const version)
Definition: elements.h:1963
oomph::FiniteElement::J_eulerian
virtual double J_eulerian(const Vector< double > &s) const
Definition: elements.cc:4103
oomph::Integral::knot
virtual double knot(const unsigned &i, const unsigned &j) const =0
Return local coordinate s[j] of i-th integration point.
oomph::Integral::nweight
virtual unsigned nweight() const =0
Return the number of integration points of the scheme.
oomph::Integral::weight
virtual double weight(const unsigned &i) const =0
Return weight of i-th integration point.
s
RealScalar s
Definition: level1_cplx_impl.h:130
ProblemParameters::exact_soln
void exact_soln(const double &time, const Vector< double > &x, Vector< double > &soln)
Definition: unstructured_two_d_curved.cc:301
UniformPSDSelfTest.r
r
Definition: UniformPSDSelfTest.py:20
calibrate.error
int error
Definition: calibrate.py:297
oomph::QuadTreeNames::W
@ W
Definition: quadtree.h:63

References calibrate::error, ProblemParameters::exact_soln(), oomph::FiniteElement::integral_pt(), interpolated_w_fvk(), oomph::FiniteElement::interpolated_x(), J, oomph::FiniteElement::J_eulerian(), oomph::Integral::knot(), oomph::FiniteElement::nnode(), oomph::Integral::nweight(), UniformPSDSelfTest::r, s, w, oomph::QuadTreeNames::W, and oomph::Integral::weight().

◆ compute_error() [2/4]

void oomph::AxisymFoepplvonKarmanEquations::compute_error ( std::ostream &  outfile,
FiniteElement::SteadyExactSolutionFctPt  exact_soln_pt,
double &  error,
double &  norm 
)
virtual

Get error against and norm of exact solution.

Reimplemented from oomph::FiniteElement.

◆ compute_error() [3/4]

void oomph::AxisymFoepplvonKarmanEquations::compute_error ( std::ostream &  outfile,
FiniteElement::UnsteadyExactSolutionFctPt  exact_soln_pt,
const double &  time,
double &  error,
double &  norm 
)
inlinevirtual

Dummy, time dependent error checker.

Reimplemented from oomph::FiniteElement.

172  {
173  throw OomphLibError(
174  "There is no time-dependent compute_error() for Foeppl von Karman"
175  "elements",
176  OOMPH_CURRENT_FUNCTION,
177  OOMPH_EXCEPTION_LOCATION);
178  }
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
OOMPH_CURRENT_FUNCTION
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86

References OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ compute_error() [4/4]

void oomph::AxisymFoepplvonKarmanEquations::compute_error ( std::ostream &  outfile,
FiniteElement::UnsteadyExactSolutionFctPt  exact_soln_pt,
const double &  time,
double &  error,
double &  norm 
)
inlinevirtual

Dummy, time dependent error checker.

Reimplemented from oomph::FiniteElement.

160  {
161  throw OomphLibError(
162  "There is no time-dependent compute_error() for Foeppl von Karman"
163  "elements",
164  OOMPH_CURRENT_FUNCTION,
165  OOMPH_EXCEPTION_LOCATION);
166  }

References OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ dshape_and_dtest_eulerian_at_knot_axisym_fvk() [1/2]

virtual double oomph::AxisymFoepplvonKarmanEquations::dshape_and_dtest_eulerian_at_knot_axisym_fvk ( const unsigned &  ipt,
Shape &  psi,
DShape &  dpsidr,
Shape &  test,
DShape &  dtestdr 
) const
protectedpure virtual

Shape/test functions and derivs w.r.t. to global coords at integration point ipt; return Jacobian of mapping

Implemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

Referenced by fill_in_contribution_to_residuals().

◆ dshape_and_dtest_eulerian_at_knot_axisym_fvk() [2/2]

virtual double oomph::AxisymFoepplvonKarmanEquations::dshape_and_dtest_eulerian_at_knot_axisym_fvk ( const unsigned &  ipt,
Shape &  psi,
DShape &  dpsidr,
Shape &  test,
DShape &  dtestdr 
) const
protectedpure virtual

Shape/test functions and derivs w.r.t. to global coords at integration point ipt; return Jacobian of mapping

Implemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

◆ dshape_and_dtest_eulerian_axisym_fvk() [1/2]

virtual double oomph::AxisymFoepplvonKarmanEquations::dshape_and_dtest_eulerian_axisym_fvk ( const Vector< double > &  s,
Shape &  psi,
DShape &  dpsidr,
Shape &  test,
DShape &  dtestdr 
) const
protectedpure virtual

Shape/test functions and derivs w.r.t. to global coords at local coord. s; return Jacobian of mapping

Implemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

◆ dshape_and_dtest_eulerian_axisym_fvk() [2/2]

virtual double oomph::AxisymFoepplvonKarmanEquations::dshape_and_dtest_eulerian_axisym_fvk ( const Vector< double > &  s,
Shape &  psi,
DShape &  dpsidr,
Shape &  test,
DShape &  dtestdr 
) const
protectedpure virtual

Shape/test functions and derivs w.r.t. to global coords at local coord. s; return Jacobian of mapping

Implemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

◆ eta() [1/2]

const double& oomph::AxisymFoepplvonKarmanEquations::eta ( ) const
inline

FvK parameter.

92  {
93  return *Eta_pt;
94  }

References Eta_pt.

Referenced by fill_in_contribution_to_residuals().

◆ eta() [2/2]

const double& oomph::AxisymFoepplvonKarmanEquations::eta ( ) const
inline

FvK parameter.

78  {
79  return *Eta_pt;
80  }

References Eta_pt.

◆ eta_pt() [1/2]

double*& oomph::AxisymFoepplvonKarmanEquations::eta_pt ( )
inline

Pointer to FvK parameter.

98  {
99  return Eta_pt;
100  }

References Eta_pt.

◆ eta_pt() [2/2]

double*& oomph::AxisymFoepplvonKarmanEquations::eta_pt ( )
inline

Pointer to FvK parameter.

84  {
85  return Eta_pt;
86  }

References Eta_pt.

◆ fill_in_contribution_to_residuals() [1/2]

void oomph::AxisymFoepplvonKarmanEquations::fill_in_contribution_to_residuals ( Vector< double > &  residuals)
virtual

Fill in the residuals with this element's contribution.

Compute contribution to element residual Vector

Pure version without hanging nodes

Reimplemented from oomph::GeneralisedElement.

47  {
48  // Find out how many nodes there are
49  const unsigned n_node = nnode();
50 
51  // Set up memory for the shape and test functions
52  Shape psi(n_node), test(n_node);
53  DShape dpsidr(n_node, 1), dtestdr(n_node, 1);
54 
55  // Set the value of n_intpt
56  const unsigned n_intpt = integral_pt()->nweight();
57 
58  // Indices at which the unknowns are stored
59  const unsigned w_nodal_index = nodal_index_fvk(0);
60  const unsigned laplacian_w_nodal_index = nodal_index_fvk(1);
61  const unsigned u_r_nodal_index = nodal_index_fvk(2);
62 
63  // Local copy of parameters
64  const double nu_local = nu();
65  const double eta_local = eta();
66 
67  // Integers to store the local equation numbers
68  int local_eqn = 0;
69 
70  // Loop over the integration points
71  for (unsigned ipt = 0; ipt < n_intpt; ipt++)
72  {
73  // Get the integral weight
74  double w = integral_pt()->weight(ipt);
75 
76  // Call the derivatives of the shape and test functions
77  double J = dshape_and_dtest_eulerian_at_knot_axisym_fvk(
78  ipt, psi, dpsidr, test, dtestdr);
79 
80  // Allocate and initialise to zero storage for the interpolated values
81  double interpolated_r = 0.0;
82 
83  double interpolated_w = 0.0;
84  double interpolated_laplacian_w = 0.0;
85  double interpolated_u_r = 0.0;
86 
87  double interpolated_dwdr = 0.0;
88  double interpolated_dlaplacian_wdr = 0.0;
89  double interpolated_du_rdr = 0.0;
90 
91  Vector<double> nodal_value(3, 0.0);
92 
93  // Calculate function values and derivatives:
94  //-----------------------------------------
95  // Loop over nodes
96  for (unsigned l = 0; l < n_node; l++)
97  {
98  // Get the nodal values
99  nodal_value[0] = raw_nodal_value(l, w_nodal_index);
100  nodal_value[1] = raw_nodal_value(l, laplacian_w_nodal_index);
101  nodal_value[2] = raw_nodal_value(l, u_r_nodal_index);
102 
103  // Add contributions from current node/shape function
104  interpolated_w += nodal_value[0] * psi(l);
105  interpolated_laplacian_w += nodal_value[1] * psi(l);
106  interpolated_u_r += nodal_value[2] * psi(l);
107 
108  interpolated_r += raw_nodal_position(l, 0) * psi(l);
109 
110  interpolated_dwdr += nodal_value[0] * dpsidr(l, 0);
111  interpolated_dlaplacian_wdr += nodal_value[1] * dpsidr(l, 0);
112  interpolated_du_rdr += nodal_value[2] * dpsidr(l, 0);
113 
114  } // End of loop over the nodes
115 
116  // Premultiply the weights and the Jacobian
117  double W = w * interpolated_r * J;
118 
119  // Get pressure function
120  //---------------------
121  double pressure = 0.0;
122  get_pressure_fvk(ipt, interpolated_r, pressure);
123 
124  // Determine the stresses
125  //-----------------------
126 
127  double sigma_r_r = 0.0;
128  double sigma_phi_phi = 0.0;
129 
130  if (!Linear_bending_model)
131  {
132  sigma_r_r =
133  1.0 / (1.0 - nu_local * nu_local) *
134  (interpolated_du_rdr + 0.5 * interpolated_dwdr * interpolated_dwdr +
135  nu_local * 1.0 / interpolated_r * interpolated_u_r);
136 
137  sigma_phi_phi =
138  1.0 / (1.0 - nu_local * nu_local) *
139  (1.0 / interpolated_r * interpolated_u_r +
140  nu_local * (interpolated_du_rdr +
141  0.5 * interpolated_dwdr * interpolated_dwdr));
142  }
143  else
144  {
145  sigma_r_r = 1.0 / (1.0 - nu_local * nu_local) *
146  (interpolated_du_rdr +
147  nu_local * 1.0 / interpolated_r * interpolated_u_r);
148 
149  sigma_phi_phi = 1.0 / (1.0 - nu_local * nu_local) *
150  (1.0 / interpolated_r * interpolated_u_r +
151  nu_local * interpolated_du_rdr);
152  }
153 
154 
155  // Assemble residuals and Jacobian:
156  //--------------------------------
157  // Loop over the test functions
158  for (unsigned l = 0; l < n_node; l++)
159  {
160  // Get the local equation
161  local_eqn = nodal_local_eqn(l, w_nodal_index);
162 
163  // IF it's not a boundary condition
164  if (local_eqn >= 0)
165  {
166  residuals[local_eqn] +=
167  (pressure * test(l) +
168  (dtestdr(l, 0)) * interpolated_dlaplacian_wdr) *
169  W;
170  if (!Linear_bending_model)
171  {
172  residuals[local_eqn] -=
173  eta_local * sigma_r_r * (dtestdr(l, 0)) * interpolated_dwdr * W;
174  }
175  }
176 
177  // Get the local equation
178  local_eqn = nodal_local_eqn(l, laplacian_w_nodal_index);
179 
180  // IF it's not a boundary condition
181  if (local_eqn >= 0)
182  {
183  residuals[local_eqn] += (test(l) * interpolated_laplacian_w +
184  (dtestdr(l, 0)) * interpolated_dwdr) *
185  W;
186  }
187 
188  // Get the local equation
189  local_eqn = nodal_local_eqn(l, u_r_nodal_index);
190 
191  // IF it's not a boundary condition
192  if (local_eqn >= 0)
193  {
194  residuals[local_eqn] +=
195  (sigma_r_r * dtestdr(l, 0) +
196  1.0 / interpolated_r * sigma_phi_phi * test(l)) *
197  W;
198  }
199 
200  } // End of loop over test functions
201  } // End of loop over integration points
202  }
oomph::AxisymFoepplvonKarmanEquations::dshape_and_dtest_eulerian_at_knot_axisym_fvk
virtual double dshape_and_dtest_eulerian_at_knot_axisym_fvk(const unsigned &ipt, Shape &psi, DShape &dpsidr, Shape &test, DShape &dtestdr) const =0
oomph::AxisymFoepplvonKarmanEquations::nodal_index_fvk
virtual unsigned nodal_index_fvk(const unsigned &i=0) const
Definition: axisym_displ_based_fvk_elements.h:111
oomph::AxisymFoepplvonKarmanEquations::eta
const double & eta() const
FvK parameter.
Definition: axisym_displ_based_fvk_elements.h:91
oomph::AxisymFoepplvonKarmanEquations::get_pressure_fvk
virtual void get_pressure_fvk(const unsigned &ipt, const double &r, double &pressure) const
Definition: axisym_displ_based_fvk_elements.h:196
oomph::AxisymFoepplvonKarmanEquations::nu
const double & nu() const
Poisson's ratio.
Definition: axisym_displ_based_fvk_elements.h:70
oomph::FiniteElement::nodal_value
double nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2593
oomph::FiniteElement::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432
oomph::FiniteElement::raw_nodal_value
double raw_nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2576
oomph::FiniteElement::raw_nodal_position
double raw_nodal_position(const unsigned &n, const unsigned &i) const
Definition: elements.cc:1686
Eigen::test
squared absolute sa ArrayBase::abs2 DOXCOMMA MatrixBase::cwiseAbs2 sa Eigen::abs2 DOXCOMMA Eigen::pow DOXCOMMA ArrayBase::square nearest sa Eigen::floor DOXCOMMA Eigen::ceil DOXCOMMA ArrayBase::round nearest integer not less than the given sa Eigen::floor DOXCOMMA ArrayBase::ceil not a number test
Definition: GlobalFunctions.h:109
test
Definition: indexed_view.cpp:20

References dshape_and_dtest_eulerian_at_knot_axisym_fvk(), eta(), get_pressure_fvk(), oomph::FiniteElement::integral_pt(), J, Linear_bending_model, oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::nodal_local_eqn(), oomph::FiniteElement::nodal_value(), nu(), oomph::Integral::nweight(), oomph::FiniteElement::raw_nodal_position(), oomph::FiniteElement::raw_nodal_value(), Eigen::test, w, oomph::QuadTreeNames::W, and oomph::Integral::weight().

◆ fill_in_contribution_to_residuals() [2/2]

void oomph::AxisymFoepplvonKarmanEquations::fill_in_contribution_to_residuals ( Vector< double > &  residuals)
virtual

Fill in the residuals with this element's contribution.

Reimplemented from oomph::GeneralisedElement.

◆ get_airy_forcing_fvk()

virtual void oomph::AxisymFoepplvonKarmanEquations::get_airy_forcing_fvk ( const unsigned &  ipt,
const double &  r,
double &  airy_forcing 
) const
inlinevirtual

Get Airy forcing term at (Eulerian) position r. This function is virtual to allow overloading in multi-physics problems where the strength of the pressure function might be determined by another system of equations.

219  {
220  // If no Airy forcing function has been set, return zero
221  if (Airy_forcing_fct_pt == 0)
222  {
223  airy_forcing = 0.0;
224  }
225  else
226  {
227  // Get Airy forcing strength
228  (*Airy_forcing_fct_pt)(r, airy_forcing);
229  }
230  }

References Airy_forcing_fct_pt, and UniformPSDSelfTest::r.

◆ get_gradient_of_deflection() [1/2]

void oomph::AxisymFoepplvonKarmanEquations::get_gradient_of_deflection ( const Vector< double > &  s,
Vector< double > &  gradient 
) const
inline

Get gradient of deflection: gradient[i] = dw/dr_i *‍/.

215  {
216  // Find out how many nodes there are in the element
217  const unsigned n_node = nnode();
218 
219  // Get the index at which the unknown is stored
220  const unsigned w_nodal_index = nodal_index_fvk(0);
221 
222  // Set up memory for the shape and test functions
223  Shape psi(n_node);
224  DShape dpsidr(n_node, 1);
225 
226  // Call the derivatives of the shape and test functions
227  dshape_eulerian(s, psi, dpsidr);
228 
229  // Initialise to zero
230  gradient[0] = 0.0;
231 
232  // Loop over nodes
233  for (unsigned l = 0; l < n_node; l++)
234  {
235  gradient[0] += this->nodal_value(l, w_nodal_index) * dpsidr(l, 0);
236  }
237  }
oomph::FiniteElement::dshape_eulerian
double dshape_eulerian(const Vector< double > &s, Shape &psi, DShape &dpsidx) const
Definition: elements.cc:3298

References oomph::FiniteElement::dshape_eulerian(), oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::nodal_value(), and s.

◆ get_gradient_of_deflection() [2/2]

void oomph::AxisymFoepplvonKarmanEquations::get_gradient_of_deflection ( const Vector< double > &  s,
Vector< double > &  gradient 
) const
inline

Get gradient of deflection: gradient[i] = dw/dr_i *‍/.

235  {
236  // Find out how many nodes there are in the element
237  const unsigned n_node = nnode();
238 
239  // Get the index at which the unknown is stored
240  const unsigned w_nodal_index = nodal_index_fvk(0);
241 
242  // Set up memory for the shape and test functions
243  Shape psi(n_node);
244  DShape dpsidr(n_node, 1);
245 
246  // Call the derivatives of the shape and test functions
247  dshape_eulerian(s, psi, dpsidr);
248 
249  // Initialise to zero
250  gradient[0] = 0.0;
251 
252  // Loop over nodes
253  for (unsigned l = 0; l < n_node; l++)
254  {
255  gradient[0] += this->nodal_value(l, w_nodal_index) * dpsidr(l, 0);
256  }
257  }

References oomph::FiniteElement::dshape_eulerian(), oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::nodal_value(), and s.

◆ get_pressure_fvk() [1/2]

virtual void oomph::AxisymFoepplvonKarmanEquations::get_pressure_fvk ( const unsigned &  ipt,
const double &  r,
double &  pressure 
) const
inlinevirtual

Get pressure term at (Eulerian) position r. This function is virtual to allow overloading in multi-physics problems where the strength of the pressure function might be determined by another system of equations.

Reimplemented in oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >.

199  {
200  // If no pressure function has been set, return zero
201  if (Pressure_fct_pt == 0)
202  {
203  pressure = 0.0;
204  }
205  else
206  {
207  // Get pressure strength
208  (*Pressure_fct_pt)(r, pressure);
209  }
210  }

References Pressure_fct_pt, and UniformPSDSelfTest::r.

Referenced by fill_in_contribution_to_residuals(), and oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >::get_pressure_fvk().

◆ get_pressure_fvk() [2/2]

virtual void oomph::AxisymFoepplvonKarmanEquations::get_pressure_fvk ( const unsigned &  ipt,
const double &  r,
double &  pressure 
) const
inlinevirtual

Get pressure term at (Eulerian) position r. This function is virtual to allow overloading in multi-physics problems where the strength of the pressure function might be determined by another system of equations.

Reimplemented in oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >.

199  {
200  // If no pressure function has been set, return zero
201  if (Pressure_fct_pt == 0)
202  {
203  pressure = 0.0;
204  }
205  else
206  {
207  // Get pressure strength
208  (*Pressure_fct_pt)(r, pressure);
209  }
210  }

References Pressure_fct_pt, and UniformPSDSelfTest::r.

◆ interpolated_stress() [1/2]

bool oomph::AxisymFoepplvonKarmanEquations::interpolated_stress ( const Vector< double > &  s,
double &  sigma_r_r,
double &  sigma_phi_phi 
)

Compute in-plane stresses. Return boolean to indicate success (false if attempt to evaluate stresses at zero radius)

280  {
281  // No in plane stresses if linear bending
282  if (Linear_bending_model)
283  {
284  sigma_r_r = 0.0;
285  sigma_phi_phi = 0.0;
286 
287  // Success!
288  return true;
289  }
290  else
291  {
292  // Get shape fcts and derivs
293  unsigned n_dim = this->dim();
294  unsigned n_node = this->nnode();
295  Shape psi(n_node);
296  DShape dpsi_dr(n_node, n_dim);
297 
298  // Check if we're dividing by zero
299  Vector<double> r(1);
300  this->interpolated_x(s, r);
301  if (r[0] == 0.0)
302  {
303  sigma_r_r = 0.0;
304  sigma_phi_phi = 0.0;
305  return false;
306  }
307 
308  // Get shape fcts and derivs
309  dshape_eulerian(s, psi, dpsi_dr);
310 
311  // Indices at which the unknowns are stored
312  const unsigned smooth_dphi_dr_nodal_index = nodal_index_fvk(5);
313 
314  // Allocate and initialise to zero storage for the interpolated values
315  double interpolated_r = 0;
316  double interpolated_dphi_dr = 0;
317  double interpolated_continuous_d2phi_dr2 = 0;
318 
319 
320  // Calculate function values and derivatives:
321  //-----------------------------------------
322  // Loop over nodes
323  for (unsigned l = 0; l < n_node; l++)
324  {
325  // Add contributions from current node/shape function
326  interpolated_r += raw_nodal_position(l, 0) * psi(l);
327  interpolated_dphi_dr +=
328  this->nodal_value(l, smooth_dphi_dr_nodal_index) * psi(l);
329  interpolated_continuous_d2phi_dr2 +=
330  this->nodal_value(l, smooth_dphi_dr_nodal_index) * dpsi_dr(l, 0);
331  } // End of loop over nodes
332 
333  // Compute stresses
334  sigma_r_r = interpolated_dphi_dr / interpolated_r;
335  sigma_phi_phi = interpolated_continuous_d2phi_dr2;
336 
337  // Success!
338  return true;
339 
340  } // End if
341 
342  } // End of interpolated_stress function
oomph::FiniteElement::dim
unsigned dim() const
Definition: elements.h:2611

References oomph::FiniteElement::dim(), oomph::FiniteElement::dshape_eulerian(), oomph::FiniteElement::interpolated_x(), Linear_bending_model, oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::nodal_value(), UniformPSDSelfTest::r, oomph::FiniteElement::raw_nodal_position(), and s.

◆ interpolated_stress() [2/2]

bool oomph::AxisymFoepplvonKarmanEquations::interpolated_stress ( const Vector< double > &  s,
double &  sigma_r_r,
double &  sigma_phi_phi 
) const

Compute in-plane stresses. Return boolean to indicate success (false if attempt to evaluate stresses at zero radius)

236  {
237  // No in plane stresses if linear bending
238  if (Linear_bending_model)
239  {
240  sigma_r_r = 0.0;
241  sigma_phi_phi = 0.0;
242 
243  // Success!
244  return true;
245  }
246  else
247  {
248  // Get shape fcts and derivs
249  unsigned n_dim = this->dim();
250  unsigned n_node = this->nnode();
251  Shape psi(n_node);
252  DShape dpsidr(n_node, n_dim);
253 
254  // Check if we're dividing by zero
255  Vector<double> r(1);
256  this->interpolated_x(s, r);
257  if (r[0] == 0.0)
258  {
259  sigma_r_r = 0.0;
260  sigma_phi_phi = 0.0;
261  return false;
262  }
263 
264  // Get shape fcts and derivs
265  dshape_eulerian(s, psi, dpsidr);
266 
267  // Allocate and initialise to zero storage for the interpolated values
268  double interpolated_r = 0.0;
269  double interpolated_u_r = 0.0;
270 
271  double interpolated_dwdr = 0.0;
272  double interpolated_du_rdr = 0.0;
273 
274  double nu_local = nu();
275 
276 
277  // Calculate function values and derivatives:
278  //-----------------------------------------
279  // Loop over nodes
280  for (unsigned l = 0; l < n_node; l++)
281  {
282  // Add contributions from current node/shape function
283  interpolated_r += raw_nodal_position(l, 0) * psi(l);
284  interpolated_u_r +=
285  this->raw_nodal_value(l, nodal_index_fvk(2)) * psi(l);
286  interpolated_dwdr +=
287  this->raw_nodal_value(l, nodal_index_fvk(0)) * dpsidr(l, 0);
288  interpolated_du_rdr +=
289  this->raw_nodal_value(l, nodal_index_fvk(2)) * dpsidr(l, 0);
290  } // End of loop over nodes
291 
292  // Compute the stresses:
293  //---------------------
294  sigma_r_r =
295  1.0 / (1.0 - nu_local * nu_local) *
296  (interpolated_du_rdr + 0.5 * interpolated_dwdr * interpolated_dwdr +
297  nu_local * 1.0 / interpolated_r * interpolated_u_r);
298 
299  sigma_phi_phi =
300  1.0 / (1.0 - nu_local * nu_local) *
301  (1.0 / interpolated_r * interpolated_u_r +
302  nu_local *
303  (interpolated_du_rdr + 0.5 * interpolated_dwdr * interpolated_dwdr));
304 
305  // Success!
306  return true;
307 
308  } // End if
309 
310  } // End of interpolated_stress function

References oomph::FiniteElement::dim(), oomph::FiniteElement::dshape_eulerian(), oomph::FiniteElement::interpolated_x(), Linear_bending_model, oomph::FiniteElement::nnode(), nodal_index_fvk(), nu(), UniformPSDSelfTest::r, oomph::FiniteElement::raw_nodal_position(), oomph::FiniteElement::raw_nodal_value(), and s.

Referenced by output(), and oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >::output().

◆ interpolated_u_fvk()

double oomph::AxisymFoepplvonKarmanEquations::interpolated_u_fvk ( const Vector< double > &  s) const
inline

Return FE representation of radial displacement.

278  {
279  // Find number of nodes
280  const unsigned n_node = nnode();
281 
282  // Get the index at which the radial displacement unknown is stored
283  const unsigned u_nodal_index = nodal_index_fvk(2);
284 
285  // Local shape function
286  Shape psi(n_node);
287 
288  // Find values of shape function
289  shape(s, psi);
290 
291  // Initialise value of u
292  double interpolated_u = 0.0;
293 
294  // Loop over the local nodes and sum
295  for (unsigned l = 0; l < n_node; l++)
296  {
297  interpolated_u += this->nodal_value(l, u_nodal_index) * psi[l];
298  }
299 
300  return (interpolated_u);
301  }
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0

References oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::nodal_value(), s, and oomph::FiniteElement::shape().

Referenced by output().

◆ interpolated_w_fvk() [1/2]

double oomph::AxisymFoepplvonKarmanEquations::interpolated_w_fvk ( const Vector< double > &  s) const
inline

Return FE representation of transverse displacement.

250  {
251  // Find number of nodes
252  const unsigned n_node = nnode();
253 
254  // Get the index at which the transverse displacement unknown is stored
255  const unsigned w_nodal_index = nodal_index_fvk(0);
256 
257  // Local shape function
258  Shape psi(n_node);
259 
260  // Find values of shape function
261  shape(s, psi);
262 
263  // Initialise value of u
264  double interpolated_w = 0.0;
265 
266  // Loop over the local nodes and sum
267  for (unsigned l = 0; l < n_node; l++)
268  {
269  interpolated_w += this->nodal_value(l, w_nodal_index) * psi[l];
270  }
271 
272  return (interpolated_w);
273  }

References oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::nodal_value(), s, and oomph::FiniteElement::shape().

Referenced by compute_error(), output(), and oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >::output().

◆ interpolated_w_fvk() [2/2]

double oomph::AxisymFoepplvonKarmanEquations::interpolated_w_fvk ( const Vector< double > &  s) const
inline

Return FE representation of vertical displacement, w_fvk(s) at local coordinate s

270  {
271  // Find number of nodes
272  const unsigned n_node = nnode();
273 
274  // Get the index at which the unknown is stored
275  const unsigned w_nodal_index = nodal_index_fvk(0);
276 
277  // Local shape function
278  Shape psi(n_node);
279 
280  // Find values of shape function
281  shape(s, psi);
282 
283  // Initialise value of u
284  double interpolated_w = 0.0;
285 
286  // Loop over the local nodes and sum
287  for (unsigned l = 0; l < n_node; l++)
288  {
289  interpolated_w += this->nodal_value(l, w_nodal_index) * psi[l];
290  }
291 
292  return (interpolated_w);
293  }

References oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::nodal_value(), s, and oomph::FiniteElement::shape().

◆ nodal_index_fvk() [1/2]

virtual unsigned oomph::AxisymFoepplvonKarmanEquations::nodal_index_fvk ( const unsigned &  i = 0) const
inlinevirtual

Return the index at which the i-th unknown value is stored. The default value, i, is appropriate for single-physics problems. By default, these are: 0: transverse displacement w 1: laplacian w 2: radial displacement u In derived multi-physics elements, this function should be overloaded to reflect the chosen storage scheme. Note that these equations require that the unknown is always stored at the same index at each node.

112  {
113  return i;
114  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9

References i.

Referenced by oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >::dposition_dt(), fill_in_contribution_to_residuals(), get_gradient_of_deflection(), interpolated_stress(), interpolated_u_fvk(), interpolated_w_fvk(), oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >::output_integration_points(), oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >::position(), and use_linear_bending_model().

◆ nodal_index_fvk() [2/2]

virtual unsigned oomph::AxisymFoepplvonKarmanEquations::nodal_index_fvk ( const unsigned &  i = 0) const
inlinevirtual

Return the index at which the i-th unknown value is stored. The default value, i, is appropriate for single-physics problems. By default, these are: 0: w 1: laplacian w 2: phi 3: laplacian phi 4-5: smooth first derivatives In derived multi-physics elements, this function should be overloaded to reflect the chosen storage scheme. Note that these equations require that the unknown is always stored at the same index at each node.

100  {
101  return i;
102  }

References i.

◆ nu()

const double& oomph::AxisymFoepplvonKarmanEquations::nu ( ) const
inline

Poisson's ratio.

71  {
72 #ifdef PARANOID
73  if (Nu_pt == 0)
74  {
75  std::stringstream error_stream;
76  error_stream << "Nu has not yet been set!" << std::endl;
77  throw OomphLibError(
78  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
79  }
80 #endif
81  return *Nu_pt;
82  }

References Nu_pt, OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

Referenced by fill_in_contribution_to_residuals(), and interpolated_stress().

◆ nu_pt()

double*& oomph::AxisymFoepplvonKarmanEquations::nu_pt ( )
inline

Pointer to Poisson's ratio.

86  {
87  return Nu_pt;
88  }

References Nu_pt.

◆ operator=() [1/2]

void oomph::AxisymFoepplvonKarmanEquations::operator= ( const AxisymFoepplvonKarmanEquations &  )
delete

Broken assignment operator.

◆ operator=() [2/2]

void oomph::AxisymFoepplvonKarmanEquations::operator= ( const AxisymFoepplvonKarmanEquations &  )
delete

Broken assignment operator.

◆ output() [1/8]

void oomph::AxisymFoepplvonKarmanEquations::output ( FILE *  file_pt)
inlinevirtual

C_style output with default number of plot points.

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

129  {
130  const unsigned n_plot = 5;
131  output(file_pt, n_plot);
132  }
oomph::AxisymFoepplvonKarmanEquations::output
void output(std::ostream &outfile)
Output with default number of plot points.
Definition: axisym_displ_based_fvk_elements.h:117

References output().

◆ output() [2/8]

void oomph::AxisymFoepplvonKarmanEquations::output ( FILE *  file_pt)
inlinevirtual

C_style output with default number of plot points.

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

117  {
118  const unsigned n_plot = 5;
119  output(file_pt, n_plot);
120  }

References output().

◆ output() [3/8]

void oomph::AxisymFoepplvonKarmanEquations::output ( FILE *  file_pt,
const unsigned &  nplot 
)
virtual

C-style output FE representation of soln: r,w at n_plot plot points

C-style output function: r,w nplot points

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

353  {
354  // Vector of local coordinates
355  Vector<double> s(1);
356 
357  // Tecplot header info
358  fprintf(file_pt, "%s", tecplot_zone_string(nplot).c_str());
359 
360  // Loop over plot points
361  unsigned num_plot_points = nplot_points(nplot);
362  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
363  {
364  // Get local coordinates of plot point
365  get_s_plot(iplot, nplot, s);
366 
367  fprintf(file_pt, "%g ", interpolated_x(s, 0));
368  fprintf(file_pt, "%g \n", interpolated_w_fvk(s));
369  fprintf(file_pt, "%g \n", interpolated_u_fvk(s));
370  }
371 
372  // Write tecplot footer (e.g. FE connectivity lists)
373  write_tecplot_zone_footer(file_pt, nplot);
374  }
oomph::AxisymFoepplvonKarmanEquations::interpolated_u_fvk
double interpolated_u_fvk(const Vector< double > &s) const
Return FE representation of radial displacement.
Definition: axisym_displ_based_fvk_elements.h:277
oomph::FiniteElement::tecplot_zone_string
virtual std::string tecplot_zone_string(const unsigned &nplot) const
Definition: elements.h:3161
oomph::FiniteElement::get_s_plot
virtual void get_s_plot(const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
Definition: elements.h:3148
oomph::FiniteElement::nplot_points
virtual unsigned nplot_points(const unsigned &nplot) const
Definition: elements.h:3186
oomph::FiniteElement::write_tecplot_zone_footer
virtual void write_tecplot_zone_footer(std::ostream &outfile, const unsigned &nplot) const
Definition: elements.h:3174

References oomph::FiniteElement::get_s_plot(), interpolated_u_fvk(), interpolated_w_fvk(), oomph::FiniteElement::interpolated_x(), oomph::FiniteElement::nplot_points(), s, oomph::FiniteElement::tecplot_zone_string(), and oomph::FiniteElement::write_tecplot_zone_footer().

◆ output() [4/8]

void oomph::AxisymFoepplvonKarmanEquations::output ( FILE *  file_pt,
const unsigned &  n_plot 
)
virtual

C-style output FE representation of soln: r,w at n_plot plot points

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

◆ output() [5/8]

void oomph::AxisymFoepplvonKarmanEquations::output ( std::ostream &  outfile)
inlinevirtual

Output with default number of plot points.

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

118  {
119  const unsigned n_plot = 5;
120  output(outfile, n_plot);
121  }

Referenced by output(), and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >::output().

◆ output() [6/8]

void oomph::AxisymFoepplvonKarmanEquations::output ( std::ostream &  outfile)
inlinevirtual

Output with default number of plot points.

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

106  {
107  const unsigned n_plot = 5;
108  output(outfile, n_plot);
109  }

References output().

◆ output() [7/8]

void oomph::AxisymFoepplvonKarmanEquations::output ( std::ostream &  outfile,
const unsigned &  nplot 
)
virtual

Output FE representation of soln: r,w,sigma_r_r,sigma_phi_phi at n_plot plot points

Output function: r, w, sigma_r_r, sigma_phi_phi nplot points

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >, oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

319  {
320  // Vector of local coordinates
321  Vector<double> s(1);
322 
323  // Tecplot header info
324  outfile << "ZONE\n";
325 
326  // Loop over plot points
327  unsigned num_plot_points = nplot_points(nplot);
328  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
329  {
330  // Get local coordinates of plot point
331  get_s_plot(iplot, nplot, s);
332 
333  // Get stress
334  double sigma_r_r = 0.0;
335  double sigma_phi_phi = 0.0;
336  bool success = interpolated_stress(s, sigma_r_r, sigma_phi_phi);
337  if (success)
338  {
339  outfile << interpolated_x(s, 0) << " " << interpolated_w_fvk(s) << " "
340  << interpolated_u_fvk(s) << " " << sigma_r_r << " "
341  << sigma_phi_phi << std::endl;
342  }
343  }
344  }
oomph::AxisymFoepplvonKarmanEquations::interpolated_stress
bool interpolated_stress(const Vector< double > &s, double &sigma_r_r, double &sigma_phi_phi) const
Definition: axisym_displ_based_fvk_elements.cc:234

References oomph::FiniteElement::get_s_plot(), interpolated_stress(), interpolated_u_fvk(), interpolated_w_fvk(), oomph::FiniteElement::interpolated_x(), oomph::FiniteElement::nplot_points(), and s.

◆ output() [8/8]

void oomph::AxisymFoepplvonKarmanEquations::output ( std::ostream &  outfile,
const unsigned &  n_plot 
)
virtual

Output FE representation of soln: r,w,sigma_r_r,sigma_phi_phi at n_plot plot points

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::FSIAxisymFoepplvonKarmanElement< NNODE_1D, FLUID_ELEMENT >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

◆ output_fct() [1/4]

virtual void oomph::AxisymFoepplvonKarmanEquations::output_fct ( std::ostream &  outfile,
const unsigned &  n_plot,
const double &  time,
FiniteElement::UnsteadyExactSolutionFctPt  exact_soln_pt 
)
inlinevirtual

Output exact soln: r,w_exact at n_plot plot points (dummy time-dependent version to keep intel compiler happy)

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

151  {
152  throw OomphLibError(
153  "There is no time-dependent output_fct() for Foeppl von Karman"
154  "elements ",
155  OOMPH_CURRENT_FUNCTION,
156  OOMPH_EXCEPTION_LOCATION);
157  }

References OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ output_fct() [2/4]

virtual void oomph::AxisymFoepplvonKarmanEquations::output_fct ( std::ostream &  outfile,
const unsigned &  n_plot,
const double &  time,
FiniteElement::UnsteadyExactSolutionFctPt  exact_soln_pt 
)
inlinevirtual

Output exact soln: r,w_exact at n_plot plot points (dummy time-dependent version to keep intel compiler happy)

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

139  {
140  throw OomphLibError(
141  "There is no time-dependent output_fct() for Foeppl von Karman"
142  "elements ",
143  OOMPH_CURRENT_FUNCTION,
144  OOMPH_EXCEPTION_LOCATION);
145  }

References OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ output_fct() [3/4]

void oomph::AxisymFoepplvonKarmanEquations::output_fct ( std::ostream &  outfile,
const unsigned &  nplot,
FiniteElement::SteadyExactSolutionFctPt  exact_soln_pt 
)
virtual

Output exact soln: r,w_exact at n_plot plot points.

Output exact solution

Solution is provided via function pointer. Plot at a given number of plot points.

r,w_exact

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >, and oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

389  {
390  // Vector of local coordinates
391  Vector<double> s(1);
392 
393  // Vector for coordinates
394  Vector<double> r(1);
395 
396  // Tecplot header info
397  outfile << tecplot_zone_string(nplot);
398 
399  // Exact solution Vector (here a scalar)
400  // Vector<double> exact_soln(1);
401  Vector<double> exact_soln(1);
402 
403  // Loop over plot points
404  unsigned num_plot_points = nplot_points(nplot);
405  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
406  {
407  // Get local coordinates of plot point
408  get_s_plot(iplot, nplot, s);
409 
410  // Get r position as Vector
411  interpolated_x(s, r);
412 
413  // Get exact solution at this point
414  (*exact_soln_pt)(r, exact_soln);
415 
416  // Output r,w_exact
417  outfile << r[0] << " ";
418  outfile << exact_soln[0] << std::endl;
419  }
420 
421  // Write tecplot footer (e.g. FE connectivity lists)
422  write_tecplot_zone_footer(outfile, nplot);
423  }

References ProblemParameters::exact_soln(), oomph::FiniteElement::get_s_plot(), oomph::FiniteElement::interpolated_x(), oomph::FiniteElement::nplot_points(), UniformPSDSelfTest::r, s, oomph::FiniteElement::tecplot_zone_string(), and oomph::FiniteElement::write_tecplot_zone_footer().

Referenced by oomph::AxisymFoepplvonKarmanElement< NNODE_1D >::output_fct().

◆ output_fct() [4/4]

void oomph::AxisymFoepplvonKarmanEquations::output_fct ( std::ostream &  outfile,
const unsigned &  n_plot,
FiniteElement::SteadyExactSolutionFctPt  exact_soln_pt 
)
virtual

Output exact soln: r,w_exact at n_plot plot points.

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::AxisymFoepplvonKarmanElement< NNODE_1D >.

◆ pressure_fct_pt() [1/4]

AxisymFoepplvonKarmanPressureFctPt& oomph::AxisymFoepplvonKarmanEquations::pressure_fct_pt ( )
inline

Access function: Pointer to pressure function.

182  {
183  return Pressure_fct_pt;
184  }

References Pressure_fct_pt.

◆ pressure_fct_pt() [2/4]

AxisymFoepplvonKarmanPressureFctPt& oomph::AxisymFoepplvonKarmanEquations::pressure_fct_pt ( )
inline

Access function: Pointer to pressure function.

170  {
171  return Pressure_fct_pt;
172  }

References Pressure_fct_pt.

◆ pressure_fct_pt() [3/4]

AxisymFoepplvonKarmanPressureFctPt oomph::AxisymFoepplvonKarmanEquations::pressure_fct_pt ( ) const
inline

Access function: Pointer to pressure function. Const version.

188  {
189  return Pressure_fct_pt;
190  }

References Pressure_fct_pt.

◆ pressure_fct_pt() [4/4]

AxisymFoepplvonKarmanPressureFctPt oomph::AxisymFoepplvonKarmanEquations::pressure_fct_pt ( ) const
inline

Access function: Pointer to pressure function. Const version.

176  {
177  return Pressure_fct_pt;
178  }

References Pressure_fct_pt.

◆ self_test() [1/2]

unsigned oomph::AxisymFoepplvonKarmanEquations::self_test ( )
virtual

Self-test: Return 0 for OK.

Reimplemented from oomph::FiniteElement.

209  {
210  bool passed = true;
211 
212  // Check lower-level stuff
213  if (FiniteElement::self_test() != 0)
214  {
215  passed = false;
216  }
217 
218  // Return verdict
219  if (passed)
220  {
221  return 0;
222  }
223  else
224  {
225  return 1;
226  }
227  }
oomph::FiniteElement::self_test
virtual unsigned self_test()
Definition: elements.cc:4440

References oomph::FiniteElement::self_test().

◆ self_test() [2/2]

unsigned oomph::AxisymFoepplvonKarmanEquations::self_test ( )
virtual

Self-test: Return 0 for OK.

Reimplemented from oomph::FiniteElement.

◆ use_linear_bending_model() [1/2]

void oomph::AxisymFoepplvonKarmanEquations::use_linear_bending_model ( )
inline

Sets a flag to signify that we are solving the linear, pure bending equations, and pin all the nodal values that will not be used in this case

322  {
323  // Set the boolean flag
324  Linear_bending_model = true;
325 
326  // Get the index of the first FvK nodal value
327  unsigned first_fvk_nodal_index = nodal_index_fvk();
328 
329  // Get the total number of FvK nodal values (assuming they are stored
330  // contiguously) at node 0 (it's the same at all nodes anyway)
331  unsigned total_fvk_nodal_indices = 3;
332 
333  // Get the number of nodes in this element
334  unsigned n_node = nnode();
335 
336  // Loop over the appropriate nodal indices
337  for (unsigned index = first_fvk_nodal_index + 2;
338  index < first_fvk_nodal_index + total_fvk_nodal_indices;
339  index++)
340  {
341  // Loop over the nodes in the element
342  for (unsigned inod = 0; inod < n_node; inod++)
343  {
344  // Pin the nodal value at the current index
345  node_pt(inod)->pin(index);
346  }
347  }
348  }
oomph::Data::pin
void pin(const unsigned &i)
Pin the i-th stored variable.
Definition: nodes.h:385
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175

References Linear_bending_model, oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::node_pt(), and oomph::Data::pin().

◆ use_linear_bending_model() [2/2]

void oomph::AxisymFoepplvonKarmanEquations::use_linear_bending_model ( )
inline

Sets a flag to signify that we are solving the linear, pure bending equations, and pin all the nodal values that will not be used in this case

309  {
310  // Set the boolean flag
311  Linear_bending_model = true;
312 
313  // Get the index of the first FvK nodal value
314  unsigned first_fvk_nodal_index = nodal_index_fvk();
315 
316  // Get the total number of FvK nodal values (assuming they are stored
317  // contiguously) at node 0 (it's the same at all nodes anyway)
318  unsigned total_fvk_nodal_indices = 6;
319 
320  // Get the number of nodes in this element
321  unsigned n_node = nnode();
322 
323  // Loop over the appropriate nodal indices
324  for (unsigned index = first_fvk_nodal_index + 2;
325  index < first_fvk_nodal_index + total_fvk_nodal_indices;
326  index++)
327  {
328  // Loop over the nodes in the element
329  for (unsigned inod = 0; inod < n_node; inod++)
330  {
331  // Pin the nodal value at the current index
332  node_pt(inod)->pin(index);
333  }
334  }
335  }

References Linear_bending_model, oomph::FiniteElement::nnode(), nodal_index_fvk(), oomph::FiniteElement::node_pt(), and oomph::Data::pin().

Member Data Documentation

◆ Airy_forcing_fct_pt

AxisymFoepplvonKarmanPressureFctPt oomph::AxisymFoepplvonKarmanEquations::Airy_forcing_fct_pt
protected

Pointer to Airy forcing function.

Referenced by airy_forcing_fct_pt(), and get_airy_forcing_fvk().

◆ Default_Physical_Constant_Value

static double oomph::AxisymFoepplvonKarmanEquations::Default_Physical_Constant_Value = 0.0
staticprivate

Default value for physical constants.

Default value physical constants.

Referenced by AxisymFoepplvonKarmanEquations().

◆ Eta_pt

double * oomph::AxisymFoepplvonKarmanEquations::Eta_pt
protected

Pointer to FvK parameter.

Referenced by AxisymFoepplvonKarmanEquations(), eta(), and eta_pt().

◆ Linear_bending_model

bool oomph::AxisymFoepplvonKarmanEquations::Linear_bending_model
protected

Flag which stores whether we are using a linear, pure bending model instead of the full non-linear Foeppl-von Karman

Referenced by AxisymFoepplvonKarmanEquations(), fill_in_contribution_to_residuals(), interpolated_stress(), and use_linear_bending_model().

◆ Nu_pt

double* oomph::AxisymFoepplvonKarmanEquations::Nu_pt
protected

Pointer to Poisson's ratio.

Referenced by nu(), and nu_pt().

◆ Pressure_fct_pt

AxisymFoepplvonKarmanPressureFctPt oomph::AxisymFoepplvonKarmanEquations::Pressure_fct_pt
protected

Pointer to pressure function:

Referenced by get_pressure_fvk(), and pressure_fct_pt().

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