Scientific Programming
Volume 2018 (2018), Article ID 2483781, 13 pages
https://doi.org/10.1155/2018/2483781

Research Article

Cylindricity Error Evaluation Based on an Improved Harmony Search Algorithm

Yang Yang,¹ Ming Li,¹ Chen Wang,² and QingYue Wei¹

¹Key Laboratory of Intelligent Manufacturing and Robotics, School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200072, China
²Department of Mechanical Engineering, Hubei University of Automotive Technology, Shiyan 442002, China

Correspondence should be addressed to Yang Yang; mryyfi@i.shu.edu.cn

Received 30 March 2018; Revised 7 June 2018; Accepted 14 June 2018; Published 19 July 2018

Academic Editor: Ricardo Soto

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

Abstract.

The cylindricity error is one of the basic form errors in mechanical parts, which greatly influences the assembly accuracy and service life of relevant parts. For the minimum zone method (MZM) in international standards, there is no specific formula to calculate the cylindricity error. Therefore, the evaluation methods of the cylindricity error under the MZM have been widely concerned by international scholars. To improve the evaluation accuracy and accelerate the iteration speed of the cylindricity, an improved harmony search (IHS) algorithm is proposed and applied to compute the cylindricity. On the basis of the standard harmony search algorithm, the logistic chaotic initialization is introduced into the generation of initial solution to improve the quality of solutions. During the iterative process, the global and local search capabilities are balanced by adopting the par and operators adaptively. After each iteration, the Cauchy mutation strategy is adopted to the best solution to further improve the calculation precision of the IHS algorithm. Finally, four test functions and three groups of cylindricity error examples were applied to validity verification of the IHS algorithm, the simulation test results show that the IHS algorithm has advantages of the computing accuracy and iteration speed compared with other traditional algorithms, and it is very effective for the application in the evaluation of the cylindricity error.

1. Introduction

In order to evaluate the manufacturing quality of mechanical parts, the inspection process has always been an important step in the whole life cycle of mechanical products. As an important parameter to evaluate the manufacturing accuracy of the internal and external cylindrical surfaces of the rotary part, the cylindricity error has become the focus of the research in the field of mechanical measurement. Coordinate measuring machine (CMM) is a common tool to inspect the cylindricity error of the mechanical parts, which is widely used because of its high precision. The spatial coordinates information of the measured point of the part is obtained by the sensor of the CMM. Then, the measurement data obtained from the sensor is analyzed through the error evaluation algorithm to calculate the cylindricity error of the measured parts. At present, in order to obtain the cylindricity error of measurement data more quickly and accurately, a lot of research results have been obtained. It is of great value in engineering application to evaluate the cylindricity error quickly and accurately, which has great influence on the function and life of the product.

For the cylindricity error evaluation, the least square method (LSM) is widely used in the field of engineering applications, but the accuracy of LSM is relatively low and cannot be applied in the high-precision evaluation field [1, 2]. Compared with the LSM, the minimum zone method (MZM) is one of the techniques that has a higher evaluation accuracy. It can meet the minimum zone principle in the international standards and obtain the accurate results of the measurement data. However, there is no specific algorithm in international standards, so it has received widely attention and research on the MZM by international scholars.

Based on the initial research results, the nonlinear optimization algorithms [3, 4] are the most important methods for cylindricity error evaluation under the MZM. Moreover, the nonlinear optimization algorithms are relatively complex and difficult in the construction of the mathematical model. However, with the continuous improvement of the artificial intelligence technology, the intelligent optimization algorithms have developed an effective method to solve such complex nonlinear problems. And some of them have been introduced into the cylindricity error evaluation because of their simplicity and effectiveness, such as GA [5] and PSO [6, 7]. In addition, compared with nonlinear optimization algorithms, the construction of the mathematical model for intelligent optimization algorithm is also very simple. But the initialization parameters of these algorithms have great influence on the calculation results. Therefore, the improvement of the accuracy and convergence speed of the cylindricity algorithm has become an important research topic at this stage.

Harmony search (HS) algorithm is an intelligent optimization algorithm, which simulates the creation process of the music [8]. It has the advantages of relatively simple computational principle and strong global search ability, so it is very suitable for engineering applications such as the sizing optimization of truss structures and reliability problems [9, 10]. But for the basic HS algorithm, there are some problems such as the premature convergence in the iteration process and low accuracy in later iteration.

To improve the evaluation accuracy and accelerate the iteration speed of the cylindricity, on the basis of the standard HS algorithm, the logistic chaotic initialization is introduced into the generation of initial solution to improve the quality of solutions. During the iterative process, the global and local search capabilities are balanced by adopting the par and operators adaptively. After each iteration, the Cauchy mutation strategy is adopted to the best solution to further improve the calculation precision of the IHS algorithm. Finally, four test functions and three groups of cylindricity error examples were applied to validity verification of the IHS algorithm, the simulation test results show that the IHS algorithm has advantages of the computing accuracy and iteration speed compared with other traditional algorithms, and it is very effective for the application in the evaluation of the cylindricity error.

The organizational structure of this paper is shown as follows: the research background and significance of the paper are introduced in Section 1. In Section 2, the present research results of this work are described. In Section 3, the mathematical model of the cylindricity error and the objective function of the problem are established. The basic HS algorithm and IHS algorithm are proposed in Sections 4 and 5, respectively. The simulation experiments of the algorithms and cylindricity error analysis are carried out in Section 6. In the last section of this paper, the conclusion of this work is summarized and the future work is prospected.

2. Related Works

In the research results of the cylindricity error evaluation under the MZM, Carr and Ferreira established the nonlinear mathematical model of the cylindricity and straightness errors and then used the linear programming method (LPM) to solve the problem, but this method requires a set of appropriate initial solutions [3]. Lai and Chen converted the spatial cylindricity into the plane problems by nonlinear transformation and used the parameters adjustment to obtain the cylindricity error [4]. Chou and Sun deduced a general mathematical model of the cylindricity error. The simulated annealing (SA) algorithm is applied to solve the model that shows the advantage in the global optimization performance [11]. Lai et al. used the genetic algorithm (GA) to solve the mathematical model of the cylindricity error under the MZM, and the experiment proves that the method has good flexibility and accuraxcy [5]. Weber et al. applied the linear approximation technique to the form error evaluation that includes the cylindricity error. The method provided an approximate solution for the form error [12]. Zhu and Ding combined the motion geometry and sequence approximation algorithm to analyze and evaluate the cylindricity error. The experiments proved the validity of the algorithm [13]. Lao et al. improved the hyperboloid technique by constructing an initial axis to establish the cylindrical error evaluation model, and the experiment showed that the method further improved the accuracy of cylindricity evaluation [14]. Cui et al. and Mao et al. applied a particle swarm optimization (PSO) algorithm to cylindricity error evaluation, and the calculation accuracy was further improved, respectively [6, 7]. Venkaiah and Shunmugam designed a cylindricity error evaluation algorithm based on the computational geometry method [15]. Bei et al. and Guo et al. improved the genetic algorithm and applied it to the cylindricity error assessment, respectively [16, 17]. Li et al. used the coordinate transformation to construct the cylindricity error mathematical model, and the model has a high evaluation accuracy [18]. Luo et al. further improved the global optimization ability of an artificial bee colony algorithm by introducing tabu strategy, thus increasing the evaluation speed of the cylindricity error [19]. Wen et al. modified the particle velocity using a contraction factor and improved the convergence ability of the PSO. The improved PSO (IPSO) was applied to the cylindricity and conicity errors which have a good flexibility [20]. Lei et al. designed a method of cylindricity error evaluation based on the geometry optimization searching algorithm. The method can obtain the cylindricity error by setting the hexagon in the measurement points [21]. Lee et al. developed a cylindricity error evaluation based on support vector machines (SVM), and the algorithm can be applied in the field of the machine vision system [22]. Wen et al. utilized a quasiparticle swarm algorithm and calculated the cylindricity error [23]. He et al. used kinematic geometry and a sequential quadratic programming algorithm to solve the cylindricity error of measured parts. The method has high stability and accuracy [24]. Hermann introduced the application of the computation geometry technique to the form errors that include cylindricity [25]. Zheng et al. adopted a linear programming method to construct the mathematical model of the cylindricity error. The model was solved by a new simplex method [26].

On the basis of the above literatures, we improved the HS algorithm through the generation of initial solution, the dynamic adjustment of operators, and the optimal solution perturbation, to further enhance the optimization ability and the iterative speed of the basic HS algorithm. It is designed to provide a reliable method for the evaluation of the cylindricity error.

3. Mathematical Model of Cylindricity Error Evaluation

According to the concepts in the relevant standard [1, 2], the cylindricity error under the MZM can be defined as follows: the area is the radius difference of two coaxial ideal cylinders that contain the measured cylinder, and when the area reaches the minimum value, the radius difference is the cylindricity error under the MZM. As shown in Figure 1, the graph shows the measurement process of the cylindricity error. For the cylindrical parts, due to the manufacturing errors of the machine tool, the actual parts will always have a certain error compared with the ideal parts, and the cylindrical surface will produce a certain deformation. Therefore, the CMM is used to obtain the coordinates information of the measured parts. Figure 2 illustrates the principle of the cylindricity error calculation, where is the axis of the two ideal cylinders and is the radius difference; when takes the minimum value, the value of is the cylindricity error of the cylinder measured under the MZM.

Figure 1: The measurement process of the cylindricity error.

Figure 2: Schematic of the cylindricity error.

As a space element, the mathematical expression of the cylindricity error requires more parameters. In order to reduce the amount of computation, the relevant literatures simplified the model by assuming conditions, which does not really match the reality [5, 17]. In order to ensure the generality of cylindricity error evaluation, the general equation of the cylinder error is established under the MZM to obtain the error. It is assumed that the points of the measured cylinder are , where is the coordinates of the measured point in the Cartesian coordinate system and is the number of measured points. is the coaxial cylindrical axis, and the point-oriented parametric equation is formula (1) [16], where , , and is the fixed point of and the , , and are the directional parameters, respectively. For the measured point , its distance equation to the space axis is shown as formula (2), where is the distance from the point to the axis and , , and are the unit vectors of the direction:

Therefore, for the intelligent optimization algorithm, the objective function of the problem needs to be established. Based on the minimum condition of the cylindricity error, it can be converted to a nonlinear minimum problem. The objective function is shown as formula (3) [16]. The target parameters are the value of , , , , , and when met the minimum value:

Formula (3) describes the mathematical model of the cylindricity error. is the value of the cylindricity error of the measured part. , , , , , and are the parameters value of the two ideal cylindrical axes when obtains the minimum value that can satisfy the MZM condition. Therefore, for this nonlinear minimization problem, it can be solved by the HS algorithm.

4. Basic Harmony Search Algorithm

The HS algorithm is a heuristic intelligent optimization algorithm proposed by Geem et al. [8] in 2001. It simulates the creation process of music to achieve the goal of optimization. The audience's evaluation of the music is the function fitness value. HS algorithm flow includes the parameters initialization, the generation of new solutions, and the updating of the harmony memory library (HMS). The HS algorithm process is shown as follows.

4.1. Parameter Initialization

In the HS algorithm, is the size of the HM, is the harmony vector of the HMS (the solution of problem), is the component of harmony vector , is the dimension of the problem, and are the lower and upper values in each dimension, the specific expression is , is the harmony memory considering the rate, par is the pitch adjusting rate, and is the bandwidth. According to the abovementioned variable meaning, the harmony vector i’s achievement is , is the function fitness value, and the optimal solution is the , and initialize the harmony memory library by

where is a random number whose value is 0 to 1.

4.2. Generation of New Harmony Vector

In this stage, a new harmonic vector is generated by three methods, which are harmony memory, tone adjustment, and random generation. The new harmony vector can be expressed as , and each component of is produced by the above three methods. Firstly, the probability of HMCR is used to control the new vector component, and the pseudocode is shown as Algorithm 1; secondly, after considering the probability of memory consideration, it needs to consider probability of the tone adjustment, and the pseudocode is shown as Algorithm 2.

If rand ≤ P then

Else

End If

Algorithm 1: Pseudocode of the first stage.

If rand < par then

End If

Algorithm 2: Pseudocode of the second stage.

4.3. Update the Harmony Memory Library

Update the HMS by comparing the new solution with the worst solution in HM. If is better than , then accept ; if is worse than , then keep .

4.4. Termination Condition

By setting the number of iterations or calculating precision, the algorithm is judged to be finished. The pseudocode of the HS algorithm is shown as Algorithm 3.

Start

Set N, HMS, D, P, par,

Generate HM harmony vector by formula (4)

Calculate

t = 0

While t ≤ T do

    For j = 1 to D do

        If rand ≤ P then

            If rand ≤ par then

            End if

        Else

        End if

    End for

    If then

    End if

    t = t + 1

End while

Output

Algorithm 3: Pseudocode of the HS algorithm.

5. Improved Harmony Search Algorithm

5.1. The Motivation of the IHS

HS algorithm is an intelligent optimization algorithm derived from the creation of the music. For the basic HS, the best solution in the HM is an individual with the best function value. Although the HS algorithm has a good optimization ability, it is still easy to fall into the local optimum. According to the principle of the basic HS algorithm, there are three main problems. Firstly, the distribution of the initial solution will greatly affect the initial computational performance of the algorithm. Secondly, the par and of the HS algorithm greatly affect the algorithm optimization efficiency. Thirdly, the HS algorithm only updates the worst solution and does not search the optimal solution. Therefore, based on the three problems existing above in the HS algorithm and combined with the problems to be solved, the initial solution generation, dynamic adjustment of the operator, and optimal solution disturbance are proposed to increase the diversity of the population and prevent the premature convergence of the algorithm.

5.2. Chaos Initialization

In the initialization process of the intelligent optimization algorithm, different generation methods will produce different initial solution sequences. The basic HS algorithm generates the initial population in a random way that has some problems such as the uneven distribution of the original individual and thus affects the optimization performance of the algorithm. To make the initial solution more uniformly distributed and improve the quality of the initial solution, the chaotic initialization is applied to the standard HS. Chaos sequence has the advantages of traversal, randomness, and regularity, and it can be obtained by a definite equation [27]. The sequence of initial solution generated by logistic chaotic [28] mapping is proposed as

where is 4, and it indicates that the system is a completely chaotic; and are the upper and lower limits of the variable in solution space, respectively; is a initialization variable generated by chaotic in (0,1]; and is the variable that is mapped to the range of value between and by .

For i = 1 to N

End

End

Algorithm 4: Pseudocode of population initialization.

5.3. New Parameters of par and

In the basic harmony search algorithm, par and are the two key control parameters. The global and local search ability of the algorithm is mainly controlled by par, and the convergence speed of the algorithm is adjusted by . However, par and are constant values, and they are impossible to adjust the global and local search ability effectively in the iterative process of the algorithm, which leads to the algorithm falling into the local optimum.

In order to improve the optimization ability and iterative speed of the basic HS, the par and are adjusted by a dynamic factor. For par, a smaller value is taken at the beginning of the iteration to increase the global exploration ability, and the local search ability is enhanced by a larger value at the end of the algorithm. For , a large value is taken at the beginning of the iteration to increase the diversity of the population, and a smaller value is taken in the late iteration inspired by a simplified artificial fish warm algorithm (SAFSA) [29]. The new par and are shown as

where is the initial value of 0.5 and is the current number of iterations.

5.4. Cauchy Mutation Strategy

In the iterative process of the harmonic search algorithm, the optimal solution usually swings between the historical optimal solutions, so it needs to further increase the search space of the optimal solution and avoid the algorithm falling into the local optimum. In order to improve the accuracy of the algorithm, the Cauchy mutation strategy is introduced into the algorithm; in each iterative process, the global optimal solution is mutated, and the mutation prevents the optimal solution oscillation and accelerates the convergence speed of the algorithm, β is defined as the random number of Cauchy (0, 1) distribution, as shown in the following formulae:

5.5. Algorithm Steps of IHS

Based on the corresponding improvement methods, the proposed IHS algorithm process is shown below:

(1)Parameters initialization: N is the size of the HM; D is the dimension of the problem; and are the upper and lower limits of the variables; P is the harmony memory considering the rate; par is the pitch adjusting rate; and is the bandwidth. Then, go to step (2).(2)Chaos initialization: As described in Section 5.2, random sequence individuals which number is N are generated, and the fitness function value of each individual is calculated; the generated N harmony vectors are deposited into the harmony memory as the initial solution vectors. Then, go to step (3).(3)Generate a new harmony vector: As described in Section 4.2, a new harmony is generated by harmony memory, tone adjustment, and random generation, respectively, where par and are replaced by Section 5.2. Then, go to step (4).(4)Update the HMS, as described in Section 4.3. Then, go to step (5).(5)Cauchy mutation strategy: According to Section 5.3, the best solution is updated, shown as formulae (7) and (8). Then, go to step (6).(6)Iteration stopping criterion: If the calculation satisfies the stopping criterion, the calculation is completed and the optimal solution is output; if not, return to step (3). The pseudocode of IHS is shown as Algorithm 5.

Start

Set N, HMS, D, P, par,

For i = 1 to N

End

End

t = 0

While t ≤ T

    For i = 1 to D

        If rand ≤ P

            If

            End

        Else

        End

    End

    If

    End

Update as formulae (7) and (8)

End while

Output

Algorithm 5: Pseudocode of IHS.

6. Experiment and Discussion

In order to verify the validity of the IHS algorithm for the calculation of the cylindricity error, this paper has selected four test functions and three groups of examples to test the performance of the IHS algorithm. The computer system in this experiment is Win 7 64 bit, the processor is i5-4200ucpu @ 1.60 GHz, the installation memory is 4.00 GB, each group experiment is calculated 30 times, experiment iterations of each group is 500, the population quantity is 50, the dimension of problem D is 6, is 0.5, P is 0.9, and the upper and lower bounds of variables are [−30, 30].

6.1. Performance Test

To test the performance of the IHS, four standard test functions are adopted. As shown in Table 1, sphere function is a simple one-peak function that mainly tests the iteration speed of the algorithm. Rosenbrock function is a relatively complex single-peak function to test the optimal performance. Rastrigin function and Griewank function are both multimodal test functions. But the Griewank function is relatively complex than the Rastrigin function, and they are commonly used to test the global search capability of the algorithm. The images of the above four test functions are shown as Figures 3–6.

Table 1: Four test functions.

Number
Name
Function
Range
Min.

f₁
Sphere

[−30, 30]
0

f₂
Rosenbrock

[−30, 30]
0

f₃
Rastrigin

[−5.12, 5.12]
0

f₄
Griewank

[−10, 10]
0

Figure 3: Sphere function.

Figure 4: Rosenbrock function.

Figure 5: Rastrigin function.

Figure 6: Griewank function.

To test the computing performance of the IHS, the HS, firefly algorithm (FA), dragonfly algorithm (DA), artificial bee colony (ABC), and modified harmony search (MHS) algorithm are applied to the test as contrast algorithms. The parameters of the IHS are shown in the following: the number of initial population is 50, the dimension of the problem is 6, the number of global iteration is 500, and the algorithm is tested 30 times.

The experimental results are shown in Table 2. For the test results of the test functions, compared with the HS, FA, DA, ABC, and MHS algorithms, the IHS algorithm has the lowest mean value that represents highest computing precision, and the lowest standard deviation that represents the high computational stability in sphere and Rosenbrock functions. For the test results of Rastrigin and Griewank, IHS can converge to the global optimal solution, while other algorithms fall into the local optimal. Therefore, IHS algorithm has good optimization ability compared with other algorithms.

Table 2: Results of test functions ().

Algorithm
Function
Minimum value
Maximum value
Mean value
Standard deviation

FA
f₁
2.40E − 06
7.46E − 06
4.79E − 06
3.41E − 11

DA
1.04E − 14
2.39E − 07
4.28E − 08
5.67E − 14

ABC
3.39E − 12
3.93E − 11
1.87E − 11
2.25E − 21

HS
1.85E − 06
3.02E − 05
1.17E − 05
7.32E − 10

MHS
1.17E − 14
4.94E − 10
5.57E − 11
2.15E − 19

IHS
1.40E − 17
4.27E − 17
2.55E − 17
1.83E − 26

FA
f₂
1.11E − 02
2.81E − 02
1.96E − 02
2.00E − 04

DA
8.36E + 00
1. 22E + 01
9.37E + 00
1.18E + 01

ABC
1.27E + 00
3.14E + 00
2.29E + 00
2.52E + 00

HS
9.56E − 01
3.80E + 00
1.48E + 00
1.13E + 00

MHS
3.45E + 00
5.51E + 00
4.65E + 00
3.89E + 00

IHS
3.17E − 03
1.71E − 02
8.04E − 03
1.75E − 05

FA
f₃
7.79E − 04
1.93E − 03
1.22E − 03
1.26E − 06

DA
9.95E + 00
1.62E + 01
1.83E + 01
3.86E + 01

ABC
3.50E + 00
7.27E + 00
6.00E + 00
9.69E + 00

HS
9.38E − 01
3.31E + 00
1.95E + 00
5.55E + 00

MHS
2.31E − 04
1.98E − 02
6.73E − 03
4.54E − 04

IHS
0
0
0
0

FA
f₄
2.46E − 02
2.10E − 01
1.04E − 01
3.27E − 02

DA
2.42E − 01
8.04E − 01
4.23E − 01
2.49E − 01

ABC
1.41E − 01
2.00E − 01
1.70E − 01
4.44E − 03

HS
1.33E − 02
1.88E − 01
4.44E − 01
9.11E − 05

MHS
1.17E − 08
2.37E − 06
8.59E − 07
6.17E − 12

IHS
0
0
0
0

Figures 7, 9, 11, and 13 describe the iterative convergence curves of FA, DA, ABC, HS, MHS, and IHS. Meanwhile, the optimal solution boxplot of the algorithms under each test function is also drawn. It can see that the IHS algorithm converges faster and converges to the optimal solution on Rastrigin and Griewank test functions. And as it also can be seen from Figures 8, 10, 12, and 14 that the IHS algorithm has a concentrated distribution of optimization results in 30 times, which shows has good stability. Based on the results of the test function of the six algorithms, it can been seen that the IHS algorithm is better than the other five algorithms from the comparison of the experimental curves in Figures 7–14.

Figure 7: Curve of sphere function.

Figure 8: Boxplot of sphere function.

Figure 9: Curve of Rosenbrock function.

Figure 10: Boxplot of Rosenbrock function.

Figure 11: Curve of Rastrigin function.

Figure 12: Boxplot of Rastrigin function.

Figure 13: Curve of Griewank function.

Figure 14: Boxplot of Griewank function.

Table 3 shows the statistical analysis results of the Wilcoxon-rank test of the algorithms. The value is the Wilcoxon test value of IHS algorithm with other five algorithms at confidence level of 0.05. According to the results in Table 3, the calculated values are all less than 0.05, and there are significant differences between IHS and the other five algorithms under these four test functions. Therefore, the conclusion is of statistical significance.

Table 3: value of the Wilcoxon-rank test.

Function
FA
DA
ABC
HS
MHS

f₁
1.82E − 04
1.82E − 04
1.82E − 04
1.82E − 04
1.82E − 04

f₂
5.85E − 04
5.85E − 04
5.85E − 04
5.85E − 04
5.85E − 04

f₃
6.39E − 05
6.39E − 05
6.39E − 05
6.39E − 05
6.39E − 05

f₄
6.38E − 05
6.38E − 05
6.38E − 05
6.38E − 05
6.38E − 05

6.2. Cylindricity Error Experiment 1

To verify the effectiveness of the IHS algorithm in the cylindricity error evaluation, the data of the measured cylindrical parts were obtained by the CMM and are shown in Table 4. HS and IHS algorithms are used to calculate the cylindricity error.

Table 4: Measurement data (mm).

1
12.6722
17.869
4.9991

2
21.9325
14.2729
5.01

3
25.8783
5.1659
5

4
22.2705
−3.9311
5

5
12.6857
−7.8611
5.0002

6
3.5711
−3.7504
5.0002

7
0.135
4.8405
5

8
4.0705
14.2751
4.9999

9
12.6777
17.8729
9.9995

10
21.9385
14.2734
10.0006

11
25.87
5.1675
10

12
22.2712
−3.9328
10

13
12.689
−7.866
9.9966

14
3.5762
−3.7559
9.9999

15
0.1314
4.8369
10.0002

16
4.0694
14.2654
10

17
12.6799
17.8697
15

18
21.9398
14.2659
15.0009

19
25.8742
5.173
15.05

20
22.2676
−3.9255
15

21
12.6891
−7.8421
14.9996

22
3.588
−3.7526
15.0005

23
0.1364
4.8404
14.9994

24
4.0696
14.2729
15

25
12.6832
17.8678
20.0001

26
21.9382
14.268
19.9995

27
25.8706
5.1725
20

28
22.2654
−3.9283
20.0001

29
12.6813
−7.8449
20.0008

30
3.5785
−3.7515
20

31
0.1296
4.8316
19.9993

32
4.073
14.2676
20.0008

The calculation results of HS and IHS algorithms are listed in Table 5. It can be seen that the calculation result of HS is 0.01941 mm. The calculation result of IHS is 0.01938 mm, and the cylinder axis coordinates are 13.004, 5.01066, 14.9803, −5.11827E − 05, −0.000733209, and −0.816712; its result is smaller than HS in the listed algorithms, so the accuracy of the cylindricity error is highest. As can be seen from Figure 15, the iteration curve of the HS and IHS algorithms are, respectively, described. As seen in Figure 8, when the number of iterations is set to 64 times, the IHS algorithm is convergence, and for HS, the number of iterations is 390 times. The IHS has a faster iterative speed than HS.

Table 5: Cylindricity error of data (mm).

Method

HS
13.0038
5.00857
12.642
−3.40985E − 05
−0.000517393
−0.576862
0.01941

IHS
13.004
5.01066
14.9803
−5.11827E − 05
−0.000733209
−0.816712
0.01938

Figure 15: Cylindricity error iteration curve of HS and IHS algorithms.

6.3. Cylindricity Error Experiment 2

To test the performance of the IHS algorithm in the cylindricity evaluation, the data obtained from correlation literature are used to experiment [21]. The evaluation results of the HS and IHS algorithms are compared with the other literatures. Table 6 is reproduced from Lei et al. [21] (under the Creative Commons Attribution License/public domain).

Table 6: Measurement data (mm) [21].

1
−7.169
−34.21
30.5

2
−14.324
−31.887
30.5

3
−25.509
−23.877
30.5

4
−33.683
−9.283
30.5

5
−34.167
7.378
30.5

6
−25.051
24.377
30.5

7
−7.1
34.238
30.5

8
8.071
34.034
30.5

9
19.797
28.86
30.5

10
28.945
19.647
30.5

11
34.891
2.595
30.5

12
34.168
−7.484
30.5

13
28.525
−20.224
30.5

14
18.113
−29.908
30.5

15
7.971
−34.044
30.5

16
−6.778
−34.290
30.5

17
−11.322
−33.073
43.3

18
−20.262
−28.475
43.3

19
−24.016
−25.389
43.3

20
−31.014
−16.098
43.3

21
−34.584
−4.991
43.3

22
−30.614
16.878
43.3

23
−19.46
29.045
43.3

24
−5.802
34.484
43.3

25
10.248
33.444
43.3

26
23.395
26.026
43.3

27
32.027
14.043
43.3

28
34.985
0.374
43.3

29
33.364
−10.449
43.3

30
27.418
−21.716
43.3

31
17.678
−30.164
43.3

32
−0.706
−34.942
43.3

33
−14.539
−31.801
60.8

34
−22.717
−26.563
60.8

35
−26.655
−22.607
60.8

36
−31.59
−14.943
60.8

37
−34.892
1.967
60.8

38
−29.347
18.995
60.8

39
−19.04
29.322
60.8

40
−6.241
34.407
60.8

41
6.174
34.448
60.8

42
14.594
31.802
60.8

43
28.945
19.636
60.8

44
34.887
2.506
60.8

45
31.458
−15.279
60.8

46
20.538
−28.299
60.8

47
6.489
−34.352
60.8

48
−7.37
−34.173
60.8

49
−15.817
−31.168
81

50
−28.573
−20.125
81

51
−33.084
−11.238
81

52
−34.611
4.838
81

53
−25.488
23.924
81

54
−16.38
30.884
81

55
−0.788
34.969
81

56
10.206
33.457
81

57
21.025
27.971
81

58
30.269
17.524
81

59
33.93
8.513
81

60
34.787
−3.557
81

61
31.626
−14.895
81

62
18.362
−29.755
81

63
3.4
−34.794
81

64
−5.318
−34.548
81

65
−8.484
−33.909
97.5

66
−18.607
−29.583
97.5

67
−26.163
−23.175
97.5

68
−33.973
−8.156
97.5

69
−34.49
5.71
97.5

70
−29.015
19.502
97.5

71
−17.564
30.231
97.5

72
−3.842
34.765
97.5

73
15.389
31.424
97.5

74
29.428
18.935
97.5

75
34.054
8.042
97.5

76
34.931
−2.026
97.5

77
33.204
−11.008
97.5

78
29.818
−18.273
97.5

79
21.451
−27.629
97.5

80
15.791
−31.216
97.5

The calculation results of the IHS and other related literature algorithms are listed in Table 7. It can be seen that the calculation results of other algorithms are less than the MIC. The results of the GOS algorithm in the paper are 0.0319 mm. The calculation results of HS and IHS are 0.03204 mm and 0.03188 mm, respectively, and the cylinder axis coordinates of the IHS are 0.0120278, 0.01206, 60.522, −0.000008223, 0.000137151, and −0.853656; their result of IHS is the smallest in the listed algorithms, so the accuracy of the cylindricity error is the highest. The iteration curves of HS and IHS are shown in Figure 16. The convergence times of IHS are at 57 iterations, and HS times is at 365 iterations. The convergence speed of the IHS algorithm is faster than that of the HS algorithm.

Table 7: Cylindricity error of data (mm).

Method

LSC [21]
0.0125
0.0123
30.5
−0.0000040172
0.000013612
1
0.0366

MZC [21]
0.0118
0.0170
30.5
−0.0000025775
−0.00016886
1
0.0324

MIC [21]
0.0121
0.0076
30.5
−0.000012014
0.000074796
1
0.0367

MCC [21]
0.0110
0.0165
30.5
0.000004866
−0.00014041
1
0.0335

GOS [21]
0.0123
0.0170
30.5
−0.000010192
−0.00016365
1
0.0319

HS
0.0121121
0.0113201
65.3792
−0.000005753
0.000109724
−0.69781
0.03204

IHS
0.0120278
0.01206
60.522
−0.000008223
0.000137151
−0.853656
0.03188

Figure 16: Cylindricity error iteration curves of the HS and IHS algorithms.

6.4. Cylindricity Error Experiment 3

To further test the performance of the IHS algorithm in cylindricity error evaluation, the data from other correlation literature are applied to the experiment [3]. The computing results of the HS and IHS algorithms are compared with the other relevant literatures. Table 8 is reproduced from Carr and Ferreira [3] (under the Creative Commons Attribution License/public domain).

Table 8: Measurement data (mm) [3].

1
60.051121
0.002953
3.946134

2
−57.932024
15.399312
15.983017

3
57.43213
17.488707
20.365942

4
55.022756
−23.936632
11.505062

5
29.1801
−52.423113
1.037163

6
−58.861558
−11.113569
20.134482

7
−44.597179
40.113733
2.005267

8
−23.247383
−55.406652
17.669299

9
34.041568
−49.309081
15.807863

10
−34.084135
−49.427745
12.479981

11
50.684216
−32.022045
22.865941

12
57.318676
17.619539
22.082457

13
−40.40813
−44.485701
22.692315

14
−39.83837
44.994386
7.411167

15
−10.261352
−59.146784
22.600675

16
53.919844
26.493193
18.949042

17
−8.540012
59.442972
13.092342

18
−59.369089
8.361285
7.133233

19
−38.029817
46.404843
4.995216

20
47.946099
−35.92538
27.276243

The calculation results of IHS and other related literature algorithms are listed in Table 9. It can be seen that the calculation results of other algorithms are less than the LSM. The calculation results of HS and IHS are 0.18748 mm and 0.183957 mm, respectively, and the cylinder axis coordinates of the IHS are 0.005953, 0.024811, 7.583571, 0.006012, 0.028294, and −9.706070; their result of IHS is the smallest in the listed algorithms, so the accuracy of the cylindricity error is the highest. The iteration curves of HS and IHS are shown in Figure 17. The convergence times of IHS are at 48 iterations, and HS times is at 285 iteration times. The convergence speed of the IHS algorithm is faster than that of the HS algorithm.

Table 9: Cylindricity error of data (mm).

Method

LSM [3]
0.08246
0.11153
0.00101
−0.00490
−0.00546
0.99997
0.21197

LPM [3]
0.01065
0.04692
0.00014
−0.00062
−0.00292
1.00000
0.18396

SA [11]
0.091257
0.096517
−
0.004542
0.004802
0.999978
0.184590

IGA [30]
0.011853
0.047689
0
−0.000674
0.002960
1
0.184274

PSO-DE [31]
0.010650
0.046918
0
−0.000619
−0.002915
1
0.18397196

HS
0.032054
0.0338191
10.996
−0.0536994
−0.0848952
15.5462
0.18748

IHS
0.005953
0.024811
7.583571
0.006012
0.028294
−9.706070
0.183957

Figure 17: Cylindricity error iteration curves of HS and IHS algorithms.

7. Conclusion and Future Work

To improve the evaluation accuracy and accelerate the iteration speed of the cylindricity, the harmony search algorithm is applied to the evaluation of the cylindricity error. An improved harmony search algorithm is proposed to solve the problems of early premature convergence of the basic harmony search algorithm and easy to get into local optimization. The logistic chaotic initialization is introduced into the generation of initial solution to improve the quality of solutions. During the iterative process, the global and local search capabilities are balanced by adopting the par and operators adaptively. After each iteration, the Cauchy mutation strategy is adopted to the best solution to further improve the calculation precision of the IHS algorithm. The results show that the improved harmony search algorithm is very suitable for the accurate evaluation of the cylindricity error.

With the continuous development of big data and artificial intelligence technology, we will carry on a depth analysis to the cylindricity error data of the parts to find the internal relation between the manufacturing information and cylindricity error, so as to further optimize the processing technology and improve the manufacturing accuracy of the parts.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This paper was supported by the National High Technology Research and Development Program of China (no. 2015AA043003).

References

ASME Y 14.5M, Dimensioning and Tolerancing, The American Society of Mechanical Engineers, New York, NY, USA, 1994.

ISO/TS 12180-1, Geometrical Product Specification GSP. Cylindricity, Part 1: Vocabulary and Parameters of Cylindricity Form, International Organization for Standardization, Geneva, Switzerland, 2011.

K. Carr and P. Ferreira, “Verification of form tolerances part II: cylindricity and straightness of a median line,” Precision Engineering, vol. 17, no. 2, pp. 144–156, 1995.

J.-Y. Lai and I.-H. Chen, “Minimum zone evaluation of circles and cylinders,” International Journal of Machine Tools and Manufacture, vol. 36, no. 4, pp. 435–451, 1996.

H.-Y. Lai, W.-Y. Jywe, C.-K. Chen, and C.-H. Liu, “Precision modeling of form errors for cylindricity evaluation using genetic algorithms,” Precision Engineering, vol. 24, no. 4, pp. 310–319, 2000.

C.-C. Cui, F.-G. Huang, R.-C. Zhang, and B. Li, “Research on cylindricity evaluation based on the particle swarm optimization algorithm,” Optics and Precision Engineering, vol. 14, no. 2, pp. 254–260, 2006.

J. Mao, H.-W. Zheng, Y.-L. Cao, and S. Xu, “Method for cylindricity errors evaluation using particle swarm optimization algorithm,” Transactions of the Chinese Society for Agricultural Machinery, vol. 38, no. 2, pp. 146–149, 2007.

Z. W. Geem, J. H. Kim, and G. V. Loganathan, “A new heuristic optimization algorithm: harmony search,” Simulation, vol. 76, no. 2, pp. 60–68, 2001.

S. O. Degertekin, “Improved harmony search algorithms for sizing optimization of truss structures,” Computers and Structures, vol. 92-93, no. 3, pp. 229–241, 2012.

D.-X. Zou, L.-Q. Gao, J.-H. Wu, S. Li, and Y. Li, “A novel global harmony search algorithm for reliability problems,” Computers and Industrial Engineering, vol. 58, no. 2, pp. 307–316, 2010.

S.-Y. Chou and C.-W. Sun, “Assessing cylindricity for oblique cylindrical features,” International Journal of Machine Tools and Manufacture, vol. 40, no. 3, pp. 327–341, 2000.

T. Weber, S. Motavalli, B. Fallahi, and S.-H. Cheraghid, “A unified approach to form error evaluation,” Precision Engineering, vol. 26, no. 3, pp. 269–278, 2002.

L.-M. Zhu and H. Ding, “Application of kinematic geometry to computational metrology: distance function based hierarchical algorithms for cylindricity evaluation,” International Journal of Machine Tools and Manufacture, vol. 43, no. 2, pp. 203–215, 2003.

Y.-Z. Lao, H.-W. Leong, F. P. Preparata, and G. Singh, “Accurate cylindricity evaluation with axis-estimation preprocessing,” Precision Engineering, vol. 27, no. 4, pp. 429–437, 2003.

N. Venkaiah and M. S. Shunmugam, “Evaluation of form data using computational geometric techniques—part II: cylindricity error,” International Journal of Machine Tools and Manufacture, vol. 47, no. 7-8, pp. 1237–1245, 2007.

G.-X. Bei, P.-H. Lou, X.-Y. Wang, H.-Y. Zhu, and H. Du, “Cylindricity error evaluation based on genetic algorithms,” Journal of Shandong University, vol. 38, no. 2, pp. 33–36, 2008.

H. Guo, D.-J. Lin, J.-Z. Pan, and S.-W. Jiang, “Cylindricity error evaluation based on multi-population genetic algorithm,” Journal of Engineering Graphics, vol. 29, no. 4, pp. 48–53, 2008.

J.-S. Li, X.-Q. Lei, Y.-J. Xue, and M.-D. Duan, “Evaluation algorithm of cylindricity error based on coordinate transformation,” China Mechanical Engineering, vol. 20, no. 16, pp. 1983–1987, 2009.

J. Luo, J.-J. Lu, W.-M. Chen et al., “Cylindricity error evaluation using artificial bee colony algorithm with tabu strategy,” Journal of Chongqing University, vol. 32, no. 12, pp. 1482–1485, 2009.

X.-L. Wen, J.-C. Huang, D.-H. Sheng, and F.-L. Wang, “Conicity and cylindricity error evaluation using particle swarm optimization,” Precision Engineering, vol. 34, no. 2, pp. 338–344, 2010.

X.-Q. Lei, H.-W. Song, Y.-J. Xue, J.-S. Li, J. Zhou, and M.-D. Duan, “Method for cylindricity error evaluation using geometry optimization searching algorithm,” Measurement, vol. 44, no. 9, pp. 1556–1563, 2011.

K. Lee, S. Cho, and S. Asfour, “Web-based algorithm for cylindricity evaluation using support vector machine learning,” Computers and Industrial Engineering, vol. 60, no. 2, pp. 228–235, 2011.

X.-L. Wen, Y.-B. Zhao, D.-X. Wang, and J. Pan, “Adaptive Monte Carlo and GUM methods for the evaluation of measurement uncertainty of cylindricity error,” Precision Engineering, vol. 37, no. 4, pp. 856–864, 2013.

G.-Y. He, P.-P. Liu, and K. Wang, “Evaluation of the cylindricity error based on the sequential quadratic programming algorithm,” Mechanical Science and Technology for Aerospace Engineering, vol. 33, no. 12, pp. 1845–1849, 2014.

G. Hermann, “Application of computational geometry in coordinate measurement,” in Proceedings of the 10th Jubilee IEEE International Symposium on Applied Computational Intelligence and Informatics, pp. 511–516, IEEE, Timisoara, Romania, May 2015.

P. Zheng, J.-Q. Wu, and L.-N. Zhang, “Research of the on-line evaluating the cylindricity error technology based on the new generation of GPS,” Procedia Engineering, vol. 174, pp. 402–409, 2017.

B. Li and W.-S. Jiang, “Chaos optimization method and its application,” Control theory and application, vol. 14, no. 4, pp. 613–615, 1997.

A. Kanso and N. Smaoui, “Logistic chaotic maps for binary numbers generations,” Chaos Solitons and Fractals, vol. 40, no. 5, pp. 2557–2568, 2009.

L.-G. Wang, Y. Hong, and D.-M. Yu, “Simplified artificial fish swarm algorithm,” Journal of Chinese Computer Systems, vol. 30, no. 8, pp. 1663–1667, 2009.

X.-L. Wen and A.-G. Song, “An improved genetic algorithm for planar and spatial straightness error evaluation,” International Journal of Machine Tools and Manufacture, vol. 43, no. 11, pp. 1157–1162, 2003.

X.-C. Zhang, X.-Q. Jiang, and P. J. Scott, “A reliable method of minimum zone evaluation of cylindricity and conicity from coordinate measurement data,” Precision Engineering, vol. 35, no. 3, pp. 484–489, 2011.