US12038413B2 - Stress gradient high-efficiency non-destructive detection system based on frequency domain calculation of broadband swept frequency signals, and detection method thereof - Google Patents
Stress gradient high-efficiency non-destructive detection system based on frequency domain calculation of broadband swept frequency signals, and detection method thereof Download PDFInfo
- Publication number
- US12038413B2 US12038413B2 US17/551,319 US202117551319A US12038413B2 US 12038413 B2 US12038413 B2 US 12038413B2 US 202117551319 A US202117551319 A US 202117551319A US 12038413 B2 US12038413 B2 US 12038413B2
- Authority
- US
- United States
- Prior art keywords
- frequency
- wedge block
- signals
- stress gradient
- measured
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active, expires
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/043—Analysing solids in the interior, e.g. by shear waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/25—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons
- G01L1/255—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons using acoustic waves, or acoustic emission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
- G01N29/075—Analysing solids by measuring propagation velocity or propagation time of acoustic waves by measuring or comparing phase angle
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/30—Arrangements for calibrating or comparing, e.g. with standard objects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/34—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor
- G01N29/348—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with frequency characteristics, e.g. single frequency signals, chirp signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
- G01N29/4436—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a reference signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0231—Composite or layered materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0234—Metals, e.g. steel
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02827—Elastic parameters, strength or force
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0421—Longitudinal waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/048—Transmission, i.e. analysed material between transmitter and receiver
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/057—Angular incidence, parallel to surface propagation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/102—Number of transducers one emitter, one receiver
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/26—Scanned objects
- G01N2291/269—Various geometry objects
- G01N2291/2694—Wings or other aircraft parts
Definitions
- the disclosure belongs to the field of ultrasonic non-destructive detection, and specifically relates to a stress gradient high-efficiency non-destructive detection system based on frequency domain calculation of broadband swept frequency signals, and a detection method thereof.
- a non-destructive detection method can perform stress measurement by characterizing a relationship between the characteristic physical quantities such as acoustic velocity and magnetic field and the stress without damaging the use performance of the measured component.
- An ultrasonic method is widely used due to unique advantages such as simple measuring device and wide application range.
- the prior art only calculates an acoustic time difference between signals under stress and signals under no stress in a time domain.
- an excitation signal thereof contains a large number of frequency components, and different frequency signals propagate at different depths inside a measured object.
- an acoustic time difference calculated by the time domain method cannot accurately reflect the influence of internal stresses of the measured object at different depths on different frequency components of ultrasonic signals.
- a stress obtained based on such acoustic time measuring method is an average stress within a certain depth range, which cannot reflect the stress gradient on the internal depth of the measured object.
- an ultrasonic stress measuring device When the internal stress of the measured object is measured by means of acoustic elasticity, the change of the ultrasonic velocity due to the stress is very small. Generally, a stress of 100 MPa causes about 1% change in an ultrasonic propagation velocity. Such a small change in velocity is difficult to measure. Therefore, actually, an ultrasonic stress measuring device usually characterizes the stress according to the change of the acoustic time in ultrasonic propagation within a fixed distance.
- the disclosure provides a stress gradient high-efficiency non-destructive detection system based on frequency domain calculation of broadband swept frequency signals, and a detection method thereof, and aims to solve the problem of low stress gradient measuring accuracy, and realize the high-efficiency characterization of the stress gradient in the depth direction.
- a stress gradient high-efficiency non-destructive detection system based on frequency domain calculation of broadband swept frequency signals is provided.
- the detection system includes an arbitrary waveform generator 1 , a power amplifier 2 , a transmitting probe 3 , a receiving transducer 4 , a high-bandwidth receiving and amplifying device 5 , a single-channel high-speed data acquisition system 6 and a PC computer 7 .
- the arbitrary waveform generator 1 is connected with the power amplifier 2 , the power amplifier 2 is connected with the transmitting probe 3 , the transmitting probe 3 and the receiving transducer 4 are both disposed on a calibration device, the receiving transducer 4 is connected with the high-bandwidth receiving and amplifying device 5 , the high-bandwidth receiving and amplifying device 5 is connected with the single-channel high-speed data acquisition system 6 , and the single-channel high-speed data acquisition system 6 is connected with the PC computer 7 .
- the calibration device includes a transmitting probe 3 , a receiving transducer 4 , a wedge block I 53 and a wedge block II 55 .
- the wedge block I 53 and the wedge block II 55 are separated at two ends of a measured object 51 , the transmitting probe 3 is disposed on the wedge block I 53 or the wedge block II 55 , the receiving transducer 4 is disposed on the wedge block II 55 or the wedge block I 53 , the transmitting probe 3 and the receiving transducer 4 are disposed oppositely, and the wedge block I 53 or the wedge block II 55 is located on a plane.
- the stress gradient high-efficiency non-destructive detection method includes the following steps:
- a calibration device for calibrating the LCR wave velocity of the object to be measured in the step 1 includes a transmitting probe 3 , a receiving transducer 4 , a wedge block I 53 and a wedge block II 55 .
- the wedge block I 53 and the wedge block II 55 are separated at two ends of a measured object 51 , the transmitting probe 3 is disposed on the wedge block I 53 or the wedge block II 55 , the receiving transducer 4 is disposed on the wedge block II 55 or the wedge block I 53 , the transmitting probe 3 and the receiving transducer 4 are disposed oppositely, and the wedge block I 53 or the wedge block II 55 is located on a plane.
- a propagation velocity of the wedge block in the calibration device in the step 1 is known, and a propagation velocity of the LCR wave of the measured object 51 is calculated by the refraction law;
- the transmitting probe 3 is excited by a 5-period sine wave modulated by a Hanning window, and at the same time, a data acquisition board is triggered to acquire ultrasonic signals of the receiving transducer 4 ; deflection angles of the receiving and transmitting probes are adjusted, and the energy of received signals is calculated; when the deflection angles of the receiving and transmitting probes meet a first critical refraction angle, the energy of the signal received by the receiving transducer is the highest; and at this time, the propagation velocity of the LCR wave of the measured object meets:
- FFT Fast Fourier Transform
- step 4 is specifically as follows:
- ⁇ 1 - 2 ⁇ 2 ⁇ ( D 2 ⁇ L ) - ⁇ 1 ⁇ ( D 1 ⁇ L ) D 2 ⁇ L - D 1 ⁇ L . ( 4 )
- a panel of ultrasonic probes with arbitrarily adjustable deflection angles are used in the measuring device of the disclosure. According to the condition that the energy of the LCR wave obtained by the receiving transducer is maximum when the deflection angles of the receiving and transmitting probes meet the first critical refraction law, the high-precision characterization of the propagation velocity of the LCR wave is realized.
- the whole measuring device is simple and easy to operate.
- the starting frequency and the cut-off frequency of broadband swept frequency signals are calculated, and corresponding linear swept frequency signals are generated as the excitation signals of the ultrasonic probes to ensure that the wavelength range of the ultrasonic waves excited by the ultrasonic probes completely covers the target detection depth.
- such excitation mode of high-bandwidth probes matched with swept frequency signals reduces the detection cost and improves the convenience of detection.
- the original excitation signals and the ultrasonic signals acquired by the receiving transducer are subjected to FFT to calculate phase delays corresponding to different frequencies, which can be converted to time delays in a time domain according to a formula, and the time delays corresponding to different frequencies correspond to the stresses of different detection depths, so as to calculate the stress gradient in the target depth range of the measured object.
- One advantage of calculating the acoustic time difference by a frequency domain method is that the effect of the Gaussian white noise in the original signals can be greatly weakened, so as to effectively improve the anti-interference ability of acoustic time difference calculation.
- FIG. 1 is a schematic diagram of a calibration device of the disclosure.
- FIG. 2 is a schematic diagram of a measuring system of the disclosure.
- FIG. 3 is a time domain waveform diagram of broadband swept frequency signals of the disclosure.
- FIG. 4 is a frequency domain waveform diagram of broadband swept frequency signals of the disclosure.
- FIG. 5 is an ultrasonic detection model diagram of the stress gradient of the disclosure.
- a stress gradient high-efficiency non-destructive detection system based on frequency domain calculation of broadband swept frequency signals is provided.
- the detection system includes an arbitrary waveform generator 1 , a power amplifier 2 , a transmitting probe 3 , a receiving transducer 4 , a high-bandwidth receiving and amplifying device 5 , a single-channel high-speed data acquisition system 6 and a PC computer 7 .
- the arbitrary waveform generator 1 is connected with the power amplifier 2 , the power amplifier 2 is connected with the transmitting probe 3 , the transmitting probe 3 and the receiving transducer 4 are both disposed on a calibration device, the receiving transducer 4 is connected with the high-bandwidth receiving and amplifying device 5 , the high-bandwidth receiving and amplifying device 5 is connected with the single-channel high-speed data acquisition system 6 , and the single-channel high-speed data acquisition system 6 is connected with the PC computer 7 .
- a positive electrode of the arbitrary waveform generator 1 is connected with a positive electrode of the power amplifier 2
- a negative electrode of the arbitrary waveform generator 1 is connected with a negative electrode of the power amplifier 2
- the positive electrode of the power amplifier 2 is connected with the transmitting probe 3
- the negative electrode of the power amplifier 2 is connected with the transmitting probe 3
- the transmitting probe 3 and the receiving transducer 4 are both disposed on a calibration device
- the receiving transducer 4 is connected with a positive electrode of the high-bandwidth receiving and amplifying device 5
- the receiving transducer 4 is connected with a negative electrode of the high-bandwidth receiving and amplifying device 5
- the positive electrode of the high-bandwidth receiving and amplifying device 5 is connected with a positive electrode of a digital-to-analog conversion chip in the single-channel high-speed data acquisition system 6
- the negative electrode of the high-bandwidth receiving and amplifying device 5 is connected with a negative electrode of the digital-to-analog conversion chip
- a signal generated by the arbitrary waveform generator 1 is a low-voltage swept frequency signal which must be amplified by the power amplifier 2 before the transmitting probe 3 can be excited.
- the LCR wave propagates along an upper surface of the measured object and is received by the receiving transducer 4 .
- the received signal is an mV-order weak signal, which is not easy to be acquired and is easily interfered by random electronic noise, so that it is necessary to perform non-distortion amplification on the original signal and then perform bandwidth filtering on the amplified signal.
- the original signal obtained by the receiving transducer 4 is input into the high-bandwidth receiving and amplifying device 5 , the signal is amplified to an input range of the digital-to-analog conversion chip through coarse gain adjustment and fine gain adjustment, then a lower cut-off frequency and an upper cut-off frequency of a low-pass filter are set according to the bandwidth of the swept frequency signal, the amplified and filtered signal is input into the single-channel high-speed data acquisition system 6 , and is encoded through the FPGA chip, and the sampling signal is transmitted into an upper computer 7 by the PXIE bus and stored for subsequent off-line processing.
- the calibration device includes a transmitting probe 3 , a receiving transducer 4 , a wedge block I 53 and a wedge block II 55 .
- the wedge block I 53 and the wedge block II 55 are separated at two ends of a measured object 51 , the transmitting probe 3 is disposed on the wedge block I 53 or the wedge block II 55 , the receiving transducer 4 is disposed on the wedge block II 55 or the wedge block I 53 , the transmitting probe 3 and the receiving transducer 4 are disposed oppositely, and the wedge block I 53 or the wedge block II 55 is located on a plane.
- the stress gradient high-efficiency non-destructive detection method includes the following steps:
- a calibration device for calibrating the LCR wave velocity of the object to be measured in the step 1 includes a transmitting probe 3 , a receiving transducer 4 , a wedge block I 53 and a wedge block II 55 .
- the wedge block I 53 and the wedge block II 55 are separated at two ends of a measured object 51 , the transmitting probe 3 is disposed on the wedge block I 53 or the wedge block II 55 , the receiving transducer 4 is disposed on the wedge block II 55 or the wedge block I 53 , the transmitting probe 3 and the receiving transducer 4 are disposed oppositely, and the wedge block I 53 or the wedge block II 55 is located on a plane.
- the propagation velocity of the wedge block in the calibration device in the step 1 is known, and the propagation velocity of the LCR wave of the measured object 51 is calculated by the refraction law; the transmitting probe 3 is excited by a 5-period sine wave modulated by a Hanning window, and at the same time, a data acquisition board is triggered to acquire ultrasonic signals of the receiving transducer 4 .
- the deflection angles of the receiving and transmitting probes are constantly adjusted, and the energy of received signals is calculated.
- the deflection angles of the receiving and transmitting probes need to be synchronously adjusted with high precision, and the adjustment accuracy of the angle has a certain influence on velocity calibration.
- the deflection angles of the receiving and transmitting probes meet a first critical refraction angle, the energy of the signal received by the receiving transducer is the highest.
- the propagation velocity of the LCR wave of the measured object meets:
- step 4 is specifically as follows:
- ⁇ 1 - 2 ⁇ 2 ⁇ ( D 2 ⁇ L ) - ⁇ 1 ⁇ ( D 1 ⁇ L ) D 2 ⁇ L - D 1 ⁇ L . ( 4 )
- the starting frequency and the cut-off frequency of swept frequency signals are set, excitation signals of the ultrasonic transducer and received pulse signals are respectively subjected to frequency domain transformation, the phase delays of different frequency components are calculated, and then, the acoustoelastic effect of the internal stress of the measured object on different frequency components in the swept frequency signals can be determined according to the corresponding relationship between the phase change and the acoustic time difference. Because the signals of different frequency components correspond to different measurement depths, the layer-by-layer calculation of the internal stress of the measured object can be realized through the frequency domain calculation of the broadband swept frequency signals, and finally, the high-efficiency characterization of the stress gradient in the depth direction can be obtained.
- a measured object is an aluminum plate which is an isotropic material
- a propagation velocity of an LCR wave inside the aluminum plate is 6300 m/s
- a target depth measuring range is 1 mm to 10 mm
- a starting frequency and a cut-off frequency of broadband swept frequency signals can be 0.6 MHz and 6 MHz as calculated.
- a broadband swept frequency signal waveform is generated by EasyWave software, and then, the data is imported into a Dingyang signal generator in a USB communication mode.
- Low-frequency signals generated by the signal generator are subjected to first-stage amplification by an Antai power amplifier, the ultrasonic probe is excited by high-voltage swept frequency signals, trigger signals are generated at the same time, and ultrasonic signals obtained by the receiving transducer are acquired by a high-speed data acquisition board.
- the signals are subjected to weak signal processing and band-pass filtering first, the signals are amplified to an input voltage range of the data acquisition board, and low-frequency noise and high-frequency noise in the signals are filtered out at the same time.
- the received signals are acquired continuously 10 times, and the acquired signals are subjected to smooth filtering to filter out a part of electronic noise to improve a signal-to-noise ratio of the received signals.
- a zero-stress measurement sample In order to obtain an absolute stress value inside the measured object, a zero-stress measurement sample must be measured first. When the zero-stress sound of different frequency signals is obtained, the actually measured sample is measured. When the actual sound of different frequency signals is obtained, stress values at different depths inside the measured object are calculated according to an acoustoelastic coefficient of the measured object and an acoustic time difference between the actually measured sample and the zero-stress sample. Finally, the stress gradient characterization in the target depth range is obtained according to the superposition relationship of the stresses.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Toxicology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
Abstract
Description
-
- step 1: calibrating an LCR wave velocity of an object to be measured;
- step 2: calculating a starting frequency and a cut-off frequency of broadband swept frequency signals based on the LCR wave velocity of the object to be measured in the
step 1 and a stress gradient measuring range in a depth direction of the object to be measured; - step 3: converting phase delay to time delay information based on the phase delay of the starting frequency and the cut-off frequency in the
step 2; and - step 4: determining stresses of the depths corresponding to different frequency components based on the time delay information in the
step 3 to finally realize layer-by-layer scanning of the stresses at different depths of the measured object.
-
- where c0 represents the propagation velocity of the wedge block, and ϑ0 represents the deflection angle of the wedge block.
S=exp(j2πBt 2 /T) (2);
-
- and the generated broadband swept frequency signal data is loaded into the arbitrary waveform generator so as to generate broadband swept frequency signals for exciting the transmitting probe.
Δφ=wΔt (3),
-
- where w represents an angular frequency of signals, and Δt represents delay time of the signal corresponding to the angular frequency.
-
- during the FFT, original signals are expanded, that is, the number of points of the FFT is increased so as to improve a frequency resolution after the FFT and correspondingly improve a stress gradient spatial resolution in the depth direction, the delay times of different frequency components in broadband swept frequency signals are calculated to obtain stresses of different internal depth layers corresponding to the measured object, and finally, the stress gradient measurement in the depth direction of the measured object is realized by means of the superposition relationship of the stresses; and
- if a measurement depth corresponding to a frequency f1 is D1 and a measured stress value thereof is σ1; and a measurement depth corresponding to a frequency f2 is D2 and a measured stress value thereof is σ2, a stress value σ1-2 of the D1-2 layer is:
-
- in
step 1, an LCR wave velocity of an object to be measured is calibrated; because a propagation velocity of ultrasonic waves and an excitation frequency jointly determine the wavelength of the ultrasonic waves and the wavelength determines a stress measurement depth of the LCR wave, it is very important to perform high-precision characterization of the propagation velocity of the LCR wave of the object to be measured; - in
step 2, a starting frequency and a cut-off frequency of broadband swept frequency signals are calculated based on the LCR wave velocity of the object to be measured in thestep 1 and a stress gradient measuring range in a depth direction of the object to be measured; a mathematical tool is used to generate corresponding waveform data which is to be loaded into a waveform generator to generate a corresponding waveform; an excitation signal generated by the waveform generator is amplified by a power amplifier and then connected to a transmitting transducer, a signal of a receiving transducer at the other end is amplified by a first-order weak signal, then received by a high-speed acquisition board and transmitted to an upper computer through a PXIE bus, and the signal is stored for off-line processing; the excitation signal and the received signal are subjected to FFT to calculate different frequency components; - in
step 3, phase delay is converted to time delay information based on the phase delay of the starting frequency and the cut-off frequency in thestep 2; and - in
step 4, stresses of the depths corresponding to different frequency components are determined based on the time delay information in thestep 3 to finally realize layer-by-layer scanning of the stresses at different depths of the measured object.
- in
-
- where c0 represents the propagation velocity of the wedge block, and ϑ0 represents the deflection angle of the wedge block.
S=exp(j2πBt 2 /T) (2);
-
- and the generated broadband swept frequency signal data is loaded into the arbitrary waveform generator so as to generate broadband swept frequency signals for exciting the transmitting probe.
Δφ=wΔt (3),
-
- where w represents an angular frequency of signals, and Δt represents delay time of the signal corresponding to the angular frequency.
-
- during the FFT, original signals are expanded, that is, the number of points of the FFT is increased so as to improve a frequency resolution after the FFT and correspondingly improve a stress gradient spatial resolution in the depth direction, the delay times of different frequency components in broadband swept frequency signals are calculated to obtain stresses of different internal depth layers corresponding to the measured object, and finally, the stress gradient measurement in the depth direction of the measured object is realized by means of the superposition relationship of the stresses; and
- if a measurement depth corresponding to a frequency f1 is D1 and a measured stress value thereof is σ1; and a measurement depth corresponding to a frequency f2 is D2 and a measured stress value thereof is σ2, a stress value σ1-2 of the D1-2 layer is:
Claims (8)
S=exp(j2πBt 2 /T) (2);
Δφ=wΔt (3),
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111062917.0A CN113758617B (en) | 2021-09-10 | 2021-09-10 | Efficient nondestructive testing system and method for stress gradient based on broadband frequency-sweep signal frequency domain calculation |
| CN2021110629170 | 2021-09-10 | ||
| CN202111062917.0 | 2021-09-10 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| US20230081998A1 US20230081998A1 (en) | 2023-03-16 |
| US12038413B2 true US12038413B2 (en) | 2024-07-16 |
Family
ID=78794787
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US17/551,319 Active 2043-03-30 US12038413B2 (en) | 2021-09-10 | 2021-12-15 | Stress gradient high-efficiency non-destructive detection system based on frequency domain calculation of broadband swept frequency signals, and detection method thereof |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US12038413B2 (en) |
| CN (1) | CN113758617B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114459648B (en) | 2022-01-19 | 2022-08-02 | 哈尔滨工业大学 | Multi-mode Lamb wave data fusion-based baseline-free stress online monitoring method |
| CN116539734B (en) * | 2023-04-06 | 2024-10-22 | 广东省智能科学与技术研究院 | Sound velocity measurement method, device and readable storage medium based on ultrasonic resonance method |
| CN120232506A (en) * | 2025-05-30 | 2025-07-01 | 零声科技(苏州)有限公司 | Ultrasonic sound velocity calibration method, thickness gauge, terminal equipment and storable medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102721747A (en) * | 2012-06-08 | 2012-10-10 | 北京航空航天大学 | Non-colinear non-linear ultrasonic nondestructive testing method |
| CN105203638A (en) * | 2015-09-18 | 2015-12-30 | 哈尔滨工业大学深圳研究生院 | Method for detecting distribution of steel member absolute stress along depth on basis of Lcr wave method |
| CN107490446A (en) * | 2017-07-04 | 2017-12-19 | 北京理工大学 | High ferro wheel tread stress supersonic damage-free detection method |
| CN112903820A (en) * | 2021-01-16 | 2021-06-04 | 大连理工大学 | Critical refraction longitudinal wave multi-material detection system based on single-angle wedge block and sound velocity measurement method thereof |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103235039A (en) * | 2013-04-25 | 2013-08-07 | 哈尔滨工业大学 | Ultrasonic TOFD (Time of Flight Diffraction) detection method for weld defects based on linear frequency modulation technology |
| CN110108402B (en) * | 2019-05-20 | 2021-08-06 | 北京工业大学 | A Nonlinear Lamb Wave Mixing Method for Stress Distribution Measurement in Thin Metal Plates |
| CN110243945B (en) * | 2019-07-23 | 2021-04-30 | 大连理工大学 | Ultrasonic TOFD Blind Spot Suppression Method Based on Synthetic Aperture Focusing and Mode-Converted Waves |
| CN112880895B (en) * | 2019-11-29 | 2022-09-20 | 哈尔滨工业大学 | Nonlinear ultrasonic wave-based large-scale high-speed rotation equipment blade residual stress measurement method |
| CN112763580B (en) * | 2020-12-22 | 2022-07-19 | 哈尔滨工业大学(深圳) | Steel member plane stress detection method based on ultrasonic transverse and longitudinal wave combination |
| CN113325075A (en) * | 2021-05-27 | 2021-08-31 | 浙江工业大学 | Nonlinear wave detection method for high-cycle fatigue damage of metal sheet |
-
2021
- 2021-09-10 CN CN202111062917.0A patent/CN113758617B/en active Active
- 2021-12-15 US US17/551,319 patent/US12038413B2/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102721747A (en) * | 2012-06-08 | 2012-10-10 | 北京航空航天大学 | Non-colinear non-linear ultrasonic nondestructive testing method |
| CN105203638A (en) * | 2015-09-18 | 2015-12-30 | 哈尔滨工业大学深圳研究生院 | Method for detecting distribution of steel member absolute stress along depth on basis of Lcr wave method |
| CN107490446A (en) * | 2017-07-04 | 2017-12-19 | 北京理工大学 | High ferro wheel tread stress supersonic damage-free detection method |
| CN107490446B (en) * | 2017-07-04 | 2021-11-02 | 北京理工大学 | Ultrasonic non-destructive testing method for tread stress of high-speed rail wheels |
| CN112903820A (en) * | 2021-01-16 | 2021-06-04 | 大连理工大学 | Critical refraction longitudinal wave multi-material detection system based on single-angle wedge block and sound velocity measurement method thereof |
Non-Patent Citations (4)
| Title |
|---|
| CN-102721747-A_trnaslated (Year: 2012). * |
| CN-105203638-A_translated (Year: 2015). * |
| CN-107490446-B_translated (Year: 2021). * |
| CN-112903820-A_translated (Year: 2021). * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20230081998A1 (en) | 2023-03-16 |
| CN113758617B (en) | 2022-10-18 |
| CN113758617A (en) | 2021-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114487109B (en) | Non-baseline data stress online monitoring method, system, equipment and medium based on single-mode multi-frequency signal fusion | |
| US12130259B2 (en) | Non-baseline on-line stress monitoring system and monitoring method based on multi-mode lamb wave data fusion | |
| CN109764986B (en) | A method for detecting plane stress of steel members based on ultrasonic shear wave phase spectrum | |
| US12203818B2 (en) | Method, system, device and medium for online monitoring of plane stress field without baseline data based on piezoelectric transducer array | |
| US12038413B2 (en) | Stress gradient high-efficiency non-destructive detection system based on frequency domain calculation of broadband swept frequency signals, and detection method thereof | |
| CN111208206B (en) | Composite material ultrasonic guided wave structure health monitoring method in temperature change environment | |
| CN106989703A (en) | A kind of gamut lubrication film thickness supersonic detection device and method | |
| CN105651215A (en) | A Coating Thickness Measuring Method Under Unknown Ultrasonic Velocity | |
| CN102042844A (en) | Sound surface wave measuring sensor and parameter analytical method | |
| CN103344709A (en) | Piezoelectric ultrasonic phased array multi-damage imaging method based on linear mapping frequency dispersion compensation | |
| CN104764803A (en) | Material stress detection technology based on changes of wavelength of supersonic waves | |
| Goujon et al. | Behaviour of acoustic emission sensors using broadband calibration techniques | |
| CN113739967B (en) | Normal stress and shear stress detection method based on acoustic elastic effect | |
| CN103075981B (en) | A kind of ultrasonic thickness test method | |
| CN115219597A (en) | An ultrasonic shear wave detection system for the stress state of thin-walled concrete structures | |
| CN113588414A (en) | Bolt axial stress detection method based on ultrasonic frequency spectrum energy attenuation | |
| US4823609A (en) | Ultrasonic method and apparatus for determining crack opening load | |
| JP4617540B2 (en) | Ultrasonic characteristic measuring method, acoustic anisotropy measuring method, and acoustic anisotropy measuring apparatus | |
| CN115508449B (en) | Defect positioning imaging method based on ultrasonic guided wave multi-frequency sparseness and application thereof | |
| JP3156012B2 (en) | Concrete structure thickness measurement method | |
| CN114562965A (en) | A Quantitative Evaluation Method of Corrosion Damage Based on Ultrasonic Guided Wave Dual-frequency Test | |
| CN110455930B (en) | A fast method for measuring the absolute wavenumber curve of ultrasonic guided waves | |
| CN119619997B (en) | Method and system for positioning damage of carbon fiber composite material plate | |
| CN112146744A (en) | Method for accurately measuring reciprocal constant of reverberation pool | |
| CN118067285A (en) | A surface wave virtual superposition interference amplitude spectrum stress measurement method based on dual-transmitter-dual-receiver sensor array |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| AS | Assignment |
Owner name: HARBIN INSTITUTE OF TECHNOLOGY, CHINA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ZHAO, BO;SHI, WEIJIA;LI, JIAXIN;AND OTHERS;REEL/FRAME:058393/0352 Effective date: 20211130 |
|
| FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
| FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: SMAL); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT RECEIVED |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
| STCF | Information on status: patent grant |
Free format text: PATENTED CASE |