Next Article in Journal
Research on Effects of Dust Removal Using Ultrasonic Vibrators
Next Article in Special Issue
Free-Form Surface Partitioning and Simulation Verification Based on Surface Curvature
Previous Article in Journal
Internet of Things-Based ECG and Vitals Healthcare Monitoring System
Previous Article in Special Issue
An Intelligent Logic-Based Mold Breakout Prediction System Algorithm for the Continuous Casting Process of Steel: A Novel Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Engineering Training Center, Changchun University of Technology, Changchun 130012, China
2
School of Mechanical and Electrical Engineering, Changchun University of Technology, Changchun 130012, China
3
Changchun Yidong Clutch Co., Ltd., Changchun 130012, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(12), 2154; https://doi.org/10.3390/mi13122154
Submission received: 14 November 2022 / Revised: 3 December 2022 / Accepted: 5 December 2022 / Published: 6 December 2022
(This article belongs to the Special Issue Advanced Manufacturing Technology and Systems, 2nd Edition)

Abstract

:
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.

1. 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 high-frequency 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. Numanoğlu [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 longitudinal–torsional 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.

2. Design of Longitudinal–Torsional Transducer

2.1. 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 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:
( S σ ) x d x = S ρ 2 ξ t 2 d x
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:
2 ξ x 2 + 1 S S x ξ x + k 2 ξ = 0
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 33 2 = 1 / ( ρ S 33 E ) . Substituting v = jωt into Equation (2), the vibration velocity for longitudinal vibration can be expressed as follows:
2 v x 2 + 1 S S x v x + k 2 v = 0
Since S(x) is constant, Equation (3) can be simplified as:
2 v x 2 + k v 2 = 0
The general solution for the vibration velocity, i.e., Equation (4), can be obtained as:
v ( x ) = A sin k x + B cos k x
Substituting Equation (5) into the stress distribution expression results in:
F ( x ) = E S j ω v x = j ρ c s ( A cos k x B sin k x ) = j Z ( A cos k x B sin k x )
where j represents an imaginary sign Z is the characteristic acoustic impedance and Z = ρ c s .
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:
v 1 ( 0 ) = v b v 1 ( l 1 ) = v 2 ( 0 ) v 2 ( l 2 ) = v 3 ( 0 ) v 3 ( l 3 ) = v 4 ( 0 ) v 4 ( l 4 ) = v 5 ( 0 ) v 5 ( l 5 ) = 0 F 1 ( 0 ) = 0 F 1 ( l 1 ) = F 2 ( 0 ) F 2 ( l 2 ) = F 3 ( 0 ) F 3 ( l 3 ) = F 4 ( 0 ) F 4 ( l 4 ) = F 5 ( 0 )
Substituting Equation (7) into Equations (5) and (6), the following results can be obtained:
A 1 = 0 A 2 = Z 1 Z 2 v b sin k 1 l 1 A 3 = Z 1 Z 3 v b sin k 1 l 1 cos k 2 l 2 Z 2 Z 3 v b cos k 1 l 1 sin k 2 l 2 A 4 = Z 1 Z 4 v b sin k 1 l 1 ( Z 3 Z 2 sin k 2 l 2 sin k 3 l 3 cos k 2 l 2 sin k 3 l 3 )                     1 Z 4 v b cos k 1 l 1 ( Z 2 sin k 2 l 2 cos k 3 l 3 + Z 3 cos k 2 l 2 sin k 3 l 3 ) A 5 = cos k 5 l 5 sin k 5 l 5 v b [ ( 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 ) ( 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 ) ]
B 1 = v b B 2 = v b cos k 1 l 1 B 3 = Z 1 Z 2 v b sin k 1 l 1 sin k 2 l 2 + v b cos k 1 l 1 cos k 2 l 2 B 4 = Z 1 v b sin k 1 l 1 ( 1 Z 3 cos k 2 l 2 sin k 3 l 3 + 1 Z 2 sin k 2 l 2 cos k 3 l 3 ) +                     v b cos k 1 l 1 ( cos k 2 l 2 cos k 3 l 3 Z 2 Z 3 sin k 2 l 2 sin k 3 l 3 ) B 5 = v b ( 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 ) v b ( 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 )
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 = ρ c s i .
The vibration velocity and stress distribution of each part can be calculated as follows:
v 1 ( x 1 ) = v b cos k 1 x 1 v 2 ( x 2 ) = Z 1 Z 2 v b sin k 1 l 1 sin k 2 x 2 + v b cos k 1 l 1 cos k 2 x 2 v 3 ( x 3 ) = v b sin k 3 x 3 ( Z 1 Z 3 sin k 1 l 1 cos k 2 l 2 Z 2 Z 3 cos k 1 l 1 sin k 2 l 2 ) +                                   v b cos k 3 x 3 ( Z 1 Z 2 sin k 1 l 1 sin k 2 l 2 + cos k 1 l 1 cos k 2 l 2 ) v 4 ( x 4 ) = sin k 4 x 4 [ Z 1 Z 4 v b sin k 1 l 1 ( Z 3 Z 2 sin k 2 l 2 sin k 3 l 3 cos k 2 l 2 cos k 3 l 3 )                                   1 Z 4 v b cos k 1 l 1 ( Z 2 sin k 2 l 2 cos k 3 l 3 + Z 3 cos k 2 l 2 sin k 3 l 3 ) ] + cos k 4 x 4                                   [ v b cos k 1 l 1 ( cos k 2 l 2 cos k 3 l 3 Z 2 Z 3 sin k 2 l 2 sin k 3 l 3 ) Z 1 v b sin k 1 l 1                                   ( 1 Z 3 cos k 2 l 2 sin k 3 l 3 + 1 Z 2 sin k 2 l 2 cos k 3 l 3 ) ] v 5 ( x 5 ) = v b ( cos k 5 l 5 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 ) ]
F 1 ( x 1 ) = j Z 1 v b sin k 1 l 1 F 2 ( x 2 ) = j Z 2 v b ( Z 1 Z 2 cos k 2 x 2 sin k 1 l 1 + sin k 2 x 2 cos k 1 l 1 ) F 3 ( x 3 ) = j Z 3 v b [ cos k 3 x 3 ( Z 1 Z 3 sin k 1 l 1 cos k 2 l 2 + Z 2 Z 3 cos k 1 l 1 sin k 2 l 2 ) +                                   sin k 3 x 3 ( Z 1 Z 2 sin k 1 l 1 sin k 2 l 2 + cos k 1 l 1 cos k 2 l 2 ) ] F 4 ( x 4 ) = j Z 4 v b cos k 4 x 4 [ Z 1 Z 4 sin k 1 l 1 ( Z 3 Z 2 sin k 2 l 2 sin k 3 l 3 cos k 2 l 2 cos k 3 l 3 ) +                                   1 Z 2 cos k 1 l 1 ( Z 2 sin k 2 l 2 cos k 3 l 3 + Z 3 cos k 2 l 2 sin k 3 l 3 ) ] j Z 4 v b sin k 4 x 4                                   [ Z 1 sin k 1 l 1 ( 1 Z 3 cos k 2 l 2 + 1 Z 2 sin k 2 l 2 cos k 3 l 3 ) cos k 1 l 1 ( cos k 2 l 2 cos k 3 l 3                                   Z 2 Z 3 sin k 2 l 2 sin k 3 l 3 ) ] F 5 ( x 5 ) = j Z 5 v b ( cos k 5 l 5 sin k 5 l 5 cos k 5 x 5 + sin 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 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 ) ]
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:
Z 1 Z 2 tan k 1 l 1 tan k 2 l 2 + Z 1 Z 3 tan k 1 l 1 tan k 3 l 3 + Z 1 Z 4 tan k 1 l 1 tan k 4 l 4 + Z 1 Z 5 tan k 1 l 1 tan k 5 l 5 + Z 2 Z 3 tan k 2 l 2 tan k 3 l 3 + Z 2 Z 4 tan k 2 l 2 tan k 4 l 4 + Z 2 Z 5 tan k 2 l 2 tan k 5 l 5 + Z 3 Z 4 tan k 3 l 3 tan k 4 l 4 + Z 3 Z 5 tan k 3 l 3 tan k 5 l 5 + Z 4 Z 5 tan k 4 l 4 tan k 5 l 5 Z 1 Z 2 Z 3 Z 4 tan k 1 l 1 tan k 2 l 2 tan k 3 l 3 tan k 4 l 4 Z 1 Z 2 Z 3 Z 5 tan k 1 l 1 tan k 2 l 2 tan k 3 l 3 tan k 5 l 5 Z 1 Z 2 Z 4 Z 5 tan k 1 l 1 tan k 2 l 2 tan k 4 l 4 tan k 5 l 5 Z 1 Z 3 Z 4 Z 5 tan k 1 l 1 tan k 3 l 3 tan k 4 l 4 tan k 5 l 5 Z 2 Z 3 Z 4 Z 5 tan k 2 l 2 tan k 3 l 3 tan k 4 l 4 tan k 5 l 5 = 1
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:
Z 1 Z 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 longitudinal– torsional 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:
tan ( k 8 l 8 ) = k 6 tan k 6 l 6 + β tan k 6 l 6 tan k 7 l 7 k 7 tan k 7 l 7 k 6 tan k 6 l 6 tan k 7 l 7 + β tan k 7 l 7 + k 7

2.2. 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 longitudinal–torsional 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 2. Where Ⅰ, Ⅱ indicates the wave propagation from both sides of the node, the arrow indicates the positive direction of the transducer, and the curves indicate the wave propagation schematic at different amplitudes.
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 longitudinal–torsional 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.
Figure 5 illustrates the front cover plate of the longitudinal–torsional transducer. 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.
The longitudinal force and tangential force can be determined by the following equation:
{ F l = F · cos φ F t = F · sin φ
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:
M = s r f d s
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 d s = 2 π r d r is the element area. Therefore, the tangential force per unit area can be obtained as:
f = F sin ( φ ) π ( r 2 2 r 1 2 )
Equation (16) can be simplified as follows:
M = F sin ( φ ) 2 ( r 2 3 r 1 3 ) 3 ( r 2 2 r 2 1 )
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 Figure 6 and Figure 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.
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 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.

3. 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.

3.1. 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:
u = u a ( x , y , z ) e [ ω t ( x , y , z ) ]
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, d p , can be expressed as the integral of the amplitude on the panel:
d p = j κ ρ 0 c 0 2 π h ( x , y , z ) d Q 0 e j [ ω t k h ( x , y , z ) α ( x , y , z ) ]
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 0 2 u a .
Total sound pressure can be obtained as follows:
p = S j k ρ 0 c 0 2 π h ( x , y , z ) e j [ ω t k h ( x , y , z ) α ( x , y , z ) ] d Q 0
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:
d p = j k ρ 0 c 0 2 π r u a d S e j ( ω t k h )
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:
p = d p = S j k ρ 0 c 0 2 π r u a e j ( ω t k r ) d S = j ω ρ 0 u a 2 π r e j ( ω t k r ) 0 a ρ d ρ 0 2 π e j k ρ sin θ cos ξ d ξ
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],
J 0 ( x ) = 1 2 π 0 2 π e j x cos ξ d ξ , x J 0 ( x ) d x = x J 1 ( x )
The integral of Equation (24) can be expressed as the following equation:
p = j ω ρ 0 u a a 2 r 2 [ 2 J 1 ( k a sin θ ) k a sin θ ] e j ( ω t k r )
The radial velocity of the mass, meanwhile, can be expressed as:
ν r = 1 j ω ρ 0 p r = 1 ρ 0 c 0 ( 1 + 1 j k r ) p
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.

3.2. 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],
J n ( x ) = k = 0 ( 1 ) k ( x 2 ) 2 k + n k ! ( n + k ) !
Therefore, the directivity of the piston sound source can be expressed as follows:
D ( θ ) = ( p a ) θ ( p a ) θ = 0 = | 2 J 1 ( k a sin θ ) k a sin θ |
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.

4. Discussion

4.1. 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.

4.2. 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.

4.3. 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.

5. Conclusions

This paper has proposed a longitudinal–torsional transducer for honeycomb aramid material during ultrasonic vibration assisted milling (UVAM). The directivity of transducer was analyzed via a simulation and an experiment, providing a new evidence for the processing of honeycomb aramid material and wide application. The major conclusions were listed as the following:
  • The longitudinal wave can be converted to a transverse wave via the longitudinal–torsional transducer, thus realizing longitudinal–torsional conversion.
  • During UVAM, the maximum longitudinal vibration can reach 19 μm, and the maximum torsional amplitude at the top of the disk milling cutter can reach 9 μm.
  • Compared with traditional milling, the cutting force along three directions during UVAM with a longitudinal–torsional transducer is significantly reduced. Moreover, the directivity is stronger and the cutting depth is greater under the same cutting speed. This shows that UVAM technology can improve the stability of the milling process.

Author Contributions

X.W. and M.Z. conceived and designed the experiments; Z.M., T.M. and X.L. performed the experiments; Z.M. and X.L. analyzed the data; X.W. and M.Z. contributed reagents/materials/analysis tools; Z.M. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brehl, D.E.; Dow, T. Review of vibration-assisted machining. Precis. Eng. 2008, 32, 153–172. [Google Scholar] [CrossRef]
  2. Kumar, M.N.; Subbu, S.K.; Krishna, P.V.; Venugopal, A. Vibration assisted conventional and advanced machining: A review. Procedia Eng. 2014, 97, 1577–1586. [Google Scholar] [CrossRef] [Green Version]
  3. Harkness, P.; Lucas, M.; Cardoni, A. Coupling and degenerating modes in longitudinal–torsional step horns. Ultrasonics 2012, 52, 980–988. [Google Scholar] [CrossRef]
  4. Celaya, A.; Lopez de Lacalle, L.N.; Campa, F.J.; Lamikiz, A. Ultrasonic Assisted Turning of mild steels. Int. J. Mater. Prod. Technol. 2010, 37, 60. [Google Scholar] [CrossRef]
  5. Liu, S.; Shan, X.; Cao, W.; Yang, Y.; Xie, T. A longitudinal-torsional composite ultrasonic vibrator with thread grooves. Ceram. Int. 2017, 43, S214–S220. [Google Scholar] [CrossRef]
  6. Radi Najafabadi, Y.; Nosouhi, R.; Sajajid, F. Design and fabrication of a longitudinal-torsional ultrasonic transducer. J. Mod. Process. Manuf. Prod. 2017, 6, 25–31. [Google Scholar]
  7. Guofu, G.; Ziwen, X.; Zhaojie, Y.; Xiang, D.; Bo, Z. Influence of longitudinal-torsional ultrasonic-assisted vibration on micro-hole drilling Ti-6Al-4V. Chin. J. Aeronaut. 2021, 34, 247–260. [Google Scholar]
  8. Songmei, Y.; Zhixiang, T.; Qi, W.; Heng, S. Design of Longitudinal Torsional Ultrasonic Transducer and Its Performance Test. J. Mech. Eng. 2019, 55, 139–148. [Google Scholar]
  9. Qiaoli, Z.; Jianfu, Z.; Pingfa, F.; Dingwen, Y.; Zhijun, W. Characteristics of the longitudinal-torsional vibration of an ultrasonic horn with slanting slots. J. Vib. Shock 2019, 38, 58–64. [Google Scholar]
  10. Wu, C.; Chen, S.; Cheng, K.; Ding, H.; Xiao, C. Innovative design and analysis of a longitudinal-torsional transducer with the shared node plane applied for ultrasonic assisted milling. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2019, 233, 4128–4139. [Google Scholar] [CrossRef]
  11. Al-Budairi, H.; Lucas, M.; Harkness, P. A design approach for longitudinal–torsional ultrasonic transducers. Sens. Actuators A Phys. 2013, 198, 99–106. [Google Scholar] [CrossRef]
  12. Weishan, C.; Yingxiang, L.; Junkao, L.; Shengjun, S. Working principle and design of a double cylinders type traveling wave ultrasonic motor using composite transducer. J. Harbin Inst. Technol. 2011, 18, 28–32. [Google Scholar]
  13. Wang, L.; Wang, J.-A.; Jin, J.-M.; Yang, L.; Wu, S.-W.; Zhou, C.C. Theoretical modeling, verification, and application study on a novel bending-bending coupled piezoelectric ultrasonic transducer. Mech. Syst. Signal Process. 2022, 168, 108644. [Google Scholar] [CrossRef]
  14. Suzuki, A.; Tsujino, J. Load characteristics of ultrasonic motors with a longitudinal-torsional converter and various nonlinear springs for inducing static pressure. Jpn. J. Appl. Phys. 2002, 41, 3267. [Google Scholar] [CrossRef]
  15. Shen, X.-H.; Zhang, J.; Xing, D.X.; Zhao, Y. A study of surface roughness variation in ultrasonic vibration-assisted milling. Int. J. Adv. Manuf. Technol. 2012, 58, 553–561. [Google Scholar] [CrossRef]
  16. Chen, Y.; Su, H.; He, J.; Qian, N.; Gu, J.; Xu, J.; Ding, K. The Effect of Torsional Vibration in Longitudinal–Torsional Coupled Ultrasonic Vibration-Assisted Grinding of Silicon Carbide Ceramics. Materials 2021, 14, 688. [Google Scholar] [CrossRef]
  17. Cleary, R.; Wallace, R.; Simpson, H.; Kontorinis, G.; Lucas, M. A longitudinal-torsional mode ultrasonic needle for deep penetration into bone. Ultrasonics 2022, 124, 106756. [Google Scholar] [CrossRef]
  18. Wang, Y.; Cai, Y.; Lee, Y.-S. Design and Development of Longitudinal and Torsional Ultrasonic Vibration-Assisted Needle Insertion Device for Medical Applications. Comput.-Aided Des. Appl. 2021, 19, 797–811. [Google Scholar] [CrossRef]
  19. Numanoğlu, H.M.; Akgöz, B.; Civalek, Ö. On dynamic analysis of nanorods. Int. J. Eng. Sci. 2018, 130, 33–50. [Google Scholar] [CrossRef]
  20. Akbaş, Ş.D.; Ersoy, H.; Akgöz, B.; Civalek, Ö. Dynamic analysis of a fiber-reinforced composite beam under a moving load by the Ritz method. Mathematics 2021, 9, 1048. [Google Scholar] [CrossRef]
  21. Kim, J.O.; Kwon, O.S. Vibration characteristics of piezoelectric torsional transducers. J. Sound Vib. 2003, 264, 453–473. [Google Scholar] [CrossRef]
  22. Mishiro, S. Torsional Vibration Apparatus. US Patent US4652786A, 1987. [Google Scholar]
  23. Harada, T.; Ishikawa, N.; Kanda, T.; Suzumori, K.; Yamada, Y.; Sotowa, K.-I. Droplet generation using a torsional Langevin-type transducer and a micropore plate. Sens. Actuators A Phys. 2009, 155, 168–174. [Google Scholar] [CrossRef]
  24. Lianjun, S.; Wenhe, L.; Zheng, K.; Wei, T.; Jinshan, L.; Jindan, F. Stability analysis of robotic longitudinal-torsional composite ultrasonic milling. Chin. J. Aeronaut. 2022, 35, 249–264. [Google Scholar]
  25. Jun, T.; Bo, Z. Stability analysis of the separated longitudinal-torsional composite ultrasonic milling. Acta Armamentarii 2015, 36, 1318. [Google Scholar]
  26. Zhou, X.; Wang, L.; Yu, D.; Zhang, C. Dynamic effective equivalent stiffness analysis on the periodical honeycomb reinforced composite laminated structure filled with viscoelastic damping material. Compos. Struct. 2018, 193, 306–320. [Google Scholar] [CrossRef]
  27. Suárez, A.; Veiga, F.; de Lacalle, L.N.L.; Polvorosa, R.; Lutze, S.; Wretland, A. Effects of ultrasonics-assisted face milling on surface integrity and fatigue life of Ni-Alloy 718. J. Mater. Eng. Perform. 2016, 25, 5076–5086. [Google Scholar] [CrossRef]
  28. Karafi, M.R.; Hojjat, Y.; Sassani, F.; Ghodsi, M. A novel magnetostrictive torsional resonant transducer. Sens. Actuators A Phys. 2013, 195, 71–78. [Google Scholar] [CrossRef]
  29. Ziomek, L.J. Fundamentals of Acoustic Field Theory and Space-Time Signal Processing; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
Figure 1. Structural diagram of longitudinal–torsional transducer.
Figure 1. Structural diagram of longitudinal–torsional transducer.
Micromachines 13 02154 g001
Figure 2. Waveform diagram of nodal plane.
Figure 2. Waveform diagram of nodal plane.
Micromachines 13 02154 g002
Figure 3. Schematic diagram of vibration propagation.
Figure 3. Schematic diagram of vibration propagation.
Micromachines 13 02154 g003
Figure 4. Schematic diagram of spiral structure.
Figure 4. Schematic diagram of spiral structure.
Micromachines 13 02154 g004
Figure 5. Front cover plate of longitudinal–torsional transducer.
Figure 5. Front cover plate of longitudinal–torsional transducer.
Micromachines 13 02154 g005
Figure 6. Effect of spiral groove parameters on longitudinal vibration.
Figure 6. Effect of spiral groove parameters on longitudinal vibration.
Micromachines 13 02154 g006
Figure 7. Effect of spiral groove parameters on torsional vibration.
Figure 7. Effect of spiral groove parameters on torsional vibration.
Micromachines 13 02154 g007
Figure 8. Modal velocity vector diagram of longitudinal–torsional transducer.
Figure 8. Modal velocity vector diagram of longitudinal–torsional transducer.
Micromachines 13 02154 g008
Figure 9. Longitudinal–torsional amplitude.
Figure 9. Longitudinal–torsional amplitude.
Micromachines 13 02154 g009
Figure 10. Distribution of radiated sound field. (a) round piston sound source, (b) rectangular piston sound source.
Figure 10. Distribution of radiated sound field. (a) round piston sound source, (b) rectangular piston sound source.
Micromachines 13 02154 g010aMicromachines 13 02154 g010b
Figure 11. Directivity pattern of the piston transducer with different radii of 30 mm (a), 40 mm (b) and 50 mm (c).
Figure 11. Directivity pattern of the piston transducer with different radii of 30 mm (a), 40 mm (b) and 50 mm (c).
Micromachines 13 02154 g011aMicromachines 13 02154 g011b
Figure 12. Experimental setup.
Figure 12. Experimental setup.
Micromachines 13 02154 g012
Figure 13. Schematic diagram of the cutting force system.
Figure 13. Schematic diagram of the cutting force system.
Micromachines 13 02154 g013
Figure 14. System for measurement cutting force.
Figure 14. System for measurement cutting force.
Micromachines 13 02154 g014
Figure 15. Cutting force under different processing methods.
Figure 15. Cutting force under different processing methods.
Micromachines 13 02154 g015
Figure 16. Directivity pattern under different processing methods.
Figure 16. Directivity pattern under different processing methods.
Micromachines 13 02154 g016aMicromachines 13 02154 g016b
Table 1. Specific parameters of transducer structure.
Table 1. Specific parameters of transducer structure.
Number of GroovesLengthDepthWidthl1Angle
622 mm1.5 mm2 mm10 mm45°
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, M.; Ma, Z.; Wang, X.; Meng, T.; Li, X. Design of Longitudinal–Torsional Transducer and Directivity Analysis during Ultrasonic Vibration-Assisted Milling of Honeycomb Aramid Material. Micromachines 2022, 13, 2154. https://doi.org/10.3390/mi13122154

AMA Style

Zhang M, Ma Z, Wang X, Meng T, Li X. Design of Longitudinal–Torsional Transducer and Directivity Analysis during Ultrasonic Vibration-Assisted Milling of Honeycomb Aramid Material. Micromachines. 2022; 13(12):2154. https://doi.org/10.3390/mi13122154

Chicago/Turabian Style

Zhang, Mingxing, Zuotian Ma, Xiaodong Wang, Ting Meng, and Xiangqun Li. 2022. "Design of Longitudinal–Torsional Transducer and Directivity Analysis during Ultrasonic Vibration-Assisted Milling of Honeycomb Aramid Material" Micromachines 13, no. 12: 2154. https://doi.org/10.3390/mi13122154

APA Style

Zhang, M., Ma, Z., Wang, X., Meng, T., & Li, X. (2022). Design of Longitudinal–Torsional Transducer and Directivity Analysis during Ultrasonic Vibration-Assisted Milling of Honeycomb Aramid Material. Micromachines, 13(12), 2154. https://doi.org/10.3390/mi13122154

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop