gmpMaterialVonMises.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_VONMISES_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_VONMISES_H_

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialVonMises : public GmpMaterialElastoplastic
{
public:

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

    virtual ~GmpMaterialVonMises() {}

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

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

    virtual bool isIsotropic() const { return true; }
    virtual double initialYieldStress(const GmElement* e, const GmVector* coord, int ip) const
    {
      return propertyAc(IYST_ID)->scalarValueAt(e, coord, ip);
    }
    virtual double plasticModulus(const GmElement* e, const GmVector* coord, int ip) const
    {
      GmCellAccessor* pmAc = propertyAc(PMD_ID);
      if (pmAc == NULL)
      {
        return 0.0;
      }
      return propertyAc(PMD_ID)->scalarValueAt(e, coord, ip);
    }
    // Elastoplastic Model Equations
    virtual double vonMisesCriterion(const GmElement*, const GmVector, const double, const GmVector*, int, unsigned, bool) const;
    virtual bool vonMisesStressGradient(const GmElement*, GmVector&, const GmVector, const double, const GmVector*, int, unsigned, bool) const;
    virtual bool vonMisesHardeningGradient(const GmElement*, GmVector&, const GmVector, const double, const GmVector*, int, unsigned, bool) const;
    virtual bool vonMisesflowVectorStressGradient(const GmElement*, GmMatrix&, const GmVector, const double, const GmVector*, int, unsigned, bool) const;
    virtual bool vonMisesflowVectorHardeningGradient(const GmElement*, GmVector&, const GmVector, const double, const GmVector*, int, unsigned, bool) const;
    virtual double hardeningLaw(const GmElement*, const GmVector, const double, const GmVector*, int, unsigned, bool) const;
    virtual bool hardeningLawStressGradient(const GmElement*, GmVector&, const GmVector, const double, const GmVector*, int, unsigned, bool) const;
    virtual bool hardeningLawHardeningGradient(const GmElement*, GmVector&, const GmVector, const double, const GmVector*, int, unsigned, bool) const;
    virtual double yieldStrengthRatio(const GmElement*, const GmVector, const double, const GmVector*, int, unsigned, bool) const;

    //
    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;
    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;
    virtual double plasticPotential(const GmElement*, const GmpMechanicPoint*, const GmVector*, int, unsigned) const;
    virtual double plasticGradient(const GmElement*, const GmpMechanicPoint*, const GmVector*, int, unsigned) const;
    virtual double plasticHessian(const GmElement*, const GmpMechanicPoint*, const GmVector*, int, unsigned) const;
    virtual double yieldStrengthRatio(const GmElement*, const GmVector&, const GmVector*, int, unsigned) const { return  0.0; };
    

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

protected:
  enum ElementPropertyIds
  {
    IYST_ID = GmpMaterialElastoplastic::NUM_PROPERTY_IDS, 
    PMD_ID,           

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

 };

#endif