An Object Finite Element Library

Beam3DL2.h

Go to the documentation of this file.

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

                                    O  F  E  L  I

                           Object  Finite  Element  Library

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

   Copyright (C) 1998 - 2004 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 Beam3DL2
                       for Beam element with 6 d.o.f. per node
                           using Mindlin - Reissner Theory

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


#ifndef __BEAM3DL2_H
#define __BEAM3DL2_H


#include "Equa_Solid.h"
#include "UserData.h"

namespace OFELI {

class Beam3DL2 : virtual public Equa_Solid<double,2,12,1,1> {

public :

    Beam3DL2() {
      _bending = _axial = _torsion = _shear = true;
      _reduced_integration = false;
#ifdef _OFELI_DEBUG
      std::clog << "An instance of class Beam3DL2 is constructed.\n";
      std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    Beam3DL2(Element *el, double A, double I1, double I2);

    Beam3DL2(Element *el, double A, double I1, double I2, const Vect<double> &u, const double &time=0);

    Beam3DL2(Mesh &ms, const NodeVect<double> &u, NodeVect<double> &d);

    ~Beam3DL2() {
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class Beam3DL2 is destructed.\n";
       std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    void LMassToLHS(double coef=1.);

    void LMassToRHS(double coef=1.);

    void MassToLHS(double coef=1.) { cerr << "Error: Beam3DL2::Mass2LHS not implemented" << endl; }

    void MassToRHS(double coef=1.) { cerr << "Error: Beam3DL2::Mass2RHS not implemented" << endl; }

    void Stiffness(double coef=1.);

// Add body right-hand side term to right hand side
//    void BodyRHS(UserData<double> &ud);

    void Load(const NodeVect<double> &f);

// Return Stress in beam
//    double Stress();

    void setBending() { _bending = true; }

    void setAxial() { _axial = true; }

    void setShear() { _shear = true; }

    void setTorsion() { _torsion = true; }

    void setNoBending() { _bending = false; }

    void setNoAxial() { _axial = false; }

    void setNoShear() { _shear = false; }

    void setNoTorsion() { _torsion = false; }

    void setReducedIntegration() { _reduced_integration = true; }

    double AxialForce() const;

    Point<double> ShearForce() const;

    Point<double> BendingMoment() const;

    double TwistingMoment() const;

private:

   double        _section, _I1, _I2, _h, _mu;
   Point<double> _x[2];
   bool          _bending, _axial, _torsion, _shear, _reduced_integration;
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani