HHTGeneralized.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.1 $
00022 // $Date: 2005/12/19 22:39:21 $
00023 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/HHTGeneralized.h,v $
00024 
00025 #ifndef HHTGeneralized_h
00026 #define HHTGeneralized_h
00027 
00028 // Written: Andreas Schellenberg (andreas.schellenberg@gmx.net)
00029 // Created: 10/05
00030 // Revision: A
00031 //
00032 // Description: This file contains the class definition for HHTGeneralized.
00033 // HHTGeneralized is an algorithmic class for performing a transient analysis
00034 // using the HHTGeneralized integration scheme.
00035 //
00036 // What: "@(#) HHTGeneralized.h, revA"
00037 
00038 #include <TransientIntegrator.h>
00039 
00040 class DOF_Group;
00041 class FE_Element;
00042 class Vector;
00043 
00044 class HHTGeneralized : public TransientIntegrator
00045 {
00046 public:
00047     // constructors
00048     HHTGeneralized();
00049     HHTGeneralized(double rhoInf);
00050     HHTGeneralized(double rhoInf,
00051         double alphaM, double betaK, double betaKi, double betaKc);
00052     HHTGeneralized(double alphaI, double alphaF, double beta, double gamma);
00053     HHTGeneralized(double alphaI, double alphaF, double beta, double gamma,
00054         double alphaM, double betaK, double betaKi, double betaKc);
00055     
00056     // destructor
00057     ~HHTGeneralized();
00058     
00059     // methods which define what the FE_Element and DOF_Groups add
00060     // to the system of equation object.
00061     int formEleTangent(FE_Element *theEle);
00062     int formNodTangent(DOF_Group *theDof);        
00063     
00064     int domainChanged(void);    
00065     int newStep(double deltaT);    
00066     int revertToLastStep(void);        
00067     int update(const Vector &deltaU);
00068     int commit(void);
00069     
00070     virtual int sendSelf(int commitTag, Channel &theChannel);
00071     virtual int recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker);
00072     
00073     void Print(OPS_Stream &s, int flag = 0);        
00074     
00075 protected:
00076     
00077 private:
00078     double alphaI;
00079     double alphaF;
00080     double beta;
00081     double gamma;
00082     double deltaT;
00083     
00084     // rayleigh damping factors
00085     double alphaM;
00086     double betaK;
00087     double betaKi;
00088     double betaKc;
00089     
00090     double c1, c2, c3;                          // some constants we need to keep
00091     Vector *Ut, *Utdot, *Utdotdot;              // response quantities at time t
00092     Vector *U, *Udot, *Udotdot;                 // response quantities at time t + deltaT
00093     Vector *Ualpha, *Ualphadot, *Ualphadotdot;  // response quantities at time t+alpha*deltaT
00094 };
00095 
00096 #endif

Generated on Mon Oct 23 15:04:58 2006 for OpenSees by doxygen 1.5.0