gmpMechanicalContact.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/

#ifndef _GEMA_PLUGIN_MECHANICAL_CONTACT_H_
#define _GEMA_PLUGIN_MECHANICAL_CONTACT_H_

#include "gmpMechanicalConfig.h"
#include <gmContactBoundaryConditionAccessor.h>
#include <gmGaussAccessor.h>

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMechanicalContact
{
public:
  //constructors
  GmpMechanicalContact(void);

  //destructor
  virtual ~GmpMechanicalContact(void);

  // Fills the node coordinates at the element border
  virtual void fillBorderNodeCoordinates(const GmElement* e, GmValueAccessor* nodeAc, GmMatrix& X, int border, int d, bool transposed = false) const;
  // 

  // Translate edges in nodes coordinates
  virtual void edgeToNodes(const GmContactBoundaryCondition* cbc, int bcIndex, int nbc, int d, GmValueAccessor* nodeAc, QList<GmMatrix>& surfaces,
          QList<GmVector>& nodesSurf) const;

  //  Fills the possible nodes in contact pair
  virtual void nodeToSurfaceSearch(QList<GmMatrix>& surfMaster, QList<GmMatrix>& surfSlave, QList<GmVector> nodesSurf1, QList<GmVector> nodesSurf2, QList<int>& nodesSlave, QList<GmMatrix>& nodesMaster) const;
  //

  // Fills actualized displacements vector for use in global coordinates actualized in each step
  virtual void coordinatesGlobal(GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, QList<double>& slaveCoordinates, QList<double>& masterCoordinates, QList<int>& nodesSlave, 
          QList<GmMatrix>& nodesMaster);
  //

  // Fills the impenetration vector for surfaces in contact
  virtual void calVectors(QList<int>& nodesSlave, QList<GmMatrix>& nodesMaster, GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, GmVector& lengthBorder, 
          GmVector& xi_c, GmVector& impVec, GmVector& slipVec, bool& hasContact);
  //

  // Fills the global force vector with contact contribution
  virtual void contactForce(GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, QList<int> nodesSlave, QList<GmMatrix> nodesMaster, GmVector lengthBorder, GmVector xi_c, GmVector impVec, 
          GmVector slipVec, double Kn, double Ks, double muAc, QList<GmMatrix>& Fc);
  //

  // Fills the global tangent matrix with contact contribution
  virtual void contactMatrix(GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, QList<int> nodesSlave, QList<GmMatrix> nodesMaster, GmVector lengthBorder, GmVector xi_c, GmVector impVec, 
          GmVector slipVec, double Kn, double Ks, double muAc, QList<GmMatrix>& Kc);
  //

  // Fills dofMap nodes to fill global tangent matrix
  virtual void contactNodes(GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, QList<int> nodesSlave, QList<GmMatrix> nodesMaster, GmVector impVec, QList<int>& nodesInContact, QList<GmVector>& dofMatrix);
  //

  // Write contact matrix to fill global tangent matrix
  virtual void translateMatrix(GmVector impVec, QList<int> nodesSlave, QList<GmMatrix> nodesMaster, QList<int> nodesInContact, QList<GmMatrix> Kc, GmMatrix& Kcc, QList<GmVector>& dofMatrix);
  //

  // Write contact forces to fill global vector internal forces
  virtual void translateContactForce(GmVector impVec, QList<int> nodesSlave, QList<GmMatrix> nodesMaster, QList<int> nodesInContact, QList<GmMatrix> Fc, GmMatrix& Fcc, QList<GmVector>& dofMatrix);
  //
};
#endif