gmpMaterialCrackingCrushingDamage.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_CRACKINGCRUSHINGDAMAGE_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_CRACKINGCRUSHINGDAMAGE_H_

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

#include "gmpMaterialElastic.h"
#include "gmpMechanicPoint.h"


class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialCrackingCrushingDamage : public GmpMechanicalMaterialElastic
{
protected:
  enum ElementPropertyIds
  {   
    FCM_ID = GmpMechanicalMaterialElastic::NUM_PROPERTY_IDS, 
    RCM_ID,      
    FT0_ID,      
    FBT_ID,      
    FTC_ID,      
    FC0_ID,      
    FBC_ID,      
    RHO_ID,      
    GF_ID,       
    GC_ID,       
    LCH_ID,      
    LCHMIN_ID,   

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

public:

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

  virtual ~GmpMaterialCrackingCrushingDamage() {}

  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialCrackingCrushingDamage(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;

  enum elasticDamageGaussAttrIds
  {
    DAM_GA_ID,         
    DAMOLD_GA_ID,      
    TISV_GA_ID,        
    TISVOLD_GA_ID,     
    CISV_GA_ID,        
    CISVOLD_GA_ID,     

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

  // 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 mechanicalConstitutiveModel(const GmElement* e, GmMatrix& Ded, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector& Time, unsigned nc, bool ips) const;

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

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

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

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

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

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

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

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

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

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

  virtual double userDefinedCharacteristicLength(const GmElement* e, const GmVector* coord, int ip) const
  {
      S_TRACE(); assert(e);
      // Get accessor for the user-defined characteristic length
      GmCellAccessor* lch = propertyAc(LCH_ID);

      if (lch == NULL)
      {
          return 0.0;
      }

      return lch->scalarValueAt(e, coord, ip);
  }

  virtual double minimumCharacteristicLength(const GmElement* e, const GmVector* coord, int ip) const
  {
      S_TRACE(); assert(e);
      // Get accessor for the minimum characteristic length
      GmCellAccessor* lchmin = propertyAc(LCHMIN_ID);

      if (lchmin == NULL)
      {
          return 0.0;
      }

      return lchmin->scalarValueAt(e, coord, ip);
  }

  // Damage criterion definition
  virtual bool damageCriterion(const GmElement*, const GmpMechanicPoint*, const GmVector*, unsigned, GmVector, double&, double&, GmVector&, GmVector&, unsigned) const;

  // Damage evolution law definition
  virtual bool damageLaw(const GmElement*, const GmVector*, const GmpMechanicPoint*, double, double, double&, double&, double& dg) const;

  // Returns the characteristic length
  virtual double characteristicLength(const GmElement*, const GmVector*, int, GmValueAccessor*, bool type = false) const;
};

#endif