Rebars made of fiber-reinforced plastic (FRP) might be the future reinforcing material, replacing mild steel rebars, which are prone to corrosion. The bond characteristics of FRP rebars differ from those of mild steel rebars due to their different stress-strain behavior than mild steel. As a result, determining the bond strength (BS) qualities of FRP rebars is critical. In this work, BS data for FRP rebars was investigated, utilizing non-linear capabilities of gene expression programming (GEP) on 273 samples. The BS of FRP and concrete was considered a function of bar surface (Bs), bar diameter (db), concrete compressive strength (fc′), concrete-cover-bar-diameter ratio (c/d), and embedment-length-bar-diameter ratio (l/d). The investigation of the variable number of genetic parameters such as number of chromosomes, head size, and number of genes was undertaken such that 11 different models (M1–M11) were created. The results of accuracy evaluation parameters, namely coefficient of determination (R2), mean absolute error (MAE), and root mean square error (RMSE) imply that the M11 model outperforms other created models for the training and testing stages, with values of (0.925, 0.751, 1.08) and (0.9285, 0.802, 1.11), respectively. The values of R2 and error indices showed that there is very close agreement between the experimental and predicted results. 30 number chromosomes, 9 head size, and 5 genes yielded the optimum model. The parametric analysis revealed that db, c/d, and l/d significantly affected the BS. The FRP rebar diameter size is greater than 10 mm, whereas a l/d ratio of more than 12 showed a considerable decrease in BS. In contrast, the rise in c/d ratio revealed second-degree increasing trend of BS.