doxygen/libmesh/linear__elasticity__with__contact_8h_source.html

#ifndef LINEAR_ELASTICITY_WITH_CONTACT_H

#define LINEAR_ELASTICITY_WITH_CONTACT_H


#include "libmesh/dof_map.h"

#include "libmesh/nonlinear_implicit_system.h"


// define the subdomain IDs

#define TOP_SUBDOMAIN 2

#define BOTTOM_SUBDOMAIN 1


// define the boundary IDs in the mesh

#define MIN_Z_BOUNDARY 1

#define MAX_Z_BOUNDARY 2

#define CONTACT_BOUNDARY_LOWER 3

#define CONTACT_BOUNDARY_UPPER 4


using libMesh::DofMap;

using libMesh::NonlinearImplicitSystem;

using libMesh::dof_id_type;

using libMesh::Point;

using libMesh::Real;

using libMesh::Number;

using libMesh::MeshBase;

using libMesh::NumericVector;

using libMesh::SparseMatrix;


class LinearElasticityWithContact :

  public NonlinearImplicitSystem::ComputeResidualandJacobian

{

private:


  NonlinearImplicitSystem & _sys;


  Real _contact_penalty;


  std::map<dof_id_type, Real> _lambda_plus_penalty_values;


  std::map<dof_id_type, Real> _lambdas;


  std::map<dof_id_type, dof_id_type> _contact_node_map;


public:


  LinearElasticityWithContact(NonlinearImplicitSystem & sys_in,

                              Real contact_penalty_in);


  void set_contact_penalty(Real contact_penalty_in);


  Real get_contact_penalty() const;


  Real kronecker_delta(unsigned int i,

                       unsigned int j);


  Real elasticity_tensor(Real young_modulus,

                         Real poisson_ratio,

                         unsigned int i,

                         unsigned int j,

                         unsigned int k,

                         unsigned int l);


  void move_mesh(MeshBase & input_mesh,

                 const NumericVector<Number> & input_solution);


  void initialize_contact_load_paths();


  void add_contact_edge_elements();


  virtual void residual_and_jacobian (const NumericVector<Number> & soln,

                                      NumericVector<Number> * residual,

                                      SparseMatrix<Number> * jacobian,

                                      NonlinearImplicitSystem & /*sys*/);


  void compute_stresses();


  std::pair<Real, Real> update_lambdas();


  std::pair<Real, Real> get_least_and_max_gap_function();

};


#endif