Next Article in Journal
Impact of Magnetic Field Environment on the EDM Performance of Al-SiC Metal Matrix Composite
Next Article in Special Issue
High-Accuracy Surface Topography Manufacturing for Continuous Phase Plates Using an Atmospheric Pressure Plasma Jet
Previous Article in Journal
A Decade of Organs-on-a-Chip Emulating Human Physiology at the Microscale: A Critical Status Report on Progress in Toxicology and Pharmacology
Previous Article in Special Issue
A New Method for Precision Measurement of Wall-Thickness of Thin-Walled Spherical Shell Parts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Elementary Approximation of Dwell Time Algorithm for Ultra-Precision Computer-Controlled Optical Surfacing

1
School of Mechanical Engineering, Dalian University of Technology, Dalian 116023, China
2
Institute of Machinery Manufacturing Technology, China Academy of Engineering Physics, Mianyang 621900, China
*
Author to whom correspondence should be addressed.
Micromachines 2021, 12(5), 471; https://doi.org/10.3390/mi12050471
Submission received: 23 March 2021 / Revised: 13 April 2021 / Accepted: 16 April 2021 / Published: 21 April 2021
(This article belongs to the Special Issue Frontiers in Ultra-Precision Machining)

Abstract

:
The dwell time algorithm is one of the key technologies that determines the accuracy of a workpiece in the field of ultra-precision computer-controlled optical surfacing. Existing algorithms mainly consider meticulous mathematics theory and high convergence rates, making the computation process more uneven, and the flatness cannot be further improved. In this paper, a reasonable elementary approximation algorithm of dwell time is proposed on the basis of the theoretical requirement of a removal function in the subaperture polishing and single-peak rotational symmetry character of its practical distribution. Then, the algorithm is well discussed with theoretical analysis and numerical simulation in cases of one-dimension and two-dimensions. In contrast to conventional dwell time algorithms, this proposed algorithm transforms superposition and coupling features of the deconvolution problem into an elementary approximation issue of function value. Compared with the conventional methods, it has obvious advantages for improving calculation efficiency and flatness, and is of great significance for the efficient computation of large-aperture optical polishing. The flatness of φ150 mm and φ100 mm workpieces have achieved PVr150 = 0.028 λ and PVcr100 = 0.014 λ respectively.

1. Introduction

With the rapid increasing requirements for the fabrication of high-precision optical elements in modern optical systems, several advanced deterministic optical surfacing technologies have been developed over the past decades, such as ultra-precision computer controlled optical surfacing (CCOS), magnetorheological finishing (MRF), ion-beam figuring (IBF), bonnet polishing (BP) [1,2]. These achieve precision material removal on certain workpiece areas by accurately controlling the dwell time on the elaborately predesigned polishing path. Therefore, the dwell time algorithm is one of the key elements in modern advanced deterministic optical surfacing technologies. Among them, most of the removal functions have the characteristics of rotational symmetry, but how to use this feature to develop a high-efficiency, high-precision dwell time algorithm is the current research focus.
In deterministic optical finishing technologies, the amount of material removal can be expressed as discrete two-dimensional convolution operations of the dwell time and removal functions. The dwell time algorithm is used to solve the deconvolution process and to ensure that the calculated dwell time not only meets the performance of machine tools, but also has a high surface error convergence efficiency. Various dwell time algorithms have earlier been proposed, including the Fourier transform method, the numerical iteration method, the matrix equation method. Ronald Aspden et al. [3] studied the polar and rectangular coordinates of the process in CCOS, and discussed the variation of the removal function with the radius of the workpiece in the gyrosymmetric correction process. Jones et al. [4] proposed an iterative method for solving the dwell time function, studied the relationship between the flatness convergence efficiency and the removal function, and pointed out that only the symmetric central single-peak removal function could converge. The convergence accuracy of these methods is not high enough, and now they are seldom applied. Carnal et al. [5] introduced the linear equation method and solved the dwell time by adopting the LSQR method. Drueding et al. [6] proposed a series expansion solution. Waluschka et al. [7] presented a one-dimensional dwell time function algorithm for cylindrical workpieces based on a graphic method, and Shanbhag et al. [8] proposed an algorithm based on wavelet analysis. Zheng et al. [9] proposed a damped iterative method for solving the CCOS dwell time function. Zhou et al. [10] used the TSVD method to solve the linear equation model, which entailed further research. Wu et al. [11] proposed a solution based on discretized matrix equations using LSQR. Jiao et al. [12] and Jiang et al. [13] improved the traditional L–R algorithm. Taking into account the scanning path of the polishing tool and the tilt angle of the workpiece, Guo [14] proposed a dwell time algorithm to achieve rapid convergence of the accuracy of the optical mold. Pan et al. [15] proposed an improved dwell time calculation algorithm to optimize tool path planning in optical figuring. Li et al. [16] developed a positive dwell time algorithm with minimum equal extra material removal to consider the machine dynamics limitations. Li and Zhou [17] gave a solution algorithm of dwell time in a slope-based figuring model. Wang et al. [18] provided a quantitative study on the performances of dwell time algorithms in ion-beam figuring. Han et al. [19,20] proposed a Gaussian mixture model to model experimental removal functions and provided the dwell time algorithm according to the dynamic characteristics of the machine tool. These methods are mainly based on matrix equations, and the computational efficiency might be much lower especially for large-aperture optical elements, so the solution might not be smooth enough.
The existing dwell time algorithms are conducted mainly based on a meticulous mathematical theory and designed to pursue high convergence rate. Nonetheless, those methods do not adequately consider the distribution characteristics of the removal function and rarely incorporate the speed-smoothing issues that are closely related to convergence efficiency and machine tool motion implementation. Actually, the convergence rate of flatness is only between 1.1X and 1.3X [21] for most CNC machines; hence, it is unnecessary to pursue high convergence rates excessively.
In this paper, an elementary approximation method for solving the dwell time algorithm based on the symmetrical distribution of single-peak rotation of removal function is proposed. The proposed method has the characteristics of clear physical meaning and was verified by simulation and experiments. By using triangular approximation of the removal function, the initial surface shape is discretized reasonably and the approximate solution of dwell time is obtained. In this paper, the performance of the algorithm is verified through residence-time mathematical modeling, accuracy analysis, simulation and experimental research. The results showed that it performs well in the profile for smoothness and convergence efficiency.

2. Dwell Time Algorithm Model

2.1. Approximation Treatment of Removal Function

According to the measured residual error, the surface is polished to achieve a theoretical profile. To eliminate the residual error, the material removal function, generated by the polishing tool in a constant time (also called the removal function) [22,23,24], and the dwell time of the polishing tool should be known first. It is generally assumed that deterministic optical surfacing technology is a linear shift-invariant system, and the mathematical model of the convolution of the dwell time and removal functions being equal to the distribution of the removal amount is generally adopted, and this is given in Equation (1) [6]:
Δ h ( x , y ) = R ( x , y ) D ( x , y )
where ∆h(x, y) is the distribution function of removal amount, R(x, y) is the removal function (also called influence function), and D(x, y) is the dwell time function.
In some deterministic optical surfacing technologies, the removal function is distributed by rotational symmetry. Suppose an ideal removal function is a two-dimensional Gaussian distribution, as shown in Figure 1.
According to the actual characteristics of its distribution [25], the removal function is approximated as shown in Figure 2.
In the one-dimensional case, the removal function can be represented by an isosceles triangle distribution with height A and bottom 2R, similar to the roof function in the one-dimensional finite element method. In the two-dimensional case, it is assumed that the kernel function can be represented by a conic distribution with a height of A and a bottom radius of R. This kind of approximation reflects the main distribution characteristics of the removal function and concentrates more than 80% of volume removal, which is completely acceptable in engineering.

2.2. Dwell Time Algorithm Model

2.2.1. One-Dimensional Analysis

The principle of an elementary approximation for a one-dimensional deconvolution is shown in Figure 3. The blue curve represents the target removal amount distribution curve H(x), and the red curve is the actual removal amount distribution curve h(x). The standard removal function has a maximum width of 2R and a height of A. The discretization distance of the nodes is L = R. For each dwell node Xi, the dwell time function is set to Di = H(Xi)/A, which eliminates superposition coupling and is only a simple elementary algebraic operation. It shows that the algorithm has the same accuracy as the trapezoidal method of the one-dimensional definite integral problem, and its approximation residual is a second-order small quantity.
When L = R, the calculated residual error is already a second-order small quantity, but it can still be seen that the actual removal curve is not very smooth, which means that the smoothness of the optical processing surface is poor. When the discretization node is doubled and the spacing L = R/2, then Di = H(Xi)/2A. The calculation principle of elementary geometric approximation for one-dimensional mesh refinement is shown in Figure 4.
Generally, let L = R/2n, and the distribution function of target removal amount is H(x) = 2n·H(x)/2n. The discretized node set {Xi} is divided into 2n groups. The node spacing in each group is L = R, and the phase difference between each group of nodes is an integer multiple of R/2n in turn. Each set of nodes after partition is decoupled according to Figure 3, and then Di = H(Xi)/2n/A. In this way, the actual processing curve gradually becomes smooth.
According to the above analysis, the basic criteria for the elementary approximation of one-dimensional deconvolution are as follows:
  • It is acceptable to use an isosceles triangle as an approximate expression of the removal function in engineering;
  • The discretization distance of the nodes should not be more than half of the width of the removal function; otherwise, the deconvolution calculation will lose the ability to approximate;
  • When the node spacing is doubled, the time weight of each node is reduced by half, so the total time remains basically unchanged. The dwell time of the subdivided nodes is not the interpolation between the original discrete nodes, but the redistribution of the dwell time. The physical meaning is that the total removal amount is constant, and the removal function is constant, so the total time is basically conserved;
  • The approximation residual of approximate solution is the same as that of definite integral trapezoid method, which is a second-order small quantity.

2.2.2. Two-Dimensional Analysis

The orthogonal grid M is divided according to the spacing R. The coordinates of each grid node are (xi, yj), and the target removal amount on the node is H (xi, yj). The plane distribution is discussed first, as shown in Figure 5.
Let the removal function be a conic distribution and the center of the conic be the origin. The expression under the rectangular coordinate system is shown as follows:
{ φ i j ( x , y ) = A R ( R ( x x i ) 2 + ( y y j ) 2 ) , ( x x i ) 2 + ( y y j ) 2 R 2 φ i j ( x , y ) = 0 , ( x x i ) 2 + ( y y j ) 2 > R 2
Where R is the radius of the circular support region of the removal function, and A is the peak value of the removal function center. The expression in polar coordinate form is shown as follows:
{ φ i j ( ρ , θ ) = A R ( R ( ρ cos θ ρ i j cos θ i j ) 2 + ( ρ sin θ ρ i j sin θ i j ) 2 ) , ( ρ cos θ ρ i j cos θ i j ) 2 + ( ρ sin θ ρ i j sin θ i j ) 2 R φ i j ( ρ , θ ) = 0 , ( ρ cos θ ρ i j cos θ i j ) 2 + ( ρ sin θ ρ i j sin θ i j ) 2 > R
On the mesh, the removal function is simplified as follows:
{ φ i ( x ) = A R ( R f a b s ( x x i ) ) , ( x x i ) 2 R 2 φ i ( x ) = 0 , ( x x i ) 2 > R 2
Or the following polar form:
{ φ j ( y ) = A R ( R f a b s ( y y j ) ) , ( y y j ) 2 R 2 φ j ( y ) = 0 , ( y y j ) 2 > R 2
Where fabs represents the function of taking absolute value.
Suppose that each node (xi, yj) of the grid is superimposed with a removal function φij of Equation (2) above, which has the same weight. Then, for the interior of the region M, according to the symmetry, only the case of the middle region R × R needs to be considered. The edge of the whole grid area is special and will not be discussed here.
For the grid area, it is customary to take the lower left corner of the grid as the origin O, so the center of the removal function is respectively located at the four corners of the grid, as shown in Figure 6. The green lines are auxiliary lines, the red lines are distances from each corner nodes, and P is an arbitrary point in the polishing area. Set the node number of O as (i, j), OP = r1, POC = θ1, CP = r2, PCB = θ2, BP = r3, PBA = θ3, AP = r4, and PAO = θ4 to give the following:
{ r i j = 0 , θ i j = 0 r i j + 1 = R , θ i j + 1 = 0 r i + 1 j + 1 = 2 R , θ i + 1 j + 1 = π 4 r i + 1 j = R , θ i + 1 j = π 2
According to Equation (3), there are the following:
{ φ i j ( r 1 , θ 1 ) = A R ( R r 1 ) φ i j + 1 ( r 2 , θ 2 ) = A R ( R r 2 ) φ i + 1 j + 1 ( r 3 , θ 3 ) = A R ( R r 3 ) φ i + 1 j ( r 4 , θ 4 ) = A R ( R r 4 )
For a discrete mesh, due to the symmetry, only one-eighth of the triangular EOD area in the grid needs to be considered. This area can be further divided into four subareas: boundary line, M1, M2, and M3. For each subregion, only the value range of any point in the region under the function of each removal function (the maximum and minimum values) can be considered for evaluating the approximation ability of the approximate solution. The points where the maximum and minimum are located are the feature points in each square. According to the symmetry, these feature points must be obtained on the symmetry axis or the boundary of the square.
Similar to the one-dimensional case, let L = R/2n, and the distribution function of target removal amount is H(x) = 2n·H(x)/2n. The discretized node set (xi, yj) is divided into 2n groups. The node spacing in each group is L = R, and the phase difference between each group of nodes is an integer multiple of R/2n in turn. Each set of nodes after partition is decoupled as above, and then there is Di = H(xi, yj)/(2nA). The basic criteria for the primary approximation of two-dimensional deconvolution are as follows:
(1)
It is acceptable to use cone as an approximate expression of the removal function in engineering.
(2)
The distance of node discretization should not be larger than the radius of the removal function support domain; otherwise, the deconvolution calculation based on this method will lose the ability to approximate.
(3)
When the node spacing is doubled, the time weight of each node is reduced by half, so the total time remains basically unchanged.
(4)
The approximation residual of the elementary geometric approximation method for two-dimensional deconvolution is completely acceptable compared with the actual polishing convergence rate.

2.3. Dwell Time Algorithm Analysis

2.3.1. Split Line Value Analysis

In the real polishing case, the workpiece surface is a two-dimensional planar. For each meridian segment or latitude segment, there are only two kernel functions. At this time, the superposition value of the removal function is always A, as shown in Figure 7.
For example, let xi + 1 > x > xi, xi + 1 = xi + R, then the value on the grid is determined by Equation (8).
f = φ i ( x ) + φ i + 1 ( x ) = A R ( R ( x x i ) ) + A R ( R ( x i + 1 x ) ) = A

2.3.2. Area M1 Value Analysis

The point P is located in the region M1 and is acted on by four removal functions as shown in Figure 8.
According to Equations (3), (6) and (7), at this time the value of any point in M1 is determined by Equation (9).
f = φ i j ( r 1 , θ 1 ) + φ i j + 1 ( r 2 , θ 2 ) + φ i + 1 j + 1 ( r 3 , θ 3 ) + φ i + 1 j ( r 4 , θ 4 ) = A R ( 4 R ( r 1 + r 2 + r 3 + r 4 ) )
The relationships of each parameter are as follows:
r 2 = ( r 1 cos θ 1 ) 2 + ( R r 1 cos θ 1 ) 2 , θ 2 = A r c c t g ( r 1 sin θ 1 R r 1 cos θ 1 ) r 3 = ( r 2 cos θ 2 ) 2 + ( R r 2 cos θ 2 ) 2 , θ 3 = A r c c t g ( r 2 sin θ 2 R r 2 cos θ 2 ) r 4 = ( r 3 cos θ 3 ) 2 + ( R r 3 cos θ 3 ) 2 , θ 4 = A r c c t g ( r 3 sin θ 3 R r 3 cos θ 3 ) ( 2 1 ) R r 1 2 2 R , 5 π 36 θ 1 π 4
According to the symmetry, when point P is located at the vertex G of EFG, that is, r3 = r4 = R, r 1 = r 2 = 2 3 R , making the sum of r1+r2+r3+r4 the maximum, then the minimum value of the function in this area is:
f min = 2 ( 1 2 3 ) A 0.9647 A
Since the sum of any two sides of a triangle is greater than the third side, it can be known that when P is located at the center of the square; that is, r 1 = r 2 = r 3 = r 4 = 2 R / 2 , the sum of r1+r2+r3+r4, is a minimum, so the maximum value in this area is
f max = 4 A R ( R 2 2 R ) 1.172 A

2.3.3. Area M2 Value Analysis

The point P is located in the region M2 and is acted on by three removal functions as shown in Figure 9.
According to Equation (3), (6) and (7), at this time, the value of any point in M2 is determined by Equation (13).
f = φ i j ( r 1 , θ 1 ) + φ i j + 1 ( r 2 , θ 2 ) + φ i + 1 j ( r 4 , θ 4 ) = A R ( 3 R ( r 1 + r 2 + r 4 ) )
The relationships of each parameter are as follows:
r 2 = ( r 1 cos θ 1 ) 2 + ( R r 1 cos θ 1 ) 2 , θ 2 = A r c c t g ( r 1 sin θ 1 R r 1 cos θ 1 ) r 3 = ( r 2 cos θ 2 ) 2 + ( R r 2 cos θ 2 ) 2 , θ 3 = A r c c t g ( r 2 sin θ 2 R r 2 cos θ 2 ) r 4 = ( r 3 cos θ 3 ) 2 + ( R r 3 cos θ 3 ) 2 , θ 4 = A r c c t g ( r 3 sin θ 3 R r 3 cos θ 3 ) 0 r 1 2 - 3 R , 0 θ 1 π 4
According to the symmetry, when point P is located at the vertex G of EFG, that is, r4 = R, r 1 = r 2 = 2 3 R , making the sum of r1+r2+r4 the maximum, then the minimum value of the function in this area is:
f min = 2 ( 1 2 3 ) A 0.9647 A
According to the symmetry, when point P is located at the vertex G of EFG, that is, r1 = 0, r2 = r4 = R, making the sum of r1 + r2 + r4 the minimum, then the maximum value in this area is:
f max = A

2.3.4. Area M3 Value Analysis

The point P is located in the region M2 and is acted by two removal functions as shown in Figure 10.
According to Equation (3), (6) and (7), at this time, the value of any point in M3 is determined by Equation (17).
f = φ i j ( r 1 , θ 1 ) + φ i j + 1 ( r 2 , θ 2 ) = A R ( 2 R ( r 1 + r 2 ) )
The relationships of each parameter are as follows:
r 2 = ( r 1 cos θ 1 ) 2 + ( R r 1 cos θ 1 ) 2 , θ 2 = A r c c t g ( r 1 sin θ 1 R r 1 cos θ 1 ) 0 r 1 2 3 R , 0 θ 1 5 π 36
According to the geometric relationship and symmetry, when point P is located at the vertex G of a curved triangle, that is, r 1 = r 2 = 2 3 R , the sum of r1 + r2 is the maximum, then the following function obtains the minimum value:
f min = 2 ( 1 2 3 ) A 0.9647 A
When point P is located on the side of the curved triangle OD, that is, r1 + r2 = R, the sum of r1 + r2 is the smallest, and the function obtains the maximum value:
f max = A

2.3.5. Numerical analysis

Figure 11 is a cloud picture of simulation calculation of an equal weight superposition distribution of removal function in the middle area. The distribution characteristics are consistent with the theoretical analysis.
According to the results of the above analysis, on the whole grid the maximum and minimum values of the weight superposition distribution of the removal function are as follows:
f min = 2 ( 1 2 3 ) A 0.9647 A f max = 4 ( 1 2 2 ) A 1.172 A
This reflects the approximation ability of the conic distribution removal function to the plane, that is, the error level based on the elementary approximation method. The convergence accuracy of the calculation is much higher than the actual polishing convergence rate.
More generally, if each removal function is weighted according to the value of the surface distribution at the center, the conical distribution removal function can approximate the general surface better. That is to say, for any initial error distribution, the dwell time function the kernel of which is a cone distribution can be approximately determined by the weight of the kernel’s center. At this time, the higher deconvolution problem can be simplified to a basic function-value calculation problem.

3. Simulations

The example of a two-dimensional deconvolution based on an elementary approximation is shown below. Three sets of simulations use the cone distribution removal function, and the removal amount distribution is a plane, sphere and arbitrary surface, respectively. The dwell time function is solved by the elementary approximation method proposed in this paper. The results are shown in Figure 12, Figure 13 and Figure 14.
Apparently, the numerical simulation results of a deconvolution calculation based on an elementary approximation are satisfactory.

4. Results and Discussion

4.1. Experiment Setup and the Parameters

To validate the effectiveness of the proposed algorithm, repetitious experimental studies on flatness figuring was carried out based on the self-developed ion-beam figuring (IBF) machine. The experiment setup is shown in Figure 15. It can process planar, spherical and aspheric parts with a maximum size of 300 mm × 300 mm, and the positioning accuracy of the linear axis is below 0.005 mm. This machine can be used for the corresponding experimental verification of the dwell time solving algorithm. The specific parameters of the polishing processing are shown in Table 1.
IBF spots were taken on the IBF machine, and the results are shown in Figure 16. The polishing result was detected by the INF300H-LP-WM interferometer made in China, with RMS repeatability of less than 0.3 nm. The removal function result of IBF, shown in Figure 16, was tested on a fixed point, and polishing time was 60 s. The peak removal rate (PRR) of influence function was 0.30614 λ/min, where λ was 658 nm, and the volume removal rate (VRR) was 0.056474 mm3/min.

4.2. Results and Discussion

4.2.1. Experiment Case 1

The workpiece used for figuring is silica flat and the diameter is φ120 mm. The scanning path is in raster with a step of 1.5 mm and spacing of 3.0 mm. Initial flatness was PV 0.259 λ, RMS 0.050 λ; Predicted flatness was PV 0.051 λ and RMS 0.004 λ, with a convergence rate of 80.0% and 92%, respectively.
Experimental result indicated that full-aperture machined flatness was PV 0.077 λ and RMS 0.013 λ with a convergence rate of 70.3% and 74.5% respectively, where the missing data in the periphery was caused by the fixture. The aperture machined flatness of 90% was PV 0.063 λ and RMS 0.012 λ with a convergence rate of 64% and 70% respectively, as shown in Figure 17 and Figure 18.

4.2.2. Experiment Case 2

The universal measuring software was employed. The PRR of influence function was 0.244 λ/min, and the VRR was 0.0098 mm3/min as shown in Figure 19.
The workpiece is fused silica with diameter of φ170 mm. Surface map is measured using the relative accuracy method. Initial flatness of φ150 mm aperture is PVφ150 = 0.048 λ, RMSφ150 = 0.008 λ, and those of φ100 mm aperture are PVφ100 = 0.030 λ, RMSφ100 = 0.007 λ respectively, as shown in Figure 20.
The scanning path was rasterized with 2 mm spacing, and the predicted machine time was 83.8 min; the predicted flatness of φ150 mm aperture was PVc150 = 0.025 λ, RMSc150 = 0.003 λ; and those of φ100 mm were PVc100 = 0.014 λ; RMSc100 = 0.002 λ, as illustrated in Figure 21.
The practical polishing time was 84 min. The surface map was also measured using the relative accuracy method. The polished flatness of the φ150 mm aperture was PVr150 = 0.028 λ, and RMSr150 = 0.005 λ, and those of φ100 mm were PVcr100 = 0.014 λ, and RMSc100 = 0.002 λ, as shown in Figure 22.
Simulated convergence rate was
η c 150 = p v o 150 p v c 150 p v o 150 × 100 % = 47.9 %
And the practical convergence rate was
η r 150 = p v r 150 p v r 150 p v r 150 × 100 % = 41.7 %
Taking the surface PV value of φ100 mm aperture as the evaluation index, the predicted convergence rate was 52.0%, and the practical convergence rate was 51.9%.

5. Conclusions

In this paper, an elementary approximation of the dwell time algorithm for single-peak rotational symmetry removal function was presented. The work showed that it is engineeringly acceptable to use the cone distribution as the removal function to approximate expressions of all the subaperture polishing with a single-peak rotational symmetry removal function, such as CCOS, IBF, or BP. The dwell time algorithm model and computation method were given. When the distance of node discretization was not more than the radius of the support domain of removal function, the coupling characteristic of the deconvolution problem could be eliminated by using the elementary approximation solution proposed in this paper. Theoretical analysis, numerical simulation and experimental results show that the proposed method had a lower calculated residual error than the initial value by one order of magnitude, and had a higher approximation ability. The flatness of φ150 mm and φ100 mm workpieces achieved PVr150 = 0.028 λ and PVcr100 = 0.014 λ, respectively.
In contrast to conventional dwell time algorithms, this work transformed the superposition and coupling features of the deconvolution problem into a simple calculation of the discretization function value. When the discrete nodes were doubled, the time weight of each node was then halved; consequently, the total time remained unchanged. The approximation ability or smoothness of the deconvolution result was greatly increased, which agreed with the engineering reality that total polishing time remains unchanged if the total removal amount and removal function are unchanged. Compared with conventional methods, the proposed algorithm has obvious advantages for improving calculation efficiency and smoothness, which is of great significance for the efficient computation of large-aperture optical polishing.
Nevertheless, the calculation accuracy of the proposed dwell time algorithm is related to the symmetry of the removal function and its approximation error. Meanwhile, the dwell time algorithm has a limited ability to correct mid-to-high frequency errors determined by the sampling characteristics of the discretization nodes and the scale of the removal function. The approximation error in the edge area of the workpiece needs further analysis.

Author Contributions

All authors have contributed to the study conception and design. Y.W. designed the experiments and wrote the paper. Y.Z. performed the experiments and analyzed the data. R.K. and F.J. supervised the whole project and gave advices on manuscript editing. All authors read and approved the final manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 52005463); the Science Challenge Project (No.TZ2016006-0504-01); High-end CNC Machine Tools and Basic Manufacturing Equipment Major National Science and Technology Project (2017ZX04022001); Laboratory of Precision Manufacturing Technology, CAEP.

Data Availability Statement

NOMENCLATURE is contained within the main text.

Acknowledgments

The authors thank the technical support from peer fellows.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

R(x, y)Removal function
D(x, y)Dwell time function
H (Xi, Yj)Target removal amount
APeak value of the removal function center
fabsFunction of taking absolute value
t(xi, yi)Dwell time at the i th path node
h(xk, yk)Desired amount of removed material at the k th figure-control node
NtTotal number of the path nodes
NhTotal number of the figure-control nodes
r(xk-xi, yk-yi)Amount of removed material at the k th figure-control node
LDiscretization distance of nodes
CCOSComputer controlled optical surfacing
MRFMagnetorheological finishing
IBFIon beam figuring
BPBonnet polishing

References

  1. Zhao, Q.Z.; Zhang, L.; Han, Y.J.; Fan, C. Polishing path generation for physical uniform coverage of the aspheric surface based on the Archimedes spiral in bonnet polishing. Proc. Inst. Mech. Eng. B J. Eng. Manuf. 2019, 233, 2251–2263. [Google Scholar] [CrossRef]
  2. Zhao, Q.Z.; Zhang, L.; Fan, C. Six-directional pseudorandom consecutive unicursal polishing path for suppressing mid-spatial frequency error and realizing consecutive uniform coverage. Appl. Opt. 2019, 58, 8529–8541. [Google Scholar] [CrossRef] [PubMed]
  3. Aspden, R.; Mcdonough, R.; Nitchie, F.R. Computer assisted optical surfacing. Appl. Opt. 1972, 11, 2739–2747. [Google Scholar] [CrossRef] [PubMed]
  4. Jones, R.A. Optimization of computer controlled polishing. Appl. Opt. 1977, 16, 218–224. [Google Scholar] [CrossRef] [PubMed]
  5. Carnal, C.L.; Egert, C.M.; Hylton, K.W. Advanced matrix-based algorithm for ion-beam milling of optical components. In Proceedings of the Current Developments in Optical Design and Optical Engineering II, SPIE, San Diego, CA, USA, 20–21 July 1992; Volume 1752, pp. 54–62. [Google Scholar]
  6. Drueding, T.W.; Fawcett, S.C.; Wilson, S.R.; Bifano, T.G. Ion beam figuring of small optical components. Opt. Eng. 1995, 34, 3565. [Google Scholar]
  7. Waluschka, E. Cylindrical optic figuring dwell time optimization. In Proceedings of the International Symposium on Optical Science and Technology, SPIE, San Diego, CA, USA, 30–31 July 2000; Volume 4138, pp. 25–32. [Google Scholar]
  8. Shanbhag, P.M.; Feinberg, M.R.; Sandri, G.; Horenstein, M.N.; Bifana, T.G. Ion-beam machining of millimeter scale optics. Appl. Opt. 2000, 39, 599–611. [Google Scholar] [CrossRef] [PubMed]
  9. Zheng, L.G.; Zhang, X.J. A novel resistance iterative algorithm for CCOS. In Proceedings of the SPIE—The International Society for Optical Engineering, San Diego, CA, USA, 14 August 2006; Volume 6288, p. 62880N-1-9. [Google Scholar]
  10. Zhou, L.; Dai, Y.F.; Xie, X.H.; Jiao, C.J.; Li, S.Y. Model and method to determine dwell time in ion beam figuring. Nanotechnol. Precis. Eng. 2007, 5, 107–112. [Google Scholar]
  11. Wu, J.F.; Lu, Z.W.; Zhang, H.X.; Wang, T.S. Dwell time algorithm in ion beam figuring. Appl. Opt. 2009, 48, 3930–3937. [Google Scholar] [CrossRef] [PubMed]
  12. Jiao, C.J.; Li, S.Y.; Xie, X. Algorithm for ion beam figuring of low-gradient mirrors. Appl. Opt. 2009, 48, 4090–4096. [Google Scholar] [CrossRef] [PubMed]
  13. Bo, J.; Zhen, M. Dwell time method based on richardson-lucy algorithm. In Proceedings of the Space Optics and Earth Imaging and Space Navigation, SPIE, Beijing, China, 4–6 June 2017; Volume 10463. [Google Scholar]
  14. Guo, J. Corrective finishing of a micro-aspheric mold made of tungsten carbide to 50 nm accuracy. Appl. Opt. 2015, 54, 5764–5770. [Google Scholar] [CrossRef] [PubMed]
  15. Pan, R.; Zhang, Y.J.; Ding, J.B.; Cao, C.; Wang, Z.Z.; Jiang, T.; Peng, Y.F. Rationality optimization of tool path spacing based on dwell time calculation algorithm. Int. J. Adv. Manuf. Technol. 2016, 84, 2055–2065. [Google Scholar] [CrossRef]
  16. Li, L.X.; Xue, D.L.; Deng, W.J.; Wang, X.; Bai, Y.; Zhang, F.; Zhang, X.J. Positive dwell time algorithm with minimum equal extra material removal in deterministic optical surfacing technology. Appl. Opt. 2017, 56, 9098–9104. [Google Scholar] [CrossRef] [PubMed]
  17. Li, Y.; Zhou, L. Solution algorithm of dwell time in slope-based figuring model. In Proceedings of the AOPC 2017: Optoelectronics and Micro/Nano-Optics, SPIE, Beijing, China, 4–6 June 2017; Volume 10460. [Google Scholar]
  18. Wang, Q.H.; Liang, Y.J.; Xu, C.Y.; Li, J.R.; Zhou, X.F. Generation of material removal map for freeform surface polishing with tilted polishing disk. Int. J. Adv. Manuf. Technol. 2019, 102, 4213–4226. [Google Scholar] [CrossRef]
  19. Wan, S.L.; Zhang, X.C.; Wang, W.; Xu, M. Effect of pad wear on tool influence function in robotic polishing of large optics. Int. J. Adv. Manuf. Technol. 2019, 102, 2521–2530. [Google Scholar] [CrossRef]
  20. Cao, Z.C.; Cheung, C.F.; Ho, L.T.; Liu, M.Y. Theoretical and experimental investigation of surface generation in swing process bonnet polishing of complex three-dimensional structured surfaces. Precis. Eng. 2017, 50, 361–371. [Google Scholar] [CrossRef]
  21. Li, L.X.; Zheng, W.J.; Deng, L.G.; Wang, X.; Wang, X.K.; Zhang, B.Z.; Bai, Y.; Hu, H.X.; Zhang, X.J. Optimized dwell time algorithm in magnetorheological finishing. Int. J. Adv. Manuf. Technol. 2015, 81, 833–841. [Google Scholar] [CrossRef]
  22. Han, Y.J.; Zhu, W.L.; Zhang, L.; Beaucamp, A. Region adaptive scheduling for time-dependent processes with optimal use of machine dynamics. Int. J. Mach. Tools Manuf. 2020, 156, 103589. [Google Scholar] [CrossRef]
  23. Han, Y.J.; Duan, F.; Zhu, W.L.; Zhang, L.; Beaucamp, A. Analytical and stochastic modeling of surface topography in time-dependent sub-aperture processing. Int. J. Mech. Sci. 2020, 175, 105575. [Google Scholar] [CrossRef]
  24. Zheng, W.M.; Cao, T.N.; Zhang, X.Z. Applications of a novel general removal function model in the CCOS. Proc. Soc. Photo Opt. Instrum. Eng. 2000, 4231, 51–58. [Google Scholar]
  25. Ghigo, M.; Canestrari, R.; Spigal, D.; Novi, A. Correction of high spatial frequency errors on optical surfaces by means of ion beam figuring. Proc. Soc. Photo Opt. Instrum. Eng. Conf. Ser. 2007, 6671, 667114. [Google Scholar]
Figure 1. The basic distribution of removal function in some deterministic optical surfacing technologies: (a) Two-dimensional distribution of ideal Gaussian removal function, (b) One-dimensional distribution of CCOS removal function, (c) Kernel function distribution of ion beam figuring (IBF).
Figure 1. The basic distribution of removal function in some deterministic optical surfacing technologies: (a) Two-dimensional distribution of ideal Gaussian removal function, (b) One-dimensional distribution of CCOS removal function, (c) Kernel function distribution of ion beam figuring (IBF).
Micromachines 12 00471 g001
Figure 2. Approximate distribution of the removal function: (a) One-dimensional approximation, (b) Two-dimensional approximation.
Figure 2. Approximate distribution of the removal function: (a) One-dimensional approximation, (b) Two-dimensional approximation.
Micromachines 12 00471 g002
Figure 3. Elementary approximation solution for one dimension.
Figure 3. Elementary approximation solution for one dimension.
Micromachines 12 00471 g003
Figure 4. Mesh refinement principle for one dimension.
Figure 4. Mesh refinement principle for one dimension.
Micromachines 12 00471 g004
Figure 5. Mesh generation.
Figure 5. Mesh generation.
Micromachines 12 00471 g005
Figure 6. Superposition distribution map in grids.
Figure 6. Superposition distribution map in grids.
Micromachines 12 00471 g006
Figure 7. Schematic diagram of value taken on longitude and latitude line.
Figure 7. Schematic diagram of value taken on longitude and latitude line.
Micromachines 12 00471 g007
Figure 8. Four kernel function scopes.
Figure 8. Four kernel function scopes.
Micromachines 12 00471 g008
Figure 9. Three removal function scopes.
Figure 9. Three removal function scopes.
Micromachines 12 00471 g009
Figure 10. Two removal function scopes.
Figure 10. Two removal function scopes.
Micromachines 12 00471 g010
Figure 11. Simulation results of removal function superposition.
Figure 11. Simulation results of removal function superposition.
Micromachines 12 00471 g011
Figure 12. Approximation results for the distribution of plane removal amount (a) Grid spacing L = R, (b) Grid spacing L = R/2, (c) Grid spacing L = R/4, (d) Ideal plan distribution.
Figure 12. Approximation results for the distribution of plane removal amount (a) Grid spacing L = R, (b) Grid spacing L = R/2, (c) Grid spacing L = R/4, (d) Ideal plan distribution.
Micromachines 12 00471 g012
Figure 13. Approximation results for the distribution of spherical removal (a) Grid spacing L = R, (b) Grid spacing L = R/2, (c) Grid spacing L = R/4, (d) Ideal hemispheric distribution.
Figure 13. Approximation results for the distribution of spherical removal (a) Grid spacing L = R, (b) Grid spacing L = R/2, (c) Grid spacing L = R/4, (d) Ideal hemispheric distribution.
Micromachines 12 00471 g013
Figure 14. Approximation results for the distribution of arbitrary removal (a) Grid spacing L = R, (b) Grid spacing L = R/2, (c) Grid spacing L = R/4, (d) Original surface distribution.
Figure 14. Approximation results for the distribution of arbitrary removal (a) Grid spacing L = R, (b) Grid spacing L = R/2, (c) Grid spacing L = R/4, (d) Original surface distribution.
Micromachines 12 00471 g014
Figure 15. The self-developed IBF machine.
Figure 15. The self-developed IBF machine.
Micromachines 12 00471 g015
Figure 16. Removal function in case 1: (a) surface map, (b) surface profile.
Figure 16. Removal function in case 1: (a) surface map, (b) surface profile.
Micromachines 12 00471 g016
Figure 17. Full aperture data in case 1: (a) Initial surface, (b) Predicted surface, (c) Machined surface.
Figure 17. Full aperture data in case 1: (a) Initial surface, (b) Predicted surface, (c) Machined surface.
Micromachines 12 00471 g017
Figure 18. Case 1 90% aperture data: (a) Initial surface, (b) Predicted surface, (c) Machined surface.
Figure 18. Case 1 90% aperture data: (a) Initial surface, (b) Predicted surface, (c) Machined surface.
Micromachines 12 00471 g018
Figure 19. Removal function in case 2.
Figure 19. Removal function in case 2.
Micromachines 12 00471 g019
Figure 20. Initial surface in case 2: (a) φ150 mm aperture, (b) φ100 mm aperture.
Figure 20. Initial surface in case 2: (a) φ150 mm aperture, (b) φ100 mm aperture.
Micromachines 12 00471 g020
Figure 21. Predicted surface in case 2 (a) φ150 mm aperture (b) φ100 mm aperture.
Figure 21. Predicted surface in case 2 (a) φ150 mm aperture (b) φ100 mm aperture.
Micromachines 12 00471 g021
Figure 22. Polished surface in case 2: (a) φ150 mm aperture, (b) φ100 mm aperture.
Figure 22. Polished surface in case 2: (a) φ150 mm aperture, (b) φ100 mm aperture.
Micromachines 12 00471 g022
Table 1. Parameters of the polishing process.
Table 1. Parameters of the polishing process.
ParameterValue
Ion beam voltage800 V
Ion beam current60 mA
Ion beam Angle
Processing distance150 mm
Grating spacing2 mm
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Zhang, Y.; Kang, R.; Ji, F. An Elementary Approximation of Dwell Time Algorithm for Ultra-Precision Computer-Controlled Optical Surfacing. Micromachines 2021, 12, 471. https://doi.org/10.3390/mi12050471

AMA Style

Wang Y, Zhang Y, Kang R, Ji F. An Elementary Approximation of Dwell Time Algorithm for Ultra-Precision Computer-Controlled Optical Surfacing. Micromachines. 2021; 12(5):471. https://doi.org/10.3390/mi12050471

Chicago/Turabian Style

Wang, Yajun, Yunfei Zhang, Renke Kang, and Fang Ji. 2021. "An Elementary Approximation of Dwell Time Algorithm for Ultra-Precision Computer-Controlled Optical Surfacing" Micromachines 12, no. 5: 471. https://doi.org/10.3390/mi12050471

APA Style

Wang, Y., Zhang, Y., Kang, R., & Ji, F. (2021). An Elementary Approximation of Dwell Time Algorithm for Ultra-Precision Computer-Controlled Optical Surfacing. Micromachines, 12(5), 471. https://doi.org/10.3390/mi12050471

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop