gmpMaterialOneDimPlasticity.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Nilthson Noreña
** 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_1DPLASTICITY_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_1DPLASTICITY_H_

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialOneDimPlasticity : public GmpMechanicalMaterialElastic
{
protected:
    enum ElementPropertyIds
    {
        IYST_ID = GmpMechanicalMaterialElastic::NUM_PROPERTY_IDS, 
        PMD_ID,           

        // --- NO ADDING BELOW THIS LINE
        NUM_PROPERTY_IDS,  
    };
    enum elemGaussAttrIds
    {
      EP_GA_ID,          
      EPOLD_GA_ID,       
      HARD_GA_ID,        
      HARDOLD_GA_ID,     
      DLAMBD_GA_ID,      

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

    GmpMaterialOneDimPlasticity(int typeIndex, QString typeName, const GmLogCategory& logger)
        : GmpMechanicalMaterialElastic(typeIndex, typeName, logger) {}

    virtual ~GmpMaterialOneDimPlasticity() {}

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

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

    virtual bool isIsotropic() const { return true; }

    virtual double initialYieldStress(const GmElement* e, const GmVector* coord, int ip) const
    {
        return propertyAc(IYST_ID)->scalarValueAt(e, coord, ip);
    }
    virtual double plasticModulus(const GmElement* e, const GmVector* coord, int ip) const
    {
        GmCellAccessor* pmAc = propertyAc(PMD_ID);
        if (pmAc == NULL)
        {
            return 0.0;
        }
        return pmAc->scalarValueAt(e, coord, ip);
    }

    // Returns the stresses according to the material behavior adopted
    virtual bool mechanicalConstitutiveModel(const GmElement*, GmMatrix&, const GmpMechanicPoint*, const GmVector*, const GmVector&, unsigned, bool) const;

 };

#endif