CN107505396A - A kind of structural damage on-line real time monitoring method and system - Google Patents
A kind of structural damage on-line real time monitoring method and system Download PDFInfo
- Publication number
- CN107505396A CN107505396A CN201710818527.9A CN201710818527A CN107505396A CN 107505396 A CN107505396 A CN 107505396A CN 201710818527 A CN201710818527 A CN 201710818527A CN 107505396 A CN107505396 A CN 107505396A
- Authority
- CN
- China
- Prior art keywords
- electrical contact
- pressing plate
- measurement point
- frequency
- energy value
- 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
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/048—Marking the faulty objects
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention discloses a kind of structural damage on-line real time monitoring method and system, this method or the system perpendicular to the microvibration response signal of structure to direction by using set Empirical Mode Decomposition algorithm, first natural frequency and second-order intrinsic frequency are obtained, and judges whether first natural frequency and second-order intrinsic frequency are equal with the intrinsic frequency of non-damaged structure；If then determine that body structure surface does not damage；If otherwise determining, body structure surface damages, and calculates the energy value of different slip number difference measurement points；Whether the energy value for comparing energy value corresponding measurement point with non-damaged structure of the different slip a certain measurement points of number is all equal；If then determine that the measurement point is not damaged；If otherwise the measurement point is damage position.Therefore, can be under high frequency, high speed background condition using method provided by the invention or system, the microlesion of real-time detection structure transient state simultaneously determines damage position.
Description
Technical field
The present invention relates to diagnosing structural damage field, more particularly to a kind of structural damage on-line real time monitoring method and is
System.
Background technology
Online structural damage detection, diagnosis are an importances of diagnosing structural damage.The method proposed in this regard
It is most of to be all based on the change of model parameter (intrinsic frequency, the vibration shape, change ratio etc.) or utilize FEM updating
Technology etc..Damage diagnosis identification essence based on vibration signal is extraction and the analysis process of a Dynamic Signal statistical property,
Identify that the theoretical research related to damage check is many on Vibrating modal parameters at present.Early stage research includes:Yuen proposes root
The method for determining damage position according to the slope variation rate of the vibration shape, carrys out the damage of location structure；Aktan et al. researchs measure soft
The method for spending matrix judges the integrality of bridge structure；Han Xi et al. carries out damage knowledge using mode experiment data to arch structure
Not etc..In recent years, with the application of Time-Frequency Analysis Method, damage is gradually applied to using the methods of wavelet analysis and wavelet packet analysis
In wound identification.S.S.Patel etc. utilizes principle of wavelet analysis, has efficiently identified six layers of reinforced concrete structure degree of impairment；
Wavelet Packet Transform Method is successfully applied in jacket type offshore platform structure non-destructive tests by Asgarian etc.；Z.D.Zheng etc.
Using power spectral-density analysis method, degree of impairment of the plane framework under single damage and multiple damages is successfully identified.But
These research methods are mostly focused on low speed, low frequency, needed in the civil engineering class formation of long term monitoring, for industrial equipment
Structural damage on-line real time monitoring under central high speed, high frequency condition is relatively fewer.Such as high-speed slide electrical contact slide rail, should
High-speed slide electrical contact slide rail is critical component in the equipment such as automation equipment, robot and communication, is entirely to set
Standby or system key structure.Under the conditions of big current-carrying, high speed, high-frequency work, high-speed slide electrical contact slide rail surface is also easy to produce
Microlesion, such as：Ablation, abrasion, planing equivalent damage phenomenon, operating efficiency is had a strong impact on, and high-speed slide electrical contact slide rail exists
Inside mechanical system, in invisible area, damage check can not be carried out by above method.Therefore, being badly in need of one kind can
Under high frequency, high speed background condition, the method for real-time detection structure degree of impairment.
The content of the invention
, can be in high frequency, high speed it is an object of the invention to provide a kind of structural damage on-line real time monitoring method and system
Under background condition, ONLINE RECOGNITION and the positioning of structure transient state microlesion are realized.
To achieve the above object, the invention provides following scheme：
A kind of structural damage on-line real time monitoring method, the structural damage on-line real time monitoring method are applied to high ski-running
Dynamic electrical contact slide rail；The high-speed slide electrical contact slide rail includes high-speed slide electrical contact pressing plate, sliding block and slide rail, the cunning
Rail is arranged in the high-speed slide electrical contact pressing plate；The sliding block slides in the slide rail；The structural damage is real online
When monitoring method include：
Obtain in high-speed slide electrical contact each measurement point of pressing plate, direction is to make electrical contact with to press perpendicular to the high-speed slide
The acceleration signal of plate；
Obtain the intrinsic frequency of the high-speed slide electrical contact slide rail when slide rail surface is not damaged and the high speed
Slide the energy value in electrical contact each measurement point of pressing plate；
Integral Processing is carried out to the acceleration signal, obtains displacement signal；
Set Empirical Mode Decomposition is carried out to institute's displacement signal, obtains gathering the first natural frequency after Empirical Mode Decomposition
With second-order intrinsic frequency；
Judge the first natural frequency whether equal with the intrinsic frequency and the second-order intrinsic frequency and institute
Whether equal state intrinsic frequency, obtain the first judged result；
If it is equal with the intrinsic frequency and described second that first judged result represents the first natural frequency
Rank intrinsic frequency is equal with the intrinsic frequency, it is determined that the slide rail surface is not damaged；
If first judged result represents that the first natural frequency and the intrinsic frequency are unequal or described
Second-order intrinsic frequency and the intrinsic frequency are unequal, it is determined that the slide rail surface is damaged, and is calculated in different cunnings
The energy value of high-speed slide electrical contact each measurement point of pressing plate described in dynamic number；
Judge that high-speed slide makes electrical contact with the energy value of each measurement point and the slide rail in pressing plate described in different slip numbers
Whether the energy value of each measurement point is equal in high-speed slide electrical contact pressing plate when surface is not damaged, obtains the second judgement knot
Fruit；
If second judged result represents the different energy values for sliding number high speed and sliding electrical contact pressing plate measurement point
The energy value of high-speed slide electrical contact pressing plate measurement point when not damaged with the slide rail surface is all equal, it is determined that described
High-speed slide electrical contact pressing plate measurement point is not damaged；
If second judged result represents any energy value for sliding number high speed and sliding electrical contact pressing plate measurement point
The energy value of high-speed slide electrical contact pressing plate measurement point when not damaged with the slide rail surface is unequal, it is determined that described
High-speed slide electrical contact pressing plate measurement point is damaged.
Optionally, described to obtain in high-speed slide electrical contact each measurement point of pressing plate, direction is perpendicular to the high speed
The acceleration signal of electrical contact pressing plate is slided, is specifically included：
Multiple piezoelectric acceleration sensors are arranged in high-speed slide electrical contact clamp surface；
Obtain the acceleration signal of the piezoelectric acceleration sensor collection；The direction of the acceleration signal is perpendicular to institute
State high-speed slide electrical contact pressing plate.
Optionally, it is described that set Empirical Mode Decomposition is carried out to institute's displacement signal, obtain gathering the after Empirical Mode Decomposition
First natural frequency and second-order intrinsic frequency, are specifically included：
On the basis of institute's displacement signal, a white Gaussian noise signal is added, obtains overall signal；
Empirical mode decomposition is carried out to the overall signal, obtains the first rank intrinsic mode function and second-order natural mode of vibration
Function；The first rank intrinsic mode function is the oscillating function that amplitude and frequency change；The second-order intrinsic mode function
The oscillating function changed for amplitude and frequency；
According to the first rank intrinsic mode function and the second-order intrinsic mode function, first natural frequency is obtained
With second-order intrinsic frequency.
Optionally, the energy for calculating high-speed slide electrical contact each measurement point of pressing plate described in different slip numbers
Value, is specifically included：
Obtain in different slip numbers, the acceleration signal in high-speed slide electrical contact each measurement point of pressing plate；
Acceleration signal in different slip each measurement points of number is handled, calculated described in different slip numbers
High-speed slide makes electrical contact with the energy value of each measurement point of pressing plate.
Optionally, the acceleration signal in different slip each measurement points of number is handled, and is calculated in different cunnings
The energy value of high-speed slide electrical contact each measurement point of pressing plate described in dynamic number, is specifically included：
According to formula (1), the acceleration signal in different slip each measurement points of number is handled, calculated in different cunnings
The acceleration autopower spectral density of high-speed slide electrical contact each measurement point of pressing plate described in dynamic number；The formula (1) is
S (f) represents the acceleration autopower spectral density of f frequencies in formula (1)；M represents acceleration signal tract hop count；
NFFTRepresent the length of Fast Fourier Transform (FFT)；Using MATLAB by the fast Fourier transform of acceleration signal and forming frequency
Domain power spectral density plot；Xi(t) Fourier transformation of the acceleration signal in i-th of data segment of t is represented；For Xi
(t) conjugate complex number；
According to formula (2), the energy that the high-speed slide described in different slip numbers makes electrical contact with each measurement point of pressing plate is calculated
Value；
The formula (2) is
E is represented from f in formula (2)1Frequency is to f2The energy value of frequency.
Present invention also offers a kind of structural damage on-line real time monitoring system, the structural damage on-line real time monitoring system
System is applied to high-speed slide electrical contact slide rail；The high-speed slide electrical contact slide rail includes high-speed slide electrical contact pressing plate, sliding block
And slide rail, the slide rail are arranged in the high-speed slide electrical contact pressing plate；The sliding block slides in the slide rail；It is described
Structural damage on-line real time monitoring system includes：
Acceleration signal acquisition module, for obtaining high-speed slide electrical contact each measurement point of pressing plate, direction is vertical
Directly in the acceleration signal of high-speed slide electrical contact pressing plate；
Intrinsic frequency and energy value acquisition module, for obtaining the high-speed slide electricity when slide rail surface is not damaged
Contact the energy value in the intrinsic frequency and high-speed slide electrical contact each measurement point of pressing plate of slide rail；
Displacement information obtains module, for carrying out Integral Processing to the acceleration signal, obtains displacement signal；
First natural frequency and second-order intrinsic frequency obtain module, for carrying out set experience to institute's displacement signal
Mould decomposes, and obtains gathering the first natural frequency and second-order intrinsic frequency after Empirical Mode Decomposition；
First judged result obtains module, for judging whether the first natural frequency is equal with the intrinsic frequency
And whether the second-order intrinsic frequency is equal with the intrinsic frequency, the first judged result is obtained；
Slide rail surface does not occur to damage determining module, for representing that first rank is intrinsic when first judged result
Frequency is equal with the intrinsic frequency and when the second-order intrinsic frequency is equal with the intrinsic frequency, determines the slide rail table
Face is not damaged；
Damage determining module occurs for slide rail surface, for representing the first natural frequency when first judged result
With the intrinsic frequency is unequal or when the second-order intrinsic frequency and the unequal intrinsic frequency, determines the slide rail
Surface is damaged, and calculates the energy value of high-speed slide electrical contact each measurement point of pressing plate described in different slip numbers；
Second judged result obtains module, for judging that high-speed slide makes electrical contact with each in pressing plate described in different slip numbers
The energy of each measurement point in the high-speed slide electrical contact pressing plate when energy value of measurement point does not damage with the slide rail surface
Whether value is equal, obtains the second judged result；
High-speed slide electrical contact pressing plate measurement point does not occur to damage determining module, for being represented when second judged result
The height when energy value that difference slides number high speed slip electrical contact pressing plate measurement point does not damage with the slide rail surface
When the energy value of the dynamic electrical contact pressing plate measurement point of ski-running is all equal, determine that the high-speed slide electrical contact pressing plate measurement point does not occur
Damage；
Damage determining module occurs for high-speed slide electrical contact pressing plate measurement point, for representing to appoint when second judged result
The high speed when energy value that one slip number high speed slides electrical contact pressing plate measurement point does not damage with the slide rail surface
When the energy value of slip electrical contact pressing plate measurement point is unequal, determine that the high-speed slide electrical contact pressing plate measurement point is damaged
Wound.
Optionally, the first natural frequency and second-order intrinsic frequency obtain module, specifically include：
Overall signal obtains unit, on the basis of institute's displacement signal, adding a white Gaussian noise signal, obtaining
To overall signal；
First rank intrinsic mode function and second-order intrinsic mode function obtain unit, for being carried out to the overall signal
Empirical mode decomposition, obtain the first rank intrinsic mode function and second-order intrinsic mode function；The first rank natural mode of vibration letter
The oscillating function that number changes for amplitude and frequency；The second-order intrinsic mode function is the vibration letter that amplitude and frequency change
Number；
First natural frequency and second-order intrinsic frequency obtain unit, for according to the first rank intrinsic mode function
With the second-order intrinsic mode function, first natural frequency and second-order intrinsic frequency are obtained.
Optionally, damage determining module occurs for the slide rail surface, specifically includes：
Acceleration signal acquiring unit, for obtaining in different slip numbers, the high-speed slide electrical contact pressing plate is respectively surveyed
Acceleration signal on amount point；
Energy value computing unit is handled the acceleration signal in different slip each measurement points of number, is calculated in difference
Slide the energy value of each measurement point of pressing plate of high-speed slide electrical contact described in number.
Optionally, the energy value computing unit, is specifically included：
Acceleration autopower spectral density computation subunit, for according to formula (1), being slided to different in each measurement points of number
Acceleration signal handled, calculate it is different slip numbers described in high-speed slide electrical contact each measurement point of pressing plate acceleration
Spend autopower spectral density；The formula (1) is
S (f) represents the acceleration autopower spectral density of f frequencies in formula (1)；M represents acceleration signal tract hop count；
NFFTRepresent the length of Fast Fourier Transform (FFT)；Using MATLAB by the fast Fourier transform of acceleration signal and forming frequency
Domain power spectral density plot；Xi(t) Fourier transformation of the acceleration signal in i-th of data segment of t is represented；For Xi
(t) conjugate complex number；
Energy value computation subunit, for according to formula (2), calculating the high-speed slide electricity described in different slip numbers and connecing
Press the energy value of each measurement point of plate；
The formula (2) is
E is represented from f in formula (2)1Frequency is to f2The energy value of frequency.
According to specific embodiment provided by the invention, the invention discloses following technique effect：The present invention's provides one
Kind structural damage on-line real time monitoring method and system, this method or system are applied to high-speed slide electrical contact slide rail；It is described
High-speed slide electrical contact slide rail includes high-speed slide electrical contact pressing plate, sliding block and slide rail, the slide rail and is arranged on the high speed
Slide in electrical contact pressing plate；The sliding block slides in the slide rail.Wherein, this method or the system are first by real-time
Each measurement point acceleration signal that direction is obtained perpendicular to high-speed slide electrical contact pressing plate carries out Integral Processing, obtains high-speed slide
The microvibration response signal (displacement signal) of pressing plate is made electrical contact with, obtains the high-speed slide electricity when slide rail surface is not damaged
Contact the intrinsic frequency of slide rail and the energy value of high-speed slide electrical contact each measurement point of slide rail；Secondly to vibration response signal
Using the set Empirical Mode Decomposition algorithm that Modal Parameter Identification can be carried out to structure microvibration, first natural frequency is obtained
With second-order intrinsic frequency, and judge whether first natural frequency and second-order intrinsic frequency are equal with intrinsic frequency；If
Then determine that the slide rail surface is not damaged, slide rail surface is damaged if otherwise determining, and is calculated in different slip numbers
Described in high-speed slide electrical contact each measurement point of pressing plate energy value, judge it is different slip numbers described in high-speed slide make electrical contact with
Respectively measured in the high-speed slide electrical contact pressing plate when energy value of each measurement point does not damage with the slide rail surface in pressing plate
Whether the energy value of point is equal；If then determine that the high-speed slide electrical contact pressing plate measurement point is not damaged；It is if otherwise true
The fixed high-speed slide electrical contact pressing plate measurement point is damaged, and determines damage position.Therefore, using method provided by the invention
Or system, can be under high frequency, high speed background condition, detection high-speed slide electrical contact slide rail transient state microlesion in real time is simultaneously true
Determine damage position, similarly push away it, also to the microlesion of other structures transient state and damage position can be determined.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to institute in embodiment
The accompanying drawing needed to use is briefly described, it should be apparent that, drawings in the following description are only some implementations of the present invention
Example, for those of ordinary skill in the art, without having to pay creative labor, can also be according to these accompanying drawings
Obtain other accompanying drawings.
Fig. 1 is the schematic flow sheet of structural damage on-line real time monitoring method of the embodiment of the present invention；
Fig. 2 is that high-speed slide of the embodiment of the present invention makes electrical contact with sled test machine structural representation；
Fig. 3 is that acceleration signal of the embodiment of the present invention changes over time schematic diagram；
Fig. 4 is that displacement signal of the embodiment of the present invention changes over time schematic diagram；
Fig. 5 is acceleration autopower spectral density comparison diagram of the embodiment of the present invention；
Fig. 6 is acceleration spectrum energy value comparison diagram of the embodiment of the present invention；
Fig. 7 is that the embodiment of the present invention damages fore-aft acceleration autopower spectral density comparison diagram；
Fig. 8 is the structural representation of structural damage on-line real time monitoring system of the embodiment of the present invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other under the premise of creative work is not made
Embodiment, belong to the scope of protection of the invention.
, can be in high frequency, high speed it is an object of the invention to provide a kind of structural damage on-line real time monitoring method and system
Under background condition, ONLINE RECOGNITION and the positioning of structure transient state microlesion are realized.
In order to facilitate the understanding of the purposes, features and advantages of the present invention, it is below in conjunction with the accompanying drawings and specific real
Applying mode, the present invention is further detailed explanation.
Embodiment one
To achieve the above object, the embodiments of the invention provide a kind of structural damage on-line real time monitoring method, the structure
Damage on-line real time monitoring method and be applied to high-speed slide electrical contact slide rail；The high-speed slide electrical contact slide rail includes high ski-running
Dynamic electrical contact pressing plate, sliding block and slide rail, the slide rail are arranged in the high-speed slide electrical contact pressing plate；The sliding block is in institute
State in slide rail and slide.
Fig. 1 is the schematic flow sheet of structural damage on-line real time monitoring method of the embodiment of the present invention, as shown in figure 1, this hair
The structural damage on-line real time monitoring method that bright embodiment provides specifically includes following steps：
Step 101：Obtain in high-speed slide electrical contact each measurement point of pressing plate, direction is perpendicular to the high-speed slide
Make electrical contact with the acceleration signal of pressing plate.
Step 102：Obtain when the slide rail surface is not damaged the high-speed slide electrical contact slide rail intrinsic frequency and
Energy value in high-speed slide electrical contact each measurement point of pressing plate.
In the present embodiment, the intrinsic frequency and the high-speed slide of high-speed slide electrical contact slide rail when not damaging
When energy value in electrical contact each measurement point of pressing plate slides for first time sliding block in the slide rail, the intrinsic frequency of the first rank of acquisition
Energy value in rate, high-speed slide electrical contact pressing plate in each measurement point.
Step 103：Integral Processing is carried out to the acceleration signal, obtains displacement signal.
Step 104：Set Empirical Mode Decomposition is carried out to institute's displacement signal, obtains gathering the first rank after Empirical Mode Decomposition
Intrinsic frequency and second-order intrinsic frequency.
Step 105：Judge whether the first natural frequency is equal with the intrinsic frequency and the second-order is intrinsic
Whether frequency is equal with the intrinsic frequency, obtains the first judged result.
If it is equal with the intrinsic frequency and described second that first judged result represents the first natural frequency
Rank intrinsic frequency is equal with the intrinsic frequency, then performs step 106.
Step 106：Determine that the slide rail surface is not damaged.
If first judged result represents that the first natural frequency and the intrinsic frequency are unequal or described
Second-order intrinsic frequency and the intrinsic frequency are unequal, then step 107.
Step 107：Determine that the slide rail surface is damaged, and calculate the high-speed slide electricity described in different slip numbers
Contact the energy value of each measurement point of pressing plate.
Step 108：Judge it is different slip numbers described in high-speed slide electrical contact pressing plate in each measurement point energy value with
Whether the energy value of each measurement point is equal in the high-speed slide electrical contact pressing plate when slide rail surface is not damaged, obtains the
Two judged results.
If second judged result represents the different energy values for sliding number high speed and sliding electrical contact pressing plate measurement point
The energy value of high-speed slide electrical contact pressing plate measurement point when not damaged with the slide rail surface is all equal, then performs step
109。
Step 109：Determine that the high-speed slide electrical contact pressing plate measurement point is not damaged.
If second judged result represents any energy value for sliding number high speed and sliding electrical contact pressing plate measurement point
The energy value of high-speed slide electrical contact pressing plate measurement point when not damaged with the slide rail surface is unequal, then performs step
110。
In order to which clear explanation step 107 is to step 110, it is illustrated with.For example, the present embodiment has carried out 5 sliding blocks
Slided in slide rail, high-speed slide electrical contact pressing plate measurement point is three, is the first measurement point, the second measurement point and the 3rd measurement
Point.
This 5 sliding blocks are calculated respectively in slide rail slip, the energy value of the first measurement point, the energy value of the second measurement point,
The energy value of 3rd measurement point.
This 5 sliding blocks are judged in slide rail slip, when the energy value of the first measurement point does not damage with the slide rail surface
Whether the energy value of the first measurement point all equal in high-speed slide electrical contact pressing plate, i.e., each sliding block in slide rail slip,
First measurement point in the high-speed slide electrical contact pressing plate when energy value of first measurement point does not damage with the slide rail surface
Energy value it is whether all equal；
If, it is determined that the first measurement point is not damaged；
In slide rail slip, if the energy value of the first measurement point does not damage wherein any sliding block with the slide rail surface
When the high-speed slide electrical contact pressing plate in the first measurement point energy value it is unequal, it is determined that the first measurement point is damaged
Wound, and can determine which time sliding block first measurement point when slide rail slides is damaged.
Similarly, the second measurement point and the 3rd measurement point judge as described above.
Step 110：Determine that the high-speed slide electrical contact pressing plate measurement point is damaged.
Wherein, step 101 specifically includes：
Multiple piezoelectric acceleration sensors are arranged in high-speed slide electrical contact clamp surface.
Obtain the acceleration signal of the piezoelectric acceleration sensor collection；The direction of the acceleration signal is perpendicular to institute
State high-speed slide electrical contact pressing plate.
Step 103 specifically includes：
Integral Processing twice is carried out to the acceleration signal of acquisition, obtains displacement signal.
Institute's displacement signal is handled using set Empirical Mode Decomposition algorithm, obtains first natural frequency and second
Rank intrinsic frequency.
The set Empirical Mode Decomposition (Ensemble Empirical Mode Decomposition, abbreviation EEMD) is calculated
Method is improved on the basis of empirical mode decomposition (Empirical Mode Decomposition, abbreviation EMD) algorithm.
EMD algorithms utilize the local feature time scale of original signal, and some rank natural mode of vibration letters are extracted from original signal
Number (Intrinsic Mode Function, abbreviation IMF) and a residual volume, each rank IMF signals decomposited highlight data
Local feature, residue signal embodies the slowly varying amount in signal.Each IMF signals are analyzed, can be more accurate
Effectively hold the characteristic information of original signal.
Wherein, following two conditions must are fulfilled for per rank IMF signals：(1) number of the zero crossing in whole IMF signals with
The number of extreme point it is equal or at most differ 1；(2) envelope that is determined by the Local modulus maxima of IMF signals and by part
The average for the envelope that minimum point determines is zero, i.e., IMF signals are on time shaft Local Symmetric.
EMD decomposition is carried out to either signal x (t) to concretely comprise the following steps：
(1) all maximum points and minimum point on primary signal x (t), are determined.
(2), all maximum points and minimum point determined are connected respectively with Cubic Spline Fitting curve, and
Using this two curves as x (t) envelope up and down.
(3) the average value curve m of upper and lower envelope, is calculated1(t).M is subtracted with x (t)1(t)：
h1(t)=x (t)-m1(t) (3)；
If h1(t) two conditions of IMF signals in the above are unsatisfactory for, it is necessary to h1(t) repeated as primary signal
Step (1)-step (3), is obtained：
h11(t)=h1(t)-m11(t) (4)；
K times is so screened until h1k(t) meet two conditions of IMF signals in the above, be changed into an IMF signal,
Referred to as the first rank IMF signals c1(t), it is designated as：
c1(t)=h1k(t) (5)；
(4) c, is subtracted from primary signal1(t) the first rank residual signal r, is obtained1(t)。
In the first rank residual signal r1(t) more macrocyclic signal is still included in, therefore, it is necessary to r1(t) step is carried out
(1) to the screening of step (4), so obtain successively second-order IMF signals ..., the i-th rank IMF signals ..., n-th order IMF signals
With second-order residual signal ..., n-th order residual signal.As n-th order residual signal rn(t) when turning into a monotonic function, screening
Terminate.Primary signal x (t) is represented by a series of intrinsic mode function ciAnd a discrepance r (t)n(t) sum, i.e. formula
(6)
From formula (6), primary signal x (t) is decomposed into the sum of each rank IMF signals by EMD algorithms according to certain rules, because
This can consider that the energy of primary signal has been assigned on each rank IMF signals, therefore, can detect letter using EMD algorithms
Minor variations in number.Overcome in the prior art using EMD algorithms because the minor variations of signal are relative to primary signal
In energy very little, the influence unobvious to whole primary signal, so be not suitable for structure microlesion detection the problem of.
Although there is spurious signal and mould in minor variations that can be in detection signal using EMD algorithms, the EMD algorithms
State Aliasing Problem.In order to solve the problems, such as spurious signal that EMD algorithms occur and modal overlap, it is necessary to add white noise to original
In signal, the yardstick of some missings is supplemented, overcomes drawbacks described above.Reach more preferable decomposition result, this namely EEMD algorithm
Basic ideas.
EEMD main algorithms structure and EMD algorithms are essentially identical, EEMD algorithm specific steps：
The first step, by being superimposed white Gaussian noise signal ω (t) to primary signal x (t), obtain an overall signal：
X (t)=x (t)+ω (t) (7)；
Second step, using EMD algorithms, X (t) is decomposed, obtain each rank IMF signals；It should be noted that divide every time
Single order solid-state mode function is solved, white Gaussian noise signal ω (t) will be superimposed again in residual signal；Obtain IMF signals
ci(t) it is represented by：
N represents overall number in formula (8), that is, adds the number of white Gaussian noise；Cj,i(t) represent to add jth time white Gaussian
The the i-th rank solid-state mode function decomposed after noise.
Due to the zero-mean principle of white Gaussian noise frequency spectrum, eliminate white Gaussian noise and brought as time domain distribution reference structure
Influence, the number of added white Gaussian noise obeys statistical law in EEMD：
In formula (9)：N is that N represents overall number, that is, adds the number of white Gaussian noise；ε is the amplitude of white Gaussian noise；
εnDecompose for primary signal and finally obtained intrinsic mode function and add error between obtained signal.That is,
In the case that noise amplitude is certain, overall number is more, finally decomposes obtained result closer to actual value.For institute's plus noise
Amplitude, if amplitude is too small, signal to noise ratio is too high, noise will be unable to have influence on the selection of limit, and then loses supplement yardstick
Effect.Generally, εn0.2, N is taken to take 100.To ensure algorithm Fast Convergent and effective detection, ε should not obtain too small.
I.e. primary signal X (t) can be analyzed to：
R in formula (10)N,n(t) the n-th order residual signal decomposed after superposition n-th white Gaussian noise is represented.
Step 104 specifically includes：
According to formula (7), on the basis of institute's displacement signal, a white Gaussian noise signal is added, obtains overall letter
Number.
The overall signal is decomposed using EMD algorithms presented hereinbefore, obtain the first rank intrinsic mode function and
Second-order intrinsic mode function；The first rank intrinsic mode function is the oscillating function that amplitude and frequency change；Described second
Rank intrinsic mode function is the oscillating function that amplitude and frequency change.
According to the first rank intrinsic mode function and the second-order intrinsic mode function, first natural frequency is obtained
With second-order intrinsic frequency.
Step 107 specifically includes：
Acceleration autopower spectral density contains the abundant information of vibration, and it shows the gross energy of vibration signal in frequency domain
Interior distribution situation.Acceleration autopower spectral density is often used to determine the self-vibration characteristic of structure, in equipment fault monitoring, also
The reason for sign may break down with searching can be occurred come failure judgement according to the change of different periods auto-power spectrum.Therefore,
Can be according to the Power spectral density of figure method average period, i.e. formula (1) calculates in different slip each measurement points of number
Acceleration signal autopower spectral density；The formula (1) is
S (f) represents the acceleration autopower spectral density of f frequencies in formula (1)；M represents acceleration signal tract hop count；
NFFTRepresent the length of Fast Fourier Transform (FFT)；Using MATLAB by the fast Fourier transform of acceleration signal and forming frequency
Domain power spectral density plot；Xi(t) Fourier transformation of the acceleration signal in i-th of data segment of t is represented；For Xi
(t) conjugate complex number；
In order to reduce random fluctuation, M.S. Charles Bartletts propose figure method average period, i.e., first signal sequence are divided into some
Section, calculates its cyclic graph to every section, then takes the average valuation as power spectrum of each cyclic graph respectively.M is represented namely
Signal sequence segments, i.e. signal sequence have several fragments to be just averaged several times.
The a certain frequency range f in area, according to formula (2), it is each to calculate the high-speed slide electrical contact pressing plate described in different slip numbers
The energy value of measurement point；
The formula (2) is
E is represented from f in formula (2)1Frequency is to f2The energy value of frequency.
In addition, method provided in an embodiment of the present invention also includes：Calculate the front and rear energy value difference of structure same position damage
Absolute value.Specially：
After structural damage produces, each point acceleration autopower spectral density changes in structure, according to each point acceleration certainly
Energy value of the power spectral density in different frequency also changes, and judges whether damage is sent out by the change size of energy value
It is raw, and determine that position occurs for damage.
Use Ei(f)、Represent that the acceleration spectrum energy value of forward and backward f frequency ranges is damaged in i-th position respectively, then i-th
The absolute value that the acceleration spectrum energy value difference of forward and backward f frequency ranges is damaged in position is：
By certain frequency Range-partition acceleration auto-power spectrum energy variation amount into m sections, f is designated as1,f2,...fm, at n
It is expressed as on position in the absolute value matrix of the acceleration spectrum energy value difference of m frequency range：
In above-mentioned matrix each column element represent same frequency, diverse location acceleration spectrum energy value difference absolute value, by
The absolute value of acceleration spectrum energy value difference is represented by caused by damage：
∑ E={ ∑sfE1(f)∑fE2(f)...∑fEn(f)} (13)。
∑ in formula (13)fEn(f) sum in representation formula (12) per a line.
If structural damage degree is lighter, the absolute value of acceleration spectrum energy value difference will be smaller.Conversely, acceleration spectrum energy
The absolute value of difference in magnitude will be larger.
Therefore, method provided in an embodiment of the present invention is applied to high-speed slide electrical contact slide rail, and is implemented by the present invention
The technical scheme that example provides also illustrate that technical scheme provided by the invention can monitor high-speed slide electrical contact slide rail table in real time
Whether face is damaged, and also can determine that damage position.But method provided by the invention is not limited to the embodiment of the present invention, the present invention
The method of offer can be used for being operated under the conditions of big load electric current (about tens of kA), superelevation sliding speed (about hundreds of m/s)
The ONLINE RECOGNITION of structure transient state microlesion and positioning.
Embodiment two
In order to verify embodiment one provide method can monitoring of structures structure transient state microlesion in real time, below by
Inventive embodiments two are verified.
Fig. 2 is that high-speed slide of the embodiment of the present invention makes electrical contact with sled test machine structural representation；Fig. 3 is the embodiment of the present invention
Acceleration signal changes over time schematic diagram；Fig. 4 is that displacement signal of the embodiment of the present invention changes over time schematic diagram；Fig. 5 is this
Inventive embodiments acceleration autopower spectral density comparison diagram；Fig. 6 is acceleration spectrum energy value comparison diagram of the embodiment of the present invention；Fig. 7
Fore-aft acceleration autopower spectral density comparison diagram is damaged for the embodiment of the present invention.
As shown in Fig. 2 the device that embodiment two uses makes electrical contact with sled test machine for high-speed slide；The high-speed slide electricity
Contacting sled test machine includes high-speed slide electrical contact pressing plate, sliding block, slide rail, 3 piezoelectric acceleration sensor (i.e. 3 measurements
Point), snap-fit connector cable, charge amplifier, oscillograph, data collecting card and host computer.The slide rail is set
In high-speed slide electrical contact pressing plate；The sliding block slides in the slide rail.
When sliding block slides one time in slide rail, the piezoelectricity for being laid in three measurement points on high-speed slide electrical contact pressing plate adds
Velocity sensor gathers acceleration signal, obtains acceleration signal as shown in Figure 3.The direction of the acceleration signal perpendicular to
The high-speed slide makes electrical contact with pressing plate, and is transmitted this acceleration signal to electric charge by snap-fit connector cable and amplified
Device is handled, the acceleration signal after then being handled by oscilloscope display, passes through adding after data collecting card acquisition process
Rate signal, and the acceleration signal after processing is transmitted to host computer.Wherein, the structural damage that the embodiment of the present invention one provides
The application hardware device of on-line real time monitoring method is host computer.
Integral Processing twice is carried out to the acceleration signal of acquisition, obtains displacement signal as shown in Figure 4.
According to the step 104 of embodiment one, by 6 experiments (i.e. sliding block has slided 6 times in slide rail), obtain such as the institute of table 1
Show first natural frequency and second-order intrinsic frequency, from table 1 it can be seen that, with the increase of experiment number, the of structure
First natural frequency, second-order intrinsic frequency are gradually reduced, and the main reason for causing this phenomenon damages for slide rail surface
Afterwards, the inner parameter of structure is changed, i.e. slide rail surface is damaged.
The rank natural frequency value of table 1 first and second-order natural frequency value
According to the step 107 of embodiment one, the acceleration signal of 5 experiments, three measurement points before acquisition, using formula
(1), each acceleration signal is handled, obtains the acceleration autopower spectral density of each measurement point, draws as shown in Figure 5 add
Speed autopower spectral density curve.From figure 5 it can be seen that frequency corresponding at each experiment acceleration autopower spectral density peak value
Rate is 9kHz, is understood using peak picking method principle：Vertical eigentone at the point is in 9kHz or so, with table 1
First natural frequency is close.After certain is damaged on slide rail, with the increase of experiment number, damage position is shown
Damage aggravation phenomenon.Can be with the degree of injury of qualitative characterization's slide rail by the peak change of acceleration autopower spectral density curve.
According to the step 107 of embodiment one, the high-speed slide electricity described in different slip numbers is calculated using formula (2) and connect
The energy value of each measurement point of plate is pressed, and takes wherein measurement point energy value curve map in 5 experiments, is drawn as shown in Figure 6
Acceleration spectrum energy value comparison diagram.From Fig. 6 this it appears that：With the increase of experiment number, when structure partial damages
When hindering, it can reduce in the acceleration spectrum energy value of the measurement point, and as experiment number increases, what acceleration spectrum energy value declined
Speed is increasing.Therefore, it may determine that damage position using the changing value of acceleration spectrum energy value curve and qualitative evaluation damage
Hinder degree；Simultaneously as damage caused by slide rail surface abrasion quality in mg levels, therefore every time experiment guide rail damage change it is smaller,
Also illustrate that this method is sensitive to Small loss, be advantageous to contact slide rail early stage micro-damage detection.
The change of fore-aft acceleration power spectral density is damaged around high-speed slide electrical contact each measurement point of pressing plate such as Fig. 7 institutes
Show, when curve understands to have near measuring point damage generation from Fig. 7, the peak value of corresponding acceleration power spectral density function curve
Diminish, when damage is reached to a certain degree, frequency can also diminish corresponding to its peak value, thus can orientated damage generation position.At this
In width figure, summit it is high be damage before, peak value it is low be damage after.
Therefore, a kind of structural damage on-line real time monitoring method is provided using the embodiment of the present invention one, can be in height
Frequently, under high speed background condition, high-speed slide electrical contact slide rail transient state microlesion is detected in real time and determines damage position, is similarly pushed away
It, also to the microlesion of other structures transient state and can determine damage position.
Embodiment three
To reach above-mentioned purpose, present invention also offers small caused by the raceway surface during big current-carrying high-speed slide
Damage, the online non-destructive tests of electrical contact raceway surface and alignment system under a kind of multi- scenarios method background in transient process, i.e., originally
Invention provides a kind of structural damage on-line real time monitoring system, can realize that structure transient state is micro- under high frequency, high speed background condition
The ONLINE RECOGNITION of Small loss and positioning.
Fig. 8 is the structural representation of structural damage on-line real time monitoring system of the embodiment of the present invention, as shown in figure 8, described
Structural damage on-line real time monitoring system includes：
Acceleration signal acquisition module 801, for obtaining high-speed slide electrical contact each measurement point of pressing plate, direction is
Perpendicular to the acceleration signal of high-speed slide electrical contact pressing plate.
Intrinsic frequency and energy value acquisition module 802, for obtaining the high ski-running when slide rail surface is not damaged
Energy value in the intrinsic frequency and high-speed slide electrical contact each measurement point of pressing plate of dynamic electrical contact slide rail.
Displacement information obtains module 803, for carrying out Integral Processing to the acceleration signal, obtains displacement signal.
First natural frequency and second-order intrinsic frequency obtain module 804, for gathering institute's displacement signal
Empirical Mode Decomposition, obtain gathering the first natural frequency and second-order intrinsic frequency after Empirical Mode Decomposition.
First judged result obtains module 805, for judging whether are the first natural frequency and the intrinsic frequency
Whether equal and described second-order intrinsic frequency is equal with the intrinsic frequency, obtains the first judged result.
Slide rail surface does not occur to damage determining module 806, for representing first rank when first judged result
Intrinsic frequency is equal with the intrinsic frequency and when the second-order intrinsic frequency is equal with the intrinsic frequency, determines the cunning
Track surface does not damage.
Damage determining module 807 occurs for slide rail surface, for representing that first rank is intrinsic when first judged result
Frequency and the intrinsic frequency be unequal or when the second-order intrinsic frequency and the unequal intrinsic frequency, it is determined that described
Slide rail surface is damaged, and calculates the energy of high-speed slide electrical contact each measurement point of pressing plate described in different slip numbers
Value.
Second judged result obtains module 808, for judging high-speed slide electrical contact pressing plate described in different slip numbers
In each measurement point in the high-speed slide electrical contact pressing plate of energy value and the slide rail surface of each measurement point when not damaging
Whether energy value is equal, obtains the second judged result.
High-speed slide electrical contact pressing plate measurement point does not occur to damage determining module 809, for when second judged result
Represent the institute when energy value that different slip number high speeds slide electrical contact pressing plate measurement point does not damage with the slide rail surface
State high-speed slide electrical contact pressing plate measurement point energy value it is all equal when, determine high-speed slide electrical contact pressing plate measurement point not
Damage.
Damage determining module 810 occurs for high-speed slide electrical contact pressing plate measurement point, for when the second judged result table
Show described when the energy value that any slip number high speed slides electrical contact pressing plate measurement point does not damage with the slide rail surface
When the energy value of high-speed slide electrical contact pressing plate measurement point is unequal, determine that the high-speed slide electrical contact pressing plate measurement point occurs
Damage.
The first natural frequency and second-order intrinsic frequency obtain module 804, specifically include：
Overall signal obtains unit, on the basis of institute's displacement signal, adding a white Gaussian noise signal, obtaining
To overall signal.
First rank intrinsic mode function and second-order intrinsic mode function obtain unit, for being carried out to the overall signal
Empirical mode decomposition, obtain the first rank intrinsic mode function and second-order intrinsic mode function；The first rank natural mode of vibration letter
The oscillating function that number changes for amplitude and frequency；The second-order intrinsic mode function is the vibration letter that amplitude and frequency change
Number.
First natural frequency and second-order intrinsic frequency obtain unit, for according to the first rank intrinsic mode function
With the second-order intrinsic mode function, first natural frequency and second-order intrinsic frequency are obtained.
Damage determining module 807 occurs for the slide rail surface, specifically includes：
Acceleration signal acquiring unit, for obtaining in different slip numbers, the high-speed slide electrical contact pressing plate is respectively surveyed
Acceleration signal on amount point.
Energy value computing unit is handled the acceleration signal in different slip each measurement points of number, is calculated in difference
Slide the energy value of each measurement point of pressing plate of high-speed slide electrical contact described in number；Wherein, the energy value computing unit, specifically
Including：
Acceleration autopower spectral density computation subunit, for according to formula (1), being slided to different in each measurement points of number
Acceleration signal handled, calculate it is different slip numbers described in high-speed slide electrical contact each measurement point of pressing plate acceleration
Spend autopower spectral density；The formula (1) is
S (f) represents the acceleration autopower spectral density of f frequencies in formula (1)；M represents acceleration signal tract hop count；
NFFTRepresent the length of Fast Fourier Transform (FFT)；Using MATLAB by the fast Fourier transform of acceleration signal and forming frequency
Domain power spectral density plot；Xi(t) Fourier transformation of the acceleration signal in i-th of data segment of t is represented；For Xi
(t) conjugate complex number.
Energy value computation subunit, for according to formula (2), calculating the high-speed slide electricity described in different slip numbers and connecing
Press the energy value of each measurement point of plate.
The formula (2) is
E is represented from f in formula (2)1Frequency is to f2The energy value of frequency.
The embodiment of the present invention three provides a kind of structural damage on-line real time monitoring system, and the system can be in high frequency, height
Under fast background condition, high-speed slide electrical contact slide rail transient state microlesion is detected in real time and determines damage position, similarly pushes away it,
To the microlesion of other structures transient state and damage position can be determined.
Each embodiment is described by the way of progressive in this specification, what each embodiment stressed be and other
The difference of embodiment, between each embodiment identical similar portion mutually referring to.For system disclosed in embodiment
For, because it is corresponded to the method disclosed in Example, so description is fairly simple, related part is said referring to method part
It is bright.
Specific case used herein is set forth to the principle and embodiment of the present invention, and above example is said
It is bright to be only intended to help the method and its core concept for understanding the present invention；Meanwhile for those of ordinary skill in the art, foundation
The thought of the present invention, in specific embodiments and applications there will be changes.In summary, this specification content is not
It is interpreted as limitation of the present invention.
Claims (9)
- A kind of 1. structural damage on-line real time monitoring method, it is characterised in that the structural damage on-line real time monitoring method should Slide rail is made electrical contact with for high-speed slide；High-speed slide electrical contact slide rail include high-speed slide electrical contact pressing plate, sliding block and Slide rail, the slide rail are arranged in the high-speed slide electrical contact pressing plate；The sliding block slides in the slide rail；The structure Damage on-line real time monitoring method includes：Obtain in high-speed slide electrical contact each measurement point of pressing plate, direction is to make electrical contact with pressing plate perpendicular to the high-speed slide Acceleration signal；Obtain the intrinsic frequency of the high-speed slide electrical contact slide rail when slide rail surface is not damaged and the high-speed slide Make electrical contact with the energy value in each measurement point of pressing plate；Integral Processing is carried out to the acceleration signal, obtains displacement signal；Set Empirical Mode Decomposition is carried out to institute's displacement signal, obtains gathering the first natural frequency and the after Empirical Mode Decomposition Second order intrinsic frequency；Judge whether the first natural frequency equal with the intrinsic frequency and the second-order intrinsic frequency with it is described solid There is frequency whether equal, obtain the first judged result；If first judged result represents that the first natural frequency is equal with the intrinsic frequency and the second-order is consolidated There is frequency equal with the intrinsic frequency, it is determined that the slide rail surface is not damaged；If it is unequal or described second that first judged result represents the first natural frequency and the intrinsic frequency Rank intrinsic frequency and the intrinsic frequency are unequal, it is determined that the slide rail surface is damaged, and is calculated in different slips time High-speed slide described in number makes electrical contact with the energy value of each measurement point of pressing plate；Judge that high-speed slide makes electrical contact with the energy value of each measurement point and the slide rail surface in pressing plate described in different slip numbers Whether the energy value of each measurement point is equal in high-speed slide electrical contact pressing plate when not damaging, obtains the second judged result；If second judged result represents different energy value and the institutes slided number high speed and slide electrical contact pressing plate measurement point The energy value for stating high-speed slide electrical contact pressing plate measurement point when slide rail surface is not damaged is all equal, it is determined that the high speed Electrical contact pressing plate measurement point is slided not damage；If second judged result represents any energy value and institute slided number high speed and slide electrical contact pressing plate measurement point The energy value for stating high-speed slide electrical contact pressing plate measurement point when slide rail surface is not damaged is unequal, it is determined that the high speed Electrical contact pressing plate measurement point is slided to damage.
- 2. structural damage on-line real time monitoring method according to claim 1, it is characterised in that described to obtain the high speed Slide in electrical contact each measurement point of pressing plate, direction is the acceleration signal that pressing plate is made electrical contact with perpendicular to the high-speed slide, specifically Including：Multiple piezoelectric acceleration sensors are arranged in high-speed slide electrical contact clamp surface；Obtain the acceleration signal of the piezoelectric acceleration sensor collection；The direction of the acceleration signal is perpendicular to the height The dynamic electrical contact pressing plate of ski-running.
- 3. structural damage on-line real time monitoring method according to claim 1, it is characterised in that described to believe the displacement Number enter row set Empirical Mode Decomposition, obtain gathering the first natural frequency and second-order intrinsic frequency after Empirical Mode Decomposition, tool Body includes：On the basis of institute's displacement signal, a white Gaussian noise signal is added, obtains overall signal；Empirical mode decomposition is carried out to the overall signal, obtains the first rank intrinsic mode function and second-order natural mode of vibration letter Number；The first rank intrinsic mode function is the oscillating function that amplitude and frequency change；The second-order intrinsic mode function is Amplitude and the oscillating function of frequency change；According to the first rank intrinsic mode function and the second-order intrinsic mode function, first natural frequency and are obtained Second order intrinsic frequency.
- 4. structural damage on-line real time monitoring method according to claim 1, it is characterised in that the calculating is in different cunnings The energy value of high-speed slide electrical contact each measurement point of pressing plate described in dynamic number, is specifically included：Obtain in different slip numbers, the acceleration signal in high-speed slide electrical contact each measurement point of pressing plate；Acceleration signal in different slip each measurement points of number is handled, calculated described in different slip numbers at a high speed Slide the energy value of electrical contact each measurement point of pressing plate.
- 5. structural damage on-line real time monitoring method according to claim 4, it is characterised in that described to different slips time Acceleration signal in each measurement point of number is handled, and it is each to calculate the high-speed slide electrical contact pressing plate described in different slip numbers The energy value of measurement point, is specifically included：According to formula (1), the acceleration signal in different slip each measurement points of number is handled, calculated in different slips time High-speed slide described in number makes electrical contact with the acceleration autopower spectral density of each measurement point of pressing plate；The formula (1) isS (f) represents the acceleration autopower spectral density of f frequencies in formula (1)；M represents acceleration signal tract hop count；NFFTTable Show the length of Fast Fourier Transform (FFT)；Using MATLAB by the fast Fourier transform of acceleration signal and forming frequency domain power Spectrum density curve；Xi(t) Fourier transformation of the acceleration signal in i-th of data segment of t is represented；For Xi(t) be total to Yoke plural number；According to formula (2), the energy value that the high-speed slide described in different slip numbers makes electrical contact with each measurement point of pressing plate is calculated；The formula (2) isE is represented from f in formula (2)1Frequency is to f2The energy value of frequency.
- 6. a kind of structural damage on-line real time monitoring system, it is characterised in that the structural damage on-line real time monitoring system should Slide rail is made electrical contact with for high-speed slide；High-speed slide electrical contact slide rail include high-speed slide electrical contact pressing plate, sliding block and Slide rail, the slide rail are arranged in the high-speed slide electrical contact pressing plate；The sliding block slides in the slide rail；The structure Damage on-line real time monitoring system includes：Acceleration signal acquisition module, for obtaining high-speed slide electrical contact each measurement point of pressing plate, direction be perpendicular to The acceleration signal of the high-speed slide electrical contact pressing plate；Intrinsic frequency and energy value acquisition module, for obtaining the high-speed slide electrical contact when slide rail surface is not damaged Energy value in the intrinsic frequency of slide rail and high-speed slide electrical contact each measurement point of pressing plate；Displacement information obtains module, for carrying out Integral Processing to the acceleration signal, obtains displacement signal；First natural frequency and second-order intrinsic frequency obtain module, for carrying out gathering Empirical Mode point to institute's displacement signal Solution, obtain gathering the first natural frequency and second-order intrinsic frequency after Empirical Mode Decomposition；First judged result obtains module, for judging whether the first natural frequency and the intrinsic frequency equal and institute It is whether equal with the intrinsic frequency to state second-order intrinsic frequency, obtains the first judged result；Slide rail surface does not occur to damage determining module, for representing the first natural frequency when first judged result When equal with the intrinsic frequency and described second-order intrinsic frequency is equal with the intrinsic frequency, the slide rail surface is determined not Damage；Damage determining module occurs for slide rail surface, for representing the first natural frequency and institute when first judged result When stating unequal intrinsic frequency or the second-order intrinsic frequency and the unequal intrinsic frequency, the slide rail surface is determined Damage, and calculate the energy value of high-speed slide electrical contact each measurement point of pressing plate described in different slip numbers；Second judged result obtains module, for judging respectively to measure in high-speed slide electrical contact pressing plate described in different slip numbers The energy value of each measurement point is in the high-speed slide electrical contact pressing plate when energy value of point does not damage with the slide rail surface It is no equal, obtain the second judged result；High-speed slide electrical contact pressing plate measurement point does not occur to damage determining module, for representing different when second judged result The high ski-running when energy value that slip number high speed slides electrical contact pressing plate measurement point does not damage with the slide rail surface When the energy value of dynamic electrical contact pressing plate measurement point is all equal, determine that the high-speed slide electrical contact pressing plate measurement point is not damaged Wound；Damage determining module occurs for high-speed slide electrical contact pressing plate measurement point, for representing any cunning when second judged result The high-speed slide when energy value that dynamic number high speed slides electrical contact pressing plate measurement point does not damage with the slide rail surface When the energy value of electrical contact pressing plate measurement point is unequal, determine that the high-speed slide electrical contact pressing plate measurement point is damaged.
- 7. structural damage on-line real time monitoring system according to claim 6, it is characterised in that the intrinsic frequency of the first rank Rate and second-order intrinsic frequency obtain module, specifically include：Overall signal obtains unit, on the basis of institute's displacement signal, adding a white Gaussian noise signal, obtains total Body signal；First rank intrinsic mode function and second-order intrinsic mode function obtain unit, for carrying out experience to the overall signal Mode decomposition, obtain the first rank intrinsic mode function and second-order intrinsic mode function；The first rank intrinsic mode function is Amplitude and the oscillating function of frequency change；The second-order intrinsic mode function is the oscillating function that amplitude and frequency change；First natural frequency and second-order intrinsic frequency obtain unit, for according to the first rank intrinsic mode function and institute Second-order intrinsic mode function is stated, obtains first natural frequency and second-order intrinsic frequency.
- 8. structural damage on-line real time monitoring system according to claim 6, it is characterised in that the slide rail surface occurs Determining module is damaged, is specifically included：Acceleration signal acquiring unit, for obtaining in different slip numbers, the high-speed slide makes electrical contact with each measurement point of pressing plate On acceleration signal；Energy value computing unit is handled the acceleration signal in different slip each measurement points of number, is calculated in different slips High-speed slide described in number makes electrical contact with the energy value of each measurement point of pressing plate.
- 9. structural damage on-line real time monitoring system according to claim 8, it is characterised in that the energy value calculates single Member, specifically include：Acceleration autopower spectral density computation subunit, for according to formula (1), adding in each measurement points of number to be slided to different Rate signal is handled, and calculates the acceleration of high-speed slide electrical contact each measurement point of pressing plate described in different slip numbers certainly Power spectral density；The formula (1) isS (f) represents the acceleration autopower spectral density of f frequencies in formula (1)；M represents acceleration signal tract hop count；NFFTTable Show the length of Fast Fourier Transform (FFT)；Using MATLAB by the fast Fourier transform of acceleration signal and forming frequency domain power Spectrum density curve；Xi(t) Fourier transformation of the acceleration signal in i-th of data segment of t is represented；For Xi(t) be total to Yoke plural number；Energy value computation subunit, for according to formula (2), calculating the high-speed slide electrical contact pressure described in different slip numbers The energy value of each measurement point of plate；The formula (2) isE is represented from f in formula (2)1Frequency is to f2The energy value of frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710818527.9A CN107505396A (en) | 2017-09-12 | 2017-09-12 | A kind of structural damage on-line real time monitoring method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710818527.9A CN107505396A (en) | 2017-09-12 | 2017-09-12 | A kind of structural damage on-line real time monitoring method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107505396A true CN107505396A (en) | 2017-12-22 |
Family
ID=60696409
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710818527.9A Pending CN107505396A (en) | 2017-09-12 | 2017-09-12 | A kind of structural damage on-line real time monitoring method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107505396A (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109472846A (en) * | 2018-12-27 | 2019-03-15 | 燕山大学 | The method for obtaining Bode diagram with MATLAB processing frequency sweep data |
| CN109521028A (en) * | 2018-12-04 | 2019-03-26 | 燕山大学 | A kind of 3-dimensional multi-layered lattice structure inside gross imperfection detection method of metal |
| CN109855770A (en) * | 2018-12-31 | 2019-06-07 | 华东交通大学 | A method of TEMPERATURE FORCE is detected based on rail transverse acceleration power spectral density |
| CN109946389A (en) * | 2019-01-31 | 2019-06-28 | 青岛理工大学 | Damage Detection of Structures based on overall experience mode decomposition and convolutional neural networks |
| CN110243465A (en) * | 2019-07-19 | 2019-09-17 | 上海市建筑科学研究院 | Bridge vibration acceleration and intrinsic frequency on line real-time monitoring device, terminal and method |
| CN110308205A (en) * | 2019-08-08 | 2019-10-08 | 陕西师范大学 | Sorting metals method based on EMD generalized phase arrangement plan method |
| CN110371162A (en) * | 2019-06-19 | 2019-10-25 | 石家庄铁道大学 | A kind of plate track of high speed railway hurt detection method based on strain mode |
| CN110657906A (en) * | 2019-10-23 | 2020-01-07 | 南京航空航天大学 | Impact monitoring method based on fiber bragg grating sensor |
| CN111323488A (en) * | 2020-04-15 | 2020-06-23 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Damage detection method, device, equipment and medium for large metal roof panel |
| CN112345184A (en) * | 2020-09-28 | 2021-02-09 | 同济大学 | Structural earthquake damage identification method based on real-time monitoring data |
| CN113702500A (en) * | 2021-09-03 | 2021-11-26 | 郑州轻工业大学 | Sandwich composite material nondestructive testing method based on two-dimensional continuous wavelet transform |
-
2017
- 2017-09-12 CN CN201710818527.9A patent/CN107505396A/en active Pending
Non-Patent Citations (2)
| Title |
|---|
| 周航: "电接触构件微振动的检测与损伤识别研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 徐佳 等: "《希尔伯特——黄变换理论及其在重大工程变形监测中的应用》", 《希尔伯特——黄变换理论及其在重大工程变形监测中的应用》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109521028B (en) * | 2018-12-04 | 2021-06-25 | 燕山大学 | Method for detecting macroscopic defects inside metal three-dimensional multilayer lattice structure |
| CN109521028A (en) * | 2018-12-04 | 2019-03-26 | 燕山大学 | A kind of 3-dimensional multi-layered lattice structure inside gross imperfection detection method of metal |
| CN109472846B (en) * | 2018-12-27 | 2022-11-29 | 燕山大学 | Method for obtaining bode diagram by processing sweep frequency data by MATLAB |
| CN109472846A (en) * | 2018-12-27 | 2019-03-15 | 燕山大学 | The method for obtaining Bode diagram with MATLAB processing frequency sweep data |
| CN109855770A (en) * | 2018-12-31 | 2019-06-07 | 华东交通大学 | A method of TEMPERATURE FORCE is detected based on rail transverse acceleration power spectral density |
| CN109946389A (en) * | 2019-01-31 | 2019-06-28 | 青岛理工大学 | Damage Detection of Structures based on overall experience mode decomposition and convolutional neural networks |
| CN110371162A (en) * | 2019-06-19 | 2019-10-25 | 石家庄铁道大学 | A kind of plate track of high speed railway hurt detection method based on strain mode |
| CN110243465A (en) * | 2019-07-19 | 2019-09-17 | 上海市建筑科学研究院 | Bridge vibration acceleration and intrinsic frequency on line real-time monitoring device, terminal and method |
| CN110308205B (en) * | 2019-08-08 | 2021-10-08 | 陕西师范大学 | Metal classification method based on EMD generalized phase permutation entropy algorithm |
| CN110308205A (en) * | 2019-08-08 | 2019-10-08 | 陕西师范大学 | Sorting metals method based on EMD generalized phase arrangement plan method |
| CN110657906B (en) * | 2019-10-23 | 2020-06-09 | 南京航空航天大学 | Impact monitoring method based on fiber bragg grating sensor |
| CN110657906A (en) * | 2019-10-23 | 2020-01-07 | 南京航空航天大学 | Impact monitoring method based on fiber bragg grating sensor |
| CN111323488A (en) * | 2020-04-15 | 2020-06-23 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Damage detection method, device, equipment and medium for large metal roof panel |
| CN111323488B (en) * | 2020-04-15 | 2022-12-13 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Damage detection method, device, equipment and medium for large metal roof panel |
| CN112345184A (en) * | 2020-09-28 | 2021-02-09 | 同济大学 | Structural earthquake damage identification method based on real-time monitoring data |
| CN113702500A (en) * | 2021-09-03 | 2021-11-26 | 郑州轻工业大学 | Sandwich composite material nondestructive testing method based on two-dimensional continuous wavelet transform |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107505396A (en) | A kind of structural damage on-line real time monitoring method and system | |
| CN103575523B (en) | The rotary machinery fault diagnosis method of kurtosis-envelope spectrum analysis is composed based on FastICA- | |
| CN104198144B (en) | Middle and small bridge fast detecting method based on long-scale-distance optical fiber strain sensor | |
| CN107014668A (en) | A kind of fatigue crack integrated monitoring based on piezoelectricity and smart coat sensor | |
| CN102175768A (en) | Method and device for detecting defects and failures of high-speed rail based on vibration signals | |
| CN103076394B (en) | Safety evaluation method for ocean platform based on integration of vibration identification frequencies and vibration mode | |
| CN102422154A (en) | System, device for structural damage detection and method for structural damage detection | |
| CN104265577B (en) | Wind generating set abnormal detection method based on compressed sensing | |
| CN108519444A (en) | A kind of accurate measurement method of contact line defective locations | |
| CN111024527B (en) | Crack propagation monitoring method based on multi-sensor data fusion | |
| CN103217248A (en) | Method for detecting tensile force of bridge steel cable | |
| CN105862935A (en) | Damage recognition method used for retaining wall structural system | |
| CN112781820B (en) | Hob performance degradation trend evaluation method | |
| CN104462862A (en) | Mechanical structure dynamic load identification method based on cubic b-spline scaling function | |
| CN104537251A (en) | Fan blade impulse load recognition method | |
| Moeinifard et al. | Lost vibration test data recovery using convolutional neural network: a case study | |
| CN107454925A (en) | For the method for the remaining life for determining wind energy plant | |
| CN111678699A (en) | Early fault monitoring and diagnosing method and system for rolling bearing | |
| CN107036808A (en) | Gearbox of wind turbine combined failure diagnostic method based on SVMs probability Estimation | |
| CN104101542A (en) | Beam structure damage monitoring device and monitoring method thereof | |
| CN110057918B (en) | Method and system for quantitatively identifying damage of composite material under strong noise background | |
| CN106932162A (en) | Track dynamic stiffness method of testing and system | |
| CN116698317A (en) | Bridge modal parameter evaluation system based on array mobile sensor | |
| CN105600682A (en) | Crane strain monitoring method based on time-to-digital conversion technology | |
| CN112305378B (en) | Photoelectric composite submarine cable fault location system and fault location 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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171222 |