Next Article in Journal
Safety Risk Assessment in Urban Public Space Using Structural Equation Modelling
Previous Article in Journal
CRISPR-Based Multi-Gene Integration Strategies to Create Saccharomyces cerevisiae Strains for Consolidated Bioprocessing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Self-Born Weighted Least Square Method for Cylindricity Error Evaluation

School of Mechanical Engineering, Shanghai Institute of Technology, Shanghai 201418, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(23), 12319; https://doi.org/10.3390/app122312319
Submission received: 24 October 2022 / Revised: 24 November 2022 / Accepted: 25 November 2022 / Published: 1 December 2022
(This article belongs to the Topic Manufacturing Metrology)

Abstract

:
In order to improve the stability of the evaluation results and the gross error resistance of the algorithm in view of the widespread gross errors in geometric error evaluation, an improved self-born weighted least square method (ISWLS) is proposed in this paper. First, the nonlinear cylindrical axial model is linearized to establish the error equation of the observed values. We use the conditional equations of the independent observations found as valid information to derive the weights of the observations. The weights of the observations are subjected to least-square iteration to calculate the error values and equation parameters. Meanwhile, the ordinal numbers of the independent sets of equations in the observed equations are updated several times. By updating the ordinal number information of the conditional equations, the influence of gross error data on the solution of the equations is minimized. Through a series of experiments, the algorithm is proved to have a strong resistance to gross differences, and operation time is shorter. According to the evaluation results of cylindricity error, the uncertainty of cylindricity error was calculated by the Guide to the expression of uncertainty in measurement method (GUM)and the Monte Carlo method (MCM). Experiments show that the uncertainty results of the MCM method can verify the results assessed by the GUM method, which proves that the results of the ISWLS method are effective and robust.

1. Introduction

Cylindricity error is an important basis for the acceptance of shaft parts. Accurate cylindricity error evaluation not only provides a reliable guarantee for improving the machining accuracy and assembly accuracy of parts but is also a prerequisite for stably improving production efficiency [1]. The spatial coordinate information of the measured point is measured, and the data are analyzed by an error evaluation algorithm to calculate the cylindricity error of the part [2]. The cylindricity error evaluation methods include the minimum zone cylinder method (MZC), least square cylinder method (LS), maximum inscribed cylinder method (MIC), and minimum circumscribed cylinder method (MCC). The minimum zone cylindrical method satisfies the minimum condition defined by the cylindricity error and is recognized as an arbitration method in case of inconsistent errors. Since it is an unconstrained nonlinear optimization problem, it cannot be solved directly by a computer. The least square method is widely used in the field of error evaluation by instruments such as the coordinate measuring machine (CMM), which has the advantages of mature theory, simple calculation, and stable evaluation results [3]. Since least square does not meet the minimum conditions defined by international standards, more in-depth research is required.
In recent years, many scholars have successfully applied the genetic algorithm [4], ant colony algorithm [5], artificial immune algorithm [6], artificial fish swarm algorithm [7], particle swarm algorithm [8], and other intelligent algorithms in geometric error assessment. In the field of geometric error evaluation, good results have been achieved. At the same time, more and more algorithms are proposed according to the characteristics of cylindricity error evaluation. He et al. [9] proposed a cylindricity error evaluation method based on the sequential quadratic programming algorithm, which uses coordinate simplification and the sequential quadratic programming algorithm to improve the accuracy of cylindricity error evaluation. Wu et al. [10] established an improved integrated learning particle swarm optimization algorithm applied to the cylindricity error evaluation problem to solve the local and globally optimal solutions of the cylindricity error evaluation results. Liu et al. [11] proposed a cylindricity error evaluation method based on incremental optimization. Li et al. [12] used the Hybrid Greedy Differential Evolution-Sine-Cosine Algorithm (HGSCADE) to solve optimization problems and evaluate cylindricity errors. Wang et al. [13] proposed a step-acceleration-based optimization algorithm to solve the efficiency problem of the crankshaft cylindricity error evaluation algorithm. The above methods for geometric error assessment have achieved good assessment accuracy, but the assessment methods are complicated, and further research and improvement are needed to improve the stability and robustness of the methods.
Based on the analysis of the observed data in production and experiments, statisticians observe that the probability of gross deviations is 1% to 10% of the total [14]. These gross errors seriously affect the error evaluation results. In the field of geodesy, to eliminate or attenuate the effect of gross errors on parameter estimation, Khaled et al. [15] proposed robust estimation to detect and eliminate gross errors at long distances. Guangfeng et al. [16] discussed the method of gross difference localization in detail: the “good” and “bad” points of gross difference are eliminated, but the “good” points are often eliminated as well. Jia et al. [17] studied the effectiveness of robust estimation in weakening and eliminating gross variances. The self-born weighted least square (SWLS) method is a robust estimation method with excellent robustness [18], which can effectively reduce the effect of gross differences, and the method is mainly applied in the field of mapping. It makes full use of the valid information provided by the conditional equations generated from independent observations to construct the weights of the observations, which can effectively eliminate or attenuate the effect of gross errors on parameter estimation. It is more effective than often other robust estimation methods [19]. At present, it has the advantages of simple theory and high accuracy of the assessment [20]. It has a good effect in reducing or eliminating gross errors, etc., but its research and application in the field of geometric error assessment have not been reported [21]. The GUM and MCM methods are widely used in uncertainty assessment. When the calculated results of the GUM method can be verified by MCM, the uncertainty evaluation results have high reliability [8]. In this paper, we consider the introduction of the SWLS method into the field of geometric error assessment and improve the adaptability of the SWLS method by using cylindricity error evaluation as an example. By combining model initialization with the screening of independent groups, the technique is extended from linear model assessment to a broader range of nonlinear error assessment, improving the adaptability, accuracy, and operational speed of the method. Meanwhile, the GUM method and MCM method were also used to calculate the cylindricity error uncertainty in this paper, which proved the reliability of the error evaluation results.

2. Model of Least Square Solution of Cylindricity Error

T Cylindricity error is the amount of deviation between the actual measured cylindrical surface and the ideal cylindrical surface, expressed as the difference between the radii of the two coaxial cylindrical surfaces. Assume that the measured coordinates in the right-angle coordinate system are P i x i , y i , z i , i = 1 , 2 , 3 , , n . The intersection point of the least square cylinder axis L and the cylinder interface is set as G a 1 , b 1 , c 1 . The direction vector is set as n = g , k , l . The equation of the axis L is:
x a 1 g = y b 1 k = z c 1 l
The axial L equation was simplified to obtain [13]:
x a p = y b q = z p = g l , q = k l a = a 1 c 1 × g l ,   b = b 1 c 1 × g l
The distance d i . from each sampling point on the cylindrical surface to the axis is:
d i = x i p × z i + a 2 + y i q × z i + b 2
The least square cylinder radius R is expressed as:
R = 1 n i = 1 n d i
Then the objective function of the least square cylindrical method is [4]:
J = min i = 1 n d i R 2  
The cylindricity error is transformed into a search for the corresponding ideal cylindrical surface axis parameters a , b , p , q , which minimize the objective function J. From this, the least square cylinder axis L equation is obtained, and the maximum distance R m a x and the minimum distance R m i n from the actual measuring point to the axis are obtained. The diagram is shown in Figure 1. The corresponding cylindricity error:
e = R m a x R m i n s

3. SWLS Solution for Cylindricity Error

3.1. Self-Born Weighted Least Square Robust Estimation Method

The SWLS method is a least square robust estimation method. It uses the observation correction number as the effective information provided by the conditional equation and constructs the weights with multiple estimates of the observation correction number. Compared with the ordinary least square method, the SWLS method uses the self-born weights of the observation value as the weights of the observation value, which effectively reduces the influence of gross errors.
SWLS assumes that the observation equation is of the linear form [18]:
L ^ = B x ^ + d + B X 0  
The error equation is:
V = B x ^ l  
B is the coefficient matrix of the observation equation, L ^ is the estimated value of the observation value L , V is the correct number of L , x ^ is the solution of the unknown, and d is the constant term matrix of the observation equation. The unknown matrix X ^ = X 0 + x ^ , where X 0 is the initial unknown value; l = d + B X 0 L .
The solution x ^ of the unknown, the observed number of corrections V , and the unit weighted variance valuation σ ^ 2 can be obtained by using the least square method. P is the weighted matrix; p ˙ j is the diagonal element in the weighted matrix with an initial value of 1 .
x ^ = B T P B 1 B T P l  
V = B T P B 1 B T P l l  
σ ^ 0 2 = V T P V r
The coefficient matrix B in Equation (8) is transformed into B t and B r , and B t is a t × t order full rank matrix; n . is the number of observations, and t is the number of parameters to be requested. Meanwhile, V and l are matrix transformed to obtain V t and l t :
V t = B t x ^ l t
V r = B r x ^ l r  
According to formulas (12) and (13), we can obtain:
x ^ = B t 1 V t + l t  
V r = B r t V t W r t
In the formula, B r t = B r B t 1 , W r t = B r B t 1 l t l r . Both B t and B r are not unique and are primarily chosen by independent groups.
Since the absolute value of the accidental error has a finite value, η σ ^ 0 is used to limit the unit weighted error to a certain estimation range, and η is the value range coefficient (take 3 according to the simulation experiment). From Equations (12) and (13), the m groups of true error estimates V 1 that satisfy Equation (16) are obtained, V 1 ,   V 2 , …, V m are:
V i = v 1 i   v 2 i     v n i T ,   i = 1 , 2 ,   , m
v j i η σ ^ 0 p j   i = 1 , 2 , , m ;   j = 1 , 2 , , n  
In the formula, when v j i V t , η σ ^ 0 p j is the initial value of v j i , η σ ^ 0 p j is the final value of v j i , and η σ ^ 0 θ p j is the step size when v j i V r , v j i is determined by Equation (15) [18]; p j is the diagonal element of the weighted matrix generated by the previous iteration.
The self-born variance is the variance of the observed value calculated from m estimates of the true error of the same observation value. The reproduction weights are the weights of the observation value calculated from the reproduction variance of the observation value. Then, the self-born variance σ ˙ j , self-born variance σ ¯ 0 , and self-born weights p ˙ j generated from different observations are calculated as follows:
σ ˙ j 2 = 1 m i = 1 m v j i 2   j = 1 , 2 , ,   n  
σ ¯ 0 2 = 1 n j = 1 n σ ˙ j 2  
p ˙ j = σ ¯ 0 2 σ ˙ j 2   j = 1 , 2 ,   , n  
The calculated self-born weight p ˙ j is taken into the weights of the observation value, and the correction number V of the observation value and the error σ ^ 0 in the unit weights are calculated by the least square method. The iteration terminates when the adjacent observation correction number iteration difference is less than the limit.

3.2. Insufficiency of Self-Born Weighted Least Square Method

SWLS is mainly applied in linear regression, and it uses analytical calculations to solve the error equation. It has not been studied for nonlinear problems and needs to be further optimized and explored. SWLS has certain requirements for the choice of independent groups in the observations and uses the number of observation corrections as valid information supplied by the conditional equation. It is necessary to artificially modify the selection of independent groups in various application settings. The outcomes will exhibit significant variances when the independent group was improperly chosen and includes serious flaws. In geometric error evaluation (e.g., CMM evaluation of geometric errors), producers need measurement methods that are more efficient, accurate and adaptable. The adjustment of the independent group will lead to an increase in production costs as well as a decrease in adaptability and accuracy, so the SWLS method is needed for further improvement.

3.3. Improved self-Born weighted Least Square Method

In response to the above shortcomings, this paper improves the SWLS method to increase the effectiveness of the method. The improved method is called ISWLS. In this paper, we take the calculation of cylindricity error as an example and apply the SWLS method to linearize the nonlinear model to calculate the cylindrical axis equation and cylindricity error. Additionally, an automatic selection step of the independent equation sets is implemented to lessen the impact of the choice of independent equation sets on the outcomes. It eliminates the manual filtering step of the independent group equations and makes the calculation simplified. The calculation steps are as follows:

3.3.1. Model Initialization

To apply the SWLS method to a nonlinear problem, the nonlinear model must first be transformed into a linear model. The model equations are assumed to be:
Z = F X  
In order to linearize the equation, a sufficient approximation X 0 of X ˜ is taken [19]:
X ˜ i = X i 0 + x ˜   ,   i = 1 , 2 , , n
Because   x ˜ is a tiny quantity, the Taylor formula can be used to omit the quadratic and more than quadratic terms [22], and then the functional equation is:
Z = F X 1 0 , X 2 0 ,   , X n 0 + F X 0 x ˜ + F X 1 x ˜ + F X 2 x ˜ + + F X n x ˜ + Δ higher   order   terms
Bringing Equations (3)–(5) into Equation (23), it is expressed in matrix form as:
J = b 1   b 2   b 3   b 4   b 5 δ a ^ δ b ^ δ p ^ δ q ^ δ R ^ D i
Equation parameters:
ε 1 = a x i + p × z i ,   ε 2 = b y i + q × z i b 1 = ε 1 ε 1 2 + ε 2 2 ,   b 2 = ε 2 ε 1 2 + ε 2 2 b 3 = ε 1 × z i ε 1 2 + ε 2 2 ,   b 4 = ε 2 × z i ε 1 2 + ε 2 2 b 5 = 1 ;   D i = ( R ε 1 2 + ε 2 2 ) 2 δ a ^ = a a 0 ,   δ b ^ = b b 0 ,   δ p ^ = p p 0 ,   δ q ^ = q q 0
The   initialization   parameters   are   a 0 ,     b 0 ,   p 0 ,   q 0 ,     δ R ^ , and they are determined by the coordinates of the measurement cylinder center position.

3.3.2. Calculation of the Self-Born Weights Function

Equation (24) is divided into t levels according to the number of unknowns, and the first set of equations in each level is selected and combined into a linearly independent combination of error equations. The error equation is expressed as:
V t = B t x ^ l t  
V r = B r x ^ l r
Equation (25) is the set of independent equations, and Equation (26) is the set of conditional equations before iteration. Then, according to Equation (9) to Equation (11), the solution x ^ of the unknown, the number of observed corrections V , and the unit weighted variance valuation σ ^ 2 can be obtained.

3.3.3. Self-Born Weight Iteration and Parameter Calculation

The self-born weights that were calculated were brought into Equation (8) to calculate the number of corrections and unit weight errors of the observations. The new self-born weights are calculated from Equations (18)–(20). The newly calculated self-born weights are brought into Equation (8) to calculate the new number of corrections V and the unit weight errors. The iteration condition is set to correct the number V before and after the difference of the iteration value is less than Δ . The gross error weights far from the column surface are reduced. The axial parameter and cylindricity error are output.

3.3.4. Updating of the Independent Equations

The independent set of equations is Equation (25), and the initial order of the independent set of equations has been obtained from step (2). The number of updates of the independent set of equations is set enough (the number of updates is generally set to the number of samples per level). Step (2) and step (3) are repeated. The set of data with the smallest error is selected. The errors and equation parameters are output. We set enough random sampling times and repeat steps (2) to (3) until a set of data with the smallest error was selected. The number of iterations is determined by the amount of data in each layer. The error and equation parameters are output. The algorithm flow chart is shown in Figure 2.

4. Cylindricity Error Uncertainty Evaluation

4.1. GUM Assessment Method

Currently, the GUM method is the most commonly used method in uncertainty evaluation. According to the GUM method, the uncertainty propagation rate model is determined, and the uncertainty source of each parameter is determined at the same time. The composite standard uncertainty and expanded uncertainty of the cylindricity error uncertainty evaluation are calculated. The sources of uncertainty are analyzed as follows:
(1) Repeat the measurement at a single point 10 times, and the introduced uncertainty component u 1 is calculated according to the Bessel formula [3]:
u 1 = i = 1 n x i x ¯ 2 n 1
(2) The uncertainty component caused by the indicated value error can be derived from the measurement experience of the instrument. The maximum allowable indication error σ of the CMM is 0.5   μ m . The coverage factor k is 3 . The indication error is expressed as u 2 [3]:
u 2 = σ k
(3) During the test, the room temperature is kept constant at 20 degrees Celsius, so the uncertainty u 3 caused by the temperature is approximately 0   μ m [3].
(4) In practical measurements, the uncertainty component u 4 is approximated as 0   μ m due to the deformation caused by the measuring force [3].
Then the single-point uncertainty u 0 is calculated as [3]:
  u 0 = u 1 2 + u 2 2 + u 3 2 + u 4 2
According to the actual situation, the conditions of a single measurement tend to be the same, and the single-point uncertainty is:
u x m a x = u y m a x = u z m a x = u x m i n = u y m i n = u z m i n = u 0
The key to calculating the uncertainty is to find the uncertainty u and ρ of each parameter as the correlation coefficient of each parameter. Uncertainties of a , b , p , q are obtained from actual data. Suppose the vector D = a , b , p , q T , then
c o v D = v a r a c o v a , b c o v a , p c o v a , q c o v b , a v a r b c o v b , p c o v b , q c o v p , a c o v q , a c o v p , b c o v q , b v a r p c o v q , p c o v p , q v a r q
The correlation coefficient can be calculated according to Equation (31) as follows:
ρ a b = u a b u a u b
Finally, the synthetic standard uncertainty u δ is calculated according to the GUM method:
u δ 2 = f x m a x u x m a x 2 + f y m a x u y m a x 2 + f z m a x u z m a x 2 + f x m i n u x m i n 2 + f y m i n u y m i n 2 + f z m i n u z m i n 2 + f p u p 2 + f q u q 2 + f a u a 2 + f b u b 2 + 2 f a f b ρ a b u a b + 2 f a f p ρ a p u a p + 2 f a f q ρ a q u a q + 2 f b f p ρ b p u b p + 2 f b f q ρ b q u b q + 2 f p f q ρ p q u p q

4.2. MCM Assessment Method

The Monte Carlo sampling method evaluation (MCM) can be used as a supplementary verification method to GUM or as a stand-alone method for uncertainty evaluation. The MCM method is mainly divided into three steps: input, propagation, and output. The evaluation steps are shown in Figure 3.

4.3. Validation of GUM Method Results by MCM Method

In order to evaluate the applicability of the GUM experimental results, the results of the GUM method need to be validated using the MCM method. The inclusion interval y ± U p of the inclusion probability p (taken as 95 % in this paper) of the output quantity is calculated using the GUM method. The standard uncertainty u y of the MCM output quantity and the endpoint values y l o w and y h i g h of the shortest inclusion interval are calculated. Finally, the absolute value deviation of the respective endpoints is found by Equations (17) and (18). When the absolute deviations d l o w and d h i g h of the respective endpoints of the two included intervals are less than or equal to the numerical tolerance δ of the uncertainty, the GUM can be verified, indicating that the GUM assessment results are consistent with the actual situation [3].

5. Experimental Verification and Literature Verification

5.1. Experimental Verification

Experimentally, a part of the shaft segment was selected for measurement with a nominal diameter of Φ = 92   mm . The measuring equipment was the German Zeiss CMM PRISMO 9/13/7 ultra-type three-coordinate measuring machine to measure the outer ring of the shaft. According to the ISO 12180-2:2011 standard cylindricity measurement guidelines, various point assignment strategies can be employed, including measurement methods such as roundness (circumference) profiles, bus bars, bird cages, spirals, and random points [22,23]. The circle profile method was selected in the experiment, 5 sections were measured, and 48 points were measured at equal angular intervals on each section. The experiment is shown in Figure 4.

5.1.1. Evaluation of Experimental Data

The least square objective function was constructed from the sampled point data. Then the least square (LS), self-born weighted least square (SWLS), improved self-born weighted least square (ISWLS), standard particle swarm algorithm (standard PSO), and standard genetic algorithm (standard GA) were used to calculate the cylindricity error, respectively. The selection of the irrelevance group of the SWLS method followed the actual situation, and the 1st, 57th, 112th, 167th, and 222nd error equations were selected as the 5 error equations that make up the maximum linear independence group. The number of updates of the independent system of equations was 30 when the experiment uses the ISWLS method. The population size of the genetic algorithm was 200, the number of evolutions was 2000, and the crossover probability and mutation probability were 0.7 and 0.2. The particle swarm optimization population size was 500, the number of evolutions was 1000, and the learning factors were c1 = c2 = 2. The results of the first group of measurement data processing are shown in Table 1. At the same time, the fitting result of the ISWLS method is shown in Figure 5.
In Table 1, the cylindricity axis parameters, cylindricity error value and calculation time were calculated, and it can be seen that the axis positions calculated by the five methods are slightly different. The cylindricity error value of the LS method was 0.1955 mm, which was significantly larger than the error of other algorithms. The standard PSO method and the standard GA method were global optimization algorithms with cylindricity error results of 0.1415 mm and 0.1514 mm, respectively, and the computing time was more than 7 s. The results were more accurate than those of the LS method, but the computation time was longer. The cylindricity error results of the ISWLS method and SWLS method were 0.1346 mm and 0.1336 mm, respectively. Their operation time were 0.0977 s and 0.08295 s, respectively. They both had higher accuracy compared to the LS method and shorter operation time than the standard PSO and standard GA methods. Since the operation time of the ISWLS method did not include the selection time of the independent set of equations, and the influence of gross errors on the selection of the independent set of equations was guaranteed to be excluded in each operation, the algorithm is generally better than the SWLS method.
Under the same conditions, the same sampling method was used to measure the same measurement object 10 times, and the error data were obtained by using the four methods of processing in Section 5.1.1, as shown in Figure 6.
As can be seen from Figure 6, the cylindricity error values calculated by the ISWLS method were consistent with those of the SWLS method, LS method, standard PSO method, and standard GA method in the evaluation of 10 measurement results, which indicated the effectiveness of the ISWLS method. The value of cylindricity error found by the LS method is larger than the value of other algorithms. The standard PSO method and the standard GA method were global optimization algorithms, and the trend of the evaluation error values for the 10 sets of measurement data was the same as that of the LS method, but the evaluation error values were significantly lower. When the SWLS method was used for cylindricity error assessment, the stability of the assessment results was poor because of the different effects of the independent groups in the selection on the assessment results, and some of the cylindricity error values were large compared with the results of the standard PSO method and the standard GA method, which was similar to the LS method. The evaluation error values of each group of the ISWLS method were smaller than those of the SWLS method and LS method. The evaluation results of the ISWLS method were similar to those of the standard PSO method and standard GA method, and the results were stable and have higher accuracy.

5.1.2. Error Evaluation after Inserting Gross Error Points

Considering the influence of gross differences on the error assessment, a numerical experiment of inserting gross differences was designed. According to the existing measurement data, the column surface points of one measurement were selected, and three gross error points were inserted randomly in the 1st, 2nd, and 3rd layer measurement sections to observe the change of the axis position. The insertion of the gross difference data is shown in Table 2. The ISWLS method was used to calculate the axial parameters and cylindricity error values, and the results are shown in Table 3.
In order to reflect the influence of the gross difference value on the error results before and after the iteration of the ISWLS method, the true error value v i of Equation (16) in the ISWLS method was selected as the variable, and the true error values v i , i = 1 , , 240 of 240 points on the column before the iteration were compared with the calculated results after the iteration, as shown in Figure 7.
According to Figure 7a, the true error value v i did not change much before and after the iteration when the gross error point was not inserted, and the result was more stable. As shown in Figure 7b, the true error value v i at the beginning of the iteration was 0.95 mm larger than that at the rest of the points after inserting a gross point, and the deviation of the true error value v i at the gross point was reduced to normal after the end of the iteration. As shown in Figure 7c, after the insertion of 2 gross points, the maximum true error value v i at the 2 gross points inserted at the beginning of the iteration was 9.5 mm compared to the remaining points, and the true error values v i at the gross points can converge to the normal value after the end of the iteration. According to Figure 7d, after the insertion of three gross points, the maximum true error value v i at the three gross points inserted at the beginning of the iteration was 9.8 mm larger than the deviation of the remaining points, and the maximum deviation of the true error value v i at the gross points returned to 0.2 mm after the end of the iteration. It can be seen that the true error values of cylindricity at each point are larger under the influence of inserting gross error points. The true error parameter values of the gross difference part tended to be smooth after iteration. It reflected that the ISWLS method was very effective in reducing the influence of gross errors.
According to Table 3, it can be seen that the radius of the cylindrical fit of the ISWLS method was 46.2953 mm for the insertion of one gross point. The radius of the cylindrical fit was 46.2966 mm for the insertion of two gross points. The radius of the cylindrical fit was 46.2998 mm when three gross points were inserted. Compared to when gross error point was not inserted, the cylindrical radius changed within 0.01 mm, and the algorithm had better resistance to gross error. Cylindricity error was 0.1336 mm when gross error point was not inserted. Cylindricity error was 0.1475 mm when one gross error point was inserted. The cylindricity error was 0.1875 mm when two gross error points were inserted. The cylindricity error was 0.2038 mm when three gross error points were inserted. With the increase of inserting gross error points, the cylindricity error changed within 0.08 mm, and the data change was small. This shows that the algorithm can reduce the influence of gross error on the evaluation results, had a certain resistance to gross error, and had good robustness.

5.1.3. Comparison and Evaluation of Algorithms after Inserting Gross Error Points

The ISWLS algorithm differed from other algorithms in terms of its ability to resist gross differences. The standard PSO, standard GA, SWLS, and ISWLS methods were used to compare the error ratings according to the number of insertion gross error point. The comparison results are shown in Table 4 and Figure 8.
As can be seen from Table 4, when one gross error point was inserted, the cylindricity error of the ISWLS method was 0.1475 mm, the cylindricity error of the SWLS method was 0.1681 mm, and the cylindricity errors of the standard PSO method, standard GA method, and LS method were 1.0698 mm, 1.4444 mm, and 1.1775 mm respectively, which showed that the five algorithms were gradually affected by gross errors. ISWLS had a better ability to resist gross error than SWLS. SWLS had a better ability to resist gross error than other algorithms. When two gross error points were inserted, the cylindricity error of the ISWLS method was 0.1875 mm, and the cylindricity error of the SWLS method was 0.2130 mm. The cylindricity errors of the standard PSO method, standard GA method, and LS method were 10.2447 mm, 9.2046 mm, and 10.6446 mm, respectively. When the number of gross error points increased, the five algorithms were affected by the number of gross error increases. As the number of gross error points increased, ISWLS outperformed SWLS, PSO, GA, and LS in resisting gross errors, and their ability to resist gross errors gradually declined. When three gross error points were inserted, the cylindricity error of the ISWLS method was 0.2038 mm, that of the SWLS method was 0.2911 mm, and that of the standard PSO method, standard GA method, and LS method were 10.8566 mm, 9.4694 mm and 19.0610 mm, respectively. The ISWLS method had less error change than the SWLS method, and its ability to resist gross error was significantly better than that of standard PSO and standard GA methods, while the cylindricity error of the LS method increased significantly with the increase of gross error points, and the ability to resist gross error declined significantly.
As can be seen from Figure 8, with the increase in inserting gross error points, the LS method was most affected by the gross error points, and the cylindricity error value increased gradually. The standard PSO and standard GA methods had a significant decrease in the error elevation of the cylindricity error value after the insertion of the 2nd error point due to the reason of the gross error point location and the advantage of the algorithm. The ISWLS method and the SWLS method had a certain resistance to gross errors. Compared with other algorithms, the cylindricity error assessment results did not increase significantly after inserting gross error points, the changes were relatively flat, and the ISWLS method was more resistant to gross errors than the SWLS method.

5.2. Literature Verification

To verify the method proposed in this paper, data samples of the literature [4] were used in this paper to evaluate the cylindricity error. The 24 data points are shown in Table 5. Reference [4] used the genetic algorithm to evaluate three cylindricity error model methods, including the minimum area method, the minimum circumscribed cylinder method, and the maximum inscribed cylinder method. The error evaluation results based on the MZC, MCC, and MIC methods in the literature [4] are shown in Table 6. The ISWLS, LS, PSO, and GA methods were used to evaluate the error of the least square cylinder method for the sampling point data, and the setting parameters were consistent with the above. The error evaluation results are shown in Table 7.
Comparing Table 6 and Table 7, it can be seen that the method proposed in this paper was consistent with the results obtained in the literature, and the evaluation error was smaller.

5.3. Cylindricity Error Uncertainty Experiment

According to the first group of cylindricity error evaluation data in Table 1, the cylindricity error uncertainty was calculated by the ISWLS method. The combined standard uncertainty of the output was calculated, and the inclusion interval of the 95% inclusion probability was obtained for the interval factor k = 1.96. The data are shown in Table 8.
The MCM method, as a supplementary method to the GUM method, can perform an uncertainty assessment of cylindricity error. The results can be compared with the GUM method. The results are shown in Table 9 and Figure 9.

5.4. Analysis of Results

The ISWLS was compared with the standard PSO, standard GA, and LS methods through experiments and literature validation. It was proved that ISWLS satisfies the definition. A more accurate cylindricity error evaluation result and a shorter running time were obtained. By inserting gross errors, it was verified that the ISWLS method is more resistant to gross errors than other algorithms, with more stable results and better robustness. The uncertainty of the cylindricity error was evaluated by the GUM method and the MCM method. The result of the GUM method was 0.002222 mm and the result of the MCM method was 0.002349 mm. The absolute deviations of d l o w and d h i g h for each endpoint of the two included intervals were 0.0002 and 0.0001, which was less than the uncertainty value tolerance of 0.0005. The uncertainty results of the MCM method can verify the results of the GUM method, which proved the validity of the results and the robustness of the ISWLS method.

6. Conclusions

Since there are cases of gross errors in cylindricity error evaluation, this paper proposes the ISWLS method to evaluate the cylindricity error. The algorithm was compared with the standard PSO, standard GA, and LS methods for error assessment. The ability of the algorithm to resist gross errors was verified by inserting gross error points, and the correctness of the algorithm was verified using literature data. The experimental results show that:
(1) The method in this paper differs from SWLS and conventional global optimization algorithms in that it does not require setting parameters such as independent group order, population size, and crossover variation, and the operation steps are simple. As a result, the algorithm is relatively simple.
(2) By randomly inserting three gross error points, it was demonstrated that the cylindrical axis parameters calculated by the ISWLS method are stable when the data contain gross error points. The variation of cylindricity error is within 0.08 mm. The variation is smaller compared with the results without inserting gross error points, so the algorithm is able to resist gross errors.
(3) In several experimental evaluations, compared with the particle swarm algorithm and genetic algorithm, the method in this paper took relatively less time and had faster computing power than the LS. In the experiment of inserting three gross error points, the LS method was most affected by the gross error points, and the standard PSO method and the standard GA method had less error increase in the cylindricity error value after the insertion of the second error point due to the advantage of the algorithm. The ISWLS method and the SWLS method have a certain resistance to gross errors. Compared with other algorithms, the cylindricity error assessment results did not increase significantly after inserting gross error points, the changes were relatively flat, and the ISWLS method was more resistant to gross errors than the SWLS method.
(4) The LS method, SWLS method, ISWLS, standard PSO method, and standard GA method were used to evaluate the model without inserting gross error points. Compared with other optimization algorithms, ISWLS has higher accuracy and stable arithmetic results with strong robustness. The results of the uncertainty evaluation experiments show that the GUM method results were verified by the MCM method, and thus the error uncertainty evaluation results are reliable. The ISWLS method has the robustness of the cylindrical error evaluation results. The algorithm can be extended to error evaluation.

Author Contributions

Methodology, Y.Y.; Formal analysis, Y.Y.; Resources, K.Z.; Data curation, Y.Y.; Writing—original draft, Y.Y.; Writing—review & editing, K.Z.; Supervision, K.Z.; Project administration, K.Z.; Funding acquisition, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The study was funded by author Y.Y.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No data was made public.

Acknowledgments

The authors would like to thank the Managing Editor and the anonymous reviewers for their pertinent comments on this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zheng, P.; Wu, J.Q.; Zhang, L.N. Research of the on-line evaluating the cylindricity error technology based on the new generation of GPS. Procedia Eng. 2017, 174, 402–409. [Google Scholar] [CrossRef]
  2. Yang, Y.; Li, M.; Wang, C.; Wei, Q. Cylindricity error evaluation based on an improved harmony search algorithm. Sci. Program. 2018, 13, 2483781. [Google Scholar] [CrossRef]
  3. Zhang, K.; Zhang, W.; Yan, W.Z.; Hou, H.S. Research on Evaluation and Uncertainty of Measurement of Circularity Errors via Neural Network Algorithm. Mech. Sci. Technol. Aerosp. Eng. 2019, 38, 428–432. [Google Scholar] [CrossRef]
  4. Wen, X.L.; Song, A. Cylindricity Error Evaluation Based on an Improved Genetic Algorithm. Acta Metrol. Sin. 2004, 25, 115–118. [Google Scholar] [CrossRef]
  5. Ambewadkar, G.N.; Gajre, S.P. Probe Path Planning for Flatness Measurement on Coordinate Measuring Machine Using Ant Colony Optimization. Adv. Eng. Forum 2021, 41, 6362. [Google Scholar] [CrossRef]
  6. Tu, J.; Chen, H.; Liu, J.; Heidari, A.A.; Zhang, X.; Wang, M.; Pham, Q.V. Evolutionary biogeography-based whale optimization methods with communication structure: Towards measuring the balance. Knowl. Based Syst. 2021, 212, 106642. [Google Scholar] [CrossRef]
  7. Zhang, Y.; Wu, M.; Fu, X.; Li, Q.; Gao, J. Reduced algorithm of laser measurement data for geometric parameters of ring forgings based on Artificial Immune Algorithm. Measurement 2019, 134, 527–538. [Google Scholar] [CrossRef]
  8. Zhang, K.; Cheng, G.; Zhang, W. Evaluation of roundness error uncertainty by a Bayesian dynamic model. Measurement 2020, 155, 107565. [Google Scholar] [CrossRef]
  9. He, G.Y.; Liu, P.P.; Wang, K. Evaluation of the cylindricity error based on the sequential quadratic programming algorithm. Mech. Sci. Technol. Aerosp. Eng. 2014, 33, 1845–1849. [Google Scholar] [CrossRef]
  10. Wu, Q.; Zhang, C.; Zhang, M.; Yang, F.; Gao, L. A modified comprehensive learning particle swarm optimizer and its application in cylindricity error evaluation problem. Math. Biosci. Eng. 2019, 16, 1190–1209. [Google Scholar] [CrossRef]
  11. Liu, D.; Zheng, P.; Wu, J.; Yin, H.; Zhang, L. A new method for cylindricity error evaluation based on increment-simplex algorithm. Sci. Prog. 2020, 103, 0036850420959878. [Google Scholar] [CrossRef]
  12. Li, Q.; Ning, H.; Gong, J.; Li, X.; Dai, B. A hybrid greedy sine cosine algorithm with differential evolution for global optimization and cylindricity error evaluation. Appl. Artif. Intell. 2021, 35, 171–191. [Google Scholar] [CrossRef]
  13. Wang, Y.X.; Wang, H.X.; Tian, H.H. Cylindricity Error Measurement and Evaluation Based on Step Acceleration Algorithm in Crankshaft Measuring Machine. MAPAN 2022, 2022, 1–10. [Google Scholar] [CrossRef]
  14. Khattab, A.A.; Abouelseoud, Y.; Younes, M.A. Effect of Sampling Rates Variation during Cylindricity Error Evaluation. Int. J. Eng. Appl. Sci. 2019, 3, 257585. [Google Scholar]
  15. Zaky, K.M.; Ghonem, A.A. Application of Robust Estimation Methods for Detecting and Removing Gross Errors from Close-Range Photogrammetric Data. Am. Acad. Sci. Res. J. Eng. Technol. Sci. 2019, 55, 111–120. [Google Scholar]
  16. Yan, G.; Cen, M.; Li, Y. Gross error detectability and identifiability analysis in track control network for high-speed railway based on GEJE. J. Surv. Eng. 2020, 146, 04019020. [Google Scholar] [CrossRef]
  17. Jia, N.N.; Ge, Y.H. Remainder reliability and robust estimation: A case study using twelve simulated leveling networks. Appl. Mech. Mater. 2013, 263, 3163–3167. [Google Scholar] [CrossRef]
  18. Ge, Y.H. Research on the method of self-born weighted least square. Bull. Surv. Mapp. 2014, 8, 36–39. [Google Scholar] [CrossRef]
  19. Survey Adjustment Subject Group, School of Surveying and Mapping, Wuhan University. Error Theory and Basis of Survey Adjustment; Wuhan University Press: Hubei, China, 2014. (In Chinese) [Google Scholar]
  20. Oliveira, E.C.; Aguiar, P.F. Data Reconciliation in the Natural Gas Industry: Analytical Applications. Energy Fuels 2009, 23, 3658–3664. [Google Scholar] [CrossRef]
  21. Ozyurt, D.B.; Pike, R.W. Theory and practice of simultaneous data reconciliation and gross error detection for chemical processes. Comput. Chem. Eng. 2004, 28, 381e402. [Google Scholar] [CrossRef]
  22. Mian, S.H.; Umer, U.; Abdulhameed, O.; Alkhalefah, H. Prognosis and Multiobjective Optimization of the Sampling Plan for Cylindricity Evaluation. Math. Probl. Eng. 2020, 2020, 8201312. [Google Scholar] [CrossRef]
  23. ISO 12180-2:2011; Geometrical Product Specifications (GPS)—Cylindricity—Part 2: Specification Operators. ISO: Geneva, Switzerland, 2021.
Figure 1. Diagram of cylindricity error.
Figure 1. Diagram of cylindricity error.
Applsci 12 12319 g001
Figure 2. ISWLS cylindricity error algorithm evaluation process.
Figure 2. ISWLS cylindricity error algorithm evaluation process.
Applsci 12 12319 g002
Figure 3. MCM process for assessing measurement uncertainty.
Figure 3. MCM process for assessing measurement uncertainty.
Applsci 12 12319 g003
Figure 4. Cylindricity error measurement.
Figure 4. Cylindricity error measurement.
Applsci 12 12319 g004
Figure 5. ISWLS method cylindrical fitting.
Figure 5. ISWLS method cylindrical fitting.
Applsci 12 12319 g005
Figure 6. Cylindricity error evaluation comparison.
Figure 6. Cylindricity error evaluation comparison.
Applsci 12 12319 g006
Figure 7. Statistics chart of true error value when inserting gross error points. (a) No gross error point insertion. (b) Insert one gross error point. (c) Insert two gross error points. (d) Insert three gross error points.
Figure 7. Statistics chart of true error value when inserting gross error points. (a) No gross error point insertion. (b) Insert one gross error point. (c) Insert two gross error points. (d) Insert three gross error points.
Applsci 12 12319 g007
Figure 8. Influence of gross error on the results of different algorithms.
Figure 8. Influence of gross error on the results of different algorithms.
Applsci 12 12319 g008
Figure 9. Statistics chart of true error value when inserting gross error point.
Figure 9. Statistics chart of true error value when inserting gross error point.
Applsci 12 12319 g009
Table 1. The first set of data evaluates the results.
Table 1. The first set of data evaluates the results.
Evaluation MethodAxis Parameters
a   m m b   m m p   m m q   m m
LS0.00040.0000−0.0009−0.0007
SWLS−0.0087−0.00100.00040.0000
ISWLS0.02670.0275−0.0010−0.0010
PSO0.00030.00100.01420.0943
GA−0.0013−0.00080.02560.0141
Evaluation MethodCylindricity Error m m Running Time s
LS0.19550.074023
SWLS0.13460.177704
ISWLS0.13360.082950
PSO0.14157.467153
GA0.151415.34568
Table 2. Insertion point data.
Table 2. Insertion point data.
No. X   m m Y   m m Z   m m
146.2840310.0000822
234.6397812.2283619
3−55.363149.2202916
Table 3. Evaluation results of ISWLS method when inserting gross error points.
Table 3. Evaluation results of ISWLS method when inserting gross error points.
Insert Gross PointAxis Parameters
a   m m b   m m p   m m q   m m
0−0.0087−0.00100.00040.0000
1−0.0467−0.00920.00330.0007
20.07860.0356−0.0090−0.0038
3−0.00270.0488−0.0090−0.0038
Insert Gross PointCylinder Radius r m m Cylindricity Error s
046.29090.1336
146.29530.1475
246.29660.1875
346.29980.2038
Table 4. Error evaluation results of different algorithms when inserting gross error points (mm).
Table 4. Error evaluation results of different algorithms when inserting gross error points (mm).
Evaluation Method1 Gross Error Point [mm]2 Gross Error Points [mm]3 Gross Error Points [mm]
LS1.177510.644619.0610
SWLS0.16810.21300.2911
ISWLS0.14750.18750.2038
PSO1.069810.244710.8566
GA1.44449.20469.4694
Table 5. Coordinates of points in reference [4]. Reprinted with permission from Ref. [4]. 2004, Wen, X.L.; Song, A.
Table 5. Coordinates of points in reference [4]. Reprinted with permission from Ref. [4]. 2004, Wen, X.L.; Song, A.
No. X   m m Y   m m Z   m m No. X   m m Y   m m Z   m m
111.09430.452265.23281310.8150.591885.2307
25.09410.84565.0765144.814810.984685.074
3−6.906310.843965.008915−7.185510.983585.0641
4−12.90650.449865.089716−13.1850.589484.8952
5−6.9063−9.942965.05417−7.1855−9.803385.0516
65.094−9.941865.2216184.8149−9.802285.2171
710.95460.52275.23161910.67540.661695.2291
84.954410.914875.0752204.675211.054495.0728
9−7.045910.913775.07721−7.325311.053394.9077
10−13.04610.519674.896422−13.32540.659295.094
11−7.0459−9.873175.052823−7.3252−9.733595.0504
124.95447−9.87275.2204244.6752−9.732395.2179
Table 6. The data evaluation results in the literature [4]. Reprinted with permission from Ref. [4]. 2004, Wen, X.L.; Song, A.
Table 6. The data evaluation results in the literature [4]. Reprinted with permission from Ref. [4]. 2004, Wen, X.L.; Song, A.
Evaluation MethodAxis ParametersCylindricity Error m m
a   m m b   m m p   m m q   m m
MZC0.0024−0.0040−0.01400.00700.00279
MCC0.0025−0.0042−0.01400.00700.00282
MIC0.0037−0.0039−0.01400.00700.00279
Table 7. ISWLS method data evaluation results.
Table 7. ISWLS method data evaluation results.
Evaluation MethodAxis Parameters
a   m m b   m m p   m m q   m m
LS−0.01400.00700.0036−0.0046
PSO0.0024−0.0250−0.01390.0069
GA0.0045−0.0068−0.01400.0070
ISWLS0.0037−0.0034−0.01400.0070
Evaluation MethodCylindricity Error m m Running Time s
LS0.002880.034277
PSO0.002866.508752
GA0.0027912.568197
ISWLS0.0027880.041420
Table 8. GUM method evaluation results of the first set of data.
Table 8. GUM method evaluation results of the first set of data.
ParameterValue μ m NameValue μ m
u 1 0.0052 u a 0.3274
u 2 1.556 u b 0.3274
u 3 0 u p 0.0012
u 4 0 u q 0.0012
u 0 1.560 u δ 2.2216
ParameterValue m m
best estimate0.1338
95% probability inclusion interval[0.1294, 0.1381]
Table 9. MCM method evaluation results of the first set of data.
Table 9. MCM method evaluation results of the first set of data.
Parameter Value   m m
best estimate0.1338
95% probability inclusion interval[0.1294, 0.1381]
Project d l o w d h i g h MCM verifies GUM results (δ ≤ 0.0005)
MCM Verification GUM0.00020.0001Yes
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yao, Y.; Zhang, K. An Improved Self-Born Weighted Least Square Method for Cylindricity Error Evaluation. Appl. Sci. 2022, 12, 12319. https://doi.org/10.3390/app122312319

AMA Style

Yao Y, Zhang K. An Improved Self-Born Weighted Least Square Method for Cylindricity Error Evaluation. Applied Sciences. 2022; 12(23):12319. https://doi.org/10.3390/app122312319

Chicago/Turabian Style

Yao, Yunhan, and Ke Zhang. 2022. "An Improved Self-Born Weighted Least Square Method for Cylindricity Error Evaluation" Applied Sciences 12, no. 23: 12319. https://doi.org/10.3390/app122312319

APA Style

Yao, Y., & Zhang, K. (2022). An Improved Self-Born Weighted Least Square Method for Cylindricity Error Evaluation. Applied Sciences, 12(23), 12319. https://doi.org/10.3390/app122312319

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