A New Approach on Vibrating Horns Design

An optimization method of the vibrating horns is presented considering the smallest action principle and the attached cutting tool mass. The model is based on Webster’s wave propagation equation and as an objective function the minimization of the volume in structural equilibrium conditions was considered. The considered input parameters were working frequency, maximum crosssectional area, magnification coefficient, and the attached mass. At the end of the study, a new shape function of the horn’s cross section is obtained. The particularity of the new obtained shape is given by the nodal point position that is the same with the position of the maximum cross-sectional area. The obtained horn was analyzed from the modal point of view using theoretical and experimental methods. As theoretical methods, both the state-space method and the finite element method were used. An experimental setup for frequency response function determination was developed using a random input signal. The verification of the magnitude value was done considering a harmonic steady-state signal. The recorded values were compared with the predicted values. The numerical simulations and tests support the validity of the assumptions used in the horns optimization design.


Introduction
Ultrasonic horns are axisymmetric devices that provide, with high efficiency, acoustic energy developed by a transducer that is working based on the magnetostriction effect, to a tool.These devices can be themselves tools or can be used with a tool attached at the small end, and thus they are tied with high-power applications of ultrasound, like cutting [1][2][3], drilling [4][5][6], turning [7], welding [8,9], and so forth, or for machining different materials as carbon fiber reinforced composites [10][11][12].
From a physical point of view, the work of such a system consists of an initial signal, emitted by a transducer at a given frequency, with an amplitude  0 , which is mechanically amplified by the horn, with an imposed magnitude coefficient .The whole system is working in a resonance regime at an imposed frequency.The mechanical magnification of the input signal amplitude is generated by the horn varying cross section from an initial value  0 , in the origin of the reference system, to an end one value   (Figure 1).Generally, the cross section variation is given by different mathematical functions (exponential [13], linear [14], tapered (or stepped) [15], catenoidal [14], Bézier [16], Gaussian, etc.).Studies concerning the wave propagation and basic concepts of horns design are presented in different papers [17][18][19][20][21].
A problem that is associated with the horns design refers to the shape optimization.The optimization procedure is connected with design parameters like frequency working values [21], amplitude of the signal [22,23], load transferred in the manufacturing area [24], the objective function [25,26], the impact of different design variables on the objective function [27], combined signals transmitted in the horn [28], and so forth.
In the present paper, we present an approach based on the principle of the least action that takes into consideration as an objective function the volume minimization.The considered input parameters are the working frequency, the maximum cross-sectional area, the magnification coefficient, and the attached mass.

Mathematical Model
In case of ultrasonic horns, a planar stress waves state is assumed.The general equation that describes one-dimensional 2 Shock and Vibration propagation, in the longitudinal direction, of ultrasonic waves along a rigid horn with a variable cross-sectional area is Webster's horn equation [29]: where (, ) is the signal,   = () represents the crosssectional area of the considered horn at the distance  from the reference (point 0), and  is the wave velocity in the horn's material.
One of the widely used techniques to find the explicit analytical solutions of the linear partial differential equations is the method of variables separation [30].The method assumes that the solution of (1) can be written as a product of two functions, one in space domain () =   and the other one in time domain () =   :  (, ) =  ()  () =     . ( Introducing solution (2) in Webster's equation (1), the new shape is obtained: where the following notations were done: This way, in (3), there were separate variables on each side of the equation.
The congruency of both sides, left and right, of (3) is achieved only if the two functions are equal with the same separation constant .It is considered as a constant value  = − 2 , where  is the angular wave frequency (rad/s).
Thus, from (3), we obtained the following ordinary differential equations: where  = / is the wave number.
The general solution of the time domain equation from ( 4) is given by with the initial conditions where  0 is the initial amplitude of the ultrasonic wave.Thus, solution (5) becomes and the total solution (2) becomes In case of a steady-state mode, considering the assumption of an infinitesimal motion in a beam made of nondissipative medium and with a variable cross-sectional area, the second equation from (4), which represents the plane wave motion, can be rewritten as It is considered that the maximum diameter of the horn should not exceed one-quarter of the wavelength  and the maximum cross section is denoted as   .Thus, we can write the condition that, on the whole length of the horn, the cross section value has to be smaller than the maximum cross section: Considering relation (10), one can write the following mathematical relation: where () is an unspecified function of the variable .In the following, the notations () =   and () =   will be used.
An optimization problem of the horn's shape refers to minimizing the horn's volume for a given initial set of design data.Based on the principle of the least action, one can write the action relation as where  is the kinetic energy and  represents the potential energy.
Considering the wave shape as for a beam of length  with attached mass at one end, the case of an ultrasonic horn with the attached tool, the total energies of the system horn tool, mentioned in (12), are given by where   is the specific mass of the horn,   is the attached mass (toll mass),  ℎ is Young's modulus of the horn material,   is Young's modulus of the tool material,   is the volume of the attached tool,  ℎ represents the horn material density, u = /, and   = /.Considering a period, action (12) becomes with Based on [31,32], the following functional is defined: where  1 and   = () are Lagrange multipliers and  is the volume of the horn, given by The optimization problem has a solution for the imposed limit values if the functional Ι given by (17) has to be stationary related to the variations of   and to the design parameter   .Thus, based on (18), it is necessary for Euler's equations [33] to be satisfied: where    =   /,    =   /, and    =   /.Considering the functional given by ( 17) and (19), one can obtain the following differential equations: (a) From the first equation in (19) results which when divided by  ℎ becomes Considering the velocity of the longitudinal wave in the horn as   = √ ℎ / ℎ and the wave number  = /  , (21) can be rewritten as (b) From the second equation, one obtains which when divided by  ℎ becomes (c) From the third equation results, it is found that Analyzing the previous equations, one can obtain the following remarks about the optimal design of the ultrasonic horn: (a) Equation ( 21) is identical to the longitudinal wave propagation equation along a horn and it is available along the whole horn.
(b) Equation ( 25) leads to two other conditions: Conditions (26) involve the following aspect: the function   is different from zero on the whole length of the horn except for one point, where   =   .
Consequently, for (25) to be satisfied for any point in the range 0 ≤  ≤ , it is necessary to have the first condition from ( 26)   = 0 in the abovementioned range.Moreover, all Lagrange's coefficients have to be constants.

The Solution of the Differential Equation
The solution of ( 28) can be found considering the initial functional conditions of the ultrasonic horn.The following are considered as initial conditions: (a) The length of the horn has to be equal to a wavelength.(b) To achieve the assumption that the propagated waves are only longitudinal ones, a maximum value is considered for the magnification coefficient  = 5 [17].(c) The initial value of the signal amplitude, at the connection point between the transducer and the horn ( = 0), is  0 equal to the emitted signal by the transducer: (d) Considering the initial signal and the magnification coefficient, the signal at the end of the horn (for  = ) is given by (e) In ( 30), the negative value (−) is considered in order to obtain a nodal point in a position along the horn length.(f) Between the two ends of the horn, there is no other extreme of the signal except for the maximum developed at the end of the horn, in the connection point with the tool.(g) The following condition is considered [32]: where  is the mass on the unit length, at the end of the horn,   = is the strain at the end of the horn,  is the attached mass, and  = is the signal at the end of the horn, given by (30).
equation ( 28) can be rewritten as and integrating it with the derivative method, one can obtain the following equation [34]: Dividing ( 34) by   ̸ = 0 results in the differential equation with the solution Introducing solution (36) in (33), after some calculations, the following is obtained: Considering (38) in solution (36) results in The constants  and  2 can be found considering the initial conditions mentioned above, for  = 0. Combining (38) and (30), the following relation is obtained: which leads to the following relation for the two values of the constant  2 as a function of : Considering ( 41) and ( 40), for both values  1,2 2 , the same dependence is obtained between the constants  and  2 : Introducing the constant , given by (42), in solution (39) leads to the following relation of the amplitude (): The constant  2 can be found considering conditions ( 39) and (40): where  ℎ is the density of the horn's material,  ℎ is the cross section at the end of the horn,   is the cross section of the tool (  =   ),   is the tool's length,   is the density of the tool's material, and  is the gain of the initial signal  0 .
From (44) results Considering dependence (42), one can define the function of the space dependent amplitude component of the signal along the horn: In the following, we use the notation and the signal function (46) becomes As it is known, an important aspect in horns design is to find the position of the neutral point   .The neutral point defines the position of the cross section where the signal is zero.Thus, the neutral point can be easily found considering relation (48).Equalizing function (48) with zero, it is obtained that ) . (49)

Cross Section Shape Definition
From relation (48), we obtained the first two derivatives of the space component of the signal: Equation ( 9) can be rewritten as Considering relations (50) in (51) results in the new formula Defining the function () =   as (52) becomes Integrating (54) results in a constant, denoted by : which depends on the initial geometrical conditions.At the end of the horn, for  = , the cross section becomes   =   and it is obtained that where   is the value of the function () =   , defined by (53), for  = .
Introducing the value of  obtained by (56) in (55) results in or and the cross section at the distance  can be calculated as An interesting design aspect is represented by the value and position of the maximum cross section.Thus, we considered the condition that the first derivative of function ( 59) is equal to zero.This condition leads to the relation with similar solution to (49).
It can be concluded that the position of the maximum cross section is in the same place as the nodal point.
The end cross section can be found as and introducing (61) in (59) results in the relation of the cross section value in a point situated at the distance  from the system origin, as a function of the maximum cross section: Taking into consideration relation (62), a script in Matlab was made to generate the shape of the horn made of steel.The input parameters were the resonance frequency   = 19900 Hz, the wave velocity in steel  = 50550 m/s, the horn's material density  ℎ = 7850 kg/m 3 , the end radius of the horn (connection surface with the tool)  ℎ = 0.005 m, the magnification coefficient  = 5, and tool radius equal to the end horn radius (  =  ℎ ).The tool was considered to be made of steel.
Based on the script written in Matlab, the shape presented in Figure 2 was obtained, and in Figure 3 the physical model of the horn is presented.
Based on the same script, the following design data were found: the radius in the origin of the system ( = 0)  0 = 29.256mm, the largest radius   = 31.721mm, and the position of the nodal point  nod = 27.734mm.
Considering relation (48), using a script written in Matlab, we obtained the magnitude variation of the spatial components () of the considered signal (13) (Figure 4).
As it can be seen from Figure 4, the nodal point is at the coordinate found by relation (49).

State-Space Approach of the Optimized Model.
The modal behavior can be described in many ways.Two of them are the state-space approach and transfer function representation.The differences between these two methods were described in [35,36].In the present paper, we considered as a theoretical modal analysis the state-space approach, and for the experimental modal analysis the FRF representation was made.
The state-space method is often used in dynamic system analysis.This method can be applied in the case of lumped systems and it is based on the concept of state vector that consists of two components: the vector of displacement and the vector of the velocities.Practically, the "" second-order coupled equations that describe the dynamic behavior of a lumped system are transformed in a set of 2 first-order coupled equations [37][38][39].
The ultrasonic horn can be considered as a linear timeinvariant (LTI) system that can be approximate with a lumped mass system.Thus, the obtained horn (Figure 3) is divided into "" elements with the same length   , where   = /, with  being the total length of the horn equal to a quarter of the wavelength.Considering the case of the free-free system, as is in the experimental modal analysis, the horn can be modeled as a lumped system made of "" masses connected with springs (Figure 5).
We considered 200 masses with a length of   = 6.344 ⋅ 10 −4 m in each part.The obtained components are, taking into consideration the shape (Figure 3), small frusta of cone elements.
The components of the modal model are established considering the following assumptions: (a) The stiffness of the springs was considered to be equal to the equivalent stiffness of each element based on the relation where   represents the average area of the region of length   : with both   and  +1 calculated using relation (62).
(b) The mass   of each small frustum of cone element was found considering the classical relation of the density using for volume the areas   and  +1 and the length   .
Considering the d' Alembert principle, the motion equations for the system from Figure 5 are given by Shock and Vibration 7 System (65) can be rewritten in matrix form as follows: or from which it is obtained that As the state vector, we considered the matrix vector which, by derivation and considering (68), leads to the equation where I is the unit matrix and 0 is the zeros matrix, both of  ×  type.The matrix A is defined as the state matrix or the transmission matrix and is a property of the system being  defined by the structure of the system.The eigenvalues of the state matrix can be found based on the relation that is developed in an equation with unknowns .
If, in relation (71), the value  =   is considered, where  = √ −1 and   is the natural frequency measured in rad/s, then the so-called characteristic equation is obtained, which gives the natural frequencies.
Thus, considering the equivalent system from Figure 5 with the associated motion equations (66), (67), and (70), one can find the natural frequencies of the system.The number and the precision of these frequencies depend on the number of the considered masses   ( = 1, ). Figure 6 presents the frequency values obtained for a model with 200 masses considering the method of eigenvalues of the transition matrix.
Figure 6 also presents some of the obtained natural frequencies.The first natural frequency found was   = 19,837 Hz, which is close to the designed frequency of   = 19.9kHz.

Finite Element Model.
In the next step, the finite element model was considered.The used specialized software was ANSYS.The 3D model (Figure 7) was made based on (59).
For good convergence, we considered a mapped mesh consisting of first-order solid hexaelements with 1 mm length.We obtained 87,076 nodes and 75,254 elements.Mode extraction is carried out in the frequency range 18-24 kHz using Block Lanczos option with double precision.Two pure longitudinal modes were obtained with the natural frequencies of  = 19,817 Hz (Figure 8) and  = 23,151 Hz (Figure 9).

Experimental Modal Analysis
The aim of the experimental modal analysis is to verify whether the considered imposed working frequency of 19.9 kHz is one of the natural frequencies of the considered horn and whether the natural frequencies found by statespace approach are close to the real natural frequencies.
The method with random signal excitation given by an electrodynamic shaker (Figure 10) was considered [40].The used equipment was a complete testing setup (Brüel & Kjaer) made of an integrated platform, PULSE 12, with a force transducer (type 8230) and an accelerometer (type 4397).The excitation signal was generated by the PULSE 12 platform and was amplified using the power amplifier type 2718 and then transmitted to the tested horn by the shaker type 4810 (Figure 11).
The random generated signal (the input signal) has a frequency span up to 25.6 kHz with an excitation time of 500 ms (Figure 12).The Fourier spectrum of the input signal is shown in Figure 13.As it is known, the frequency response function estimation from the measured data represents the main step of the experimental modal analysis.The analysis is based on nonparametric estimators.In the present study, we considered the often used estimator  1 .The obtained frequency response function, for the estimator  1 , is presented in Figure 14.
The range of the measured frequencies was up to 20,000 Hz and the values are presented in Figure 15.
At the same time, based on the experimental modal analysis, we obtained the damping factors  for each natural frequency found.The distribution of these values is presented in Figure 16.
Considering the facilities of the Matlab software, a power function was defined, which can approximate the damping factor for different values of the frequency :  () = 321.9−0.7427 − 0.07775.
(72) Function ( 72) can be used to study the influence of the damping in case of horn vibration at different frequencies.

Testing the Designed Horn
The test was designed to check whether the magnification of the designed horn is the initial considered  = 5 for the resonance frequency of 19.9 kHz.The used setup is made of the same components as the setup used in the experimental modal analysis with the exception of the accelerometers.For testing, two Brüel & Kjaer accelerometers type 4517-002 were used.The signal was generated by a magnetostrictive transducer for a harmonic signal () =  sin  with a frequency  = 19.9kHz.The input signal, represented in time domain, is shown in Figure 17, and the output signal (response) also in time domain is designed in Figure 18.At the same time, for frequency control, we performed a check in frequency domain for both signals: the input (Figure 19) and the output (Figure 20).
Considering the values of the input and output signals (Figures 17 and 18), the magnification obtained by experiment is given by the following relation: which represents an error of about 8.2% between the considered value  = 5 and the obtained value  real .

Conclusions
Based on the principle of the smallest action, an optimization of the shape of an ultrasonic horn was performed considering Webster's equation.The design was made considering a working resonance frequency of   = 19.9kHz and a magnification coefficient  = 5, and the main assumption was that, through the horn and all attached coupling devices, only longitudinal plane waves propagate.The considered boundary conditions were the classical ones referring to the amplitude of the signal, relations (29) and (30), to which a new one was added, relation (31), which highlights the influence of the tool mass, added at the end of the horn.This condition (31) is a particularity of the study.
The optimization procedure led to a new shape function (62), with the space dependent signal component function (48) and the position of the neutral point   (49) being determined at the same time.Comparing the obtained shape with the classical shape, it can be seen that it is a different one, with the largest diameter in the neutral point.In the case of the other horns' shape, the largest diameter is at one end.The natural frequencies were determined both in simulation (theoretically) and experimentally.In the case of theoretical modal analysis, the state-space representation was used.The horn was divided into 200 masses connected by springs (Figure 5) and the eigenvalues of the transmission matrix of system (71) were calculated.The first obtained value of the natural frequency was the imposed resonant frequency  1, = 19,982 Hz.
The second considered theoretical approach in theoretical modal analysis was the finite element method.A model with 75,254 elements was designed, and based on ANSYS software facilities, two natural frequencies were found around the designed resonance frequency (Figures 8 and 9), both for the longitudinal mode shape.In order to provide clarity to the results, a run was made in double precision mode.
The experimental modal analysis offered the possibility of finding a set of natural frequencies (Figure 16) and the corresponding values of the damping factor.Considering the facilities of the Matlab tool, a function of damping factor variation was obtained according to the frequencies values (72) useful in the dynamic analysis.As it can be seen, a natural frequency around the resonance designed frequency was found.
At the end, the horn was tested to check whether the magnification coefficient is the same as the designed one.A harmonic signal with the frequency value equal to the resonance frequency was considered as the input, and the output signal was measured.The found magnification coefficient is close to the designed one (72).
In the frame of the test, the attached tool was not considered.It is important for the assembly transducer-horntool to work in resonance regime and each part to work within the parameters for which it was designed.
The aim of the test was to check whether the considered horn, designed in the proposed set of boundary conditions ( 29), (30), and (31), generates the initial designed mechanical magnification of the input signal.
The influence of the attached tool mass is the particular design element that was taken into consideration compared with the other methods of optimization and design.
Considering the abovementioned aspects, it can be concluded that the numerical simulations and tests that were done support the validity of the assumptions used in the horn's optimization design.
The proposed model offers as an advantage a more convenient placement of the nodal point that is situated in the same place as the largest cross section.
The described design and optimization method offer the possibility of obtaining the shape based on the main working parameters.

Figure 2 :
Figure 2: The half-horn shapes along the longitudinal axis.

Figure 3 :Figure 4 :
Figure 3: The physical model of the horn.

Figure 5 :
Figure 5: The lumped model of the horn.

6 Figure 6 :
Figure 6: The natural frequencies obtained for the lumped model.

Figure 7 :
Figure 7: The finite element model of the horn.

Freq. = 19817 Figure 8 :
Figure 8: The mode shape for the natural frequency of  = 19,817 Hz.

Figure 9 :
Figure 9: The mode shape for the natural frequency of  = 23,151 Hz.

Figure 18 :
Figure 18: The output signal in time domain.

Figure 20 :
Figure 20: The output signal in frequency domain.