xfemMaterial.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2017 by Geomechanic group - Tecgraf/Puc-Rio
** 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_XFEM_MATERIAL_H_
#define _GEMA_XFEM_MATERIAL_H_

#include <gmpHydraulicMaterial.h>
#include <gmpMechanicalMaterial.h>
#include <gmpInterfaceMaterial.h>

class XfemMaterial : public GmpFemPhysicsCommonMaterial
{
public:

  XfemMaterial(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpFemPhysicsCommonMaterial(typeIndex, typeName, logger)
  {
    _hydroMaterial = NULL;
    _mechanoMaterial = NULL;
        _mechDiscontinuityMaterial = NULL;
        _hydroDiscontinuityMaterial = NULL;

  }

  // destructor
  virtual ~XfemMaterial(){};

  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {    
    Q_UNUSED(simulation);
    return new XfemMaterial(typeIndex, typeName, logger);
  }

  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap();

  // Returns the constitutive tangent matrix
  virtual void tangentModulus(const GmElement* e, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned nc, unsigned ips) const;

  // Returns the stresses according to the material behavior adopted
  virtual void returnMapping(const GmElement* e, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned nc, unsigned ips) const;

  void setPhysicsMaterials(GmpFemPhysicsCommonMaterial* hydro, GmpFemPhysicsCommonMaterial* mechano)
  {
    _hydroMaterial = hydro;
    _mechanoMaterial = mechano;
  }

  void setPhysicsDiscontinuityMaterials(GmpFemPhysicsCommonMaterial* hydro, GmpFemPhysicsCommonMaterial* mechano)
  {
          _hydroDiscontinuityMaterial = hydro;
          _mechDiscontinuityMaterial  = mechano;
  }

  GmpHydraulicMaterial* hydro() const { return dynamic_cast<GmpHydraulicMaterial*>(_hydroMaterial); }

  GmpMechanicalMaterial* mechanical() const { return dynamic_cast<GmpMechanicalMaterial*>(_mechanoMaterial); }
  
  GmpInterfaceMaterial* mechanicalDiscontinuity() const { return dynamic_cast<GmpInterfaceMaterial*>(_mechDiscontinuityMaterial); }

  virtual double maximumPrincipalStress(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(SPMax_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double fractureOpening(const GmElement* e, const GmVector* coord, int ip) const
  {
          S_TRACE();
          assert(e);
          return propertyAc(GAP_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double normalElasticStiffness(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(Kni_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double shearElasticStiffness(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(Ksi_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double shearElasticStiffnessTwo(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(Kti_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double faultCohesion(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(CF_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double faultFrictionAngle(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(PHIF_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double faultDilationAngle(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(PSIF_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double cutOffTension(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(TCUT_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double fractureEnergy(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(Gf_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double shearModulus(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(G_ID)->scalarValueAt(e, coord, ip);
  }

  virtual double gravitationalAcceleration(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(GRAV_ID)->scalarValueAt(e, coord, ip);
  }

protected:
  enum ElemXfemPropertyIds
  {
    SPMax_ID,      
        GAP_ID,        
        Kni_ID,        
        Ksi_ID,        
        Kti_ID,        
        CF_ID,         
        PHIF_ID,       
        PSIF_ID,       
        TCUT_ID,       
        Gf_ID,         
        G_ID,         
        GRAV_ID,         

    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };

private:
  GmpFemPhysicsCommonMaterial*  _hydroMaterial;   
  GmpFemPhysicsCommonMaterial*  _mechanoMaterial; 
  GmpFemPhysicsCommonMaterial*  _mechDiscontinuityMaterial;  // TODO: Try to avoid this
  GmpFemPhysicsCommonMaterial*  _hydroDiscontinuityMaterial; // TODO: Try to avoid this
};
#endif