gmpMaterialEquivalentContinuum.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2016 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_EQUIVALENTCONTINUUM_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_EQUIVALENTCONTINUUM_H_

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

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMechanicalMaterialEquivalentContinuum : public GmpMechanicalMaterialElastic
{
protected:
  enum FracRockElementPropertyIds
  { 
    Knf_ID = GmpMechanicalMaterialElastic::NUM_PROPERTY_IDS,  
    Ksf_ID,        
    Sfr_ID,        
    DIP_ID,        
    STRIKE_ID,     
    Bfr_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:

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

  virtual ~GmpMechanicalMaterialEquivalentContinuum() {}

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

  enum fracturedRockGaussAttrIds
  {
    VMFR_GA_ID,   
    BFR_GA_ID,    

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

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

  // Sets the initial conditions required by fractured rock 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& Dep, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector&, unsigned nc, bool ips) const;

  // calculates derived results
  virtual bool calcDerivedResults(const GmElement* e, GmpMechanicPoint* mp, const GmVector* coord, unsigned sc) const;
  
  virtual GmVector normalElasticStiffnessFracture(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector Kn;
    GmCellAccessor* knAc = propertyAc(Knf_ID);
    // Fills property values
    vectorPropertyValue(e, knAc, Kn, coord, ip);
    return Kn;
  }

  virtual GmVector shearElasticStiffnessFracture(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector Ks;
    GmCellAccessor* ksAc = propertyAc(Ksf_ID);
    // Fills property values
    vectorPropertyValue(e, ksAc, Ks, coord, ip);
    return Ks;
  }

  virtual GmVector fractureSpacing(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector Sf;
    GmCellAccessor* sfAc = propertyAc(Sfr_ID);
    // Fills property values
    vectorPropertyValue(e, sfAc, Sf, coord, ip);
    return Sf;
  }

  virtual GmVector fracturedDip(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector Dip;
    GmCellAccessor* dipAc = propertyAc(DIP_ID);
    // Fills property values
    vectorPropertyValue(e, dipAc, Dip, coord, ip);
    return Dip;
  }

  virtual GmVector fractureStrike(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector Strike;
    GmCellAccessor* strikeAc = propertyAc(STRIKE_ID);
    if (!vectorPropertyValue(e, strikeAc, Strike, coord, ip)) {
      Strike.set_size(1);
      return Strike.fill(0);
    }
    return Strike;
  }

  virtual GmVector fractureAperture(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector Aperture;
    GmCellAccessor* apertureAc = propertyAc(Bfr_ID);
    // Fills property values
    vectorPropertyValue(e, apertureAc, Aperture, coord, ip);

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

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

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

    return type;
  }
  
  virtual GmVector updatedFractureAperture(const GmElement* e, const GmVector* coord, int ip) const;

  virtual void stiffness(const GmElement*, unsigned, GmMatrix&, GmMatrix& Cf, const GmVector*, unsigned, int) const;

  virtual void factureFlexibility(GmMatrix& Cf, double Kn, double Ks, double Sf, int ns) const;

};

#endif