The transfer matrix method (TMM) was used for description of harmonic vibrations of piping with transported medium. Apart from 12 well-known mechanical parameters which characterize the state of piping system in each cross section two additional parameters that characterize the vibration of the medium, namely its translation and pressure pulsations were considered. The solution of these equations, which take into account the Poisson contraction of the pipe wall, in the form suitable for the transfer matrix method application was derived. The biggest uncertainty in the analytical modeling is to adopt the boundary conditions for above mentioned 2 parameters for the considered piping section. To solve this problem of identification of the most probable induced frequency we developed the technique of choosing such boundary conditions at which the maximum of energy is confined within the considered piping section. The validity of the approach was tested on some analytical examples. This method was used to analyze the forced vibration of the second circuit loop of unit 1 Zaporizhia Nuclear Power Plant (ZNPP) with VVER-1000 (from Russian: Vodo-Vodyanoi Energetichesky Reactor; Water-Water Power Reactor) arising from turbulent eddies in the flow of steam. Natural frequencies and forms of mechanical, hydrodynamic, and related hydro-mechanical vibration were found, a number of recommendations were given to reduce the vibration levels.