gmpMaterialFluid.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_FLUID_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_FLUID_H_

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

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

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialFluid : public GmpMechanicalMaterial
{
protected:
  enum ElementPropertyIds
  {
    Bf_ID,          

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

public:

  GmpMaterialFluid(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpMechanicalMaterial(typeIndex, typeName, logger) {}

  virtual ~GmpMaterialFluid() {}

  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap()
  {
    S_TRACE();
    static QVariantMap m;
    if (m.isEmpty())
    {
      m["properties"] = GmpFemPhysicsCommon::ValueList()
        << GmpFemPhysicsCommon::ScalarValue(Bf_ID, "Bf", QObject::tr("Fluid compressibility"), "1/kPa", true);
    }
    return &m;
  }

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

  // Returns the stresses according to the material behavior adopted
  virtual bool mechanicalConstitutiveModel(const GmElement* e, GmMatrix& De, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector& Time, unsigned sc, bool ips) const
  {
    S_TRACE();
    Q_UNUSED(Time);
    assert(e);

    //stress components
    const unsigned ns = GmMathUtils::bit_enum_count(sc, (unsigned)GmpMechanicPoint::stress::last);
    // Integration point
    int ip = mp->_index;
    // Defines scalar varibles
    GmVector s(ns), dStrain(ns);

    //old strain
    const GmVector eu_old(mp->_oldStrain->valueAt(e, ip, coord), ns);
    // old stress
    const GmVector s_old(mp->_oldStress->valueAt(e, ip, coord), ns);
    //residual stress
    const GmVector sr(mp->_residualStress->valueAt(e, ip, coord), ns);

    //new strain
    const GmVector eu(mp->_newStrain->valueAt(e,ip, coord), ns);

    // Fluid Compressibility
    double Cf = fluidCompressibility(e, coord, ip);
    // Bulk modulus
    double K = 1/Cf;

    // Fills the elastic matrix;
    GmpMechanicUtils::elasticStiffnessKG(De, K, 0.0, sc, ips);

    // Strain increment
    dStrain = eu - eu_old;

    //Updates of stresses
    s = s_old + De*dStrain + sr;

    //Updates state variables
    mp->_newStress->setValue(e, ip, s.mem);
    mp->_newYSR->setScalarValue(e, ip, 1.0);
    return true;
  }

  virtual bool isIsotropic() const { return true; }

  virtual double fluidCompressibility(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    double Cf = propertyAc(Bf_ID)->scalarValueAt(e, coord, ip);

    // Computes Bulk modulus
    if (Cf == 0.0)
    {
      return 1.0e-25;
    }

    return Cf;
  }

};

#endif