gmpMaterialModifiedCamClay.h

Go to the documentation of this file.

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

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialModifiedCamClay : public GmpMaterialElastoplastic
{
protected:
  enum ElementPropertyIds
  { 
    M_ID = GmpMaterialElastoplastic::NUM_PROPERTY_IDS,          
    l_ID,               
    k_ID,               
    e0_ID,              
        beta_ID,        

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

  enum GaussAttributeIds
  {
          void_ID = GmpMaterialElastoplastic::NUM_GA_IDS,       

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

public:

  GmpMaterialModifiedCamClay(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpMaterialElastoplastic(typeIndex, typeName, logger) {}

  virtual ~GmpMaterialModifiedCamClay() {}

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


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

  virtual bool isIsotropic() const { return true; }

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

  double BetaParameter(const GmElement* e, const GmVector& s, const GmVector* coord, int ip, unsigned st) const
  {
          //Hidrostatic stress
          const double p_ = -GmpMechanicUtils::stressInvariantI1(s, st) / 3;
          //const double p_ = GmpMechanicUtils::hydrostaticStress(s, st);
          //von Mises Stress
          const double q_ = sqrt(GmpMechanicUtils::deviatoricInvariantJ2(s, st) * 3);
          //const double q_ = GmpMechanicUtils::vonMisesStress(s, st);

          //Critical line slope
          const double M_ = CriticalFrictionSlope(e, coord, ip);

          //Beta parameter
          double beta_ = 1.0;
          GmCellAccessor* B_ = propertyAc(beta_ID);
          if (B_ != NULL && p_ > -1 * M_*q_)
                  beta_ = propertyAc(beta_ID)->scalarValueAt(e, coord, ip);

          return beta_;
  }

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

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

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

  //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 double yielddP(const GmElement*, const GmVector&, const GmVector*, int, unsigned) const;

  //plastic potential
  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;
  virtual bool plasticdsdp(const GmElement*, GmVector&, const GmVector&, const GmVector*, int, unsigned) const;

  //internal variable
  virtual double CamClayHardeningLaw(const GmElement*, const GmVector&, const GmVector*, int, unsigned, bool) const;
  virtual bool CamClayHardeningGradient(const GmElement*, GmVector&, const GmVector&, const GmVector*, int, unsigned, bool) const;
  double CamClayHardeningGradientP(const GmElement*, const GmVector&, const GmVector*, int, unsigned, bool) const;

  //Returns the stresses according to the material behavior adopted
  virtual bool mechanicalConstitutiveModel(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, const GmVector&, unsigned, bool) const;
  
  virtual double yieldStrengthRatio(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const;

  //Returns the current void ratio scalar
  void voidratio(const GmElement*, const GmpMechanicPoint*, const GmVector*, unsigned) const;

  //Returns the constitutive elastic matrix
  void ElasticModulus(const GmElement*, const GmpMechanicPoint*, GmMatrix&, const GmVector*, unsigned) const;

  //Returns the second order unit tensor
  GmVector I2(unsigned) const;

  //Returns the symmetric fourth order unit tensor
  GmMatrix I4symm(unsigned) const;

};

#endif