gmpHydraulicFlow.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_HYDRAULICFLOW_H_
#define _GEMA_PLUGIN_HYDRAULICFLOW_H_
 
#include "gmpHydraulic.h"
#include <gmTrace.h>
#include <math.h>
 
 
template <class T> class GmpHydraulicFlow : public T
{
public:
  GmpHydraulicFlow(const char* pluginType, GmSimulationData* simulation, QString id, QString description, 
                   const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger)
    : T(pluginType, simulation, id, description, matFactory, logger)
  {}
 
protected:
  // Comments on the base class
  virtual double fillBpMatrix(const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bp)
  {
    S_TRACE();
    Q_UNUSED(X);  Q_UNUSED(N);
 
    return shape->shapeCartesianPartialsFromJacobian(ncoord, J, Bp, true);
  }
};
 
template <class T> class GmpHydraulicAxFlow : public T
{
public:
    GmpHydraulicAxFlow(const char* pluginType, GmSimulationData* simulation, QString id, QString description, 
                       const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger)
      : T(pluginType, simulation, id, description, matFactory, logger)
    {}
 
protected:
  virtual bool isAxisymmetric() { return true; }
 
  // Comments on the base class
  virtual double fillBpMatrix(const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bp)
  {
    S_TRACE();
    Q_UNUSED(X);  Q_UNUSED(N);
    return shape->shapeCartesianPartialsFromJacobian(ncoord, J, Bp, true);
  }
 
  // Comments on the base class
  virtual double axisymmetricFactor(const GmElement* e, const GmMatrix& X, const GmVector& N) const
  {
    S_TRACE();
    int n = e->numNodes();
    // evaluate radio at integration points
    double rGauss = 0.0;
    for (int j = 0; j < n; ++j)
      rGauss += N(j)*X(j, 0);
 
    return 2 * arma::datum::pi * rGauss;
  }
};
 
template <class T> class GmpHydraulicBrinkman : public T
{
public:
  GmpHydraulicBrinkman(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                       const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger)
    : T(pluginType, simulation, id, description, matFactory, logger)
  {}
 
protected:
  // Comments on the base class
  virtual double fillBpMatrix(const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bp)
  {
    S_TRACE();
    Q_UNUSED(X);  Q_UNUSED(N);
 
    // fill the elemental shape cartesian partials
    double detJ = shape->shapeCartesianPartialsFromJacobian(ncoord, J, Bp, true);
 
    return detJ;
  }
};
 
#endif