Evolutionary algorithms (EAs) have received increasing interests both in the academy and industry. One main difficulty in applying EAs to real-world applications is that EAs usually need a large number of fitness evaluations before a satisfying result can be obtained. However, fitness evaluations are not always straightforward in many real-world applications. Either an explicit fitness function does not exist, or the evaluation of the fitness is computationally very expensive. In both cases, it is necessary to estimate the fitness function by constructing an approximate model. In this paper, a comprehensive survey of the research on fitness approximation in evolutionary computation is presented. Main issues like approximation levels, approximate model management schemes, model construction techniques are reviewed. To conclude, open questions and interesting issues in the field are discussed.Evolutionary computation has found a wide range of applications in various fields of science and engineering. Among others, evolutionary algorithms have been proved to be powerful global optimizers. Generally, evolutionary algorithms outperform conventional optimization algorithms for problems which are discontinuous, non-differential, multi-modal, noisy and not well-defined problems, such as art design, music composition and experimental designs [76]. Besides, evolutionary algorithms are also well suitable for multi-criteria problems.Despite the great successes achieved in real-world applications, evolutionary algorithms have also encountered many challenges. For most evolutionary algorithms, a large number of fitness evaluations (performance calculations) are needed before a well acceptable solution can be found. In many real-world applications, fitness evaluation is not trivial. There are several situations in which fitness evaluation becomes difficult and computationally efficient approximations of the fitness function have to be adopted.Several issues need to be addressed in employing fitness approximations in evolutionary computation. First, which levels of the fitness approximation should be used. While an experimental verification can be seen as the true fitness value of a given solution, fully computational simulations, simplified computational simulations as well as functional approximations (meta-models) can be used for fitness calculation. So far, several models have been used for fitness approximation. The most popular ones are polynomials (often known as response surface methodology), the kriging model, most popular in design and analysis of computer experiments (DACE), the feedforward neural networks, including multi-layer perceptrons and radialbasis-function networks and the support vector machines. Due to the lack of data and the high dimensionality of input space, it is very difficult to obtain a perfect global functional approximation of the original fitness function. To tackle this problem, two main measures can be taken.Firstly, the approximate model should be used together with the original fitness functi...