FiberSection2d.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.12 $
00022 // $Date: 2006/09/05 23:29:17 $
00023 // $Source: /usr/local/cvs/OpenSees/SRC/material/section/FiberSection2d.h,v $
00024                                                                         
00025 // Written: fmk
00026 // Created: 04/01
00027 //
00028 // Description: This file contains the class definition for 
00029 // FiberSection2d.h. FiberSection2d provides the abstraction of a 
00030 // 2d beam section discretized by fibers. The section stiffness and
00031 // stress resultants are obtained by summing fiber contributions.
00032 
00033 #ifndef FiberSection2d_h
00034 #define FiberSection2d_h
00035 
00036 #include <SectionForceDeformation.h>
00037 #include <Vector.h>
00038 #include <Matrix.h>
00039 
00040 class UniaxialMaterial;
00041 class Fiber;
00042 class Response;
00043 class SectionIntegration;
00044 
00045 class FiberSection2d : public SectionForceDeformation
00046 {
00047   public:
00048     FiberSection2d(); 
00049     FiberSection2d(int tag, int numFibers, Fiber **fibers);
00050     FiberSection2d(int tag, int numFibers, UniaxialMaterial **mats,
00051                    SectionIntegration &si);
00052     ~FiberSection2d();
00053 
00054     const char *getClassType(void) const {return "FiberSection2d";};
00055 
00056     int   setTrialSectionDeformation(const Vector &deforms); 
00057     const Vector &getSectionDeformation(void);
00058 
00059     const Vector &getStressResultant(void);
00060     const Matrix &getSectionTangent(void);
00061     const Matrix &getInitialTangent(void);
00062 
00063     int   commitState(void);
00064     int   revertToLastCommit(void);    
00065     int   revertToStart(void);
00066  
00067     SectionForceDeformation *getCopy(void);
00068     const ID &getType (void);
00069     int getOrder (void) const;
00070     
00071     int sendSelf(int cTag, Channel &theChannel);
00072     int recvSelf(int cTag, Channel &theChannel, 
00073                  FEM_ObjectBroker &theBroker);
00074     void Print(OPS_Stream &s, int flag = 0);
00075             
00076     Response *setResponse(const char **argv, int argc, 
00077                           Information &info, 
00078                           OPS_Stream &s);
00079     int getResponse(int responseID, Information &info);
00080 
00081     int addFiber(Fiber &theFiber);
00082 
00083     // AddingSensitivity:BEGIN //////////////////////////////////////////
00084     int setParameter(const char **argv, int argc, Parameter &param);
00085     const Vector& getStressResultantSensitivity(int gradNumber,
00086                                                 bool conditional);
00087     const Vector& getSectionDeformationSensitivity(int gradNumber);
00088     const Matrix& getInitialTangentSensitivity(int gradNumber);
00089     int commitSensitivity(const Vector& sectionDeformationGradient,
00090                           int gradNumber, int numGrads);
00091     // AddingSensitivity:END ///////////////////////////////////////////
00092 
00093   protected:
00094     
00095     //  private:
00096     int numFibers;                   // number of fibers in the section
00097     UniaxialMaterial **theMaterials; // array of pointers to materials
00098     double   *matData;               // data for the materials [yloc and area]
00099     double   kData[4];               // data for ks matrix 
00100     double   sData[2];               // data for s vector 
00101     
00102     double yBar;       // Section centroid
00103   
00104     SectionIntegration *sectionIntegr;
00105 
00106     static ID code;
00107 
00108     Vector e;          // trial section deformations 
00109     Vector eCommit;    // committed section deformations 
00110     Vector *s;         // section resisting forces  (axial force, bending moment)
00111     Matrix *ks;        // section stiffness
00112 
00113 // AddingSensitivity:BEGIN //////////////////////////////////////////
00114     Vector dedh; // MHS hack
00115 // AddingSensitivity:END ///////////////////////////////////////////
00116 };
00117 
00118 #endif

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