**Stress Estimation Using the Acoustoelastic E**ff**ect of Surface Waves in Weak Anisotropic Materials**

**Jihyun Jun <sup>1</sup> , Young-Dae Shim <sup>1</sup> and Kyung-Young Jhang 2,\***


Received: 20 November 2019; Accepted: 18 December 2019; Published: 24 December 2019 -

**Abstract:** This paper proposes a novel stress measurement method using the acoustoelastic effect of surface wave to estimate the stress of a homogeneous material plate with orthogonal anisotropy, in which the surface wave velocities are measured in three different directions before and after loading stress. The effectiveness of the proposed method was verified by numerical simulations and experiments. For the simulations, the surface wave velocities in three directions were obtained from a conventional perturbation model for weak anisotropic materials. The simulation results showed that the stress estimation error was less than 3% for an anisotropic rate up to 2% under stress conditions up to 90 MPa. Two specimens were prepared for the experiments, one was almost isotropic and another that had a relatively larger anisotropy rate of 2.6%. Then, the stresses loaded by a tensile test machine were estimated. The results showed good agreement with the given stresses for both specimens. These results confirm that the proposed method can be applied to estimate the surface stress state in anisotropic material plates. The proposed method is simple, practical, and is expected to be useful for monitoring changes of surface stress before and after machining such as the punching or bending of plate.

**Keywords:** acoustoelastic effect; stress estimation; surface waves; weak anisotropy material

#### **1. Introduction**

Monitoring of applied or residual stress is important for quality control of industrial parts manufacturing and health management of structures. The conventional methods of measuring stress include mechanical methods such as hole-drilling and nondestructive techniques such as X-ray diffraction (XRD) and acoustoelastic method. However, mechanical methods are complicated and not suitable for online monitoring. X-ray diffraction has the disadvantage of being harmful to the human body and requires separate pretreatment of the specimen. Unlike these approaches, the acoustoelastic method is harmless and has advantages that it can be applied to real products without pre-treatment and it is suitable for online monitoring.

The acoustoelastic method utilizes the propagation velocity of an elastic wave, which changes with the stress state of a medium [1]. The acoustoelastic effect appears in all kinds of ultrasonic waves including longitudinal, transverse and surface waves. For many applications, however, surface stresses are needed and in this case, surface waves are suitable since they only penetrate to a depth of approximately one wavelength. In addition, the transmitter and receiver of the surface waves are located on the same plane so that it is possible to perform inspection using only one side. Furthermore, it is not necessary to know the exact thickness of the applied object to measure the wave velocity. These features enhance the application of surface waves in the field. Accordingly, various studies on stress estimation using acoustoelastic effect of surface wave have been conducted. In particular, many

studies have been carried out to analyze stresses in weld [2] or rail specimen [3,4]. In those studies, however, the direction of the dominant stress is already known, so it is only possible to determine the presence of the stress or to estimate the stress in the known direction of the dominant stress.

Meanwhile, there is also a method using longitudinal critical refraction (Lcr) waves, which have advantage of high sensitivity to stress [2,5]. However, the Lcr wave is only applicable in contact manner because it uses the critical refraction at the contact interface between the wedge and the test object. On the other hand, surface waves are easy to extend in a non-contact manner such as laser ultrasonic technology.

Early studies regarding the acoustoelastic effect of surface wave were conducted to identify the linear relationships between stress and surface wave velocity in isotropic material [6,7]. However, even if a material is initially isotropic, changes of the micro-structure occur during material production and processing, which may result in anisotropic properties. In this case, no matter how weak the anisotropy is, the use of the isotropic theory causes a large error in stress estimation and therefore cannot be applied.

Only a very limited number of studies have investigated on the acoustoelastic effect in anisotropic material. Delsanto et al. [8] and Mase and Delsanto [9] proposed a model that separates the contribution of the anisotropic effect and the contributions of the acoustoelastic effect to evaluate how the surface wave velocity changes from the isotropic state when a stress is applied in a weakly anisotropic material. Although, this method has not been experimentally verified, it is worthwhile to analyze the effect of anisotropy. However, it is very inconvenient for actual application, because it is necessary to know in advance how the actual elastic modulus of the material differs from the isotropic elastic modulus before the stress is applied. Another study by Thompson et al. [10] proposed a general solution for plane wave propagation in a symmetry plane of an orthorhombic, biaxial stressed anisotropic material. However, it does not provide a direct solution to the case of using surface waves. Additionally, it is necessary to know the direction of anisotropy in advance, which leads to an inconvenient step of measuring the shear wave velocity in many different polarization directions.

In this study, we propose a simpler and more practical method for estimating surface stress using the acoustoelastic effect of surface waves in a plate that has orthogonal anisotropy on its surface. In the proposed method, we modified the isotropic theory, in which the rate of change of surface wave velocity according to the stress is expressed using two acoustoelastic coefficients, two principal stresses and the angle between the principal stress and the wave propagation direction. The difference between the proposed method and isotropic theory is that in isotropic theory, the surface wave velocity at the unstressed state is constant regardless of the propagation direction, while in the proposed method it is dependent on the propagation direction due to anisotropy.

Proposed method can estimate not only magnitude but also direction of principal stress by measuring the rate of change of the surface wave velocity in any three directions before and after stress is applied, which does not require any information of stress direction in advance. The requirement to measure the initial surface wave velocities in three directions before stress is applied also will not be a problem at all in the comparison of the stress states before and after the machining such as punching or bending. Furthermore, this method is very simple and practical compared to the previous methods of Delsanto et al. [8] and Thompson et al. [10]. Since there is no need to pre-measure the directional elasticity of the material, or pre-check the direction of anisotropy. In addition, this approach requires measuring only the surface wave velocities in any three directions before and after the stress is applied.

For verification of the proposed method, numerical simulations were performed. For the simulation, we generated the surface wave velocity data with respect to the propagation direction using Delsanto's model for several typical cases when both anisotropy and stress exist. We then analyzed the error that occurs when the stress is estimated by the proposed method using this data. In the simulation, three different combinations of the anisotropy direction, principal stress direction and measurement direction, and two different stress states, uniaxial stress and biaxial stress, were evaluated.

To verify the proposed method experimentally, stresses applied by a tensile testing machine were estimated by measuring the surface wave velocities in three directions. The experiments were carried out for two types of aluminum plates, one nearly isotropic and one with relatively greater anisotropy. Two specimens were prepared for each type where one was used for the measurement of the acoustoelastic coefficients and the other was used for stress estimation. Stress up to 90 MPa was applied and the performance of stress estimation was verified by comparing the estimated stress with the applied stress.

#### **2. Theory**

#### *2.1. Stress Estimation Method in Isotropic Material*

The acoustoelastic effect of surface wave in isotropic materials is expressed as follows,

$$\frac{V\_{\theta} - V\_{R}}{V\_{R}} = \frac{K\_{1} + K\_{2}}{2}(\sigma\_{1} + \sigma\_{2}) + \frac{K\_{1} - K\_{2}}{2}(\sigma\_{1} - \sigma\_{2})\cos 2\theta \tag{1}$$

which represents the relationship between the surface wave velocity change and the principal stresses [9]. Here, *V<sup>R</sup>* is the reference velocity of the surface wave in the unstressed state, *V*<sup>θ</sup> is the surface wave velocity in the stressed state in the θ direction, and θ is the angle between the surface wave propagation direction and the principal stress σ1. σ<sup>1</sup> and σ<sup>2</sup> are principal stresses, and *K*<sup>1</sup> and *K*<sup>2</sup> are the acoustoelastic coefficients.

The acoustoelastic coefficients *K*<sup>1</sup> and *K*<sup>2</sup> represent the linear proportional coefficients between the velocity change rate and stress when the surface wave propagates in the stress direction and in the direction perpendicular to the stress direction under the uniaxial stress condition, respectively. That is, when σ<sup>1</sup> = σ and σ<sup>2</sup> = 0, the coefficient can be expressed as follows:

$$K\_1 = \frac{V\_0 - V\_R}{V\_R} \frac{1}{\sigma} \tag{2}$$

$$K\_2 = \frac{V\_{90} - V\_R}{V\_R} \frac{1}{\sigma} \tag{3}$$

where, *V*<sup>0</sup> is the surface wave velocity in the direction of stress and *V*<sup>90</sup> is the velocity in the direction perpendicular to the stress.

To determine the three terms (σ1, σ2, θ) in Equation (1), three-directional surface wave velocities can be used. As an example, using the velocities in the three-directions with a 45◦ difference, as shown in Figure 1, the following three equations can be obtained.

$$\frac{V\_{\theta} - V\_{R}}{V\_{R}} = \frac{K\_{1} + K\_{2}}{2}(\sigma\_{1} + \sigma\_{2}) + \frac{K\_{1} - K\_{2}}{2}(\sigma\_{1} - \sigma\_{2})\cos 2\theta \tag{4}$$

$$\frac{V\_{\theta-45} - V\_R}{V\_R} = \frac{K\_1 + K\_2}{2}(\sigma\_1 + \sigma\_2) + \frac{K\_1 - K\_2}{2}(\sigma\_1 - \sigma\_2)\cos 2(\theta - 45) \tag{5}$$

$$\frac{V\_{\theta-90} - V\_R}{V\_R} = \frac{K\_1 + K\_2}{2}(\sigma\_1 + \sigma\_2) + \frac{K\_1 - K\_2}{2}(\sigma\_1 - \sigma\_2)\cos 2(\theta - 90) \tag{6}$$

In the above equations, *V*θ, *V*θ−45, and *V*θ−<sup>90</sup> are surface wave velocities propagating in the θ, θ − 45◦ , and θ − 90◦ directions, respectively, in the stressed material, and *V<sup>R</sup>* is the reference velocity in the unstressed state. From Equations (4)–(6), the principal stresses, σ<sup>1</sup> and σ2, and the angle θ can be determined as follows:

$$\sigma\_1 = \frac{1}{2} \left[ \frac{1}{K\_1 + K\_2} \left( \frac{V\_\theta}{V\_R} + \frac{V\_{\theta-90}}{V\_R} - 2 \right) + \frac{1}{K\_1 - K\_2} \left( \frac{V\_\theta}{V\_R} - \frac{V\_{\theta-90}}{V\_R} \right) \cdot \frac{1}{\cos 2\theta} \right] \tag{7}$$

$$
\sigma\_2 = \frac{1}{2} \left[ \frac{1}{K\_1 + K\_2} \left( \frac{V\_\theta}{V\_R} + \frac{V\_{\theta-90}}{V\_R} - 2 \right) - \frac{1}{K\_1 - K\_2} \left( \frac{V\_\theta}{V\_R} - \frac{V\_{\theta-90}}{V\_R} \right) \cdot \frac{1}{\cos 2\theta} \right] \tag{8}
$$

1

$$\theta = \frac{1}{2} \left[ \tan^{-1} \left( \frac{\frac{V\_{\theta}}{V\_{R}} + \frac{V\_{\theta - 90}}{V\_{R}} - 2 \frac{V\_{\theta - 45}}{V\_{R}}}{\frac{V\_{\theta}}{V\_{R}} - \frac{V\_{\theta - 90}}{V\_{R}}} \right) \right] \tag{9}$$

*θ*

Using these three equations, when the acoustoelastic coefficients (*K*1, *K*2) and the reference surface wave velocity (*VR*) in the unstressed condition are known, the principal stresses (σ1, σ2) and the angle (θ) can be estimated using the three-directional measurements of the surface wave velocity. *σ σ θ*

0

1

1 2 12

2 2


1 2 12

() ( - )cos2

*θ θ σ σ σ*

*θ* **Figure 1.** Schematic of the relationships of the angles between the principal stress direction and three-measuring directions in 45◦ increments, in which θ is the angle from a measuring direction to the principal stress (σ<sup>1</sup> ) direction, and *V*θ, *V*θ−45, *V*θ−<sup>90</sup> are surface wave velocities propagating in the θ, θ − 45◦ , and θ − 90◦ directions, respectively.

#### *2.2. Stress Estimation Method in Weakly Anisotropic Material*

The method mentioned above is based on the theory for an ideal isotropic material, in which the surface wave velocity in the unstressed condition is constant regardless of the direction of propagation. However, the actual material is not always produced in a perfectly isotropic condition and has anisotropic properties. In this case the surface wave velocity is directionally dependent, and the isotropic method cannot be directly applied. As will be shown later in the numerical simulations, a very large error occurs even in the case of a very small anisotropy.

To compensate for this anisotropic effect, in this study, Equation (1) for the isotropic case was modified as follows,

$$\frac{V\_{\theta} - V\_{\theta, R}}{V\_{\theta, R}} = \frac{K\_1 + K\_2}{2}(\sigma\_1 + \sigma\_2) + \frac{K\_1 - K\_2}{2}(\sigma\_1 - \sigma\_2)\cos 2\theta \tag{10}$$

where *V<sup>R</sup>* in Equation (1) is replaced with *V*θ,*<sup>R</sup>* which is the surface wave velocity in the θ direction in the unstressed state. This equation represents the rate of change of the surface wave velocity before and after applying stress in the θ direction.

Then, the acoustoelastic coefficients *K*<sup>1</sup> and *K*<sup>2</sup> can be rederived as follows under the uniaxial stress condition (σ<sup>1</sup> = σ and σ<sup>2</sup> = 0):

$$K\_1 = \frac{V\_0 - V\_{0,R}}{V\_{0,R}} \frac{1}{\sigma} \tag{11}$$

$$K\_2 = \frac{V\_{90} - V\_{90,R}}{V\_{90,R}} \frac{1}{\sigma} \tag{12}$$

where *V*0,*<sup>R</sup>* and *V*90,*<sup>R</sup>* are the reference surface wave velocities parallel and perpendicular to the stress direction in the unstressed state, respectively. As a result, Equations (11) and (12) replace *V<sup>R</sup>* in Equations (2) and (3) with the initial directional velocities, *V*0,*<sup>R</sup>* and *V*90,*R*, respectively.

Next, for stress estimation using three-directional measurements of the surface wave velocity in the same way as isotropic theory, we derive the following three equations from Equation (10),

$$\sigma\_1 = \frac{1}{2} \left[ \frac{1}{K\_1 + K\_2} \left( \frac{V\_\theta}{V\_{\theta, R}} + \frac{V\_{\theta-90}}{V\_{\theta-90, R}} - 2 \right) + \frac{1}{K\_1 - K\_2} \left( \frac{V\_\theta}{V\_{\theta, R}} - \frac{V\_{\theta-90}}{V\_{\theta-90, R}} \right) \cdot \frac{1}{\cos 2\theta} \right] \tag{13}$$

$$
\sigma\_2 = \frac{1}{2} \left[ \frac{1}{K\_1 + K\_2} \left( \frac{V\_\theta}{V\_{\theta, R}} + \frac{V\_{\theta-90}}{V\_{\theta-90, R}} - 2 \right) - \frac{1}{K\_1 - K\_2} \left( \frac{V\_\theta}{V\_{\theta, R}} - \frac{V\_{\theta-90}}{V\_{\theta-90, R}} \right) \cdot \frac{1}{\cos 2\theta} \right] \tag{14}
$$

$$\theta = \frac{1}{2} \left[ \tan^{-1} \left( \frac{\frac{V\_{\theta}}{V\_{\theta,R}} + \frac{V\_{\theta-90}}{V\_{\theta+90,R}} - 2 \frac{V\_{\theta-45}}{V\_{\theta-45,R}}}{\frac{V\_{\theta}}{V\_{\theta,R}} - \frac{V\_{\theta-90}}{V\_{\theta-90,R}}} \right) \right] \tag{15}$$

where, *V*θ*,R*, *V*θ−45,*R*, and *V*θ−90,*<sup>R</sup>* are the reference surface wave velocities in the θ, θ − 45◦ and θ − 90◦ directions in the unstressed state, respectively, and *V*θ, *V*θ−45, and *V*θ−<sup>90</sup> are the surface wave velocities in the stressed state. As a results, Equations (13)–(15) replace the reference velocity V<sup>R</sup> in the unstressed condition in Equations (7)–(9) with the initial velocity in each measurement direction.

#### **3. Numerical Simulations**

Numerical simulations were performed to verify the stress estimation performance in typical stress states when the proposed method was applied to anisotropic materials. For this, the anisotropy of a material is assumed to be orthogonal, and the acoustic anisotropy rate (η) is defined as the ratio of the difference between surface wave velocity in the direction in which the wave velocity becomes maximum (*Vmax*) and the surface wave velocity in the direction in which the wave velocity becomes minimum (*Vmin*) to the minimum surface wave velocity (*Vmin*) as follows:

$$
\eta = \frac{V\_{\text{max}} - V\_{\text{min}}}{V\_{\text{min}}} \tag{16}
$$

The numerical simulation requires three directional surface wave velocities in unstressed and stressed states. For this, the model suggested by Delsanto et al. [8] and Mase and Delsanto [9] was used. For the simple analysis applied in the study of Mase, only one elastic constant, *C*11, was changed by anisotropy, and thus Equation (17) can be applied

$$\begin{split} \frac{V\_{\theta} - V\_{R}}{V\_{R}} &= \frac{1}{2\mu} A\_{2222} \cos^{4}\phi \mathbf{C}'\_{11} + \quad (a\_{\text{o0}} + \sum\_{k=1}^{3} \frac{a\_{\text{ok}} p\_{k}}{\mu}) \frac{\sigma\_{1} + \sigma\_{2}}{\mu} \\ &+ (a\_{\text{A0}} + \sum\_{k=1}^{3} \frac{a\_{\text{Ak}} p\_{k}}{\mu}) \cos 2\theta \frac{\sigma\_{1} - \sigma\_{2}}{\mu} \end{split} \tag{17}$$

Here, *V<sup>R</sup>* is the surface wave velocity in the unstressed isotropic medium, and *V*<sup>θ</sup> is the wave velocity affected by anisotropy and the acoustoelastic effect in the θ direction. µ, *A*2222, ασ*<sup>k</sup>* , α∆*<sup>k</sup>* (*k*: 0, 1, 2, 3), and ν*<sup>k</sup>* (*k*: 1, 2, 3) are material constants of the isotropic material, and *C'*<sup>11</sup> is the anisotropic variation of *C*11. The left side of Equation (17) corresponds to the rate of change of the surface wave velocity, the first term on the right side reflects the effect of anisotropy on the change of velocity, and the other terms refer to the contribution of the acoustoelastic effect on the change of velocity In this model, φ is the angle between the acoustic anisotropy direction and the surface wave propagation direction, and θ is the angle between the principal stress direction and the surface wave propagation direction, as shown in Figure 2. The acoustic anisotropy direction is defined as the direction in which the velocity of the surface wave is maximum.

*θ*

*θ*

*θ*

*θ*

*σ θ σ ϕ <sup>θ</sup> θ− θ− θ θ*− *θ*− **Figure 2.** Principal stress (σ<sup>1</sup> ) direction, and three-measuring directions in 45◦ increments, in which θ is the angle from a measuring direction to the principal stress (σ<sup>1</sup> ) direction and φ is the angle to the anisotropy direction, and *V*θ, *V*θ−45, *V*θ−<sup>90</sup> are surface wave velocities propagating in the θ, θ − 45◦ , and θ − 90◦ directions, respectively. *σ θ σ ϕ <sup>θ</sup> θ− θ− θ θ*− *θ*− *σ θ σ ϕ <sup>θ</sup> θ− θ− θ θ*− *θ*− *<sup>σ</sup> θσ <sup>ϕ</sup> θ θ*−*θ*−

*θ* − − *θ σ* The three-directional measurement technique was applied using the velocities of the surface waves propagating towards the θ direction and −45◦ and −90◦ to the θ direction. Simulations were performed for the six cases of the combinations of the different angle conditions shown in Table 1 and different stress conditions shown in Table 2. Table 1 shows the three cases for setting the anisotropy direction, principal stress direction and propagation direction. Case 1 is the case where all of the anisotropy, principal stress (σ1), and surface wave propagation directions are the same. In Case 2, the anisotropy direction and the principal stress direction are coincident, but the surface wave propagation direction is different. Case 3 is a general case where all directions are inconsistent. *θ* − − *θ σ θ* − − *θ σ*

*σ σ σ σ σ* Table 2 shows the given stresses for the two groups. Group A is the case where only uniaxial stress (σ<sup>2</sup> = 0) is applied and σ<sup>1</sup> is varied among 30, 60 and 90 MPa. Group B is the case where biaxial stress is applied while fixing σ<sup>1</sup> at 90 MPa and varying σ<sup>2</sup> among three stresses (σ<sup>2</sup> = 30, 60, 90 MPa). *σ σ σ σ σ σ σ σ σ σ σσσ <sup>σ</sup>*


**Table 1.** Directions of the anisotropy, principal stress, and the surface wave propagation for numerical simulation.

**Table 2.** Two groups of applied stress setting for numerical simulation.


The isotropic elastic constants of aluminum in the COMSOL database, which are listed in Table 3, were used to calculate the surface wave velocity, where λ and µ are Lame constants and *l*, *m* and *n* are Murnaghan constants.


**Table 3.** Material constant used to numerical verification.

The anisotropy rate η changed from 0 to 0.02 in 15 steps by varying the value of *C*11. That is, the anisotropy was varied up to 2%. Figure 3 shows an example of the calculation results for η = 0.02, and principal stress direction is coincidence with the anisotropy direction and principal stress magnitudes are σ<sup>1</sup> = 90 MPa and σ<sup>2</sup> = 0. The velocity distribution of the unstressed state has maximum values in the directions of 0◦ and 180◦ , and a minimum value in the directions perpendicular to them. Also, it is possible to verify that the velocity distribution is slightly changed by the acoustoelastic effect when stress is applied. *λ λ μ*

To demonstrate the difficulty of applying the isotropic theory to the anisotropic material first, we applied the surface wave velocity data set obtained for the Case 1A to the isotropic theory of Equations (7)–(9), prior to applying the proposed method. Figure 4 shows the stress estimation error of σ1, in which surface wave velocities in three directions of 0◦ , 45◦ , and 90◦ were used. As can be seen from the results, the stress estimation error is ridiculously large even for very weak anisotropy. Therefore, it is very difficult to estimate the stress of anisotropic materials using the isotropic theory. − − − η *η σ σ*

Next, to verify the performance of the proposed method, the surface wave velocity data set for all six conditions shown in Tables 1 and 2 was generated. Stress estimation was performed by using surface wave velocities in three directions of 0◦ , 45◦ , and 90◦ . The errors between the estimated stress and given stress of σ<sup>1</sup> were then calculated, where the results are shown in Figure 5. In all cases tested, the errors were less than 3%. In particular, in Case 1A, where all directions of anisotropy, principal stress (σ1) and surface wave propagation are identical, no errors are generated in the stress estimation regardless of the stress. This is because Delsanto's model in Equation (17) used in this numerical simulation is consistent with our modified Equation (10). Actually, by rearranging Equation (17) for the uniaxial stress and angle conditions of Case 1A, it becomes equal to Equation (10) (see the Appendix A). *σ*

**Figure 3.** Surface wave velocity according to the propagation direction for unstressed and stressed states in an anisotropic material.

1 estimation error [%]

*σ*

**Figure 4.** Stress estimation errors of the isotropic theory with respect to the anisotropy rate (Case 1A).

 **Figure 5.** Stress estimation errors of the proposed method with respect to the anisotropy rate (**a**) Case 1A; (**b**) Case 1B; (**c**) Case 2A; (**d**) Case 2B; (**e**) Case 3A; (**f**) Case 3B.

Meanwhile, in the case of A with uniaxial stress, the stress estimation errors are independent of the given stress, and only increase with the anisotropy rate. In the case of B with biaxial stress, the stress estimation errors increase with both anisotropy and stress. However, it must be noted again that the estimation error in all cases are very small. Based on these results, it was confirmed that the proposed method for stress estimation can properly estimate stress in a weakly anisotropic material.

#### **4. Specimens and Experiments**

#### *4.1. Specimens and Measurement of the Anisotropy Rate*

To verify the proposed method experimentally, stress estimation experiments were performed. Two kinds of Al6061 material with different acoustic anisotropy rates were used as the test specimens. One (SA) was produced by a hot rolling process and the other (SB) was produced by an extrusion process. Two test specimens were prepared for each material as shown in Figure 6. One specimen was used to check the anisotropy and acoustoelastic coefficient, and the other was used to verify the stress estimation method. The size of the hot rolled specimen is 350 mm × 62 mm × 4 mm (length × width × thickness) and the size of the extruded specimen is 350 mm × 60 mm × 5 mm.

**Figure 6.** Aluminum plate specimens used in experiments.

To verify the surface wave velocity distribution and anisotropy rate of each specimen, the velocity was measured in 24 directions with an interval of 15◦ . The experimental setup is shown in Figure 7. PZT transducers with a center frequency of 2.25 MHz and wedges (ABWML-7T-90, Olympus) were used for surface wave excitation and reception. A pulser-receiver (Olympus 5077PR) was used for the excitation with a single pulse and a high-resolution digital oscilloscope (Lecroy HDO4034A) with a 10 GS/s sampling rate was used for precise measurement of small velocity changes. The surface wave propagation distance was fixed at 10 mm using a jig. The time of flight (TOF) was used to measure the change rate of the surface wave velocities or surface wave velocity ratios used in Equations (11)–(15). The TOF was measured based on the arrival time of the peak of the pulse signal. In addition, the delay time irrespective of surface wave propagation on the specimens, such as the delay in the wedge, was compensated by calibration in advance.

**Figure 7.** Surface wave velocity measurement set-up, in which the distance fixing zig holds the wedges to keep the surface wave propagation distance constant.

**–**

The measurement results of the surface wave velocity distribution are shown in Figure 8, in which the velocity was normalized with respect to the maximum velocity. Maximum velocity (*Vmax*) and minimum velocity (*Vmin*) in the hot-rolled specimen are 2967.1 m/s and 2957.5 m/s, respectively, and they are 2951.8 m/s and 2876.2 m/s, respectively, in extruded specimen. Figure 8a display the results of the hot-rolled aluminum specimen (SA.1) and Figure 8b shows the results of the extruded aluminum specimen (SB.1). The acoustic anisotropy rate of the hot-rolled aluminum specimen (SA.1) was 0.005 (0.5%) and that of the extruded aluminum specimen was 0.026 (2.6%). This demonstrates that specimen SA is nearly isotropic while specimen SB has a weak anisotropic feature.

**Figure 8.** Surface wave velocity distributions of the two specimens: (**a**) Hot rolled specimen (SA.1); (**b**) Extruded specimen (SB.1).

#### *4.2. Measurement of Acoustoelastic Coe*ffi*cients*

**–** The acoustoelastic coefficients were obtained prior to the stress estimation experiments. The measurement set-up is shown in Figure 9. Based on the definition of the acoustoelastic coefficient, the change rates of the surface wave velocity were measured with stress in the stress direction (*K*1) and in the perpendicular direction (*K*2) in the uniaxial stress state. A tensile tester (Zwick-Roell) was used to apply the stress. The loading force was varied between 2500 N–22,500 N with an interval of 2500 N. The loading interval of 2500 N corresponds to a stress of about 10 MPa. The velocity measurement set-up is identical to the set-up described in Figure 7.

Figure 10 shows the change of the received pulsed surface wave signal according to the applied stress. Figure 10a corresponds to the *K*<sup>1</sup> measurement and Figure 10b reflects *K*2. As shown in the zoomed-in box, the peak points of the pulse signal shift depend on the stress increment. That is, the surface wave velocity in the stress direction decreases with the stress increment, while it increases in the direction perpendicular to the stress.

Figures 11 and 12 show the change rates of the time of flight (TOF) with respect to the applied stress in each direction for the SA.1 and SB.1 specimens, respectively. The results were linearly fitted and the slopes of the fitted lines are the acoustoelastic coefficients defined in Equations (11) and (12). The obtained acoustoelastic coefficients are listed in Table 4.

**Figure 9.** Set-up for acoustoelastic coefficient measurement. The change in surface wave velocity under tensile stress loading is measured. For the measurement of *K1*, transducers (blue colored) are placed in the loading direction, while for the measurement of *K2*, transducers (orange colored) are placed in the direction perpendicular to the loading direction.

 **Figure 10.** Received signals obtained in the acoustoelastic coefficient measurements: (**a**) Measurement of *K*1—The peak moves to the right because the velocity of surface wave in the direction of tensile stress decreases with increasing stress; (**b**) Measurement of *K*2—The peak moves to the left because the velocity of surface wave in the direction perpendicular to the tensile stress increases with increasing stress.

 **Figure 11.** Change rates of the TOF with respect to the applied stress in specimen SA.1 for measurements of (**a**) *K*<sup>1</sup> and (**b**) *K*<sup>2</sup> .

−

−

− −

− −

θ

 **Figure 12.** Change rates of the TOF with respect to the applied stress in specimen SB.1 for measurements of (**a**) *K*<sup>1</sup> and (**b**) *K*<sup>2</sup> .

**Table 4.** Values of the acoustoelastic coefficients for specimens SA.1 and SB.1.


− −

−

#### *4.3. Experiments for Stress Estimation*

θ The stress estimation set-up is shown in Figure 13. The experimental set-up is similar to the acoustoelastic coefficient measurement except for the velocity measurement direction. The surface wave velocities were measured in three directions with difference of 45◦ , where the angle θ was set to 30◦ . This measurement condition corresponds to the Case 2B in the numerical simulation. Stress was applied by tensile loading from 7500 N to 22,500 N at an increment of 2500 N. Specimens SA.2 and SB.2 were used for the stress estimation.

**Figure 13.** Experimental set-up for the proposed stress estimation.

The stress estimation results of the hot-rolled (SA.2 specimen) and extruded (SB.2 specimen) aluminum plates are plotted in Figure 14 and the detailed values are given in Tables 5 and 6, respectively. For the hot-rolled aluminum specimen, the estimated results match the given stresses and angle well. Similarly, the stress and angle estimations for the extruded specimens also exhibit a good agreement with the applied values even though this specimen has a larger anisotropy rate than the hot-rolled

Estimated Stress [MPa]

−

Estimated stress [MPa]

specimen. The maximum error was 5.14 MPa at a given stress of 90.92 MPa in the hot-rolled specimen, and the other errors were less than 5.00 MPa for both specimens. As a result, it was confirmed that stresses can be estimated in weakly anisotropic materials using the proposed method. Note that the experiments were carried out for the case where the stress in the thickness direction was uniform only to verify the proposed method. If the stress varies with depth, surface waves of different wavelengths can be used to probe the stress field at different depths.

**Figure 14.** Stress estimation results (**a**) Hot-rolled specimen (SA.2) and (**b**) Extruded specimen (SB.2).

**Table 5.** Estimated results and errors for the hot-rolled specimen (SA.2).


**Table 6.** Estimated results and errors for the hot-rolled specimen (SB.2).


#### **5. Conclusions**

This study proposed a method to estimate stress in anisotropic materials, which involved modification of the equation for the isotropic theory. In this method, the ratios of the surface wave velocities in the stressed state to those in the stress-free state in three directions were utilized, enabling simultaneous estimation of the principal stresses and their directions.

Numerical simulations and experiments were conducted for verification. The results of the numerical simulation showed that the stress estimation error occurred at a very small order when the acousitc anisotropy rate was less than 2%. Experiments were performed on hot-rolled and extruded aluminum plates that had different anisotropy rates of, 0.05% and 2.6%, respectively. The stresses applied by a tensile tester were estimated by measuring the ratios of the surface wave velocities in the stressed state to those in the unstressed state in three-directions. The estimated results showed good agreement with the given stresses in the range from 30 to 90 MPa. These results demonstrate that the stress estimation method proposed in this study can effectively estimate the stress in weakly anisotropic materials. The proposed method is simple, practical and is expected to be useful for monitoring changes in surface stress before and after machining such as punching or bending of plates. **Author Contributions:** Conceptualization, J.J. and Y.-D.S.; Methodology, J.J. and Y.-D.S.; Software, Y.-D.S.; Validation, J.J., Y.-D.S. and K.-Y.J.; Formal Analysis, J.J.; Investigation, J.J. and Y.-D.S.; Resources, K.-Y.J.; Data Curation, Y.-D.S.; Writing-Original Draft Preparation, J.J.; Writing-Review and Editing, K.-Y.J.; Visualization, Y.-D.S.; Supervision, K.-Y.J.; Project Administration, K.-Y.J.; Funding Acquisition, K.-Y.J. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was supported by the Nuclear Power Research and Development Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (NRF-2013M2A2A9043241).

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Nomenclature**


#### **Appendix A**

Considering the uniaxial stress (σ<sup>1</sup> = σ, σ<sup>2</sup> = 0) and angle condition of Φ = θ, Equation (17) can be rearranged as follows,

$$\frac{V\_{\partial} - V\_{R}}{V\_{R}} = \frac{1}{2\mu} A\_{2222} \cos^{4}\theta C\_{11}' + (a\_{\circ 0} + \sum\_{k=1}^{3} \frac{a\_{0k} v\_{k}}{\mu}) \frac{\sigma}{\mu} - (a\_{\Delta 0} + \sum\_{k=1}^{3} \frac{a\_{\Delta k} v\_{k}}{\mu}) \cos 2\theta \frac{\sigma}{\mu} \tag{A1}$$

By substituting σ = 0, the velocity in the unstressed state, in the θ direction can be obtained as follows:

$$V\_{\theta,R} = (\frac{1}{2\mu} A\_{2222} \cos^4 \theta C\_{11}' + 1) V\_R \tag{A2}$$

On the other hand, the velocity in the stressed state in the θ direction is represented as follows,

$$V\_{\theta} = \left[\frac{1}{2\mu} A\_{2222} \cos^4 \theta C\_{11}' + (a\_{\theta 0} + \sum\_{k=1}^{3} \frac{a\_{0k} p\_k}{\mu}) \frac{\sigma}{\mu} + (a\_{\Delta 0} + \sum\_{k=1}^{3} \frac{a\_{\Delta k} p\_k}{\mu}) \frac{\sigma}{\mu} \cos 2\theta + 1\right] V\_R \tag{A3}$$

Considering the condition of Case 1A shown in Tables 1 and 2, for θ = 0 ◦ , 45◦ and 90◦ , three-directional surface wave velocities in the unstressed state can be expressed by Equations (A4)–(A6)

$$V\_{0,R} = (\frac{1}{2\mu} A\_{2222} \cos^4 0 \mathbf{C}'\_{11} + 1) V\_R = (\frac{1}{2\mu} A\_{2222} \mathbf{C}'\_{11} + 1) V\_R \tag{A4}$$

$$V\_{45,R} = (\frac{1}{2\mu} A\_{2222} \cos^4 45\mathbf{C}'\_{11} + 1)V\_R \tag{A5}$$

$$V\_{90,R} = (\frac{1}{2\mu} A\_{2222} \cos^4 90 \mathcal{C}'\_{11} + 1) V\_R = V\_R \tag{A6}$$

The three-directional velocities in the stressed state are represented by Equations (A7)–(A9).

$$V\_0 = \left[\frac{1}{2\mu} A\_{2222} \mathcal{C}\_{11}' + (a\_{o0} + \sum\_{k=1}^3 \frac{a\_{0k} \nu\_k}{\mu}) \frac{\sigma}{\mu} + (a\_{\Delta 0} + \sum\_{k=1}^3 \frac{a\_{\Delta k} \nu\_k}{\mu}) \frac{\sigma}{\mu} + 1\right] V\_R \tag{A7}$$

$$V\_{45} = \left[\frac{1}{2\mu} \frac{1}{4} A\_{2222} C\_{11} + (a\_{o0} + \sum\_{k=1}^{3} \frac{a\_{0k} \nu\_k}{\mu}) \frac{\sigma}{\mu} + 1\right] V\_R \tag{A8}$$

$$V\_{90} = \left[ (a\_{o0} + \sum\_{k=1}^{3} \frac{a\_{0k} v\_k}{\mu}) \frac{\sigma}{\mu} + 1 \right] V\_R \tag{A9}$$

The acoustoelastic coefficients defined in Equations (11) and (12) then become the following:

$$K\_{1} = \frac{(a\_{\circ 0} + \sum\_{k=1}^{3} \frac{a\_{k0} v\_{k}}{\mu})\frac{1}{\mu} + (a\_{\Delta 0} + \sum\_{k=1}^{3} \frac{a\_{\Delta k} v\_{k}}{\mu})\frac{1}{\mu}}{(\frac{1}{2\mu} A\_{2222} C\_{11}' + 1)}\tag{A10}$$

$$K\_2 = (a\_{o0} + \sum\_{k=1}^{3} \frac{a\_{0k} v\_k}{\mu}) \frac{1}{\mu} \tag{A11}$$

Finally, Equations (A4)–(A11) are substituted into the stress estimation equations shown in Equations (13)–(15) to obtain σ<sup>1</sup> = σ, σ<sup>2</sup> = 0 and θ = 0, which are identical to the given conditions.

#### **References**


© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
