gmpMaterialElastoplasticMultiSurface.h

/************************************************************************
**
** Copyright (C) 2014 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_ELASTOPLASTIC_MULTIPLESURFACE_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_ELASTOPLASTIC_MULTIPLESURFACE_H_

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

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

    //destructor
    virtual ~GmpMaterialElastoplasticMultiSurface() {}

    // 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 closet point return algorithm
    virtual bool implicitReturnAlgorithm(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;

    // Yield criterion
    virtual double yieldCriterion(const GmElement*, const GmVector&, const GmVector*, int, unsigned) const = 0;
    // Stress gradient of Yield function
    virtual bool yieldStressGradient(const GmElement*, GmVector&, const GmVector&, const GmVector*, int, unsigned, bool) const = 0;

    // Flow vector (Stress gradient of potencial function)
    virtual bool flowVector(const GmElement*, GmVector&, const GmVector&, const GmVector*, int, unsigned, bool) const = 0;
    // Stress gradient of Flow vector
    virtual bool flowVectorStressGradient(const GmElement*, GmMatrix&, const GmVector&, const GmVector*, int, unsigned, bool) 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 hardemningHessian(const GmElement*, const GmpMechanicPoint*, const GmVector*, int, unsigned) const { return  0.0; };

};

#endif