gmpPorousElasticMaterial.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_POROUS_ELASTIC_MATERIAL_H_
#define _GEMA_PLUGIN_POROUS_ELASTIC_MATERIAL_H_
 
#include <gmpFemPhysicsCommonMaterial.h>
#include <gmpCoupledHMMaterial.h>
#include <gmpHydraulicMaterial.h>
 
#include "gmpMechanicalMaterial.h"
#include "gmpMaterialElastic.h"
#include "gmpMechanicPoint.h"
 
class GmpPorousElasticMaterial : public GmpCoupledHMMaterial
{
protected:
  enum ElementPropertyIds
  {
    KSS_ID,         
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
public:
 
  GmpPorousElasticMaterial(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpCoupledHMMaterial(typeIndex, typeName, logger)
  {
  }
 
  virtual ~GmpPorousElasticMaterial() {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpPorousElasticMaterial(typeIndex, typeName, logger);
  }
 
  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap();
 
  virtual double solidBulkModulus(const GmElement* e, const GmVector* coord, int ip) const
  {
    if (propertyAc(KSS_ID) == NULL)
    {
      return 1.0e25;
    }
        // read solid 
        double Kss = propertyAc(KSS_ID)->scalarValueAt(e, coord, ip);
        if (Kss <= 0.0) { Kss = 1.0e25; }
 
        return Kss;
  }
 
  // Returns the bulk modulus of the porous skeleton
  virtual double porousBulkModulus(const GmElement* e, const GmVector* coord, int ip) const
  {
      const GmpMechanicalMaterialElastic* mmat = dynamic_cast<const GmpMechanicalMaterialElastic*>(mechanical());
      assert(mmat);
      return mmat->bulkModulus(e, coord, ip);
  }
 
  // Returns the bulk modulus of the porous skeleton
  virtual double biotCoeffient(const GmElement* e, const GmVector* coord, int ip) const
  {
      const GmpMechanicalMaterialElastic* mmat = dynamic_cast<const GmpMechanicalMaterialElastic*>(mechanical());
      assert(mmat);
      double Kb = mmat->bulkModulus(e, coord, ip);
          double Ks = solidBulkModulus(e, coord, ip);
          double alpha = 1.0 - (Kb / Ks);
      return alpha;
  }
  
//private:
  
};
 
#endif