Shock and Vibration
Volume 2017 (2017), Article ID 8532021, 12 pages
https://doi.org/10.1155/2017/8532021

Research Article

A New Approach on Vibrating Horns Design

Maria Violeta Guiman and Ioan Călin Roșca

Department of Mechanical Engineering, Faculty of Mechanical Engineering, Transilvania University of Brașov, 29 Eroilor Av., 500036 Brașov, Romania

Correspondence should be addressed to Ioan Călin Roșca; icrosca@unitbv.ro

Received 5 February 2017; Revised 1 June 2017; Accepted 19 June 2017; Published 7 August 2017

Academic Editor: Tai Thai

Copyright © 2017 Maria Violeta Guiman and Ioan Călin Roșca. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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 cross-sectional 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.

1. 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–3], drilling [4–6], turning [7], welding [8, 9], and so forth, or for machining different materials as carbon fiber reinforced composites [10–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 , 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 , in the origin of the reference system, to an end one value (Figure 1).

Figure 1: The solid horn geometry.

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

2. Mathematical Model

In case of ultrasonic horns, a planar stress waves state is assumed. The general equation that describes one-dimensional 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 cross-sectional 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: , , and .

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 , 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 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, , and .

Considering a period, action (12) becomes

with

Based on [31, 32], the following functional is defined:

where 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 , it is necessary to have the first condition from (26) in the abovementioned range. Moreover, all Lagrange’s coefficients have to be constants.

(a)From (24) and (26), the following is obtained:

or

3. 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 [17].(c)The initial value of the signal amplitude, at the connection point between the transducer and the horn , is 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).

Denoting

equation (28) can be rewritten as

and integrating it with the derivative method, one can obtain the following equation [34]:

Dividing (34) by 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 can be found considering the initial conditions mentioned above, for . Combining (38) and (30), the following relation is obtained:

which leads to the following relation for the two values of the constant as a function of :

Considering (41) and (40), for both values , the same dependence is obtained between the constants and :

Introducing the constant , given by (42), in solution (39) leads to the following relation of the amplitude :

The constant 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 .

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

4. 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  Hz, the wave velocity in steel  m/s, the horn’s material density  kg/m³, the end radius of the horn (connection surface with the tool)  m, the magnification coefficient , 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.

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

Figure 3: The physical model of the horn.

Based on the same script, the following design data were found: the radius in the origin of the system ()  mm, the largest radius  mm, and the position of the nodal point  mm.

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

Figure 4: The signal magnitude along the horn.

As it can be seen from Figure 4, the nodal point is at the coordinate found by relation (49).

5. Theoretical Modal Analysis

5.1. 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 first-order coupled equations [37–39].

The ultrasonic horn can be considered as a linear time-invariant (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).

Figure 5: The lumped model of the horn.

We considered 200 masses with a length of  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 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 and the length .

Considering the d’Alembert principle, the motion equations for the system from Figure 5 are given by

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 is the unit matrix and is the zeros matrix, both of type.

The matrix 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 and is the natural frequency measured in , 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 . Figure 6 presents the frequency values obtained for a model with 200 masses considering the method of eigenvalues of the transition matrix.

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

Figure 6 also presents some of the obtained natural frequencies. The first natural frequency found was  Hz, which is close to the designed frequency of  kHz.

5.2. 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).

Figure 7: The finite element model of the horn.

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  Hz (Figure 8) and  Hz (Figure 9).

Figure 8: The mode shape for the natural frequency of  Hz.

Figure 9: The mode shape for the natural frequency of  Hz.

6. 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 state-space 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 & Kjær) 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).

Figure 10: The modal testing setup.

Figure 11: The detailed connection shaker horn.

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.

Figure 12: The input random signal in time domain.

Figure 13: The Fourier spectrum of the input signal.

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 . The obtained frequency response function, for the estimator , is presented in Figure 14.

Figure 14: The frequency response function of the horn.

The range of the measured frequencies was up to 20,000 Hz and the values are presented in Figure 15.

Figure 15: The frequencies obtained by experimental data.

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.

Figure 16: The damping factor values distribution versus natural frequencies.

Considering the facilities of the Matlab software, a power function was defined, which can approximate the damping factor for different values of the frequency :

Function (72) can be used to study the influence of the damping in case of horn vibration at different frequencies.

7. Testing the Designed Horn

The test was designed to check whether the magnification of the designed horn is the initial considered 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 & Kjær accelerometers type 4517-002 were used. The signal was generated by a magnetostrictive transducer for a harmonic signal with a frequency  kHz. 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.

Figure 17: The input signal in time domain.

Figure 18: The output signal in time domain.

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

Figure 19: The input signal in frequency domain.

Figure 20: The output signal in frequency domain.

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 and the obtained value .

8. 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  kHz and a magnification coefficient , 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  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-horn-tool 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.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

V. K. Astashev and V. I. Babitsky, “Ultrasonic cutting as a nonlinear (vibro-impact) process,” Ultrasonics, vol. 36, no. 1-5, pp. 89–96, 1998.

R. Singh and J. S. Khamba, “Ultrasonic machining of titanium and its alloys: a review,” Journal of Materials Processing Technology, vol. 173, no. 2, pp. 125–135, 2006.

T. Asami and H. Miura, “Study of ultrasonic machining by longitudinal-torsional vibration for processing brittle materials-observation of machining marks,” Physics Procedia, vol. 70, pp. 118–121, 2015.

V. I. Babitsky, V. K. Astashev, and A. Meadows, “Vibration excitation and energy transfer during ultrasonically assisted drilling,” Journal of Sound and Vibration, vol. 308, no. 3-5, pp. 805–814, 2007.

L. DeFu, W. L. Cong, Z. J. Pei, and Y. JunTang, “A cutting force model for rotary ultrasonic machining of brittle materials,” International Journal of Machine Tools & Manufacture, vol. 52, pp. 77–84, 2012.

J. Wang, P. Feng, J. Zhang, W. Cai, and H. Shen, “Investigations on the critical feed rate guaranteeing the effectiveness of rotary ultrasonic machining,” Ultrasonics, vol. 74, pp. 81–88, 2017.

V. Sharma and P. M. Pandey, “Optimization of machining and vibration parameters for residual stresses minimization in ultrasonic assisted turning of 4340 hardened steel,” Ultrasonics, vol. 70, pp. 172–182, 2016.

M. Roopa Rani, K. Prakasan, and R. Rudramoorthy, “Studies on thermo-elastic heating of horns used in ultrasonic plastic welding,” Ultrasonics, vol. 55, no. 1, pp. 123–132, 2015.

T. V. da Cunha and C. E. N. Bohórquez, “Ultrasound in arc welding: a review,” Ultrasonics, vol. 56, pp. 201–209, 2015.

W. L. Cong, Z. J. Pei, T. W. Deines, A. Srivastava, L. Riley, and C. Treadwell, “Rotary ultrasonic machining of CFRP composites: a study on power consumption,” Ultrasonics, vol. 52, no. 8, pp. 1030–1037, 2012.

P. Feng, J. Wang, J. Zhang, and J. Zheng, “Drilling induced tearing defects in rotary ultrasonic machining of C/SiC composites,” Ceramics International, vol. 43, no. 1, pp. 791–799, 2017.

K. Ding, Y. Fu, H. Su, Y. Chen, X. Yu, and G. Ding, “Experimental studies on drilling tool load and machining quality of C/SiC composites in rotary ultrasonic machining,” Journal of Materials Processing Technology, vol. 214, no. 12, pp. 2900–2907, 2014.

E. Eisner, “Design of sonic amplitude transformers for large magnification,” in Proceedings of the IEEE, vol. 51, pp. 512–512, IEEE.

K. Graff, Wave Motion in Elastic Solids, The Clarendon Press, Oxford, UK, 1975.

D. Sindayihebura, L. Bolle, A. Cornet, and L. Joannes, “Theoretical and experimental study of transducers aimed at low-frequency ultrasonic atomization of liquids,” Journal of the Acoustical Society of America, vol. 103, no. 3, pp. 1442–1448, 1998.

D.-A. Wang, W.-Y. Chuang, K. Hsu, and H.-T. Pham, “Design of a Bézier-profile horn for high displacement amplification,” Ultrasonics, vol. 51, no. 29, pp. 148–156, 2011.

G. Amza and D. Drimer, “The design and construction of solid concentrators for ultrasonic energy,” Ultrasonics, vol. 14, no. 5, pp. 223–226, 1976.

J. J. Coy and F. S. Tse, “Synthesis of solid elastic horns, transactions of ASME,” Journal of Engineering for Industry, vol. 96, no. 2, pp. 627–632, 1974.

X.-P. He and J. Gao, “A review of ultrasonic solid horn design,” Technical Acoustics, vol. 25, no. 1, pp. 82–86, 2006.

V. Salmon, “Generalized plane wave horn theory,” Journal of the Acoustical Society of America, vol. 17, no. 3, pp. 199–211, 1946.

M. P. Kamat, V. B. Venkayya, and N. S. Khot, “Optimization with frequency constraints-limitations,” Journal of Sound and Vibration, vol. 91, no. 1, pp. 147–154, 1983.

J.-C. Hung, Y.-P. Tsai, and C. Hung, “Optimization of ultrasonic plastic welding horns on amplitude uniformity,” Applied Mechanics and Materials, vol. 121, no. 126, pp. 278–282, 2012.

T. He, X.-Q. Ye, and Y. Zhao, “Optimization design for ultrasonic horn with large amplitude based on genetic algorithm,” Journal of Vibroengineering, vol. 17, no. 3, pp. 1157–1168, 2015.

P. Harkness, A. Mathieson, C. Murray, and M. Lucas, “Optimization of ultrasonic horns for momentum transfer and survivability in high-frequency/low frequency planetary drill tools,” in Proceedings of the AIAA SPACE Conference and Exposition 2011, Long Beach, Calif., USA, September 2011.

I. C. Rosca, S. T. Chiriacescu, and N. C. Cretu, “Ultrasonic horns optimization,” Physics Procedia, vol. 3, pp. 1033–1040, 2010.

I.-C. Rosca, M.-I. Pop, and N. Cretu, “Experimental and numerical study on an ultrasonic horn with shape designed with an optimization algorithm,” Applied Acoustics, vol. 95, pp. 60–69, 2015.

K.-M. Shu, Y.-G. Li, C.-C. Chan, and J.-B. Kuan, “Optimized design of the horn of ultrasonic roll welding,” Advanced Materials Research, vol. 482-484, pp. 2223–2226, 2012.

L. Xu, S. Lin, and W. Hu, “Optimization design of high power ultrasonic circular ring radiator in coupled vibration,” Ultrasonics, vol. 51, no. 7, pp. 815–823, 2011.

A. G. Webster, “Acoustical impedance, and the theory of horns and the phonograph,” Proceedings of the National Academy of Science of the United States of America, vol. 5, pp. 275–282, 1919.

E. Eisner, “Complete Solutions of the “Webster” Horn Equation,” Journal of the Acoustical Society of America, vol. 41, no. 4B, pp. 1126–1146, 1967.

H. Popescu and V. Chiroiu, “Calculul structurilor optimale (Calculation of optimal structures – in romanian),” Editura Academiei Române, Bucuresti (Romanian Academy Printhouse), 1981.

M. J. Turner, “Design of minimum mass structures with specified natural frequencies,” AIAA Journal, vol. 5, no. 3, pp. 406–412, 1967.

M. J. Forray, Variational Calculus in Science and Engineering, McGraw Hill Inc., New York, NY, USA, 1968.

V. A. Ilyin and E. G. Poznyak, “Fundamentals of Mathematical Analysis,” Mir Publishers, vol. I & II, 1982.

C. T. Chen, Linear System Theory and Design, Holt, Rinehart and Winston, New York, NY, USA, 1984.

L. M. Silverman, “Realization of linear dynamical systems,” IEEE Transactions on Automatic Control, vol. 16, no. 6, pp. 554–567, 1971.

B. De Schutter, “Minimal state-space realization in linear system theory: an overview,” Journal of Computational and Applied Mathematics, vol. 121, no. 1-2, pp. 331–354, 2000.

D. E. Newland, Mechanical Vibration Analysis and Computation, Longman, Harlow and John Wiley, New York, NY, USA, 1989.

R. Fazlollah, Linear Spaces in Engineering, Ginn & Company, Waltham, Mass, USA, 1971.

J. He and Z.-F. Fu, Modal Analysis, Butterworth-Heinemann, 2001.