An Object Finite Element Library

Elas2DT3.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 Elas2DT3
             for Two-Dimensional Linearized Elasticity equations using
                           3-node triangular finite element

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


#ifndef __ELAS2DT3_H
#define __ELAS2DT3_H


#include "Equa_Solid.h"
#include "Triang3.h"
#include "Line2.h"

namespace OFELI {

class Elas2DT3 : public Equa_Solid<double,3,6,2,4>
{

 public :

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

    Elas2DT3(const Element *el)
    {
       ElementSet(el);
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class Elas2DT3 is constructed.\n";
       std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    Elas2DT3(const Side *sd)
    {
       SideSet(sd);
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class Elas2DT3 is constructed.\n";
       std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

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

    Elas2DT3(const Element *el, const Vect<double> &u, double time, double deltat, int scheme);

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

    Elas2DT3(const Side *sd, const Vect<double> &u, double time, double deltat, int scheme);

#ifndef DOXYGEN_SHOULD_SKIP_THIS
/*    Elas2DT3(Mesh &mesh, SkSMatrix<double> &A, Vect<double> &b) : Equation<double,3,6,2,4>(mesh, A, b) { }
    Elas2DT3(Mesh &mesh, SkMatrix<double>  &A, Vect<double> &b) : Equation<double,3,6,2,4>(mesh, A, b) { }
    Elas2DT3(Mesh &mesh, SpMatrix<double>  &A, Vect<double> &b) : Equation<double,3,6,2,4>(mesh, A, b) { }*/

// Constructors for a mesh, matrix, right-hand side and previous solution
/*   Elas2DT3(Mesh &mesh, SkSMatrix<double> &A, Vect<double> &b, Vect<double> &u) :
            Equation<double,3,6,2,4>(mesh, A, b, u) { }
            Elas2DT3(Mesh &mesh, SkMatrix<double>  &A, Vect<double> &b, Vect<double> &u) :
            Equation<double,3,6,2,4>(mesh, A, b, u) { }
   Elas2DT3(Mesh &mesh, SpMatrix<double>  &A, Vect<double> &b, Vect<double> &u) :
            Equation<double,3,6,2,4>(mesh, A, b, u) { }*/
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

    ~Elas2DT3();

    void Media(double E, double nu, double rho);

    void PlaneStrain();

    void PlaneStrain(double E, double nu);

    void PlaneStress();

    void PlaneStress(double E, double nu);

    void ElementSet(const Element *el);

    void SideSet(const Side *sd);

#ifndef DOXYGEN_SHOULD_SKIP_THIS
    void Build();
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

    void LMassToLHS(double coef=1.);

    void LMassToRHS(double coef=1.);

    void LMass(double coef=1.) { LMassToLHS(coef); LMassToRHS(coef); };

    void MassToLHS(double coef=1.);

    void MassToRHS(double coef=1.);

    void Mass(double coef=1.) { MassToLHS(coef); MassToRHS(coef); };

    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 NodeVect<double> &f, int opt=LOCAL_ARRAY);

    void BoundaryRHS(UserData<double> &ud);

    void BoundaryRHS(SideVect<double> &f);

    int SignoriniContact(UserData<double> &ud, double coef=1.e07);

    int SignoriniContact(SideVect<double> &f, double coef=1.e07);

    void Reaction(SideVect<double> &r);

    void ContactPressure(const SideVect<double> &f, double penal, Point<double> &p);

    void Strain(Vect<double> &eps);

    void Stress(Vect<double> &s, double &vm);

    void Periodic(double coef=1.e20);

 private :

   double   _E1, _E2, _E3, _E6;
   Triang3  *_tr;
   Line2    *_ln;
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani