Modeling the Radiation Field of Normal-Beam Ultrasonic Transducers by Multi-Gaussian Beam Technique

Similar documents
UNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

Analysis of composites with multiple rigid-line reinforcements by the BEM

INF-GEO Solutions, Geometrical Optics, Part 1

Solutions to Final Exam Review Problems

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Types of Waves Transverse Shear. Waves. The Wave Equation

The Method of Least Squares. To understand least squares fitting of data.

Numerical Methods in Fourier Series Applications

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

Warped, Chirp Z-Transform: Radar Signal Processing

Ma 530 Introduction to Power Series

Curve Sketching Handout #5 Topic Interpretation Rational Functions

Chapter 9: Numerical Differentiation

MATH 31B: MIDTERM 2 REVIEW

Deterministic Model of Multipath Fading for Circular and Parabolic Reflector Patterns

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

Math 116 Second Exam

Chapter 2 Feedback Control Theory Continued

The Relative Angle Distribution Function in the Langevin Theory of Dilute Dipoles. Robert D. Nielsen

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

Review Problems for the Final

Formation of A Supergain Array and Its Application in Radar

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =

There is no straightforward approach for choosing the warmup period l.

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Math 113, Calculus II Winter 2007 Final Exam Solutions

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

GROUND MOTION OF NON-CIRCULAR ALLUVIAL VALLEY FOR INCIDENT PLANE SH-WAVE. Hui QI, Yong SHI, Jingfu NAN

Free Space Optical Wireless Communications under Turbulence Channel Effect

CALCULATION IN THE FIELD OF SEGMENTAL ROTOR MACHINES TAKING INTO ACCOUNT WINDING HARMONICS AND ROTOR AIRGAP IRREGULARITIES

Characterization of anisotropic acoustic metamaterials

FIR Filter Design: Part II

DISTRIBUTION LAW Okunev I.V.

THE LEVEL SET METHOD APPLIED TO THREE-DIMENSIONAL DETONATION WAVE PROPAGATION. Wen Shanggang, Sun Chengwei, Zhao Feng, Chen Jun

Notes The Incremental Motion Model:

Rotationally invariant integrals of arbitrary dimensions

Random Variables, Sampling and Estimation

6.3 Testing Series With Positive Terms

U8L1: Sec Equations of Lines in R 2

EXPERIMENT OF SIMPLE VIBRATION

COMM 602: Digital Signal Processing

PROPOSING INPUT-DEPENDENT MODE CONTRIBUTION FACTORS FOR SIMPLIFIED SEISMIC RESPONSE ANALYSIS OF BUILDING SYSTEMS

EVALUATION OF GLASS FIBER/EPOXY INTERFACIAL STRENGTH BY THE CRUCIFORM SPECIMEN METHOD

CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS

Section 11.8: Power Series

THE KALMAN FILTER RAUL ROJAS

Chapter 6 Sampling Distributions

SOLID MECHANICS TUTORIAL BALANCING OF RECIPROCATING MACHINERY

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

1 6 = 1 6 = + Factorials and Euler s Gamma function

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

Vector Quantization: a Limiting Case of EM

Probability, Expectation Value and Uncertainty

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Solving equations (incl. radical equations) involving these skills, but ultimately solvable by factoring/quadratic formula (no complex roots)

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

577. Estimation of surface roughness using high frequency vibrations

A numerical Technique Finite Volume Method for Solving Diffusion 2D Problem

Research on real time compensation of thermal errors of CNC lathe based on linear regression theory Qiu Yongliang

Obtaining Constants of Johnson-Cook Material Model Using a Combined Experimental, Numerical Simulation and Optimization Method

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

MATH 2411 Spring 2011 Practice Exam #1 Tuesday, March 1 st Sections: Sections ; 6.8; Instructions:

Math 113 Exam 3 Practice

Quantum Annealing for Heisenberg Spin Chains

REFLECTION AND REFRACTION

Monte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem

ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION

Chapter 4. Fourier Series

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Sequences, Series, and All That

Chapter 4 : Laplace Transform

Chapter 7 z-transform

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem

Optics. n n. sin. 1. law of rectilinear propagation 2. law of reflection = 3. law of refraction

[ ] sin ( ) ( ) = 2 2 ( ) ( ) ( ) ˆ Mechanical Spectroscopy II

9.3 The INTEGRAL TEST; p-series

Streamfunction-Vorticity Formulation

Outline. Linear regression. Regularization functions. Polynomial curve fitting. Stochastic gradient descent for regression. MLE for regression

Balancing. Rotating Components Examples of rotating components in a mechanism or a machine. (a)

6.1 Analysis of frequency selective surfaces

On Random Line Segments in the Unit Square

Orthogonal Gaussian Filters for Signal Processing

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

mx bx kx F t. dt IR I LI V t, Q LQ RQ V t,

CATHOLIC JUNIOR COLLEGE General Certificate of Education Advanced Level Higher 2 JC2 Preliminary Examination MATHEMATICS 9740/01

PHYS-3301 Lecture 9. CHAPTER 5 Wave Properties of Matter and Quantum Mechanics I. 5.3: Electron Scattering. Bohr s Quantization Condition

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

Math 116 Practice for Exam 3

MATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS

Probability of error for LDPC OC with one co-channel Interferer over i.i.d Rayleigh Fading

Information-based Feature Selection

The Advection-Diffusion equation!

Characterization and Simulation of HF Propagation in a realistic disturbed and horizontally inhomogeneous Ionosphere

Physics 116A Solutions to Homework Set #9 Winter 2012

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

PILOT STUDY ON THE HORIZONTAL SHEAR BEHAVIOUR OF FRP RUBBER ISOLATORS

The Scattering Matrix

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

Transcription:

Modelig the Radiatio Field of Normal-Beam Ultrasoic Trasducers by Multi-Gaussia Beam Techique Ami Yaghootia, Farhag Hoarvar, ad Mehdi Zeighami NDE Lab., Faculty of Mechaical Egieerig, K. N. Toosi Uiversity of Techology, ardis St., Molasadra Ave., Vaak Sq., Tehra, Ira. Abstract The modelig of ultrasoic tests ca help i gaiig a better uderstadig of the test process. The model ca also predict the test results ad provide a basis for choosig the optimum parameter settigs. Oe of the first steps i modelig a ultrasoic test is to uderstad the characteristics of the ultrasoic wave field. I this paper, the multi-gaussia beam () techique is used for modelig the radiatio field of a ormal-beam ultrasoic trasducer. By usig this approach, oe ca avoid the computatioal difficulties ecoutered i other modelig techiques such as Rayleigh-Sommerfeld ad Gree fuctio itegrals. By modelig the radiatio field of a MHz trasducer, the beam profile of the reflected sigal have bee calculated. Experimets cofirm the computatioal results. Keywords: Ultrasoic Modelig, Multi-Gaussia Beam, Wave Field. Itroductio Modelig of wave field geerated by a trasducer ca help i odestructive evaluatio of materials by predictio of echoes, simulatio of defects, ad also i sigal processig. Moreover, a good model ca be very helpful i the desig ad optimizatio of a testig procedure ad i the iterpretatio of experimetal results. It is also much simpler ad cheaper to perform parametric studies with by usig a good model compared to coductig experimets. By usig the model, the probability of detectio of various types of defects ca be assessed. I the earliest works of wave field modelig the trasducer surface was cosidered to act as a source of poit forces (Cherry 96). This method had computatio difficulties. By usig Johso's Gree fuctio, Tag et. al. (99) derived itegral expressios of the elastic displacemets for both compressio ad shear pisto sources actig o the surface of a elastic half-space. The shortcomig of this method is its iability to model the ear-field of the wave field as well as oblique icidece modelig. Harris (98) reviewed three approximatios icludig Rayleigh surface itegral, Kig itegral, ad Schoch solutio to fid the pressure field of a plaar probe o a ifiite half-space. Other umerical methods were also developed to show the displacemet field i a solid media or the pressure field i a fluid due to wave radiated from ormal or agular trasducers (Alia et. al. 4). Based o a elastodyamics solutio, Bostrom ad Wirdelius (995) modeled the wave field for elliptical ad rectagular trasducers i a ormal or agular positio, emittig wave ito a solid half space. Although the exact o-axis field ad a approximatio to the far field ca be expressed i a closed form, describig the ear field, the trasiet zoe, ad off-axis beams are difficult ad time cosumig i umerical methods by itegratio over poits of the probe surface.

We ad Breazeale (988) proposed a alterative approach. They computed the total field by superimposig a umber of Gaussia beam solutios. They have show that by superimpositio of oly Gaussia solutios, the field radiated by a circular pisto trasducer ca be accurately modeled. I this paper, the multi-gaussia beam techique is used for modelig the radiatio field of a ormal-beam ultrasoic trasducer. By usig this approach, oe ca avoid the computatio difficulties ecoutered i other techiques icludig the Rayleigh-Sommerfeld ad Gree fuctio itegrals. I some ultrasoic tests, it is required to use ultrasoic beams with large spread agles. Based o the locatio of defects i the material, the wave ca be reflected from ay poit o the beam cross sectio. By modelig the radiatio field of a MHz trasducer, the beam profile ad frequecy spectrum of the reflected sigal have bee calculated. Multi-Gaussia Beam Modelig Modelig the fields radiated by ultrasoic trasducers is a challegig task because of the large umber of possible trasducer types, sizes, ad cofiguratios that are used i practice. The multi-gaussia beam model ca describe the wave field of a circular pisto trasducer by superpositio of a umber of idividual Gaussia beams with a proper set of weightig factors. Although ultrasoic trasducers do ot geerate Gaussia beams but this modelig works well i cosiderig the trasducer to act as a pisto source. Oe of the advatages of multi-gaussia beam modelig is sigificat reductio i computatio time as well as ability i modelig the wave field i aisotropic materials. I Gaussia beam modelig of circular trasducers, it is first assumed that the trasducer is a time harmoic Gaussia source with time depedecy ad waves are emitted perpedicular to the trasducer surface ( x directio) as a quasi-plae wave. The displacemet field i Gaussia beam is obtaied from the followig formulatio (Kim et. al. 4), u ( x, x, x [ M ( x )] T det ( ), ω iv M x ik x x ) = exp( ik ) exp x ρc det M () ω () where v is the velocity o the trasducer surface, k = ω / c is the wave umber, ω is frequecy i radias per secod, c is logitudial wave speed i specime, ad ρ is the T desity of the meduium i which the wave radiates. Moreover, x = [ x, x ] ad M is a complex-valued symmetric matrix. The multi-gaussia beam modelig ca fially be writte as (Schmerr ad Sog 7), u ( x, x, x [ M ( x )] N T ia iωx x, ) = exp( )exp ω ik x + () = ib x / DR where N is the umber of Gaussia beams, DR = ka / is the Rayleigh distace, a is the radius of the trasducer ad

ib / c DR + ib = x / DR M () ib / cdr + ib x / DR [ ( x )] A ad B are complex-valued expasio coefficiets that eed to be determied to match the velocity field o the face of the trasducer. We ad Breazeale (988) foud te coefficiets by a optimizatio method for circular plaar pisto trasducers. Experimet The results of multi-gaussia beam model have to be verified by compariso with other models or with experimets. A experimet was carried o a 69 mm thick alumium block. The compressioal ad shear velocity i alumium are 6 m/s ad 5 m/s, respectively. The wave was geerated by a 6.5 mm diameter trasducer havig a ceter frequecy of MHz. The test was performed by ultrasoic pulse-echo techique i which the sigle trasducer acted as both the trasmitter ad receiver. Results ad Discussio I the first step, to show the ability of multi-guassia beam for modelig both the ear field ad far field of the trasducer, the computed pressure field of a circular trasducer with?? mm diameter ad 5 MHz cetral frequecy is illustrated i Fig.. Fig. - Modelig of the trasducer wave field wave field The backwall echo obtaied from the alumium test block is show i Fig a. The frequecy spectrum of this echo is also show i Fig. b. Both the multi-gaussia model ad the Gree fuctio mdoel were used for modelig the backwall echo. The Gree fuctio model was developed followig the approach used by Tag et. al. (994). The modeled echo obtaied from this approach is show i Fig. c alog with its frequecy spectrum show i Fig. d. The echo obtaied from the multi-gaussia model is show i Fig. e ad its reletive spectrum frequecy is show i Fig. f. Both modeled echoes are i good agreemet with the experimetal result.

(a) (b) 5 4 5 (c) (d) 5 4 5 (e) (f) 5 Time (µs) 4 5 Frequecy (MHz) Fig. - Comparsio of experimetal sigal (a), Gree fuctio model (c) ad multi-gaussia beam modelig (e). Frequecy spectra of three echo is show i (b), (d) ad (f) respectively. Next, we compare the capability of the two models i modelig the radiated ultrasoic wave at various agles. We cosider a 5 MHz trasduce with 6.5 mm diameter to be trasmittig the ultrasoic wave ito a semi-circular specime havig a radius of 75 mm. The wave is the received by a idetical trasducer o the cured surface of the specime as show i Fig.. Fig. - Cofiguratio of experimetal set-up 4

.8 θ = 5 o Tag(994).8 θ = o Tag(994).6.6.4.4.. 5 5 frequecy (MHz) 5 5 frequecy (MHz).8 θ = 5 o Tag(994).8 θ = o Tag(994).6.6.4.4.. 5 5 frequecy (MHz) 5 5 frequecy (MHz) Fig. 4- Frequecy spectrum comparssio of (solid lie) ad Gree's fuctio (dashed lie) As show i Fig. 4, by icreasig i agle with respect to the ceteral axis of the trasmittig trasducer, the frequecy spectrum of the received echoes chage. This chage of frequecy spectra is due to the destructive iterferece of waves origiatig from various poits o the trasducer surface. Agai the two models are i good agreemet. Experimeal verificatio of the two models is curretly uderway. Coclusio I this paper, the multi-gaussia beam () approach was used to model the wave field of a ormal beam trasducer. The advatage of usig this method is modelig the ear-field as well as the far-field of the trasducer i both fluids ad solids. It also avoids the computatioal difficulties ecoutered with other modelig techiques. To examie the capability of this model, a backwall echo was modeled by the multi-gaussia beam approach ad compared with the experimetal as well as Gree fuctio modelig results. The result obtaied from the model was i good agreemet with the experimetal result. Moreover, the frequecy spectra of beam radiated from a ormal trasducer was modeled by the approach ad compared with Gree's fuctio results. Agai, the two models were i good agreemet. The techique ca be used i modelig the refractio ad reflectio from iterfaces without ay computatio difficulties. 5

Referece [] J. T. Cherry, Jr., "The Azimuthal ad olar Radiatio atters Obtaied From a Horizotal Stress Applied at the Surface of Elastic Half Space", Bull. Seismol. Soc. Am. 5, 7-6 (96). [] X. M. Tag, M. N. Toksoz, ad C. H. Cheg, "Elastic wave radiatio ad diffractio of a pisto source", J. Acoust. Soc. Am. 87, 894-9, (99). [] G. R. Harris, "Review of trasiet field theory for a bared plaar pisto", J. Acoust. Soc. Am. 7, - (98). [4] A. Alia, H. Djelouah, N. Bouaoua, "Fiite differece modelig of the ultrasoic field radiated by circular trasducers", J. Comput. Acoust. (4), 475 499, (4). [5] A. Bostrom ad H. Wirdelius, "Ultrasoic probe modelig ad odestructive crack detectio", J. Acoust. Soc. Am. 97(5), 86-848, May (995). [6] J.J. We ad M.A. Breazeale, "A diffractio beam field expressed as the superpositio of Gaussia beams", J. Acoust. Soc. Am. 8 (5), 75-756, (988). [7] H. J. Kim, J. S. ark, S. J. Sog, ad L.W. Schmerr Jr, "Modelig Agel Beam Ultrasoic Testig Usig Multi-Gaussia Beam", Joural of Nodestructive Evaluatio, 4,(),8-9. [8] L.W. Schmerr, S.J. Sog, "Ultrasoic Nodestructive Evaluatio Systems: Models ad Measuremets", Spriger, 7. 6