gmpMaterialCohesiveExponentialSoftening.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_MECHANICALMATERIAL_COHESIVEEXPONENTIALSOFTENING_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_COHESIVEEXPONENTIALSOFTENING_H_


#include "gmpMaterialCohesiveLinearSoftening.h"
#include "gmpMaterialElasticInterface.h"
#include "gmpMechanicPoint.h"
#include "gmpFemPhysics.h"
#include "gmMathUtils.h"
#include <math.h>

class GmpCohesiveExponentialSoftening : public GmpCohesiveLinearSoftening
{
protected:
  enum CohesiveExpElemPropertyIds
  {              
    ALPHA_ID = GmpCohesiveLinearSoftening::NUM_PROPERTY_IDS,

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

  // returns the scalar damage value "Sdvg"
  virtual double scalarDamageVariable(const GmElement* e, const GmVector* coord, int ip, double de, double di, double df) const
  {
    S_TRACE(); 
    double alpha = exponentialSofteningParameter(e, coord, ip);
    double Sdvg = 1.0 - di / de * (1.0 - (1.0 - exp(-alpha*((de - di) / (df - di)))) / (1.0 - exp(-alpha)));
    return Sdvg;
  }

  // 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
  {
    S_TRACE();
    double alpha = exponentialSofteningParameter(e, coord, ip);
    double Hd = di*(1 - 1 / (1 - exp(-alpha))) / pow(de, 2.0) + (di / pow(de, 2.0) + di*alpha / (de*(df - di)))*
                exp(-alpha*(de - di) / (df - di)) / (1 - exp(-alpha));
    return Hd;
  }

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

  virtual ~GmpCohesiveExponentialSoftening() {}

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

  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap()
  {
    S_TRACE();
    static QVariantMap m;
    if (m.isEmpty())
    {
      // read the material attribute map from elastic interface material
      m = *GmpCohesiveLinearSoftening::materialMetaDataMap();
      // Adds material properties
      GmpFemPhysicsCommon::ValueList matExpProp = m["properties"].value<GmpFemPhysicsCommon::ValueList>();
      matExpProp.append(GmpFemPhysicsCommon::ScalarValue(ALPHA_ID, "alpha", QObject::tr("Exponential material parameter"), "", true));

      m["properties"] = matExpProp;
    }   

    return &m;
  }

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

};
#endif