Fiber reinforced polymer (FRP) serves as a prospective alternative to reinforcement in concrete slabs. However, similarly to traditional reinforced concrete slabs, FRP reinforced concrete slabs are susceptible to punching shear failure. Accounts of the insufficient consideration of impact factors, existing empirical models and design provisions for punching strength of FRP reinforced concrete slabs have some problems such as high bias and variance. This study established machine learning-based models to accurately predict the punching shear strength of FRP reinforced concrete slabs. A database of 121 groups of experimental results of FRP reinforced concrete slabs are collected from a literature review. Several machine learning algorithms, such as artificial neural network, support vector machine, decision tree, and adaptive boosting, are selected to build models and compare the performance between them. To demonstrate the predicted accuracy of machine learning, this paper also introduces 6 empirical models and design codes for comparative analysis. The comparative results demonstrate that adaptive boosting has the highest predicted precision, in which the root mean squared error, mean absolute error and coefficient of determination of which are 29.83, 23.00 and 0.99, respectively. GB 50010-2010 (2015) has the best predicted performance among these empirical models and design codes, and ACI 318-19 has the similar result. In addition, among these empirical models, the model proposed by El-Ghandour et al. (1999) has the highest predicted accuracy. According to the results obtained above, SHapley Additive exPlanation (SHAP) is adopted to illustrate the predicted process of AdaBoost. SHAP not only provides global and individual interpretations, but also carries out feature dependency analysis for each input variable. The interpretation results of the model reflect the importance and contribution of the factors that influence the punching shear strength in the machine learning model.