CentralDifferenceNoDamping.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/02/22 22:21:10 $
00023 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/CentralDifferenceNoDamping.h,v $
00024                                                                         
00025 #ifndef CentralDifferenceNoDamping_h
00026 #define CentralDifferenceNoDamping_h
00027 
00028 // Written: fmk 
00029 // Created: 02/05
00030 // Revision: A
00031 //
00032 // Description: This file contains the class definition for CentralDifferenceNoDamping.
00033 // CentralDifferenceNoDamping is an algorithmic class for performing a transient 
00034 // analysis using the Central Difference Scheme as implemented in Dyna
00035 //       An = M(-1) (Pn - Fn)
00036 //       Vn+1/2 = Vn-1/2 + dT * An
00037 //       Dn+1   = Dn + deltaT * Vn+1/2
00038 // which is an explicit direct integration scheme as outlined in the paper:
00039 // Goudreau, G.L. and J.O. Hallquist, "Recent Developments in Large Scale Finite Element Lagrangian 
00040 // Hydrocode Technology", Journal of Computer Methods in Applied Mechanics and Engineering, 30, 1982.
00041 //
00042 // What: "@(#) CentralDifferenceNoDamping.h, revA"
00043 
00044 #include <TransientIntegrator.h>
00045 
00046 class DOF_Group;
00047 class FE_Element;
00048 class Vector;
00049 
00050 class CentralDifferenceNoDamping : public TransientIntegrator
00051 {
00052   public:
00053     CentralDifferenceNoDamping();
00054     ~CentralDifferenceNoDamping();
00055 
00056     // methods which define what the FE_Element and DOF_Groups add
00057     // to the system of equation object.
00058     int formEleTangent(FE_Element *theEle);
00059     int formNodTangent(DOF_Group *theDof);        
00060     int formEleResidual(FE_Element *theEle);
00061     int formNodUnbalance(DOF_Group *theDof);    
00062 
00063     int domainChanged(void);    
00064     int newStep(double deltaT);    
00065     int update(const Vector &deltaU);
00066 
00067     int commit(void);
00068 
00069     virtual int sendSelf(int commitTag, Channel &theChannel);
00070     virtual int recvSelf(int commitTag, Channel &theChannel, 
00071                          FEM_ObjectBroker &theBroker);
00072 
00073     void Print(OPS_Stream &s, int flag =0);        
00074     
00075   protected:
00076     
00077   private:
00078     int updateCount;    // method should only have one update per step
00079     Vector *U;          // disp response quantities at time t + deltaT
00080     Vector *Udot;       // vel response quantity at time t-1/2 delta t
00081     Vector *Udotdot;    // accel response at time t
00082     double deltaT;
00083 };
00084 
00085 #endif
00086 

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