gmpMaterialCohesivePotentialSoftening.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_COHESIVEPOTENTIALSOFTENING_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_COHESIVEPOTENTIALSOFTENING_H_

#include <gmpFemPhysics.h>
#include <gmTrace.h>

#include "gmpInterfaceMaterial.h"
#include "gmpMechanicPoint.h"


class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialCohesivePotentialSoftening : public GmpInterfaceMaterial
{
public:
  GmpMaterialCohesivePotentialSoftening(int typeIndex, QString typeName, const GmLogCategory& logger)
    :GmpInterfaceMaterial(typeIndex, typeName, logger) {}

  virtual ~GmpMaterialCohesivePotentialSoftening() {}

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

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

  virtual bool mechanicalConstitutiveModel(const GmElement* e, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, unsigned nc) const;

protected:
  enum CohesiveLinearPropertyIds
  {              
    GN_ID,      
    GS_ID,      
    Sni_ID,     
    Tsi_ID,     
    Tti_ID,     
    ALPHA_ID,   
    BETA_ID,    
    Ln_ID,      
    Ls_ID,      

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

  enum cohesiveGaussAttributeIds
  {
    DELMAX_GA_ID,          
    DELMAXOLD_GA_ID,       


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

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

  virtual double shearEnergy(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(GS_ID)->scalarValueAt(e, coord, ip);
  }
  
  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 normalShape(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(ALPHA_ID)->scalarValueAt(e, coord, ip);
  }

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

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

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

#endif