Starting with the diagonal spacetime metric tensor, the Einstein gravitational field equation is solved, and a set of exact (3 + 1) dimensional cylindrically symmetric wave solutions with two arbitrary functions are found. In these solutions all nonvanishing components of spacetime metric tensor are varying with the same propagating factor (ct-z) while the waves are travelling along z axis. The physical picture and the condition of positive energy density of the wave solutions are discussed.