QuadBase9N.h

// -*-c++-*-
//----------------------------------------------------------------------------
//  XC program; finite element analysis code
//  for structural analysis and design.
//
//  Copyright (C)  Luis C. Pérez Tato
//
//  This program derives from OpenSees <http://opensees.berkeley.edu>
//  developed by the  «Pacific earthquake engineering research center».
//
//  Except for the restrictions that may arise from the copyright
//  of the original program (see copyright_opensees.txt)
//  XC is free software: you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or 
//  (at your option) any later version.
//
//  This software is distributed in the hope that it will be useful, but 
//  WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details. 
//
//
// You should have received a copy of the GNU General Public License 
// along with this program.
// If not, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------------
//QuadBase9N.h
                                                                        
#include "PlaneElement.h"

#ifndef QuadBase9N_h
#define QuadBase9N_h

#include "preprocessor/multi_block_topology/matrices/ElemPtrArray3d.h"
#include "preprocessor/multi_block_topology/aux_meshing.h"
#include "preprocessor/prep_handlers/LoadHandler.h"
#include "domain/load/plane/BidimStrainLoad.h"
#include "vtkCellType.h"

namespace XC {
template <class PhysProp>
class QuadBase9N: public PlaneElement<9,PhysProp>
  {
  protected:
    ElemPtrArray3d put_on_mesh(const NodePtrArray3d &,meshing_dir dm) const;

  public:
    QuadBase9N(int classTag);
    QuadBase9N(int tag, int classTag,const PhysProp &);

    BoolArray3d getNodePattern(void) const;
    Element::NodesEdge getNodesEdge(const size_t &i) const;
    ID getLocalIndexNodesEdge(const size_t &i) const;
    int getEdgeNodes(const Node *,const Node *) const;
    Polygon3d getPolygon(bool initialGeometry= true) const;
    Segment3d getSide(const size_t &i,bool initialGeometry= true) const;

    int getVtkCellType(void) const;
    std::vector<int> getIdxNodes(void) const;

    void zeroLoad(void);    
    int addLoad(ElementalLoad *theLoad, double loadFactor);
  };

template <class PhysProp>
XC::QuadBase9N<PhysProp>::QuadBase9N(int classTag)
  : PlaneElement<9,PhysProp>(0,classTag) {}

template <class PhysProp>
XC::QuadBase9N<PhysProp>::QuadBase9N(int tag,int classTag,const PhysProp &pp)
  :PlaneElement<9,PhysProp>(tag,classTag,pp) {}

template <class PhysProp>
BoolArray3d XC::QuadBase9N<PhysProp>::getNodePattern(void) const
  {
    BoolArray3d retval(1,3,3,true); //One layer, three rows, three columns.
    return retval;
  }

template <class PhysProp>
XC::ElemPtrArray3d XC::QuadBase9N<PhysProp>::put_on_mesh(const NodePtrArray3d &nodes,meshing_dir dm) const
  { return put_quad9N_on_mesh(*this,nodes,dm); }

//
//   3   6    2
//   +---+---+
//   |       |
// 7 +   +8  + 5
//   |       |
//   +---+---+
//   0   4   1
template <class PhysProp>
Polygon3d XC::QuadBase9N<PhysProp>::getPolygon(bool initialGeometry) const
  {
    const std::deque<Pos3d> positions= this->getPosNodes(initialGeometry);
    //Not collinear nodes:
    Polygon3d retval(positions[7], positions[0], positions[4]);
    //rest of the nodes:
    retval.push_back(positions[1]);
    retval.push_back(positions[5]);
    retval.push_back(positions[2]);
    retval.push_back(positions[6]);
    retval.push_back(positions[3]);
    return retval;
  }

template <class PhysProp>
std::vector<int> XC::QuadBase9N<PhysProp>::getIdxNodes(void) const
  {
    std::vector<int> retval= this->getNodePtrs().getIdx();
    retval.pop_back(); //Node at center removed.
    return retval; 
  }

 
template <class PhysProp>
Segment3d XC::QuadBase9N<PhysProp>::getSide(const size_t &i,bool initialGeometry) const
  {
    Segment3d retval;
    std::cerr << this->getClassName() << "::" << __FUNCTION__
          << "; not implemented." << std::endl;
    return retval;
  }

template <class PhysProp>
XC::Element::NodesEdge XC::QuadBase9N<PhysProp>::getNodesEdge(const size_t &i) const
  {
    Element::NodesEdge retval(3,static_cast<Node *>(nullptr));
    std::cerr << this->getClassName() << "::" << __FUNCTION__
          << "; not implemented." << std::endl;
    return retval;    
  }

template <class PhysProp>
int XC::QuadBase9N<PhysProp>::getEdgeNodes(const Node *n1,const Node *n2) const
  {
    int retval= -1;
    std::cerr << this->getClassName() << "::" << __FUNCTION__
              << "; not implemented." << std::endl;
    return retval;
  }

template <class PhysProp>
ID XC::QuadBase9N<PhysProp>::getLocalIndexNodesEdge(const size_t &i) const
  {
    ID retval(2);
    std::cerr << this->getClassName() << "::" << __FUNCTION__
          << "; not implemented." << std::endl;
    return retval;
  }

template <class PhysProp>
void XC::QuadBase9N<PhysProp>::zeroLoad(void)
  {
    PlaneElement<9,PhysProp>::zeroLoad();
    this->physicalProperties.getMaterialsVector().zeroInitialGeneralizedStrains(); //Removes initial deformations.
    return;
  }
 
template <class PhysProp>
int XC::QuadBase9N<PhysProp>::addLoad(ElementalLoad *theLoad, double loadFactor)
  {
    if(this->isDead())
      std::cerr << this->getClassName() 
                << "; load over inactive elements: "
                << this->getTag() << std::endl;
    else
      {
        if(const BidimStrainLoad *strainLoad= dynamic_cast<const BidimStrainLoad *>(theLoad)) //Prescribed strains.
          {
            static std::vector<Vector> initStrains;
            initStrains= strainLoad->getStrains();
            for(std::vector<Vector>::iterator i= initStrains.begin();i!=initStrains.end();i++)
              (*i)*= loadFactor;
            this->physicalProperties.getMaterialsVector().incrementInitialGeneralizedStrains(initStrains);
          }
        else
          {
            std::cerr << this->getClassName() << "::" << __FUNCTION__
              << "; load type unknown for element with tag: " <<
            this->getTag() << std::endl;
            return -1;
          }
      }
    return 0;
  }

template <class PhysProp>
int XC::QuadBase9N<PhysProp>::getVtkCellType(void) const
  { return VTK_QUADRATIC_QUAD; }

} // end of XC namespace
#endif