Analyses of the properties of conical piezoelectric acoustic emission transducer are presented here. The conical piezoelectric transducer is used to scanning normal displacement on a surface of testing constructions in NDT methods, such as acoustic emission. The simplified one-dimensional analysis based on the equivalent circuit of conical waveguide is recapitulated. FEM analysis on this transducer is presented. Both the frequency response characteristic of the transducer and his time response on pulse excitation are results of this analysis. All of computations are done for various cone angles of piezoceramic (PTZ-H5) element and various sizes of the cylindrical backing of the sensor

Hora P.

Červená O.

Directory of Open Access Journals

Applied and Computational Mechanics 2 (2008) 13–24Analysis of the conical piezoelectric acoustic emission transducerO. ˇ Cerven´ aa,∗,P .H o r a aaInstitute of Thermomechanics of the ASCR, v.v.i., Veleslav´ ınova 11, 301 14 Plzeˇ n, Czech RepublicReceived 28 August 2008; received in revised form 2 October 2008AbstractAnalyses of the properties of conical piezoelectric acoustic emission transducer are presented here. The conicalpiezoelectric transducer is used to scanning normal displacement on a surface of testing constructions in NDTmethods, such as acoustic emission. The simpliﬁed one-dimensional analysis based on the equivalent circuitof conical waveguide is recapitulated. FEM analysis on this transducer is presented. Both the frequency responsecharacteristic of the transducer and his time response on pulse excitation are results of this analysis. All of com-putations are done for various cone angles of piezoceramic (PTZ-H5) element and various sizes of the cylindricalbacking of the sensor.c   2008 University of West Bohemia in Pilsen. All rights reserved.Keywords: acoustic emission, conical transducer, FEM1. IntroductionA conical transducer is a highly sensitive wide-band device (operating up to several MHz)for measuring a vertical component of displacement of a small area on the surface of a body.This transducer is designed for a wide region of applications, e.g., for testing by means ofacoustic emission, to be used as a standard transducer, and the like. The scheme of a coni-cal transducer is shown in ﬁg. 1. The basic characteristics of the transducer are as follows:∅ 38 mm25mm∅ 1m mFig. 1. A conical transducer1) The active element is of piezoelectric ceramics and is conicaland its polarization is parallel with the axis of the cone;2) The greater base of the cone with scalded silver or gold elec-trode is soldered on a relatively great cylindrical backing blockmade usually from brass;3) The smaller base of the cone with scalded silver or gold elec-trode is connected through the medium of as thin binding layeras possible with the part of relatively large surface on which ismeasured the vertical component of displacement;4) The output voltage is measured between the brass cylindricalbacking block and the surface of the active element.5) The output voltage signal of a conical transducer is directlyproportional to the normal surface displacement at the contact area with frequency up to fewMHz.The conical transducers were developed by Proctor [7] in the 1980s for quantitativeacousticemission at the National Bureau of Standards (NBS), now the National Institute of Standardsand Technology (NIST), USA.∗Corresponding author. Tel.: +420 377 236 415, e-mail: cervena@cdm.it.cas.cz.13O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–24Fig. 2. Various kinds of conical transducersThere are many various transducers based on the NBS conical transducer. Some of them areshown in ﬁg. 2. The top left image shows the unique embedded acoustic emission sensor usedby Glaser [2]. The top right image shows the robust impact-echo sensor used by Sebastian [9].The transducer was designed for continuous use up to at least 600 ◦C. The bottom left imageshows the conical transducer used by Theobald [11] as an energy source for system calibra-tion. The last image shows the new type miniature–conical transducer for acoustic emissionmeasurements used by Lee [4].1D analysis based on the model of the conical transducer developed by Greenspan [3] ispresented in the second section. Greenspan extended the well-known Mason model [5] to ac-count for the conical piezoelectric element. The third section deals with FEM analysis of theconical transducer.14O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–242. One-dimensional analysisNow we are going to determine an equivalent circuit of the conical transducer, and to calculateits frequency characteristic from geometrical dimensions and from material constants of anactive element and terminal impedances.By designing an equivalent circuit we will use as a basis an impedance matrix of a conicaldivergent waveguide, see ﬁg. 3 — left in solid phase [5]Zm =⎡⎢⎢⎣ρcS1 cot(kl)i+1ikx1 −ρc S1S21isin(kl)ρc S1S21isin(kl)−ρcS2 cot(kl)i−1ikx2 ⎤⎥⎥⎦,where ρ is a density, c is a phase velocity, i is imaginary unit, S1,S 2 are areas of bases of a conewaveguide, k is wave number, l is the length and x1,x 2 are parameters of a cone waveguide,see ﬁg. 3 — left. This impedance matrix was derived on the assumption that the sections per-pendicular to the axis of waveguide remain by deformation planar, the axial stress is uniformlydistributed all over the surface and radial displacements can be neglected. Comparing the aboverelation with the matrix of T-element shown in ﬁg. 3 — right,0 x1 x2 xr1r2lαF1 F2ZA ZBZCv1 v2Fig. 3. A conical divergent waveguide (left), an equivalent circuit (right)ZT = ZA + ZC −ZCZC −ZB − ZC (1)we get for ZA, ZB and ZC the following relationsZA = ρcS1 cot(kl)i+1ikx1−qisin(kl) (2)ZB = ρcS2 cot(kl)i−1ikx2−1/qisin(kl) (3)ZC = ρc S1S21isin(kl)(4)where q = x2/x1 = r2/r1.15O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–24F1 F2ZA ZBZCv1 v2U-C0C01:kbFig. 4. Equivalent circuit of a conical piezoelectric elementIf the T-element is supplemented by Mason’s model, we get an equivalent circuit of a piezo-electric conical element, see ﬁg. 4. Here kb is a transformation ratio and C0 is a transducercapacity.Transferring the electric part to the mechanical side and on supplementing the circuit bysource impedance Z1 and load impedance Z2 (see Norton’s theorem — velocity short-circui-ted+parallel impedance Z1), we get an equivalent circuit of a conic transducer, indicated inﬁg. 5. It should be remarked that v, v1 and v2 are loop velocities (not branch velocities), forinstance, the velocity on Zc is v1 − v2.For the circuit shown in ﬁg.5 it holdsv1 − v2 = vF, (5)whereF =Z1 (ZB + Z2)(ZA + Z1 + ZC)(ZB + Z2 + ZC) − Z2C(6)Z1 Z2ZA ZBZCv1 v2vUkb-C0/k2bC0/k2bFig. 5. Equivalent circuit of a conical transducer16O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–24andUkb =v1 − v2iωC0/k2b. (7)On substituting for v1 − v2 from equation (5) into equation (7), we obtainUkb =viωk2bC0F. (8)Further, substituting ξ = v/iω and −p33 = kb/C0 we getU/ξ = −p33F, (9)which is the resultant relation for sensitivity of a conical transducer.The value of p33 lies between h33 (disc) and g33/sD33 (rod). h33 is piezoelectric constant (forcalculation we used the value 21.5 · 108 V/m). g33 is piezoelectric constant (for calculation weused the value 24.9 · 10−3Vm/N) and sD33 is the yielding constant at D =0(for calculation weused the value 9.46 · 10−12 m2/N).Substituting relations (2)–(4) into relation (6) for individual impedances, we get the numer-ator of fraction FNF = z1 z2 + i 1qsin(kl)+q − 1qkl− cot(kl)  (10)and denominator of fraction FDF =1+z1z2 +(q − 1)2q(kl)2 −(q − 1)2qklcotkl++ i z1 q − 1qkl− cot(kl) − z2 q − 1kl+c o t ( kl)  .(11)Here z1 and z2 are scaled speciﬁc impedances Z1/ρcS1 and Z2/ρcS2 (the scaling is carried outto ρc of the cone!).In the case of the cylindrical transducer, q =1 ,w eg e tf o rFF =z1 [z2 + itan(kl/2)]1+z1z2 − i(z1 + z2)cot(kl). (12)For kl an odd multiple of π, the numerator and denominator in equation (12) are both inﬁ-nite, but the ratio is ﬁnite; however, for kl a multiple of 2π, i.e., if length l is equal to integralmultiple of wave length, the denominator is inﬁnite for any z1 and z2, hence, both F and sensi-tivity U/ξ are zero. For the case of a cone, the expression 1/(qsin(kl)) for q>1 preserves theﬁnite value of fraction F in equation (6). Instead of zeros, in this case, there are only minima inthe characteristic. This is the reason why a conical transducer yields a suitably even wide-bandcharacteristic (which cannot be achieved with a cylindrical transducer), see ﬁg. 6.In practical applications the signals detected by a conic transducer are so short that both thetestmediumand thebacking may besupposedto behalf–spaces. Theproblemof theimpedancein a half-space (values of Z1 and Z2), considered in a circular area with radius a on the surfaceof a half-space, was solved in [3].17O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–240 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 288.599.510f [MHz]log10|U/ξ|[V/m]cylinderconeFig. 6. Frequency dependencies of sensitivity for cylindrical (r1 = r2 =0 .5,l=4 mm) and conicaltransducer (r1 =0 .5, r2 =3 , l =4mm)3. FEM analysisFEM calculations were performed in the commercial environment COMSOL Multiphysics [1].Because of the axial symmetry of the geometry, an axisymmetric 2D application mode wasused. Elements were the Lagrange-Quadratic type. It was used the linear system solver, direct(UMFPACK).Computations were performed for nine various piezoelectric element sizes, see tab. 1–left.The parameters in this table correspond to labelling in ﬁg. 3 — left and q is a ratio of radii ofthe cone. The inﬂuence of the backing size was studied on four various examples, see tab.1 —right.Table 1. Transducer parameters: element (left), backing (right)elementType r1 [mm] r2 [mm] q[–] l[mm] α[deg]cylinder 0.5 0.5 1 4.0 0.00cone 10 0.5 1.0 2 4.0 7.13cone 15 0.5 1.5 3 4.0 14.04cone 20 0.5 2.0 4 4.0 20.56cone 25 0.5 2.5 5 4.0 26.57cone 30 0.5 3.0 6 4.0 32.01cone 35 0.5 3.5 7 4.0 36.87cone 40 0.5 4.0 8 4.0 41.19cone 45 0.5 4.5 9 4.0 45.00backingType height radius[mm] [mm]A 25.0 3.0B 25.0 11.0C 12.5 19.0D 25.0 19.018O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–24Material of piezoelectric elements was lead zirconate titanate (PZT-5H) with parameters:elasticity matrixc =⎡⎢⎢⎢⎢⎢ ⎢⎢⎢⎢ ⎢⎣127.205 80.2122 84.6702 0 0 080.2122 127.205 84.6702 0 0 084.6702 84.6702 117.436 0 0 0000 2 2 .9885 0 00000 2 2 .9885 000000 2 3 .4742⎤⎥⎥⎥⎥⎥ ⎥⎥⎥⎥ ⎥⎦GPa,coupling matrix (stress-charge form)e =⎡⎢⎢⎣00 0 0 1 7 .0345 000 0 1 7 .0345 0 0−6.62281 −6.62281 23.2403 0 0 0⎤⎥⎥⎦C/m2,permittivity matrix (stress-charge form)εS =⎡⎢⎢⎣1704.40 00 1704.400 0 1433.6⎤⎥⎥⎦,and density ρ = 7500kg/m3.Thetransducerbackingwas createdbybrasswithYoung’smodulusE = 110GPa, Poisson’sratio ν =0 .35 and density ρ = 8700kg/m3.Twokindsofproblemswererealized. Firstly,thefrequency responseanalysis, andsecondly,the time dependent analysis. A mapped mesh consisting of quadrilateral elements was createdin both cases. But size of elements was different.12345678Fig. 7. BoundaryThe following boundary conditions were set for both of cases.The numbering of boundaries is presented in ﬁg. 7. The constraintwas given by prescribed displacement on boundary2 and on the otherboundaries it was free. The electric boundary conditions were ax-ial symmetry on boundary1, ground on boundary2 and zero chargeon boundary6.The boundary integration variable Vint was deﬁned on boundary4as Vint = V 2πr,h e r eV is a voltage and r is a radius.3.1. Frequency response analysisThe maximal size of one mesh element was 0.5 × 0.5mm. The amplitude of the prescribeddisplacement(in axis ofsymmetry-direction)onboundary2 (see ﬁg. 7) was 1pm. TheRayleighdamping was used with mass damping parameter αdM = 209440s−1 and stiffness dampingparameter βdK =1 .06 · 10−8s. It was solved frequency response by parametric solver. Theexcitation frequency was in range from 103 to 106 with step 103 Hz. We studied the responseof voltage U (U = Vint/(πr22)) on boundary4,s e eﬁ g .7 .19O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–24The analysis was performed both for whole transducers and for separate piezoelectric ele-ments. It was always done for nine shapes of cone (see tab. 1 — left). Because some verysimilar jobs were performed, we used advantageously MATLAB [10] for scripting these tasksof COMSOL. We created three scripts. The ﬁrst one was for whole transducers, the second onefor separate elements and the third one was for transducer with the “cone 30” element and fourvarious sizes of backing (tab. 1). We deﬁned in those scripts mapped mesh, materials, boundaryconditions, boundary integration variable Vint, damping and excitation frequency. The problemwas solved for particular shapes of geometry, which we changed in loop.The ﬁg. 8 — top displays the frequency responses of transducers with backing D for par-ticular piezoelectric element types (tab. 1). The curves get smoother the bigger a ratio of theelement radii is. The ﬁg. 8 — bottom shows the frequency response of transducer elementsonly (without backing). By comparing these ﬁgures it is evident that using of the backing intransducers is necessary.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1012345U[mV]cylindercone 10cone 15cone 20cone 25cone 30cone 35cone 40cone 450 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 105101520f [MHz]U[mV]cylindercone 10cone 15cone 20cone 25cone 30cone 35cone 40cone 45Fig. 8. Frequency responses of transducers (top) and elements (bottom)20O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–240 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 11.41.51.61.71.81.922.12.2f [MHz]U[mV]ABCDABCDFig. 9. Inﬂuence of backing geometry on the frequency responseThe ﬁg. 9 shows inﬂuence of backing shape on the frequency response. The backing geom-etry A, B, C and D is given in tab. 1 — right. The piezoelectric element with the “cone 30”(tab. 1 — left) was used in all cases. Note, the radius of backing inﬂuences more the frequencyresponse of transducer than the height of backing.3.2. Time dependence analysisThe maximal size of one mesh element was 0.25 × 0.25mm. No damping and backing D wasused in the time dependence problems. The time stepping parameters were: times (from 0 to50 with step 0.01μs), relative tolerance (10−5) and absolute tolerance (10−10). We studied theresponse of voltage U (U = Vint/(πr22)) on boundary4 (see ﬁg. 7) as well as for the frequencyresponse analysis.1pmtime0.5μsdisplacementFig. 10. A sinus — pulse functionThe prescribed displacement on boundary2 (seeﬁg.7) was given by a function or a data ﬁle. The sinus-pulse-displacement function had amplitude 1pm, fre-quency 2MHz and width 0.5μs, see ﬁg. 10.A relative frequency response can be determinedby dividing the Fourier transform of the time wave-form produced by the transducer by the Fourier trans-form of the prescribed time waveform, see ﬁg. 11.The frequency response obtained by the Fourier trans-form (ﬁg. 11 — bottom) is in a good agreement withthat obtained by 1D analysis (ﬁg. 6).The data ﬁle with the prescribed displacement onboundary2 (see ﬁg. 7) included the vertical displacement on surface of half-space in distance2cm from the step-point-force source, [6]. The force amplitude was 10N.21O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–24-2-1.5-1-0.500.511.52U[mV]cylinder0 5 10 15 20 25 30 35 40 45 50-2-1.5-1-0.500.511.52t [μs]U[mV]cone 300 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 28.599.51010.5f [MHz]log10|U/ξ|[V/m]cylindercone 30Fig. 11. The pulse time response of cylinder (top), cone (middle) transducer and their frequency re-sponses (bottom). Thin line in bottom image notes mean value of frequency response for the “cone 30”transducer22O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–242.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5-12-10-8-6-4-20246t [μs]U[V]theory cylinder cone 10 cone 30P S R15 xFig. 12. The theoretical waveform and the pulse time response of cylinder, “cone 10”, “cone 30” trans-ducerThe time response of cylinder, “cone 10” and “cone 30” transducer on the step-point-forcesourceareshowninﬁg.12. Thetheoreticaldisplacementwaveformmultipliedbythetransducersensitivity (1.9 · 109V/m) is added to these time responses. The value 1.9 · 109 was taken fromﬁg. 11 — bottom as a mean value of sensitivity for the “cone 30” transducer (log10 1.9 · 109 . =9.28). As shown, the “cone 30” transducer signal is almost exactly as predicted by theory.The letters P, S and R in the ﬁgure noted the arrival time of longitudinal (primary), transverse(secondary) and Rayleigh wave, respectively.4. ConclusionAnalyses of the properties of conical piezoelectric acousticemission transducer were presented.The conical piezoelectric transducer is used to scanning normal displacement on a surface oftesting constructions in NDT methods, such as acoustic emission. The simpliﬁed 1D analysisbased on the equivalent circuit of conical waveguide was recapitulated. FEM analysis of thistransducer was presented. Both the frequency response characteristic of the transducer and histime response on pulse excitation were studied. All of computations were done for variouscone angles of piezoceramic (PTZ-H5) element and various sizes of the cylindrical backing oftransducer.FEM calculations were performed in the environment COMSOL Multiphysics [1] withMATLAB [10] (for scripting of COMSOL tasks).The 1D analysis explained the normal mode degeneracies of the usual cylindrical piezoelec-tric disk element. The FEM analysis showed necessity of using the large transducer backing.The radius of backing inﬂuences more the frequency response of transducer than the heightof backing. It was shown in the time dependence analysis that the ”cone 30” transducer signalis almost exactly as predicted by theory.23O. ˇ Cerven´ a et al. / Applied and Computational Mechanics 2 (2008) 13–24A transducer for measuring tangential dynamic displacement of a transient nature at a pointlocation on the surface of a mechanical body was patented by Proctor [8]. Our next goal will bethe FEM analysis of this transducer. Further we consider supplementing FEM analysis of theconical transducer by horizontal excitation. Both of these tasks require full 3D analysis, whichposes considerable claims of the computer technique parameters.AcknowledgementsThe work was supported by the Czech Science Foundation under the grant 101/06/1689 andthe research project AV0Z20760514 of AS CR.References[1] COMSOL, Inc., http://www.comsol.com[2] S. D. Glaser, G. G. Weiss, L. R. Johnson, Body waves recorded inside an elastic half-space by anembedded, wideband velocity sensor, J. Acoust. Soc. Am. 104 (3) (1998) 1404–1412.[3] M. Greenspan, TheNBSconical transducer: Analysis, J. Acoust. Soc. Am.81(1) (1987) 173–183.[4] Y. C. Lee, Z. Lin, Miniature piezoelectric conical transducer: Fabrication, evaluation and applica-tion, Ultrasonics 44 (2006) e693–e697.[5] A. Mohamed, Equivalent circuits of solid horns undergoing longitudinal vibration, J. Acoust. Soc.Am. 38 (1965) 862–866.[6] H. M. Mooney, Some numerical solutions for Lamb’s problem, B. Seism. Soc. Am. 64 (1974)473–491.[7] T. M. Proctor Jr., An improved piezoelectric acoustic emission transducer, J. Acoust. Soc. Am.71 (5) (1982) 1163–1168.[8] T. M. Proctor Jr., Transducer for measuring transient tangential motion, United States Patent,Patent Number: 4782701.[9] J. Sebastian, Monitoring of refractory wall recession using high-temperature impact-echo instru-mentation, UDRI University of Dayton, Research institute, Dayton, 2004.[10] The MathWorks, Inc., http://www.mathworks.com[11] P. D. Theobald, Towards traceable calibration of acoustic emission measurement systems: devel-opment of a reference source at the UK’s National Physical Laboratory, DGZfP-Proceedings BB90-CD, DGZfP-Proceedings BB 90-CD, EWGAE, 72, 2004, pp. 683–690.24

Analysis of the conical piezoelectric acoustic emission transducer

Abstract

Similar works

Full text

Available Versions

Directory of Open Access Journals