gmInt2DShape.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_INT2D_SHAPE_H_
#define _GEMA_INT2D_SHAPE_H_
 
#include "gmShape.h"
 
//--------------------------------------------------------------------------------------------------------------------------
//      2D Interface Elements
//--------------------------------------------------------------------------------------------------------------------------
 
class GMC_API_EXPORT GmInt2DShape : public GmShape
{
public:
  // See comments on the base class
  virtual int numNaturalCoord() const { return 1; }
 
  // See comments on the base class
  virtual int numCartesianCoord() const { return 2; }
 
  // See comments on the base class
  virtual void naturalCoordLimits(int coord, double* min, double* max) const { Q_UNUSED(coord); *min = -1.0; *max = 1.0; }
 
  // See comments on the base class
  virtual void naturalCenter(GmVector& coord) const { coord.zeros(1); }
 
  // See comments on the base class. 
  virtual bool translateEdgePoint(int edge, const GmVector& srcEdgeCoord, GmVector& elementCoord) const
  {
    Q_UNUSED(edge); 
    elementCoord = srcEdgeCoord;   // edge and element coordinates are the same!!!
    return !(edge % 2);  // Returns true for edges 0 & 2 and false for edges 1 & 3, since those should not be used for integration
  }
 
  // See comments on the base class. Undefined for 2D interface elements.
  virtual bool translateFacePoint(int face, const GmVector& srcFaceCoord, GmVector& elementCoord) const
  {
    Q_UNUSED(face); Q_UNUSED(srcFaceCoord); Q_UNUSED(elementCoord); return false;
  }
 
  virtual double shapeCartesianPartialsFromPJ(const GmMatrix& P, const GmMatrix& J, GmMatrix& dN, bool transposed = false) const;
 
  // See comments on the base class
  virtual double scaledJacobianDet(const GmMatrix& J) const 
  { 
    assert(J.n_elem == 2);
    // Get dx/dxi and dy/dxi derivate. Remember that J might be transposed
    double dx = J(0, 0); 
    double dy = (J.n_cols == 2) ? J(0, 1) : J(1, 0); 
    return sqrt(dx*dx + dy*dy);
  }
 
  // See comments on the base class
  virtual void jacobianAndPartials(const GmVector& ncoord, const GmMatrix& X, GmMatrix& J,
                                   GmMatrix& P, bool transposed = false) const;
  
  // See comments on the base class
  virtual void jacobianRotation(const GmMatrix& J, GmMatrix& R) const { jacRotation2D(J, R);}
 
  // Undefined
  virtual double borderScalingFactor(int border, const GmVector& borderCoord, const GmVector& elementCoord,
                                     const GmMatrix& X, bool transposed = false) const
  {
    Q_UNUSED(border); Q_UNUSED(borderCoord); Q_UNUSED(elementCoord); Q_UNUSED(X); Q_UNUSED(transposed);
    return 0.0;
  }
 
  // Undefined
  virtual double edgeScalingFactor(int border, const GmVector& borderCoord, const GmVector& elementCoord,
                                   const GmMatrix& X, bool transposed = false) const
  {
    Q_UNUSED(border); Q_UNUSED(borderCoord); Q_UNUSED(elementCoord); Q_UNUSED(X); Q_UNUSED(transposed);
    return 0.0;
  }
 
  // Undefined
  virtual double faceScalingFactor(int border, const GmVector& borderCoord, const GmVector& elementCoord,
                                   const GmMatrix& X, bool transposed = false) const
  {
    Q_UNUSED(border); Q_UNUSED(borderCoord); Q_UNUSED(elementCoord); Q_UNUSED(X); Q_UNUSED(transposed);
    return 0.0;
  }
 
  // See comments on the base class
  void naturalToCartesian(const GmVector& ncoord, const GmMatrix& X, GmVector& coord) const;
 
protected:
  // See comments on the base class. At the moment we don't have a working version for interface elements
  // Can we just reuse the Bar implementation for all elements?
  virtual bool gradientBasedCartesianToNatural(const GmVector& coord, const GmMatrix& X, double tol,
                                               int maxIter, double natTol, GmVector& ncoord, bool* inside) const
  {
    Q_UNUSED(coord); Q_UNUSED(X); Q_UNUSED(tol); Q_UNUSED(maxIter); Q_UNUSED(natTol); Q_UNUSED(ncoord); Q_UNUSED(inside);
    return false;
  }
 
  virtual void midPlaneMatrix(const GmMatrix& X, GmMatrix& MX, bool transposed) const = 0;
};
 
class GMC_API_EXPORT GmInt2DL4Shape : public GmInt2DShape
{
public:
  // See comments on the base class
  virtual GmCellType elemType() const { return GM_INT2DL4; }
 
  // Returns the number of shape functions of this cohesive element (equal to the node number at middle plane)
  virtual int numFunctions() const { return 2; }
 
  // See comments on the base class
  virtual void nodeNaturalCoord(int node, GmVector& coord) const;
 
  virtual void shapeValues(const GmVector& ncoord, GmVector& N) const;
  virtual void shapePartials(const GmVector& ncoord, GmMatrix& dN, bool transposed = false) const;
 
protected:
  // See comments on the base class
  virtual bool hasGeometryBasedCartesianToNatural() const { return true; }
 
  virtual bool geometryBasedCartesianToNatural(const GmVector& coord, const GmMatrix& X, double natTol,
                                               GmVector& ncoord, bool* inside) const;
 
  virtual void midPlaneMatrix(const GmMatrix& X, GmMatrix& MX, bool transposed) const;
};
 
class GMC_API_EXPORT GmInt2DL6Shape : public GmInt2DL4Shape
{
public:
  // See comments on the base class
  virtual GmCellType elemType() const { return GM_INT2DL6; }
};
 
class GMC_API_EXPORT GmInt2DQ6Shape : public GmInt2DShape
{
public:
  // See comments on the base class
  virtual GmCellType elemType() const { return GM_INT2DQ6; }
 
  // Returns the number of shape functions of this interface element (equal to the node number at middle plane)
  virtual int numFunctions() const { return 3; }
 
  // See comments on the base class
  virtual void nodeNaturalCoord(int node, GmVector& coord) const;
 
  virtual void shapeValues(const GmVector& ncoord, GmVector& N) const;
  virtual void shapePartials(const GmVector& ncoord, GmMatrix& dN, bool transposed = false) const;
protected:
   virtual void midPlaneMatrix(const GmMatrix& X, GmMatrix& MX, bool transposed) const;
};
 
class GMC_API_EXPORT GmInt2DQ8Shape : public GmInt2DQ6Shape
{
public:
  // See comments on the base class
  virtual GmCellType elemType() const { return GM_INT2DQ8; }
};
 
#endif