ElasticBeam3dBase.h

// -*-c++-*-
//----------------------------------------------------------------------------
//  XC program; finite element analysis code
//  for structural analysis and design.
//
//  Copyright (C)  Luis C. Pérez Tato
//
//  This program derives from OpenSees <http://opensees.berkeley.edu>
//  developed by the  «Pacific earthquake engineering research center».
//
//  Except for the restrictions that may arise from the copyright
//  of the original program (see copyright_opensees.txt)
//  XC 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 3 of the License, or 
//  (at your option) any later version.
//
//  This software 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, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------------

#ifndef ElasticBeam3dBase_h
#define ElasticBeam3dBase_h

#include "domain/mesh/element/truss_beam_column/ProtoBeam3d.h"

namespace XC {
class Channel;
class Information;
class Response;
class SectionForceDeformation;

class ElasticBeam3dBase: public ProtoBeam3d
  {
  protected: 
    CrdTransf3d *theCoordTransf; 

    void set_transf(const CrdTransf *trf);

    int sendData(Communicator &comm);
    int recvData(const Communicator &comm);
  public:
    ElasticBeam3dBase(int tag, int classTag);
    ElasticBeam3dBase(int tag, int classTag,const Material *m,const CrdTransf *trf);
    ElasticBeam3dBase(int tag, int classTag, double A, double E, double G, double Jx, double Iy, double Iz, int Nd1, int Nd2, CrdTransf3d &theTransf, double rho = 0.0);
    ElasticBeam3dBase(int tag, int classTag, double A, double alpha_y, double alpha_z, double E, double G, double Jx, double Iy, double Iz, int Nd1, int Nd2, CrdTransf3d &theTransf, double rho = 0.0);

    ElasticBeam3dBase(int tag, int classTag, int Nd1, int Nd2, SectionForceDeformation *section, CrdTransf3d &theTransf, double rho = 0.0);
    ElasticBeam3dBase(const ElasticBeam3dBase &);
    ElasticBeam3dBase &operator=(const ElasticBeam3dBase &);
    ~ElasticBeam3dBase(void);

    void setDomain(Domain *theDomain);
    
    virtual CrdTransf *getCoordTransf(void);
    virtual const CrdTransf *getCoordTransf(void) const;

    const Vector &getVDirStrongAxisGlobalCoord(bool initialGeometry) const;
    const Vector &getVDirWeakAxisGlobalCoord(bool initialGeometry) const;    

    virtual double getN1(void) const= 0;
    virtual double getN2(void) const= 0;
    inline virtual double getN(void) const //Average axial force.
      { return (getN1()+getN2())/2.0; }
    
    virtual double getMz1(void) const= 0;
    virtual double getMz2(void) const= 0;
    inline virtual double getMz(void) const
      { return (this->getMz1()+this->getMz2())/2.0; }
    
    virtual double getVy1(void) const= 0;
    virtual double getVy2(void) const= 0;
    inline virtual double getVy(void) const
      { return (getVy1()+getVy2())/2.0; }
    
    virtual double getVz1(void) const= 0;
    virtual double getVz2(void) const= 0;
    inline virtual double getVz(void) const
      { return (getVz1()+getVz2())/2.0; }
    
    virtual double getMy1(void) const= 0;
    virtual double getMy2(void) const= 0;    
    inline virtual double getMy(void) const
      { return (this->getMy1()+this->getMy2())/2.0; }
    
    virtual double getT1(void) const= 0;
    virtual double getT2(void) const= 0;
    inline virtual double getT(void) const
      { return (getT1()+getT2())/2.0; }
    
    boost::python::list getValuesAtNodes(const std::string &, bool silent= false) const;
  };
} // end of XC namespace

#endif