Method getResponse

SRC/material/section/FiberSection3d.cpp:1344–1422

int 
FiberSection3d::getResponse(int responseID, Information &sectInfo)
{
  static double yLocs[10000];
  static double zLocs[10000];
  static double fiberArea[10000];
  
  if (sectionIntegr != 0) {
    sectionIntegr->getFiberLocations(numFibers, yLocs, zLocs);
  }  
  else {
    for (int i = 0; i < numFibers; i++) {
      yLocs[i] = matData[3*i];
      zLocs[i] = matData[3*i+1];
      fiberArea[i] = matData[3*i+2];      
    }
  }
  
  if (responseID == 5) {
    int numData = 5*numFibers;
    Vector data(numData);
    int count = 0;
    for (int j = 0; j < numFibers; j++) {
      data(count)   = yLocs[j]; // y
      data(count+1) = zLocs[j]; // z
      data(count+2) = fiberArea[j]; // A
      data(count+3) = theMaterials[j]->getStress();
      data(count+4) = theMaterials[j]->getStrain();
      count += 5;
    }
    return sectInfo.setVector(data);	
  } else if (responseID == 55) {
    int numData = 6*numFibers;
    Vector data(numData);
    int count = 0;
    for (int j = 0; j < numFibers; j++) {
      data(count)   = yLocs[j]; // y
      data(count+1) = zLocs[j]; // z
      data(count+2) = fiberArea[j]; // A
      data(count+3) = (double)theMaterials[j]->getTag();
      data(count+4) = theMaterials[j]->getStress();
      data(count+5) = theMaterials[j]->getStrain();	    
      count += 6;
    }
    return sectInfo.setVector(data);		  
  } else  if (responseID == 6) {
    int count = 0;
    for (int j = 0; j < numFibers; j++) {    
      if (theMaterials[j]->hasFailed() == true)
	count++;
    }
    return sectInfo.setInt(count);
  } else  if (responseID == 7) {
    int count = 0;
    for (int j = 0; j < numFibers; j++) {    
      if (theMaterials[j]->hasFailed() == true) {
	count+=1;
      }