gmpMaterialBreslerPister.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Nilthson Noreña
** 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_BRESLER_PISTER_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_BRESLER_PISTER_H_
 
#include "gmpMaterialElastoplastic.h"
#include "gmpMechanicPoint.h"
#include "gmpFemPhysics.h"
#include "gmMathUtils.h"
 
class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialBreslerPister : public GmpMaterialElastoplastic
{
protected:
  enum ElementPropertyIds
  { 
    ST_ID = GmpMaterialElastoplastic::NUM_PROPERTY_IDS, 
    SC_ID,         
    SB_ID,         
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
public:
 
  GmpMaterialBreslerPister(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpMaterialElastoplastic(typeIndex, typeName, logger) {}
 
  virtual ~GmpMaterialBreslerPister() {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialBreslerPister(typeIndex, typeName, logger);
  }
 
  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap();
 
  virtual bool isIsotropic() const { return true; }
 
  virtual double cohesion(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(ST_ID)->scalarValueAt(e, coord, ip);
  }
  virtual double frictionAngle(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(SC_ID)->scalarValueAt(e, coord, ip);
  }
  virtual double dilationAngle(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(SB_ID)->scalarValueAt(e, coord, ip);
  }
 
  // Yield criterion
  virtual double yieldCriterion(const GmElement* e, const GmVector& s, const GmVector* coord, int ip, unsigned sc) const;
  virtual bool yieldStressGradient(const GmElement* e, GmVector& dfs, const GmVector& s, const GmVector* coord, int ip, unsigned sc, bool ips) const;
  virtual bool yieldHessian(const GmElement* e, GmMatrix& dfss, const GmVector& s, const GmVector* coord, int ip, unsigned sc, bool ips) const;
 
  // Plastic potential
  virtual double plasticFPotential(const GmElement* e, const GmVector& s, const GmVector* coord, int ip, unsigned sc) const;
  virtual bool flowVector(const GmElement* e, GmVector& dgs, const GmVector& s, const GmVector* coord, int ip, unsigned sc, bool ips) const;
  virtual bool flowVectorStressGradient(const GmElement* e, GmMatrix& dgss, const GmVector& s, const GmVector* coord, int ip, unsigned sc, bool ips) const;
 
  virtual double yieldStrengthRatio(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const;
 
};
 
#endif