Plate end (PE) debonding and intermediate crack (IC) debonding are the two main failure modes of beams strengthened with fiber-reinforced polymer (FRP) in flexure. Therefore, it is essential to clarify the force state of the structure when debonding occurs in strengthened beams. This paper collected 229 beams with debonding failure as the database, of which 128 were PE debonding and 101 were IC debonding. Correlation and grey correlation analysis were used to establish the indicator systems for predicting PE and IC debonding and to identify the critical indicators among them. Five machine learning models, linear regression, ridge regression, decision trees, random forests, and back propagation (BP) neural networks, were used to build the two debonding prediction models. Optimization of the best prediction among the five machine learning models took place using the Dung Beetle Optimizer (DBO) algorithm, which has competitive performance with state-of-the-art optimization approaches in terms of convergence rate, solution accuracy, and stability. Finally, the optimal prediction model was compared with the models suggested by codes, and it was found that the established model can well predict PE and IC debonding.