Mathematical Problems in Engineering
Volume 2021 (2021), Article ID 5363992, 10 pages
https://doi.org/10.1155/2021/5363992

Research Article

Simulation of Constant-Volume Removal Rate Machining of Middle-Convex and Varying Ellipse Piston

Yanwei Xu,1,2 Yinhao Wang,¹ and Tancheng Xie1,2

¹School of Mechatronics Engineering, Henan University of Science and Technology, Luoyang 471003, China
²Intelligent Numerical Control Equipment Engineering Laboratory of Henan Province, Luoyang 471003, China

Correspondence should be addressed to Yanwei Xu; xuyanweiluoyang@163.com

Received 7 May 2021; Revised 19 May 2021; Accepted 21 May 2021; Published 30 May 2021

Academic Editor: Qiang Guo

Copyright © 2021 Yanwei Xu 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.

One kind of constant-volume removal rate machining method of the middle-convex and varying ellipse piston is proposed in this paper. By analyzing the structure and movement relationship of the middle-convex and varying ellipse piston machine, the NC machining model is built. And, the constant-volume removal rate machining model is also built by superposing the variable rotation satisfying the dynamic performance constraints on the uniform rotation of the spindle of the CNC piston lathe. Then, the instantaneous position parameters of each axis of the CNC piston lathe are obtained and turned into NC code. The functional feasibility of the method finally is verified by simulation machining.

1. Introduction

In the reciprocating process of crankshaft output torque, the piston transmits the gas pressure in the cylinder to this process through the piston pin and connecting rod, but the clearance between the piston and the cylinder sleeve should be as small and uniform as possible to prevent the piston from being “roughened” or “bitten” in the cylinder. However, in the coupled working environment of high temperature, high pressure, and alternating mechanical and thermal loads, the force acting on the top of the piston leads to the deformation of the piston skirt along the axis of the pin seat, the deformation of the side of the piston skirt, and the uneven thermal expansion. These causes lead to serious and irregular deformation of the piston, which makes the cross section of the piston skirt that is originally a round cross section cannot keep round in the working environment, and then become oval, resulting in that the “cylindrical” piston processed at room temperature presents the “elliptical cylinder” shape. This deformation directly affects the uniformity of fit clearance between piston and cylinder liner. The piston with middle-convex and varying ellipse piston skirt is designed, which can keep the ideal geometry under working conditions [1], ensure the good fit between the piston and the cylinder wall, reduce the cylinder clearance and the impact of the piston on the cylinder wall, and reduce the specific pressure and noise of the piston skirt. Moreover, a wedge-shaped oil gap is formed between the middle-convex skirt of the piston and the cylinder wall, which ensures the good lubrication of the piston and reduces the wear of parts.

For a long time, many scholars have deeply studied the machining method of middle-convex and varying ellipse piston [2–12].These achievements have effectively promoted the progress of machining technology of middle-convex and varying ellipse piston, but the noncircular section machining technology of piston skirt studied by these achievements is basically based on numerical control turning method. NC turning is an important method for machining the piston with middle-convex and varying ellipse piston, which can effectively improve the machining accuracy and efficiency of the piston with middle-convex and varying ellipse piston. However, in the process of cutting a circular cross-sectional workpiece into an elliptical cross-sectional workpiece on a CNC piston lathe, the workpiece will rotate at a constant speed, and the cutting tool will move with high frequency and reciprocation. Moreover, the cutting depth of the cutting tool will change continuously with the rotation position of the workpiece, and the cutting area of the cutting tool in the same cross section will change with the processing position. The material removal rate and cutting force of the workpiece change periodically, which will form dynamic excitation force and affect the machining accuracy of piston skirt to a certain extent.

Based on the basic principle that the cutting force is approximately proportional to the cutting area, this paper studies the numerical control machining technology and simulation of the equal volume cutting rate of the middle-convex and varying ellipse piston. On the premise that the rotation position of the shaft section of the piston skirt and the cutting depth of the tool tip conform to the ellipse trajectory, the variable-speed rotation meeting the dynamic performance constraint is superimposed on the uniform speed rotation of the main shaft of the piston machine tool. The rotation speed of the main shaft of the machine tool is continuously regulated, and the material removal rate and cutting force of the middle-convex and varying ellipse piston skirt are basically constant.

2. Processing Principle of Middle-Convex and Varying Ellipse Piston Skirt

The profile characteristics of the middle-convex and varying ellipse piston skirt are as follows: the generatrix of the piston skirt is a middle-convex curve which represents the variation of the diameter of the long axis of the elliptic section along the direction of the piston axis. The shape of cross section of piston skirt is similar to an ellipse, and the ovality of different sections is different. In any section, the diameter of long axis of ellipse section is the largest diameter, the long axis is along the axis of pin seat, the diameter of short axis of this section is the smallest diameter, and the short axis is perpendicular to the axis of pin seat. The ellipticity of the elliptic section is the difference between the diameter of major axis and that of minor axis.

2.1. Elliptical Profile of Skirt Cross Section of Middle-Convex and Varying Ellipse Piston

The geometric characteristics of the skirt cross section of the piston with middle-convex and varying ellipse can be described as follows [2]:

In the formula, is polar radius value;is ellipticity; and is elliptical long half-axis.

In the formula, is long half-axis of section with maximum ellipticity.

In the formula,is the maximum ellipticity; isthe relative rotation angle of polar radius to long axis; is the skirt height; is the skirt height of section with maximum ellipticity; is the dimensionless coefficient ( is an ellipse; is a quadratic ellipse); is the profile shape characteristic coefficient of longitudinal profile; is the dimensionless coefficient; and is the ellipticity change rate.

In the formula, is the dimensionless coefficient (when , the cross section is circular; when , the cross section is elliptical).

Formula (1) is the variation of polar radius, and formula (2) is the radial variation of the middle-convex profile along the piston axis. If the piston is an elliptical piston, the formula of any elliptical cross section of the piston skirt can be obtained by substituting formula (4) into (1).

among which, , is angular velocity, and

2.2. Middle-Convex Profile of Middle-Convex and Varying Ellipse Piston

The middle-convex profile of the middle-convex and varying ellipse piston skirt is usually given discrete points in design. Figure 1 is the design parameters of the middle-convex profile of the Perkins 240 piston skirt.

Figure 1: Geometric parameters of the Perkins 240 piston skirt.

Usually, the discrete points of the convex profile are fitted to smooth curves by cubic spline interpolation, and then the equation of the middle-convex profile is obtained. In the process of fitting the discrete points of middle-convex profile in piston skirt shown in Figure 2, N discrete points are put into the XOZ coordinate system in which the X-axis is parallel to the piston cross section and the Z-axis is the piston axis (X represents the long half-axis of elliptical cross section and Z represents the piston skirt height). If the first derivatives of the curve composed of discrete points are and at the beginning and end points, respectively, the function value x (z) of any point on the Z-axis can be expressed as follows [13]:

Figure 2: Fitting of discrete points of middle-convex profile in piston section.

In the formula, satisfies the equation

among which, , , and . Solution , tdhe fitting equation x (z) of middle-convex profile in piston skirt can be obtained by substitution formula (1).

2.3. Analysis of Processing Process of Middle-Convex and Varying Ellipse Piston Skirt

The structure of CNC piston machine tool is shown in Figure 3. Machine motion consists of workpiece (spindle C-axis) rotary motion, tool holder slider linear motion along the Z-axis (parallel to workpiece rotary axis), tool holder linear motion along the X axis (perpendicular to workpiece rotary axis), and tool holder U-axis reciprocating linear motion (control tool high-frequency reciprocating microdisplacement linear motion), X-axis. The cutting depth is controlled by U-axis, and the travel of the tool in high-speed reciprocating linear motion depends on the ellipticity of elliptical cross section at different skirt heights of piston skirt.

Figure 3: Structure sketch of CNC piston machine tool.

The process of turning skirt profile of middle-convex and varying ellipse piston by motion synthesis method can be divided into two independent motions: (1) the motion of tool relative to workpiece forming middle-convex profile is realized by the joint motion of X-axis and Z-axis driven by servo motor; (2) the motion of tool relative to workpiece forming elliptical profile, the reciprocating linear motion of U-axis driven by linear motor, and the motion of spindle.

2.3.1. Middle-Convex Profile Processing

The motion of forming the middle-convex profile is synthesized by the motion of the tool in the piston axis (Z-axis) and the radial direction (X-axis). After the fitting equation of discrete points on the middle-convex profile is obtained, the middle-convex profile is interpolated, and the cubic spline curve is approximated and fitted by micro line segments. According to the requirements of piston skirt surface processing accuracy and tool feed, the number of interpolation points is determined, and the middle-convex profile is interpolated in the piston axis direction. In this paper, the interpolation points are divided by equal interval method. Each step of interpolation, the workpiece rotates one week to complete a micro-short elliptical cylinder processing. Then, the Z-axis position is calculated, and the Z-axis servo motor drives the trawler to the next interpolation point. Figure 4 is a schematic diagram of NC machining of middle-convex profile in piston skirt.

Figure 4: NC machining schematic diagram of middle-convex profile.

2.3.2. Elliptical Profile Processing

The forming motion of the elliptical section of the piston skirt can be decomposed into a high-speed rotational motion of the workpiece and a high-speed reciprocating linear feed motion of the tool in the radial direction of the piston. During the forming process, the piston rotates with the spindle every one revolution, and the tool sequentially processes the long axis of the elliptical section ⟶ the short axis ⟶ the long axis ⟶ the short axis ⟶ the long axis, and the tool feeds twice in rapid reciprocating direction. The higher the spindle speed, the higher the tool feed frequency; the greater the ellipticity of the elliptical section, the greater the turning radius change, the greater the displacement of the tool’s fast reciprocating linear feed, and the greater the speed and acceleration of the tool. The key to the elliptical section turning of the piston skirt is the control of the tool path: (1) the high-frequency reciprocating linear motion of the tool; (2) the microdisplacement of the tool reciprocating linear feed motion and the angular displacement of the workpiece high-speed rotation maintain a strict one-to-one correspondence. The schematic diagram of the elliptical section machining process of the piston is shown in Figure 5.

Figure 5: Schematic diagram of elliptical section processing.

Assuming that the starting position of the tool tip is located at the apex B of the ellipse long axis, the expression of the motion displacement of the tool tip vertex can be described as follows:

In the formula, is the rotation frequency of the lathe spindle and t is time.

When the workpiece rotates with the spindle at a constant speed to turn the piston skirt contour with middle-convex and varying ellipse, the displacement relationship of the four moving axes of the machine tool can be expressed as follows:

This is the NC machining model of the middle-convex and varying ellipse piston skirt in the uniform turning of the workpiece.

2.4. Distribution Law of Cutting Area of Elliptical Section

In the NC machining process of the middle-convex and varying ellipse piston, the rough machining of the piston is completed on an ordinary CNC lathe, and the final forming finishing is completed by a numerical control piston lathe. During the conventional forming and finishing of the piston skirt, the workpiece rotates at a constant speed. The area cut-off in the same elliptical cross section varies with the machining position.

In Figure 5, it is assumed that the cutting depth is changed from the workpiece diameter dimension A of the previous process to the long axis B of the elliptical cross section, and the cutting depth is ap from 0° to 90°, and the tool is continuously advanced, 90°. When the time is fed to the short axis K of the elliptical cross section, the path of the cutting edge along the surface of the workpiece is substantially like a smooth curve BDGK. During the process of moving the workpiece from 0° to 90°, the area cut by the workpiece through the same corner gradually increases. When the workpiece is turned from 0° to angle, the area to be cut is

In the formula, the elliptical sector area is calculated as follows:

When the workpiece is turned from 1 to 2, the area to be cut is

Select the Perkins 240 piston skirt height H = 20 elliptical section in Figure 1 (diameter d = Φ92, corresponding to ellipticity G = 0.20). Take the cutting depth ap = 0.1. Calculate the workpiece in the first quadrant every revolution. The depth of cut ap and its variation ap and the area cut within the same corner. The results are shown in Table 1.

Table 1: Cutting depth and its variation and the area cut in the same corner.

Angle (°)
ap (mm)
ap (µm)
Equal-angle cut area (mm²)

3
0.10027
0.27
0.24134

6
0.10109
0.82
0.24266

9
0.10244
1.35
0.24529

12
0.10432
1.87
0.24919

15
0.10669
2.37
0.25433

18
0.10954
2.85
0.26065

21
0.11283
3.29
0.26807

24
0.11653
3.70
0.27652

27
0.12059
4.06
0.28590

30
0.12498
4.39
0.29611

33
0.12964
4.66
0.30703

36
0.13452
4.88
0.31855

39
0.13958
5.05
0.33054

42
0.14475
5.17
0.34286

45
0.14997
5.23
0.35538

48
0.15520
5.23
0.36797

51
0.16037
5.17
0.38048

54
0.16543
5.06
0.39278

57
0.17031
4.89
0.40473

60
0.17498
4.67
0.41621

63
0.17937
4.39
0.42709

66
0.18344
4.07
0.43724

69
0.18715
3.70
0.44657

72
0.19044
3.30
0.45496

75
0.19329
2.85
0.46233

78
0.19567
2.38
0.46860

81
0.19755
1.88
0.47370

84
0.19891
1.36
0.47758

87
0.19973
0.82
0.48018

90
0.2
0.27
0.48149

It can be seen from Table 1 that in the first quadrant, the cutter starts to cut from the long semiaxis of the elliptical section, and the cutting depth ap gradually increases as the workpiece rotation angle changes; the amount of change in the cutting depth gradually increases at the beginning of the angle ap at 45°. The vicinity ap increases to the maximum value and then gradually decreases; the area cut by the workpiece rotating through the same angle (3°) gradually increases. In the first quadrant, the cutting depth ap, the cutting depth variation ap, and the change trend of the cutting area in the workpiece rotating through the same angle (3°) are shown in Figures 6–8, respectively.

Figure 6: Cutting depth ap change trend.

Figure 7: Change in cutting depth ap change trend.

Figure 8: Area cut in the same corner.

3. Processing Principle of Equal-Volume Excision Rate

During the process of uniform turning of the elliptical cross section of the middle-convex and varying ellipse piston skirt during uniform turning, the resection area changes periodically with the workpiece rotation angle per unit time, resulting in periodic changes in material removal rate and cutting force, and processing of the piston skirt to some extent. Accuracy has an impact. In order to reduce this effect, this paper proposes the NC machining concept of equal-volume resection rate of middle-convex and varying ellipse piston—the volume ρV of the material removed per unit time is equal during the elliptical cross section of the middle-convex and varying ellipse piston skirt during processing.

The projected area of the volume on the elliptical section is

The projected area of the maximum volume allowed to be cut per unit time in the elliptical section is

In the schematic diagram of the elliptical cross section machining of the middle-convex and varying ellipse piston shown in Figure 4, the total area cut-off in the first quadrant is equal to the difference between the quarter circle area SAOI and the quarter ellipse area SBOK. Then, there is

The total area cut out in the first quadrant is divided into n equal parts for processing, so that the time taken to cut each aliquot area is the same.

In the formula, n must meet the conditions

In the first quadrant, the area cut-off from is equal to. Then, there is

Combining (17) and (19) gives

Solving the equation yields angle values of

When the elliptical section of the middle-convex and varying ellipse piston skirt is machined by the equal-volume resection rate method, the displacement relationship of the four axes of motion of the machine tool can be expressed as follows:

This is the numerical control machining model when machining the middle-convex and varying ellipse piston skirt in the same volume resection rate method.

4. Equal-Volume Resection Rate Simulation Processing and Experiment

The skirt height H of the middle-convex and varying ellipse piston in Perkins 240 is shown in Figure 1 and its corresponding elliptical cross-sectional long axis value and ellipticity value are shown in Tables 2 and 3, respectively.

Table 2: Piston skirt height H and corresponding elliptical section long axis value.

Skirt height H (mm)
Long axis (mm)
Skirt height H (mm)
Long axis (mm)

4
91.9835
36
91.969

8
91.9911
40
91.948

12
91.9951
44
91.922

16
91.9976
48
91.8924

20
92
52
91.859

24
91.9972
56
91.8191

28
91.9925
60
91.7693

32
91.9829
62
91.737

Table 3: Piston skirt height H and corresponding cross-sectional ellipticity values.

Skirt height H
Cross-sectional ellipticity

0–41
0.2

41–62
Linear gradient to 0.3

According to Table 2, using cubic spline interpolation, the fitting equation and fitting curve of the ellipse long axis value of the Perkins 240 piston skirt can be obtained by MATLAB program fitting (Figure 9).

Figure 9: Contour line fitting curve in the Perkins 240 piston skirt.

According to Table 3, using linear interpolation, the equation for the ellipticity of the elliptical section G with the height of the skirt is

Take the elliptical cross-sectional long axis maximum value plus twice the depth of cut ap as the cylinder workpiece diameter d₀ before the piston skirt forming process, starting from the skirt height H = 0, according to the feed amount f = 0.001 [14], calculate the different ellipse values of the long semiaxis a, the ellipticity G, and the short semiaxis b of the section, and calculate the feed of the tool from the skirt height H = 0, the tool is cut from d₀/2 to the long semiaxis of each different elliptical section. For each different elliptical section, divide the area cut-off in the first quadrant into n equal parts, and obtain the angle corresponding to each aliquot and the corresponding tool feed amount; then, according to the principle of symmetry, the workpiece is rotated within one week. The angle corresponding to each aliquot and the tool feed. Then, the numerical values of workpiece rotation angle and tool feed per equal part are transformed into corresponding NC machining program. According to the NC machining model of equal volume removal rate for middle-convex and varying ellipse piston, the piston skirt can be machined with equal-volume removal rate.

Taking the maximum elliptical section of the Perkins 240 piston as an example, the skirt height H = 20, the long axis diameter, the ellipticity G = 0.20, the long half axis a = 46, the short half-axis b = 45.9, and the first quadrant 0-degree cutting allowance ap = 0.1. According to these parameters, the angle corresponding to the area of each aliquot and the actual rotation angle of the workpiece and the corresponding depth of cut ap and its variation ap can be obtained. Take n = 30 aliquots, and the calculation results are shown in Table 4.

Table 4: Equal-volume resection rate processing calculation results.

No.
ap (mm)
ap (µm)
Every aliquot (°)
Actual corner (°)

1
0.1006
0.612
4.4897
4.4897

2
0.1024
1.793
4.4358
8.9255

3
0.1053
2.852
4.3361
13.2617

4
0.1090
3.74
4.2015
17.4632

5
0.1134
4.434
4.0463
21.5094

6
0.1184
4.939
3.8818
25.3912

7
0.1236
5.275
3.7174
29.1085

8
0.1291
5.467
3.5592
32.6678

9
0.1347
5.542
3.4102
36.0780

10
0.1402
5.525
3.2733
39.3513

11
0.1456
5.435
3.1478
42.4991

12
0.1509
5.289
3.0338
45.5330

13
0.1560
5.102
2.9313
48.4642

14
0.1609
4.882
2.8390
51.3032

15
0.1655
4.639
2.7565
54.0597

16
0.1699
4.377
2.6826
56.7423

17
0.1740
4.102
2.6167
59.3590

18
0.1778
3.817
2.5583
61.9173

19
0.1813
3.524
2.5061
64.4234

20
0.1846
3.227
2.4603
66.8837

21
0.1875
2.925
2.4196
69.3033

22
0.1901
2.621
2.3841
71.6873

23
0.1924
2.315
2.3537
74.0410

24
0.1944
2.007
2.3274
76.3684

25
0.1961
1.699
2.3056
78.6740

26
0.1975
1.391
2.2878
80.9618

27
0.1986
1.081
2.2735
83.2353

28
0.1994
0.772
2.2632
85.4985

29
0.1998
0.463
2.2557
87.7542

30
0.2
0.154
2.2458
90.0000

According to Table 4, when the maximum elliptical cross section of the middle-convex and varying ellipse piston skirt in Perkins 240 is processed by equal-volume resection; in the first quadrant, the change trend of the workpiece turning angle of each aliquot is as shown in Figure 10. The trend of the actual turning angle is shown in Figure 11; the change trend of each cutting area corresponding to the cutting depth is shown in Figure 12, and the corresponding cutting depth is shown in Figure 13.

Figure 10: Each resection area corresponds to the workpiece rotation angle.

Figure 11: The cut area corresponds to the actual rotation angle of the workpiece.

Figure 12: Cut area corresponds to the depth of cut.

Figure 13: Cut area corresponds to the change in the depth of cut.

The obtained actual cutting angle and the tool feed amount corresponding to each aliquot area and the positional parameters of the corresponding X-axis and Z-axis are converted into a numerical control machining program, and the elliptical section of the middle-convex and varying ellipse piston skirt can be obtained. Perform equal-volume resection rate processing. Skirt height H = 20 elliptical section first image inner limit equal volume resection rate machine tool motion parameters when middle-convex and varying ellipse piston skirt is processed as shown in Table 5.

Table 5: Machine tool motion parameters for equal-volume resection rate machining in the first quadrant.

No.
Z (mm)
X (mm)
(°)
U (mm)

N 3310
20.00
46.100
4.4897
0.1006

N 3320
20.00
46.100
8.9256
0.1024

N 3330
20.00
46.100
13.2617
0.1053

N 3340
20.00
46.100
17.4632
0.1090

N 3350
20.00
46.100
21.5094
0.1134

N 3360
20.00
46.100
25.3912
0.1184

N 3370
20.00
46.100
29.1086
0.1236

N 3380
20.00
46.100
32.6678
0.1291

N 3390
20.00
46.100
36.0780
0.1347

N 3400
20.00
46.100
39.3513
0.1402

N 3410
20.00
46.100
42.4991
0.1456

N 3420
20.00
46.100
45.5330
0.1509

N 3430
20.00
46.100
48.4642
0.1560

N 3440
20.00
46.100
51.3032
0.1609

N 3450
20.00
46.100
54.0597
0.1655

N 3460
20.00
46.100
56.7423
0.1699

N 3470
20.00
46.100
59.3590
0.1740

N 3480
20.00
46.100
61.9173
0.1778

N 3490
20.00
46.100
64.4234
0.1813

N 3500
20.00
46.100
66.8837
0.1846

N 3510
20.00
46.100
69.3033
0.1875

N 3520
20.00
46.100
71.6873
0.1901

N 3530
20.00
46.100
74.0410
0.1924

N 3540
20.00
46.100
76.3684
0.1944

N 3550
20.00
46.100
78.6740
0.1961

N 3560
20.00
46.100
80.9618
0.1975

N 3570
20.00
46.100
83.2353
0.1986

N 3580
20.00
46.100
85.4985
0.1994

N 3590
20.00
46.100
87.7542
0.1998

N 3600
20.00
46.100
90.0000
0.2000

According to the machine motion parameters, based on the VERICUT CNC machining simulation platform, the equal-volume resection rate simulation was performed on the Perkins 240 piston skirt in Figure 1 (as shown in Figure 14).

Figure 14: Perkins 240 piston skirt equal volume resection rate simulation processing.

The ellipticity value of the geometric parameters of the middle-convex and varying ellipse piston skirt in Perkins 240 is magnified 30 times, the remaining parameters are unchanged, recalculated, and simulated, and the simulation processing effect of the middle-convex line and the variable elliptic section can be clearly seen (as shown in Figure 15).

Figure 15: Perkins 240 piston skirt ovality magnification 30 times simulation processing.

The simulation processing of Perkins 240 piston meets the functional verification requirements of the middle-convex and varying ellipse piston machine. The results show that the simulation machining can simulate the actual machining process of the medium-convex elliptical piston machine tool and also verify the volume of the middle-convex and varying ellipse piston skirt.

By the maximum speed 2500 r/min of the workpiece and the feed rate 0.001 mm/r, the equal-volume cutting rate machining test was carried out on Perkins 240 piston based on the middle-convex and varying ellipse piston NC machining test platform. The machining test results show that the simulation can meet the requirements of the actual machining process of the medium-convex elliptical piston machine tool, and the method is feasible. The experimental process is shown in Figure 16.

Figure 16: The experimental process of the Perkins 240 piston skirt.

5. Conclusion

This paper analyzes the forming principle of the oval section profile of the middle-convex and varying ellipse piston skirt. The middle-convex profile and ellipse profile of the middle-convex and varying ellipse piston are described mathematically, and the NC machining model of turning piston skirt is established. The distribution of the cutting area in the process of the ellipse section is analyzed. The numerical control machining model of equal-volume cutting rate for middle-convex and varying ellipse piston is established, and the machining method of equal volume cutting rate for middle-convex and varying ellipse piston is verified by simulation.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors are grateful for the financial support of the National Natural Science Foundation of China under Grant no. 51805151 and the Key Scientific Research Project of the University of Henan Province of China under Grant no. 21B460004.

References

A. Li, J. Zhao, Z. Gong et al., “Optimal selection of cutting tool materials based on multi-criteria decision-making methods in machining Al-Si piston alloy,” International Journal of Advanced Manufacturing Technology, vol. 86, no. 1-4, pp. 1055–1062, 2016.

S. Jiang, X. Wang, Z. Shi et al., “Study on the high performance linear servo system for middle-converx and varying ellipse piston machining,” Chinese Journal of Mechanical Engineering, vol. 37, no. 9, pp. 58–61, 2001.

S. Pan, R. Zhang, M. Cheng et al., “The shaping method of middle convex and varied elliptic pattern with NC turning,” China Mechinacl Eegineering, vol. 12, no. 8, pp. 911–914, 2001.

X. Wang, S. Yong, C. Liu et al., “The analysis for real-time performance and realization of CNC system on middle-convex and varying ellipse piston machining,” Manufacturing Technology and Machine Tool, no. 3, pp. 15–18, 2004.

Y. Qin, X. Song, Q. Wang et al., “Research on error control method of noncircular piston turning,” China Mechinacl Eegineering, vol. 18, no. 5, pp. 517–519, 2007.

B. Zhang and T. Zhu, “Research on CNC lathe for middle-convex and varying ellipse piston,” Machine Tool and Hydraulics, vol. 38, no. 20, pp. 15–17, 2010.

X. Ma, B. Zhang, and H. Xu, “NC turning processing principle and realization for convex oval piston,” Machine Building and Automation, vol. 41, no. 5, pp. 63–65, 2012.

K. Liu, Z. Wang, and G. Gui, “Processing technology of middle-convex and varying ellipse piston turning,” Journal of Hefei University of Technology (Natural Science), vol. 35, no. 10, pp. 1316–1319, 2012.

C. Gong, “The research and development of CNC system for varying ellipse piston lathe,” Huazhong University of Science and Technology, Wuhan, China, 2012, Master Thesis.

X. Tian, “Variable oval convex surface of the piston-type CNC turning technology innovation into research,” Shandong University, Jinan, China, 2014, Master Thesis.

H. Wu, T. Jiang, X. Chang et al., “Optimizing process algorithm and calculating data points of the skirt part of piston that combination of elliptical part and eccentric circular,” Journal of Xiamen University (Natural Science), vol. 53, no. 2, pp. 212–216, 2014.

S. Yun, “Research on the key technology of middle-convex and varying ellipse piston machining,” Jiangsu University of Science and Technology, Zhen Jiang, China, 2013, Thesis.

Y. Zhang, Y. Huang, W. J. Shao et al., “Research on MCVE piston machining and process parameter optimization,” International Journal of Advanced Manufacturing Technology, vol. 93, no. 9–12, pp. 3955–3966, 2017.

L. Peng and H. Zhang, Quick Search Manual for Machining Parameters, Chemical Industry Press, Beijing, China, 2010.