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 int ns = GmpMechanicUtils::getNumStressComp(sc);
    // Integration point
    int ip = mp->_index.localData();
    // 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