Modeling and analysis of the material removal rate for ultrasonic vibration-assisted polishing of optical glass BK7

The emergence of ultrasonic vibration-assisted polishing technology has effectively improved the machining accuracy and efficiency of hard and brittle materials in modern optical industry, however, the material removal mechanism of ultrasonic vibration-assisted polishing (UVAP) still needs to be further revealed. This paper focuses on the material removal mechanism of ultrasonic vibration-assisted polishing of optical glass (BK7), the application of ultrasonic vibration to axial vibration and the atomization of polishing slurry, the material removal model was established. Based on the analysis of the relationship between the nominal distance d of the polishing pad and the actual contact area distribution, the prediction of the material removal profile is realized. In addition, the effects of different parameters on the material removal rate (MRR) were analyzed, including polishing force, spindle speed, abrasive particle size, ultrasonic amplitude, feed rate, and flow-rate of polishing slurry. Based on the motion equation of abrasive particles, the trajectory of abrasive particles in the polishing slurry was simulated, and the simulation results show that the introduction of the ultrasonic vibration field changes the motion state and trajectory of embedded and free abrasive particles. The new model can not only qualitatively analyze the influence of different process parameters on MRR, but also predict the material removal depth and MRR, providing a possibility for deterministic material removal and a theoretical basis for subsequent polishing of complex curved surfaces of optical glass.

efficiency for hard and brittle materials. Ultrasonic machining was first proposed by G August, s Roland [16] in 1927, and was applied to actual machining in 1920s. In recent decades, scholars have done a lot of research on UVAP. Sanjay Agarwal [17] analyze the material removal mechanism and removal rate of UVAP of glass with a mild steel tool using boron carbide abrasive in water as slurry. The propagation and intersection of radial and median cracks that are induced due to repeated impacts of abrasive particles, which will lead to material removal. The models of shocking force and MRR are derived, which take into account the elastic properties and non-uniformity of abrasive particles, and the mechanical properties of workpiece and abrasive particle Zhang [18] considered the lateral extension of the contact area, the periodic changes of polishing force and contact radius, modeled the local surface profile and material removal distribution function of UVAP, and obtained that larger axial ultrasonic amplitude can improve the Preston coefficient through experiments. The model has good prediction results for local surface profile, maximum contact radius, material removal depth and MRR. Han [19] fixed the workpiece on a rectangle hexahedron ultrasonic sonotrode platform with optimized slots platform to realize the assisted ultrasonic vibration of workpiece, and carried out experimental research on austenitic stainless steel materials. It is concluded that the transverse ultrasonic vibration applied to the workpiece can reduce polishing force, improve machining accuracy, and enhance the proportion of plastic shear effect in material removal. When the vibration direction of the workpiece is consistent with the polishing path direction, the polishing force distribution is more uniform and the surface roughness is lower. When the included angle is 45°, the average polishing force can be reduced by 75.2%. Wang [20] studied the mechanism of ultrasonic vibration assisted grinding (UAG) of hard and brittle materials, established the mathematical model of UAG of brittle materials, conducted in-depth research on the influence of input variables on grinding force and prediction of surface roughness, and obtained the advantages of UAG through theoretical analysis and experimental verification.
In addition to the above studies of UVAP, ultrasonic has also been used for the ultrasonic atomization of polishing slurry. Ultrasonic atomization can produce droplets less than 100 nm, which have small temperature change [21]. Zhang [22] analyzed the effects of operating ultrasonic parameters such as input power and flow rate, and concluded that equipment parameters and liquid physicochemical properties have a great relationship with the mean particle size and particle size distribution of droplets. An effective ultrasonic atomization method was obtained through optimizing the parameters, and the ultrasonic atomization mechanism was established. Kazuhiko Sekiguchi [23] determined the relation between the ultrasonic frequency and the size distribution of the droplets, proposed the generation mechanism of ultrasonic atomization, and concluded that ultrasonic frequency has a negative correlation with atomized droplet diameter of atomized droplet, and is positively correlated with the number and concentration. It is considered that the atomized droplet diameter can be controlled by controlling ultrasonic frequency, power intensity and density of ultrasonic energy transducer.
In conclusion, ultrasonic vibration-assisted polishing is proved to be an effective method for machining hard and brittle materials. However, the theory of ultrasonic vibration-assisted polishing of optical glass is still rarely studied, and the mechanism has not been fully explained so far. Therefore, this paper focuses on the mechanism of ultrasonic vibration-assisted polishing of optical glass (BK7), a predict model of the MRR based on ultrasonic atomization is proposed. The motion of abrasive particles of UVAP was discussed to understand the complexity polishing process. The actual contact area between the polishing pad and the workpiece and its influence on the polishing accuracy are analyzed. At the same time, the influence of different process parameters on the MRR is analyzed, including not only the process parameters, which refer to the polishing force, spindle speed and feed rate, but also the abrasive particle size, as well as ultrasonic amplitude and flow-rate of polishing slurry. This study can provide a theoretical basis for complex surfaces with UVAP.

Motion analysis of ultrasonic vibration
In this study, the cylindrical polishing tool is used to polish the optical glass (BK7). The polishing tool will leave a circular material removal area on the workpiece to be processed. In the polishing process, if the load is constant, the contact area is also constant [15]. Because the polishing head is directly driven by the ultrasonic electric spindle, the equation of motion of the polishing head is： where, 0 z is initial position of polishing head, f and A is ultrasonic frequency and ultrasonic amplitude, respectively.

Contact between the polishing pad and workpiece surface
It is known that the polishing pad surface is much rougher than the workpiece surface, so when microscopic analysis is performed, the workpiece surface is often simplified to a flat surface and the polishing pad surface is considered as a rough surface, as shown in Figure 1. When the polishing pad touches the workpiece surface under a constant load, the distance between the reference surface and the workpiece surface is the nominal distance (d). The polishing pad microscopic surface is randomly distributed, so a Gaussian distribution is used to describe the polishing pad microscopic surface height distribution [23], as shown in Eq (2).
where, h is the average distance between polishing pad and workpiece,  is standard deviation of polishing pad surface height distribution, and = 2 The actual contact area between the polishing pad and the workpiece is given by Where,  is the convexity density, 0 A is the apparent area and 0 R is the average convexity radius. The pressure between pad and workpiece is where, f is stress of single convex peak, t N is the total number of abrasive particles involved in polishing. pw E is the contact modulus of the workpiece and the polishing pad given by

Total number of abrasive particles involved in polishing
In the UVAP, not all the abrasive particles distributed in the polishing liquid participate in the polishing process, and the abrasive particles participating in the polishing are called effective abrasive particles. In this paper, polishing fluid is sprayed by ultrasonic peristaltic pump [16], a calculation method of effective abrasive particles based on ultrasonic atomization is proposed. If  is the utilization rate of effective abrasive particles and N is the total number of abrasive particles ejected by the slurry atomizing nozzle, the effective abrasive number e N can be expressed as: It is known that the mass of a single abrasive is So, the total number of abrasive particles is: Assuming that the ratio of micro convex peaks on the surface of polishing head is  , the effective number of abrasive particles s N participating in the polishing process is:

Relative sliding speed of abrasive particles
As shown in Fig 2, for any point M (x, y) on the polishing path, the speed r V generated by the rotation of the polishing tool can be expressed as: where, p  is spindle speed. When the feed speed of polishing tool is V j , the relative sliding speed of abrasive particles is: The ultrasonic polishing process can be abstracted as a large number of abrasive particles high-frequency impact on brittle material surface. As shown in the Fig 3, based on the indentation fracture model [17], when the abrasive particles impact the workpiece, radial cracks appear on the surface of the workpiece, and transverse cracks appear when the abrasive particles are separated from the workpiece. Assuming that only the broken volume of the abrasive particles is considered in the ultrasonic polishing period, the material removal is determined by the length of the transverse cracks and the depth of the longitudinal cracks. In the ultrasonic vibration period, the material removal amount is approximately equal to the elliptic cone volume shown in Fig. 3, which is the sum of the crushing volumes from the point when the abrasive particles directly contact the workpiece to the point when the abrasive particles just leave the workpiece [20].
where, 1 Then the material removal volume of a single abrasive particle in a single vibration period T is Set a correction factor K (K =1.2), the material removal volume of single abrasive particle on the micro convex peak and in the pore is respectively, 2.5 The material removal rate of numerous abrasive particles The total material removal volume in the UVAP is equal to the sum of the material removal volume of the abrasive particles in the pores and the abrasive particles removal volume in the micro convex peak of the polishing head. Among them, the removal volume of abrasive particles in pores and micro convex peaks in a single cycle are expressed as p V and v V respectively.
Total material removal volume in a single vibration period V is: So, the per unit time MRR in the polishing process can be expressed as: The material removal profile prediction is shown in

3.Analysis and discussions
Theories analyses are carried out in this section for the MRR prediction model, the abrasive motion characteristics in the process of UVAP is analyzed, and the relationship between different process parameters (spindle speed, feed rate, abrasive particle diameter, ultrasonic amplitude and the flow-rate of polishing slurry) and MRR are analyzed in detail respectively to provide a theoretical basis for subsequent research.

The simulation parameters conditions
The detailed parameter is explained in this section. The associated material properties are listed in Table 1 and Table 2. The detailed simulation parameters are listed in Table 3.

Motion analyses of the abrasive particles
From the above analysis, it can be seen the part of the abrasive particles in the polishing liquid embedded the polishing pad micro-raised part with the polishing pad movement. Another part of the abrasive particles exists in the polishing pad pores, forming the free abrasive particles. The two kinds of abrasive particles movement form are slightly different.

Motion analyses of the embedded abrasive particles
From the above formula, the instantaneous velocity V at point M is:   It can be seen from Eq. 20 that the factor that has the greatest influence on the law of abrasive particles motion is the spindle speed. The trajectories of abrasive particles at different spindle speeds are shown in Fig 5. As can be seen from Fig 6, the distribution of abrasive particles motion trajectory under different rotational speeds is very different. When the spindle speed of is higher, the horizontal relative velocity is higher, the instantaneous velocity is higher, the change frequency of abrasive trajectory is faster, and the number of scratches on the workpiece per unit time is increased, which helps to improve the polishing accuracy and efficiency.

Motion analyses of the free abrasive particles
The free abrasive particles are only affected by the axial ultrasonic vibration, so they do periodic simple harmonic vibration in the axial direction. The free abrasive particles motion trajectories under different amplitudes as shown in Fig 8 and the free abrasive particle velocity curves at different frequencies are shown in the Fig 9. From  Fig. 8 and Fig. 9 that the greater the ultrasonic amplitude is, the greater the motion distance of abrasive particles is; the higher the ultrasonic frequency is, the more times of scratching with the workpiece in unit time is, and the higher the MRR is.

The relationship between nominal distance and material removal rate
It can be seen from the above analysis that the actual contact area between the polishing pad and the workpiece is directly related to the nominal diameter. Only the micro convex peak of the polishing pad higher than the nominal diameter can contact the workpiece, and this parameter directly affects the polishing accuracy and efficiency.

The relationship between nominal distance and actual contact area
Under the same conditions, the larger the nominal distance is, the smaller the actual contact area is, and the values of nominal distance and actual contact area can be fitted into a curve by equation, as shown in Fig 10. It can be seen from figure 10 that there is a negative correlation between the nominal distance and the actual contact area. The curve shows that the attenuation rate of the actual contact area and the nominal distance is getting lower and lower, and when the nominal distance is 5 6.25 10   m, the polishing pad is disconnected from the workpiece. It can be seen from the pressure formula / P F S  , when F is constant, the smaller the actual contact area, the greater the pressure between the contact surfaces, and as shown in Fig 11, the greater the nominal distance d, the smaller the actual contact area, the greater the pressure, and the greater the MRR. Fig 11 plots the relationship between the nominal distance d and the MRR, and it can be seen from the figure that the value of d is negatively correlated with the MRR, and when d is greater than 22μm, the MRR decreases rapidly and approaches zero. The MRR at different spindle speeds are shown in Fig 12. In the simulation analysis, the spindle speed varies from 2000 rpm to 8000 rpm, while the other process parameters are fixed as d=0.5μm, F=3N, Vj=400mm/min, f=25kHz, A=5μm, Q=12ml•min -1 . It can be seen the faster spindle speed the greater the MRR, and the two are proportional to each other. It can be seen from Eq 9, 10 and 12, that the spindle speed directly affects the relative speed s V . The faster spindle speed is, the faster the relative speed is, the more times the workpiece surface is scratched in unit time, the larger the removal volume of single abrasive particle is, and the greater the MRR is. However, in the actual polishing, the spindle speed is limited by the stiffness of the machine tool, and can not be increased indefinitely. Only when the processing conditions are met, the spindle speed can be increased as much as possible under the processing conditions to improve the MRR.  F=3N, ωp=8000rpm, f=25kHz, A=5μm, Q=12ml•min -1 .It can be seen that the faster the feed speed the greater the MRR. It can be seen from Eq 10, and 12 that the speed of the spindle directly affects the relative speed j V .The higher the feed speed, the higher the relative speed, but the difference between the feed speed and the spindle speed is very large, and the feed speed has limited influence on the relative speed. In the actual polishing, the feed rate should be selected appropriately to prevent too high feed rate leading to low polishing accuracy or too low feed rate leading to lower machining efficiency. 6.5237x10 -6 6.5237x10 -6 6.5237x10 -6 6.5237x10 -6 6.5237x10 -6 6.5237x10 -6 6.5237x10 -6 6.5238x10 -6 6.5238x10 -6 Material removal rate (m/s) Feed rate (mm/min) 3.4.3 The relationship between abrasive particle size and material removal rate The relationship between different abrasive particle diameters and MRR is shown in Fig 14, which shows that the larger the abrasive particle diameter, the lower the MRR. when the abrasive particle diameter increases from 0.5μm to 1μm, the MRR decreases rapidly. The simulation settings are F=3N, ωp=8000rpm, Vj=400mm/min, f=25kHz, A=5μm, Q=12ml•min -1 . When the abrasive particle size increases, the removal volume of a single abrasive particle increases, but the linear density of the abrasive decreases (when the mass concentration or volume concentration of the polishing slurry is a fixed value), and the number of abrasive particles embedded in the polishing tool per unit area decreases. That is to say, the number of effective abrasive particles decreases greatly, resulting in the decrease of MRR. The relationship between different polishing slurry flow rates and MRR is shown in Fig 16, simulation settings include d=0.5μm, F=3N, ωp=8000rpm, Vj=400mm/min, f=25kHz, A=5μm. It can be seen the flow rate of polishing slurry and MRR are approximately proportional, the more polishing slurry flow rate the greater the MRR. The effective abrasive particles increase with the increase of the slurry flow rate, which can improve the MRR. In the actual polishing, too little polishing slurry flow will result in insufficient abrasive particles, affecting polishing accuracy. However, excessive use of polishing slurry wastes resources, polishing fluid affects the environment, and does not conform to the green processing trend.

4.Conclusions
UVAP can effectively improve the polishing efficiency. In this paper, the process parameters are qualitatively analyzed, and the MRR model is established. In order to better apply the UVAP technology to actual production, this paper draws the following conclusions: 1. The MRR model of UVAP of optical glass (BK7) based on ultrasonic atomization was established by analyzing the abrasive particle motion during the polishing, which can be used to predict the material removal profile and analyze the effect of different process parameters on the MRR.
2. In this paper, the influence of technological parameters (spindle speed, feed speed, abrasive diameter, ultrasonic amplitude, polishing slurry flow rate) on MRR in the polishing process was analyzed qualitatively. Through a detailed analysis of a single variable, it is concluded that the MRR is approximately proportional to the spindle speed and the polishing slurry flow rate. The feed rate and ultrasonic amplitude are positively related to the MRR while negatively related to the diameter of the abrasive particles. Among them, the spindle speed has the greatest influence on MRR, followed by the abrasive diameter and polishing slurry flow rate, while the ultrasonic amplitude and feed speed have less influence on MRR.
3.The characteristics of abrasive particle motion are changed by the introduction of ultrasonic vibration. The trajectory of the embedded abrasive particles is a space sine curve, and the free abrasive particles do a simple harmonic motion in the vertical direction. Each abrasive particle makes intermittent contact with the BK7 specimen only the first half vibration period, which can reduce micro-polishing force and improve machining accuracy. Ultrasonic vibration can greatly improve the impact between the abrasive particles and the workpiece, which can improve the polishing efficiency. Ultrasonic vibration and polishing efficiency are directly affected by ultrasonic frequency and amplitude.