CN104792364A - Dynamic bridge parameter extracting system and dynamic bridge parameter extracting method based on laser Doppler - Google Patents

Dynamic bridge parameter extracting system and dynamic bridge parameter extracting method based on laser Doppler Download PDF

Info

Publication number
CN104792364A
CN104792364A CN201510173926.5A CN201510173926A CN104792364A CN 104792364 A CN104792364 A CN 104792364A CN 201510173926 A CN201510173926 A CN 201510173926A CN 104792364 A CN104792364 A CN 104792364A
Authority
CN
China
Prior art keywords
frequency
photosignal
sequence
vibration
bridge
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510173926.5A
Other languages
Chinese (zh)
Other versions
CN104792364B (en
Inventor
王翔
汪正兴
汪双炎
王波
柴小鹏
伊建军
马长飞
荆国强
刘鹏飞
王艳芬
高阳
张华兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Bridge Science Research Institute Ltd of MBEC
China Railway Major Bridge Engineering Group Co Ltd MBEC
Original Assignee
China Railway Major Bridge Engineering Group Co Ltd MBEC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Railway Major Bridge Engineering Group Co Ltd MBEC filed Critical China Railway Major Bridge Engineering Group Co Ltd MBEC
Priority to CN201510173926.5A priority Critical patent/CN104792364B/en
Publication of CN104792364A publication Critical patent/CN104792364A/en
Application granted granted Critical
Publication of CN104792364B publication Critical patent/CN104792364B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses a dynamic bridge parameter extracting system and a dynamic bridge parameter extracting method based on laser Doppler and relates to the technical field of bridge detection. The dynamic bridge parameter extracting method includes the steps of S1, collecting photoelectric signals of a photoelectric detector by a high-frequency signal collector, and converting the photoelectric signals into voltage values to obtain a photoelectric signal discrete sequence ; S2, by the aid of a computer, subjecting the photoelectric signal discrete sequence to fast Fourier transform so as to obtain a frequency movement sequence, seeking the specific position of the photoelectric signal corresponding to a peak value of the frequency movement sequence, subjecting a frequency-subdivided interval to further subdivision to compute a frequency spectrum within the frequency-subdivided interval, seeking the specific position of the photoelectric signal corresponding to a peak value of a frequency spectrum sequence, and computing a Doppler frequency shift; S3, computing moving speed of an object; S4, judging whether vibrating frequency of the object needs to be computed or not; S5, computing the vibrating frequency of the object according to a speed sequence. The dynamic bridge parameter extracting system and the dynamic bridge parameter extracting method based on the laser Doppler are high in extracting precision and used for measuring the low moving speed and the vibrating frequency of an engineering structure like a bridge.

Description

The dynamic parameter extraction system and method for bridge based on laser-Doppler
Technical field
The present invention relates to Bridge Inspection field, be specifically related to the dynamic parameter extraction system and method for a kind of bridge based on laser-Doppler.
Background technology
In recent years, the sensor technology for bridge detection of dynamic is flourish, comprises the acceleration transducer of the principles such as based semiconductor material piezoresistive effect, material piezoelectric effect, but is substantially all touch sensor.Touch sensor drawback is very large, and installing touch sensor wastes time and energy, and can cause damage in various degree to bridge structure, can hinder site traffic after installation.Utilize contactless bridge dynamic detection technology can solve above-mentioned drawback, wherein, the knotmeter based on laser Doppler principle is a Typical Representative of contactless technique of dynamic measurement.Knotmeter based on laser Doppler principle has the excellent properties such as rate accuracy is high, dynamic range large, response of testing the speed is fast, be widely used in measuring atomic flowing velocity in fluid, and the flowing velocity etc. of various types of cells in measurement biological blood, be also applied to the dynamic parameter measuring bridge.
Shown in Figure 1, a kind of bridge vibration detection device based on laser, comprise interconnective test box 1 and computing machine 14, the inside of described test box 1 is provided with optical-fiber laser transmitter 5, first fiber coupler 4, first fiber optic collimator mirror 2, second fiber optic collimator mirror 9, second fiber coupler 11, photodetector 12, high-frequency signal Acquisition Instrument 13 and laser pen 10; Described laser pen 10 is between the first fiber optic collimator mirror 2 and the second fiber optic collimator mirror 9; Described optical-fiber laser transmitter 5, first fiber coupler 4 is connected in turn with the first fiber optic collimator mirror 2, and described second fiber optic collimator mirror 9, second fiber coupler 11, photodetector 12 are connected in turn with high-frequency signal Acquisition Instrument 13; Described first fiber coupler 4 is connected with the second fiber coupler 11, and described high-frequency signal Acquisition Instrument 13 is connected with the computing machine 14 being provided with laser Doppler signal disposal system.The principle of work of this pick-up unit: optical-fiber laser transmitter 5 Emission Lasers, laser is by being divided into the first laser and the second laser after the first fiber coupler 4, first laser can change frequency because of Doppler effect through the body surface of vibration, therefore, the frequency diffused and the first laser, the frequency of the second laser is different, vibration velocity due to incident light and radiative difference on the frequency and object has a fixing analytic relationship, laser Doppler signal disposal system is provided with in computing machine 14, can by diffusing and the difference on the frequency of the first laser, calculate the vibration information of determinand.
Different target movement velocity scope to the requirement very different of laser Doppler signal disposal system, because bridge structure movement velocity is less, therefore flashlight and reference light difference on the frequency very little, signal to noise ratio (S/N ratio) decline, the extraction accuracy of light frequency is lower.For the engineering structure of the low-speed motions such as bridge, current laser Doppler signal disposal system is difficult to the dynamic parameter in measuring-signal to extract, and therefore, is difficult to the low-speed motion measuring bridge structure, such as speed is 100m/s ~ 1mm/s, even the measurement of 100m/s ~ 0.1mm/s.
Summary of the invention
For the defect existed in prior art, the object of the present invention is to provide the dynamic parameter extraction system and method for a kind of bridge based on laser-Doppler, the extraction accuracy of light frequency is high, for measuring low-speed motion speed and the vibration frequency of the engineering structures such as bridge.
For reaching above object, the technical scheme that the present invention takes is: the dynamic parameter extraction system of a kind of bridge based on laser-Doppler, and this system comprises:
Photodetector, for detecting photosignal;
High-frequency signal Acquisition Instrument, for gathering the photosignal of certain length in photodetector, and being converted into magnitude of voltage, obtaining photosignal discrete series;
Computing machine, for carrying out Fast Fourier Transform (FFT) to photosignal discrete series, obtains frequency movement sequence; Find the particular location of the photosignal corresponding to frequency movement sequence peaks, obtain frequency fine by stages; Segmented further frequency fine by stages, the frequency spectrum on calculated rate subdivided interval, obtains spectrum sequence; Find the particular location of the photosignal corresponding to peak value of spectrum sequence, calculate Doppler shift; According to Doppler shift, calculate target speed; Judge whether the vibration frequency needing to calculate target, judge whether the vibration frequency needing to calculate target, if not, then calculate movement velocity corresponding to different length photosignal; If so, then formed according to movement velocity and obtain the vibration frequency that velocity series calculates target; Fast Fourier Transform (FFT) is carried out to velocity series, obtains frequency vibration sequence; Find the particular location of the photosignal corresponding to frequency vibration sequence peaks, obtain vibration frequency subdivided interval; Frequency fine by stages is segmented further, calculates the rumble spectrum on vibration frequency subdivided interval, obtain rumble spectrum sequence, find the particular location of the photosignal corresponding to rumble spectrum sequence peaks, calculate the vibration frequency of target;
The dynamic parameter extracting method of the bridge based on laser-Doppler for described system, comprises the following steps:
S1. high-frequency signal Acquisition Instrument gathers the photosignal of certain length in photodetector, and is converted into magnitude of voltage, obtains photosignal discrete series;
S2. computing machine carries out Fast Fourier Transform (FFT) to photosignal discrete series, obtains frequency movement sequence; Find the particular location of the photosignal corresponding to frequency movement sequence peaks, obtain frequency fine by stages; Segmented further frequency fine by stages, the frequency spectrum on calculated rate subdivided interval, obtains spectrum sequence; Find the particular location of the photosignal corresponding to peak value of spectrum sequence, calculate Doppler shift; According to Doppler shift, calculate target speed;
S3. judge whether the vibration frequency needing to calculate target, if not, then repeat step S1 ~ S2, calculate the movement velocity that different length photosignal is corresponding; If so, then step S4 is forwarded to;
S4. corresponding according to different length photosignal movement velocity, obtains velocity series;
S5. Fast Fourier Transform (FFT) is carried out to velocity series, obtain frequency vibration sequence; Find the particular location of the photosignal corresponding to frequency vibration sequence peaks, obtain vibration frequency subdivided interval; Frequency fine by stages is segmented further, calculates the rumble spectrum on vibration frequency subdivided interval, obtain rumble spectrum sequence, find the particular location of the photosignal corresponding to rumble spectrum sequence peaks, calculate the vibration frequency of target.
On the basis of technique scheme, the detailed process of step S1 is: in a period of time △ t, high-frequency signal Acquisition Instrument is according to voltage sample rate FS, gather the photosignal in photodetector, and being converted into magnitude of voltage, the length of described photosignal is N, N is positive integer, photosignal is at the corresponding magnitude of voltage x in each position of length direction, obtain photosignal discrete series X={x (n), n=0,1 ... N-1}, the resolution of frequency is Δ f
On the basis of technique scheme, step S2 specifically comprises the following steps:
S201. computing machine is to photosignal discrete series X={x (n), n=0, and 1 ... N-1} carries out Fast Fourier Transform (FFT), obtain frequency movement sequence Z=T (X)={ z (n), n=0,1 ... N-1}, wherein, T is Fourier transform operator;
S202. photosignal discrete series X={x (n) that collects at every batch of frequency movement sequence Z=T (X), n=0,1 ... simple spike value is there is in the scope of N-1}, find the photosignal particular location k in the longitudinal direction corresponding to frequency movement sequence Z=T (X) peak value, k ∈ [1, N], obtaining frequency fine by stages is
S203. according to accuracy class M 1, by frequency fine by stages be further subdivided into 2M 1equal portions, M 1for positive integer, further frequency analysis be spaced apart △, adopt the frequency spectrum on following formula calculated rate subdivided interval, obtain spectrum sequence Y (f):
Y ( f ) = Σ n = 0 N - 1 x ( n ) e - j 2 πfn ,
Wherein, Doppler shift e is natural constant, and j is imaginary unit;
S204. there is simple spike value in spectrum sequence Y (f) within the scope of f, finds the position p in the longitudinal direction of the photosignal corresponding to peak value of spectrum sequence Y (f), p ∈ [1,2M 1], adopt following formula to calculate photosignal discrete series X={x (n), n=0,1 ... the Doppler shift f of N-1}:
f = FS N × ( k - 1 ) + FS N × 2 M 1 p ;
S205. according to following formula, target speed v is calculated:
v = 1 2 f λ 0 ,
Wherein, λ 0for the wavelength of transmitted light of optical-fiber laser transmitter.
On the basis of technique scheme, in step S3, if not, then repeat step S1 ~ S2, computational length is the movement velocity that the photosignal of l is corresponding respectively, wherein, l=0,1 ..., L-1, L be positive integer.
On the basis of technique scheme, the detailed process of step S4 is: be the movement speed v that the photosignal of l is corresponding according to length, forms velocity series V={v (l), l=0,1 ... L-1}, wherein, the sampling rate of speed is Fs, the resolution of speed in frequency is new △ f, newly Δf = Fs L .
On the basis of technique scheme, step S5 specifically comprises the following steps:
S501. to velocity series V={v (l), l=0,1 ... L-1} carries out Fast Fourier Transform (FFT), obtain frequency vibration sequence W=T (V)=w (l), l=0,1...L-1}, wherein T is Fourier transform operator;
S502. frequency vibration sequence W=T (V) is at velocity series V={v (l), l=0,1 ... there is simple spike value in the scope of L-1}, find the particular location k of the photosignal corresponding to frequency vibration sequence W=T (V) peak value at length direction i, k i∈ [1, L], obtains vibration frequency subdivided interval: wherein, m is vibration exponent number, and m is positive integer;
S503. according to accuracy class M 2, by vibration frequency subdivided interval be further subdivided into 2M 2equal portions, M 2for positive integer, further vibration Frequency Analysis be spaced apart Δ 2, adopt the frequency spectrum on following formula calculating vibration frequency subdivided interval, obtain rumble spectrum sequence U (f i):
U ( f i ) = Σ l = 0 L - 1 v ( l ) e - j 2 π f i l ,
Wherein e is natural constant, and j is imaginary unit;
S504. rumble spectrum sequence U (f i) at f ithere is simple spike value in scope, find rumble spectrum sequence U (f i) the position q in the longitudinal direction of the photosignal corresponding to peak value, q ∈ [1,2M 2], adopt following formula to calculate every rank vibration frequency f of target i:
f i = Fs L × ( k i - 1 ) + Fs L × 2 M 2 q .
On the basis of technique scheme, described voltage sample rate FS=625KHz, accuracy class M 1=20, M 2=40.
On the basis of technique scheme, when the movement velocity estimating bridge test position is 0 ~ 15mm/s, the length N=65536 of setting photosignal, the sampling rate Fs=5Hz of speed; When to estimate bridge test position movement velocity be 16 ~ 50mm/s, the length N=32768 of setting photosignal, the sampling rate Fs=10Hz of speed; When to estimate bridge test position movement velocity be 51 ~ 100mm/s, the length N=16384 of setting photosignal, the sampling rate Fs=15Hz of speed; When to estimate bridge test position movement velocity be 101 ~ 250mm/s, the length N=8192 of setting photosignal, the sampling rate Fs=20Hz of speed.
Beneficial effect of the present invention is:
The present invention adds the further high precision computation of frequency on the basis of conventional fast Fourier transform algorithm, the extraction accuracy of light frequency is high, can when laser Doppler frequencies signal and noise ratio lower accurately extract frequency values, finally can meet the demand of bridge structure dynamic test; For the accuracy class value of low-speed motion different sample frequency, sampling number and the optimization by test target classification setting further again of bridge, ensure that the precision of the dynamic parameter extraction of bridge and the computing velocity of on-line measurement simultaneously, reduce the requirement to hardware computing power, reduce cost, for measuring low-speed motion speed and the vibration frequency of the engineering structures such as bridge.
Accompanying drawing explanation
Fig. 1 is the structural representation based on the bridge vibration detection device of laser in background technology of the present invention;
Fig. 2 is the schematic flow sheet of the dynamic parameter extracting method of bridge that the present invention is based on laser-Doppler.
In figure: 1-test box, 2-first fiber optic collimator mirror, 3-adjustable type clamper, 4-first fiber coupler, 5-optical-fiber laser transmitter, 6-first laser port, 7-second laser port, 8-the 3rd laser port, 9-second fiber optic collimator mirror, 10-laser pen, 11-second fiber coupler, 12-photodetector, 13-high-frequency signal Acquisition Instrument, 14-computing machine.
Embodiment
Below in conjunction with accompanying drawing, embodiments of the invention are described in further detail.
The embodiment of the present invention provides a kind of bridge based on laser-Doppler dynamic parameter extraction system, and this system comprises:
Photodetector 12, for detecting photosignal.
High-frequency signal Acquisition Instrument 13, for gathering the photosignal of certain length in photodetector, and being converted into magnitude of voltage, obtaining photosignal discrete series.
Computing module, for carrying out Fast Fourier Transform (FFT) to photosignal discrete series, obtains frequency movement sequence; Find the particular location of the photosignal corresponding to frequency movement sequence peaks, obtain frequency fine by stages; Segmented further frequency fine by stages, the frequency spectrum on calculated rate subdivided interval, obtains spectrum sequence; Find the particular location of the photosignal corresponding to peak value of spectrum sequence, calculate Doppler shift; According to Doppler shift, calculate target speed; Judge whether the vibration frequency needing to calculate target, judge whether the vibration frequency needing to calculate target, if not, then calculate movement velocity corresponding to different length photosignal; If so, then formed according to movement velocity and obtain the vibration frequency that velocity series calculates target; Fast Fourier Transform (FFT) is carried out to velocity series, obtains frequency vibration sequence; Find the particular location of the photosignal corresponding to frequency vibration sequence peaks, obtain vibration frequency subdivided interval; Frequency fine by stages is segmented further, calculates the rumble spectrum on vibration frequency subdivided interval, obtain rumble spectrum sequence, find the particular location of the photosignal corresponding to rumble spectrum sequence peaks, calculate the vibration frequency of target.
Shown in Figure 2, the embodiment of the present invention also provides a kind of bridge based on laser-Doppler for said system dynamic parameter extracting method, comprises the following steps:
S1. high-frequency signal Acquisition Instrument 13 gather certain length in photodetector 12 photosignal, and be converted into magnitude of voltage, obtain photosignal discrete series, detailed process is: in a period of time △ t, high-frequency signal Acquisition Instrument 13 is according to voltage sample rate FS, gather the photosignal in photodetector 12, and be converted into magnitude of voltage, the length of described photosignal is N, photosignal is at the corresponding magnitude of voltage in each position of length direction, obtain photosignal discrete series X={x (n), n=0, 1 ... N-1}, N is positive integer, the resolution of frequency is Δ f,
S2. computing machine carries out Fast Fourier Transform (FFT) to photosignal discrete series, obtains frequency movement sequence; Find the particular location of the photosignal corresponding to frequency movement sequence peaks, obtain frequency fine by stages; Segmented further frequency fine by stages, the frequency spectrum on calculated rate subdivided interval, obtains spectrum sequence; Find the particular location of the photosignal corresponding to peak value of spectrum sequence, calculate Doppler shift; According to Doppler shift, calculate target speed.
Step S2 specifically comprises the following steps:
S201. computing machine is to photosignal discrete series X={x (n), n=0, and 1 ... N-1} carries out Fast Fourier Transform (FFT), obtain frequency movement sequence Z=T (X)={ z (n), n=0,1 ... N-1}, wherein, T is Fourier transform operator;
S202. photosignal discrete series X={x (n) that collects at every batch of frequency movement sequence Z=T (X), n=0,1 ... simple spike value is there is in the scope of N-1}, find the photosignal particular location k in the longitudinal direction corresponding to frequency movement sequence Z=T (X) peak value, k ∈ [1, N], obtaining frequency fine by stages is
S203. according to accuracy class M 1, by frequency fine by stages be further subdivided into 2M 1equal portions, M 1for positive integer, further frequency analysis be spaced apart △, adopt the frequency spectrum on following formula calculated rate subdivided interval, obtain spectrum sequence Y (f):
Y ( f ) = Σ n = 0 N - 1 x ( n ) e - j 2 πfn ,
Wherein, Doppler shift e is natural constant, and j is imaginary unit;
S204. there is simple spike value in spectrum sequence Y (f) within the scope of f, finds the position p in the longitudinal direction of the photosignal corresponding to peak value of spectrum sequence Y (f), p ∈ [1,2M 1], adopt following formula to calculate photosignal discrete series X={x (n), n=0,1 ... the Doppler shift f of N-1}:
f = FS N × ( k - 1 ) + FS N × 2 M 1 p ;
S205. according to following formula, target speed v is calculated:
v = 1 2 f λ 0 ,
Wherein, λ 0for the wavelength of transmitted light of optical-fiber laser transmitter 5.
S3. judge whether the vibration frequency needing to calculate target, if not, then repeat step S1 ~ S2, calculate the movement velocity that different length photosignal is corresponding, that is, computational length is the movement velocity that the photosignal of l is corresponding respectively, wherein, l=0,1 ..., L-1, L be positive integer; If so, then step S4 is forwarded to.
S4. corresponding according to different length photosignal movement velocity, obtains velocity series, and detailed process is: be the movement speed v that the photosignal of l is corresponding according to length, forms velocity series V={v (l), l=0,1 ... L-1}, wherein, the sampling rate of speed is Fs the resolution of speed in frequency is new △ f, newly
S5. Fast Fourier Transform (FFT) is carried out to velocity series, obtain frequency vibration sequence; Find the particular location of the photosignal corresponding to frequency vibration sequence peaks, obtain vibration frequency subdivided interval; Frequency fine by stages is segmented further, calculates the rumble spectrum on vibration frequency subdivided interval, obtain rumble spectrum sequence, find the particular location of the photosignal corresponding to rumble spectrum sequence peaks, calculate the vibration frequency of target.
Step S5 specifically comprises the following steps:
S501. to velocity series V={v (l), l=0,1 ... L-1} carries out Fast Fourier Transform (FFT), obtain frequency vibration sequence W=T (V)=w (l), l=0,1...L-1}, wherein T is Fourier transform operator;
S502. frequency vibration sequence W=T (V) is at velocity series V={v (l), l=0,1 ... there is simple spike value in the scope of L-1}, find the particular location k of the photosignal corresponding to frequency vibration sequence W=T (V) peak value at length direction i, k i∈ [1, L], obtains vibration frequency subdivided interval: wherein, m is vibration exponent number, and m is positive integer;
S503. according to accuracy class M 2, by vibration frequency subdivided interval be further subdivided into 2M 2equal portions, M 2for positive integer, further vibration Frequency Analysis be spaced apart Δ 2, adopt the frequency spectrum on following formula calculating vibration frequency subdivided interval, obtain rumble spectrum sequence U (f i):
U ( f i ) = Σ l = 0 L - 1 v ( l ) e - j 2 π f i l ,
Wherein e is natural constant, and j is imaginary unit;
S504. rumble spectrum sequence U (f i) at f ithere is simple spike value in scope, find rumble spectrum sequence U (f i) the position q in the longitudinal direction of the photosignal corresponding to peak value, q ∈ [1,2M 2], adopt following formula to calculate every rank vibration frequency f of target i:
f i = Fs L × ( k i - 1 ) + Fs L × 2 M 2 q .
In said process, voltage sample rate FS=625KHz, accuracy class M 1=20, M 2=40.
When the movement velocity estimating bridge test position is 0 ~ 15mm/s, the length N=65536 of setting photosignal, the sampling rate Fs=5Hz of speed; When to estimate bridge test position movement velocity be 16 ~ 50mm/s, the length N=32768 of setting photosignal, the sampling rate Fs=10Hz of speed; When to estimate bridge test position movement velocity be 51 ~ 100mm/s, the length N=16384 of setting photosignal, the sampling rate Fs=15Hz of speed; When to estimate bridge test position movement velocity be 101 ~ 250mm/s, the length N=8192 of setting photosignal, the sampling rate Fs=20Hz of speed.
Below in conjunction with embodiment, embodiments of the invention are described in further detail.
Example is measured as, sampling rate FS=625KHz, accuracy class M with the span centre vertical displacement of certain bridge 1=20, M 2=40, according to design feature, the vertical displacement speed in estimation bridge span, at 101 ~ 250mm/s, sets N=8192, Fs=20Hz.The dynamic parameter extracting method of bridge based on laser-Doppler comprises the following steps:
S1. every △ t=50ms, signal in high-frequency signal Acquisition Instrument 13 pairs of photodetectors 12 carries out once batch sampling, in 50ms at front 13ms according to sampling rate FS=625KHz, gather N=8192 data, obtain photosignal discrete series X={x (1), x (2) ... x (8192) }, the resolution of frequency
S2. to X={x (1), x (2) ... x (8192) } carry out Fast Fourier Transform (FFT) (FFT conversion), calculate frequency movement sequence Z=T (X)={ z (1), z (2) ... z (8192) }, the peak k=3277 of Z=T (X), obtain the subdivided interval of frequency for [249052Hz, 255606Hz]; Accuracy class M 1=20, interval △=3.8Hz that cline frequency is analyzed, adopt the frequency spectrum on formulae discovery frequency fine by stages: wherein f ∈ [249052Hz, 255606Hz], goes forward one by one and is spaced apart △=3.8Hz, and e is natural constant, and j is imaginary unit; Find the peak value to Y (f) sequence, find the p=15 that crest frequency is corresponding; Calculate Doppler shift f = FS N × ( k - 1 ) + FS N × 2 M 1 p = 249052 + 3.8 × 15 = 249109 Hz ; Calculate target speed v = 1 2 f λ 0 = 1 2 × 249109 Hz × 1550 nm = 190 mm / s , Vertical displacement speed so within the sampling time of △ t=50ms in bridge span is 190mm/s.
S3. judge whether the vibration frequency needing to calculate target, if not, then repeat step S1 ~ S2, calculate the movement velocity that different length photosignal is corresponding, that is, computational length is the movement velocity that the photosignal of l is corresponding respectively, wherein, l=0,1 ..., L-1, L be positive integer; If so, then step S4 is forwarded to.
S4. corresponding according to different length photosignal movement velocity, obtains velocity series V={v (l), l=0,1 ... L-1}, wherein, the sampling rate of speed the resolution of speed in frequency is
S5. Fast Fourier Transform (FFT) is carried out to velocity series, obtain frequency vibration sequence; Find the particular location of the photosignal corresponding to frequency vibration sequence peaks, obtain vibration frequency subdivided interval; Frequency fine by stages is segmented further, calculates the rumble spectrum on vibration frequency subdivided interval, obtain rumble spectrum sequence, find the particular location of the photosignal corresponding to rumble spectrum sequence peaks, calculate the vibration frequency of target.
The present invention is not limited to above-mentioned embodiment, and for those skilled in the art, under the premise without departing from the principles of the invention, can also make some improvements and modifications, these improvements and modifications are also considered as within protection scope of the present invention.The content be not described in detail in this instructions belongs to the known prior art of professional and technical personnel in the field.

Claims (9)

1., based on the dynamic parameter extraction system of bridge of laser-Doppler, it is characterized in that, this system comprises:
Photodetector, for detecting photosignal;
High-frequency signal Acquisition Instrument, for gathering the photosignal of certain length in photodetector, and being converted into magnitude of voltage, obtaining photosignal discrete series;
Computing machine, for carrying out Fast Fourier Transform (FFT) to photosignal discrete series, obtains frequency movement sequence; Find the particular location of the photosignal corresponding to frequency movement sequence peaks, obtain frequency fine by stages; Segmented further frequency fine by stages, the frequency spectrum on calculated rate subdivided interval, obtains spectrum sequence; Find the particular location of the photosignal corresponding to peak value of spectrum sequence, calculate Doppler shift; According to Doppler shift, calculate target speed; Judge whether the vibration frequency needing to calculate target, judge whether the vibration frequency needing to calculate target, if not, then calculate movement velocity corresponding to different length photosignal; If so, then formed according to movement velocity and obtain the vibration frequency that velocity series calculates target; Fast Fourier Transform (FFT) is carried out to velocity series, obtains frequency vibration sequence; Find the particular location of the photosignal corresponding to frequency vibration sequence peaks, obtain vibration frequency subdivided interval; Frequency fine by stages is segmented further, calculates the rumble spectrum on vibration frequency subdivided interval, obtain rumble spectrum sequence, find the particular location of the photosignal corresponding to rumble spectrum sequence peaks, calculate the vibration frequency of target.
2., for the dynamic parameter extracting method of the bridge based on laser-Doppler of system described in claim 1, it is characterized in that, comprise the following steps:
S1. high-frequency signal Acquisition Instrument gathers the photosignal of certain length in photodetector, and is converted into magnitude of voltage, obtains photosignal discrete series;
S2. computing machine carries out Fast Fourier Transform (FFT) to photosignal discrete series, obtains frequency movement sequence; Find the particular location of the photosignal corresponding to frequency movement sequence peaks, obtain frequency fine by stages; Segmented further frequency fine by stages, the frequency spectrum on calculated rate subdivided interval, obtains spectrum sequence; Find the particular location of the photosignal corresponding to peak value of spectrum sequence, calculate Doppler shift; According to Doppler shift, calculate target speed;
S3. judge whether the vibration frequency needing to calculate target, if not, then repeat step S1 ~ S2, calculate the movement velocity that different length photosignal is corresponding; If so, then step S4 is forwarded to;
S4. corresponding according to different length photosignal movement velocity, obtains velocity series;
S5. Fast Fourier Transform (FFT) is carried out to velocity series, obtain frequency vibration sequence; Find the particular location of the photosignal corresponding to frequency vibration sequence peaks, obtain vibration frequency subdivided interval; Frequency fine by stages is segmented further, calculates the rumble spectrum on vibration frequency subdivided interval, obtain rumble spectrum sequence, find the particular location of the photosignal corresponding to rumble spectrum sequence peaks, calculate the vibration frequency of target.
3. as claimed in claim 2 based on the dynamic parameter extracting method of bridge of laser-Doppler, it is characterized in that, the detailed process of step S1 is: in a period of time Δ t, high-frequency signal Acquisition Instrument is according to voltage sample rate FS, gather the photosignal in photodetector, and be converted into magnitude of voltage, the length of described photosignal is N, N is positive integer, and photosignal, at the corresponding magnitude of voltage x in each position of length direction, obtains photosignal discrete series X={x (n), n=0,1 ... N-1}, the resolution of frequency is Δ f
4., as claimed in claim 3 based on the dynamic parameter extracting method of bridge of laser-Doppler, it is characterized in that, step S2 specifically comprises the following steps:
S201. computing machine is to photosignal discrete series X={x (n), n=0, and 1 ... N-1} carries out Fast Fourier Transform (FFT), obtain frequency movement sequence Z=T (X)={ z (n), n=0,1 ... N-1}, wherein, T is Fourier transform operator;
S202. photosignal discrete series X={x (n) that collects at every batch of frequency movement sequence Z=T (X), n=0,1 ... simple spike value is there is in the scope of N-1}, find the photosignal particular location k in the longitudinal direction corresponding to frequency movement sequence Z=T (X) peak value, k ∈ [1, N], obtaining frequency fine by stages is
S203. according to accuracy class M 1, by frequency fine by stages be further subdivided into 2M 1equal portions, M 1for positive integer, further frequency analysis be spaced apart Δ, adopt the frequency spectrum on following formula calculated rate subdivided interval, obtain spectrum sequence Y (f):
Wherein, Doppler shift e is natural constant, and j is imaginary unit;
S204. there is simple spike value in spectrum sequence Y (f) within the scope of f, finds the position p in the longitudinal direction of the photosignal corresponding to peak value of spectrum sequence Y (f), p ∈ [1,2M 1], adopt following formula to calculate photosignal discrete series X={x (n), n=0,1 ... the Doppler shift f of N-1}:
S205. according to following formula, target speed v is calculated:
Wherein, λ 0for the wavelength of transmitted light of optical-fiber laser transmitter.
5., as claimed in claim 4 based on the dynamic parameter extracting method of bridge of laser-Doppler, it is characterized in that: in step S3, if not, then repeat step S1 ~ S2, computational length is the movement velocity that the photosignal of l is corresponding respectively, wherein, l=0,1 ..., L-1, L be positive integer.
6. as claimed in claim 5 based on the dynamic parameter extracting method of bridge of laser-Doppler, it is characterized in that, the detailed process of step S4 is: be the movement speed v that the photosignal of l is corresponding according to length, form velocity series V={v (l), l=0,1, ... L-1}, wherein, the sampling rate of speed is Fs the resolution of speed in frequency is new Δ f, newly
7., as claimed in claim 6 based on the dynamic parameter extracting method of bridge of laser-Doppler, it is characterized in that, step S5 specifically comprises the following steps:
S501. to velocity series V={v (l), l=0,1 ... L-1} carries out Fast Fourier Transform (FFT), obtain frequency vibration sequence W=T (V)=w (l), l=0,1...L-1}, wherein T is Fourier transform operator;
S502. frequency vibration sequence W=T (V) is at velocity series V={v (l), l=0,1 ... there is simple spike value in the scope of L-1}, find the particular location k of the photosignal corresponding to frequency vibration sequence W=T (V) peak value at length direction i, k i∈ [1, L], obtains vibration frequency subdivided interval: wherein, m is vibration exponent number, and m is positive integer;
S503. according to accuracy class M 2, by vibration frequency subdivided interval be further subdivided into 2M 2equal portions, M 2for positive integer, further vibration Frequency Analysis be spaced apart Δ 2, adopt the frequency spectrum on following formula calculating vibration frequency subdivided interval, obtain rumble spectrum sequence U (f i):
Wherein e is natural constant, and j is imaginary unit;
S504. rumble spectrum sequence U (f i) at f ithere is simple spike value in scope, find rumble spectrum sequence U (f i) the position q in the longitudinal direction of the photosignal corresponding to peak value, q ∈ [1,2M 2], adopt following formula to calculate every rank vibration frequency f of target i:
8., as claimed in claim 7 based on the dynamic parameter extracting method of bridge of laser-Doppler, it is characterized in that: described voltage sample rate FS=625KHz, accuracy class M 1=20, M 2=40.
9. as claimed in claim 7 based on the dynamic parameter extracting method of bridge of laser-Doppler, it is characterized in that: when the movement velocity estimating bridge test position is 0 ~ 15mm/s, the length N=65536 of setting photosignal, the sampling rate Fs=5Hz of speed; When to estimate bridge test position movement velocity be 16 ~ 50mm/s, the length N=32768 of setting photosignal, the sampling rate Fs=10Hz of speed; When to estimate bridge test position movement velocity be 51 ~ 100mm/s, the length N=16384 of setting photosignal, the sampling rate Fs=15Hz of speed; When to estimate bridge test position movement velocity be 101 ~ 250mm/s, the length N=8192 of setting photosignal, the sampling rate Fs=20Hz of speed.
CN201510173926.5A 2015-04-10 2015-04-10 Dynamic bridge parameter extracting system and dynamic bridge parameter extracting method based on laser Doppler Active CN104792364B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510173926.5A CN104792364B (en) 2015-04-10 2015-04-10 Dynamic bridge parameter extracting system and dynamic bridge parameter extracting method based on laser Doppler

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510173926.5A CN104792364B (en) 2015-04-10 2015-04-10 Dynamic bridge parameter extracting system and dynamic bridge parameter extracting method based on laser Doppler

Publications (2)

Publication Number Publication Date
CN104792364A true CN104792364A (en) 2015-07-22
CN104792364B CN104792364B (en) 2017-04-12

Family

ID=53557361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510173926.5A Active CN104792364B (en) 2015-04-10 2015-04-10 Dynamic bridge parameter extracting system and dynamic bridge parameter extracting method based on laser Doppler

Country Status (1)

Country Link
CN (1) CN104792364B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932135A (en) * 2017-05-17 2017-07-07 苏交科集团股份有限公司 The flexible cable cable tension test method that peak recognizes vibration frequency is searched based on weighting arrowband
CN107036751A (en) * 2017-05-17 2017-08-11 苏交科集团股份有限公司 The flexible strand cable force computational methods that peak recognizes vibration frequency are searched by weighting broadband
CN107192448A (en) * 2017-05-17 2017-09-22 苏交科集团股份有限公司 Peak method is searched in a kind of broadband for recognizing flexible rope vibration frequency
CN107607928A (en) * 2017-08-14 2018-01-19 北京理工大学 A kind of rotor laser-Doppler and micro-doppler composite signal emulator
CN110346035A (en) * 2019-06-28 2019-10-18 中铁大桥科学研究院有限公司 Bridge real-time frequency test method and system
CN110361723A (en) * 2019-07-22 2019-10-22 深圳锐越微技术有限公司 The time-frequency characteristics extracting method of Doppler radar motion target
CN117168604A (en) * 2023-09-04 2023-12-05 中冶建筑研究总院有限公司 Doppler vectorization test method for structural vibration frequency

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6240783B1 (en) * 1998-05-02 2001-06-05 Usbi, Co Bridge monitoring system
CN101315274A (en) * 2008-06-19 2008-12-03 西安交通大学 Monitoring device and real-time monitoring method for bridge vibration deformation
CN202101754U (en) * 2011-05-25 2012-01-04 黑龙江工程学院 Bridge vibration detection device
CN102564569A (en) * 2011-12-26 2012-07-11 西北工业大学 Bridge vibration frequency detection method based on deeply-optimized particle filter
CN103335858A (en) * 2013-06-06 2013-10-02 湖南大学 Method for measuring bridge structure dynamic displacement and vibration frequency
JP2014109536A (en) * 2012-12-04 2014-06-12 Nihon Univ Nondestructive inspection system for concrete structure using distortion of concrete, monitoring system and nondestructive inspection method
CN203929213U (en) * 2014-07-11 2014-11-05 公安部第一研究所 A kind of multipoint mode laser-Doppler vibration measuring system
CN104132693A (en) * 2014-08-06 2014-11-05 电子科技大学 Method for simultaneously extracting position and frequency of vibration signal in phase OTDR system
CN204255494U (en) * 2014-12-18 2015-04-08 中国神华能源股份有限公司 Bridge vibration monitoring device

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6240783B1 (en) * 1998-05-02 2001-06-05 Usbi, Co Bridge monitoring system
CN101315274A (en) * 2008-06-19 2008-12-03 西安交通大学 Monitoring device and real-time monitoring method for bridge vibration deformation
CN202101754U (en) * 2011-05-25 2012-01-04 黑龙江工程学院 Bridge vibration detection device
CN102564569A (en) * 2011-12-26 2012-07-11 西北工业大学 Bridge vibration frequency detection method based on deeply-optimized particle filter
JP2014109536A (en) * 2012-12-04 2014-06-12 Nihon Univ Nondestructive inspection system for concrete structure using distortion of concrete, monitoring system and nondestructive inspection method
CN103335858A (en) * 2013-06-06 2013-10-02 湖南大学 Method for measuring bridge structure dynamic displacement and vibration frequency
CN203929213U (en) * 2014-07-11 2014-11-05 公安部第一研究所 A kind of multipoint mode laser-Doppler vibration measuring system
CN104132693A (en) * 2014-08-06 2014-11-05 电子科技大学 Method for simultaneously extracting position and frequency of vibration signal in phase OTDR system
CN204255494U (en) * 2014-12-18 2015-04-08 中国神华能源股份有限公司 Bridge vibration monitoring device

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
何鹍: "时域模态分析方法在桥梁动态无损检测中的应用研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *
朱靖等: "桥梁振动的激光干涉测量方法", 《物理与工程》 *
李妍: "关头坝大桥振动监测系统研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *
杨焕龙: "基于环境振动下模态试验分析方法在桥梁结构中的应用", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932135A (en) * 2017-05-17 2017-07-07 苏交科集团股份有限公司 The flexible cable cable tension test method that peak recognizes vibration frequency is searched based on weighting arrowband
CN107036751A (en) * 2017-05-17 2017-08-11 苏交科集团股份有限公司 The flexible strand cable force computational methods that peak recognizes vibration frequency are searched by weighting broadband
CN107192448A (en) * 2017-05-17 2017-09-22 苏交科集团股份有限公司 Peak method is searched in a kind of broadband for recognizing flexible rope vibration frequency
CN106932135B (en) * 2017-05-17 2022-09-23 苏交科集团股份有限公司 Flexible inhaul cable force testing method for identifying vibration frequency based on weighted narrow-band peak searching method
CN107607928A (en) * 2017-08-14 2018-01-19 北京理工大学 A kind of rotor laser-Doppler and micro-doppler composite signal emulator
CN107607928B (en) * 2017-08-14 2020-06-12 北京理工大学 Rotor wing laser Doppler and micro Doppler composite signal simulator
CN110346035A (en) * 2019-06-28 2019-10-18 中铁大桥科学研究院有限公司 Bridge real-time frequency test method and system
CN110361723A (en) * 2019-07-22 2019-10-22 深圳锐越微技术有限公司 The time-frequency characteristics extracting method of Doppler radar motion target
CN110361723B (en) * 2019-07-22 2021-11-30 深圳锐越微技术有限公司 Time-frequency feature extraction method for Doppler radar moving target
CN117168604A (en) * 2023-09-04 2023-12-05 中冶建筑研究总院有限公司 Doppler vectorization test method for structural vibration frequency

Also Published As

Publication number Publication date
CN104792364B (en) 2017-04-12

Similar Documents

Publication Publication Date Title
CN104792364A (en) Dynamic bridge parameter extracting system and dynamic bridge parameter extracting method based on laser Doppler
Grare et al. Growth and dissipation of wind-forced, deep-water waves
CN101586997A (en) Method for calculating guy cable vibrating base frequency
CN103196407A (en) Method, device and system for vibration displacement measurement of cantilever crane of pump truck and engineering machinery device
CN101183050B (en) Electrohydraulic servo valve dynamic performance testing method for measuring displacement
CN103808406A (en) Oil-gas pipeline vibration monitoring method and device based on vibration wire type sensor
CN1831485A (en) Cavity length demodulating algorithm of fibre-optical F-P sensor
CN103412137B (en) With speed-measuring method and device in twiddle factor
CN102955004B (en) Subway tunnel segment service performance detection method based on wave velocity determination
CN104316160A (en) Underwater sound signal instantaneous frequency demodulation method based on wavelet ridges
Soulier et al. Low-Reynolds-number investigations on the ability of the strip of e-TellTale sensor to detect the flow features over wind turbine blade section: flow stall and reattachment dynamics
CN102269803A (en) Method for correcting low-frequency components in discrete spectrum based on time delay
CN102298072B (en) High precision wind measuring device with micro-differential pressure type and method thereof
CN101718798B (en) Debris flow velocity measurement method and implementing system thereof
CN103630604B (en) The recognition methods of centrifugal compressor half-opened impeller crack fault
CN101334325A (en) Cable force vibration detection method and its detecting apparatus
CN106199060B (en) Dual-Phrase Distribution of Gas olid speed measurement method based on sliding average and capacitance sensor
Zhang et al. The method for determining optimal analysis length of vibration data based on improved multiscale permutation entropy
CN201083518Y (en) Waveform height measuring systems
CN202150036U (en) Measurement and control system facing solid-liquid double-phase soft abrasive particle flow precision processing
CN103884862A (en) Secondary correlation time delay estimation method for monitoring wind speed of supersonic wave of wind power station
Fairhurst Modelling and design of an oscillating wave energy converter
CN102539825A (en) Wind speed spectrum acquisition method based on wind speed re-sampling technology
CN109870389A (en) Based on magnetostrictive displacement sensor irrigation water coefficient of viscosity detection device and method
Chen et al. Automatic transformation between water level meter measurement data and altitude on the large-scale physical model

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 103 No. 430000 Hubei Province in Qiaokou District of Wuhan city road construction

Applicant after: BRIDGE SCIENCE AND TECHNOLOGY INSTITUTE CO., LTD., CHINA RAILWAY ENGINEERING GROUP

Applicant after: Group Co., Ltd., Zhongtie Daqiao Bureau

Address before: 103 No. 430000 Hubei Province in Qiaokou District of Wuhan city road construction

Applicant before: China Zhongtie Major Bridge Engineering Group Co., Ltd.

Applicant before: Group Co., Ltd., Zhongtie Daqiao Bureau

COR Change of bibliographic data
GR01 Patent grant
GR01 Patent grant