AlphaOS.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/AlphaOS.h,v $
00024 
00025 
00026 #ifndef AlphaOS_h
00027 #define AlphaOS_h
00028 
00029 // Written: Andreas Schellenberg (andreas.schellenberg@gmx.net)
00030 // Created: 02/05
00031 // Revision: A
00032 //
00033 // Description: This file contains the class definition for AlphaOS.
00034 // AlphaOS is an algorithmic class for performing a transient analysis
00035 // using the Alpha-Operator-Splitting integration scheme.
00036 // The parameter alpha corresponds to 1+alpha_{HHT}.
00037 //
00038 // What: "@(#) AlphaOS.h, revA"
00039 
00040 #include <TransientIntegrator.h>
00041 
00042 class DOF_Group;
00043 class FE_Element;
00044 class Vector;
00045 
00046 class AlphaOS : public TransientIntegrator
00047 {
00048 public:
00049     // constructors
00050     AlphaOS();
00051     AlphaOS(double alpha);
00052     AlphaOS(double alpha, 
00053         double alphaM, double betaK, double betaKi, double betaKc);
00054     AlphaOS(double alpha, double beta, double gamma);
00055     AlphaOS(double alpha, double beta, double gamma,
00056         double alphaM, double betaK, double betaKi, double betaKc);
00057 
00058     // destructor
00059     ~AlphaOS();
00060     
00061     // methods which define what the FE_Element and DOF_Groups add
00062     // to the system of equation object.
00063     int formEleTangent(FE_Element *theEle);
00064     int formNodTangent(DOF_Group *theDof);       
00065     
00066     int domainChanged(void);
00067     int newStep(double deltaT);    
00068     int revertToLastStep(void);
00069     int update(const Vector &deltaU);
00070     int commit(void);
00071     
00072     virtual int sendSelf(int commitTag, Channel &theChannel);
00073     virtual int recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker);
00074     
00075     void Print(OPS_Stream &s, int flag = 0);        
00076     
00077 protected:
00078     virtual int formElementResidual(void);
00079     
00080 private:
00081     double alpha;
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;     // 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