AlphaOSGeneralized.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.2 $
00022 // $Date: 2005/12/21 00:31:57 $
00023 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/AlphaOSGeneralized.h,v $
00024 
00025 #ifndef AlphaOSGeneralized_h
00026 #define AlphaOSGeneralized_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 AlphaOSGeneralized.
00033 // AlphaOSGeneralized is an algorithmic class for performing a transient analysis
00034 // using the generalized Alpha-Operator-Splitting integration scheme.
00035 // The parameters alpha correspond to 1+alpha_{HHT}.
00036 //
00037 // What: "@(#) AlphaOSGeneralized.h, revA"
00038 
00039 #include <TransientIntegrator.h>
00040 
00041 class DOF_Group;
00042 class FE_Element;
00043 class Vector;
00044 
00045 class AlphaOSGeneralized : public TransientIntegrator
00046 {
00047 public:
00048     // constructors
00049     AlphaOSGeneralized();
00050     AlphaOSGeneralized(double rhoInf);
00051     AlphaOSGeneralized(double rhoInf, 
00052         double alphaM, double betaK, double betaKi, double betaKc);
00053     AlphaOSGeneralized(double alphaI, double alphaF, double beta, double gamma);
00054     AlphaOSGeneralized(double alphaI, double alphaF, double beta, double gamma,
00055         double alphaM, double betaK, double betaKi, double betaKc);
00056 
00057     // destructor
00058     ~AlphaOSGeneralized();
00059     
00060     // methods which define what the FE_Element and DOF_Groups add
00061     // to the system of equation object.
00062     int formEleTangent(FE_Element *theEle);
00063     int formNodTangent(DOF_Group *theDof);       
00064     
00065     int domainChanged(void);
00066     int newStep(double deltaT);    
00067     int revertToLastStep(void);
00068     int update(const Vector &deltaU);
00069     int commit(void);
00070     
00071     virtual int sendSelf(int commitTag, Channel &theChannel);
00072     virtual int recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker);
00073     
00074     void Print(OPS_Stream &s, int flag = 0);        
00075     
00076 protected:
00077     virtual int formElementResidual(void);
00078     
00079 private:
00080     double alphaI;
00081     double alphaF;
00082     double beta;
00083     double gamma;
00084     double deltaT;
00085     
00086     // rayleigh damping factors
00087     double alphaM;
00088     double betaK;
00089     double betaKi;
00090     double betaKc;
00091     
00092     int updateCount;                            // method should only have one update per step
00093     double c1, c2, c3;                          // some constants we need to keep
00094     Vector *Ut, *Utdot, *Utdotdot;              // response quantities at time t
00095     Vector *U, *Udot, *Udotdot;                 // response quantities at time t+deltaT
00096     Vector *Ualpha, *Ualphadot, *Ualphadotdot;  // response quantities at time t+alpha*deltaT
00097     Vector *Upt, *Uptdot;                       // predictor quantities at time t
00098 };
00099 
00100 #endif

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