The multiple exp-function technique, as an extension to Hirota's perturbation algorithm, is employed to find N-soliton solutions of the 3D Vakhnenko-Parkes (3D-VP) equation. This is a straightforward and effective alternative method. Consequently, one- two- and three-soliton solutions besides the front wave solutions and their wave numbers are obtained, which are completely distinct from the other published literature. 3D and density plots were investigated via the help of the Maple package program by suitably choosing the free parameters in the solutions. Further, the physical meaning for these solutions is clarified. We mention that the propagation of high-frequency waves in a relaxing medium can be investigated by discussing the dynamical behavior of the obtained solutions for the 3D-VP equation.