Quick Links Home Download Documentation Message Board Source Code Class Interface Bugs Tasks Search

EightNodeBrick.h Go to the documentation of this file. // // COPYRIGHT (C): :-)) // PROJECT: Object Oriented Finite Element Program // FILE: EightNodeBrick.h // CLASS: EightNodeBrick // MEMBER FUNCTIONS: // // MEMBER VARIABLES // // PURPOSE: Finite Element Class // RETURN: // VERSION: // LANGUAGE: C++.ver >= 3.0 // TARGET OS: DOS || UNIX || . . . // DESIGNER: Boris Jeremic, Zhaohui Yang and Xiaoyan Wu // PROGRAMMER: Boris Jeremic, Zhaohui Yang and Xiaoyan Wu // DATE: Aug. 2000 // UPDATE HISTORY: Modified from Brick3D and FourNodeQuad.hh 07/06/00 // Sept. - Oct 2000 connected to OpenSees by Zhaohui // // #ifndef EIGHTNODEBRICK_H #define EIGHTNODEBRICK_H #ifndef _bool_h #include "bool.h" #endif #include <Element.h> #include <Node.h> // Commented by Xiaoyan. Use ~/fem/node.hh 08/04/00 // Released Node.h now. Wu use Opensees's Node.09/27/00 #include <ID.h> #include <Renderer.h> #include <Domain.h> #include <string.h> #include <GaussQuadRule1d.h> #include <G3Globals.h> #include <Matrix.h> #include <Vector.h> #include <ID.h> // Xiaoyan added from brick3d.hh 07/11/00 #include <basics.h> #include <BJtensor.h> #include <nDarray.h> #include <stresst.h> #include <straint.h> //#include <node.h> //#include <mmodel.h> //#include <fe.h> #include <MatPoint3D.h> //#include <NDMaterial.h> #include <Template3Dep.h> #include <Node.h> //#include <NDMaterial.h> class EightNodeBrick: public Element { public: EightNodeBrick(int element_number, int node_numb_1, int node_numb_2, int node_numb_3, int node_numb_4, int node_numb_5, int node_numb_6, int node_numb_7, int node_numb_8, NDMaterial * Globalmmodel, double b1, double b2, double b3, double r, double p); // int dir, double surflevel); EightNodeBrick (); ~EightNodeBrick(); //Not needed Zhaohui //int Initialize(int element_number, // int node_numb_1, int node_numb_2, int node_numb_3, int node_numb_4, // int node_numb_5, int node_numb_6, int node_numb_7, int node_numb_8, // NDMaterial * Globalmmodel, double b1, double b2, double b3, // double p, double r); int getNumExternalNodes () const; const ID &getExternalNodes (); int getNumDOF (); void setDomain(Domain *theDomain); int commitState (); int revertToLastCommit (); int revertToStart (); // public methods to obtain stiffness, mass, damping and residual information // We haven't build the following functions. // All the value of K M Dmp and F are nothing. just const Matrix &getTangentStiff (); const Matrix &getSecantStiff (); const Matrix &getDamp (); const Matrix &getMass (); const Matrix &getConsMass (); void zeroLoad (); int addLoad(const Vector &addP); int addInertiaLoadToUnbalance(const Vector &accel); const Vector FormEquiBodyForce(void); const Vector &getResistingForce (); const Vector &getResistingForceIncInertia (); int sendSelf (int commitTag, Channel &theChannel); int recvSelf (int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker); int displaySelf (Renderer &theViewer, int displayMode, float fact); void Print(ostream &s, int flag =0); // Do nothing with void Print (ostream &s, int flag =0); Response *setResponse (char **argv, int argc, Information &eleInformation); int getResponse (int responseID, Information &eleInformation); protected: private: // private attributes - a copy for each object of the class int numDOF; ID connectedExternalNodes; Node *nd1Ptr; Node *nd2Ptr; Node *nd3Ptr; Node *nd4Ptr; Node *nd5Ptr; Node *nd6Ptr; Node *nd7Ptr; Node *nd8Ptr; Matrix K; Matrix C; Matrix M; Vector P; Vector Q; Vector bf; // Body forces double rho; double pressure; int order; // Order of the quadrature rule //int dir; // Direction of vertial coord. //double surflevel; // free surface level above or below this element Matrix J; Matrix L; Matrix B; // Strain interpolation matrix // // static data - single copy for all objects of the class // static G3Matrix N; // Displacement interpolation matrix // // private member functions - only objects of this class can call these // void setJacobian (double r, double s, double t); // Xiaoyan changed // double formDetJ (double r, double s, double t); // xi, eta to r and s // void formBMatrix (double r, double s, double t); // and added t // static void formNMatrix (double r, double s, double t); // 07/06/00 private: // element number (tag) double determinant_of_Jacobian; int nodes_in_brick; // number of nodes ( from 8-20 //8 now Zhaohui) NDMaterial * mmodel; int r_integration_order; int s_integration_order; int t_integration_order; // Gauss-Legendre integration order in t direction // Now I want 3D array of Material points! // MatPoint3D[r_integration_order][s_integration_order][t_integration_order] MatPoint3D ** matpoint; // pointer to array of Material Points // 3D array of material models for each Material points // NDMaterial *GPmmodel; // pointer to array of material models for Material Points // Do we need this one? //..NDMaterial *MatPoint; // Zhaohui 10-01-2000 // this tensor is tangent constitutive tensor MS definition. // It represents the state at Material point before applying strain increment // that was produced by incremental displacements // Zhaohui 10-01-2000 //..tensor * GPtangent_E; // pointer to array of constitutive tensors for Material Points // this stress tensor is start_stress from my MS definition. // It represents the state at Material point before applying strain increment // that was produced by incremental displacements // Zhaohui 10-01-2000 //..stresstensor * GPstress; // pointer to array of stresstensors for Material Points // 3D array of stresstensors models for each Material points // this stress tensor is iterative stress. // It represents the state at Material point during iterative procedure on FEM level. // Zhaohui 10-01-2000 //stresstensor * GPiterative_stress; // pointer to array of stresstensors for Material Points //double * GPq_ast_iterative; // pointer to array of iterative values of internal variable // 3D array of straintensors models for each Material points // this strain tensor is strain_increment from my MS definition. // It represents the additional strains that are to be // integrated. After numerical integration is done, the return value // ( from one of the numerical integration procedures ) // is stresstensor that should then be put in GPstress place! // Zhaohui 10-01-2000 //straintensor * GPstrain; // pointer to array of straintensors for Material Points int LM[24]; // for 8noded x 3 = 24 public: //CONSTRUCTOR //EightNodeBrick(int tag = 0, // default constructor // //int r_int_order = 0, //int s_int_order = 0, //int t_int_order = 0, // int node_numb_1 = 0, // int node_numb_2 = 0, // int node_numb_3 = 0, // int node_numb_4 = 0, // int node_numb_5 = 0, // int node_numb_6 = 0, // int node_numb_7 = 0, // int node_numb_8 = 0, // NDMaterial * Globalmmodel = 0, // //Node * GlobalNodes = 0, // maybe put NULL // //tensor * IN_tangent_E = 0, //stresstensor * INstress = 0, //stresstensor * INiterative_stress = 0, //double * INq_ast_iterative =0, //straintensor * INstrain = 0 // EPState *eps =0 // ); //void Initialize(int element_numb, // Initialize function // //short int r_int_order, // //short int s_int_order, // //short int t_int_order, // int node_numb_1, // int node_numb_2, // int node_numb_3, // int node_numb_4, // int node_numb_5, // int node_numb_6, // int node_numb_7, // int node_numb_8, // //Node * GlobalNodes, // NDMaterial * Globalmmodel, // EPState *eps ); void incremental_Update(void); tensor H_3D(double r1, double r2, double r3); tensor interp_poli_at(double r, double s, double t); tensor dh_drst_at(double r, double s, double t); //CE Dynamic Allocation for for brick3d s. EightNodeBrick & operator[](int subscript); //Finite_Element & operator[](int subscript); tensor getStiffnessTensor(void); void set_strain_stress_tensor(FILE *fp, double * u); tensor getMassTensor(void); tensor Jacobian_3D(tensor dh); tensor Jacobian_3Dinv(tensor dh); tensor Nodal_Coordinates(void); tensor incr_disp(void); tensor total_disp(void); tensor total_disp(FILE *fp, double * u); tensor stiffness_matrix(const tensor & K); tensor mass_matrix(const tensor & M); int get_global_number_of_node(int local_node_number); int get_Brick_Number(void); int * get_LM(void); //void set_LM(Node * node); // commented out temporarily 09-27-2000 Zhaohui double get_Gauss_p_c(short order, short point_numb); double get_Gauss_p_w(short order, short point_numb); tensor nodal_forces(void); tensor iterative_nodal_forces(void); tensor nodal_forces_from_stress(stresstensor & ); // returns nodal forces for given incremental strain field in an element tensor linearized_nodal_forces(void); tensor update_stress_strain(tensor & disp); void report(char *); void reportshort(char *); void reportPAK(char *); void reportpqtheta(int); void reportLM(char *); void reportTensor(char *); void reportCIPIC(char *); void reportTensorF(FILE *); // Setting initial E according to the initial pressure //void setInitE(void); //void reportStressTensorF(FILE *); }; #endif ÿ

Copyright

Contact Us