Method getResponse

SRC/element/mixedBeamColumn/MixedBeamColumn2d.cpp:1384–1516

int MixedBeamColumn2d::getResponse(int responseID, Information &eleInfo) {
  if (responseID == 1) { // global forces
    return eleInfo.setVector(this->getResistingForce());

  } else if (responseID == 2) { // local forces
    // Axial
    double N = internalForce(0);
    theVector(3) =  N;
    theVector(0) = -N+p0[0];

    // Moments about z and shears along y
    double M1 = internalForce(1);
    double M2 = internalForce(2);
    theVector(2)  = M1;
    theVector(5) = M2;
    double L = crdTransf->getInitialLength();
    double V = (M1+M2)/L;
    theVector(1) =  V+p0[1];
    theVector(4) = -V+p0[2];

    return eleInfo.setVector(theVector);

  } else if (responseID == 3) { // basic forces
    return eleInfo.setVector(internalForce);

  } else if (responseID == 4) { // section deformation (from forces)

    int i;
    Vector tempVector(2*numSections);
    tempVector.Zero();
    for ( i = 0; i < numSections; i++ ){
      tempVector(2*i)   = sectionDefFibers[i](0);
      tempVector(2*i+1) = sectionDefFibers[i](1);
    }

    return eleInfo.setVector(tempVector);

  } else if (responseID == 5) { // plastic section deformation (from forces)

    int i;
    Vector tempVector(2*numSections);
    Vector sectionForce(NDM_SECTION);
    Vector plasticSectionDef(NDM_SECTION);
    tempVector.Zero();
    for ( i = 0; i < numSections; i++ ){

      getSectionStress(i,sectionForce);
      getSectionTangent(i,2,ks);
      invertMatrix(NDM_SECTION,ks,fs);

      plasticSectionDef = sectionDefFibers[i] - fs*sectionForce;

      tempVector(2*i)   = plasticSectionDef(0);
      tempVector(2*i+1) = plasticSectionDef(1);
    }

    return eleInfo.setVector(tempVector);