Nowadays, elliptical sections are among the geometric shapes that are becoming more and more popular in the architecture world. Finite element analysis (FEA) software, such as ABAQUS, is used to simulate elliptical concrete columns confined with fiber reinforced polymer (FRP) jackets based on experimental data published in the literature. The validity of the finite element modeling (FEM) approach was confirmed by comparing it to experimental data obtained from 45 specimens in existing studies, demonstrating a favorable agreement. Additionally, a parametric study was also carried out, which included simulating 40 more specimens, resulting in a total of 85 specimens for analysis. The effect of various test variables such as sectional aspect ratio, the amount of FRP layers, and concrete strength on the elliptical column’s behavior is investigated. The axial compressive strength is predicted using current models from the literature. The outcomes revealed that the higher unconfined concrete strength decreased FRP confinement efficacy. Insufficient confinement with post-peak softening is more likely in FRP-confined high strength concrete columns, especially if the confinement was not stiff enough. Also, by adding more FRP layers, the stress and strain capabilities of elliptical concrete columns confined with FRP composites are improved. Physical experiments require a significant investment of time and resources, but they can produce valuable results. Finite element analysis, on the other hand, depends heavily on the modeler's competence and can minimize the number of experiments required by employing computer simulation, but it also requires high computer settings. In recent years, there has been a major advancement in the usage of machine learning (ML) algorithms in the applications of FRP composites with different concrete components. Taking this into consideration, this investigation is extended to estimate the compressive strength of elliptical FRP-confined columns using four tree-based ML algorithms including Decision Tree, Random Forest, Gradient Boosting and XGBoost. The XGBoost model achieved the highest predictive accuracy with the R2 values of 0.95 for both the compressive strength and axial strain targets.