AlphaOS.hGo 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 |