A ship–ice–water interaction model is established using smoothed-particle hydrodynamics (SPH) to predict the ice breaking resistance of the icebreaker in the Yellow River effectively. This method includes the numerical process of the constitutive equation, yield criterion, and the coupling model in SPH. The ice breaking resistance is determined under different conditions. The numerical results of the ice breaking resistance agree with the empirical formula results. Results show that the prediction accuracy of ice resistance is less than 17.6%compared with the empirical formula in the level ice. The method can also be extended to predict the floe motion and ice breaking resistance in actual river channels. The validation against the empirical formula indicates that the proposed ship–ice–water SPH method can predict the ice breaking resistance of icebreakers in actual rivers effectively. The predicted ice breaking resistance is analyzed under different conditions. The ice breaking resistance increases with increasing bending strength and ice thickness, and the latter is the most important factor influencing ice resistance.