TAPsFEMTetrahedron.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsFEMTetrahedron.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (12/15/2009)
UPDATE          (05/01/2010)
******************************************************************************/
#ifndef TAPs_FEM_TETRAHEDRON_HPP
#define TAPs_FEM_TETRAHEDRON_HPP

#include "TAPsFEMElement.hpp"

//  Node Numbering Convention
//
//             3
//           / |\
//         /   | \
//      /      |  \
//   /        |    \
//  2---------|-----1
//   \        |    /
//    \      |    /
//     \     |   /
//      \    |  /
//       \  |  /
//        \ | /
//         \|/
//          0

BEGIN_NAMESPACE_TAPs__FEM
//=============================================================================
// Forward declaration
template <typename T> class MeshTetrahedra;

template <typename T>
class Tetrahedron : public Element<T> {
//=============================================================================
    friend  class MeshTetrahedra<T>;
public:
    // Member Functions -------------------------------------------------------
    friend std::ostream & operator<< ( std::ostream &output, Tetrahedron<T> const &obj )
    {
        output << obj.StrInfo();
        return output;
    }
    //-------------------------------------------------------------------------
    Tetrahedron ( 
        std::vector< Vector3<T> > * deformedNodes,      
        std::vector< Vector3<T> > * undeformedNodes,    
        T YoungModulus,     
        T PoissionRatio,    
        unsigned int globalNode0,   
        unsigned int globalNode1,   
        unsigned int globalNode2,   
        unsigned int globalNode3    
    );
    //Tetrahedron ( Tetrahedron<T> const &orig );
    virtual ~Tetrahedron ();
    //-------------------------------------------------------------------------
    virtual std::string StrInfo () const;
    //-------------------------------------------------------------------------
    //Tetrahedron<T> & operator= ( Tetrahedron<T> const &orig );

    //-------------------------------------------------------------------------
    // Get/Set Functions

    inline Matrix3x3<T> const & GetStiffnessMatrix ( unsigned int n, unsigned int m ) const;

    //-------------------------------------------------------------------------
    // Operations

    virtual inline Vector3<T> * const DeformedNode ( unsigned int i );

    virtual inline Vector3<T> * const UndeformedNode ( unsigned int i );

    inline Vector3<T> ComputeDeformedPositionBasedOnCoordinates ( T coord0, T coord1, T coord2, T coord3 );

    inline Vector3<T> ComputeUndeformedPositionBasedOnCoordinates ( T coord0, T coord1, T coord2, T coord3 );

    //---------------------------------------------------------------
    // Inherited from FEMElement class
    //---------------------------------------------------------------
    inline void SetNode ( unsigned int i, unsigned int globalNodeNumber );

    //inline virtual void ResetToDeformedState ();

    inline virtual void ResetToUndeformedState ();

    inline virtual T    CalDeformedVolume () const;

    inline virtual T    CalUndeformedVolume ();

    virtual void UpdateStrainDisplacementMatrixElements ();

    virtual void UpdateStiffnessMatrixElements ();

    inline unsigned int GetGlobalNodeNumberFromLocalNode ( unsigned int i ) const
    {
        assert( 0 <= i && i <= 3 );
        return ids[i];
    }
    //---------------------------------------------------------------
    // Calculations

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

    // Data Members -----------------------------------------------------------

    unsigned int    ids[4]; 

    // Element's properties
    T   b[4];   
    T   c[4];   
    T   d[4];   

    // Stiffness matrix
    // The tetrahedron stiffness matrix is a 12x12 matrix which can be broken into 4x4 of 3x3 matrices
    Matrix3x3<T>    Ke[4][4];   

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

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

    virtual void DrawUndeformedMesh ();

    virtual void DrawDeformedMesh ();
#endif

//=============================================================================
}; // END CLASS Tetrahedron
//=============================================================================
END_NAMESPACE_TAPs__FEM
//-----------------------------------------------------------------------------
#include "TAPsFEMTetrahedron.cpp"

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