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

EigenAnalysis.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:16 $ // $Source: /usr/local/cvs/OpenSees/SRC/analysis/analysis/EigenAnalysis.cpp,v $ // File: ~/analysis/analysis/eigenAnalysis/EigenAnalysis.C // // Written: Jun Peng // Created: Wed Jan 27, 1999 // Revision: A // // Description: This file contains the class definition of EigenAnalysis. // EigenAnalysis is a subclass of Analysis, it is used to perform the // eigen vlaue analysis on the FE_Model. // // This class is inheritanted from the base class of Analysis // which was created by fmk (Frank). #include <EigenAnalysis.h> #include <EigenAlgorithm.h> #include <AnalysisModel.h> #include <EigenSOE.h> #include <DOF_Numberer.h> #include <ConstraintHandler.h> #include <EigenIntegrator.h> #include <Domain.h> #include <Timer.h> #include <FE_Element.h> #include <DOF_Group.h> #include <FE_EleIter.h> #include <DOF_GrpIter.h> #include <Matrix.h> #include <ID.h> #include <Graph.h> EigenAnalysis::EigenAnalysis(Domain &the_Domain, ConstraintHandler &theHandler, DOF_Numberer &theNumberer, AnalysisModel &theModel, EigenAlgorithm &theAlgo, EigenSOE &theEigenSOE, EigenIntegrator &theEigenIntegrator) :Analysis(the_Domain), theConstraintHandler(&theHandler), theDOF_Numberer(&theNumberer), theAnalysisModel(&theModel), theAlgorithm(&theAlgo), theSOE(&theEigenSOE), theIntegrator(&theEigenIntegrator), domainStamp(0) { // first set up the links needed by the elements in the aggregation. theAnalysisModel->setLinks(the_Domain); theConstraintHandler->setLinks(the_Domain, theModel, theEigenIntegrator); theDOF_Numberer->setLinks(theModel); theIntegrator->setLinks(theModel, theEigenSOE); theAlgorithm->setLinks(theModel, theEigenIntegrator, theEigenSOE); } EigenAnalysis::~EigenAnalysis() { // do nothing now. this->clearAll(); } void EigenAnalysis::clearAll(void) { // invoke the destructor on all the objects in the aggregation delete theAnalysisModel; delete theConstraintHandler; delete theDOF_Numberer; delete theIntegrator; delete theAlgorithm; delete theSOE; } int EigenAnalysis::analyze(int numModes) { int result = 0; Domain *the_Domain = this->getDomainPtr(); // check for change in Domain since last step. As a change can // occur in a commit() in a domaindecomp with load balancing // this must now be inside the loop int stamp = the_Domain->hasDomainChanged(); if (stamp != domainStamp) { domainStamp = stamp; result = this->domainChanged(); if (result < 0) { cerr << "EigenAnalysis::analyze() - domainChanged failed\n"; return -1; } } result = theIntegrator->newStep(); if (result < 0) { cerr << "EigenAnalysis::analyze() - integrator failed\n"; return -2; } result = theAlgorithm->solveCurrentStep(numModes); if (result < 0) { cerr << "EigenAnalysis::analyze() - algorithm failed\n"; return -3; } return 0; } int EigenAnalysis::domainChanged() { theAnalysisModel->clearAll(); theConstraintHandler->clearAll(); theConstraintHandler->handle(); theDOF_Numberer->numberDOF(); Graph &theGraph = theAnalysisModel->getDOFGraph(); theSOE->setSize(theGraph); theIntegrator->domainChanged(); theAlgorithm->domainChanged(); return 0; } int EigenAnalysis::setAlgorithm(EigenAlgorithm &theAlgo) { cerr << "EigenAnalysis::setAlgorithm() - does nothing yet\n"; return 0; } int EigenAnalysis::setIntegrator(EigenIntegrator &theIntegrator) { cerr << "EigenAnalysis::setIntegrator() - does nothing yet\n"; return 0; } int EigenAnalysis::setEigenSOE(EigenSOE &theSOE) { cerr << "EigenAnalysis::setEigenSOE() - does nothing yet\n"; return 0; } ConstraintHandler * EigenAnalysis::getConstraintHandlerPtr() const { return theConstraintHandler; } DOF_Numberer * EigenAnalysis::getDOF_NumbererPtr() const { return theDOF_Numberer; } AnalysisModel * EigenAnalysis::getAnalysisModelPtr() const { return theAnalysisModel; } EigenAlgorithm * EigenAnalysis::getEigenAlgorithm() const { return theAlgorithm; } EigenSOE * EigenAnalysis::getEigenSOE() const { return theSOE; } EigenIntegrator * EigenAnalysis::getEigenIntegrator() const { return theIntegrator; }

Copyright

Contact Us