gmpMaterialCapModel.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_MODEL_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_CAP_MODEL_H_
 
#include "gmpMaterialElastoplasticMultiSurface.h"
#include "gmpMechanicPoint.h"
#include "gmpFemPhysics.h"
#include "gmMathUtils.h"
 
class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialCapModel : public GmpMaterialElastoplasticMultiSurface
{
protected:
        enum ElementPropertyIds
        {
                ALPHA_ID = GmpMaterialElastoplasticMultiSurface::NUM_PROPERTY_IDS, 
                BETA_ID,             
                LAMBDA_ID,           
                THETA_ID,            
                INITIALCENTER_ID,    
                RATIO_ID,            
                CUTOFF_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:
 
        GmpMaterialCapModel(int typeIndex, QString typeName, const GmLogCategory& logger)
                : GmpMaterialElastoplasticMultiSurface(typeIndex, typeName, logger) {}
 
        virtual ~GmpMaterialCapModel() {}
 
        static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
        {
                Q_UNUSED(simulation);
                return new GmpMaterialCapModel(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;
 
        // Returns the stresses according to the implicit return algorithm
        virtual bool implicitReturnAlgorithm(const GmElement* c, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned sc, bool ips) const;
 
        // Returns the stresses according to a Newton-Raphson implicit state-update algorithm based on complementary functions
        virtual bool cFunctionsReturnAlgorithm(const GmElement* c, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned sc, bool ips) const;
 
        // Returns the stresses according to a Newton-Krylov implicit state-update algorithm based on complementary functions
        virtual bool newtonKrylovReturnAlgorithm(const GmElement* c, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, 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 double yieldStrengthRatio(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const;
 
        virtual bool isIsotropic() const { return true; }
 
        virtual double alphaParameter(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(ALPHA_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double betaParameter(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(BETA_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double lambdaParameter(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(LAMBDA_ID)->scalarValueAt(e, coord, ip);
        }
        virtual double thetaParameter(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(THETA_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 tensionCutOff(const GmElement* e, const GmVector* coord, int ip) const
        {
                return propertyAc(CUTOFF_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;
 
        // Yield criterion - Tension Cut-off
        virtual double tensionCutOffCriterion(const GmElement*, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool tensionCutOffGradient(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool tensionCutOffGradient_Hardening(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool tensionCutOffHessian(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;
 
        // Plastic potential - Tension Cut-off
        virtual double tensionCutOffPlasticFPotential(const GmElement*, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool tensionCutOffPlasticFGradient(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool tensionCutOffPlasticFHessian(const GmElement*, GmMatrix&, const GmVector&, const double, const GmVector*, int, unsigned, bool) const;
        virtual bool tensionCutOffFlowVectorDerivative_Hardening(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned) 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;
 
        // Hardening Evolution - Tension Cut-off
        virtual bool tensionCutOffHardeningEvolution(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool tensionCutOffHardeningEvolutionDerivative_Stress(const GmElement*, GmVector&, const GmVector&, const double, const GmVector*, int, unsigned) const;
        virtual bool tensionCutOffHardeningEvolutionDerivative_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;
 
        // Method for the Newton-Krylov algorithm
        // Apply Givens Rotations
        virtual bool applyGivensRotation(const GmVector c, const GmVector s, GmVector& vrot, const int k) const;
        // Vector function to be solved by the Newton-Krylov method
        virtual bool vectorFunction(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, GmVector& r, 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 cphi, const double cdg, const double delta, unsigned sc, bool ips) const;
        // GMRES method
        virtual bool generalizedMinimumResidual(const GmVector fi, const GmMatrix M_inv, double tol, const int max_iter, const GmElement* e, const GmpMechanicPoint* mp,
                const GmVector* coord, const GmMatrix C, const GmVector de, GmVector& d, double& res, int& counter, const GmVector s_current, const GmVector s_old,
                const double q_current, const double q_old, const GmVector dg_current, const double cphi, const double cdg, const double delta, unsigned sc, bool ips) const;
        // Numerical Jacobian
        virtual bool numericalJacobian(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, GmMatrix& J, 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 cphi, const double cdg, const double delta, unsigned sc, bool ips) const;
 
};
 
#endif