/Users/jack/Code/basso_dev/inc/ElasticityFormulation.h

Go to the documentation of this file.

#ifndef _ELASTICITY_FORMULATION_BASSO_H_
#define _ELASTICITY_FORMULATION_BASSO_H_

// std includes

// Basso includes
#include "basso.h"
#include "StructuralFormulation.h"
#include "MatElastic.h"
#include "basic_fem_operations.h"

namespace basso {


class ElasticityFormulation : public StructuralFormulation {

public:
        
        // PUBLIC DEFS
        
        // CONSTRUCTORS 
        ElasticityFormulation( list< Element > &elem, const Array< Node > &nodes, 
                const Material *material, Numeric t=1.0, StressStateType ss=kSTRESS_3D ) 
                : StructuralFormulation( elem, nodes, ss )
        { 
                this->fMat=material; 
                Initialize(t);
                CheckElements();
        }
        
        ~ElasticityFormulation() { }
        
        // ACCESSORS
        
        // MEMBER FUNCTIONS
        
protected:
        virtual void SetCMatrix()  
                { this->fMat->ComputeTangentStiffness( this->mCmat, this->fStressState );       }
        
        void CheckElements() const;
        
        void Initialize( Numeric t );
        
private:
        Numeric mThickness;

};
        
void ElasticityFormulation::Initialize( Numeric t ) 
{ 
        this->mSdim = this->fStressState->Dimension();
        mThickness=t; 
        resize( this->activeLocalDofs, this->mSdim );
        this->activeLocalDofs[0]=kDISPx; 
        if ( this->mSdim>1 ) this->activeLocalDofs[1]=kDISPy;
        if ( this->mSdim>2 ) this->activeLocalDofs[2]=kDISPz;
}

void ElasticityFormulation::CheckElements() const
{
        list< Element >::const_iterator eItr;
        for ( eItr=this->fElement->begin(); eItr!=this->fElement->end(); ++eItr ) 
                if ( eItr->Dimension()!=this->mSdim )
                        WarningMessage("ElasticityFormulation",
                                "Elements must be of proper dimension; other elements will be ignored");
}


} // end namespace

#endif

---