gmpHydraulicMaterial.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_HYDRAULICMATERIAL_H_
#define _GEMA_PLUGIN_HYDRAULICMATERIAL_H_
 
#include "gmpHydraulic.h"  
 
#include <gmpFemPhysicsCommonMaterial.h>
#include <gmTrace.h>
#include <assert.h>
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpHydraulicMaterial : public GmpFemPhysicsCommonMaterial
{
//public:
protected:
  enum ElementPropertyIds
  {
    K_ID,          
    KWW_ID,        
    GW_ID,         
    Ki_ID,         
    rhof_ID,       
    Uf_ID,         
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
public:
 
  GmpHydraulicMaterial(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpFemPhysicsCommonMaterial(typeIndex, typeName, logger) {}
 
  virtual ~GmpHydraulicMaterial() {}
 
  virtual double waterBulkModulus(const GmElement* e, const GmVector* coord, int ip) const
  {
    if (propertyAc(KWW_ID) == NULL)
    {
      return 1.0e25;
    }
    return propertyAc(KWW_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double gammaWater(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(GW_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fluidDensity(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(rhof_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fluidViscosity(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(Uf_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double porosity(const GmElement* e, const GmVector* coord, int ip) const
  {
    Q_UNUSED(e); Q_UNUSED(coord); Q_UNUSED(ip);
    return 0.0; 
  }
 
  virtual double poreCompressibility(const GmElement* e, const GmVector* coord, int ip) const
  {
    Q_UNUSED(e); Q_UNUSED(coord); Q_UNUSED(ip);
    return 0.0;
  }
 
  virtual double referencePorePressure(const GmElement* e, const GmVector* coord, int ip) const = 0;
  
  virtual double relativePermeability (const GmElement* e, const GmVector* coord, int ip) const = 0;
  
  virtual void hydraulicConstitutiveModel(const GmElement* e, const GmVector* coord, int ip, const GmpHydraulic::HMaterialPointAccessors& mp, GmMatrix& K) const = 0;
  
  virtual void permeabilityStiffness(const GmElement* e, const GmVector* coord, int ip, GmMatrix& K) const;
  
  virtual bool permeabilityStiffnessIsSymmetric() const;
 
  virtual const QVariantMap* materialMetaDataMap();
};
#endif