#### Abstract

The improvement of ultrasonic phased array detection technology is a major concern of engineering community. Orthotropic piezoelectric fiber composite (OPFC) can be constructed to multielement linear array which may be applied conveniently to actuators and sensors. The phased array transducers can generate special directional strong actuator power and high sensitivity for its orthotropic performance. Focusing beam of the linear phased array transducer is obtained simply only by adjusting a parabolic time delay. In this work, the distributed point source method (DPSM) is used to model the ultrasonic field. DPSM is a newly developed mesh-free numerical technique that has been developed for solving a variety of engineering problems. This work gives the basic theory of this method and solves the problems from the application of new OPFC phased array transducer. Compared with traditional transducer, the interaction effect of two OPFC linear phased array transducers is also modeled in the same medium, which shows that the pressure beam produced by the new transducer is narrower or more collimated than that produced by the conventional transducer at different angles. DPSM can be used to analyze and optimally design the OPFC linear phased array transducer.

#### 1. Introduction

The ultrasonic phased array detection technology which was proposed by Crawford of America Conoco Company originated in the 1960s [1]. The phased array transducer is formed by connecting different components and space distribution of piezoelectric elements. The structure can be detected with rapid multi-angles by controlling emitting order and time delay. These arrays have been also utilized in the field of nondestructive evaluation (NDE) of materials [2, 3]. But conventional phased array transducers have some deficiencies such as brittle fracture, different design ability, and illegitimate mechanical-electrical parameters [4] with piezoelectric material as core element, and the final resolution will be affected. Therefore, the ultrasonic phased array detecting technology can be further improved by improving the quality of key materials and developing new phased array transducer. Compared with traditional piezoceramics, orthotropic piezoelectric composite materials (OPFC) which are used as actuator/sensor elements in damage detection for engineering structures show the obvious advantages for their excellent orthotropic characteristic [5, 6].

The group of Luo [7] has developed models of actuators and sensors with OPFC to satisfy the detection of different structures. These models of sensors will be used to fabricate new phased array transducers in many different engineering structures inspection. OPFC phased array is treated as a linear array consisting of single element sources. The new structure has orthotropic performance. This special performance can be used to actuate and receive the stress wave in different directions and reduce the interfering signals.

In this paper, new OPFC phased array transducers are semianalytically modeled by the distributed point source method (DPSM). The DPSM is used to model ultrasonic transducers for computing pressure and velocity fields generated by the new transducers. The group of Placko [8] first developed the new method for ultrasonic transducer simulation. In the present paper, an OPFC phased array transducer is modeled using DPSM to obtain the ultrasonic field in a homogeneous medium. The results compare the phased array transducer and conventional planar transducer generated ultrasonic beam in different directions. The importance of the interaction effect between two phased arrays is also considered by computing the ultrasonic field.

#### 2. DPSM Theory Analysis

Based on a spatial distribution of point source, DPSM modeling can be used to compute the pressure field generated by an ultrasonic transducer which is placed in homogeneous or nonhomogeneous fluid media. A traditional transducer composed of many point sources which oscillate at the same time. Equation (1) shows the total field at a point generated by number of point sources, but the point sources are excited at different times for a phased array transducer. Let the th point source of strength be excited at time; the pressure at the distance from a specific point source is given in (2): where is the angular frequency and is the wave number. For the phased array transducer, the point sources on the surface are assumed to be excited after different time intervals and can be expressed as where is a reference time . If there are point sources distributed over the element surface, as shown in Figure 1, the total pressure at point can be obtained: where is the distance of the th point source from . From the pressure-velocity relation [9], the velocity in direction which is generated by the th point source is given by

After adding the contributions of sources, the velocity components in directions , and (see Figure 1) are obtained:

Because the transducer surface vibrates with a nonzero velocity in direction only, on the surface, and . For traditional transducers, and for phased array transducers, . If the observation point is placed on the surface of transducer, the three expressions of equation (6) should be equal to 0, 0, and , respectively. Thus, equations will be obtained from points on the transducer surface. We get more equations than unknowns. To get the same number of unknowns as equations, the number of unknowns can be increased from to by replacing each point source by a triplet source. All sources are placed on the same plane at parallel to the transducer surface. Random orientations of triplet sources are necessary to preserve the isotropic material properties. By solving linear equations, the source strength of all point sources can be obtained. After getting , the pressure can be calculated at any point from (1). The velocity components parallel to the transducer face are not necessarily continuous because slippage may occur between transducer face and the fluid. Equation (6) can be combined into the following matrix equation:
where is the vector of the velocity components at number of surface points *. * is the vector of the strengths of number of point sources. is the matrix relating the two vectors and . So can be written as

Vector of the source strengths is given by

Finally, the square matrix can be obtained:where

In the above equations stands for the distance between the th point source and th aim point of the transducer surface. The first subscript of of can take 1, 2, or 3 and indicate the , , or directions. The subscript of and can take values from 1 to depending on the considered point source and the superscript can take values from 1 to *.* As mentioned formulation from point sources, three boundary conditions of velocity are satisfied to every point of boundary points. From the point source strengths by inverting matrix ,

Then, the pressure and velocity components , at the observation points, are obtained from the following matrix relation:

If the observation points are identical to the transducer surface points, matrix will be the same as of (11). However, by computing the velocity field at different points, will be as shown below:

When two transducers are placed in the same fluid, the pressure field can be modeled by superimposing two simpler parts, as shown in Figure 2. The total field at any point can be looked as part 1 and part 2. The right face is replaced by an imaginary scattering surface in part 1 and the left face is replaced by the imaginary scatterer. The two problems can be modeled as distributed point sources. The total ultrasonic field at point generated by the two OPFC phased array transducers can be obtained by combining the above equations: where and are the ultrasonic fields generated by part 1 and part 2, respectively. is the total effect for the two transducers. Note that takes into account the interaction effect. is the pressure field by the original transducer.

#### 3. The Pressure Field of OPFC Element in a Homogeneous Fluid

The DPSM theory presented above is used to model ultrasonic field generated in a homogeneous fluid (water) by OPFC ultrasonic phased array transducer. The near field length and the divergence angle can be calculated in the method. Compared with the closed form analytical values, we can validate the accuracy of DPSM. The near field length and divergence angle equation can be shown as follows [10]: where is the wave length, is the half length of section, and is the rectangle area of the transducer, as shown in Figure 3. The analytical expression of the pressure in a homogeneous fluid along the central axis of a transducer is given in (17):

Figure 4 shows the variations of acoustic pressure along -axis by DPSM and analytical expression. The signal frequency MHz. The dimension of the transducer is 21.6 mm × 3 mm 0.8 mm, which is the optimal dimension to design the OPFC element. For the DPSM modeling 80 points are placed at , while the transducer surface is at . According to the condition of DPSM geometry problem, . From (16) is found to be 0.51 mm. As shown in Figure 3, the theoretical value of matches very well with the transition point.

Figure 6 shows the corresponding acoustic pressure variation in plane in DPSM technology. The theoretical divergence angle is calculated from (17). The measured divergence angle from Figure 5 is about 30°, which is very close to the theoretical value. The two results validate the accuracy of DPSM.

#### 4. Numerical Simulation Results of OPFC Phased Array Transducer by DPSM

##### 4.1. Dynamic Focusing and Time Delay Calculation

Focusing beam is accomplished by combining a spherical timing relationship to produce a beam, which is focused at a given range and propagates at a specific angle. The focusing delays can be calculated by the following traditional formula [11]: where is the required delay for the th element, is the focal length, and is a constant to keep the delays positive. The focusing has two limitations as it stands. Firstly, it is only valid for an odd number of elements. However, most phased arrays have an even number of elements. Secondly, the constant to keep the delays positive is very clumsy to utilize. In order to overcome these deficiencies, a generalized formula for handling both even and odd number of elements was derived. As shown in Figure 7, the following geometric relationship can be attained: for the element , where is the total number of elements. is the actuate time of the th element. Solving for , we obtain

Consider Taylor expression,

Substituting (20) into (21), the focusing delay of adjacent arrays can be obtained: and .

So,

The delay of the th element and the first element can be written as

OPFC phased array transducer is composed of many elements arranged in a certain pattern that emit acoustic energy at different times. The elements can be pulsed in certain sequence to control the beam angle depending on the time delay.

##### 4.2. Numerical Simulation Results and Discussion

OPFC phased array transducer can easily steer an ultrasonic beam in different directions without rotating or moving. But the traditional transducers are rotated to obtain the desired steering direction. Figure 8 shows the steering sketch of phased array transducer and conventional transducer rotated.

**(a)**

**(b)**

Figure 9 shows the acoustic pressure field generated by an OPFC phased array transducer in water. The length and width of the transducer side are 5 mm and 0.8 mm, respectively. And the exciting frequency is 2 MHz. The acoustic pressure fields are modeled for steering angles equal to 0°, 15°, 30°, and 45°, respectively. Figures 10(a)–10(d) show the pressure field generated by eight single OPFC sensor elements when the axis of the transducer is rotated by 0°, 15°, 30°, and 45°, respectively. For both of these results, -axis stands for the coordinate in length direction of OPFC phased array transducer and -axis for the coordinate in width direction. The surface area is equal to 4.0 mm^{2} and the number of point sources is 80. Each OPFC element is discredited into five point sources distributed along the direction, giving rise to a total of 60 point sources. In Figure 9, the distance of the focal point is set at 5 mm in the direction. The normal distance in direction kees stable at the location of 5 mm as the steering angles changes. From Figures 9 and 10, it can be seen that the maximum pressure intensities are and units for phased array and conventional transducers, respectively. It is also observed that the pressure beam produced by the OPFC phased array transducer is narrower or more collimated than that produced by the traditional transducer. Similar conclusions can be drawn from Figures 9(b) and 10(b) for 15° steering angle, Figures 9(c) and 10(c) for 30° steering angle, and Figures 9(d) and 10(d) for 45° steering angle.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

Figure 11 is the acoustic pressure along the axis by phased array and traditional transducers for steering angles equal to 0°, 15°, 30°, and 45°, respectively. From the pictures, the OPFC phased array transducer produces stronger pressure beam than the traditional transducer in the near field along steering direction. The simulation shows that the pressure excited by OPFC phased array transducer is stronger than that by traditional transducer between the near field distance and the focal point .

**(a)**

**(b)**

**(c)**

**(d)**

The interaction effect between two OPFC phased array transducers placed in the same medium is studied. Two transducers faces are placed parallel to each other and the central axes are collinear. The acoustic field generated by the one end as well as the scattered fields by the other end transducers is calculated. The surrounding medium is water. Figures 12(a) and 12(b) show the pressure fields generated by the phased array transducers at the steering angle 0° when the interaction effect is ignored. Figure 13 shows the pressure field by the same OPFC phased array transducers when the interaction effect is included. The pressure beam generated by the bottom transducer is partly reflected or scattered by the top transducer and affects the total acoustic field. The interaction effect can be seen clearly. The reflected beam from the top transducer caused peak value is much greater than that observed in Figure 12. In Figure 13, the beam at near mm is slightly wider than that in Figure 12(b).

**(a)**

**(b)**

#### 5. Conclusion

The fabrication principle of OPFC linear phased array is demonstrated for the special orthotropic performance. An OPFC linear phased array transducer with dynamic focusing is modeled by using DPSM. DPSM technology was utilized to compare the behavior of focusing in different angles of OPFC linear phased array and traditional transducer. Two linear phased array transducers placed face to face in the same medium are also modeled to study modeled to study the interaction effect. This demonstrated the importance of focusing in the near fields, while the directivity for focusing is very well defined. The OPFC linear phased array transducer can produce stronger and better collimated beams to detect the damages.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors thank the support of “863 Project” (2009AA03Z107), Natural Science Foundation of China (11272138), Jiangsu University Foundation (14JDG022), Doctor Point Foundation Project (20123227130002), and their colleagues for their contribution to the work.