gmpMaterialIsotropicTMDamage.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_COUPLED_TM_MATERIALISOTROPICTMDAMAGE_H_
#define _GEMA_PLUGIN_COUPLED_TM_MATERIALISOTROPICTMDAMAGE_H_
 
#include "gmpCoupledTMConfig.h"
#include "gmpCoupledTMMaterial.h"
 
#include <gmpFemPhysics.h>
#include <gmMathUtils.h>
#include <gmpMechanicPoint.h>
#include <gmpFemPhysicsCommonMaterial.h>
 
 
class GMP_COUPLED_TM_PHYSICS_API_EXPORT GmpMaterialIsotropicTMDamage : public GmpCoupledTMMaterial
{
protected:
  enum ElementPropertyIds
  {   
    TDAMLAW_ID, 
    ALPHA_ID,   
    TD0_ID,     
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
  enum DamageGaussAttrIds
  {
    TDAM_GA_ID,         
    TDAMOLD_GA_ID,      
    TISV_GA_ID,         
    TISVOLD_GA_ID,      
 
    // ------    NO ADDING BELOW THIS LINE
    NUM_GA_IDS,          
  };
 
  enum ThermalDamageLaw
  {
    pearce,     
    luccioni,   
 
    // --- NO ADDING BELOW THIS LINE
    NUM_THERMALDAMAGELAW_IDS,  
  };
 
public:
 
  GmpMaterialIsotropicTMDamage(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpCoupledTMMaterial(typeIndex, typeName, logger) {}
 
  virtual ~GmpMaterialIsotropicTMDamage() {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialIsotropicTMDamage(typeIndex, typeName, logger);
  }
 
  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap();
  
  // Checks material loaded data. 
  virtual bool checkLoadedData(const GmElement* e) const;
 
  // Sets the initial conditions required by isotropic damage material
  virtual bool setInitialConditions(const GmElement* e, GmpMechanicPoint* mp, const GmVector* coord, unsigned sc) const;
 
  // Returns the stresses according to the material behavior adopted
  virtual bool thermoMechanicalConstitutiveModel(const GmElement* e, GmMatrix& Ded, const GmpMechanicPoint* mp, const GmVector* coord, 
          const GmVector& Time, double Tg, unsigned nc, bool ips) const;
  
  virtual double thermalExpansionCoeff(const GmElement* e, const GmVector* coord, int ip) const
  {
    return  propertyAc(ALPHA_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double thermalDamageThreshold(const GmElement* e, const GmVector* coord, int ip) const
  {
    return  propertyAc(TD0_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual int thermalDamageLaw(const GmElement* e, const GmVector* coord, int ip) const //alterei
  {
    S_TRACE(); assert(e);
    GmCellAccessor* TDLawAcc = propertyAc(TDAMLAW_ID);
    if (TDLawAcc == NULL)
    {
      return 0;
    }
    return TDLawAcc->scalarValueAt(e, coord, ip);
  }
 
  // Thermal damage Law definition
  virtual bool thermalDamageCriterion(const GmElement*, const GmpMechanicPoint*, const GmVector*, unsigned, double, double& rtd, unsigned) const;
  //virtual bool thermalDamageCriterion(const GmElement*, const GmpMechanicPoint*, const GmVector*, unsigned, double, double, double, GmVector&, unsigned) const;
  
 
  // Thermal damage Law definition
  virtual bool thermalDamageLaw(const GmElement* e, const GmVector* coord, const GmpMechanicPoint* mp, unsigned sc, double& g, double rtd, double& dgt, double Tg) const;
  
};
 
#endif