Public Types | Public Member Functions | Static Public Member Functions | Protected Attributes | Static Protected Attributes | List of all members
oomph::LinearElasticityEquationsBase< DIM > Class Template Referenceabstract

#include <linear_elasticity_elements.h>

+ Inheritance diagram for oomph::LinearElasticityEquationsBase< DIM >:

Public Types

typedef void(* BodyForceFctPt) (const double &t, const Vector< double > &x, Vector< double > &b)
 
- 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

virtual unsigned u_index_linear_elasticity (const unsigned i) const
 
double d2u_dt2_linear_elasticity (const unsigned &n, const unsigned &i) const
 d^2u/dt^2 at local node n More...
 
void interpolated_u_linear_elasticity (const Vector< double > &s, Vector< double > &disp) const
 Compute vector of FE interpolated displacement u at local coordinate s. More...
 
double interpolated_u_linear_elasticity (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated displacement u[i] at local coordinate s. More...
 
 LinearElasticityEquationsBase ()
 
ElasticityTensor *& elasticity_tensor_pt ()
 Return the pointer to the elasticity_tensor. More...
 
double E (const unsigned &i, const unsigned &j, const unsigned &k, const unsigned &l) const
 Access function to the entries in the elasticity tensor. More...
 
const double & lambda_sq () const
 Access function for timescale ratio (nondim density) More...
 
double *& lambda_sq_pt ()
 Access function for pointer to timescale ratio (nondim density) More...
 
BodyForceFctPt & body_force_fct_pt ()
 Access function: Pointer to body force function. More...
 
BodyForceFctPt body_force_fct_pt () const
 Access function: Pointer to body force function (const version) More...
 
void enable_inertia ()
 Switch on solid inertia. More...
 
void disable_inertia ()
 Switch off solid inertia. More...
 
bool is_inertia_enabled () const
 Access function to flag that switches inertia on/off (const version) More...
 
virtual void pin_elemental_redundant_nodal_solid_pressures ()
 Pin the element's redundant solid pressures (needed for refinement) More...
 
virtual void get_stress (const Vector< double > &s, DenseMatrix< double > &sigma) const =0
 
void get_strain (const Vector< double > &s, DenseMatrix< double > &strain) const
 Return the strain tensor. More...
 
void body_force (const Vector< double > &x, Vector< double > &b) const
 
unsigned ndof_types () const
 
void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- 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 (std::ostream &outfile)
 
virtual void output (std::ostream &outfile, const unsigned &n_plot)
 
virtual void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output (FILE *file_pt)
 
virtual void output (FILE *file_pt, const unsigned &n_plot)
 
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)
 
- 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
 

Static Public Member Functions

static void pin_redundant_nodal_solid_pressures (const Vector< GeneralisedElement * > &element_pt)
 

Protected Attributes

ElasticityTensor * Elasticity_tensor_pt
 Pointer to the elasticity tensor. More...
 
double * Lambda_sq_pt
 Timescale ratio (non-dim. density) More...
 
bool Unsteady
 Flag that switches inertia on/off. More...
 
BodyForceFctPt Body_force_fct_pt
 Pointer to body force 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 Protected Attributes

static double Default_lambda_sq_value = 1.0
 Static default value for timescale ratio (1.0 – for natural scaling) 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
 

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
 
- 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)
 
virtual void fill_in_contribution_to_residuals (Vector< double > &residuals)
 
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)
 

Detailed Description

template<unsigned DIM>
class oomph::LinearElasticityEquationsBase< DIM >

A base class for elements that solve the equations of linear elasticity in Cartesian coordinates. Combines a few generic functions that are shared by LinearElasticityEquations and LinearElasticityEquationsWithPressure (hierher: The latter don't exist yet but will be written as soon as somebody needs them...)

Member Typedef Documentation

◆ BodyForceFctPt

template<unsigned DIM>
typedef void(* oomph::LinearElasticityEquationsBase< DIM >::BodyForceFctPt) (const double &t, const Vector< double > &x, Vector< double > &b)

Function pointer to function that specifies the body force as a function of the Cartesian coordinates and time FCT(t,x,b) – x and b are Vectors!

Constructor & Destructor Documentation

◆ LinearElasticityEquationsBase()

template<unsigned DIM>
oomph::LinearElasticityEquationsBase< DIM >::LinearElasticityEquationsBase ( )
inline

Constructor: Set null pointers for constitutive law and for isotropic growth function. Set physical parameter values to default values, switch on inertia and set body force to zero.

166  : Elasticity_tensor_pt(0),
167  Lambda_sq_pt(&Default_lambda_sq_value),
168  Unsteady(true),
169  Body_force_fct_pt(0)
170  {
171  }
oomph::LinearElasticityEquationsBase::Elasticity_tensor_pt
ElasticityTensor * Elasticity_tensor_pt
Pointer to the elasticity tensor.
Definition: linear_elasticity_elements.h:352
oomph::LinearElasticityEquationsBase::Default_lambda_sq_value
static double Default_lambda_sq_value
Static default value for timescale ratio (1.0 – for natural scaling)
Definition: linear_elasticity_elements.h:364
oomph::LinearElasticityEquationsBase::Lambda_sq_pt
double * Lambda_sq_pt
Timescale ratio (non-dim. density)
Definition: linear_elasticity_elements.h:355
oomph::LinearElasticityEquationsBase::Body_force_fct_pt
BodyForceFctPt Body_force_fct_pt
Pointer to body force function.
Definition: linear_elasticity_elements.h:361
oomph::LinearElasticityEquationsBase::Unsteady
bool Unsteady
Flag that switches inertia on/off.
Definition: linear_elasticity_elements.h:358

Member Function Documentation

◆ body_force()

template<unsigned DIM>
void oomph::LinearElasticityEquationsBase< DIM >::body_force ( const Vector< double > &  x,
Vector< double > &  b 
) const
inline

Evaluate body force at Eulerian coordinate x at present time (returns zero vector if no body force function pointer has been set)

271  {
272  // If no function has been set, return zero vector
273  if (Body_force_fct_pt == 0)
274  {
275  // Get spatial dimension of element
276  unsigned n = dim();
277  for (unsigned i = 0; i < n; i++)
278  {
279  b[i] = 0.0;
280  }
281  }
282  else
283  {
284  // Get time from timestepper of first node (note that this must
285  // work -- body force only makes sense for elements that can be
286  // deformed and given that the deformation of solid finite elements
287  // is controlled by their nodes, nodes must exist!)
288  double time = node_pt(0)->time_stepper_pt()->time_pt()->time();
289 
290  // Now evaluate the body force
291  (*Body_force_fct_pt)(time, x, b);
292  }
293  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
n
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
b
Scalar * b
Definition: benchVecAdd.cpp:17
oomph::Data::time_stepper_pt
TimeStepper *& time_stepper_pt()
Return the pointer to the timestepper.
Definition: nodes.h:238
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::dim
unsigned dim() const
Definition: elements.h:2611
oomph::TimeStepper::time_pt
Time *const & time_pt() const
Access function for the pointer to time (const version)
Definition: timesteppers.h:572
oomph::Time::time
double & time()
Return the current value of the continuous time.
Definition: timesteppers.h:123
plotDoE.x
list x
Definition: plotDoE.py:28

References b, oomph::LinearElasticityEquationsBase< DIM >::Body_force_fct_pt, oomph::FiniteElement::dim(), i, n, oomph::FiniteElement::node_pt(), oomph::Time::time(), oomph::TimeStepper::time_pt(), oomph::Data::time_stepper_pt(), and plotDoE::x.

◆ body_force_fct_pt() [1/2]

template<unsigned DIM>
BodyForceFctPt& oomph::LinearElasticityEquationsBase< DIM >::body_force_fct_pt ( )
inline

Access function: Pointer to body force function.

202  {
203  return Body_force_fct_pt;
204  }

References oomph::LinearElasticityEquationsBase< DIM >::Body_force_fct_pt.

Referenced by oomph::RefineableLinearElasticityEquations< DIM >::further_build().

◆ body_force_fct_pt() [2/2]

template<unsigned DIM>
BodyForceFctPt oomph::LinearElasticityEquationsBase< DIM >::body_force_fct_pt ( ) const
inline

Access function: Pointer to body force function (const version)

208  {
209  return Body_force_fct_pt;
210  }

References oomph::LinearElasticityEquationsBase< DIM >::Body_force_fct_pt.

◆ d2u_dt2_linear_elasticity()

template<unsigned DIM>
double oomph::LinearElasticityEquationsBase< DIM >::d2u_dt2_linear_elasticity ( const unsigned &  n,
const unsigned &  i 
) const
inline

d^2u/dt^2 at local node n

69  {
70  // Get the timestepper
71  TimeStepper* time_stepper_pt = node_pt(n)->time_stepper_pt();
72 
73  // Storage for the derivative - initialise to 0
74  double d2u_dt2 = 0.0;
75 
76  // If we are doing an unsteady solve then calculate the derivative
77  if (!time_stepper_pt->is_steady())
78  {
79  // Get the nodal index
80  const unsigned u_nodal_index = u_index_linear_elasticity(i);
81 
82  // Get the number of values required to represent history
83  const unsigned n_time = time_stepper_pt->ntstorage();
84 
85  // Loop over history values
86  for (unsigned t = 0; t < n_time; t++)
87  {
88  // Add the contribution to the derivative
89  d2u_dt2 +=
90  time_stepper_pt->weight(2, t) * nodal_value(t, n, u_nodal_index);
91  }
92  }
93 
94  return d2u_dt2;
95  }
oomph::FiniteElement::nodal_value
double nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2593
oomph::GeomObject::time_stepper_pt
TimeStepper *& time_stepper_pt()
Definition: geom_objects.h:192
oomph::LinearElasticityEquationsBase::u_index_linear_elasticity
virtual unsigned u_index_linear_elasticity(const unsigned i) const
Definition: linear_elasticity_elements.h:62
oomph::TimeStepper::ntstorage
unsigned ntstorage() const
Definition: timesteppers.h:601
oomph::TimeStepper::weight
virtual double weight(const unsigned &i, const unsigned &j) const
Access function for j-th weight for the i-th derivative.
Definition: timesteppers.h:594
oomph::TimeStepper::is_steady
bool is_steady() const
Definition: timesteppers.h:389
plotPSD.t
t
Definition: plotPSD.py:36

References i, oomph::TimeStepper::is_steady(), n, oomph::FiniteElement::nodal_value(), oomph::FiniteElement::node_pt(), oomph::TimeStepper::ntstorage(), plotPSD::t, oomph::GeomObject::time_stepper_pt(), oomph::Data::time_stepper_pt(), oomph::LinearElasticityEquationsBase< DIM >::u_index_linear_elasticity(), and oomph::TimeStepper::weight().

◆ disable_inertia()

template<unsigned DIM>
void oomph::LinearElasticityEquationsBase< DIM >::disable_inertia ( )
inline

Switch off solid inertia.

221  {
222  Unsteady = false;
223  }

References oomph::LinearElasticityEquationsBase< DIM >::Unsteady.

◆ E()

template<unsigned DIM>
double oomph::LinearElasticityEquationsBase< DIM >::E ( const unsigned &  i,
const unsigned &  j,
const unsigned &  k,
const unsigned &  l 
) const
inline

Access function to the entries in the elasticity tensor.

184  {
185  return (*Elasticity_tensor_pt)(i, j, k, l);
186  }
k
char char char int int * k
Definition: level2_impl.h:374
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References oomph::LinearElasticityEquationsBase< DIM >::Elasticity_tensor_pt, i, j, and k.

◆ elasticity_tensor_pt()

template<unsigned DIM>
ElasticityTensor*& oomph::LinearElasticityEquationsBase< DIM >::elasticity_tensor_pt ( )
inline

Return the pointer to the elasticity_tensor.

175  {
176  return Elasticity_tensor_pt;
177  }

References oomph::LinearElasticityEquationsBase< DIM >::Elasticity_tensor_pt.

Referenced by oomph::RefineableLinearElasticityEquations< DIM >::further_build().

◆ enable_inertia()

template<unsigned DIM>
void oomph::LinearElasticityEquationsBase< DIM >::enable_inertia ( )
inline

Switch on solid inertia.

215  {
216  Unsteady = true;
217  }

References oomph::LinearElasticityEquationsBase< DIM >::Unsteady.

◆ get_dof_numbers_for_unknowns()

template<unsigned DIM>
void oomph::LinearElasticityEquationsBase< DIM >::get_dof_numbers_for_unknowns ( std::list< std::pair< unsigned long, unsigned >> &  dof_lookup_list) const
inlinevirtual

Create a list of pairs for all unknowns in this element, so that the first entry in each pair contains the global equation number of the unknown, while the second one contains the number of the "DOF types" that this unknown is associated with. (Function can obviously only be called if the equation numbering scheme has been set up.)

Reimplemented from oomph::GeneralisedElement.

313  {
314  // temporary pair (used to store dof lookup prior to being added
315  // to list)
316  std::pair<unsigned long, unsigned> dof_lookup;
317 
318  // number of nodes
319  const unsigned n_node = this->nnode();
320 
321  // Integer storage for local unknown
322  int local_unknown = 0;
323 
324  // Loop over the nodes
325  for (unsigned n = 0; n < n_node; n++)
326  {
327  // Loop over dimension
328  for (unsigned i = 0; i < DIM; i++)
329  {
330  // If the variable is free
331  local_unknown = nodal_local_eqn(n, i);
332 
333  // ignore pinned values
334  if (local_unknown >= 0)
335  {
336  // store dof lookup in temporary pair: First entry in pair
337  // is global equation number; second entry is dof type
338  dof_lookup.first = this->eqn_number(local_unknown);
339  dof_lookup.second = i;
340  // dof_lookup.second = DIM;
341 
342  // add to list
343  dof_lookup_list.push_front(dof_lookup);
344  }
345  }
346  }
347  }
oomph::FiniteElement::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
oomph::GeneralisedElement::eqn_number
unsigned long eqn_number(const unsigned &ieqn_local) const
Definition: elements.h:704
DIM
#define DIM
Definition: linearised_navier_stokes_elements.h:44

References DIM, oomph::GeneralisedElement::eqn_number(), i, n, oomph::FiniteElement::nnode(), and oomph::FiniteElement::nodal_local_eqn().

◆ get_strain()

template<unsigned DIM>
void oomph::LinearElasticityEquationsBase< DIM >::get_strain ( const Vector< double > &  s,
DenseMatrix< double > &  strain 
) const

Return the strain tensor.

Compute the strain tensor at local coordinate s.

/////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////

Now fill in the entries of the strain tensor

49  {
50 #ifdef PARANOID
51  if ((strain.ncol() != DIM) || (strain.nrow() != DIM))
52  {
53  std::ostringstream error_message;
54  error_message << "Strain matrix is " << strain.ncol() << " x "
55  << strain.nrow() << ", but dimension of the equations is "
56  << DIM << std::endl;
57  throw OomphLibError(
58  error_message.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
59  }
60 
61  // Find out how many position types there are
62  unsigned n_position_type = this->nnodal_position_type();
63 
64  if (n_position_type != 1)
65  {
66  throw OomphLibError("LinearElasticity is not yet implemented for more "
67  "than one position type",
68  OOMPH_CURRENT_FUNCTION,
69  OOMPH_EXCEPTION_LOCATION);
70  }
71 #endif
72 
73 
74  // Find out how many nodes there are in the element
75  unsigned n_node = nnode();
76 
77  // Find the indices at which the local velocities are stored
78  unsigned u_nodal_index[DIM];
79  for (unsigned i = 0; i < DIM; i++)
80  {
81  u_nodal_index[i] = u_index_linear_elasticity(i);
82  }
83 
84  // Set up memory for the shape and derivative functions
85  Shape psi(n_node);
86  DShape dpsidx(n_node, DIM);
87 
88  // Call the derivatives of the shape functions
89  (void)dshape_eulerian(s, psi, dpsidx);
90 
91  // Calculate interpolated values of the derivative of global position
92  DenseMatrix<double> interpolated_dudx(DIM, DIM, 0.0);
93 
94  // Loop over nodes
95  for (unsigned l = 0; l < n_node; l++)
96  {
97  // Loop over velocity components
98  for (unsigned i = 0; i < DIM; i++)
99  {
100  // Get the nodal value
101  const double u_value = this->nodal_value(l, u_nodal_index[i]);
102 
103  // Loop over derivative directions
104  for (unsigned j = 0; j < DIM; j++)
105  {
106  interpolated_dudx(i, j) += u_value * dpsidx(l, j);
107  }
108  }
109  }
110 
112  for (unsigned i = 0; i < DIM; i++)
113  {
114  // Do upper half of matrix
115  // Note that j must be signed here for the comparison test to work
116  // Also i must be cast to an int
117  for (int j = (DIM - 1); j >= static_cast<int>(i); j--)
118  {
119  // Off diagonal terms
120  if (static_cast<int>(i) != j)
121  {
122  strain(i, j) =
123  0.5 * (interpolated_dudx(i, j) + interpolated_dudx(j, i));
124  }
125  // Diagonal terms will including growth factor when it comes back in
126  else
127  {
128  strain(i, i) = interpolated_dudx(i, i);
129  }
130  }
131  // Matrix is symmetric so just copy lower half
132  for (int j = (i - 1); j >= 0; j--)
133  {
134  strain(i, j) = strain(j, i);
135  }
136  }
137  }
oomph::DenseMatrix::nrow
unsigned long nrow() const
Return the number of rows of the matrix.
Definition: matrices.h:485
oomph::DenseMatrix::ncol
unsigned long ncol() const
Return the number of columns of the matrix.
Definition: matrices.h:491
oomph::FiniteElement::nnodal_position_type
unsigned nnodal_position_type() const
Definition: elements.h:2463
oomph::FiniteElement::dshape_eulerian
double dshape_eulerian(const Vector< double > &s, Shape &psi, DShape &dpsidx) const
Definition: elements.cc:3298
s
RealScalar s
Definition: level1_cplx_impl.h:130
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 DIM, i, j, oomph::DenseMatrix< T >::ncol(), oomph::DenseMatrix< T >::nrow(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, and s.

Referenced by oomph::RefineableLinearElasticityEquations< DIM >::get_Z2_flux(), and oomph::TLinearElasticityElement< DIM, NNODE_1D >::get_Z2_flux().

◆ get_stress()

template<unsigned DIM>
virtual void oomph::LinearElasticityEquationsBase< DIM >::get_stress ( const Vector< double > &  s,
DenseMatrix< double > &  sigma 
) const
pure virtual

Return the Cauchy stress tensor, as calculated from the elasticity tensor at specified local coordinate Virtual so separaete versions can (and must!) be provided for displacement and pressure-displacement formulations.

Implemented in oomph::LinearElasticityEquations< DIM >.

◆ interpolated_u_linear_elasticity() [1/2]

template<unsigned DIM>
double oomph::LinearElasticityEquationsBase< DIM >::interpolated_u_linear_elasticity ( const Vector< double > &  s,
const unsigned &  i 
) const
inline

Return FE interpolated displacement u[i] at local coordinate s.

129  {
130  // Find number of nodes
131  unsigned n_node = nnode();
132 
133  // Local shape function
134  Shape psi(n_node);
135 
136  // Find values of shape function
137  shape(s, psi);
138 
139  // Get nodal index at which i-th velocity is stored
140  unsigned u_nodal_index = u_index_linear_elasticity(i);
141 
142  // Initialise value of u
143  double interpolated_u = 0.0;
144 
145  // Loop over the local nodes and sum
146  for (unsigned l = 0; l < n_node; l++)
147  {
148  interpolated_u += nodal_value(l, u_nodal_index) * psi[l];
149  }
150 
151  return (interpolated_u);
152  }
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0

References i, oomph::FiniteElement::nnode(), oomph::FiniteElement::nodal_value(), s, oomph::FiniteElement::shape(), and oomph::LinearElasticityEquationsBase< DIM >::u_index_linear_elasticity().

◆ interpolated_u_linear_elasticity() [2/2]

template<unsigned DIM>
void oomph::LinearElasticityEquationsBase< DIM >::interpolated_u_linear_elasticity ( const Vector< double > &  s,
Vector< double > &  disp 
) const
inline

Compute vector of FE interpolated displacement u at local coordinate s.

100  {
101  // Find number of nodes
102  unsigned n_node = nnode();
103 
104  // Local shape function
105  Shape psi(n_node);
106 
107  // Find values of shape function
108  shape(s, psi);
109 
110  for (unsigned i = 0; i < DIM; i++)
111  {
112  // Index at which the nodal value is stored
113  unsigned u_nodal_index = u_index_linear_elasticity(i);
114 
115  // Initialise value of u
116  disp[i] = 0.0;
117 
118  // Loop over the local nodes and sum
119  for (unsigned l = 0; l < n_node; l++)
120  {
121  disp[i] += nodal_value(l, u_nodal_index) * psi[l];
122  }
123  }
124  }

References DIM, i, oomph::FiniteElement::nnode(), oomph::FiniteElement::nodal_value(), s, oomph::FiniteElement::shape(), and oomph::LinearElasticityEquationsBase< DIM >::u_index_linear_elasticity().

Referenced by oomph::RefineableLinearElasticityEquations< DIM >::get_interpolated_values().

◆ is_inertia_enabled()

template<unsigned DIM>
bool oomph::LinearElasticityEquationsBase< DIM >::is_inertia_enabled ( ) const
inline

Access function to flag that switches inertia on/off (const version)

227  {
228  return Unsteady;
229  }

References oomph::LinearElasticityEquationsBase< DIM >::Unsteady.

Referenced by oomph::RefineableLinearElasticityEquations< DIM >::further_build().

◆ lambda_sq()

template<unsigned DIM>
const double& oomph::LinearElasticityEquationsBase< DIM >::lambda_sq ( ) const
inline

Access function for timescale ratio (nondim density)

190  {
191  return *Lambda_sq_pt;
192  }

References oomph::LinearElasticityEquationsBase< DIM >::Lambda_sq_pt.

◆ lambda_sq_pt()

template<unsigned DIM>
double*& oomph::LinearElasticityEquationsBase< DIM >::lambda_sq_pt ( )
inline

Access function for pointer to timescale ratio (nondim density)

196  {
197  return Lambda_sq_pt;
198  }

References oomph::LinearElasticityEquationsBase< DIM >::Lambda_sq_pt.

Referenced by oomph::RefineableLinearElasticityEquations< DIM >::further_build().

◆ ndof_types()

template<unsigned DIM>
unsigned oomph::LinearElasticityEquationsBase< DIM >::ndof_types ( ) const
inlinevirtual

The number of "DOF types" that degrees of freedom in this element are sub-divided into: for now lump them all into one DOF type. Can be adjusted later

Reimplemented from oomph::GeneralisedElement.

300  {
301  return 2;
302  // return 1;
303  }

◆ pin_elemental_redundant_nodal_solid_pressures()

template<unsigned DIM>
virtual void oomph::LinearElasticityEquationsBase< DIM >::pin_elemental_redundant_nodal_solid_pressures ( )
inlinevirtual

Pin the element's redundant solid pressures (needed for refinement)

232 {}

Referenced by oomph::LinearElasticityEquationsBase< DIM >::pin_redundant_nodal_solid_pressures().

◆ pin_redundant_nodal_solid_pressures()

template<unsigned DIM>
static void oomph::LinearElasticityEquationsBase< DIM >::pin_redundant_nodal_solid_pressures ( const Vector< GeneralisedElement * > &  element_pt)
inlinestatic

Loop over all elements in Vector (which typically contains all the elements in a refineable solid mesh) and pin the nodal solid pressure degrees of freedom that are not being used. Function uses the member function

  • LinearElasticityEquationsBase<DIM>:: pin_elemental_redundant_nodal_pressure_dofs()

which is empty by default and should be implemented for elements with nodal solid pressure degrees of freedom (e.g. linear elasticity elements with continuous pressure interpolation.)

247  {
248  // Loop over all elements and call the function that pins their
249  // unused nodal solid pressure data
250  unsigned n_element = element_pt.size();
251  for (unsigned e = 0; e < n_element; e++)
252  {
253  dynamic_cast<LinearElasticityEquationsBase<DIM>*>(element_pt[e])
254  ->pin_elemental_redundant_nodal_solid_pressures();
255  }
256  }
e
Array< double, 1, 3 > e(1./3., 0.5, 2.)
oomph::LinearElasticityEquationsBase::pin_elemental_redundant_nodal_solid_pressures
virtual void pin_elemental_redundant_nodal_solid_pressures()
Pin the element's redundant solid pressures (needed for refinement)
Definition: linear_elasticity_elements.h:232

References e(), and oomph::LinearElasticityEquationsBase< DIM >::pin_elemental_redundant_nodal_solid_pressures().

◆ u_index_linear_elasticity()

template<unsigned DIM>
virtual unsigned oomph::LinearElasticityEquationsBase< DIM >::u_index_linear_elasticity ( const unsigned  i) const
inlinevirtual

Return the index at which the i-th unknown displacement component is stored. The default value, i, is appropriate for single-physics problems.

63  {
64  return i;
65  }

References i.

Referenced by oomph::LinearElasticityEquationsBase< DIM >::d2u_dt2_linear_elasticity(), oomph::RefineableLinearElasticityEquations< DIM >::get_interpolated_values(), and oomph::LinearElasticityEquationsBase< DIM >::interpolated_u_linear_elasticity().

Member Data Documentation

◆ Body_force_fct_pt

template<unsigned DIM>
BodyForceFctPt oomph::LinearElasticityEquationsBase< DIM >::Body_force_fct_pt
protected

Pointer to body force function.

Referenced by oomph::LinearElasticityEquationsBase< DIM >::body_force(), oomph::LinearElasticityEquationsBase< DIM >::body_force_fct_pt(), and oomph::RefineableLinearElasticityEquations< DIM >::further_build().

◆ Default_lambda_sq_value

template<unsigned DIM>
double oomph::LinearElasticityEquationsBase< DIM >::Default_lambda_sq_value = 1.0
staticprotected

Static default value for timescale ratio (1.0 – for natural scaling)

◆ Elasticity_tensor_pt

template<unsigned DIM>
ElasticityTensor* oomph::LinearElasticityEquationsBase< DIM >::Elasticity_tensor_pt
protected

Pointer to the elasticity tensor.

Referenced by oomph::LinearElasticityEquationsBase< DIM >::E(), oomph::LinearElasticityEquationsBase< DIM >::elasticity_tensor_pt(), and oomph::RefineableLinearElasticityEquations< DIM >::further_build().

◆ Lambda_sq_pt

template<unsigned DIM>
double* oomph::LinearElasticityEquationsBase< DIM >::Lambda_sq_pt
protected

Timescale ratio (non-dim. density)

Referenced by oomph::RefineableLinearElasticityEquations< DIM >::further_build(), oomph::LinearElasticityEquationsBase< DIM >::lambda_sq(), and oomph::LinearElasticityEquationsBase< DIM >::lambda_sq_pt().

◆ Unsteady

template<unsigned DIM>
bool oomph::LinearElasticityEquationsBase< DIM >::Unsteady
protected

Flag that switches inertia on/off.

Referenced by oomph::LinearElasticityEquationsBase< DIM >::disable_inertia(), oomph::LinearElasticityEquationsBase< DIM >::enable_inertia(), oomph::RefineableLinearElasticityEquations< DIM >::further_build(), and oomph::LinearElasticityEquationsBase< DIM >::is_inertia_enabled().

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