CentralDifferenceNoDamping.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.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 |