xfemMechanic.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_MECHANIC_H_
#define _GEMA_XFEM_MECHANIC_H_

#include <gmpFemPhysicsCoupled.h>
#include <gmpMechanicalPhysics.h>
#include <gmpMechanicPoint.h>

#include "xfemMesh.h"

class GmpMechanicalPhysics;
class GmpHydraulic;
class GmGaussAccessor;

template <class T, class F> class XfemPhysicsWrapper: public T
{
public:
  XfemPhysicsWrapper(const char* physicsType, const char* objType, GmSimulationData* simulation, QString id, const GmLogCategory& logger)
    : T(objType, simulation, id, "",
        dynamic_cast<F*>(simulation->loader()->objectFactory(physicsType))->matFactory(simulation),
        logger) {}

protected:
  virtual GmGaussAccessor* createGaussAttributeAccessor(GmElementMesh* mesh, QString id, int snum, bool locked, 
                                                        Unit desiredUnit, const GmLogCategory& logger) const
  {
    XfemMesh* xmesh = dynamic_cast<XfemMesh*>(mesh);
    assert(xmesh);
    return xmesh->enrichedGaussAttributeAccessor(id, snum, locked, desiredUnit, logger);
  }
};

class XfemMechanic : public GmpFemPhysicsCoupled
{
public:
  // constructor
  XfemMechanic(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
    const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
  // destructor
  virtual ~XfemMechanic();

  // Comments on the base class
  virtual const char* pluginName() const { return "Xfem"; }

  virtual const char* pluginType() const { return _pluginType; }
  
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();

  // returns the list of elements that had their dofs changed since the physics was create
  virtual QList<int> changedElements() const;

  // clear the internal list used to track dof changes
  virtual void clearChangedElements();
  // clear the internal list used to track dof changes
  virtual bool fillFixedEnrichedBC(QList<int>& nodes, QList<int>&dof, QList<double>& values) const;

  // Comments on the base class
  virtual bool dofByElement(bool* fixed, bool* addOnly, bool* trackChanges) const;

  virtual const GmElementDof* dofMapping(const GmElement* e) const;
  
  virtual bool fixedNodalDofsBc(QList<int>& nodes, QList<int>&dof, QList<double>& values, bool* constantValues) const;

  virtual bool fixedNodalForcesBc(QList<int>& nodes, QList<int>&dof, QList<double>& values) const;

  virtual FemResultType fillElementDataForBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
                                             int border, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);

  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);

protected:

  enum xfemStateVarIds
  { 
    U_ID,         
    A_ID,         

    // ------      NO ADDING BELOW THIS LINE
    NUM_STATEVAR_IDS,  
  };

  enum xfemPhysicsAttributeIds
  {
    ISOPARAMETRIC_ID,        
        FIXEDENRICHEDDOF_ID,      


    // ------         NO ADDING BELOW THIS LINE
    NUM_ATTRI_IDS,    
  };


  enum xfemBoundaryConditionIds
  {                  
    FIXED_A_BC_ID, 


    // ------    NO ADDING BELOW THIS LINE
    NUM_BC_IDS,  
  };

  enum xfemBoundaryConditionValueIds
  {               
    BC_AX_ID,     
    BC_AY_ID,     
    BC_AZ_ID,     

    // ------     NO ADDING BELOW THIS LINE
    NUM_BCV_IDS,  
  };

  virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
  virtual bool checkAndLoadPrivateData(LuaTable& table);
  virtual void fillElementDisplacements(const GmElement* e, GmVector& ue, const int& nEnrich);
  void fillEnhancedNuJumpMatrix(const GmElement* e, GmMatrix& Nu, const GmVector& N, const QList<int>& BottomLevelSet, const QList<int>& TopLevelSet, const GmVector& nodes, GmElementDof* dofMap, XfemEnrichedElementData* exData, const int& iSubCrack);
  void fillEnhancedBuMatrix(const GmElement* e, const GmShape* shape, const GmVector& ncoord, const GmMatrix& X,
          const GmMatrix& J, GmMatrix& Bu, const QList<int>& gaussLevelSet, const int& nEnrich, const GmVector& nodes, GmElementDof* dofMap, XfemEnrichedElementData* exData, const int& gaussIP);
  virtual void fillEnhancedDof(const XfemElement* xe, QList<int>& iUe, const int& np);
  void tipDOFTreat(const int& nelem, GmVector& nodes);
  void fillNodalLevelSetMatrix(GmMatrix X, GmMatrix Xg, GmMatrix& H);
  void fillGaussLevelSetMatrix(GmMatrix X, GmMatrix Xg, GmMatrix& H);
  void fillInteractingCraksNodalLevelSets(const GmElement* e, XfemEnrichedElementData* exData, int nEnrich, int nCracks, QList<double>crackIndex, QList<GmMatrix>& nodalLevelSet);
  void fillInteractingCraksGaussLevelSets(const GmElement* e, XfemEnrichedElementData* exData, int nEnrich, int nCracks, QList<double>crackIndex, QList<GmMatrix>& gaussLevelSet);
  // computes the fracture rotation matrix
  virtual void fillFractureRotationMatrix(const GmElement* e, const GmShape* shp, GmMatrix Je, GmMatrix &Rot);
  //Return a list with the pair of ghost nodes of each fracture segment contained in this element
  virtual bool getCrackNodePair(XfemEnrichedElementData* exData, QList<QList<int>>& pfNodes);
  // remove repeted values
  virtual bool insertionSortWORepeatedint(QList<int>& a, double tol);
  // Sort vector values
  virtual bool insertionSortint(QList<int>& tVec);
  //Returns cossines of perpendicular unit vector of segment a-b
  void unitVectSegment(const GmMatrix& Xg, GmMatrix& uVec, GmMatrix& uVec2, GmMatrix& uVec3, GmMatrix& ROT);
  // Returns the elastic constitutive matrix for the fracture segment
  virtual bool fillNaturalCrackStiffnessMatrix(const GmElement* e, GmMatrix& Dep, GmVector coord, const int& ip, const GmVector& S_tn, const GmVector& Etn_old, const GmVector& v_w_tn, const GmMatrix& ROT, GmMatrix& D_xy, bool open);
  // Returns the Cohesive constitutive matrix for the fracture segment
  virtual bool fillCohesiveCrackStiffnessMatrix(const GmElement* e, GmMatrix& Dep, GmVector coord, const int& ip, const GmVector& nrd, const GmVector& trd, const GmVector& Etn_old, const GmVector& v_w_tn, const GmVector& sold,
          const GmMatrix& ROT, GmMatrix& D_xy, const double& rockStrength, bool open);
  // Returns fracture segment stresses
  virtual bool naturalFractureReturnMapping(const GmElement* e, XfemEnrichedElementData* exData, const GmVector& sold, const GmMatrix& Dep, const GmVector& enew, const GmVector& eold, GmVector& Stn_new, const GmMatrix& ROT, GmVector& S_w_xy, bool open, const int& iFracSubElement, const GmVector& ip, const int& k);
  // Returns fracture cohesive segment stresses
  virtual bool cohesiveFractureReturnMapping(const GmElement* e, const GmVector& sold, const GmMatrix& Dep, const GmVector& enew, const GmVector& eold, GmVector& Stn_new, const GmMatrix& ROT, GmVector& S_w_xy, const double& rockStrength, bool open, const GmVector& ip, const int& k);
  virtual bool assemblyCrackStiffnessMatrix(const GmElement* e, GmMatrix& elemK, const GmMatrix& Kfr, const QList<int>& iUe, const QList<int>& pfList, const int& numDof);
  virtual void assemblyCrackInternalForceVector(const GmElement* e, GmMatrix& elemFi, const GmVector& Fint_St, const QList<int>& iUe, const QList<int>& pfList, const int& numDof);
  // Split tensor to interface or Join interface to tensor
  virtual bool splitJoinTensorToInterface(GmVector& Tensor, GmVector& Cohesive, bool type);


  virtual GmGaussAccessor* createGaussAttributeAccessor(GmElementMesh* mesh, QString id, int snum, bool locked, 
                                                        Unit desiredUnit, const GmLogCategory& logger) const
  {
    XfemMesh* xmesh = dynamic_cast<XfemMesh*>(mesh);
    assert(xmesh);
    return xmesh->enrichedGaussAttributeAccessor(id, snum, locked, desiredUnit, logger);
  }

  const char*           _pluginType;  
  GmMatrix              _Dsn;         
  GmpMechanicalPhysics* _mechPhys;    
  GmpHydraulic*         _hydroPhys;   
  //GmSimulationData*   _SimulationData;   //!< Simulation data
};
#endif