An Object Finite Element Library

Elas3DH8.h

Go to the documentation of this file.

/*==============================================================================

                                    O  F  E  L  I

                           Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2008 Rachid Touzani

   This program 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; Version 2 of the License.

   This program 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 this program; if not, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================

                         Definition of class Elas3DH8
             for Three-Dimensional Linearized Elasticity equations using
                       8-node hexahedral finite element

  ==============================================================================*/


#ifndef __ELAS3DH8_H
#define __ELAS3DH8_H


#include "Equa_Solid.h"
#include "Hexa8.h"
#include "Quad4.h"

namespace OFELI {

class Elas3DH8 : virtual public Equa_Solid<double,8,24,4,12>
{

public :

    Elas3DH8()
    {
       _hexa = NULL;
       _quad = NULL;
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class Elas3DH8 is constructed.\n";
       std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    Elas3DH8(const Element *el);

    Elas3DH8(const Side *sd);

    Elas3DH8(const Element *element, const Vect<double>  &u, const double &time=0.);

    Elas3DH8(const Side *side, const Vect<double> &u, const double &time=0.);

    ~Elas3DH8();

    void LMassToLHS(double coef=1.);

    void LMassToRHS(double coef=1.);

    void LMass(double coef) { LMassToLHS(coef); LMassToRHS(coef); }

    void Mass(double coef=1.) { coef=1; cout << "Sorry, consistent mass matrix is not implemented !\n"; }

    void Deviator(double coef=1.);

    void DeviatorToRHS(double coef=1.);

    void Dilatation(double coef=1.);

    void DilatationToRHS(double coef=1.);

    void BodyRHS(UserData<double> &ud);

    void BodyRHS(const Vect<double> &bf, int opt=LOCAL_ARRAY);

private :

   Hexa8                           *_hexa;
   Quad4                           *_quad;
   LocalVect<double,8>             _w;
   double                          _xg[2], _wg[2];
   Point<double>                   _cg, _xl[8];
   LocalMatrix<Point<double> ,8,8> _dsh;
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani