PEER NEES NEESit

• HOME
• USER
• DEVELOPER
• PROJECTS
• SUPPORT
• SITE MAP

• Dev Doc
• API
• Source
• Download
• Builds
• Bug Reports
• Message Board

DistributedDisplacementControl.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 $ // $Date: 2005/11/29 21:59:49 $ // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/DistributedDisplacementControl.cpp,v $ // Written: fmk // Created: 07/98 // // Description: This file contains the class definition for DistributedDisplacementControl. // DistributedDisplacementControl is an algorithmic class for perfroming a static analysis // using the arc length scheme, that is within a load step the follwing // constraint is enforced: dU^TdU + alpha^2*dLambda^2 = DistributedDisplacementControl^2 // where dU is change in nodal displacements for step, dLambda is // change in applied load and DistributedDisplacementControl is a control parameter. // // What: "@(#) DistributedDisplacementControl.C, revA" #include <DistributedDisplacementControl.h> #include <AnalysisModel.h> #include <LinearSOE.h> #include <Vector.h> #include <Channel.h> #include <math.h> #include <Domain.h> #include <Node.h> #include <DOF_Group.h> #include <ID.h> DistributedDisplacementControl::DistributedDisplacementControl(int node, int dof, double increment, int numIncr, double min, double max) :StaticIntegrator(INTEGRATOR_TAGS_DistributedDisplacementControl), processID(0), theChannels(0), numChannels(0), theNode(node), theDof(dof), theIncrement(increment), theDofID(0), deltaUhat(0), deltaUbar(0), deltaU(0), deltaUstep(0), phat(0), deltaLambdaStep(0.0), currentLambda(0.0), specNumIncrStep(numIncr), numIncrLastStep(numIncr), minIncrement(min), maxIncrement(max) { // to avoid divide-by-zero error on first update() ensure numIncr != 0 if (numIncr == 0) { opserr << "WARNING DistributedDisplacementControl::DistributedDisplacementControl() -"; opserr << " numIncr set to 0, 1 assumed\n"; specNumIncrStep = 1.0; numIncrLastStep = 1.0; } } DistributedDisplacementControl::DistributedDisplacementControl() :StaticIntegrator(INTEGRATOR_TAGS_DistributedDisplacementControl), processID(0), theChannels(0), numChannels(0), theNode(0), theDof(0), theIncrement(0), theDofID(0), deltaUhat(0), deltaUbar(0), deltaU(0), deltaUstep(0), phat(0), deltaLambdaStep(0.0), currentLambda(0.0), specNumIncrStep(0), numIncrLastStep(0), minIncrement(0), maxIncrement(0) { } DistributedDisplacementControl::~DistributedDisplacementControl() { // delete any vector object created if (deltaUhat != 0) delete deltaUhat; if (deltaU != 0) delete deltaU; if (deltaUstep != 0) delete deltaUstep; if (deltaUbar != 0) delete deltaUbar; if (phat != 0) delete phat; if (theChannels != 0) delete [] theChannels; } int DistributedDisplacementControl::newStep(void) { // get pointers to AnalysisModel and LinearSOE AnalysisModel *theModel = this->getAnalysisModelPtr(); LinearSOE *theLinSOE = this->getLinearSOEPtr(); if (theModel == 0 || theLinSOE == 0) { opserr << "WARNING DistributedDisplacementControl::newStep() "; opserr << "No AnalysisModel or LinearSOE has been set\n"; return -1; } // determine increment for this iteration double factor = specNumIncrStep/numIncrLastStep; theIncrement *=factor; if (theIncrement < minIncrement) theIncrement = minIncrement; else if (theIncrement > maxIncrement) theIncrement = maxIncrement; // get the current load factor currentLambda = theModel->getCurrentDomainTime(); // determine dUhat this->formTangent(); if (processID == 0) theLinSOE->setB(*phat); else theLinSOE->zeroB(); theLinSOE->solve(); (*deltaUhat) = theLinSOE->getX(); Vector &dUhat = *deltaUhat; double dUahat = dUhat(theDofID); if (dUahat == 0.0) { opserr << "WARNING DistributedDisplacementControl::newStep() "; opserr << "dUahat is zero -- zero reference displacement at control node DOF\n"; return -1; } // determine delta lambda(1) == dlambda double dLambda = theIncrement/dUahat; deltaLambdaStep = dLambda; currentLambda += dLambda; // opserr << "DistributedDisplacementControl: " << dUahat << " " << theDofID << endln; // opserr << "DistributedDisplacementControl::newStep() : " << deltaLambdaStep << endln; // determine delta U(1) == dU (*deltaU) = dUhat; (*deltaU) *= dLambda; (*deltaUstep) = (*deltaU); // update model with delta lambda and delta U theModel->incrDisp(*deltaU); theModel->applyLoadDomain(currentLambda); if (theModel->updateDomain() < 0) { opserr << "DistributedDisplacementControl::newStep - model failed to update for new dU\n"; return -1; } numIncrLastStep = 0; return 0; } int DistributedDisplacementControl::update(const Vector &dU) { AnalysisModel *theModel = this->getAnalysisModelPtr(); LinearSOE *theLinSOE = this->getLinearSOEPtr(); if (theModel == 0 || theLinSOE == 0) { opserr << "WARNING DistributedDisplacementControl::update() "; opserr << "No AnalysisModel or LinearSOE has been set\n"; return -1; } (*deltaUbar) = dU; // have to do this as the SOE is gonna change double dUabar = (*deltaUbar)(theDofID); // determine dUhat if (processID == 0) theLinSOE->setB(*phat); else theLinSOE->zeroB(); theLinSOE->solve(); (*deltaUhat) = theLinSOE->getX(); double dUahat = (*deltaUhat)(theDofID); if (dUahat == 0.0) { opserr << "WARNING DistributedDisplacementControl::update() "; opserr << "dUahat is zero -- zero reference displacement at control node DOF\n"; return -1; } // determine delta lambda(1) == dlambda double dLambda = -dUabar/dUahat; // determine delta U(i) (*deltaU) = (*deltaUbar); deltaU->addVector(1.0, *deltaUhat,dLambda); // update dU and dlambda (*deltaUstep) += *deltaU; deltaLambdaStep += dLambda; currentLambda += dLambda; // update the model theModel->incrDisp(*deltaU); theModel->applyLoadDomain(currentLambda); if (theModel->updateDomain() < 0) { opserr << "DistributedDisplacementControl::update - model failed to update for new dU\n"; return -1; } // set the X soln in linearSOE to be deltaU for convergence Test theLinSOE->setX(*deltaU); numIncrLastStep++; return 0; } int DistributedDisplacementControl::domainChanged(void) { // we first create the Vectors needed AnalysisModel *theModel = this->getAnalysisModelPtr(); LinearSOE *theLinSOE = this->getLinearSOEPtr(); if (theModel == 0 || theLinSOE == 0) { opserr << "WARNING DistributedDisplacementControl::update() "; opserr << "No AnalysisModel or LinearSOE has been set\n"; return -1; } //int size = theModel->getNumEqn(); // ask model in case N+1 space const Vector &b = theLinSOE->getB(); // ask model in case N+1 space int size = b.Size(); // first we determine the id of the nodal dof Domain *theDomain = theModel->getDomainPtr(); if (theDomain == 0) { opserr << "BUG WARNING DistributedDisplacementControl::domainChanged() - no Domain associated!!"; return -1; } theDofID = -1; Node *theNodePtr = theDomain->getNode(theNode); if (theNodePtr != 0) { DOF_Group *theGroup = theNodePtr->getDOF_GroupPtr(); if (theGroup == 0) { opserr << "BUG DistributedDisplacementControl::domainChanged() - no DOF_Group associated with the node!!\n"; return -1; } const ID &theID = theGroup->getID(); if (theDof < 0 || theDof >= theID.Size()) { opserr << "DistributedDisplacementControl::domainChanged() - not a valid dof " << theDof << endln; return -1; } theDofID = theID(theDof); if (theDofID < 0) { opserr << "DistributedDisplacementControl::domainChanged() - constrained dof not a valid a dof\n";; return -1; } } static ID data(1); if (processID != 0) { Channel *theChannel = theChannels[0]; data(0) = theDofID; theChannel->sendID(0, 0, data); theChannel->recvID(0, 0, data); theDofID = data(0); } // if main domain, collect all theDofID if not -1 then this is the value & send value out else { for (int j=0; j<numChannels; j++) { Channel *theChannel = theChannels[j]; theChannel->recvID(0, 0, data); if (data(0) != -1) theDofID = data(0); } for (int k=0; k<numChannels; k++) { Channel *theChannel = theChannels[k]; data(0) = theDofID; theChannel->sendID(0, 0, data); } } if (deltaUhat == 0 || deltaUhat->Size() != size) { // create new Vector if (deltaUhat != 0) delete deltaUhat; // delete the old deltaUhat = new Vector(size); if (deltaUhat == 0 || deltaUhat->Size() != size) { // check got it opserr << "FATAL DistributedDisplacementControl::domainChanged() - ran out of memory for"; opserr << " deltaUhat Vector of size " << size << endln; exit(-1); } } if (deltaUbar == 0 || deltaUbar->Size() != size) { // create new Vector if (deltaUbar != 0) delete deltaUbar; // delete the old deltaUbar = new Vector(size); if (deltaUbar == 0 || deltaUbar->Size() != size) { // check got it opserr << "FATAL DistributedDisplacementControl::domainChanged() - ran out of memory for"; opserr << " deltaUbar Vector of size " << size << endln; exit(-1); } } if (deltaU == 0 || deltaU->Size() != size) { // create new Vector if (deltaU != 0) delete deltaU; // delete the old deltaU = new Vector(size); if (deltaU == 0 || deltaU->Size() != size) { // check got it opserr << "FATAL DistributedDisplacementControl::domainChanged() - ran out of memory for"; opserr << " deltaU Vector of size " << size << endln; exit(-1); } } if (deltaUstep == 0 || deltaUstep->Size() != size) { if (deltaUstep != 0) delete deltaUstep; deltaUstep = new Vector(size); if (deltaUstep == 0 || deltaUstep->Size() != size) { opserr << "FATAL DistributedDisplacementControl::domainChanged() - ran out of memory for"; opserr << " deltaUstep Vector of size " << size << endln; exit(-1); } } if (phat == 0 || phat->Size() != size) { if (phat != 0) delete phat; phat = new Vector(size); if (phat == 0 || phat->Size() != size) { opserr << "FATAL DistributedDisplacementControl::domainChanged() - ran out of memory for"; opserr << " phat Vector of size " << size << endln; exit(-1); } } // now we have to determine phat // do this by incrementing lambda by 1, applying load // and getting phat from unbalance. currentLambda = theModel->getCurrentDomainTime(); currentLambda += 1.0; theModel->applyLoadDomain(currentLambda); this->formUnbalance(); // NOTE: this assumes unbalance at last was 0 (*phat) = theLinSOE->getB(); currentLambda -= 1.0; theModel->setCurrentDomainTime(currentLambda); // check there is a reference load int haveLoad = 0; for (int i=0; i<size; i++) if ( (*phat)(i) != 0.0 ) { haveLoad = 1; i = size; } if (haveLoad == 0) { opserr << "WARNING DistributedDisplacementControl::domainChanged() - zero reference load"; return -1; } if (theDofID == -1) { opserr << "DistributedDisplacementControl::setSize() - failed to find valid dof - are the node tag and dof values correct?\n"; return -1; } return 0; } int DistributedDisplacementControl::sendSelf(int cTag, Channel &theChannel) { int sendID =0; // if P0 check if already sent. If already sent use old processID; if not allocate a new process // id for remote part of object, enlarge channel * to hold a channel * for this remote object. // if not P0, send current processID if (processID == 0) { // check if already using this object bool found = false; for (int i=0; i<numChannels; i++) if (theChannels[i] == &theChannel) { sendID = i+1; found = true; } // if new object, enlarge Channel pointers to hold new channel * & allocate new ID if (found == false) { int nextNumChannels = numChannels + 1; Channel **nextChannels = new Channel *[nextNumChannels]; if (nextChannels == 0) { opserr << "DistributedDisplacementControl::sendSelf() - failed to allocate channel array of size: " << nextNumChannels << endln; return -1; } for (int i=0; i<numChannels; i++) nextChannels[i] = theChannels[i]; nextChannels[numChannels] = &theChannel; numChannels = nextNumChannels; if (theChannels != 0) delete [] theChannels; theChannels = nextChannels; // allocate new processID for remote object sendID = numChannels; } } else sendID = processID; // send remotes processID & info about node, dof and numIncr static ID idData(3); idData(0) = sendID; idData(1) = theNode; idData(2) = theDof; int res = theChannel.sendID(0, cTag, idData); if (res < 0) { opserr <<"WARNING DistributedDisplacementControl::sendSelf() - failed to send data\n"; return -1; } static Vector dData(5); dData(0) = theIncrement; dData(1) = minIncrement; dData(2) = maxIncrement; dData(3) = specNumIncrStep; dData(4) = numIncrLastStep; res = theChannel.sendVector(0, cTag, dData); if (res < 0) { opserr <<"WARNING DistributedDisplacementControl::recvSelf() - failed to recv vector data\n"; return -1; } return 0; } int DistributedDisplacementControl::recvSelf(int cTag, Channel &theChannel, FEM_ObjectBroker &theBroker) { static ID idData(3); int res = theChannel.recvID(0, cTag, idData); if (res < 0) { opserr <<"WARNING DistributedDisplacementControl::recvSelf() - failed to recv id data\n"; return -1; } processID = idData(0); theNode = idData(1); theDof = idData(2); static Vector dData(5); res = theChannel.recvVector(0, cTag, dData); if (res < 0) { opserr <<"WARNING DistributedDisplacementControl::recvSelf() - failed to recv vector data\n"; return -1; } theIncrement = dData(0); minIncrement = dData(1); maxIncrement = dData(2); specNumIncrStep = dData(3); numIncrLastStep = dData(4); // set the Channel & processID numChannels = 1; theChannels = new Channel *[1]; theChannels[0] = &theChannel; return 0; } void DistributedDisplacementControl::Print(OPS_Stream &s, int flag) { // TO FINISH } Generated on Mon Oct 23 15:04:58 2006 for OpenSees by  1.5.0

---

opensees-support @ berkeley.edu

Supported by the National Science Foundation

©2006, UC Regents