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