CN112147225A - Nonlinear wave detection method for underwater gate - Google Patents
Nonlinear wave detection method for underwater gate Download PDFInfo
- Publication number
- CN112147225A CN112147225A CN202010795436.XA CN202010795436A CN112147225A CN 112147225 A CN112147225 A CN 112147225A CN 202010795436 A CN202010795436 A CN 202010795436A CN 112147225 A CN112147225 A CN 112147225A
- Authority
- CN
- China
- Prior art keywords
- wave
- frequency
- waves
- gate
- transducer
- 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.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 40
- 230000006378 damage Effects 0.000 claims abstract description 47
- 235000019687 Lamb Nutrition 0.000 claims abstract description 38
- 230000005284 excitation Effects 0.000 claims abstract description 33
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 31
- 239000002184 metal Substances 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims abstract description 10
- 238000012800 visualization Methods 0.000 claims abstract description 7
- 238000005457 optimization Methods 0.000 claims abstract description 5
- 230000008878 coupling Effects 0.000 claims abstract description 4
- 238000010168 coupling process Methods 0.000 claims abstract description 4
- 238000005859 coupling reaction Methods 0.000 claims abstract description 4
- 238000006073 displacement reaction Methods 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 14
- 239000006185 dispersion Substances 0.000 claims description 6
- 230000003902 lesion Effects 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000035882 stress Effects 0.000 description 5
- 230000007797 corrosion Effects 0.000 description 4
- 238000005260 corrosion Methods 0.000 description 4
- 238000009659 non-destructive testing Methods 0.000 description 3
- 230000000903 blocking effect Effects 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 239000004567 concrete Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000002905 metal composite material Substances 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010257 thawing Methods 0.000 description 1
- 238000001363 water suppression through gradient tailored excitation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/069—Defect imaging, localisation and sizing using, e.g. time of flight diffraction [TOFD], synthetic aperture focusing technique [SAFT], Amplituden-Laufzeit-Ortskurven [ALOK] technique
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/015—Attenuation, scattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0234—Metals, e.g. steel
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/044—Internal reflections (echoes), e.g. on walls or defects
Abstract
The invention discloses a nonlinear wave detection method of an underwater gate, which comprises the steps of firstly determining the relation of the group velocity of lamb waves and the product of the frequency plate thickness through wave velocity calculation, then carrying out coupling optimization, reducing the energy of lamb waves leaking into water, increasing the energy of lamb waves in a gate plate, adopting 7 transducers to excite the waves to simultaneously reach a focusing point, generating the maximum scattered waves at the focusing point, improving the signal-to-noise ratio, obtaining a detection signal, extracting integral frequency multiplication components of the excitation waves from a frequency spectrum, adding the frequency multiplication components to be used as the damage degree of the focusing point, and sequentially transforming the focusing point to form a damage image of the whole gate; the invention removes the linear wave scattered by the structure, and only adopts the nonlinear wave scattered by the damage to detect the gate, so that the lamb wave can detect the metal plate with a complex structure and can carry out detection without influencing the water retaining of the gate; the damage visualization of the whole gate can be realized, and the damage degree of each point can be quantitatively represented.
Description
Technical Field
The invention belongs to the technical field of nondestructive testing, and particularly relates to a nonlinear wave detection method for an underwater gate.
Background
The hydraulic gate is a thin-wall metal structure and is easily damaged by aging, freeze thawing, corrosion and cavitation. The method is characterized in that nondestructive testing is regularly carried out on the hydraulic gate under the condition of not influencing water retaining, is the key for guaranteeing engineering safety and realizing scientific maintenance, and therefore an integral, rapid and nondestructive testing method for the underwater gate needs to be developed.
The existing hydraulic gate detection mainly comprises coating and corrosion amount detection, and the gate safety is evaluated through the combination of the erosion residual thickness of the gate and stress calculation. The existing detection method only can detect point by point, cannot evaluate the overall safety of the gate, and has hidden danger on the safety of the gate. Particularly for deep water gates of high dam and large-span gates with span of nearly one hundred meters, the risk of structural failure is higher, and an integral detection method is more needed.
The current mature overall detection methods comprise A scanning, B scanning and C scanning. A, exciting and receiving ultrasonic waves by a point of a sweep aiming at a plate, wherein the abscissa is time, and the ordinate is an ultrasonic signal, so that the defect in the plate thickness direction can be detected; b, scanning a vertical section of the plate along the thickness direction, wherein the probe moves along a line, the abscissa is the displacement of the probe, and the ordinate is the thickness of the plate; c-scan a horizontal section to a fixed depth of the plate, as shown in fig. 13, the probe is moved across the surface of the plate to image a signal at a particular depth. Therefore, the scanning process requires moving the probe, which is time-consuming and labor-consuming.
For metal plates, one of the more common methods of bulk inspection is guided wave inspection. When the guided wave encounters damage, scattering occurs. By collecting and analyzing the scattered waves, the information such as whether the damage occurs, the position of the damage, the size of the damage and the like can be reversely calculated. The guided wave detection has the advantages that: during detection, the position of the probe does not need to be moved, and large-area integral detection can be realized through the propagation, scattering and return of the guided wave in the plate. Meanwhile, when water exists on one side of the metal plate, the wave guide method still has good detection capability.
However, the gate is not a smooth metal plate, and structures such as stiffening, riveting and the like are often provided on the surface of the gate, and since the structures also cause scattering of guided waves, the scattering of the structures may submerge and damage the scattering, so that damage cannot be identified. Indoor experiments show that only linear scattering waves of damage can be used to detect damage far away from a structure (in the moment, linear scattering of the damage is not submerged by structural scattering), a detection blind area exists, and real damage is likely to occur near the structure (such as fatigue cracks at an opening). Therefore, although lamb waves can well detect smooth metal plates, it is difficult to apply the lamb waves to gate detection having a complicated configuration. Meanwhile, common damages of the gate, such as corrosion, fatigue, bolt looseness and the like, are difficult to detect through linear scattered waves.
In essence, the damage scattering removes linear components, and also includes nonlinear components. Linear scattering of the lesion, whose frequency is the same as the frequency of the constructive scattered wave, and nonlinear scattering of the lesion, whose frequency is an integer multiple of the constructive scattered wave. By time-frequency analysis, linear scattering and nonlinear scattering can be separated. Because the nonlinear scattering can only come from damage, the damage can be detected, positioned and quantified only by considering nonlinear scattering waves, the interference of the constructed linear scattering waves can be avoided, and the lamb waves can be applied to the detection of components with complex structures, so that the problem of integral nondestructive detection of the underwater gate is solved.
Disclosure of Invention
In order to solve the problems, the invention discloses a nonlinear wave detection method for an underwater gate, which can realize the whole damage visualization of the gate, quantitatively represent the damage degree of each point, is not interfered by the complex structure of the gate, has no blind area in detection and improves the detection accuracy.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a nonlinear wave detection method for an underwater gate comprises the following steps:
(1) wave velocity calculation
Determining the length, width, thickness, elastic modulus and Poisson ratio of the gate panel by looking up engineering data, field measurement and inversion by combining monitoring data; solving a wave equation of lamb waves:
where ω is the angular frequency, cLIs the transverse wave velocity, cTDrawing a frequency dispersion curve of lamb waves, and determining the relation between the phase velocity of the lamb waves and the product of the frequency plate thickness;
according to the relationship between phase velocity and group velocity:
determining the relation of the group velocity of lamb waves and the product of frequency and plate thickness;
(2) coupling optimization
Modulating the sine function by using a window function m (t) to generate an excitation signal:
the window function modulation is to ensure that the excitation signal is at the central frequency omega0Peaks with energy, in the dispersion curve of lamb waves, by looking up ω0The propagation velocity of the excitation signal is obtained.
Decomposing displacement of gate panel into scalar potential function by Helmholtz decompositionAnd the sum of the vector potential functions Ψ, assuming that water cannot transfer shear, its displacement is represented by a scalar potential only, resulting in the following equation of motion:
the waterless side of the gate panel is a free surface, the stress is 0, and the boundary condition is obtained:
the water retaining side of the gate panel has the conditions of stress balance and off-plane displacement continuity:
the united vertical type (4) - (6) solves the displacement of the water body, and represents the energy leaked into the water body by the lamb wave; the wave number in the metal plate and the water body has the following geometrical relationship:
the geometric relation (7) is used as a limiting condition, the conditional extreme value of the water displacement is solved through a Lagrange method, the optimal wave number kx and ky are found, the energy of lamb waves leaking into water is reduced, and the energy of the lamb waves in the gate plate is increased;
(3) single point focusing
On the water blocking surface, a linear array consisting of 7 transducers is pasted on the gate panel, and the serial numbers are 1-7 in sequence. The distance between every two transducers is equal, and the transducers can independently excite and receive lamb wave signals; selecting the center of the array as a coordinate origin O, and recording the coordinate of the focusing point as O', wherein
Where Δ t is the time for transducer i to excite earlier than transducer 4, and if Δ t is negative, it indicates that transducer i is delayed by Δ t than transducer 4. Setting the excitation time of each transducer according to the criterion, ensuring that excitation waves of each transducer simultaneously reach a focusing point, and generating maximum scattered waves at the focusing point;
(4) signal-to-noise ratio enhancement
After the excitation is completed, the transducers acquire the excitation of the adjacent 6 transducers; taking two adjacent transducers as an example, the distance between the two transducers is recorded as deltax; the waveform of the excitation signal after Δ x has been delivered is as follows:
when k is0,k1First two-stage Taylor expansion of A-mode wavenumberWhen the coefficient is calculated, the A mode direct wave excited by the transducer a is obtained;
when k is0,k1When the first two-order Taylor expansion coefficients of the S-mode wave number are taken, the S-mode direct wave excited by the transducer a is obtained;
removing direct waves of A and S modes from a signal received by the transducer b to obtain a detection signal of the transducer b;
(5) calculation of damage product
The detection signal contains two components, one is a linear wave, the frequency is equal to the excitation frequency and comes from damage linear scattering, structural scattering and edge scattering, and the other is a nonlinear wave, the frequency is an integral multiple of the excitation frequency and comes from damage nonlinear scattering; extracting frequency components of integral multiple frequency from the frequency spectrum, and adding the frequency spectrum energy of nonlinear waves to be used as the damage degree of the focusing point;
(6) lesion visualization
And (5) continuously changing the focus point, repeating the step 3-5 to obtain the damage degree of each point, and forming a damage image of the whole gate.
The invention has the beneficial effects that:
(1) detection can be carried out under the condition that the water retaining of the gate is not influenced;
(2) the device is not interfered by the complex structure of the gate, and has no blind area in detection;
(3) the damage visualization of the whole gate can be realized, and the damage degree of each point can be quantitatively represented;
(4) contact damage such as corrosion, bolt loosening, fatigue, etc. can be detected;
(5) the probe does not need to be moved during detection.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a graph showing the product of the phase velocity of lamb wave and the frequency plate thickness in step (1) of the present invention;
FIG. 3 is a graph showing the product of the group velocity of lamb waves and the frequency plate thickness in step (1) of the present invention;
FIG. 4 is a schematic diagram of the optimization of the excitation wavenumber in step (2) of the present invention;
FIG. 5 is a phase difference calculation chart of single-point focusing in step (3) of the present invention;
FIG. 6 is a schematic diagram of single point focusing in step (3) of the present invention;
FIG. 7 is a diagram of the original signal received by transducer b in step (4) of the present invention;
FIG. 8 shows the direct wave of the A-mode emitted from the transducer a in step (4) of the present invention;
FIG. 9 shows the S-mode direct wave emitted from the transducer a in step (4) of the present invention;
FIG. 10 is a diagram of the detection signal of the transducer b after the signal-to-noise ratio is increased in step (4) of the present invention;
FIG. 11 is a frequency spectrum of the detection signal of the transducer b in step (5) of the present invention;
FIG. 12 is an image of the damage of the entire gate in step (6) of the present invention;
FIG. 13 is a schematic diagram of the scanning methods A, B, and C in the background art.
Detailed Description
The present invention will be further illustrated with reference to the accompanying drawings and specific embodiments, which are to be understood as merely illustrative of the invention and not as limiting the scope of the invention.
As shown in the figure, the nonlinear wave detection method of the underwater gate comprises the following steps:
(1) wave velocity calculation
The length, width, thickness, elastic modulus and Poisson ratio of the gate panel are determined by means of looking up engineering data, field measurement, monitoring data inversion and the like. Solving a wave equation of lamb waves:
where ω is the angular frequency, cLIs the transverse wave velocity, cTIs the longitudinal wave velocity, h is the half plate thickness, and k is the wave number. Drawing a frequency dispersion curve of the lamb wave, and determining the relation between the phase velocity of the lamb wave and the product of the frequency and the plate thickness, as shown in FIG. 2:
according to the relationship between phase velocity and group velocity:
the relationship between the group velocity of lamb waves and the product of frequency and plate thickness is determined, as shown in FIG. 3:
(2) coupling optimization
Modulating the sine function by using a window function m (t) to generate an excitation signal:
the purpose of the window function modulation is to ensure that the excitation signal has a peak with energy at the center frequency ω 0, and in the dispersion curve of the lamb wave, the propagation velocity of the excitation signal can be obtained by querying ω 0, as shown in fig. 4.
Decomposing displacement of gate panel into scalar potential function by Helmholtz decompositionAnd the sum of the vector potential functions Ψ, assuming that water cannot transfer shear, its displacement is represented by a scalar potential only, resulting in the following equation of motion:
the waterless side of the gate panel is a free surface, the stress is 0, and the boundary condition is obtained:
the water retaining side of the gate panel has the conditions of stress balance and off-plane displacement continuity:
and (4) to (6) in a joint type, solving the displacement of the water body, and representing the energy of lamb wave leaking into the water body, wherein as shown in fig. 4, the metal plate and the wave number in the water body have the following geometrical relationship:
and (3) taking the geometric relation (7) as a limiting condition, solving the conditional extreme value of the water displacement by a Lagrange method, finding the optimal wave number kx and ky, reducing the energy leaked into the water by the lamb wave, and increasing the energy of the lamb wave in the gate plate.
(3) Single point focusing
On the water blocking surface, a linear array composed of 7 transducers (such as piezoelectric transducer and electromagnetic transducer) is pasted on the gate panel, and the serial numbers are 1-7. The distance between every two transducers is equal, and the transducers can independently excite and receive lamb wave signals. Selecting the center of the array as a coordinate origin O, and recording the coordinate of the focusing point as O', wherein
As shown in fig. 5, where Δ t is the time for transducer No. i to fire earlier than transducer No. 4, if Δ t is negative, it indicates that transducer No. i is delayed by Δ t than transducer No. 4. Setting the excitation time of each transducer according to the criterion can ensure that the excitation waves of each transducer simultaneously reach the focus point, and the maximum scattered waves are generated at the focus point, as shown in fig. 6.
(4) Signal-to-noise ratio enhancement
After the excitation is completed, the transducer will acquire the excitation of the adjacent 6 transducers. Due to the multi-modal phenomenon of lamb waves, the excitation of one transducer contains two wave packets of A and S, and the detection signal is very complex. Taking two adjacent transducers as an example, the received signal of one transducer is shown in FIG. 7
Let us note the separation of the two transducers as Δ x. The waveform of the excitation signal after Δ x has been delivered is as follows:
when k is0,k1When the first two-order Taylor expansion coefficients of the A-mode wave number are taken, the A-mode direct wave excited by the transducer a is obtained;
when k is0,k1When the first two-order Taylor expansion coefficients of the S-mode wave number are taken, the S-mode direct wave excited by the transducer a is obtained; in the signal received by the transducer b, the direct waves of the a and S modes are removed, and a detection signal of the transducer b is obtained, as shown in fig. 10.
Since the direct wave (linear wave) is removed, the energy ratio of the nonlinear wave is further improved in the detection signal. In the existing nonlinear wave extraction process, a great difficulty is that in a frequency spectrum, the peak value of a linear wave is too large, the energy is too high, so that the frequency spectrum peak of the nonlinear wave is not obvious, and the nonlinear wave cannot be extracted. The invention can highlight the frequency spectrum wave crest of the nonlinear wave through the operation of the step.
(5) Calculation of damage product
The detection signal contains two components, one is a linear wave with equal frequency and excitation frequency and comes from damage linear scattering, structural scattering and edge scattering, and the other is a nonlinear wave with integral multiple of the excitation frequency and comes from damage nonlinear scattering. Frequency components of an integral multiple of the frequency are extracted from the spectrum, and the spectral energy of the nonlinear wave is added as the degree of damage at the focusing point, as shown in fig. 11.
(6) Lesion visualization
And (5) continuously changing the focus point, and repeating the step 3-5 to obtain the damage degree of each point. The damage image shown in fig. 12 is formed, and the point with a high pixel is the damage position.
The invention removes the linear wave scattered by the structure, and only adopts the nonlinear wave scattered by the damage to detect the gate, so that the lamb wave can detect the metal plate with a complex structure and can carry out detection when the water retaining of the gate is not influenced; the damage visualization of the whole gate can be realized, and the damage degree of each point can be quantitatively represented.
The description of the points is as follows:
(1) the phased array may be not only linear, but also circular, annular, rectangular, and the like. The directional detection can be realized not only by a phase modulation method, but also by other methods such as frequency modulation, or by using a transducer capable of being directionally excited in hardware.
(2) By replacing Δ x in step 4 with the path length of the scattered waves through the edge of the shutter plate, the scattered waves at the edge of the shutter plate can also be removed, thereby further improving the composition ratio of the nonlinear waves. The method for improving the signal-to-noise ratio also comprises a plurality of signal processing methods such as wavelet transformation, Fourier transformation, duffing oscillator and the like, and can replace the method in the step 4.
(3) The invention is not only suitable for gates, but also suitable for thin-wall structures such as airplanes, ship shells, bridge decks, gas storage tanks, pipelines and the like, the structures can be made of metal materials or composite materials, and the structures can be filled with liquid and can also be wrapped by soil, concrete and water (such as seabed petroleum pipelines and underground buried pipes).
The technical means disclosed in the invention scheme are not limited to the technical means disclosed in the above embodiments, but also include the technical scheme formed by any combination of the above technical features.
Claims (1)
1. A nonlinear wave detection method of an underwater gate is characterized by comprising the following steps: the method comprises the following steps:
(1) wave velocity calculation
Determining the length, width, thickness, elastic modulus and Poisson ratio of the gate panel by looking up engineering data, field measurement and inversion by combining monitoring data; solving a wave equation of lamb waves:
whereinOmega is the angular frequency, cLIs the transverse wave velocity, cTDrawing a frequency dispersion curve of lamb waves, and determining the relation between the phase velocity of the lamb waves and the product of the frequency plate thickness;
according to the relationship between phase velocity and group velocity:
determining the relation of the group velocity of lamb waves and the product of frequency and plate thickness;
(2) coupling optimization
Modulating the sine function by using a window function m (t) to generate an excitation signal:
the window function modulation is to ensure that the excitation signal is at the central frequency omega0Peaks with energy, in the dispersion curve of lamb waves, by looking up ω0Obtaining the propagation speed of the excitation signal;
decomposing displacement of gate panel into scalar potential function by Helmholtz decompositionAnd the sum of the vector potential functions Ψ, assuming that water cannot transfer shear, its displacement is represented by a scalar potential only, resulting in the following equation of motion:
the waterless side of the gate panel is a free surface, the stress is 0, and the boundary condition is obtained:
the water retaining side of the gate panel has the conditions of stress balance and off-plane displacement continuity:
the united vertical type (4) - (6) solves the displacement of the water body, and represents the energy leaked into the water body by the lamb wave; the wave number in the metal plate and the water body has the following geometrical relationship:
the geometric relation (7) is used as a limiting condition, the conditional extreme value of the water displacement is solved through a Lagrange method, the optimal wave number kx and ky are found, the energy of lamb waves leaking into water is reduced, and the energy of the lamb waves in the gate plate is increased;
(3) single point focusing
On the water retaining surface, a linear array consisting of 7 transducers is stuck on the gate panel, and the serial numbers are 1 to 7 in sequence; the distance between every two transducers is equal, and the transducers can independently excite and receive lamb wave signals; selecting the center of the array as a coordinate origin O, and recording the coordinate of the focusing point as O', wherein
Wherein, delta t is the time that the transducer I needs to excite earlier than the transducer 4, if delta t is a negative value, the transducer I needs to delay delta t than the transducer 4; setting the excitation time of each transducer according to the criterion, ensuring that excitation waves of each transducer simultaneously reach a focusing point, and generating maximum scattered waves at the focusing point;
(4) signal-to-noise ratio enhancement
After the excitation is completed, the transducers acquire the excitation of the adjacent 6 transducers; taking two adjacent transducers as an example, the distance between the two transducers is recorded as deltax; the waveform of the excitation signal after Δ x has been delivered is as follows:
when k is0,k1When the first two-order Taylor expansion coefficients of the A-mode wave number are taken, the A-mode direct wave excited by the transducer a is obtained;
when k is0,k1When the first two-order Taylor expansion coefficients of the S-mode wave number are taken, the S-mode direct wave excited by the transducer a is obtained;
removing direct waves of A and S modes from a signal received by the transducer b to obtain a detection signal of the transducer b;
(5) calculation of damage product
The detection signal contains two components, one is a linear wave, the frequency is equal to the excitation frequency and comes from damage linear scattering, structural scattering and edge scattering, and the other is a nonlinear wave, the frequency is an integral multiple of the excitation frequency and comes from damage nonlinear scattering; extracting frequency components of integral multiple frequency from the frequency spectrum, and adding the frequency spectrum energy of nonlinear waves to be used as the damage degree of the focusing point;
(6) lesion visualization
And (5) continuously changing the focus point, repeating the step 3-5 to obtain the damage degree of each point, and forming a damage image of the whole gate.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010795436.XA CN112147225A (en) | 2020-08-10 | 2020-08-10 | Nonlinear wave detection method for underwater gate |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010795436.XA CN112147225A (en) | 2020-08-10 | 2020-08-10 | Nonlinear wave detection method for underwater gate |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112147225A true CN112147225A (en) | 2020-12-29 |
Family
ID=73887857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010795436.XA Pending CN112147225A (en) | 2020-08-10 | 2020-08-10 | Nonlinear wave detection method for underwater gate |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112147225A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114166942A (en) * | 2021-11-25 | 2022-03-11 | 武汉大学 | Laser ultrasonic-based method for measuring defect burying depth between metal additive manufacturing layers |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040255678A1 (en) * | 2003-06-20 | 2004-12-23 | Yoshiaki Nagashima | Nondestructive inspection apparatus and nondestructive inspection method using elastic guided wave |
US20130279561A1 (en) * | 2012-03-06 | 2013-10-24 | Carnegie Mellon University | System and method for time reversal data communications on pipes using guided elastic waves |
CN103969339A (en) * | 2014-05-15 | 2014-08-06 | 厦门大学 | Nonlinear ultrasonic guided wave detection method and device for pipeline micro damage |
CN104297346A (en) * | 2014-09-11 | 2015-01-21 | 天津大学 | Nondestructive detection system of sheet metal by ultrasonic planar guided-wave and detection method thereof |
CN104407054A (en) * | 2014-12-17 | 2015-03-11 | 厦门大学 | Ultrasonic micro-damage location detection method and device based on collinear aliasing of Lamb waves |
CN105372330A (en) * | 2015-11-09 | 2016-03-02 | 北京工业大学 | Non-linear Lamb wave frequency mixing method for detecting microcrack in plate |
CN105842348A (en) * | 2016-04-07 | 2016-08-10 | 北京大学 | Piezoelectric transducer used for exciting and receiving non-dispersive ultrasonic guided waves, and production method thereof |
CN107290746A (en) * | 2017-06-08 | 2017-10-24 | 北京小米移动软件有限公司 | Material tests method and device |
CN110338846A (en) * | 2019-07-19 | 2019-10-18 | 河南科技大学第一附属医院 | Long bone cortex bone microcrack zone system and method based on non-linear ultrasonic guided wave |
CN210037711U (en) * | 2019-03-29 | 2020-02-07 | 广东省特种设备检测研究院珠海检测院 | Fatigue crack detection device for gantry crane boom structure |
CN110967400A (en) * | 2019-12-25 | 2020-04-07 | 电子科技大学 | Lamb wave tomography method for defects of large storage tank bottom plate |
CN111044613A (en) * | 2019-12-26 | 2020-04-21 | 武汉工程大学 | Metal plate micro-defect detection method based on nonlinear Lamb wave |
CN111103153A (en) * | 2018-10-25 | 2020-05-05 | 上海铁路通信有限公司 | Bogie structure safety monitoring device |
-
2020
- 2020-08-10 CN CN202010795436.XA patent/CN112147225A/en active Pending
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040255678A1 (en) * | 2003-06-20 | 2004-12-23 | Yoshiaki Nagashima | Nondestructive inspection apparatus and nondestructive inspection method using elastic guided wave |
US20130279561A1 (en) * | 2012-03-06 | 2013-10-24 | Carnegie Mellon University | System and method for time reversal data communications on pipes using guided elastic waves |
CN103969339A (en) * | 2014-05-15 | 2014-08-06 | 厦门大学 | Nonlinear ultrasonic guided wave detection method and device for pipeline micro damage |
CN104297346A (en) * | 2014-09-11 | 2015-01-21 | 天津大学 | Nondestructive detection system of sheet metal by ultrasonic planar guided-wave and detection method thereof |
CN104407054A (en) * | 2014-12-17 | 2015-03-11 | 厦门大学 | Ultrasonic micro-damage location detection method and device based on collinear aliasing of Lamb waves |
CN105372330A (en) * | 2015-11-09 | 2016-03-02 | 北京工业大学 | Non-linear Lamb wave frequency mixing method for detecting microcrack in plate |
CN105842348A (en) * | 2016-04-07 | 2016-08-10 | 北京大学 | Piezoelectric transducer used for exciting and receiving non-dispersive ultrasonic guided waves, and production method thereof |
CN107290746A (en) * | 2017-06-08 | 2017-10-24 | 北京小米移动软件有限公司 | Material tests method and device |
CN111103153A (en) * | 2018-10-25 | 2020-05-05 | 上海铁路通信有限公司 | Bogie structure safety monitoring device |
CN210037711U (en) * | 2019-03-29 | 2020-02-07 | 广东省特种设备检测研究院珠海检测院 | Fatigue crack detection device for gantry crane boom structure |
CN110338846A (en) * | 2019-07-19 | 2019-10-18 | 河南科技大学第一附属医院 | Long bone cortex bone microcrack zone system and method based on non-linear ultrasonic guided wave |
CN110967400A (en) * | 2019-12-25 | 2020-04-07 | 电子科技大学 | Lamb wave tomography method for defects of large storage tank bottom plate |
CN111044613A (en) * | 2019-12-26 | 2020-04-21 | 武汉工程大学 | Metal plate micro-defect detection method based on nonlinear Lamb wave |
Non-Patent Citations (9)
Title |
---|
F. CIAMPA等: "Nonlinear imaging of damage in composite structures using sparse ultrasonic sensor arrays", STRUCTURAL CONTROL AND HEALTH MONITORING, 19 July 2016 (2016-07-19), pages 1 - 13 * |
JAN HETTLER等: "Linear and Nonlinear Guided Wave Imaging of Impact Damage in CFRP Using a Probabilistic Approach", MATERIALS, 7 November 2016 (2016-11-07), pages 1 - 17 * |
刘镇清: "无损检测中的超声兰姆波", 无损检测, no. 09, 31 December 1999 (1999-12-31), pages 28 - 42 * |
张锐等: "低频多模式超声兰姆波超薄弹性层定征方法_", 计量学报, vol. 21, no. 1, 31 January 2000 (2000-01-31), pages 59 - 67 * |
李仕俊;廖非易;高敏;张胤;林媛;: "金属铝板弯折形变积累的超声兰姆波无损检测", 电子元件与材料, no. 06, 31 December 2015 (2015-12-31), pages 82 - 85 * |
王帮峰;李迎;施益峰;: "复合材料结构健康主动监测中激励信号的优化", 南京航空航天大学学报, no. 05, 20 October 2006 (2006-10-20), pages 89 - 94 * |
王高平;石政;吴明丽;魏鸣;: "利用Lamb波的铝板损伤识别实验研究", 内燃机与配件, no. 03, 26 January 2018 (2018-01-26), pages 56 - 58 * |
郑阳等: "采用超声导波阵列技术研究板类结构大面积检测", 声学学报, 31 January 2023 (2023-01-31), pages 71 - 79 * |
郭鹏;田振华;李鸿源;徐鸿;: "板中Lamb波与管中纵向模态的传播特性的对比分析", 无损检测, no. 04, 10 April 2017 (2017-04-10), pages 48 - 54 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114166942A (en) * | 2021-11-25 | 2022-03-11 | 武汉大学 | Laser ultrasonic-based method for measuring defect burying depth between metal additive manufacturing layers |
CN114166942B (en) * | 2021-11-25 | 2023-08-15 | 武汉大学 | Method for measuring interlayer defect burial depth in metal additive manufacturing based on laser ultrasonic |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ostachowicz et al. | Damage localisation in plate-like structures based on PZT sensors | |
Cawley et al. | Practical long range guided wave inspection-applications to pipes and rail | |
US20060203086A1 (en) | Signal processing arrangement | |
Cawley et al. | Practical long range guided wave inspection-managing complexity | |
US20050268720A1 (en) | Matrix switched phased array ultrasonic guided wave system | |
Wu et al. | Ultrasonic guided wave-based switch rail monitoring using independent component analysis | |
Hayashi et al. | Single mode extraction from multiple modes of Lamb wave and its application to defect detection | |
Angulo et al. | Mooring integrity management: Novel approaches towards in situ monitoring | |
CN112147225A (en) | Nonlinear wave detection method for underwater gate | |
Zhang et al. | Detection of damages in underwater metal plate using acoustic inverse scattering and image processing methods | |
Cawley | Practical guided wave inspection and applications to structural health monitoring | |
Naqiuddin et al. | Ultrasonic signal processing techniques for Pipeline: A review | |
Robert et al. | Assessment of real-time techniques for ultrasonic non-destructive testing | |
Mukherjee et al. | NDE of composite structures using microwave time reversal imaging | |
CN113406202B (en) | Structural surface defect detection method based on high-frequency Lamb wave frequency domain information | |
CN113298805B (en) | Structure surface defect detection method based on active Lamb wave acoustic emission | |
Dhutti et al. | Monitoring of critical metallic assets in oil and gas industry using ultrasonic guided waves | |
CN110568084B (en) | Method for extracting low signal-to-noise ratio guided wave signal reaching time suitable for guided wave transducer array | |
EP3983790B1 (en) | A method for detecting faults in plates using guided lamb waves | |
El Mountassir et al. | An ultrasonic guided waves based prognostic approach for predictive maintenance: Experimental study cases | |
Rose et al. | An introduction to ultrasonic guided waves | |
Norli et al. | Ultrasonic detection of crack defects in pipe samples with a 132-channel test scanner in gas | |
Philtron et al. | Guided wave phased array sensor tuning for improved defect detection and characterization | |
Fromme | Monitoring of corrosion damage using high-frequency guided ultrasonic waves | |
CN109342570B (en) | Silt submarine channel instability and slump prediction and early warning method |
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 |