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 #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
00054 #include <Response.h>
00055
00057 class NLBeamColumn3d: public Element
00058 {
00059 public:
00061 NLBeamColumn3d ();
00063 NLBeamColumn3d (int tag, int nodeI, int nodeJ,
00064 int numSections, SectionForceDeformation *sectionPtrs[],
00065 CrdTransf3d &coordTransf, double massDensPerUnitLength = 0.0,
00066 int maxNumIters = 1, double tolerance = 1e-10);
00067
00069 ~NLBeamColumn3d();
00070
00072 int getNumExternalNodes(void) const;
00074 const ID &getExternalNodes(void);
00075
00077 int getNumDOF(void);
00078
00080 void setDomain(Domain *theDomain);
00082 int commitState(void);
00084 int revertToLastCommit(void);
00086 int revertToStart(void);
00088 int update(void);
00089
00091 const Matrix &getTangentStiff(void);
00093 const Matrix &getSecantStiff(void);
00095 const Matrix &getDamp(void);
00097 const Matrix &getMass(void);
00098
00100 void zeroLoad(void);
00102 int addLoad(const Vector &load);
00104 const Vector &getResistingForce(void);
00106 const Vector &getResistingForceIncInertia(void);
00107
00109 bool isSubdomain(void);
00110
00112 int sendSelf(int cTag, Channel &theChannel);
00114 int recvSelf(int cTag, Channel &theChannel, FEM_ObjectBroker &theBroker);
00116 int displaySelf(Renderer &theViewer, int displayMode, float fact);
00117
00118
00120 friend ostream &operator<<(ostream &s, NLBeamColumn3d &E);
00122 void Print(ostream &s, int flag =0);
00123
00125 Response *setResponse(char **argv, int argc, Information &eleInformation);
00127 int getResponse(int responseID, Information &eleInformation);
00128
00129 private:
00131 void getGlobalDispls(Vector &dg) const;
00133 void getGlobalAccels(Vector &ag) const;
00135 void getForceInterpolatMatrix(double xi, Matrix &b, const ID &code);
00137 void getDistrLoadInterpolatMatrix(double xi, Matrix &bp, const ID &code);
00139 void compSectionDisplacements(Vector sectionCoords[], Vector sectionDispls[]) const;
00141 void setSectionInterpolation (void);
00143 void initializeSectionHistoryVariables (void);
00144
00145
00147
00148 ID connectedExternalNodes;
00149 int nSections;
00151
00152 CrdTransf3d *crdTransf;
00154
00155 int maxIters;
00156 double tol;
00157
00158 double L;
00159 int initialFlag;
00160 bool isTorsion;
00161
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 #endif
00194
00195 ÿ