tensors.h

/*************************************************************************************
 * MechSys - A C++ library to simulate (Continuum) Mechanical Systems                *
 * Copyright (C) 2005 Dorival de Moraes Pedroso <dorival.pedroso at gmail.com>       *
 * Copyright (C) 2005 Raul Dario Durand Farfan  <raul.durand at gmail.com>           *
 *                                                                                   *
 * This file is part of MechSys.                                                     *
 *                                                                                   *
 * MechSys is free software; you can redistribute it and/or modify it under the      *
 * terms of the GNU General Public License as published by the Free Software         *
 * Foundation; either version 2 of the License, or (at your option) any later        *
 * version.                                                                          *
 *                                                                                   *
 * MechSys is distributed in the hope that it will be useful, but WITHOUT ANY        *
 * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A   *
 * PARTICULAR PURPOSE. See the GNU General Public License for more details.          *
 *                                                                                   *
 * You should have received a copy of the GNU General Public License along with      *
 * MechSys; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, *
 * Fifth Floor, Boston, MA 02110-1301, USA                                           *
 *************************************************************************************/

#ifndef MECHSYS_TENSORS_H
#define MECHSYS_TENSORS_H

#ifdef HAVE_CONFIG_H
  #include "config.h"
#else
  #ifndef REAL
    #define REAL double
  #endif
#endif

#include <blitz/tinyvec-et.h>
#include <blitz/tinymat.h>

namespace Tensors
{



typedef blitz::TinyVector<REAL,6> Tensor2;


typedef blitz::TinyMatrix<REAL,6,6> Tensor4;

// Three-dimensional vector == Tensor1
typedef blitz::TinyVector<REAL,3> Tensor1;
 // end of group Tensor2nd4th

Tensor2 I;

Tensor4 IIsym;

Tensor4 IdyI;

Tensor4 Psd;

Tensor4 Piso;




// Initialize (global) second order identity tensor
int __initialize_the_I(Tensor2 & theI) // {{{
{
    theI = 1.0, 1.0, 1.0, 0.0, 0.0, 0.0;
    return 0;
} // }}}

// Initialize (global) fourth order identity tensor
int __initialize_the_IIsym(Tensor4 & theIIsym) // {{{
{
    theIIsym = 1.0, 0.0, 0.0, 0.0, 0.0, 0.0,
               0.0, 1.0, 0.0, 0.0, 0.0, 0.0,
               0.0, 0.0, 1.0, 0.0, 0.0, 0.0,
               0.0, 0.0, 0.0, 1.0, 0.0, 0.0,
               0.0, 0.0, 0.0, 0.0, 1.0, 0.0,
               0.0, 0.0, 0.0, 0.0, 0.0, 1.0;
    return 0;
} // }}}

// Initialize (global) IdyI fourth order tensor
int __initialize_the_IdyI(Tensor4 & theIdyI) // {{{
{
    theIdyI = 1.0, 1.0, 1.0, 0.0, 0.0, 0.0,
              1.0, 1.0, 1.0, 0.0, 0.0, 0.0,
              1.0, 1.0, 1.0, 0.0, 0.0, 0.0,
              0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
              0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
              0.0, 0.0, 0.0, 0.0, 0.0, 0.0;
    return 0;
} // }}}

// Initialize (global) fourth order symmetric-deviatoric tensor
int __initialize_the_Psd(Tensor4 & thePsd) // {{{
{
    thePsd =  2.0/3.0, -1.0/3.0, -1.0/3.0, 0.0, 0.0, 0.0,
             -1.0/3.0,  2.0/3.0, -1.0/3.0, 0.0, 0.0, 0.0,
             -1.0/3.0, -1.0/3.0,  2.0/3.0, 0.0, 0.0, 0.0,
                  0.0,      0.0,      0.0, 1.0, 0.0, 0.0,
                  0.0,      0.0,      0.0, 0.0, 1.0, 0.0,
                  0.0,      0.0,      0.0, 0.0, 0.0, 1.0;
    return 0;
} // }}}

// Initialize (global) fourth order isotropic tensor
int __initialize_the_Piso(Tensor4 & thePiso) // {{{
{
    thePiso = 1.0/3.0, 1.0/3.0, 1.0/3.0, 0.0, 0.0, 0.0,
              1.0/3.0, 1.0/3.0, 1.0/3.0, 0.0, 0.0, 0.0,
              1.0/3.0, 1.0/3.0, 1.0/3.0, 0.0, 0.0, 0.0,
                  0.0,     0.0,     0.0, 0.0, 0.0, 0.0,
                  0.0,     0.0,     0.0, 0.0, 0.0, 0.0,
                  0.0,     0.0,     0.0, 0.0, 0.0, 0.0;
    return 0;
} // }}}

int __dummy1=__initialize_the_I(I);
int __dummy2=__initialize_the_IIsym(IIsym);
int __dummy3=__initialize_the_IdyI(IdyI);
int __dummy4=__initialize_the_Psd(Psd);
int __dummy5=__initialize_the_Piso(Piso);

}; // namespace Tensors

#endif // MECHSYS_TENSORS_H

// vim:fdm=marker

---