gmpHMEmbeddedFractureConductive.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_HYDROMECHANICAL_EMBEDDED_FRACTURE_CONDUCTIVE_H_
#define _GEMA_PLUGIN_HYDROMECHANICAL_EMBEDDED_FRACTURE_CONDUCTIVE_H_
 
#include "gmpHydroMechanical.h"
#include "gmpCoupledHMFemPhysics.h"
#include <gmpFemPhysicsCoupled.h>
#include <gmpHydraulicEmbeddedFracture.h>
#include <gmpMechanicalSolid.h>
#include <gmpMechanicalEmbeddedFracture.h>
#include <gmTrace.h>
#include <gmDiscontinuitySet.h>
#include <gmDiscontinuity.h>
#include "gmDiscontinuityAccessor.h"
 
class GmGaussAccessor;
 
class GmpHMEmbeddedFractureConductive : public GmpHydroMechanical
{
public:
  // constructor
    GmpHMEmbeddedFractureConductive(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                     const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
 
  // destructor
  virtual ~GmpHMEmbeddedFractureConductive();
 
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
 
  // Comments on the base class
  bool checkAndLoadPrivateData(LuaTable& table);
 
  // Comments on the base class
  bool checkAndLoadAttributeAccessors(LuaTable& nodeTable, LuaTable& cellTable, LuaTable& gaussTable);
 
  // Fill the element old displacement vector
  void fillElementOldDisplacement(const GmElement* e, GmVector& ur);
 
  // Fill the element old pressure vector
  void fillElementOldPorePressure(const GmElement* e, GmVector& pe);
 
  // Comments on the base class
  virtual bool checkAndLoadDiscontinuitySetAccessors(LuaTable& table);
 
  // Sets the initial conditions for embedded fracture models.
  virtual bool setInitialConditions();
 
  // Comments on the base class
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // Discretizes the system of equations in time using a fully implicit scheme.
  virtual FemResultType fullyImplicitTimeDiscretization(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors, GmMatrix& elemK, GmMatrix& elemC, GmVector& elemFi);
 
  // Compute the element matrices and vectors for the elements with embedded discontinuities
  virtual FemResultType fillEnrichedElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // Assemble the element matrix and vector 
  void  assembleElementData(const GmElement* e, GmMatrix& A, GmVector& b, GmMatrix& Al, GmVector& bl, bool needsLHS, bool needsRHS);
  void  assembleElementDofVector(const GmElement* e, GmVector& dof, GmMatrix& dof_u, GmVector& dof_p);
 
  // Solves the local equations associated with the enriched dofs
  FemResultType solveLocalEquations(const GmElement* e, GmMatrix& A_rr, GmMatrix& A_re, GmMatrix& A_er, GmMatrix& A_ee, GmVector& b_r, GmVector& b_e, GmVector& dof_e);
 
  // Fill the system sub-matrices and vectors for the enriched element
  FemResultType fillSystemSubData(const GmElement* e, GmMatrix& A_rr, GmMatrix& A_re, GmMatrix& A_er, GmMatrix& A_ee, GmVector& b_r, GmVector& b_e, GmVector& dof_e, GmMatrix& stressCoeffs, QMap<int, int> dNodesGlobalToLocal);
 
  // Fill the system sub-matrices and vectors for the enriched element
  FemResultType fillElementSubData(const GmElement* e, QMap<QString, GmMatrix>& sub_matrices, QMap<QString, GmVector>& sub_vectors, GmVector& dof_ur, GmVector& dof_ue, GmVector& dof_pc, GmMatrix& stressCoeffs, QMap<int, int> dNodesGlobalToLocal);
 
  // Initialize the element sub-matrices and sub-vectors
  void initializeSubMatrices(QMap<QString, GmMatrix>& sub_matrices, QMap<QString, GmVector>& sub_vectors, int ndofs_ur, int ndofs_ue, int ndofs_pc, int ndofs_pd);
 
  // Update the discontinuity elements hydromechanical coupling parameters
  FemResultType updateDiscontinuitySetHMCouplingParameters(const GmElement* e);
  bool updateHMCouplingParameters(const GmElement* e, const GmVector* coord, int ip);
  bool updateAperture(const GmElement* e, const GmVector* coord, int ip);
 
  // Fill the coupling hydromechanical sub-matrices
  FemResultType fillHMCouplingData(const GmElement* e, GmMatrix& Quc, GmMatrix& Qac, GmMatrix& Qad, GmMatrix& Qca, GmMatrix& stressCoeffs, QMap<int, int> dNodesGlobalToLocal);
  FemResultType fillContinuumHMCouplingData(const GmElement* e, GmMatrix& Quc, GmMatrix& Qac, GmMatrix& Qca, GmMatrix& stressCoeffs);
  FemResultType fillDiscontinuitySetHMCouplingData(const GmElement* e, GmMatrix& Qad, QMap<int, int> dNodesGlobalToLocal);
 
  enum embeddedPhysicsAttributeIds
  {
          UPDATE_DSET_APERTURE = GmpHydroMechanical::NUM_ATTRI_IDS,       
 
          // ------         NO ADDING BELOW THIS LINE
          NUM_ATTRI_IDS,    
  };
protected:
    GmpHydraulicEmbeddedFractureConductive*        _embeddedHydroPhys;      
    GmpMechanicalEmbeddedFracture*                 _embeddedMechPhys;       
private:
    GmValueAccessor* _uOldAcc;           
    GmValueAccessor* _pOldAcc;           
};
 
#endif