CN114184682A - Dual-chaotic system detection method for weak ultrasonic guided wave signals - Google Patents

Dual-chaotic system detection method for weak ultrasonic guided wave signals Download PDF

Info

Publication number
CN114184682A
CN114184682A CN202111505537.XA CN202111505537A CN114184682A CN 114184682 A CN114184682 A CN 114184682A CN 202111505537 A CN202111505537 A CN 202111505537A CN 114184682 A CN114184682 A CN 114184682A
Authority
CN
China
Prior art keywords
chaotic
target signal
driving force
chaotic system
signal
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.)
Granted
Application number
CN202111505537.XA
Other languages
Chinese (zh)
Other versions
CN114184682B (en
Inventor
张伟伟
成梦菲
武静
马宏伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dongguan University of Technology
Original Assignee
Dongguan University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dongguan University of Technology filed Critical Dongguan University of Technology
Priority to CN202111505537.XA priority Critical patent/CN114184682B/en
Publication of CN114184682A publication Critical patent/CN114184682A/en
Priority to JP2022131323A priority patent/JP7309236B2/en
Application granted granted Critical
Publication of CN114184682B publication Critical patent/CN114184682B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/04Analysing solids
    • G01N29/07Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4409Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
    • G01N29/4418Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a model, e.g. best-fit, regression analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids

Landscapes

  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Acoustics & Sound (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The invention provides a dual chaotic system detection method of weak ultrasonic guided wave signals, which comprises the following steps: determining a target signal, constructing a first chaotic system, determining a first chaotic system threshold value, constructing a first chaotic detection system, identifying a target signal amplitude value, determining the reduction of the first chaotic system threshold value when considering phase influence, constructing a second chaotic system, constructing a second chaotic detection system, detecting and positioning the target signal by combining the first chaotic detection system, and quantitatively detecting the pipeline defect to determine the position and the size of the pipeline defect. The invention utilizes the double-chaos detection system to quickly and accurately detect the occurrence time of the defect echo target signal through the phase difference, saves the solving resource, improves the analysis efficiency and provides a direction for detecting small defects or long distance.

Description

Dual-chaotic system detection method for weak ultrasonic guided wave signals
Technical Field
The invention belongs to the field of ultrasonic guided wave nondestructive detection methods, and particularly relates to a dual chaotic system detection method for weak ultrasonic guided wave signals.
Background
The guided wave is considered to have unique advantages in detecting buried pipelines and pipelines with wrapping layers because of the advantages of long propagation distance, capability of simultaneously detecting defects on the surface and the inside of the pipeline and the like, however, for small-defect or long-distance defect detection, a defect echo obtained by testing is often very weak and even completely submerged in a noise signal.
In order to overcome the difficulty, researchers utilize the small defect of the sensitivity enhancement of the chaotic system to detect. However, signal identification using the chaos detection system can only determine whether the signal to be detected contains a defect echo, but cannot determine the occurrence time and size of the defect echo, which becomes a difficult problem to overcome in the field, so that the chaos detection method has very limited application in ultrasonic guided wave signal detection.
The invention discloses a method for identifying ultrasonic guided waves of an inclined crack pipeline by using an improved duffing chaotic system, which comprises the following steps: 1) exciting an ultrasonic guided wave signal on one side end face of the pipeline to enable the ultrasonic guided wave to traverse all positions of the pipeline; 2) receiving an ultrasonic guided wave echo signal at a receiving end near the excitation end and recording a time history curve propagated in an ultrasonic guided wave pipeline, wherein the ultrasonic guided wave echo signal comprises an end face echo signal, a noise signal and a bevel crack defect signal submerged in noise; 3) and detecting and analyzing the ultrasonic guided wave echo signals by using the improved duffing chaotic system, extracting and identifying the inclined crack defect information submerged in the noise, and obtaining the inclined crack defect condition of the whole pipeline. Although the method can detect weak defect echo signals, the amplitude and the occurrence time of the signals cannot be positioned, and quantitative analysis of defects cannot be realized.
The invention discloses a method for positioning and detecting small defects of ultrasonic guided waves by combining a TR system and a Duffing system, which is applied to CN 105954358B in China, and comprises the following steps: 1) exciting an ultrasonic guided wave signal on one end face of the pipeline to enable the ultrasonic guided wave to traverse all positions of the pipeline; 2) receiving an ultrasonic guided wave echo signal at the other end surface near the excitation end and recording a time history curve propagated in the ultrasonic guided wave pipeline; 3) detecting and analyzing the ultrasonic guided wave echo signals by using an improved duffing chaotic system, and extracting and identifying small defect information submerged in noise; 4) carrying out time reversal on the obtained signal with the small defect information and then carrying out secondary excitation; 5) an inverted signal is received. The chaos detecting system is an auxiliary means to judge whether the detecting signal contains defect signal to decide whether the time reversal detecting is necessary.
Therefore, in order to realize the detection, quantification and positioning technologies of the defect echo signals at the same time, it is very urgent and necessary to find a dual-chaos system detection method of the weak ultrasonic guided wave signals to solve the problem that the quantitative analysis of the defect signals cannot be performed in the field.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a method for detecting a double-chaotic system of a weak ultrasonic guided wave signal. The method comprises the steps of determining a target signal, constructing a first chaotic system, determining a first chaotic system threshold value, constructing a first chaotic detection system, identifying a target signal amplitude value, determining the reduction amount of the first chaotic system threshold value when the phase influence is considered, constructing a second chaotic system, constructing a second chaotic detection system, detecting and positioning the target signal, and quantitatively detecting the pipeline defect to determine the position and the size of the pipeline defect. The invention utilizes the double-chaos detection system to quickly and accurately detect the occurrence time of the defect echo target signal through phase difference compensation, saves solving resources, improves analysis efficiency, and provides direction for small defect or long-distance detection.
The invention provides a method for detecting a double-chaotic system of weak ultrasonic guided wave signals, which comprises the following steps:
s1, determining a target signal: setting a target signal as a modulated single-frequency signal s (t);
s2, constructing a first chaotic system and determining a first chaotic system threshold value Fc1
S21, constructing a first chaotic system by a Duffing equation:
Figure BDA0003404172460000021
wherein ω represents the first outer driving force circle frequency; c represents a first damping and is constant; x represents a displacement variable of the chaotic system; f represents the magnitude of the external driving force; t represents the evolution time of the chaotic system;
s22, under the condition of a first parameter, converting the external driving force amplitude F to obtain a first chaotic system threshold F converted from a periodic state to a chaotic statec1Taking the parameter as a first chaotic system threshold value Fc1A system determined by the first damping c and the first external driving force circular frequency omega is set as a first chaotic system;
s3, constructing a first chaotic detection system and identifying a target signal amplitude: a target signal is superposed on an external driving force item on the right side of the formula (2) to construct a first chaotic detection system, and a target signal part is superposed to realize the identification of the amplitude of the target signal;
s4, determining the threshold reduction amount of the first chaotic system when the phase influence is considered;
s41, assuming that the target signal in the step S1 is out of phase with the right external driving force in the formula (2) in the step S2
Figure BDA0003404172460000022
When in use
Figure BDA0003404172460000023
In time, the first chaotic detecting system in step S3 becomes:
Figure BDA0003404172460000024
wherein, ω is1Represents the target signal circular frequency, and12 pi f, f denotes the target signal frequency; n represents the number of troughs of the target signal; k is an integer; a represents the target signal amplitude; Δ t represents a target signal length, and Δ t ═ 2 π n/ω;
s42, when the phase influence is considered, the threshold decrease amount Δ F of the first chaotic system becomes:
Figure BDA0003404172460000031
when in use
Figure BDA0003404172460000032
When the target signal amplitude A is obtained, the threshold reduction quantity delta F of the first chaotic system obtains the maximum value;
s5, constructing a second chaotic system and constructing a second chaotic detection system, and performing target signal positioning by combining the second chaotic detection system and the first chaotic detection system;
s51, adding an initial phase to the external driving force of the first chaotic system
Figure BDA00034041724600000315
The phase compensation is performed, and equation (7) becomes:
Figure BDA0003404172460000033
s52, transforming the initial phase
Figure BDA0003404172460000034
Superposing the target signals, reducing the driving force amplitude value after the system enters the chaotic state, enabling the system to return to the periodic state again, recording the driving force amplitude value reduction quantity as a threshold value reduction quantity delta F of the first chaotic system, recording delta F as 0 if the system is still in the periodic state, and drawing
Figure BDA0003404172460000035
Curve when
Figure BDA0003404172460000036
In time, the threshold reduction amount Δ F of the first chaotic system takes a maximum value of (Δ F)maxWhen the phase is A, the initial phase is recorded
Figure BDA0003404172460000037
Is composed of
Figure BDA0003404172460000038
S53, constructing a second chaotic system: under the condition of a second parameter, the amplitude F of the external driving force is converted to obtain the second parameter of the system from the periodic state to the chaotic stateTwo chaos system threshold Fc2(ii) a Taking the parameter as a second chaotic system threshold value Fc2A second damping c2Second external driving force circular frequency omega2The determined system is set as a second chaotic system;
s54, constructing a mathematical model of the second chaotic detection system:
Figure BDA0003404172460000039
wherein,
Figure BDA00034041724600000310
a second external driving force compensation initial phase representing a second chaotic detection system;
s55, converting the second external driving force initial phase of the second chaotic detection system
Figure BDA00034041724600000311
Superposing the target signals, reducing the driving force amplitude value after the system enters the chaotic state, enabling the system to return to the periodic state again, recording the driving force amplitude value reduction quantity as a threshold value reduction quantity delta F of the second chaotic system, recording delta F as 0 if the system is still in the periodic state, and drawing
Figure BDA00034041724600000312
The curve marks the maximum value of the threshold value reduction quantity delta F of the second chaotic system
Figure BDA00034041724600000313
Is composed of
Figure BDA00034041724600000314
S56 peak generation time t of target signalxComprises the following steps:
Figure BDA0003404172460000041
wherein T represents the length of a signal to be detected;
initial time t of superposition of signals to be detected0Obtaining the real peak occurrence time t of the target signaldComprises the following steps:
Figure BDA0003404172460000042
and completing the detection and positioning of the target signal.
Further, the step S3 specifically includes the following steps:
s31, superimposing the target signal onto the external driving force on the right side of the Duffing equation, and the first chaotic detection system is:
Figure BDA0003404172460000043
s32, examining the superimposed target signal portion, including:
Figure BDA0003404172460000044
s33, wherein n in the formula (4) is not more than 20, and meets the limiting condition, the latter two terms are negligible at this time, and are as follows:
Figure BDA0003404172460000045
s34, when the Duffing equation is in the critical state of cycle-to-chaos transition, the guided wave signals are superposed, the system changes into the chaos state, and F in the formula (5) is reducedc1So that it enters the cycle again, the threshold decrease amount Δ F of the first chaotic system:
ΔF=A (6);
the identification of the target signal amplitude is achieved by using equation (6).
Preferably, in the step S1, the target signal S (t) satisfies a functional relationship:
Figure BDA0003404172460000046
preferably, the first parameter condition in step S2 is a first damping c-c of the first chaotic system1And the first external driving force circular frequency omega of the first chaotic system and the target signal circular frequency omega1Unity, i.e. ω ═ ω1
Preferably, the second parameter condition in step S53 is that the second damping of the second chaotic system is the same as the first damping of the first chaotic system, i.e. c2=c=c1And a second external driving force circular frequency omega of the second chaotic system2Satisfies the relationship:
12)T=π (10)。
preferably, the limiting condition in step S33 is that when the frequency of the signal to be detected is different from the driving force frequency of the Duffing equation by more than 3%, the signal to be detected will not cause the solution of the Duffing equation to change.
Preferably, in the step S41
Figure BDA0003404172460000055
The steps S51 and S54 perform target signal generation timing detection by driving force phase compensation of the dual chaotic detection system in relation to the target signal generation timing.
In a preferred embodiment of the present invention, the method of the present invention can be applied to the quantitative detection of pipe defects. It also includes the following steps:
s6, carrying out quantitative detection on the pipeline defects: determining the position and size of the pipeline defect based on the propagation speed of the ultrasonic guided wave and the reflection law of the ultrasonic guided wave vertically incident in the structure;
s61, when the propagation velocity of the ultrasonic guided wave is known as c, the real peak occurrence time t of the target signal is obtained in step S5dI.e. the moment of occurrence of the defect echo, the distance L from the defect to the signal excitation, i.e. the reception position, is calculatedx
Figure BDA0003404172460000051
S62, according to the reflection law of the ultrasonic guided wave vertically incident in the structure, obtaining a reflection coefficient R as follows:
Figure BDA0003404172460000052
wherein β represents a structural section loss rate;
s63, the reflection coefficient R is the defect echo amplitude AdA and incident wave amplitude A0The ratio of:
Figure BDA0003404172460000053
wherein K represents an attenuation correction coefficient, 0< K < 1;
s64, obtaining structural section loss ratio β from equations (15) and (16):
Figure BDA0003404172460000054
s65, based on the structure section loss rate beta and the distance L between the defect and the signal excitation or receiving positionxAnd determining the position and the size of the pipeline defect.
Compared with the prior art, the invention has the technical effects that:
1. the invention provides a double-chaos system detection method of weak ultrasonic guided wave signals, which takes a Duffing equation critical state converted from a period to a chaos state as a chaos detection system, superposes a target signal on a driving force item, the system is converted into the chaos state, and then the amplitude of the driving force is adjusted to return the system to the period state, and the amplitude of the target signal is determined by the amplitude adjustment quantity of the driving force; the amplitude of the weak signal is detected by using the chaotic system, and a basis is provided for the evaluation of the size of the target signal.
2. The invention provides a method for detecting a dual chaotic system of weak ultrasonic guided wave signals, which constructs two driving force with different frequenciesAnd the Duffing chaotic system meeting a certain relation carries out initial phase compensation on the driving force item constructed by the double chaotic system, simultaneously carries out weak signal detection and respectively draws the two systems
Figure BDA0003404172460000061
Determining the occurrence time of a target signal through the phase difference corresponding to the peak values of the two curves so as to determine the defect position; by using the double-chaos detection system, the occurrence time of the target signal is quickly and accurately detected through the phase difference, the solving resource is saved, the analysis efficiency is improved, and the proposed utilization
Figure BDA0003404172460000062
The curve is used for analyzing the signal amplitude and the occurrence time, and has the advantages of intuition and quickness.
3. According to the double-chaos system detection method for the weak ultrasonic guided wave signals, on the aspect of quantitative detection of the pipeline defects, weak defect echoes are caused by small defects, and meanwhile, the ultrasonic guided wave signals are propagated for a long distance, so that the defect echoes are very weak.
Drawings
Other features, objects and advantages of the present application will become more apparent upon reading of the following detailed description of non-limiting embodiments thereof, made with reference to the accompanying drawings.
FIG. 1 is a flow chart of a dual chaotic system detection method of weak ultrasonic guided wave signals of the present invention;
FIG. 2 is a graph of a target signal to be identified in accordance with the present invention;
FIG. 3 is a schematic diagram of the pipeline ultrasonic guided wave defect detection of the present invention;
FIG. 4a is a diagram of the experimental setup of an embodiment of the present invention;
FIG. 4b is an experimental schematic of one embodiment of the present invention;
FIG. 5a is a graph of an ultrasonic guided wave test signal for an intact pipe according to an embodiment of the present invention;
FIG. 5b is a diagram of an ultrasonic guided wave test signal for a damaged pipe according to an embodiment of the present invention;
FIG. 6a is a graph showing the effect of a perfect pipeline signal after being filtered by a Butterworth band-pass filter according to the present invention;
FIG. 6b is a diagram showing the effect of the signal of the damaged pipeline after being filtered by the Butterworth band-pass filter according to the present invention;
FIG. 7a is a graph of the first chaotic detection system for intact pipes, in accordance with the present invention;
FIG. 7b is a diagram of the second chaotic detection system for intact pipes, in accordance with the present invention;
FIG. 8a is a diagram of the first chaotic detection system for damaged pipelines according to the present invention;
FIG. 8b is a diagram of the identification result of the second chaotic detection system for damaged pipelines in accordance with the present invention.
Detailed Description
The present application will be described in further detail with reference to the following drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the relevant invention and not restrictive of the invention. It should be noted that, for convenience of description, only the portions related to the related invention are shown in the drawings.
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments with reference to the attached drawings.
Fig. 1 shows a dual chaotic system detection method of a weak ultrasonic guided wave signal of the present invention, which comprises the following steps:
s1, determining a target signal: setting the target signal as a single-frequency signal s (t) modulated by a Hanning window, and satisfying the functional relation:
Figure BDA0003404172460000071
wherein, ω is1Represents the target signal circular frequency, and12 pi f, f denotes the target signal frequency; n represents the number of troughs of the target signal; a represents the target signal amplitude; Δ t represents the target signal length, and Δ t ═ 2 π n/ω.
In one embodiment, the frequency f is 70kHz, and in order to match the chaotic system solution, unit transformation is performed, where 70kHz is 0.07(1/μ s), ω 12 pi · 0.07rad/μ s ≈ 0.439823rad/μ s. Taking a complete wave packet length as 142 μ s, a ═ 0.1mV, n ═ 10, Δ t ≈ 142(μ s), the target signal shape to be identified is shown in fig. 2.
S2, constructing a first chaotic system and determining a first chaotic system threshold value Fc1
S21, constructing a first chaotic system by a Duffing equation:
Figure BDA0003404172460000072
wherein ω represents the first outer driving force circle frequency; c represents a first damping and is constant; x represents a displacement variable of the chaotic system; f represents the magnitude of the external driving force; and t represents the evolution time of the chaotic system.
S22, under the condition of a first parameter, converting the amplitude F of the external driving force to obtain a first chaotic system threshold value F of the system from the periodic state to the chaotic statec1Taking the parameter as a first chaotic system threshold value Fc1A system determined by the first damping c and the first external driving force circular frequency omega is set as a first chaotic system; the first parameter condition is that a first damping c of the first chaotic system is c1And the first external driving force circular frequency omega of the first chaotic system and the target signal circular frequency omega1Unity, i.e. ω ═ ω1
In a specific embodiment, the first damping c is 0.4, and the critical value of the transition of the system from the periodic state to the chaotic state is obtained as Fc10.45781, the parameter is Fc1=0.45781,c=0.4,ω=ω1The determined system is set as a first chaotic system.
S3, constructing a first chaotic detection system and identifying a target signal amplitude: and (3) superposing the target signal to the external driving force on the right side of the formula (2) to construct a first chaotic detection system, superposing the target signal part, and realizing the identification of the target signal amplitude.
S31, superimposing the target signal onto the external driving force on the right side of the Duffing equation, and the first chaotic detection system is:
Figure BDA0003404172460000073
s32, examining the superimposed target signal portion, including:
Figure BDA0003404172460000081
s33, wherein n in the formula (4) is not more than 20, and the limiting condition is satisfied: when the frequency of the signal to be detected is different from the frequency of the driving force of the Duffing equation by more than 3 percent, the signal to be detected can not cause the change of the solution of the Duffing equation. At this time, the last two terms are negligible, among them:
Figure BDA0003404172460000082
s34, when the Duffing equation is in the critical state of cycle-to-chaos transition, the guided wave signals are superposed, the system changes into the chaos state, and F in the formula (5) is reducedc1So that it enters the cycle again, the threshold decrease amount Δ F of the first chaotic system:
ΔF=A (6);
the identification of the target signal amplitude is achieved by using equation (6).
And S4, determining the threshold reduction amount of the first chaotic system when the phase influence is considered.
S41, assuming that the target signal in the step S1 has a phase difference with the right side external driving force item in the formula (2) in the step S2
Figure BDA0003404172460000083
When in use
Figure BDA0003404172460000084
In time, the first chaotic detecting system in step S3 becomes:
Figure BDA0003404172460000085
wherein k is an integer;
Figure BDA0003404172460000086
related to the moment of occurrence of the target signal.
S42, when the phase influence is considered, the threshold decrease amount Δ F of the first chaotic system becomes:
Figure BDA0003404172460000087
when in use
Figure BDA0003404172460000088
And obtaining the maximum value of the threshold reduction quantity delta F of the first chaotic system, wherein the maximum value is the target signal amplitude A.
And S5, constructing a second chaotic system and constructing a second chaotic detection system, and performing target signal positioning by combining the second chaotic detection system and the first chaotic detection system.
S51, adding an initial phase to the external driving force of the first chaotic system
Figure BDA0003404172460000089
The phase compensation is performed, and equation (7) becomes:
Figure BDA00034041724600000810
s52, transforming the initial phase
Figure BDA0003404172460000091
Superposing the target signal, if the system enters the chaotic state, reducing the driving force amplitude to make the system return to the periodic state again, and recording the driving force amplitudeThe reduction is the threshold reduction amount delta F of the first chaotic system, if the system still keeps a periodic state, the delta F is recorded as 0, and the graph is drawn
Figure BDA0003404172460000092
Curve when
Figure BDA0003404172460000093
In time, the threshold reduction amount Δ F of the first chaotic system takes a maximum value of (Δ F)maxWhen the phase is A, the initial phase is recorded
Figure BDA0003404172460000094
Is composed of
Figure BDA0003404172460000095
S53, constructing a second chaotic system: under the condition of a second parameter, the external driving force amplitude F is converted to obtain a second chaotic system threshold value F for the system to be converted from the periodic state to the chaotic statec2(ii) a Taking the parameter as a second chaotic system threshold value Fc2A second damping c2Second external driving force circular frequency omega2The determined system is set as a second chaotic system; the second parameter condition in step S53 is that the second damping of the second chaotic system is the same as the first damping of the first chaotic system, i.e. c2=c=c1And a second external driving force circular frequency omega of the second chaotic system2Satisfies the relationship:
12)T=π (10)
wherein T represents the length of the signal to be detected.
In a specific embodiment, the second damping c2Obtaining a critical value F of the system from the periodic state to the chaotic state as 0.4c20.45609, the parameter is Fc2=0.45609,c2=0.4,
Figure BDA0003404172460000096
The determined system is set as a second chaotic system.
S54, constructing a mathematical model of the second chaotic detection system:
Figure BDA0003404172460000097
wherein,
Figure BDA0003404172460000098
and a second external driving force compensation initial phase of the second chaotic detecting system is shown.
S55, converting the second external driving force initial phase of the second chaotic detection system
Figure BDA0003404172460000099
Superposing the target signals, reducing the driving force amplitude value after the system enters the chaotic state, enabling the system to return to the periodic state again, recording the driving force amplitude value reduction quantity as a threshold value reduction quantity delta F of the second chaotic system, recording delta F as 0 if the system is still in the periodic state, and drawing
Figure BDA00034041724600000910
The curve marks the maximum value of the threshold value reduction quantity delta F of the second chaotic system
Figure BDA00034041724600000911
Is composed of
Figure BDA00034041724600000912
S56 peak generation time t of target signalxComprises the following steps:
Figure BDA00034041724600000913
initial time t of superposition of signals to be detected0Obtaining the real peak occurrence time t of the target signaldComprises the following steps:
Figure BDA0003404172460000101
and completing the detection and positioning of the target signal.
S6, carrying out quantitative detection on the pipeline defects: the position and the size of the pipeline defect are determined based on the propagation speed of the ultrasonic guided wave and the reflection law of the ultrasonic guided wave vertically incident in the structure, and a pipeline ultrasonic guided wave defect detection schematic diagram is shown in fig. 3.
S61, when the propagation velocity of the ultrasonic guided wave is known as c, the real peak occurrence time t of the target signal is obtained in step S5dI.e. the moment of occurrence of the defect echo, the distance L from the defect to the signal excitation, i.e. the reception position, is calculatedx
Figure BDA0003404172460000102
Since the reflection coefficient of the defect echo is in direct proportion to the size of the defect, the reflection coefficient can be used for solving the reflection coefficient after the defect echo amplitude is identified by using the first chaotic detection system, so that the size of the defect is obtained.
S62, according to the reflection law of the ultrasonic guided wave vertically incident in the structure, obtaining a reflection coefficient R as follows:
Figure BDA0003404172460000103
where β represents a structural section loss rate.
S63, the reflection coefficient R is the defect echo amplitude AdA and incident wave amplitude A0The ratio of:
Figure BDA0003404172460000104
where K represents the attenuation correction coefficient, 0< K < 1.
S64, obtaining structural section loss ratio β from equations (15) and (16):
Figure BDA0003404172460000105
s65, based on the structure section loss rate beta and the distance L between the defect and the signal excitation or receiving positionxAnd determining the position and the size of the pipeline defect.
In order to verify the effectiveness of the method, experimental studies were carried out on a seamless steel tube of 5m length, the equipment used is shown in fig. 4a and mainly comprises: the system comprises an arbitrary signal generator, a low-frequency amplifier and an oscilloscope, and adopts a piezoelectric transducer to excite and collect ultrasonic guided wave signals. The test system was set up as shown in fig. 4b and used a Hanning window modulated signal with n-10 and a center frequency of 70kHz as the excitation signal.
3 working conditions shown in table 1 are set, wherein the working condition 1 is a perfect pipeline, then a cutting machine is used for machining crack simulation defects at a position 3m away from an excitation end, and the cross section is reduced by 3.2% and 6.4% respectively to be the working condition 2 and the working condition 3.
Figure BDA0003404172460000111
TABLE 1
Fig. 5a and 5b show the test signal for a 3.2% loss in section of a sound pipe compared to a sound pipe, respectively, and it can be seen that the defect echo is not clearly seen in the loss pipe test signal compared to the sound pipe. On the left in FIGS. 5a and 5b is the incident wave ω1Right side is end face echo omega2And the middle part is a signal Cs to be detected.
In the conventional detection method, the influence of noise can be reduced in a filtering mode, so that defect echo information is obtained. According to the method, the test signal in the figure 4 is filtered by using a Butterworth band-pass filter, wherein the center frequency is 70kHz, the low-frequency cut-off frequency is 60kHz, and the high-frequency cut-off frequency is 80 kHz. The effect of the filtered intact pipe and the loss of 3.2% in section is shown in fig. 6a and 6b, respectively. In the fig. 6a and 6B diagrams the upper line is the original signal a and the lower line is the filtered signal B.
Obviously, the defect echo is still difficult to observe under the damage working condition after filtering, which mainly lies in that a small defect causes a weak defect echo, and meanwhile, the ultrasonic guided wave signal is propagated for a long distance, so that the defect echo is very weak.
The double-chaos detection system can effectively identify the small defect echo information:
when a signal between an incident wave under a perfect pipeline working condition and an end face echo (0.65 ms-1.65 ms) is intercepted as a signal to be detected (see fig. 5a), T is 1ms, because the signal does not contain a defect echo and only contains a noise signal, the signal is identified by using a double-chaos detection system, the results of a first chaos detection system and a second chaos detection system are respectively shown in fig. 7a and 7b, and therefore, no matter the first chaos detection system or the second chaos detection system, the value of delta F is very small and is distributed randomly, and the situation that no defect echo (target signal) exists in the signal to be detected is shown.
Similarly, signals between damage conditions (0.65ms to 1.65ms) are intercepted (see fig. 5b), and are identified by using a dual-chaos detection system, and the results of the first chaos detection system and the second chaos detection system are respectively shown in fig. 8a and 8b, from which it can be clearly seen that
Figure BDA0003404172460000112
The peak on the curve. Note that each chaotic detection system has multiple peaks, which are periodically related to the identification result and the phase, using equation (12), as long as the guarantee is made
Figure BDA0003404172460000113
And
Figure BDA0003404172460000114
satisfy the relationship
Figure BDA0003404172460000115
And (4) finishing. It can be seen that the requirement can be satisfied as long as the first peak is selected in both chaotic detection systems.
From the peaks in FIG. 8a, a defect echo peak of 0.04mV was identified, and the defect size was identified as a cross-sectional loss of 3.45% according to equation (17) with a relative error of 7.81%.
By using fig. 8a and 8b for defect localization, the wave velocity can be determined by the time difference between the incident wave and the end echo, and their propagation distance,
Figure BDA0003404172460000121
the initial moment of intercepting the signal is: t is t0=0.65ms。
The defect echo occurrence time is as follows:
Figure BDA0003404172460000122
the defect locations are:
Figure BDA0003404172460000123
the relative error is 5.63%.
The same method is adopted to identify another working condition, the defect size identification result is that the section loss is 6.61%, the relative error is 3.28%, the defect position identification result is 3.203m, the relative error is 6.7%, and the identification precision can meet the requirement.
The invention designs a dual chaotic system detection method of weak ultrasonic guided wave signals, which takes a Duffing equation critical state converted from a period to a chaotic state as a chaotic detection system, superposes a target signal on a driving force item, converts the system into the chaotic state, returns the system to the period state by adjusting the amplitude of the driving force, and determines the amplitude of the target signal by the amplitude adjustment amount of the driving force; the amplitude of the weak signal is detected by using the chaotic system, and a basis is provided for the evaluation of the size of the target signal; two Duffing chaotic systems with different driving force frequencies and meeting a certain relation are constructed, the initial phase of the driving force item of the double chaotic system is constructed, weak signal detection is simultaneously carried out, and the initial phase and the weak signal detection are respectively drawn for the two systems
Figure BDA0003404172460000124
Determining the occurrence time of a target signal through the phase difference corresponding to the peak values of the two curves so as to determine the defect position; by using the double-chaos detection system, the occurrence time of the target signal is quickly and accurately detected through the phase difference, thereby savingSolving resources, improving analysis efficiency, and the proposed utilization
Figure BDA0003404172460000125
The curve is used for analyzing the signal amplitude and the occurrence time, and has the advantages of intuition and quickness; in the quantitative detection of the pipeline defects, because the small defects cause weak defect echoes, and meanwhile, ultrasonic guided wave signals are propagated for a long distance, the defect echoes are very weak, the method can still accurately identify and is obviously superior to the traditional method, the detection efficiency of the pipeline ultrasonic guided wave detection on the small defects or the long distance is improved, and the method has important significance for improving the ultrasonic guided wave detection efficiency.
Finally, it should be noted that: although the present invention has been described in detail with reference to the above embodiments, it should be understood by those skilled in the art that: modifications and equivalents may be made thereto without departing from the spirit and scope of the invention and it is intended to cover in the claims the invention as defined in the appended claims.

Claims (8)

1. A dual chaotic system detection method of weak ultrasonic guided wave signals is characterized by comprising the following steps:
s1, determining a target signal: setting a target signal as a modulated single-frequency signal s (t);
s2, constructing a first chaotic system and determining a first chaotic system threshold value Fc1
S21, constructing a first chaotic system by a Duffing equation:
Figure FDA0003404172450000011
wherein ω represents the first outer driving force circle frequency; c represents a first damping and is constant; x represents a displacement variable of the chaotic system; f represents the magnitude of the external driving force; t represents the evolution time of the chaotic system;
s22, under the condition of the first parameterAnd transforming the amplitude F of the external driving force to obtain a first chaotic system threshold value F converted from a periodic state to a chaotic statec1Taking the parameter as a first chaotic system threshold value Fc1A system determined by the first damping c and the first external driving force circular frequency omega is set as a first chaotic system;
s3, constructing a first chaotic detection system and identifying a target signal amplitude: a target signal is superposed on an external driving force item on the right side of the formula (2) to construct a first chaotic detection system, and a target signal part is superposed to realize the identification of the amplitude of the target signal;
s4, determining the threshold reduction amount of the first chaotic system when the phase influence is considered;
s41, assuming that the target signal in the step S1 is out of phase with the right external driving force in the formula (2) in the step S2
Figure FDA0003404172450000012
When in use
Figure FDA0003404172450000013
In time, the first chaotic detecting system in step S3 becomes:
Figure FDA0003404172450000014
wherein, ω is1Represents the target signal circular frequency, and12 pi f, f denotes the target signal frequency; n represents the number of troughs of the target signal; k is an integer; a represents the target signal amplitude; Δ t represents a target signal length, and Δ t ═ 2 π n/ω;
s42, when the phase influence is considered, the threshold decrease amount Δ F of the first chaotic system becomes:
Figure FDA0003404172450000015
when in use
Figure FDA0003404172450000016
When the temperature of the water is higher than the set temperature,obtaining the maximum value of the threshold value reduction delta F of the first chaotic system, wherein the maximum value is a target signal amplitude A;
s5, constructing a second chaotic system and constructing a second chaotic detection system, and performing target signal positioning by combining the second chaotic detection system and the first chaotic detection system;
s51, adding an initial phase to the external driving force of the first chaotic system
Figure FDA0003404172450000017
The phase compensation is performed, and equation (7) becomes:
Figure FDA0003404172450000018
s52, transforming the initial phase
Figure FDA0003404172450000021
Superposing the target signals, reducing the driving force amplitude value after the system enters the chaotic state, enabling the system to return to the periodic state again, recording the driving force amplitude value reduction quantity as a threshold value reduction quantity delta F of the first chaotic system, recording delta F as 0 if the system is still in the periodic state, and drawing
Figure FDA0003404172450000022
Curve when
Figure FDA0003404172450000023
In time, the threshold reduction amount Δ F of the first chaotic system takes a maximum value of (Δ F)maxWhen the phase is A, the initial phase is recorded
Figure FDA0003404172450000024
Is composed of
Figure FDA0003404172450000025
S53, constructing a second chaotic system: under the condition of a second parameter, the external driving force amplitude F is converted to obtain a second chaos of the system from the periodic state to the chaos stateSystem threshold Fc2(ii) a Taking the parameter as a second chaotic system threshold value Fc2A second damping c2Second external driving force circular frequency omega2The determined system is set as a second chaotic system;
s54, constructing a mathematical model of the second chaotic detection system:
Figure FDA0003404172450000026
wherein,
Figure FDA0003404172450000027
a second external driving force compensation initial phase representing a second chaotic detection system;
s55, converting the second external driving force initial phase of the second chaotic detection system
Figure FDA0003404172450000028
Superposing the target signals, reducing the driving force amplitude value after the system enters the chaotic state, enabling the system to return to the periodic state again, recording the driving force amplitude value reduction quantity as a threshold value reduction quantity delta F of the second chaotic system, recording delta F as 0 if the system is still in the periodic state, and drawing
Figure FDA0003404172450000029
The curve marks the maximum value of the threshold value reduction quantity delta F of the second chaotic system
Figure FDA00034041724500000210
Is composed of
Figure FDA00034041724500000211
S56 peak generation time t of target signalxComprises the following steps:
Figure FDA00034041724500000212
wherein T represents the length of a signal to be detected;
initial time t of superposition of signals to be detected0Obtaining the real peak occurrence time t of the target signaldComprises the following steps:
Figure FDA00034041724500000213
and completing the detection and positioning of the target signal.
2. The dual chaotic system detection method for the weak ultrasonic guided wave signal according to claim 1, wherein the step S3 specifically includes the following steps:
s31, superimposing the target signal onto the external driving force on the right side of the Duffing equation, and the first chaotic detection system is:
Figure FDA0003404172450000031
s32, examining the superimposed target signal portion, including:
Figure FDA0003404172450000032
s33, wherein n in the formula (4) is not more than 20, and meets the limiting condition, the latter two terms are negligible at this time, and are as follows:
Figure FDA0003404172450000033
s34, when the Duffing equation is in the critical state of cycle-to-chaos transition, the guided wave signals are superposed, the system changes into the chaos state, and F in the formula (5) is reducedc1So that it enters the cycle again, the threshold decrease amount Δ F of the first chaotic system:
ΔF=A (6);
the identification of the target signal amplitude is achieved by using equation (6).
3. The method for detecting the dual chaotic system of the weak ultrasonic guided wave signal according to claim 1, wherein the target signal S (t) in the step S1 satisfies a functional relationship:
Figure FDA0003404172450000034
4. the method for detecting a dual chaotic system in a weak ultrasonic guided wave signal according to claim 1, wherein the first parameter condition in step S2 is a first damping c-c of a first chaotic system1And the first external driving force circular frequency omega of the first chaotic system and the target signal circular frequency omega1Unity, i.e. ω ═ ω1
5. The dual chaotic system detection method of the weak ultrasonic guided wave signal of claim 1, wherein the second parameter condition in the step S53 is that a second damping of the second chaotic system is the same as a first damping of the first chaotic system, i.e., c2=c=c1And a second external driving force circular frequency omega of the second chaotic system2Satisfies the relationship:
12)T=π (10)。
6. the dual chaotic system detection method of a weak ultrasonic guided wave signal according to claim 1, wherein the limiting condition in step S33 is that when the frequency of the signal to be detected is different from the driving force frequency of Duffing equation by more than 3%, the signal to be detected will not cause the solution of Duffing equation to change.
7. The method for detecting the dual chaotic system in the weak ultrasonic guided wave signal according to claim 1, wherein the step S41 is performed in
Figure FDA0003404172450000035
The steps S51 and S54 perform target signal generation timing detection by driving force phase compensation of the dual chaotic detection system in relation to the target signal generation timing.
8. The dual chaotic system detection method of weak ultrasonic guided wave signals according to claim 1, further comprising the steps of:
s6, carrying out quantitative detection on the pipeline defects: determining the position and size of the pipeline defect based on the propagation speed of the ultrasonic guided wave and the reflection law of the ultrasonic guided wave vertically incident in the structure;
s61, when the propagation velocity of the ultrasonic guided wave is known as c, the real peak occurrence time t of the target signal is obtained in step S5dI.e. the moment of occurrence of the defect echo, the distance L from the defect to the signal excitation, i.e. the reception position, is calculatedx
Figure FDA0003404172450000041
S62, according to the reflection law of the ultrasonic guided wave vertically incident in the structure, obtaining a reflection coefficient R as follows:
Figure FDA0003404172450000042
wherein β represents a structural section loss rate;
s63, the reflection coefficient R is the defect echo amplitude AdA and incident wave amplitude A0The ratio of:
Figure FDA0003404172450000043
wherein K represents an attenuation correction coefficient, 0< K < 1;
s64, obtaining structural section loss ratio β from equations (15) and (16):
Figure FDA0003404172450000044
s65, based on the structure section loss rate beta and the distance L between the defect and the signal excitation or receiving positionxAnd determining the position and the size of the pipeline defect.
CN202111505537.XA 2021-12-10 2021-12-10 Double-chaotic system detection method for weak ultrasonic guided wave signals Active CN114184682B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202111505537.XA CN114184682B (en) 2021-12-10 2021-12-10 Double-chaotic system detection method for weak ultrasonic guided wave signals
JP2022131323A JP7309236B2 (en) 2021-12-10 2022-08-19 Dual Chaos Analysis System Detection Method for Weak Ultrasonic Guided Wave Signals

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111505537.XA CN114184682B (en) 2021-12-10 2021-12-10 Double-chaotic system detection method for weak ultrasonic guided wave signals

Publications (2)

Publication Number Publication Date
CN114184682A true CN114184682A (en) 2022-03-15
CN114184682B CN114184682B (en) 2024-04-30

Family

ID=80604333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111505537.XA Active CN114184682B (en) 2021-12-10 2021-12-10 Double-chaotic system detection method for weak ultrasonic guided wave signals

Country Status (2)

Country Link
JP (1) JP7309236B2 (en)
CN (1) CN114184682B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115097009A (en) * 2022-06-30 2022-09-23 东莞理工学院 Frequency domain regulation-based ultrasonic guided wave phased array CFRP defect detection method
CN115713844A (en) * 2022-10-18 2023-02-24 兰州交通大学 Alarm method and system

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117168665B (en) * 2023-08-14 2024-08-06 哈尔滨工业大学 Lyapunov index dictionary-based space-coupled ultrasonic Lamb wave acoustic time feature extraction method

Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1324069A1 (en) * 2001-12-21 2003-07-02 STMicroelectronics S.r.l. A system for detecting distances using chaotic signals
CN101256230A (en) * 2007-11-21 2008-09-03 北京理工大学 Continuous wave chaos phase-modulation radio fuse detection system and demodulation method
CN101650428A (en) * 2009-09-04 2010-02-17 西北工业大学 Method for detecting chaotic oscillator of submarine weak target signal
GB201012597D0 (en) * 2010-07-28 2010-09-08 Guided Ultrasonics Ltd Processing signals acquired during guided wave testing
JP2012149980A (en) * 2011-01-19 2012-08-09 Hitachi Engineering & Services Co Ltd Method and device for inspecting guide wave
KR20120122440A (en) * 2011-04-29 2012-11-07 원광대학교산학협력단 Ultrasonic nondestructive inspection device and ultrasonic nondestructive inspection method
CN103323529A (en) * 2013-05-23 2013-09-25 暨南大学 Method for identifying ultrasonic guided wave of oblique crack pipeline through utilizing improved Duffing chaotic system
CN104155364A (en) * 2014-07-31 2014-11-19 暨南大学 Method and device for identifying defect position of pipeline based on phase locus partitioned by round window
CN104462695A (en) * 2014-12-12 2015-03-25 燕山大学 Weak signal detection method based on double-coupling Duffing vibrators and scale varying
CN105675718A (en) * 2016-02-11 2016-06-15 广东工业大学 Ultrasonic guided-wave method for identifying pipeline cracks by utilizing cubic van der pol chaos system
CN105954358A (en) * 2016-05-04 2016-09-21 中国计量大学 TR and Duffing system combined ultrasonic guided wave small defect positioning and detecting method
JP2017106828A (en) * 2015-12-10 2017-06-15 日本電信電話株式会社 Wave reception method and wave reception device
CN107228905A (en) * 2017-06-02 2017-10-03 东莞理工学院 Ultrasonic guided wave signals detection method based on bistable system
CN109521099A (en) * 2018-01-25 2019-03-26 四川升拓检测技术股份有限公司 Mobile lossless detection method based on sound pick-up phase difference phase battle array
CN110274956A (en) * 2019-07-10 2019-09-24 广东海洋大学 A kind of pipe ultrasonic guided wave detecting method based on Kolmogorov entropy
CN110412132A (en) * 2019-07-01 2019-11-05 合肥通用机械研究院有限公司 Based on the multi-modal ultrasonic phase array detection system of multi-mode and nonlinear imaging method
US20200191754A1 (en) * 2015-08-12 2020-06-18 Chevron U.S.A. Inc. Detection, monitoring, and determination of location of changes in metallic structures using multimode acoustic signals
CN111537607A (en) * 2020-05-07 2020-08-14 太原科技大学 Ultrasonic guided wave defect detection method based on Poincare section mutation characteristic
CN111781275A (en) * 2020-06-16 2020-10-16 东莞理工学院 Steel rail ultrasonic guided wave defect identification and positioning method and device based on Lyapunov index
CN111781276A (en) * 2020-06-16 2020-10-16 东莞理工学院 Steel rail ultrasonic guided wave defect identification and positioning method and device based on fractal dimension
CN112697881A (en) * 2020-12-09 2021-04-23 东莞理工学院 Steel rail guided wave defect identification and positioning method, device and system based on K-S entropy
CN113567560A (en) * 2021-08-06 2021-10-29 东南大学 Ultrasonic guided wave-based damage detection method for pipeline with auxiliary structure

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104155368A (en) * 2014-07-31 2014-11-19 暨南大学 Method and device for identifying cracking angle of pipeline based on phase locus partitioned by round window
CN104777222A (en) * 2015-03-30 2015-07-15 暨南大学 Pipeline defect identification and visualization method based on three-dimensional phase trajectory of Duffing system
CN104792865B (en) * 2015-03-30 2017-07-04 暨南大学 The small defect recognition of pipe ultrasonic guided wave and localization method based on dimension
CN108680643A (en) * 2018-03-29 2018-10-19 东莞理工学院 A kind of new-type ultrasonic guided wave detecting method
CN112505155B (en) * 2020-12-09 2022-04-19 暨南大学 Pipeline guided wave damage identification and positioning method, device and system based on information entropy

Patent Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1324069A1 (en) * 2001-12-21 2003-07-02 STMicroelectronics S.r.l. A system for detecting distances using chaotic signals
CN101256230A (en) * 2007-11-21 2008-09-03 北京理工大学 Continuous wave chaos phase-modulation radio fuse detection system and demodulation method
CN101650428A (en) * 2009-09-04 2010-02-17 西北工业大学 Method for detecting chaotic oscillator of submarine weak target signal
GB201012597D0 (en) * 2010-07-28 2010-09-08 Guided Ultrasonics Ltd Processing signals acquired during guided wave testing
JP2012149980A (en) * 2011-01-19 2012-08-09 Hitachi Engineering & Services Co Ltd Method and device for inspecting guide wave
KR20120122440A (en) * 2011-04-29 2012-11-07 원광대학교산학협력단 Ultrasonic nondestructive inspection device and ultrasonic nondestructive inspection method
CN103323529A (en) * 2013-05-23 2013-09-25 暨南大学 Method for identifying ultrasonic guided wave of oblique crack pipeline through utilizing improved Duffing chaotic system
CN104155364A (en) * 2014-07-31 2014-11-19 暨南大学 Method and device for identifying defect position of pipeline based on phase locus partitioned by round window
CN104462695A (en) * 2014-12-12 2015-03-25 燕山大学 Weak signal detection method based on double-coupling Duffing vibrators and scale varying
US20200191754A1 (en) * 2015-08-12 2020-06-18 Chevron U.S.A. Inc. Detection, monitoring, and determination of location of changes in metallic structures using multimode acoustic signals
JP2017106828A (en) * 2015-12-10 2017-06-15 日本電信電話株式会社 Wave reception method and wave reception device
CN105675718A (en) * 2016-02-11 2016-06-15 广东工业大学 Ultrasonic guided-wave method for identifying pipeline cracks by utilizing cubic van der pol chaos system
CN105954358A (en) * 2016-05-04 2016-09-21 中国计量大学 TR and Duffing system combined ultrasonic guided wave small defect positioning and detecting method
CN107228905A (en) * 2017-06-02 2017-10-03 东莞理工学院 Ultrasonic guided wave signals detection method based on bistable system
CN109521099A (en) * 2018-01-25 2019-03-26 四川升拓检测技术股份有限公司 Mobile lossless detection method based on sound pick-up phase difference phase battle array
CN110412132A (en) * 2019-07-01 2019-11-05 合肥通用机械研究院有限公司 Based on the multi-modal ultrasonic phase array detection system of multi-mode and nonlinear imaging method
CN110274956A (en) * 2019-07-10 2019-09-24 广东海洋大学 A kind of pipe ultrasonic guided wave detecting method based on Kolmogorov entropy
CN111537607A (en) * 2020-05-07 2020-08-14 太原科技大学 Ultrasonic guided wave defect detection method based on Poincare section mutation characteristic
CN111781275A (en) * 2020-06-16 2020-10-16 东莞理工学院 Steel rail ultrasonic guided wave defect identification and positioning method and device based on Lyapunov index
CN111781276A (en) * 2020-06-16 2020-10-16 东莞理工学院 Steel rail ultrasonic guided wave defect identification and positioning method and device based on fractal dimension
CN112697881A (en) * 2020-12-09 2021-04-23 东莞理工学院 Steel rail guided wave defect identification and positioning method, device and system based on K-S entropy
CN113567560A (en) * 2021-08-06 2021-10-29 东南大学 Ultrasonic guided wave-based damage detection method for pipeline with auxiliary structure

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
WU JING ET. AL.: "Detection of Arbitrary Frequency Ultrasonic Guided Wave Signals Based on the Time-Shift Duffing Oscillator", IEEE ACCESS, 21 July 2021 (2021-07-21), pages 95823 - 95831, XP011865868, DOI: 10.1109/ACCESS.2021.3094561 *
张伟伟;成梦菲;赵子龙: "利用Duffing***周期跳跃的弱超声导波识别", 声学学报, 15 July 2021 (2021-07-15) *
张伟伟;王志华;马宏伟;: "含缺陷管道超声导波检测信号的相关性分析", 暨南大学学报(自然科学与医学版), no. 03, 15 June 2009 (2009-06-15) *
戴华平;胡红亮;王旭;王玉涛;钱嘉伟;: "基于耦合同步混沌***的超声波测距方法", 计算机工程, no. 05, 15 May 2013 (2013-05-15) *
时培明;孙彦龙;韩东颖;: "基于双耦合混沌振子变尺度微弱信号检测方法研究", 计量学报, no. 03, 22 May 2016 (2016-05-22) *
武静: "基于杜芬方程的超声导波检测方法研究", 中国优秀硕士学位论文全文数据库 基础科学辑, 15 August 2014 (2014-08-15) *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115097009A (en) * 2022-06-30 2022-09-23 东莞理工学院 Frequency domain regulation-based ultrasonic guided wave phased array CFRP defect detection method
CN115097009B (en) * 2022-06-30 2024-05-24 东莞理工学院 Ultrasonic guided wave phased array CFRP defect detection method based on frequency domain regulation and control
CN115713844A (en) * 2022-10-18 2023-02-24 兰州交通大学 Alarm method and system

Also Published As

Publication number Publication date
JP7309236B2 (en) 2023-07-18
CN114184682B (en) 2024-04-30
JP2023086662A (en) 2023-06-22

Similar Documents

Publication Publication Date Title
CN114184682A (en) Dual-chaotic system detection method for weak ultrasonic guided wave signals
CN106770664A (en) A kind of method that edge defect detection is improved based on total focus imaging algorithm
CN102043015B (en) Ultrasonic guided wave device and method for detecting defect at rail bottom of steel rail at long distance
CN105954358B (en) A kind of small defect location detection method of supersonic guide-wave that TR is combined with Duffing systems
CN109958432B (en) Method and device for evaluating cementing quality of well cementation II interface by utilizing ultrasonic echo logging
CN102980942B (en) Metal pipeline detection method
CN108802203B (en) rod-shaped member internal defect positioning method based on multi-mode technology
CN107121497A (en) The ultrasonic guided wave detecting method of accidental resonance characteristic based on Duffing systems
US20020194916A1 (en) Method for inspecting clad pipe
CN107153095A (en) The health detecting system and its method of work of a kind of composite plate structure
CN106770668B (en) Method for detecting quality of single-hole foundation pile by acoustic transmission method
CN113533510B (en) Rail fatigue micro-crack identification method and device
He et al. Quantitative detection of surface defect using laser-generated Rayleigh wave with broadband local wavenumber estimation
CN107576726B (en) Damage distinguishing and damage extension identifying method for guided wave detection
CN108195934B (en) Ultrasonic guided wave detection frequency optimization method based on time-frequency analysis
CN204575596U (en) A kind of spiral welded pipe twist and warping guided wave detection device adopting spiral transduction device
CN113406202B (en) Structural surface defect detection method based on high-frequency Lamb wave frequency domain information
CN114755300A (en) Defect positioning quantitative detection method based on ultrasonic nondestructive detection
Masserey et al. IN‐SITU MONITORING OF FATIGUE CRACK GROWTH AT FASTENER HOLES USING RAYLEIGH‐LIKE WAVES
CN109085238B (en) Method for identifying welding seam and hoop reflection signals in torsional mode guided wave pipeline detection
CN113466336A (en) Weld defect detection method and device, electronic equipment and storage medium
Hesse et al. A single probe spatial averaging technique for guided waves and its application to surface wave rail inspection
Gauthier et al. EMAT generation of horizontally polarized guided shear waves for ultrasonic pipe inspection
CN117761165B (en) Pipeline crack positioning method based on electromagnetic ultrasonic array torsion guided wave
CN114689689B (en) Ultrasonic guided wave long-distance detection method based on pipeline axial local scanning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant