Open access peer-reviewed chapter

Modeling and Identification of Parameters the Piezoelectric Transducers in Ultrasonic Systems

By Pawel Fabijanski and Ryszard Lagoda

Submitted: November 15th 2010Reviewed: April 13th 2011Published: September 6th 2011

DOI: 10.5772/21619

Downloaded: 6073

1. Introduction

This chapter is dedicated to ultrasonic piezoelectric ceramic power transducers. These elements are now the most popular source of high power ultrasound and is used in many industrial applications. High power ultrasonic waves are generally used in such industrial processes as welding, acceleration of chemical reactions, scavenging in gas medium, echo sounding and underwater communication (sonar systems), picture transmission, and, above all, ultrasonic cleaning. In practice is now the most widely used the sandwich type power transducers.

Stage design power converters high power ultrasonic devices usually preceded by computer analysis of currents and voltages waveforms the elements of the system, particularly in semiconductor instruments of power. Competent representation requires the use of these waveforms of electrical models of piezoelectric ceramic transducers under the parameters of line with reality and allows to calculate the electrical operating parameters used in the layout of semiconductor switches, capacitors and reactors. Application to simulation circuit of the main generators of ultrasonic piezoelectric ceramic transducers correct model also allows analysis of different variants of control systems and regulation of voltage-frequency converters

For example the standard ultrasonic system for cleaning technology (Fig. 1) includes:

  1. ultrasonic generator,

  2. transducer or set of transducers,

  3. cleaning tank.

Piezoelectric ceramic transducers placed in the tub generate ultrasonic waves that pass through the liquid and reach the element immersed in the tank. As a result, created in the liquid, with very high frequency, alternating areas of high and low pressure.In areas, where low pressure is forming millions of bubbles of vacuum. When the pressure in the alveoli increases and is high enough, bubbles implode, releasing enormous energy at the same time. This phenomenon is called cavitation. Emerging implosions work as a whole series of small cleaning brush. The phenomenon is spreading in all directions and causes intense but controlled detachment of particles of pollutants on the entire surface of cleaning detail. Washed away dirt particles collect on the surface of the cleaning solution from where they are blown into a nearby basin, and then be filtered and recycled.

Figure 1.

Ultrasonic cleaning system

Ultrasonic cleaning is more effective in cleaning hard materials, than the cleaning of soft or porous materials. It was found that, the harder the surface, including the operation of ultrasound is more efficient. Hence, metals, glass, hard plastics well led by ultrasound and are ideally suited for ultrasonic cleaning.

2. Sandwich type transducer

In the technological equipment for cleaning, welding, etc. are generated ultrasound with high intensity and frequency from 20 kHz to 100 kHz. Currently, most teams or a single power are a source of ultrasonic piezoelectric ceramic transducers. The construction of such a transducer is shown in Fig. 2.

Figure 2.

Construction ofsandwich type piezoelectric ceramic transducer; 1 – screw or pin settings (gripping), transmitter, 2, 5 – blocks of metal (eg. aluminum, iron, brass), 3, 4 – piezoelectric ceramic plates (cylindrical, annular)

Such transducers consists of two metal blocks (2, 5), between which are clamped to the material of the piezoelectric ceramic plate (3, 4). Metal blocks and plates are twisted with one or more screws (1). This construction has a much lower own resonant frequency compared to the same frequency of vibration plates, and what is more important allows you to generate high intensity ultrasound. Characteristic of ultrasonic power converters is that they work in a state of mechanical resonance. Thus, in this case the wave frequency of the supply voltage must be equal to the natural frequency of the transducer.

2.1. Mechanical model

Piezoelectric ceramic power converter in resonance state is a mechanically vibrating block, which can model the system with one degree of freedom shown in Fig. 3.

This model consists of massM, which represents the mass of the whole converter, a damper with a coefficient of frictionRand spring with a coefficient of mechanical sensitivityK.

Figure 3.

The model of the mechanical oscillating system with one degree of freedom

In this system, there are four forces: the external forcefZ, the force of elasticity (Hook)fK, the force of friction fRand force of inertiafM, which satisfy the equation:

fZ=fM+fR+fKE1

Assuming: K=1S, where Sthe coefficient of elasticity,

x- deviation from the equilibrium position, v- linear velocity of particles (acoustic) dependence (Eq. 1) can be written as:
fZ=Md2xdt2+Rdxdt+1Kx=Mdvdt+Rv+1K-tvdϢE2

Assuming that the vibration exciting force FZis sinusoidal variable

fZ=FZMsinϧtE3

where: FZMamplitude, and ϧpulse of this force, the relationship (Eq. 2) can be written as:

FZ¯=FZejϧt=jϧMV¯+RV¯+1jϧKV¯E4

Complex mechanical impedance of the transducer is thus equal to:

Zmech¯=FZ¯V¯=R+jϧM+1jϧK=R+jϧM-1ϧK=ZmecheE5

The impedance module is equal to:

Zmech=R2+ϧM-1ϧK2=R2+X2E6
Ϥ=arctgXRE7

Finding a model similar to the mechanical model of electrical converter provides digital modeling of complete systems of generators supplying power electronic converters, the analysis found their current and voltage waveforms and to verify the different concepts of control algorithms of such systems.

The relationship (Eq. 5) describing the complex mechanical impedance of the transducer is similar to the relationship describing the Zm¯impedance of the electrical serial circuitRmLmCmshown on Fig. 4.

Zm¯=Rm+Lm+1Cm=Rm+jϧLm-1ϧCmE8

Figure 4.

Electrical serial circuit RmLmCm

In this circuit the resistance Rmis equal toRm=Rs+Ra, where Rsrepresents the mechanical losses of the converter, which are practically constant, and Rathe acoustic resistance, which is inversely proportional to the intensity Iaproduced by the ultrasound transducer.

2.2. Electric model

Equation (Eq. 5) shows that the linear velocity (molecular), vibrating transducer has the greatest amplitude when the mechanical impedance value Zmechis the smallest. This condition occurs when the elastic force fKbalance the power of inertiafM,

fM+fK=0E9

and pulse of external forcefZ, to compensate for the frictional force fRwill be equal to the pulsation resonant converterϧm, whereϧm2=1MK.

In the piezoelectric phenomenon, the relationship between external force fZcausing the deformation of the piezoelectric ceramic element and the intensity of the electric field E generated inside, is linear.

As for the specific thickness piezoelectric ceramic plate l we can write:

u=ElE10

thenassuming a coefficient of proportionality kpcan be said that an external force kpisproportional to the voltage that was applied to the plate.

fZ=kpuE11

Comparing the equation describing a serial electrical circuitRmLmCm,

u=Lmd2qdt2+Rmdqdt+qCm=Lmdidt+Rmi+1Cm-tidϢE12

of equation (Eq. 2) can be seen that supplied to the circuit an electric charge is proportional to the deformation:

q=kpxE13

Therefore meets the current relationship:

i=dqdt=kpdxdt=kpvv=1kpiE14

After inserting equation (Eq. 1) and (Eq. 14) to (Eq. 2) is:

kpu=Mkpdidt+Rmkpi+1kpK-tiϢE15
u=Mkp2didt+Rkp2i+1kp2K-tiϢE16

Since the components of the sum in equation (Eq. 16) May dimension of tension is present in them can be replaced by permanent mechanical respective electrical components:

Lm=Mkp2,Rm=Rkp2,Cm=kp2KE17

2.2.1. Linearmodel

On the basis of equation (Eq. 16) can be drawn, well-known and frequently used, an electric equivalent circuit oscillating piezoelectric ceramic transducer shown in Fig. 5.

Additional components Ceand Rerepresent the electrical capacitance and resistance of the piezoelectric ceramic plate transducer.

Figure 5.

Linear, a replacement wiring diagram piezoelectric ceramic transducer,ZPT - Electrical impedance of the transmitter: a) a detailed diagram, b) the distribution of the branch of electrical and mechanical

Figure 6.

Frequency characteristics of impedance modulus normalized and the phase of the mechanical industry Zm and the entire transducerZPT

Examples of the frequency characteristics of a standard impedance modulus and phase of mechanical industry and the entire transducer is shown in Fig. 6.

When the frequency f1reaches a minimum impedance ZPTof the module. Near this frequency is the frequency of mechanical resonancefr.

Resistance Rmcan be divided into two parts:

Rm=Rs+RaE18

whereRsrepresents the mechanical losses of the converter, which are practically constant, and Rathe acoustic resistance.

Presented in this chapter, a linear model of the piezoelectric ceramic transducer is not mapped correctly the actual frequency characteristics of the transmitter. Frequency response of linear model is much wider than in the actual transmitter. This is important in the process of creating digital models and design of ultrasonic generators.

2.2.2. Non-linear model

As mentioned above, the linear model presented in Fig. 5 does not correctly reproduces the frequency characteristics of the transmitter. This feature can be eliminated by supplementing the model with two non-linear electrical resistanceRm1fand Rm2fthe accompanying in mechanical industries, as shown in Fig. 7.

The resultant mechanical resistance Rmwfin the industry is thus equal to:

Rmwf=Rm+Rm1fRm2fRm+Rm1f+Rm2fE19

Assuming that

Rm1f=Rm1fr-ffrnE20
Rm2f=Rm21+fr-ffrmE21

the current state of the branch of mechanical resonance, which is proportional to the speed of molecular vstill is determined by the resistanceRmbecause, Rm1ffr=0andRm2f(fr).

Figure 7.

Non-linear electrical model of piezoelectric ceramic transducer

3. Admittance characteristics of sandwich type transducers

The electrical characteristics of ceramic transducer resistance Re and electrical capacitance Cecan be determined with high accuracy on the basis of electrical measurement. Changes in these parameters while the transducer can be considered to be negligible because of the transmitter power is assumed that the CeandReare fixed. Parameters of dynamic branch RmCmLmof power piezoelectric ceramic transducer model electric are not physically measurable. We can only designate empirically. Identification algorithm is therefore based on numerical calculations. Input parameter for these calculations is the image of the actual characteristics of the transducer admittance.

Example characteristics of the actual and approximated admittance converter a resonance frequency around 25 kHz are shown in Fig. 8.

Assuming a digital model of the transmitter as shown in Fig. 5 the susceptanceB, conductanceG, admittance Yof the transducers in frequency ϧfunction may be expressed by the following equations:

Bϧ=ImYϧ=ϧCe+ϧCm(1-ϧ2LmCm)ϧ2Cm2Rm2+(1-ϧ2LmCm)2E22
Gϧ=ReYϧ=Ge+ϧ2Cm2Rmϧ2Cm2Rm2+(1-ϧ2LmCm)2E23
Yϧ=G(ϧ)2+B(ϧ)2E24

whereGe=1Re. For the resonant pulsation ϧm=1LmCmare:

Bϧm=CeLmCm0E25
Gϧm=1Re+1Rm1RmE26

Figure 8.

Exemplary admittance characteristics of sandwich type transducers1. Real characteristic,2. Equivalent characteristic

4. The genetic algorithm in use to identifyRm, Lm, Cmparameters of the mechanical branch

As mentioned earlier parameters of dynamic branch RmLmCmof power piezoelectric ceramic transducer model electric are not physically measurable.

We can only designate empirically. Identification algorithm is therefore based on numerical calculations. Input parameter for these calculations is the image of the actual characteristics of the transducer admittance.

An example measurements of results for the transducer mechanical resonance frequency equal to about 43 kHz is locate in appendix in Table 1 and one-to-one correspondence real admittance characteristics of this transduceris shown in Fig. 9. The graph consists of points whose coordinates correspond to the conductance G=Re(Y)and susceptanceB=Im(Y)of the converter, measured at a certain frequency.

Analysing the shape of this characteristic can be clearly observed that with increasing frequency (Figure accordance with the frequency increasing clockwise), outlines the main loop of the graph. In the general case, the image of the curve is an ellipse. This ellipse is interpolated electrical admittance characteristics of an ideal replacement transducer schedule shown in Fig. 9. It may be noted that the actual characteristics, in addition to the main loop also includes many smaller "loops" that testify to the presence of additional resonances in the transducer side. These resonances, however, will not occur in the adopted system replacement transmitter.

Figure 9.

Real admittance characteristics of the piezoelectric ceramic transducer

After eliminating these "loops" with the actual picture of the actual characteristics of the admittance characteristics of the transducer will be the figure presented in Fig. 10.

In Table 1 are marked in bold points, which missed the chart.

Figure 10.

The actual characteristics of the transducer after eliminating resonances fringe

The set of coordinates of the actual approximation characteristics transducer (Fig. 9) is a database for further calculations.

4.1. Calculationalgorithm

Designation of alternative modes of dynamic parameters of the transmitter requires the implementation of numerical calculations and find such a set of RmLmCmparameters, which in a given error will allow mapping of the main loop, the actual characteristics of the transducer admittance presented Fig. 10.

One element of the identification of these parameters is a genetic algorithm.Generally one can say that thegenetic algorithm is a regula - system learner who makes a certain number of iteration steps. It is a group of strict security procedures that are based on the fundamental mechanisms of biological evolution such as natural selection and inheritance. It works interactively with the environment in discrete time. The algorithm of this type of reproduction may take place subject to the diversity of the population:

  • model with preload,

  • measures niche

A special feature of this type is that the algorithm is not seeking a single optimal solution, but a group of cooperating the best solutions. At any time, the algorithm works evolutionary principle of survival, which is always available some of the best solutions at the moment. As the proceedings algorithm solution to optimize and adapt to the conditions in which the algorithm works. For further calculations is always the best solutions are selected and rejected solutions are worse. Here there is a process of succession. In order to obtain optimal solutions group for further reproduction of the best solutions are selected at the time of the algorithm. It is a natural selection process occurs in nature. It is known that the probability of obtaining better result is greater if we use it to generate the best available solution than the use of inferior results. According to the law of nature and genetics "survive" the best and strongest.

Figure 11.

Simply genetic algorithm

The genetic algorithms is a natural process of mating genes. It is no different genetic algorithm. There also are choosing the best "genes" of each solution and verify their combinations. Alongside the cross as a natural evolution is mutation, a random change in the gene. Both these processes are the values of the genetic operator.

Using the genetic algorithm, remember to keep the best balance between the transfer of genes to the next generation, and a draw solution space. Too broad conditions imposed solutions may give erroneous results in spite of every generation the best available solution at the moment. Genetic algorithm is an excellent tool to monitor and maintain the balance between these two dependencies.

The overall pattern of genetic algorithm is illustrated in Fig. 11.

Using a genetic algorithm to identify the parameters of dynamic model of the electrical branch of the transducer sandwich working near mechanical resonance is shown in Fig. 12.

Figure 12.

Algorithm for parameter identification of dynamic model of the branch of electrical power ultrasonic transducer

The initial value of the dynamic resistance Rmis calculated by selecting the characteristics shown in Fig. 9 point with coordinates(Gmax,Be). In this point there is a mechanical resonance of the transducer and is a good approximation condition.

ϧm2LmCm-1=0E27

The dependence (Eq. 26) shows that the resonance:

Gϧm=Gmax=Ge+1RmE28

Therefore:

Rm=1Gmax-GeE29

Figure 13.

Iterative improvement of the results obtained in numerical calculationsand a computer simulation electric model of an sandwich type transducers

Choosing the actual characteristics of the admittance of any two points with coordinates(G1,B1), (G2,B2), which lie beyond the point of mechanical resonance(ϧ1ϧm,ϧ2ϧm)can calculate the initial value of the replacement of mechanical capacityCm.

Transforming the system of equations:

G1-Ge=ϧ12Cm2Rm(1-ϧ12LmCm)2+ϧ12Rm2Cm2E30
G2-Ge=ϧ22Cm2Rm(1-ϧ22LmCm)2+ϧ22Rm2Cm2E31

and taking into account that for Gmaxis equal to the pulsation vibrations ϧmapproximate value of Cmis described by the formula:

Cm'=1-ϧ12ϧm2ϧ1G1-Ge(Gmax-Ge)1-G1-GeGmax-GeE32
Cm''=1-ϧ22ϧm2ϧ2G2-Ge(Gmax-Ge)1-G2-GeGmax-GeE33
Cm=Cm'+Cm''2E34

The initial value of the replacement of mechanical inductance Lmdetermined from the relationship:

Lm=1ϧm2CmE35

After substituting the calculated value of the RmLmCmand the measurementCe,Rethe relationship (Eq. 22) and (Eq. 23) is determined by numerical coordinates of the points of the electrical characteristics of the admittance model for the same pulse, for which measurements of actual performance. If the error resulting from a comparison of the approximation characteristic with real characteristic (Fig. 10) is greater than the accepted values, generate a new population ofCe,Re, RmLmCmparameters, repeat the calculation of coordinates of the electrical characteristics of the model and then analyze the resulting error. Calculations should be continued until the error resulting from a comparison of the characteristics shown in Fig. 10 with the characteristics of the accepted model of the electrical transducer is smaller than the set value. Measurement error should be performed in the pulsation ϧchanging between points Bmaxand Bmin(Fig. 10).

Procedure described above can be used and introduces the appropriate algorithm in the DSP program simulation.

//Program piezoelectric ceramic transducer

input {Circuit parameters [Re,Ce,R,L,C,u(t)];

Simulation parameters (tp,tk,Δ,u(0))

Evolution parameters (Num,Size, NumM,Initial,NumG,stop);)

Output (t, u;tp; uo=u(o);};

for (t,tk, t++){New function form;

start initial population;}

for (integer i=1,Num,i++){New population,

Select the best}

where:

Genotype dimension: Num,

Size of the generation created by mutation: Size,

Number of the best selected genotypes for mutation: NumM,

Number of the best selected genotypes for crossing: NumC,

initial parameter mutation range: Initial

4.1.2. Description of algorithm parameters

The algorithmbuildsan initialpopulation from completelyrandomvalues. Eachalgorithm result isthereforepartlyrandom. Valuesimplicitlyincludedinthe programwhenitstartsin mostcasesproduce the desiredresult. In orderto increasethe effectiveness ofthe algorithm we canincrease thepopulation sizeorquantity. This implies, however, withprolongedtime neededforcalculations, andthus thewaiting timeforresults. Increasingthe likelihoodmay leadtoan erroneousfindingandlosing theoptimal solution. In thecase ofreduction inpopulation size, increasedthe likelihood ofcrossinghas a better chanceofgettingan individualwitha desiredadjustmenteven thoughhisspeechislesscertain.

Similaralgorithmbehavesina situationof increasingthe probabilityof mutation, but inthiscase, the resultis morerandom.

4.1.3. Description of program code

The main application file is interface.m. He is responsible for building the interface, retrieve data from the controls and calling the functions of the buttons.

The first step is to create an initial population for genetic algorithm.

// INITIAL POPULATION

startpop = zeros (popsize, 2);

fori = 1: popsize

startpop (i,:) = individual (2, max_val); / of the initial random population

end

The creation of the initial population is to create a matrix of population size 'number of individuals' x 'dimension of the individual', where the number of individuals is the value declared in the program window, and the dimension is equal to the number of individual factors.

The population is filled with random individuals.

The next step is to assess the population. It is used here for the sum of squared differences of the values obtained and references. Since the genetic algorithm seeks to maximize, a function of adaptation used in this case is the inverse of the sum of the squares, in addition to the square have raised to accelerate the convergence of the algorithm.

Based on the resulting evaluation, subjects were assigned to be the appropriate probability of the next population. The draw of the population makes the next feature:

function [new] = sel (evaluation, current) / disproportionate Roulette

globalpopsize;

for j = 1: popsize

number (j, 1) = round (100 * (score (j, 3)));

end

roulette = zeros (100.2);

fori = 1:100

Roulette (i, 1) = i;

end

// 100-point roulette wheel

all = 100;

item = popsize;

while an element> 0; / complement vector numbers in accordance with the adaptation of individual

for a = 1: all

Roulette (a, 2) = element;

end

all = all-number (element, 1);

element = element-1;

end

// Draw elements

items = Random (popsize, 100); / draws n = popsize position

for j = 1: popsize

number (j, 1) = roulette ((positions (j, 1)), 2); / individual in this position

end

fori = 1: popsize

new (i,:) = current ((number (i, 1 )),:); / choose a subject to reproduce

end

After drawing the new population are made on the operations of crossover and mutation.

// CROSSING

xn = fix (popsize / 2); / number of crossover operations

remainder = mod (popsize, 2); / population parity check

steam = select (popsize); / choice pairs of interbreeding

fori = 1: xn

[Child1 Child2] = cross ((new (pair (2 * i-1, 1 ),:)),( new (pair (2 * i, a cross ),:)));

pox (i,:) = Child1;% first child

pox (i+ xn,:) = Child2;% second child

end

pox (2 * xn +1,:) = new (pair (xn, 1 2 * ),:); / last rewritten without change

end

Crossing begins by selecting pairs of individuals, in such a way that a given individual could occur in only one hand. If the population is odd that an individual will be rewritten without change. Then it is drawn for each pair whether the operation will cross (the probability of randomly selected cross <put the probability of crossover):

- if so, part of the code is converted into fish,

- if not, individuals are no changes.

According to the function below:

function [Child1, Child2] = cross (A, B); / crucifixion simple

globalpx;

or = rand (1);

if and>px / cross check whether

Child1 = A;

Child2 = B;

else / crossing straight

Child1 (1) = A (1);

Child1 (2) = B (2);

Child2 (1) = B (1);

Child2 (2) = A (2);

end

After the crossover operation for each individual of the population by analogy, we check if there is a mutation. If so, an individual code rate is changed to a random value according to the following function:

function [X] = mut (A) / mutation (replacing a randomly selected for another random value)

globalmax_val;

index = Random (1,2); / random mutation alleles

new = individual (1, max_val); / new value allele

A (index) = new;

X = A;

After mutation we get a new population, which also is subjected to evaluation of the adjustment operations, reproduction, crossover and mutation according to the number of iterations the algorithm set by the user.

Examples of results obtained in subsequent iterative steps for the transducer mechanical resonance frequency of 42.9 kHz is shown in Fig. 10.

When the next iterative step will be the condition for the completion of the calculations should be considered that the characteristics of the model shown in Fig. 5 parameters Re, Ceand RmLmCmgenerated in the genetic algorithm in the last population of values, coincides with the actual characteristics of the interpolated transducer. The shape of this characteristic is shown in Fig. 14.

Figure 14.

Circle admittance obtained using genetic algorithm

4.2. Effect of the electrical equivalent circuit parameters on the shape characteristics of the transducer admittance

Location of figures which creates a characteristic admittance of the transducer coordinate system Bϧ=f[Gϧ]can be determined by changing the value of capacity Ceand the loss of electrical resistance of ceramicsRe. The higher value of Rethe greater displacement figures to the right people along the conductanceG(ϧ). If you forget in the pattern of replacement sensor resistance Re(Re)will move to the left approximation of thecharacteristics. In this case, the graph will be almost tangential to the axisB(ϧ). Reduction in the capacity of Cewill move down the sheet, along the axisB(ϧ). Conductance Gmaxvalue determines the width of the main loop of the graph. From relation (28) shows that the value of Gmaxdecreases when the resistance Rmincreases.

The impact of dynamic capacity Cmcan analyze the shape characteristics of saving the relationship (Eq. 22) in the form

B(ϧ)ϧCm(1-ϧ2ϧm2)ϧ2Cm2Rm2+(1-ϧ2ϧm2)E36

For small, ± 2% of the pulsation untuning mechanical resonance, it can be assumed that:

B(ϧ)1-ϧ2ϧm2ϧCmRm2E37

From relation (Eq. 30) shows that reducing the value of the dynamic capacity Cmwill increase susceptance system. This also applies to the maximum Bmaxand minimumBmin. In this case, it means that the values of Bmaxand │Bmin│ increased. The difference BminBminis equal to the amount of main loop of the graph.

These rule changes in the electrical equivalent circuit parameters of the transducer should be used during the generation of successive populations of a set of values ofRe, Ceand RmLmCmgenetic algorithm.

5. Conclusion

Linear model presented in Fig. 5 does not correctly reproduces the frequency characteristics of the transmitter. This feature can be eliminated by complementing the electrical model of two nonlinear resistances Rm1fand included Rm2fin the mechanical industry, as shown in Fig. 7 This is important in the process of creating digital models and design of ultrasonic generators.

The actual characteristics of the admittance is a source of transparent, condensed information about the properties and parameters of the piezoelectric ceramic transducer stimulated by the mechanical vibrations. Used to identify the dynamic parameters of RmLmCmand optimization of electrical parametersRe, Ceoccurring electric transducer model, genetic algorithm uses these data and represents one of many possible methods to use here. It has the advantage that it can be realized in automatic cycle.

Based on preliminary data obtained experimentally from measurements of electrical parameters(Re, Ce)and measurements of voltage and current waveforms created image of the actual transducer admittance characteristics ofBϧ=f[Gϧ]. Genetic algorithm can find the optimal approximation of the characteristics, correction of the value ofRe, andCeand calculation of parameters RmLmCmmodel of the electrical industry dynamic.

An important part of the described method is the proper preparation of input data. They should be given in the form of an ordered table of coordinates of points lying on the actual characteristics of the transducer according to those points that do not belong to the main loop of the graph. The process of eliminating points of the loop parasitic resonances can be implemented step by step analysis of the data included in the table of measurement results or by analyzing the image created by a set of characteristic points ofBϧ=f[Gϧ].

Attempts to identify the parameters of a single transducer mechanical resonance frequency of 23.8 kHz and 43.2 kHz have confirmed the effectiveness of the genetic algorithm. Identification algorithm described above can be particularly useful in studies of larger groups of ultrasound transducers in high-power washing.

By introducing the necessary modifications can also be used to identify the parameters of transducers working in real conditions and to implement control systems that track the frequency of mechanical resonance.

6. Appendix

Table 1 shows selected results of detailed studies of the piezoelectric ceramic transducer mechanical resonance frequency of 42.9 kHz. Measurements were performed impedance meter HP 4192 type IMPEDANCE Analyzer. Indicated in bold letter omitted in the genetic algorithm performance admittance points beyond the main loop of the graph. The contents of Table 1 refers to Fig. 9, and Fig. 10.

f
kHz
R
X
tg(φ)
φ
rad
cos φ
sin φ
G
mS
B
mS
42,000,0721-0,23383,24261,271660,29470,95561,2043204,650189
42,040,0674-0,21323,16591,264850,30120,95361,3469224,223453
42,080,0703-0,19252,73711,220520,34320,93931,6747034,456680
42,120,0841-0,17452,07501,121710,43410,90082,2410504,963926
42,160,1006-0,18091,79881,063420,48590,87402,3479064,296001
42,200,0921-0,17631,91381,089300,46310,88632,3286933,289479
42,240,0888-0,14831,67051,031380,51360,85802,9715943,192256
42,280,1148-0,13581,18340,869190,64540,76383,6303033,907909
42,320,1515-0,14030,92550,746750,73390,67933,5540973,487659
42,360,1501-0,20151,34230,930510,59740,80192,3781893,725892
42,400,1047-0,20141,92311,091270,46140,88722,0321283,739428
42,440,1432-0,15021,04920,809420,68990,72393,3239914,162891
42,480,1212-0,19181,58311,007420,53400,84552,3535074,431376
42,520,1227-0,18681,52280,989730,54890,83592,4556374,700036
42,560,1041-0,18001,72841,046290,50080,86562,4084784,873619
42,600,0992-0,16661,67911,033650,51170,85922,6391184,997220
42,640,0980-0,14781,50790,985220,55270,83343,1169304,959969
42,680,0973-0,13511,38840,946600,58440,81143,5102814,685777
42,720,0974-0,12291,26100,900340,62130,78353,9628064,035905
42,760,1005-0,10811,07520,821620,68100,73224,6130982,849480
42,800,1057-0,09240,87380,718140,75300,65805,3625971,944590
42,840,1139-0,07520,65990,583310,83460,55086,1158132,348010
42,880,1295-0,05710,44080,415170,91500,40336,4645112,602020
42,920,1542-0,05140,33310,321520,94880,31605,8383411,397040
42,960,1582-0,07040,44510,418800,91360,40675,274907-0,305750
43,000,1374-0,05780,42070,398190,92180,38786,185462-0,238180
43,040,1365-0,02710,19820,195690,98090,19447,047725-0,525000
43,080,15090,0070-0,0463-0,046220,9989-0,04626,610636-1,104520
43,120,18800,0084-0,0449-0,044830,9990-0,04485,309238-1,992950
43,150,15850,0133-0,0838-0,083620,9965-0,08356,263456-1,955190
43,170,15770,0283-0,1798-0,177890,9842-0,17706,143423-1,362000
43,200,16970,0661-0,3894-0,371350,9318-0,36295,117705-1,591100
43,240,20420,1018-0,4985-0,462450,8950-0,44613,922151-1,756280
43,280,25790,1058-0,4104-0,389400,9251-0,37963,319092-1,271460
43,320,23950,1108-0,4627-0,433390,9075-0,42003,438489-0,718400
43,360,26380,1777-0,6736-0,592800,8294-0,55872,6072434,650189
43,400,35300,2199-0,6231-0,557200,8487-0,52882,0406784,223453

Table 1.

Results of laboratory tests by the ultrasonic transducer sandwich-type of the resonance frequency around 43 kHz

© 2011 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike-3.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited and derivative works building on this content are distributed under the same license.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Pawel Fabijanski and Ryszard Lagoda (September 6th 2011). Modeling and Identification of Parameters the Piezoelectric Transducers in Ultrasonic Systems, Advances in Ceramics - Electric and Magnetic Ceramics, Bioceramics, Ceramics and Environment, Costas Sikalidis, IntechOpen, DOI: 10.5772/21619. Available from:

chapter statistics

6073total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Influence of Dopants, Temperature and Atmosphere of Sintered on the Microstructure and Behavior of Lead Free Ceramics

By Maria A. Zaghete, Francisco Moura, Alexandre Z. Simões, José A.Varela and Elson Longo

Related Book

First chapter

Advanced Ceramic Target Materials Produced by Self-Propagating High-Temperature Synthesis for Deposition of Functional Nanostructured Coatings - Part 1: Four Elements and Less Systems

By Evgeny A. Levashov, Yury S. Pogozhev and Victoria V. Kurbatkina

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us