gmpMaterialInterfaceMohrCoulomb.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_MECHANICALMATERIAL_INTERFACEMOHRCOULOMB_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_INTERFACEMOHRCOULOMB_H_


#include "gmpInterfaceMaterial.h"
#include "gmpMaterialElasticInterface.h"
#include "gmpMechanicPoint.h"
#include <gmpFemPhysics.h>
#include <gmMathUtils.h>

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpInterfaceMohrCoulomb: public GmpElasticInterface
{
protected:
  enum InterfaceMCPropertyIds
  {                 
    CF_ID = GmpElasticInterface::NUM_PROPERTY_IDS,        
    PHIF_ID,        
    PSIF_ID,        
    TCUT_ID,        

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


public:
  GmpInterfaceMohrCoulomb(int typeIndex, QString typeName, const GmLogCategory& logger)
    :GmpElasticInterface(typeIndex, typeName, logger){}

  virtual ~GmpInterfaceMohrCoulomb() {}

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

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

  enum GaussAttributeIds
  {
    EP_GA_ID = GmpElasticInterface::NUM_GA_IDS, 
    EPOLD_GA_ID,       
    //STR_GA_ID,         //!< Base Id for retrieving the slip tendency ratio accessor at the current state
    //STROLD_GA_ID,      //!< Id for retrieving the internal variable accessor at the previous state


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

  virtual bool mechanicalConstitutiveModel(const GmElement* e, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned nc) const;

  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);
  }
};
#endif