Air-coupled Ultrasonic Detection Method and Device Based on Defect Probability Reconstruction Algorithm

20230061816 · 2023-03-02

    Inventors

    Cpc classification

    International classification

    Abstract

    The disclosure discloses an air-coupled ultrasonic detection method and device based on a defect probability reconstruction algorithm. The method includes the following steps: determining the excitation frequency of a transmitting air-coupled transducer according to a frequency dispersion curve of guided waves and the thickness of a to-be-detected piece; determining the group velocity of an antisymmetric mode according to the excitation frequency, and determining the inclination angle of the transmitting/receiving air-coupled transducer according to the Snell law; obtaining an initial waveform of a defect-free test piece as reference data by adopting a same-side penetration method, then rotating the transmitting/receiving transducer by 360 degrees by taking the Z direction as an axis at preset angle intervals by adopting a rotary scanning method, collecting N groups of signal data of the to-be-detected piece again, comparing the N groups of signal data with the reference data to determine whether the signal characteristics have great changes or not, calculating the defect distribution probability on the to-be-detected piece, and carrying out defect imaging on a rotating coverage area of the transmitting/receiving air-coupled transducer according to the defect distribution probability. According to the method, the precision of traditional air-coupled ultrasonic X and Y scanning detection is improved, and compared with a complex imaging technology, the air-coupled ultrasonic detection method consumes less time.

    Claims

    1. An air-coupled ultrasonic detection method based on a defect probability reconstruction algorithm, comprising the following steps: step S1, determining excitation frequency of a transmitting air-coupled transducer according to a frequency dispersion curve of guided waves and a thickness of a to-be-detected piece; step S2, determining group velocity of an antisymmetric mode according to the excitation frequency, and determining an inclination angle of the transmitting air-coupled transducer and a receiving air-coupled transducer according to the Snell law; step S3, placing a transmitting air-coupled transducer and a receiving air-coupled transducer on a same-side upper surface of a defect-free test piece according to the inclination angle, and obtaining an initial waveform of the defect-free test piece as reference data by adopting a same-side penetration method; step S4, placing the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece according to the inclination angle, and obtaining N groups of signal data of the to-be-detected piece by adopting a rotary scanning method; and step S5, comparing the reference data with the N groups of signal data, solving N signal change correlation coefficients, processing the N signal change correlation coefficients based on the defect probability reconstruction algorithm to obtain defect distribution probability on the to-be-detected piece, and carrying out defect imaging on a rotating coverage area of the transmitting air-coupled transducer and the receiving air-coupled transducer according to the defect distribution probability.

    2. The air-coupled ultrasonic detection method based on the defect probability reconstruction algorithm according to claim 1, wherein the step S4 comprises: step S401, placing the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece according to the inclination angle; and step S402, rotating the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece by adopting the rotary scanning method at present angle intervals to collect new signals, and obtaining the N groups of signal data when the transmitting air-coupled transducer and the receiving air-coupled transducer rotate by 360 degrees.

    3. The air-coupled ultrasonic detection method based on the defect probability reconstruction algorithm according to claim 2, wherein the signal change correlation coefficient is as follows: ρ = C XY σ X σ Y wherein ρ is the signal change correlation coefficient, C.sub.XY is covariance of X and Y, X is a reference data set, Y is signal data after a period of service time, and σ.sub.X and σ.sub.Y are standard deviations of X and Y.

    4. The air-coupled ultrasonic detection method based on the defect probability reconstruction algorithm according to claim 1, wherein the defect distribution probability is the sum of all signal change effects of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer.

    5. The air-coupled ultrasonic detection method based on the defect probability reconstruction algorithm according to claim 4, wherein the defect distribution probability is as follows: P ( x , y ) = .Math. i = 1 N P i ( x , y ) = .Math. i = 1 N A i [ β - R i ( x , y ) β - 1 ] wherein P.sub.i(x,y) is defect distribution probability estimation of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in a multi-order symmetric mode S.sub.i, A.sub.i(x,y)=1−ρ.sub.i is a signal difference coefficient of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in the multi-order symmetric mode S.sub.i, (β−R.sub.i(x,y))/(β−1) is a non-negative space distribution function of the multi-order symmetric mode S.sub.i, and the outline thereof is a group of ellipses.

    6. An air-coupled ultrasonic detection device based on a defect probability reconstruction algorithm, comprising: an excitation frequency determination module, configured to determine the excitation frequency of a transmitting air-coupled transducer according to a frequency dispersion curve of guided waves and the thickness of a to-be-detected piece; an inclination angle determination module, configured to determine the group velocity of an antisymmetric mode according to the excitation frequency, and determine the inclination angle of the transmitting air-coupled transducer and a receiving air-coupled transducer according to the Snell law; a reference data acquisition module, configured to place a transmitting air-coupled transducer and a receiving air-coupled transducer on the same-side upper surface of a defect-free test piece according to the inclination angle, and obtain an initial waveform of the defect-free test piece as reference data by adopting a same-side penetration method; a signal data acquisition module, configured to compare the reference data with N groups of signal data, solve N signal change correlation coefficients, place the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece according to the inclination angle, and obtain the N groups of signal data of the to-be-detected piece by adopting the rotary scanning method; and a defect imaging module, configured to process the N signal change correlation coefficients based on the defect probability reconstruction algorithm to obtain the defect distribution probability on the to-be-detected piece, and carry out defect imaging on a rotating coverage area of the transmitting air-coupled transducer and the receiving air-coupled transducer according to the defect distribution probability.

    7. The air-coupled ultrasonic detection device based on the defect probability reconstruction algorithm according to claim 6, wherein the signal data acquisition module further comprises: a placement unit, configured to place the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece according to the inclination angle; and an acquisition unit, configured to rotate the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece by adopting the rotary scanning method at preset angle intervals to collect new signals, and obtain the N groups of signal data when the transmitting air-coupled transducer and the receiving air-coupled transducer rotate by 360 degrees.

    8. The air-coupled ultrasonic detection device based on the defect probability reconstruction algorithm according to claim 7, wherein the signal change correlation coefficient is as follows: ρ = C XY σ X σ Y wherein ρ is the signal change correlation coefficient, C.sub.XY is a covariance of X and Y, X is a reference data set, Y is signal data after a period of service time, and σ.sub.X and σ.sub.Y are standard deviations of X and Y.

    9. The air-coupled ultrasonic detection device based on the defect probability reconstruction algorithm according to claim 6, wherein the defect distribution probability is the sum of all signal change effects of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer.

    10. The air-coupled ultrasonic detection device based on the defect probability reconstruction algorithm according to claim 9, wherein the defect distribution probability is as follows: P ( x , y ) = .Math. i = 1 N P i ( x , y ) = .Math. i = 1 N A i [ β - R i ( x , y ) β - 1 ] wherein P.sub.i(x,y) is defect distribution probability estimation of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in a multi-order symmetric mode S.sub.i, A.sub.i(x,y)=1−ρ.sub.i is a signal difference coefficient of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in the multi-order symmetric mode S.sub.i, (β−R.sub.i(x,y))/(β−1) is a non-negative space distribution function of the multi-order symmetric mode S.sub.i, and the outline thereof is a group of ellipses.

    Description

    BRIEF DESCRIPTION OF FIGURES

    [0028] The above and/or additional aspects and advantages of the disclosure will become obvious and easy to understand from the description of the embodiments in conjunction with the following drawings, wherein:

    [0029] FIG. 1 is a flowchart of an air-coupled ultrasonic detection method based on a defect probability reconstruction algorithm according to an embodiment of the disclosure;

    [0030] FIG. 2 is a graph of variation of a guided wave dispersion characteristic curve according to an embodiment of the disclosure;

    [0031] FIG. 3 is a schematic diagram of an air-coupled ultrasonic same-side penetration method for detection according to an embodiment of the disclosure;

    [0032] FIG. 4 is a schematic diagram of an air-coupled ultrasonic rotary scanning method for detection according to an embodiment of the disclosure;

    [0033] FIG. 5 is a schematic diagram of an ellipse distribution function of a defect probability reconstruction algorithm according to an embodiment of the disclosure; and

    [0034] FIG. 6 is a schematic structural diagram of an air-coupled ultrasonic detection device based on a defect probability reconstruction algorithm according to an embodiment of the disclosure.

    [0035] Description of reference numerals: 1-transmitting air-coupled transducer, 2-receiving air-coupled transducer, 3-to-be-detected piece, 10-air-coupled ultrasonic detection device based on defect probability reconstruction algorithm, 100-excitation frequency determination module, 200-inclination angle determination module, 300-reference data acquisition module, 400-signal data acquisition module and 500-defect imaging module.

    DETAILED DESCRIPTION

    [0036] Embodiments of the disclosure are described in detail below, and examples of the embodiments are illustrated in the accompanying drawings. The same or similar numbers from the beginning to the end denote the same or similar elements or the elements having the same or similar functions. The embodiments described below with reference to the drawings are illustrative and intended to be explanatory of the disclosure and are not to be construed as limiting the disclosure.

    [0037] An air-coupled ultrasonic detection method and device based on a defect probability reconstruction algorithm, provided by the embodiments of the disclosure, are described below with reference to the accompanying drawings, and the air-coupled ultrasonic detection method based on the defect probability reconstruction algorithm, provided by the embodiment of the disclosure, is described first with reference to the accompanying drawings.

    [0038] FIG. 1 is a flowchart of an air-coupled ultrasonic detection method based on a defect probability reconstruction algorithm according to an embodiment of the disclosure.

    [0039] As shown in FIG. 1, the air-coupled ultrasonic detection method based on the defect probability reconstruction algorithm includes the following steps.

    [0040] In step S1, the excitation frequency of a transmitting air-coupled transducer is determined according to a frequency dispersion curve of guided waves and the thickness of a to-be-detected piece.

    [0041] Specifically, as shown in FIG. 2, the guided waves have a symmetric mode and an antisymmetric mode as well as dispersion characteristics, and a multi-order symmetric mode (S0, S1, . . . , Si) and an antisymmetric mode (A0, A1, . . . , Ai) may be excited at the same excitation frequency. In order to excite a relatively pure mode in the to-be-detected piece by the air-coupled transducers, according to the frequency dispersion curve of the guided waves and the thickness of the to-be-detected piece, it may be known that the excitation frequency of the transmitting air-coupled transducer is smaller than a certain upper limit value, and then the excitation frequency is determined according to the actual performance of the air-coupled transducers.

    [0042] In step S2, the group velocity of the antisymmetric mode is determined according to the excitation frequency, and the inclination angle of the transmitting air-coupled transducer and the receiving air-coupled transducer is determined according to the Snell law.

    [0043] Specifically, those skilled in the art know from research analysis that the in-plane displacement of the symmetric mode is relatively large, and the out-of-plane displacement of the antisymmetric mode is relatively large, so that the antisymmetric mode is adopted to carry out air-coupled ultrasonic detection. After the product of frequency and thickness (the product of the frequency and the thickness of the to-be-detected piece) is determined, the group velocity of the antisymmetric mode A0 is also known, and then the inclination angle of the air-coupled transducers is determined according to the first critical refraction angle of the Snell law and the propagation velocity in air.

    [0044] In step S3, a transmitting air-coupled transducer and a receiving air-coupled transducer are placed on the same-side upper surface of a defect-free test piece according to the inclination angle, and an initial waveform of the defect-free test piece is obtained as reference data by adopting a same-side penetration method.

    [0045] In step S4, the transmitting air-coupled transducer and the receiving air-coupled transducer are placed on the same-side upper surface of the to-be-detected piece according to the inclination angle, and N groups of signal data of the to-be-detected piece are obtained by adopting a rotary scanning method.

    [0046] Further, in an embodiment of the disclosure, step S4 specifically includes: [0047] step S401, the transmitting air-coupled transducer and the receiving air-coupled transducer are placed on the same-side upper surface of the to-be-detected piece according to the inclination angle; and [0048] step S402, the transmitting air-coupled transducer and the receiving air-coupled transducer are rotated on the same-side upper surface of the to-be-detected piece by adopting the rotary scanning method at present angle intervals to collect new signals, and the N groups of signal data are obtained when the transmitting air-coupled transducer and the receiving air-coupled transducer rotate by 360 degrees.

    [0049] Specifically, as shown in FIG. 3 and FIG. 4, firstly, the initial waveform of the defect-free test piece is obtained as the reference data by adopting the same-side penetration method, then the transmitting air-coupled transducer and the receiving air-coupled transducer are placed on the same-side upper surface of the to-be-detected piece according to the inclination angle, and the transmitting air-coupled transducer and the receiving air-coupled transducer are rotated by taking a Z direction as an axis at preset angle intervals by adopting the rotary scanning method to collect new signals again, the rest can be done in the same manner, when the transmitting air-coupled transducer and the receiving air-coupled transducer are rotated by 360 degrees, the N groups of signal data are obtained, then the N groups of signal data are compared with the reference data, and N signal change correlation coefficients are solved to determine whether signal characteristics have great change or not.

    [0050] Signal change may be represented by a signal change correlation coefficient p as follows:

    [00005] ρ = C XY σ X σ Y

    [0051] In the formula, C.sub.XY is covariance of X and Y,

    [00006] C XY = .Math. k = 1 K ( X k - μ x ) ( Y k - μ y )

    [0052] In the formula, μ is an average value of various data sets, K is the length of a data set, X is a reference data set, Y is signal data after a period of service time, σ.sub.X and σ.sub.Y are standard deviations of X and Y.

    [00007] σ X σ Y = .Math. k = 1 K ( X k - μ x ) 2 .Math. k = 1 K ( Y k - μ y ) 2 .

    [0053] In step S5, the reference data are compared with the N groups of signal data, N signal change correlation coefficients are solved, the N signal change correlation coefficients are processed based on the defect probability reconstruction algorithm to obtain defect distribution probability on the to-be-detected piece, and defect imaging is carried out on a rotating coverage area of the transmitting air-coupled transducer and the receiving air-coupled transducer according to the defect distribution probability.

    [0054] That is, according to the embodiment of the disclosure, small defects are detected by measuring the difference between guided wave signals under normal conditions and guided wave signals under fault conditions based on the defect probability reconstruction algorithm.

    [0055] Specifically, as shown in FIG. 5, in order to determine the location of a defect, according to an embodiment of the disclosure, it assumes that the probability of the defect occurring at a certain point may be estimated from the severity of signal change of different sensor pairs and the location of the defect relative to the sensor pairs. This means that the defect will cause the most significant signal change in the direct wave path, and if the defect is far from the direct path of the sensor pairs, the signal change effect will be reduced. The defect distribution probability in a sensor network may be expressed as the sum of all signal change effects for each pair of the transmitting air-coupled transducer and the receiving air-coupled transducer. It assumes here that in simple linearly decreased ellipse distribution, the transmitting air-coupled transducer and the receiving air-coupled transducer are arranged at the focus of an ellipse. It assumes that in a network of N total numbers of pairs of air-coupled transducers, the defect probability estimation at a location (x, y) within the reconstruction region may be written as:

    [00008] P ( x , y ) = .Math. i = 1 N P i ( x , y ) = .Math. i = 1 N A i [ β - R i ( x , y ) β - 1 ]

    [0056] In the formula, P.sub.i(x,y) is defect distribution probability estimation from the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in a multi-order symmetric mode S.sub.i, A.sub.i(x,y) is a signal difference coefficient of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in the multi-order symmetric mode S.sub.i, and A.sub.i(x,y)=1−ρ.sub.i.Math.(β−R.sub.i(x,y))/(β−1) is a non-negative space distribution function of the multi-order symmetric mode S.sub.i, and the outline thereof is a group of ellipses.

    [00009] R i ( x , y ) = { R D i ( x , y ) when R D i ( x , y ) < β β when R D i ( x , y ) β

    [0057] In the formula:

    [00010] R D i ( x , y ) = ( x - x i 1 ) 2 + ( y - y i 1 ) 2 + ( x - x i 2 ) 2 + ( y - y i 2 ) 2 ( x i 2 - x i 1 ) 2 + ( y i 2 - y i 1 ) 2

    is the ratio of the sum of the distance from a point (x, y) to the transmitting air-coupled transducer and the distance from the point (x, y) to the receiving air-coupled transducer in the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in the multi-order symmetric mode S.sub.i to the distance between the transmitting air-coupled transducer and the receiving air-coupled transducer. β is a scale parameter for controlling the size of an effective ellipse distribution region, and β>1. When RD.sub.i, (x,y)=1 means that the point (x,y) is located on a straight line of the pair of the transmitting/receiving air-coupled transducers in the multi-order symmetric mode S.sub.i, at the moment, P.sub.i(x,y)=A.sub.i(x,y); and when RD.sub.i(x,y)=β means that the point (x,y) is located on the boundary of the effective distribution region, at the moment, P.sub.i(x,y)=β. In general, β is selected to be about 1.05. If β is too small, pseudomorphism may be introduced, and if β is too large, resolution may be lost. Usually, if defects occur, a group of signals of the pair of the transmitting/receiving air-coupled transducers will be affected. Consequently, in a defect distribution probability image, compared with other points, the point where the defect is located has obvious larger probability. Therefore, the defect distribution probability is treated by using an image processing technology, for example, a threshold value of a defect estimation image is selected, a defect position may be estimated, and then, defects of the to-be-detected piece are imaged.

    [0058] According to the air-coupled ultrasonic detection method based on the defect probability reconstruction algorithm provided by the embodiment of the disclosure, the accuracy of traditional air-coupled ultrasonic X and Y scanning detection is improved; compared with a complex imaging technology, the air-coupled ultrasonic detection method consumes less time; and under the same environment and measurement conditions, the change is caused by the generation of defects. According to the embodiment of the disclosure, the defect growth may be monitored by tracking increase of the signal difference relative to the normal condition, and the fact health monitoring of the to-be-detected piece is realized.

    [0059] Secondly, the air-coupled ultrasonic detection device based on the defect probability reconstruction algorithm is described according to the embodiment of the disclosure by referring to the accompanying drawings.

    [0060] FIG. 6 is a schematic structural diagram of an air-coupled ultrasonic detection device based on the defect probability reconstruction algorithm according to an embodiment of the disclosure.

    [0061] As shown in FIG. 6, the device 10 includes: an excitation frequency determination module 100, an inclination angle determination module 200, a reference data acquisition module 300, a signal data acquisition module 400 and a defect imaging module 500.

    [0062] The excitation frequency determination module 100 is configured to determine the excitation frequency of a transmitting air-coupled transducer according to a frequency dispersion curve of guided waves and the thickness of a to-be-detected piece. The inclination angle determination module 200 is configured to determine the group velocity of an antisymmetric mode according to the excitation frequency, and determine the inclination angle of the transmitting air-coupled transducer and a receiving air-coupled transducer according to the Snell law. The reference data acquisition module 300 is configured to place a transmitting air-coupled transducer and a receiving air-coupled transducer on the same-side upper surface of a defect-free test piece according to the inclination angle, and obtain an initial waveform of the defect-free test piece as reference data by adopting a same-side penetration method. The signal data acquisition module 400 is configured to place the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece according to the inclination angle, and obtain N groups of signal data of the to-be-detected piece by adopting the rotary scanning method. The defect imaging module 500 is configured to compare the reference data with the N groups of signal data, solve N signal change correlation coefficients, process the N signal change correlation coefficients based on the defect probability reconstruction algorithm to obtain the defect distribution probability on the to-be-detected piece, and carry out defect imaging on a rotating coverage area of the transmitting air-coupled transducer and the receiving air-coupled transducer according to the defect distribution probability.

    [0063] Further, in an embodiment of the disclosure, the signal data acquisition module 400 further includes: a placement unit, configured to place the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece according to the inclination angle; and an acquisition unit, configured to rotate the transmitting air-coupled transducer and the receiving air-coupled transducer on the same-side upper surface of the to-be-detected piece by adopting the rotary scanning method at preset angle intervals to collect new signals, and obtain the N groups of signal data when the transmitting air-coupled transducer and the receiving air-coupled transducer rotate by 360 degrees.

    [0064] Further, in an embodiment of the disclosure, the signal change correlation coefficient is as follows:

    [00011] ρ = C XY σ X σ Y

    [0065] ρ is the signal change correlation coefficient, C.sub.XY is covariance of X and Y, X is a reference data set, Y is signal data after a period of service time, and σ.sub.X and σ.sub.Y are standard deviations of X and Y.

    [0066] Further, in an embodiment of the disclosure, the defect distribution probability is the sum of all signal change effects of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer.

    [0067] Further, in an embodiment of the disclosure, the defect distribution probability is as follows:

    [00012] P ( x , y ) = .Math. i = 1 N P i ( x , y ) = .Math. i = 1 N A i [ β - R i ( x , y ) β - 1 ]

    [0068] P.sub.i(x,y) is defect distribution probability estimation of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in a multi-order symmetric mode S.sub.i, A.sub.i(x,y)=1−ρ.sub.i is a signal difference coefficient of the pair of the transmitting air-coupled transducer and the receiving air-coupled transducer in the multi-order symmetric mode S.sub.i, (β−R.sub.i(x,y))/(β−1) is a non-negative space distribution function of the multi-order symmetric mode S.sub.i, and the outline thereof is a group of ellipses.

    [0069] It should be noted that the explanation of the embodiment of the air-coupled ultrasonic detection method based on the defect probability reconstruction algorithm is also applicable to the device and is not repeated here.

    [0070] According to the air-coupled ultrasonic detection device based on the defect probability reconstruction algorithm provided by the embodiment of the disclosure, the accuracy of traditional air-coupled ultrasonic X and Y scanning detection is improved; compared with a complex imaging technology, the air-coupled ultrasonic detection device consumes less time; and under the same environment and measurement conditions, the change is caused by the generation of defects. According to the embodiment of the disclosure, the defect growth may be monitored by tracking increase of the signal difference relative to the normal condition, and the fact health monitoring of the to-be-detected piece is realized.

    [0071] In addition, terms “first” and “second” are only intended for description, but cannot be construed as indicating or implying relative importance or implicitly indicating the number of the specified technical features. Thus, the features defined with “first” and “second” may explicitly or implicitly include at least one of the features. In the description of the disclosure, the meaning of “a plurality of” is at least two, such as two and three, unless otherwise specifically defined.

    [0072] In the description of the present description, descriptions with reference to terms “one embodiment”, “some embodiments”, “example”, “specific example”, or “some examples” and the like mean that specific features, structures, materials, or characteristics described in combination with the embodiments or examples are included in at least one embodiment or example of the disclosure. In the present description, the schematic representations of the foregoing terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in a suitable manner in any one or more embodiments or examples. Moreover, those skilled in the art connect and combine different embodiments or examples described in the present description and features of different embodiments or examples without mutual conflict.

    [0073] Although the embodiments of the disclosure have been shown and described above, it can be understood that the foregoing embodiments are illustrative and are not intended to be understood as limiting the disclosure. A person of ordinary skill in the art may make changes, modifications, replacements and variations to the described embodiments without departing from the scope of the disclosure.