Method addLoad

SRC/element/pipe/Pipe.cpp:547–618

}

int Pipe::addLoad(ElementalLoad *theLoad, double loadFactor) {
    if (theCoordTransf == 0) {
        return 0;
    }

    int type;
    const Vector &data = theLoad->getData(type, loadFactor);
    double L = theCoordTransf->getInitialLength();

    if (type == LOAD_TAG_Beam3dUniformLoad) {
        // transformation matrix
        Matrix T(3, 3);
        Vector xAxis(3), yAxis(3), zAxis(3);
        Vector global_w(3);
        if (theCoordTransf->getLocalAxes(xAxis, yAxis, zAxis) < 0) {
            return -1;
        }
        xAxis.Normalize();
        yAxis.Normalize();
        zAxis.Normalize();
        for (int j = 0; j < 3; ++j) {
            T(0, j) = xAxis(j);
            T(1, j) = yAxis(j);
            T(2, j) = zAxis(j);
        }

        // global wy, wz, wx
        global_w(1) = data(0) * loadFactor;
        global_w(2) = data(1) * loadFactor;
        global_w(0) = data(2) * loadFactor;

        // local
        Vector local_w(3);
        local_w.addMatrixVector(0.0, T, global_w, 1.0);

        double wx = local_w(0);  // Transverse
        double wy = local_w(1);  // Transverse
        double wz = local_w(2);  // Axial (+ve from node I to J)

        this->wx += wx;
        this->wy += wy;
        this->wz += wz;

        double Vy = 0.5 * wy * L;
        double Mz = Vy * L / 6.0;  // wy*L*L/12
        double Vz = 0.5 * wz * L;
        double My = Vz * L / 6.0;  // wz*L*L/12
        double Pw = wx * L;

        // Reactions in basic system
        p0[0] -= Pw;  // plw1
        p0[1] -= Vy;  // plw2
        p0[2] -= Vy;  // plw8
        p0[3] -= Vz;  // plw3
        p0[4] -= Vz;  // plw9

        // Fixed end forces in basic system
        q0[0] -= 0.5 * Pw;  // pb0_1