xc
HHT.h
1 // -*-c++-*-
2 //----------------------------------------------------------------------------
3 // XC program; finite element analysis code
4 // for structural analysis and design.
5 //
6 // Copyright (C) Luis C. Pérez Tato
7 //
8 // This program derives from OpenSees <http://opensees.berkeley.edu>
9 // developed by the «Pacific earthquake engineering research center».
10 //
11 // Except for the restrictions that may arise from the copyright
12 // of the original program (see copyright_opensees.txt)
13 // XC is free software: you can redistribute it and/or modify
14 // it under the terms of the GNU General Public License as published by
15 // the Free Software Foundation, either version 3 of the License, or
16 // (at your option) any later version.
17 //
18 // This software is distributed in the hope that it will be useful, but
19 // WITHOUT ANY WARRANTY; without even the implied warranty of
20 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
21 // GNU General Public License for more details.
22 //
23 //
24 // You should have received a copy of the GNU General Public License
25 // along with this program.
26 // If not, see <http://www.gnu.org/licenses/>.
27 //----------------------------------------------------------------------------
28 /* ****************************************************************** **
29 ** OpenSees - Open System for Earthquake Engineering Simulation **
30 ** Pacific Earthquake Engineering Research Center **
31 ** **
32 ** **
33 ** (C) Copyright 1999, The Regents of the University of California **
34 ** All Rights Reserved. **
35 ** **
36 ** Commercial use of this program without express permission of the **
37 ** University of California, Berkeley, is strictly prohibited. See **
38 ** file 'COPYRIGHT' in main directory for information on usage and **
39 ** redistribution, and for a DISCLAIMER OF ALL WARRANTIES. **
40 ** **
41 ** Developed by: **
42 ** Frank McKenna (fmckenna@ce.berkeley.edu) **
43 ** Gregory L. Fenves (fenves@ce.berkeley.edu) **
44 ** Filip C. Filippou (filippou@ce.berkeley.edu) **
45 ** **
46 ** ****************************************************************** */
47 
48 // $Revision: 1.7 $
49 // $Date: 2005/12/19 22:43:36 $
50 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/HHT.h,v $
51 
52 
53 #ifndef HHT_h
54 #define HHT_h
55 
56 // Written: Andreas Schellenberg (andreas.schellenberg@gmx.net)
57 // Created: 02/05
58 // Revision: A
59 //
60 // Description: This file contains the class definition for HHT.
61 // HHT is an algorithmic class for performing a transient analysis
62 // using the HHT integration scheme.
63 //
64 // What: "@(#) HHT.h, revA"
65 
66 #include "solution/analysis/integrator/transient/rayleigh/HHTBase.h"
67 #include "solution/analysis/integrator/transient/ResponseQuantities.h"
68 
69 namespace XC {
70 
72 //
106 class HHT: public HHTBase
107  {
108  protected:
109  int sendData(Communicator &);
110  int recvData(const Communicator &);
111 
112  friend class SolutionStrategy;
113  friend class FEM_ObjectBroker;
115  HHT(SolutionStrategy *,double alpha);
116  HHT(SolutionStrategy *,double alpha,const RayleighDampingFactors &rF);
117  HHT(SolutionStrategy *,double alpha, double beta, double gamma);
118  HHT(SolutionStrategy *,double alpha, double beta, double gamma,const RayleighDampingFactors &);
119  Integrator *getCopy(void) const;
120  public:
121  // methods which define what the FE_Element and DOF_Groups add
122  // to the system of equation object.
123  int formEleTangent(FE_Element *theEle);
124  int formNodTangent(DOF_Group *theDof);
125 
126  int domainChanged(void);
127  int newStep(double deltaT);
128  int revertToLastStep(void);
129  int update(const Vector &deltaU);
130  int commit(void);
131 
132  virtual int sendSelf(Communicator &);
133  virtual int recvSelf(const Communicator &);
134 
135  void Print(std::ostream &s, int flag = 0) const;
136  };
137 } // end of XC namespace
138 
139 #endif
int formNodTangent(DOF_Group *theDof)
This performs the following: {tabbing} while \= + while \= while \= if (RayleighDamping == false) + ...
Definition: HHT.cpp:290
void Print(std::ostream &s, int flag=0) const
The object sends to its type, the current time, , and .
Definition: HHT.cpp:525
int newStep(double deltaT)
The following are performed when this method is invoked: {enumerate} First sets the values of the th...
Definition: HHT.cpp:160
Float vector abstraction.
Definition: Vector.h:94
double alpha
should be between 0.67 and 1.0 (alpha factor corresponds to Newmark method.
Definition: HHTRayleighBase.h:47
HHT is an algorithmic class for performing a transient analysis using the Hilber-Hughes-Taylor integr...
Definition: HHT.h:106
Communication parameters between processes.
Definition: Communicator.h:66
HHT(SolutionStrategy *)
Constructor.
Definition: HHT.cpp:71
virtual int sendSelf(Communicator &)
Send object through the communicator argument.
Definition: HHT.cpp:486
int domainChanged(void)
If the size of the LinearSOE has changed, the object deletes any old Vectors created and then creates...
Definition: HHT.cpp:312
FEM_ObjectBroker is is an object broker class for the finite element method.
Definition: FEM_ObjectBroker.h:151
int sendData(Communicator &)
Send object members through the communicator argument.
Definition: HHT.cpp:472
int commit(void)
First the response quantities at the DOF_Group objects are updated with the new approximations by inv...
Definition: HHT.cpp:445
Base class for the objects that performs the integration of physical properties over the domain to fo...
Definition: Integrator.h:100
int recvData(const Communicator &)
Receives object members through the communicator argument.
Definition: HHT.cpp:479
A DOF_Group object is instantiated by the ConstraintHandler for every unconstrained node in the domai...
Definition: DOF_Group.h:107
HHTBase is an algorithmic class for performing a transient analysis using the HHTBase integration sch...
Definition: HHTBase.h:53
double gamma
factor.
Definition: HHTRayleighBase.h:49
Rayleigh damping factors.
Definition: RayleighDampingFactors.h:59
Solution strategy for the finite element problem.
Definition: SolutionStrategy.h:94
int update(const Vector &deltaU)
Invoked this first causes the object to increment the DOF_Group response quantities at time ...
Definition: HHT.cpp:387
Finite element as seen by analysis.
Definition: FE_Element.h:108
int formEleTangent(FE_Element *theEle)
This tangent for each FE_Element is defined to be , where c1,c2 and c3 were determined in the last in...
Definition: HHT.cpp:252
Open source finite element program for structural analysis.
Definition: ContinuaReprComponent.h:35
virtual int recvSelf(const Communicator &)
Receives object through the communicator argument.
Definition: HHT.cpp:501