gmpHydraulicEmbeddedFracture.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_HYDRAULICS_EMBEDDED_FRACTURE_H_
#define _GEMA_PLUGIN_HYDRAULICS_EMBEDDED_FRACTURE_H_
 
 
#include "gmpHydraulicConfig.h"
#include "gmpHydraulicEmbeddedFractureConductive.h"
 
#include <gmpFemPhysicsCommon.h>
#include <gmpFemPhysicsData.h>
#include <gmpHydraulic.h>
#include <gmGaussAccessor.h>
#include <gmMathUtils.h>
#include <gmDiscontinuitySet.h>
#include <gmDiscontinuity.h>
#include "gmDiscontinuityAccessor.h"
 
class GmGaussAccessor;
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpHydraulicEmbeddedFracture : public GmpHydraulicEmbeddedFractureConductive
{
public:
  // constructor
  GmpHydraulicEmbeddedFracture(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
  // destructor
  virtual ~GmpHydraulicEmbeddedFracture();
 
  // Comments on the base class
  virtual const char* pluginType() const { return _pluginType; }
 
  // Comments on the base class
  virtual const QVariantMap* physicsMetaDataMap();
 
  // Comments on the base class
  bool checkAndLoadAttributeAccessors(LuaTable& nodeTable, LuaTable& cellTable, LuaTable& gaussTable);
 
  // Comments on the base class
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // Fills the element matrices and vectors for a enriched element.
  virtual FemResultType fillEnrichedElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors, int nDiscontinuities);
 
  // Discretizes the system of equations in time using a fully implicit scheme.
  virtual FemResultType fullyImplicitTimeDiscretization(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // Fills nodal pressure of the discontinuity elements inside element e
  virtual bool fillDiscontinuityInternalPressureNodalValues(GmElement* e, GmValueAccessor* nAc);
 
  // Fills degrees of freedom
  void fillElementOldPorePressure(const GmElement* e, GmVector& pe);
  void fillElementPorePressureJump(const GmElement* e, GmVector& Dp);
  void fillElementOldPorePressureJump(const GmElement* e, GmVector& dpe);
  void setElementPorePressureJump(const GmElement* e, GmVector& dpe);
  void fillDiscontinuitySetInternalPressure(const GmElement* e, GmVector& pd, QMap<int, int> dNodesGlobalToLocal);
  void fillDiscontinuitySetOldInternalPressure(const GmElement* e, GmVector& pdOld, QMap<int, int> dNodesGlobalToLocal);
  void setDiscontinuitySetInternalPressure(const GmElement* e, GmVector& pd, QMap<int, int> dNodesGlobalToLocal);
  
  // Fills the element matrices and vectors of the enriched element that considers a discontinuous pressure field.
  FemResultType  fillDiscontinuousPressureElementSubData(const GmElement* e, bool needsK, bool needsC, bool needsFi, GmMatrix& Hcc, GmMatrix& Hcj, GmMatrix& Hjc, GmMatrix& Hjj, GmMatrix& Scc, GmMatrix& Scj, GmMatrix& Sjj, QMap<int, int> dNodesGlobalToLocal, int nDeltaPDofs);
  
  // Fills the element compressibility matrices for a enriched element with discontinuous pressure field.
  void fillCompressibilityMatrices(const GmElement* e, GmMatrix& Xe, const GmVector* coord, int ip, GmMatrix& Scc, GmMatrix& Scj, GmMatrix& Sjj, const GmVector& N, double c, int nDeltaPDofs);
 
  // Fills the element fluid-flow matrices for a enriched element with discontinuous pressure field.
  void fillFluidFlowMatrices(const GmElement* e, GmMatrix& Xe, const GmVector* coord, int ip, GmMatrix& Hcc, GmMatrix& Hcj, GmMatrix& Hjc, GmMatrix& Hjj, const GmMatrix& Bp, double c, int nDeltaPDofs);
  
  // Enriched shape function associated with the discontinuity idDiscontinuityElem
  void fillEnrichedShapeValues(const GmElement* e, GmMatrix& Xe, GmVector& N, GmVector& X, int idDiscontinuityElem, GmVector& Nenr);
 
  // Enriched shape function matrix
  void fillNenrMatrix(const GmElement* e, GmMatrix& Xe, GmVector& N, GmVector& X, GmVector& Nenr);
  
  // Enriched gradient pressure matrix
  void fillGMatrix(const GmElement* e, GmMatrix& Xe, const GmMatrix& B, GmMatrix& G);
 
  // Gets the number of enriched dofs (pressure jumps) at the element
  int getNumberEnrichedDofs(const GmElement* e);
 
  // Fills the element matrices and vectors of the discontinuity elements.
  FemResultType fillDiscontinuitySetElementData(const GmElement* e, GmMatrix& Hd, GmMatrix& Sd, GmMatrix& Lcc, GmMatrix& Lcj, GmMatrix& Lcd, GmMatrix& Ljj, GmMatrix& Ljd, GmMatrix& Ldd, GmMatrix& Tcd, GmMatrix& Tjd, QMap<int, int> globalToLocal);
 
  enum embeddedPhysicsAttributeIds
  {
          DP_ID,                                        // Current pressure jump
          DPOLD_ID,                                     // Old pressure jump
 
          // ------         NO ADDING BELOW THIS LINE
          NUM_ATTRI_IDS,    
  };
 
private:
  const char*       _pluginType;     
  GmValueAccessor* _pOldAcc;             
};
 
#endif