NDT.net • Aug 2004 • Vol. 9 No.08 
Tomographic Imaging of an Ultrasonic Field in a Plane using a Linear ArrayPatrick H. JohnstonNASA Langley Research Center, Hampton, VA 23681 Kendall R. Waters Department of Physics, Washington University, St. Louis, MO 63130 (Present affiliation NIST, Boulder, CO.w v) ABSTRACTA method is presented based on tomographic reconstruction techniques for the visualization of an ultrasonic field. The tall, narrow piezoelectric elements of a onedimensional receiving array, rotated about its central axis, provide parallel projection input data for tomographic reconstruction of the field. The technique is demonstrated using simulated ultrasonic beams and experimental measurements of an ultrasonic beam propagated through a stitched woven composite.INTRODUCTIONResearch is producing increasingly complex materials for use in aerospace and other applications. These advanced material architectures include such developments as multidimensional fiberreinforced plastics incorporating textiles and stitching [1], materials with embedded sensors and actuators [2], and other multifunctional materials concepts. A feature common to all of these “smart materials” is that they exhibit increased spatial variability of material properties, requiring an increased number of parameters to describe their state. These spatial variations occur on a scale significant to ultrasonic waves at the frequencies used for NDE. Because of this informationrich propagation environment, an interrogating ultrasonic field becomes phasedistorted as it gathers material information along its path. Conventional singleelement ultrasonic sensors phasesensitively average the received ultrasonic field. The phasedistortions within the field thus manifest themselves as a kind of coherent “noise” on the received signal. A conclusion can be reached that increasingly complex materials require increasingly sophisticated NDE. Multipleelement ultrasonic arrays can be used as a means to provide additional degrees of measurement freedom. By exploiting the additional degrees of freedom, it may be possible to measure and interpret the distorted ultrasonic field emerging from a complex material, effectively converting the “noise” into “signal” for more complete quantitative NDE of the material. Within the context of exploiting these additional degrees of freedom, the concept was identified of using a rotated linear array receiver to perform tomographic reconstruction of an ultrasonic beam.
CONCEPTConsider the linear array depicted in Fig. 1A. The elements of the linear array have height h and width w, with centertocenter spacing p. The ultrasonic field incident on the array has magnitude and phase characteristics, and can be expressed in the form
Here, the complex spatial variation has been separated into real and imaginary parts. The voltage V(u_{n}, t) produced by the n^{th} element, centered at u_{n}, is proportional to the integral of the pressure over the active area of the element,
If the pressure variation in the u direction over a distance w is small, the u integral can be performed trivially, yielding a line integral of pressure over v,
This can be viewed as a projection integral of the pressure along v at the horizontal location u n . The physical array elements are parallel and evenly spaced along u, and thus produce a set of voltages which represent parallel projection integrals of the incident pressure. If the field is now considered to be fixed in a reference frame (x, y) and the linear array lies in a reference frame (u, v) rotated by angle m relative to (x, y), as depicted in Fig. 1B, the coordinate frames are related according to
The set of voltages {V(u_{n}(q_{m}), t); n = 1, 2, … , N; m = 1, 2, …, M}, acquired at each of the N array elements at each of M angles, equally distributed between 0 and , represent a complete set of data input for tomographic reconstruction by filtered backprojection of parallel projection data ‡ The filtered backprojection algorithm can thus be applied to {V(u_{n}(q_{m}),t)} in order to reconstruct the pressure field [p_{R} (x, y) + i p_{I} (x, y)]. TOMOGRAPHIC RECONSTRUCTION  SIMULATIONSimulated Ultrasonic Beam in WaterSimulated ultrasonic beams were computed using the angular spectrum decomposition technique [3]. Fig. 2 shows the imaginary part of a simulated 2.25 MHz ultrasonic beam produced in water by a 12.7 mm diameter, 50 mm focal length circular transducer, at a distance of 50 mm. Superimposed on the image are the outlines of the elements of a simulated linear array transducer. (For the sake of visualization, the figure contains only 32 elements which are twice the width of the actually simulated ones.) The field was calculated over a 25.6 mm square area, at a pixel resolution of 0.1 mm (256 x 256 pixels). The 64 simulated array elements were 10.00 mm (100 pixels) tall by 0.20 mm (2 pixels) wide, with a gap of 0.10 mm (1 pixel) between elements. The signal from each element was computed by summing all the pixels within the element. The process was repeated for the real and imaginary parts of the simulated beam. Because of symmetry, the resulting set of data for all rotation angles would be identical, so the data were reproduced to simulate measurement at 64 angles, spaced by 5.625° increments. Filtered backprojection was performed to reconstruct the real and imaginary parts of the simulated beam. From the real and imaginary parts, the magnitude and phase were computed, and a twodimensional phase unwrapping algorithm was applied.
‡ The filtered backprojection algorithm is also known as the Radon transform. Some commercial math packages include an implementation of the Radon transform [4,5]. A comprehensive development of the technique is given by Kak [6].
Reconstructed FieldThe log magnitude and the unwrapped phase [7] of the simulated ultrasonic beam is presented on the lefthand panels of Fig. 3. On the righthand side are the log magnitude and unwrapped phase of the reconstructed ultrasonic beam. The field is not reconstructed beyond the area swept out by the rotating array. The reconstructed data match very well with the original field.TOMOGRAPHIC RECONSTRUCTION  EXPERIMENTData acquisitionThe data acquisition configuration is presented in Fig. 4. In a water tank, a 64 element linear array (2.25 MHz, 10 mm by 19.2 mm area) was mounted vertically on a turntable, centered on the axis of rotation. Coaxial with the array, a singleelement circular transmitting transducer (2.25 MHz, 12.7 mm diameter, 50 mm focal length) was mounted.
The two were separated by 55 mm to allow the placement of a specimen at the focal distance of the transmitting transducer. The transmitting transducer and the receiving array were connected to the NASA Langley Ultrasonic Phased Array Testbed System (UPATS). The transmitter was driven by a 10cycle tone burst at 2.25 MHz. The signal received by each of the 64 array elements was digitized (50 MHz sampling rate, 1024 points) and a discrete Fourier transform was performed to extract the real and imaginary parts at 2.25 MHz. The measurement was repeated for each of 64 rotational positions of the turntable separated by 5.625°. This measurement procedure was followed first with only water between the transducers. Subsequently, a stitched composite specimen was positioned between the transducers, and the measurement procedure was repeated at each of two locations on the specimen, to be described below. Stitched composite specimenFig. 5A is a crosssectional sketch of the fiber architecture in a stitched fiberreinforced composite material. The composite is built up of layers of fibers, depicted here as straight, but may be mats of woven, braided, or knitted fiber bundles. The layers are held together by a sewing stitch, which consists of a doubled length of fibers penetrating through the thickness of the composite. Under a typical standard ultrasonic inspection, an ultrasonic beam would be incident perpendicular to the surface of the part. The stitch fibers present a locally increased stiffness, and hence wave speed, which results in wavefront distortion, and the stitch fibers also provide additional scattering and other loss mechanisms.In Fig. 5B is a transmission Cscan of a stitched composite material, 10 cm on each side. The textile layers are woven carbon fibers, the stitch material is Kevlar, and the matrix is epoxy. At the locations of the throughthickness stitches, the signal loss is greatest, producing dark spots in the image. Also, the woven fabric produces variations in the background image. For this study, an ultrasonic beam was propagated through a composite location between stitches, where the woven textile architecture alone was interrogated, and through a location directly over a stitch. As a dimensional reference, the outline of a reconstructed beam area is illustrated by a broken square in Fig. 5B.
RESULTSReconstructed beam through waterFig. 6 presents the reconstructed field derived from data taken through a wateronly path. From left to right, the images are of the real part and the imaginary part of the pressure, followed by the log magnitude and the unwrapped phase derived from the first two quantities. The real and imaginary part data have both positive and negative values. Contours are superposed to indicate the zero crossings, and + and – symbols indicate the sense of the data between the contours. Unwrapped phase has contours at multiples of . In this, as in subsequent images of unwrapped phase, there are some phase artifacts around the outer portion of the image, where the signaltonoise ratio is smaller. While there is evidence in these results of an angular misalignment between the transmitting transducer and the plane of the linear array, the error is small and the circular symmetry of the beam is predominantly preserved.Reconstructed beam through woven compositeIn Fig. 7 are presented the reconstructed images for a beam transmitted through a woven composite between the throughthethickness stitches. As before, contours are shown marking zero in the real and imaginary parts, as well as changes of in the phase image. Comparison of these images with those in Fig. 6 shows that the beam has been distorted by propagating through the woven composite, not simply attenuated uniformly.Reconstructed beam through a stitch in a woven compositeFig. 8 shows the reconstructed images for a beam transmitted through a woven composite at the location of a throughthethickness stitch. Contours are indicated as in the preceding two figures. Comparing these results with those in Figs. 6 and 7, one can observe a data “dimple” near the center of the main lobe of the beam, which should correspond to effects caused by the Kevlar stitch. In addition to the localized “dimple”, there is a general distortion similar to that seen in Fig. 7, caused by the woven textile fibers.
DISCUSSIONModern aerospace materials and emerging experimental materials possess higher spatial variability in their elastic properties. Because these variations occur on a spatial scale similar to the wavelengths used for ultrasonic NDE, the interaction between an interrogating ultrasonic beam and these material inhomogeneities produces distortions in the beam. In a conventional ultrasonic inspection, the diameter of the receiving aperture might typically be just less than the areas reconstructed above. These receiving apertures phasesensitively integrate the real and imaginary parts of the field over their areas, and the beam distortions, such as those displayed in Figs. 89, are lost. Consider the log magnitude results from Figs. 68, plotted again in Fig. 9 with the grayscale compressed to highlight only the largest 15 dB of the beam, which is concentrated in the main sound beam where 7080% of the energy is found. There is an apparent structure to the distortion caused by the woven textile fibers, in the center image. In the righthand image, one can interpret the beam as having a similar distortion as in the center image, rotated slightly counterclockwise, with an additional localized loss caused by the stitch. In a conventional ultrasonic inspection, these three images would be reduced to three values of a single parameter, representing the three integrated amplitudes.
It clearly takes more than one parameter to characterize the actual differences among these three ultrasonic beams, and this additional information may be valuable for the quantitative characterization of the composite. In this paper, a novel method has been introduced which enables the highresolution imaging of a compact ultrasonic field, such as a transmitted beam. This method provides a means to characterize the spatial information carried by the beam after passing through a complex material, such as a stitched, woven composite. ACKNOWLEDGMENTThis work was supported in part by a NASA Graduate Student Research Fellowship NGT152197. The authors thank Dr. Michael Johnson of Lawrence Livermore National Laboratory for the use of his code for unwrapping of the twodimensional phase. We thank Jeffrey Boghosian for technical assistance with the UPATS data acquisition.REFERENCES

© NDT. net  Top 