print out element data
| 326 | |
| 327 | //print out element data |
| 328 | void BrickUP::Print( OPS_Stream &s, int flag ) |
| 329 | { |
| 330 | |
| 331 | if (flag == 2) { |
| 332 | |
| 333 | s << "#Brick\n"; |
| 334 | |
| 335 | int i; |
| 336 | const int numNodes = 8; |
| 337 | const int nstress = 6 ; |
| 338 | |
| 339 | for (i=0; i<numNodes; i++) { |
| 340 | const Vector &nodeCrd = nodePointers[i]->getCrds(); |
| 341 | const Vector &nodeDisp = nodePointers[i]->getDisp(); |
| 342 | s << "#NODE " << nodeCrd(0) << " " << nodeCrd(1) << " " << nodeCrd(2) |
| 343 | << " " << nodeDisp(0) << " " << nodeDisp(1) << " " << nodeDisp(2) << endln; |
| 344 | } |
| 345 | |
| 346 | // spit out the section location & invoke print on the scetion |
| 347 | const int numMaterials = 8; |
| 348 | |
| 349 | static Vector avgStress(7); |
| 350 | static Vector avgStrain(nstress); |
| 351 | avgStress.Zero(); |
| 352 | avgStrain.Zero(); |
| 353 | for (i=0; i<numMaterials; i++) { |
| 354 | avgStress += materialPointers[i]->getStress(); |
| 355 | avgStrain += materialPointers[i]->getStrain(); |
| 356 | } |
| 357 | avgStress /= numMaterials; |
| 358 | avgStrain /= numMaterials; |
| 359 | |
| 360 | s << "#AVERAGE_STRESS "; |
| 361 | for (i=0; i<7; i++) |
| 362 | s << avgStress(i) << " " ; |
| 363 | s << endln; |
| 364 | |
| 365 | s << "#AVERAGE_STRAIN "; |
| 366 | for (i=0; i<nstress; i++) |
| 367 | s << avgStrain(i) << " " ; |
| 368 | s << endln; |
| 369 | |
| 370 | /* |
| 371 | for (i=0; i<numMaterials; i++) { |
| 372 | s << "#MATERIAL\n"; |
| 373 | // materialPointers[i]->Print(s, flag); |
| 374 | s << materialPointers[i]->getStress(); |
| 375 | } |
| 376 | */ |
| 377 | |
| 378 | } else { |
| 379 | |
| 380 | s << endln ; |
| 381 | s << "Eight Node BrickUP \n" ; |
| 382 | s << "Element Number: " << this->getTag() << endln ; |
| 383 | s << "Node 1 : " << connectedExternalNodes(0) << endln ; |
| 384 | s << "Node 2 : " << connectedExternalNodes(1) << endln ; |
| 385 | s << "Node 3 : " << connectedExternalNodes(2) << endln ; |