In this paper, the threshold of chaos in a forced vibration system with a square nonlinear term is given by using the Melnikov method. The effects of parameters on the threshold are analyzed. The known results obtained respectively by approximate 2-D iterated map and by numerically integrating the differential equation are compared with the results here.