J:/Tecgraf/prj/femoop/trunk/elc/elcplfrm.h

// -------------------------------------------------------------------------
// ElcPlFrame.h - This header file contains the definition of the plane
//                frame element class. 
// -------------------------------------------------------------------------
// Protected methods:
// -------------------------------------------------------------------------
//
// void DisplToArray( sNodeDispl *displ, double *u )
//                       
//      displ      -  nodal displacements (u,v,w)                (  in )
//      u          -  nodal displacements (array)                ( out )
//
// Transforms displacements from structure (u,v,w) form to array form.
// Obs: displ[0]..displ[1] => u[0]..u[5].
// -------------------------------------------------------------------------
//
// void CalcConfig( sNodeCoord *coord, double *l, double *a )
//                       
//      coord      -  nodal coordinates                          (  in )
//      l          -  element straight length                    ( out )
//      a          -  element angle in global system             ( out )
//
// Evaluates element straight length and global angle from given end 
// coordinates.
// -------------------------------------------------------------------------
//
// void UpdateGeom( sNodeCoord *coord, sNodeDispl *displ )
//                       
//      coord      -  nodal coordinates                       ( in/out )
//      displ      -  nodal displacements                        (  in )
//
// Updates the coordinates of end nodes.
// -------------------------------------------------------------------------
//
// void GetTrnVec( double c, double s, double *r, double *z )
//                       
//      c          -  element cosine                             (  in )
//      s          -  element sine                               (  in )
//      r          -  r vector                                   ( out )
//      z          -  z vector                                   ( out )
//
// Evaluates the transformation vectors ({r} and {z}). These vectors are
// used in the computation of element geomtric stiffness matrix.
// -------------------------------------------------------------------------
//
// void GetTrnMat( double l, double c, double s, double **T )
//                       
//      l          -  element length                             (  in )
//      c          -  element cosine                             (  in )
//      s          -  element sine                               (  in )
//      T          -  tranformation matrix                       ( out )
//
// Evaluates the transformation matrix between local (natural) and global 
// systems. In linear analysis, the given data should correspond to initial
// configuration. In nonlinear analysis, these data should correspond to
// new configuration.
// -------------------------------------------------------------------------
//
// void GetElastMat( double l, double **C )
//                       
//      l          -  element length                             (  in )
//      C          -  elastic constitutive matrix                ( out )
//
// Evaluates the elstic constitutive matrix of local (natural) system.
// -------------------------------------------------------------------------
//
// double GetAtan(double y, double x)
//                       
//      y,x        -  projections                                (  in )
//
// Returns the atan of a given angle in the range [0, 2*PI].
// -------------------------------------------------------------------------
//
// double GetIncrRot( sNodeCoord *coord, sNodeDispl *displ )
//                       
//      coord      -  nodal coordinates                          (  in )
//      displ      -  nodal displacements                        (  in )
//
// Returns the incremental rigid body rotation from the current to the new
// configuration. The element rigid body rotation is updated in the function
// InternalForce.
// -------------------------------------------------------------------------
//
// $Author: evandro $
// $Revision: 1.1 $
// $Date: 1999/04/12 14:48:47 $
// $State: Exp $
//
// -------------------------------------------------------------------------

#ifndef _ELCPLFRAME_H
#define _ELCPLFRAME_H

#include "gbldefs.h"
#include "elcbar.h"

// -------------------------------------------------------------------------
// Definition of the Plane Frame Element Class:
//
class cElcPlFrame : public cElcBar
{
 protected:
          double  _dForce[6];    // Frame forces

 public:
                  cElcPlFrame    ( void );
  virtual        ~cElcPlFrame    ( void );
          int     ThermalEffects ( void ) { return 1; }
          void    Stiffness      ( double * );
          void    GeoStiffness   ( double * );
          int     InternalForce  ( double * );
          void    CalcIntPtStress( void );
          void    NodeStress     ( sTensor * );

 protected:
          int     GetStrnTemp    ( double *, double * );
          void    CompRotMat     ( void );
          void    GetRMat        ( double ** );
          void    ElmVector      ( sNodeCoord *, double *);
          void    DisplToArray   ( sNodeDispl *, double *);
          void    CalcConfig     ( sNodeCoord * , double *, double * );
          void    GetTrnMat      ( double, double, double, double ** );
          void    GetProp        ( double *, double * );
          void    GetElastMat    ( double, double ** );
};

// -------------------------------------------------------------------------
// Definition of the Elastic Foundation (for Plane Frame) Element Class:
//
class cElcPlFrWink : public cElcPlFrame
{
 public:
                  cElcPlFrWink   ( void );
  virtual        ~cElcPlFrWink   ( void ) { }
          void    Stiffness      ( double * );
          int     InternalForce  ( double * );
          void    CalcIntPtStress( void );
};

// -------------------------------------------------------------------------
// Definition of the Corotational Plane Frame Element Class:
//
class cElcCoPlFrame : public cElcPlFrame
{
 protected:
          double  _dRot;         // Rigig body rotation 

 public:
                  cElcCoPlFrame  ( void );
  virtual        ~cElcCoPlFrame  ( void ) { }
          int     GeoNonLinear   ( void ) { return 1; }
          void    Stiffness      ( double * );
          int     InternalForce  ( double * );
          void    CalcIntPtStress( void );
          void    UpdateState    ( void );

 protected:
          void    BMatrix        ( double, double, double, double ** );
          void    GetTrnVec      ( double, double, double *, double * );
          void    UpdateGeom     ( sNodeCoord *, sNodeDispl * );
          double  GetCurrRot     ( sNodeCoord *, sNodeDispl * );
};

// -------------------------------------------------------------------------
// Definition of the Enhanced Corotational Plane Frame Element Class:
//
class cElcEnCoPlFrame : public cElcCoPlFrame
{
 public:
                  cElcEnCoPlFrame( void );
  virtual        ~cElcEnCoPlFrame( void ) { }
          void    Stiffness      ( double * );
          int     InternalForce  ( double * );
          void    CalcIntPtStress( void );

 protected:
          void    BMatrix        ( double, double, double, double, double,
                                   double, double ** );
};

// -------------------------------------------------------------------------
// Definition of the Constante Curvature Corotational Plane Frame Element 
// Class:
//
class cElcCcCoPlFrame : public cElcCoPlFrame
{
 public:
                  cElcCcCoPlFrame( void );
  virtual        ~cElcCcCoPlFrame( void ) { }
          void    Stiffness      ( double * );
          int     InternalForce  ( double * );
          void    CalcIntPtStress( void );

 protected:
          void    BMatrix        ( double, double, double, double, double,
                                   double, double ** );
};

// -------------------------------------------------------------------------
// Definition of the Corotational Plane Frame (Timoshenko Theory) Element
// Class:
//
class cElcTCPlFrame : public cElcCoPlFrame
{
 public:
                  cElcTCPlFrame  ( void );
  virtual        ~cElcTCPlFrame  ( void ) { }
          void    Stiffness      ( double * );
          int     InternalForce  ( double * );
          void    CalcIntPtStress( void );

 protected:
          void    GetProp        ( double *, double *, double * );
          void    BMatrix        ( double, double, double, double,
                                   double ** );
          void    CMatrix        ( double ** );
};
#endif

---