gmpMaterialCrackingCrushingPlasticDamage.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_CRACKINGCRUSHINGPLASTICDAMAGE_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_CRACKINGCRUSHINGPLASTICDAMAGE_H_

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialCrackingCrushingPlasticDamage : public GmpMaterialElastoplasticMultiSurface
{
protected:
        enum ElementPropertyIds
        {
                FCM_ID = GmpMaterialElastoplasticMultiSurface::NUM_PROPERTY_IDS, 
                RCM_ID,              
                FT0_ID,              
                FBT_ID,              
                FTC_ID,              
                FC0_ID,              
                FBC_ID,              
                RHO_ID,              
                GF_ID,               
                GC_ID,               
                LCH_ID,              
                LCHMIN_ID,           

                COH_ID,              
                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,      
                DAM_GA_ID,           
                DAMOLD_GA_ID,        
                TISV_GA_ID,          
                TISVOLD_GA_ID,       
                CISV_GA_ID,          
                CISVOLD_GA_ID,       

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

public:

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

        virtual ~GmpMaterialCrackingCrushingPlasticDamage() {}

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

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

        // Returns the stresses according to the return mapping algorithm
        virtual bool returnMappingAlgorithm(const GmElement* c, GmVector& s, GmVector de, GmVector& ep, double& q, GmVector& DLamb, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned sc, bool ips, bool stress_state) const;

        // 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;

        // Checks material loaded data. 
        // virtual bool checkLoadedData(const GmElement* e) 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 ultimateUniaxialCompressiveStrength(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(FCM_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double strainAtUltimateCompressiveStrength(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(RCM_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double uniaxialTensileStrength(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(FT0_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double biaxialTensileStrength(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(FBT_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double tensileCriterionShapeParameter(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(FTC_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double elasticUniaxialCompressiveLimit(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(FC0_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double biaxialCompressiveStrength(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(FBC_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double parameterRho(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(RHO_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double fractureEnergy(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(GF_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double crushingEnergy(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(GC_ID)->scalarValueAt(e, coord, ip);
        }

        virtual double userDefinedCharacteristicLength(const GmElement* e, const GmVector* coord, int ip) const
        {
                S_TRACE(); assert(e);
                // Get accessor for the user-defined characteristic length
                GmCellAccessor* lch = propertyAc(LCH_ID);

                if (lch == NULL)
                {
                        return 0.0;
                }

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

        virtual double minimumCharacteristicLength(const GmElement* e, const GmVector* coord, int ip) const
        {
                S_TRACE(); assert(e);
                // Get accessor for the minimum characteristic length
                GmCellAccessor* lchmin = propertyAc(LCHMIN_ID);

                if (lchmin == NULL)
                {
                        return 0.0;
                }

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

        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);
        }

        // Damage criterion definition
        virtual bool damageCriterion(const GmElement*, const GmpMechanicPoint*, const GmVector*, unsigned, GmVector, double&, double&, GmVector&, GmVector&, unsigned) const;

        // Damage evolution law definition
        virtual bool damageLaw(const GmElement*, const GmVector*, const GmpMechanicPoint*, double, double, double&, double&, double& dg) const;

        // Returns the characteristic length
        virtual double characteristicLength(const GmElement*, const GmVector*, int, GmValueAccessor*, bool type = false) const;

        // 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;

        // Line Search
        virtual double lineSearch(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, GmMatrix C, GmVector de, const double dde_zz,
                GmVector s, GmVector s_old, GmVector ds, const double q, const double q_old, const double dq, GmVector dg, GmVector ddg, const double cdg, const double delta,
                unsigned sc, bool ips) const;

        // Golden Section Method
        virtual double goldenSectionMethod(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, GmMatrix C, GmVector de, const double dde_zz,
                GmVector s, GmVector s_old, GmVector ds, const double q, const double q_old, const double dq, GmVector dg, GmVector ddg, const double cdg, const double delta,
                unsigned sc, bool ips) const;

        // Quadratic Interpolation
        virtual double quadraticInterpolation(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, GmMatrix C, GmVector de, const double dde_zz,
                GmVector s, GmVector s_old, GmVector ds, const double q, const double q_old, const double dq, GmVector dg, GmVector ddg, const double cdg, const double delta,
                unsigned sc, bool ips) const;

        // Cubic Interpolation
        virtual double cubicInterpolation(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, GmMatrix C, GmVector de, const double dde_zz,
                GmVector s, GmVector s_old, GmVector ds, const double q, const double q_old, const double dq, GmVector dg, GmVector ddg, const double cdg, const double delta,
                unsigned sc, bool ips) const;

        // Objective Function to be minimized by the Newton-Raphson method
        virtual double objectiveFunction(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, const GmMatrix C, const GmVector de,
                const GmVector s, const GmVector s_old, const double q, const double q_old, const GmVector dg, const double cdg, const double delta, unsigned sc, bool ips) const;

};

#endif