PEER NEES NEESit

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

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

ElasticSection2d.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.10 $ // $Date: 2006/09/06 20:17:34 $ // $Source: /usr/local/cvs/OpenSees/SRC/material/section/ElasticSection2d.cpp,v $ #include <ElasticSection2d.h> #include <Matrix.h> #include <Vector.h> #include <Channel.h> #include <FEM_ObjectBroker.h> #include <MatrixUtil.h> #include <stdlib.h> #include <Information.h> #include <Parameter.h> #include <classTags.h> Vector ElasticSection2d::s(2); Matrix ElasticSection2d::ks(2,2); ID ElasticSection2d::code(2); ElasticSection2d::ElasticSection2d(void) :SectionForceDeformation(0, SEC_TAG_Elastic2d), E(0.0), A(0.0), I(0.0), e(2), eCommit(2) { if (code(0) != SECTION_RESPONSE_P) { code(0) = SECTION_RESPONSE_P; // P is the first quantity code(1) = SECTION_RESPONSE_MZ; // Mz is the second } } ElasticSection2d::ElasticSection2d (int tag, double E_in, double A_in, double I_in) :SectionForceDeformation(tag, SEC_TAG_Elastic2d), E(E_in), A(A_in), I(I_in), e(2), eCommit(2) { if (E <= 0.0) { opserr << "ElasticSection2d::ElasticSection2d -- Input E <= 0.0 ... setting E to 1.0\n"; E = 1.0; } if (A <= 0.0) { opserr << "ElasticSection2d::ElasticSection2d -- Input A <= 0.0 ... setting A to 1.0\n"; A = 1.0; } if (I <= 0.0) { opserr << "ElasticSection2d::ElasticSection2d -- Input I <= 0.0 ... setting I to 1.0\n"; I = 1.0; } if (code(0) != SECTION_RESPONSE_P) { code(0) = SECTION_RESPONSE_P; // P is the first quantity code(1) = SECTION_RESPONSE_MZ; // Mz is the second } } ElasticSection2d::~ElasticSection2d(void) { } int ElasticSection2d::commitState(void) { eCommit = e; return 0; } int ElasticSection2d::revertToLastCommit(void) { e = eCommit; return 0; } int ElasticSection2d::revertToStart(void) { eCommit.Zero(); return 0; } int ElasticSection2d::setTrialSectionDeformation(const Vector &def) { e = def; return 0; } const Vector & ElasticSection2d::getSectionDeformation(void) { return e; } const Vector & ElasticSection2d::getStressResultant(void) { s(0) = E*A*e(0); s(1) = E*I*e(1); return s; } const Matrix & ElasticSection2d::getSectionTangent(void) { ks(0,0) = E*A; ks(1,1) = E*I; return ks; } const Matrix & ElasticSection2d::getInitialTangent(void) { ks(0,0) = E*A; ks(1,1) = E*I; return ks; } const Matrix & ElasticSection2d::getSectionFlexibility(void) { ks(0,0) = 1.0/(E*A); ks(1,1) = 1.0/(E*I); return ks; } const Matrix & ElasticSection2d::getInitialFlexibility(void) { ks(0,0) = 1.0/(E*A); ks(1,1) = 1.0/(E*I); return ks; } SectionForceDeformation* ElasticSection2d::getCopy(void) { // Make a copy of the hinge ElasticSection2d *theCopy = new ElasticSection2d (this->getTag(), E, A, I); theCopy->eCommit = eCommit; return theCopy; } const ID& ElasticSection2d::getType(void) { return code; } int ElasticSection2d::getOrder(void) const { return 2; } int ElasticSection2d::sendSelf(int commitTag, Channel &theChannel) { int res = 0; static Vector data(6); int dataTag = this->getDbTag(); data(0) = this->getTag(); data(1) = E; data(2) = A; data(3) = I; data(4) = eCommit(0); data(5) = eCommit(1); res += theChannel.sendVector(dataTag, commitTag, data); if (res<0) { opserr << "ElasticSection2d::sendSelf -- failed to send data\n"; return res; } return res; } int ElasticSection2d::recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker) { int res = 0; static Vector data(6); int dataTag = this->getDbTag(); res += theChannel.recvVector(dataTag, commitTag, data); if(res < 0) { opserr << "ElasticSection2d::recvSelf -- failed to receive data\n"; return res; } this->setTag((int)data(0)); E = data(1); A = data(2); I = data(3); eCommit(0) = data(4); eCommit(1) = data(5); return res; } void ElasticSection2d::Print(OPS_Stream &s, int flag) { s << "ElasticSection2d, tag: " << this->getTag() << endln; s << "\tE: " << E << endln; s << "\tA: " << A << endln; s << "\tI: " << I << endln; } int ElasticSection2d::setParameter(const char **argv, int argc, Parameter &param) { if (argc < 1) return -1; if (strcmp(argv[0],"E") == 0) return param.addObject(1, this); if (strcmp(argv[0],"A") == 0) return param.addObject(2, this); if (strcmp(argv[0],"I") == 0) return param.addObject(3, this); return -1; } int ElasticSection2d::updateParameter(int paramID, Information &info) { if (paramID == 1) E = info.theDouble; if (paramID == 2) A = info.theDouble; if (paramID == 3) I = info.theDouble; return 0; } int ElasticSection2d::activateParameter(int paramID) { parameterID = paramID; return 0; } const Vector& ElasticSection2d::getStressResultantSensitivity(int gradNumber, bool conditional) { s.Zero(); if (parameterID == 1) { // E s(0) = A*e(0); s(1) = I*e(1); } if (parameterID == 2) // A s(0) = E*e(0); if (parameterID == 3) // I s(1) = E*e(1); return s; } const Vector& ElasticSection2d::getSectionDeformationSensitivity(int gradNumber) { s.Zero(); return s; } const Matrix& ElasticSection2d::getInitialTangentSensitivity(int gradNumber) { ks.Zero(); return ks; } int ElasticSection2d::commitSensitivity(const Vector& sectionDeformationGradient, int gradNumber, int numGrads) { // Nothing to commit, path independent return 0; } Generated on Mon Oct 23 15:05:17 2006 for OpenSees by  1.5.0

---

opensees-support @ berkeley.edu

Supported by the National Science Foundation

©2006, UC Regents