gmpHydraulicEmbeddedFractureConductive.h

/************************************************************************
**
** 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_CONDUCTIVE_H_
#define _GEMA_PLUGIN_HYDRAULICS_EMBEDDED_FRACTURE_CONDUCTIVE_H_
 
 
#include "gmpHydraulicConfig.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 GmpHydraulicEmbeddedFractureConductive : public GmpHydraulic
{
public:
  // constructor
        GmpHydraulicEmbeddedFractureConductive(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
  // destructor
  virtual ~GmpHydraulicEmbeddedFractureConductive();
 
  // 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 bool checkAndLoadDiscontinuitySetAccessors(LuaTable& table);
 
  // Comments on the base class
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // Fills nodal pressure of the discontinuity elements inside element e
  virtual bool fillDiscontinuityInternalPressureNodalValues(GmElement* e, GmValueAccessor* nAc);
 
  // Sets the initial conditions for hydraulic embedded fracture models.
  virtual bool setInitialConditions();
 
  // calculate derived results
  virtual bool calcDerivedResults(bool nonLinear);
 
  // Assemble the contribution of the discontinuity element
  void assembleDiscontinuityElementMatrix(GmMatrix& Kg, const GmMatrix& Kl, const  QVector<int>& dof) const;
 
  // Fills the element matrices and vectors for a enriched element.
  virtual FemResultType fillEnrichedElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors, int nDiscontinuities);
  
  // Fills the element matrices and vectors of the discontinuity elements.
  FemResultType fillDiscontinuitySetElementData(const GmElement* e, GmMatrix& Hd, GmMatrix& Sd, GmMatrix& Tcd, const QMap<int, int>& dNodesGlobalToLocal);
 
  // Fills the element matrices and vectors of the continuum (Used by the hydro-mechanical physics).
  FemResultType fillContinuumElementData(const GmElement* e, GmMatrix& Hcc, GmMatrix& Scc);
private:
  const char*       _pluginType;     
};
 
#endif