xfemMechanicEmbeddedFracture.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_EMBEDDEDFRACTURE_H_
#define _GEMA_XFEM_MECHANIC_EMBEDDEDFRACTURE_H_
 
#include <gmpFemPhysicsCoupled.h>
#include <gmpMechanicalPhysics.h>
#include <gmpMechanicPoint.h>
#include <gmDiscontinuitySet.h>
#include <gmDiscontinuity.h>
#include "gmDiscontinuityAccessor.h"
#include "gmPolylineDiscontinuity.h"
 
#include "xfemMesh.h"
 
class GmpMechanicalPhysics;
class GmGaussAccessor;
 
template <class T, class F> class embeddedPhysicsWrapper : public T
{
public:
  embeddedPhysicsWrapper(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 XfemMechanicEmbeddedFracture : public GmpFemPhysicsCoupled
{
public:
  // constructor
  XfemMechanicEmbeddedFracture(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
    const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
  // destructor
  virtual ~XfemMechanicEmbeddedFracture();
 
  // 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(QVector<int>& nodes, QVector<int>&dof, QVector<double>& values, bool* constantValues) const;
 
  virtual bool fixedNodalForcesBc(QVector<int>& nodes, QVector<int>&dof, QVector<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& ue, GmVector& we, int nRegDofs, int nEnrDof, GmMatrix& Kuu, GmMatrix& Kuw, GmMatrix& Kww, GmVector& fiu, GmVector& fiw);
 
  virtual FemResultType solveDiscontinuityEquilibriumEq(const GmElement* e, GmVector& ue, GmVector& we, int nRegDofs, int nEnrDof, GmMatrix& Kuu, GmMatrix& Kuw, GmMatrix& Kww, GmVector& fiu, GmVector& fiw);
 
protected:
 
  enum embeddedStateVarIds
  { 
    U_ID,         
    W_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,  
  };
 
  // See comments on the base class
  virtual bool loadPrivateData(LuaTable& table);
 
  virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
  virtual bool checkAndLoadPrivateData(LuaTable& table);
  virtual void fillElementDisplacementJumps(const GmElement* e, GmVector& we);
  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 displacementJumpShapeFunction(GmVector& N, GmMatrix& Nd);
  virtual void discontinuityCohesiveLaw(double ks, double kn, GmVector& w, GmVector& td, GmMatrix& Td);
  virtual void addTranslationDeformationMode(GmMatrix& Mu);
  virtual void addRelativeRotationDeformationMode(GmMatrix& Mu, GmVector& X, GmVector& Xr, GmVector& n, GmVector& m);
  virtual void transformSpaceDisplacementJumpMappingMatrix(GmMatrix& Mu, const GmCompact2DSegmentCellIntersection* discontinuitySeg);
  virtual void fillDisplacementJumpMappingMatrix(const GmElement* e, GmMatrix& Xe, GmMatrix& Mu);
  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);
 
  const char*           _pluginType;  
  GmMatrix              _Dsn;         
  GmpMechanicalPhysics* _mechPhys;    
  //GmSimulationData*   _SimulationData;   //!< Simulation data
 
private:
  GmDiscontinuitySet*   _discontinuitySet;
};
#endif