The ship shows strong nonlinear characteristics when surf-riding/ broaching occurs, which is a failure mode with the most complex mechanism in the second generation intact stability. The wave drift force generated by the ship in low-frequency motion will affect the motion state of the ship. The wave drift force calculated by the empirical formula is therefore introduced into the wave force formula to study the numerical calculation method of the surf-riding/broaching with consideration of the influence of the wave drift force based on the MMG model of four degrees of freedom maneuvering motion equation. The variation law of the ship motion under different wavelength ratios and different rudder angles is then analyzed. In addition, the critical dangerous region of the ship is also predicted to judge the safe speed and stable heading angle during the sailing of the ship. The simulation results are illustrated and compared with the experimental data of Japan ship Mutual Insurance Association and the limits of the standard formula, thereby providing suggestions on avoiding dangerous regions.