Method getResponse

SRC/element/pipe/Pipe.cpp:894–996

}

int Pipe::getResponse(int responseID, Information &info) {
    if (theCoordTransf == 0) {
        return 0;
    }

    // update section data
    double A, Jx, Iy, Iz, rho, alphaV;
    if (updateSectionData(A, Jx, Iy, Iz, rho, alphaV) < 0) {
        opserr << "Pipe::getResistingForce failed to update section "
                  "data\n";
        return -1;
    }

    // update material data
    double E, nu, G, alp;
    if (updateMaterialData(E, nu, G, alp) < 0) {
        opserr << "Pipe::getResistingForce failed to update material "
                  "data\n";
        return -1;
    }

    double L = theCoordTransf->getInitialLength();
    double oneOverL = 1.0 / L;

    switch (responseID) {
        case 1:
            // stiffness
            return info.setMatrix(this->getTangentStiff());
        case 2:
            // global forces
            return info.setVector(this->getResistingForce());
        case 3:
            // local forces
            P(0) = -q(0) + p0[0];  // -pb1+plw1
            P(1) =
                (q(1) + q(2)) * oneOverL + p0[1];  // (pb2+pb3)/L+plw2
            P(2) = -(q(3) + q(4)) * oneOverL +
                   p0[3];  // -(pb4+pb5)/L+plw3
            P(3) = -q(5);  // -pb6
            P(4) = q(3);   // pb4
            P(5) = q(1);   // pb2
            P(6) = q(0);   // pb1
            P(7) = -(q(1) + q(2)) * oneOverL +
                   p0[2];  // -(pb2+pb3)/L+plw8
            P(8) =
                (q(3) + q(4)) * oneOverL + p0[4];  // (pb4+pb5)/L+plw9
            P(9) = q(5);                           // pb6
            P(10) = q(4);                          // pb5
            P(11) = q(2);                          // pb3
            return info.setVector(P);
        case 4:
            // basic forces
            return info.setVector(q);
        case 5:
            // basic deformations
            return info.setVector(
                theCoordTransf->getBasicTrialDisp());
        case 6: {