This study is primarily aimed at creating three machine learning models: artificial neural network (ANN), random forest (RF), and k-nearest neighbour (KNN), so as to predict the crippling load (CL) of I-shaped steel columns. Five input parameters, namely length of column (L), width of flange (bf), flange thickness (tf), web thickness (tw) and height of column (H), are used to compute the crippling load (CL). A range of performance indicators, including the coefficient of determination (R2), variance account factor (VAF), a-10 index, root mean square error (RMSE), mean absolute error (MAE) and mean absolute deviation (MAD), are used to assess the effectiveness of the established machine learning models. The results show that all of the three ML (machine learning) models can accurately predict the crippling load, but the performance of ANN is superior: it delivers the highest value of R2 = 0.998 and the lowest value of RMSE = 0.008 in the training phase, as well as the highest value of R2 = 0.996 and the smaller value of RMSE = 0.012 in the testing phase. Additional methods, including rank analysis, reliability analysis, regression plot, Taylor diagram and error matrix plot, are employed to assess the models’ performance. The reliability index (β) of the models is calculated by using the first-order second moment (FOSM) technique, and the result is compared with the actual value. Additionally, sensitivity analysis is performed to check the impact of the input variables on the output (CL), finding that bf has the greatest impact on the crippling load, followed by tf, tw, H and L, in that order. This study demonstrates that ML techniques are useful for developing a reliable numerical tool for measuring the crippling load of I-shaped steel columns. It is found that the proposed techniques can also be used to predict other kinds of failures as well as different kinds of perforated columns.