00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046 #ifndef NLBeamColumn2d_h
00047 #define NLBeamColumn2d_h
00048
00049 #include <Element.h>
00050 #include <Node.h>
00051 #include <Matrix.h>
00052 #include <Vector.h>
00053 #include <Channel.h>
00054 #include <SectionForceDeformation.h>
00055 #include <CrdTransf2d.h>
00056
00057 #include <Response.h>
00058
00060 class NLBeamColumn2d: public Element
00061 {
00062 public:
00064 NLBeamColumn2d ();
00066 NLBeamColumn2d (int tag, int nodeI, int nodeJ,
00067 int numSections, SectionForceDeformation *sectionPtrs[],
00068 CrdTransf2d &coordTransf, double massDensPerUnitLength = 0.0,
00069 int maxNumIters = 1, double tolerance = 1e-10);
00070
00072 ~NLBeamColumn2d();
00073
00075 int getNumExternalNodes(void) const;
00077 const ID &getExternalNodes(void);
00078
00080 int getNumDOF(void);
00081
00083 void setDomain(Domain *theDomain);
00085 int commitState(void);
00087 int revertToLastCommit(void);
00089 int revertToStart(void);
00091 int update(void);
00092
00094 const Matrix &getTangentStiff(void);
00096 const Matrix &getSecantStiff(void);
00098 const Matrix &getDamp(void);
00100 const Matrix &getMass(void);
00101
00103 void zeroLoad(void);
00105 int addLoad(const Vector &load);
00107 const Vector &getResistingForce(void);
00109 const Vector &getResistingForceIncInertia(void);
00110
00112 bool isSubdomain(void);
00113
00115 int sendSelf(int cTag, Channel &theChannel);
00117 int recvSelf(int cTag, Channel &theChannel, FEM_ObjectBroker &theBroker);
00119 int displaySelf(Renderer &theViewer, int displayMode, float fact);
00120
00122 friend ostream &operator<<(ostream &s, NLBeamColumn2d &E);
00124 void Print(ostream &s, int flag =0);
00125
00127 Response *setResponse(char **argv, int argc, Information &eleInformation);
00129 int getResponse(int responseID, Information &eleInformation);
00130
00131 private:
00133 void getGlobalDispls(Vector &dg) const;
00135 void getGlobalAccels(Vector &ag) const;
00137 void getForceInterpolatMatrix(double xi, Matrix &b, const ID &code);
00139 void getDistrLoadInterpolatMatrix(double xi, Matrix &bp, const ID &code);
00141 void compSectionDisplacements(Vector sectionCoords[], Vector sectionDispls[]) const;
00143 void setSectionInterpolation (void);
00145 void initializeSectionHistoryVariables (void);
00146
00147
00149 ID connectedExternalNodes;
00150 int nSections;
00152
00153 CrdTransf2d *crdTransf;
00155
00156 int maxIters;
00157 double tol;
00158
00159 double L;
00160 double cosTheta, sinTheta;
00161 int initialFlag;
00162
00163 Node *node1Ptr, *node2Ptr;
00164
00165 Matrix K;
00166 Matrix m;
00167 Matrix d;
00168
00169 Vector P;
00170 Vector Pinert;
00171 Vector load;
00172
00173 Vector prevDistrLoad;
00174 Vector Uepr;
00175 Matrix kv;
00176 Vector Se;
00177
00178 Vector distrLoadcommit;
00179 Vector Uecommit;
00180 Matrix kvcommit;
00181 Vector Secommit;
00182
00183 Matrix *b;
00184 Matrix *bp;
00185
00186 Matrix *fs;
00187 Vector *vs;
00188 Vector *Ssr;
00189
00190 Vector *vscommit;
00191
00192 };
00193
00194 #endif
00195
00196
00197 ÿ