PEER NEES NEESit

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

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

BeamFiberMaterial.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.4 $ // $Date: 2003/02/14 23:01:24 $ // $Source: /usr/local/cvs/OpenSees/SRC/material/nD/BeamFiberMaterial.cpp,v $ // Written: MHS // Created: Aug 2001 // // Description: This file contains the class definition of BeamFiberMaterial. // The BeamFiberMaterial class is a wrapper class that performs static // condensation on a three-dimensional material model to give the 11, 12, and 13 // stress components which can then be integrated over an area to model a // shear flexible 3D beam. #include <BeamFiberMaterial.h> #include <Channel.h> #include <FEM_ObjectBroker.h> Vector BeamFiberMaterial::stress(3); Matrix BeamFiberMaterial::tangent(3,3); BeamFiberMaterial::BeamFiberMaterial(void) : NDMaterial(0, ND_TAG_BeamFiberMaterial), Tstrain22(0.0), Tstrain33(0.0), Tgamma23(0.0), Cstrain22(0.0), Cstrain33(0.0), Cgamma23(0.0), theMaterial(0), strain(3) { // Nothing to do } BeamFiberMaterial::BeamFiberMaterial(int tag, NDMaterial &theMat) : NDMaterial(tag, ND_TAG_BeamFiberMaterial), Tstrain22(0.0), Tstrain33(0.0), Tgamma23(0.0), Cstrain22(0.0), Cstrain33(0.0), Cgamma23(0.0), theMaterial(0), strain(3) { // Get a copy of the material theMaterial = theMat.getCopy("ThreeDimensional"); if (theMaterial == 0) { opserr << "BeamFiberMaterial::BeamFiberMaterial -- failed to get copy of material\n"; exit(-1); } } BeamFiberMaterial::~BeamFiberMaterial(void) { if (theMaterial != 0) delete theMaterial; } NDMaterial* BeamFiberMaterial::getCopy(void) { BeamFiberMaterial *theCopy = new BeamFiberMaterial(this->getTag(), *theMaterial); theCopy->Tstrain22 = this->Tstrain22; theCopy->Tstrain33 = this->Tstrain33; theCopy->Tgamma23 = this->Tgamma23; theCopy->Cstrain22 = this->Cstrain22; theCopy->Cstrain33 = this->Cstrain33; theCopy->Cgamma23 = this->Cgamma23; return theCopy; } NDMaterial* BeamFiberMaterial::getCopy(const char *type) { if (strcmp(type, "BeamFiber") == 0) return this->getCopy(); else return 0; } int BeamFiberMaterial::getOrder(void) const { return 3; } const char* BeamFiberMaterial::getType(void) const { return "BeamFiber"; } int BeamFiberMaterial::commitState(void) { Cstrain22 = Tstrain22; Cstrain33 = Tstrain33; Cgamma23 = Tgamma23; return theMaterial->commitState(); } int BeamFiberMaterial::revertToLastCommit(void) { Tstrain22 = Cstrain22; Tstrain33 = Cstrain33; Tgamma23 = Cgamma23; return theMaterial->revertToLastCommit(); } int BeamFiberMaterial::revertToStart() { this->Tstrain22 = 0.0; this->Tstrain33 = 0.0; this->Tgamma23 = 0.0; this->Cstrain22 = 0.0; this->Cstrain33 = 0.0; this->Cgamma23 = 0.0; return theMaterial->revertToStart(); } double BeamFiberMaterial::getRho(void) { return theMaterial->getRho(); } //receive the strain //NDmaterial strain order = 11, 22, 33, 12, 23, 31 //BeamFiberMaterial strain order = 11, 12, 31, 22, 33, 23 int BeamFiberMaterial::setTrialStrain(const Vector &strainFromElement) { static const double tolerance = 1.0e-08; this->strain(0) = strainFromElement(0); this->strain(1) = strainFromElement(1); this->strain(2) = strainFromElement(2); //newton loop to solve for out-of-plane strains double norm; static Vector condensedStress(3); static Vector strainIncrement(3); static Vector threeDstress(6); static Vector threeDstrain(6); static Matrix threeDtangent(6,6); static Vector threeDstressCopy(6); static Matrix threeDtangentCopy(6,6); static Matrix dd22(3,3); int i, j; int ii, jj; do { //set three dimensional strain threeDstrain(0) = this->strain(0); threeDstrain(1) = this->Tstrain22; threeDstrain(2) = this->Tstrain33; threeDstrain(3) = this->strain(1); threeDstrain(4) = this->Tgamma23; threeDstrain(5) = this->strain(2); if (theMaterial->setTrialStrain(threeDstrain) < 0) { opserr << "BeamFiberMaterial::setTrialStrain - setStrain failed in material with strain " << threeDstrain; return -1; } //three dimensional stress threeDstress = theMaterial->getStress(); //three dimensional tangent threeDtangent = theMaterial->getTangent(); //NDmaterial strain order = 11, 22, 33, 12, 23, 31 //BeamFiberMaterial strain order = 11, 12, 31, 22, 33, 23 //swap matrix indices to sort out-of-plane components for (i=0; i<6; i++) { ii = this->indexMap(i); threeDstressCopy(ii) = threeDstress(i); for (j=0; j<6; j++) { jj = this->indexMap(j); threeDtangentCopy(ii,jj) = threeDtangent(i,j); }//end for j }//end for i //out of plane stress and tangents for (i=0; i<3; i++) { condensedStress(i) = threeDstressCopy(i+3); for (j=0; j<3; j++) dd22(i,j) = threeDtangentCopy(i+3,j+3); }//end for i //set norm norm = condensedStress.Norm(); //condensation dd22.Solve(condensedStress, strainIncrement); //update out of plane strains this->Tstrain22 -= strainIncrement(0); this->Tstrain33 -= strainIncrement(1); this->Tgamma23 -= strainIncrement(2); } while (norm > tolerance); return 0; } const Vector& BeamFiberMaterial::getStrain(void) { return this->strain; } const Vector& BeamFiberMaterial::getStress() { const Vector &threeDstress = theMaterial->getStress(); static Vector threeDstressCopy(6); int i, ii; //swap matrix indices to sort out-of-plane components for (i=0; i<6; i++) { ii = this->indexMap(i); threeDstressCopy(ii) = threeDstress(i); } for (i=0; i<3; i++) this->stress(i) = threeDstressCopy(i); return this->stress; } const Matrix& BeamFiberMaterial::getTangent() { static Matrix dd11(3,3); static Matrix dd12(3,3); static Matrix dd21(3,3); static Matrix dd22(3,3); static Matrix dd22invdd21(3,3); static Matrix threeDtangentCopy(6,6); const Matrix &threeDtangent = theMaterial->getTangent(); //swap matrix indices to sort out-of-plane components int i, j , ii, jj; for (i=0; i<6; i++) { ii = this->indexMap(i); for (j=0; j<6; j++) { jj = this->indexMap(j); threeDtangentCopy(ii,jj) = threeDtangent(i,j); }//end for j }//end for i for (i=0; i<3; i++) { for (j=0; j<3; j++) { dd11(i,j) = threeDtangentCopy(i, j ); dd12(i,j) = threeDtangentCopy(i, j+3); dd21(i,j) = threeDtangentCopy(i+3,j ); dd22(i,j) = threeDtangentCopy(i+3,j+3); } } //int Solve(const Vector &V, Vector &res) const; //int Solve(const Matrix &M, Matrix &res) const; //condensation dd22.Solve(dd21, dd22invdd21); this->tangent = dd11; this->tangent -= (dd12*dd22invdd21); return this->tangent; } const Matrix& BeamFiberMaterial::getInitialTangent() { static Matrix dd11(3,3); static Matrix dd12(3,3); static Matrix dd21(3,3); static Matrix dd22(3,3); static Matrix dd22invdd21(3,3); static Matrix threeDtangentCopy(6,6); const Matrix &threeDtangent = theMaterial->getInitialTangent(); //swap matrix indices to sort out-of-plane components int i, j , ii, jj; for (i=0; i<6; i++) { ii = this->indexMap(i); for (j=0; j<6; j++) { jj = this->indexMap(j); threeDtangentCopy(ii,jj) = threeDtangent(i,j); }//end for j }//end for i for (i=0; i<3; i++) { for (j=0; j<3; j++) { dd11(i,j) = threeDtangentCopy(i, j ); dd12(i,j) = threeDtangentCopy(i, j+3); dd21(i,j) = threeDtangentCopy(i+3,j ); dd22(i,j) = threeDtangentCopy(i+3,j+3); } } //int Solve(const Vector &V, Vector &res) const; //int Solve(const Matrix &M, Matrix &res) const; //condensation dd22.Solve(dd21, dd22invdd21); this->tangent = dd11; this->tangent -= (dd12*dd22invdd21); return this->tangent; } //NDmaterial strain order = 11, 22, 33, 12, 23, 31 //BeamFiberMaterial strain order = 11, 12, 31, 22, 33, 23 int BeamFiberMaterial::indexMap(int i) { int ii; if (i == 3) ii = 1; else if (i == 5) ii = 2; else if (i == 1) ii = 3; else if (i == 2) ii = 4; else if (i == 4) ii = 5; else ii = i; return ii; } void BeamFiberMaterial::Print(OPS_Stream &s, int flag) { s << "BeamFiberMaterial, tag: " << this->getTag() << endln; s << "\tWrapped material: "<< theMaterial->getTag() << endln; theMaterial->Print(s, flag); } int BeamFiberMaterial::sendSelf(int commitTag, Channel &theChannel) { int res = 0; // put tag and assocaited materials class and database tags into an id and send it static ID idData(3); idData(0) = this->getTag(); idData(1) = theMaterial->getClassTag(); int matDbTag = theMaterial->getDbTag(); if (matDbTag == 0) { matDbTag = theChannel.getDbTag(); theMaterial->setDbTag(matDbTag); } idData(2) = matDbTag; res = theChannel.sendID(this->getDbTag(), commitTag, idData); if (res < 0) { opserr << "BeamFiberMaterial::sendSelf() - failed to send id data\n"; return res; } // put the strains in a vector and send it static Vector vecData(3); vecData(0) = Cstrain22; vecData(1) = Cstrain33; vecData(2) = Cgamma23; res = theChannel.sendVector(this->getDbTag(), commitTag, vecData); if (res < 0) { opserr << "BeamFiberMaterial::sendSelf() - failed to send vector data\n"; return res; } // now send the materials data res = theMaterial->sendSelf(commitTag, theChannel); if (res < 0) opserr << "BeamFiberMaterial::sendSelf() - failed to send vector material\n"; return res; } int BeamFiberMaterial::recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker) { int res = 0; // recv an id containg the tag and associated materials class and db tags static ID idData(3); res = theChannel.sendID(this->getDbTag(), commitTag, idData); if (res < 0) { opserr << "BeamFiberMaterial::sendSelf() - failed to send id data\n"; return res; } this->setTag(idData(0)); int matClassTag = idData(1); // if the associated material has not yet been created or is of the wrong type // create a new material for recvSelf later if (theMaterial == 0 || theMaterial->getClassTag() != matClassTag) { if (theMaterial != 0) delete theMaterial; theMaterial = theBroker.getNewNDMaterial(matClassTag); if (theMaterial == 0) { opserr << "BeamFiberMaterial::recvSelf() - failed to get a material of type: " << matClassTag << endln; return -1; } } theMaterial->setDbTag(idData(2)); // recv a vector containing strains and set the strains static Vector vecData(3); res = theChannel.recvVector(this->getDbTag(), commitTag, vecData); if (res < 0) { opserr << "BeamFiberMaterial::sendSelf() - failed to send vector data\n"; return res; } Cstrain22 = vecData(0); Cstrain33 = vecData(1); Cgamma23 = vecData(2); Tstrain22 = Cstrain22; Tstrain33 = Cstrain33; Tgamma23 = Cgamma23; // now receive the materials data res = theMaterial->recvSelf(commitTag, theChannel, theBroker); if (res < 0) opserr << "BeamFiberMaterial::sendSelf() - failed to send vector material\n"; return res; } Generated on Mon Oct 23 15:05:13 2006 for OpenSees by  1.5.0

---

opensees-support @ berkeley.edu

Supported by the National Science Foundation

©2006, UC Regents