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

Go to the documentation of this file.

#ifndef _MAT_ELASTIC_BASSO_H_
#define _MAT_ELASTIC_BASSO_H_

// std includes

// Basso includes
#include "basso.h"
#include "Material.h"

namespace basso {


class MatElastic : public Material  
{

public:
        
        // PUBLIC DEFS
        
        // CONSTRUCTORS 
        MatElastic( Numeric E, Numeric nu, Numeric rho=1.0, Numeric alpha=0.0, Numeric diff=0.0 ) : Material(), C(6,6)
                { 
                        fE=E; fNu=nu; fRho=rho; 
                        fThermalExpansion=alpha; 
                        fDiffusivity=diff;
                        FormC(); 
                }
        
        virtual ~MatElastic( ) { }
        
        // ACCESSORS
        
        // MEMBER FUNCTIONS
        
        virtual Numeric Density(  const MaterialState sv=MaterialState() ) const { return fRho; }
        
        virtual Numeric ThermalExpansion( ) const       { return fThermalExpansion; }
                
        virtual Numeric Diffusivity( ) const    { return fDiffusivity; }
        
        virtual void ComputeStress( nArray &stress, const StressState *ssptr, const nArray &strain  ) 
                {
//                      this->CheckVoigtInput(stress);
//                      mult( C, strain, stress );  
                }
                
        virtual void ComputeStress( nArray &stress, const StressState *ssptr, const MaterialState &sv  ) 
                {
                        // ComputeStress( stress, sv.Strain() ); 
                }

        virtual void ComputeStress( nMatrix &stress, const StressState *ssptr, const MaterialState &sv ) 
                {   
//                      this->CheckInput(stress);
/*                      nArray vstress( 6 );
                        ComputeStress( vstress, sv.Strain() );
                        voigt2tensor( vstress, stress ); 
*/
                }
        
        
        virtual void ComputeTangentStiffness( nMatrix &cmat, const StressState *ssptr, const MaterialState sv=MaterialState() ) const
                { 
                        ssptr->Convert( C, cmat );
//                      this->CheckVoigtInput(cmat);
//                      copy( C, cmat ); 
                }
        
protected:
        void SetProperties( Numeric E, Numeric nu, Numeric rho=1.0 ) { fE=E; fNu=nu; fRho=rho; FormC(); }
        void FormC();
        
protected:
        Numeric fE, fNu, fRho, fDiffusivity, fThermalExpansion;
        nMatrix C;
};

void MatElastic::FormC()  
{
        Numeric Cv=fE*(1.0-fNu)/(1.0+fNu)/(1.0-2*fNu), 
                Cp=Cv*fNu/(1.0-fNu), 
                Cs=fE/(1+fNu)/2.0;
        C(0,0)=Cv; C(0,1)=Cp; C(0,2)=Cp; C(0,3)=0.; C(0,4)=0.; C(0,5)=0.;                               
        C(1,0)=Cp; C(1,1)=Cv; C(1,2)=Cp; C(0,3)=0.; C(0,4)=0.; C(0,5)=0.;
        C(2,0)=Cp; C(2,1)=Cp; C(2,2)=Cv; C(0,3)=0.; C(0,4)=0.; C(0,5)=0.;
        C(3,0)=0.; C(3,1)=0.; C(3,2)=0.; C(3,3)=Cs; C(3,4)=0.; C(3,5)=0.; 
        C(4,0)=0.; C(4,1)=0.; C(4,2)=0.; C(4,3)=0.; C(4,4)=Cs; C(4,5)=0.;
        C(5,0)=0.; C(5,1)=0.; C(5,2)=0.; C(5,3)=0.; C(5,4)=0.; C(5,5)=Cs; 
        
}

} // end namespace


#endif

---