gmpDualPorosityDualPermeability.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_DUALPOROSITYDUALPERMEABILITY_H_
#define _GEMA_PLUGIN_DUALPOROSITYDUALPERMEABILITY_H_
 
 
#include "gmpCoupledHMFemPhysics.h"
#include <gmpFemPhysicsCoupled.h>
#include <gmTrace.h>
 
class GmGaussAccessor;
 
 
class GmpDualPorosityDualPermeability : public GmpCoupledHMFemPhysics
{
public:
  // constructor
  GmpDualPorosityDualPermeability(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                                  const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
 
  // destructor
  virtual ~GmpDualPorosityDualPermeability();
 
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  virtual bool fillCoupledMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& R1, const GmVector& Np, const GmMatrix& Dmfr, const GmMatrix& Cm, const GmMatrix& Bu, double c);
  virtual bool fillCoupledFracture(const GmElement* e, const GmVector* coord, int ip, GmMatrix& R2, const GmVector& Np, const GmMatrix& Dmfr, const GmMatrix& Cm, const GmMatrix& Bu, double c);
  virtual void fillMechanicalStiffnessMatrix(const GmElement* e, GmMatrix& Ke, const GmMatrix& Dmfr, const GmMatrix& B, double c);
  virtual void fillInternalForceVector(const GmElement* e, GmVector& Fi, const GmMatrix& Kc, const GmMatrix& Kfc, const GmMatrix& Wfm, const GmMatrix& R1, const GmMatrix& R2, const GmVector& pe, const GmVector& pfe, const GmVector& pfu, const QVector<int>& iUe, const QVector<int>& iPe, const QVector<int>& iPfe);
  virtual void fillDampingMatrix(const GmElement* e, GmMatrix& C, const GmMatrix& Sc, const GmMatrix& Sfc, const GmMatrix& R1, const GmMatrix& R2, const QVector<int>& iUe, const QVector<int>& iPe, const QVector<int>& iPfe);
  virtual void assemblyStiffnessMatrix(const GmElement*e, GmMatrix& elemK, const GmMatrix& Ke, const GmMatrix& R1, const GmMatrix& R2, const GmMatrix& Wfm, const GmMatrix& Kc, const GmMatrix& Kfc, const QVector<int>& iUe, const QVector<int>& iPe, const QVector<int>& iPfe);
  virtual void FillDofIndexPosition(const GmElement* e, QVector<int>& iUe, QVector<int>& iPe, QVector<int>& iPfe);
  virtual bool calcDerivedResults(bool nonLinearSolver);
  virtual FemResultType fillElementDataForBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex, int border, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // See comments on the base class
  virtual bool addStateItemsToGroup(GmStateDump* state, bool fixedHint, int groupId);
 
protected:
  // See comments on the base class
  virtual const GmElementDof* dofMapping(GmCellType type) const;
 
  enum HMDPDKVarIds
  {                    
    PF_ID = GmpCoupledHMFemPhysics::NUM_STATEVAR_IDS,
 
    // ------     NO ADDING BELOW THIS LINE
    NUM_STATEVAR_IDS,  
  };
 
};
 
#endif