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: 2001/03/29 05:23:32 $
00023 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/StaticIntegrator.cpp,v $
00024
00025
00026 // File: ~/analysis/integrator/StaticIntegrator.C
00027 //
00028 // Written: fmk
00029 // Created: 11/96
00030 // Revision: A
00031 //
00032 // Description: This file contains the class definition for StaticIntegrator.
00033 // StaticIntegrator is an algorithmic class for setting up the finite element
00034 // equations for a static analysis and for Incrementing the nodal displacements
00035 // with the values in the soln vector to the LinearSOE object.
00036 //
00037 // What: "@(#) StaticIntegrator.C, revA"
00038
00039 #include <StaticIntegrator.h>
00040 #include <FE_Element.h>
00041 #include <LinearSOE.h>
00042 #include <AnalysisModel.h>
00043 #include <Vector.h>
00044 #include <DOF_Group.h>
00045 #include <Channel.h>
00046 #include <FEM_ObjectBroker.h>
00047
00048StaticIntegrator::StaticIntegrator(int clasTag)
00049 :IncrementalIntegrator(clasTag)
00050 {
00051 // for subclasses
00052 }
00053
00054StaticIntegrator::~StaticIntegrator()
00055 {
00056 }
00057
00058 int
00059StaticIntegrator::formEleTangent(FE_Element *theEle)
00060 {
00061 // only elements stiffness matrix needed
00062
00063 if (statusFlag == CURRENT_TANGENT) {
00064 theEle->zeroTangent();
00065 theEle->addKtToTang();
00066 } else if (statusFlag == INITIAL_TANGENT) {
00067 theEle->zeroTangent();
00068 theEle->addKiToTang();
00069 } else if (statusFlag == CURRENT_SECANT) {
00070 theEle->zeroTangent();
00071 theEle->addKsToTang();
00072 }
00073
00074 return 0;
00075 }
00076
00077 int
00078StaticIntegrator::formEleResidual(FE_Element *theEle)
00079 {
00080 // only elements residual needed
00081 theEle->zeroResidual();
00082 theEle->addRtoResidual();
00083 return 0;
00084 }
00085
00086 int
00087StaticIntegrator::formNodTangent(DOF_Group *theDof)
00088 {
00089 // should never be called
00090 cerr << "StaticIntegrator::formNodTangent() -";
00091 cerr << " this method should never have been called!\n";
00092 return -1;
00093 }
00094
00095 int
00096StaticIntegrator::formNodUnbalance(DOF_Group *theDof)
00097 {
00098 // only nodes unbalance need be added
00099 theDof->zeroUnbalance();
00100 theDof->addPtoUnbalance();
00101 return 0;
00102 }
00103
Copyright