We present an iterative method for simulating beam propagation in nonlinear media using Hamiltonian ray tracing. The Wigner distribution function of the input beam is computed at the entrance plane and is used as the initial condition for solving the Hamiltonian equations. Examples are given for the study of periodic self-focusing, spatial solitons, and Gaussian-Schell model in Kerr-effect media. Simulation results show good agreement with the split-step beam propagation method. The main advantage of ray tracing, even in the nonlinear case, is that ray diagrams are intuitive and easy to interpret in terms of traditional optical engineering terms, such as aberrations, ray-intercept plots, etc.