The safety and economy of an infrastructure project depends on the material and design equations used to simulate the performance of a particular member. A variety of materials can be used in conjunction to achieve a composite action, such as a hollow steel section filled with concrete, which can be successfully utilized in the form of an axially loaded member. This study aims to model the ultimate compressive strength (Pu) of concrete-filled hollow steel sections (CFSS) by formulating a mathematical expression using gene expression programming (GEP). A total of 149 datapoints were obtained from the literature, considering ten input parameters, including the outer diameter of steel tube (D), wall thickness of steel tube, compressive strength of concrete (fc’), elastic modulus of concrete (Ec), yield strength of steel (fv), elastic modulus of steel (Es), length of the column (L), confinement factor (ζ), ratio of D to thickness of column, and the ratio of length to D of column. The performance of the developed models was assessed using coefficient of regression R2, root mean squared error RMSE, mean absolute error MAE and comparison of regression slopes. It was found that the optimal GEP Model T3, having number of chromosomes Nc = 100, head size Hs = 8 and number of genes Ng = 3, outperformed all the other models. For this particular model, R2overall equaled 0.99, RMSE values were 133.4 and 162.2, and MAE = 92.4 and 108.7, for training (TR) and testing (TS) phases, respectively. Similarly, the comparison of regression slopes analysis revealed that the Model T3 exhibited the highest R2 of 0.99 with m = 1, in both the TR and TS stages, respectively. Finally, parametric analysis showed that the Pu of composite steel columns increased linearly with the value of D, t and fy.