gmpMaterialCapMohrCoulomb.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2016 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_CAPMOHRCOULOMB_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_CAPMOHRCOULOMB_H_

#include "gmpMaterialElastoplastic.h"
#include "gmpMaterialModifMohrCoulomb.h"
#include "gmpMechanicPoint.h"
#include "gmpFemPhysics.h"
#include "gmMathUtils.h"

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialCapMohrCoulomb : public GmpMaterialModifiedMohrCoulomb
{
protected:
  enum ElementPropertyIds
  {
    PA_ID = GmpMaterialModifiedMohrCoulomb::NUM_PROPERTY_IDS, 
        R_ID,   

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

public:

  GmpMaterialCapMohrCoulomb(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpMaterialModifiedMohrCoulomb(typeIndex, typeName, logger) {}

  virtual ~GmpMaterialCapMohrCoulomb() {}

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

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

  // Returns the stresses according to the material behavior adopted
  virtual bool closestPointReturnAlgorithm(const GmElement* c, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned sc, bool ips) const;

          //yield criterion
  virtual double yieldCriterion(const GmElement*, const GmVector&, const GmVector*, int, unsigned) const;
  virtual bool yieldStressGradient(const GmElement*, GmVector&, const GmVector&, const GmVector*, int, unsigned, bool) const;
  virtual bool yieldHessian(const GmElement*, GmMatrix&, const GmVector&, const GmVector*, int, unsigned, bool) const;


  //plastic potential
  virtual double plasticFPotential(const GmElement*, const GmVector&, const GmVector*, int, unsigned) const;
  virtual bool flowVector(const GmElement*, GmVector&, const GmVector&, const GmVector*, int, unsigned, bool) const;
  virtual bool flowVectorStressGradient(const GmElement*, GmMatrix&, const GmVector&, const GmVector*, int, unsigned, bool) const;

  //Returns the stresses according to the material behavior adopted
  //virtual void returnMapping(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, unsigned, unsigned) const;
  // Line Search
  virtual double lineSearch(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, GmVector, double, unsigned, bool) const;
  // Residual Function
  virtual double residualFunction(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, GmVector S, double dLamb, unsigned, bool) const;

  virtual double yieldStrengthRatio(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const;

};

#endif