gmpMaterialWDCreep.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_WDCREEP_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_WDCREEP_H_
 
#include "gmpMechanicalConfig.h"
#include <gmpFemPhysics.h>
#include <gmMathUtils.h>
 
#include <gmpMaterialElastic.h>
#include "gmpMechanicPoint.h"
 
class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialWDCreep : public GmpMechanicalMaterialElastic
{
protected:
  enum ElementPropertyIds
  {
    A_ID = GmpMechanicalMaterialElastic::NUM_PROPERTY_IDS,  
    B_ID,         
    C1_ID,      
    C2_ID,      
    N_ID,        
    F_ID,         
    R_ID,         
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
  enum GaussAttributeIds
  {
    EC_GA_ID,      
    ECOLD_GA_ID,   
    SDV_GA_ID,     
    SDVOLD_GA_ID,  
    SDV2_GA_ID,    
    SDV3_GA_ID,    
    SDV4_GA_ID,    
    SDV5_GA_ID,    
    SDV6_GA_ID,    
    SDV7_GA_ID,    
    SDV8_GA_ID,    
    SDV9_GA_ID,    
    SDV10_GA_ID,   
    SDV11_GA_ID,   
    SDV12_GA_ID,   
    SDV13_GA_ID,   
    SDV14_GA_ID,   
    SDV15_GA_ID,   
    SDV16_GA_ID,   
    SDV17_GA_ID,   
    SDV18_GA_ID,   
    SDV19_GA_ID,   
    SDV20_GA_ID,   
    SDV21_GA_ID,   
    SDV22_GA_ID,   
    SDV23_GA_ID,   
    SDV24_GA_ID,   
    
                   
    // ------    NO ADDING BELOW THIS LINE
    NUM_GA_IDS,  
  };
 
public:
 
  GmpMaterialWDCreep(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpMechanicalMaterialElastic(typeIndex, typeName, logger) {}
 
  virtual ~GmpMaterialWDCreep() {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialWDCreep(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; }
 
  virtual double firstprimarycreepConstant(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(A_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double secondprimarycreepConstant(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(B_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double secondarycreepConstant(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(C1_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double tertiarycreepConstant(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(C2_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double stressExponential(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(N_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double firstdamageConstant(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(F_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double seconddamageConstant(const GmElement* e, const GmVector* coord, int ip) const
  {
          return propertyAc(R_ID)->scalarValueAt(e, coord, ip);
  }
 
  // Computes the elastic stiffness using lame and G
  void elasticStiffness(const GmElement*, unsigned, GmMatrix&, const GmVector*, unsigned, int) const;
  
};
 
#endif