NLBeamColumn3d.h

Go to the documentation of this file.
00001 /* ****************************************************************** **
00002 **    OpenSees - Open System for Earthquake Engineering Simulation    **
00003 **          Pacific Earthquake Engineering Research Center            **
00004 **                                                                    **
00005 **                                                                    **
00006 ** (C) Copyright 1999, The Regents of the University of California    **
00007 ** All Rights Reserved.                                               **
00008 **                                                                    **
00009 ** Commercial use of this program without express permission of the   **
00010 ** University of California, Berkeley, is strictly prohibited.  See   **
00011 ** file 'COPYRIGHT'  in main directory for information on usage and   **
00012 ** redistribution,  and for a DISCLAIMER OF ALL WARRANTIES.           **
00013 **                                                                    **
00014 ** Developed by:                                                      **
00015 **   Frank McKenna (fmckenna@ce.berkeley.edu)                         **
00016 **   Gregory L. Fenves (fenves@ce.berkeley.edu)                       **
00017 **   Filip C. Filippou (filippou@ce.berkeley.edu)                     **
00018 **                                                                    **
00019 ** ****************************************************************** */
00020                                                                         
00021 // $Revision: 1.15 $
00022 // $Date: 2003/02/27 17:15:16 $
00023 // $Source: /usr/local/cvs/OpenSees/SRC/element/nonlinearBeamColumn/element/NLBeamColumn3d.h,v $
00024                                                                         
00025                                                                         
00026 // File: ~/model/element/NLBeamColumn3d.h
00027 //
00028 // Written: Remo Magalhaes de Souza on 03/99 
00029 // Revised: rms 06/99 (mass matrix)
00030 //          rms 07/99 (using setDomain)
00031 //          rms 08/99 (included P-Delta effect)
00032 //          fmk 10/99 setResponse() & getResponse()
00033 //          rms 04/00 (using CrdTransf class)
00034 //          mhs 06/00 (using new section class)
00035 //          mhs 06/00 (using new section class w/ variable dimensions)
00036 //          rms 06/00 (torsional stiffness considered at the section level)
00037 //          rms 06/00 (making copy of the sections)
00038 //          rms 06/00 (storing section history variables at the element level)
00039 //            
00040 // Purpose: This file contains the class definition for NLBeamColumn3d.
00041 // NLBeamColumn3d is a materially nonlinear flexibility based frame element.
00042 
00043 #ifndef NLBeamColumn3d_h
00044 #define NLBeamColumn3d_h
00045 
00046 #include <Element.h>
00047 #include <Node.h>
00048 #include <Matrix.h>
00049 #include <Vector.h>
00050 #include <Channel.h>
00051 #include <SectionForceDeformation.h>
00052 #include <CrdTransf3d.h>
00053 #include <GaussLobattoQuadRule1d01.h>
00054 
00055 class Response;
00056 
00057 class NLBeamColumn3d: public Element
00058 {
00059   public:
00060     NLBeamColumn3d ();
00061     NLBeamColumn3d (int tag, int nodeI, int nodeJ,  
00062                     int numSections, SectionForceDeformation *sectionPtrs[], 
00063                     CrdTransf3d &coordTransf, double massDensPerUnitLength = 0.0, 
00064                     int maxNumIters = 10, double tolerance = 1e-12);
00065     
00066     ~NLBeamColumn3d();
00067 
00068     int getNumExternalNodes(void) const;
00069     const ID &getExternalNodes(void);
00070     Node **getNodePtrs(void);
00071 
00072     int getNumDOF(void);
00073     void setDomain(Domain *theDomain);
00074 
00075     int commitState(void);
00076     int revertToLastCommit(void);        
00077     int revertToStart(void);
00078     int update(void);    
00079     
00080     const Matrix &getTangentStiff(void);
00081     const Matrix &getInitialStiff(void);
00082     const Matrix &getMass(void);
00083 
00084     void zeroLoad(void);        
00085     int addLoad(ElementalLoad *theLoad, double loadFactor);
00086     int addInertiaLoadToUnbalance(const Vector &accel);
00087 
00088     const Vector &getResistingForce(void);
00089     const Vector &getResistingForceIncInertia(void);            
00090     
00091     int sendSelf(int cTag, Channel &theChannel);
00092     int recvSelf(int cTag, Channel &theChannel, FEM_ObjectBroker &theBroker);
00093     int displaySelf(Renderer &theViewer, int displayMode, float fact);
00094 
00095     
00096     friend OPS_Stream &operator<<(OPS_Stream &s, NLBeamColumn3d &E);        
00097     void Print(OPS_Stream &s, int flag =0);    
00098 
00099     Response *setResponse(const char **argv, int argc, Information &eleInformation);
00100     int getResponse(int responseID, Information &eleInformation);
00101     
00102     int setParameter(const char **argv, int argc, Information &info);
00103     int updateParameter(int parameterID, Information &info);
00104 
00105   private:
00106     void getGlobalDispls(Vector &dg) const;
00107     void getGlobalAccels(Vector &ag) const;             
00108     void getForceInterpolatMatrix(double xi, Matrix &b, const ID &code);
00109     void getDistrLoadInterpolatMatrix(double xi, Matrix &bp, const ID &code);
00110     void compSectionDisplacements(Vector sectionCoords[], Vector sectionDispls[]) const;
00111     void initializeSectionHistoryVariables (void);
00112     
00113     
00114     // internal data
00115         
00116     ID     connectedExternalNodes; // tags of the end nodes
00117     int    nSections;              // number of sections (integration 
00118                                    // points) along the element
00119     SectionForceDeformation **sections;          // array of pointers to sections
00120     CrdTransf3d *crdTransf;        // pointer to coordinate tranformation object 
00121                                    // (performs the transformation between the global and basic system)
00122     Node *theNodes[2];
00123 
00124     double rho;                    // mass density per unit length
00125     int    maxIters;               // maximum number of local iterations
00126     double tol;                    // tolerance for relative energy norm for local iterations
00127 
00128     int    initialFlag;            // indicates if the element has been initialized
00129     bool isTorsion;
00130         
00131     Vector load;                   // equivalent nodal loads ????
00132 
00133     Matrix kv;                     // stiffness matrix in the basic system 
00134     Vector Se;                     // element resisting forces in the basic system
00135 
00136     Matrix kvcommit;               // commited stiffness matrix in the basic system
00137     Vector Secommit;               // commited element end forces in the basic system
00138 
00139     Matrix *fs;                    // array of section flexibility matrices
00140     Vector *vs;                    // array of section deformation vectors
00141     Vector *Ssr;                   // array of section resisting force vectors
00142  
00143     Vector *vscommit;              // array of commited section deformation vectors
00144 
00145     Matrix *sp;  // Applied section forces due to element loads, 5 x nSections
00146     double p0[5]; // Reactions in the basic system due to element loads
00147     Matrix *Ki;
00148 
00149     static Matrix theMatrix;
00150     static Vector theVector;
00151     static GaussLobattoQuadRule1d01 quadRule;
00152     static double workArea[];
00153 };
00154 
00155 #endif
00156 

Generated on Mon Oct 23 15:05:09 2006 for OpenSees by doxygen 1.5.0