gmpMaterialElastoplastic.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_ELASTOPLASTIC_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_ELASTOPLASTIC_H_

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialElastoplastic : public GmpMechanicalMaterialElastic
{
public:
  //constructors
  GmpMaterialElastoplastic(int typeIndex, QString typeName, const GmLogCategory& logger)
  : GmpMechanicalMaterialElastic(typeIndex, typeName, logger) {}

  //destructor
  virtual ~GmpMaterialElastoplastic() {}

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

  enum returnMappingStrategy
  {
          Explicit,            
          SemiImplicit,        
          Implicit,            
          NewtonKrylov,        
          CFunctions,          
          Alternative,         
  };

  enum lineSearchStrategy
  {
      none,            
      goldenSection,   
      quadratic,       
      cubic,           
  };

  enum ElementPropertyIds
  {
    RMA_ID = GmpMechanicalMaterialElastic::NUM_PROPERTY_IDS, 
    SUBSTEP_ID,    
    LINESEARCH_ID, 
    FTOL_ID,       
    STOL_ID,       
    HTOL_ID,       

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

  enum elastoplasticGaussAttrIds
  {
    EP_GA_ID,          
    EPOLD_GA_ID,       
    HARD_GA_ID,        
    HARDOLD_GA_ID,     
    DLAMBD_GA_ID,      

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

  //Returns the stresses according to the material behavior adopted
  virtual bool mechanicalConstitutiveModel(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, const GmVector&, unsigned, bool) const;

  // Returns the stresses according to the explicit return algorithm
  virtual bool explicitReturnAlgorithm(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const ;

  // Returns the stresses according to the cutting plane return algorithm
  virtual bool semiImplicitReturnAlgorithm(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;

  // Returns the stresses according to the implicit return algorithm
  virtual bool implicitReturnAlgorithm(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const ;

  // Returns the stresses according to a Newton-Krylov return mapping algorithm
  virtual bool newtonKrylovReturnAlgorithm(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;

  // Returns the stresses according to a Newton-Raphson implicit state-update algorithm based on complementary functions
  virtual bool cFunctionsReturnAlgorithm(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;

  // Returns the stresses according to an alternative return mapping algorithm
  virtual bool alternativeReturnAlgorithm(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;

  // Yield criterion
  virtual double yieldCriterion(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const = 0;
  // Stress gradient of Yield function
  virtual bool yieldStressGradient(const GmElement* e, GmVector& dfs, const GmVector& S, const GmVector* coord, int ip, unsigned sc, bool ips) const = 0;

  // Flow vector (Stress gradient of potencial function)
  virtual bool flowVector(const GmElement* e, GmVector& dgs, const GmVector& s, const GmVector* coord, int ip, unsigned sc, bool ips) const = 0;
  // Stress gradient of Flow vector
  virtual bool flowVectorStressGradient(const GmElement* e, GmMatrix& dgss, const GmVector& s, const GmVector* coord, int ip, unsigned sc, bool ips) const = 0;

  // Hardening law
  virtual double hardeningLaw(const GmElement* , const GmpMechanicPoint*, const GmVector*, int, unsigned) const { return  0.0; }
  // Stress gradient of Hardening law
  virtual double hardeningStressGradient(const GmElement* , const GmpMechanicPoint*, const GmVector*, int, unsigned) const { return  0.0; };
  // hardening gradient of hardening law
  virtual double hardeningHessian(const GmElement* , const GmpMechanicPoint*, const GmVector*, int, unsigned) const { return  0.0; };
  
  virtual double yieldStrengthRatio(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const = 0;

  // Line Search for the optimized return algorithm
  virtual double lineSearch(const GmElement*, GmMatrix&, GmVector&, GmVector&, const GmpMechanicPoint*, const GmVector*, 
              double, double, unsigned, bool) const;
  // Golden section method
  virtual double goldenSectionMethod(const GmElement*, GmMatrix&, GmVector&, GmVector&, const GmpMechanicPoint*, const GmVector*, 
              double, double, unsigned, bool) const;
  // Quadratic interpolation
  virtual double quadraticInterpolation(const GmElement*, GmMatrix&, GmVector&, GmVector&, const GmpMechanicPoint*, const GmVector*,
      double, double, unsigned, bool) const;
  // Quadratic interpolation
  virtual double cubicInterpolation(const GmElement*, GmMatrix&, GmVector&, GmVector&, const GmpMechanicPoint*, const GmVector*,
      double, double, unsigned, bool) const;

  // Residual Function for the optimized return algorithm
  virtual double residualFunction(const GmElement*, GmMatrix&, GmVector&, const GmpMechanicPoint*, const GmVector*, double, unsigned, bool) const;

  virtual double returnAlgorithm(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    GmCellAccessor* rmAcc = propertyAc(RMA_ID);

    // Read return mapping option 
    double value = 2;
    if (rmAcc != NULL)
    {
      value = rmAcc->scalarValueAt(e, coord, ip);
    }

    return value;
  }

  virtual bool substepping(const GmElement* e, const GmVector* coord, int ip) const
  {
      S_TRACE();
      assert(e);
      GmCellAccessor* subStepAcc = propertyAc(SUBSTEP_ID);

      // Read substepping option 
      bool value = false;
      if (subStepAcc != NULL)
      {
          value = (subStepAcc->scalarValueAt(e, coord, ip) > 0)? true : false;
      }

      return value;
  }

  virtual double lineSearchStrategy(const GmElement* e, const GmVector* coord, int ip) const
  {
      S_TRACE();
      assert(e);
      GmCellAccessor* lineSearchAcc = propertyAc(LINESEARCH_ID);

      // Read line search option 
      double value = 1;
      if (lineSearchAcc != NULL)
      {
          value = lineSearchAcc->scalarValueAt(e, coord, ip);
      }

      return value;
  }

  virtual double yieldTolerance(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    // Reads accessor of Yield function tolerance
    GmCellAccessor* fTolAcc = propertyAc(FTOL_ID);
    if (fTolAcc == NULL)
    {
      return 1.0e-6;
    }
    return fTolAcc->scalarValueAt(e, coord, ip);
  }

  virtual double stressTolerance(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    GmCellAccessor* sTolAcc = propertyAc(STOL_ID);
    // Reads accessor of Stress tolerance
    if (sTolAcc == NULL)
    {
      return 1.0e-6;
    }
    return sTolAcc->scalarValueAt(e, coord, ip);
  }

  virtual double hardeningTolerance(const GmElement* e, const GmVector* coord, int ip) const
  {
        S_TRACE();
        assert(e);
    // Reads accessor of hardening tolerance
        GmCellAccessor* hTolAcc = propertyAc(HTOL_ID);
        if (hTolAcc == NULL)
        {
         return 1.0e-6;
        }
        return hTolAcc->scalarValueAt(e, coord, ip);
  }
};

#endif