This paper presents a novel application of machine learning models to clarify the intricate behaviors of expansive soils, focusing on the impact of sand content, saturation level, and dry density. Departing from conventional methods, this research utilizes a data-centric approach, employing a suite of sophisticated machine learning models to predict soil properties with remarkable precision. The inclusion of a 30% sand mixture is identified as a critical threshold for optimizing soil strength and stiffness, a finding that underscores the transformative potential of sand amendment in soil engineering. In a significant advancement, the study benchmarks the predictive power of several models including extreme gradient boosting (XGBoost), gradient boosting regression (GBR), random forest regression (RFR), decision tree regression (DTR), support vector regression (SVR), symbolic regression (SR), and artificial neural networks (ANNs and proposed ANN-GMDH). Symbolic regression equations have been developed to predict the elasticity modulus and unconfined compressive strength of the investigated expansive soil. Despite the complex behaviors of expansive soil, the trained models allow for optimally predicting the values of unconfined compressive parameters. As a result, this paper provides for the first time a reliable and simply applicable approach for estimating the unconfined compressive parameters of expansive soils. The proposed ANN-GMDH model emerges as the pre-eminent model, demonstrating exceptional accuracy with the best metrics. These results not only highlight the ANN’s superior performance but also mark this study as a groundbreaking endeavor in the application of machine learning to soil behavior prediction, setting a new benchmark in the field.