TPZMatPorous.h

Go to the documentation of this file.

#ifndef PZPORO_H
#define PZPORO_H


#include "TPZMaterial.h"
#include "TPZMatTemporal.h" 
#include "pzporoelastoplasticmem.h"
#include "TPZMatElastoPlastic.h"
#include "pzbndcond.h"

#define TBASEPOROUS(T, TMEM) TPZMatElastoPlastic< T, TMEM >

template <class T, class TMEM = TPZPoroElastoPlasticMem >
class TPZMatPorous : public TPZMatTemporal, public TPZMatElastoPlastic< T, TMEM >
{
   public:

      enum SOLUTIONVARS{ENone = -1,
                                                                                                EPorePressure = 90};

      TPZMatPorous();
                                
      TPZMatPorous(int id);

      TPZMatPorous(const TPZMatPorous<T, TMEM > &mat);
      virtual ~TPZMatPorous();
                
      virtual std::string Name() override;
                
                  void SetUp(const REAL &k, const REAL &Mu, 
                                                                                                                                const REAL &StorageEps,
                                                                                                                                const REAL &Alpha,
                                                                                                                                const REAL &Rhof);

      virtual int Dimension() const override { return TBASEPOROUS(T, TMEM)::Dimension(); }

      virtual int NStateVariables() const override{ return TBASEPOROUS(T, TMEM)::NStateVariables() + 1; }

      virtual void Print(std::ostream &out = std::cout, const int memory = 0) override;

      virtual int VariableIndex(const std::string &name) override;

      virtual int NSolutionVariables(int var) override;

      virtual void Solution(TPZMaterialData &data, int var, TPZVec<REAL> &Solout) override;
                
      virtual int NFluxes()  override
                  {
         PZError << "TPZMatPorous<TBASEPOROUS(T, TMEM)>::NFluxes() - Method not implemented\n";
         return 0;
      }

      virtual void Flux(TPZVec<REAL> &x, TPZVec<REAL> &Sol, TPZFMatrix<REAL> &DSol, TPZFMatrix<REAL> &axes, TPZVec<REAL> &flux) override
                  {
         PZError << "TPZMatPorous<TBASEPOROUS(T, TMEM)>::Flux - Method not implemented\n";
      }

      virtual void Errors(TPZVec<REAL> &x,TPZVec<REAL> &u, TPZFMatrix<REAL> &dudx,
                          TPZFMatrix<REAL> &axes, TPZVec<REAL> &flux,
                          TPZVec<REAL> &u_exact,TPZFMatrix<REAL> &du_exact,TPZVec<REAL> &values) override;
      virtual int NEvalErrors() override {return NStateVariables();}

      virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<REAL> &ek, TPZFMatrix<REAL> &ef) override;

      virtual void ContributeBC(TPZMaterialData &data, REAL weight, TPZFMatrix<REAL> &ek, TPZFMatrix<REAL> &ef, TPZBndCond &bc) override;

      virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<REAL> &ef) override;

      virtual void ContributeBC(TPZMaterialData &data, REAL weight, TPZFMatrix<REAL> &ef, TPZBndCond &bc) override;

      virtual TPZMaterial * NewMaterial() override;

      public:
virtual int ClassId() const override;


      virtual void Write(TPZStream &buf, int withclassid) const override;

      virtual void Read(TPZStream &buf, void *context) override;
                
                  virtual void FillDataRequirements(TPZMaterialData &data) override;
                
                  void ComputePorePressure(TPZMaterialData & data, REAL & Pp, TPZVec<REAL> & dPp);
                
                  void UpdatePorePressure(TPZMaterialData & data);
                
                  void SetPorePressure(const REAL Pp);
                
protected:
                
                  REAL fk;
                
                  REAL fMu;
                
                  REAL fStorageEps;
                
                  REAL m_hardening;
                
                  REAL fRhof;
                
};

template <class T, class TMEM>
int TPZMatPorous<T, TMEM >::ClassId() const{
    return Hash("TPZMatPorous") ^ TPZMatTemporal::ClassId() << 1 ^ TPZMatElastoPlastic<T, TMEM>::ClassId() << 2;
}
#endif