TAPsModelDefBasedOnFEMTet.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsModelDefBasedOnFEMTet.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (03/16/2010)
UPDATE          (10/27/2010)
******************************************************************************/
#ifndef TAPs_MODEL_DEF_BASED_ON_FEM_TET_HPP
#define TAPs_MODEL_DEF_BASED_ON_FEM_TET_HPP

#include "ModelDefBasedOnFEMSupport/TAPsModelDefBasedOnFEM.hpp"
#include "../FEM/TAPsFEMMeshTetrahedraGen.hpp"
#include "ModelDefBasedOnFEMSupport/TAPsModelDefBasedOnFEMTet_CD.hpp"

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

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

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

    virtual bool    InitCD ();

    inline virtual FEM::ModelDefBasedOnFEMTet_CD<T> const & CDUnit () const { return m_CD; }
    inline virtual FEM::ModelDefBasedOnFEMTet_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 ();

#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                
    );

    // 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 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 ();

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

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

    inline virtual void UpdateSurfaceMesh ();

//=============================================================================
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 ModelDefBasedOnFEMTet
//=============================================================================
END_NAMESPACE_TAPs
//-----------------------------------------------------------------------------
#include "TAPsModelDefBasedOnFEMTet.cpp"

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