The paper illustrates the development of metamodels of the response of steel piles driven in sand and subjected to pull-out. The metamodels are created for the prediction of the pile tensile capacity and secant stiffness. They were developed using the results of finite element analyses, which made use of finite element models of robustness assessed employing a selection of available data from large-scale model pile tests. Four hundred finite element analyses allowed for the calibration of very accurate metamodels, which were also demonstrated to closely track the outputs of the experimental results. Once calibrated, the metamodels can be used independently from the finite element models they stemmed from. The outcomes of the study show that metamodels of piles response can yield very accurate results within a wide and realistic range of soil-pile configuration, avoiding the laborious implementation and computational cost which underpins the use of finite element models. As the use of metamodels in this context is new, the paper relies on particularly simplified problem, but the procedure could be extended to accommodate modelling features of higher complexity.