gmpMechanicalEmbeddedFracture.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** 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_PLUGIN_MECHANICAL_EMBEDDED_FRACTURE_H_
#define _GEMA_PLUGIN_MECHANICAL_EMBEDDED_FRACTURE_H_

#include "gmpMechanicalPhysics.h"

#include <gmpFemPhysicsCommon.h>
#include <gmpFemPhysicsData.h>
#include <gmDiscontinuitySet.h>
#include <gmDiscontinuity.h>
#include "gmDiscontinuityAccessor.h"
#include "gmPolylineDiscontinuity.h"
#include <gmTrace.h>

class GmGaussAccessor;

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMechanicalEmbeddedFracture: public GmpMechanicalSolid
{
public:

  // constructor
  GmpMechanicalEmbeddedFracture(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                     const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);

  // destructor
  virtual ~GmpMechanicalEmbeddedFracture();

  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
  bool numberEnhancedDofs(const GmElement* e, int & ndofw);
  //virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  // Fills fracture contribution
  virtual FemResultType fillEnrichedElementData(const GmElement* e, GmMatrix& K, GmVector& Fi, bool needsK, bool needsFi);
  void initializeFractureGeometry(double & ld, GmVector & m1, GmVector & m2, GmVector & nd, GmVector & Xref, GmMatrix & Xf, int nfn);
  void elementStiffnessSubMtrcs(const GmElement * e, const GmElement * ef, GmVector & we, GmMatrix & Kaa, GmMatrix & Kaw, GmMatrix & Kwa, GmVector & Fi, bool needsFi);
  void fillBulkContribution(const GmElement * e, const GmElement * ef, GmVector & we, GmMatrix & Kww, GmVector & Fint);
  void fillGrMtrx(const GmElement * e, GmMatrix & Gr, GmMatrix & Bu, GmVector & Xip);
  void elementStiffnessSubMtrx_Kaw(GmMatrix & Kaa, GmMatrix & Kaw, GmMatrix & Xe, const GmElement * e, const GmElement * ef);
  //void elementStiffnessSubMtrx_Kaw(GmMatrix & Kaa, GmMatrix & Kaw, GmMatrix & Xe, GmVector & m1, GmVector & nd, GmVector & Xref, double ld, int d, int ne, int nfn);
  void elementStiffnessSubMtrx_Kwa(GmMatrix & Kaa, GmMatrix & Kwa, GmMatrix & Xe, GmVector & m1, GmVector & nd, GmVector & Xref, double ld, int d, int ne, int nfn);
  void elementStiffnessSubMtrx_Kww(GmMatrix & Kaa, GmMatrix & Kww, GmMatrix & Xe, GmVector & m1, GmVector & nd, GmVector & Xref, double ld, int d, int ne, int nfn);
  void computeJumpDisplacements(const GmElement * e, GmVector & we, GmMatrix & Kww);
  void fillFractureElementData(const GmElement * ef, GmVector & strain, GmMatrix & Kd, GmVector & Fdi);
  void updateFractureDisplacementJumps(const GmElement * ef, GmVector & we);
  void getNodalDisplacementJumps(const GmElement * ef, GmVector & we);
  double heavisideFnc(GmVector & X, GmVector & Xref, GmVector & nd);
  void elementHeavisideMtrx(GmMatrix & Hd, GmVector & Xref, GmMatrix & Xe, GmVector & nd, int ne);
  void rotationMtrx(const GmElement * ef, GmMatrix & Rot);
  void mappingTranslationMtrx(GmMatrix & Mtr);
  void mappingRelativeRotationMtrx(GmMatrix & Mrot, GmVector & X, GmVector & m1, GmVector & Xref, double ld);
  void mappingRigidBodyMtrx(GmMatrix & Mrb, GmVector & X, GmVector & m1, GmVector & Xref, double ld, int nfn);
  void elementMappingMtrx(GmMatrix & Me, GmVector & m1, GmMatrix & Xe, GmVector & Xref, double ld, int nfn, int ne);
  void fillMiddleCoordinates(const GmElement * e, GmMatrix & X, GmMatrix & MX);
  void constructShapeFncMtrx(GmMatrix & Nw, GmVector & N);
  // See comments on the base class
  virtual bool calcDerivedResults(bool nonLinearSolver);
  // Comments on the base class
  //virtual bool dofByElement(bool* fixed, bool* addOnly, bool* trackChanges) const;
  // Comments on the base class
  virtual const GmElementDof* dofMapping(const GmElement* e) const;

protected:
  enum StateVarIds
  { 
    W_ID= GmpMechanicalSolid::NUM_STATEVAR_IDS, 

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

  enum embeddedPhysicsAttributeIds
  {
    CONDENSATION_ID = GmpMechanicalSolid::NUM_ATTRI_IDS,  

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

  // See comments on the base class
  virtual bool loadPrivateData(LuaTable& table);

private:
  const char*       _pluginType;     
  GmpFemPhysicsData _fractureData;   
};

#endif