TAPsFEMElement.cpp

Go to the documentation of this file.

/******************************************************************************
TAPsFEMElement.cpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (12/15/2009)
UPDATE          (03/19/2010)
******************************************************************************/
#include "TAPsFEMElement.hpp"
// Using Inclusion Model (i.e. definitions are included in declarations)
//                       (this name.cpp is included in name.hpp)
// Each friend is defined directly inside its declaration.

BEGIN_NAMESPACE_TAPs__FEM
//=============================================================================
// Constructors
//-----------------------------------------------------------------------------
template <typename T>
Element<T>::Element ( 
        std::vector< Vector3<T> > * deformedNodes,      
        std::vector< Vector3<T> > * undeformedNodes,    
        T _YoungModulus,    
        T _PoissonRatio     
) :
    pU( deformedNodes ),
    pX( undeformedNodes ),
    YoungModulus(_YoungModulus),
    PoissonRatio(_PoissonRatio),
    Volume(0),
    Ee(0), Ef(0), Eg(0)
{
}
//-----------------------------------------------------------------------------
//template <typename T>
//Element<T>::Element ( Element<T> const &orig )
//{}
//-----------------------------------------------------------------------------
template <typename T>
Element<T>::~Element ()
{}
//-----------------------------------------------------------------------------
//template <typename T>
//std::string Element<T>::StrInfo () const
//{}
//-----------------------------------------------------------------------------
//=============================================================================
// Assignment Operator
//-----------------------------------------------------------------------------
//template <typename T>
//Element<T> & Element<T>::operator= ( Element<T> const &orig )
//{}
//-----------------------------------------------------------------------------


//-----------------------------------------------------------------------------
// Set Young's modulus
template <typename T>
void Element<T>::SetYoungModulus ( T t )
{
    YoungModulus = t;
    UpdateElasticMatrixElements();
}

//-----------------------------------------------------------------------------
// Set Poisson's Ratio
template <typename T>
void Element<T>::SetPoissonRatio ( T t )
{
    PoissonRatio = t;
    UpdateElasticMatrixElements();
}

//-----------------------------------------------------------------------------
template <typename T>
void Element<T>::UpdateElasticMatrixElements ()
{
    //     | e f f 0 0 0 |
    //     | f e f 0 0 0 |
    // E = | f f e 0 0 0 |
    //     | 0 0 0 g 0 0 |
    //     | 0 0 0 0 g 0 |
    //     | 0 0 0 0 0 g |

    T VE = Volume*YoungModulus;
    T scale = VE / ((1+PoissonRatio)*(1-2*PoissonRatio));
    Ee = (1-PoissonRatio) * scale;
    Ef = PoissonRatio * scale;
    Eg = (0.5-PoissonRatio) * scale;

    // DEBUG
    //std::cout << "Volume, Young Modulus, Poisson Ratio: " << Volume << ", " << YoungModulus << ", " << PoissonRatio << "\n";
    //std::cout << "Ee, Ef, Eg: " << Ee << ", " << Ef << ", " << Eg << "\n";
}

//=============================================================================
END_NAMESPACE_TAPs__FEM
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----