CN108910079A - A kind of rotor model.test system dynamic balancing adjustment phase recognition methods - Google Patents
A kind of rotor model.test system dynamic balancing adjustment phase recognition methods Download PDFInfo
- Publication number
- CN108910079A CN108910079A CN201810668288.8A CN201810668288A CN108910079A CN 108910079 A CN108910079 A CN 108910079A CN 201810668288 A CN201810668288 A CN 201810668288A CN 108910079 A CN108910079 A CN 108910079A
- Authority
- CN
- China
- Prior art keywords
- array
- data
- phase
- dynamic balancing
- rotor
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
A kind of rotor model.test system dynamic balancing adjustment phase recognition methods of the present invention, belong to helicopter model Rotor Test technical field, include the following steps:Step 1: according to rotor blade the piece number n, test speed Speed, distinguishable minimum phase angle P0 and percentage error revolving speed Rt, calculate sample frequency Fn as follows:Step 2: under the liter rotary regimes of rotor platform tach signal and vibration signal synchronize acquisition;Step 3: carrying out location finding to revolving speed peak value present in rotary speed data array;Step 4: carrying out the interception of vibration signal complete cycle using the peak position that step 3 obtains;Step 5: to vibration signal order analysis complete cycle of interception;Step 6: carrying out identification judgement to rotor model testing stand dynamic balancing adjustment phase.The present invention can quickly obtain the corresponding phase of high-precision reference frequency, at the same can the rotor model test dynamic balancing adjustments of adaptive various different rotating speeds need.
Description
Technical field
The invention belongs to helicopter model Rotor Test technical fields, and in particular to a kind of rotor model.test system dynamic balancing adjustment
Phase recognition methods.
Background technique
Rig for model test dynamic balancing adjustment is the work must carry out before model rotor test, only by testing stand
Dynamic balancing adjustment just can be carried out subsequent section's purpose test to perfect condition.
It is all to carry out spectrum analysis using acquisition set time length vibration signal to obtain dynamic balancing adjustment both at home and abroad at present
The mode of required phase carries out rig for model test dynamic balancing adjustment, and this method needs long data acquisition time, while by
The trueness error for the reference frequency that the limitation of sampling time length is calculated is larger, causes to correspond to phase to the reference frequency of selection
Position and actual phase value deviation are larger, so that testing stand dynamic balancing adjustment efficiency is lower.
Summary of the invention
The purpose of the present invention:To solve the above-mentioned problems, the invention proposes a kind of rotor model.test system dynamic balancing adjustment phases
Position recognition methods, by carrying out interception complete cycle to the rotor model testing stand vibration signal that collects, then to intercepting
The data in the multiple periods arrived carry out order analysis, can quickly obtain the corresponding phase of high-precision reference frequency, this method is not
Only can quickly obtain needing the dynamic balancing phase that is adjusted, at the same can adaptive various different rotating speeds rotor model test
Dynamic balancing adjustment needs.
Technical solution of the present invention:A kind of rotor model.test system dynamic balancing adjustment phase recognition methods are suitable for rank complete cycle
The rotor model of secondary analysis is tested, and is included the following steps:
Step 1: according to rotor blade the piece number n, test speed Speed, distinguishable minimum phase angle P0 and percentage
Error revolving speed Rt calculates sample frequency Fn as follows:
Fn=max (int (Speed × n × 6/ (P0 × 1000))+1) × 1000, (int (1250/ (Rt × 30000))+
1)×1000}
Wherein, int is downward round numbers, and max is to take maximum number from two numbers in bracket;
Step 2: under the liter rotary regimes of rotor platform tach signal and vibration signal synchronize acquisition;
The tach signal of rotor platform and vibration signal are accessed simultaneously in the data collection system with synchronous acquisition function,
Sample frequency Fn is provided in the data collection system;
By test speed Speed as defined in test bench running to test mission book, then carry out always away from be incremented by be manipulated to
Carry out dynamic balancing adjustment always away from;
It is to be tested it is in stable condition after, based on data collection system setting sample frequency Fn obtain rotary speed data array
[Spa1, Spa2 ..., Spax] and vibration data number [Vi1, Vi2 ..., Vix];
Wherein, X is number of sampling points;
Step 3: carrying out location finding to revolving speed peak value present in rotary speed data array;
A), the numerical value in the rotary speed data array collected is taken absolute value;
B), be arranged a revolving speed peak threshold Psp, using threshold value Psp to rotary speed data array [Spa1, Spa2 ...,
Spax] carry out peak value searching;
C), after having searched for rotary speed data, one group of peak position array [La1, La2 ..., Lay] can be obtained;
Wherein y is the circle number of peak position number or revolving speed and vibration signal;
Step 4: carrying out the interception of vibration signal complete cycle using the peak position that step 3 obtains;
D), using peak position array [La1, La2 ..., Lay] be calculated y-1 peak value starting point array [La1,
La2 ..., La (y-1)] and corresponding data length number [La2-La1, La3-La2 ..., Lay-La (y-1)];
E), calculated peak value starting point array and data length number group, by vibration data array [Vi1,
Vi2 ..., Vix] interception at the whole circle vibration data array of y-1 [Vi11, Vi12 ... Vi0 (La2-La1)], [Vi21,
Vi22 ... Vi2 (La2-La1)] ... [2 ... Vi (y-2) of Vi (y-1) 1, Vi (y-2) (Lay-La (y-1))], and calculate separately
The corresponding reference frequency Fb of the whole circle vibration data array of y-1:
Fb (y-1)=Fn/ (Lay-La (y-1))
Step 5: to vibration signal order analysis complete cycle of interception;
Using be calculated reference frequency Fb0, Fb1 ... Fb (y-1) is respectively to its corresponding y-1 whole circle vibration datas
Array carries out FFT spectrum analysis, to obtain the order amplitude number and order phase data of y-1 circle;
Step 6: carrying out identification judgement to rotor model testing stand dynamic balancing adjustment phase;
Choose the y-1 amplitude composition amplitude array and y-1 phase composition phase array of the first rank;
It is ranked up to obtain the sequence amplitude array [A after one group of sequence by sizen0,An1,…An9] and the sorting position
Corresponding sequence phase array [Pn0,Pn1,…Pn(n-1)];
Intermediate n-4 value is taken averagely to be obtained average amplitude At=(A from sequence amplitude arrayn2+An3+An4+An5+…
+A(n-3)+A(n-2))/(n-4),
Take phase array [Pn0,Pn1,…Pn(n-1)] intermediate n-4 value averagely obtained average phase Pt=(Pn2+Pn3+
Pn4+Pn5+…+P(n-3)+P(n-2))/(n-4);
If average amplitude At is greater than the vibration values of setting, need to carry out dynamic balancing adjustment, and needs to increase weight stack
Phase position is Pt+180.
Preferably, in the step 3, first by rotary speed data array [Spa1, Spa2 ..., Spax] from current search position
It sets and takes out odd number t (t >=3) a point,
This t point if (t-1)/2 all bigger than threshold value Psp and preceding point is gradually increased, at the same rear (t-1)/2 point by
Decrescence small, then the position of (t-1)/2+1 point in this section of rotary speed data had both been a peak point;
Next records position of the peak point in rotary speed data array [Spa1, Spa2 ..., Spax], while will work as
Preceding searching position skips s (s >=t) a data point backward, repeats step 1;Otherwise searching position is moved 1 from current location backward
A position, repeats the above steps.
Preferably, in the step 2, it is to be tested it is in stable condition after, based on data collection system setting sample frequency Fn
Obtain the rotary speed data array [Spa1, Spa2 ..., Spax] and vibration data number [Vi1, Vi2 ..., Vix] of at least 10 circles.
The beneficial effect of technical solution of the present invention:The present invention passes through to the rotor model testing stand vibration signal collected
Complete cycle (i.e. whole circle) interception is carried out, the data in the multiple periods then obtained to interception carry out order analysis, can quickly obtain
The corresponding phase of high-precision reference frequency, i.e. dynamic balancing adjustment phase.This method not only can quickly obtain needing to be adjusted
Dynamic balancing phase, while can the rotor model test dynamic balancing adjustments of adaptive various different rotating speeds need.
Detailed description of the invention
Fig. 1 is the process signal of a preferred embodiment of rotor model.test system dynamic balancing adjustment phase recognition methods of the present invention
Figure.
Specific embodiment
To keep the purposes, technical schemes and advantages of the invention implemented clearer, below in conjunction in the embodiment of the present invention
Attached drawing, technical solution in the embodiment of the present invention is further described in more detail.In the accompanying drawings, identical from beginning to end or class
As label indicate same or similar element or element with the same or similar functions.Described embodiment is the present invention
A part of the embodiment, instead of all the embodiments.The embodiments described below with reference to the accompanying drawings are exemplary, it is intended to use
It is of the invention in explaining, and be not considered as limiting the invention.Based on the embodiments of the present invention, ordinary skill people
Member's every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.Under
Face is described in detail the embodiment of the present invention in conjunction with attached drawing.
As shown in Figure 1, a kind of rotor model.test system dynamic balancing adjustment phase recognition methods of the present invention, are suitable for order complete cycle
The rotor model of analysis is tested, and is included the following steps:
1) sample frequency calculates
According to rotor blade the piece number n (unit:Piece), test speed Speed (unit:Rev/min), distinguishable minimum phase
Parallactic angle P0 (unit:Degree) and percentage error revolving speed Rt (unit:Percentage %), sample frequency Fn is calculated as follows
(unit:Sampling number/second):
Fn=max (int (Speed × n × 6/ (P0 × 1000))+1) × 1000, (int (1250/ (Rt × 30000))+
1)×1000} (1)
Wherein, int () is downward round numbers, and max { } is to take maximum number from two numbers in bracket.
2) revolving speed and Vibration Synchronization acquisition
Tach signal and vibration signal are accessed simultaneously in the data collection system with synchronous acquisition function, be powered preheating
30 minutes.By test speed Speed as defined in test bench running to test mission book, then carry out always away from be incremented by be manipulated to into
Action balance adjustment always away from.Wait 10 seconds it is to be tested it is in stable condition after, by the result calculated in (1), sample frequency is set,
Then acquire at least n circle (n >=10) rotary speed data array [Spa1, Spa2 ..., Spax] and vibration data array [Vi1,
Vi2 ..., Vix], x is the number of sampled point.
3) revolving speed peak position is searched for
By the numerical value of the collected n rotary speed data enclosed take absolute value to obtain new absolute value rotary speed data array [Spb1,
Spb2,…,Spbx];
One revolving speed peak threshold Psp is set, using threshold value Psp to rotary speed data array [Spa1, Spa2 ...,
Spax] peak value searching is carried out, search step is as follows:
Step 1, by rotary speed data array [Spa1, Spa2 ..., Spax] from current search position (initiating searches first
A data point Spa1) take out odd number t (t >=3) a point, if this t point (t-1)/2 all bigger than threshold value Psp and preceding point by
It is cumulative big, while rear (t-1)/2 point is gradually reduced, then the position of (t-1)/2+1 point in this section of rotary speed data was both
One peak point;
Step 2, if finding peak point, record the peak point rotary speed data array [Spa1, Spa2 ...,
Spax] in position, while s (s >=t) a data point is skipped into current search position backward, repeats step 1;It otherwise will search
1 position is moved from current location in position backward, repeats step 1;
Step 3 if the peak position number found is greater than n, is needed revolving speed after having searched for rotary speed data
The numerical value of peak value threshold values increases (new revolving speed peak value threshold values=original revolving speed peak value threshold values × 110%), then by current search position
It installs and is set to first data point Spa1, repeat step 1 and arrive step 3;If the peak position number found is less than n,
It needs the numerical value of revolving speed peak value threshold values reducing (new revolving speed peak value threshold values=original revolving speed peak value threshold values × 90%), then will work as
Preceding searching position is set as first data point Spa1, repeats step 1 and arrives step 3;If the peak position number found is equal to n
When a, then terminate this search and peak position array [La0, La1 ..., La10].
4) vibration signal complete cycle intercepts
Using peak position array [La0, La1 ..., La10] be calculated 10 peak value starting point arrays [La0, La2 ...,
La9] and corresponding data length array [La1-La0, La2-La1 ..., La10-La9];
Calculated peak value starting point array and data length number group, by vibration data array [Vi1, Vi2 ...,
Vix] interception at the whole circle vibration data array of n-1 [Vi01, Vi02 ... Vi0 (La1-La0)], [Vi11, Vi12 ... Vi1
(La2-La1)],…[Vi91,Vi(n-1)2,…Vi(n-1)(La10-La9)].It calculates separately to obtain n-1 whole circle vibrations simultaneously
The dynamic corresponding reference frequency Fb0, Fb1 of array of data ... Fb (n-1), calculation formula is as follows:
Fb0=Fn/ (La1-La0), Fb2=Fn/ (La2-La1) ... Fb (n-1)=Fn/ (Lan-La (n-1))
5) vibration signal complete cycle order analysis
Using be calculated reference frequency Fb0, Fb1 ... Fb (n-1) is respectively to its corresponding n-1 whole circle vibration datas
Array carries out FFT spectrum analysis.Since data are in the data for intercepting into whole circle before, spectrum analysis here is exactly order
Analysis, to obtain the order amplitude array [Am01, Am02 ... Am0 ((La1-La0)/2)] of n-1 circle, [Am11, Am12 ...
Am1 ((La2-La1)/2)] ... [Am (n-1) 1, Am (n-1) 2 ... Am (n-1) ((Lan-La (n-1))/2)] and order phase
Data [Ph01, Ph02 ... Ph0 ((La1-La0)/2)], [Ph11, Ph12 ... Ph1 ((La2-La1)/2)] ... [Ph (n-1)
1,Ph(n-1)2,…Ph(n-1)((Lan-La(n-1))/2)]。
6) dynamic balancing adjustment phase identification
Extract 10 amplitude groups of the 1st rank respectively from obtained n-1 circle order amplitude array and order phase data
At amplitude array [Am02, Am12 ... Am (n-1) 2] and n-1 phase composition phase array [Ph02, Ph12 ... Ph (n-1)
2].Amplitude array [Am02, Am12 ... Am (n-1) 2] is ranked up to obtain the sequence amplitude array after one group of sequence by size
[An0, An1 ... An9] and the corresponding sequence phase array of the sorting position [Pn0, Pn1 ... Pn (n-1)], then from sequence
Intermediate n-4 value is taken averagely to be obtained average amplitude At=(An2+An3+An4+An5+ ...+A (n-3)+A in amplitude array
(n-2))/(n-4), while taking phase array [Pn0, Pn1 ... Pn (n-1)]
Intermediate n-4 value is averagely obtained average phase Pt=(Pn2+Pn3+Pn4+Pn5+ ...+P (n-3)+P (n-
2))/(n-4).If average amplitude At is greater than the vibration values of setting, need to carry out dynamic balancing adjustment, and needs to increase counterweight
The phase position of piece is Pt+180.
Specific example is as follows:
Sample frequency calculates:
Rotor blade the piece number n=4, test speed Speed=960, distinguishable minimum phase angle P0=5 and percentage
Error revolving speed Rt=0.2% calculates sample frequency Fn by above-mentioned formula (1):
Fn=max { 5000,21000 }=21000
Revolving speed and Vibration Synchronization acquisition:
Tach signal and vibration signal are accessed simultaneously in the data collection system with synchronous acquisition function, be powered preheating
30 minutes.Then test speed Speed=960 as defined in test bench running to test mission book is carried out always away from incremental manipulation
To pending dynamic balancing adjustment always away from 7 °.Wait 10 seconds it is to be tested it is in stable condition after, by calculating sample frequency setting adopt
Sample frequency, then acquire 11 circle rotary speed data array [Spa1, Spa2 ..., Spax] and vibration data array [Vi1,
Vi2 ..., Vix], x is the number of sampled point.
The search of revolving speed peak position:
The numerical value for the rotary speed data that collected 11 enclose is taken absolute value to obtain new absolute value rotary speed data array
[Spb1,Spb2,…,Spbx];One revolving speed peak threshold Psp=1.5 is set, using threshold value Psp to rotary speed data array
[Spa1, Spa2 ..., Spax] progress peak value searching, the peak position array searched for [37,1348,2661,3971,
5284,6594,7907,9218,10529]。
The interception of vibration signal complete cycle:
It is calculated using peak position array [37,1348,2661,3971,5284,6594,7907,9218,10529]
10 peak value starting point arrays [37,1348,2661,3971,5284,6594,7907,9218] and corresponding data length array
[1311,1313,1310,1313,1310,1313,1311,1311];
Calculated peak value starting point array and data length number group, by vibration data array [Vi1, Vi2 ...,
Vix] interception at 10 whole circle vibration data arrays (i.e. 10 enclose) [Vi01, Vi02 ... Vi0 (La1-La0)], [Vi11,
Vi12,…Vi1(La2-La1)],…[Vi91,Vi92,…Vi9(La10-La9)].It calculates separately to obtain 10 whole circle vibrations simultaneously
The dynamic corresponding reference frequency (unit of array of data:Hertz):[16.0134,15.9932,16.0285,15.9832,16.0293,
15.9917、16.0156、16.0212]。
Complete cycle vibration signal order analysis:
FFT spectrum point is carried out to its corresponding 10 whole circle vibration data array respectively using reference frequency is calculated
Analysis.Since data are in the data for intercepting into whole circle before, spectrum analysis here is exactly order analysis.To obtain 10 circles
Order amplitude array [Am01, Am02 ... Am0 ((La1-La0)/2)], [Am11, Am12 ... Am1 ((La2-La1)/
2)] ... [Am91, Am92 ... Am9 ((La11-La10)/2)] and order phase data [Ph01, Ph02 ... Ph0 ((La1-
La0)/2)]、[Ph11,Ph12,…Ph1((La2-La1)/2)]、…[Ph91,Ph92,…Ph9((La11-La10)/2)]。
Dynamic balancing adjustment phase identification
Extract 10 amplitudes composition of the 1st rank respectively from obtained 10 circle order amplitude arrays and order phase data
Amplitude array
[0.2502,0.2399,0.2450,0.2531,0.2489,0.2479,0.2392,0.2522,0.2512,
0.2488] and 10 phases form phase array
[1.888,2.8449,3.035,2.5879,2.5755,3,4231,3,9742,4.0746,2.9350,
3.4603]。
To amplitude array by the sequence amplitude array after being ranked up to obtain one group of sequence from small to large
[0.2392,0.2399,0.2450,002479,0.2488,0.2489,0.2502,0.2512,0.2522,
0.2531] and the corresponding cenotype bit array of the sorting position
[3.9742,2.8449,3.0350,3.4231,3.4603,2.5755,1.8888,2.9350,4.0746,
2.5876], then intermediate 6 values are taken averagely to be obtained average amplitude from sequence amplitude array
At=(0.2450+002479+0.2488+0.2489+0.2502+0.2512)/6=0.2487, while taking cenotype
6 values are averagely obtained average phase among bit array
Pt=(3.0350+3.4231+3.4603+2.5755+1.8888+2.9350)/6=2.8863.Due to average width
Value At is greater than the vibration values 0.1 of setting, then needs to carry out dynamic balancing adjustment, and the phase position for needing to increase weight stack is Pt+
180=182.8863.
It is last it is to be noted that:The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations.To the greatest extent
Present invention has been described in detail with reference to the aforementioned embodiments for pipe, those skilled in the art should understand that:It is still
It is possible to modify the technical solutions described in the foregoing embodiments, or part of technical characteristic is equally replaced
It changes;And these are modified or replaceed, the essence for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution
Mind and range.
Claims (3)
1. a kind of rotor model.test system dynamic balancing adjustment phase recognition methods, the rotor model suitable for order analysis complete cycle is tried
It tests, which is characterized in that include the following steps:
Step 1: according to rotor blade the piece number n, test speed Speed, distinguishable minimum phase angle P0 and percentage error
Revolving speed Rt calculates sample frequency Fn as follows:
Fn=max (int (Speed × n × 6/ (P0 × 1000))+1) × 1000, (int (1250/ (Rt × 30000))+1) ×
1000}
Wherein, int is downward round numbers, and max is to take maximum number from two numbers in bracket;
Step 2: under the liter rotary regimes of rotor platform tach signal and vibration signal synchronize acquisition;
The tach signal of rotor platform and vibration signal are accessed simultaneously in the data collection system with synchronous acquisition function, it is described
Sample frequency Fn is provided in data collection system;
By test speed Speed as defined in test bench running to test mission book, then carry out always away from be incremented by be manipulated to it is pending
Dynamic balancing adjustment always away from;
It is to be tested it is in stable condition after, based on data collection system setting sample frequency Fn obtain rotary speed data array [Spa1,
Spa2 ..., Spax] and vibration data number [Vi1, Vi2 ..., Vix];
Wherein, X is number of sampling points;
Step 3: carrying out location finding to revolving speed peak value present in rotary speed data array;
A), the numerical value in the rotary speed data array collected is taken absolute value;
B), a revolving speed peak threshold Psp is set, using threshold value Psp to rotary speed data array [Spa1, Spa2 ..., Spax]
Carry out peak value searching;
C), after having searched for rotary speed data, one group of peak position array [La1, La2 ..., Lay] can be obtained;
Wherein y is the circle number of peak position number or revolving speed and vibration signal;
Step 4: carrying out the interception of vibration signal complete cycle using the peak position that step 3 obtains;
D), using peak position array [La1, La2 ..., Lay] be calculated y-1 peak value starting point array [La1, La2 ...,
La (y-1)] and corresponding data length number [La2-La1, La3-La2 ..., Lay-La (y-1)];
E), calculated peak value starting point array and data length number group, by vibration data array [Vi1, Vi2 ...,
Vix] interception at the whole circle vibration data array of y-1 [Vi11, Vi12 ... Vi0 (La2-La1)], [Vi21, Vi22 ... Vi2
(La2-La1)] ... [2 ... Vi (y-2) of Vi (y-1) 1, Vi (y-2) (Lay-La (y-1))], and calculate separately y-1 whole circle vibrations
The dynamic corresponding reference frequency Fb of array of data:
Fb (y-1)=Fn/ (Lay-La (y-1))
Step 5: to vibration signal order analysis complete cycle of interception;
Using be calculated reference frequency Fb0, Fb1 ... Fb (y-1) is respectively to its corresponding y-1 whole circle vibration data arrays
FFT spectrum analysis is carried out, to obtain the order amplitude number and order phase data of y-1 circle;
Step 6: carrying out identification judgement to rotor model testing stand dynamic balancing adjustment phase;
Choose the y-1 amplitude composition amplitude array and y-1 phase composition phase array of the first rank;
It is ranked up to obtain the sequence amplitude array [A after one group of sequence by sizen0,An1,…An9] and sorting position correspondence
Sequence phase array [Pn0,Pn1,…Pn(n-1)];
Intermediate n-4 value is taken averagely to be obtained average amplitude At=(A from sequence amplitude arrayn2+An3+An4+An5+…+
A(n-3)+A(n-2))/(n-4),
Take phase array [Pn0,Pn1,…Pn(n-1)] intermediate n-4 value averagely obtained average phase Pt=(Pn2+Pn3+Pn4+
Pn5+…+P(n-3)+P(n-2))/(n-4);
If average amplitude At is greater than the vibration values of setting, need to carry out dynamic balancing adjustment, and needs to increase the phase of weight stack
Position is Pt+180.
2. rotor model.test system dynamic balancing adjustment phase recognition methods according to claim 1, it is characterised in that:
In the step 3,
Rotary speed data array [Spa1, Spa2 ..., Spax] is taken out into odd number t (t >=3) a point from current search position first,
If this t point (t-1)/2 all bigger than threshold value Psp and preceding point is gradually increased, while rear (t-1)/2 point gradually subtracts
Small, then the position of (t-1)/2+1 point in this section of rotary speed data had both been a peak point;
Next records position of the peak point in rotary speed data array [Spa1, Spa2 ..., Spax], while will currently search
S (s >=t) a data point is skipped in rope position backward, repeats step 1;Otherwise searching position is moved into 1 position from current location backward
It sets, repeats the above steps.
3. rotor model.test system dynamic balancing adjustment phase recognition methods according to claim 1, it is characterised in that:The step
In two,
It is to be tested it is in stable condition after, based on data collection system setting sample frequency Fn obtain at least 10 circle rotary speed data numbers
Group [Spa1, Spa2 ..., Spax] and vibration data number [Vi1, Vi2 ..., Vix].
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810668288.8A CN108910079B (en) | 2018-06-26 | 2018-06-26 | Rotor wing test bed dynamic balance adjustment phase identification method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810668288.8A CN108910079B (en) | 2018-06-26 | 2018-06-26 | Rotor wing test bed dynamic balance adjustment phase identification method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108910079A true CN108910079A (en) | 2018-11-30 |
CN108910079B CN108910079B (en) | 2020-07-17 |
Family
ID=64422467
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810668288.8A Active CN108910079B (en) | 2018-06-26 | 2018-06-26 | Rotor wing test bed dynamic balance adjustment phase identification method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108910079B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110794367A (en) * | 2019-10-12 | 2020-02-14 | 中国直升机设计研究所 | Method for positioning sound source interfered by propeller vortexes |
CN110844115A (en) * | 2019-10-18 | 2020-02-28 | 中国直升机设计研究所 | Method for judging effectiveness of data of propeller vortex interference noise and blade flapping load |
CN110844110A (en) * | 2019-10-11 | 2020-02-28 | 中国直升机设计研究所 | Method for determining phase of blade motion parameter |
CN112326164A (en) * | 2020-09-25 | 2021-02-05 | 中国航空工业集团公司上海航空测控技术研究所 | Coaxial forward and reverse rotation dual-rotor vibration source separation method |
CN112468050A (en) * | 2020-11-03 | 2021-03-09 | 中国直升机设计研究所 | Rotating speed control method capable of controlling motor phase |
CN112487360A (en) * | 2020-11-20 | 2021-03-12 | 中国直升机设计研究所 | Rotor blade vortex interference noise whole-period averaging method based on statistical analysis |
CN112504426A (en) * | 2020-11-20 | 2021-03-16 | 中国直升机设计研究所 | Peak search-based rotor blade vortex interference noise whole-period averaging method |
CN112556930A (en) * | 2020-12-25 | 2021-03-26 | 中国人民解放军32382部队 | Helicopter moving part vibration signal data quality calculation method |
CN113720614A (en) * | 2021-07-19 | 2021-11-30 | 北京化工大学 | Fundamental frequency vibration signal reduction method for propeller |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3802273A (en) * | 1972-10-17 | 1974-04-09 | Chadwick Elect Inc H | Helicopter rotor balancing method and system |
CN104316266A (en) * | 2014-08-26 | 2015-01-28 | 中国直升机设计研究所 | Correction-function-contained dynamic balance adjustment phase calculating method of helicopter model |
CN106289645A (en) * | 2016-11-07 | 2017-01-04 | 北京东方振动和噪声技术研究所 | Rotor dynamic balancing assay method based on steady stimulation method and device |
CN106768792A (en) * | 2016-11-29 | 2017-05-31 | 中国直升机设计研究所 | A kind of method for carrying out rotor model testing stand dynamic balancing adjustment |
CN108088546A (en) * | 2017-12-03 | 2018-05-29 | 中国直升机设计研究所 | Amplifier terminals loosen recognition methods under a kind of rotor platform working speed state |
CN108106827A (en) * | 2017-12-03 | 2018-06-01 | 中国直升机设计研究所 | A kind of rotor platform rises amplifier terminals under rotary regimes and loosens recognition methods |
CN108119384A (en) * | 2017-10-31 | 2018-06-05 | 沈阳透平机械股份有限公司 | The rotor dynamic balancing adjusting method and device of centrifugal compressor |
-
2018
- 2018-06-26 CN CN201810668288.8A patent/CN108910079B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3802273A (en) * | 1972-10-17 | 1974-04-09 | Chadwick Elect Inc H | Helicopter rotor balancing method and system |
CN104316266A (en) * | 2014-08-26 | 2015-01-28 | 中国直升机设计研究所 | Correction-function-contained dynamic balance adjustment phase calculating method of helicopter model |
CN106289645A (en) * | 2016-11-07 | 2017-01-04 | 北京东方振动和噪声技术研究所 | Rotor dynamic balancing assay method based on steady stimulation method and device |
CN106768792A (en) * | 2016-11-29 | 2017-05-31 | 中国直升机设计研究所 | A kind of method for carrying out rotor model testing stand dynamic balancing adjustment |
CN108119384A (en) * | 2017-10-31 | 2018-06-05 | 沈阳透平机械股份有限公司 | The rotor dynamic balancing adjusting method and device of centrifugal compressor |
CN108088546A (en) * | 2017-12-03 | 2018-05-29 | 中国直升机设计研究所 | Amplifier terminals loosen recognition methods under a kind of rotor platform working speed state |
CN108106827A (en) * | 2017-12-03 | 2018-06-01 | 中国直升机设计研究所 | A kind of rotor platform rises amplifier terminals under rotary regimes and loosens recognition methods |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110844110A (en) * | 2019-10-11 | 2020-02-28 | 中国直升机设计研究所 | Method for determining phase of blade motion parameter |
CN110794367A (en) * | 2019-10-12 | 2020-02-14 | 中国直升机设计研究所 | Method for positioning sound source interfered by propeller vortexes |
CN110794367B (en) * | 2019-10-12 | 2022-06-21 | 中国直升机设计研究所 | Method for positioning sound source interfered by propeller vortexes |
CN110844115B (en) * | 2019-10-18 | 2022-04-12 | 中国直升机设计研究所 | Method for judging effectiveness of data of propeller vortex interference noise and blade flapping load |
CN110844115A (en) * | 2019-10-18 | 2020-02-28 | 中国直升机设计研究所 | Method for judging effectiveness of data of propeller vortex interference noise and blade flapping load |
CN112326164A (en) * | 2020-09-25 | 2021-02-05 | 中国航空工业集团公司上海航空测控技术研究所 | Coaxial forward and reverse rotation dual-rotor vibration source separation method |
CN112468050A (en) * | 2020-11-03 | 2021-03-09 | 中国直升机设计研究所 | Rotating speed control method capable of controlling motor phase |
CN112468050B (en) * | 2020-11-03 | 2023-09-01 | 中国直升机设计研究所 | Rotating speed control method capable of controlling motor phase |
CN112487360B (en) * | 2020-11-20 | 2022-12-30 | 中国直升机设计研究所 | Rotor wing propeller vortex interference noise whole-period averaging method based on statistical analysis |
CN112504426B (en) * | 2020-11-20 | 2022-10-18 | 中国直升机设计研究所 | Peak search-based rotor blade vortex interference noise whole-period averaging method |
CN112504426A (en) * | 2020-11-20 | 2021-03-16 | 中国直升机设计研究所 | Peak search-based rotor blade vortex interference noise whole-period averaging method |
CN112487360A (en) * | 2020-11-20 | 2021-03-12 | 中国直升机设计研究所 | Rotor blade vortex interference noise whole-period averaging method based on statistical analysis |
CN112556930A (en) * | 2020-12-25 | 2021-03-26 | 中国人民解放军32382部队 | Helicopter moving part vibration signal data quality calculation method |
CN112556930B (en) * | 2020-12-25 | 2023-03-03 | 中国人民解放军32382部队 | Helicopter movable part vibration signal data quality calculation method |
CN113720614A (en) * | 2021-07-19 | 2021-11-30 | 北京化工大学 | Fundamental frequency vibration signal reduction method for propeller |
CN113720614B (en) * | 2021-07-19 | 2024-05-28 | 北京化工大学 | Fundamental frequency vibration signal restoration method for propeller |
Also Published As
Publication number | Publication date |
---|---|
CN108910079B (en) | 2020-07-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108910079A (en) | A kind of rotor model.test system dynamic balancing adjustment phase recognition methods | |
CN203053919U (en) | Testing device for vibrating dynamic characteristic transfer tests of trees | |
US7961550B2 (en) | Systems and methods of processing impulses including bullet pulses and/or muzzle pulses in association with time domain representations | |
CN107478990B (en) | A kind of generator electromechanical transient process method for dynamic estimation | |
CN103558003B (en) | A kind of rotor torsion oscillation excitation and vibration analysis experimental system | |
CN106802630B (en) | The Parallel frequency domain method of milling parameter stability prediction based on thin-wall part | |
CN108092304B (en) | Electric power system transient stability discrimination method based on phase track and relative kinetic energy change rate | |
CN104698325B (en) | A kind of method of discrimination of the low-frequency oscillation of power system negative damping mechanism and forced oscillation | |
CN109228886B (en) | Parameter matching method, control device and equipment for electric vehicle power assembly | |
CN108763040A (en) | Collecting method and data acquisition device | |
CN111126633B (en) | Automatic fault diagnosis method for mining multistage centrifugal pump | |
CN105865793B (en) | A kind of method for improving more rotor vibration monitoring of aero-engine precision | |
CN105629060B (en) | Power grid frequency measurement method and device based on optimal baseband filtering | |
CN104977505B (en) | A kind of grid disturbance source localization method based on comprehensive oscillatior | |
CN103335785A (en) | Unsteady dynamic balance method based on order domain-holographic spectrum principle | |
CN104316768A (en) | Negative sequence impedance parameter estimation method for locating three-phase unbalanced disturbance source | |
CN110837003A (en) | Double-window full-phase DFT (discrete Fourier transform) synchronous phasor measurement method and system based on triangular window | |
CN112906210A (en) | Wind turbine generator test bed time delay identification method and system based on instability feature extraction | |
CN108458781A (en) | A kind of remote oscillation analysis system and its analysis method | |
CN106053000A (en) | Missile and carrier rocket modal test method | |
CN108106827A (en) | A kind of rotor platform rises amplifier terminals under rotary regimes and loosens recognition methods | |
CN106991076A (en) | Multimodal extremum search method based on gradient method extremum search | |
CN105510621B (en) | A kind of centrifugal compressor critical speed acquisition methods and device | |
CN107894279B (en) | Synchronous detection method and device for vibration data and working condition data of wind generating set | |
CN106452621B (en) | A kind of peak power, the determination method and device of peak to average power ratio |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |