BrickBase.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/>.
//----------------------------------------------------------------------------
//BrickBase.h
                                                                        
#ifndef BrickBase_h
#define BrickBase_h

#include <domain/mesh/element/ElemWithMaterial.h>
#include "domain/mesh/element/utils/physical_properties/NDMaterialPhysicalProperties.h"

namespace XC {
class NDMaterial;
class ParticlePos3d;
class BrickRawLoad;
  
//              hexahedron
//
//                        s              
//                       /               
//            t         /                
//            | 5---------4             
//            |/.        /|             
//            / .       / |             
//           /  .      /  |             
//          6---------7   |             
//          |   .     |   |             
//          |   1.....|...0            
//          |  .      |  /              
//          | .       | /--- r                
//          |.        |/                
//          2---------3            
//
// Inverse mapping code taken from the file trilinearintepolator.c
// created by Dirk-Philip van Herwaarden on 4/21/17.

 
//
// 
class BrickBase: public ElemWithMaterial<8,NDMaterialPhysicalProperties>
  {
  public:
    static const int numberNodes= 8; 
    static const int ndm= 3; 
    static const int ndf= 3; 
    static const double mNodesR[]; 
    static const double mNodesS[]; 
    static const double mNodesT[]; 
  private:
    static Matrix &compute_extrapolation_matrix(void);
  protected:
    //local nodal coordinates, three coordinates for each of eight nodes
    static double xl[ndm][numberNodes];
    ElemPtrArray3d put_on_mesh(const NodePtrArray3d &,meshing_dir dm) const;
    void computeBasis(void) const; //compute coordinate system
  public:
    BrickBase(int classTag);
    BrickBase(int tag,int classTag,const NDMaterialPhysicalProperties &);
    BrickBase(int tag, int classTag,int nd1, int nd2, int nd3, int nd4,int nd5,int nd6,int nd7,int nd8, const NDMaterialPhysicalProperties &);
    size_t getDimension(void) const;
    double getVolume(bool initialGeometry= true) const;
    BoolArray3d getNodePattern(void) const;
    int getVtkCellType(void) const;

    const BrickRawLoad *vector3dRawLoadLocal(const std::vector<Vector> &);
    const BrickRawLoad *vector3dRawLoadGlobal(const std::vector<Vector> &);
    void zeroLoad(void);    
    int addLoad(ElementalLoad *theLoad, double loadFactor);
    
    Matrix getLocalAxes(bool initialGeometry= true) const;
    Pos3d getGlobalCoordinates(const double &r, const double &s, const double &t) const;  
    ParticlePos3d getNaturalCoordinates(const Pos3d &, bool initialGeometry= true) const;
    const Matrix &getExtrapolationMatrix(void) const;
  };

} // end of XC namespace
#endif