gmpMaterialElasticInterface.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_ELASTICINTERFACE_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_ELASTICINTERFACE_H_


#include "gmpInterfaceMaterial.h"
#include "gmpMechanicPoint.h"
#include <gmpFemPhysics.h>
#include <gmMathUtils.h>
#include <gmTrace.h>

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpElasticInterface: public GmpInterfaceMaterial
{
protected:
  enum InterfaceElementPropertyIds
  {
    Kni_ID,        
    Ksi_ID,        
    Kti_ID,        
    GAP_ID,        
    Closure_ID,    

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

  enum closureModel
  {
    normal,          
    hardStiffness,   
    bartonBandis,    


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

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

  virtual ~GmpElasticInterface() {}

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

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

  enum elasticInterfaceGaussAttrIds
  {
    VM_GA_ID,    

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

  virtual bool fillElasticStiffnessMatrix(const GmElement* e, GmMatrix& Dep, const GmVector* coord, int ip, double penalty = 1.0) const;
  
  virtual bool fillBartonBandisStiffnessMatrix(const GmElement* e, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord)const;
  
  virtual bool fillBartonBandisStresses(const GmElement* e, GmVector& Snew, const GmpMechanicPoint* mp, const GmVector* coord, GmVector& eNew, GmVector& eOld, GmVector& sOld)const;

  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 normalElasticStiffness(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(Kni_ID)->scalarValueAt(e, coord, ip);
  }

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

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

  virtual double fractureOpening(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    if (propertyAc(GAP_ID) == NULL)
    {
      return 0.0;
    }

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

  virtual closureModel fractureClosureModel(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    GmpElasticInterface::closureModel type;

    if (propertyAc(Closure_ID) == NULL)
    {
      return GmpElasticInterface::closureModel::normal;
    }

    int closureType = propertyAc(Closure_ID)->scalarValueAt(e, coord, ip);
    //
    if (closureType == 1)
    {
      type = GmpElasticInterface::closureModel::hardStiffness;
    }
    else if (closureType == 2)
    {
      type = GmpElasticInterface::closureModel::bartonBandis;
    }
    else
    {
      type = GmpElasticInterface::closureModel::normal;
    } 

    return type;
  }
};
#endif