Function TclCommand_addUCFiberSection

SRC/runtime/commands/modeling/section.cpp:1630–1730

int
TclCommand_addUCFiberSection(ClientData clientData, Tcl_Interp *interp,
                             int argc, TCL_Char ** const argv)
{
  assert(clientData != nullptr);
  BasicModelBuilder* builder = static_cast<BasicModelBuilder*>(clientData);
  G3_Runtime *rt = G3_getRuntime(interp);
  int secTag;

  if (argc < 4)
    return TCL_ERROR;

  if (Tcl_GetInt(interp, argv[2], &secTag) != TCL_OK) {
    opserr << "could not read section tag\n";
    return TCL_ERROR;
  }

  builder->setCurrentSectionBuilder(secTag);

  // first create an empty FiberSection

  SectionForceDeformation *section = nullptr;
  FiberSection2d *section2d = nullptr;
  FiberSection3d *section3d = nullptr;
  bool computeCentroid = false;

  int NDM = builder->getNDM();
  if (NDM == 2) {
    section2d = new FiberSection2d(secTag, 30, computeCentroid);
    section = section2d;
    // SectionForceDeformation *section = new FiberSection(secTag, 0, 0);
  } else if (NDM == 3) {
    UniaxialMaterial *theGJ = new ElasticMaterial(0, 1e10);
    section3d =
        new FiberSection3d(secTag, 30, *theGJ, computeCentroid);
    section = section3d;
    delete theGJ;
  }

  if (section == nullptr) {
    return TCL_ERROR;
  }

  //
  // Now parse the output file containing the fiber data,
  // create fibers and add them to the section
  //

  // open the file
  TCL_Char *fileName = argv[3];
  std::ifstream theFile;
  theFile.open(fileName, std::ios::in);
  if (!theFile) {
    opserr << "section UCFiber - could not open file named " << fileName;
    return TCL_ERROR;
  } else {
    int foundStart = 0;
    static char garbage[100];