For the vibration analysis of 3-dimensional piping system containing fluid flow, a transfer matrix method is presented. The fluid velocity and pressure were considered, that coupled to longitudinal and flexural vibrations. Transfer matrices and point matrices were derived from direct solutions of the differential equations of motion of pipe conveying fluids, and the variations of natural frequency with flow velocity for 3-dimensional piping system were investigated.