Quick Links Home Download Documentation Message Board Source Code Class Interface Bugs Tasks Search

EigenIntegrator.cpp Go to the documentation of this file. /* ****************************************************************** ** ** OpenSees - Open System for Earthquake Engineering Simulation ** ** Pacific Earthquake Engineering Research Center ** ** ** ** ** ** (C) Copyright 1999, The Regents of the University of California ** ** All Rights Reserved. ** ** ** ** Commercial use of this program without express permission of the ** ** University of California, Berkeley, is strictly prohibited. See ** ** file 'COPYRIGHT' in main directory for information on usage and ** ** redistribution, and for a DISCLAIMER OF ALL WARRANTIES. ** ** ** ** Developed by: ** ** Frank McKenna (fmckenna@ce.berkeley.edu) ** ** Gregory L. Fenves (fenves@ce.berkeley.edu) ** ** Filip C. Filippou (filippou@ce.berkeley.edu) ** ** ** ** ****************************************************************** */ // $Revision: 1.1.1.1 $ // $Date: 2000/09/15 08:23:17 $ // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/EigenIntegrator.cpp,v $ // File: ~/analysis/integrator/eigenIntegrator/EigenIntegrator.C // // Written: Jun Peng // Created: Wed Jan 27, 1999 // Revision: A // // Description: This file contains the class definition of EigenIntegrator. // EigenIntegrator is an algorithmic class for setting up the finite element // equations for a eigen problem analysis. // // This class is inheritanted from the base class of Integrator which was // created by fmk (Frank). #include "EigenIntegrator.h" #include <FE_Element.h> #include <AnalysisModel.h> #include <EigenSOE.h> #include <Vector.h> #include <DOF_Group.h> #include <FE_EleIter.h> #include <DOF_GrpIter.h> EigenIntegrator::EigenIntegrator() :Integrator(EigenINTEGRATOR_TAGS_Eigen), theSOE(0), theAnalysisModel(0) { } EigenIntegrator::~EigenIntegrator() { } void EigenIntegrator::setLinks(AnalysisModel &theModel, EigenSOE &theSysOE) { theAnalysisModel = &theModel; theSOE = &theSysOE; } int EigenIntegrator::formEleTangent(FE_Element *theEle) { if (flagK == 0) return this->formEleTangK(theEle); else return this->formEleTangM(theEle); } int EigenIntegrator::formNodTangent(DOF_Group *theDof) { return this->formNodTangM(theDof); } int EigenIntegrator::formEleResidual(FE_Element *theEle) { return 0; } int EigenIntegrator::formNodUnbalance(DOF_Group *theDof) { return 0; } int EigenIntegrator::newStep() { return 0; } int EigenIntegrator::getLastResponse(Vector &result, const ID &id) { return 0; } int EigenIntegrator::formK() { if (theAnalysisModel == 0 || theSOE == 0) { cerr << "WARNING EigenIntegrator::formK -"; cerr << " no AnalysisModel or EigenSOE has been set\n"; return -1; } // the loops to form and add the tangents are broken into two for // efficiency when performing parallel computations // loop through the FE_Elements getting them to form the tangent // FE_EleIter &theEles1 = theAnalysisModel->getFEs(); FE_Element *elePtr; flagK = 0; theSOE->zeroA(); //while((elePtr = theEles1()) != 0) // elePtr->formTangent(this); // loop through the FE_Elements getting them to add the tangent int result = 0; FE_EleIter &theEles2 = theAnalysisModel->getFEs(); while((elePtr = theEles2()) != 0) { if (theSOE->addA(elePtr->getTangent(this), elePtr->getID()) < 0) { cerr << "WARNING EigenIntegrator::formK -"; cerr << " failed in addA for ID " << elePtr->getID(); result = -2; } } return result; } int EigenIntegrator::formM() { if (theAnalysisModel == 0 || theSOE == 0) { cerr << "WARNING EigenIntegrator::formK -"; cerr << " no AnalysisModel or EigenSOE has been set\n"; return -1; } // the loops to form and add the tangents are broken into two for // efficiency when performing parallel computations // loop through the FE_Elements getting them to form the tangent // FE_EleIter &theEles1 = theAnalysisModel->getFEs(); FE_Element *elePtr; flagK = 1; theSOE->zeroM(); // while((elePtr = theEles1()) != 0) // elePtr->formTangent(this); // loop through the FE_Elements getting them to add the tangent int result = 0; FE_EleIter &theEles2 = theAnalysisModel->getFEs(); while((elePtr = theEles2()) != 0) { if (theSOE->addM(elePtr->getTangent(this), elePtr->getID()) < 0) { cerr << "WARNING EigenIntegrator::formK -"; cerr << " failed in addA for ID " << elePtr->getID(); result = -2; } } DOF_Group *dofPtr; DOF_GrpIter &theDofs = theAnalysisModel->getDOFs(); while((dofPtr = theDofs()) != 0) { // dofPtr->formTangent(this); if (theSOE->addM(dofPtr->getTangent(this),dofPtr->getID()) < 0) { cerr << "WARNING EigenIntegrator::formM -"; cerr << " failed in addM for ID " << dofPtr->getID(); result = -3; } } return result; } int EigenIntegrator::formEleTangK(FE_Element *theEle) { theEle->zeroTangent(); theEle->addKtToTang(1.0); return 0; } int EigenIntegrator::formEleTangM(FE_Element *theEle) { theEle->zeroTangent(); theEle->addMtoTang(1.0); return 0; } int EigenIntegrator::formNodTangM(DOF_Group *theDof) { theDof->zeroTangent(); theDof->addMtoTang(1.0); return 0; } int EigenIntegrator::update(const Vector &deltaU) { return 0; } EigenSOE * EigenIntegrator::getEigenSOEPtr() const { return theSOE; } AnalysisModel * EigenIntegrator::getAnalysisModelPtr() const { return theAnalysisModel; } int EigenIntegrator::sendSelf(int commitTag, Channel &theChannel) { return 0; } int EigenIntegrator::recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker) { return 0; } void EigenIntegrator::Print(ostream &s, int flag) { s << "\t EigenIntegrator: \n"; }

Copyright

Contact Us