Path planning of mechanical polishing process for freeform surface with a small polishing tool

Products with freeform surface are widely applied in industries, and the surface quality plays an important role in order to fulfill the targeted functions. As polishing path of small polishing tool affects the polishing removal function considerably, it is highly necessary to study the polishing path of freeform surface for obtaining good polishing efficiency and well-proportioned surface quality. By combining the Preston polishing removal function, the material removal model of small polishing tool under the control of constant polishing force and pressure is established. Based on this model, the material removal functions of scan line, Archimedean spiral, and Hilbert fractal polishing path are derived. The simulation results show that the Hilbert fractal polishing path has the best comprehensive performance. By using the projection relation of differential geometry, the optimal path generation algorithm of the Bézier surface based on Hilbert fractal polishing path is established. The polishing experiments are conducted on a self-developed polishing machine which is based on a parallel manipulator. The experimental results demonstrate that the surface roughness is improved from level 9 to level 11.


Background
Product structure can be optimized by using freeform surfaces, which opens the door for solutions with improved performance, reduced complexity, lower mass, and smaller size. However, due to the geometric particularity, freeform surfaces face more challenges and difficulties in precision manufacturing. Polishing is usually one of the final processing steps of precision manufacturing, and the results directly affect the appearance and longevity of parts [1]. It is fundamentally different from other precision manufacturing technologies. Removal of polishing does not only depend on the position of the tool orthogonal to the workpiece, as for grinding and cutting processes but also proportional to the product of local pressure and relative-speed between tool and workpiece and the dwell time. In recent years, there were efforts to develop versatile polishing processes in order to achieve high accuracy. The current polishing methods mainly include electrochemical polishing [2], magnetorheological fluid polishing [3,4], plasma polishing [5,6], ultrasonic polishing [7,8] and *Correspondence: libing.sgs@hit.edu.cn 1 Shenzhen Graduate School, Harbin Institute of Technology, Shenzhen University Town, Xili, 518055 Shenzhen, China Full list of author information is available at the end of the article computer-controlled mechanical polishing [9,10]. Nowadays, a lot of freeform components are still produced by final manual polishing. It not only heavily relies on the know-how and experience of technicians but also needs much attention for processing and testing. To achieve a given level of precision with high efficiency and reliability, process automation is clearly the way forward. As the computer-controlled mechanical polishing has high efficiency and can be controlled easily [11], it is the focus of this study. The mechanical polishing is a statistical 'rubbing' process that the microscopic loose-abrasive particles in the polishing liquid which is driven by high-speed rotational polishing tool could produce friction with the part surface. Protruding portions of the surface are removed to meet the roughness requirement. However, due to the various factors, the polishing process heavily relies on trial and experience, which leads to a slow development in the mechanical polishing process.
Compared with the large polishing tool, the small polishing tool used in this study has many advantages [12]. The small tool can follow the freeform surface with relatively large curvature, while large tool cannot polish the freeform surface with curvature smaller than its radius. The small tool can redress the error of local surface, while a large tool may polish the nearby surface when it polishes a local surface. The small can also operate at greater pressure and velocity than a large tool; therefore, it can remove the material in a rapid manner. In the case of manufacturing one single piece, small polishing tool has a higher polishing efficiency.
As the polishing surface is highly nondeterministic, one of the reasons is that the polishing path of the tool affects the removal considerably. A basic requirement for polishing paths is that the surface can be completely and uniformly covered during a polishing cycle. Therefore, the study of removal function of different polishing path has a very important significance. A lots of polishing paths in mechanical polishing process can be found in literatures [13][14][15][16][17][18]. The ultimate goal of the research about different polishing removal function is to select an optimal polishing path so as to avoid any under-polishing and over-polishing phenomenon and ensure the uniformity of the polishing process. Meanwhile, it can also present a clear understanding of the various factors that influence the polishing quality. The core issue of this research is to study the polishing removal function and the polishing performance in a given surface shape and a specified path, so that it can facilitate the subsequent polishing path control.
The rest of the paper is organized as follows. In the 'Theoretical background' section, some basic polishing assumptions and fundamental theories are introduced. In the 'Methods' section, the scan line path, Archimedean spiral path, and Hilbert fractal path are described and analyzed in details, while the modeling method of polishing path from plane surface to freeform surface is also given. The simulation and experiment are conducted, and the results are depicted in the 'Results and discussion' section. Conclusions are drawn in the last section.

Removal rate of polishing
Removal rate of polishing, which depicts the properties of part surface's removal amount at a local polishing position during every unit of time, is an important part of study about polishing process [19,20]. Many factors affect the polishing removal rate, such as polishing tool and polishing environment [21][22][23]. For the polishing tool, the main factors are material property, shape, and surface roughness. While for the polishing part, the results are relevant to material property and friction coefficient between the part and the polishing tool. Additionally, for the polishing environment, the main factors include the particle size and viscosity of polishing paste and polishing liquid, the temperature and pressure during polishing and the chemical correction of part in the polishing liquid.
Preston presented theoretical basis for the prediction of material removal for mechanical polishing [18]. Although there are many other polishing formulations in polishing literatures, most of them are based on the Preston equation. Tseng and Wang [24] proposed a modification equation based on the Preston equation in chemicalmechanical polishing; Nanz [25] considered the bending of pad and flow of slurry into the Preston equation to obtain a new polishing formula. The formula of Preston polishing removal function is displayed as below: where f ( q, t) is the material removal in unit time; q is the position vector of polishing point; t is the time; k is the Preston coefficient which is related to the material, polishing liquid, and environment; p( q, t) is the pressure at the polishing point; and v( q, t) is the instantaneous polishing velocity at the polishing point. Actually, the problem of polishing removal rate is a very complex nonlinear problem. Preston polishing removal function is a linearization assumption, but it is approximately valid in many practical polishing experiments.

Removal amount of polishing path
In terms of mechanical polishing, the polishing velocity and polishing pressure can be controlled by computer numerical control (CNC) machine precisely. Therefore, during the polishing process, the Preston coefficient k is regarded as constant. Moreover, the polishing tool is usually made of relatively soft materials. At the micro level, the soft tool is adaptive to the shape of surface. It can keep the polishing tool and word-piece contact with each other throughout the polishing process. Thus, the polishing pressure applied by small polishing tool on freeform surface is assumed to be uniform; that is to say, the pressure of contact surface is equal everywhere at the same time.
At time τ and point q in the part coordinate system, the polishing velocity and pressure are set as v( q, τ ) and n(τ ), respectively. When the radius of polishing tool is R, the polishing removal amount is defined as follow: The polishing pressure during the polishing process can be divided into the constant polishing pressure n s and the time-varying polishing pressure n v (τ ), namely n(τ ) = n s + n v (τ ). However, in order to simplify the subsequent discussion, we assume that the polishing force remains constant during the polishing process. Then, the polishing removal function can be rewritten as: Equation 3 is a general polishing removal function. Given different polishing paths, the polishing velocity in Equation 3 will change correspondingly, and this will result in different integration results. Figure 1 is the diagram of polishing model at a fixed point. When the polishing tool is stationary and rotates with angular velocity ω t , the polishing removal function of the point with radius r to the center of polishing tool is described as:

Removal function at a fixed point
The Equation 4 is a conical surface as shown in Figure 2. The velocity of polishing tool center equals zero which results in zero removal amount at this point.
In practice, the polishing velocity is the sum of feed velocity and rotation velocity; that is to say, v = v t + r · ω t . However, the angular velocity ω t of the polishing tool is very large which can reach about thousands of revolutions per minute. As a result, v t is very small and r · ω t v t , so it can be considered that v ≈ r · ω t .

Methods
Tool paths are one of the key factors for the automatic mechanical surface polishing. Three polishing paths are covered in this study: scan line path, Archimedean spiral path, and Hilbert fractal path. How they affect the polishing results are investigated here.

Scan line polishing path
The scan line path, a conventional machining path, is frequently used in CNC machining. In the scan line polishing  process, the polishing tool offsets an interval at the boundary and continues the machining until the processing is done.

Removal function of single line polishing path
In the modeling of single line polishing path along the Xaxis, the removal function Equation 3 at point (x p , y p ) on plane can be rewritten as: where The coefficient term of Equation 5 is proportional to the removal function, and it does not impact the shape of the removal function curve. Thus, let I represent the coefficient term: Figure 3 is the removal function of single line polishing path in the XZ plane. The maximum removal value is not at the polishing tool center, although the center point has a longest valid polishing length.

Removal function of scan line polishing path
Let Z y 0 (x, y) = Z s (x, y − y 0 ) denote removal function of the polishing line y = y 0 . Calculating the sum of all the polishing lines with interval L, the removal function of scan line is:

Archimedean spiral polishing path
Archimedean spiral is an arithmetic spiral, which is one of the most useful machining paths. Due to the rotational motion of workpiece, with no need of reversing, the polishing machine can achieve smooth and steady motion. Figure 5 is the model of Archimedean spiral polishing path. Let L be the constant separation distance and v t be the feed velocity of the polishing tool. The Archimedean spiral formula in the polar coordinate system is ρ = (L · θ)/2π. The feed velocity of polishing tool is v t = ρ 2 + ρ 2 ω 2 . Therefore, the angular velocity of the Archimedean spiral ω =θ can be obtained:

Modeling of Archimedean spiral polishing path
With the initial condition θ(0) = 0, the integral for Equation 9 is: where arcsinh(θ) = ln θ + √ 1 + θ 2 is the inverse hyperbolic function. There is no analytical expression for θ in Equation 10. The value of θ should be calculated by numerical method and then the value of ρ can be obtained. Figure 6 is the model of Archimedean spiral polishing path. R is the radius of the polishing tool and r is the distance between point (ρ 0 , θ 0 ). (ρ, θ) denotes the polishing tool center. According to the cosine law r(ρ, θ, ρ 0 , θ 0 ) = where the piecewise function σ R is defined as follows: For a given velocity v t and a separation distance L, the variables ρ and θ are the function of t. With Equations 3 Figure 6 Polishing model of Archimedean spiral path. and 11, the integral of Archimedean spiral path can be obtained as: The polishing process of Archimedean spiral is the uniform motion polishing process, the natural coordinate is s = v t · t. Thus, Equations 10 and 13 can be rewritten as: In the natural coordinate system, Archimedean spiral path is expressed by natural coordinate s rather than time t.

Removal function of Archimedean spiral polishing path
The expression of natural coordinate s can be obtained from Equation 14: Set the radius of the polishing tool R as a reference unit 1. The polar coordinate of polishing tool center is (ρ, θ).
Using numerical method and discrete method, the removal function of Equation 17 can be rewritten approximately as: The value of Equation 18 at polishing point (ρ 0 , θ 0 ) can be calculated by Newton's method.

Hilbert fractal polishing path
The scan line and Archimedean spiral polishing path discussed above are very directional. Therefore, the removal functions of them have strong anisotropy property. The anisotropy property destroys the polishing isotropy of the polished surface. In order to obtain high polishing quality, the removal function of designed polishing path with isotropy property is expected to apply into freeform surface polishing task.

Fractal curve
The fractal curve refers to the curve with fractal dimension. The local simple structure of fractal curve is suitable for CNC programming and controlling, while the global randomness property shows a global isotropy characteristic to remove the streaks of polished freeform surface and improves the polishing quality.
Fractal curves include snowflake fractal curve, triangle fractal curve, Hilbert fractal curve, and so on, but not every fractal curve is qualified for surface polishing. The qualified fractal curve needs to satisfy the following conditions: surface ergodic property, easy motion control for CNC, and strong isotropy property.
The surface ergodic property requires the fractal path with dimension greater than one, and the fractal curve can cover the polished surface when the order of the fractal curve is large enough. In mathematics, it means that the closure of the fractal curve is homeomorphism to two-dimensional plane.
The surface ergodic property requires that fractal curve can cover the whole higher-order polished surface. The snowflake fractal curve does not fulfill the surface ergodic property. For example, the cubic Koch snowflake cannot cover the two-dimensional surface. On the other hand, the easy motion control property requires the simple curve. The isometric segment is a good choice for this requirement because of easy planning, programming, and CNC machining. The snowflake fractal curve does not have well isotropic property globally, and it has a strong global orderliness with the high orders. Hence, it is not suitable for polishing purpose. Fulfilling the three requirements for freeform surface polishing, the Hilbert fractal curve is selected for freeform surface polishing here.
As shown in Figure 7, the Hilbert fractal path is organized with segments and traverses the whole twodimensional plane. The segments of the Hilbert curve are parallel to X axis or Y axis. They can be easily controlled by CNC machine. Also, Hilbert fractal path has a good isotropy and no orderliness in global. The basic unit of the Hilbert curve is segment. Therefore, Equation 3 with segment path can be applied for Hilbert polishing path. The segments of Hilbert curve are all parallel to X axis or Y axis.

Modeling of Hilbert fractal polishing path
For a given segment along X axis direction, the motion of polishing tool center is uniform linear motion in segment C x (x s , y) − (x e , y) . The removal function is nonnegative. Hence, changing the integral variable from dτ to dx needs to use absolute function to guarantee the nonnegative value of Equation 3. The removal function on point (x p , y p ) generated by the polishing segment C x can be obtained: Similarly, C y [(x, y s ) − (x, y e )] denotes the uniform linear segment of polishing tool center along Y axis direction, the removal function on point (x p , y p ) generated by the polishing segment C y can be obtained as: The σ R in Equations 19 and 20 has a mutation at the point where |r| = R. Therefore, the mutation is not convenient for solving the integration of Equations 19 and 20. x − and x + are the x-coordinate of intersection points between the line and the circle. Let X min = min(x s , x e ) and The valid lower and upper limits x s ≤ x e of integral function along X axis direction the polishing segment can be defined as: x e = ⎧ ⎪ ⎪ ⎨ ⎪ ⎪ ⎩ x + , X min ≤ x + ≤ X max , and |y − y p | ≤ R X max , x − ≤ X max ≤ x + , and |y − y p | ≤ R not defined , other Similarly, let Y min = min(y s , y e ) and Y max = max(y s , y e ); let y − = y p − R 2 − (x − x p ) 2 and y + = y p + lower and upper limits y s ≤ y e of integral function along Y axis direction polishing segment can be defined as: The integrands of Equations 19 and 20 are given, respectively: f x p ,y p (x, y) = δx δx 2 + δy 2 + δy 2 ln δx + δx 2 + δy 2 2 (25) g x p ,y p (x, y) = δy δx 2 + δy 2 + δx 2 ln δy + δx 2 + δy 2 2 (26) where δx = x − x p , δy = y − y p . When x in f x p ,y p (x, y) is not defined, f x p ,y p (x, y) = 0; when y in g x p ,y p (x, y) is not defined, g x p ,y p (x, y) = 0. With new defined lower and upper limits, the absolute function can be removed from the integration. The single uniform linear polishing segment can be rewritten as: Z Y x p , y p , x, y s , y e = I · g x p ,y p x, y s − g x p ,y p x, y s where x s , x e , y s , and y e are defined in Equations 21, 22, 23, and 24. Applying uniform linear motion segment, the removal functions of Equations 27 and 28 with the Hilbert fractal curve are the segments set H, the removal function of Hilbert fractal curve at point (x p , y p ) is the sum of all the segments removal functions:

Freeform surface polishing modeling
The previous polishing paths are all discussed in the planer plane. Referring to differential geometry method, the planer plane can be mapped on a freeform surface. The following discusses the general mapping method.

Natural coordinate expression of freeform surface path
The freeform curve expressed by natural coordinate s can be directly used in path planning. Freeform surface refers to a smooth two-dimensional surface with finite or infinite degrees of freedom. Because the order of freeform surface is unknown, it can hardly use an unified expression to describe the surface exactly. The freeform surface is expressed by linear combining of several basic functions to approach the freeform surface itself. The polynomial is one of the most useful basic functions due to its simpleness and easy calculation. According to the approximation theory, piecewise polynomial surfaces can approximate a smooth freeform surface in arbitrary precision.
Let C uv be polishing path in parameter u-v plane and t be the parameter of path C uv . Mapping the plane to the freeform surface obtains the project path C f . The space path C f can be expressed as follows: In application, x(t) and y(t) are selected as the parameter variables of u(t) and v(t). The path C f on freeform surface is described as z = z (x(t), y(t)).
The natural coordinate parameter s which is the length of curve is used to replace the parameter variable t to describe the space path C f . Referring to formula ds = dx 2 + dy 2 + dz 2 , natural coordinate parameter s can be written as: where z x = dz/dx. With Equation 31, the interpolation algorithm of the polishing path on freeform surface can be easily designed.

Modeling of polishing path on the Bézier freeform surface
The traditional surface fitting methods include the Hermite method, Bézier method, B-spline method, and  NURBS method. The Bézier surface can approximate any freeform surface and has fewer parameters which make it easier to use in practice. Thus, the Bézier method is studied in this paper for freeform surface polishing.
Bézier method is polynomial-based approximation algorithm. The degree n of the polynomial is related to the number of the knots. Polynomial can approximate a smooth function in arbitrary precision. Hence, the polynomial-based Bézier surface can approximate the smooth freeform surface as well in arbitrary precision.
Given m rows and n columns Bézier control points P ij = , P x mn = P x ij , P y mn = P y ij , and P z mn = P z ij . Thus, the Bézier surface can be expressed as: where u ∈[ 0, 1] , v ∈[ 0, 1]. The Bernstein basic function B i n (α) is defined as:   The surface is evenly divided into m × n matrix points. Therefore, P x ij = X s + i−1 m−1 · (X e − X s ) and P y ij = Y s + j−1 n−1 · (Y e − Y s ). Bernstein polynomial is endowed with a binomial expansion n i=0 B i n (α) = 1. The identical relation for the x(u, v) and y (u, v) in Equation 32 is: The Bézier surface can be simplified as follows: Thus, z x and z y in Equation 31 are expressed as: Figure 15 Removal difference ratio function curve of Archimedean spiral.

The first order derivative of Bernstein function is
can be calculated by a series of simple algebraic operation. With the natural coordinate parameter s, the interpolation algorithm of freeform surface curve can be directly obtained.

Results and discussion
Different parameters for the same type of polishing path may cause different quality. The following will discuss removal function of the different polishing paths and optimize the parameters of the polishing path. Figure 9 shows the removal function curves of the scan line polishing with different spaces L. Ignoring the boundary effect, the points on the surface have no more than two times of repeatedly polishing along the feed direction when R < L < 2R. However, when L < R, there are several times of repeatedly polishing of the surface. If the equivalent coefficient and the ratio between the rotational speed and the feed speed of polishing tool is fixed, the polishing  removal amount and the uniformity will increase with the decrease of the polishing space L < R. As the main role of polishing is to remove the microscopic irregularities of the surface, the smaller removal amounts the better the surface uniformity based on the condition that the accuracy of polishing is ensured.

Scan line polishing path
Set the maximum removal amount and the minimum removal amount at a specific polishing point during a polishing cycle as Z max and Z min , respectively. The difference of the removal amount Z = Z max − Z min which is varied with polishing space is studied. By numerical calculation, the curve of the Z changed with polishing space can be obtained as shown in Figure 10. It displays that the removal amount increases with the decrease of the polishing space according to a zigzag trend. By only considering the difference of removal amount and eliminating the effect of the removal amount, it can be seen that the smaller polishing space has a better polishing results. Of  course, this value needs to be selected from the trough point on the zigzag curve. The uniformity of the polishing cannot be fully seen if only the difference of the polishing removal amount Z is discussed, because different polishing space has different difference maximum polishing amount Z max . Thus, the difference ratio coefficient R Z = Z/Z max is introduced to eliminate the effect of Z max on Z. Figure 11 shows that the ratio coefficient R Z varies with the polishing space. Although the trend of the results is similar to the curve line of the polishing removal amount during a cycle, Figure 11 can better reflect the polishing amount uniformity of the polishing path.
In order to ensure the polishing efficiency and avoid excessive repeatedly polishing of the surface, the polishing space should be set within the interval of [ R, 2R]. From Figures 10 and 11, it can be found that the optimal value appear in the vicinity of 1.9R. The planar graph and the three-dimensional graph of the removal function at this value is shown in Figure 9b and Figure 12, respectively. Figure 13 shows that the curves of the polishing removal function vary with ρ 0 when θ 0 = 0. When L > 2R, the spiral polishing path do not interfere with each other. Thus, the curve of removal function is an independent waveform. Similar to the scan line polishing path, the polishing amount increased with the decrease of the polishing space. What different to the scan line polishing path is that the Archimedean spiral polishing path displays quasi-periodicity characteristic. The periodicity becomes increasingly stronger with the distance from the center of Archimedean spiral path and the removal amount near the center of the Archimedean spiral path varies considerably.

Archimedean spiral polishing path
In order to depict the uniformity of the Archimedean spiral polishing, Z max is set as the maximum removal amount and Z min is set as the minimum removal amount. The difference of the removal amount is represented as Z = Z max − Z min . By numerical calculation, Z changing with the space can be derived as shown in Figure 14. Different from the scan line polishing path, the smaller polishing space does not display a better polishing results.
The difference ratio coefficient R Z = Z/Z max is introduced, and the curve line that it varies with the polishing space is shown in Figure 15. It can be seen from the figure that with the decrease of the polishing space, the uniformity of the surface becomes better and better.
In practical, the polishing space of the Archimedean spiral cannot tend to zero. In order to improve the polishing efficiency, the polishing space usually set within the interval of [ R, 2R]. Figure 16 demonstrates the threedimensional graph of removal function at the optimal polishing space L = 1.3R which can be obtained from  the Figure 14. Although this is the optimized result, the removed amount near the center of Archimedean spiral is still highly asymmetric.

Hilbert fractal polishing path
The quintic Hilbert fractal curve is projected on a plane with the range [ −2R, 2R] ×[ −2R, 2R]. In order to ignore the boundary effect, only the results in the area within the range [ −R, R] ×[ −R, R] are studied. The polishing space of small segment line in the quintic Hilbert fractal curve is 0.125R. Figure 17 is the polishing removal function curve of the quintic Hilbert fractal curve along the X-axis. It can be seen from the figure that the removal amount is symmetrical to the Y-axis. In addition, the periodic of the polishing removal function is not obvious, which means it satisfies the overall randomness requirement of polishing.
The overall randomness characteristic of the fractal curve polishing path can be seen from the threedimensional graph of the removal function which is shown in Figure 18. This kind of characteristic is very suitable for surface polishing, and it can make the polishing results achieve good performance.

Experiment
The features for the three kinds of polishing path are listed in Table 1. Thus, the polishing path selected in this study is the Hilbert fractal polishing path which has geometric ergodicity and affine projection invariance characteristics.
In order to verify the performance of the presented polishing path, a part with Bézier freeform surface is selected. First, the discrete point positions and their normal vectors of the surface are saved. Next, by importing them into the fractal path generation module developed with Equation 31, the Hilbert fractal path is obtained. Then, the data file of the Hilbert fractal path is imported into the polishing simulation software of the self-developed polishing machine and the path on the freeform surface is obtained, which is shown in Figure 19.
Based on the freeform surface polishing requirement, a self-developed automatic polishing machine is proposed, which mainly includes a five-DOF parallel manipulator, a single DOF rotary table, and a linear motion platform with a force feedback system. It has been described extensively elsewhere in the literatures [26][27][28]. The parallel manipulator is a closed mechanical structure, and it has high stiffness, low inertia, and high dynamic performance. In the polishing process, the high stiffness helps to reduce the vibration generated by high speed rotation of polishing spindle, and the low inertia feature allows high polishing speed. A redundant translation is provided by a linear motion actuator on the moving platform to realize a constant polishing force. The rotary table is used to enlarge the polishing workspace of the automatic polishing machine. Moreover, with the help of rotary table, the path generation module could generate a singularfree and optimized dexterity polishing path for the parallel manipulator. Wool wheel polishing tool and polishing paste are used to polish the freeform surface part. During the experiment, the objective constant control force is set as 6N, and the polishing process is shown in Figure 20.
The roughness of the part surface is detected by the Taylor Hobson roughmeter. The original surface is only dealt with polishing cloth, and the roughness S a equals 173.5 nm and the maximum height of the profile S z equals 383.9 nm. This belongs to level 9 precision. After polishing, the roughness S a decrease to 43.1 nm and the maximum height of the profile S z decrease to 119.4 nm. This belongs to level 11 precision. The measuring roughness of the freeform surface after polishing is demonstrated in Figure 21.

Conclusions
Path planning is one of the key issues of the polishing process for freeform surface. Using the classical Preston polishing function, the mathematical model of small polishing tool is derived. With the assumption of a constant polishing pressure, the removal functions of scan line polishing path, Archimedean spiral polishing path, and Hilbert fractal polishing path are derived. The modeling method for polishing paths on a planar surface projected on a freeform surface is also discussed.
As different parameters for the same type polishing path may cause different polishing quality, the optimize parameters of polishing path are analyzed. The optimal polishing space for the scan line path and the Archimedean spiral polishing path is 1.9R and 1.3R, respectively. By comparing the removal functions of different polishing path, the reasonable polishing path and parameters are chosen. In order to verify the effectiveness of the selected Hilbert fractal polishing path, a polishing experiment is conducted on a self-developed mechanical polishing machine and the sequence of polishing process can be automatically scheduled. The results demonstrate that the roughness of freeform surface is improved from level 9 to level 11.