{'fraction_non_alphanumeric': 0.030696827083649614, 'fraction_numerical': 0.012691665401548504, 'mean_word_length': 4.547582954354051, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Analysis and Optimization of the Machining Characteristics of High-Volume Content SiCp/Al Composite in Wire Electrical Discharge Machining

: With the properties of high speciﬁc strength, small thermal expansion and good abrasive resistance, the particle-reinforced aluminum matrix composite is widely used in the ﬁelds of aerospace, automobile and electronic communications, etc. However, the cutting performance of the particle-reinforced aluminum matrix composite is very poor due to severe tool wear and low machining efﬁciency. Wire electrical discharge machining has been proven to be a good machining method for conductive material with any hardness. Even so, the high-volume SiCp/Al content composite is still a difﬁcult-to-machine material in wire electrical discharge machining due to the inﬂuence of insulative the SiC particle. The goal of this paper is to analyze the machining characteristics and ﬁnd the optimal process parameters for the high-volume content (65 vol.%) SiCp/Al composite in wire electrical discharge machining. Experimental results show that the material removal method of the SiCp/Al composite includes sublimating, decomposing and particle shedding. The material removal rate is found to increase with the increasing pulse-on time, ﬁrst increasing and then decreasing with the increasing pulse-off time, servo voltage, wire feed and wire tension. Pulse-on time and servo voltage are the dominant factors for surface roughness. In addition, the multi-objective optimization method of the nondominated neighbor immune algorithm is presented to optimize the process parameters for a fast material removal rate and low surface roughness. The optimized process parameters can increase the material removal rate by 34% and reduce the surface roughness by 6%. Furthermore, the effectiveness of the Pareto optimal solution is proven by the veriﬁed experiment.


Introduction
The particle-reinforced aluminum matrix composite is a material that is prepared by adding reinforcement to the aluminum matrix, such as carbide, nitride or graphite. Compared with the aluminum matrix, the particle-reinforced aluminum matrix composite has better physical and chemical properties, such as low density, high specific strength, excellent high-temperature properties, high wear resistance and excellent stability dimensional [1][2][3]. The SiCp/Al composite is one of the most common particle-reinforced aluminum matrix composites, which is widely used in the fields of aerospace, automobiles and electronic communications, etc. Due to the non-uniform distribution of super-hard SiC particles, SiCp/Al is a difficult-to-machine material in the traditional cutting method. The major displays of machining difficulties are severe tool wear, low machining efficiency and surface defects [4]. With an increasing volume content of the SiC particle, the machining process becomes more and more difficult. This fact severely limits the application and extension of the particle-reinforced aluminum matrix composite.
In wire electrical discharge machining (EDM/WEDM), a good deal of pulse sparks occurs between the electrode and the workpiece. Every pulse spark can produce a small discharge crater (diameter of 1-100 µm) due to melting or vaporizing from high-density thermal energy (1-10 × 10 6 J/m 2 ) [5,6]. Then, continuous pulse sparks can cause considerable material removal efficiency. Because the maximum instantaneous temperature between the electrode and the workpiece can be up to 10,000 • C, EDM/WEDM can process various conductive materials regardless of hardness [7][8][9]. Hence, EDM/WEDM is an alternative method for SiCp/Al composites.
In recent years, much research has been carried out to investigate the machining characteristics of the particle-reinforced aluminum matrix composite in EDM/WEDM. Balasubramaniam V. et al. [10], Gu L. [11], Dey A. [12], Daneshmand S. [13] and Shelvaraj S.G. et al. [14] analyzed the effects of process parameters on the material removal rate (MRR), surface roughness (SR) and electrode wear rate (TWR) in EDM of aluminum matrix composites with a particle content of 7.5-20%. It was pointed out that the method of process parameter optimization could evidently improve the machining characteristics. Singh B [15] compared TWR in traditional EDM and powder-mixed EDM of aluminum matrix composites with a particle content of 10%. It was found that tungsten powder could effectively reduce TWR. Pramanik A [16] studied the effect of wire tension and discharge current on the MRR and surface quality in WEDM of aluminum matrix composites with a particle content of 10%. It was pointed out that surface roughness first decreased and then increased with an increasing discharge current. Kumar N.M. et al. [17] investigated the influence of particle content on the machining characteristics of aluminum matrix composites with particle contents of 0-8% in EDM. It was pointed out that the performance of EDM decreased with an increasing particle content. Bains P.S. [18] employed the magnetic field method to improve surface properties of aluminum matrix composites with a particle content of 37-50% in EDM. It was pointed out that this method could significantly reduce surface microhardness and the thickness of the recast layer. Kumar T.T.S. et al. [19] adopted a response surface methodology to determine the optimal process parameters for aluminum matrix composites with a particle content of 20% in WEDM. Uthayakumar M. [20] analyzed the effects of process parameters on the machining speed and surface roughness of aluminum hybrid composites with a particle content of 20% in EDM. Besides, the gray relational analysis method was adopted to obtain the optimal process parameters for aluminum hybrid composites. Senthilkumar T.S. [21] investigated the effect of particle content on the surface topography in EDM of aluminum hybrid composites with a particle content of 5-8%. It was found that, with an increase in particle content, MRR decreased, while the surface hardness and the diameter of the craters increased. Paswan K. et al. [22] utilized steam as a dielectric medium for machining metal matrix composites with a particle content of 10% in EDM. Compared with the traditional kerosene medium, steam could significantly improve machining efficiency, surface quality and economic benefit. Devi M.B. et al. [23] completed an experimental study to determine the optimal process parameters for aluminum hybrid composites with a particle content of 6% in EDM. It was pointed out that the optimal process parameters for aluminum hybrid composites changed with the content of the reinforced particle.
From abovementioned research, we can find that EDM/WEDM has been proven to be a good machining method for particle-reinforced aluminum matrix composites. Besides, the process parameters are key factors for the machining characteristics of particlereinforced aluminum matrix composites. However, the particle contents of aluminum matrix composites in the abovementioned research are relatively low. As pointed out in reference [17], with an increasing of particle content, the aluminum matrix composite becomes more and more difficult to machine. The optimal process parameters for aluminum matrix composites with different particle contents are also different.
The research object of this paper is the aluminum matrix composite with a highvolume content of reinforced particles (65 vol.% SiCp/Al composite). A set of discharge cutting experiments is carried out to investigate the effects of process parameters on the MRR and SR of the SiCp/Al composite. The machining mechanism of the SiCp/Al composite is revealed through a scanning electron microscope (SEM). In addition, the multi-objective optimization method of the nondominated neighbor immune algorithm (NNIA) is presented to optimize the process parameters for fast MRR and low SR. The feasibility and precision of the optimal process parameters are evaluated by a verified experiment.

Material
The mechanical and physical properties of the 65 vol.% SiCp/Al composite are excellent, such as high thermal conductivity, high specific strength and good abrasive resistance. The specific stiffness of the 65 vol.% SiCp/Al composite is three times higher than the aluminum matrix and 25 times higher than copper. This material is praised as a thirdgeneration electronic packaging material, which is widely used in civil electronic equipment, IGBT plate substrates and wireless base stations. The 65 vol.% SiCp/Al composite in this study is from Xi'An Fadi Technology Co., Ltd. (Xi'An, China) The material properties of the 65 vol.% SiCp/Al composite are listed in Table 1. In addition, the material properties of the SiCp/Al composite can be obtained according to the theory in Ref. [24]. The geometric dimension of the SiCp/Al composite basal plate is 150 mm × 50 mm × 4 mm.

Machine Tools
All discharge cutting experiments are carried out on a wire EDM machine (ACCUTX EZ-43 SA) from ACCUTEX technologies Co., Ltd. (Taiwan, China), as shown in Figure 1. It mainly consists of a workbench, a motion platform, a wire-moving system, a cooling system, CNC and a high-frequency pulse electrical source (the peak discharge voltage is 80 V). The dielectric is deionized water. The wire electrode is copper wire with a diameter of 0.25 mm. The workpiece is completely submerged in deionized water during the discharge process.

The Experiment Design
A set of discharge cutting experiments is implemented to investigate the effects of process parameters on machining characteristics of the 65 vol.% SiCp/Al composite. Consequently, five important process parameters are selected as input factors, which include the pulse-on time (Ton), pulse-off time (Toff), servo voltage (SV), wire feed (WF) and wire

The Experiment Design
A set of discharge cutting experiments is implemented to investigate the effects of process parameters on machining characteristics of the 65 vol.% SiCp/Al composite. Consequently, five important process parameters are selected as input factors, which include the pulse-on time (T on ), pulse-off time (T off ), servo voltage (SV), wire feed (WF) and wire tension (WT). Each process parameter has five levels, as shown in Table 2. Besides this, the material removal rate (MRR) and surface roughness (SR) are chosen as output factors. The calculating formula of MRR is as shown in Equation (1). The arithmetical mean deviation of the profile (Ra) is selected to represent the surface roughness (SR), which is measured by an optical profilometer (WYKO NT9100). The mean value of three measured data is treated as the final value of Ra. The design of discharge cutting experiments is as shown in Table 3. In this study, the discharge current is a constant value, because the machining efficiency is too slow if the discharge current is lower than 10 A, and frequent wire breakages will happen if the discharge current is higher than 10 A.
Here, H is the thickness of the workpiece (mm), L is the cutting length (mm) and t is the cutting time recorded by a stopwatch. The cutting length is set as 10 mm.
The machined surface of the SiCp/Al composite is characterized by a scanning electron microscope (SEM, MIRA 3 LMU) under an acceleration voltage of 20.0 kV and magnification of 1000×.

Experiment Result
On the basis of the design of discharge cutting experiments, the results of discharge cutting experiments can be obtained, as shown in Table 4. The relative deviation of surface roughness is about 0.1-0.5 µm due to the instrumental error and different measure position.  Figure 2 shows the machined surface of the SiCp/Al composite characterized by SEM. Combined with the XRD results of our previous research [25], a large quantity of microspheres is found on the machined surface. This is because the aluminum matrix can be sublimated under an ultrahigh temperature field (up to 10,000 • C) [8,[26][27][28][29] due to discharge sparks. The sublimated aluminum matrix can become solid due to the cooling effect of the dielectric during the pulse-off time. Then, this material can adhere to the machined surface again in the form of a sphere. Besides, many microspheres can accumulate and form a blocky solid metal with many tentacles. This solid metal is called the recast layer. In addition, a large number of micropores are found on the machined surface. A part of the micropores is produced as a result of the gas entering the sublimated aluminum matrix during the recrystallizing process [30]. The other part of the micropores is produced in the preparation process of the SiCp/Al composite. Moreover, microcracks are found on the machined surface. This is the result of the non-uniform temperature field and rapid cooling [31]. Furthermore, many SiC particles and SiC shedding pits are found on the machined surface. As we know, the decomposition point of the SiC particle is higher than the boiling point of the aluminum matrix. Then, it is difficult to remove the SiC particle. When a part of the aluminum matrix around the SiC particle is sublimated, this SiC particle will be exposed. When the aluminum matrix around the SiC particle is completely sublimated, this SiC particle will be shed. Then, the shedding pits will be formed. This is consistent with the perspective in Ref. [32]. In addition, the influence of the direct sublimating of the aluminum matrix from a solid to a gas (thermal dissociation), the thermochemical interaction between ions and the deposition of more-complex secondary compounds of the second order may also contribute to the method of removal of SiCp/Al in EDM/WEDM [33].   Table 4. Figure 3 shows the results of EDS measurement on the machined su Table 4. Table 5 shows the element composition on the machined surface 4. In region A, the contents of C, Si and O elements are obviously high other elements. Besides, in region B, the contents of C, Al and O elemen higher than those of other elements. Hence, it can be inferred the main m A and region B are SiC particle and Al substrate, respectively. Moreover O element means the redox reaction occurs during the machining process on the machined surface is transferred from the wire electrode due to the  Table 4.   Table 4. Table 5 shows the element composition on the machined surface of No.4 in Table 4. In region A, the contents of C, Si and O elements are obviously higher than those of other elements. Besides, in region B, the contents of C, Al and O elements are obviously higher than those of other elements. Hence, it can be inferred the main material in region A and region B are SiC particle and Al substrate, respectively. Moreover, the existence of O element means the redox reaction occurs during the machining process. The Cu element on the machined surface is transferred from the wire electrode due to the violent collision between electron and ion. In addition, in region A, the volume content of the Si element is significantly lower than that of the C element. This may have resulted from the thermal decomposition of the SiC particle. Figure 3 shows the results of EDS measurement on the machined surface of No.4 in Table 4. Table 5 shows the element composition on the machined surface of No.4 in Table  4. In region A, the contents of C, Si and O elements are obviously higher than those of other elements. Besides, in region B, the contents of C, Al and O elements are obviously higher than those of other elements. Hence, it can be inferred the main material in region A and region B are SiC particle and Al substrate, respectively. Moreover, the existence of O element means the redox reaction occurs during the machining process. The Cu element on the machined surface is transferred from the wire electrode due to the violent collision between electron and ion. In addition, in region A, the volume content of the Si element is significantly lower than that of the C element. This may have resulted from the thermal decomposition of the SiC particle.

The Effects of Process Parameters on MRR and SR
According to Table 4, the effects of process parameters on MRR and SR can be acquired, as shown in Figures 4-8. The degree of influence for MRR from high to low in order is T on , SV, Toff, WS and WF. T on and SV are the dominant factors for surface roughness (SR). The other three process parameters have a small effect on SR. Figure 4 shows the effect of T on on MRR and SR. MRR and SR are found to increase with the increase in T on . The growth rate of MRR decreases with the increase in T on . This is because a longer T on can produce larger discharge energy in the single-pulse discharge process. Then, more material can be sublimated and decomposed, which can result in a fast machining speed. A larger discharge crater can be formed, which can lead to a rougher workpiece surface. In addition, the discharge debris between the wire electrode and the workpiece will be greater and greater corresponding to T on . The probability of an arc discharge or short circuit will increase alongside T on , which is harmful to the material removal. Hence, with the increasing of T on , the growth rate of MRR becomes slower and slower.
process. Then, more material can be sublimated and decomposed, which can result in a fast machining speed. A larger discharge crater can be formed, which can lead to a rougher workpiece surface. In addition, the discharge debris between the wire electrode and the workpiece will be greater and greater corresponding to Ton. The probability of an arc discharge or short circuit will increase alongside Ton, which is harmful to the material removal. Hence, with the increasing of Ton, the growth rate of MRR becomes slower and slower.  Figure 5 shows the effect of Toff on MRR and SR. MRR is found to first increase and then decrease with an increasing Toff. On the one hand, an increasing pulse-off time means there is more time to flush the sublimated and decomposed material, which is beneficial to the material removal. On the other hand, a longer Toff can lead to a smaller discharge energy produced in the continuous discharge process. Then, less material can be sublimated and decomposed, which can result in a slow machining speed. In addition, Toff has a small effect on SR. This is because shedding is one form of material removal of the SiCp/Al composite. The shedding pit is a key factor affecting SR. The geometric dimensioning of the shedding pit is decided by the size of the SiC particle.  Figure 5 shows the effect of T off on MRR and SR. MRR is found to first increase and then decrease with an increasing T off . On the one hand, an increasing pulse-off time means there is more time to flush the sublimated and decomposed material, which is beneficial to the material removal. On the other hand, a longer T off can lead to a smaller discharge energy produced in the continuous discharge process. Then, less material can be sublimated and decomposed, which can result in a slow machining speed. In addition, T off has a small effect on SR. This is because shedding is one form of material removal of the SiCp/Al composite. The shedding pit is a key factor affecting SR. The geometric dimensioning of the shedding pit is decided by the size of the SiC particle.  Figure 6 shows the effect of SV on MRR and SR. MRR is found to first increase and then decrease with an increasing SV. This is because increasing SV means increasing the discharge energy for removing the material, which is beneficial to material removal. It will result in the wire frequently drawing back if SV exceeds the critical value, which is harmful to material removal [26]. SR is found to first increase and then decrease with an increasing SV. This is because, on the one hand, an increasing SV can increase the discharge energy in the single-pulse discharge process, which can result in a rough workpiece surface. On the one hand, increasing SV can increase the discharge gap between the wire electrode and the workpiece. Then, more discharge debris can be expelled, which can lead to a smoother workpiece surface.  Figure 6 shows the effect of SV on MRR and SR. MRR is found to first increase and then decrease with an increasing SV. This is because increasing SV means increasing the discharge energy for removing the material, which is beneficial to material removal. It will result in the wire frequently drawing back if SV exceeds the critical value, which is harmful to material removal [26]. SR is found to first increase and then decrease with an increasing SV. This is because, on the one hand, an increasing SV can increase the discharge energy in the single-pulse discharge process, which can result in a rough workpiece surface. On the one hand, increasing SV can increase the discharge gap between the wire electrode and the workpiece. Then, more discharge debris can be expelled, which can lead to a smoother workpiece surface. ful to material removal [26]. SR is found to first increase and then decrease with an in creasing SV. This is because, on the one hand, an increasing SV can increase the discharg energy in the single-pulse discharge process, which can result in a rough workpiece sur face. On the one hand, increasing SV can increase the discharge gap between the wir electrode and the workpiece. Then, more discharge debris can be expelled, which can lea to a smoother workpiece surface.  Figure 7 shows the effect of WF on MRR and SR. MRR is found to first increase an then decrease with an increasing WF. When WF is relatively low, increasing WF can en hance the flow of the dielectric, which is beneficial to the discharge debris being expelled When WF exceeds the critical value, increasing the WF can result in obvious wire vibra tion, which is harmful to the stability of the discharge process. In addition, WF does no have a significant effect on SR.  Figure 7 shows the effect of WF on MRR and SR. MRR is found to first increase and then decrease with an increasing WF. When WF is relatively low, increasing WF can enhance the flow of the dielectric, which is beneficial to the discharge debris being expelled. When WF exceeds the critical value, increasing the WF can result in obvious wire vibration, which is harmful to the stability of the discharge process. In addition, WF does not have a significant effect on SR.  Figure 8 shows the effect of WT on MRR and SR. MRR is found to first increase an then decrease with an increasing WT. When WT is relatively low, increasing WF can re duce the deflection of the wire electrode, which is beneficial to the discharge debris bein expelled. When WT exceeds the critical value, increasing WF can result in wire electrod plastic deformation so as to enhance the wire vibration, which is harmful to the stabilit of the discharge process. In addition, WT does not have a significant effect on SR.   Figure 8 shows the effect of WT on MRR and SR. MRR is found to first increase and then decrease with an increasing WT. When WT is relatively low, increasing WF can reduce the deflection of the wire electrode, which is beneficial to the discharge debris being expelled. When WT exceeds the critical value, increasing WF can result in wire electrode plastic deformation so as to enhance the wire vibration, which is harmful to the stability of the discharge process. In addition, WT does not have a significant effect on SR. Figure 8 shows the effect of WT on MRR and SR. MRR is found to first incr then decrease with an increasing WT. When WT is relatively low, increasing WF duce the deflection of the wire electrode, which is beneficial to the discharge deb expelled. When WT exceeds the critical value, increasing WF can result in wire e plastic deformation so as to enhance the wire vibration, which is harmful to the of the discharge process. In addition, WT does not have a significant effect on SR

The Numerical Relationship between Process Parameters on MRR/SR
Based on the experimental data in Table 4, the numerical relationship betw cess parameters on MRR/SR can be obtained through the method of nonlinear re fitting, as shown in Equations (2) and (3). The numerical analysis software of Min used to obtain the nonlinear regression fitting equation. The nonlinear regressi rithm is Gauss-Newton regression, whereby the maximum number of iteration and the convergence tolerance is 0.00001. Figures 9 and 10 show the residual MRR and SR, respectively. The fitting residuals of MRR and SR essentially obey a distribution. In addition, Table 6 shows the comparative results of experimental d fitting data. The relative errors between experimental data and fitting data are l ±8%. The obtained nonlinear regression fitting equations of MRR and SR can be optimize the multi-objective process parameters.

The Numerical Relationship between Process Parameters on MRR/SR
Based on the experimental data in Table 4, the numerical relationship between process parameters on MRR/SR can be obtained through the method of nonlinear regression fitting, as shown in Equations (2) and (3). The numerical analysis software of Minitab was used to obtain the nonlinear regression fitting equation. The nonlinear regression algorithm is Gauss-Newton regression, whereby the maximum number of iterations is 200 and the convergence tolerance is 0.00001. Figures 9 and 10 show the residual plots for MRR and SR, respectively. The fitting residuals of MRR and SR essentially obey a normal distribution.
In addition, Table 6 shows the comparative results of experimental data and fitting data. The relative errors between experimental data and fitting data are less than ±8%. The obtained nonlinear regression fitting equations of MRR and SR can be used to optimize the multi-objective process parameters.
where T 0 on is 1 ns, T 0 off is 1 µs, SV 0 is 1 V, WF 0 is mm/s, WT 0 is 1 N, MRR 0 is 1 mm 2 /s and SR 0 is 1 µm. The units of MRR and Ra are mm 2 /s and µm, respectively. where T 0 on is 1 ns, T 0 off is 1 μs, SV 0 is 1 V, WF 0 is mm/s, WT 0 is 1 N, MRR 0 is 1 mm 2 /s and S is 1 μm. The units of MRR and Ra are mm 2 /s and μm, respectively.

NNIA
As pointed out in Section 3.3, the effect degree and impact trend of process parameters on MRR and SR are different. In the practical machining process, it is desired that the workpiece is quickly removed with low surface roughness. Hence, the method of multiobjective process parameter optimization is suitable for the above issue.
In this study, the multi-objective optimization method of the nondominated neighbor immune algorithm is presented to optimize the process parameters for fast MRR and low SR. NNIA is a multi-objective optimization algorithm, which simulates the natural immune function. This algorithm is inspired by immunology, which simulates the phenomena of the commensalism of various antibodies and the activation of a small number of antibodies during the immunologic process. This small number of relatively independent nondominated individuals is treated as active antibodies. According to the degree of crowdedness, the active antibodies can clone, recombine and hyper mutate through the selection of a nondominated domain. NNIA has an obvious advantage in the high-dimensional multiobjective optimization problem because it pays more attention to the region with a low degree of crowdedness. Besides, NNIA is a multi-objective optimization algorithm on the basis of the Pareto optimal solution. Figure 11 shows the flow chart of NNIA, and the main procedures of optimization are as follows: (1) Initialization The primary antibody group (B 0 ), dominated antibody group, activity antibody group and clone antibody group are generated in this procedure, where the size of the primary antibody group is n D .
(2) Update dominant groups The dominant antibodies (B t ) are recognized in this procedure. All dominant antibodies are copied to form the temporary dominant antibody group (DT t+1 ).

(3) Select based on nondominated neighbor
If DT t+1 is not more than n D , DT t+1 is set as D t+1 . Otherwise, the crowding distance between all individuals in the DT t+1 is calculated to arrange individuals in descending order. The top-n D individuals in the first group form D t+1 according to the crowding distance in descending order. If D t is not more than n A , A t is set as D t . Otherwise, the top-n D individuals in the first group form A t according to the crowding distance in descending order.
(4) Proportional clone Clone group (C t ) is obtained through applying the proportional clone on A t .

(5) Recombination and hypermutation
Clone group (C t ) is reorganized and hyper mutated. C is set as a new clone group (C t ) and proceeds to step 2.

(6) End
If t is more than G max , D t+1 is exported as the result of the multi-objective optimization algorithm. Otherwise, t is set as t + 1. According to the experience and configuration of the WEDM machine tool, the multiobjective optimization model is developed to obtain high machining efficiency and good surface quality, as shown in Equation (4).  Figure 12 shows the partial solution set of the multi-objective optimization algorithm. Table 7 shows the partial Pareto optimal solution of MRR and SR. In the Pareto optimal solution, MRR is found to be negatively correlated with SR. This means that there is no process parameter combination that can simultaneously obtain the highest MRR and lowest SR. Besides, when a single objective is taken into account, the maximum MRR and the According to the experience and configuration of the WEDM machine tool, the multiobjective optimization model is developed to obtain high machining efficiency and good surface quality, as shown in Equation (4). Figure 12 shows the partial solution set of the multi-objective optimization algorithm. Table 7 shows the partial Pareto optimal solution of MRR and SR. In the Pareto optimal Crystals 2021, 11,1342 14 of 17 solution, MRR is found to be negatively correlated with SR. This means that there is no process parameter combination that can simultaneously obtain the highest MRR and lowest SR. Besides, when a single objective is taken into account, the maximum MRR and the minimum SR can reach 0.501 mm 2 /s and 4.32 µm, respectively. Moreover, this Pareto optimal solution of MRR and SR can be utilized for selecting process parameters in different machining conditions.    Comparing Tables 4 and 7, the comparative results of MRR and SR under the optimized and original process parameters can be obtained, as shown in Table 8. It can be found that for No. 1-2 in Table 4, MRR with the original process parameters is almost the same as that using the optimized process parameters, and SR can be reduced by nearly 6.4%. For No. 3-4 in Table 4, SR with the original process parameters is almost the same as that using the optimized process parameters, and MRR can be increased by 28-34%. This proves that the proposed multi-objective optimization method of NNIA can effectively improve the machining characteristics of the SiCp/Al composite in WEDM.

Verified Experiment
To evaluate the reliability and precision of the Pareto optimal solution, a set of verified experiments is conducted. Table 9 shows the comparison of verified experimental data and predicted data. The relative error between the verified experimental data and predicted data in the Pareto optimal solution ranges from 3.14% to 10.61%. This means that the Pareto optimal solution with NNIA has high reliability and precision.

Conclusions
(1) The methods of material removal of the SiCp/Al composite include sublimating, decomposing and particle shedding. The shedding pit is the primary cause of high surface roughness on the machined surface.
(2) The material removal rate (MRR) is found to increase with an increasing pulse-on time (from 0.265 mm 2 /s to 0.465 mm 2 /s), which first increases and then decreases with an increasing pulse-off time (from 0.374 mm 2 /s to 0.404 mm 2 /s, and to 0.315 mm 2 /s), servo voltage (from 0.408 mm 2 /s to 0.430 mm 2 /s, and to 0.308 mm 2 /s), wire feed (from 0.364 mm 2 /s to 0.404 mm 2 /s, and to 0.351 mm 2 /s) and wire tension (from 0.348 mm 2 /s to 0.404 mm 2 /s, and to 0.364 mm 2 /s). The pulse-on time (the maximum difference up to 0.74 µm) and servo voltage (the maximum difference up to 0.45 µm) are the dominant factors for surface roughness (SR).
(3) The proposed multi-objective optimization method of NNIA can increase the machining speed and reduce the surface roughness of the SiCp/Al composite in WEDM. Specifically, NNIA can increase MRR by 34% and reduce SR by 6.4%.
The Pareto optimal solution by NNIA is proved to possess high reliability and precision, which can be utilized for selecting process parameters in different machining conditions. In future work, we will adopt more direct methods to reveal the machining mechanism of SiCp/Al in EDM/WEDM, such as thermal FEM, molecular dynamics simulation and high-speed observation.