Design of Longitudinal–Torsional Transducer and Directivity Analysis during Ultrasonic Vibration-Assisted Milling of Honeycomb Aramid Material

This paper presents a longitudinal–torsional transducer for use during the ultrasonic vibration-assisted milling (UVAM) of honeycomb aramid material. The mechanism of longitudinal–torsional conversion was analyzed to guide the design of a vibration transducer. The transducer features five spiral grooves around the front cover plate, which function under the excitation of a group of longitudinal piezoelectric ceramics. A portion of the longitudinal vibration was successfully converted into torsional vibration. The resonant frequency, longitudinal vibration displacement and torsional amplitude at the top of the disk milling cutter were 24,609 Hz, 19 μm and 9 μm, respectively. In addition, the directivity of the longitudinal–torsional transducer was theoretically analyzed. Compared with conventional milling, UVAM with the longitudinal-torsional could significantly reduce the cutting force (40–50%) and improve the machining stability.


Introduction
Honeycomb aramid materials have been widely used in the field of aerospace because of their excellent structural and mechanical properties. However, honeycomb core deformation, the risk of hole grid collapse, and a low yield are among challenges faced by traditional machining methods in the machining of honeycomb aramid materials. Based on the resonance effect, ultrasonic vibration-assisted milling (UVAM) can generate a highfrequency relative motion between the tool and the workpiece, which is conducive to improving processing stability and surface quality.
At present, most ultrasonic transducers used in UVAM vibrate along the longitudinal direction. However, longitudinal vibration has a limited effect on the improvement in terms of workpiece machinability [1,2]. In order to improve the machining stability of honeycomb aramid materials during UVAM, this paper presents a longitudinal-torsional transducer, which can excite the milling cutter to vibrate in both longitudinal and torsional directions simultaneously. The longitudinal and torsional vibrations of the longitudinal-torsional transducer can be realized in two ways: the vibration mode conversion of the transducer and the tangential polarization of a piezoelectric ceramic stack.
In the process of vibration mode conversion, torsional vibrations in the transmission rod can be obtained based on the coupling and degradation of horn form longitudinal vibration. Patrick Harkness [3] proposed two design methods for longitudinal-torsional mode transitions. Celaya [4] proposed a booster that can improve the surface quality but with a low amplitude. Liu [5] demonstrated two forms of thread grooves which can degrade longitudinal vibrations into longitudinal-torsional vibrations. Yahya [6], Gao [7], Yuan [8], Zhang [9] and Wu [10] developed L-T transducers based on theoretical and simulation analysis and applied them to the processing of Ti-6Al-4V. Budairi [11] used numerical, analytical and experimental methods to evaluate transducer performance under different excitations based on the L-T response. In addition, by setting different degradation modes on the cylinder surface, the same method can be applied to ultrasonic motors [12][13][14]. Shen [15] conducted an assisted ultrasonic vibration micro-milling process on aluminum alloy material and used a combination of machining parameters to demonstrate that ultrasonic vibration results in a significant improvement in terms of surface roughness; Chen [16] came to the same conclusion for silicon carbide ceramic grinding. Cleary [17] and Wang [18] proposed an L-T ultrasound needle device for medical treatment and implantation applications which provides a potential application scenario for L-T mode devices. Numanoglu [19] and Akbaş [20] determined the vibrations of composite rods and beams for different boundary conditions using different theoretical methods, providing a theoretical foundation for the derivation of an L-T transducer. However, in the processing of honeycomb aramid materials by UVAM, there is a lack of research on longitudinaltorsional transducers, which limits the possible improvements to be made regarding the machinability of honeycomb aramid materials.
The tangential polarization of a piezoelectric ceramic stack can be realized by combining axial piezoelectric ceramics and tangential piezoelectric ceramics. Kim [21] proposed a transducer with two frequencies by combining two piezoelectric torsional discs [22] designed a torsional device. Harada [23] used a torsional transducer and a microporous plate to generate droplets. However, due to the complex fabrication process associated with tangentially polarized piezoelectric wafers, problems such as electrical breakdown and incomplete polarization will occur when the chip size is large. Therefore, it is difficult to develop torsional vibration transducers with higher frequencies using the current level of ceramic production technology.
In addition, the stability lobes diagram is of great significance for the optimization of equipment structure and process parameters, with it being affected by the directivity of a transducer during UVAM. At present, the relevant research primarily focuses on analysis methods for stability [24], the correctness of the stability lobes diagram [25] and increasing the equivalent stiffness by filling the honeycomb material with viscoelastic damping material [26]. Suárez [27] showed that the application of ultrasonic milling can improve the form of tool wear patter and render the evolution more stable. However, the research on the directivity of transducers during ultrasonic vibration-assisted milling (UVAM) of honeycomb aramid material is relatively sparse. This paper presents a longitudinal-torsional transducer for use during the ultrasonic vibration-assisted milling (UVAM) of honeycomb aramid material. To avoid energy loss, the front cover plate and the horn of transducer were designed as a whole, and spiral grooves were machined on the front cover plate. Based on the theory of longitudinal-torsional transformation, the spiral groove parameters on the front cover plate were analyzed, and the directivity of the longitudinal-torsional transducer was studied theoretically. Finally, the performance of the longitudinal-torsional transducer and its processing stability were verified by experiments.

Theoretical Analysis of Longitudinal-Torsional Transducer
Ultrasonic transducers can be divided into two types according to different internal working mediums: sandwich piezoelectric ceramic transducers and magnetostrictive transducers [28]. Sandwich piezoelectric ceramic transducers are widely used because of the advantages of uniform vibration, a simple structure andease in terms of forming and processing, etc. This type of transducer, which is composed of two quarter wavelengths in length, as shown in Figure 1, was selected for this study. Where Z i is the characteristic acoustic impedance of each part of the material, X i is the boundary of each component, L i is the length of each component (i = 1,2, . . . , 8).
in length, as shown in Figure 1, was selected for this study. Where Zi is the characteristic acoustic impedance of each part of the material, Xi is the boundary of each component, Li is the length of each component(i=1,2,…,8). According to the theory related to Newtonian dynamics in the principle of vibration transmission, the equation for vibration at any section can be expressed as: where S is area function of the rod at any cross section, S = S( ),  is the displacement function of the mass, is the stress function of the rod, which can be represented as and c represent the circular frequency and propagation of the speed of sound wave, respectively. Therefore, the vibration equation for the transducer at any section can be written as follows: where k represents the number of circular waves and k = /c. For metallic materials, / cE  = ; for piezoelectric ceramic materials, the equivalent sound velocity of the longitudinal wave is approximately expressed as: (2), the vibration velocity for longitudinal vibration can be expressed as follows: Since S(x) is constant, Equation (3) can be simplified as: The general solution for the vibration velocity, i.e., Equation (4), can be obtained as: According to the theory related to Newtonian dynamics in the principle of vibration transmission, the equation for vibration at any section can be expressed as: where S is area function of the rod at any cross section, S = S(x), ξ is the displacement function of the mass, σ is the stress function of the rod, which can be represented as σ = σ(x) = E ∂ξ ∂X and ρ and E are the density and Young's modulus of horn, respectively. In the process of UVAM, the vibration form of the transducer is simple harmonic vibration. ξ = ξ · e jωt , κ = ω c and c represent the circular frequency and propagation of the speed of sound wave, respectively. Therefore, the vibration equation for the transducer at any section can be written as follows: where k represents the number of circular waves and k = ω/c. For metallic materials, c = E/ρ; for piezoelectric ceramic materials, the equivalent sound velocity of the longitudinal wave is approximately expressed as: c ≈ c 1 − k 2 33 = 1/(ρS E 33 ). Substituting v = jωt into Equation (2), the vibration velocity for longitudinal vibration can be expressed as follows: Since S(x) is constant, Equation (3) can be simplified as: The general solution for the vibration velocity, i.e., Equation (4), can be obtained as: Substituting Equation (5) into the stress distribution expression results in: where j represents an imaginary sign Z is the characteristic acoustic impedance and Z = ρcs.
During the theoretical analysis of ultrasonic transducers, the system is usually divided into two parts using the nodal plane. In order to simplify the calculation, the nodal plane of an ultrasonic transducer was designed on the joint surfaces of the piezoelectric crystal stack and front cover plate. One part is composed of the rear cover plate, electrode plate and piezoelectric ceramic plate, and the other part is composed of the front cover plate and amplitude transformer, as shown in Figure 1. Vf and Vb are the vibration velocities of the end face of the rear cover plate and the front face of the piezoelectric ceramic plate, respectively. When the ultrasonic waves propagate at the medium partition interface, the boundary conditions on the left side of the nodal plane can be introduced as follows: Substituting Equation (7) into Equations (5) and (6), the following results can be obtained: 1 Z 4 cos k 3 l 3 sin k 4 l 4 + 1 Z 3 sin k 3 l 3 cos k 4 l 4 (Z 1 sin k 1 l 1 cos k 2 l 2 + Z 2 cos k 1 l 1 sin k 2 l 2 ) (9) where Z i is the characteristic acoustic impedance of each part of the material on the left side of the nodal plane and Z i = ρcs i .
The vibration velocity and stress distribution of each part can be calculated as follows: sin k 5 l 5 sin k 5 x 5 + cos k 5 x 5 Z 3 Z 4 sin k 3 l 3 sin k 4 l 4 − cos k 3 l 3 cos k 4 l 4 ) Z 1 Z 2 sin k 1 l 1 sin k 2 l 2 − cos k 1 l 1 cos k 2 l 2 ) − 1 Z 2 cos k 3 l 3 sin k 4 l 4 + 1 Z 3 sin k 3 l 3 cos k 4 l 4 (Z 1 sin k 1 l 1 cos k 2 l 2 + Z 2 cos k 1 l 1 sin k 2 l 2 )] Micromachines 2022, 13, 2154 5 of 18 Substituting the coefficients of Equations (8) and (9) into the equations of Equation (7), the frequency equation on the left of the ultrasonic transducer can be written as: There are only a rear cover plate and four piezoelectric ceramics on the left side of transducer, and the relevant parameters for piezoelectric ceramics are known. Therefore, Equation (12) can be simplified as: tan k 1 l 1 tan k l 5 − 4 tan 2 k l + 6 tan 2 k l − tan 4 k l = 1 The right side of the section surface, AB, is the front cover plate of the longitudinaltorsional transducer and the index cylinder composite variable amplitude rod. According to the same method and results can be derived from the frequency equation, the vibration speed equation and the general solution. As a result, the frequency equation on the right side of the nodal plane is:

Structural Design of Longitudinal-Torsional Transducer
As shown in Figure 1, 1 and 3 denote piezoelectric ceramics, 2 denotes electrode sheets, 4 represents the rear cover plate and 5 and 6 represent the front cover plates. The vibration velocity of the outer surface of the front cover plate is V b , and the elastic force can be written as F = Z w V b , where Z w represents the input impedance. The vibration speed of the outer surface of the rear cover plate is V f . Due to the front and rear ends of the longitudinaltorsional transducer being exposed to air, it can be considered that the elastic force is zero. The waveforms of each part under different input frequencies on both sides of the transducer node planes are shown in Figure  In order to generate longitudinal-torsional vibration, the shape and section of the horn needs to be changed to realize the transformation of vibration modes. In this paper, a spiral groove sandwich longitudinal-torsional transducer is presented. The spiral groove structure, which can convert the longitudinal vibration into longitudinal vibration and torsional vibration, was machined at the front end of the sandwich longitudinal-torsional transducer. As a result, there will be torsional and longitudinal vibration components at the tip of tool, as shown in Figure 3. Where A-B denotes the nodal surface, C-D denotes the helix start line, F denotes the total force, Ft denotes the torsional force, Fl denotes the longitudinal force, and  denotes the helix angle. In Figure 3, the solid line between A-B shows the node plane in the longitudinaltorsional vibration mode, and the dashed line between C-D shows the partition interface between the longitudinal vibration part and longitudinal-torsional vibration part. At the left end of the partition interface, the piezoelectric ceramics generate axial vibration under the excitation of a sinusoidal signal generated by an ultrasonic generator, and the force is a longitudinal force, F. As shown in Figure 4, when the longitudinal wave moves forward along the axis of the cylinder, it will produce a transverse wave and a longitudinal wave when encountering the spiral structure. Due to the effect of the groove, the longitudinal force is separated into a longitudinal force, l F , along the axis direction and a tangential force, t F , perpendicular to the radius. In order to generate longitudinal-torsional vibration, the shape and section of the horn needs to be changed to realize the transformation of vibration modes. In this paper, a spiral groove sandwich longitudinal-torsional transducer is presented. The spiral groove structure, which can convert the longitudinal vibration into longitudinal vibration and torsional vibration, was machined at the front end of the sandwich longitudinal-torsional transducer. As a result, there will be torsional and longitudinal vibration components at the tip of tool, as shown in Figure 3. Where A-B denotes the nodal surface, C-D denotes the helix start line, F denotes the total force, F t denotes the torsional force, F l denotes the longitudinal force, and ϕ denotes the helix angle. In order to generate longitudinal-torsional vibration, the shape and section of the horn needs to be changed to realize the transformation of vibration modes. In this paper, a spiral groove sandwich longitudinal-torsional transducer is presented. The spiral groove structure, which can convert the longitudinal vibration into longitudinal vibration and torsional vibration, was machined at the front end of the sandwich longitudinal-torsional transducer. As a result, there will be torsional and longitudinal vibration components at the tip of tool, as shown in Figure 3. Where A-B denotes the nodal surface, C-D denotes the helix start line, F denotes the total force, Ft denotes the torsional force, Fl denotes the longitudinal force, and  denotes the helix angle. In Figure 3, the solid line between A-B shows the node plane in the longitudinaltorsional vibration mode, and the dashed line between C-D shows the partition interface between the longitudinal vibration part and longitudinal-torsional vibration part. At the left end of the partition interface, the piezoelectric ceramics generate axial vibration under the excitation of a sinusoidal signal generated by an ultrasonic generator, and the force is a longitudinal force, F. As shown in Figure 4, when the longitudinal wave moves forward along the axis of the cylinder, it will produce a transverse wave and a longitudinal wave when encountering the spiral structure. Due to the effect of the groove, the longitudinal force is separated into a longitudinal force, l F , along the axis direction and a tangential force, t F , perpendicular to the radius. In Figure 3, the solid line between A-B shows the node plane in the longitudinaltorsional vibration mode, and the dashed line between C-D shows the partition interface between the longitudinal vibration part and longitudinal-torsional vibration part. At the left end of the partition interface, the piezoelectric ceramics generate axial vibration under the excitation of a sinusoidal signal generated by an ultrasonic generator, and the force is a longitudinal force, F. As shown in Figure 4, when the longitudinal wave moves forward along the axis of the cylinder, it will produce a transverse wave and a longitudinal wave when encountering the spiral structure. Due to the effect of the groove, the longitudinal force is separated into a longitudinal force, F l , along the axis direction and a tangential force, F t , perpendicular to the radius.      Where C-D denotes the helix start line, r is the radius of the circle formed by the helix. The inner diameter and outer diameter of the front cover plate are r 1 and r 2 , respectively. ϕ is the angle between the spiral groove and the axis of front cover plate. The length of the cylinder in front of the groove is l 1 , and the length of rest is l 2 . Due to the torsional vibration generated by the spiral structure, it can be considered that some sections (l 1 ) do not produce torsional vibration, while the rest produces not only longitudinal vibration but also torsional vibration. To simplify the analysis, it is assumed that the groove is an ideal geometric segment, i.e., the width of the groove is infinitely small.   According to vibration theory, all tangential force components will produce The longitudinal force and tangential force can be determined by the following equation: where ϕ is the inclination angle of the helical groove. According to vibration theory, all tangential force components will produce a torsional movement for any section of a spiral structure, and the total torque can be expressed as: where s = π r 2 2 − r 1 2 represents the cross-sectional area of the cylinder, r is the section radius at any position, f is the tangential force per unit area and ds = 2πrdr is the element area. Therefore, the tangential force per unit area can be obtained as: Equation (16) can be simplified as follows: According to Equation (18), the torque, M, increases with the increases in rotation angle, ϕ, and notch radius, r 1 , while the sectional area and polar moment of inertia decrease with the increase in notch radius, r 1 . Therefore, the torsional vibration amplitude increases with the increase in rotation angle, ϕ, and notch radius, r 1 . Due to the existence of the spiral structure, in addition to longitudinal vibration, there is also a torsional vibration in the front cover plate. Therefore, the vibration mode of the front cover plate of the transducer is a longitudinal-torsional composite vibration.
In order to avoid the influence of spiral groove parameters on the output amplitude, it was necessary to simulate the spiral structure with different parameters and select the optimal chute parameters. In this study, the single factor analysis method was used for analysis by changing the spiral parameters, such as the inclination angle, width, depth, length and distance from the small end, etc. The specific simulation results are shown in Figures 6 and 7.  r . Due to the existence of the spiral structure, in addition to longitudinal vibration, there is also a torsional vibration in the front cover plate. Therefore, the vibration mode of the front cover plate of the transducer is a longitudinal-torsional composite vibration. In order to avoid the influence of spiral groove parameters on the output amplitude, it was necessary to simulate the spiral structure with different parameters and select the optimal chute parameters. In this study, the single factor analysis method was used for analysis by changing the spiral parameters, such as the inclination angle, width, depth, length and distance from the small end, etc. The specific simulation results are shown in Figures 6 and 7.  As shown in Figure 6, the groove depth has a great influence on the longitudinal amplitude, which changes as a quadratic function. As the groove width becomes wider, the longitudinal amplitude first decreases, then increases and finally decreases. When the groove width is 0.5 mm and 2 mm, the longitudinal amplitude reaches the maximum value. The longitudinal amplitude value fluctuates between 17 mm-20 mm. The effect of spiral groove parameters on the torsional vibration is illustrated in Figure 7. As shown in Figure 7, the influence of spiral groove structural parameters is greater than the longitudinal amplitude and natural frequency. As shown in Figure 6, the groove depth has a great influence on the longitudinal amplitude, which changes as a quadratic function. As the groove width becomes wider, the longitudinal amplitude first decreases, then increases and finally decreases. When the groove width is 0.5 mm and 2 mm, the longitudinal amplitude reaches the maximum value. The longitudinal amplitude value fluctuates between 17 mm-20 mm. The effect of spiral groove parameters on the torsional vibration is illustrated in Figure 7. As shown in Figure 7, the influence of spiral groove structural parameters is greater than the longitudinal amplitude and natural frequency.
To sum up, the structural parameters of spiral grooves have a significant influence on the torsional amplitude. The reason for this may be that when the width, length and depth of the spiral groove change, more energy is converted into torsional vibration, which makes the longitudinal vibration decrease and torsional vibration increase. Therefore, when considering the influence of spiral groove structural parameters on the longitudinal-torsional vibration, the length and width of the spiral groove should be considered. In order to make the longitudinal-torsional composite vibration system output a larger longitudinal torsional amplitude towards the design frequency of 25 kHz, the specific parameters of the spiral groove are shown in Table 1. The vibration modes and resonant frequencies of the longitudinal-torsional transducer were confirmed using finite element method (FEM). Figure 8 illustrates the modal velocity vector diagram of the longitudinal-torsional transducer. As shown in Figure 8, the resonance frequency of the transducer which corresponds to the longitudinaltorsional vibration is 24609 Hz. This indicates that the longitudinal vibration that was To sum up, the structural parameters of spiral grooves have a significant influence on the torsional amplitude. The reason for this may be that when the width, length and depth of the spiral groove change, more energy is converted into torsional vibration, which makes the longitudinal vibration decrease and torsional vibration increase. Therefore, when considering the influence of spiral groove structural parameters on the longitudinaltorsional vibration, the length and width of the spiral groove should be considered. In order to make the longitudinal-torsional composite vibration system output a larger longitudinal torsional amplitude towards the design frequency of 25 kHz, the specific parameters of the spiral groove are shown in Table 1. The vibration modes and resonant frequencies of the longitudinal-torsional transducer were confirmed using finite element method (FEM). Figure 8 illustrates the modal velocity vector diagram of the longitudinal-torsional transducer. As shown in Figure 8, the resonance frequency of the transducer which corresponds to the longitudinal-torsional vibration is 24,609 Hz. This indicates that the longitudinal vibration that was input was successfully converted into a longitudinal-torsional composite vibration. Celaya [4] measurements with a vibrometer showed that the bending vibration mode was obtained at 20.1 kHz, with the maximum vibration amplitude obtained at 8~10 µm at the tip of the tool. In contrast to the UVAM in this paper, the maximum longitudinal vibration can reach 19 µm, and the torsional amplitude can reach 9 µm, as shown in Figure 9.

Directivity of Longitudinal-Torsional Transducer
During ultrasonic propagation, the conversion of longitudinal and transverse waves in transducers is essentially a sound field problem. Sound fields are generated by transducers, and most transducers generate a finite beam sound field. This sound field is closely related to the shape, size, vibration mode and working parameters of the transducer. Therefore, it was necessary to analyze the directivity of the transducer and evaluate the directivity of the finite beam sound field.

Distribution of Radiated Sound Field
The transducer source selected in this study was a piston type sound source, which generally refers to a plane shaped vibrator. When the sound source vibrates along the normal direction of the plane, the amplitude and phase of the vibration velocity at any point are the same. Based on the theory of sound field calculation [29], the sound field radiated by the transducer within a finite size can be analyzed according to the linear superposition principle. The effective radiation surface of the transducer is regarded as the combination of countless point sound sources. The sound pressure at a certain point in the radiation field is the result of the superposition of sound pressure generated at that point by all point sources. If there is a surface sound source of an arbitrary shape, the amplitude and phase of vibration at each point on the surface may be different. The

Directivity of Longitudinal-Torsional Transducer
During ultrasonic propagation, the conversion of longitudinal and transverse waves in transducers is essentially a sound field problem. Sound fields are generated by transducers, and most transducers generate a finite beam sound field. This sound field is closely related to the shape, size, vibration mode and working parameters of the transducer. Therefore, it was necessary to analyze the directivity of the transducer and evaluate the directivity of the finite beam sound field.

Distribution of Radiated Sound Field
The transducer source selected in this study was a piston type sound source, which generally refers to a plane shaped vibrator. When the sound source vibrates along the normal direction of the plane, the amplitude and phase of the vibration velocity at any point are the same. Based on the theory of sound field calculation [29], the sound field radiated by the transducer within a finite size can be analyzed according to the linear superposition principle. The effective radiation surface of the transducer is regarded as the combination of countless point sound sources. The sound pressure at a certain point in the radiation field is the result of the superposition of sound pressure generated at that point by all point sources. If there is a surface sound source of an arbitrary shape, the amplitude and phase of vibration at each point on the surface may be different. The

Directivity of Longitudinal-Torsional Transducer
During ultrasonic propagation, the conversion of longitudinal and transverse waves in transducers is essentially a sound field problem. Sound fields are generated by transducers, and most transducers generate a finite beam sound field. This sound field is closely related to the shape, size, vibration mode and working parameters of the transducer. Therefore, it was necessary to analyze the directivity of the transducer and evaluate the directivity of the finite beam sound field.

Distribution of Radiated Sound Field
The transducer source selected in this study was a piston type sound source, which generally refers to a plane shaped vibrator. When the sound source vibrates along the normal direction of the plane, the amplitude and phase of the vibration velocity at any point are the same. Based on the theory of sound field calculation [29], the sound field radiated by the transducer within a finite size can be analyzed according to the linear superposition principle. The effective radiation surface of the transducer is regarded as the combination of countless point sound sources. The sound pressure at a certain point in the radiation field is the result of the superposition of sound pressure generated at that point by all point sources. If there is a surface sound source of an arbitrary shape, the amplitude and phase of vibration at each point on the surface may be different. The sound source surface, S, is decomposed into an infinite number of small panels, dS. On each panel, dS, the vibration of each point can be seen as uniform, so these panels, dS, are considered to be point sound sources. The vibration law of the point source located in space (x, y, z) can be assumed to be: where u a is the amplitude of the surface source, α is a position function, which is related to the position of the point source in space. The point source intensity, dp, can be expressed as the integral of the amplitude on the panel: where h(x, y, z) is the radial propagation distance function from the center of the point source; κ represents the wave number; κ = 2π f /c 0 ;Q 0 refers to the point source intensity and represents the volumetric velocity amplitude of the spherical source formed by the point source; and Q 0 = 4πr 2 0 u a . Total sound pressure can be obtained as follows: Figure 10 illustrates the sound field distributions of two common piston radiation sources, where r is the radius of the circular piston (the rectangular piston length and width are a and b). The center of a circular piston transducer is set as coordinate origin o, and the plane of piston is located at the xoy plane. The sound field is rotationally symmetric with respect to the z axis passing through the center of the piston. Therefore, for an observation point p, in the sound field, the distance from the origin, o, is r, and the included angle between the position vector, r, and z axis is θ.
Based on the principle of point source combination, the radiated sound field of the circular piston was analyzed. The piston surface can be divided into an infinite number of small panels. The sound pressure generated at the observation point, p, by the panel with polar diameter, r 0 , and polar angle, ξ, is expressed as follows: dp = j kρ 0 c 0 2πr u a dSe j(ωt−kh) For circular pistons, each face element vibrates in the same phase. In order to facilitate calculations, the radiation sound pressure of the whole piston can be obtained by superimposing the face elements: Where the initial position function can be expressed as α(x, y, z) = 0. For the region of r > a, the amplitude can be approximated as the distance from the piston center to the observation point, expressed as r.
Based on the Bessel function [29], The integral of Equation (24) can be expressed as the following equation: The radial velocity of the mass, meanwhile, can be expressed as: As can be seen from Equation (26), in the area far away from the piston, the sound pressure decreases inversely with distance. However, the sound pressure is uneven along different directions at the same distance. Because the sound waves from different positions have different phases when reaching the observation point, the interference results in the directivity of the sound field.  (20) where ( , , ) h x y z is the radial propagation distance function from the center of the point source;  represents the wave number;  Figure 10 illustrates the sound field distributions of two common piston radiation sources, where r is the radius of the circular piston (the rectangular piston length and width are a and b). The center of a circular piston transducer is set as coordinate origin o, and the plane of piston is located at the xoy plane. The sound field is rotationally symmetric with respect to the z axis passing through the center of the piston. Therefore, for an observation point p , in the sound field, the distance from the origin, o, is r, and the included angle between the position vector, r, and z axis is . kc dp j u dSe r (22) For circular pistons, each face element vibrates in the same phase. In order to facilitate calculations, the radiation sound pressure of the whole piston can be obtained by superimposing the face elements: Where the initial position function can be expressed as ( ) , , 0 xyz  = . For the region of r ＞ a, the amplitude can be approximated as the distance from the piston center to

Directivity Functions and Properties
The directivity function, i.e., the sensitivity of the directivity of the transducer, is a spatial distribution function used to describe the sound field radiated by a sound source (free far field) or the resolution and anti-interference ability of the receiver. The directivity response pattern is used to characterize the energy loss and reflect the stability during processing.
Based on the Bessel function [29], Therefore, the directivity of the piston sound source can be expressed as follows: It can be seen that the directivity is related to the relative ratio of the size of the piston to the wavelength. Under certain conditions in terms of the sound field environment (a certain ultrasonic frequency, f , and a constant acoustic propagation medium), it is possible to analyze the directivity by changing the radius of the piston transducer. Figure 11 shows the directivity pattern of the piston transducer with different radii of 30 mm, 40 mm and 50 mm. As shown in Figure 11, no matter what the radius of the piston transducer is, side lobes will appear on both sides of the main lobe. As the radius of the piston transducer increases, the width of the main lobe narrows, the number of side lobes increase and the directivity of the transducer becomes stronger. With the increase in radius, the position of the maximum sidelobe moves to the middle. When θ is close to 90 • , the directivity tends to be 0. (c) a = 50mm

Experimental Setup
During UVAM of honeycomb aramid material, the wireless power supply system for the longitudinal-torsional transducer was pre-installed on the machine spindle in order to ensure that there was no effect on the machine's electric spindle, as shown in Figure 12. The high-frequency power of the wireless power supply system was generated by an ultrasonic generator, which was switched on and off by a control valve or control switch, and the interface was generally guaranteed for three to four years of service life. The force date was recorded by a three-way force sensor, which was mounted on the acrylic fixing plate of the worktable.

Experimental Setup
During UVAM of honeycomb aramid material, the wireless power supply system for the longitudinal-torsional transducer was pre-installed on the machine spindle in order to ensure that there was no effect on the machine's electric spindle, as shown in Figure 12. The high-frequency power of the wireless power supply system was generated by an ultrasonic generator, which was switched on and off by a control valve or control switch, and the interface was generally guaranteed for three to four years of service life. The force date was recorded by a three-way force sensor, which was mounted on the acrylic fixing plate of the worktable.

Experimental Setup
During UVAM of honeycomb aramid material, the wireless power supply system for the longitudinal-torsional transducer was pre-installed on the machine spindle in order to ensure that there was no effect on the machine's electric spindle, as shown in Figure 12. The high-frequency power of the wireless power supply system was generated by an ultrasonic generator, which was switched on and off by a control valve or control switch, and the interface was generally guaranteed for three to four years of service life. The force date was recorded by a three-way force sensor, which was mounted on the acrylic fixing plate of the worktable.  It is difficult to measure the cutting force in x, y and z directions simultaneously, so a three-way force sensor was arranged on the acrylic plate. First, the honeycomb aramid material was fixed on the acrylic plate, and when the disc cutter cut the material, the three component forces generated were transmitted to the three-way force sensor through the acrylic plate. The sensor amplified the received signals through a charge amplifier, stored them in a data acquisition card, and finally presented the data through a data acquisition device (computer). Figure 13 shows a schematic diagram of the cutting force system. The test field installation is shown in Figure 14.
It is difficult to measure the cutting force in x, y and z directions simultaneously, so a three-way force sensor was arranged on the acrylic plate. First, the honeycomb aramid material was fixed on the acrylic plate, and when the disc cutter cut the material, the three component forces generated were transmitted to the three-way force sensor through the acrylic plate. The sensor amplified the received signals through a charge amplifier, stored them in a data acquisition card, and finally presented the data through a data acquisition device (computer). Figure 13 shows a schematic diagram of the cutting force system. The test field installation is shown in Figure 14.

Cutting Force
Compared with traditional milling, when the spindle speed is 1500 r/min and the cutting depth is 2 mm, the cutting forces along three directions are reduced by approximately 40-50%. The measurement results for the cutting force during UVAM are shown in Figure 15. It is difficult to measure the cutting force in x, y and z directions simultaneously, so a three-way force sensor was arranged on the acrylic plate. First, the honeycomb aramid material was fixed on the acrylic plate, and when the disc cutter cut the material, the three component forces generated were transmitted to the three-way force sensor through the acrylic plate. The sensor amplified the received signals through a charge amplifier, stored them in a data acquisition card, and finally presented the data through a data acquisition device (computer). Figure 13 shows a schematic diagram of the cutting force system. The test field installation is shown in Figure 14.

Cutting Force
Compared with traditional milling, when the spindle speed is 1500 r/min and the cutting depth is 2 mm, the cutting forces along three directions are reduced by approximately 40-50%. The measurement results for the cutting force during UVAM are shown in Figure 15.

Cutting Force
Compared with traditional milling, when the spindle speed is 1500 r/min and the cutting depth is 2 mm, the cutting forces along three directions are reduced by approximately 40-50%. The measurement results for the cutting force during UVAM are shown in Figure 15.

Directivity Evaluation
Concerning machining in the context of UVAM, machining stability can be expressed in terms of the maximum depth of the cut produced by the milling tool at different cutting speeds. The maximum axial cutting depth under different processing methods when machining honeycomb aramid material with a length, width and thickness of 100 mm, 100 mm and 12 mm, respectively, is illustrated in Figure 12. As shown in Figure 16, when the spindle speed is low, the maximum cutting depth during UVAM is higher than that of traditional machining, and the directivity of the transducer directivity function is also stronger, which indicates that UVAM can significantly improve the stability of the milling process.
(a) Traditional milling

Directivity Evaluation
Concerning machining in the context of UVAM, machining stability can be expressed in terms of the maximum depth of the cut produced by the milling tool at different cutting speeds. The maximum axial cutting depth under different processing methods when machining honeycomb aramid material with a length, width and thickness of 100 mm, 100 mm and 12 mm, respectively, is illustrated in Figure 12. As shown in Figure 16, when the spindle speed is low, the maximum cutting depth during UVAM is higher than that of traditional machining, and the directivity of the transducer directivity function is also stronger, which indicates that UVAM can significantly improve the stability of the milling process.

Directivity Evaluation
Concerning machining in the context of UVAM, machining stability can be expressed in terms of the maximum depth of the cut produced by the milling tool at different cutting speeds. The maximum axial cutting depth under different processing methods when machining honeycomb aramid material with a length, width and thickness of 100 mm, 100 mm and 12 mm, respectively, is illustrated in Figure 12. As shown in Figure 16, when the spindle speed is low, the maximum cutting depth during UVAM is higher than that of traditional machining, and the directivity of the transducer directivity function is also stronger, which indicates that UVAM can significantly improve the stability of the milling process.
(a) Traditional milling