The Melnikov function for prediction of Smale horseshoe chaos is applied to the rf-driven Josephson junction. Linear and quadratic damping resistors are considered. In the latter case the analytic solution including damping and dc bias is used to obtain an improved threshold curve for the onset of chaos. The prediction is compared to new computational solutions. The Melnikov technique provides a good, but slightly low, estimate of the chaos threshold.