TAPsModelDefBasedOnFEM.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsModelDefBasedOnFEM.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (03/31/2010)
UPDATE          (11/29/2010)
******************************************************************************/
#ifndef TAPs_MODEL_DEF_BASED_ON_FEM_HPP
#define TAPs_MODEL_DEF_BASED_ON_FEM_HPP

#ifndef TAPs_ADD_DEFORMABLE_MODEL_BASED_ON_FEM
    #define TAPs_ADD_DEFORMABLE_MODEL_BASED_ON_FEM
#endif

//#define   TAPs_ADD_INTERACTIONS
#ifdef  TAPs_ADD_INTERACTIONS
    #ifndef TAPs_ADVANCED_SIMULATION
        //#define   TAPs_ADVANCED_SIMULATION
    #endif//TAPs_ADVANCED_SIMULATION
    //#include "../../Support/TAPsAdvSimCtrl.hpp"
#endif//TAPs_ADD_INTERACTIONS

#include "../../OpenGLModelFromFile/TAPsReadModels.hpp"
#include "../../FEM/TAPsFEMElement.hpp"
#include "../../FEM/TAPsFEMMesh.hpp"

// For collision detection
#include "TAPsModelDefBasedOnFEM_CD.hpp"

BEGIN_NAMESPACE_TAPs
//=============================================================================
template <typename T, typename DATA>
class ModelDefBasedOnFEM {
//=============================================================================
public:
    // Member Functions -------------------------------------------------------
    friend std::ostream & operator<< ( std::ostream &output, ModelDefBasedOnFEM<T,DATA> const &obj )
    {
        output << obj.StrInfo();
        return output;
    }
    //-------------------------------------------------------------------------
    ModelDefBasedOnFEM ();
    ModelDefBasedOnFEM ( std::string const & modelFile );
    //ModelDefBasedOnFEM ( ModelDefBasedOnFEM<T,DATA> const &orig );
    virtual ~ModelDefBasedOnFEM ();
    //-------------------------------------------------------------------------
    virtual std::string StrInfo () const;
    //-------------------------------------------------------------------------
    //inline ModelDefBasedOnFEM<T,DATA> & operator= ( ModelDefBasedOnFEM<T,DATA> const &orig );
    //-------------------------------------------------------------------------

    OpenGL::OpenGLHalfEdgeTrigonalModel<T> const &  RefToSurfaceMesh () const   { return *m_pHEModel; }
    OpenGL::OpenGLHalfEdgeTrigonalModel<T> &        RefToSurfaceMesh ()         { return *m_pHEModel; }

    virtual FEM::Mesh<T> *  PtrToFEMMesh ()         = 0;

    BVHTree<T> const *  PtrToBVHTreeOfSurfaceMesh ()    const   { return RefToSurfaceMesh().GetBVHTree(); }
    BVHTree<T> *        PtrToBVHTreeOfSurfaceMesh ()            { return RefToSurfaceMesh().GetBVHTree(); }

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

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

    void ClearExtraVertex ( unsigned int i );

    void ClearAllExtraVertices ();

    virtual bool    InitCD ();

    inline virtual FEM::ModelDefBasedOnFEM_CD<T> const &    CDUnit () const = 0;
    inline virtual FEM::ModelDefBasedOnFEM_CD<T>    &       CDUnit ()       = 0;

    inline virtual void AddConnectionForceInLocalCoordinates ( PointForce<T> * const ptrConnectionForce ) = 0;
    inline virtual void RemoveConnectionForce ( PointForce<T> * const ptrConnectionForce ) = 0;
    inline virtual void ClearAllConnectionLoads () = 0;

    inline virtual void UpdateSurfaceMeshNormals () = 0;

    virtual void AdvanceSimulationStatic () = 0;

    virtual void AdvanceSimulationStatic_woFixedDisplacementNodes () = 0;

#ifdef  TAPs_SIM_DYNAMIC_SYSTEM

    virtual void AdvanceSimulationDynamic () = 0;

    virtual void AdvanceSimulationDynamic_woFixedDisplacementNodes () = 0;
#endif//TAPs_SIM_DYNAMIC_SYSTEM

    inline virtual T GetTimeStep () const = 0;
    inline virtual void SetTimeStep ( T timestep ) = 0;

    inline virtual TransformationSupport<T> const & RefToTransformationSupport () const = 0;
    inline virtual TransformationSupport<T> &       RefToTransformationSupport ()       = 0;

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

    // Data Members -----------------------------------------------------------
    OpenGL::OpenGLHalfEdgeTrigonalModel<T> *    m_pHEModel; 

    // Extra (HE-based) vertices (not nodes) are designed to be inside the FEM mesh.
    // Each extra vertices must be accompanied by a BaseBarycentricCoordsForVertexRef<T,DATA>.
    std::vector< HEVertex<T> * >    m_SetOfExtraVertices;   
    std::vector< BaseBarycentricCoordsForVertexRef<T> * > m_SetOfExtraVertexCoords; 

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

    inline virtual void UpdateSurfaceMesh ();

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

//=============================================================================
#if defined(__gl_h_) || defined(__GL_H__)
//-----------------------------------------------------------------------------
public:
    virtual void Draw ( OpenGL::Enum::DrawMode DM = OpenGL::Enum::POLYGON );

    virtual void DrawUndeformedFEMMesh () = 0;

    virtual void DrawDeformedFEMMesh () = 0;

    virtual void DrawCollidedMeshParts () = 0;

    virtual void DrawAllExtraVertices ();

protected:
    TAPs::OpenGL::OpenGLUsefulObj<T>    m_HelpDrawingUsefulObj;

//-----------------------------------------------------------------------------
#endif
//=============================================================================

//=============================================================================
}; // END CLASS ModelDefBasedOnFEM
//=============================================================================
END_NAMESPACE_TAPs
//-----------------------------------------------------------------------------
#include "TAPsModelDefBasedOnFEM.cpp"

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