xfemHydraulicEmbeddedFracture.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_HYDRAULIC_EMBEDDEDFRACTURE_H_
#define _GEMA_XFEM_HYDRAULIC_EMBEDDEDFRACTURE_H_
 
#include <gmpFemPhysicsCoupled.h>
#include <gmpHydraulic.h>
#include <gmDiscontinuitySet.h>
#include <gmDiscontinuity.h>
#include "gmDiscontinuityAccessor.h"
#include "gmPolylineDiscontinuity.h"
 
#include "xfemMesh.h"
 
class GmpHydraulicPhysics;
class GmGaussAccessor;
 
class XfemHydraulicEmbeddedFracture : public GmpFemPhysicsCoupled
{
public:
  // constructor
  XfemHydraulicEmbeddedFracture(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
    const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
  // destructor
  virtual ~XfemHydraulicEmbeddedFracture();
 
  // 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();
 
  // 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);
 
  virtual FemResultType fillEnrichedElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  virtual FemResultType fillEnrichedElementSubData(const GmElement* e, GmVector& pe, GmVector& Dpe, int nRegDofs, int nEnrDof, GmMatrix& Scc, GmMatrix& Scd, GmMatrix& Sdd, GmMatrix& Hcc, GmMatrix& Hcd, GmMatrix& Hdd, GmVector& qic, GmVector& qid);
 
protected:
 
  enum embeddedStateVarIds
  { 
    P_ID,         
    Pd_ID,        
 
    // ------      NO ADDING BELOW THIS LINE
    NUM_STATEVAR_IDS,  
  };
 
  enum embeddedPhysicsAttributeIds
  {
    ISOPARAMETRIC_ID,     
    CONDENSATION_ID,      
        CONDENSATION_TOL,         
 
    // ------         NO ADDING BELOW THIS LINE
    NUM_ATTRI_IDS,    
  };
 
 
  enum embeddedBC_Ids
  {                  
    FIXED_W_BC_ID, 
 
 
    // ------    NO ADDING BELOW THIS LINE
    NUM_BC_IDS,  
  };
 
  enum embeddedBCValueIds
  {               
    BC_WX_ID,     
    BC_WY_ID,     
    BC_WZ_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 fillElementPressureJumps(const GmElement* e, GmVector& Dpe);
  virtual void oldElementDisplacementJumps(const GmElement* e, GmVector& we);
  virtual void fillDiscontinuitySetElementData(const GmElement* e, GmVector& we, GmMatrix& Kd, GmVector& fid);
  virtual void assembleElementData(bool needsK, bool needsFi, GmMatrix& elemK, GmMatrix& Kuu, GmMatrix& Kuw, GmMatrix& Kww, GmVector& elemFi, GmVector& fiu, GmVector& fiw);
  virtual void get2DDiscontinuityRotationMatrix(const GmCompact2DSegmentCellIntersection* discontinuitySeg, GmMatrix& R);
  virtual void discontinuityCohesiveLaw(double ks, double kn, GmVector& w, GmVector& td, GmMatrix& Td);
  virtual void fillPressureJumpMappingMatrix(const GmElement* e, GmMatrix& Xe, GmMatrix& Mp);
  virtual void get2DDiscontinuityGeometry(const GmCompact2DSegmentCellIntersection* discontinuitySeg, GmVector& Xref, GmVector& m, GmVector& nd);
  virtual void fillEnrichedBuMatrix(const GmElement* e, const GmShape* elemShape, const GmVector& Xn, GmMatrix& Xe, GmMatrix& Mu, GmMatrix& Bu, GmMatrix& Buenr);
  double heavisideFunction(GmVector& X, GmVector& Xref, GmVector& nd);
  int getNumberEnrichedDofs(const GmElement* e);
  virtual void fillHeavisideMatrix(GmMatrix& Hd, GmMatrix& Xe, GmVector& Xref, GmVector& nd);
 
  virtual double fillBpMatrix(const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bp)
  {
          S_TRACE();
          Q_UNUSED(X);  Q_UNUSED(N);
          return shape->shapeCartesianPartialsFromJacobian(ncoord, J, Bp, true);
  }
 
  const char*           _pluginType;   
  GmMatrix              _Dsn;          
  GmpHydraulic*             _hydroPhys;    
  //GmSimulationData*   _SimulationData;   //!< Simulation data
 
private:
  GmDiscontinuitySet*   _discontinuitySet;
};
#endif