PEER NEES NEESit

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

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

NodeRecorder.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.26 $ // $Date: 2006/08/07 22:18:06 $ // $Source: /usr/local/cvs/OpenSees/SRC/recorder/NodeRecorder.cpp,v $ // Written: fmk // // Description: This file contains the class definition for NodeRecorder. // A NodeRecorder is used to record the specified dof responses // at a collection of nodes over an analysis. (between commitTag of 0 and // last commitTag). // // What: "@(#) NodeRecorder.C, revA" #include <NodeRecorder.h> #include <Domain.h> #include <Node.h> #include <Vector.h> #include <ID.h> #include <Matrix.h> #include <FE_Datastore.h> #include <FEM_ObjectBroker.h> #include <string.h> NodeRecorder::NodeRecorder() :Recorder(RECORDER_TAGS_NodeRecorder), theDofs(0), theNodalTags(0), theNodes(0), response(0), theDomain(0), theOutputHandler(0), echoTimeFlag(true), dataFlag(0), deltaT(0), nextTimeStampToRecord(0.0), sensitivity(0), initializationDone(false), numValidNodes(0) { } NodeRecorder::NodeRecorder(const ID &dofs, const ID &nodes, int psensitivity, const char *dataToStore, Domain &theDom, OPS_Stream &theOutputHandler, double dT, bool timeFlag) :Recorder(RECORDER_TAGS_NodeRecorder), theDofs(0), theNodalTags(0), theNodes(0), response(1 + nodes.Size()*dofs.Size()), theDomain(&theDom), theOutputHandler(&theOutputHandler), echoTimeFlag(timeFlag), dataFlag(0), deltaT(dT), nextTimeStampToRecord(0.0), sensitivity(psensitivity), initializationDone(false), numValidNodes(0) { // // store copy of dof's to be recorder, verifying dof are valid, i.e. >= 0 // int numDOF = dofs.Size(); if (numDOF != 0) { theDofs = new ID(numDOF); int count = 0; int i; for (i=0; i<numDOF; i++) { int dof = dofs(i); if (dof >= 0) { (*theDofs)[count] = dof; count++; } else { opserr << "NodeRecorder::NodeRecorder - invalid dof " << dof; opserr << " will be ignored\n"; } } } // // create memory to hold nodal ID's (neeed parallel) // int numNode = nodes.Size(); if (numNode != 0) { theNodalTags = new ID(nodes); if (theNodalTags == 0) { opserr << "NodeRecorder::NodeRecorder - out of memory\n"; } } // // set the data flag used as a switch to get the response in a record // if (dataToStore == 0 || (strcmp(dataToStore, "disp") == 0)) { dataFlag = 0; } else if ((strcmp(dataToStore, "vel") == 0)) { dataFlag = 1; } else if ((strcmp(dataToStore, "accel") == 0)) { dataFlag = 2; } else if ((strcmp(dataToStore, "incrDisp") == 0)) { dataFlag = 3; } else if ((strcmp(dataToStore, "incrDeltaDisp") == 0)) { dataFlag = 4; } else if ((strcmp(dataToStore, "unbalance") == 0)) { dataFlag = 5; } else if ((strcmp(dataToStore, "unbalanceInclInertia") == 0)) { dataFlag = 6; } else if ((strcmp(dataToStore, "reaction") == 0)) { dataFlag = 7; } else if (((strcmp(dataToStore, "reactionIncInertia") == 0)) || ((strcmp(dataToStore, "reactionIncludingInertia") == 0))) { dataFlag = 8; } else if ((strncmp(dataToStore, "eigen",5) == 0)) { int mode = atoi(&(dataToStore[5])); if (mode > 0) dataFlag = 10 + mode; else dataFlag = 6; } else if ((strncmp(dataToStore, "sensitivity",11) == 0)) { int grad = atoi(&(dataToStore[11])); if (grad > 0) dataFlag = 1000 + grad; else dataFlag = 6; } else if ((strncmp(dataToStore, "velSensitivity",14) == 0)) { int grad = atoi(&(dataToStore[14])); if (grad > 0) dataFlag = 2000 + grad; else dataFlag = 6; } else if ((strncmp(dataToStore, "accSensitivity",14) == 0)) { int grad = atoi(&(dataToStore[14])); if (grad > 0) dataFlag = 3000 + grad; else dataFlag = 6; } else { dataFlag = 6; opserr << "NodeRecorder::NodeRecorder - dataToStore " << dataToStore; opserr << "not recognized (disp, vel, accel, incrDisp, incrDeltaDisp)\n"; } } NodeRecorder::~NodeRecorder() { theOutputHandler->endTag(); // Data theOutputHandler->endTag(); // OpenSeesOutput if (theOutputHandler != 0) delete theOutputHandler; if (theDofs != 0) delete theDofs; if (theNodalTags != 0) delete theNodalTags; if (theNodes != 0) delete [] theNodes; } int NodeRecorder::record(int commitTag, double timeStamp) { if (theDomain == 0 || theNodalTags == 0 || theDofs == 0) { return 0; } if (theOutputHandler == 0) { opserr << "NodeRecorder::record() - no DataOutputHandler has been set\n"; return -1; } if (initializationDone != true) if (this->initialize() != 0) { opserr << "NodeRecorder::record() - failed in initialize()\n"; return -1; } int numDOF = theDofs->Size(); if (deltaT == 0.0 || timeStamp >= nextTimeStampToRecord) { if (deltaT != 0.0) nextTimeStampToRecord = timeStamp + deltaT; // // if need nodal reactions get the domain to calculate them // before we iterate over the nodes // if (dataFlag == 7) theDomain->calculateNodalReactions(false); else if (dataFlag == 8) theDomain->calculateNodalReactions(true); // // add time information if requested // int timeOffset = 0; if (echoTimeFlag == true) { timeOffset = 1; response(0) = timeStamp; } // // now we go get the responses from the nodes & place them in disp vector // for (int i=0; i<numValidNodes; i++) { int cnt = i*numDOF + timeOffset; Node *theNode = theNodes[i]; if (dataFlag == 0) { // AddingSensitivity:BEGIN /////////////////////////////////// if (sensitivity==0) { const Vector &theResponse = theNode->getTrialDisp(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else { response(cnt) = 0.0; } cnt++; } } else { for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); response(cnt) = theNode->getDispSensitivity(dof+1, sensitivity); cnt++; } } // AddingSensitivity:END ///////////////////////////////////// } else if (dataFlag == 1) { const Vector &theResponse = theNode->getTrialVel(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else response(cnt) = 0.0; cnt++; } } else if (dataFlag == 2) { const Vector &theResponse = theNode->getTrialAccel(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else response(cnt) = 0.0; cnt++; } } else if (dataFlag == 3) { const Vector &theResponse = theNode->getIncrDisp(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else response(cnt) = 0.0; cnt++; } } else if (dataFlag == 4) { const Vector &theResponse = theNode->getIncrDeltaDisp(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else response(cnt) = 0.0; cnt++; } } else if (dataFlag == 5) { const Vector &theResponse = theNode->getUnbalancedLoad(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else response(cnt) = 0.0; cnt++; } } else if (dataFlag == 6) { const Vector &theResponse = theNode->getUnbalancedLoadIncInertia(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else response(cnt) = 0.0; cnt++; } } else if (dataFlag == 7) { const Vector &theResponse = theNode->getReaction(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else response(cnt) = 0.0; cnt++; } } else if (dataFlag == 8) { const Vector &theResponse = theNode->getReaction(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (theResponse.Size() > dof) { response(cnt) = theResponse(dof); } else response(cnt) = 0.0; cnt++; } } else if (10 <= dataFlag && dataFlag < 1000) { int mode = dataFlag - 10; int column = mode - 1; const Matrix &theEigenvectors = theNode->getEigenvectors(); if (theEigenvectors.noCols() > column) { int noRows = theEigenvectors.noRows(); for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); if (noRows > dof) { response(cnt) = theEigenvectors(dof,column); } else response(cnt) = 0.0; cnt++; } } } else if (dataFlag >= 1000 && dataFlag < 2000) { int grad = dataFlag - 1000; for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); dof += 1; // Terje uses 1 through DOF for the dof indexing; the fool then subtracts 1 // his code!! response(cnt) = theNode->getDispSensitivity(dof, grad); cnt++; } } else if (dataFlag >= 2000 && dataFlag < 3000) { int grad = dataFlag - 2000; for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); dof += 1; // Terje uses 1 through DOF for the dof indexing; the fool then subtracts 1 // his code!! response(cnt) = theNode->getVelSensitivity(dof, grad); cnt++; } } else if (dataFlag >= 3000) { int grad = dataFlag - 3000; for (int j=0; j<numDOF; j++) { int dof = (*theDofs)(j); dof += 1; // Terje uses 1 through DOF for the dof indexing; the fool then subtracts 1 // his code!! response(cnt) = theNode->getAccSensitivity(dof, grad); cnt++; } } else { // unknown response for (int j=0; j<numDOF; j++) { response(cnt) = 0.0; } } } // insert the data into the database theOutputHandler->write(response); } return 0; } int NodeRecorder::setDomain(Domain &theDom) { theDomain = &theDom; return 0; } int NodeRecorder::sendSelf(int commitTag, Channel &theChannel) { if (theChannel.isDatastore() == 1) { opserr << "NodeRecorder::sendSelf() - does not send data to a datastore\n"; return -1; } static ID idData(6); idData.Zero(); if (theDofs != 0) idData(0) = theDofs->Size(); if (theNodalTags != 0) idData(1) = theNodalTags->Size(); if (theOutputHandler != 0) { idData(2) = theOutputHandler->getClassTag(); } if (echoTimeFlag == true) idData(3) = 1; else idData(3) = 0; idData(4) = dataFlag; idData(5) = sensitivity; if (theChannel.sendID(0, commitTag, idData) < 0) { opserr << "NodeRecorder::sendSelf() - failed to send idData\n"; return -1; } if (theDofs != 0) if (theChannel.sendID(0, commitTag, *theDofs) < 0) { opserr << "NodeRecorder::sendSelf() - failed to send dof id's\n"; return -1; } if (theNodalTags != 0) if (theChannel.sendID(0, commitTag, *theNodalTags) < 0) { opserr << "NodeRecorder::sendSelf() - failed to send nodal tags\n"; return -1; } static Vector data(2); data(0) = deltaT; data(1) = nextTimeStampToRecord; if (theChannel.sendVector(0, commitTag, data) < 0) { opserr << "NodeRecorder::sendSelf() - failed to send data\n"; return -1; } if (theOutputHandler->sendSelf(commitTag, theChannel) < 0) { opserr << "NodeRecorder::sendSelf() - failed to send the DataOutputHandler\n"; return -1; } return 0; } int NodeRecorder::recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker) { if (theChannel.isDatastore() == 1) { opserr << "NodeRecorder::sendSelf() - does not send data to a datastore\n"; return -1; } static ID idData(6); if (theChannel.recvID(0, commitTag, idData) < 0) { opserr << "NodeRecorder::recvSelf() - failed to send idData\n"; return -1; } int numDOFs = idData(0); int numNodes = idData(1); if (idData(3) == 1) echoTimeFlag = true; else echoTimeFlag = false; dataFlag = idData(4); sensitivity = idData(5); // // get the DOF ID data // if (theDofs == 0 || theDofs->Size() != numDOFs) { if (theDofs != 0) delete theDofs; if (numDOFs != 0) { theDofs = new ID(numDOFs); if (theDofs == 0 || theDofs->Size() != numDOFs) { opserr << "NodeRecorder::recvSelf() - out of memory\n"; return -1; } } } if (theDofs != 0) if (theChannel.recvID(0, commitTag, *theDofs) < 0) { opserr << "NodeRecorder::recvSelf() - failed to recv dof data\n"; return -1; } // // get the NODAL tag data // if (theNodalTags == 0 || theNodalTags->Size() != numNodes) { if (theNodalTags != 0) delete theNodalTags; if (numNodes != 0) { theNodalTags = new ID(numNodes); if (theNodalTags == 0 || theNodalTags->Size() != numNodes) { opserr << "NodeRecorder::recvSelf() - out of memory\n"; return -1; } } } if (theNodalTags != 0) if (theChannel.recvID(0, commitTag, *theNodalTags) < 0) { opserr << "NodeRecorder::recvSelf() - failed to recv dof data\n"; return -1; } static Vector data(2); data(0) = deltaT; data(1) = nextTimeStampToRecord; if (theChannel.recvVector(0, commitTag, data) < 0) { opserr << "NodeRecorder::sendSelf() - failed to receive data\n"; return -1; } if (theOutputHandler != 0) delete theOutputHandler; theOutputHandler = theBroker.getPtrNewStream(idData(2)); if (theOutputHandler == 0) { opserr << "NodeRecorder::sendSelf() - failed to get a data output handler\n"; return -1; } if (theOutputHandler->recvSelf(commitTag, theChannel, theBroker) < 0) { opserr << "NodeRecorder::sendSelf() - failed to send the DataOutputHandler\n"; return -1; } return 0; } int NodeRecorder::initialize(void) { if (theDofs == 0 || theNodalTags == 0 || theDomain == 0) { opserr << "NodeRecorder::initialize() - either nodes, dofs or domain has not been set\n"; return -1; } theOutputHandler->tag("OpenSeesOutput"); if (echoTimeFlag == true) { theOutputHandler->tag("TimeOutput"); theOutputHandler->tag("ResponseType", "time"); theOutputHandler->endTag(); } // // create & set nodal array pointer // if (theNodes != 0) delete [] theNodes; numValidNodes = 0; int i; int numNode = theNodalTags->Size(); for (i=0; i<numNode; i++) { int nodeTag = (*theNodalTags)(i); Node *theNode = theDomain->getNode(nodeTag); if (theNode != 0) { numValidNodes++; } } theNodes = new Node *[numValidNodes]; if (theNodes == 0) { opserr << "NodeRecorder::domainChanged - out of memory\n"; return -1; } int count = 0; for (i=0; i<numNode; i++) { int nodeTag = (*theNodalTags)(i); Node *theNode = theDomain->getNode(nodeTag); if (theNode != 0) { theNodes[count] = theNode; count++; } } // // resize the response vector // int timeOffset = 0; if (echoTimeFlag == true) timeOffset = 1; int numValidResponse = numValidNodes*theDofs->Size() + timeOffset; response.resize(numValidResponse); response.Zero(); // // need to create the data description, i.e. what each column of data is // char outputData[32]; char dataType[10]; if (dataFlag == 0) { strcpy(dataType,"D"); } else if (dataFlag == 1) { strcpy(dataType,"V"); } else if (dataFlag == 2) { strcpy(dataType,"A"); } else if (dataFlag == 3) { strcpy(dataType,"dD"); } else if (dataFlag == 4) { strcpy(dataType,"ddD"); } else if (dataFlag == 5) { strcpy(dataType,"U"); } else if (dataFlag == 6) { strcpy(dataType,"U"); } else if (dataFlag == 7) { strcpy(dataType,"R"); } else if (dataFlag == 8) { strcpy(dataType,"R"); } else if (dataFlag > 10) { sprintf(dataType,"E%d", dataFlag-10); } else strcpy(dataType,"Unknown"); /************************************************************ } else if ((strncmp(dataToStore, "sensitivity",11) == 0)) { int grad = atoi(&(dataToStore[11])); if (grad > 0) dataFlag = 1000 + grad; else dataFlag = 6; } else if ((strncmp(dataToStore, "velSensitivity",14) == 0)) { int grad = atoi(&(dataToStore[14])); if (grad > 0) dataFlag = 2000 + grad; else dataFlag = 6; } else if ((strncmp(dataToStore, "accSensitivity",14) == 0)) { int grad = atoi(&(dataToStore[14])); if (grad > 0) dataFlag = 3000 + grad; else dataFlag = 6; ***********************************************************/ char nodeCrdData[20]; sprintf(nodeCrdData,"coord"); for (i=0; i<numValidNodes; i++) { int nodeTag = theNodes[i]->getTag(); const Vector &nodeCrd = theNodes[i]->getCrds(); int numCoord = nodeCrd.Size(); theOutputHandler->tag("NodeOutput"); theOutputHandler->attr("nodeTag", nodeTag); for (int j=0; j<3; j++) { sprintf(nodeCrdData,"coord%d",j+1); if (j < numCoord) theOutputHandler->attr(nodeCrdData, nodeCrd(j)); else theOutputHandler->attr(nodeCrdData, 0.0); } for (int k=0; k<theDofs->Size(); k++) { sprintf(outputData, "%s%d", dataType, k+1); theOutputHandler->tag("ResponseType",outputData); } theOutputHandler->endTag(); } theOutputHandler->tag("Data"); initializationDone = true; return 0; } Generated on Mon Oct 23 15:05:25 2006 for OpenSees by  1.5.0

---

opensees-support @ berkeley.edu

Supported by the National Science Foundation

©2006, UC Regents