gmpInterfaceFlow.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_HYDRAULIC_INTERFACEFLOW_H_
#define _GEMA_PLUGIN_HYDRAULIC_INTERFACEFLOW_H_
 
#include "gmpHydraulicInterface.h"
#include <gmTrace.h>
#include <math.h>
 
 
template <class T> class GmpInterfaceFlow : public T
{
public:
  GmpInterfaceFlow(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 GmpInterfaceAxFlow : public T
{
public:
  GmpInterfaceAxFlow(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);
    //int    n      = Bp.n_cols;
    //double rGauss = 0.0;
    //for(int j = 0; j < n; ++j)
    //  rGauss += N(j)*X(0, j);
 
    double detJ = shape->shapeCartesianPartialsFromJacobian(ncoord, J, Bp, true);
 
    //for(int k = 0; k < n; ++k)
    //  Bp(0, k) = Bp(0, k)/rGauss;
 
    return detJ;
  }
 
  // Comments on the base class
  virtual double axisymmetricFactor(const GmElement* e, const GmMatrix& MX, const GmVector& N) const
  { 
    S_TRACE();
    assert(e);
    // number of nodes
    int n = e->numNodes();
    int d = GmpFemPhysicsCommon::nodeDim();
    // number of nodes for tri-nodded interface elements
    if (e->type() == GM_INT2DL6)
    {
      n = n - pow(2, d - 1);
    }
    assert(MX.n_rows == n / 2 && MX.n_cols == d);
 
    double rGauss = 0.0;
    for (int j = 0; j < n / 2; ++j)
      rGauss += N(j)*MX(j, 0);
 
    return 2 * arma::datum::pi * rGauss;
  }
};
 
#endif