CN105844250B - A method of maximum pressure rate of rise is recognized based on vibration acceleration signal - Google Patents
A method of maximum pressure rate of rise is recognized based on vibration acceleration signal Download PDFInfo
- Publication number
- CN105844250B CN105844250B CN201610195062.1A CN201610195062A CN105844250B CN 105844250 B CN105844250 B CN 105844250B CN 201610195062 A CN201610195062 A CN 201610195062A CN 105844250 B CN105844250 B CN 105844250B
- Authority
- CN
- China
- Prior art keywords
- rise
- maximum pressure
- acceleration signal
- vibration acceleration
- pressure rate
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Abstract
The invention discloses a kind of methods based on vibration acceleration signal identification maximum pressure rate of rise, comprising the following steps: step 1 filters out vibration acceleration signal medium-high frequency using chebyshev low-pass filter and interferes;Step 2 utilizes the low-frequency disturbance in empirical mode decomposition algorithm removal vibration acceleration signal;Step 3 characterizes maximum pressure rate of rise using the area in current moment is gone out from timing of combustion to maximum pressure rate of rise on processing curve.The High-frequency Interference that actual measureed value of acceleration signal is removed by low-pass filter removes low-frequency disturbance therein using EMD method, effectively raises the signal-to-noise ratio of measured signal.The present invention characterizes maximum pressure rate of rise using the area gone out current moment on processing curve from timing of combustion to maximum pressure rate of rise.The invention proposes the linear relationships between the maximum pressure rate of rise of foundation and IMF area, are modified based on the linear relationship to its identification result.
Description
Technical field
The present invention relates to signal processing technology fields, more particularly to a kind of vibration acceleration signal that is based on to recognize maximum pressure
The method of rate of rise.
Background technique
Maximum pressure rate of rise is as the important feature parameter in Combustion Process Control, generally by adding to the vibration of actual measurement
Speed signal is integrated to obtain.Common time-domain integration method has trapezoid formula integration method and Simpson integration method, trapezoidal
Formula integration method algorithm is simple, and the required calculating time is shorter, but precision is poor, and convergence rate is slow;Simpson integration method phase
For trapezoid formula, accurate height is calculated, however there are still interference in pretreated vibration acceleration signal, especially fire
It burns near initial point, leading to the initial value of integral, there are deviations, influence subsequent integrated value, and the calculating effect of Simpson integration method
Rate is lower.
Summary of the invention
The purpose of the present invention is to solve the above-mentioned problems, provides a kind of based on the maximum pressure of vibration acceleration signal identification
The method of power rate of rise recognizes maximum pressure rate of rise using cylinder lid surface vibration acceleration speed signal, passes through low-pass filtering
Device removes the High-frequency Interference of actual measureed value of acceleration signal, removes low-frequency disturbance therein using EMD method, effectively raises actual measurement
The signal-to-noise ratio of signal.
To achieve the goals above, the present invention adopts the following technical scheme:
A method of maximum pressure rate of rise is recognized based on vibration acceleration signal, comprising the following steps:
Step 1 filters out vibration acceleration signal medium-high frequency using chebyshev low-pass filter and interferes;
Step 2 utilizes the low-frequency disturbance in empirical mode decomposition (EMD) algorithm removal vibration acceleration signal;
Step 3 goes out the characterization of the area in current moment most from timing of combustion to maximum pressure rate of rise using on processing curve
The big rate of pressure rise.
Chebyshev low-pass filter ripple factor is 0.001 in the step 1, order 10.
It is filtered out in the step 1 more than 150 times of fundamental frequencies unrelated with burning exciter response signal in vibration acceleration signal
High-frequency Interference.
Processing curve in the step 3 is to calculate each rank intrinsic mode function and cylinder pressure signal second derivative
Correlation, by composition processing curve after the summation of preceding two ranks intrinsic mode function.
Beneficial effects of the present invention:
The High-frequency Interference that actual measureed value of acceleration signal is removed by low-pass filter removes low frequency therein using EMD method
Interference, effectively raises the signal-to-noise ratio of measured signal.
The present invention is characterized most using the area gone out current moment on processing curve from timing of combustion to maximum pressure rate of rise
The big rate of pressure rise.
The invention proposes the linear relationships between the maximum pressure rate of rise of foundation and IMF area, are based on the linear relationship
Its identification result is modified.
The present invention evaluates in-cylinder combustion state to the research using vibration signal to in-cylinder combustion process feature parameter, right
Combustion process of the internal-combustion engine on-line monitoring, closed-loop control and fault diagnosis have important theory significance and practical value.
Detailed description of the invention
Fig. 1 is that cylinder pressure, the rate of pressure rise and cylinder press second derivative curve;
Fig. 2 is that processing curve and cylinder press second derivative curve;
Fig. 3 (a) is 800r/min recognition result;Fig. 3 (b) is 1000r/min recognition result;
Fig. 3 (c) is 1200r/min recognition result;Fig. 3 (d) is 1400r/min recognition result.
Specific embodiment
The invention will be further described with embodiment with reference to the accompanying drawing.
A method of maximum pressure rate of rise is recognized based on vibration acceleration signal, utilizes chebyshev low-pass filter
High-frequency Interference is filtered out, using the low-frequency disturbance in EMD empirical mode decomposition removal vibration acceleration signal, using on processing curve
Go out the characterization maximum pressure rate of rise of the area in current moment from timing of combustion to maximum pressure rate of rise.
Above-mentioned chebyshev low-pass filter ripple factor is 0.001, and order 10 filters out in vibration acceleration signal
The High-frequency Interferences unrelated with burning exciter response signal more than 150 times of fundamental frequencies.
Above-mentioned EMD is empirical mode decomposition (Empirical Mode Decomposition, EMD), and EMD is according to non-thread
Property, the characteristic time scale of non-stationary signal itself by signal decomposition at one group of stable state and linear data sequence collection, i.e., it is intrinsic
Mode function (IMF).IMF reflects that the internal feature of signal, residual components indicate the trend of signal, need to meet following two item
Part: zero point number is equal with pole number in entire signal or at most differs 1;For any point on signal, by local maximum
The mean value for the envelope that determining envelope and local minimum determine is 0.It has the advantage of identification local feature.
Above-mentioned processing curve is the correlation for calculating each rank IMF and cylinder pressure signal second derivative, preceding two ranks IMF and flat
Equal cylinder pressure signal second derivative correlation is higher, by composition processing curve after preceding two ranks IMF summation.
Fig. 1 shows cylinder pressure, rate of pressure rise curve and cylinder pressure signal second derivative curve, often will pressure in internal combustion engine
Break of the power rate of rise in compression stroke is as timing of combustion, and in cylinder pressure second derivative curve, timing of combustion corresponds to peak value
First preceding zero crossing, i.e. A point, B corresponding is that maximum pressure rate of rise goes out current moment.Fig. 2 shows the cylinder of a circulation
The vibration acceleration curve of interior pressure signal second derivative and IMF1 and the sum of the IMF2 obtained after EMD is decomposed (are defined as
Handle curve), the circle of the first two is A, B two o'clock in Fig. 1.It is corresponding, handle the first two triangle in curve
Shape is respectively used to characterization, and timing of combustion goes out current moment, maximum pressure rate of rise goes out current moment.
Therefore the present invention utilizes and goes out the face in current moment from timing of combustion to maximum pressure rate of rise on vibration acceleration curve
Product, i.e. shaded area shown in Fig. 2 (being defined as IMF area) characterize maximum pressure rate of rise.Fig. 3 (a)-Fig. 3 (d) gives
Maximum pressure rate of rise and IMF area under different operating conditions give simultaneously in figure for the ease of analysis and utilize linear regression method
Linear relationship between the maximum pressure rate of rise established offline and IMF area.As seen from the figure, under different operating conditions, degree of fitting difference
It is 0.97,0.97,0.91,0.91, there is good linear relationship therebetween.Therefore in actual application, after processing being based on
Vibration acceleration signal calculate corresponding IMF area, in conjunction with off-line calibration different operating conditions under relationship between the two, into
And estimate the maximum pressure rate of rise under actual condition.Wherein, 800r/min, 1000r/min, 1200r/min, 10Nm and
The maximum deviation of 1400r/min, 10Nm operating condition is respectively 2.2%, 0.72%, 3.41% and 3.97%, shows that IMF area can
For characterizing maximum pressure rate of rise.
Above-mentioned, although the foregoing specific embodiments of the present invention is described with reference to the accompanying drawings, not protects model to the present invention
The limitation enclosed, those skilled in the art should understand that, based on the technical solutions of the present invention, those skilled in the art are not
Need to make the creative labor the various modifications or changes that can be made still within protection scope of the present invention.
Claims (1)
1. a kind of method based on vibration acceleration signal identification maximum pressure rate of rise, characterized in that the following steps are included:
Step 1 filters out vibration acceleration signal medium-high frequency using chebyshev low-pass filter and interferes;
Step 2 utilizes the low-frequency disturbance in empirical mode decomposition algorithm removal vibration acceleration signal;
Step 3 characterizes maximum pressure using the area in current moment is gone out from timing of combustion to maximum pressure rate of rise on processing curve
Power rate of rise;
Chebyshev low-pass filter ripple factor is 0.001 in the step 1, and order 10 filters out vibration acceleration signal
In the more than 150 times of fundamental frequencies High-frequency Interference unrelated with burning exciter response signal;
Processing curve in the step 3 is that each rank intrinsic mode function of calculating is related to cylinder pressure signal second derivative
Property, by composition processing curve after the summation of preceding two ranks intrinsic mode function.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610195062.1A CN105844250B (en) | 2016-03-31 | 2016-03-31 | A method of maximum pressure rate of rise is recognized based on vibration acceleration signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610195062.1A CN105844250B (en) | 2016-03-31 | 2016-03-31 | A method of maximum pressure rate of rise is recognized based on vibration acceleration signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105844250A CN105844250A (en) | 2016-08-10 |
CN105844250B true CN105844250B (en) | 2019-12-03 |
Family
ID=56596223
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610195062.1A Active CN105844250B (en) | 2016-03-31 | 2016-03-31 | A method of maximum pressure rate of rise is recognized based on vibration acceleration signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105844250B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108644023A (en) * | 2018-04-25 | 2018-10-12 | 江苏大学 | Internal combustion engine cylinder pressure signal high frequency harmonic components characterizing method |
CN110567727B (en) * | 2019-08-13 | 2021-10-01 | 北京化工大学 | Diesel engine output power evaluation method based on vibration acceleration signal |
CN113532711B (en) * | 2021-06-03 | 2024-04-05 | 西安理工大学 | Printing pressure identification method for central embossing cylinder of satellite type flexographic printing machine |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008037260A3 (en) * | 2006-09-26 | 2008-05-15 | Morpheus Medical | Methods for a movement and vibration analyzer (mva) |
CN102494626A (en) * | 2011-11-18 | 2012-06-13 | 中国船舶重工集团公司第七0四研究所 | Method for eliminating vibration noise in shaft torsional deformation test |
CN104866671A (en) * | 2015-05-26 | 2015-08-26 | 山东大学 | Vibration speed signal-based method for recognizing maximal pressure rise rate |
-
2016
- 2016-03-31 CN CN201610195062.1A patent/CN105844250B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008037260A3 (en) * | 2006-09-26 | 2008-05-15 | Morpheus Medical | Methods for a movement and vibration analyzer (mva) |
CN102494626A (en) * | 2011-11-18 | 2012-06-13 | 中国船舶重工集团公司第七0四研究所 | Method for eliminating vibration noise in shaft torsional deformation test |
CN104866671A (en) * | 2015-05-26 | 2015-08-26 | 山东大学 | Vibration speed signal-based method for recognizing maximal pressure rise rate |
Non-Patent Citations (2)
Title |
---|
"EMD算法研究及其在信号去噪中的应用;王婷;《中国博士学位论文全文数据库 信息科技辑》;20110715;第I136-26页 * |
基于振动信号评价柴油机缸内燃烧状态的研究;唐娟;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20110715;第C039-22页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105844250A (en) | 2016-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105844250B (en) | A method of maximum pressure rate of rise is recognized based on vibration acceleration signal | |
CN105893673B (en) | A method of combustion characteristic parameter is recognized based on vibration acceleration characteristic point | |
CN105067966B (en) | The low-voltage alternating-current fault arc detection method of feature based modal components energy spectrometer | |
CN103216461B (en) | The stall inception identification method of axial-flow compressor | |
CN102401866B (en) | Detecting method of detuning fault element of triple-tuned DC filter | |
CN106203355A (en) | The detection method of a kind of low-frequency oscillation of electric power system and device | |
CN110726875B (en) | New energy flexible direct-current grid-connected transient harmonic detection method and system | |
Pillay et al. | A data driven approach to performance assessment of PID controllers for setpoint tracking | |
CN105003373A (en) | Fault diagnosing method of observing injection timing based on characteristic parameters of waveform of rail pressure | |
Broatch et al. | A direct transform for determining the trapped mass on an internal combustion engine based on the in-cylinder pressure resonance phenomenon | |
CN102928225A (en) | Self-adaptive resonance demodulation method for fault diagnosis of rolling bearing | |
JPWO2012035647A1 (en) | Apparatus for estimating heat release rate of internal combustion engine | |
CN106226605A (en) | A kind of PWM self-sensing method | |
CN105181120A (en) | High-sensitivity transformer winding loosening determination method | |
CN103995180A (en) | Power system frequency estimation method taking inequality constraints into consideration | |
CN110567727B (en) | Diesel engine output power evaluation method based on vibration acceleration signal | |
Dambros et al. | Signal preprocessing for stiction detection methods | |
Ou et al. | Identification and reconstruction of anomalous sensing data for combustion analysis of marine diesel engines | |
CN104849548B (en) | A kind of electric system instantaneous frequency monitoring method and system | |
CN108875581A (en) | A kind of description method at delayed phase angle between Internal Combustion Engine vibration velocity signal and cylinder pressure signal | |
CN110542855B (en) | Load switch event detection method and system based on discrete cosine transform | |
CN103631991A (en) | Parameter identification system and method of prime mover speed regulating system | |
CN103632017A (en) | Pattern recognition-based method for improving signal-to-noise ratio of vibration signals on surface of internal combustion engine | |
CN113820004A (en) | Robust vibration signal initial phase estimation method | |
Zhang et al. | Knock threshold detection in turbocharged gasoline engine using EEMD and bispectrum |
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 |