WebSVN - OpenSees - Diff - Rev 1231 and 1271 - /trunk/SRC/element/truss/Truss.cpp

 **   Gregory L. Fenves (fenves@ce.berkeley.edu)                       **
 **   Filip C. Filippou (filippou@ce.berkeley.edu)                     **
 **                                                                    **
 ** ****************************************************************** */
-// $Revision: 1.14 $
-// $Date: 2002-12-16 21:10:08 $
+// $Revision: 1.15 $
+// $Date: 2003-02-14 23:01:18 $
 // $Source: /usr/local/cvs/OpenSees/SRC/element/truss/Truss.cpp,v $
 // File: ~/element/truss/Truss.C
 //
 #include <Node.h>
 #include <Channel.h>
 #include <FEM_ObjectBroker.h>
 #include <UniaxialMaterial.h>
 #include <Renderer.h>
-#include <G3Globals.h>
 #include <math.h>
 #include <stdlib.h>
 #include <string.h>
   theMatrix(0), theVector(0), t(0),
   L(0.0), A(a), M(rho)
+{
     // get a copy of the material and check we obtained a valid copy
     theMaterial = theMat.getCopy();
-    if (theMaterial == 0)
-      g3ErrorHandler->fatal("FATAL Truss::Truss - %d %s %d\n", tag,
-                            "failed to get a copy of material with tag ",
-                            theMat.getTag());
+    if (theMaterial == 0) {
+      opserr << "FATAL Truss::Truss - " << tag <<
+        "failed to get a copy of material with tag " << theMat.getTag() << endln;
+      exit(-1);
+    }
     // ensure the connectedExternalNode ID is of correct size & set values
-    if (connectedExternalNodes.Size() != 2)
-      g3ErrorHandler->fatal("FATAL Truss::Truss - %d %s\n", tag,
-                            "failed to create an ID of size 2");
+    if (connectedExternalNodes.Size() != 2) {
+      opserr << "FATAL Truss::Truss - " <<  tag << "failed to create an ID of size 2\n";
+      exit(-1);
+    }
     connectedExternalNodes(0) = Nd1;
     connectedExternalNodes(1) = Nd2;
     // set node pointers to NULL
     for (int i=0; i<2; i++)
  dimension(0), numDOF(0),
  theMatrix(0), theVector(0), t(0),
   L(0.0), A(0.0), M(0.0)
+{
     // ensure the connectedExternalNode ID is of correct size
-    if (connectedExternalNodes.Size() != 2)
-      g3ErrorHandler->fatal("FATAL Truss::Truss - %s\n",
-                            "failed to create an ID of size 2");
+  if (connectedExternalNodes.Size() != 2) {
+      opserr << "FATAL Truss::Truss - failed to create an ID of size 2\n";
+      exit(-1);
+  }
-    for (int i=0; i<2; i++)
-      theNodes[i] = 0;
+  for (int i=0; i<2; i++)
+    theNodes[i] = 0;
 // AddingSensitivity:BEGIN /////////////////////////////////////
         gradientIdentifier = 0;
 // AddingSensitivity:END //////////////////////////////////////
+}
     theNodes[1] = theDomain->getNode(Nd2);
     // if can't find both - send a warning message
     if ((theNodes[0] == 0) || (theNodes[1] == 0)) {
       if (theNodes[0] == 0)
-        g3ErrorHandler->warning("Truss::setDomain() - truss %d node %d %s\n",
-                                this->getTag(), Nd1,
-                                "does not exist in the model");
+        opserr <<"Truss::setDomain() - truss" << this->getTag() << " node " << Nd1 <<
+          "does not exist in the model\n";
       else
-        g3ErrorHandler->warning("Truss::setDomain() - truss %d node %d %s\n",
-                                this->getTag(), Nd2,
-                                "does not exist in the model");
+        opserr <<"Truss::setDomain() - truss" << this->getTag() << " node " << Nd2 <<
+          "does not exist in the model\n";
       // fill this in so don't segment fault later
       numDOF = 2;
       theMatrix = &trussM2;
       theVector = &trussV2;
     int dofNd1 = theNodes[0]->getNumberDOF();
     int dofNd2 = theNodes[1]->getNumberDOF();
     // if differing dof at the ends - print a warning message
     if (dofNd1 != dofNd2) {
-      g3ErrorHandler->warning("WARNING Truss::setDomain(): nodes %d and %d %s %d\n",Nd1, Nd2,
-                              "have differing dof at ends for truss",this->getTag());
+      opserr <<"WARNING Truss::setDomain(): nodes " << Nd1 << " and " << Nd2 <<
+        "have differing dof at ends for truss " << this->getTag() << endln;
       // fill this in so don't segment fault later
       numDOF = 2;
       theMatrix = &trussM2;
       theVector = &trussV2;
         numDOF = 12;
         theMatrix = &trussM12;
         theVector = &trussV12;
+    }
     else {
-      g3ErrorHandler->warning("WARNING Truss::setDomain cannot handle %d dofs at nodes in %d d problem\n",
-                              dimension,dofNd1);
+      opserr <<"WARNING Truss::setDomain cannot handle " << dimension << " dofs at nodes in " <<
+        dofNd1  << " problem\n";
       numDOF = 2;
       theMatrix = &trussM2;
       theVector = &trussV2;
       return;
     // create a transformation matrix for the element
     t = new Matrix(1,numDOF);
     theLoad = new Vector(numDOF);
     if (t == 0 || (t->noCols() != numDOF)) {
-      g3ErrorHandler->fatal("Truss::setDomain - truss %d %s %d\n",
-                            this->getTag(),
-                            "out of memory creating T matrix of size 1 x",
-                            numDOF);
+      opserr << "Truss::setDomain - truss " << this->getTag() <<
+        "out of memory creating T matrix of size 1 x" << numDOF << endln;
+      exit(-1);
       return;
+    }
     if (theLoad == 0) {
-        g3ErrorHandler->fatal("Truss::setDomain - truss %d %s %d\n",
-                              this->getTag(),
-                              "out of memory creating vector of size",
-                              numDOF);
+      opserr << "Truss::setDomain - truss " << this->getTag() <<
+        "out of memory creating vector of size" << numDOF << endln;
+      exit(-1);
       return;
+    }
     // now determine the length, cosines and fill in the transformation
     // NOTE t = -t(every one else uses for residual calc)
         trans(0,1) = 1;
         double dx = end2Crd(0)-end1Crd(0);
         L = sqrt(dx*dx);
         if (L == 0.0) {
-          g3ErrorHandler->warning("WARNING Truss::setDomain() - truss %d has zero length\n",
-                                  this->getTag());
+          opserr <<"WARNING Truss::setDomain() - truss " << this->getTag() << " has zero length\n";
           return;
+        }
     } else if (dimension == 2) {
         double dx = end2Crd(0)-end1Crd(0);
         double dy = end2Crd(1)-end1Crd(1);
         L = sqrt(dx*dx + dy*dy);
         if (L == 0.0) {
-          g3ErrorHandler->warning("WARNING Truss::setDomain() - truss %d has zero length\n",
-                                  this->getTag());
+          opserr <<"WARNING Truss::setDomain() - truss " << this->getTag() << " has zero length\n";
           return;
+        }
         double cs = dx/L;
         double sn = dy/L;
         double dz = end2Crd(2)-end1Crd(2);
         L = sqrt(dx*dx + dy*dy + dz*dz);
         if (L == 0.0) {
-          g3ErrorHandler->warning("WARNING Truss::setDomain() - truss %d has zero length\n",
-                                  this->getTag());
+          opserr <<"WARNING Truss::setDomain() - truss " << this->getTag() << " has zero length\n";
           return;
+        }
         double cx = dx/L;
         double cy = dy/L;
 Truss::commitState()
+{
   int retVal = 0;
   // call element commitState to do any base class stuff
   if ((retVal = this->Element::commitState()) != 0) {
-    cerr << "Truss::commitState () - failed in base class";
+    opserr << "Truss::commitState () - failed in base class";
+  }
   retVal = theMaterial->commitState();
   return retVal;
+}
 int
 Truss::addLoad(ElementalLoad *theLoad, double loadFactor)
+{
-  g3ErrorHandler->warning("Truss::addLoad - load type unknown for truss with tag: %d\n",
-                          this->getTag());
+  opserr <<"Truss::addLoad - load type unknown for truss with tag: " << this->getTag() << endln;
   return -1;
+}
 int
   int nodalDOF = numDOF/2;
 #ifdef _G3DEBUG
   if (nodalDOF != Raccel1.Size() || nodalDOF != Raccel2.Size()) {
-    g3ErrorHandler->warning("Truss::addInertiaLoadToUnbalance %s\n",
-                            "matrix and vector sizes are incompatable");
+    opserr <<"Truss::addInertiaLoadToUnbalance " <<
+      "matrix and vector sizes are incompatable\n";
     return -1;
+  }
 #endif
     // want to add ( - fact * M R * accel ) to unbalance
+  }
   data(5) = matDbTag;
   res = theChannel.sendVector(dataTag, commitTag, data);
   if (res < 0) {
-    g3ErrorHandler->warning("WARNING Truss::sendSelf() - %d failed to send Vector\n",this->getTag());
+    opserr <<"WARNING Truss::sendSelf() - " << this->getTag() << " failed to send Vector\n";
     return -1;
+  }
   // truss then sends the tags of it's two end nodes
   res = theChannel.sendID(dataTag, commitTag, connectedExternalNodes);
   if (res < 0) {
-    g3ErrorHandler->warning("WARNING Truss::sendSelf() - %d failed to send Vector\n",this->getTag());
+    opserr <<"WARNING Truss::sendSelf() - " << this->getTag() << " failed to send Vector\n";
     return -2;
+  }
   // finally truss asks it's material object to send itself
   res = theMaterial->sendSelf(commitTag, theChannel);
   if (res < 0) {
-    g3ErrorHandler->warning("WARNING Truss::sendSelf() - %d failed to send its Material\n",this->getTag());
+    opserr <<"WARNING Truss::sendSelf() - " << this->getTag() << " failed to send its Material\n";
     return -3;
+  }
   return 0;
+}
   // internal data with the data in the Vector
   static Vector data(7);
   res = theChannel.recvVector(dataTag, commitTag, data);
   if (res < 0) {
-    g3ErrorHandler->warning("WARNING Truss::recvSelf() - failed to receive Vector\n");
+    opserr <<"WARNING Truss::recvSelf() - failed to receive Vector\n";
     return -1;
+  }
   this->setTag((int)data(0));
   dimension = data(1);
   M = data(6);
   // truss now receives the tags of it's two external nodes
   res = theChannel.recvID(dataTag, commitTag, connectedExternalNodes);
   if (res < 0) {
-    g3ErrorHandler->warning("WARNING Truss::recvSelf() - %d failed to receive ID\n", this->getTag());
+    opserr <<"WARNING Truss::recvSelf() - " << this->getTag() << " failed to receive ID\n";
     return -2;
+  }
   // finally truss creates a material object of the correct type,
   // sets its database tag and asks this new object to recveive itself.
       delete theMaterial;
     // create a new material object
     theMaterial = theBroker.getNewUniaxialMaterial(matClass);
     if (theMaterial == 0) {
-      g3ErrorHandler->warning("WARNING Truss::recvSelf() - %d failed to get a blank Material of type %d\n",
-                              this->getTag(), matClass);
+      opserr <<"WARNING Truss::recvSelf() - " << this->getTag()
+        << " failed to get a blank Material of type " << matClass << endln;
       return -3;
+    }
+  }
   theMaterial->setDbTag(matDb); // note: we set the dbTag before we receive the material
   res = theMaterial->recvSelf(commitTag, theChannel, theBroker);
   if (res < 0) {
-    g3ErrorHandler->warning("WARNING Truss::recvSelf() - %d failed to receive its Material\n", this->getTag());
+    opserr <<"WARNING Truss::recvSelf() - "<< this->getTag() << "failed to receive its Material\n";
     return -3;
+  }
   return 0;
+}
     return 0;
+}
 void
-Truss::Print(ostream &s, int flag)
+Truss::Print(OPS_Stream &s, int flag)
+{
     // compute the strain and axial force in the member
     double strain, force;
     strain = theMaterial->getStrain();
     force = A * theMaterial->getStress();
             (*theVector)(i) = (*t)(0,i)*force;
           s << " \n\t unbalanced load: " << *theVector;
+        }
         s << " \t Material: " << *theMaterial;
-        s << endl;
+        s << endln;
     } else if (flag == 1) {
         s << this->getTag() << "  " << strain << "  ";
-        s << force << endl;
+        s << force << endln;
+    }
+}
 double
 Truss::computeCurrentStrain(void) const

Subversion Repositories OpenSees

(root)/trunk/SRC/element/truss/Truss.cpp @ 1271 - Rev 1231 → 1271