In this work, propagation of harmonic waves in initially stressed cylindrical viscoelastic thick tubes filled with a Newtonian fluid is studied. The tube, subjected to a static inner pressure P(i) and a positive axial stretch 2, will be considered as an incompressible viscoelastic and fibrous material. The fluid is assumed as an incompressible, viscous and dusty fluid. The field equations for the fluid are obtained in the cylindrical coordinates. The governing differential equations of the tube's viscoelastic material are obtained also in the cylindrical coordinates utilizing the theory of small deformations superimposed on large initial static deformations. For the axially symmetric motion the field equations are solved by assuming harmonic wave solutions. A closed form solution can be obtained for equations governing the fluid body, but due to the variability of the coefficients of resulting differential equations of the solid body, such a closed form solution is not possible to obtain. For that reason, equations for the solid body and the boundary conditions are treated numerically by the finite-difference method to obtain the effects of the thickness of the tube on the wave characteristics. Dispersion relation is obtained using the long wave approximation and, the wave velocities and the transmission coefficients are computed. (C) 2009 Elsevier Inc. All rights reserved.