gmpMaterialCreep.h

/************************************************************************
**
** 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_CREEP_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_CREEP_H_

#include "gmpMechanicalConfig.h"
#include <gmpFemPhysics.h>
#include <gmMathUtils.h>

#include <gmpMaterialElastic.h>
#include "gmpMechanicPoint.h"

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialCreep : public GmpMechanicalMaterialElastic
{
protected:
  enum ElementPropertyIds
  {
    T_ID = GmpMechanicalMaterialElastic::NUM_PROPERTY_IDS, 
    STOL_ID,           
    DILATCRIT_ID,      
    DILATPARAM_ID,     

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

  enum GaussAttributeIds
  {
    EC_GA_ID,       
    ECOLD_GA_ID,    

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

  enum NodeAttributeIds
  {
          T_NA_ID,           

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

  enum dilatancyCriterion
  {
    vanSambeek,      
    DeVries,         
    HunscheMod,      

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

public:

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

  virtual ~GmpMaterialCreep() {}

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

  // Checks material loaded data. 
  virtual bool checkLoadedData(const GmElement* e) const;

  // 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 double yieldStrengthRatio(const GmElement* e, const GmVector& S, const GmVector* coord, int ip, unsigned sc) const;
  
  virtual double fillTemperatureFromNodalAttr(const GmElement* e, const GmVector* coord) const;

  virtual double fillCreepStrainRate(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector& Time, double Smises) const = 0;

  virtual void evaluatesCreepStrain(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector& Stress, const GmVector& Time, GmVector& dStrainCreep, unsigned sc) const;

  virtual bool isIsotropic() const { return true; }

  virtual double elasticModulus(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(E_ID)->scalarValueAt(e, coord, ip);
  }
  virtual double poissonRatio(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(NU_ID)->scalarValueAt(e, coord, ip);
  }
  virtual double temperatureProperty(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(T_ID)->scalarValueAt(e, coord, ip);
  }
  virtual double stressTolerance(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    GmCellAccessor* sTolAcc = propertyAc(STOL_ID);
    // Reads accessor of Stress tolerance
    if (sTolAcc == NULL)
    {
      return 1.0e-4;
    }
    return sTolAcc->scalarValueAt(e, coord, ip);
  }
  virtual int dilatancyCriterion(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    GmCellAccessor* DilCrtAcc = propertyAc(DILATCRIT_ID);
    if (DilCrtAcc == NULL)
    {
      return 0;
    }
    return DilCrtAcc->scalarValueAt(e, coord, ip);
  }
  virtual GmVector dilatancyParameters(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector data;
    GmCellAccessor* DPAcc = propertyAc(DILATPARAM_ID);
    if (DPAcc != NULL) {
      vectorPropertyValue(e, DPAcc, data, coord, ip);
    }
    return data;
  }
};

#endif