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

Go to the documentation of this file.

#ifndef _STRESS_STATE_BASSO_H_
#define _STRESS_STATE_BASSO_H_

// std includes

// Basso includes
#include "basso.h"

namespace basso {

        
void voigt2tensor( const nArray &vstress, nMatrix &tstress ) 
{
        tstress(0,0)=vstress[0]; tstress(0,1)=vstress[5]; tstress(0,2)=vstress[4];
        tstress(1,0)=vstress[5]; tstress(1,1)=vstress[1]; tstress(1,2)=vstress[3];
        tstress(2,0)=vstress[4]; tstress(2,1)=vstress[3]; tstress(2,2)=vstress[2];
}

void tensor2voigt( const nMatrix &tstress, nArray &vstress ) 
{
        vstress[0]=tstress(0,0);
        vstress[1]=tstress(1,1);
        vstress[2]=tstress(2,2);
        vstress[3]=tstress(1,2);
        vstress[4]=tstress(2,0);
        vstress[5]=tstress(0,1);
}

class VoigtStress
{

public:
        VoigtStress() : vs(6) { }
        
        Numeric &operator() ( int i ) { return vs[i]; }
        Numeric &operator() ( int i, int j ) 
        {
                if ( i==j )
                        return vs[i];
                else if ( ( i==1 && j==0 ) || ( i==0 && j==1 ) )
                        return vs[5];
                else if ( ( i==2 && j==1 ) || ( i==1 && j==2 ) )
                        return vs[3];
                else  // if ( ( i==2 && j==0 ) || ( i==0 && j==2 ) )
                        return vs[4];
        }
        
        void FillTensor( nMatrix &T ) const
        {
                T(0,0)=vs[0]; T(0,1)=vs[5]; T(0,2)=vs[4]; 
                T(1,0)=vs[5]; T(1,1)=vs[1]; T(1,2)=vs[3];
                T(1,0)=vs[4]; T(2,1)=vs[3]; T(2,2)=vs[2];
        }
        
private:
        nArray vs;
        
};

class VoigtStrain 
{
public:
        VoigtStrain() : vs(6) { }
        
        Numeric operator() ( int i ) { return vs[i]; }
        Numeric operator() ( int i, int j ) 
        {
                if ( i==j )
                        return vs[i];
                else if ( ( i==1 && j==0 ) || ( i==0 && j==1 ) )
                        return 0.5*vs[5];
                else if ( ( i==2 && j==1 ) || ( i==1 && j==2 ) )
                        return 0.5*vs[3];
                else  // if ( ( i==2 && j==0 ) || ( i==0 && j==2 ) )
                        return 0.5*vs[4];
        }
        
        void FillTensor( nMatrix &T ) const
        {
                T(0,0)=vs[0]; T(0,1)=0.5*vs[5]; T(0,2)=0.5*vs[4]; 
                T(1,0)=0.5*vs[5]; T(1,1)=vs[1]; T(1,2)=0.5*vs[3];
                T(1,0)=vs[4]; T(2,1)=0.5*vs[3]; T(2,2)=0.5*vs[2];
        }
        
private:
        nArray vs;
        
};


class StressState
{
public:
        StressState() { }
        
        virtual ~StressState() { }
        
        virtual void Convert( const nMatrix &C, nMatrix &Cmat ) const = 0;
        
        virtual int VoigtDimension() const = 0;
        
        virtual int Dimension() const = 0;
        
};

class PlaneStrainState : public StressState
{

public:
        PlaneStrainState(   ) : StressState() { }
        
        virtual void Convert( const nMatrix &C, nMatrix &Cmat ) const 
        {
#ifdef ALLOW_DYNAMIC_RESIZE
                if ( mat_nrows(Cmat)!=3 || mat_ncols(Cmat)!=3 )
                        resize(Cmat,3,3);
#elif EXPLICIT_BOUNDS_CHECK
                if (  ( mat_nrows(Cmat)!=3 || mat_ncols(Cmat)!=3 ) ||
                        ( mat_nrows(C)!=6 || mat_ncols(C)!=6 ) )
                        ErrorMessage("PlaneStrainState.Convert()","incorrct matrix size");
#endif
                Cmat(0,0)=C(0,0); Cmat(0,1)=C(0,1); Cmat(0,2)=C(0,5);
                Cmat(1,0)=C(1,0); Cmat(1,1)=C(1,1); Cmat(1,2)=C(1,5);
                Cmat(2,0)=C(5,0); Cmat(2,1)=C(5,1); Cmat(2,2)=C(5,5);           
        }
        
        virtual int VoigtDimension() const { return 3; }
        
        virtual int Dimension() const { return 2; }
         
};

class PlaneStressState : public StressState
{

public:
        PlaneStressState( ) : StressState() { }
        
        virtual void Convert( const nMatrix &C, nMatrix &Cmat ) const 
        {
#ifdef ALLOW_DYNAMIC_RESIZE
                if ( mat_nrows(Cmat)!=3 || mat_ncols(Cmat)!=3 )
                        resize(Cmat,3,3);
#elif EXPLICIT_BOUNDS_CHECK
                if (  ( mat_nrows(Cmat)!=3 || mat_ncols(Cmat)!=3 ) ||
                        ( mat_nrows(C)!=6 || mat_ncols(C)!=6 ) )
                        ErrorMessage("PlaneStressState.Convert()","incorrct matrix size");
#endif
                Numeric a=1.0/C(2,2);
                Cmat(0,0)=C(0,0)-C(0,2)*C(2,0)*a; Cmat(0,1)=C(0,1)-C(0,2)*C(2,1)*a; Cmat(0,2)=C(0,5);
                Cmat(1,0)=C(1,0)-C(1,2)*C(2,0)*a; Cmat(1,1)=C(1,1)-C(1,2)*C(2,1)*a; Cmat(1,2)=C(1,5);
                Cmat(2,0)=C(5,0); Cmat(2,1)=C(5,1); Cmat(2,2)=C(5,5);           
        }
        
        virtual int VoigtDimension() const { return 3; }
        
        virtual int Dimension() const { return 2; }
        
private:
         
};

class AxisymmetricState : public StressState
{

public:
        AxisymmetricState(   ) : StressState() {  }
        
        virtual void Convert( const nMatrix &C, nMatrix &Cmat ) const 
        {
#ifdef ALLOW_DYNAMIC_RESIZE
                if ( mat_nrows(Cmat)!=4 || mat_ncols(Cmat)!=4 )
                        resize(Cmat,4,4);
#elif EXPLICIT_BOUNDS_CHECK
                if (  ( mat_nrows(Cmat)!=4 || mat_ncols(Cmat)!=4 ) ||
                        ( mat_nrows(C)!=6 || mat_ncols(C)!=6 ) )
                        ErrorMessage("PlaneStrainState.Convert()","incorrct matrix size");
#endif
                Cmat(0,0)=C(0,0); Cmat(0,1)=C(0,1); Cmat(0,2)=0.0;    Cmat(0,3)=C(0,5);
                Cmat(1,0)=C(1,0); Cmat(1,1)=C(1,1); Cmat(0,2)=0.0;    Cmat(0,3)=C(0,5);
                Cmat(2,0)=0.0;    Cmat(2,1)=0.0;    Cmat(2,2)=C(5,5); Cmat(2,3)=0.0;
                Cmat(3,0)=C(5,0); Cmat(3,1)=C(5,1); Cmat(3,2)=0.0;    Cmat(3,3)=1.0;    
        }
        
        virtual int VoigtDimension() const { return 4; }
        
        virtual int Dimension() const { return 2; }
        
private:

         
};

class StressState3D : public StressState
{

public:
        StressState3D(   ) : StressState() { }

        virtual void Convert( const nMatrix &C, nMatrix &Cmat ) const 
        { 
#ifdef ALLOW_DYNAMIC_RESIZE
                if ( mat_nrows(Cmat)!=6 || mat_ncols(Cmat)!=6 )
                        resize(Cmat,6,6);
#elif EXPLICIT_BOUNDS_CHECK
                if (  ( mat_nrows(Cmat)!=6 || mat_ncols(Cmat)!=6 ) ||
                        ( mat_nrows(C)!=6 || mat_ncols(C)!=6 ) )
                        ErrorMessage("StressState3D.Convert()","incorrct matrix size");
#endif
                copy(C,Cmat);   
        }
        
        virtual int VoigtDimension() const { return 6; }
        
        virtual int Dimension() const { return 3; }
        
private:
         
};

class StressState1D : public StressState
{

public:
        StressState1D(   ) : StressState() { }

        virtual void Convert( const nMatrix &C, nMatrix &Cmat ) const 
        { 
#ifdef ALLOW_DYNAMIC_RESIZE
                if ( mat_nrows(Cmat)!=1 || mat_ncols(Cmat)!=1 )
                        resize(Cmat,1,1);
#elif EXPLICIT_BOUNDS_CHECK
                if (  ( mat_nrows(Cmat)!=1 || mat_ncols(Cmat)!=1 ) ||
                        ( mat_nrows(C)!=1 || mat_ncols(C)!=1 ) )
                        ErrorMessage("StressState1D.Convert()","incorrct matrix size");
#endif
                Cmat(0,0)=C(0,0)+2*C(0,1)*C(0,1)/(C(1,1)+C(1,2));       
        }
        
        virtual int VoigtDimension() const { return 1; }
        
        virtual int Dimension() const { return 1; }
        
private:
         
};

} // end namespace

#endif

---