gmpMaterialCapGeneralYieldSurface.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_CAP_GENERAL_YIELD_SURFACE_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_CAP_GENERAL_YIELD_SURFACE_H_

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialCapGeneralYieldSurface : public GmpMaterialElastoplasticMultiSurface
{
protected:
        enum ElementPropertyIds
        {
                COH_ID = GmpMaterialElastoplasticMultiSurface::NUM_PROPERTY_IDS, 
                PHI_ID, 
                PSI_ID, 
                ALPHA_ID, 
                BETA_ID, 
                GAMMA_ID, 
                A_ID, 
                INITIALCENTER_ID,    
                RATIO_ID,            
                MAXVSTRAIN_ID,       
                VSTRAINRATE_ID,      

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

        enum GaussAttributeIds
        {
                MDLAMBD_GA_ID = GmpMaterialElastoplastic::NUM_GA_IDS,      

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

public:

        GmpMaterialCapGeneralYieldSurface(int typeIndex, QString typeName, const GmLogCategory& logger)
                : GmpMaterialElastoplasticMultiSurface(typeIndex, typeName, logger) {}

        virtual ~GmpMaterialCapGeneralYieldSurface() {}

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

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

        // Returns the stresses according to the material behavior adopted
        virtual bool mechanicalConstitutiveModel(const GmElement* c, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector& Time, unsigned sc, bool ips) const;
  
  // Sets the initial conditions required by Cap model material
  virtual bool setInitialConditions(const GmElement* e, GmpMechanicPoint* mp, const GmVector* coord, unsigned sc) const;

        virtual bool isIsotropic() const { return true; }

  virtual double yieldStrengthRatio(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const;

        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);
        }
        virtual double initialCenterEllipse(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(INITIALCENTER_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double ratioSemiaxes(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(RATIO_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double maxVolumetricStrain(const GmElement* e, const GmVector* coord, int ip) const
        {
                S_TRACE(); assert(e);
                // Get accessor for the ultimate compressive volumetric strain
                GmCellAccessor* wAccc = propertyAc(MAXVSTRAIN_ID);

                if (wAccc == NULL)
                {
                        return 1.0e10;
                }

                return wAccc->scalarValueAt(e, coord, ip);
        }

        virtual double volumetricStrainRate(const GmElement* e, const GmVector* coord, int ip) const
        {
                S_TRACE(); assert(e);
                // Get accessor for the volumetric strain rate with respect to the compressive hydrostatic strain
                GmCellAccessor* dAccc = propertyAc(VSTRAINRATE_ID);

                if (dAccc == NULL)
                {
                        return 1.0e-3;
                }

                return dAccc->scalarValueAt(e, coord, ip);
        }

        // Yield criterion - Conical Surface
        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 yieldGradient_Hardening(const GmElement*, GmVector&, const GmVector&, const GmVector*, int, unsigned) const;
        virtual bool yieldHessian(const GmElement*, GmMatrix&, const GmVector&, const GmVector*, int, unsigned, bool) const;

        // Yield criterion - Cap Surface
        virtual double capSurfaceCriterion(const GmElement*, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool capSurfaceGradient(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool capSurfaceGradient_Hardening(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool capSurfaceHessian(const GmElement*, GmMatrix&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;

        // Plastic potential - Conical Surface
        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 bool flowVectorDerivative_Hardening(const GmElement*, GmVector&, const GmVector&, const GmVector*, int, unsigned) const;

        // Plastic potential - Cap Surface
        virtual double capSurfacePlasticFPotential(const GmElement*, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool capSurfacePlasticFGradient(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool capSurfacePlasticFHessian(const GmElement*, GmMatrix&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool capSurfaceFlowVectorDerivative_Hardening(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;

        // Hardening Evolution - Conical Surface
        virtual bool hardeningEvolution(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool hardeningEvolutionDerivative_Stress(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool hardeningEvolutionDerivative_Hardening(const GmElement*, GmMatrix&, const GmVector&, const double, const GmVector*, int, unsigned) const;

        // Hardening Evolution - Cap Surface
        virtual bool capSurfaceHardeningEvolution(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool capSurfaceHardeningEvolutionDerivative_Stress(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool capSurfaceHardeningEvolutionDerivative_Hardening(const GmElement*, GmMatrix&, const GmVector&, const double, const GmVector*, int, unsigned) const;

        // Complementary Methods
        virtual bool yieldSurfaces(const GmElement*, GmVector&, GmVector, const double, const GmVector*, int, unsigned) const;
        virtual bool flowVectors(const GmElement*, GmVector&, int, GmVector, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool hardeningLaws(const GmElement*, GmVector&, int, GmVector, const double, const GmVector*, int, unsigned) const;
        virtual bool gradientVectors_Stress(const GmElement*, GmVector&, int, GmVector, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool gradientVectors_Hardening(const GmElement*, GmVector&, int, GmVector, const double, const GmVector*, int, unsigned) const;
        virtual bool flowVectorDerivatives_Stress(const GmElement*, GmMatrix&, int, GmVector, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool flowVectorDerivatives_Hardening(const GmElement*, GmVector&, int, GmVector, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool hardeningLawDerivatives_Stress(const GmElement*, GmVector&, int, GmVector, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool hardeningLawDerivatives_Hardening(const GmElement*, GmMatrix&, int, GmVector, const double q, const GmVector*, int, unsigned) const;
        void updateDLamb(GmVector&, GmVector, GmVector) const;
        virtual int numberActiveSurface(GmVector, const double) const;
        void indexActiveSurfaces(GmVector&, GmVector, int, const double) const;
        void deactivateSurfaces(GmVector&, GmVector, int) const;

        // Line Search
        virtual double lineSearch(const GmElement*, GmMatrix, GmVector, GmVector, GmVector, GmVector, const double,
                const double, const double, GmVector, GmVector, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;
  virtual double goldenSectionMethod(const GmElement*, GmMatrix, GmVector, GmVector, GmVector, GmVector, const double,
        const double, const double, GmVector, GmVector, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;
  virtual double quadraticInterpolation(const GmElement*, GmMatrix, GmVector, GmVector, GmVector, GmVector, const double,
        const double, const double, GmVector, GmVector, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;
  virtual double cubicInterpolation(const GmElement*, GmMatrix, GmVector, GmVector, GmVector, GmVector, const double,
        const double, const double, GmVector, GmVector, const GmpMechanicPoint*, const GmVector*, unsigned, bool) const;

        // Residual Function
        virtual double residualFunction(const GmElement*, GmMatrix, const GmpMechanicPoint*, const GmVector*, GmVector,
                GmVector, GmVector, const double, const double, GmVector, unsigned, bool) const;

};

#endif