Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
oomph::SurfaceContactElementBase< ELEMENT > Class Template Referenceabstract

#include <contact_elements.h>

+ Inheritance diagram for oomph::SurfaceContactElementBase< ELEMENT >:

Public Member Functions

 SurfaceContactElementBase (FiniteElement *const &element_pt, const int &face_index, const unsigned &id=0, const bool &called_from_refineable_constructor=false)
 value of the index and its limit More...
 
 SurfaceContactElementBase ()
 Default constructor. More...
 
void enable_stick ()
 
double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double s_min () const
 Min value of local coordinate. More...
 
double s_max () const
 Max. value of local coordinate. More...
 
void disable_stick ()
 
bool is_stick_enabled ()
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Return the residuals. More...
 
Penetrator * penetrator_pt () const
 Pointer to penetrator. More...
 
void set_penetrator_pt (Penetrator *penetrator_pt)
 Set pointer to penetrator. More...
 
void output (FILE *file_pt)
 C_style output function. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 C-style output function. More...
 
void output (std::ostream &outfile)
 Output function. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 Output function. More...
 
void shape_p (const Vector< double > &s, Shape &psi) const
 Shape fct for lagrange multiplier. More...
 
void shape_i (const Vector< double > &s, Shape &psi) const
 
- Public Member Functions inherited from oomph::FaceElement
 FaceElement ()
 Constructor: Initialise all appropriate member data. More...
 
virtual ~FaceElement ()
 Empty virtual destructor. More...
 
 FaceElement (const FaceElement &)=delete
 Broken copy constructor. More...
 
const unsigned & boundary_number_in_bulk_mesh () const
 Broken assignment operator. More...
 
void set_boundary_number_in_bulk_mesh (const unsigned &b)
 Set function for the boundary number in bulk mesh. More...
 
double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double J_eulerian (const Vector< double > &s) const
 
double J_eulerian_at_knot (const unsigned &ipt) const
 
void check_J_eulerian_at_knots (bool &passed) const
 
double interpolated_x (const Vector< double > &s, const unsigned &i) const
 
double interpolated_x (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
void interpolated_x (const Vector< double > &s, Vector< double > &x) const
 
void interpolated_x (const unsigned &t, const Vector< double > &s, Vector< double > &x) const
 
double interpolated_dxdt (const Vector< double > &s, const unsigned &i, const unsigned &t)
 
void interpolated_dxdt (const Vector< double > &s, const unsigned &t, Vector< double > &dxdt)
 
int & normal_sign ()
 
int normal_sign () const
 
int & face_index ()
 
int face_index () const
 
const Vector< double > * tangent_direction_pt () const
 Public access function for the tangent direction pointer. More...
 
void set_tangent_direction (Vector< double > *tangent_direction_pt)
 Set the tangent direction vector. More...
 
void turn_on_warning_for_discontinuous_tangent ()
 
void turn_off_warning_for_discontinuous_tangent ()
 
void continuous_tangent_and_outer_unit_normal (const Vector< double > &s, Vector< Vector< double >> &tang_vec, Vector< double > &unit_normal) const
 
void continuous_tangent_and_outer_unit_normal (const unsigned &ipt, Vector< Vector< double >> &tang_vec, Vector< double > &unit_normal) const
 
void outer_unit_normal (const Vector< double > &s, Vector< double > &unit_normal) const
 Compute outer unit normal at the specified local coordinate. More...
 
void outer_unit_normal (const unsigned &ipt, Vector< double > &unit_normal) const
 Compute outer unit normal at ipt-th integration point. More...
 
FiniteElement *& bulk_element_pt ()
 Pointer to higher-dimensional "bulk" element. More...
 
FiniteElement * bulk_element_pt () const
 Pointer to higher-dimensional "bulk" element (const version) More...
 
CoordinateMappingFctPt & face_to_bulk_coordinate_fct_pt ()
 
CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt () const
 
BulkCoordinateDerivativesFctPt & bulk_coordinate_derivatives_fct_pt ()
 
BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt () const
 
Vector< double > local_coordinate_in_bulk (const Vector< double > &s) const
 
void get_local_coordinate_in_bulk (const Vector< double > &s, Vector< double > &s_bulk) const
 
void get_ds_bulk_ds_face (const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction) const
 
unsigned & bulk_position_type (const unsigned &i)
 
const unsigned & bulk_position_type (const unsigned &i) const
 
void bulk_node_number_resize (const unsigned &i)
 Resize the storage for the bulk node numbers. More...
 
unsigned & bulk_node_number (const unsigned &n)
 
const unsigned & bulk_node_number (const unsigned &n) const
 
void bulk_position_type_resize (const unsigned &i)
 Resize the storage for bulk_position_type to i entries. More...
 
unsigned & nbulk_value (const unsigned &n)
 
unsigned nbulk_value (const unsigned &n) const
 
void nbulk_value_resize (const unsigned &i)
 
void resize_nodes (Vector< unsigned > &nadditional_data_values)
 
void output_zeta (std::ostream &outfile, const unsigned &nplot)
 Output boundary coordinate zeta. More...
 
- 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
 
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...
 
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)
 
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)
 
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, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output an exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 Output a time-dependent exact solution over the element. More...
 
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, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, 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 self_test ()
 
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
 
- Public Member Functions inherited from oomph::TemplateFreeContactElementBase
 TemplateFreeContactElementBase ()
 
virtual ~TemplateFreeContactElementBase ()
 Destructor. More...
 
virtual void resulting_contact_force (Vector< double > &contact_force)=0
 Resulting contact force. More...
 
void traction_fct (const Vector< double > &x, Vector< double > &p)
 
bool use_isoparametric_flag () const
 
bool use_collocated_penetration_flag ()
 
bool use_collocated_contact_pressure_flag ()
 
TractionFctPt & traction_fct_pt ()
 Access function: Pointer to body force function. More...
 
TractionFctPt traction_fct_pt () const
 Access function: Pointer to body force function (const version) More...
 
bool *& use_isoparametric_flag_pt ()
 Access function: Pointer to flag to use isoparametric. More...
 
bool * use_isoparametric_flag_pt () const
 
bool *& use_collocated_penetration_flag_pt ()
 Access function: Pointer to flag to use collocated penetration. More...
 
bool * use_collocated_penetration_flag_pt () const
 
bool *& use_collocated_contact_pressure_flag_pt ()
 Access function: Pointer to flag to use collocated contact pressure. More...
 
bool * use_collocated_contact_pressure_flag_pt () const
 

Protected Member Functions

double get_interpolated_lagrange_p (const Vector< double > &s)
 
virtual void fill_in_contribution_to_residuals_surface_contact (Vector< double > &residuals)=0
 
void penetration (const Vector< double > &x, const Vector< double > &n, double &d, bool &intersection) const
 Work out penetration of point. More...
 
- Protected Member Functions inherited from oomph::FaceElement
void add_additional_values (const Vector< unsigned > &nadditional_values, const unsigned &id)
 
- 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

Penetrator * Penetrator_pt
 Pointer to penetrator. More...
 
unsigned Contact_id
 
bool Enable_stick
 
Vector< int > Penetrator_eq_data_type
 
Vector< int > Penetrator_eq_data_index
 
Vector< int > Penetrator_eq_data_data_index
 
- Protected Attributes inherited from oomph::FaceElement
unsigned Boundary_number_in_bulk_mesh
 The boundary number in the bulk mesh to which this element is attached. More...
 
FiniteElement * Bulk_element_pt
 Pointer to the associated higher-dimensional "bulk" element. More...
 
Vector< unsigned > Bulk_node_number
 
Vector< unsigned > Nbulk_value
 
Vector< double > * Tangent_direction_pt
 
- 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
 
- Protected Attributes inherited from oomph::TemplateFreeContactElementBase
TractionFctPt Traction_fct_pt
 
bool * Use_isoparametric_flag_pt
 Set whether or not to use isoparametric basis function for pressure. More...
 
bool * Use_collocated_penetration_flag_pt
 Set options for basis/test functions for penetration and pressure. More...
 
bool * Use_collocated_contact_pressure_flag_pt
 Set options for basis/test functions for penetration and pressure. More...
 

Additional Inherited Members

- 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 Types inherited from oomph::TemplateFreeContactElementBase
typedef void(* TractionFctPt) (const double &t, const Vector< double > &x, Vector< double > &p)
 
- 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
 
- Protected Types inherited from oomph::TemplateFreeContactElementBase
enum  { Nodal_data , Nodal_position_data , External_data }
 Enumeration for type of penetrator data. More...
 
- 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

template<class ELEMENT>
class oomph::SurfaceContactElementBase< ELEMENT >

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// Base class for elements that impose contact boundary conditions either enforcing non-penetration (but "without stick"; default) or permament contact (in which case contact force can be positive or negative. Uses Lagrange-multiplier-like pressure to enforce contact/non-penetration. Almost certainly works only with 2D Penetrator at the moment. Specific implementation required for linear or nonlinear elasticity.

Constructor & Destructor Documentation

◆ SurfaceContactElementBase() [1/2]

template<class ELEMENT >
oomph::SurfaceContactElementBase< ELEMENT >::SurfaceContactElementBase ( FiniteElement *const &  element_pt,
const int &  face_index,
const unsigned &  id = 0,
const bool &  called_from_refineable_constructor = false 
)
inline

value of the index and its limit

Constructor, which takes a "bulk" element and the

685  :
686  FaceGeometry<ELEMENT>(), FaceElement()
687  {
688 
689  // By default we only to proper non-penetration (without "stick", i.e.
690  // we don't allow negative contact pressures)
691  Enable_stick=false;
692 
693  // Initialise pointer to penetrator
694  Penetrator_pt=0;
695 
696  //Attach the geometrical information to the element. N.B. This function
697  //also assigns nbulk_value from the required_nvalue of the bulk element
698  element_pt->build_face_element(face_index,this);
699 
700 
701 #ifdef PARANOID
702  {
703  //Check that the bulk element is not a refineable 3d element
704  if (!called_from_refineable_constructor)
705  {
706  if(element_pt->dim()==3)
707  {
708  //Is it refineable
709  RefineableElement* ref_el_pt=
710  dynamic_cast<RefineableElement*>(element_pt);
711  if(ref_el_pt!=0)
712  {
713  if (this->has_hanging_nodes())
714  {
715  throw OomphLibError(
716  "This face element will not work correctly if nodes are hanging.\nUse the refineable version instead. ",
717  OOMPH_CURRENT_FUNCTION,
718  OOMPH_EXCEPTION_LOCATION);
719  }
720  }
721  }
722  }
723  }
724 #endif
725 
726  // Store the ID of the FaceElement -- this is used to distinguish
727  // it from any others
728  Contact_id=id;
729 
730  // We need one additional value for each FaceElement node:
731  // the normal traction (Lagrange multiplier) to be
732  // exerted onto the solid
733  unsigned n_nod=nnode();
734  Vector<unsigned> n_additional_values(n_nod,1);
735 
736  // Now add storage for Lagrange multipliers and set the map containing
737  // the position of the first entry of this face element's
738  // additional values.
739  add_additional_values(n_additional_values,id);
740 
741 #ifdef PARANOID
742  // Check spatial dimension
743  if (element_pt->dim()!=2)
744  {
745  //Issue a warning
746  throw OomphLibError(
747  "This element will almost certainly not work in non-2D problems, though it should be easy enough to upgrade... Volunteers?\n",
748  OOMPH_CURRENT_FUNCTION,
749  OOMPH_EXCEPTION_LOCATION);
750  }
751 #endif
752 
753  // Always use piecewise Gauss -- it "over-integrates" for isoparametric formulations
754  // but is likely to be better for discontinuous basis fcts.
755  //Only implemented for 1,3 (1 is from surface of 2D bulk) (3 for quadratic elements)
756 #ifdef PARANOID
757  // Bit hacky... Only really works for three-noded 1D elements
758  if (n_nod!=3)
759  {
760  //Issue a warning
761  throw OomphLibError(
762  "Piecwise Gauss used here isn't appropriate for non-3-node elements.\n",
763  OOMPH_CURRENT_FUNCTION,
764  OOMPH_EXCEPTION_LOCATION);
765  }
766 #endif
767  set_integration_scheme(new PiecewiseGauss<1,3>(this->s_min(),this->s_max()));
768  }
oomph::FaceElement::face_index
int & face_index()
Definition: elements.h:4626
oomph::FaceElement::FaceElement
FaceElement()
Constructor: Initialise all appropriate member data.
Definition: elements.h:4446
oomph::FaceElement::add_additional_values
void add_additional_values(const Vector< unsigned > &nadditional_values, const unsigned &id)
Definition: elements.h:4428
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
oomph::FiniteElement::set_integration_scheme
virtual void set_integration_scheme(Integral *const &integral_pt)
Set the spatial integration scheme.
Definition: elements.cc:3210
oomph::FiniteElement::has_hanging_nodes
bool has_hanging_nodes() const
Definition: elements.h:2470
oomph::SurfaceContactElementBase::Contact_id
unsigned Contact_id
Definition: contact_elements.h:1035
oomph::SurfaceContactElementBase::Enable_stick
bool Enable_stick
Definition: contact_elements.h:1039
oomph::SurfaceContactElementBase::s_max
double s_max() const
Max. value of local coordinate.
Definition: contact_elements.h:793
oomph::SurfaceContactElementBase::Penetrator_pt
Penetrator * Penetrator_pt
Pointer to penetrator.
Definition: contact_elements.h:1030
oomph::SurfaceContactElementBase::s_min
double s_min() const
Min value of local coordinate.
Definition: contact_elements.h:788
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::FaceElement::add_additional_values(), oomph::FiniteElement::build_face_element(), oomph::SurfaceContactElementBase< ELEMENT >::Contact_id, oomph::FiniteElement::dim(), oomph::SurfaceContactElementBase< ELEMENT >::Enable_stick, oomph::FaceElement::face_index(), oomph::FiniteElement::has_hanging_nodes(), oomph::FiniteElement::nnode(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_pt, oomph::SurfaceContactElementBase< ELEMENT >::s_max(), oomph::SurfaceContactElementBase< ELEMENT >::s_min(), and oomph::FiniteElement::set_integration_scheme().

◆ SurfaceContactElementBase() [2/2]

template<class ELEMENT >
oomph::SurfaceContactElementBase< ELEMENT >::SurfaceContactElementBase ( )
inline

Default constructor.

772 {}

Member Function Documentation

◆ disable_stick()

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::disable_stick ( )
inline

Allow only proper non-penetration (without "stick", i.e. we don't allow negative contact pressures)

801  {
802  Enable_stick=false;
803  }

References oomph::SurfaceContactElementBase< ELEMENT >::Enable_stick.

◆ enable_stick()

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::enable_stick ( )
inline

Enforce permanent contact with penetrator, allowing for negative contact pressures.

777  {
778  Enable_stick=true;
779  }

References oomph::SurfaceContactElementBase< ELEMENT >::Enable_stick.

Referenced by MeltContactProblem< ELEMENT >::complete_problem_setup(), and SolarRadiationProblem< ELEMENT >::complete_problem_setup().

◆ fill_in_contribution_to_residuals()

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::fill_in_contribution_to_residuals ( Vector< double > &  residuals)
inlinevirtual

Return the residuals.

Reimplemented from oomph::GeneralisedElement.

815  {
816  fill_in_contribution_to_residuals_surface_contact(residuals);
817  }
oomph::SurfaceContactElementBase::fill_in_contribution_to_residuals_surface_contact
virtual void fill_in_contribution_to_residuals_surface_contact(Vector< double > &residuals)=0

References oomph::SurfaceContactElementBase< ELEMENT >::fill_in_contribution_to_residuals_surface_contact().

◆ fill_in_contribution_to_residuals_surface_contact()

template<class ELEMENT >
virtual void oomph::SurfaceContactElementBase< ELEMENT >::fill_in_contribution_to_residuals_surface_contact ( Vector< double > &  residuals)
protectedpure virtual

Helper function that actually calculates the residuals This small level of indirection is required to avoid calling fill_in_contribution_to_residuals in fill_in_contribution_to_jacobian which causes all kinds of pain if overloading later on

Implemented in oomph::HeatedLinearSurfaceContactElement< ELEMENT >, oomph::LinearSurfaceContactElement< ELEMENT >, and oomph::NonlinearSurfaceContactElement< ELEMENT >.

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::fill_in_contribution_to_residuals().

◆ get_interpolated_lagrange_p()

template<class ELEMENT >
double oomph::SurfaceContactElementBase< ELEMENT >::get_interpolated_lagrange_p ( const Vector< double > &  s)
inlineprotected

Get interpolated pressure (essentially a Lagrange multiplier that enforces the imposed boundary motion to ensure non-penetration or contact)

981  {
982  // Initialise pressure
983  double p=0;
984 
985  //Find out how many nodes there are
986  unsigned n_node = nnode();
987 
988  //Set up memory for the shape functions
989  Shape psi(n_node);
990 
991  // Evaluate shape function
992  shape_p(s,psi);
993 
994  // Build up Lagrange multiplier (pressure)
995  for (unsigned j=0;j<n_node;j++)
996  {
997  // Cast to a boundary node
998  BoundaryNodeBase *bnod_pt =
999  dynamic_cast<BoundaryNodeBase*>(node_pt(j));
1000 
1001  // Get the index of the first nodal value associated with
1002  // this FaceElement
1003  unsigned first_index=
1004  bnod_pt->index_of_first_value_assigned_by_face_element(Contact_id);
1005 
1006  // Pressure (Lagrange multiplier) is the first (and only) additional
1007  // value created by this face element
1008  p+=node_pt(j)->value(first_index)*psi[j];
1009  }
1010  return p;
1011  }
p
float * p
Definition: Tutorial_Map_using.cpp:9
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::Node::value
double value(const unsigned &i) const
Definition: nodes.cc:2408
oomph::SurfaceContactElementBase::shape_p
void shape_p(const Vector< double > &s, Shape &psi) const
Shape fct for lagrange multiplier.
Definition: contact_elements.h:913
s
RealScalar s
Definition: level1_cplx_impl.h:130
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References oomph::SurfaceContactElementBase< ELEMENT >::Contact_id, oomph::BoundaryNodeBase::index_of_first_value_assigned_by_face_element(), j, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), p, s, oomph::SurfaceContactElementBase< ELEMENT >::shape_p(), and oomph::Node::value().

Referenced by oomph::NonlinearSurfaceContactElement< ELEMENT >::output(), and oomph::LinearSurfaceContactElement< ELEMENT >::output().

◆ is_stick_enabled()

template<class ELEMENT >
bool oomph::SurfaceContactElementBase< ELEMENT >::is_stick_enabled ( )
inline

Do we allow only proper non-penetration (without "stick", i.e. no negative contact pressures)?

808  {
809  return Enable_stick;
810  }

References oomph::SurfaceContactElementBase< ELEMENT >::Enable_stick.

◆ output() [1/4]

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::output ( FILE *  file_pt)
inlinevirtual

C_style output function.

Reimplemented from oomph::FiniteElement.

893  {FiniteElement::output(file_pt);}
oomph::FiniteElement::output
virtual void output(std::ostream &outfile)
Definition: elements.h:3050

References oomph::FiniteElement::output().

◆ output() [2/4]

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::output ( FILE *  file_pt,
const unsigned &  n_plot 
)
inlinevirtual

C-style output function.

Reimplemented from oomph::FiniteElement.

897  {FiniteElement::output(file_pt,n_plot);}

References oomph::FiniteElement::output().

◆ output() [3/4]

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::output ( std::ostream &  outfile)
inlinevirtual

Output function.

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::LinearSurfaceContactElement< ELEMENT >, and oomph::NonlinearSurfaceContactElement< ELEMENT >.

901  {
902  unsigned n_plot=5;
903  FiniteElement::output(outfile,n_plot);
904  }

References oomph::FiniteElement::output().

◆ output() [4/4]

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::output ( std::ostream &  outfile,
const unsigned &  n_plot 
)
inlinevirtual

Output function.

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::HeatedLinearSurfaceContactElement< ELEMENT >, oomph::LinearSurfaceContactElement< ELEMENT >, and oomph::NonlinearSurfaceContactElement< ELEMENT >.

908  {
909  FiniteElement::output(outfile,n_plot);
910  }

References oomph::FiniteElement::output().

◆ penetration()

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::penetration ( const Vector< double > &  x,
const Vector< double > &  n,
double &  d,
bool &  intersection 
) const
inlineprotected

Work out penetration of point.

1025  {
1026  Penetrator_pt->penetration(x,n,d,intersection);
1027  }
n
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
oomph::Penetrator::penetration
virtual void penetration(const Vector< double > &x, const Vector< double > &n, double &d, bool &intersection) const =0
plotDoE.x
list x
Definition: plotDoE.py:28

References n, oomph::Penetrator::penetration(), oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_pt, and plotDoE::x.

Referenced by oomph::NonlinearSurfaceContactElement< ELEMENT >::output(), and oomph::LinearSurfaceContactElement< ELEMENT >::output().

◆ penetrator_pt()

template<class ELEMENT >
Penetrator* oomph::SurfaceContactElementBase< ELEMENT >::penetrator_pt ( ) const
inline

Pointer to penetrator.

839  {
840  return Penetrator_pt;
841  }

References oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_pt.

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::set_penetrator_pt().

◆ s_max()

template<class ELEMENT >
double oomph::SurfaceContactElementBase< ELEMENT >::s_max ( ) const
inlinevirtual

Max. value of local coordinate.

Reimplemented from oomph::FiniteElement.

794  {
795  return FaceGeometry<ELEMENT>::s_max();
796  }

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::shape_i(), oomph::SurfaceContactElementBase< ELEMENT >::shape_p(), and oomph::SurfaceContactElementBase< ELEMENT >::SurfaceContactElementBase().

◆ s_min()

template<class ELEMENT >
double oomph::SurfaceContactElementBase< ELEMENT >::s_min ( ) const
inlinevirtual

Min value of local coordinate.

Reimplemented from oomph::FiniteElement.

789  {
790  return FaceGeometry<ELEMENT>::s_min();
791  }

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::shape_i(), oomph::SurfaceContactElementBase< ELEMENT >::shape_p(), and oomph::SurfaceContactElementBase< ELEMENT >::SurfaceContactElementBase().

◆ set_penetrator_pt()

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::set_penetrator_pt ( Penetrator *  penetrator_pt)
inline

Set pointer to penetrator.

845  {
846  Penetrator_pt=penetrator_pt;
847  Vector<std::pair<Data*,unsigned> >
848  eq_data(Penetrator_pt->equilibrium_data());
849  unsigned n=eq_data.size();
850  Penetrator_eq_data_data_index.resize(n,-1);
851  Penetrator_eq_data_index.resize(n,-1);
852  Penetrator_eq_data_type.resize(n,-1);
853  for (unsigned i=0;i<n;i++)
854  {
855  if (eq_data[i].first!=0)
856  {
857  bool is_duplicate=false;
858  unsigned nnod=nnode();
859  for (unsigned j=0;j<nnod;j++)
860  {
861  if (eq_data[i].first==node_pt(j))
862  {
863  Penetrator_eq_data_type[i]=Nodal_data;
864  Penetrator_eq_data_data_index[i]=j;
865  Penetrator_eq_data_index[i]=eq_data[i].second;
866  is_duplicate=true;
867  break;
868  }
869  if (eq_data[i].first==
870  dynamic_cast<SolidNode*>(node_pt(j))->variable_position_pt())
871  {
872  Penetrator_eq_data_type[i]=Nodal_position_data;
873  Penetrator_eq_data_data_index[i]=j;
874  Penetrator_eq_data_index[i]=eq_data[i].second;
875  is_duplicate=true;
876  break;
877  }
878  }
879 
880  if (!is_duplicate)
881  {
882  Penetrator_eq_data_type[i]=External_data;
883  Penetrator_eq_data_data_index[i]=
884  this->add_external_data(eq_data[i].first);
885  Penetrator_eq_data_index[i]=eq_data[i].second;
886  }
887  }
888  }
889  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
oomph::GeneralisedElement::add_external_data
unsigned add_external_data(Data *const &data_pt, const bool &fd=true)
Definition: elements.cc:307
oomph::Penetrator::equilibrium_data
virtual Vector< std::pair< Data *, unsigned > > equilibrium_data()
Definition: contact_elements.h:195
oomph::SurfaceContactElementBase::Penetrator_eq_data_type
Vector< int > Penetrator_eq_data_type
Definition: contact_elements.h:1047
oomph::SurfaceContactElementBase::Penetrator_eq_data_index
Vector< int > Penetrator_eq_data_index
Definition: contact_elements.h:1055
oomph::SurfaceContactElementBase::penetrator_pt
Penetrator * penetrator_pt() const
Pointer to penetrator.
Definition: contact_elements.h:838
oomph::SurfaceContactElementBase::Penetrator_eq_data_data_index
Vector< int > Penetrator_eq_data_data_index
Definition: contact_elements.h:1063
oomph::TemplateFreeContactElementBase::Nodal_position_data
@ Nodal_position_data
Definition: contact_elements.h:637
oomph::TemplateFreeContactElementBase::External_data
@ External_data
Definition: contact_elements.h:638
oomph::TemplateFreeContactElementBase::Nodal_data
@ Nodal_data
Definition: contact_elements.h:636

References oomph::GeneralisedElement::add_external_data(), oomph::Penetrator::equilibrium_data(), oomph::TemplateFreeContactElementBase::External_data, i, j, n, oomph::FiniteElement::nnode(), oomph::TemplateFreeContactElementBase::Nodal_data, oomph::TemplateFreeContactElementBase::Nodal_position_data, oomph::FiniteElement::node_pt(), oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_eq_data_data_index, oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_eq_data_index, oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_eq_data_type, oomph::SurfaceContactElementBase< ELEMENT >::penetrator_pt(), and oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_pt.

Referenced by MeltContactProblem< ELEMENT >::complete_problem_setup(), SolarRadiationProblem< ELEMENT >::complete_problem_setup(), ContactProblem< ELEMENT >::complete_problem_setup(), ContactProblem< ELEMENT >::create_contact_elements(), and ContactProblem< ELEMENT >::switch_to_force_control().

◆ shape_i()

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::shape_i ( const Vector< double > &  s,
Shape &  psi 
) const
inline

Top hat function used in discretising either penetration or the contact pressure

949  {
950  const double smin=this->s_min();
951  const double smax=this->s_max();
952  const double sl=smin+0.25*(smax-smin);
953  const double sr=smin+0.75*(smax-smin);
954  if (s[0]<=sl)
955  {
956  psi[0]=1.0;
957  psi[1]=0.0;
958  psi[2]=0.0;
959  }
960  else if (s[0]<=sr)
961  {
962  psi[0]=0.0;
963  psi[1]=1.0;
964  psi[2]=0.0;
965  }
966  else
967  {
968  psi[0]=0.0;
969  psi[1]=0.0;
970  psi[2]=1.0;
971  }
972  }
smax
float BLASFUNC() smax(int *, float *, int *)
smin
float BLASFUNC() smin(int *, float *, int *)

References s, oomph::SurfaceContactElementBase< ELEMENT >::s_max(), oomph::SurfaceContactElementBase< ELEMENT >::s_min(), smax(), and smin().

◆ shape_p()

template<class ELEMENT >
void oomph::SurfaceContactElementBase< ELEMENT >::shape_p ( const Vector< double > &  s,
Shape &  psi 
) const
inline

Shape fct for lagrange multiplier.

914  {
915  if (this->use_isoparametric_flag())
916  {
917  FaceGeometry<ELEMENT>::shape(s,psi);
918  }
919  else
920  {
921  const double smin=this->s_min();
922  const double smax=this->s_max();
923  const double sl=smin+0.25*(smax-smin);
924  const double sr=smin+0.75*(smax-smin);
925  if (s[0]<=sl)
926  {
927  psi[0]=1.0;
928  psi[1]=0.0;
929  psi[2]=0.0;
930  }
931  else if (s[0]<=sr)
932  {
933  psi[0]=0.0;
934  psi[1]=1.0;
935  psi[2]=0.0;
936  }
937  else
938  {
939  psi[0]=0.0;
940  psi[1]=0.0;
941  psi[2]=1.0;
942  }
943  }
944  }
oomph::TemplateFreeContactElementBase::use_isoparametric_flag
bool use_isoparametric_flag() const
Definition: contact_elements.h:561
oomph::OneDimLagrange::shape
void shape(const double &s, double *Psi)
Definition: shape.h:564

References s, oomph::SurfaceContactElementBase< ELEMENT >::s_max(), oomph::SurfaceContactElementBase< ELEMENT >::s_min(), oomph::OneDimLagrange::shape(), smax(), smin(), and oomph::TemplateFreeContactElementBase::use_isoparametric_flag().

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::get_interpolated_lagrange_p().

◆ zeta_nodal()

template<class ELEMENT >
double oomph::SurfaceContactElementBase< ELEMENT >::zeta_nodal ( const unsigned &  n,
const unsigned &  k,
const unsigned &  i 
) const
inlinevirtual

Specify the values of the "global" intrinsic coordinate, zeta, of a compound geometric object (a mesh of elements) when the element is viewied as a sub-geometric object. The default assumption is that the element will be treated as a sub-geometric object in a bulk Mesh of other elements (geometric objects). The "global" coordinate of the compound geometric object is simply the Eulerian coordinate, x. The second default assumption is that the coordinate zeta will be stored at the nodes and interpolated using the shape functions of the element. This function returns the value of zeta stored at local node n, where k is the type of coordinate and i is the coordinate direction. The function is virtual so that it can be overloaded by different types of element: FaceElements and SolidFiniteElements

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::HeatedLinearSurfaceContactElement< ELEMENT >.

783  {
784  return oomph::FiniteElement::zeta_nodal(n,k,i);
785  }
oomph::FiniteElement::zeta_nodal
virtual double zeta_nodal(const unsigned &n, const unsigned &k, const unsigned &i) const
Definition: elements.h:2722
k
char char char int int * k
Definition: level2_impl.h:374

References i, k, and oomph::FiniteElement::zeta_nodal().

Member Data Documentation

◆ Contact_id

template<class ELEMENT >
unsigned oomph::SurfaceContactElementBase< ELEMENT >::Contact_id
protected

ID of the contact constraint (used for the identification of the nodal value that corresponds to the pressure-like Lagrange multipliers.

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::get_interpolated_lagrange_p(), and oomph::SurfaceContactElementBase< ELEMENT >::SurfaceContactElementBase().

◆ Enable_stick

template<class ELEMENT >
bool oomph::SurfaceContactElementBase< ELEMENT >::Enable_stick
protected

Do we allow only proper non-penetration (without "stick", i.e. no negative contact pressures)?

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::disable_stick(), oomph::SurfaceContactElementBase< ELEMENT >::enable_stick(), oomph::SurfaceContactElementBase< ELEMENT >::is_stick_enabled(), and oomph::SurfaceContactElementBase< ELEMENT >::SurfaceContactElementBase().

◆ Penetrator_eq_data_data_index

template<class ELEMENT >
Vector<int> oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_eq_data_data_index
protected

Vector containing the index of the value in the Data object (e.g. the node number in the current element) that stores the unknown that is determined from the Penetrator's i-th equilibrium of forces equation. -1 if the relevant equilibrium equation is not employed to determine an unknown that's associated with the current element.

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::set_penetrator_pt().

◆ Penetrator_eq_data_index

template<class ELEMENT >
Vector<int> oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_eq_data_index
protected

Vector containing the index of the Data object (e.g. the node number in the current element) that stores the unknown that is determined from the Penetrator's i-th equilibrium of forces equation. -1 if the relevant equilibrium equation is not employed to determine an unknown that's associated with the current element.

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::set_penetrator_pt().

◆ Penetrator_eq_data_type

template<class ELEMENT >
Vector<int> oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_eq_data_type
protected

Vector containing the type of Data (labeled by the class' private enumeration) that is determined (via pseudo-hijacking) from the Penetrator's i-th equilibrium of forces equation. -1 if this relevant equilibrium equation is not employed to determine an unknown that's associated with the current element.

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::set_penetrator_pt().

◆ Penetrator_pt

template<class ELEMENT >
Penetrator* oomph::SurfaceContactElementBase< ELEMENT >::Penetrator_pt
protected

Pointer to penetrator.

Referenced by oomph::SurfaceContactElementBase< ELEMENT >::penetration(), oomph::SurfaceContactElementBase< ELEMENT >::penetrator_pt(), oomph::SurfaceContactElementBase< ELEMENT >::set_penetrator_pt(), and oomph::SurfaceContactElementBase< ELEMENT >::SurfaceContactElementBase().

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