PEER NEES NEESit

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

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

WideFlangeSectionIntegration.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: 2006/08/11 18:32:56 $ // $Source: /usr/local/cvs/OpenSees/SRC/material/section/integration/WideFlangeSectionIntegration.cpp,v $ #include <WideFlangeSectionIntegration.h> #include <Matrix.h> #include <Vector.h> #include <Channel.h> #include <FEM_ObjectBroker.h> #include <Information.h> #include <Parameter.h> WideFlangeSectionIntegration::WideFlangeSectionIntegration(double D, double TW, double BF, double TF, int NFDW, int NFTF): SectionIntegration(SECTION_INTEGRATION_TAG_WideFlange), d(D), tw(TW), bf(BF), tf(TF), Nfdw(NFDW), Nftf(NFTF), parameterID(0) { } WideFlangeSectionIntegration::WideFlangeSectionIntegration(): SectionIntegration(SECTION_INTEGRATION_TAG_WideFlange), d(0.0), tw(0.0), bf(0.0), tf(0.0), Nfdw(0), Nftf(0), parameterID(0) { } WideFlangeSectionIntegration::~WideFlangeSectionIntegration() { } int WideFlangeSectionIntegration::getNumFibers(void) { return Nfdw + 2*Nftf; } int WideFlangeSectionIntegration::arrangeFibers(UniaxialMaterial **theMaterials, UniaxialMaterial *theSteel) { int numFibers = this->getNumFibers(); for (int i = 0; i < numFibers; i++) theMaterials[i] = theSteel; return 0; } void WideFlangeSectionIntegration::getFiberLocations(int nFibers, double *xi) { double dw = d-2*tf; int loc; double yIncr = tf/Nftf; double yStart = 0.5 * (d-yIncr); for (loc = 0; loc < Nftf; loc++) { xi[loc] = yStart - yIncr*loc; xi[nFibers-loc-1] = -xi[loc]; } yIncr = dw/Nfdw; yStart = 0.5 * (dw-yIncr); for (int count = 0; loc < nFibers-Nftf; loc++, count++) { xi[loc] = yStart - yIncr*count; } return; } void WideFlangeSectionIntegration::getFiberWeights(int nFibers, double *wt) { double dw = d-2*tf; double a_f = bf*tf/Nftf; double a_w = dw*tw/Nfdw; int loc = 0; for (loc = 0; loc < Nftf; loc++) { wt[loc] = a_f; wt[nFibers-loc-1] = a_f; } for ( ; loc < nFibers-Nftf; loc++) { wt[loc] = a_w; } return; } SectionIntegration* WideFlangeSectionIntegration::getCopy(void) { return new WideFlangeSectionIntegration(d, tw, bf, tf, Nfdw, Nftf); } int WideFlangeSectionIntegration::setParameter(const char **argv, int argc, Parameter &param) { if (argc < 1) return -1; if (strcmp(argv[0],"d") == 0 || strcmp(argv[0],"db") == 0) return param.addObject(1, this); if (strcmp(argv[0],"tw") == 0) return param.addObject(2, this); if (strcmp(argv[0],"bf") == 0) return param.addObject(3, this); if (strcmp(argv[0],"tf") == 0) return param.addObject(4, this); return -1; } int WideFlangeSectionIntegration::updateParameter(int parameterID, Information &info) { switch (parameterID) { case 1: d = info.theDouble; return 0; case 2: tw = info.theDouble; return 0; case 3: bf = info.theDouble; return 0; case 4: tf = info.theDouble; return 0; default: return -1; } } int WideFlangeSectionIntegration::activateParameter(int paramID) { parameterID = paramID; return 0; } void WideFlangeSectionIntegration::getLocationsDeriv(int nFibers, double *dptsdh) { double dw = d-2*tf; double dddh = 0.0; double ddwdh = 0.0; //double dtwdh = 0.0; //double dbfdh = 0.0; double dtfdh = 0.0; if (parameterID == 1) { // d dddh = 1.0; ddwdh = 1.0; } //if (parameterID == 2) // tw // dtwdh = 1.0; //if (parameterID == 3) // bf // dbfdh = 1.0; if (parameterID == 4) { // tf dtfdh = 1.0; ddwdh = -2.0; } //double yIncr = tf/Nftf; //double yStart = 0.5*d - 0.5*yIncr; double dyIncrdh = dtfdh/Nftf; double dyStartdh = 0.5 * (dddh-dyIncrdh); int loc; for (loc = 0; loc < Nftf; loc++) { //xi[loc] = yStart - yIncr*loc; //xi[nFibers-loc-1] = -xi[loc]; dptsdh[loc] = dyStartdh - dyIncrdh*loc; dptsdh[nFibers-loc-1] = -dptsdh[loc]; } //yIncr = dw/Nfdw; //yStart = 0.5*dw - 0.5*yIncr; dyIncrdh = ddwdh/Nfdw; dyStartdh = 0.5 * (ddwdh-dyIncrdh); for (int count = 0; loc < nFibers-Nftf; loc++, count++) { //xi[loc] = yStart - yIncr*count; dptsdh[loc] = dyStartdh - dyIncrdh*count; } //for (int i = 0; i < nFibers; i++) // opserr << dptsdh[i] << ' '; //opserr << endln; return; } void WideFlangeSectionIntegration::getWeightsDeriv(int nFibers, double *dwtsdh) { double dw = d-2*tf; double ddwdh = 0.0; double dtwdh = 0.0; double dbfdh = 0.0; double dtfdh = 0.0; if (parameterID == 1) // d ddwdh = 1.0; if (parameterID == 2) // tw dtwdh = 1.0; if (parameterID == 3) // bf dbfdh = 1.0; if (parameterID == 4) { // tf dtfdh = 1.0; ddwdh = -2.0; } double dAfdh = (bf*dtfdh + dbfdh*tf) / Nftf; double dAwdh = (dw*dtwdh + ddwdh*tw) / Nfdw; int loc = 0; for (loc = 0; loc < Nftf; loc++) { dwtsdh[loc] = dAfdh; dwtsdh[nFibers-loc-1] = dAfdh; } for ( ; loc < nFibers-Nftf; loc++) dwtsdh[loc] = dAwdh; //for (int i = 0; i < nFibers; i++) // opserr << dwtsdh[i] << ' '; //opserr << endln; return; } void WideFlangeSectionIntegration::Print(OPS_Stream &s, int flag) { s << "WideFlange" << endln; s << " d = " << d; s << " tw = " << tw; s << " bf = " << bf; s << " tf = " << tf << endln; s << " Nfdw = " << Nfdw; s << " Nftf = " << Nftf << endln; return; } int WideFlangeSectionIntegration::sendSelf(int cTag, Channel &theChannel) { static Vector data(6); data(0) = d; data(1) = tw; data(2) = bf; data(3) = tf; data(4) = Nfdw; data(5) = Nftf; int dbTag = this->getDbTag(); if (theChannel.sendVector(dbTag, cTag, data) < 0) { opserr << "WideFlangeSectionIntegration::sendSelf() - failed to send Vector data\n"; return -1; } return 0; } int WideFlangeSectionIntegration::recvSelf(int cTag, Channel &theChannel, FEM_ObjectBroker &theBroker) { static Vector data(6); int dbTag = this->getDbTag(); if (theChannel.recvVector(dbTag, cTag, data) < 0) { opserr << "WideFlangeSectionIntegration::recvSelf() - failed to receive Vector data\n"; return -1; } d = data(0); tw = data(1); bf = data(2); tf = data(3); Nfdw = (int)data(4); Nftf = (int)data(5); return 0; } Generated on Mon Oct 23 15:05:18 2006 for OpenSees by  1.5.0

---

opensees-support @ berkeley.edu

Supported by the National Science Foundation

©2006, UC Regents