gmpCoupledDualPorosityMaterial.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_COUPLED_DUALPOROSITY_MATERIAL_H_
#define _GEMA_PLUGIN_COUPLED_DUALPOROSITY_MATERIAL_H_
 
#include <gmpFemPhysicsCommonMaterial.h>
 
#include <gmpMechanicalMaterial.h>
#include <gmpCoupledHMMaterial.h>
#include <gmpHydraulicDualFlow.h>
#include <gmpHydraulicMaterial.h>
#include "gmpMechanicPoint.h"
 
class GmpCoupledDualPorosityMaterial : public GmpCoupledHMMaterial
{
protected:
  enum DualPorosityElementPropertyIds
  {
    Knf_ID,        
    Ksf_ID,        
    Ktf_ID,        
    Closure_ID,    
  };
 
  enum closureModel
  {
    normal,          
    hardStiffness,   
    bartonBandis,    
 
  };
 
public:
 
  GmpCoupledDualPorosityMaterial(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpCoupledHMMaterial(typeIndex, typeName, logger)
  {
  }
 
  virtual ~GmpCoupledDualPorosityMaterial() {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpCoupledDualPorosityMaterial(typeIndex, typeName, logger);
  }
 
  virtual const QVariantMap* materialMetaDataMap();
 
  virtual double normalElasticStiffness(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(Knf_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double shearElasticStiffness(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(Ksf_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double shearElasticStiffnessTwo(const GmElement* e, const GmVector* coord, int ip) const
  {
    if (propertyAc(Ktf_ID) == NULL)
    {
      return 0;
    }
    return propertyAc(Ktf_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual closureModel fractureClosureModel(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    GmpCoupledDualPorosityMaterial::closureModel type;
 
    if (propertyAc(Closure_ID) == NULL)
    {
      return GmpCoupledDualPorosityMaterial::closureModel::normal;
    }
 
    int closureType = propertyAc(Closure_ID)->scalarValueAt(e, coord, ip);
    //
    if (closureType == 1)
    {
      type = GmpCoupledDualPorosityMaterial::closureModel::hardStiffness;
    }
    //else if (closureType == 2)
    //{
    //  type = GmpCoupledDualPorosityMaterial::closureModel::bartonBandis;
    //}
    else
    {
      type = GmpCoupledDualPorosityMaterial::closureModel::normal;
    }
 
    return type;
  }
  
  virtual bool fillElasticStiffnessMatrix(const GmElement* e, GmMatrix& Dmfr, GmMatrix& Cm, GmMatrix& Cfr, const GmVector* coord, int ip, const GmpMechanicPoint* mp, unsigned sc) const;
  
  virtual void updateElasticStiffnessMatrix(const GmElement* e, GmMatrix& Dmfr, GmMatrix& Cm, GmMatrix& Cfr, const GmpMechanicPoint* mp, const GmVector* coord, int ip, unsigned sc) const;
  virtual void penalty(const GmElement* e, GmVector& b,  const GmVector* coord, int ip, const GmpMechanicPoint* mp, unsigned sc)const;
  
  virtual bool mechanicalConstitutiveModel(const GmElement* e, GmMatrix& Dmfr, GmMatrix& Cm, GmMatrix& Cfr, const GmpMechanicPoint* mp, const GmVector* coord, unsigned sc, unsigned ips) const;
  virtual void tangentModulus(const GmElement* e, GmMatrix& Dmfr, GmMatrix& Cm, GmMatrix& Cfr, const GmpMechanicPoint* mp, const GmVector* coord, unsigned nc) const;
 
  virtual bool fillFluidFlowFracture(const GmElement* e, GmMatrix& Kfc, GmMatrix B, double c, int d, const GmpMechanicPoint* mp, const GmVector* coord, unsigned nc) const;
  
  virtual bool fillFractureStorageCompMatrix(const GmElement* e, GmpHydraulicDualFlow* hp, GmMatrix& Sfc, const GmVector& N, double c, int d, const GmpMechanicPoint* mp, const GmVector* coord, unsigned sc)const;
 
  virtual bool fillTransferMatrix(const GmElement* e, GmpHydraulicDualFlow* hp, GmMatrix& Wfm, const GmVector& N, double c, int d, const GmpMechanicPoint* mp, const GmVector* coord, unsigned sc)const;
 
};
 
#endif