The Scientific World Journal

Volume 2013, Article ID 641269, 12 pages

http://dx.doi.org/10.1155/2013/641269

## A Quasiphysics Intelligent Model for a Long Range Fast Tool Servo

^{1}College of Mechanical Science and Engineering, Jilin University, Changchun 130022, China^{2}College of Electromechanical Engineering, Changchun University of Technology, Changchun 130012, China

Received 18 July 2013; Accepted 13 August 2013

Academic Editors: C. Bao and A. Szekrenyes

Copyright © 2013 Qiang Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Accurately modeling the dynamic behaviors of fast tool servo (FTS) is one of the key issues in the ultraprecision positioning of the cutting tool. Herein, a quasiphysics intelligent model (QPIM) integrating a linear physics model (LPM) and a radial basis function (RBF) based neural model (NM) is developed to accurately describe the dynamic behaviors of a voice coil motor (VCM) actuated long range fast tool servo (LFTS). To identify the parameters of the LPM, a novel *Opposition-based Self-adaptive Replacement Differential Evolution (OSaRDE)* algorithm is proposed which has been proved to have a faster convergence mechanism without compromising with the quality of solution and outperform than similar evolution algorithms taken for consideration. The modeling errors of the LPM and the QPIM are investigated by experiments. The modeling error of the LPM presents an obvious trend component which is about ±1.15% of the full span range verifying the efficiency of the proposed OSaRDE algorithm for system identification. As for the QPIM, the trend component in the residual error of LPM can be well suppressed, and the error of the QPIM maintains noise level. All the results verify the efficiency and superiority of the proposed modeling and identification approaches.

#### 1. Introduction

In recent years, the fast tool servo (FTS) based single point diamond turning (SPDT) has been gaining increasing applications in extensive areas such as dynamic compensation of machining errors, generation of nonaxisymmetric optical and components, fabrication of microstructured functional surfaces [1–3]. The determination of trajectory tracking strategy for the cutting tool is an ongoing topic for the FTS system to achieve excellent control performance with high accuracy and high robustness, accordingly to achieve fine profile accuracy and high quality of finished surface.

Up to now, various trajectory tracking strategies have been developed. They could be mainly categorized into model-based control approaches relying on the models of the controlled systems and adaptive control approaches with no dependences on the models. As for the adaptive control approaches, the slide mode controller (SMC) has been mainly implemented for the tool positioning of FTS [4, 5]. The SMC does not need accurate model of controlled system and has relatively high robustness to external disturbances and model uncertainties or variations, but only at the expense of positioning accuracy. The model-based controllers which have been applied to FTS generally comprise of proportional integral derivate (PID) controller [6–8], adaptive feedforward cancellation (AFC) controller [3, 9, 10], repetitive controller (RC) [11, 12], and active disturbance rejection controller (ADRC) [13, 14]. To further enhance the robustness and the stability of FTS system, extended state observer (ESO) based compensation approach serving as one of the key components of the ADRC has been extensively embedded in both model-based controllers and adaptive controllers [10, 11, 13–15]. As for the optimal design process of the PID controller, the feedforward component and phase advance parameters of AFC controller, and the extensively employed ESO compensator, accurate model for describing the physical process of FTS should be critically established. Unfortunately, almost all of the FTSs are approximately equivalent to linear second-order dynamic systems [8–10, 12–14], even for the piezoelectric actuated FTS with inherent hysteresis and creep nonlinearity (HCN) effects [4, 15]. It is evident that there should exist inherent nonlinearities in any practical FTS system, the linear approximations of these servo systems would lead to significant modeling errors and consequently deteriorate control performances of FTS systems. Hence, considerable work should be further done to improve the modeling accuracy.

Generally, the system modeling approaches can be mainly divided into two categories, which are physical approach and phenomenological approach, respectively. On one hand, the description capability of physical model is limited by its assumptions. On the other hand, phenomenological model has its own drawbacks, such as outlier and overfitting. Such problems probably lead to conclusions against physical laws, especially when phenomenological model is used for extrapolation [16–19]. So it is believed that integrating these models has the potential to overcome the drawbacks of both physical and phenomenological models [16–19].

Motivated by this, a quasiphysics intelligent model (QPIM) integrating a linear physics model (LPM) and a radial basis function (RBF) based neural model (NM) is proposed to accurately describe the dynamic behaviors of a long range FTS (LFTS) driven by a voice coil motor (VCM), and a novel *Opposition-based Self-adaptive Replacement Differential Evolution (OSaRDE)* algorithm is developed for the parameters identification of the LPM governed by physics laws. The remainder of this paper is organised as follows: the mechanical structure and the modeling procedure of the LFTS is introduced in Section 2; the preliminary of the differential evolution (DE) algorithm and the proposed OSaRDE algorithm is introduced in Section 3; the experiment results and corresponding discussions are carried out in Section 4; the main conclusions of this paper are drawn in Section 5.

#### 2. Modeling and Identification Procedure of the LFTS

##### 2.1. Mechanical Structure of the LFTS

The authors of this paper have developed a LFTS for ultraprecision diamond turning; the mechanical structure of the LFTS is illustrated in Figure 1(a). The LFTS mainly contains two parts: the VCM actuator and the flexure-based guide mechanism. A BEI Kimco VCM is employed as the driving element. It can provide maximum continuous stall forces of 82.74 N (18.6 lbs) with a stroke of ±5.72 mm. Being benefited from the long stroke of the actuator, the stroke of the LFTS could reach up to 1 mm. Two groups of cross shape flexure hinges are specially designed for the guide of the moving component of the LFTS. The two hinges are connected by a linking member crossing the VCM stator. The mechanical design process and performance evaluations of this LFTS are detailed in [20].

##### 2.2. The Quasiphysics Intelligent Model of the LFTS

From an electromagnetic point of view, the VCM can be considered as an inductance component with a resistance [21, 22]. Meanwhile, from the dynamic point of view, the flexure-based mechanism can be equivalent as a damped mass-spring system [8, 23]. With consideration of the mass of the moving part including the flexure hinges and the bobbin of the VCM, the dynamic model of the LFTS mechanism can be obtained, as illustrated in Figure 1(b). The electromagnetic dynamic balance of the VCM in the LFTS, as equivalently shown in Figure 1(b), can be derived based on the Kirchhoff’s law, yielding [21, 22] where is the voltage applied to the VCM, is the back electromotive force (EMF), and is the equivalent back EMF constant. , , and represent the equivalent resistance of voice coils, inductance and applied current of coil, respectively.

With the electrical dynamics in (1) derived, the Flemining’s left-hand rule is then employed to derive the electromagnetic forces for actuation, which can be expressed as [21, 22] where is the number of coil loops, is the magnetic flux densities within the air gap between bobbin and magnets, and is the effective coil lengths considering all loops. Generally, the term can be considered negligible due to the fact that the electrical dynamics of the conducted current is much faster than mechanical behaviors; thus

Based on Newton’s second law of motion, the differential equation of dynamic motion for the LFTS mechanism can be given as follows: where and represent the equivalent moving mass of the flexure-based mechanism and the bobbin of the VCM, respectively. and represent the equivalent damping coefficient and the equivalent stiffness of the flexure-based mechanism, respectively.

Submitting (3) into (4) yields

If this mechanical system is a linear time-invariant (LTI) system, the transfer function of this LFTS system can be obtained as follows by applying the Laplace transform: where , , .

Unfortunately, due to manufacturing tolerances on various crucial dimensions of both bobbin and magnets and possible misalignments of magnetic components, there often exist an uneven magnetic filed within the air gap, as a result, it would lead to a more complex and position-varying magnetic flux density [21, 22]. On the other hand, the temperature in the VCM would increase with the motion ongoing, consequently leads to the variations of the equivalent resistance and inductance of the VCM. Meanwhile, the equivalent damping coefficient and stiffness of the flexure-based mechanism could not be permanent constants due to the inherent properties of the employed material for the mechanical structure. Let be the real output of the LFTS; it would be difficult for the description of due to the presence of noise (uncontrollable) factors, measurement error, and inherent nonlinearities of system. Hence, all the factors making the system to be nonlinear and random would significantly deteriorate the accuracy of the LTI model.

For nonlinear dynamic modelling, neural model has been an extremely popular choice due to its high adaptability to various nonlinear systems [24–27]. In the feedforward neural networks, RBF neural networks are local neural networks with series advantages of the faster training speed, small computing, and simple structure [27]. To enhance the accuracy of the pure physics model as obtained in (6), an RBF neural network (RBFNN) based NM is further employed to model the components with nonlinearities and uncertainties. By introducing this NM, the governing equation of the QPIM for the LFTS mechanism in the time domain can be modified as where represent the inverse Laplace transform operation and represent the RBF based NM. The command signal and its differentiation are employed as the input of the NM, which are three-layer feedforward networks consisting of an input layer, a hidden layer, and an output layer. Generally, the Gaussian function is employed as the basis function of the hidden layer, given by [27] where , , and are the input vector, centre vector, and width of the Gaussian, respectively. denotes the Euclidean norm operation.

The output of the NM is a linear combiner, which is defined by [27] where is the weight between the hidden layer and the output layer; is the number of nodes in the hidden layer.

##### 2.3. Parameters Identification Procedure

The basic idea in the LFTS system identification is to compare the input dependent response of the system and a corresponding mathematic model by certain performance criterions giving a measure to how well the model response fits the system response. The principle of the identification procedure for the parameter identification of the LFTS system is illustrated in Figure 2, where denotes the response of the LFTS system; and denote the response of the LPM and the NM, respectively. As discussed above, the LPM depends on a set of three parameters, that is, . Hence, the objective of the LPM identification is to find a set of parameters that minimize the difference between the response of the real LFTS system and the output of the LPM at each discrete sampling point. Therefore, the identification process lies on minimizing the following criterion: where represents the index of the th discrete sampling point.

The criterion is commonly called a fitness function or objective function reflecting the goodness of solution in the parameter tuning process. The identification problem thus can be treated as a linearly constrained multidimensional optimization problem and formalized as where is the 3-dimensional feasible search space; and denote the upper bounds and the lower bounds of the th parameter , respectively.

serving as objective function maps decision variables into the objective space. Obviously, the fitness landscape of this identification task may have many local optima and a highly complex topology. Hence, an improved OSaRDE algorithm, which is a simple yet efficient heuristic approach, is employed as the searching tool and will be detailed in the following section.

As for the NM, the similar minimization criterion is defined as

For the training process, the command signal and its differentiation are employed as the input signal, while the residual error of the LPM is employed as the target signal. To obtain the optimal weights of the RBFNN, the developed OSaRDE algorithm could also be implemented as shown in [26]. However, the training of the NM to obtain the optimal weights and structure is executed through the NEURB function in Matlab just for simplicity. Following the scheme provided by Matlab, new neuron is created at each iteration, and the criterion of the new network is checked. If it is low enough, the training is finished; otherwise, the next neuron is added. This procedure is repeated until the criterion goal is met or the maximum number of iteration is reached [27].

#### 3. The Opposition-Based Self-Adaptive Replacement Differential Evolution Algorithm

The DE algorithm proposed by Storn and Price (1997) is a heuristic approach for minimizing or maximizing possibly nonlinear and nondifferentiable continuous space functions [28]. Due to some attractive characteristics such as simple differential operator, one-to-one competition scheme, and constructive cooperation between current individuals and memory of the best individuals, DE has been extensively implemented for the identification of both linear and nonlinear systems [29–31]. The obtained results demonstrate that DE is a very promising tool for optimal parameters determination and outperforms other evolutionary technologies [30–32].

##### 3.1. The Classical DE Algorithm

The DE algorithm is a population-based algorithm like general evolution algorithms using the similar operators: crossover, mutation, and selection. In DE, a population of potential solution vectors is initialized covering the entire parameter space at the start, which is then evolved to find optimal solutions through the mutation, crossover, and selecting operation procedures [28, 31]. More specifically DE’s basic strategy can be formalized as follows.

*(S1) Initialization.* Assuming that an optimization task consists of parameters, let be the search space of the problem under consideration; the initial vector population randomly covering the entire parameter space can be expressed as
where denotes the population size; denotes a potential solution for the task at start.

*(S2) Mutation.* For each target vector at the th generation, an associated mutant vector should be generated via a certain mutation operator. The objective of mutation is to enable search diversity in the parameter space as well as to direct the existing object vectors to better results at a suitable time. Motivated by this, various mutation strategies have been proposed, and certain typical strategies can be summarized as [33, 34]
where , , and represent the random and mutually different integers generated within the range , and also different from index ; is a mutation scale factor within the range , usually less than 1; denotes the best individual in generation .

*(S3) Crossover.* After the mutation phase, the crossover operation is applied to each pair of the generated mutant vector and its corresponding target vector to generate a trial vector as follows:
where ; ; is the th independent random number uniformly distributed in the range of . is a randomly chosen index from the set ; CR is a user-specified crossover factor within that controls the diversity of the population [28, 31].

*(S4) Selection.* DE employs a greedy selection process that the better one of new offspring and its parent win the competition. If denotes the objective function under consideration, the next generation can be determined by
For the minimum problem, the current best individual at the th generation is defined as [35]
which denotes the best individual in the previous entire generations is defined as [35]

##### 3.2. The Improved OSaRDE Algorithm

However, the optimum performances of different problems are not only highly dependent on the control parameters and learning strategies involved in DE, but also dependent on the characteristics of the region of the search landscape being explored [36, 37]. For a given task, we should spend a huge amount of time to try through various strategies and fine tune the corresponding parameters to achieve the best optimum performance. Motivated by this dilemma, large number of researchers have devoted themselves to develope adaptive DE algorithm to solve general problems more efficiently. Up to now, a classification, into three macrogroups of the main adaptive strategies employed in the DE scheme can be proposed here. Self-adapting parameter setting of and CR, including the pseudorandom based updating approach [34, 36, 38] and the local optimum search based updating approach [39, 40]. Adaptive selection of mutation strategies according to their capacities of generating promising solutions [36, 37, 41]. Dynamic population sizing strategy based on self-adaptation [42] or progressive reduction [43, 44].

The classical DE algorithm typically converges at a rapid pace in the initial stages of the search procedure and then gradually slows down or keeps unchanged as it approaches global optimum. In case there is no improvement in the fitness in successive generations, it would be an indication that the individuals are clustered together in a local region and the operations are very insignificant to allow any improvement of the evolution process [45]. As is evident from this phenomenon, the information of the evolutionary process contained in these individuals would overlap and become extensively redundant. However, we do not suggest that we should discard the redundant individuals by the dynamic population sizing strategy. If the redundant individuals were simply discarded, the information contained in these individuals would be totally lost. Moreover, the population sizing strategy also has no assistance to the performance of escaping from local attractors.

In order to avoid the redundancy and enhance the population diversity, Ali and Pant (2011) take advantage of Cauchy mutation to perturb the clustered individuals and force them to move to someother location, thereby providing them an opportunity to improve their performance [45]. Being totally different from Cauchy mutation based and other adaptive strategies based DE algorithm mentioned above, a novel *Opposition-based Self-adaptive Replacement Differential Evolution (OSaRDE)* algorithm is proposed in this research to enhance the performance of DE. In the OSaRDE scheme, the redundant inferior individuals are self-adaptively selected according to a criterion defined by Euclidean distance after a specified number of generations, and then each of them is replaced by its opposite based on the opposition-based learning (OBL) strategy. After a period of regular evolutions, the process will be repeated. More specifically OSaRDE’s basic strategy considering the constraint task will be formalized as follows.

###### 3.2.1. The OSaRDE’s Basic Strategy

For the replacement operation in the generation, the current worst individual can be defined as

The redundant inferior individuals are defined as where denotes the number of inferior individuals, which is self-determined according to (22); representing the replacement ratio is a user-defined constant; denotes the operator for Euclidean distance; denotes the total number of the replacement operations during the evolution process. Let , , and represent the number of total generations, the specified number of generations from which the replacement operation begins, and the evolution period during two adjacent replacement operations. Considering that the replacement operation should not be executed at the last generation, can be obtained as

The opposite of each selected inferior individual can be defined as [26, 46] where and represent the predefined boundaries of the optimum task.

By applying the replacement operation, the population in the th generation can be obtained as

By applying the replacement approach to the self-adaptively selected inferior individuals, the evolutionary process dynamically calculates the corresponding opposite according to the current location of each selected individual, thereby forcing them to jump to new solution candidates. The evolution process would possess much higher diversity and accordingly escape from the problem of premature convergence. Then the evolution process could have more opportunities to find better individuals during a period of regular evolutions after a replacement operation.

###### 3.2.2. The Constraint Handling Technique

Generally, most of the optimum problems in engineering application are constrained. Unfortunately, the constraint of the searching space is not considered in the conventional DE scheme. A simple and popular repair operation of the individuals beyond the boundary works as follows [47]: if the th element of the mutant vector is out of the search region , then is reset as

For sake of clarity, the framework of the proposed OSaRDE is emphasized in Algorithm 1.

##### 3.3. A Brief Testing on Benchmark Functions

A comprehensive set of benchmark functions, including eight typical global optimization problems, has been used only for performance verification and comparison of the proposed OSaRDE. The definition of the selected benchmark functions and their global optimum(s) are listed in the appendix [45, 48].

The eight algorithms described in the appendix have been tested and compared with standard DE and Cauchy mutation based DE (MDE) proposed by Ali and Pant (2011) [45]. Experiments related to DE and MDE in [45] have been employed as a reference. In order to have a fair comparison, the control parameters of OSaRDE are set as constant, where and , which are the same as employed in [45]. The number of function evaluations (NFEs) of OSaRDE are kept as 1 × 10^{5}, which is half of that employed in [45]. The DE/cur-to-best/1 is chosen as the mutation strategy. In order to have statistically significant results, each test has been run 30 times. The average absolute error and standard deviation are recorded for comparisons. The corresponding results are given in Table 1, and the best results are highlighted in bold face. As evidence from the results in Table 1, it can be clearly observed that both DE and MDE algorithms give the same results for lower dimension problem, while almost all of the absolute errors obtained by OSaRDE are significantly smaller than that obtained by DE and MDE except for . As for , the average absolute error obtained by OSaRDE is larger than that obtained by MDE, while it is quite smaller than those obtained by DE. The results demonstrate that the proposed OSaRDE has a faster convergence mechanism without compromising with the quality of solution, significantly outperforming DE and MDE for both low and high dimension optimum tasks. It could be a fine choice for the OSaRDE to be applied to the engineering practice.

#### 4. Experiment Results and Discussion

Figure 3 illustrates the testing condition of the experiment part, where the LFTS is installed on a precision machine tool. The command signal to be applied to the LFTS is generated by a computer and converted through a data-acquisition card from ADLINK. The displacement of the LFTS is measured by a Renishaw linear encoder with a resolution of 10 nm. The measured displacement which is a digital signal is then directly gathered and stored in the computer through the data-acquisition card with a sample time of 0.45 ms for further analysis.

##### 4.1. Identification Results of the LPM

To catch some prior knowledge of the LFTS system and accordingly help achieve more accurate parameters of the system, the sweep check is carried out to estimate the damped natural frequency. Thus, a sweep harmonic signal with constant amplitude and varying frequency ranged from 0 Hz to 400 Hz is applied to the VCM actuator. Fourier transform is then applied to the corresponding response, and the result is illustrated in Figure 4. As shown in Figure 4, the first damped natural frequency is about 172.6 Hz. Generally, the damping factor is about 0.05 and less than 0.1. With the roughly estimated first damped natural frequency and damping factor, the proportionality factor can be investigated as 5 × 10^{6}. Thus, the feasible search space of the physical model can be defined as
where and are the scale factors of the first damped natural frequency and damping factor, respectively.

To identify the parameters of the LPM, the step command with 1 V voltage and the corresponding response of the LFTS are employed for the parameter identification due to the reason that the step response contains abundant information of the system dynamics, such as the natural frequency and the damping behavior. To guarantee that the identified results have reached up to the maximum description capacity of the LPM, the optimum procedure is carried out by applying OSaRDE with different population sizes and generations, and the obtained best fitness values are further compared with them obtained by standard DE. The best fitness values are summarized in Table 2, and the convergence processes are illustrated in Figure 5. As given in Table 2, the proposed OSaRDE outperforms DE for low NFEs and gradually tends to be the same for large NFEs from the perspective of best fitness values. It also should be noticed that the best fitness values obtained by OSaRDE with different NFEs remain the same. It means that the OSaRDE is of high robustness for optimum tasks and could afford a much faster convergence rate than DE, thereby effectively reducing the computational time.

The transfer function of the LFTS can be obtained by introducing the optimum parameters corresponding to the best fitness value, that is

##### 4.2. Comprehensive Results of the QPIM

In this experiment, an input excitation scheme with varying frequency and amplitude as shown in Figure 6(a) is applied to the VCM, and the corresponding output of the LFTS is illustrated in Figure 6(b). The 13600 sampling data are divided into two parts in which the former 6000 groups of data are used for training the NM and the latter 7600 groups of data are used for prediction. The modeling and prediction performance of the trained RBF NM is illustrated in Figure 7. As is evident from Figure 7, the NM could well approximate to the residual error of the LPM and present a good prediction capacity. To give a comprehensive comparison of the performances between the LPM and the proposed QPIM with neural model calibration, the modeling errors of the two models are illustrated in Figure 8. It can be seen that the maximum modeling error of the LPM is about ±1.15% of the full span range; meanwhile, obvious trend component attributing to the unmodeled parts of the physical process which could not be described by the LPM can be observed. The relatively small error verifies the efficiency of the proposed OSaRDE algorithm for system identification. As for the QPIM, the trend component in the residual modeling error of the LPM is well suppressed and the corresponding error maintains noise level. The result demonstrates that the QPIM is able to coincide more effectively with the response of the LFTS system; namely, the proposed modeling and identification method is of excellent performance.

#### 5. Conclusion

A quasiphysics intelligent model (QPIM) integrating a linear physics model (LPM) and a radial basis function neural network (RBFNN) based neural model (NM) is developed in this paper to accurately describe the dynamic behaviors of a voice coil motor (VCM) actuated long range fast tool servo (LFTS). A linear second-order differential equation of the LFTS system is established according to physics laws, and a novel Opposition-based Self-adaptive Replacement Differential Evolution (OSaRDE) algorithm is further developed for the parameter identification of the LPM. The NM is developed for the calibration of the unmodeled parts of the physics process attributing to the nonlinearity effect of the system. Both numerical and experiment examinations are carried out to respectively assess the performance of the OSaRDE algorithm and the QPIM. Certain conclusions can be drawn as follows.(1)A comprehensive set of benchmark functions, including eight typical global optimization problems, has been used for performance verification and comparison of the proposed OSaRDE. In comparison with both standard DE and a modified DE (MDE), it is observed that the OSaRDE performs remarkably better for both low and high dimension optimum tasks. The observation demonstrates that the proposed OSaRDE has a faster convergence mechanism without compromising with the quality of solution and outperforms DE and MDE taken for consideration. We may suggest that it should be a fine choice for the OSaRDE to be implemented for engineering practice.(2)As is evident from the identification results of the LPM and the QPIM, the modeling error of the LPM presents an obvious trend component which is about ±1.15% of the full span range. The relatively small error verifies the efficiency of the proposed OSaRDE algorithm for system identification. The residual modeling error is further reduced by NM calibration; the trend component in the residual error of the LPM is well suppressed, and the error of the QPIM maintains noise level. The result demonstrates that the QPIM is able to coincide more effectively with the response of the LFTS system, namely, the proposed modeling and identification method is of excellent performance.

#### Appendix

#### Benchmark Problems

* Eason Function (EP):*

with and .

It is a multimodal, nonseparable function having several optima.

* Sphere Function (SP):*

with and .

It is a simple, continuous unimodal, separable, and highly convex function. It serves as test case for validating the convergence speed of an algorithm.

* Ackley’s Function (ACK):*

with and .

It has numerous local minima, but its complexity is moderated. Any search strategies that analyzes a wider region will be able to cross the valley among the optima and achieve the global optimum.

* Griewank Function (GW):*

with and .

It is a highly multimodal, nonseparable function. It has many regularly distributed local minima. It tests both convergence speed and the stability to escape from a shallow local minimum.

* Levy and Montalvo 2 Problem (LM2):*

with and .

It is a multimodal, nonseparable function having several local optima.

* Schwefel’s Problem (SWF):*

with and .

It is a multimodal function with very deep sinusoidal interactions. It is generally considered to be difficult for optimization algorithms.

* Rosenbrock Problem (RB):*

with and .

As the number of dimensions increase, it ceases to be unimodal. The optimum of RB function lies inside a long, narrow, parabolic-shaped float valley. It helps in testing the ability of an algorithm to prevent premature convergence.

* Rastrigin’s Function (RG):*

with and .

It is a multimodal, nonseparable function having large numbers of local optima.

#### Acknowledgments

The authors are grateful to the financial support from the NSF of China (51175221; 51075041; 50995077), the Ministry of Science and Technology (MoSt) of China (2008AA04Z125), the Ministry of Education (MoE) of China (20070183104), and the Department of Science and Technology (DoST) of Jilin Province, China (20090337; 20100359; 20130522155JH).

#### References

- W. Gao, M. Tano, T. Araki, S. Kiyono, and C. H. Park, “Measurement and compensation of error motions of a diamond turning machine,”
*Precision Engineering*, vol. 31, no. 3, pp. 310–316, 2007. View at Publisher · View at Google Scholar · View at Scopus - D. P. Yu, Y. S. Wong, and G. S. Hong, “Ultraprecision machining of micro-structured functional surfaces on brittle materials,”
*Journal of Micromechanics and Microengineering*, vol. 21, no. 9, Article ID 095011, 2011. View at Publisher · View at Google Scholar · View at Scopus - H.-S. Kim, K.-I. Lee, K.-M. Lee, and Y.-B. Bang, “Fabrication of free-form surfaces using a long-stroke fast tool servo and corrective figuring with on-machine measurement,”
*International Journal of Machine Tools and Manufacture*, vol. 49, no. 12-13, pp. 991–997, 2009. View at Publisher · View at Google Scholar · View at Scopus - Y. Altintas and A. Woronko, “A piezo tool actuator for precision turning of hardened shafts,”
*CIRP Annals*, vol. 51, no. 1, pp. 303–306, 2002. View at Google Scholar · View at Scopus - D. P. Yu, G. S. Hong, and Y. S. Wong, “Profile error compensation in fast tool servo diamond turning of micro-structured surfaces,”
*International Journal of Machine Tools and Manufacture*, vol. 52, no. 1, pp. 13–23, 2012. View at Publisher · View at Google Scholar · View at Scopus - O. Sosnicki, A. Pages, C. Pacheco, and T. Maillard, “Servo piezo tool SPT400MML for the fast and precise machining of free forms,”
*International Journal of Advanced Manufacturing Technology*, vol. 47, no. 9–12, pp. 903–910, 2010. View at Publisher · View at Google Scholar · View at Scopus - Q. Liu, X. Zhou, P. Xu, Q. Zou, and C. Lin, “A flexure-based long-stroke fast tool servo for diamond turning,”
*International Journal of Advanced Manufacturing Technology*, vol. 59, no. 9–12, pp. 859–867, 2012. View at Publisher · View at Google Scholar · View at Scopus - S. Rakuff and J. F. Cuttino, “Design and testing of a long-range, precision fast tool servo system for diamond turning,”
*Precision Engineering*, vol. 33, no. 1, pp. 18–25, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. F. Byl, S. J. Ludwick, and D. L. Trumper, “A loop shaping perspective for tuning controllers with adaptive feedforward cancellation,”
*Precision Engineering*, vol. 29, no. 1, pp. 27–40, 2005. View at Publisher · View at Google Scholar · View at Scopus - M. F. Byl,
*Design and control of a long stroke fast tool servo [Ph.D. thesis]*, Massachusetts Institute of Technology, 2005. - M. Crudele and T. R. Kurfess, “Implementation of a fast tool servo with repetitive control for diamond turning,”
*Mechatronics*, vol. 13, no. 3, pp. 243–257, 2003. View at Publisher · View at Google Scholar · View at Scopus - H. Wang and S. Yang, “Design and control of a fast tool servo used in noncircular piston turning process,”
*Mechanical Systems and Signal Processing*, vol. 36, no. 1, pp. 87–94, 2013. View at Publisher · View at Google Scholar · View at Scopus - D. Wu and K. Chen, “Design and analysis of precision active disturbance rejection control for noncircular turning process,”
*IEEE Transactions on Industrial Electronics*, vol. 56, no. 7, pp. 2746–2753, 2009. View at Publisher · View at Google Scholar · View at Scopus - D. Wu, K. Chen, and X. Wang, “Tracking control and active disturbance rejection with application to noncircular machining,”
*International Journal of Machine Tools and Manufacture*, vol. 47, no. 15, pp. 2207–2217, 2007. View at Publisher · View at Google Scholar · View at Scopus - D. P. Yu, G. S. Hong, and Y. S. Wong, “Integral Sliding Mode Control for Fast Tool Servo diamond turning of micro-structured surfaces,”
*International Journal of Automation Technology*, vol. 5, no. 1, pp. 4–10, 2011. View at Google Scholar · View at Scopus - V. R. Joseph and S. N. Melkote, “Statistical adjustments to engineering models,”
*Journal of Quality Technology*, vol. 41, no. 4, pp. 362–375, 2009. View at Google Scholar · View at Scopus - X. Wang, B. Li, K. Wang, and S. Wu, “An engineering-statistical model for synthesis process of nanomaterials,” in
*Proceedings of the International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering (ICQR2MSE '11)*, pp. 72–76, Xi'an, China, June 2011. View at Publisher · View at Google Scholar · View at Scopus - R. K. Singh, V. R. Joseph, and S. N. Melkote, “A statistical approach to the optimization of a laser-assisted micromachining process,”
*International Journal of Advanced Manufacturing Technology*, vol. 53, no. 1–4, pp. 221–230, 2011. View at Publisher · View at Google Scholar · View at Scopus - X. Deng, V. R. Joseph, W. Mai, L. W. Zhong, and C. F. J. Wu, “Statistical approach to quantifying the elastic deformation of nanomaterials,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 106, no. 29, pp. 11845–11850, 2009. View at Publisher · View at Google Scholar · View at Scopus - Q. Liu,
*Study on freeform optical surfaces diamond turning assisted by fast tool servo [Ph.D. thesis]*, Jilin University, 2012. - P. C.-P. Chao, C.-L. Lai, and J.-S. Huang, “Intelligent actuation strategy for a three-DOF four-wire type optical pickup,”
*Sensors and Actuators A*, vol. 117, no. 1, pp. 28–40, 2005. View at Publisher · View at Google Scholar · View at Scopus - P. C.-P. Chao, C.-L. Lai, and J.-S. Huang, “Nonlinear dynamic analysis and actuation strategy for a three-DOF four-wire type optical pickup,”
*Sensors and Actuators A*, vol. 105, no. 2, pp. 171–182, 2003. View at Publisher · View at Google Scholar · View at Scopus - Y. Tian, B. Shirinzadeh, and D. Zhang, “A flexure-based mechanism and control methodology for ultra-precision turning operation,”
*Precision Engineering*, vol. 33, no. 2, pp. 160–166, 2009. View at Publisher · View at Google Scholar · View at Scopus - C. X. Wong and K. Worden, “Generalised NARX shunting neural network modelling of friction,”
*Mechanical Systems and Signal Processing*, vol. 21, no. 1, pp. 553–572, 2007. View at Publisher · View at Google Scholar · View at Scopus - V. Rouss, W. Charon, and G. Cirrincione, “Neural model of the dynamic behaviour of a non-linear mechanical system,”
*Mechanical Systems and Signal Processing*, vol. 23, no. 4, pp. 1145–1159, 2009. View at Publisher · View at Google Scholar · View at Scopus - B. Subudhi and D. Jena, “A differential evolution based neural network approach to nonlinear system identification,”
*Applied Soft Computing Journal*, vol. 11, no. 1, pp. 861–871, 2011. View at Publisher · View at Google Scholar · View at Scopus - X. Dang and Y. Tan, “RBF neural networks hysteresis modelling for piezoceramic actuator using hybrid model,”
*Mechanical Systems and Signal Processing*, vol. 21, no. 1, pp. 430–440, 2007. View at Publisher · View at Google Scholar · View at Scopus - R. Storn and K. Price, “Differential evolution—a simple and efficient heuristic for global optimization over continuous spaces,”
*Journal of Global Optimization*, vol. 11, no. 4, pp. 341–359, 1997. View at Google Scholar · View at Scopus - H. Yousefi, H. Handroos, and A. Soleymani, “Application of differential evolution in system identification of a servo-hydraulic system with a flexible load,”
*Mechatronics*, vol. 18, no. 9, pp. 513–528, 2008. View at Publisher · View at Google Scholar · View at Scopus - G. Quaranta, G. Monti, and G. C. Marano, “Parameters identification of Van der Pol-Duffing oscillators via particle swarm optimization and differential evolution,”
*Mechanical Systems and Signal Processing*, vol. 24, no. 7, pp. 2076–2095, 2010. View at Publisher · View at Google Scholar · View at Scopus - H. Tang, S. Xue, and C. Fan, “Differential evolution strategy for structural system identification,”
*Computers and Structures*, vol. 86, no. 21-22, pp. 2004–2012, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. Biswas, S. Das, A. Abraham, and S. Dasgupta, “Design of fractional-order PI
*λ*D*μ*controllers with an improved differential evolution,”*Engineering Applications of Artificial Intelligence*, vol. 22, no. 2, pp. 343–350, 2009. View at Publisher · View at Google Scholar · View at Scopus - F. Neri and V. Tirronen, “Recent advances in differential evolution: a survey and experimental analysis,”
*Artificial Intelligence Review*, vol. 33, no. 1-2, pp. 61–106, 2010. View at Publisher · View at Google Scholar · View at Scopus - J. Brest, B. Boškovć, S. Greiner, V. Žumer, and M. S. Maučec, “Performance comparison of self-adaptive and adaptive differential evolution algorithms,”
*Soft Computing*, vol. 11, no. 7, pp. 617–629, 2007. View at Publisher · View at Google Scholar · View at Scopus - Z. Zhu, X. Zhou, Q. Liu, and S. Zhao, “Multi-objective optimum design of fast tool servo based on improved differential evolution algorithm,”
*Journal of Mechanical Science and Technology*, vol. 25, no. 12, pp. 3141–3149, 2011. View at Publisher · View at Google Scholar · View at Scopus - A. K. Qin and P. N. Suganthan, “Self-adaptive differential evolution algorithm for numerical optimization,” in
*Proceedings of the IEEE Congress on Evolutionary Computation (IEEE CEC '05)*, pp. 1785–1791, September 2005. View at Scopus - W. Gong, Á. Fialho, Z. Cai, and H. Li, “Adaptive strategy selection in differential evolution for numerical optimization: an empirical study,”
*Information Sciences*, vol. 181, no. 24, pp. 5364–5386, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. Brest, S. Greiner, B. Bošković, M. Mernik, and V. Zumer, “Self-adapting control parameters in differential evolution: a comparative study on numerical benchmark problems,”
*IEEE Transactions on Evolutionary Computation*, vol. 10, no. 6, pp. 646–657, 2006. View at Publisher · View at Google Scholar · View at Scopus - V. Tirronen, F. Neri, and T. Rossi, “Enhancing differential evolution frameworks by scale factor local search—part I,” in
*Proceedings of the IEEE Congress on Evolutionary Computation (CEC '09)*, pp. 94–101, Trondheim, Norway, May 2009. View at Publisher · View at Google Scholar · View at Scopus - F. Neri, V. Tirronen, and T. Kärkkäinen, “Enhancing differential evolution frameworks by scale factor local search—part II,” in
*Proceedings of the IEEE Congress on Evolutionary Computation (CEC '09)*, pp. 118–125, Trondheim, Norway, May 2009. View at Publisher · View at Google Scholar · View at Scopus - A. K. Qin, V. L. Huang, and P. N. Suganthan, “Differential evolution algorithm with strategy adaptation for global numerical optimization,”
*IEEE Transactions on Evolutionary Computation*, vol. 13, no. 2, pp. 398–417, 2009. View at Publisher · View at Google Scholar · View at Scopus - N. S. Teng, J. Teo, and M. H. A. Hijazi, “Self-adaptive population sizing for a tune-free differential evolution,”
*Soft Computing*, vol. 13, no. 7, pp. 709–724, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. Brest and M. Sepesy Maučec, “Population size reduction for the differential evolution algorithm,”
*Applied Intelligence*, vol. 29, no. 3, pp. 228–247, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. Brest and M. S. Maučec, “Self-adaptive differential evolution algorithm using population size reduction and three strategies,”
*Soft Computing*, vol. 15, no. 11, pp. 2157–2174, 2011. View at Publisher · View at Google Scholar · View at Scopus - M. Ali and M. Pant, “Improving the performance of differential evolution algorithm using Cauchy mutation,”
*Soft Computing*, vol. 15, no. 5, pp. 991–1007, 2011. View at Publisher · View at Google Scholar · View at Scopus - R. S. Rahnamayan, H. R. Tizhoosh, and M. M. A. Salama, “Opposition-based differential evolution,”
*IEEE Transactions on Evolutionary Computation*, vol. 12, no. 1, pp. 64–79, 2008. View at Publisher · View at Google Scholar · View at Scopus - Y. Wang, Z. Cai, and Q. Zhang, “Enhancing the search ability of differential evolution through orthogonal crossover,”
*Information Sciences*, vol. 185, no. 1, pp. 153–177, 2012. View at Publisher · View at Google Scholar - M. M. Ali, C. Khompatraporn, and Z. B. Zabinsky, “A numerical evaluation of several stochastic algorithms on selected continuous global optimization test problems,”
*Journal of Global Optimization*, vol. 31, no. 4, pp. 635–672, 2005. View at Publisher · View at Google Scholar · View at Scopus