linear_elasticity_with_contact.h

Go to the documentation of this file.

#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