MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2016/6762076 6762076 Research Article The Elastic Constants Measurement of Metal Alloy by Using Ultrasonic Nondestructive Method at Different Temperature Hu Eryi 1 Wang Wenjin 1 Merodio Jose School of Mechatronic Engineering China University of Mining & Technology 1 Daxue Road Xuzhou Jiangsu 221116 China cumt.edu.cn 2016 612016 2016 08 10 2015 13 12 2015 612016 2016 Copyright © 2016 Eryi Hu and Wenjin Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The ultrasonic nondestructive method is introduced into the elastic constants measurement of metal material. The extraction principle of Poisson’s ratio, elastic modulus, and shear modulus is deduced from the ultrasonic propagating equations with two kinds of vibration model of the elastic medium named ultrasonic longitudinal wave and transverse wave, respectively. The ultrasonic propagating velocity is measured by using the digital correlation technique between the ultrasonic original signal and the echo signal from the bottom surface, and then the elastic constants of the metal material are calculated. The feasibility of the correlation algorithm is verified by a simulation procedure. Finally, in order to obtain the stability of the elastic properties of different metal materials in a variable engineering application environment, the elastic constants of two kinds of metal materials in different temperature environment are measured by the proposed ultrasonic method.

1. Introduction

The elastic constants are significant properties of some new metal alloy materials. Furthermore, the stability of the elastic properties at different temperature will affect the engineering application of the material in a terrible environment condition, such as the aerospace spacecraft design and deep ocean resources development. For instance, the basic vibration frequency of some special construction will change corresponding with the environment temperature. Thus, the metal material with a much more steady elastic characteristic is more popular in the engineering design. On the other hand, sometimes the engineers are interested in the elastic constants of a new metal alloy. Particularly in the reverse engineering area, the material property coefficients must be measured without the destruction of the rare sample components. Thus, the nondestructive evaluation approach is very useful in these material characteristic researches.

The ultrasonic wave method can be introduced into the elastic constants measurement of the metal material. Nowadays, many researchers are using ultrasonic nondestructive testing to evaluate the mechanical characteristics on metallic materials, mainly, due to its accuracy, besides being fast and very reliable. As we know, the ultrasonic detection methods are very popular in the medical inspection of the human body for a long time. Several researchers have focused on Young’s modulus of the bone material or carotid artery measurement by using the ultrasonic technique . And then this method is also used in the elastic modulus and other geometrical coefficients detection of other materials, such as alumina ceramics . The ultrasonic method is also introduced into the civil engineering area to obtain the rock’s material properties . Tavares et al.  proposed a nondestructive determination method of the elastic properties in shape memory alloy, stainless steel, and plain carbon steel by using ultrasonic measurements. From the previous useful research works, it is found that the detection of the propagating velocity of the ultrasonic wave is the key point in the elastic constants measurement of these different materials . The digital correlation and Fourier transform methods can both be used in the ultrasonic wave flight time extraction; however the precision should be increased by some novel algorithms .

In this paper, we focus on the elastic constant measurement of metal alloy materials especially in a temperature changing environment. Poisson’s ratio, elastic modulus, and shear modulus will be measured by the ultrasonic wave propagating velocity detection. And then temperature stability of the elastic properties can be provided for the further engineering design application.

2. Principle of the Ultrasonic Measurement Method

As we know, the ultrasonic wave is a kind of elastic waves propagating in the medium with higher frequency. Thus, when the elastic wave is propagating in the isotropic medium without the influence of the volume stress, the wave propagating equation can be expressed as follows: (1) λ + G θ + G 2 U = ρ 2 U t 2 , where θ is the volume strain and U = u , v , w denotes the deformation of one point in the elastic medium. When the propagating elastic wave is an ultrasonic longitudinal wave, the deformation in the elastic medium can be defined as u = x , t , v = 0 , and w = 0 . It is meant that the wave propagating equation can be simplified as follows: (2) 2 u t 2 = V L 2 2 u x 2 . Thus, the propagating velocity of the ultrasonic longitudinal wave denoted as V L is expressed as (3) V L 2 = λ + 2 G ρ = E 1 - μ 1 + μ 1 - 2 μ ρ .

On the other hand, when the propagating elastic wave is an ultrasonic transverse wave, the deformation in the elastic medium can be defined as u = 0 , v = v x , t , and w = 0 . The wave propagating equation can be simplified similarly: (4) 2 v t 2 = V T 2 2 v x 2 . Following the same process, the propagating velocity of the ultrasonic transverse wave denoted as V T is expressed as (5) V T 2 = G ρ = E 2 1 + μ ρ . Combining (3) and (5), the elastic constants of the detected material can be deduced: (6) μ = 1 - 2 V T / V L 2 2 - 2 V T / V L 2 , E = V L 2 ρ 1 + μ 1 - 2 μ 1 - μ , G = V T 2 ρ , where ρ is the density of the material, V L is the ultrasonic longitudinal wave propagating velocity, V T is the ultrasonic transverse wave propagating velocity, μ is Poisson’s ratio, E is the elastic modulus, and G is the shear modulus. From (6), it is found that the elastic constants can be calculated after getting the elastic wave propagating velocities and the material density in a certain experimental environment temperature.

3. Ultrasonic Experimental Measurement System

The schematic diagram of the ultrasonic measurement experimental system is illustrated in Figure 1. The experimental system is constructed by the ultrasonic pulse-receiver, ultrasonic wave transducer, digital oscilloscope, personal computer, detected specimen, and temperature control box. The ultrasonic pulse-receiver is used to excite the ultrasonic transducer and receives the signals reflected from the bottom surface of the detected metal specimen. The type of the ultrasonic pulse-receiver manufactured by Panametrics-NDT is Model 5072PR, and the bandwidth is 35 MHz which is sufficient for this experimental research project. And then, the digital oscilloscope is used for sampling ultrasonic waveform, and the digitized waveforms are stored in the memory of the personal computer. In this experiment, the InfiniiVision 4000X digital oscilloscope is manufactured by Agilent Technologies with 5 GHz maximum sample frequency, so that the single point sampling time period is 0.2 ns. The temperature control box called ZTH100UK(S) is manufactured by Shanghai Zundar. The experimental temperature is controlled by a computer automatically in the range of −70~150°C, and the temperature controlling precision is about ±0.5°C.

Sketch map of the experimental system.

The single ultrasonic wave transducer reflect model is used to transmit and receive the ultrasonic waves. In this experiment, two kinds of ultrasonic transducers which can transmit ultrasonic longitudinal wave and transverse wave are fixed on the detected metal specimen, respectively. The main frequencies of the ultrasonic transducers are both 5 MHz. Particularly, the ultrasonic transverse wave transducer is one kind of normal incidence transverse wave transducer, and this single element contact transducer introduces transverse waves directly into the test piece. The special transducer incorporates a transverse wave crystal in a contact transducer case. Rather than using the principles of refraction by the longitudinal wave, as with the angle beam transducers, to produce transverse waves in a material, the crystal itself produces the transverse wave. However, because the transverse waves do not propagate in liquids, it is necessary to use a very viscous couplant when making measurements with the normal incidence transverse wave transducer.

4. Digital Correlation Method for Measuring the Ultrasonic Flight Time 4.1. Basic Principle of the Algorithm

As we know, the ultrasonic transducer is excited by the high frequency electrical pulse, and then the ultrasonic wave propagates in the detected medium. The ultrasonic echo is reflected back by the bottom surface of the medium, and it can be detected by the transducer at the same time. It is known that there is an obvious correlation between the ultrasonic original signal and the echo signal. By using the digital signal analysis method, the similarity of the two signal waveforms can be expressed as the correlation coefficient. Thus, the ultrasonic wave flight time along the propagating direction in the medium can be extracted by this correlation method. In detail, the ultrasonic wave signal sampled and stored in the computer is denoted as s ( t ) . The sections of the original signal and the echo signal are separated and denoted as s 1 ( t ) and s 2 ( t ) , respectively, and then the correlation coefficient R 12 can be expressed as follows: (7) R 12 = - + s 1 t s 2 t d t - + s 1 2 t d t - + s 2 2 t d t 1 / 2 . However, the ultrasonic signals recorded by the digital oscilloscope and computer are two series of discrete signals arrays. Thus, the normalized correlation coefficient with the discrete signals is written as (8) R 12 = s 1 i s 2 i - s 1 i s 2 i / n s 1 2 i - s 1 i 2 / n s 2 2 i - s 2 i 2 / n 1 / 2 , i = 1 , , n , where n is the calculation length of the signal array and i is the location in the signal array.

In the time domain of signal s ( t ) , there is time delay between the ultrasonic original signal s 1 ( t ) and the echo signal s 2 ( t ) . Actually, the time delay is just the fight time of the ultrasonic propagating in the thickness direction for a round trip. In order to obtain the flight time, a circle calculation is applied to get the array of the correlations. In the circle calculation, the maximum correlation coefficient between the discrete data array s 1 ( 1 ) ~ s 1 ( n ) of s 1 ( t ) and s 2 ( 1 + m ) ~ s 2 ( n + m ) of s 2 ( t ) should be found to track the time position denoted as m in the signal array. And then if the time sampling frequency is denoted as f , the corresponding flight time T is equal to m / f . Finally, the ultrasonic velocity can be obtained by the following formula: (9) V = 2 h f m , where h is the thickness of the detected specimen.

4.2. Calculation Example by the Simulation Signal

Firstly, the reliability of the proposed algorithm should be proved by simulation, and the simulated ultrasonic signal is a declined Sinusoidal function (10) S t = e - k t sin 2 π f 0 t , where k is the declining coefficient of the Sinusoidal function set as 9 × 1 0 - 6 and f 0 is the prime frequency of the ultrasonic transducer set as 5 × 1 0 6  Hz. The signal sample frequency f is set as 5 × 1 0 9  Hz, and the flight time T between the original and the echo signals is set as 4.0 × 1 0 - 6  s. The waveform of the simulated signal is shown in Figure 2, and the calculation result of the correlation distribution is shown in Figure 3.

The waveform of the simulated signal.

The result of the correlation calculation.

In Figure 3, a significant peak of the correlation coefficients obtained in the circle calculation is found in the time axis, and the flight time corresponding to the correlation peak point is 4.0 × 1 0 - 6  s. The calculated result fits well with the theoretical setting value; thus the proposed algorithm is suitable for the ultrasonic signal process in the elastic constants experimental measurement.

4.3. Calculation Example with the Ultrasonic Signal

Following the previous simulation, now the fight time calculation of the practical ultrasonic signal can be achieved. For instance, the signal shown in Figure 4 is a series of ultrasonic transverse waves with the sampling frequency 5 GHz. It is found that there are high frequency noises in the practical ultrasonic signal. But because the correlation of the noise signal is so small, the effect of this type of noise can be removed in the correlation calculation procedure. By using the proposed method the flight time calculation result is shown in Figure 5. Similarly, there is also a peak value in the correlation coefficient distribution map. Furthermore, the flight time corresponding to this peak is T = 1.0061 × 1 0 - 5  s.

The ultrasonic echo signal.

The correlation coefficient distribution.

5. Experimental Results and Analysis

The experimental specimen is manufactured as a circular column shape, and the cross-sectional area is bigger enough than that of the ultrasonic transducers. In the beginning of the experiment, the couplant is smeared uniformly onto the interface between the specimen and the transducer. The couplant corresponding to the ultrasonic longitudinal wave is one kind of normal lubricating oil, but the couplant corresponding to the transverse wave must be the special viscous media because of the normal incidence transverse transducer’s characteristic. The ultrasonic transducers fixed on the experimental specimen are put into the temperature control box, and the experimental environment temperature is increased smoothly in the range of −20~70°C. The ultrasonic echo signals, which are displayed on the screen of the digital oscilloscope, are manually sampled and stored by the computer with a temperature interval of 10°C. In order to decrease the detection errors, the echo signals are sampled three times in each temperature point. By using the proposed correlation algorithm method, the ultrasonic echo flight time can be extracted, and the average values are used for the elastic constants calculation of the detected metal materials. However, the changing of the volume and density of the experimental specimen is not considered in the experiment.

In this study, the ultrasonic nondestructive method is used to measure two kinds of different material specimens’ elastic constants. One is named 3J53 in Chinese standard of metal materials, and the other is a new kind of metal alloy which is manufactured by ourselves. And the elastic constants stability of the metal alloys will be researched in different environmental temperature. The parameters of the experimental specimen, such as diameter, thickness, and the density, are shown in Table 1.

The parameters of the experimental specimen.

Specimen material Diameter Thickness Density
(mm) (mm) (kg/m3)
3J53 metal alloy 24.02 19.35 8000.0
New metal alloy 24.04 20.28 8166.1

After the proposed experimental measurement and data analysis, the calculated Poisson ratio, elastic modulus, and shear modulus of the 3J53 metal alloy material are shown in Figures 6, 7, and 8, respectively. And the calculated results of the new metal alloy material are shown in Figures 9, 10, and 11, respectively.

Poisson’s ratio of 3J53 alloy corresponding to different temperature.

The elastic modulus of 3J53 alloy corresponding to different temperature.

The shear modulus of 3J53 alloy corresponding to different temperature.

Poisson’s ratio of the new metal alloy corresponding to different temperature.

The elastic modulus of the new metal alloy corresponding to different temperature.

The shear modulus of the new metal alloy corresponding to different temperature.

When the metal materials are used in the variable temperature environment, the presence of metal inside the thermal activation leads to dislocation movement and would inevitably lead to dislocation motion and local small plastic deformation. Thus the metal material intensity and elastic parameters will be variable. Furthermore, the heat-activated fan-shaped deformation rate caused by the temperature can affect the propagation speed of the ultrasound in the metal materials. However, because it is different to describe the behavior of this kind of temperature influence with a simple formula, the complex environmental temperature effect can be studied by this special ultrasonic nondestructive method. For different metal materials, the changing of the elastic parameters is different because of the material internal crystal structure characteristics and the manufacturing procedure. From the calculation results, it is found that Poisson’s ratio of the new metal is bigger than that of the 3J53 alloy, and Poisson’s ratios of the two different materials are both decreasing with the increase of the experimental temperature. In the opposite side, the elastic modulus and the shear modulus are increasing with the increase of the experimental temperature. But in the new metal, the elastic modulus and the shear modulus first decrease and then start increasing from 5 to 10 degrees.

The two different materials tested in this study will be used to manufacture primary components in a special frequency detection sensor, so the stability of the elastic parameters plays an important role in the sensor design. In order to give a much more clear difference between the elastic constants of the two kinds of materials with the increase of temperature, a series of new coefficients are defined as follows: (11) β μ = Δ μ μ · Δ T , β E = Δ E E · Δ T , β G = Δ G G · Δ T .

The change in slope of Poisson’s ratio, elastic modulus, and shear modulus is shown in Figures 12, 13, and 14. The maximum changing range of the measured elastic parameters of the β μ , β E , and β G is shown in Table 2.

The maximum changing range of the elastic parameters.

Specimen material β μ β E β G
3J53 metal alloy 2.206 × 1 0 - 4 1.283 × 1 0 - 4 1.417 × 1 0 - 4
New metal alloy 4.645 × 1 0 - 4 2.774 × 1 0 - 4 3.431 × 1 0 - 4

The change in slope of Poisson’s ratio.

The change in slope of elastic modulus.

The change in slope of shear modulus.

It indicates that the change of the elastic constants of the 3J53 metal material is smaller than that of the new metal material. In other words, the elastic property of the 3J53 metal material in different environmental temperature is much more steady than that of the new metal material. Thus, the 3J53 metal material is much more appropriate for the engineering application in the temperature changing environment.

6. Conclusion

Through the simulation and experimental analysis, the following conclusions can be drawn. The proposed ultrasonic nondestructive method can be used for the metal material’s elastic constants measurement. Based on the velocity extraction of the ultrasonic longitudinal wave and transverse wave by the digital correlation method, the elastic constants can be calculated conveniently. The feasibility of the correlation algorithm is verified well by the simulation and experimental approach. The application of the ultrasonic nondestructive method for two kinds of metal materials indicates that Poisson’s ratio is decreasing with the increase of the experimental temperature, and the elastic modulus and the shear modulus are increasing with the increase of the experimental temperature. Furthermore, the elastic modulus and the shear modulus of the new metal first decrease and then start increasing from 5 to 10 degrees. However, the change of the elastic constants of the 3J53 metal material is smaller than that of the new metal material. Thus, it is much more appropriate for the temperature changing environment.

Conflict of Interests

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

Acknowledgments

The authors gratefully acknowledge support from National Natural Science Foundation of China under Grant no. 51205396, a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions, and the Young Teacher Overseas Training Program of CUMT.

Rho J. Y. Ashman R. B. Turner H. Young's modulus of trabecular and cortical bone material: ultrasonic and microtensile measurements Journal of Biomechanics 1993 26 2 111 119 10.1016/0021-9290(93)90042-d 2-s2.0-0027551874 Riley W. A. Barnes R. W. Evans G. W. Burke G. L. Ultrasonic measurement of the elastic modulus of the common carotid artery: The atherosclerosis risk in communities (ARIC) study Stroke 1992 23 7 952 956 10.1161/01.str.23.7.952 2-s2.0-0026690858 Töyräs J. Lyyra-Laitinen T. Niinimäki M. Lindgren R. Nieminen M. T. Kiviranta I. Jurvelin J. S. Estimation of the Young's modulus of articular cartilage using an arthroscopic indentation instrument and ultrasonic measurement of tissue thickness Journal of Biomechanics 2001 34 2 251 256 10.1016/s0021-9290(00)00189-5 2-s2.0-0035141953 Asmani M. Kermel C. Leriche A. Ourak M. Influence of porosity on Young's modulus and poisson's ratio in alumina ceramics Journal of the European Ceramic Society 2001 21 8 1081 1086 10.1016/s0955-2219(00)00314-9 2-s2.0-0035427418 Moradian Z. A. Behnia M. Predicting the uniaxial compressive strength and static Young's modulus of intact sedimentary rocks using the ultrasonic test International Journal of Geomechanics 2009 9 1 14 19 10.1061/(asce)1532-3641(2009)9:1(14) 2-s2.0-58549083569 Khandelwal M. Singh T. N. Correlating static properties of coal measures rocks with P-wave velocity International Journal of Coal Geology 2009 79 1-2 55 60 10.1016/j.coal.2009.01.004 2-s2.0-67349242287 Verma D. Kainthola A. Singh R. Singh T. N. Assessment of Geo-mechanical properties of some Gondwana Coal using P-Wave Velocity International Research Journal of Geology and Mining 2012 2 9 2276 6618 Kainthola A. Singh P. K. Verma D. Singh R. Sarkar K. Singh T. N. Prediction of strength parameters of himalayan rocks: a statistical and ANFIS approach Geotechnical and Geological Engineering 2015 33 5 1255 1278 10.1007/s10706-015-9899-z de Albuquerque V. H. C. Melo T. A. D. A. de Oliveira D. F. Gomes R. M. Tavares J. M. R. S. Evaluation of grain refiners influence on the mechanical properties in a CuAlBe shape memory alloy by ultrasonic and mechanical tensile testing Materials & Design 2010 31 7 3275 3281 10.1016/j.matdes.2010.02.010 2-s2.0-77952744982 Freitas V. L. D. A. Albuquerque V. H. C. D. Silva E. D. M. Silva A. A. Tavares J. M. R. S. Nondestructive characterization of microstructures and determination of elastic properties in plain carbon steel using ultrasonic measurements Materials Science and Engineering A 2010 527 16-17 4431 4437 10.1016/j.msea.2010.03.090 2-s2.0-77952883515 de Albuquerque V. H. C. de Macedo Silva E. Pereira Leite J. de Moura E. P. de Araújo Freitas V. L. Tavares J. M. R. S. Spinodal decomposition mechanism study on the duplex stainless steel UNS S31803 using ultrasonic speed measurements Materials & Design 2010 31 4 2147 2150 10.1016/j.matdes.2009.11.010 2-s2.0-73049112709 Freitas V. L. D. A. Normando P. G. de Albuquerque V. H. C. Silva E. D. M. Silva A. A. Tavares J. M. R. S. Nondestructive characterization and evaluation of embrittlement kinetics and elastic constants of duplex stainless steel SAF 2205 for different aging times at 425°C and 475°C Journal of Nondestructive Evaluation 2011 30 3 130 136 10.1007/s10921-011-0100-1 2-s2.0-80053384835 Hu E. He Y. Chen Y. Experimental study on the surface stress measurement with Rayleigh wave detection technique Applied Acoustics 2009 70 2 356 360 10.1016/j.apacoust.2008.03.002 2-s2.0-55949094616 Friman O. Borga M. Lundberg P. Knutsson H. Detection of neural activity in fMRI using maximum correlation modeling NeuroImage 2002 15 2 386 395 10.1006/nimg.2001.0972 2-s2.0-0036335151 Kim T. J. Grill W. Determination of the velocity of ultrasound by short pulse switched sinusoidal excitation and phase-sensitive detection by a computer-controlled pulse-echo system Ultrasonics 1998 36 1–5 233 238 10.1016/s0041-624x(97)00059-0 2-s2.0-0031996051 Tsai D.-M. Lin C.-T. Fast normalized cross correlation for defect detection Pattern Recognition Letters 2003 24 15 2625 2631 10.1016/S0167-8655(03)00106-5 2-s2.0-0041328321 Yue H.-M. Su X.-Y. Liu Y.-Z. Fourier transform profilometry based on composite structured light pattern Optics & Laser Technology 2007 39 6 1170 1175 10.1016/j.optlastec.2006.08.014 2-s2.0-33847354311