gmpMaterialDPCreep.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_DPCREEP_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_DPCREEP_H_

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

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

// Double power creep material
class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialDPCreep : public GmpMaterialCreep
{
protected:
  enum ElementPropertyIds
  {
    A1_ID = GmpMaterialCreep::ElementPropertyIds::NUM_PROPERTY_IDS, 
    AN1_ID,        
    Q1_ID,         
    A2_ID,         
    AN2_ID,        
    Q2_ID,         
    So_ID,         
    R_ID,          

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

public:

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

  virtual ~GmpMaterialDPCreep() {}

  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialDPCreep(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
      m["properties"] = m["properties"].value<GmpFemPhysicsCommon::ValueList>()
        << GmpFemPhysicsCommon::ScalarValue(A1_ID,  "A1", QObject::tr("Structure factor for DCL"), "1/s", true)
        << GmpFemPhysicsCommon::ScalarValue(AN1_ID, "n1", QObject::tr("Stress power for DCL"), "", true)
        << GmpFemPhysicsCommon::ScalarValue(Q1_ID,  "Q1", QObject::tr("Thermal activation energy for DCL"), "J/mol", true)
        << GmpFemPhysicsCommon::ScalarValue(A2_ID,  "A2", QObject::tr("Structure factor for UMC"), "1/s", true)
        << GmpFemPhysicsCommon::ScalarValue(AN2_ID, "n2", QObject::tr("Stress power for UMC"), "", true)
        << GmpFemPhysicsCommon::ScalarValue(Q2_ID,  "Q2", QObject::tr("Thermal activation energy for UMC"), "J/mol", true)
        << GmpFemPhysicsCommon::ScalarValue(So_ID,  "So", QObject::tr("Threshold deviatoric stress"), "kPa", true)
        << GmpFemPhysicsCommon::ScalarValue(R_ID,   "R",  QObject::tr("Universal gas constant"), "J/(mol*K)", true);
    }
    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, A1, Q1, N1, A2, Q2, N2, So, R, TEMP;

    // Reads Double Power creep Properties 
   
     A1  = structureFactorDCL(e, coord, mp->_index);
     Q1  = thermalActivEnergyDCL(e, coord, mp->_index);
     N1 = stressPowerDCL(e, coord, mp->_index);
     A2  = structureFactorUMC(e, coord, mp->_index);
     Q2  = thermalActivEnergyUMC(e, coord, mp->_index);
     N2 = stressPowerUMC(e, coord, mp->_index);
     So  = thresholdDevStress(e, coord, mp->_index);
     R   = universalGasConstant(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);
        }
  // Computes the creep strain rate 
  eCreepRate = A1 * exp(-Q1 / (R*TEMP))*pow(MISES / So, N1) + A2 * exp(-Q2 / (R*TEMP))*pow(MISES / So, N2);

  return eCreepRate;
  }
    
  double structureFactorDCL(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(A1_ID)->scalarValueAt(e, coord, ip);
  }
  
  double thermalActivEnergyDCL(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(Q1_ID)->scalarValueAt(e, coord, ip);
  }
  
  double stressPowerDCL(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(AN1_ID)->scalarValueAt(e, coord, ip);
  }
  
  double structureFactorUMC(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(A2_ID)->scalarValueAt(e, coord, ip);
  }

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

  double stressPowerUMC(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(AN2_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 universalGasConstant(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(R_ID)->scalarValueAt(e, coord, ip);
  }

};

#endif