xfemPos.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2017 by Geomechanic group - Tecgraf/Puc-Rio
** 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_XFEM_POS_H_
#define _GEMA_XFEM_POS_H_
 
#include <assert.h>
#include "xfemTools.h"
 
class XfemPos : public XfemTools
{
public:
        //XfemPos(XfemMesh* _xmesh, XfemOptions& xfemOptions, GmLogCategory& logger);
        XfemPos(XfemMesh* _xmesh, XfemOptions& xfemOptions, GmSimulationData* simulation, QVector<GmpFemPhysics*>& physics, const GmLogCategory& logger);
 
  virtual ~XfemPos();
 
  virtual bool checkFractureTips(XfemMesh* xmesh) const;
  //set regions to check fracture propagation
  virtual bool setRegionsToCheck(XfemMesh* xmesh) const;
  //set regions to check fracture propagation
  virtual bool setRegionsToCheck3D(XfemMesh* xmesh) const;
  //check fracture propagation at the crack front
  bool setFrontToCheck3D(XfemMesh * xmesh) const;
  //Treatment for No inital cracks in the model
  virtual bool noInitialCracks(XfemMesh* xmesh) const;
  //set nonLocal Quad region
  virtual bool nonLocalQuad(XfemMesh* xmesh, XfemFracture* fracture, QList<GmVector>& headIp, QList<GmVector>& tailIp, int idCrack) const;
  //set nonLocal Circle region
  virtual bool nonLocalCircle(XfemMesh* xmesh, XfemFracture* fracture, QList<GmVector>& headIp, QList<GmVector>& tailIp) const;
  //get polygon coords
  virtual bool createPolygon(XfemMesh* xmesh, GmVector Crk, GmMatrix& pCoords) const;
  //get centroid 
  //virtual bool createCentroid(XfemMesh * xmesh, GmVector Crk, GmMatrix & polygonVert) const;
  //get circle coords
  //virtual bool createCircle(XfemMesh* xmesh, GmVector Crk, GmMatrix& polygonVert) const;
  //get polygon Area
  virtual bool polygonArea(GmMatrix Xnodes, double& S) const;
  //get circle Area
  //virtual bool circleArea(GmMatrix Xnodes, double& S) const;
  //check propagation of fracture in integration points inside a region
  virtual bool checkPropagation(XfemMesh* xmesh) const;
  //check propagation of fracture in integration points inside a region
  virtual bool checkPropagation3D(XfemMesh* xmesh) const;
  virtual bool testPropagation3D(XfemMesh* xmesh) const;
  //check propagation of fracture in integration points inside this element
  virtual bool checkCrackReactivation(XfemMesh* xmesh, GmVector crackPath, XfemEnrichedElementData* exData) const;
  //get average stress or strain field according to propagation criteria and weight function
  virtual bool averageStressStrainField(XfemMesh* xmesh, XfemFracture* fracture, GmVector& fieldValue, double& strength, bool isHead) const;
  //get average stress or strain field according to propagation criteria and weight function
  virtual bool averageStressStrainField3D(XfemMesh* xmesh, XfemFracture* fracture, int index, GmMatrix ROT2, GmVector& fieldValue, GmMatrix& fieldDir, bool& next) const;
  virtual bool testAverageStressStrainField3D(XfemMesh* xmesh, XfemFracture* fracture, GmCell* e, int index, GmMatrix ROT2, GmVector& fieldValue, GmMatrix& fieldDir, bool& next) const;
  //get average stress or strain field according to propagation criteria and weight function
  virtual bool averageUnitVectors3D(XfemMesh* xmesh, XfemFracture* fracture, int index, GmMatrix& ROT2) const;
  //extend Crack Plane 3D
  virtual bool extendCrackPlane3D(XfemMesh* xmesh, XfemFracture* fracture, XfemElement* xe, GmCell* e, int index, GmMatrix ROT2, GmVector& fieldValue, GmMatrix& fieldDir, GmMatrix& Xg, GmVector& v_n, GmMatrix& planeCoord, QList<int>& crackNodeId) const;
  //extend Crack Plane 3D
  virtual bool testExtendCrackPlane3D(XfemMesh* xmesh, XfemFracture* fracture, XfemElement* xe, GmCell* e, int index, GmMatrix ROT2, GmVector& fieldValue, GmMatrix& fieldDir, GmMatrix& Xg, GmVector& v_n, GmMatrix& planeCoord, QList<int>& crackNodeId, bool & next) const;
  //get stress or strain field according to propagation criteria and weight function
  virtual bool stressStrainField(XfemMesh* xmesh, XfemFracture* fracture, GmVector& fieldValue, double& strength, bool isHead) const;
  //average stress/strain field for uniform weight function
  virtual bool weightFunction(XfemMesh* xmesh, XfemFracture* fracture, GmVector Xip, double& w, bool isHead) const;
  //Allow to update crack geometry after propagation
  virtual bool updateGeomEnrich(XfemMesh* xmesh) const;
  //Allow to update crack geometry after propagation
  virtual bool updateGeomEnrich3D(XfemMesh* xmesh) const;
  //Allow to update crack geometry after propagation
  virtual bool updateCrackData3D(XfemMesh* xmesh, GmMatrix& planeCoord, QList<GmVector> aCrack, int index, int iCr, QList<int>& crackNodeId, int& nodeIndex, GmVector& pathData, bool& next, QList<GmVector>& ghostList) const;
  //treat propagation of fractures 
  virtual bool crackPropagation(XfemMesh* xmesh, XfemFracture* Fr, int index, int secondaryCrack, bool isHead, GmVector crackInfo, GmVector stress, GmVector B, QList<GmVector>& aCrack, int& nInterElem, bool& noCrossingNFReactivation) const;
  //check if this crack segment cuts other crack to restrict propagation until this other crack is touched not crossed
  virtual bool updateElmentIntersection(XfemMesh* xmesh, GmCell* e, XfemElement* xe, GmVector Xg, XfemEnrichedElementData* exData, QList<GmVector>& aCrack, int secondaryCrack) const;
  //Get Crack propagation direction 
  virtual bool crackPropDirection(XfemMesh* xmesh, XfemElement* xe, GmVector stress, GmVector aPStress, int nPS, bool isHead, GmVector crackInfo, GmMatrix& XGhost, GmVector& uPs, int index, bool& noCrossingNFReactivation) const;
  //update Crack Geometry coordinates
  virtual bool updateCrackGeom(XfemMesh* xmesh, XfemFracture* fracture, int index, int iNodePath, GmVector B, GmVector Xhtcoord, GmVector& mElemCutPos) const;
  //Update ghost nodes in Xmesh
  virtual bool updateGhostNodes(XfemMesh* xmesh, XfemFracture* fracture, int nodeIndex, GmVector coord, GmVector pathData, int iCr) const;
  //Update ghost nodes in Xmesh
  virtual bool updateGhostNodes3D(XfemMesh* xmesh, int nodeIndex, QList<GmVector> aCrack, QList<GmVector>& ghostList, GmVector pathData, int iCr) const;
  //Update subelements division in Xmesh
  virtual bool updatesubElemDiv(XfemMesh* xmesh, QList<int> enrichList) const;
  //Update subelements division in Xmesh
  virtual bool updatesubElemDiv3D(XfemMesh* xmesh, QList<int> enrichList) const;
  //treat propagation of a fracture in a cross type intersection
  virtual bool crossTypeCrackProp(XfemMesh* xmesh, XfemFracture* fracture, GmVector crackInfo, GmVector B, QList<GmVector> aCrack, bool isHead, int& iCr, int& nCracks) const;
  virtual bool setUpdateAssembler(XfemMesh* xmesh) const;
  // update Enrichment Functions
  virtual bool updateEnrichmentFunctions(XfemMesh* xmesh)const;
  //upadate elemCutInfo based on the number of fractures
  virtual bool upadateElemCutInfo(XfemMesh* xmesh, QList<int> enrichList) const;
  // Define DOF sets for intersected fractures
  virtual bool upadateDofSets(XfemMesh* xmesh) const;
  //update Crack Front Data
  bool updateCrackFrontData(XfemMesh * xmesh, QList<int> enrichList) const;
  // set data for integration points for each subelement
  virtual bool setInterEnrichedData(XfemMesh* xmesh, GmCell* e, int numRegularSubElements, int numFractureSubElements) const;
};
 
#endif