CN106404159B - A kind of determining method of Rocket sled test continuous shaking transmitting spectrum - Google Patents
A kind of determining method of Rocket sled test continuous shaking transmitting spectrum Download PDFInfo
- Publication number
- CN106404159B CN106404159B CN201610423010.5A CN201610423010A CN106404159B CN 106404159 B CN106404159 B CN 106404159B CN 201610423010 A CN201610423010 A CN 201610423010A CN 106404159 B CN106404159 B CN 106404159B
- Authority
- CN
- China
- Prior art keywords
- equipment
- sequence
- vibration
- oscillation power
- spectrum
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
A kind of determining method of Rocket sled test continuous shaking transmitting spectrum, belongs to technical field of mechanical environmental test.Steps are as follows: (1) going out the acceleration value sequence a of equipment X and equipment Y at runtime by sensor measurementIx(k)、aIy(k)；(2) aIx(k)、aIy(k) it is used as input quantity, vibration value sequence a is acquired after high-pass filteringdx(k), ady(k)；(3) it is based on AR model solution adx(k) and ady(k) the continuous function P of oscillation power spectrumx(f)、Py(f)；(4) oscillation power spectral continuous function P is correctedx(f)、Py(f) amplitude, obtains P 'x(f)、P′y(f)；(5) spectrum is transmitted in the continuous shaking for calculating equipment X to equipment YThis method can directly acquire high-frequency vibration information, effectively remove noise, realize the quantificational description of the vibration transmissibilities at different levels with frequency consecutive variations, and precision is high.
Description
Technical field
It is transmitted the present invention relates to a kind of Rocket sled test continuous shaking and composes the method for determination, belong to dynamic environmental test technology neck
Domain.
Background technique
Rocket sledge is to use rocket engine as thrust, the sled of high-speed cruising in orbit.Since rocket sledge is being run
It can produce big overload in the process, therefore, different types of product, which is all mounted in skid body, carries out Rocket sled test.
During Rocket sled test, for subject product according to test requirements document there are two types of mounting means, first way is subject
Product is directly mounted in skid body；The second way is that subject product is installed on vibration damping plate, and vibration damping plate passes through metal
Damper is installed in skid body.By vibration test system real-time measurement and record skid body and subject product (second way
Words further include vibration damping plate) vibration and impact signal, it is to be tested after read record storage data and carry out at data
Reason.
Vibration test system is made of sensor, data acquisition logging system, trigger device and signal-transmitting cable etc..It surveys
Sensor, data collection system and trigger device are connected by signal-transmitting cable in test system, peace of each component in skid body
The mode of press strip, bolt fastening can be used in dress, and specific installation position includes:
A) installation is no less than the orthogonal sensor of 3 input shafts respectively on skid body, vibration damping plate and subject product shell,
Course acceleration, vertical acceleration and the side acceleration of installed position are measured respectively；
B) data collection system, trigger device are directly mounted in skid body；
C) signal-transmitting cable is connected with skid body nearby using bunchy is tied up according to wire laying mode.
Currently, having the following problems in data processing:
(1) since the acceleration information of sensor measurement includes low frequency overload messages and high-frequency vibration information, to obtain
Vibration information is filtered acceleration information to obtain overload messages using low-pass filter first, and then, acceleration information subtracts
Overload messages are gone to obtain vibration information.And lacks and realize that high-frequency vibration information directly acquires method；
(2) Cyclical power spectrum analysis has only been carried out to vibration informations at different levels, feature is that noise is big, is lacked at effective noise reduction
Reason method；
It (3) is qualitative analysis vibration level size at different levels, whether there is or not resonance peak etc., it cannot be flat to vibration damping to skid body
Platform, vibration reduction platform to subject product, skid body to subject product vibration transmissibility analyzed；
(4) vibration transmitting can only realize discrete point in frequency, and can not achieve the analysis of cline frequency section, therefore, resolution ratio
It is lower.
Summary of the invention
Technical problem solved by the present invention is overcoming the deficiencies of the prior art and provide a kind of Rocket sled test continuous shaking
The method of determination is composed in transmitting, can be directly acquired high-frequency vibration information, effectively removed noise, realize vibration transmissibilities at different levels
Continuous quantificational description, precision is high.
The technical scheme is that the method for determination is composed in a kind of Rocket sled test continuous shaking transmitting, steps are as follows:
(1) during Rocket sled test, acceleration value when equipment X and equipment Y is run is acquired by sensor respectively,
Obtain the acceleration value sequence a of equipment XIx(k) and the acceleration value sequence a of equipment YIy(k), k is sampling instant, k=1,
The acceleration value of 2 ..., N, the equipment X and equipment Y include random vibration information and overload messages, the sensor sample
Frequency fsMeet 4KHz < fs< 10KHz, sampling period Ts=1/fs；
(2) using high-pass filter to acceleration value sequence aIx(k) it is filtered, obtains the vibration value sequence a of equipment Xdx
(k)；Using high-pass filter to acceleration value sequence aIy(k) it is filtered, obtains the vibration value sequence a of equipment Ydy(k)；
(3) a is calculated based on vibration AR modeldx(k) oscillation power spectral continuous function Px(f) and ady(k) vibration function
Rate spectral continuous function Py(f)；
(4) amplitude correction is carried out to the oscillation power spectrum that step (3) obtains respectively, obtains revised oscillation power spectrum and connects
Continuous function P 'x(f) and P 'y(f)；
(5) formula is utilizedDetermine the continuous shaking transmitting spectrum Φ of equipment X to equipment Yxy(f), described to set
The continuous shaking transmitting spectrum of standby X to equipment Y is the company that skid body transmits spectrum, skid body to vibration reduction platform to the continuous shaking of subject product
Spectrum is transmitted in the continuous shaking of continuous vibration transmitting spectrum or vibration reduction platform to subject product.
High-pass filter is filtered acceleration value sequence to obtain the implementation method of vibration value sequence in the step (2)
Are as follows:
(2.1) corner frequency for defining high-pass filter is fz, damped coefficient ξ, then the coefficient of high-pass filter be
(2.2) when k=1, the state variable x of high-pass filter1(1)=0, x2(1)=0；
(2.3) high-pass filter utilizes formula ad(k)=(n1-d1)x1(k)+(n2-d2)x2(k)+b0aI(k) to k-th when
The acceleration value a at quarterI(k) it is filtered to obtain vibration values ad(k), (2.4) are entered step；
(2.4) value of k adds 1, when judging whether k is less than or equal to N, if it is, utilizing formulaUpdate x1(k) and x2(k), return step (2.3)；Otherwise, high
Bandpass filter filtering terminates, and obtains vibration value sequence.
The implementation of the step (3) are as follows:
(3.1) a is setdx(k) autocorrelation sequence is Rx(i), ady(k) autocorrelation sequence is Ry(i), i=0,1,2 ...,
N-1, Rx(i) and Ry(i) calculation formula is as follows:
(3.2) q < N is taken, according to autocorrelation sequence Rx(i) a is calculateddx(k) each term coefficient a of q rank AR model1、a2、…、
aq, according to autocorrelation sequence Ry(i) a is calculateddy(k) each term coefficient b of q rank AR model1、b2、…、bq, calculation formula is as follows:
σx 2=Rx(0)
σy 2=Ry(0)；
(3.3) a is setdxIt (k) is the white noise sequence u for being 1 by variance1(n) motivate after as a result, then based on vibration AR model
Calculate adx(k) oscillation power spectral continuous function Px(f) formula are as follows:
Wherein q is the order for vibrating AR model, and f is frequency, and j is real number unit, j2=-1；
(3.4) a is setdyIt (k) is the white noise sequence u for being 1 by variance2(n) motivate after as a result, then based on vibration AR model
Calculate ady(k) oscillation power spectral continuous function Py(f) formula are as follows:
The implementation of the step (4) are as follows:
(4.1) f is takenp=(p-1)/T, wherein p=1,2 ..., N/2, substitution step (3.3) and (3.4) can acquire a respectivelydx
(k) the power spectrum sequence P of oscillation power spectrumx(fp) and ady(k) oscillation power spectral sequence Py(fp), T is N number of sampling week
The duration of phase；
(4.2) according to adx(k) the oscillation power spectrum density sequence Φ of equipment X is obtainedx(fp), according to ady(k) equipment Y is obtained
Oscillation power spectrum density sequence Φy(fp)；
(4.3) formula is utilizedCalculate the amplitude correction system of the oscillation power spectral continuous function of equipment X
Number cx, utilize formulaCalculate the amplitude correction coefficient c of the oscillation power spectral continuous function of equipment Yy；
(4.4) the revised oscillation power spectral continuous function P of equipment Xx' (f)=cxPx(f), the revised vibration of equipment Y
Power spectral continuous function Py' (f)=cyPy(f)。
According to a in the step (4.2)dx(k) the oscillation power spectrum density sequence Φ of equipment X is obtainedx(fp) realization side
Formula are as follows:
(i) it is based on positive frequency domain, the line spectrum sequence c of equipment X vibration data is calculated using following formulax(fp):
Wherein fpFor frequency, fp=(p-1)/T, p=1,2 ..., N/2；
(ii) the oscillation power spectrum density sequence Φ of equipment X is calculated according to line spectrum sequencex(fp):
Wherein F (p) is Discrete Fourier Transform；
According to ady(k) the oscillation power spectrum density sequence Φ of equipment Y is obtainedy(fp) implementation are as follows:
(i) the line spectrum sequence c of the vibration data of equipment Y is calculated using following formulay(fp):
(ii) the oscillation power spectrum density sequence Φ of equipment Y is calculated according to line spectrum sequencey(fp):
Compared with the prior art, the invention has the advantages that:
(1) present invention uses the method based on AR model spectra estimation, and it is continuous to give a kind of accurate description Rocket sled test
The calculation method of vibration transmitting spectrum overcomes the vibration transmitting qualitative description between at different levels in the past or only realizes discrete point in frequency,
And can not achieve the deficiency of cline frequency piecewise analysis, it realizes to skid body to vibration reduction platform, vibration reduction platform to subject product, skid body
To subject product with the quantificational description of the vibration transmissibility of frequency consecutive variations, precision is high.
(2) acceleration value of inventive sensor acquisition is by realizing and believing high-frequency vibration after high pass filter filters
The direct extraction of breath improves the transfer efficiency of data relative to the method overloaded and seek difference again is first extracted.
(3) present invention obtains power spectral continuous function using vibration AR model and its amplitude correction method, both can be accurate
Power spectrum is provided in the peak value and valley information of each Frequency point, and can accurately provide amplitude information, is realized in Rocket sled test
The accurate description of vibration transmitting spectrum of no resolution ratio limitation.
(4) relatively traditional Cyclical power spectrum has the characteristics that noise is big, and the present invention passes through vibration AR model and its amplitude
Modification method carries out noise reduction process to vibration data, effectively reduces the noise of vibration.
Detailed description of the invention
The flow chart of Fig. 1 this method；
Fig. 2 power spectrum of test item based on cyclic spectrum in Rocket sled test；
Fig. 3 test item in Rocket sled test is composed based on the initial power before AR Modifying model；
Fig. 4 test item in Rocket sled test is composed based on the accurate power after AR Modifying model；
Fig. 5 is the oscillation power spectrum density of skid body in Rocket sled test；
Fig. 6 is the oscillation power spectrum density of vibration damping plate in Rocket sled test；
Fig. 7 is the oscillation power spectrum density in Rocket sled test in subject product；
Fig. 8 is that Spectral structure is transmitted in vibration of the vibration damping plate under skid body vibrational excitation in Rocket sled test；
Fig. 9 transmits Spectral structure with respect to the vibration under the excitation of vibration damping plate vibration for subject product in Rocket sled test；
Figure 10 is that Spectral structure is transmitted in vibration of the subject product under skid body vibrational excitation in Rocket sled test.
Specific embodiment
1 couple of present invention elaborates with reference to the accompanying drawing.
Rocket sled test continuous shaking transmitting spectrum includes continuous shaking transmitting spectrum, skid body to vibration damping of the skid body to subject product
Spectrum is transmitted in the continuous shaking of the continuous shaking transmitting spectrum or vibration reduction platform to subject product of platform.Therefore, as shown in Figure 1, this hair
It is bright propose a kind of Rocket sled test during, the method for determination is composed in the continuous shaking transmitting of equipment X to equipment Y, and steps are as follows:
1, during Rocket sled test, acceleration value when equipment X and equipment Y is run is acquired by sensor respectively,
Obtain the acceleration value sequence a of equipment XIx(k) and the acceleration value sequence a of equipment YIy(k), k is sampling instant, k=1,
The acceleration value of 2 ..., N, the equipment X and equipment Y include random vibration information and overload messages, the sensor sample
Frequency fsMeet 4KHz < fs< 10KHz, sampling period Ts=1/fs。
By taking the vibration of skid body to test item transmitting spectrum as an example, the sensor for being installed on skid body samples sledge in each sampling period
Exciting source acceleration value when running body obtains exciting source acceleration value sequence a when skid body operationIx(k)；Meanwhile it being installed on
Acceleration value of the sensor of test item in sampling of each sampling period test item operation, obtains acceleration when test item operation
Angle value sequence aIy(k)。
2, using high-pass filter to acceleration value sequence aIx(k) it is filtered, obtains the vibration value sequence a of equipment Xdx
(k)；Using high-pass filter to acceleration value sequence aIy(k) it is filtered, obtains the vibration value sequence a of equipment Ydy(k), real
Existing method are as follows:
(2.1) corner frequency for defining high-pass filter is fz, damped coefficient ξ, then the coefficient of high-pass filter be
(2.2) when k=1, the state variable x of high-pass filter1(1)=0, x2(1)=0；
(2.3) high-pass filter utilizes formula ad(k)=(n1-d1)x1(k)+(n2-d2)x2(k)+b0aI(k) to k-th when
The acceleration value a at quarterI(k) it is filtered to obtain vibration values ad(k), (2.4) are entered step；
(2.4) value of k adds 1, when judging whether k is less than or equal to N, if it is, utilizing formulaUpdate x1(k) and x2(k), return step (2.3)；Otherwise, high
Bandpass filter filtering terminates, and obtains vibration value sequence.
3, a is calculateddx(k) the oscillation power spectral continuous function P based on AR modelx(f) and ady(k) vibration based on AR model
Dynamic power spectral continuous function Py(f), method are as follows:
(3.1) a is setdx(k) autocorrelation sequence is Rx(i), ady(k) autocorrelation sequence is Ry(i), i=0,1,2 ...,
N-1, Rx(i) and Ry(i) calculation formula is as follows:
(3.2) q < N is taken, according to autocorrelation sequence Rx(i) a is calculateddx(k) each term coefficient a of q rank AR model1、a2、…、
aq, according to autocorrelation sequence Ry(i) a is calculateddy(k) each term coefficient b of q rank AR model1、b2、…、bq, calculation formula is as follows:
σx 2=Rx(0)
σy 2=Ry(0)；
(3.3) a is setdxIt (k) is the white noise sequence u for being 1 by variance1(n) motivate after as a result, then based on vibration AR model
(autoregression model, auto-regressive) calculates adx(k) oscillation power spectral continuous function Px(f) formula are as follows:
Wherein q is the order that equipment X vibrates AR model, and f is frequency, and j is real number unit, j2=-1.
(3.4) a is setdyIt (k) is the white noise sequence u for being 1 by variance2(n) motivate after as a result, then based on vibration AR model
Calculate ady(k) oscillation power spectral continuous function Py(f) formula are as follows:
4, amplitude correction is carried out to the power spectrum that step (3) obtains, obtains revised power spectral continuous function P 'x(f)
With P 'y(f), method are as follows:
(4.1) f is takenp=(p-1)/T, wherein p=1,2 ..., N/2, can acquire the power of equipment X vibration AR model respectively
Spectral sequence Px(fp) and Y vibration AR model power spectrum sequence Py(fp)；
(4.2) according to adx(k) the oscillation power spectrum density sequence Φ of equipment X is obtainedx(fp), according to ady(k) equipment Y is obtained
Vibration values power spectrum degree series Φy(fp)；
According to adx(k) the oscillation power spectrum density sequence Φ of equipment X is obtainedx(fp) implementation are as follows:
(a) it is based on positive frequency domain, the line spectrum sequence c of equipment X vibration data is calculated using following formulax(fp):
Wherein fpFor frequency, fp=(p-1)/T, p=1,2 ..., N/2；
(b) the oscillation power spectrum density sequence Φ of equipment X is calculated according to line spectrum sequencex(fp):
Wherein T is the duration in N number of sampling period；F (p) is Discrete Fourier Transform；
According to ady(k) the oscillation power spectrum density sequence Φ of equipment Y is obtainedy(fp) implementation are as follows:
(a) the line spectrum sequence c of the vibration data of equipment Y is calculated using following formulay(fp):
(b) the oscillation power spectrum density sequence Φ of equipment Y is calculated according to line spectrum sequencey(fp):
(4.3) formula is utilizedCalculate the amplitude correction system of the oscillation power spectral continuous function of equipment X
Number cx, utilize formulaCalculate the amplitude correction coefficient c of the oscillation power spectral continuous function of equipment Yy；
(4.4) the revised oscillation power spectral continuous function P of equipment Xx' (f)=cxPx(f), the revised vibration of equipment Y
Power spectral continuous function Py' (f)=cyPy(f)。
5, the power spectral continuous function P ' obtained according to step (4)x(f)、P′y(f), formula is utilized
Determine the continuous shaking transmitting spectrum Φ of equipment X to equipment Yxy(f), wherein the continuous shaking transmitting spectrum of equipment X to equipment Y is
The continuous shaking transmitting spectrum, the continuous shaking transmitting spectrum of skid body to vibration reduction platform or vibration reduction platform of skid body to subject product are extremely tested
Spectrum is transmitted in the continuous shaking of product.
The present invention in order to obtain continuously compose by vibration transmitting, uses AR model spectra estimation in Rocket sled test.Due to AR mould
Type is a rational fraction, can realize the analysis to cline frequency section using smoothness properties, while when order is higher, each frequency
The wave crest and trough of point are also more accurate.But the shortcomings that AR model is whole amplitude, and by whether random signal steady etc., factors are influenced
And opposite true value has deviation, therefore, present invention further proposes the methods for carrying out amplitude correction to AR model, to adapt to fire
The Non-stationary Data processing requirement of arrow sledge test.Believe to both can accurately provide power spectrum in the peak value and valley of each Frequency point
Breath, and can accurately provide amplitude information, realize in Rocket sled test, skid body to vibration reduction platform, vibration reduction platform to subject product,
The accurate description of vibration transmitting spectrum that skid body is limited to subject product without resolution ratio, precision are high.
Embodiment 1:
By taking test item as an example, the validity of power spectrum amplitude correction method proposed by the present invention is verified.
Wherein impose a condition as follows: in Rocket sled test, sensor sample frequency is 9.6kHz, and test item is based on the period
The oscillation power of spectrum compose as shown in Fig. 2, the initial power spectrum based on AR model fig. 3, it is shown that the two amplitude some
Perhaps difference.Using the method for the present invention, power spectrum amplitude correction coefficient is c=0.723, after which is substituted into correction formula, subject
Accurate power spectrum of the product based on AR model is as shown in Figure 4.Comparing Fig. 3 and Fig. 4 can be seen that the power spectrum curve of Fig. 4 with respect to Fig. 3
Power spectrum curve it is whole toward moving down, amplitude is corrected, and peak valley shape remains unchanged.Compare Fig. 4 and Fig. 2 and can be seen that and adopts
With the power spectrum relative smooth of AR model, noise is decayed.
Embodiment 2:
The transmitting spectrum of the continuous shaking between at different levels is calculated from the vibration information in Rocket sled test according to multisensor, is realized
Vibration transmissibility analysis.
It imposes a condition as follows: in Rocket sled test, using multistage mounting means, it may be assumed that it is flat that subject product is installed on vibration damping
On plate, then by vibration damping plate skid body is installed on by vibration isolator.Sensor sample frequency is 9.6kHz.It is bringing into operation to stopping
Process, by skid body installation course transmitter, vibration damping plate installation course transmitter and inertial platform (subject product)
The corresponding acceleration value information of course transmitter difference acquisition and recording of installation, obtains vibration data by high-pass filter, and
Further obtain power spectral density.
The oscillation power spectrum density information of skid body as shown in Figure 5, oscillation power spectrum density information such as Fig. 6 institute of vibration damping plate
Show, the oscillation power spectrum density information of inertial platform (subject product) it is as shown in Figure 7.Continuous shaking transmitting spectrum according to the present invention
Calculation method can calculate separately out:
(1) Spectral structure, such as Fig. 8 are transmitted in the continuous shaking of vibration damping plate to skid body (vibration damping plate is under skid body vibrational excitation)
It is shown；
(2) to vibration damping plate, (inertial platform (subject product) motivates inertial platform (subject product) in vibration damping plate vibration
Under) continuous shaking transmit Spectral structure, as shown in Figure 9；
(3) inertial platform (subject product) is continuous to skid body (inertial platform (subject product) is under skid body vibrational excitation)
Vibration transmitting Spectral structure, as shown in Figure 10.
It can be seen from the figure that inertial platform (subject product) amplifies relative to the vibration level of vibration damping plate, due to subtracting
The attenuation of vibration plate is limited, so that the vibration level of inertial platform (subject product) when skid body vibrational excitation be caused also to put
Greatly.
The above, a specific embodiment only of the invention, but scope of protection of the present invention is not limited thereto, appoints
In the technical scope disclosed by the present invention, any changes or substitutions that can be easily thought of, all by what those familiar with the art
It is covered by the protection scope of the present invention.
Unspecified part of the present invention belongs to common sense well known to those skilled in the art.
Claims (5)
1. the method for determination is composed in a kind of Rocket sled test continuous shaking transmitting, it is characterised in that steps are as follows:
(1) during Rocket sled test, acceleration value when equipment X and equipment Y is run is acquired by sensor respectively, is obtained
The acceleration value sequence a of equipment XIx(k) and the acceleration value sequence a of equipment YIy(k), k is sampling instant, k=1,2 ..., N,
The acceleration value of the equipment X and equipment Y includes random vibration information and overload messages, the sensor sample frequency fsIt is full
Sufficient 4KHz < fs< 10KHz, sampling period Ts=1/fs；
(2) using high-pass filter to acceleration value sequence aIx(k) it is filtered, obtains the vibration value sequence a of equipment Xdx(k)；
Using high-pass filter to acceleration value sequence aIy(k) it is filtered, obtains the vibration value sequence a of equipment Ydy(k)；
(3) a is calculated based on vibration AR modeldx(k) oscillation power spectral continuous function Px(f) and ady(k) oscillation power spectrum
Continuous function Py(f)；
(4) amplitude correction is carried out to the oscillation power spectrum that step (3) obtains respectively, obtains revised oscillation power and composes continuous letter
Number P 'x(f) and P 'y(f)；
(5) formula is utilizedDetermine the continuous shaking transmitting spectrum Φ of equipment X to equipment Yxy(f), the equipment X
To equipment Y continuous shaking transmitting spectrum be skid body to subject product continuous shaking transmitting spectrum, skid body to vibration reduction platform it is continuous
Spectrum is transmitted in the continuous shaking of vibration transmitting spectrum or vibration reduction platform to subject product.
2. the method for determination is composed in a kind of Rocket sled test continuous shaking transmitting according to claim 1, it is characterised in that: described
High-pass filter is filtered acceleration value sequence to obtain the implementation method of vibration value sequence in step (2) are as follows:
(2.1) corner frequency for defining high-pass filter is fz, damped coefficient ξ, then the coefficient of high-pass filter be
(2.2) when k=1, the state variable x of high-pass filter1(1)=0, x2(1)=0；
(2.3) high-pass filter utilizes formula ad(k)=(n1-d1)x1(k)+(n2-d2)x2(k)+b0aI(k) to k-th moment
Acceleration value aI(k) it is filtered to obtain vibration values ad(k), (2.4) are entered step；
(2.4) value of k adds 1, when judging whether k is less than or equal to N, if it is, utilizing formulaUpdate x1(k) and x2(k), return step (2.3)；Otherwise, high
Bandpass filter filtering terminates, and obtains vibration value sequence.
3. the method for determination is composed in a kind of Rocket sled test continuous shaking transmitting according to claim 1, it is characterised in that: described
The implementation of step (3) are as follows:
(3.1) a is setdx(k) autocorrelation sequence is Rx(i), ady(k) autocorrelation sequence is Ry(i), i=0,1,2 ..., N-1,
Rx(i) and Ry(i) calculation formula is as follows:
(3.2) q < N is taken, according to autocorrelation sequence Rx(i) a is calculateddx(k) each term coefficient a of q rank AR model1、a2、…、aq, root
According to autocorrelation sequence Ry(i) a is calculateddy(k) each term coefficient b of q rank AR model1、b2、…、bQ,Calculation formula is as follows:
σx 2=Rx(0)
σy 2=Ry(0)；
(3.3) a is setdxIt (k) is the white noise sequence u for being 1 by variance1(n) after motivating to be calculated as a result, being then based on vibration AR model
adx(k) oscillation power spectral continuous function Px(f) formula are as follows:
Wherein q is the order for vibrating AR model, and f is frequency, and j is imaginary unit, j2=-1；
(3.4) a is setdyIt (k) is the white noise sequence u for being 1 by variance2(n) after motivating to be calculated as a result, being then based on vibration AR model
ady(k) oscillation power spectral continuous function Py(f) formula are as follows:
4. the method for determination is composed in a kind of Rocket sled test continuous shaking transmitting according to claim 3, it is characterised in that: described
The implementation of step (4) are as follows:
(4.1) f is takenp=(p-1)/T, wherein p=1,2 ..., N/2, substitution step (3.3) and (3.4) can acquire a respectivelydx(k)
Oscillation power spectrum power spectrum sequence Px(fp) and ady(k) oscillation power spectral sequence Py(fp), T is N number of sampling period
Duration；
(4.2) according to adx(k) the oscillation power spectrum density sequence Φ of equipment X is obtainedx(fp), according to ady(k) vibration of equipment Y is obtained
Dynamic power spectrum degree series Φy(fp)；
(4.3) formula is utilizedCalculate the amplitude correction coefficient c of the oscillation power spectral continuous function of equipment Xx,
Utilize formulaCalculate the amplitude correction coefficient c of the oscillation power spectral continuous function of equipment Yy；
(4.4) the revised oscillation power spectral continuous function P ' of equipment Xx(f)=cxPx(f), the revised oscillation power of equipment Y
Spectral continuous function P 'y(f)=cyPy(f)。
5. the method for determination is composed in a kind of Rocket sled test continuous shaking transmitting according to claim 4, it is characterised in that: described
According to a in step (4.2)dx(k) the oscillation power spectrum density sequence Φ of equipment X is obtainedx(fp) implementation are as follows:
(i) it is based on positive frequency domain, the line spectrum sequence c of equipment X vibration data is calculated using following formulax(fp):
Wherein fpFor frequency, fp=(p-1)/T, p=1,2 ..., N/2；
(ii) the oscillation power spectrum density sequence Φ of equipment X is calculated according to line spectrum sequencex(fp):
Wherein F (p) is Discrete Fourier Transform；
According to ady(k) the oscillation power spectrum density sequence Φ of equipment Y is obtainedy(fp) implementation are as follows:
(i) the line spectrum sequence c of the vibration data of equipment Y is calculated using following formulay(fp):
(ii) the oscillation power spectrum density sequence Φ of equipment Y is calculated according to line spectrum sequencey(fp):
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610423010.5A CN106404159B (en) | 2016-06-14 | 2016-06-14 | A kind of determining method of Rocket sled test continuous shaking transmitting spectrum |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610423010.5A CN106404159B (en) | 2016-06-14 | 2016-06-14 | A kind of determining method of Rocket sled test continuous shaking transmitting spectrum |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106404159A CN106404159A (en) | 2017-02-15 |
| CN106404159B true CN106404159B (en) | 2019-03-12 |
Family
ID=58006057
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610423010.5A Active CN106404159B (en) | 2016-06-14 | 2016-06-14 | A kind of determining method of Rocket sled test continuous shaking transmitting spectrum |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106404159B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108982140B (en) * | 2018-07-02 | 2019-10-11 | 西安交通大学 | A kind of rocket sledge track road spectrum reproducing apparatus and its method |
| CN109508504B (en) * | 2018-11-23 | 2023-08-04 | 南京航空航天大学 | Aerodynamic characteristic correction method for rocket sled pulley system test |
| CN110758394B (en) * | 2019-09-29 | 2020-10-16 | 一汽解放青岛汽车有限公司 | Engine starting control method and system under traffic jam working condition, vehicle and storage medium |
| CN110696807B (en) * | 2019-09-29 | 2021-04-09 | 一汽解放青岛汽车有限公司 | Engine shutdown control method under traffic jam condition, vehicle and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0318739A (en) * | 1989-06-16 | 1991-01-28 | Tech Res Assoc Openair Coal Min Mach | Trouble diagnostic method for speed reducer gear |
| CN102735267A (en) * | 2012-06-20 | 2012-10-17 | 北京航天控制仪器研究所 | Measuring method for inertial measurement device in sled testing |
| CN102937522A (en) * | 2012-08-30 | 2013-02-20 | 桂林电子科技大学 | Composite fault diagnosis method and system of gear case |
| CN103033329A (en) * | 2012-12-24 | 2013-04-10 | 北京工业大学 | Space grid structure frequency identification method of improved power spectrum peak method |
| CN104699980A (en) * | 2015-03-24 | 2015-06-10 | 江苏科技大学 | Method for realizing IRI value soft measurement based on acceleration sensor |
-
2016
- 2016-06-14 CN CN201610423010.5A patent/CN106404159B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0318739A (en) * | 1989-06-16 | 1991-01-28 | Tech Res Assoc Openair Coal Min Mach | Trouble diagnostic method for speed reducer gear |
| CN102735267A (en) * | 2012-06-20 | 2012-10-17 | 北京航天控制仪器研究所 | Measuring method for inertial measurement device in sled testing |
| CN102937522A (en) * | 2012-08-30 | 2013-02-20 | 桂林电子科技大学 | Composite fault diagnosis method and system of gear case |
| CN103033329A (en) * | 2012-12-24 | 2013-04-10 | 北京工业大学 | Space grid structure frequency identification method of improved power spectrum peak method |
| CN104699980A (en) * | 2015-03-24 | 2015-06-10 | 江苏科技大学 | Method for realizing IRI value soft measurement based on acceleration sensor |
Non-Patent Citations (1)
| Title |
|---|
| 一种惯性测量装置火箭橇试验误差分离方法;刘璠等;《中国惯性技术学报》;20140228;第22卷(第1期);第131-134页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106404159A (en) | 2017-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106404159B (en) | A kind of determining method of Rocket sled test continuous shaking transmitting spectrum | |
| Lin et al. | A new method for measuring engine rotational speed based on the vibration and discrete spectrum correction technique | |
| CN102095885B (en) | Vibration measuring method of engine speed by applying four-point energy centrobaric correction method | |
| CN105426644B (en) | Modal damping recognition methods, device and system | |
| US9798838B2 (en) | Method and system for determining friction coefficient μ for an aircraft landing event | |
| CN110530509B (en) | High-speed motor train unit axle box vibration dominant frequency prediction method based on maximum entropy spectrum analysis | |
| US20090193900A1 (en) | Order tracking method and system | |
| CN104050147B (en) | Time-domain signal is converted into into the method and system of frequency-region signal | |
| CN109116048A (en) | A kind of rotating speed of automobile engine detection and control method | |
| CN111487318B (en) | Time-varying structure instantaneous frequency extraction method | |
| CN104155054A (en) | Rotational inertia frequency domain detection method based on floating torsion pendulum table | |
| CN110398343A (en) | Utilize the model ship drag measurement system of active vibration control technology | |
| CN111780866B (en) | Diesel engine test base installation state natural frequency testing method and device | |
| CN109540482A (en) | A kind of turbine blade is without the mutually high frequency multiplication synchronous vibration parameters analysis method of key and analytical equipment | |
| CN104090126B (en) | A kind of method of testing of accelerometer bandwidth | |
| CN108196313B (en) | The compensation method of shake table guide rail irregularity in a kind of very broadband seismometer dynamic calibration | |
| Rizzoni et al. | Identification of a non-linear internal combustion engine model for on-line indicated torque estimation | |
| CN109635399A (en) | A kind of adding window Integral Transformation method of vibration acceleration signal | |
| CN110133665B (en) | Doppler measurement method of ultrasonic transducer based on envelope curve | |
| CN110094194A (en) | Electric work figure calculates oilwell produced fluid amount method | |
| CN106404322B (en) | A kind of determining method of Rocket sled test vibration transmitting spectrum | |
| CN110514861A (en) | It is a kind of to utilize the autocorrelative high-speed rail speed of service estimation method of amplitude spectrum | |
| CN103913271B (en) | Method for extracting dynamic unbalance signals of rotor at non-stable rotational speed | |
| CN103712794B (en) | Gear train assembly impact shock energy slippages rate detection system and method | |
| CN110395351A (en) | Utilize the model ship resistance measurement method of active vibration control technology |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |