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

StaticAnalysis.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.3 $ // $Date: 2001/05/03 06:15:34 $ // $Source: /usr/local/cvs/OpenSees/SRC/analysis/analysis/StaticAnalysis.cpp,v $ // File: ~/analysis/analysis/StaticAnalysis.C // // Written: fmk // Created: 11/96 // Revision: A // // Description: This file contains the implementation of StaticAnalysis. // // What: "@(#) StaticAnalysis.C, revA" #include <StaticAnalysis.h> #include <EquiSolnAlgo.h> #include <AnalysisModel.h> #include <LinearSOE.h> #include <DOF_Numberer.h> #include <ConstraintHandler.h> #include <StaticIntegrator.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> // Constructor // sets theModel and theSysOFEqn to 0 and the Algorithm to the one supplied StaticAnalysis::StaticAnalysis(Domain &the_Domain, ConstraintHandler &theHandler, DOF_Numberer &theNumberer, AnalysisModel &theModel, EquiSolnAlgo &theSolnAlgo, LinearSOE &theLinSOE, StaticIntegrator &theStaticIntegrator) :Analysis(the_Domain), theConstraintHandler(&theHandler), theDOF_Numberer(&theNumberer), theAnalysisModel(&theModel), theAlgorithm(&theSolnAlgo), theSOE(&theLinSOE), theIntegrator(&theStaticIntegrator), domainStamp(0) { // first we set up the links needed by the elements in the // aggregation theAnalysisModel->setLinks(the_Domain); theConstraintHandler->setLinks(the_Domain,theModel,theStaticIntegrator); theDOF_Numberer->setLinks(theModel); theIntegrator->setLinks(theModel,theLinSOE); theAlgorithm->setLinks(theModel,theStaticIntegrator,theLinSOE); } StaticAnalysis::~StaticAnalysis() { } void StaticAnalysis::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 StaticAnalysis::analyze(int numSteps) { int result = 0; Domain *the_Domain = this->getDomainPtr(); for (int i=0; i<numSteps; i++) { // 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 << "StaticAnalysis::analyze() - domainChanged failed"; cerr << " at step " << i << " of " << numSteps << endl; return -1; } } result = theIntegrator->newStep(); if (result < 0) { cerr << "StaticAnalysis::analyze() - the Integrator failed"; cerr << " at iteration: " << i << " with domain at load factor "; cerr << the_Domain->getCurrentTime() << endl; the_Domain->revertToLastCommit(); return -2; } result = theAlgorithm->solveCurrentStep(); if (result < 0) { cerr << "StaticAnalysis::analyze() - the Algorithm failed"; cerr << " at iteration: " << i << " with domain at load factor "; cerr << the_Domain->getCurrentTime() << endl; the_Domain->revertToLastCommit(); theIntegrator->revertToLastStep(); return -3; } result = theIntegrator->commit(); if (result < 0) { cerr << "StaticAnalysis::analyze() - "; cerr << "the Integrator failed to commit"; cerr << " at iteration: " << i << " with domain at load factor "; cerr << the_Domain->getCurrentTime() << endl; the_Domain->revertToLastCommit(); theIntegrator->revertToLastStep(); return -4; } } return 0; } int StaticAnalysis::initialize(void) { Domain *the_Domain = this->getDomainPtr(); // check if domain has undergone change int stamp = the_Domain->hasDomainChanged(); if (stamp != domainStamp) { domainStamp = stamp; if (this->domainChanged() < 0) { cerr << "DirectIntegrationAnalysis::initialize() - domainChanged() failed\n"; return -1; } } if (theIntegrator->initialize() < 0) { cerr << "DirectIntegrationAnalysis::initialize() - integrator initialize() failed\n"; return -2; } else theIntegrator->commit(); return 0; } int StaticAnalysis::domainChanged(void) { int result = 0; // Timer theTimer; theTimer.start(); theAnalysisModel->clearAll(); theConstraintHandler->clearAll(); // theTimer.pause(); // cout << "StaticAnalysis::clearAll() " << theTimer.getReal(); // cout << theTimer.getCPU() << endl; // theTimer.start(); // now we invoke handle() on the constraint handler which // causes the creation of FE_Element and DOF_Group objects // and their addition to the AnalysisModel. result = theConstraintHandler->handle(); if (result < 0) { cerr << "StaticAnalysis::handle() - "; cerr << "ConstraintHandler::handle() failed"; return -1; } // we now invoke number() on the numberer which causes // equation numbers to be assigned to all the DOFs in the // AnalysisModel. result = theDOF_Numberer->numberDOF(); if (result < 0) { cerr << "StaticAnalysis::handle() - "; cerr << "DOF_Numberer::numberDOF() failed"; return -2; } // we invoke setSize() on the LinearSOE which // causes that object to determine its size Graph &theGraph = theAnalysisModel->getDOFGraph(); result = theSOE->setSize(theGraph); if (result < 0) { cerr << "StaticAnalysis::handle() - "; cerr << "LinearSOE::setSize() failed"; return -3; } // finally we invoke domainChanged on the Integrator and Algorithm // objects .. informing them that the model has changed result = theIntegrator->domainChanged(); if (result < 0) { cerr << "StaticAnalysis::setAlgorithm() - "; cerr << "Integrator::domainChanged() failed"; return -4; } result = theAlgorithm->domainChanged(); if (result < 0) { cerr << "StaticAnalysis::setAlgorithm() - "; cerr << "Algorithm::domainChanged() failed"; return -5; } // if get here successfull return 0; } int StaticAnalysis::setAlgorithm(EquiSolnAlgo &theNewAlgorithm) { int result = 0; // invoke the destructor on the old one if (theAlgorithm != 0) delete theAlgorithm; // first set the links needed by the Algorithm theAlgorithm = &theNewAlgorithm; theAlgorithm->setLinks(*theAnalysisModel,*theIntegrator,*theSOE); // invoke domainChanged() either indirectly or directly Domain *the_Domain = this->getDomainPtr(); int stamp = the_Domain->hasDomainChanged(); if (stamp != domainStamp) { domainStamp = stamp; result = this->domainChanged(); if (result < 0) { cerr << "StaticAnalysis::setAlgorithm() - domainChanged() failed"; return -1; } } else { result = theAlgorithm->domainChanged(); if (result < 0) { cerr << "StaticAnalysis::setAlgorithm() - "; cerr << "algorithm::domainChanged() failed"; return -2; } } return 0; } int StaticAnalysis::setIntegrator(StaticIntegrator &theNewIntegrator) { int result = 0; // invoke the destructor on the old one if (theIntegrator != 0) { delete theIntegrator; } // set the links needed by the other objects in the aggregation Domain *the_Domain = this->getDomainPtr(); theIntegrator = &theNewIntegrator; theIntegrator->setLinks(*theAnalysisModel,*theSOE); theConstraintHandler->setLinks(*the_Domain,*theAnalysisModel,*theIntegrator); theAlgorithm->setLinks(*theAnalysisModel,*theIntegrator,*theSOE); // invoke domainChanged() either indirectly or directly int stamp = the_Domain->hasDomainChanged(); if (stamp != domainStamp) { domainStamp = stamp; result = this->domainChanged(); if (result < 0) { cerr << "StaticAnalysis::setAlgorithm() - domainChanged() failed"; return -1; } } else { result = theIntegrator->domainChanged(); if (result < 0) { cerr << "StaticAnalysis::setAlgorithm() - "; cerr << "Integrator::domainChanged() failed"; return -2; } } return 0; } int StaticAnalysis::setLinearSOE(LinearSOE &theNewSOE) { int result = 0; // invoke the destructor on the old one if (theSOE != 0) delete theSOE; // set the links needed by the other objects in the aggregation theSOE = &theNewSOE; theIntegrator->setLinks(*theAnalysisModel,*theSOE); theAlgorithm->setLinks(*theAnalysisModel,*theIntegrator,*theSOE); // set the size either indirectly or directly Domain *the_Domain = this->getDomainPtr(); int stamp = the_Domain->hasDomainChanged(); if (stamp != domainStamp) { domainStamp = stamp; result = this->domainChanged(); if (result < 0) { cerr << "StaticAnalysis::setAlgorithm() - domainChanged failed"; return -1; } } else { Graph &theGraph = theAnalysisModel->getDOFGraph(); theSOE->setSize(theGraph); if (result < 0) { cerr << "StaticAnalysis::setAlgorithm() - "; cerr << "LinearSOE::setSize() failed\n"; return -2; } } return 0; }

Copyright

Contact Us