CN101975634A - Engine excitation force measurement method employing window length varying phase difference correction method - Google Patents
Engine excitation force measurement method employing window length varying phase difference correction method Download PDFInfo
- Publication number
- CN101975634A CN101975634A CN 201010288084 CN201010288084A CN101975634A CN 101975634 A CN101975634 A CN 101975634A CN 201010288084 CN201010288084 CN 201010288084 CN 201010288084 A CN201010288084 A CN 201010288084A CN 101975634 A CN101975634 A CN 101975634A
- Authority
- CN
- China
- Prior art keywords
- engine
- centerdot
- frequency
- point
- axis
- 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
Links
Landscapes
- Testing Of Engines (AREA)
Abstract
The invention discloses an engine excitation force measurement method employing a window length varying phase difference correction method. In the method, signal frequency error of an engine excitation frequency interval is solved by utilizing phase difference produced by adding Hanning with different time lengths during signal processing after a vibration acceleration signal on the surface or a suspension point of an engine cylinder is discretely sampled, and the signal frequency error is used for solving actual excitation frequency f, amplitude information and phase information of an engine in a vibration displacement frequency spectrum, so that the vibration displacement vector value in an excitation force measurement formula is accurately obtained and combined with other known engine parameters to calculate the excitation force F (f). The method of the invention can solve the problem that due to the lack of accurate phase information during the engine excitation force measurement, a complex non-linear equation needs to be solved.
Description
Technical field
This aspect relates to power machine and signal Processing field, specifically is a kind of frequency, amplitude, phase method for position of accurate identification of Engine vibration signal, is applied to the identification of engine exciting force.
Background technology
The exciting force of engine mainly is to obtain by the rigid dynamics Model Calculation at present, the excitation that the principal element of considering in calculating is a gas burst power in the cylinder, rotating inertia force and reciprocal inertia force produced.The major defect of this method is need to obtain precise parameters, and the accurate acquisition of some parameter is very difficult, the moment resulting from sidesway that causes as gas burst etc.On the other hand, in actual conditions, engine always links to each other with many annexes, and theoretical calculation formula is not considered the coupling of engine and these annexes usually, this will inevitably cause result of calculation and engineering is actual that certain deviation is arranged, and deviation can be very big in some cases.
In application of practical project, method is by experiment obtained the engine exciting force and is studied by numerous scholars, and this is the problem of typically inverting, and normally adopts the frequency response function matrix method of inverting.But in the vibration control problem of engine, what need is the equivalent force and the equivalent moment at Motor Mass Centre place, owing to can't directly apply excitation at the barycenter place, also the reciprocity method is used in the response that can't directly obtain the barycenter place, therefore can't record transport function, cause the frequency response function matrix method of inverting in the engine exciting force is measured, can't directly use.On the other hand, concerning the exciting force at Motor Mass Centre place was measured, the frequency response function matrix method of inverting was too complicated comparatively speaking.
External scholar J.S.TAO, G.R.LIU and K.Y.LAM etc. has proposed a kind of method (J.S.TAO that measures Motor Mass Centre place exciting force, G.R.LIU and K.Y.LAM is in 2001 articles of delivering on periodical " Journal of Sound and Vibration " " Excitation force identification of an engine with velocity data at mounting points "), point out the amplitude and the phase information of the three-way vibration velocity spectrum at each suspension point place of engine of measuring by extraction, can rebuild the exciting force of engine accurately, the exciting force measure equation
But this method is when extracting phase information, usually can't obtain absolute phase information, therefore proposed only to utilize phase differential to obtain the method for exciting force, and this processing has been introduced and the measurement identical phase variant of counting that suspends, the non-linear overdetermined equation that makes finding the solution of problem become to find the solution a complexity has increased the complicacy of finding the solution.
Summary of the invention
The objective of the invention is defective at existing recognition methods, the engine exciting force measuring method that a kind of application changes the correction method of window appearance potential difference is proposed, solve for want of accurate phase information in the measurement of engine exciting force, and need find the solution the problem of complex nonlinear equation.
The object of the invention is achieved through the following technical solutions:
A kind of application changes the engine exciting force measuring method of window appearance potential difference correction method, comprises the steps:
(1) coordinate system is set up with correlation parameter and is collected: utilize the engine moment inertia experiment table, test out engine quality m, Motor Mass Centre O, with O is that initial point is set up coordinate system O-XYZ, the Y-axis forward points to crankshaft free-end, and Z axle forward is determined the X-axis forward by the right-hand rule straight up; Utilize the engine moment inertia experiment table to test out the moment of inertia J of engine again around the coordinate system X-axis
x, Y-axis moment of inertia J
y, the Z axle moment of inertia J
z, engine is to the product of inertia J of X-axis and Y-axis
Xy, to the inertia J of Y-axis and Z axle
Yz, to the inertia J of Z axle and X-axis
ZxForm the engine quality matrix M;
(2) obtain the multiple stiffness matrix that suspends: with engine by mount supports at threst stand; Determine engine mounting number h, h=3 or h=3; In the coordinate system O-XYZ of engine, give the label that suspends according to the descending order of directions X, be defined as first suspend, second suspend ... h suspends, and with the elastic body test macro the multiple stiffness characteristics that respectively suspends is tested, and obtains the first multiple rigidity that suspends
The second multiple rigidity that suspends
The multiple rigidity that h suspends
Form the multiple stiffness matrix K that suspends
*=[K
* 1K
* 2K
* h];
(3) sensor installation: on the different test points of engine cylinder surface or suspension point, L three-dimensional acceleration transducer is installed, L 〉=3, each test point is installed a three-dimensional acceleration transducer; Sensor connects data acquisition unit, and data acquisition unit connects portable computer; In the coordinate system O-XYZ of engine, determine that according to the descending order of directions X the coordinate of three-dimensional acceleration transducer is respectively [x
1y
1z
1], [x
2y
2z
2] ..., [x
Ly
Lz
L];
(4) on threst stand, setting engine moves with operate as normal rotating speed w, after treating that operating mode is stable, the acquisition time in 10s~30s, the vibration acceleration signal X of X-axis forward, Y-axis forward and Z axle forward by data acquisition unit collection and each test point of synchronous recording
1n, X
2nX
LnY
1n, Y
2nY
LnZ
1n, Z
2nZ
LnX wherein
1nBe the vibration acceleration signal sequence with n point of the first test point X-axis forward collection, X
LnBe n the vibration acceleration signal sequence of putting that L test point X-axis forward is gathered, Y
1nBe n the vibration acceleration signal sequence of putting that the first test point Y-axis forward is gathered, X
LnThe vibration acceleration signal sequence of n the point that L test point Y-axis forward is gathered; Z
1nBe n the vibration acceleration signal sequence of putting that the first test point Z axle forward is gathered, Z
LnThe vibration acceleration signal sequence of n the point that L test point Z axle forward is gathered; Sample frequency f
s, sampling number N represents the data point number that discrete series comprises, wherein n=0,1 ..., N-1, require f
s=N=2
p, P=9 or 10, frequency resolution Δ f=f
s/ N=1;
(6) according to formula
Measure the engine exciting force, wherein:
Exciting force F (f) be mean engine under humorous excitation frequency f of rotating speed w η, the exciting force F of engine X-axis forward
x(f), around the excitation torque M of X-axis
x(f), the exciting force F of Y direction
y(f), around the excitation torque M of Y-axis
y(f), the exciting force F of Z-direction
z(f), around the excitation torque M of Z axle
z(f);
K
*=[K
* 1?K
* 2…K
* h]
Humorous excitation frequency f of η and its corresponding transposed matrix
Determine by the following method:
Humorous the η that selects the engine exciting force to analyze, η=0.5,1,1.5,2,2.5 or 3; Behind humorous time of the η of designated analysis engine speed w, carry out following computing:
Vibration acceleration signal sequence X to the first test point X-axis forward collection with n point
1n, add Hanning window (Hanning) w (t) 0≤t≤T of time length T=1, carry out N point quick Fourier conversion (FFT) according to formula (1-1), obtain amplitude spectrum X
1(f); To X
1nAdd time length
Hanning window (Hanning) w (t) 0≤t≤T carries out N/2 point quick Fourier conversion (FFT) according to formula (1-2), obtains amplitude spectrum X
2(f);
According to engine speed w, determine emending frequency interval [f
1, f
2], wherein
Round numbers;
Round numbers;
At emending frequency interval [f
1, f
2] in, X
1(f) peak-peak respective frequencies f
i, f
iReal part is R
1(f
i), imaginary part is I
1(f
i), phase place is Φ
1(f
i), amplitude is y
1(f
i); X
2(f) peak-peak frequency f
j, phase place is Φ
2(f
j); If emending frequency interval [f
1, f
2] interior peak-peak actual frequency is f
0, f
0Corresponding true phase place
True amplitude A
0(f
0); If f
iFrequency error be
f
iWith actual frequency f
0Difference
For people's formula (1-3), obtain frequency error
The spectral function of Hanning window (Hanning)
Calculate amplitude A after the correction according to formula (1-6)
0(f
0), form actual frequency f
0Acceleration with correspondence
To sequence X
2nX
LnY
1n, Y
2nY
LnZ
1n, Z
2nZ
LnWith the vibration acceleration signal sequence X
1nCarry out same processing, obtain the humorous inferior excitation frequency f of engine speed w η respectively
2xF
Lxf
1y, f
2yF
Lyf
1z, f
2zF
Lz
(1≤H≤L), the engine excitation frequency that records its directions X is f for any sensor H in L the sensor test point
Hx, corresponding acceleration is
The engine excitation frequency that records the Y direction is f
Hy, corresponding acceleration is
The engine excitation frequency that records the Z direction is f
Hz, corresponding acceleration is
The directions X displacement of sensor H is dX
H, the displacement of Y direction is dY
H, the displacement of Z direction is dZ
H,
Obtain the transposed matrix of sensor H thus:
For realizing that further the object of the invention, described crankshaft free-end are the end that bent axle connects belt pulley.
Described engine is defined as suspension point with the contact point that suspends that supports it.
Described operate as normal rotating speed w is preferably 750r/min to 5500r/min; Described operating mode is stable to be meant that the engine speed fluctuation is in 20r/min.
Described engine moment inertia experiment table preferably adopts three string pendulum method of testing experiment tablees.
The present invention only needs the value by conversion η, η=0.5,1,1.5,2,2.5,3 to the exciting force test of different humorous the η of engine; Calculate different frequency refinement interval [f
1, f
2], recalculate this excitation frequency f and corresponding transposed matrix of humorous time
Solve humorous the exciting force F of η (f) of engine speed w.
With respect to prior art, the present invention has following advantage:
The present invention is by after sampling to the vibration displacement signal discrete on engine cylinder surface or the suspension point, the phase differential that adds the Hanning window (Hanning) of different time length and produce, utilize precise frequency, amplitude and the phase information of each response point vibration signal of this phase difference correction, thereby improve the measuring accuracy of exciting force.The present invention is in conjunction with J.S.TAO, the exciting force measure equation that G.R.LIU and K.Y.LAM proposes
Middle K
*, M,
Method for solving proposes to use the correction method of change window appearance potential difference and accurately extracts transposed matrix
In frequency, amplitude and phase information, improved and utilized the phase differential method to calculate
Method, avoided finding the solution of problem having been become the non-linear overdetermined equation of finding the solution a complexity because of utilizing the phase differential method, improved and found the solution
Precision and efficient, and improved the simplicity that F (f) finds the solution.
Embodiment:
The invention will be further described below in conjunction with embodiment, need to prove, embodiment does not constitute the qualification to the claimed scope of the present invention.
Embodiment 1
(1) utilize three-way pendulum engine rotary inertia experimental bench to test out: certain 1.8L in-line four cylinder engine quality m, Motor Mass Centre O is that initial point is set up coordinate system O-XYZ with barycenter O; And then test out the moment of inertia J of engine around the coordinate system X-axis
x, Y-axis moment of inertia J
y, the Z axle moment of inertia J
z, engine is to the product of inertia J of X-axis and Y-axis
Xy, to the inertia J of Y-axis and Z axle
Yz, to the inertia J of Z axle and X-axis
Zx, form the engine quality matrix M;
(2) determine engine mounting number h, h=3; In the coordinate system O-XYZ of engine, give the label that suspends according to the descending order of directions X, be defined as and suspend 1, suspend 2 ... the h that suspends adopts German MTS831 type elastic body test system and test to go out the multiple rigidity that respectively suspends, and forms the multiple stiffness matrix K that suspends
*=[k
* 1K
* 2K
* 3]=[be N/m (70000+7000j) N/m (80000+8000j) N/m (75000+7500j)];
(3) 1.8L in-line four cylinder engine is installed on the German FEV threst stand, adopts CPB three-dimensional acceleration transducer, number L=3 is arranged in the suspension point place.Measure the coordinate of 3 sensor test points by three-coordinates measuring machine, form the test point coordinate matrix
(4) input of CPB three-dimensional acceleration transducer is the data acquisition system (DAS) of German Miller shellfish nurse (BBM) testing apparatus, and data acquisition signal input portable computer is by the PAK test analysis software of installing on the portable computer, real-time monitored and tracer signal;
(5) on threst stand, set engine with certain operate as normal rotating speed w=1815rpm operation, treat the stable back of operating mode (being that the fluctuation of speed is in 20r/min), the beginning image data, the data acquisition system (DAS) of German Miller shellfish nurse (BBM) testing apparatus imported signal by CPB three-dimensional acceleration transducer, by PAK test analysis software, real-time monitored and the tracer signal of installing on the portable computer.Acquisition time length is 30s; Sample frequency f
s=512, sampling number N=512, then frequency resolution Δ f=f
s/ N=1, the vibration acceleration signal X of X-axis forward, Y-axis forward and Z axle forward by data acquisition unit collection and each test point of synchronous recording
1n, X
2n, X
3nY
1n, Y
2n, Y
3nZ
1n, Z
2n, Z
3nWherein n=0,1 ..., N-1; To vibration acceleration signal X
1n, X
2n, X
3nY
1n, Y
2n, Y
3nZ
1n, Z
2n, Z
3nCarry out the conversion of N=512 point fast Fourier.
According to engine speed w=1815rpm, determine frequency separation [f
1, f
2], the refinement interval of 0.5 humorous excitation frequency f of engine speed w is f
1=10HZ and f
2=20HZ; 1 humorous time refinement interval is f
1=25HZ and f
2=35HZ; 1.5 humorous time refinement interval is f
1=40HZ and f
2=50HZ; 2 humorous times refinement interval is f
1=55HZ and f
2=65HZ; 2.5 humorous time refinement interval is f
1=71HZ and f
2=81HZ; 3 humorous times refinement interval is f
1=86HZ and f
2=96HZ; To each frequency separation [f
1, f
2], to use to change window appearance potential difference method, concrete steps are carried out according to formula (1-1) to (1-7).
3), only list 1 humorous time and 2 humorous times the excitation frequency f of engine speed w here, the 1 humorous secondary frequencies interval [f of engine speed w owing to the length relation
1, f
2]=[25,35], use to change window appearance potential difference method, obtain 3 sensor correspondences the excitation frequency f and the acceleration of test point separately
As shown in table 1; The 2 humorous secondary frequencies interval [f of engine speed w
1, f
2]=[55,65], use to change window appearance potential difference method, obtain 3 sensor correspondences the excitation frequency f and the acceleration of test point separately
As shown in table 2;
1 humorous secondary frequencies and the acceleration of table 1 engine speed w
2 humorous secondary frequencies and the accelerations of table 2 engine speed w
Transposed matrix
Wherein:
(7) fundamental formular of measuring according to the engine exciting force
Wherein:
For the mass matrix M among the people (1), the multiple stiffness matrix K that suspends in (2)
*, the measuring point coordinates matrix in (3)
When asking 1 humorous exciting force of engine speed w, for people f=30.25HZ and
1 humorous exciting force obtaining engine start machine rotating speed w is as shown in table 3; When asking 2 humorous exciting forces of engine speed w, for people f=60.5HZ and
2 humorous exciting forces obtaining engine start machine rotating speed w are as shown in table 4;
1 humorous exciting force of table 3 engine speed w
2 humorous exciting forces of table 4 engine speed w
Embodiment 1 has illustrated that the present invention is by after sampling to the vibration displacement signal discrete on engine cylinder surface or the suspension point, the phase differential that adds the Hanning window (Hanning) of different time length and produce, utilize precise frequency, amplitude and the phase information of each response point vibration signal of this phase difference correction, thereby improve the measuring accuracy of exciting force.
Claims (5)
1. use the engine exciting force measuring method that changes the correction method of window appearance potential difference, it is characterized in that comprising the steps:
(1) coordinate system is set up with correlation parameter and is collected: utilize the engine moment inertia experiment table, test out engine quality m, Motor Mass Centre O, with O is that initial point is set up coordinate system O-XYZ, the Y-axis forward points to crankshaft free-end, and Z axle forward is determined the X-axis forward by the right-hand rule straight up; Utilize the engine moment inertia experiment table to test out the moment of inertia J of engine again around the coordinate system X-axis
x, Y-axis moment of inertia J
y, the Z axle moment of inertia J
z, engine is to the product of inertia J of X-axis and Y-axis
Xy, to the inertia J of Y-axis and Z axle
Yz, to the inertia J of Z axle and X-axis
Zx, form the engine quality matrix M;
(2) obtain the multiple stiffness matrix that suspends: with engine by mount supports at threst stand; Determine engine mounting number h, h=3 or h=3; In the coordinate system O-XYZ of engine, give the label that suspends according to the descending order of directions X, be defined as first suspend, second suspend ... h suspends, and with the elastic body test macro the multiple stiffness characteristics that respectively suspends is tested, and obtains the first multiple rigidity that suspends
The second multiple rigidity that suspends
The multiple rigidity that h suspends
Form the multiple stiffness matrix K that suspends
*=[K
* 1K
* 2K
* h];
(3) sensor installation: on the different test points of engine cylinder surface or suspension point, L three-dimensional acceleration transducer is installed, L 〉=3, each test point is installed a three-dimensional acceleration transducer; Sensor connects data acquisition unit, and data acquisition unit connects portable computer; In the coordinate system O-XYZ of engine, determine that according to the descending order of directions X the coordinate of three-dimensional acceleration transducer is respectively [x
1y
1z
1], [x
2y
2z
2] ..., [x
Ly
Lz
L];
(4) on threst stand, setting engine moves with operate as normal rotating speed w, after treating that operating mode is stable, the acquisition time in 10s~30s, the vibration acceleration signal X of X-axis forward, Y-axis forward and Z axle forward by data acquisition unit collection and each test point of synchronous recording
1n, X
2nX
LnY
1n, Y
2nY
LnZ
1n, Z
2nZ
LnX wherein
1nBe the vibration acceleration signal sequence with n point of the first test point X-axis forward collection, X
LnBe n the vibration acceleration signal sequence of putting that L test point X-axis forward is gathered, Y
1nBe n the vibration acceleration signal sequence of putting that the first test point Y-axis forward is gathered, X
LnThe vibration acceleration signal sequence of n the point that L test point Y-axis forward is gathered; Z
1nBe n the vibration acceleration signal sequence of putting that the first test point Z axle forward is gathered, Z
LnThe vibration acceleration signal sequence of n the point that L test point Z axle forward is gathered; Sample frequency f
s, sampling number N represents the data point number that discrete series comprises, wherein n=0,1 ..., N-1, require f
s=N=2
p, P=9 or 10, frequency resolution Δ f=f
s/ N=1;
(5) according to formula
Measure the engine exciting force, wherein:
Exciting force F (f) be mean engine under humorous excitation frequency f of rotating speed w η, the exciting force F of engine X-axis forward
x(f), around the excitation torque M of X-axis
x(f), the exciting force F of Y direction
y(f), around the excitation torque M of Y-axis
y(f), the exciting force F of Z-direction
z(f), around the excitation torque M of Z axle
z(f);
K
*=[K
* 1?K
* 2?…K
* h]
Humorous excitation frequency f of η and its corresponding transposed matrix
Determine by the following method:
Humorous the η that selects the engine exciting force to analyze, η=0.5,1,1.5,2,2.5 or 3; Behind humorous time of the η of designated analysis engine speed w, carry out following computing:
Vibration acceleration signal sequence X to the first test point X-axis forward collection with n point
1n, add Hanning window w (t) 0≤t≤T of time length T=1, carry out the conversion of N point quick Fourier according to formula (1-1), obtain amplitude spectrum X
1(f); To X
1nAdd time length
Hanning window w (t) 0≤t≤T carries out N/2 point quick Fourier conversion (FFT) according to formula (1-2), obtains amplitude spectrum X
2(f);
According to engine speed w, determine emending frequency interval [f
1, f
2], wherein
Round numbers;
Round numbers;
At emending frequency interval [f
1, f
2] in, X
1(f) peak-peak respective frequencies f
i, f
iReal part is R
1(f
i), imaginary part is I
1(f
i), phase place is Φ
1(f
i), amplitude is y
1(f
i); X
2(f) peak-peak frequency f
j, phase place is Φ
2(f
j); If emending frequency interval [f
1, f
2] interior peak-peak actual frequency is f
0, f
0Corresponding true phase place
True amplitude A
0(f
0); If f
iFrequency error be
f
iWith actual frequency f
0Difference
For people's formula (1-3), obtain frequency error
The spectral function of Hanning window
Calculate amplitude A after the correction according to formula (1-6)
0(f
0), form actual frequency f
0Acceleration with correspondence
To sequence X
2nX
LnY
1n, Y
2nY
LnZ
1n, Z
2nZ
LnWith the vibration acceleration signal sequence X
1nCarry out same processing, obtain the humorous inferior excitation frequency f of engine speed w η respectively
2xF
Lxf
1y, f
2yF
Lyf
1z, f
2zF
Lz
(1≤H≤L), the engine excitation frequency that records its directions X is f for any sensor H in L the sensor test point
Hx, corresponding acceleration is
The engine excitation frequency that records the Y direction is f
Hy, corresponding acceleration is
The engine excitation frequency that records the Z direction is f
Hz, corresponding acceleration is
The directions X displacement of sensor H is dX
H, the displacement of Y direction is dY
H, the displacement of Z direction is dZ
H, obtain the transposed matrix of sensor H thus:
2. application according to claim 1 changes the engine exciting force measuring method of window appearance potential difference correction method, it is characterized in that: described crankshaft free-end is the end that bent axle connects belt pulley.
3. application according to claim 1 changes the engine exciting force measuring method of window appearance potential difference correction method, it is characterized in that: described engine is defined as suspension point with the contact point that suspends that supports it.
4. application according to claim 1 changes the engine exciting force measuring method of window appearance potential difference correction method, and it is characterized in that: described operate as normal rotating speed w is 750r/min to 5500r/min; Described operating mode is stable to be meant that the engine speed fluctuation is in 20r/min.
5. application according to claim 1 changes the engine exciting force measuring method of window appearance potential difference correction method, and it is characterized in that: described engine moment inertia experiment table adopts three string pendulum method of testing experiment tablees.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010288084 CN101975634B (en) | 2010-09-17 | 2010-09-17 | Engine excitation force measurement method employing window length varying phase difference correction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010288084 CN101975634B (en) | 2010-09-17 | 2010-09-17 | Engine excitation force measurement method employing window length varying phase difference correction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101975634A true CN101975634A (en) | 2011-02-16 |
CN101975634B CN101975634B (en) | 2012-12-05 |
Family
ID=43575537
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010288084 Expired - Fee Related CN101975634B (en) | 2010-09-17 | 2010-09-17 | Engine excitation force measurement method employing window length varying phase difference correction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101975634B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680270A (en) * | 2012-06-07 | 2012-09-19 | 哈尔滨工程大学 | Indirect measurement method for ship excitation load caused by equipment under combined effect of unbalanced excitation force and unbalanced bending moment |
CN105403344A (en) * | 2015-12-16 | 2016-03-16 | 浙江大学 | Pipeline real-time stress obtaining method |
CN106855119A (en) * | 2015-12-08 | 2017-06-16 | 通用汽车环球科技运作有限责任公司 | Control using Fast Fourier Transform (FFT) to buncher |
CN106970264A (en) * | 2017-03-02 | 2017-07-21 | 浙江大学 | A kind of improvement phase difference correction method for considering mains frequency rate of change |
CN109507495A (en) * | 2018-10-17 | 2019-03-22 | 华北水利水电大学 | It is a kind of to become the long quasi- simultaneous interconnecting measurement method of parameters of window |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998041835A1 (en) * | 1997-03-17 | 1998-09-24 | Hitachi, Ltd. | Vibration exciting apparatus and vibration testing apparatus for structure using same |
CN201408105Y (en) * | 2009-05-15 | 2010-02-17 | 北京工业大学 | Time-sharing fast steady-state sine sweep frequency shock excitation frequency response function measuring system |
-
2010
- 2010-09-17 CN CN 201010288084 patent/CN101975634B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998041835A1 (en) * | 1997-03-17 | 1998-09-24 | Hitachi, Ltd. | Vibration exciting apparatus and vibration testing apparatus for structure using same |
CN201408105Y (en) * | 2009-05-15 | 2010-02-17 | 北京工业大学 | Time-sharing fast steady-state sine sweep frequency shock excitation frequency response function measuring system |
Non-Patent Citations (3)
Title |
---|
《Journal of Sound and Vibration》 20011230 J.S TAO,G.R.LIU,K.Y.LAM EXCITATION FORCE IDENTIFICATION OF AN ENGINE WITH VELOCITY DATA AT MOUNTING POINTS 第242卷, 第2期 2 * |
《华南理工大学学报》 20071031 丁康,杨志坚 改变窗长相位差校正法的改进与抗噪性能分析 第35卷, 第10期 2 * |
《振动工程学报》 20091230 杨志坚,丁康 调制FFT及其在离散频谱校正技术中的应用 第23卷, 第6期 2 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680270A (en) * | 2012-06-07 | 2012-09-19 | 哈尔滨工程大学 | Indirect measurement method for ship excitation load caused by equipment under combined effect of unbalanced excitation force and unbalanced bending moment |
CN102680270B (en) * | 2012-06-07 | 2014-08-06 | 哈尔滨工程大学 | Indirect measurement method for ship excitation load caused by equipment under combined effect of unbalanced excitation force and unbalanced bending moment |
CN106855119A (en) * | 2015-12-08 | 2017-06-16 | 通用汽车环球科技运作有限责任公司 | Control using Fast Fourier Transform (FFT) to buncher |
US10082206B2 (en) | 2015-12-08 | 2018-09-25 | GM Global Technology Operations LLC | Control of a continuously variable transmission using fast fourier transformation |
CN105403344A (en) * | 2015-12-16 | 2016-03-16 | 浙江大学 | Pipeline real-time stress obtaining method |
CN105403344B (en) * | 2015-12-16 | 2018-03-16 | 浙江大学 | The acquisition methods of the real-time stress of pipeline |
CN106970264A (en) * | 2017-03-02 | 2017-07-21 | 浙江大学 | A kind of improvement phase difference correction method for considering mains frequency rate of change |
CN106970264B (en) * | 2017-03-02 | 2020-02-21 | 浙江大学 | Improved phase difference correction method considering power grid frequency change rate |
CN109507495A (en) * | 2018-10-17 | 2019-03-22 | 华北水利水电大学 | It is a kind of to become the long quasi- simultaneous interconnecting measurement method of parameters of window |
CN109507495B (en) * | 2018-10-17 | 2020-12-15 | 华北水利水电大学 | Variable-window-length quasi-synchronization grid-connected parameter measurement method |
Also Published As
Publication number | Publication date |
---|---|
CN101975634B (en) | 2012-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101975633B (en) | Method for measuring energizing force of engine by continuous refinement analytical Fourier transform method | |
CN105424160A (en) | Method for realizing blade synchronous vibration parameter identification | |
CN101975634B (en) | Engine excitation force measurement method employing window length varying phase difference correction method | |
CN102183363B (en) | Dynamic characteristic parameter test device and method of sliding guide rail junction surface | |
Xiang et al. | Torsional vibration measurements on rotating shaft system using laser doppler vibrometer | |
CN103175602B (en) | Modal testing system and modal testing method on basis of single-point laser continuous plane-scanning vibration measurement | |
Di Maio et al. | Continuous Scan, a method for performing modal testing using meaningful measurement parameters; Part I | |
Koene et al. | Internet of things based monitoring of large rotor vibration with a microelectromechanical systems accelerometer | |
Cheng et al. | Envelope deformation in computed order tracking and error in order analysis | |
CN102661782B (en) | Rapid measurement method and device of super-low-frequency vibration parameter quantity value | |
CN110108430B (en) | High-precision damage positioning method for fiber reinforced composite beam | |
JPS60214228A (en) | Vibration analyzing method | |
Sun et al. | Modal identification of vibrating structures using singular value decomposition and nonlinear iteration based on high-speed digital image correlation | |
CN102878989A (en) | Triaxial angular vibration measuring method through adopting satellite-borne linear accelerometers | |
CN104236588A (en) | Determination method of directions and frequency cracking of stiff shafts of cylindrical vibratory gyro harmonic oscillator | |
CN107064559A (en) | A kind of SINS accelerometer frequency characteristic test methods based on angle oscillating motion | |
CN111537170A (en) | Dynamic stiffness testing method for servo actuator | |
RU2013113939A (en) | METHOD FOR RECONSTRUCTION OF THE THREE-DIMENSIONAL PHYSICAL STATE OF THE MONITORING OBJECT IN THE MEASURING POINT | |
CN111487318A (en) | Time-varying structure instantaneous frequency extraction method | |
Abd_Elsalam et al. | Modal analysis on tire with respect to different parameters | |
CN101368870B (en) | Preparation method of amplitude frequency spectrum used for mechanical rotor single cross section shaft vibration analysis | |
CN102322942B (en) | Vibration signal vector ellipse spectrum analyzing method of rotating machinery | |
Kaynardag et al. | Identification of bending modes of vibration in rails by a laser Doppler vibrometer on a moving platform | |
CN112327957A (en) | Method and device for controlling low-frequency vibration multi-order line spectrum of rotor | |
Badri et al. | A method to calibrate the measured responses by MEMS accelerometers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121205 Termination date: 20150917 |
|
EXPY | Termination of patent right or utility model |