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(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& 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