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 <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 checkAndLoadAttributeAccessors(LuaTable& nodeTable, LuaTable& cellTable, LuaTable& gaussTable);
  //virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  virtual FemResultType fillEnrichedElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  int getNumberEnrichedDofs(const GmElement* e);
  int getNumberEnrichedDofsPerDiscontinuity() const;
  int getDisplacementJumpOrder() const;
  FemResultType solveDiscontinuityEquilibriumEq(const GmElement* e, int nRegDofs, int nEnrDof, GmMatrix& Kuu, GmMatrix& Kuw, GmMatrix& Kwu, GmMatrix& Kww, GmVector& fiu, GmVector& fia, GmVector& ue, GmVector& ae);
  FemResultType fillEnrichedElementSubData(const GmElement* e, const GmVector& ue, const GmVector& we, const GmMatrix& stressCoeffs, GmMatrix& Kuu, GmMatrix& Kuw, GmMatrix& Kwu, GmMatrix& Kww, GmVector& fiu, GmVector& fiw);
  FemResultType fillDiscontinuitySetElementData(const GmElement* e, const GmVector& ue, const GmVector& ae, GmMatrix& Kd, GmVector& fid, GmMatrix& Kstab) const;
  FemResultType  fillInternalBorderElementData(const GmElement* e, GmVector& ue, GmVector& ae, GmMatrix& Kb_uu, GmMatrix& Kb_ua, GmMatrix& Kb_au, GmMatrix& Kb_aa, GmVector& fbu, GmVector& fba);
  void getDiscontinuitySubElementIds(const GmElement* e, int k, int& idDiscontinuityElem, int& idDiscontinuityMesh);
  void displacementJumpShapeFunction(const GmMatrix& DX, const GmMatrix& mn, GmMatrix& Nd) const;
  void fillGrMatrix(const GmElement* e, const GmMatrix& Xe, const GmMatrix& Bu, GmMatrix& Gc) const;
  void fillGvMatrix(const GmElement* e, const GmShape* elemShape, const GmVector& Xn, const GmMatrix& Xe, const GmMatrix& C, GmMatrix& G) const;
  void projectionMatrix(const GmVector& nd, GmMatrix& P) const;
  int dimPolynomialStressInterpolation(const GmElement* e) const;
  void polynomialStressInterpolation(int etype, const GmVector& X, GmMatrix& p) const;
  void GramMatrixPolynomialStress(const GmElement* e, int stressOrder, GmMatrix& H) const;
  void discontinuityIntegralPolynomialStress(const GmElement* ef, int displacementJumpOrder, int stressOrder, GmMatrix& S) const;
  double lineSearch(const GmElement* e, GmMatrix& stressCoeffs, GmVector& ue, GmVector& x0, GmVector& d, GmVector& r);
  FemResultType fillEnrichedElementLocalEquation(const GmElement* e, const GmMatrix& stressCoeffs, const GmVector& ue, const GmVector& ae, GmVector& fia);
  void fillElementOldEnrichedDof(const GmElement* e, GmVector& ae);
  void setElementEnrichedDof(const GmElement* e, GmVector& ae);
  void fillEnrichedShapeFncMatrix(const GmElement* e, const GmMatrix& Xe, const GmVector& N, GmMatrix& Nenr) const;
  double auxiliaryFunctionPhi(const GmVector& N, const GmMatrix& Xe, const GmVector& Xref, const GmVector& nd) const;
  void fillDisplacementShapeFncMatrix(const GmElement* e, const GmVector& Nv, GmMatrix& Nm) const;
  void fillStressInterpolationCoeffs(const GmElement* e, GmMatrix& stressCoeffs) const;
  FemResultType fillElementDataForLoads(const GmElement* e, int loadId, GmpFemVectorSet& elemVectors);
  void fillStabilizationMatrix(const GmElement* e, const GmElement* ef, const GmMatrix& Td, const GmMatrix& Dep, GmMatrix& S) const;
 
  enum embeddedPhysicsAttributeIds
  {
          LOCALEQ_TOL = GmpMechanicalSolid::NUM_ATTRI_IDS,        
          LINESEARCH,         // Flag to use the line search method during the local equilibrium equation solution
          RELROTATIONMODE,    // Flag to add the relative rotation modes
          STRETCHINGMODE,     // Flag to add the stretching mode
          SYMFORMULATION,     // Flag to use the symmetric formulation
 
          // ------         NO ADDING BELOW THIS LINE
          NUM_ATTRI_IDS,    
  };
 
  enum embeddedCellAttributeIds
  {
          A_ID,                                         // Current enrichment dof vector associated with a discontinuity
          AOLD_ID,                                      // Old enrichment dof vector associated with a discontinuity            
 
          // ------         NO ADDING BELOW THIS LINE
          NUM_CELLATTRI_IDS,    
  };
 
private:
  const char*       _pluginType;     
};
 
#endif