gmpMaterialHencky.h

Go to the documentation of this file.

/************************************************************************
**
** 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_HENCKY_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_HENCKY_H_
 
#include "gmpMechanicalConfig.h"
#include "gmpFemPhysics.h"
#include "gmMathUtils.h"
 
#include "gmpMechanicalMaterial.h"
#include "gmpMechanicPoint.h"
#include "gmpMaterialHyperelastic.h"
 
class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMechanicalMaterialHencky : public GmpMechanicalMaterialHyperelastic
{
protected:
 
        // TO DO: how to implement an internal variable, free from concurrence?
        // principal streches and directions in the 
        // reference (Total Lagrangian) or current (Updated Lagrangian) configuration
        //GmVector* _pStrains;
        //GmMatrix* _pDirections;
 
public:
 
        GmpMechanicalMaterialHencky(int typeIndex, QString typeName, const GmLogCategory& logger)
                : GmpMechanicalMaterialHyperelastic(typeIndex, typeName, logger) {}
 
        virtual ~GmpMechanicalMaterialHencky() {}
 
        static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
        {
    Q_UNUSED(simulation);
                return new GmpMechanicalMaterialHencky(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* e, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector& Time, unsigned nc, bool ips) const;
 
        virtual bool isIsotropic() const { return true; }
 
        // Updates deformation gradient F 
        virtual bool updateDeformationGradient(GmMatrix& F, const GmElement* e, const GmVector* coord, int ip, bool ips = false) const;
 
        //
        virtual void potentialdLambda(const GmElement* e, const GmVector* coord, int ip, GmVector& dL, const GmVector& L) const;
 
        //
        virtual void potentialdLambdaSquare(const GmElement* e, const GmVector* coord, int ip, GmVector& dL2, const GmVector& L) const;
 
        //
        virtual void potentiald2LambdaSquare(const GmElement* e, const GmVector* coord, int ip, GmMatrix& d2L2, const GmVector& L) const;
 
 
        virtual void materialC(const GmVector& L, const GmMatrix& N, const GmVector& ss, const GmElement* e, unsigned sc, GmMatrix& D, const GmVector* coord, int ip) const;
        //virtual void dEGdPK2(const GmElement*, unsigned, GmMatrix&, const GmVector*, unsigned, int) const;
 
        //virtual void dsdea(const GmMatrix& C, const GmMatrix& F, GmMatrix& c, int d) const;
        virtual void spatialC(const GmVector& L, const GmMatrix& N, const GmVector& ss, const GmElement* e, unsigned sc, GmMatrix& D, const GmMatrix& F, const GmVector* coord, int ip) const;
 
        //GmVector PrincipalStrains;
        //GmMatrix PrincipalDirections;
};
 
#endif