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 #ifndef BeamColumnJoint3d_h
00036 #define BeamColumnJoint3d_h
00037
00038 #include <Element.h>
00039 #include <ID.h>
00040 #include <Matrix.h>
00041 #include <Vector.h>
00042 #include <FileStream.h>
00043 #include <OPS_Stream.h>
00044
00045 class Node;
00046 class Channel;
00047 class FEM_ObjectBroker;
00048 class Response;
00049 class Renderer;
00050 class UniaxialMaterial;
00051
00052
00053 class BeamColumnJoint3d : public Element
00054 {
00055 public:
00056
00057 BeamColumnJoint3d();
00058
00059
00060 BeamColumnJoint3d(int tag,int Nd1, int Nd2, int Nd3, int Nd4,
00061 UniaxialMaterial& theMat1, UniaxialMaterial& theMat2,
00062 UniaxialMaterial& theMat3, UniaxialMaterial& theMat4,
00063 UniaxialMaterial& theMat5, UniaxialMaterial& theMat6,
00064 UniaxialMaterial& theMat7, UniaxialMaterial& theMat8,
00065 UniaxialMaterial& theMat9, UniaxialMaterial& theMat10,
00066 UniaxialMaterial& theMat11, UniaxialMaterial& theMat12,
00067 UniaxialMaterial& theMat13);
00068
00069 BeamColumnJoint3d(int tag,int Nd1, int Nd2, int Nd3, int Nd4,
00070 UniaxialMaterial& theMat1, UniaxialMaterial& theMat2,
00071 UniaxialMaterial& theMat3, UniaxialMaterial& theMat4,
00072 UniaxialMaterial& theMat5, UniaxialMaterial& theMat6,
00073 UniaxialMaterial& theMat7, UniaxialMaterial& theMat8,
00074 UniaxialMaterial& theMat9, UniaxialMaterial& theMat10,
00075 UniaxialMaterial& theMat11, UniaxialMaterial& theMat12,
00076 UniaxialMaterial& theMat13, double Hgtfac, double Wdtfac);
00077
00078
00079 ~BeamColumnJoint3d();
00080
00081 const char *getClassType(void) const {return "BeamColumnJoint3d";};
00082
00084 bool isSubdomain(void) { return false; } ;
00085
00086
00087 int getNumExternalNodes(void) const;
00088
00089
00090 const ID &getExternalNodes(void);
00091 Node **getNodePtrs(void);
00092
00093
00094 int getNumDOF(void);
00095
00096
00097 void setDomain(Domain *theDomain);
00098
00100
00101
00102 int commitState(void);
00103
00104
00105 int revertToLastCommit(void);
00106
00107
00108 int revertToStart(void);
00109
00110
00111 int update(void);
00112
00115
00116
00117 const Matrix &getTangentStiff(void);
00118 const Matrix &getInitialStiff(void);
00119
00120
00121 const Matrix &getDamp(void);
00122 const Matrix &getMass(void);
00123
00124
00125 void zeroLoad(void);
00126 int addLoad(ElementalLoad *theLoad, double loadFactor);
00127 int addInertiaLoadToUnbalance(const Vector &accel);
00128
00129
00130 const Vector &getResistingForce(void);
00131
00132
00133 const Vector &getResistingForceIncInertia(void);
00134
00135
00136 int sendSelf(int commitTag, Channel &theChannel);
00137 int recvSelf(int commitTag, Channel &theChannel,
00138 FEM_ObjectBroker &theBroker);
00139
00140
00141 int displaySelf(Renderer &theViewer, int displayMode, float fact);
00142
00143
00144 void Print(OPS_Stream &s, int flag =0);
00145
00146
00147 Response *setResponse(const char **argv, int argc, Information &eleInfo, OPS_Stream &s);
00148 int getResponse(int responseID, Information &eleInformation);
00149
00150 int setParameter (char **argv, int argc, Information &info);
00151 int updateParameter (int parameterID, Information &info);
00152
00153
00154 protected:
00155
00156 private:
00157
00158
00159 void getGlobalDispls(Vector&) ;
00160 void getBCJoint(void);
00161 void getdg_df(void);
00162 void getdDef_du(void);
00163 void matDiag(Vector, Matrix&);
00164 void getMatResponse(Vector, Vector&, Vector&);
00165 void formR(Vector);
00166 void formK(Vector);
00167 void formTransfMat();
00168 double getStepSize(double,double,Vector,Vector,Vector,Vector,double);
00169
00170
00171 UniaxialMaterial **MaterialPtr;
00172
00173
00174 ID connectedExternalNodes;
00175 Node* nodePtr[4];
00176
00177 int nodeDbTag, dofDbTag;
00178
00179
00180 Vector Node1; Vector Node2; Vector Node3; Vector Node4;
00181 double elemActHeight;
00182 double elemActWidth;
00183 double elemWidth;
00184 double elemHeight;
00185 double HgtFac;
00186 double WdtFac;
00187
00188 Vector Uecommit;
00189 Vector UeIntcommit;
00190 Vector UeprCommit;
00191 Vector UeprIntCommit;
00192 Matrix BCJoint;
00193 Matrix dg_df;
00194 Matrix dDef_du;
00195
00196 Matrix K;
00197 Vector R;
00198
00199
00200 static Matrix Transf;
00201 static Matrix Tran;
00202
00203 };
00204
00205 #endif