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:
 
  // Fills the node coordinates at the element border
  static void fillBorderNodeCoordinates(const GmElement* e, GmValueAccessor* nodeAc, GmMatrix& X, int border, int d, bool transposed = false);
  // 
 
  // Translate edges in nodes coordinates
  static void edgeToNodes(const GmContactBoundaryCondition* cbc, int bcIndex, int nbc, int d, GmValueAccessor* nodeAc, QVector<GmMatrix>& surfaces,
      QVector<GmVector>& nodesSurf);
 
  //  Fills the possible nodes in contact pair
  static void nodeToSurfaceSearch(const QVector<GmMatrix>& surfMaster, const QVector<GmMatrix>& surfSlave, const QVector<GmVector>& nodesSurf1, const QVector<GmVector>& nodesSurf2, QVector<int>& nodesSlave, QVector<GmMatrix>& nodesMaster);
  //
 
  // Fills actualized displacements vector for use in global coordinates actualized in each step
  static void coordinatesGlobal(GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, QVector<double>& slaveCoordinates, QVector<double>& masterCoordinates, const QVector<int>& nodesSlave,
          const QVector<GmMatrix>& nodesMaster);
  //
 
  // Fills the impenetration vector for surfaces in contact
  static void calVectors(const QVector<int>& nodesSlave, const QVector<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
  static void contactForce(GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, const QVector<int>& nodesSlave, const QVector<GmMatrix>& nodesMaster, const GmVector& lengthBorder, const GmVector& xi_c, const GmVector& impVec,
          const GmVector& slipVec, double Kn, double Ks, double muAc, QVector<GmMatrix>& Fc);
  //
 
  // Fills the global tangent matrix with contact contribution
  static void contactMatrix(GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, const QVector<int>& nodesSlave, const QVector<GmMatrix>& nodesMaster, const GmVector& lengthBorder, const GmVector& xi_c, const GmVector& impVec,
          const GmVector& slipVec, double Kn, double Ks, double muAc, QVector<GmMatrix>& Kc);
  //
 
  // Fills dofMap nodes to fill global tangent matrix
  static void contactNodes(GmValueAccessor* nodeAcc, GmValueAccessor*uAcc, const QVector<int>& nodesSlave, const QVector<GmMatrix>& nodesMaster, const GmVector& impVec, QVector<int>& nodesInContact, QVector<GmVector>& dofMatrix);
  //
 
  // Write contact matrix to fill global tangent matrix
  static void translateMatrix(const GmVector& impVec, const QVector<int>& nodesSlave, const QVector<GmMatrix>& nodesMaster, const QVector<int>& nodesInContact, const QVector<GmMatrix>& Kc, GmMatrix& Kcc, const QVector<GmVector>& dofMatrix);
  //
 
  // Write contact forces to fill global vector internal forces
  static void translateContactForce(const GmVector& impVec, const QVector<int>& nodesSlave, const QVector<GmMatrix>& nodesMaster, const QVector<int>& nodesInContact, const QVector<GmMatrix>& Fc, GmMatrix& Fcc, const QVector<GmVector>& dofMatrix);
  //
};
#endif