gmpMaterialGeneralYieldSurface.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2016 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_GENERAL_YIELD_SURFACE_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_GENERAL_YIELD_SURFACE_H_
 
#include "gmpMaterialElastoplastic.h"
#include "gmpMechanicPoint.h"
#include "gmpFemPhysics.h"
#include "gmMathUtils.h"
 
class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialGeneralYieldSurface : public GmpMaterialElastoplastic
{
protected:
        enum ElementPropertyIds
        {
                COH_ID = GmpMaterialElastoplastic::NUM_PROPERTY_IDS, 
                PHI_ID, 
                PSI_ID, 
                ALPHA_ID, 
                BETA_ID, 
                GAMMA_ID, 
                A_ID, 
 
                // --- NO ADDING BELOW THIS LINE
                NUM_PROPERTY_IDS,  
        };
 
public:
 
        GmpMaterialGeneralYieldSurface(int typeIndex, QString typeName, const GmLogCategory& logger)
                : GmpMaterialElastoplastic(typeIndex, typeName, logger) {}
 
        virtual ~GmpMaterialGeneralYieldSurface() {}
 
        static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
        {
    Q_UNUSED(simulation);
                return new GmpMaterialGeneralYieldSurface(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(COH_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double frictionAngle(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(PHI_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double dilationAngle(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(PSI_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double firstShapeFactor(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(ALPHA_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double secondShapeFactor(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(BETA_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double thirdShapeFactor(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(GAMMA_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double fourthShapeFactor(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(A_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 bool yieldHessian(const GmElement*, GmMatrix&, const GmVector&, const GmVector*, int, unsigned, bool) const;
 
        // Plastic potential
        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 yieldStrengthRatio(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const;
 
};
 
#endif