Next Article in Journal
A Novel Approach for Target Attraction and Obstacle Avoidance of a Mobile Robot in Unknown Environments Using a Customized Spiking Neural Network
Previous Article in Journal
Enhancing Privacy in Large Language Model with Homomorphic Encryption and Sparse Attention
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Approach for High-Precision Evaluation of Sphericity Errors Using Computational Geometric Method and Differential Evolution Algorithm

1
Center of Ultra-precision Optoelectronic Instrument engineering, Harbin Institute of Technology, Harbin 150080, China
2
Key Lab of Ultra-precision Intelligent Instrumentation (Harbin Institute of Technology), Ministry of Industry and Information Technology, Harbin 150080, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(24), 13144; https://doi.org/10.3390/app132413144
Submission received: 31 October 2023 / Revised: 2 December 2023 / Accepted: 6 December 2023 / Published: 11 December 2023

Abstract

:
The sphericity error is a critical form and position tolerance for spheres. We explored the distribution of sphericity errors within the solution space to achieve a high-precision evaluation using the minimum zone criteria. Within local solution spaces, we propose treating the evaluation of sphericity errors as a unimodal function optimization task. And computational geometric methods are employed to achieve highly accurate solutions within the local solution spaces. Subsequently, we integrated the computational geometric method with the differential evolution algorithm (DE algorithm). By centering on individual population members of the DE algorithm, we partitioned the local solution spaces and utilized the best solutions within them to optimize the population. With the gradual convergence of the DE algorithm, we successfully achieved the high-precision resolution of sphericity errors. The experimental results demonstrate a significant order-of-magnitude improvement in precision compared to traditional algorithms in the field of sphericity error evaluation, with uncertainty levels reaching magnitudes of 10−14 mm. Moreover, this method enhances the accuracy of sphericity error evaluation by approximately 10% for three-coordinate measuring machines. Additionally, we conducted ablation experiments to validate the effectiveness of the proposed computational geometric method. In summary, this approach enables the high-precision evaluation of sphericity errors and provides a practical methodology for applying ultra-precision spheres in precision engineering.

1. Introduction

The high-precision sphere constitutes a pivotal component within the realm of precision instruments. A sphere’s primary form tolerance characteristic manifests as a sphericity error. This inherent sphericity error, in turn, introduces unwelcome vibrations and noise, directly exerting a tangible influence on these precision instruments’ precision and service longevity [1,2]. The measurement and evaluation of the sphericity error are primarily conducted using three-coordinate measuring machines (CMMs) [3,4]. This advanced technology enables the acquisition of extensive contour data for the spheres under examination, which are then meticulously stored in the XYZ coordinate format, forming a substantial dataset [5]. It is crucial to underscore that this dataset’s precise and efficient processing stands as a critical imperative for the accurate evaluation of the sphericity error in spherical objects.
According to the standards set by the American Society of Mechanical Engineers (ASME) and the International Organization for Standardization (ISO) [6,7], the processing of datasets can be executed by applying four distinct criteria: the least squares criterion, the minimum circumscribed criterion, the maximum inscribed criterion, and the minimum zone criterion. The minimum zone criterion is the general default criterion for evaluating the sphericity error [8,9,10]. The aforementioned criterion enables the acquisition of a unique minimum sphericity error, and it is also the arbitration method when there is inconsistency in the evaluation of the sphericity error. The minimum zone criterion is an unconstrained and non-differentiable optimization problem, which poses significant challenges in terms of finding a viable solution [11,12,13].
Although the ASME and ISO standards stipulate that the evaluation of sampling point determination should be based on the minimum zone criterion, these standards do not provide a specific computational method for evaluating the spherical error based on the minimum zone criterion. Therefore, the development of a highly accurate and robust method for evaluating sphericity errors is of utmost urgency. Many scholars have extensively examined this issue by presenting a range of computational geometric methods [8,14,15,16,17,18], region search techniques [19], heuristic algorithms [20,21,22], multi-objecitve optimization methods [23,24,25], and bio-inspired intelligent optimization algorithms [26,27,28,29,30] to address the problem.
Among these approaches, region search methods are subject to inherent limitations that hinder their ability to accurately determine solutions for sphericity errors. Computational geometric methods for sampling point filtration exhibit enhanced search efficiency and improved results when dealing with fewer sampling points. Nevertheless, as the number of sampling points increases, the search cost escalates rapidly, resulting in unpredictable and fluctuating results. Heuristic algorithms are prone to being influenced by local optima and often struggle to converge effectively. Bio-inspired intelligent optimization algorithms yield more accurate results. However, they involve complex computations and numerous parameters that require optimization. Moreover, a sphericity error assessment, which is carried out to determine the optimal sphere center, minimizing the radius difference of concentric spheres that encompass the sampled contour, is fundamentally a single-objective optimization problem. Despite the potential advantages of applying multi-objective optimization methods to this evaluation, there is a paucity of scholars incorporating them into sphericity error assessment.
To address the problem of inadequate precision in sphericity error evaluation methods, we conducted an in-depth examination of the spatial distribution characteristics of sphericity errors. Based on our research, it is recommended to approach the issue of sphericity errors in the local solution space by treating it as a problem of optimizing a unimodal function. By utilizing computational geometric techniques, it is possible to achieve highly accurate solutions within the confines of the local solution space. Subsequently, by leveraging the advantages of heuristic algorithms for global space exploration, we successfully achieved a high-precision evaluation of sphericity errors. This method presents a novel research avenue for sphericity error evaluation methodologies.

2. Methods

2.1. Mathematical Model of the Minimum Zone Sphericity Error

According to the definitions of the minimum zone spherical error in ASME and ISO, for a given set of sampled contour points, if all data points fall within two concentric spheres, the difference in the radii of these two spheres is referred to as the minimum zone spherical error. Therefore, when evaluating spherical error based on the minimum zone criterion, the central task is to find an ideal sphere center p 0 that minimizes the difference in radii between the two concentric spheres (where the radius difference represents the sphericity error, denoted as e) [27,28,31,32,33]. This task is expressed by Equation (1).
e = min P 0 Q max i 1 , n p i p 0 2 min i 1 , n p i p 0 2
In the equation, e represents the spherical error, where p i d , i = 1, 2, …; n denotes the measurement points in a d-dimensional space; and p 0 is the center of the concentric sphere. Q is a cubic region containing the globally optimal solution, which we call the global solution space.
Figure 1 illustrates the spatial geometric structure of spherical error evaluation based on the minimum zone criterion and the various parameters involved. The construction of the global solution space Q requires the computation of the least squares sphere center and the least squares spherical error based on the least square criterion. Subsequently, a cube with a side length equal to the least squares spherical error is constructed around the least squares sphere center, allowing us to infer that the ideal sphere center must lie within this search region.

2.2. Sphericity Error Spatial Distribution Characteristics

From Equation (1) and Figure 1, it can be observed that the problem of evaluating the sphericity error based on the minimum zone criterion is non-differentiable and unconstrained. Finding an ideal sphere center p 0 within three-dimensional space, such that the difference in the radii is minimized when the concentric spheres encompass the spherical contour, poses a highly challenging problem. Meanwhile, existing studies suggest that the ideal sphere center p 0 is not unique [10,21]. There are several local optimal solutions near the optimal solution of the sphericity error, and the sphericity error decreases rapidly with a large gradient at the location far from the optimal solution [21].
Therefore, we have undertaken a study on the spatial distribution characteristics of spherical error. We conducted an analysis of the data sourced from reference [15]. Subsequently, we performed 27,000 samplings from the global solution space, calculating the spherical errors for each sampled point. Figure 2a illustrates the distribution of these sampled points, with the color gradient representing the magnitudes of the associated spherical error values. Following this, within this array of 27,000 points, we scrutinized the variations in spherical error values along the directions parallel to the XYZ axes. This scrutiny was carried out with the center OLS (X, Y, Z) of the least squares sphere as the central point, independently assessing the variations along each axis.
Meanwhile, Figure 2b–d show the trend of the sphericity error near the center OLS (X, Y, Z) of the least squares sphere, where the three lines passing through OLS and running parallel to the X-axis, Y-axis, and Z-axis, respectively, are constructed. It is clearly discernible from Figure 2b–d that, at positions considerably distant from the ideal solution, there is a significant gradient in the variation of the spherical error. Conversely, as the position approaches the ideal solution, the evolution of the spherical error gradually becomes more intricate.
Subsequently, we arbitrarily elected several localized solution spaces from the global solution space, analyzing the spherical error variation of the sampling points within these specific solution spaces. Figure 2e–g illustrate the distribution of sampling points within these localized solution spaces. It can be observed from Figure 2e–g that within smaller regions, the variation of the spherical error manifests a relatively uncomplicated nature, approximating the characteristics of a unimodal function. If the global solution space is divided into different size subregions, then each subregion can be approximated as a single-peaked function for processing. In local spatial domains, unimodal functions can be effectively solved using computational geometry methods, yielding high-precision results.
The algorithm presented in this paper primarily leverages the characteristics of spherical error gradient variations and the approximate unimodal nature of spherical error changes within local spaces for its design. This is the foundation of the local enhancement principle in this paper.

2.3. Principle of the Proposed Method

As demonstrated in Section 2, a comprehensive investigation was carried out to analyze the spatial distribution characteristics of the sphericity error. Given the non-uniqueness of the ideal solution p 0 and the potential existence of local optima, the high-precision sphericity error evaluation directly from the global solution space presents a considerable challenge. However, through the adept partitioning of the global solution space into multiple local solutions, it becomes possible to approach the evaluation of the sphericity error as a unimodal function optimization problem within each local space. By employing computational geometry techniques, it is possible to attain high-precision solutions within the confines of these local solution spaces.
Subsequently, the integration of the computational geometry method employed with the DE algorithm is performed. We establish the positions of local solution spaces by utilizing the individuals that are present in the population of differential evolution. During each iteration, the sizes of the local solution spaces are adaptively adjusted based on the distances between the individuals in the population and the global optimum. Additionally, incorporating Lévy flights and simulated annealing techniques allows for the perturbation of the positions of individuals within the population, thereby augmenting the diversity of local solution spaces. After attaining high-precision solutions within the designated local spaces, the DE algorithm individuals are relocated to their respective optimal positions within these local solution spaces. With each successive iteration of the DE algorithm, the population progressively approaches convergence towards the global optimum. As a result, an optimal solution for minimizing the sphericity error is eventually obtained. The overall algorithmic process is depicted in Figure 3. We denote the computational geometric method integrated with the DE algorithm as “Adaptive Local Enhancement”. The detailed steps of this approach are outlined in Section 3.2.

2.4. Adaptive Local Enhancement Method

2.4.1. Adaptive Partitioning of Local Solution Spaces

Considering the gradually complex variations in the sphericity error when the candidate sphere center approaches the ideal sphere center p 0 , we devised an adaptive method for partitioning the local solution space. By leveraging the individuals of differential evolution, this method divides the neighborhood of the population into local solution spaces of varying sizes.
During each iteration of the DE algorithm, we calculate the distance denoted as d i   , i = 1,2 , , n between an individual’s optimal solution and the population. The determination of the distance is contingent upon the expression d i = d i / m a x ( d i ) . Subsequently, we determine the size of the individual’s neighborhood space according to L i = K d i + e p s . Here, K is a scaling factor, and eps is a constant that prevents L from becoming zero. The neighborhood search range for each individual, represented as x i , y i , z i , falls within the interval of x i L i , x i + L i , y i L i , y i + L i and z i L i , z i + L i . Figure 4 shows the sizes and distributions of these local solution spaces.

2.4.2. Computational Geometric Method for Local Solution Spaces

Once the local solution spaces are determined based on the positions of population individuals, a precise solution need to be conducted for each local solution space. Since we regard the problem of solving the spherical deviation error under the constraint of local solution space as a unimodal function problem, we employ the computational geometry approach shown below.
In the three-dimensional space, the local enhancement process of individuals is shown in Figure 5. According to the predetermined search interval, the initial Z-directional search interval is determined as Z L i = z i L i , Z R i = z i + L i , and then the search is performed on two planes    Z 1 i , Z 2 i , where Z 1 i = Z R i + 2 Z L i / 3 , Z 2 i = 2 Z R i + Z L i / 3 . The initial search interval in the Y direction is Y L i = y i L i , Y R i = y i + L i on two planes   Z 1 i , Z 2 i , and a search for the local optimal solution in the one-dimensional direction, Y 1 i = Y R i + 2 Y L i / 3 and Y 1 i = Y R i + 2 Y L i / 3 , is performed on both planes. Shrinking the search interval according to strategy 1, four positions, Y 1 Z 1 X best , Y 2 Z 1 X best , Y 2 Z 1 X best , and Y 2 Z 2 X best , of the local optimal solution can be obtained when satisfying the following condition: XRXL < 1 × 10n. The local optimal solutions on the Z 1 i , Z 2 i planes are compared, respectively, and the search interval in the Y-direction on the two planes, Z 1 i , Z 2 i , is shrunken according to strategy 2. After satisfying the condition YRYL < 1 × 10n, the optimal solutions f X , Y , Z 1 i and f X , Y , Z 2 i on the Z 1 i , Z 2 i planes are obtained. Afterward, the Z-directional search space is shrunken according to strategy 3. While satisfying Z R Z L < 1 × 10 n , the local optimal solution is obtained.
Strategy   1 :   X L i = X 1 i 1 , i f : f X 1 i 1 , Y , Z > f X 2 i 1 , Y , Z X R i = X 2 i 1 , e l s e
Strategy   2 :   Y L i = Y 1 i 1 , i f : f Y 1 Z X best , Y 1 i 1 , Z > f X best Y 2 Z , Y 2 i 1 , Z Y R i = Y 2 i 1 , e l s e
Strategy   3 :   Z L i = Z 1 i 1 , i f : f Z b e s t , Y b e s t , Z 1 i 1 > f X best , Y b e s t , Z 2 i 1 Z R i = Z 2 i 1 , e l s e
Furthermore, for the one-dimensional search area corresponding to Y 1 i on the Z 1 i plane in Figure 5, the detailed search process is illustrated in Figure 6. The middle points, X 1 i 1 = ( X R i 1 + 2 X L i 1 ) / 3 and X 2 i 1 = ( 2 X R i 1 + X L i 1 ) / 3   , are determined based on the neighborhood searching range of the individuals, X L i 1 = x i 1 L i 1 , X R i 1 = x i 1 + L i 1 , after which the magnitudes of f ( X 1 i 1 ) and f ( X 2 i 1 ) are compared. If f X 1 i 1 > f X 2 i 1 , then set X L i = X 1 i 1 ; otherwise, set X R i = X 2 i 1 . Through the strategy above, used to shrink the search interval, the optimal solution is considered as the sphericity error at the X R i position when X R i X L i < 1 × 10 n is satisfied.

2.5. Algorithm Flow

2.5.1. Initialization

Firstly, the population and the global temperature of the simulated annealing are initialized, and the initial least square spherical center x 0 , y 0 , z 0 and least square sphericity error S 0 are solved based on the least square solution. The following is the approach for solving the least squares sphericity error S0 and the least squares sphere center OLS based on the least squares criterion. According to the coordinate dataset of the measurement points, p i = [ x i , y i , z i   ] T , i = 1 , 2 , , N , establish the equation set A x = b , where A = [ p 1 , p 2 , , p N , I ] T R m × n , b = [ p 1 2 2 , p 2 2 2 , , p N 2 2 ] T , and I = [ 1 , 1 , , 1 ] T . By solving the equation set based on the least square solution, the coordinates of the sphere center OLS  x 0 , y 0 , z 0 and the least square sphericity S0 are obtained. Then, the initial population of the algorithm can be generated in the ranges of x 0 S 0 , x 0 + S 0 , y 0 S 0 , y 0 + S 0 , and z 0 S 0 , z 0 + S 0 , denoted as X.

2.5.2. Mutation

The mutation strategy is an important part of the DE algorithm. A suitable mutation strategy can promote the convergence of the DE algorithm to the global optimal solution. This algorithm adopts the DE/current-to-best/1/bin strategy, which uses the elite solution to guide the mutation of the existing solutions to approach the optimal one. This mutation strategy is shown in Equation (2).
V i t + 1 = X i t + F X b e s t t X i t + F X k 1 t X k 2 t
where V i t + 1 is the i-th element of the vectors of mutation in the (t + 1)th generation; F is the scaling factor and a random number in (0, 0.5); X b e s t t is the optimal solution in the t-th generation; X i t is the i-th element among the target vectors of the t-th generation; and X k 1 t , X k 2 t is the random point ( k 1 k 2 ).

2.5.3. Lévy Flight Based on Simulated Annealing Strategy

A Lévy flight perturbed population based on simulated annealing is introduced to enhance the diversity of the local solution space distribution. This strategy enables individuals to explore the global solution space with occasionally long trajectories and short random movements. During the exploration process of the population individuals, their corresponding local solution spaces also move within the solution space. The formula for the Lévy flight in this paper is given by s = u / | v | 1 / β , where   u ~ N 0 , σ 2 ,   v ~ N ( 0 ,   1 ) , and σ = Γ 1 + β sin π β / 2 β Γ 1 + β / 2 2 β / 2 1 / β .
The steps of the Lévy flight based on the simulated annealing strategy are as follows:
Step 1: Set an initial temperature T to a sufficiently high value.
Step 2: Consider the original solution of the individuals as S1 and the new solution after Lévy flight perturbation as S2. Calculate the difference df = f (S2) − f (S1) for each individual in the population, where f represents the object function.
Step 3: If df < 0, accept the new solution S2.
Step 4: If df ≥ 0, generate a random number R uniformly distributed in the interval (0, 1) and calculate the acceptance probability exp(−df/T) for S2. If exp(−df/T) > R, accept S2 as the new solution; otherwise, keep the solution unchanged as S1.
Step 5: Reduce the global temperature by setting T = q × T, where q < 1.
As the DE algorithm progresses through iterations, the global temperature gradually decreases, resulting in a reduced probability of accepting inferior solutions. When the temperature reaches a certain threshold, the search process terminates.

2.5.4. Adaptive Local Enhancement

For a comprehensive understanding of the adaptive local enhancement process, please consult Section 2.4. Figure 5 illustrates the process of local enhancement.
Once the local solution space’s search range is determined, a local optimal solution would be searched within the local solution spaces using the computational geometric method mentioned earlier. Following a greedy strategy, individuals move to the position of the local solution space’s optimal solution if it is superior to their current position; otherwise, they remain at their initial positions. Then, an enhanced population is generated due to the local enhancement process, denoted as V.

2.5.5. Crossover and Selection

By mixing the new population V and the initial population X with a certain probability, a new population U is generated. Compare each individual in the populations U and X, and retain the ones with smaller sphericity errors for the next iteration.
The loop is ended until the conditions of convergence are satisfied, and the global optimal solution by the minimum zone criterion for evaluating the sphericity error is obtained.

3. Results

3.1. Verification with Datasets in the Literature

Due to the limited availability of publicly accessible datasets in the field of sphericity error research, current research papers often refer to the same set of publicly available spherical datasets for comparative analyses. Notably, the four datasets [15,16,31,33] used in the experiments of this section have been widely cited in numerous high-quality publications in the field of sphericity error research. As a result, these datasets hold considerable authority. By conducting experiments using these four datasets, we can effectively compare our approach with the leading algorithms in the field of sphericity error evaluation under identical experimental conditions. This reasoning highlights one of the main reasons for selecting these four datasets for experimentation in our study.
The experiments were conducted under identical conditions. The number of populations (N) was set to 30, the local enhancement parameter (n) was set to 12, the initial temperature (T) was set to 100, the annealing rate (q) was set to 0.9, and the temperature threshold was set to 0.001. Using the aforementioned settings, each dataset was independently calculated ten times. The mean and standard deviation of the results were calculated and compared with the results obtained using state-of-the-art algorithms.
Dataset 1, obtained from reference [34], comprises 100 randomly generated coordinate points within two concentric spheres with radii of 10.0 and 9.0. The data distribution is illustrated in Figure 7, where the sphericity error is determined as 1.0 based on five data points from the inner and outer spheres. The sphericity error of dataset 1 has been validated in several studies [21,34], and the algorithm’s accuracy in this paper was initially assessed using dataset 1.
Ten sets of replication experiments were performed for dataset 1. The mean experimental result is 1.00000000000003 mm, with a standard deviation of 2.0 × 10−14 mm. Figure 8 displays one convergence curve from the ten experiments. Table 1 presents a comparison with the results reported in the literature.
Dataset 2, consisting of 384 coordinate points, is provided in reference [16]. The 384 coordinate points were obtained through real spherical sampling using the birdcage method and were evenly distributed along 12 lines on the sphere. Figure 9 illustrates the distribution of dataset 2. The reference provides the optimal solution as 0.015385 mm. Ten sets of replication experiments were conducted for dataset 2. Table 2 presents the results of the ten experiments, with a mean experimental result of 0.01538487054930 mm and a standard deviation of 3.0 × 10−14 mm.
Figure 10 displays one convergence curve from the ten experiments. The least square solution obtained by our algorithm is 0.01639862 mm, and the final convergence solution is 0.01538487054931 mm. Table 3 presents a comparison with the results reported in the literature.
Dataset 3, consisting of 25 coordinate points, is provided in reference [36]. Figure 11 illustrates the distribution of dataset 3, with the coordinate points uniformly distributed on the sphere. Ten sets of replication experiments were conducted for dataset 3. Table 4 presents the results of the ten experiments, with a mean experimental result of 3.33251761616355 mm and a standard deviation of 0.8 × 10−14 mm.
Figure 12 displays a convergence curve from one of the ten experiments. The least square solution obtained by our algorithm is 3.735706 mm, and the final converged solution is 3.33251761616355 mm. Table 5 presents a comparison with the results reported in the literature.
Dataset 4, consisting of 50 coordinate points, is provided in reference [15]. Figure 13 illustrates the distribution of dataset 4, with the coordinate points randomly distributed on the sphere. Ten sets of replication experiments were conducted for dataset 4. Table 6 presents the results of the ten experiments, with a mean experimental result of 0.00766019493788 mm and a standard deviation of 0.2 × 10−14 mm.
Figure 14 displays a convergence curve from one of the ten experiments. The least square solution obtained by our algorithm is 0.008486342 mm, and the final converged solution is 0.00766019493785 mm. Table 7 presents a comparison with the results reported in the literature.

3.2. Simulated Experiments

To further assess the algorithm’s accuracy, we employed the methodology described in reference [34] to generate simulated data, creating four distinct datasets. Initially, two concentric spherical surfaces, S 1 and S 2 , were established with radii r s 1 and r s 2 , where r s 1 = r s 2 + e , and e represents the simulated spherical error.
The points on these concentric spheres were identified as control points, with ten control points generated on each inner and outer sphere. Utilizing Formula (3), by randomly generating values for θ and ϕ , we generated 100 random sampling points within the concentric sphere, culminating in a dataset comprising 120 simulated sampling points. The 20 control points on both the inner and outer spheres effectively constrained the spherical error of the simulated dataset to the difference in the radii between the inner and outer spheres. Table 8 provides comprehensive information of the spherical error, center coordinate, radii, and the number of simulated sampling points for each of the four datasets. Additionally, Figure 15 visually portrays the data distribution of these four simulated datasets.
X = R cos θ sin ϕ ,       Y = R sin θ cos ϕ ,       Z = R cos ϕ
We conducted ten experiments for each dataset, showcasing the optimal estimates and standard deviations in Table 9. Furthermore, detailed simulated data are provided in Appendix B, Table A2, Table A3, Table A4 and Table A5. Table 9 shows that the algorithm’s solutions for spherical error exhibit a deviation within the magnitude of 10−14 mm from the ideal spherical error. Moreover, the standard deviation of the results from ten experiments also lies within the magnitude of 10−14 mm. The results attest to the algorithm’s successful identification of the ideal sphere center and spherical error.

3.3. Ablation Experiment

From the perspective of the algorithmic framework, this algorithm is based on the differential evolution algorithm, thus sharing evident similarities with traditional single-objective search algorithms. The algorithm is specifically designed to address the distribution characteristics of the sphericity error. It proposes that the resolution of the sphericity error can be treated as a unimodal function in the local solution space. Subsequently, the local enhancement method for solving this unimodal function is synergistically integrated with the differential evolution algorithm, demonstrating higher precision and robustness.
To delve into the pivotal aspects of the algorithm in this paper, we conducted ablation experiments, employing the same parameter configuration outlined in Section 3.1. We designated the algorithm with the removal of local enhancement steps as algorithm A and designated the algorithm retaining these steps as algorithm B. Subsequently, ten experiments were performed for each algorithm, with the standard deviation of the results computed. The experimental findings are documented in Table 10, from which it is evident that removing local enhancement methods significantly diminishes the algorithm’s robustness, reducing the standard deviation from a magnitude of 10−14 mm to a magnitude of 10−2~10−7 mm. This experiment’s results substantiate the favorable outcome achieved by amalgamating the proposed local enhancement method and differential evolution algorithm in this paper.

3.4. Experiments with Datasets Obtained from CMM

To assess the proposed method’s practicality, we measured a metal hemispherical shell resonator, as depicted in Figure 16b. The measurements were carried out in a controlled environment, as shown in Figure 16a, using a Hexagon CMM.
For the measurement of the specimen, we selected five cross-sections on the surface. As detailed in Table A1, one hundred and thirty measurement points were taken. The iterative process convergence, reconstruction model of the dataset, and evaluation results are presented in Figure 17 and Figure 18.
Ten replication experiments were conducted for the measurement data. The sphericity error obtained by the CMM is 0.008396638 mm. The results of these ten experiments are presented in Table 11, indicating a mean result of 0.00747188442066 mm, with a standard deviation of 2.3 × 10−14 mm. The spherical center obtained by our algorithm is (0.00218039158405 mm, −0.00226432610011 mm, −23.9943439405 mm). Our algorithm demonstrated an 11% improvement in accuracy compared to the CMM result.

4. Discussion

Evaluating the sphericity error based on the minimum zone criterion is an indispensable and challenging task in metrology. Thus, this study aims to achieve a high-precision evaluation of spherical contour sampling data. Consequently, we conducted an in-depth investigation into the spatial distribution characteristics of the sphericity error. We propose treating the local space of the sphericity error as a unimodal function and utilizing computational geometry methods to obtain high-precision solutions for the local space. Subsequently, combined with the advantages of heuristic algorithms in searching the global space, we successfully achieved exact solutions for sphericity errors.
To validate the efficacy of the proposed approach, we conducted experiments using datasets that are widely cited in the realm of spherical error research, which are featured in four high-level scholarly papers. The experimental results demonstrate a significant improvement in precision, surpassing the existing achievements in the field of spherical error research by several orders of magnitude. Meanwhile, the uncertainty of the algorithm’s solution reaches 10−14 orders of magnitude, which far exceeds that of the traditional heuristic algorithm. Furthermore, when processing the distribution of data with specific sphericity or dealing with a large number of data points, the global optimal solution can be searched out accurately using the proposed algorithm, which shows stronger robustness.
Subsequently, we conducted ablation experiments on the four datasets, as mentioned earlier. The experimental results show that, after removing the local enhancement method, the algorithm’s standard deviation decreases from the order of 10−14 mm to the order of 10−2~10−7 mm. The ablation experiments demonstrate the effectiveness of the proposed local enhancement method in this paper.
Finally, we apply the algorithm to evaluate the sphericity errors of the hemispherical shell resonator and obtain its sphericity error of 0.00747188442066 mm with a standard deviation of 2.3 × 10−14 mm for ten experiments. Compared to the results obtained from the CMM, this algorithm can improve the accuracy of sphericity error evaluation by approximately 10%.
In this study, each member of the population in the DE algorithm is considered the center point for partitioning local solution spaces. During each iteration of the algorithm, computational geometric techniques are employed to solve for these individual local solution spaces. This approach has yielded outstanding precision in evaluating sphericity errors. However, it has also resulted in an algorithmic time complexity of O(n4), leading to elevated time costs for each evaluation. Furthermore, we conducted assessments of the average time required for each evaluation. Under the specified parameter settings, completing a single evaluation takes approximately 300 s. Compared to the time costs of 0.1 to 1 s in existing studies, the average time cost of our algorithm has increased several-fold. It must be acknowledged that this algorithm is not suitable for components with spherical geometry that require large-scale production and have relatively low precision requirements, such as ball bearings. In these contexts, the speed of the sphericity error assessment takes precedence over precision. However, for crucial spherical components in the entire realm of research that necessitate only one or two ultra-high precision devices, such as the ball bearings in astronomical telescopes and the electrostatic rotor in Gravity Probe B, the precision and reliability of the sphericity error assessment outweighs speed. Nevertheless, pursuing high speed and precision remains a common goal across all research domains. Therefore, reducing time costs while maintaining the existing precision level is a critical direction for future research improvements.
In summary, this study proposes an algorithm that is aimed at the high-precision evaluation of sphericity errors and that is capable of performing a highly accurate evaluation of sphericity errors. However, to achieve this high precision in sphericity error evaluation, the algorithm has made certain sacrifices in terms of the time cost. This method provides practical tools and techniques for applying ultra-precision spheres in precision engineering. It introduces new perspectives for sphericity error evaluation based on the minimum zone criterion.

Author Contributions

Conceptualization, D.Z. and J.C.; methodology, D.Z.; software, D.Z.; validation, J.C.; formal analysis, D.Z.; investigation, D.Z. and J.C.; resources, J.C., Z.W. and Y.S.; data curation, D.Z.; writing—original draft preparation, D.Z.; writing—review and editing, D.Z.; visualization, D.Z., Z.W. and Y.S.; supervision, J.C.; project administration, J.C.; funding acquisition, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China of No. 52075133.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data that support the findings of this study are included within this article.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

Table A1. CMM raw coordinate indication point dataset from a hemispherical shell resonator measurement (units: mm).
Table A1. CMM raw coordinate indication point dataset from a hemispherical shell resonator measurement (units: mm).
No.XYZNo.XYZ
18.296414−3.02143−11.32676611.612795.831947−15.6435
26.763515−5.67803−11.32676712.645082.993801−15.6456
34.413842−7.64931−11.32666812.994360.00061−15.6446
41.53153−8.6978−11.32646914.364860.00079−18.3304
5−1.53198−8.69712−11.32627014.050592.984275−18.3304
6−4.41293−7.6497−11.32717113.123125.840188−18.3305
7−6.76176−5.67733−11.32647211.621048.44174−18.3301
8−8.2954−3.02146−11.3276739.61038510.67402−18.33
9−8.827620.00175−11.3272747.18064612.43983−18.33
10−8.295493.019205−11.3266754.43613413.66083−18.3294
11−6.762265.675262−11.3273761.49960914.28478−18.3294
12−4.412117.645938−11.327777−1.5003414.28458−18.3299
13−1.531428.694617−11.327878−4.4363713.66036−18.3297
141.532118.694537−11.327179−7.1794312.43851−18.33
154.4130927.645988−11.326980−9.6091910.673−18.3294
166.7631655.674052−11.326781−11.61988.44048−18.3295
178.2961843.017535−11.327482−13.12195.839438−18.3305
188.8287840.00144−11.327583−14.04922.983165−18.3305
1911.124240.0013−13.281384−14.36220.00134−18.3307
2010.712112.999729−13.281685−14.0472−2.98888−18.3309
219.5055385.779181−13.280386−13.1193−5.84373−18.331
227.5924888.12935−13.279987−11.6175−8.44454−18.3309
235.1155969.876324−13.279988−9.60736−10.6763−18.3306
242.26061110.89058−13.281489−7.17731−12.4413−18.3302
25−0.7588711.09661−13.281490−4.43455−13.6626−18.3298
26−3.7231810.48024−13.282191−1.49886−14.287−18.33
27−6.412679.08655−13.2818921.500999−14.2875−18.3296
28−8.627047.018315−13.2813934.437544−13.6644−18.3289
29−10.20194.429719−13.2817947.181086−12.4437−18.3295
30−11.01871.511798−13.2827959.611795−10.6784−18.3291
31−11.0178−1.51657−13.28319611.62196−8.44553−18.3291
32−10.2016−4.43394−13.28359713.12441−5.84411−18.33
33−8.62646−7.02263−13.28159814.05164−2.98806−18.33
34−6.41128−9.08998−13.28149914.90208−2.96403−21.2295
35−3.72154−10.4842−13.281210014.03837−5.81562−21.2292
36−0.75777−11.0995−13.280810112.63434−8.44279−21.2289
372.261661−10.8914−13.278210210.74393−10.7462−21.2288
385.116226−9.87919−13.27981038.43965−12.6363−21.2285
397.592868−8.13236−13.27971045.8119−14.04−21.2286
409.505198−5.78223−13.28071052.961514−14.9037−21.2287
4110.71221−3.00317−13.2805106−0.00016−15.1941−21.2291
4212.64529−2.99836−15.6442107−2.96109−14.9016−21.229
4311.61302−5.83438−15.6435108−5.81005−14.0377−21.2287
449.954915−8.35568−15.6437109−8.43761−12.633−21.2283
457.759662−10.4272−15.6439110−10.7396−10.7434−21.229
465.14489−11.9356−15.6438111−12.6302−8.44041−21.2291
472.254347−12.8003−15.6441112−14.0335−5.81371−21.2291
48−0.75498−12.9743−15.6445113−14.8978−2.96381−21.2286
49−3.72379−12.4499−15.6451114−15.1907−0.00053−21.2288
50−6.49358−11.2552−15.6455115−14.92.961734−21.2291
51−8.91403−9.45329−15.6456116−14.0365.81165−21.2281
52−10.8542−7.14139−15.6459117−12.63188.43835−21.2283
53−12.2073−4.44475−15.6464118−10.74110.74128−21.2287
54−12.9032−1.50946−15.647119−8.4370212.63058−21.2273
55−12.90381.506108−15.6464120−5.8101514.0345−21.2297
56−12.20854.442633−15.6461121−2.960514.89967−21.229
57−10.85447.138399−15.64621220.0003515.19031−21.2285
58−8.914639.450499−15.64661232.96275414.8987−21.2293
59−6.4932111.25221−15.64571245.8122914.03412−21.2292
60−3.7233612.44723−15.64561258.4385612.63031−21.2276
61−0.7542612.97149−15.645912610.7423110.73945−21.2293
622.25487712.79544−15.644112712.631958.4362−21.2295
635.1452411.93096−15.644712814.035555.80935−21.2291
647.76003210.42313−15.646512914.899842.960004−21.229
659.9566958.350528−15.644913015.192340.00022−21.229

Appendix B

Table A2. Simulation dataset A’s coordinate data (units: mm).
Table A2. Simulation dataset A’s coordinate data (units: mm).
No.XYZNo.XYZ
12.707165893120495.507393568470157.8955727410491761−2.89280574054528−3.656410165450428.84661176098279
2−1.15214718780568−9.200813715486163.7440243058442562−1.51828118667085−9.13434218820972−3.77606589087291
3−1.01772769992569−1.844851049122499.77552546524520632.24430605422733−4.05526495095270−8.86103883157432
41.097457488867578.990301412068314.23912436513771641.006554784039672.749901943114569.56164599503323
5−1.28865793782677−2.147129614924069.6813849315151465−5.910083905502146.654267861560214.55980563803232
62.23712088095506−5.642802602635157.94695444739222661.056499865364443.12247361232877−9.44108128748343
70.963893786892309.38875252004527−3.30491072270252677.855010716126622.19242740865312−5.78724691162584
80.095160368306670.57085048124596−9.9832419151937268−7.39966754477964−4.849651414565394.66109444057282
93.858225913350034.61648956587795−7.98762771847785691.09631832016297−7.49379191536427−6.53003590114144
10−1.067211595385231.93885104125363−9.75203135750703703.53745810772145−1.78161148450475−9.18217909565151
116.26235696234287−7.591961875686091.77341482860745716.76213768840946−0.595111447414857.34298594770830
12−2.624997454589199.430011714067832.0455481993694872−7.89611121553721−3.79056887081821−4.82525144178618
130.291521555976138.309811991311325.5553691854505373−3.45995327527011−1.48232328353041−9.26452594664461
140.206305226061672.260365061277489.7390057254142974−0.44629221645521−7.21658669463085−6.90809564438316
152.99706563681856−2.10102303257860−9.30609523486725750.285201942022450.22255761270425−9.99346401882636
169.43834225665788−2.593744626098332.0470408294385676−7.976894768515462.822849261803815.32924080069324
17−1.606167491266341.85271355326272−9.69472740038624777.27215932101236−5.155586764188424.53163160673812
18−7.53700381107871−6.54949836084845−0.54573054674074787.56765130979587−1.134144082415976.43773397992703
199.990105250440840.161151671858430.41461656648064799.10595015430749−1.992758386136163.62086572061669
203.077943685955638.161928177949684.8897177929932180−1.65318551090508−8.646786929635124.74342379395656
21−8.881382197549103.814495336192212.5633428646052081−2.160289799959842.194621863242159.51403989099589
22−1.505284920914936.149811630761817.74041825385070823.712147142738651.380354751019599.18230164358843
23−2.633008373295042.66539640573202−9.2716212115014283−2.815831226182871.25258183256518−9.51326201418868
248.56581407890788−4.96012947251450−1.42269021052671840.25318226559820−3.54288213674587−9.34794068182515
251.82718891044369−7.94962037867871−5.78490312181898859.509594692478500.374316563755983.07042275507506
268.24289514570036−0.72547048296934−5.6150210875828486−0.30358937469407−0.37790575801959−9.98825311874762
27−7.700922920498965.216345343575033.67226461861028876.57243761907439−4.871885268688665.75047796586977
287.127078051268943.010621109054726.3356970025725188−6.014829235024795.099788241901126.14931759713210
29−0.589990657775370.821563780095479.9487207576771289−5.38098973200806−0.365722879882448.42088488112053
30−2.567160433330881.390561475089629.56431299758929909.64669012283484−2.633195622767270.08860297026600
318.22189230883354−2.82653115285224−4.94078328448101918.330831016592330.34658097412821−5.52061770906743
325.61203365939748−4.28243254741348−7.08278544660527921.16667377294868−5.985855510845737.92518310069517
33−8.30864658564425−2.911337117442284.74241852182297935.09725538078174−8.205832035287792.58502319927636
340.66411528936322−9.27973866238164−3.6668074905605994−7.94577182052644−0.15895145001070−6.06955829228184
355.595558088126514.52368413040483−6.9445014263983895−2.81033726124221−0.34145895649012−9.59090915571642
364.318218673580962.98193413994706−8.5124135271346796−8.84178799766724−1.30105332857164−4.48666162729917
37−7.27229290635095−0.94799136393752−6.7981707952144197−5.13240086305896−2.53115221949677−8.20072497861750
381.76614778493561−2.673175917623939.4728556908584898−3.26949290444136−7.24375242567628−6.06947007104547
398.904882203984933.82356465101649−2.4664859145076099−4.21044225736679−0.918016415191929.02383533011225
40−1.342163644617930.195807219696379.907595563406611000.972073480912521.35811285724233−9.85954941070868
41−2.144138332492056.54013975411457−7.254615216131381019.97467105975975−0.621340899713150.34633679834256
42−0.146415301362900.00329534577709−9.998934010711811023.07164278927927−7.79824967120529−5.45458638585138
439.33663475800817−2.262522302120202.776381432231371034.214240904615274.13157893765339−8.07281448962633
448.606153476163993.257784877801423.91419979493624104−8.209548626202172.492603050369185.13715697784912
450.92038470926800−8.232207130763795.60211794144169105−6.056467849467833.333396735207367.22549856527915
460.19039170935725−0.095598772820609.997732605039881069.290848095070471.345423569466773.44530959845342
474.630416991711348.84538152765776−0.564500385274651073.101061004908869.506145193798270.12969193913494
484.909406211640296.172215985803186.148295820985681088.25810929801433−2.842410877681324.87077201262455
49−1.356018129824572.67900851431207−9.538563880885401094.82089142825648−1.244325077204038.67241820026406
500.106976805159440.161044169835999.998136418364901108.08804152765836−4.99592166582711−3.10234592604136
510.54543741944347−2.743527033802539.600813061465511110.292299836573360.547274160586849.98074000592601
52−5.88145323430974−3.426093249214177.326015849529491122.353652950759061.64164975916963−9.57942345351573
531.983735888666362.69419177896352−9.423710658879681130.041984898792160.07803540541192−9.99961698811101
54−5.114691444164840.775700605079268.557939242570761146.311614333264790.52599244332027−7.73867005704883
55−1.616076899802905.507776474546118.188574449607061155.53492638583934−8.32813270602568−0.08363931052334
56−6.05720577992347−7.779578833681811.66986302157556116−1.94264324774318−1.185824954297039.73756059512558
571.04302816178437−3.270358363334189.39239093167803117−2.813715601087291.337730549191869.50229259300673
582.81296261497125−0.15442235185868−9.594973179666771186.371648021799890.02784080933777−7.70723405788024
590.75873469126528−0.27199953420838−9.96747399905116119−9.21326646969875−0.100052155232933.88660655643943
60−0.28319142384936−4.90514063534900−8.709741460284891206.25059255805694−7.800835365985160.27782060448905
Table A3. Simulation dataset B’s coordinate data (units: mm).
Table A3. Simulation dataset B’s coordinate data (units: mm).
No.XYZNo.XYZ
11.34866310044556−3.429644341977359.3726736443017361−6.03813011492763−7.76545546035172−0.57679931500770
26.28762586455096−5.810931588952135.2743801548729562−8.37817918397411−3.427349223696004.05958662316408
30.583982719366051.859559265503859.93207905167647630.522336956464313.252693968828149.58062372503846
4−6.79335005575682−6.713322094892382.5617921987511264−2.090266916102443.655294550109749.18639845896016
54.554466186821850.04056169255594−8.8528876949322365−9.81740328667709−0.736943093915161.07199737940335
63.79620842970985−7.535886161391825.39445802461468660.154455046924980.1262714899129110.09982417701270
7−9.836742007465241.066757296533960.67151129024692670.501013343107980.2275285412921410.09115232442710
8−6.61532798308528−3.65557088498325−6.28750106664630689.740768047983520.732312511585272.67987835617744
9−7.179850713523431.349742363526296.84106883141816690.10312758013623−9.700024955492562.08989240315612
107.85510652795512−5.947337236259451.9132982533428070−8.825752641630783.184054144062313.38931694468529
11−5.081748275717288.287468856296982.57282723147979716.95710495543139−1.95896587491257−6.88146426353406
12−4.17476966096297−1.389277223470429.01674947155353728.14147891565005−0.62754515133087−5.79960757253585
137.39523184514386−1.64735595814567−6.51258975851859734.73021563852859−8.34043159678178−2.60562618722498
141.757914251293149.35998423612448−3.29177150326352740.14885226061423−2.552333879475839.74172384433377
159.530374131005022.96783323871382−1.58629448613162750.241230849989833.35834622385500−9.35321815255340
161.06836312520230−9.740812530521691.5901954533232776−9.176012524170741.72155553651928−3.26545594466951
17−1.359614924449538.85325919098271−4.50977755465138771.24207417896070−9.833198022715710.26564340074693
18−2.54886701298226−1.16362897579127−9.45964508557108783.95019276558823−0.43913982247625−9.11333454632599
198.34553977725865−3.82986427025728−3.97030760282703791.01943786643079−8.650116784710334.85290752422392
20−9.067266783458752.76301281032927−2.8781987119394380−1.37244324949174−0.44407560653978−9.77602856563217
21−2.701206962042157.930418627385585.6531875634250681−3.01546757149110−8.628539949422543.85585826136877
220.40484779625019−2.400890292169869.77743521432033820.898645701248042.04062555358302−9.67733289275361
238.344191834031230.27817441678260−5.55699346407081831.042011855547989.23553399933964−3.85660415833079
24−9.002143968240302.942617187790473.11174899152034847.94146024657417−1.16315196293439−5.97586605514242
254.332159246809728.94616766180038−1.8581171846721185−0.10384830725595−2.048507366202209.86434844241350
26−8.593448948880644.593427810586632.1574381847950286−2.22329768993876−1.155008663514929.74507341073813
27−4.960573645916418.470293371153012.18057426683214872.454916663529245.78849624684043−7.78007284896142
28−3.22871530473500−0.06432575276105−9.3282980270683888−5.687749040808771.63805685551568−7.90852031663867
29−0.41110361489163−2.679928727983309.69222845496206895.590406684628465.338756293138616.61236943075913
30−4.689038998874977.258896447856515.18088675419929900.438187793809638.97196368206258−4.50152918880770
319.412910699909090.66348567366363−3.49891347368823912.03348649684924−0.050042587833919.91015953014887
320.94818881511365−6.65575236452663−7.22396030057950920.422406707996931.11071074115523−9.84357670621854
33−1.67590490013474−9.354558041289372.83086859870856938.59284720462050−4.75746898722982−1.96799907697002
34−4.169891829167017.860702841380044.7410778593883594−7.756459647999962.114083930738905.94974621604521
353.25227913344670−6.587395337392496.8336460762108795−6.573870311070057.250239130102772.18173076947982
362.794969594260608.194660924076495.31668384613840960.269166520600420.2558960452081510.09736113827280
378.93129544156749−0.38889049269008−4.5657555613609197−3.46440861622413−5.006147499072897.92447755025208
383.18101226342569−1.52942164932200−9.27296903448642987.75906351118220−1.287798774990836.37796001575860
397.042424717129402.41100805551266−6.7163138900011099−9.50971105265897−2.34766336070543−1.18936949313789
40−4.846198560080498.30681576997033−2.76067603756258100−7.842484255640784.566846611867314.21877363426257
41−4.85601964526380−5.430987532808626.796724190907491010.545936318116520.03153202952144−9.88981746265683
42−8.124209493316874.47991361320555−3.530270842368351020.116573537632620.5799327039053810.08847086418930
438.47199163087032−5.08576199441910−1.637164688891581034.283769531673797.700726957579265.07243383167886
44−1.490643008956799.97071885482340−0.09724232947774104−1.842550951464510.60951826428967−9.69627313630365
452.608427722563433.60002146556946−8.925399719753591051.714285022939595.137231788932668.58649396213649
46−3.83790555020136−9.074138493539800.672792905830421067.51446382502934−6.59902638575811−0.28577127558744
47−0.02360421418595−8.77097210922236−4.51418638304886107−7.041813536465215.453514380972044.60937233528599
48−9.33039977030179−1.145754920212223.184751785064611086.124873585267073.77291668319337−6.98594726658885
496.95165799200241−6.229004170997283.705367942718001094.537456584573340.852098465156109.02991192334709
503.218240469951967.570460713903435.97101496322205110−1.12776216963448−1.821495943472249.83655759421291
51−4.226340345629554.975677521365817.68357899479138111−3.63485132545628−2.96062615963476−8.65691411102701
523.133664137053994.08476978070495−8.555555710927441120.15207573590189−9.50293289978279−2.68947447382110
53−8.92065921311551−4.00273826814577−1.23991189847763113−5.67454541558127−5.354414909870126.17486513203551
54−1.278263436217825.652585779822058.301779139802711142.462582579266359.580637476534752.22976234541775
553.42467444593053−5.91207099422011−7.166480729289371154.151214288254380.208956630273399.24198564645230
564.421585682244051.50997493058889−8.807074024069291165.016260773640174.647503594670247.52635676060415
57−2.520274700994935.12888276079702−8.136777201641701177.145533502138797.113307225695211.18355483513316
58−7.66578663482778−4.23480394172058−4.47188682700721118−3.658142818454176.93857742396894−6.15382574716082
592.107660418546427.07201098260663−6.781894862101311196.912172391256440.228703360106117.41968187538919
603.626753323281404.249150248483788.48729188985141120−4.21831012911677−0.73580262084969−8.88074058906206
Table A4. Simulation dataset C’s coordinate data (units: mm).
Table A4. Simulation dataset C’s coordinate data (units: mm).
No.XYZNo.XYZ
15.95996262889717−2.53087159418200−7.6206281090066261−3.51246706823379−9.341825058616130.62737974656962
2−9.85107314551679−1.30444876074718−1.12023703431943624.50473221625104−6.052034990457176.56362954642608
3−0.425620095836071.60551591345338−9.86110403170941633.941835364309417.279245188170525.61030276751094
47.15248741603897−5.79696049895956−3.90354930806442641.105470408576357.49865115145544−6.52295453806957
5−8.798100845682034.724709327226560.52086306669802650.288265245957275.08233406156181−8.60741492035918
62.28942095507661−3.32448286731950−9.1491182829470166−9.33209675238265−1.88494780973294−3.05933742850091
7−0.236485938173690.53665373133653−9.98283913404341671.24171798290095−4.99923896512844−8.57125197473802
81.209685611716752.773181212662799.5313314525438068−1.22194940877846−1.79255182349007−9.76189481363779
95.93420067374265−3.08518708773780−7.4342339931756869−2.336589862929160.408349115635869.71465342116742
10−1.83054232077633−0.60978790824329−9.8121168734577970−2.507914521602769.322682332967922.60746212216540
11−6.42127372632781−5.52810041703587−5.31108701348930717.43197381777530−4.30470977913549−5.12213226028162
121.198489173514916.75231597828767−7.27810776457924724.236797188036026.885749582028495.88530467085766
13−0.425831973379066.133956109929087.8862826966555173−3.353336708143019.24977931999824−1.78798670354391
146.57382699453694−4.698703395788705.89130064077983749.71840199364162−1.82443538824880−1.49161565565493
15−0.672056022615953.154166288649229.4657121369517775−8.191935519956401.67782531605284−5.48429567206496
160.271554343874970.27118450797139−9.99266903431794760.40140787070271−1.33286253492707−9.90266817867122
170.25308083307023−0.04010002834001−9.99676658135817774.537260532605640.94151974265130−8.86154770747293
18−1.17785831771628−0.24510239526266−9.92740700432568781.615168133500646.96306357529998−6.99336339363929
197.065363922071052.215102206239396.72121934542273799.25819273048234−1.16503736461130−3.59564896620229
20−5.311281654706661.05688392418140−8.4067847756139580−6.28584646497884−2.28309780959583−7.43482055476037
21−2.11065116617518−4.16503976045407−8.8429577213881381−0.191789745868900.09521260975715−9.99775058306743
222.15886526922128−8.547799925826574.71964413912382828.50066661341672−2.401792197865824.68722320418105
235.88541066182448−1.27357658674910−7.9837649019267783−6.41886291974595−6.690209249917843.74696167129575
24−1.853898229824941.07643711617595−9.7675320001769484−2.17639276771440−9.70204539665501−1.06518995694937
257.271170647784136.83385042030464−0.6553375905380785−1.249020714423489.915816139867890.34144617943944
26−5.719546214365504.100209775916567.10462799391528864.009101138085789.088623619575771.15069225193287
275.08401851378667−0.029606723361428.6111765490215887−3.79323507426728−8.78984680257001−2.88972600375784
28−0.26443730535060−0.030538473300099.99647246809361887.372070445484416.549965488819431.65869014847620
299.07084930001875−3.86540812389358−1.66690047914021894.37274176835174−0.22139490848028−8.99057566771884
302.86805392843574−5.57345457408119−7.7917477926348390−9.20458184798071−2.74927686394480−2.77804077503713
31−3.95616029511580−1.76938209469346−9.01216304362207913.84083333935019−1.51569516697420−9.10772569965450
322.368185127141666.745258082133026.9923924296263892−2.59334387393980−0.262391267247169.65431239809196
336.08354542163128−5.08223798323700−6.09601467195711930.78419338417084−3.04987739990699−9.49124240726619
348.51756401769817−2.59707326287336−4.55050613105236949.006765797471053.265128395185452.86659322301019
35−0.001508491377140.0017661746306010.0000199739768095−2.929026317778723.69839622648143−8.81720186526334
36−8.87466688583470−4.100786072000662.10329295561852960.98664286932946−8.90409024626217−4.44345701534720
37−2.035490572947966.117017736945907.64453217229997977.86678077709667−6.13458879926826−0.69374280606344
388.537059445814644.35922961278074−2.8489880050908098−1.157258009797532.631288959767839.57794718658249
39−9.115281661035954.06506417766708−0.6227956137797099−0.14579857774420−0.21121056560729−9.99673047090975
40−0.581021841655769.61081481618440−2.70102749001126100−1.25580657367142−0.10475745442277−9.92030984831879
410.353210823774308.17644980500951−5.746456501084881017.01763907424743−0.69202363480080−7.09047566210989
425.67461369627328−4.99499972058684−6.54595569871264102−0.653886179863572.57232761613757−9.64137918690304
432.795532003967059.411251527211341.901022782297291036.043985300352037.30171439964982−3.18676666787718
44−5.761169668670948.168290198585500.29803366260499104−5.39008334649265−2.990990114411927.87410240772113
450.81231109499248−1.761305917244669.810105615188311056.003759433053084.59441140239776−6.54573283811845
466.78433265308395−2.49420282986743−6.910321291060411069.044075507723261.533040349178103.98189176899673
471.864202861587899.164311006808143.54120762753932107−4.97523299033828−2.57612725498324−8.28317226619606
480.463102930390699.736968128662902.231068681023051082.276480315340532.92945614989988−9.28636253779010
491.88555388454493−2.55997340072303−9.481136948224411091.720837372865126.52616757239851−7.37887901737959
50−6.56595707040798−4.98992000387135−5.65595158951884110−2.78810706129170−1.41914467735619−9.49804109497103
510.609309580041520.31233873762694−9.976561495260661116.13684273879640−3.07613391500037−7.27163312280270
524.81820410624501−5.81186464235642−6.55805905523815112−2.652065200946118.059988973511775.29183351617852
53−6.92416985159192−6.38950809436772−3.351224426753021137.65293149411519−6.408403657059490.60418841047588
54−4.01794200780133−5.597757537768357.247154791385581140.564865451001153.08894380518371−9.49422379061154
55−9.79487189140082−1.04454957796157−1.72335488751927115−1.04130724007516−1.53503208605850−9.82651289771900
562.52865790796924−5.60942584080151−7.88293580244636116−4.235617519075410.832626795155599.02037958499682
57−7.466476407114605.76777831438991−3.31443258802386117−1.361568736098513.967803529243889.07759140577493
581.60575217539826−4.852362051718498.595148018577701183.01460153226631−2.371591029893489.23514540246705
590.45562100548285−1.85585394689852−9.81574065321849119−9.085387753841363.961727350362341.32686651808811
600.608986779435334.94763038814187−8.66894807615183120−1.13534715705872−5.15041230845435−8.49618460850138
Table A5. Simulation dataset D’s coordinate data (units: mm).
Table A5. Simulation dataset D’s coordinate data (units: mm).
No.XYZNo.XYZ
11.01827397046852−16.58223949132140−10.89363718109050616.0638323945852418.620295491528004.72943319369740
2−13.34805033351610−7.65151186352272−12.5120736460994062−17.4926872840853−9.37509388839011−0.74855543825462
3−10.85895498342310−7.82787277400729−14.6326215912328063−2.9505082162349013.32859826311390−14.5866933276464
4−7.8268820469204018.45532417824060−0.3967957359506364−1.46204046352872−0.8133028842110920.01798450174830
510.30900200779320−10.09332954684140−13.75181254401260650.6169412017550914.37587094279790−13.8975848364538
6−14.43968045809420−2.42937674159383−13.39815874219520667.22220058167960−1.5111052769033118.71931724681080
7−8.7062653311866315.720649586150808.9569404727011467−6.388496996666280.99329051206819−18.7971296197525
8−9.46047435993298−4.3311846480302817.0988803648163068−11.10926788135752.70467302295429−16.2575098727216
9−8.55574557041805−3.6181382714638617.742389759172806919.176247882951503.180084007788015.25848048259532
10−4.4757535874048412.0239404834543015.49097535747460700.574003322982363.18735409862370−19.6545834053648
111.18157843626710−19.84688771532510−0.8758353009275571−8.84843005167697−10.9152758645158−13.9921756195152
1212.72462790755350−1.79655645813156−15.2955263880729072−3.001773439171240.17099678521819−19.6578936374463
13−7.356433505547956.18067041315640−17.4335976543294073−12.789645282877610.09334240088980−11.4754176273606
141.16280745257179−0.7207894883590920.054870575952207417.90839450288090−3.541566581351828.44268197993650
15−0.2081310483326414.5820863299503013.89037089100960752.9610659620201616.529812923765611.13973450479910
1618.253677358890107.116080710636074.7063683729416176−13.14656178906257.62914121276286−12.8553476744905
1711.49696314639260−15.441658035262205.4447301899657377−8.7464606559119111.632294657624013.83850132902650
18−2.163924690890010.83039496171992−19.75803535584460782.7996815357340319.3874358147267−4.45045599529497
191.683556330128951.8862005487455919.9570392097982079−1.484964042859760.5938852930200420.0309817973941
2019.643948305824601.215318891976134.19761492813433801.86825643905109−4.75766614049884−19.2203666266285
213.911444260053753.7946234859560419.3827033879984081−0.83656074379294−0.96746618861219−19.8495205433453
2218.554850353350306.58707162584348−4.0637055097497582−5.4949264466146910.73418866681580−15.8878462959520
23−14.97456160827260−1.11066147552205−12.9878643226808083−3.800932202007703.13541687496168−19.2796042424803
243.445858016758399.94215773899642−16.9861805320120084−15.58808235190499.40106327444076−8.10819659849956
25−6.09016836774572−2.68963529620687−18.7122246049292085−18.25765470319611.128956404882747.97006084712438
26−13.04807418539690−12.549476682658408.292617046095418610.4953657178619013.79659705098480−10.1146838692175
272.09518671398468−1.2722456951722219.95286733485780878.80543241338193−17.6458457320953−2.94971131041900
284.7545518948881216.99344821169060−9.54088810370740889.25865041278559−1.10120558793086−17.6391233771429
2919.475230879797400.94602374245893−4.7872353759850889−17.85055620228704.57435019426862−7.49985287327967
30−9.719394839687721.8099935622106317.4394220555959090−10.25239541889853.66547938695731−16.6366468132108
314.23624746745487−16.37834724788520−10.452512917124909111.56502410835130−11.1667931585712−11.8001152889516
321.9412209156761013.9845728829226014.3768539657455092−11.4073320118934−13.9933247567537−8.20420717407658
33−5.03249637040109−6.5126803071895618.263977139149809310.68001801215560−15.4856240086703−6.61951623427540
34−1.151799703496946.29548586292803−18.87496615486860945.90864697584834−17.34228352775267.97570731093767
351.648980275939599.77027856807626−17.33808038472760950.4666772296990320.02206521434100−1.62541336683918
36−2.056293561249589.4820996348847617.6307354191773096−10.96159421603574.4719192940741316.17879762919220
373.854113797783981.69588783533586−19.47958316498940972.28810905920513−3.44111150716674−19.4620277927884
3814.15697204354040−13.90269899135140−2.415225223480139816.884017412668710.08724836843570−4.20718583643146
39−0.204269920026110.5926367329620720.0916167409768099−11.2759557549732−5.1843391253262215.6776614308102
40−2.525445442604644.92288722155417−19.13141167663400100−5.8232810574598513.770783830111013.4425890803288
41−4.81520178654368−2.5833188493722319.300025160486201011.81833972088083−11.227039842772216.4934610087666
42−0.89198658551585−10.55279166077970−16.79776256611030102−11.316955006585415.5247657046385−5.53294932892596
43−11.52970092308250−12.92085836386540−9.6574357442062910318.497176362380005.51582670251050−5.57565972709238
44−4.3068791107027719.03743999244430−4.58540004091944104−17.5915674679496−9.22534652387289−0.11610300895466
45−1.79046902193470−14.74656897912440−13.1667214708170010515.51425837266370−4.5332252312352311.97156707998060
46−0.82428299505741−1.48334775795102−19.815800196845101061.8060231613819516.5792733261551011.30372090674270
4717.979581171109706.266417548701416.6035364250342610718.880282532438106.779162761414181.74012279350019
48−11.55695697104530−16.104529861126901.336487726417151082.44185338612605−11.5128302064991−16.0139102367006
49−6.51472496606937−9.82199878493186−15.9561400324750010910.96228303494210−14.4365979646038−8.30823285377158
50−10.628437628916305.6569058585682616.0380525605356011014.934181888356107.4999044689972711.28877720486080
511.49454173751660−2.95598487258305−19.61588724356240111−2.00751746160698−4.96094780607129−19.1339669360993
526.235608286731154.92507746739824−18.31393326764820112−19.8774126579963−0.34144282531075−0.74157929219974
5318.506885448227507.85024987758855−0.958559261928271130.5978161433241417.79372679846470−9.21044241612647
5417.03503404541520−1.2920086282884310.64832253787790114−6.7284839648613714.21862165584770−12.3111457743369
55−4.99871766746520−16.35724069453180−10.0568959065283011518.403469334180506.73345294480625−4.48043076896942
5617.46517828652930−9.81754829167634−0.204860781865851169.6839095094262313.1301170821182011.86285452827500
57−12.386158980153105.4236998373396114.78857190219600117−9.07980486558268−2.47679609248717−17.4809966968707
5814.0156120246696012.06231488177630−7.8535505162328611815.532689896697705.5674083581117111.58649332044970
59−1.087527961699600.04350469111739−19.86463520577670119−6.313954384983430.32145413074370−18.8423392515690
6017.268575071638909.805991818610713.42176555650988120−1.1038210982803015.73986163459290−12.3074793094377

References

  1. Everitt, C.W.F.; Lipa, J.A.; Siddall, G.J. Precision engineering and Einstein: The Relativity Gyroscope experiment. Precis. Eng. 1979, 1, 5–11. [Google Scholar] [CrossRef]
  2. Lipa, J.A.; Siddall, G.J. High precision measurement of gyro rotor sphericity. Precis. Eng. 1980, 2, 123–128. [Google Scholar] [CrossRef]
  3. Stojadinovi, S.M.; Majstorovi, V.D. An Intelligent Inspection Planning System for Prismatic Parts on CMMs; Springer International Publishing: Cham, Switzerland, 2019. [Google Scholar]
  4. Pathak, V.K.; Singh, R. A Comprehensive Review on Computational Techniques for Form Error Evaluation. Arch. Comput. Method Eng. 2022, 29, 1199–1228. [Google Scholar] [CrossRef]
  5. Acko, B.; Mccarthy, M.; Haertig, F.; Buchmeister, B. Standards for testing freeform measurement capability of optical and tactile coordinate measuring machines. Meas. Sci. Technol. 2012, 23, 94013–94025. [Google Scholar] [CrossRef]
  6. ISO 1101:2017; Geometrical Product Specifications (GPS), Geometrical Tolerancing, Tolerances of Form, Orientation, Location and Run-out. ISO: Geneva, Switzerland, 2004.
  7. Cogorno, G.R. Geometric Dimensioning and Tolerancing for Mechanical Design; McGraw-Hill Education: New York, NY, USA, 2020. [Google Scholar]
  8. Samuel, G.L.; Shunmugam, M.S. Evaluation of sphericity error from coordinate measurement data using computational geometric techniques. Comput. Methods Appl. Mech. Eng. 2001, 190, 6765–6781. [Google Scholar] [CrossRef]
  9. Liu, F.; Xu, G.; Zhang, Q.; Liang, L.; Liu, D. An intersecting chord method for minimum circumscribed sphere and maximum inscribed sphere evaluations of sphericity error. Meas. Sci. Technol. 2015, 26, 115005. [Google Scholar] [CrossRef]
  10. Zheng, Y. A simple unified branch-and-bound algorithm for minimum zone circularity and sphericity errors. Meas. Sci. Technol. 2020, 31, 045005. [Google Scholar] [CrossRef]
  11. Fanwua, M.; Chunguang, X.U.; Haiming, L.I.; Juan, H. Sphericity evaluation using maximum inscribed sphere method. Procedia Eng. 2011, 24, 737–742. [Google Scholar] [CrossRef]
  12. Meng, F.W.; Xu, C.G.; Li, H.M.; Hao, J. The Relationship between Minimum Zone Sphere and Minimum Circumscribed Sphere and Maximum Inscribed Sphere. Appl. Mech. Mater. 2012, 157–158, 658–665. [Google Scholar] [CrossRef]
  13. Meng, F.W.; Xu, C.G.; Hao, J.; Xiao, D.G. Sphericity evaluation based on minimum circumscribed sphere method. Adv. Mater. Res. 2012, 433, 3146–3151. [Google Scholar] [CrossRef]
  14. Fei, L.; Guanghua, X.U.; Lin, L.; Qing, Z.; Dan, L. Intersecting Chords Method in Minimum Zone Evaluation of Sphericity Deviation. J. Mech. Eng. 2016, 52, 137–143. [Google Scholar]
  15. Fana, K.C.; Lee, J.C. Analysis of minimum zone sphericity error using minimum potential energy theory. Precis. Eng. 1999, 23, 65–72. [Google Scholar] [CrossRef]
  16. Liu, F.; Xu, G.; Liang, L.; Zhang, Q.; Liu, D. Minimum zone evaluation of sphericity deviation based on the intersecting chord method in Cartesian coordinate system. Precis. Eng. 2016, 45, 216–229. [Google Scholar] [CrossRef]
  17. Janecki, D.; Stępień, K.; Adamczak, S. Sphericity measurements by the radial method: I. Mathematical fundamentals. Meas. Sci. Technol. 2015, 27, 015005. [Google Scholar] [CrossRef]
  18. Janecki, D.; Stępień, K.; Adamczak, S. Sphericity measurements by the radial method: II. Experimental verification. Meas. Sci. Technol. 2015, 27, 015006. [Google Scholar] [CrossRef]
  19. Lei, X.; Gao, Z.; Duan, M.; Pan, W. Method for sphericity error evaluation using geometry optimization searching algorithm. Precis. Eng. 2015, 42, 101–112. [Google Scholar]
  20. Kanada, T. Evaluation of spherical form errors—Computation of sphericity by means of minimum zone method and some examinations with using simulated data. J. Jpn. Soc. Precis. Eng. 1995, 17, 281–289. [Google Scholar] [CrossRef]
  21. Huang, J.; Jiang, L.; Zheng, H.; Xie, M.; Tan, J. Evaluation of minimum zone sphericity by combining single-space contraction strategy with multi-directional adaptive search algorithm. Precis. Eng. 2019, 55, 189–216. [Google Scholar] [CrossRef]
  22. Mao, J.; Zhao, M. An Approach for the Evaluation of Sphericity Error and Its Uncertainty. Adv. Mech. Eng. 2013, 2013, 67–77. [Google Scholar] [CrossRef]
  23. Cha, Y.; Agrawal, A.K.; Kim, Y.; Raich, A.M. Multi-objective genetic algorithms for cost-effective distributions of actuators and sensors in large structures. Expert Syst. Appl. 2012, 39, 7822–7833. [Google Scholar] [CrossRef]
  24. Cha, Y.; Kim, Y.; Raich, A.M.; Agrawal, A.K. Multi-objective optimization for actuator and sensor layouts of actively controlled 3D buildings. J. Vib. Control 2013, 19, 942–960. [Google Scholar] [CrossRef]
  25. Cha, Y.J.; Raich, A.; Barroso, L.; Agrawal, A. Optimal placement of active control devices and sensors in frame structures using multi-objective genetic algorithms. Struct. Control. Health Monit. 2013, 20, 16–44. [Google Scholar] [CrossRef]
  26. Cui, C.C.; Che, R.S.; Ye, D. Sphericity error evaluation using the genetic algorithm. Opt. Precis. Eng. 2002, 10, 333–339. [Google Scholar]
  27. Wen, X.; Li, H.; Wang, F.; Wang, D. Sphericity error united evaluation using particle swarm optimization technique. In Proceedings of the 2009 9th International Conference on Electronic Measurement & Instruments, Beijing, China, 16–19 August 2009; IEEE: New York, NY, USA, 2009; pp. 1-156–1-161. [Google Scholar]
  28. Wen, X.; Song, A. An immune evolutionary algorithm for sphericity error evaluation. Int. J. Mach. Tools Manuf. 2004, 44, 1077–1084. [Google Scholar] [CrossRef]
  29. Jie, H.; Changcai, C.; Fugui, H. Sphericity error evaluation based on a modified particle swarm optimizer using dynamic inertia weight. J. Graph. 2012, 33, 99–103. [Google Scholar]
  30. Yu, X.; Huang, M.; Xia, P. Sphericity Error Evaluation Based on an Improved Particle Swarm Optimization. Comput. Syst. Appl. 2009, 12, 201–203. [Google Scholar]
  31. Jiang, L.; Huang, J.; Ding, X.; Chao, X. Method for spherical form error evaluation using cuckoo search algorithm. In Proceedings of the 10th International Symposium on Precision Engineering Measurements and Instrumentation (ISPEMI 2018), Kunming, China, 8–10 August 2018. [Google Scholar]
  32. Xiulan, W.; Aiguo, S. An improved genetic algorithm for sphericity error evaluation. In Proceedings of the International Conference on Neural Networks and Signal Processing, 2003. Proceedings of the 2003, Nanjing, China, 14–17 December 2003; IEEE: New York, NY, USA, 2003; pp. 549–553. [Google Scholar]
  33. Chen, M. Analysis of spherical form errors to coordinate measuring machine data. JSME Int. J. Ser. C 2002, 45, 647–656. [Google Scholar] [CrossRef]
  34. Huang, J. An exact minimum zone solution for sphericity evaluation. Comput.-Aided Des. 1999, 31, 845–853. [Google Scholar] [CrossRef]
  35. Calvo, R. Sphericity measurement through a new minimum zone algorithm with error compensation of point coordinates. Measurement 2019, 138, 291–304. [Google Scholar] [CrossRef]
  36. Wang, M.; Cheraghi, S.H.; Masud, A.S. Sphericity error evaluation: Theoretical derivation and algorithm development. IIE Trans. 2001, 33, 281–292. [Google Scholar] [CrossRef]
  37. Cui, C.; Fan, W.; Huang, F.; Li, B. An Iterative Method for Minimum Zone Sphericity Error Evaluation. Iran. J. Mater. Sci. Eng. 2011, 1, 367. [Google Scholar]
Figure 1. Geometric structure of sphericity error evaluation based on the minimum zone criterion.
Figure 1. Geometric structure of sphericity error evaluation based on the minimum zone criterion.
Applsci 13 13144 g001
Figure 2. Schematic representation of spherical error distribution characteristics. (a) Spherical error distribution of sampling points in the global space; (b) change in sphericity along the X-axis; (c) change in sphericity along the Y-axis; (d) change in sphericity along the Z-axis; (e) schematic distribution of sampling points in Local Space A; (f) schematic distribution of sampling points in Local Space B; (g) schematic distribution of sampling points in Local Space C.
Figure 2. Schematic representation of spherical error distribution characteristics. (a) Spherical error distribution of sampling points in the global space; (b) change in sphericity along the X-axis; (c) change in sphericity along the Y-axis; (d) change in sphericity along the Z-axis; (e) schematic distribution of sampling points in Local Space A; (f) schematic distribution of sampling points in Local Space B; (g) schematic distribution of sampling points in Local Space C.
Applsci 13 13144 g002
Figure 3. Flowchart of the algorithm.
Figure 3. Flowchart of the algorithm.
Applsci 13 13144 g003
Figure 4. Local solution spaces’ sizes and distribution.
Figure 4. Local solution spaces’ sizes and distribution.
Applsci 13 13144 g004
Figure 5. Schematic representation of the computational geometric method for local solution spaces.
Figure 5. Schematic representation of the computational geometric method for local solution spaces.
Applsci 13 13144 g005
Figure 6. Schematic diagram of the search principle in a one-dimensional space.
Figure 6. Schematic diagram of the search principle in a one-dimensional space.
Applsci 13 13144 g006
Figure 7. The distribution of dataset 1.
Figure 7. The distribution of dataset 1.
Applsci 13 13144 g007
Figure 8. Algorithm optimization process (corresponding experiments based on dataset 1).
Figure 8. Algorithm optimization process (corresponding experiments based on dataset 1).
Applsci 13 13144 g008
Figure 9. The distribution of dataset 2.
Figure 9. The distribution of dataset 2.
Applsci 13 13144 g009
Figure 10. Algorithm optimization process(corresponding experiments based on dataset 2).
Figure 10. Algorithm optimization process(corresponding experiments based on dataset 2).
Applsci 13 13144 g010
Figure 11. The distribution of dataset 3.
Figure 11. The distribution of dataset 3.
Applsci 13 13144 g011
Figure 12. Algorithm optimization process (corresponding experiments based on dataset 3).
Figure 12. Algorithm optimization process (corresponding experiments based on dataset 3).
Applsci 13 13144 g012
Figure 13. The distribution of dataset 4.
Figure 13. The distribution of dataset 4.
Applsci 13 13144 g013
Figure 14. Algorithm optimization process (corresponding experiments based on dataset 4).
Figure 14. Algorithm optimization process (corresponding experiments based on dataset 4).
Applsci 13 13144 g014
Figure 15. The distribution of simulated datasets. (a) Simulated dataset A; (b) simulated dataset B; (c) simulated dataset C; (d) simulated dataset D.
Figure 15. The distribution of simulated datasets. (a) Simulated dataset A; (b) simulated dataset B; (c) simulated dataset C; (d) simulated dataset D.
Applsci 13 13144 g015
Figure 16. Measurement of a metal hemispherical shell resonator surface with a Hexagon CMM. (a) The Hexagon CMM and environment. (b) A metal hemispherical shell resonator sample was used for the spherical surface measurement.
Figure 16. Measurement of a metal hemispherical shell resonator surface with a Hexagon CMM. (a) The Hexagon CMM and environment. (b) A metal hemispherical shell resonator sample was used for the spherical surface measurement.
Applsci 13 13144 g016
Figure 17. The dataset distribution.
Figure 17. The dataset distribution.
Applsci 13 13144 g017
Figure 18. Algorithm optimization process.
Figure 18. Algorithm optimization process.
Applsci 13 13144 g018
Table 1. Comparison with results based on dataset 1 in the literature (rounding is as given in the respective papers; units: mm).
Table 1. Comparison with results based on dataset 1 in the literature (rounding is as given in the respective papers; units: mm).
AlgorithmSpherical Center Coordinate (X, Y, Z)Sphericity Error
Least squares(−0.01913, 0.08935, 0.01762)1.15224
[34](0.00000, 0.00000, 0.00000)1.00000
[21](0.00000000, 0.00000000, 0.00000000)1.00000000
Ours(−5.3 × 10−14, 0.1 × 10−14, 1.6 × 10−14)1.00000000000003
The number of zeros after the decimal points represents the precision provided in the original data from the cited references.
Table 2. Results of ten experiments based on dataset 2 (units: mm).
Table 2. Results of ten experiments based on dataset 2 (units: mm).
Experiment1234
Sphericity Error0.015384870549310.015384870549270.015384870549340.01538487054926
Experiment6789
Sphericity Error0.015384870549280.015384870549300.015384870549350.01538487054926
Experiment510
Sphericity Error0.015384870549290.01538487054932
Table 3. Comparison with results based on dataset 2 in the literature (rounding is as given in the respective papers; units: mm).
Table 3. Comparison with results based on dataset 2 in the literature (rounding is as given in the respective papers; units: mm).
AlgorithmSpherical Center Coordinate (X, Y, Z)Sphericity Error
Least squares(−0.0001146, −0.0002871, 0.0078932)0.01639862
[35](−0.0011426, 0.0000027, 0.0107779)0.0154077
[16](0.000179, −0.000332, 0.011747)0.015385
Ours(0.0002125983629424, −0.0003506987099171, 0.01174819302036)0.01538487054931
Table 4. Results of ten experiments based on dataset 3 (units: mm).
Table 4. Results of ten experiments based on dataset 3 (units: mm).
Experiment1234
Sphericity Error3.332517616163553.332517616163543.332517616163553.33251761616355
Experiment5678
Sphericity Error3.332517616163573.332517616163553.332517616163563.33251761616356
Experiment910
Sphericity Error3.332517616163553.33251761616355
Table 5. Comparison with results based on dataset 3 in the literature (rounding is as given in the respective papers; units: mm).
Table 5. Comparison with results based on dataset 3 in the literature (rounding is as given in the respective papers; units: mm).
AlgorithmSpherical Center Coordinate (X, Y, Z)Sphericity Error
Least squares(−0.2545316, −0.6600826, 0.0380931)3.735706
[27](−0.388729, −0.355488, −0.299887)3.332518
[37](−0.38873, −0.355488, 0.299888)3.33252
[26](−0.412356, −0.335014, −0.326140)3.351375
[36](−0.38872967, −0.35548811, −0.29988752)3.33251813257
[35](−0.3887296, −0.3554883, −0.2998876)3.33251762
Ours(−0.38872964944657, −0.35548834466725, −0.29988759232766)3.33251761616355
Table 6. Results of ten experiments based on dataset 4 (units: mm).
Table 6. Results of ten experiments based on dataset 4 (units: mm).
Experiment1234
Sphericity Error0.007660194937850.007660194937910.007660194937910.00766019493788
Experiment6789
Sphericity Error0.007660194937890.007660194937870.007660194937880.00766019493784
Experiment510
Sphericity Error0.007660194937870.00766019493786
Table 7. Comparison with results based on dataset 4 in the literature (rounding is as given in the respective papers; units: mm).
Table 7. Comparison with results based on dataset 4 in the literature (rounding is as given in the respective papers; units: mm).
AlgorithmSpherical Center Coordinate (X, Y, Z)Sphericity Error
Least squares(0.0008555, −0.0003732, 0.0007853)0.00848634
[26](0.000993, 0.000024, 0.000058)0.007928
[28](0.002495, −0.000097, 0.000479)0.007660
[10](0.002504156298, −0.000096126920, 0.000481529278)0.007660194938
Ours(0.00250415629838, −0.00009612691940, 0.00048152927799)0.00766019493788
Table 8. Detailed information on the simulated datasets (units: mm).
Table 8. Detailed information on the simulated datasets (units: mm).
Simulation DatasetSpherical Center Coordinate (X, Y, Z)Sphericity ErrorRadiusNumber of Coordinates
A(0, 0, 0)0.0000110120
B(0.1, 0.1, 0.1)0.0000110120
C(0, 0, 0)0.0000510120
D(0, 0, 0)0.0000120120
Table 9. Experimental results based on simulated datasets (units: mm).
Table 9. Experimental results based on simulated datasets (units: mm).
Simulation DatasetSpherical Center Coordinate
(X, Y, Z)
Sphericity ErrorStandard Deviation of ResultsDeviation from the Ideal Sphericity Error
A(6.3 × 10−16, −1.5 × 10−15, −2.8 × 10−15)0.000010000000010.2 × 10−141 × 10−14
B(0.100000000000001, 0.099999999999994, 0.100000000000000)0.000010000000091.0 × 10−149 × 10−14
C(−3.5 × 10−15, −3.9 × 10−15, −1.6 × 10−15)0.000050000000091.8 × 10−149 × 10−14
D(5.8 × 10−14, −3.7 × 10−14, −6.1 × 10−14)0.000010000000050.9 × 10−145 × 10−14
Table 10. Ablation experiment results (units: mm).
Table 10. Ablation experiment results (units: mm).
DatasetAlgorithmMean Value of ResultsStandard Deviation of Results
1A1.000577194293391.7 × 10−3
B1.000000000000032.0 × 10−14
2A0.015395049196697.6 × 10−6
B0.015384870549313.0 × 10−14
3A3.395530118622619.5 × 10−2
B3.332517616163550.8 × 10−14
4A0.007660805880455.5 × 10−7
B0.007660194937880.2 × 10−14
Table 11. Results of ten experiments based on measurement data (units: mm).
Table 11. Results of ten experiments based on measurement data (units: mm).
Experiment1234
Sphericity Error0.007471884420660.007471884420720.007471884420670.00747188442064
Experiment6789
Sphericity Error0.007471884420670.007471884420660.007471884420640.00747188442065
Experiment510
Sphericity Error0.007471884420650.00747188442064
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, D.; Cui, J.; Wang, Z.; Sun, Y. A Novel Approach for High-Precision Evaluation of Sphericity Errors Using Computational Geometric Method and Differential Evolution Algorithm. Appl. Sci. 2023, 13, 13144. https://doi.org/10.3390/app132413144

AMA Style

Zhao D, Cui J, Wang Z, Sun Y. A Novel Approach for High-Precision Evaluation of Sphericity Errors Using Computational Geometric Method and Differential Evolution Algorithm. Applied Sciences. 2023; 13(24):13144. https://doi.org/10.3390/app132413144

Chicago/Turabian Style

Zhao, Dongfang, Junning Cui, Zhisheng Wang, and Yanxu Sun. 2023. "A Novel Approach for High-Precision Evaluation of Sphericity Errors Using Computational Geometric Method and Differential Evolution Algorithm" Applied Sciences 13, no. 24: 13144. https://doi.org/10.3390/app132413144

APA Style

Zhao, D., Cui, J., Wang, Z., & Sun, Y. (2023). A Novel Approach for High-Precision Evaluation of Sphericity Errors Using Computational Geometric Method and Differential Evolution Algorithm. Applied Sciences, 13(24), 13144. https://doi.org/10.3390/app132413144

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