gmpMaterialCohesiveLinearSoftening.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_COHESIVELINEARSOFTENING_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_COHESIVELINEARSOFTENING_H_


#include "gmpMaterialElasticInterface.h"
#include "gmpMechanicPoint.h"
#include "gmpFemPhysics.h"
#include "gmMathUtils.h"

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpCohesiveLinearSoftening: public GmpElasticInterface
{
protected:
  enum CohesiveLinearPropertyIds
  {              
    Ed_ID = GmpElasticInterface::NUM_PROPERTY_IDS, 
    Sni_ID,      
    Tsi_ID,      
    Tti_ID,      
    DamIni_ID,   


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

  enum damageCriteria
  {
    quadS,        
    maxS,         


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

  // Check the damage initiation criteria
  virtual double damageInitiationCriteria(const GmElement* e, GmVector stress, const GmVector* coord, int ip, int nc) const;

  // returns the scalar damage value "Sdvg"
  virtual double scalarDamageVariable(const GmElement* e, const GmVector* coord, int ip, double de, double di, double df) const;

  // returns the gradient of the damage variable in relation to the equivalente displacement
  // Hd = d(Sdvg)/d(del_e)
  virtual double deriveDamage(const GmElement* e, const GmVector* coord, int ip, double de, double di, double df) const;

public:
  GmpCohesiveLinearSoftening(int typeIndex, QString typeName, const GmLogCategory& logger)
    :GmpElasticInterface(typeIndex, typeName, logger){}

  virtual ~GmpCohesiveLinearSoftening() {}

  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpCohesiveLinearSoftening(typeIndex, typeName, logger);
  }

  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap();

  enum cohesiveGaussAttributeIds
  {
    SDV_GA_ID = GmpElasticInterface::NUM_GA_IDS, 
    SDVOLD_GA_ID,        
    IV_GA_ID,            
    IVOLD_GA_ID,         
    AFI_GA_ID,           


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

  // Returns the stresses according to the material behavior adopted
  virtual bool mechanicalConstitutiveModel(const GmElement* e, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned nc) const;

  //Sets the initial conditions required by interface materials
  virtual bool setInitialConditions(const GmElement* e, GmpMechanicPoint* mp, const GmVector* coord, unsigned nc) const;

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

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

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

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

  virtual damageCriteria damageCriterion(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    GmpCohesiveLinearSoftening::damageCriteria type;

    if (propertyAc(DamIni_ID) == NULL)
    {
      return GmpCohesiveLinearSoftening::damageCriteria::quadS;
    }

    int criterionType = propertyAc(DamIni_ID)->scalarValueAt(e, coord, ip);
    //
    if (criterionType == 1)
    {
      type = GmpCohesiveLinearSoftening::damageCriteria::maxS;
    }
    else
    {
      type = GmpCohesiveLinearSoftening::damageCriteria::quadS;
    }

    return type;
  }

};
#endif