gmpMaterialDMCreep.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_DMCREEP_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_DMCREEP_H_

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

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialDMCreep : public GmpMaterialCreep
{
protected:
  enum ElementPropertyIds
  {
    QT_ID = GmpMaterialCreep::ElementPropertyIds::NUM_PROPERTY_IDS,  
    R_ID,         
    DEDT_ID,      
    So_ID,        
    To_ID,        
    AN1_ID,       
    AN2_ID,       

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

public:

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

  virtual ~GmpMaterialDMCreep() {}

  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialDMCreep(typeIndex, typeName, logger);
  }
 
  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap()
  {
    S_TRACE();
    static QVariantMap m;
    if (m.isEmpty())
    {
      // read the material attribute map from elastic interface material
      m = *GmpMaterialCreep::materialMetaDataMap();

      // Adds material properties
      GmpFemPhysicsCommon::ValueList matDMProp = m["properties"].value<GmpFemPhysicsCommon::ValueList>();
      matDMProp.append(GmpFemPhysicsCommon::ScalarValue(QT_ID, "QT", QObject::tr("Thermal activation energy"), "J/mol", true));
      matDMProp.append(GmpFemPhysicsCommon::ScalarValue(R_ID,  "R", QObject::tr("Universal gas constant"), "J/(mol*K)", true));
      matDMProp.append(GmpFemPhysicsCommon::ScalarValue(DEDT_ID, "DEDT", QObject::tr("Threshold creep rate"), "1/s", true));
      matDMProp.append(GmpFemPhysicsCommon::ScalarValue(So_ID,  "So", QObject::tr("Threshold deviatoric stress"), "kPa", true));
      matDMProp.append(GmpFemPhysicsCommon::ScalarValue(To_ID,  "To", QObject::tr("Threshold temperature"), "K", true));
      matDMProp.append(GmpFemPhysicsCommon::ScalarValue(AN1_ID, "AN1", QObject::tr("Dislocation creep"), "", true));
      matDMProp.append(GmpFemPhysicsCommon::ScalarValue(AN2_ID, "AN2", QObject::tr("Steady state cracking"), "", true));
      m["properties"] = matDMProp;
    }
    return &m;
  }

  // Returns the stresses according to the material behavior adopted
  virtual double fillCreepStrainRate(const GmElement* e, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector& Time, double MISES) const
  {
    S_TRACE();
    assert(e); Q_UNUSED(Time);

    double eCreepRate, NP, So, Eo, QT, Rg, To, TEMP;

    // Reads DM creep Properties 
    So = thresholdDevStress(e, coord, mp->_index);
    Eo = thresholdCreepRate(e, coord, mp->_index);
    QT = thermalActEnergy(e, coord, mp->_index);
    Rg = universalGasConstant(e, coord, mp->_index);
    To = thresholdTemperature(e, coord, mp->_index);

    // Read the Temperature at integration point
        if (nodeAttrAc(T_NA_ID)!= NULL)
        {
                TEMP = fillTemperatureFromNodalAttr(e, coord);
        }
        else
        {
                TEMP = temperatureProperty(e, coord, mp->_index);
        }

    // Set the suitable stress power
    if (MISES < So)
    {
      NP = disCreep(e, coord, mp->_index);
    }
    else
    {
      NP = steadyStateCracking(e, coord, mp->_index);
    }

    // Computes the creep strain rate 
    eCreepRate = Eo*exp(QT / (Rg*To))*exp((-1.0)*QT / (Rg*TEMP))* pow(MISES / So, NP);

    return eCreepRate;
  }

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

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

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

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

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

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

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

#endif