PEER NEES NEESit

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

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

FiberSection2d.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/09/05 23:29:17 $ // $Source: /usr/local/cvs/OpenSees/SRC/material/section/FiberSection2d.cpp,v $ // Written: fmk // Created: 04/04 // // Description: This file contains the class implementation of FiberSection2d. #include <stdlib.h> #include <Channel.h> #include <Vector.h> #include <Matrix.h> #include <MatrixUtil.h> #include <Fiber.h> #include <classTags.h> #include <FiberSection2d.h> #include <ID.h> #include <FEM_ObjectBroker.h> #include <Information.h> #include <MaterialResponse.h> #include <UniaxialMaterial.h> #include <SectionIntegration.h> ID FiberSection2d::code(2); // constructors: FiberSection2d::FiberSection2d(int tag, int num, Fiber **fibers): SectionForceDeformation(tag, SEC_TAG_FiberSection2d), numFibers(num), theMaterials(0), matData(0), yBar(0.0), sectionIntegr(0), e(2), eCommit(2), s(0), ks(0), dedh(2) { if (numFibers != 0) { theMaterials = new UniaxialMaterial *[numFibers]; if (theMaterials == 0) { opserr << "FiberSection2d::FiberSection2d -- failed to allocate Material pointers"; exit(-1); } matData = new double [numFibers*2]; if (matData == 0) { opserr << "FiberSection2d::FiberSection2d -- failed to allocate double array for material data\n"; exit(-1); } double Qz = 0.0; double A = 0.0; for (int i = 0; i < numFibers; i++) { Fiber *theFiber = fibers[i]; double yLoc, zLoc, Area; theFiber->getFiberLocation(yLoc, zLoc); Area = theFiber->getArea(); A += Area; Qz += yLoc*Area; matData[i*2] = yLoc; matData[i*2+1] = Area; UniaxialMaterial *theMat = theFiber->getMaterial(); theMaterials[i] = theMat->getCopy(); if (theMaterials[i] == 0) { opserr << "FiberSection2d::FiberSection2d -- failed to get copy of a Material\n"; exit(-1); } } yBar = Qz/A; } s = new Vector(sData, 2); ks = new Matrix(kData, 2, 2); sData[0] = 0.0; sData[1] = 0.0; kData[0] = 0.0; kData[1] = 0.0; kData[2] = 0.0; kData[3] = 0.0; code(0) = SECTION_RESPONSE_P; code(1) = SECTION_RESPONSE_MZ; } FiberSection2d::FiberSection2d(int tag, int num, UniaxialMaterial **mats, SectionIntegration &si): SectionForceDeformation(tag, SEC_TAG_FiberSection2d), numFibers(num), theMaterials(0), matData(0), yBar(0.0), sectionIntegr(0), e(2), eCommit(2), s(0), ks(0), dedh(2) { if (numFibers != 0) { theMaterials = new UniaxialMaterial *[numFibers]; if (theMaterials == 0) { opserr << "FiberSection2d::FiberSection2d -- failed to allocate Material pointers"; exit(-1); } matData = new double [numFibers*2]; if (matData == 0) { opserr << "FiberSection2d::FiberSection2d -- failed to allocate double array for material data\n"; exit(-1); } } sectionIntegr = si.getCopy(); if (sectionIntegr == 0) { opserr << "Error: FiberSection2d::FiberSection2d: could not create copy of section integration object" << endln; exit(-1); } double fiberLocs[10000]; sectionIntegr->getFiberLocations(numFibers, fiberLocs); double fiberArea[10000]; sectionIntegr->getFiberWeights(numFibers, fiberArea); double Qz = 0.0; double A = 0.0; for (int i = 0; i < numFibers; i++) { A += fiberArea[i]; Qz += fiberLocs[i]*fiberArea[i]; theMaterials[i] = mats[i]->getCopy(); if (theMaterials[i] == 0) { opserr << "FiberSection2d::FiberSection2d -- failed to get copy of a Material\n"; exit(-1); } } yBar = Qz/A; s = new Vector(sData, 2); ks = new Matrix(kData, 2, 2); sData[0] = 0.0; sData[1] = 0.0; kData[0] = 0.0; kData[1] = 0.0; kData[2] = 0.0; kData[3] = 0.0; code(0) = SECTION_RESPONSE_P; code(1) = SECTION_RESPONSE_MZ; } // constructor for blank object that recvSelf needs to be invoked upon FiberSection2d::FiberSection2d(): SectionForceDeformation(0, SEC_TAG_FiberSection2d), numFibers(0), theMaterials(0), matData(0), yBar(0.0), sectionIntegr(0), e(2), eCommit(2), s(0), ks(0), dedh(2) { s = new Vector(sData, 2); ks = new Matrix(kData, 2, 2); sData[0] = 0.0; sData[1] = 0.0; kData[0] = 0.0; kData[1] = 0.0; kData[2] = 0.0; kData[3] = 0.0; code(0) = SECTION_RESPONSE_P; code(1) = SECTION_RESPONSE_MZ; } int FiberSection2d::addFiber(Fiber &newFiber) { // need to create larger arrays int newSize = numFibers+1; UniaxialMaterial **newArray = new UniaxialMaterial *[newSize]; double *newMatData = new double [2 * newSize]; if (newArray == 0 || newMatData == 0) { opserr <<"FiberSection2d::addFiber -- failed to allocate Fiber pointers\n"; return -1; } // copy the old pointers and data int i; for (i = 0; i < numFibers; i++) { newArray[i] = theMaterials[i]; newMatData[2*i] = matData[2*i]; newMatData[2*i+1] = matData[2*i+1]; } // set the new pointers and data double yLoc, zLoc, Area; newFiber.getFiberLocation(yLoc, zLoc); Area = newFiber.getArea(); newMatData[numFibers*2] = yLoc; newMatData[numFibers*2+1] = Area; UniaxialMaterial *theMat = newFiber.getMaterial(); newArray[numFibers] = theMat->getCopy(); if (newArray[numFibers] == 0) { opserr <<"FiberSection2d::addFiber -- failed to get copy of a Material\n"; delete [] newMatData; return -1; } numFibers++; if (theMaterials != 0) { delete [] theMaterials; delete [] matData; } theMaterials = newArray; matData = newMatData; double Qz = 0.0; double A = 0.0; // Recompute centroid for (i = 0; i < numFibers; i++) { yLoc = -matData[2*i]; Area = matData[2*i+1]; A += Area; Qz += yLoc*Area; } yBar = Qz/A; return 0; } // destructor: FiberSection2d::~FiberSection2d() { if (theMaterials != 0) { for (int i = 0; i < numFibers; i++) if (theMaterials[i] != 0) delete theMaterials[i]; delete [] theMaterials; } if (matData != 0) delete [] matData; if (s != 0) delete s; if (ks != 0) delete ks; if (sectionIntegr != 0) delete sectionIntegr; } int FiberSection2d::setTrialSectionDeformation (const Vector &deforms) { int res = 0; e = deforms; kData[0] = 0.0; kData[1] = 0.0; kData[2] = 0.0; kData[3] = 0.0; sData[0] = 0.0; sData[1] = 0.0; double d0 = deforms(0); double d1 = deforms(1); double fiberLocs[10000]; double fiberArea[10000]; if (sectionIntegr != 0) { sectionIntegr->getFiberLocations(numFibers, fiberLocs); sectionIntegr->getFiberWeights(numFibers, fiberArea); } else { for (int i = 0; i < numFibers; i++) { fiberLocs[i] = matData[2*i]; fiberArea[i] = matData[2*i+1]; } } for (int i = 0; i < numFibers; i++) { UniaxialMaterial *theMat = theMaterials[i]; double y = fiberLocs[i] - yBar; double A = fiberArea[i]; // determine material strain and set it double strain = d0 - y*d1; double tangent, stress; res = theMat->setTrial(strain, stress, tangent); double ks0 = tangent * A; double ks1 = ks0 * -y; kData[0] += ks0; kData[1] += ks1; kData[3] += ks1 * -y; double fs0 = stress * A; sData[0] += fs0; sData[1] += fs0 * -y; } kData[2] = kData[1]; return res; } const Vector& FiberSection2d::getSectionDeformation(void) { return e; } const Matrix& FiberSection2d::getInitialTangent(void) { static double kInitial[4]; static Matrix kInitialMatrix(kInitial, 2, 2); kInitial[0] = 0.0; kInitial[1] = 0.0; kInitial[2] = 0.0; kInitial[3] = 0.0; double fiberLocs[10000]; double fiberArea[10000]; if (sectionIntegr != 0) { sectionIntegr->getFiberLocations(numFibers, fiberLocs); sectionIntegr->getFiberWeights(numFibers, fiberArea); } else { for (int i = 0; i < numFibers; i++) { fiberLocs[i] = matData[2*i]; fiberArea[i] = matData[2*i+1]; } } for (int i = 0; i < numFibers; i++) { UniaxialMaterial *theMat = theMaterials[i]; double y = fiberLocs[i] - yBar; double A = fiberArea[i]; double tangent = theMat->getInitialTangent(); double ks0 = tangent * A; double ks1 = ks0 * -y; kInitial[0] += ks0; kInitial[1] += ks1; kInitial[3] += ks1 * -y; } kInitial[2] = kInitial[1]; return kInitialMatrix; } const Matrix& FiberSection2d::getSectionTangent(void) { return *ks; } const Vector& FiberSection2d::getStressResultant(void) { return *s; } SectionForceDeformation* FiberSection2d::getCopy(void) { FiberSection2d *theCopy = new FiberSection2d (); theCopy->setTag(this->getTag()); theCopy->numFibers = numFibers; if (numFibers != 0) { theCopy->theMaterials = new UniaxialMaterial *[numFibers]; if (theCopy->theMaterials == 0) { opserr <<"FiberSection2d::getCopy -- failed to allocate Material pointers\n"; exit(-1); } theCopy->matData = new double [numFibers*2]; if (theCopy->matData == 0) { opserr << "FiberSection2d::getCopy -- failed to allocate double array for material data\n"; exit(-1); } for (int i = 0; i < numFibers; i++) { theCopy->matData[i*2] = matData[i*2]; theCopy->matData[i*2+1] = matData[i*2+1]; theCopy->theMaterials[i] = theMaterials[i]->getCopy(); if (theCopy->theMaterials[i] == 0) { opserr <<"FiberSection2d::getCopy -- failed to get copy of a Material"; exit(-1); } } } theCopy->eCommit = eCommit; theCopy->e = e; theCopy->yBar = yBar; theCopy->kData[0] = kData[0]; theCopy->kData[1] = kData[1]; theCopy->kData[2] = kData[2]; theCopy->kData[3] = kData[3]; theCopy->sData[0] = sData[0]; theCopy->sData[1] = sData[1]; if (sectionIntegr != 0) theCopy->sectionIntegr = sectionIntegr->getCopy(); else theCopy->sectionIntegr = 0; return theCopy; } const ID& FiberSection2d::getType () { return code; } int FiberSection2d::getOrder () const { return 2; } int FiberSection2d::commitState(void) { int err = 0; for (int i = 0; i < numFibers; i++) err += theMaterials[i]->commitState(); eCommit = e; return err; } int FiberSection2d::revertToLastCommit(void) { int err = 0; // Last committed section deformations e = eCommit; kData[0] = 0.0; kData[1] = 0.0; kData[2] = 0.0; kData[3] = 0.0; sData[0] = 0.0; sData[1] = 0.0; double fiberLocs[10000]; double fiberArea[10000]; if (sectionIntegr != 0) { sectionIntegr->getFiberLocations(numFibers, fiberLocs); sectionIntegr->getFiberWeights(numFibers, fiberArea); } else { for (int i = 0; i < numFibers; i++) { fiberLocs[i] = matData[2*i]; fiberArea[i] = matData[2*i+1]; } } for (int i = 0; i < numFibers; i++) { UniaxialMaterial *theMat = theMaterials[i]; double y = fiberLocs[i] - yBar; double A = fiberArea[i]; // invoke revertToLast on the material err += theMat->revertToLastCommit(); // get material stress & tangent for this strain and determine ks and fs double tangent = theMat->getTangent(); double stress = theMat->getStress(); double ks0 = tangent * A; double ks1 = ks0 * -y; kData[0] += ks0; kData[1] += ks1; kData[3] += ks1 * -y; double fs0 = stress * A; sData[0] = fs0; sData[1] = fs0 * -y; } kData[2] = kData[1]; return err; } int FiberSection2d::revertToStart(void) { // revert the fibers to start int err = 0; kData[0] = 0.0; kData[1] = 0.0; kData[2] = 0.0; kData[3] = 0.0; sData[0] = 0.0; sData[1] = 0.0; double fiberLocs[10000]; double fiberArea[10000]; if (sectionIntegr != 0) { sectionIntegr->getFiberLocations(numFibers, fiberLocs); sectionIntegr->getFiberWeights(numFibers, fiberArea); } else { for (int i = 0; i < numFibers; i++) { fiberLocs[i] = matData[2*i]; fiberArea[i] = matData[2*i+1]; } } for (int i = 0; i < numFibers; i++) { UniaxialMaterial *theMat = theMaterials[i]; double y = fiberLocs[i] - yBar; double A = fiberArea[i]; // invoke revertToLast on the material err += theMat->revertToStart(); // get material stress & tangent for this strain and determine ks and fs double tangent = theMat->getTangent(); double stress = theMat->getStress(); double ks0 = tangent * A; double ks1 = ks0 * -y; kData[0] += ks0; kData[1] += ks1; kData[3] += ks1 * -y; double fs0 = stress * A; sData[0] = fs0; sData[1] = fs0 * -y; } kData[2] = kData[1]; return err; } int FiberSection2d::sendSelf(int commitTag, Channel &theChannel) { int res = 0; // create an id to send objects tag and numFibers, // size 3 so no conflict with matData below if just 1 fiber static ID data(3); data(0) = this->getTag(); data(1) = numFibers; int dbTag = this->getDbTag(); res += theChannel.sendID(dbTag, commitTag, data); if (res < 0) { opserr << "FiberSection2d::sendSelf - failed to send ID data\n"; return res; } if (numFibers != 0) { // create an id containingg classTag and dbTag for each material & send it ID materialData(2*numFibers); for (int i=0; i<numFibers; i++) { UniaxialMaterial *theMat = theMaterials[i]; materialData(2*i) = theMat->getClassTag(); int matDbTag = theMat->getDbTag(); if (matDbTag == 0) { matDbTag = theChannel.getDbTag(); if (matDbTag != 0) theMat->setDbTag(matDbTag); } materialData(2*i+1) = matDbTag; } res += theChannel.sendID(dbTag, commitTag, materialData); if (res < 0) { opserr << "FiberSection2d::sendSelf - failed to send material data\n"; return res; } // send the fiber data, i.e. area and loc Vector fiberData(matData, 2*numFibers); res += theChannel.sendVector(dbTag, commitTag, fiberData); if (res < 0) { opserr << "FiberSection2d::sendSelf - failed to send material data\n"; return res; } // now invoke send(0 on all the materials for (int j=0; j<numFibers; j++) theMaterials[j]->sendSelf(commitTag, theChannel); } return res; } int FiberSection2d::recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker) { int res = 0; static ID data(3); int dbTag = this->getDbTag(); res += theChannel.recvID(dbTag, commitTag, data); if (res < 0) { opserr << "FiberSection2d::recvSelf - failed to recv ID data\n"; return res; } this->setTag(data(0)); // recv data about materials objects, classTag and dbTag if (data(1) != 0) { ID materialData(2*data(1)); res += theChannel.recvID(dbTag, commitTag, materialData); if (res < 0) { opserr << "FiberSection2d::recvSelf - failed to recv material data\n"; return res; } // if current arrays not of correct size, release old and resize if (theMaterials == 0 || numFibers != data(1)) { // delete old stuff if outa date if (theMaterials != 0) { for (int i=0; i<numFibers; i++) delete theMaterials[i]; delete [] theMaterials; if (matData != 0) delete [] matData; matData = 0; theMaterials = 0; } // create memory to hold material pointers and fiber data numFibers = data(1); if (numFibers != 0) { theMaterials = new UniaxialMaterial *[numFibers]; if (theMaterials == 0) { opserr <<"FiberSection2d::recvSelf -- failed to allocate Material pointers\n"; exit(-1); } for (int j=0; j<numFibers; j++) theMaterials[j] = 0; matData = new double [numFibers*2]; if (matData == 0) { opserr <<"FiberSection2d::recvSelf -- failed to allocate double array for material data\n"; exit(-1); } } } Vector fiberData(matData, 2*numFibers); res += theChannel.recvVector(dbTag, commitTag, fiberData); if (res < 0) { opserr << "FiberSection2d::recvSelf - failed to recv material data\n"; return res; } int i; for (i=0; i<numFibers; i++) { int classTag = materialData(2*i); int dbTag = materialData(2*i+1); // if material pointed to is blank or not of corrcet type, // release old and create a new one if (theMaterials[i] == 0) theMaterials[i] = theBroker.getNewUniaxialMaterial(classTag); else if (theMaterials[i]->getClassTag() != classTag) { delete theMaterials[i]; theMaterials[i] = theBroker.getNewUniaxialMaterial(classTag); } if (theMaterials[i] == 0) { opserr <<"FiberSection2d::recvSelf -- failed to allocate double array for material data\n"; exit(-1); } theMaterials[i]->setDbTag(dbTag); res += theMaterials[i]->recvSelf(commitTag, theChannel, theBroker); } double Qz = 0.0; double A = 0.0; double yLoc, Area; // Recompute centroid for (i = 0; i < numFibers; i++) { yLoc = matData[2*i]; Area = matData[2*i+1]; A += Area; Qz += yLoc*Area; } yBar = Qz/A; } return res; } void FiberSection2d::Print(OPS_Stream &s, int flag) { s << "\nFiberSection2d, tag: " << this->getTag() << endln; s << "\tSection code: " << code; s << "\tNumber of Fibers: " << numFibers << endln; s << "\tCentroid: " << yBar << endln; if (flag == 1) { for (int i = 0; i < numFibers; i++) { s << "\nLocation (y) = (" << matData[2*i] << ")"; s << "\nArea = " << matData[2*i+1] << endln; theMaterials[i]->Print(s, flag); } } } Response* FiberSection2d::setResponse(const char **argv, int argc, Information &sectInfo, OPS_Stream &output) { const ID &type = this->getType(); int typeSize = this->getOrder(); Response *theResponse =0; output.tag("SectionOutput"); output.attr("secType", this->getClassType()); output.attr("secTag", this->getTag()); // deformations if (strcmp(argv[0],"deformations") == 0 || strcmp(argv[0],"deformation") == 0) { for (int i=0; i<typeSize; i++) { int code = type(i); switch (code){ case SECTION_RESPONSE_MZ: output.tag("ResponseType","kappaZ"); break; case SECTION_RESPONSE_P: output.tag("ResponseType","eps"); break; case SECTION_RESPONSE_VY: output.tag("ResponseType","gammaY"); break; case SECTION_RESPONSE_MY: output.tag("ResponseType","kappaY"); break; case SECTION_RESPONSE_VZ: output.tag("ResponseType","gammaZ"); break; case SECTION_RESPONSE_T: output.tag("ResponseType","theta"); break; default: output.tag("ResponseType","Unknown"); } } theResponse = new MaterialResponse(this, 1, this->getSectionDeformation()); // forces } else if (strcmp(argv[0],"forces") == 0 || strcmp(argv[0],"force") == 0) { for (int i=0; i<typeSize; i++) { int code = type(i); switch (code){ case SECTION_RESPONSE_MZ: output.tag("ResponseType","Mz"); break; case SECTION_RESPONSE_P: output.tag("ResponseType","P"); break; case SECTION_RESPONSE_VY: output.tag("ResponseType","Vy"); break; case SECTION_RESPONSE_MY: output.tag("ResponseType","My"); break; case SECTION_RESPONSE_VZ: output.tag("ResponseType","Vz"); break; case SECTION_RESPONSE_T: output.tag("ResponseType","T"); break; default: output.tag("ResponseType","Unknown"); } } theResponse = new MaterialResponse(this, 2, this->getStressResultant()); // force and deformation } else if (strcmp(argv[0],"forceAndDeformation") == 0) { for (int j=0; j<typeSize; j++) { int code = type(j); switch (code){ case SECTION_RESPONSE_MZ: output.tag("ResponseType","kappaZ"); break; case SECTION_RESPONSE_P: output.tag("ResponseType","eps"); break; case SECTION_RESPONSE_VY: output.tag("ResponseType","gammaY"); break; case SECTION_RESPONSE_MY: output.tag("ResponseType","kappaY"); break; case SECTION_RESPONSE_VZ: output.tag("ResponseType","gammaZ"); break; case SECTION_RESPONSE_T: output.tag("ResponseType","theta"); break; default: output.tag("ResponseType","Unknown"); } } for (int i=0; i<typeSize; i++) { int code = type(i); switch (code){ case SECTION_RESPONSE_MZ: output.tag("ResponseType","Mz"); break; case SECTION_RESPONSE_P: output.tag("ResponseType","P"); break; case SECTION_RESPONSE_VY: output.tag("ResponseType","Vy"); break; case SECTION_RESPONSE_MY: output.tag("ResponseType","My"); break; case SECTION_RESPONSE_VZ: output.tag("ResponseType","Vz"); break; case SECTION_RESPONSE_T: output.tag("ResponseType","T"); break; default: output.tag("ResponseType","Unknown"); } } theResponse = new MaterialResponse(this, 4, Vector(2*this->getOrder())); } else { if (argc > 2 || strcmp(argv[0],"fiber") == 0) { int key = numFibers; int passarg = 2; if (argc <= 3) { // fiber number was input directly key = atoi(argv[1]); } else if (argc > 4) { // find fiber closest to coord. with mat tag int matTag = atoi(argv[3]); double yCoord = atof(argv[1]); double closestDist; double ySearch, dy; double distance; int j; // Find first fiber with specified material tag for (j = 0; j < numFibers; j++) { if (matTag == theMaterials[j]->getTag()) { ySearch = -matData[2*j]; dy = ySearch-yCoord; closestDist = fabs(dy); key = j; break; } } // Search the remaining fibers for ( ; j < numFibers; j++) { if (matTag == theMaterials[j]->getTag()) { ySearch = -matData[2*j]; dy = ySearch-yCoord; distance = fabs(dy); if (distance < closestDist) { closestDist = distance; key = j; } } } passarg = 4; } else { // fiber near-to coordinate specified double yCoord = atof(argv[1]); double closestDist; double ySearch, dy; double distance; ySearch = -matData[0]; dy = ySearch-yCoord; closestDist = fabs(dy); key = 0; for (int j = 1; j < numFibers; j++) { ySearch = -matData[2*j]; dy = ySearch-yCoord; distance = fabs(dy); if (distance < closestDist) { closestDist = distance; key = j; } } passarg = 3; } if (key < numFibers && key >= 0) { output.tag("FiberOutput"); output.attr("yLoc",-matData[2*key]); output.attr("zLoc",0.0); output.attr("area",matData[2*key+1]); theResponse = theMaterials[key]->setResponse(&argv[passarg], argc-passarg, sectInfo, output); output.endTag(); } } } output.endTag(); return theResponse; } int FiberSection2d::getResponse(int responseID, Information &sectInfo) { // Just call the base class method ... don't need to define // this function, but keeping it here just for clarity return SectionForceDeformation::getResponse(responseID, sectInfo); } // AddingSensitivity:BEGIN //////////////////////////////////// int FiberSection2d::setParameter(const char **argv, int argc, Parameter &param) { if (argc < 1) return -1; // Check if the parameter belongs to the material (only option for now) if (strstr(argv[0],"material") != 0) { if (argc < 3) return 0; // Get the tag of the material int materialTag = atoi(argv[1]); // Loop over fibers to find the right material int ok = 0; for (int i = 0; i < numFibers; i++) if (materialTag == theMaterials[i]->getTag()) ok += theMaterials[i]->setParameter(&argv[2], argc-2, param); return ok; } else return sectionIntegr->setParameter(argv, argc, param); } const Vector & FiberSection2d::getSectionDeformationSensitivity(int gradNumber) { static Vector dummy(2); return dummy; } const Vector & FiberSection2d::getStressResultantSensitivity(int gradNumber, bool conditional) { static Vector ds(2); ds.Zero(); double y, A, stressGradient, stress, tangent, sig_dAdh; double fiberLocs[10000]; double fiberArea[10000]; if (sectionIntegr != 0) { sectionIntegr->getFiberLocations(numFibers, fiberLocs); sectionIntegr->getFiberWeights(numFibers, fiberArea); } else { for (int i = 0; i < numFibers; i++) { fiberLocs[i] = matData[2*i]; fiberArea[i] = matData[2*i+1]; } } double locsDeriv[10000]; double areaDeriv[10000]; if (sectionIntegr != 0) { sectionIntegr->getLocationsDeriv(numFibers, locsDeriv); sectionIntegr->getWeightsDeriv(numFibers, areaDeriv); } else { for (int i = 0; i < numFibers; i++) { locsDeriv[i] = 0.0; areaDeriv[i] = 0.0; } } for (int i = 0; i < numFibers; i++) { y = fiberLocs[i] - yBar; A = fiberArea[i]; stressGradient = theMaterials[i]->getStressSensitivity(gradNumber,true); stressGradient = stressGradient * A; ds(0) += stressGradient; ds(1) += stressGradient * -y; stress = theMaterials[i]->getStress(); sig_dAdh = stress*areaDeriv[i]; ds(0) += sig_dAdh; ds(1) += sig_dAdh * -y; //ds(0) += 0.0; ds(1) += (stress*A) * -locsDeriv[i]; tangent = theMaterials[i]->getTangent(); tangent = tangent * A * e(1); ds(0) += -locsDeriv[i]*tangent; ds(1) += fiberLocs[i]*locsDeriv[i]*tangent; //opserr << locsDeriv[i] << ' ' << areaDeriv[i] << endln; } return ds; } const Matrix & FiberSection2d::getInitialTangentSensitivity(int gradNumber) { static Matrix dksdh(2,2); dksdh.Zero(); double y, A, dydh, dAdh, tangent, dtangentdh; double fiberLocs[10000]; double fiberArea[10000]; if (sectionIntegr != 0) { sectionIntegr->getFiberLocations(numFibers, fiberLocs); sectionIntegr->getFiberWeights(numFibers, fiberArea); } else { for (int i = 0; i < numFibers; i++) { fiberLocs[i] = matData[2*i]; fiberArea[i] = matData[2*i+1]; } } double locsDeriv[10000]; double areaDeriv[10000]; if (sectionIntegr != 0) { sectionIntegr->getLocationsDeriv(numFibers, locsDeriv); sectionIntegr->getWeightsDeriv(numFibers, areaDeriv); } else { for (int i = 0; i < numFibers; i++) { locsDeriv[i] = 0.0; areaDeriv[i] = 0.0; } } for (int i = 0; i < numFibers; i++) { y = fiberLocs[i] - yBar; A = fiberArea[i]; dydh = locsDeriv[i]; dAdh = areaDeriv[i]; tangent = theMaterials[i]->getInitialTangent(); dtangentdh = theMaterials[i]->getInitialTangentSensitivity(gradNumber); dksdh(0,0) += dtangentdh*A + tangent*dAdh; dksdh(0,1) += -y*(dtangentdh*A+tangent*dAdh) - dydh*(tangent*A); dksdh(1,1) += 2*(y*dydh*tangent*A) + y*y*(dtangentdh*A+tangent*dAdh); } dksdh(1,0) = dksdh(0,1); return dksdh; } int FiberSection2d::commitSensitivity(const Vector& defSens, int gradNumber, int numGrads) { double d0 = defSens(0); double d1 = defSens(1); dedh = defSens; double fiberLocs[10000]; if (sectionIntegr != 0) sectionIntegr->getFiberLocations(numFibers, fiberLocs); else { for (int i = 0; i < numFibers; i++) fiberLocs[i] = matData[2*i]; } double locsDeriv[10000]; double areaDeriv[10000]; if (sectionIntegr != 0) { sectionIntegr->getLocationsDeriv(numFibers, locsDeriv); sectionIntegr->getWeightsDeriv(numFibers, areaDeriv); } else { for (int i = 0; i < numFibers; i++) { locsDeriv[i] = 0.0; areaDeriv[i] = 0.0; } } double y; double kappa = e(1); for (int i = 0; i < numFibers; i++) { UniaxialMaterial *theMat = theMaterials[i]; y = fiberLocs[i] - yBar; // determine material strain and set it double strainSens = d0 - y*d1 - locsDeriv[i]*kappa; theMat->commitSensitivity(strainSens,gradNumber,numGrads); } return 0; } // AddingSensitivity:END /////////////////////////////////// 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