TAPsModelDefBasedOnFEMHex.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsModelDefBasedOnFEMHex.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (02/26/2010)
UPDATE          (10/27/2010)
******************************************************************************/
#ifndef TAPs_MODEL_DEF_BASED_ON_FEM_HEX_HPP
#define TAPs_MODEL_DEF_BASED_ON_FEM_HEX_HPP

#include "ModelDefBasedOnFEMSupport/TAPsModelDefBasedOnFEM.hpp"
#include "../FEM/TAPsFEMMeshHexahedraGen.hpp"
#include "ModelDefBasedOnFEMSupport/TAPsModelDefBasedOnFEMHex_CD.hpp"
#include "../Core/TAPsException.hpp"

BEGIN_NAMESPACE_TAPs
//=============================================================================
template <typename T, typename DATA>
class ModelDefBasedOnFEMHex : public ModelDefBasedOnFEM<T,DATA> {
//=============================================================================
public:
    // Member Functions -------------------------------------------------------
    friend std::ostream & operator<< ( std::ostream &output, ModelDefBasedOnFEMHex<T,DATA> const &obj )
    {
        output << obj.StrInfo();
        return output;
    }
    //-------------------------------------------------------------------------
    ModelDefBasedOnFEMHex ( std::string const & modelFile );
    //ModelDefBasedOnFEMHex ( ModelDefBasedOnFEMHex<T,DATA> const &orig );
    virtual ~ModelDefBasedOnFEMHex ();
    //-------------------------------------------------------------------------
    virtual std::string StrInfo () const;
    //-------------------------------------------------------------------------
    bool GenerateGridsWithDimension (
        T iGridWidth,   
        T iGridHeight,  
        T iGridDepth,   
        T   mass = 1,                   
        T   YoungsModulus = 96,         
        T   PoissonsRatio = T(1)/T(3),  
        bool bGenElementsInsideVolume = true,       
        bool bRecordsElementsGridLocations = true   
    );
    bool GenerateGridsWithSize (
        int iNumOfXSlices,  
        int iNumOfYSlices,  
        int iNumOfZSlices,  
        T   mass = 1,                   
        T   YoungsModulus = 96,         
        T   PoissonsRatio = T(1)/T(3),  
        bool bGenElementsInsideVolume = true,       
        bool bRecordsElementsGridLocations = true   
    );
    //-------------------------------------------------------------------------
    //inline ModelDefBasedOnFEMHex<T,DATA> & operator= ( ModelDefBasedOnFEMHex<T,DATA> const &orig );
    //-------------------------------------------------------------------------

    unsigned int GetNumberOfElementsPerGridSize_X ()    { return m_NumOfElementsPerGridSize[0]; }
    unsigned int GetNumberOfElementsPerGridSize_Y ()    { return m_NumOfElementsPerGridSize[1]; }
    unsigned int GetNumberOfElementsPerGridSize_Z ()    { return m_NumOfElementsPerGridSize[2]; }

    FEM::MeshHexahedra<T> const &   RefToFEMMesh () const   { return m_FEMMesh; }
    FEM::MeshHexahedra<T> &         RefToFEMMesh ()         { return m_FEMMesh; }

    virtual FEM::Mesh<T> *  PtrToFEMMesh ()         { return &m_FEMMesh; }

    virtual bool    InitCD ();

    inline virtual FEM::ModelDefBasedOnFEMHex_CD<T> const & CDUnit () const { return m_CD; }
    inline virtual FEM::ModelDefBasedOnFEMHex_CD<T> &       CDUnit ()       { return m_CD; }

    inline virtual void AddConnectionForceInLocalCoordinates ( PointForce<T> * const ptrConnectionForce ) { m_FEMMesh.AddConnectionForceInLocalCoordinates( ptrConnectionForce ); }
    inline virtual void RemoveConnectionForce ( PointForce<T> * const ptrConnectionForce )  { m_FEMMesh.RemoveConnectionForce( ptrConnectionForce ); }
    inline virtual void ClearAllConnectionLoads ()  { m_FEMMesh.ClearAllConnectionLoads(); }

    inline virtual void UpdateSurfaceMeshNormals () { m_pHEModel->CalAndSetNormals(); }

    virtual void AdvanceSimulationStatic ();

    virtual void AdvanceSimulationStatic_woFixedDisplacementNodes ();

    void AllocateElementsGridLocations ( unsigned int xWidth, unsigned int yHeight, unsigned int zDepth ) throw(...);
    void DeallocateElementsGridLocations ();
    inline int GetElementIDFromGridLocation_wErrChk ( unsigned int xWidth, unsigned int yHeight, unsigned int zDepth ) const;
    inline int GetElementIDFromGridLocation ( unsigned int xWidth, unsigned int yHeight, unsigned int zDepth ) const;
    void ListAllElementIDsInGridLocationFormat_zYx () const;

#ifdef  TAPs_SIM_DYNAMIC_SYSTEM

    virtual void AdvanceSimulationDynamic ();

    virtual void AdvanceSimulationDynamic_woFixedDisplacementNodes ();
#endif//TAPs_SIM_DYNAMIC_SYSTEM

    inline virtual T GetTimeStep () const   { return m_FEMMesh.GetTimeStep(); }
    inline virtual void SetTimeStep ( T timestep )  { return m_FEMMesh.SetTimeStep( timestep ); }

    inline virtual TransformationSupport<T> const & RefToTransformationSupport () const { return m_FEMMesh.RefToTransformationSupport(); }
    inline virtual TransformationSupport<T> &       RefToTransformationSupport ()       { return m_FEMMesh.RefToTransformationSupport(); }

    virtual HEVertex<T> * AddExtraVertexInLocalSpace (
        Vector3<T> location,                    
        int elementID,                          
        FEM::Element<T> * elementPtr,           
        Vector3<unsigned int> gridLocations,    
        Vector3<T> coordinates                  
    );

    virtual bool FindElementThatContainsPointInLocalSpace (
        Vector3<T> const & point,               
        unsigned int & elementID,               
        Vector3<unsigned int> & gridLocations,  
        Vector3<T> & coordinates                
    );

    //-------------------------------------------------------------------------
    // FIX/UNFIX

    void FixOrUnfixNodesByRowColumnDepthRange (
        unsigned int x0,    
        unsigned int x1,    
        unsigned int y0,    
        unsigned int y1,    
        unsigned int z0,    
        unsigned int z1,    
        bool    doFix       
    );

    void FixNodesByRowColumnDepthRange (
        unsigned int x0,    
        unsigned int x1,    
        unsigned int y0,    
        unsigned int y1,    
        unsigned int z0,    
        unsigned int z1     
    )   { FixOrUnfixNodesByRowColumnDepthRange( x0, x1, y0, y1, z0, z1, true ); }

    void UnfixNodesByRowColumnDepthRange (
        unsigned int x0,    
        unsigned int x1,    
        unsigned int y0,    
        unsigned int y1,    
        unsigned int z0,    
        unsigned int z1     
    )   { FixOrUnfixNodesByRowColumnDepthRange( x0, x1, y0, y1, z0, z1, false ); }

    //-------------------------------------------------------------------------
    // Data Members -----------------------------------------------------------
//=============================================================================
protected:
    // Member Functions -------------------------------------------------------

    bool GenerateMeshFromGridGenerator (
        GridGenerator<T,DATA> & GridGen,    
        T   mass = 1,                       
        T   YoungsModulus = 96,             
        T   PoissonsRatio = T(1)/T(3),      
        bool bGenElementsInsideVolume = true,       
        bool bRecordsElementsGridLocations = true   
    );

    //bool GenerateBaryCoordsFromGridGenerator ( GridGenerator<T,DATA> & GridGen );

    //void RecordBarycentric ();
    //void RecordBarycentric ( OpenGL::XPolygonalModel<T> * pXPModel );
    //void RecordBarycentric ( OpenGL::HalfEdgeModel<T> * pHEModel );

    //-------------------------------------------------------------------------

    inline void LocateGrid ( T x, T y, T z, unsigned int &i, unsigned int &j, unsigned int &k );

    void SetHEVertexIDToElementLocation ();

    //-------------------------------------------------------------------------
    // Simulation

    inline virtual void UpdateSurfaceMesh ();

    // Data Members -----------------------------------------------------------
    FEM::MeshHexahedra<T>   m_FEMMesh;  
    std::vector< HexBarycentricCoordsForVertexRef<T,DATA> >     m_BaryCoords;   
    //std::vector< unsigned int >   m_ElementID;    //!< identify the element that the barycentric coordinates are embedded in
    GridGenerator<T,DATA>   m_GridGen;  
    FEM::ModelDefBasedOnFEMHex_CD<T>    m_CD;   

    // The element ID is ranging from 0 to the number of elements minus one.
    unsigned int m_NumOfElementsPerGridSize[3]; 
    int *** m_ElementsGridLocs; 

//=============================================================================
private:
    // Member Functions -------------------------------------------------------
    // Data Members -----------------------------------------------------------
//=============================================================================

#if defined(__gl_h_) || defined(__GL_H__)
public:
    virtual void DrawUndeformedFEMMesh ();

    virtual void DrawDeformedFEMMesh ();

    virtual void DrawCollidedMeshParts ();
#endif

//=============================================================================
}; // END CLASS ModelDefBasedOnFEMHex
//=============================================================================
END_NAMESPACE_TAPs
//-----------------------------------------------------------------------------
#include "TAPsModelDefBasedOnFEMHex.cpp"

// Include definition if TAPs_USE_EXPORT is not defined
//#if !defined( TAPs_USE_EXPORT )
//  #include "TAPsModelDefBasedOnFEMHex.cpp"
//#endif
//-----------------------------------------------------------------------------
#endif
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----