My title

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: 2000/12/13 08:27:10 $
00023 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/Newmark.h,v $
00024                                                                         
00025                                                                         
00026 #ifndef Newmark_h
00027 #define Newmark_h
00028 
00029 // File: ~/analysis/integrator/Newmark.h
00030 // 
00031 // Written: fmk 
00032 // Created: 11/98
00033 // Revision: A
00034 //
00035 // Description: This file contains the class definition for Newmark.
00036 // Newmark is an algorithmic class for performing a transient analysis
00037 // using the Newmark integration scheme.
00038 //
00039 // What: "@(#) Newmark.h, revA"
00040 
00041 #include <TransientIntegrator.h>
00042 
00043 #include <DOF_Group.h>
00044 #include <FE_Element.h>
00045 #include <Vector.h>
00046 
00118 class Newmark : public TransientIntegrator
00119 {
00120   public:
00128     Newmark();
00129 
00136     Newmark(double gamma, double beta, bool disp = true);
00137 
00147     Newmark(double gamma, double beta, double alphaM, double betaKcurrent,
00148 
00149      double betaKinit, double betaKlastCommit, bool disp = true);
00153     ~Newmark();
00154 
00155 
00157     virtual int formEleResidual(FE_Element *theEle);
00159     virtual int formNodUnbalance(DOF_Group *theDof);    
00160 
00161     
00162     // methods which define what the FE_Element and DOF_Groups add
00164     int formEleTangent(FE_Element *theEle);
00166     int formNodTangent(DOF_Group *theDof);        
00167 
00169     int domainChanged(void);    
00171     int initialize(void);
00173     int newStep(double deltaT);    
00175     int revertToLastStep(void);        
00177     int update(const Vector &deltaU);
00178 
00180     virtual int sendSelf(int commitTag, Channel &theChannel);
00182     virtual int recvSelf(int commitTag, Channel &theChannel, 
00183     FEM_ObjectBroker &theBroker);
00184 
00186     void Print(ostream &s, int flag =0);        
00187     
00188   protected:
00189     
00190   private:
00192     bool displ;      
00193     double gamma;
00195     double beta;
00197     bool rayleighDamping; // flag indicating wheter rayleigh damping being used
00198 
00200     double alphaM, betaK, betaKi, betaKc;      
00201 
00202     double c1, c2, c3;  
00203     Vector *Ut, *Utdot, *Utdotdot; 
00204     Vector *U, *Udot, *Udotdot; 
00205     bool determiningMass; // flag to check if just want the mass contribution
00206 };
00207 
00208 #endif
00209 
00210 
00211