Next Article in Journal
A Digital Microfluidic RT-qPCR Platform for Multiple Detections of Respiratory Pathogens
Next Article in Special Issue
Research on Continuous Error Compensation of a Sub-Arc-Second Macro/Micro Dual-Drive Rotary System
Previous Article in Journal
Microfluidic Synthesis, Doping Strategy, and Optoelectronic Applications of Nanostructured Halide Perovskite Materials
Previous Article in Special Issue
Design and Performance Research of a Precision Micro-Drive Reduction System without Additional Motion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Surface Roughness in Gas-Solid Two-Phase Abrasive Flow Machining Based on Multivariate Linear Equation

1
School of Mechanical Engineering, Shandong University of Technology, Zibo 255000, China
2
Tianrun Industrial Technology Co., Ltd., Weihai 264200, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(10), 1649; https://doi.org/10.3390/mi13101649
Submission received: 6 September 2022 / Revised: 25 September 2022 / Accepted: 28 September 2022 / Published: 30 September 2022
(This article belongs to the Special Issue Ultra-Precision Manufacturing Technology)

Abstract

:
The main purpose of this study is to explore a surface roughness prediction model of Gas-Solid Two-Phase Abrasive Flow Machining. In order to achieve the above purpose, an orthogonal experiment was carried out. Q235 steel as processing material and white corundum with different particle sizes as abrasive particles were used in the experiment. Shape and spindle speed were the main reference factors. The range method and factor trend graph are used to comprehensively analyze the experimental results of different processing stages of the detection point, and the optimal parameter combination of A3B2C1D2 was obtained. According to the experimental results, a multiple linear regression equation was established to predict the surface roughness, and the experimental results were solved and significantly analyzed by software to obtain a highly reliable prediction model. Through experiments, modeling and verification, it is known that the maximum error between the obtained model and the actual value is 0.339 μm and the average error is 0.00844 μm, which can better predict the surface roughness of the gas-solid two-phase flow abrasive pool.

1. Introduction

Abrasive Flow Machining (AFM) is a new finishing method for deburring and rounding workpieces by extruding fluid abrasives. The polishing technology was originally developed according to the deburring needs of aerospace parts in the US aerospace industry, mainly for the finishing of narrow slits, tiny holes and special-shaped holes. After successful research and development, its application range has gradually expanded [1]. Since the abrasive medium is fluid, AFM has a good application in micropores, complex pores and other pore structures [2,3,4,5]. In the process of AFM, the abrasive particles are relatively small, meaning less damage to the surface of the workpiece. The residual stress on the surface of the workpiece can be reduced, and the concentrated stress on the surface of the workpiece can be reduced while the machining quality is improved, and the wear resistance and corrosion resistance of the workpiece can be increased [6,7]. However, the fluid modulus used in some AFM technologies is usually viscous abrasives, which will encounter greater resistance when the abrasives move relative to each other, and the plasticity and cohesion of some polymer media will cause the abrasives to avoid the processing area, resulting in the surface processing quality of the material being poor. Therefore, in the field of abrasive pools, new grinding media can be developed, and the reduction in media cost, media sustainability, and waste generation still needs to be continuously explored [8].
In response to the exploration of new processing methods of AFM, Barletta [9] innovatively proposed the fluidized bed to assist the spindle-type tumbling finishing and found that the fluidized abrasive grains can reduce the energy required to move the workpiece, and the fluidized abrasive grains can machine the ductile metal with a better mirror finish; for the study of new abrasive media for AFM, Pham et al. [10] explored an abrasive air jet technology to process alloy materials, which reduced the surface roughness Ra by 40 times. The gas-solid two-phase flow abrasive flow processing proposed in this paper is to use high-speed airflow to accelerate the solid-phase abrasive particles to realize the relative motion between the gas-solid two-phase flow and the workpiece, so as to achieve the purpose of finishing (and solid phase abrasive particles are innovative).
However, the flow state of gas-solid two-phase flow, workpiece shape and abrasive grain shape will affect the machining quality of the workpiece surface, and the surface roughness depends on many machining parameters. Therefore, many researchers try to analyze and predict the roughness of the machined surface through modeling. For the prediction of the surface roughness of machined workpieces, Artificial Neural Network (ANN) [11,12,13,14], Adaptive Neuro-Fuzzy Inference System (ANFIS) [15,16], Multi-Objective Genetic Algorithm (MOGA) [17], quadratic regression models [18], and other common techniques are used [19]. Mirko Ficko et al. predicted the surface roughness of abrasive water jet machining based on artificial neural network, and used k-fold cross-validation to verify the ANN model, which greatly shortened the time of experimental verification [14]; Jain et al. used the probabilistic simulation method to calculate the number of active abrasives in the AFM process [19]. Liu et al. based their work on the crystal plasticity finite element model and the coupled Eulerian–Lagrangian method to analyze the roughness evolution of constrained surfaces [20]. However, there is currently no research on the application of multiple linear regression to predict surface roughness. The advantages of multiple linear regression are that the calculation is simple, the output is a linear relationship, the model is stable and intuitive, and the correlation between factors can be accurately evaluated.
The novelty of this study is to analyze the influence of different parameters (such as workpiece shape, abrasive flow state, rotational speed and abrasive quality) on the surface roughness of the AFM-processed workpiece, and determine the priority of the effect of machining parameters on the surface roughness. The optimal parameter combination is obtained by designing an orthogonal test, and a prediction model based on the linear regression equation is developed using the orthogonal test data. Finally, the model is tested and analyzed to verify the reliability of the model.

2. Materials and Methods

2.1. Description of the Test Device

The experiment was conducted by using a self-made gas-solid two-phase flow abrasive pool processing experimental platform, as shown in Figure 1. The processing principle of the self-made abrasive pool processing platform used in this experiment is to mix solid-phase abrasives with gas to obtain gas-solid two-phase flows with various flow states. The abrasive contacts, rubs and collides with the surface of the workpiece, and micro-cuts the surface of the workpiece with the abrasive to realize the finishing of the surface of the workpiece by the gas-solid two-phase flow. Figure 2 is a schematic diagram of the gas-solid two-phase flow abrasive pool finishing process.
Figure 3 shows the cloud diagram of contact force between abrasive particles and workpiece and between abrasive particles simulated by using EDEM software. Figure 3a is the force when the abrasive grains contact the workpiece when no wind force is applied, and the force of the abrasive grains is concentrated on the workpiece. Figure 3b is the force diagram of the abrasive grains when the wind force is applied; it is seen that the force of the abrasive particles on the workpiece gradually decreases. Figure 3c is the force diagram of the abrasive particles after the stable gas-solid two-phase flow is formed. The result shows that the force of the abrasive particles is no longer concentrated around the workpiece, the polishing is more disordered, which can effectively reduce the stress concentration on the surface of the workpiece, prevent the workpiece from rubbing for a long time at a certain point during processing, and generate furrows or pits, so that the processing quality of the workpiece surface is improved.
The workpiece material processed during the experiment was Q235 steel, and the chemical composition of the material is shown in Table 1. The round tubular, square tubular and cylindrical Q235 steel materials were selected for the experiment. The abrasive grains used in the processing of the abrasive pool were white corundum with different particle sizes (as shown in Figure 4). White corundum is suitable for grinding and polishing high carbon steel, high-speed steel and hardened steel, and has excellent grinding and polishing performance. The specific parameters are shown in Table 2.
In the test, five factors, such as workpiece shape, abrasive grain size, gas-solid two-phase flow state, abrasive grain shape, and rotating speed were selected as the main parameters. The parameter ranges are set in Table 3. The minimum air pressure required by the air compressor is 0.735 MPa.
Since there were many factors involved, to reduce the number of experiments without reducing the experimental effect, this paper adopts the orthogonal test to obtain the optimal combination of experiments. The level of each factor was quite different; it was necessary to select the orthogonal experiment method at different levels. Orthogonal experiments at different levels are divided into two methods: the mixed orthogonal experiment method and the quasi-horizontal method. Mixed orthogonal experiments are mostly suitable for situations where there are many levels of a single factor and the levels of other factors are the same, and its expression is as follows:
L n m 1 k 1 m 2 k 2 ,     n = k 1 × m 1 1 + k 2 × m 2 1 + k x × m x 1 + 1
where L is an orthogonal table, n is the number of experiments, k is a factor, and m is the number of levels.
After substituting the data, the number of mixed orthogonal experiments was calculated as 15 times. The quasi-horizontal rule is to complete the factors with fewer levels through repeated levels or other methods to ensure that each factor can meet the same level, and then design a standard orthogonal experiment. After the orthogonal table was quasi-horizontal transformation, generally the table that has been transformed into a quasi-horizontal column is listed as a quasi-horizontal column. According to Table 3, after the shape of the abrasive grains was filled by the pseudo-horizontal method, the standard orthogonal experiment table of four factors and three levels L9 (43) can be obtained. The orthogonal experiment table is shown in Table 4.
The gas-solid two-phase flow state, workpiece shape and abrasive particle shape are all dimensionless factors. To simplify the orthogonal experiment, the above three parameters were digitized, and the gas-solid two-phase flow state can be replaced by the Reynolds number; the value range was set to 2000, 4000, and 6000. The shape of the workpiece can be replaced by the contact surface between the workpiece and the abrasive flow; the contact between the workpiece and the abrasive flow is generally the entire curved surface, the contact of the square tube is four planes, and the contact of the cylinder is the bottom plane. Process these three values into data and use 1/3, 4/3, and 1/2 instead. There were two shapes of abrasives, spherical and irregular, and numbers 1 and 2 can be used to represent spherical and irregular shapes. For testing the overall effect of the abrasive pool processing, it was necessary to ensure the continuity of the experiment. Therefore, the same workpiece was divided into three stages for experimental testing, that is, the orthogonal experiment of the workpiece under 24, 80 and 120 mesh abrasives, respectively. The roughness of the workpiece will be measured in each stage of the experiment, and the experimental workpiece will not be replaced during the period, but the experimental measurement was still required for each stage of processing. After the experiment, the surface roughness of the workpiece was measured with a roughness meter, and the measurement points are shown in Figure 5.

2.2. Analysis of Orthogonal Experiment Results

Firstly, the surface roughness of each marked point of the round tube workpiece under 24-mesh abrasive was analyzed by the intuitive analysis method. Figure 6 shows the average of five measurement roughness values of the same marked points of the round tube under the orthogonal experimental conditions.
To fully analyze the polishing effect of the workpiece, it is necessary to analyze the results of each point. Take the analysis of point 1 as an example. Table 5 shows the average and range of surface roughness changes at point 1. K represents the sum of surface roughness changes measured by different factors at the same level, K1–K3 are rust removal experiments, K4–K6 are rough polishing experiments, and K7–K9 are fine polishing experiments. k represents the average value of different factors at the same level. R represents the extreme difference, R = Kmax − Kmin, the larger the extreme difference, the greater the influence of this factor on finishing. A, B, C and D are used to replace the spindle speed, gas-solid two-phase flow, workpiece shape and abrasive particle shape in the processing factors, respectively. The subscripts 1, 2, and 3 are the level 1, level 2, and level 3 corresponding to the factor, respectively.
It is seen in Table 5 that the range of the D factor is the smallest when the abrasive pool is processed, indicating that the shape of the abrasive grains has the smallest influence on the surface roughness; the range of C is the largest, indicating that the shape of the workpiece has the greatest influence on the surface roughness of the workpiece. With the reduction in the abrasive and the improvement of the polishing quality, the influence of the gas-solid two-phase flow on the workpiece gradually increases, and the range value gradually increases from 0.401 μm at the beginning to 0.9397 μm at the end. With the improvement of the processing quality, the influence of the particle shape of the abrasive on the workpiece gradually increases. From the range, it is seen that the four factors act on the workpiece at the same time, the optimal combination of roughness is always A3B2C1D2, that is, the spindle speed is 1200 rpm, the gas-solid two-phase flow state is turbulent, the shape of the workpiece is a circular tube, and the shape of the abrasive particles is irregular, which shows that the workpiece shape has the greatest influence on the machining effect until the final polishing experiment.
However, it is impossible to quantitatively judge which factor has the greatest impact on the surface processing effect of the workpiece at which level by only using the average value and the range. To more intuitively determine which factor has a comparatively large influence on the roughness, it is necessary to make a factor trend graph for the orthogonal result. Continue to take point 1 as an example, and evaluate it according to the average value of roughness. Under the same factor, which type of surface roughness is the smallest, indicating that this level is the optimal level. The abscissa represents the different levels of different factors, the ordinate represents the surface roughness value during processing, and the final trend graph of the factors is formed, as shown in Figure 7.
It is presented in Figure 7 that the trend of the C factor in the same experimental stage changes greatly, indicating that the shape of the workpiece has the greatest influence on the processing of the abrasive pool. With the gradual precision of the experimental processing and the decrease in the number of abrasive grains, the influence of the gas-solid two-phase flow on the processing effect was more and more obvious. The influence of the abrasive shape on the roughness of the machined workpiece was gradually increasing. The slope of factors B1-B3 gradually increases, indicating that 120-mesh abrasive grains may still be relatively large abrasive grains for abrasive pool processing, and were not the most suitable abrasive grains for processing. It can be speculated that the most suitable abrasive pool processing may be a more precise finishing. It can be seen from Figure 7 that the optimal combination of the rust removal experiment at point 1 is A3B1C1D2, the optimal combination of the rough polishing experiment is A2B2C1D3, and the optimal combination of the fine polishing experiment is A1B2C1D2.
By analogy, the three-stage experiments of the remaining detection points were comprehensively analyzed. Table 6 shows the influence degree of each factor in all points and the optimal parameter combination table [11].
The analysis of the results in Table 6 shows that the more factors and levels in the optimal combination, the greater the impact on the surface machining quality of the workpiece:
As the first influencing factor, A appears 0 times, B is 0 times, C is 15 times, and D is 0 times; as the second influencing factor, A appeared 3 times, B was 11 times, C was 0 times, and D was 1 time; as the third influencing factor, A appeared 10 times, B was 3 times, C was 0 times, and D was 2 times; as the fourth influencing factor, A appears 2 times, B is 1 time, C is 0 times, and D is 12 times.
From the results of all stages, the primary and secondary degrees of each factor affecting the workpiece are workpiece shape, abrasive flow state, rotating speed, and abrasive shape. In the same way, under the same factor, the one with the highest number of levels is the optimal condition, and the results are specified in Table 7.
Therefore, the comprehensive optimal parameter combination for abrasive pool processing is A3B2C1D2.

3. Roughness Prediction Model

The optimal parameter combination within a certain range was obtained through the orthogonal experiment, and a surface roughness model was established on this basis to predict the surface roughness of the workpiece, reduce the number of experiments and save the experiment time. Based on the orthogonal experimental data, the multiple linear regression analysis methods were used to establish the surface roughness prediction model of abrasive pool machining as the basis of the surface roughness prediction of abrasive pool machining.
Prediction requirements for the surface roughness of abrasive pool processing: while the surface quality of the workpiece meets the requirements, the processing efficiency was maximized as much as possible. The processing parameters of the experiment, namely, spindle speed, gas-solid two-phase flow state, workpiece shape, and abrasive shape are used as the input end, and the surface roughness of the workpiece was used as the output end, and the results should be more accurate, which can provide convenience for the subsequent experimental progress. It can provide an important reference for subsequent abrasive pool processing.

3.1. Multiple Linear Regression Prediction Models

Suppose there is a linear correlation between the random variable y and the p independent variables x1, x2, …, xp, the actual sample size is n, and the observed value is xi1, xi2, …, xip, yi (i = 1, 2,…, n). Then, the n observations can be expressed as Equation (2):
y 1 = β 0 + β 1 x 11 + β 2 x 12 + + β p x 1 p + ε 1 y 2 = β 0 + β 1 x 21 + β 2 x 22 + + β p x 2 p + ε 2 y n = β 0 + β 1 x n 1 + β 2 x n 2 + + β p x n p + ε n
Among them, β0, β1, …, βn are position parameters, x1, x2, …, xp are p general variables that can be accurately measured and controlled, and ε1, ε2, …, εp are random errors. It is assumed that εi are independent random variables that obey the same normal distribution N(0, σ). Equation (3) can be represented by a matrix:
y = X β + ε
where
y = y 1 y 2 y n   X = 1 x 11 x 11 x 11 1 x 11 x 11 x 11 1 1 x 11 x 11 x 11   β = β 1 β 2 β n     ε = ε 1 ε 2 ε n
Solving the multiple linear regression equation is to establish the multiple linear regression equation by solving the estimated value b of β:
y ^ = b 0 + b 1 x 1 + b 2 x 2 + + b p
There is a linear relationship between milling and the surface roughness of the workpiece. Through the similarity principle, it can be inferred that the abrasive pool machining also has this relationship with the workpiece. Tipnis et al. established an empirical model of surface roughness for cutting speed, feed and depth of cut. Based on this, scholars established a general model of milling and surface roughness:
R a = c a p b 1 a e b 2 n b 3 f z b 4
where c is the correction coefficient of the workpiece material, ap is the cutting depth, ae is the cutting width, f is the feed amount, and b1, b2, b3, and b4 are to be estimated. By the second law of similarity, for similar physical quantities, the similarity criterion should be the same. After substituting the abrasive pool processing parameters into the above equation, we can get:
R a = c m b 1 R e p m f b 2 n b 3 M b 4
where m represents the shape of the workpiece, Repmf is the Reynolds number of the abrasive in the abrasive pool, that is, the flow state after the gas-solid two-phase flow is formed, M is the shape of the abrasive, and n is the spindle speed. Since the output value is a linear function, it is necessary to linearize the nonlinear function of the above equation, take the logarithm of both sides, and Equation (8) can be obtained:
l g R a = l g c + b 1 l g m + b 2 l g R e p m f + b 3 l g n + b 4 l g M
Suppose that:
l g R a = y , l g c = b 0 , l g m = x 1 , l g R e p m f = x 2 , l g n = x 3 , l g M = x 4
transform it into:
y = b 0 + b 1 x 1 + b 2 x 2 + b 3 x 3 + b 4 x 4
In Equation (10), y is a statistical variable. The regression coefficient b can be obtained using Equation (11):
b = X X 1 X Y
First, solve the model. The least-squares method was used to find b, and the surface roughness prediction model of the rust removal stage was first calculated. To optimize the data, the surface roughness measurements were performed on multiple points of different workpieces in the previous test, each factor level corresponds to five roughness results. The points with the same factor and level were taken as the average surface roughness. It can be obtained by calculation:
X = 2.778 3.301 0.476 0 2.778 3.602 0.125 0.301 2.778 3.778 0.301 0.301 2.954 3.301 0.125 0.301 2.954 3.602 0.301 0 2.954 3.778 0.476 0.301 3.079 3.301 0.301 0.301 3.079 3.602 0.476 0.301 3.079 3.778 0.125 0 Y = 0.6784 0.7118 0.7436 0.7218 0.7023 0.6890 0.6958 0.6627 0.7197
The multiple linear regression equation can be solved by using the regress command of Matlab. The regression coefficient b, residual r, and F value can be obtained by the code:
[b, bint, r, rint, stats] = regress (Y, X).
b =   1.6487 0.0265 0.0208 0.0562 0.0083 T
Then, conduct residual analysis. Randomness and unpredictability are key components of all regression models, and random error needs to be random and unpredictable. After obtaining the regression model, it is necessary to ensure that the residuals are untraceable, usually based on residual plots to analyze the reliability of the data and whether the model is correct.
As can be seen from Figure 8, the residual distribution of the model is close to the normal distribution and belongs to the standardized residual. It shows that the model is not disturbed by specific factors, which ensures the correctness of the model. Therefore, the surface roughness prediction model in the rust removal stage is:
R a = e 1.7071 m 0.0348 R e p m f 0.0230 n 0.0898 M 0.0146
In the same way, the residual analysis of the surface roughness prediction model in the rough polishing stage and the fine polishing stage was conducted. The residual distribution diagram of the rough throwing stage is shown in Figure 9.
As can be seen from Figure 9 the residual distribution in the rough throwing stage is a normal distribution, so the roughness prediction model in the rough throwing stage is:
R a = e 0.0599 m 0.0275 R e p m f 0.0424 n 0.1429 M 0.0074
As can be seen from Figure 10, the residual distribution is a normal distribution, so the roughness prediction model in the fine polishing stage is:
R a = e 1.3465 m 0.0808 R e p m f 0.1019 n 0.6812 M 0.0597
y = β 0 + β 1 x 1 + β 2 x 2 + + β p x p

3.2. Significance Test of Prediction Model

Firstly, use the F-test to test the overall significance level of the model. The significance test for the multiple linear regression equation is to test whether the parameters in Equation (18) are significantly different from 0.
y = b 0 + b 1 x 1 + b 2 x 2 + b 3 x 3 + + b n x n
suppose that:
H 0 : b 1 = b 2 = = b k = 0   F F α H 1 : b j j = 1 , 2 , , k   n o t   a l l   0   F > F α
Since y obeys the normal distribution, the sum of the squares of the samples corresponding to y still obeys the square of the normal distribution, so:
E S S = 1 n y ^ t y ¯ 2 / σ 2 ~ x 2 k
R S S = 1 n y ¯ y ^ t 2 / σ 2 ~ x 2 n k 1
The equation for the F-test is:
F = E S S / k R S S / n k 1
Among them, RSS is derived from the regression sum of squares, and ESS is derived from the residual sum of squares, both of them follow the square of the normal distribution, n represents the amount of data, and k is the number of calculation parameters. Calculated by Matlab, the F value of the rust removal stage is 15.352 greater than 15.22 [F (4,5) (α = 0.01)], that is, the assumption H1 is accepted, which means that bj (j = 1, 2,…, k) is not all 0, indicating that the prediction model of the rust removal stage is significant.
T-test uses t-distribution theory to infer the probability of differences, that is, to test the significance of the independent variables in the equation. Provided the independent variables in each group show a normal distribution, it means that there are consistent differences between groups, that is, the rejection H0 can be maximized. The Equation T-test is:
T = X ¯ μ σ X N
where X ¯ is the mean value of each factor, N is the number of samples, and T needs to follow a normal distribution. The significance test is conducted between the factors. If T < , hypothesis H0 is accepted, indicating that the two data are similar and there is no discrimination; if T > , hypothesis H1 is accepted, indicating that the two dates are from different distribution data, and differentiated. Check the T value distribution table, when α = 0.01, is 2.8214. After substituting the orthogonal experimental data into Equation (23), the T value of each coefficient is obtained as shown in Table 8. It can be seen that the absolute value of T of each factor is greater than , so the assumption H0 is accepted, and there is a clear distinction between independent variables. It shows that each factor in the prediction model has a significant effect on the surface roughness.
Similarly, the F value of the rough throwing stage is 32.156 greater than 15.22 [F (4,5) (α = 0.01)], the absolute value of T is greater than (2.8214), the F value of the fine throwing stage is 50.317 greater than 15.22 [F (4,5) (α = 0.01)], and the absolute value of T is greater than (2.8214).
In summary, the surface roughness models of each polishing stage have good significance, high reliability and data fit.

3.3. Significance Test of the Prediction Model

Abrasive pool processing is most effective in fine polishing. To verify the accuracy of the surface roughness prediction model, 120-mesh abrasive particles were selected to conduct the precision polishing experiment on the cylindrical workpiece. Five combinations were randomly selected for the experiment. The feature points were measured, and the above experimental conditions were substituted into the prediction model. The experimental results and the predicted results are shown in Table 9. Figure 11 illustrated a comparison diagram of the predicted value and the measured value:
The average error ratio calculation equation is:
δ = 1 n 1 n M m M f M m × 100 %
where Mm is the measured value and Mf is the predicted value. It can be seen from Table 10 that the maximum error between the predicted value and the actual value of the surface roughness processed by the abrasive pool is 0.068 μm, the minimum error is 0.0094 μm, the error ratio is basically below 10%, and the average error is only 6.44%, indicating that in the fine polishing experiment stage, the roughness prediction model is basically consistent with the actual value, and the established surface roughness prediction model can predict the surface roughness of the abrasive pool well.
The roughness prediction model of the rust removal and rough polishing experiments is the same as the roughness prediction model of the fine polishing experiment. The results are shown in Table 11. Experiments 1–5 are rust removal experiments, and 6–10 are rough polishing experiments. R1 is the mean value of the rust removal stage, and R2 is the mean value of the rough throwing stage.
From the perspective of error ratio, the prediction model of workpiece surface roughness during rust removal is more accurate, with an error ratio of only 1.4%, followed by the fine polishing experiment with an error ratio of 6.44%, and the last rough polishing experiment with an error ratio of 6.8%. Observing the error value, it can be found that the error value during fine polishing is only 0.00844 μm. From the point of view of the error value, the average error of the surface roughness prediction model during fine polishing is the smallest, and the surface roughness model is more accurate. In general, the errors of the surface roughness prediction models for rust removal, rough polishing and fine polishing are all within an appropriate range, and these three surface roughness models can predict the surface roughness of the actual abrasive pool machining.

4. Conclusions

The influence of the parameters involved in the abrasive pool machining on the workpiece machining quality was evaluated by designing an orthogonal test, and the optimal parameter combination was determined by using the range method and factor trend diagram. The surface roughness prediction model was established in three stages, and the accuracy of the model was verified. The specific conclusions are as follows:
(1)
The results of the orthogonal test show that the abrasive pool machining has the advantages of high machining efficiency and small surface loss, and the optimal machining parameter combination of A3B2C1D2 was obtained.
(2)
The surface roughness prediction model of abrasive pool machining was obtained through multiple linear regression equations, residual analysis and significance analysis were carried out on the model, and MATLAB was used for solution analysis to determine the randomness and the randomness of the surface roughness prediction model.
(3)
The actual surface roughness was obtained through the experiment and compared with the predicted value of the model. As a result, the maximum value of the average error between the predicted value of the surface model and the actual value was 0.339 μm, and the minimum value was 0.008 μm. Predict surface roughness for abrasive pool machining.
In future research, the influence of other process parameters on the grinding and polishing accuracy should be discussed in more detail, and the flow regime of the abrasive pool should be systematically analyzed through hydrodynamic simulation to better predict the surface roughness of the machined workpiece.

Author Contributions

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

Funding

This research received no external funding.

Acknowledgments

This work has been partly supported by National Natural Science Foundation of China (No. 51,805,299 and No. 52075306), Key R&D project of Shandong Province (2019GGX104081, 2019GGX104033), Shandong Province’s Key Support Regions Introducing Urgently Needed Talent Projects, and Young Innovative Talents Introduction and Training Program Project of Shandong Provincial Department of Education.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiao, L.; Wu, Y.; Wang, X.; Guo, H.; Liang, Z. Fundamental performance of Magnetic Compound Fluid (MCF) wheel in the ultra-fine surface finishing of optical glass. Int. J. Mach. Tools Manuf. 2013, 75, 109–118. [Google Scholar] [CrossRef]
  2. Dixit, N.; Sharma, V.; Kumar, P. Research trends in abrasive flow machining: A systematic review. J. Manuf. Process. 2021, 64, 1434–1461. [Google Scholar] [CrossRef]
  3. Kenda, J.; Pusavec, F.; Kermouche, G.; Kopac, J. Surface integrity in abrasive flow machining of hardened tool steel AISI D2. Procedia Eng. 2011, 19, 172–177. [Google Scholar] [CrossRef]
  4. Li, J.Y.; Liu, W.N.; Yang, L.F.; Liu, B.; Zhao, L.; Li, Z. The development of nozzle micro-hole abrasive flow machining equipment. Appl. Mech. Mater. 2011, 44, 251–255. [Google Scholar] [CrossRef]
  5. Wang, A.C.; Chen, K.Y.; Cheng, K.C.; Chiu, H. Elucidating the Effects of Helical Passageways in Abrasive Flow Machining. Adv. Mater. Res. 2011, 264, 1862–1867. [Google Scholar] [CrossRef]
  6. Singh, S.; Shan, H.S.; Kumar, P. Experimental studies on mechanism of material removal in abrasive flow machining process. Mater. Manuf. Process. 2008, 23, 714–718. [Google Scholar] [CrossRef]
  7. Xie, W.B.; Zhang, K.H.; Zhang, S.W.; Xu, B. Research on Abrasive Flow Machining for the outer rotor of cycloidal pump. Key Eng. Mater. 2013, 546, 50–54. [Google Scholar] [CrossRef]
  8. Wang, A.C.; Weng, S.H. Developing the polymer abrasive gels in AFM processs. Mater. Process. Technol. 2007, 192–193, 486–490. [Google Scholar] [CrossRef]
  9. Barletta, M.; Rubino, G.; Valentini, P.P. Experimental investigation and modeling of fluidized bed assisted drag finishing according to the theory of localization of plastic deformation and energy absorption. Int. J. Adv. Manuf. Technol. 2015, 77, 2165–2180. [Google Scholar] [CrossRef]
  10. Loc, P.H.; Shiou, F.J.; Yu, Z.R.; Hsu, W.Y. Investigation of Optimal Air-Driving Fluid Jet Polishing Parameters for the Surface Finish of N-BK7 Optical Glass. J. Manuf. Sci. Eng. 2013, 135, 011015. [Google Scholar] [CrossRef]
  11. Boga, C.; Koroglu, T. Proper estimation of surface roughness using hybrid intelligence based on artificial neural network and genetic algorithm. J. Manuf. Process. 2021, 70, 560–569. [Google Scholar] [CrossRef]
  12. Ficko, M.; Begic-Hajdarevic, D.; Husic, M.C.; Berus, L.; Cekic, A.; Klancnik, S. Prediction of Surface Roughness of an Abrasive Water Jet Cut Using an Artificial Neural Network. Materials 2021, 14, 3108. [Google Scholar] [CrossRef] [PubMed]
  13. Spaic, O.; Krivokapic, Z.; Kramar, D. Development of family of artificial neural networks for the prediction of cutting tool condition. Adv. Prod. Eng. Manag. 2020, 15, 164–178. [Google Scholar] [CrossRef]
  14. Karabulut, S. Optimization of surface roughness and cutting force during AA7039/Al2O3 metal matrix composites milling using neural networks and Taguchi method. Measurement 2015, 66, 139–149. [Google Scholar] [CrossRef]
  15. Alajmi, M.S.; Almeshal, A.M. Prediction and Optimization of Surface Roughness in a Turning Process Using the ANFIS-QPSO Method. Materials 2020, 13, 2986. [Google Scholar] [CrossRef] [PubMed]
  16. Hribersek, M.; Berus, L.; Pusavec, F.; Klancnik, S. Empirical Modeling of Liquefied Nitrogen Cooling Impact during Machining Inconel 718. Appl. Sci. 2020, 10, 3603. [Google Scholar] [CrossRef]
  17. Radovanovic, M. Multi-Objective Optimization of Abrasive Water Jet Cutting Using MOGA. Procedia Manuf. 2020, 47, 781–787. [Google Scholar] [CrossRef]
  18. Liu, D.; Huang, C.Z.; Wang, J.; Zhu, H.T.; Yao, P.; Liu, Z.W. Modeling and optimization of operating parameters for abrasive waterjet turning alumina ceramics using response surface methodology combined with Box-Behnken design. Ceram. Int. 2014, 40, 7899–7908. [Google Scholar] [CrossRef]
  19. Jain, R.K.; Jain, V.K. Stochastic simulation of active grain density in abrasive flow machining. J. Mater. Process. Technol. 2004, 152, 17–22. [Google Scholar] [CrossRef]
  20. Liu, C.S.; Xu, W.J.; Niu, T.H.; Chen, Y.T. Roughness evolution of constrained surface based on crystal plasticity finite element model and coupled Eulerian-Lagrangian method. Comp. Mater. Sci. 2022, 201, 110900. [Google Scholar] [CrossRef]
Figure 1. Three-dimensional diagram of abrasive pool test bench.
Figure 1. Three-dimensional diagram of abrasive pool test bench.
Micromachines 13 01649 g001
Figure 2. The principle of gas-solid two-phase flow finishing.
Figure 2. The principle of gas-solid two-phase flow finishing.
Micromachines 13 01649 g002
Figure 3. Abrasive contact with workpiece diagram. (a); not applied wind force (b) applied wind power; (c) stable gas-solid two-phase flow is formed.
Figure 3. Abrasive contact with workpiece diagram. (a); not applied wind force (b) applied wind power; (c) stable gas-solid two-phase flow is formed.
Micromachines 13 01649 g003
Figure 4. Abrasive.
Figure 4. Abrasive.
Micromachines 13 01649 g004
Figure 5. Workpiece diagram: (a) round tube; (b) square tube; (c) cylinder.
Figure 5. Workpiece diagram: (a) round tube; (b) square tube; (c) cylinder.
Micromachines 13 01649 g005
Figure 6. The roughness height of circular tube marking point by orthogonal experiment: (a) 120 mesh abrasive; (b) 80 mesh abrasive; (c) 24 mesh abrasive.
Figure 6. The roughness height of circular tube marking point by orthogonal experiment: (a) 120 mesh abrasive; (b) 80 mesh abrasive; (c) 24 mesh abrasive.
Micromachines 13 01649 g006
Figure 7. Factor trend of point 1 by test.
Figure 7. Factor trend of point 1 by test.
Micromachines 13 01649 g007
Figure 8. Distribution of surface roughness residuals in the rust removal stage.
Figure 8. Distribution of surface roughness residuals in the rust removal stage.
Micromachines 13 01649 g008
Figure 9. Distribution of surface roughness residuals in the rough polishing stage.
Figure 9. Distribution of surface roughness residuals in the rough polishing stage.
Micromachines 13 01649 g009
Figure 10. The residual distribution of the fine throwing stage.
Figure 10. The residual distribution of the fine throwing stage.
Micromachines 13 01649 g010
Figure 11. The measured and predicted values of the finishing polish.
Figure 11. The measured and predicted values of the finishing polish.
Micromachines 13 01649 g011
Table 1. Workpiece material properties.
Table 1. Workpiece material properties.
MaterialYield Strength
/MPa
Hardness
/HB
Melting Point/°CC/%Mn/%Si/%S/%
Q23523516516000.12–0.200.3–0.650.30.05
Table 2. Physical properties of white corundum.
Table 2. Physical properties of white corundum.
MaterialHardnessShapeDensity (Bulk) g/cm3
Al2O39.0spherical/irregular1.53–1.99
Table 3. Orthogonal experiment parameter range setting.
Table 3. Orthogonal experiment parameter range setting.
Processing Parameters
Workpiece ShapeGas-solid Two-Phase FlowAbrasive ShapeSpindle Speed/rpm
1round tubebulk fluidized bedspherical600
2square tubeturbulent fluidized bedirregular900
3cylinderspouted fluidized bed/1200
Table 4. Orthogonal experiment table.
Table 4. Orthogonal experiment table.
Experiment NumberFactors
Spindle Speed/rpmGas-Solid Two-Phase FlowWorkpiece ShapeAbrasive Shape
160020001/31
260040004/32
360060001/22
490020004/32
590040001/21
690060001/32
7120020001/22
8120040001/32
9120060004/31
Table 5. Range analysis results at point 1.
Table 5. Range analysis results at point 1.
Spindle Speed (A)Gas-Solid Two-Phase Flow (B)Workpiece Shape (C)Abrasive Shape(D)
K115.33615.03614.28815.105
K215.24815.15616.25115.419
K314.97915.43715.0915.105
K45.07545.03774.63185.4249
K54.98964.94915.98255.1204
K65.49385.5724.94455.0135
K72.06462.20131.49462.3848
K82.29041.77293.4822.0438
K92.3322.71281.71042.2584
k15.1225.0124.7635.035
k25.0955.0525.4175.140
k34.9335.1465.035.035
k41.69181.67921.54391.8083
k51.66321.64971.99411.7068
k61.83121.85731.64811.6711
k70.68820.73370.49820.7949
k80.76340.59091.16060.6812
k90.77730.90420.57010.7528
R10.3870.4011.9630.314
R20.41840.62291.35070.3045
R30.26720.93971.98720.3408
The biggest influencing factor in the rust removal stageCBAD
The biggest influencing factor in rough throwing stageCBAD
The biggest influencing factor in the polishing stageCBDA
Rust removal experimentRough throwing experimentFine polishing experiment/
The best combination of processing qualityA3B2C1D2A3B2C1D2A3B2C1D2/
Table 6. The influence degree of each factor and the optimal parameter combination table of each experiment.
Table 6. The influence degree of each factor and the optimal parameter combination table of each experiment.
Marked Point12345
Rust removal experimentInfluence levelCBADCABDCABDCBADCBAD
Optimal
combination
A3B1C1D2A3B2C1D2A1B1C1D1A1B2C2D2A3B2C1D2
Rough
throwing
experiment
Influence levelCBADCDBACADBCBADCBAD
Optimal
combination
A2B2C1D3A1B1C1D1A3B2C1D2A1B1C1D1A1B1C1D1
Fine polishing experimentInfluence levelCBADCBADCBDACBADCBAD
Optimal
combination
A1B2C1D2A3B2C1D2A2B2C2D1A3B2C1D2A3B2C1D2
Table 7. Optimal level.
Table 7. Optimal level.
FactorsLevel 1Level 2Level 3
A519
B4110
C1320
D5100
Table 8. Distribution of T value in each stage of abrasive pool processing.
Table 8. Distribution of T value in each stage of abrasive pool processing.
Regression Coefficientsb1b2b3b4
Derusting stage T value3.424.283.466.84
Rough throwing stage T value5.68−3.519.67−4.85
Table 9. Parameter settings of surface roughness prediction experiment.
Table 9. Parameter settings of surface roughness prediction experiment.
Experiment
Number
Experimental Parameters
ABCD
160060001/22
290040001/21
390060001/32
4120020001/22
5120060004/31
Table 10. The measured and predicted values of the finishing polish.
Table 10. The measured and predicted values of the finishing polish.
ExperimentMeasured ValuePredictive ValueDifferenceError Ratio
10.70130.63330.0689%
20.59120.6545−0.063310%
30.55250.49610.066411%
40.58950.5989−0.00941.5%
51.35211.3616−0.00950.7%
R0.757320.748880.008446.44%
Table 11. The measured and predicted values of rust removal and rough polishing.
Table 11. The measured and predicted values of rust removal and rough polishing.
ExperimentMeasured ValuePredictive ValueDifferenceError Ratio
15.1355.0630.0721.4%
25.0364.9670.0691.3%
34.874.80.071.4%
44.9554.8330.1222.4%
55.3655.2560.1092%
61.6451.6580.0130.7%
71.9871.6480.33917%
81.6121.5740.0382.3%
91.6931.6050.0882.1%
102.1151.9440.1718.0%
R15.07364.9978−0.07581.4%
R21.81041.6858−0.12466.8%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, W.; Yuan, W.; Yu, J.; Guo, Q.; Chen, S.; Yang, X.; Cong, J. Prediction of Surface Roughness in Gas-Solid Two-Phase Abrasive Flow Machining Based on Multivariate Linear Equation. Micromachines 2022, 13, 1649. https://doi.org/10.3390/mi13101649

AMA Style

Wang W, Yuan W, Yu J, Guo Q, Chen S, Yang X, Cong J. Prediction of Surface Roughness in Gas-Solid Two-Phase Abrasive Flow Machining Based on Multivariate Linear Equation. Micromachines. 2022; 13(10):1649. https://doi.org/10.3390/mi13101649

Chicago/Turabian Style

Wang, Wenhua, Wei Yuan, Jie Yu, Qianjian Guo, Shutong Chen, Xianhai Yang, and Jianchen Cong. 2022. "Prediction of Surface Roughness in Gas-Solid Two-Phase Abrasive Flow Machining Based on Multivariate Linear Equation" Micromachines 13, no. 10: 1649. https://doi.org/10.3390/mi13101649

APA Style

Wang, W., Yuan, W., Yu, J., Guo, Q., Chen, S., Yang, X., & Cong, J. (2022). Prediction of Surface Roughness in Gas-Solid Two-Phase Abrasive Flow Machining Based on Multivariate Linear Equation. Micromachines, 13(10), 1649. https://doi.org/10.3390/mi13101649

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