gmpHydraulicMaterialEquivalentContinuum.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_HYDRAULIC_MATERIAL_EQUIVALENTCONTINUUM_H_
#define _GEMA_PLUGIN_HYDRAULIC_MATERIAL_EQUIVALENTCONTINUUM_H_
 
#include "gmpMaterialSaturated.h"
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpHydraulicMaterialEquivalentContinuum : public GmpHydraulicMaterialSaturated
{
protected:
  enum EquivalentElementPropertyIds
  {
    Sfr_ID = GmpHydraulicMaterialSaturated::NUM_PROPERTY_IDS,    
    Bfr_ID,                                                      
    DIP_ID,                                                      
    STRIKE_ID,                                                   
    Uf_ID,                                                       
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
public:
 
  GmpHydraulicMaterialEquivalentContinuum(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpHydraulicMaterialSaturated(typeIndex, typeName, logger) {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpHydraulicMaterialEquivalentContinuum(typeIndex, typeName, logger);
  }
 
  // Fills equivalent continuum material map
  virtual const QVariantMap* materialMetaDataMap();
 
  virtual GmVector fractureSpacing(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    GmVector Sf;
    GmCellAccessor* sfAc = propertyAc(Sfr_ID);
    // fill property values
    if (!vectorPropertyValue(e, sfAc, Sf, coord, ip)) {
      Sf.set_size(1);
      return Sf.fill(0);
    }
    return Sf;
  }
 
  virtual GmVector fractureAperture(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    GmVector Bfr;
    GmCellAccessor* bfAc = propertyAc(Bfr_ID);
    // Fills property values
    vectorPropertyValue(e, bfAc, Bfr, coord, ip);
 
    return Bfr;
  }
 
  virtual GmVector fracturedDip(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    GmVector Dip;
    GmCellAccessor* dipAc = propertyAc(DIP_ID);
    // Fills property values
    vectorPropertyValue(e, dipAc, Dip, coord, ip);
 
    return Dip;
  }
 
  virtual GmVector fractureStrike(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    GmVector Strike;
    GmCellAccessor* strikeAc = propertyAc(STRIKE_ID);
    // fill property values
    if (!vectorPropertyValue(e, strikeAc, Strike, coord, ip)) {
      Strike.set_size(1);
      return Strike.fill(0);
    }
    return Strike;
  }
 
  virtual double fractureViscosity(const GmElement* e, const GmVector* coord, int ip) const
  {
          S_TRACE(); assert(e);
          return propertyAc(Uf_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual void hydraulicConstitutiveModel(const GmElement* e, const GmVector* coord, int ip, const GmpHydraulic::HMaterialPointAccessors& mp, GmMatrix& K) const;
  virtual bool fillFracturePermeability(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Kper) const;
  virtual GmMatrix fillFractureLocalPermeability(int dim, double bi, double si) const;
  virtual GmMatrix fillRotationMatrix(int dim, double dip, double strike) const;
};
 
#endif