WO2022007382A1 - Method for solving degree of contribution of non-combustion excitation signal on basis of partial coherence function - Google Patents

Method for solving degree of contribution of non-combustion excitation signal on basis of partial coherence function Download PDF

Info

Publication number
WO2022007382A1
WO2022007382A1 PCT/CN2021/071201 CN2021071201W WO2022007382A1 WO 2022007382 A1 WO2022007382 A1 WO 2022007382A1 CN 2021071201 W CN2021071201 W CN 2021071201W WO 2022007382 A1 WO2022007382 A1 WO 2022007382A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
excitation
vibration
combustion
contribution
Prior art date
Application number
PCT/CN2021/071201
Other languages
French (fr)
Chinese (zh)
Inventor
赵秀亮
王丽梅
潘邦雄
汪若尘
杨政宇
赵明明
Original Assignee
江苏大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 江苏大学 filed Critical 江苏大学
Publication of WO2022007382A1 publication Critical patent/WO2022007382A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • G01M15/04Testing internal-combustion engines
    • G01M15/12Testing internal-combustion engines by monitoring vibrations

Definitions

  • the invention relates to a method for calculating a contribution degree, in particular to a method for calculating the contribution degree of a non-combustion excitation signal based on a partial coherence function, and belongs to the technical field of internal combustion engines.
  • the vibration signal of the internal combustion engine is the result of the combined action of combustion excitation and non-combustion excitation such as reciprocating inertial force, valve seating impact force, piston reversing impact force, piston side pressure and crankshaft main journal load.
  • combustion excitation and non-combustion excitation such as reciprocating inertial force, valve seating impact force, piston reversing impact force, piston side pressure and crankshaft main journal load.
  • non-combustion excitation is coupled with the combustion excitation in the time and frequency domains, which makes it difficult to eliminate the non-combustion excitation response signal, and even affects the integrity of the combustion excitation response signal.
  • Contribution degree refers to the ratio of the spectral amplitude of the vibration velocity signal when the main excitation is applied alone to the amplitude of the vibration velocity signal when all excitations are applied at a certain frequency.
  • the characteristics of the excitation and its response signals are comprehensively analyzed, the coupling relationship between the excitation and response signals is studied, and the contribution of the main excitation response signals in the vibration signal is verified, which will eliminate the non-combustion excitation response signals and realize the acquisition of combustion information. and combustion process control to provide a theoretical basis.
  • the present invention provides a method for solving the contribution degree of the non-combustion excitation signal based on the partial coherence function. Analyze the characteristics and intensity of non-combustion excitation sources to determine the main non-combustion excitation sources; determine the non-combustion excitation that plays a primary role in the vibration signal by analyzing the characteristics of the vibration signal; use the partial coherence function to calculate the contribution and variation of the combustion excitation .
  • Step 1 Determine the main non-combustion excitation source, analyze the time domain characteristics of valve opening and closing excitation and adjacent cylinder combustion excitation, solve the excitation source strength of reciprocating inertial force and piston commutation impact force, and obtain excitation sources sensitive to the influence of vibration signals ;
  • Step 2 Determine the non-combustion excitation that plays a primary role in the vibration signal, simplify the actual internal combustion engine structure, retain the body and cylinder head, and establish an engine simulation model; select the measuring point on the model to apply a force hammer, and compare the simulation vibration signal and Measured vibration signal to verify the validity of the model;
  • each excitation vibration response signal is analyzed by compressive wavelet transform technology, and the non-combustion excitation which plays a primary role in the vibration signal is obtained;
  • Step 3 According to the primary non-combustion excitation determined in Step 2, use the partial coherence function to calculate the contribution of the excitation signal to the vibration signal at a specific frequency.
  • the method for solving the excitation source strength of the reciprocating inertial force and the reversing impact force of the piston in the step 1 is as follows: according to the measured pressure signal in the cylinder, the excitation curve of the reciprocating inertial force and the excitation curve of the reversing impact force of the piston, respectively obtain their maximum value.
  • a 10-node tetrahedral mesh is used when establishing the engine simulation model.
  • the body is fixed on the test bench through 4 brackets, and the brackets are restricted to be rigid bodies without deformation. degrees of freedom of orientation.
  • the method for verifying the validity of the simulated vibration signal and the measured vibration signal model in the second step is:
  • the method of analyzing the time-frequency characteristics of each excitation vibration response signal through the compressive wavelet transform technology is as follows:
  • the time-frequency analysis of the excitation response signal is performed to obtain the time-frequency distribution of the vibration displacement, vibration velocity and vibration acceleration signals, and the vibration signal characteristics caused by the reciprocating inertial force excitation and the piston commutation excitation alone are obtained;
  • Two inputs x 1 (t), x 2 (t), the partial coherence function of a single output signal y(t) can be defined as:
  • the present invention determines the contribution of non-combustion excitation signals to the vibration signal, which is the basis for research to eliminate interference in the vibration signal, and helps to repair part of the combustion information in the vibration signal, which is helpful for realizing the extraction of combustion information by using the vibration signal, and the internal combustion engine.
  • Working condition monitoring, closed-loop control of combustion process and fault diagnosis have important theoretical significance and practical value.
  • Fig. 1 is the working phase diagram of valve opening and closing of an internal combustion engine of the present invention
  • Fig. 2 is the measured vibration acceleration signal diagram of the present invention
  • Fig. 3 is the pressure signal in the cylinder of the present invention and the reciprocating inertial force excitation curve obtained by calculation;
  • Fig. 4 is the motive simulation model diagram of the present invention.
  • Fig. 5 is the graph of the vibration acceleration signal obtained by the actual measurement vibration acceleration signal of the present invention and simulation
  • Fig. 6 is the time frequency spectrum diagram of the reciprocating inertial force signal of the present invention.
  • Fig. 7 is the time-frequency distribution diagram of the reciprocating inertial force excitation response signal of the present invention; wherein, Fig. 7 (a), (b) are the time-frequency distribution diagrams of vibration displacement and vibration velocity signals respectively, Fig. 7 (c) is The time-frequency distribution diagram of the vibration acceleration signal, Fig. 7(d) is a partial enlarged view of Fig. 5(c);
  • Fig. 8 is a statistical diagram of the contribution degree of the second-order reciprocating inertial force of the present invention
  • Fig. 8(a) is a statistical diagram of the contribution degree when the torque is 40N m and the rotational speeds are respectively 800rpm, 1200rpm, 1650rpm and 2200rpm
  • Fig. 8( b) shows the contribution statistics when the rotational speed is 1650rpm and the torque is 20N ⁇ m, 40N ⁇ m, 60N ⁇ m, 80N ⁇ m, and 100N ⁇ m, respectively;
  • valve opening and closing excitations are also important excitation sources of vibration response signals.
  • FIG. 2 it shows the measured vibration acceleration signal at 2200 r/min and 60 N m.
  • the working process of the two-cylinder engine is that the two cylinders are alternately performed.
  • the firing time of cylinder 1 lags behind the firing time of cylinder 2 by 180 °C A
  • the ignition time of cylinder 2 lags behind cylinder 1 by 540°CA. It takes a certain amount of time for the vibration signal to recover from excitation. The higher the rotational speed, the shorter the interval between the occurrences of each excitation.
  • the measurement is performed at the rated speed.
  • the vibration acceleration signal of the cylinder head surface is taken and analyzed.
  • A1 is the surface vibration acceleration of the 1-cylinder cylinder head
  • A2 is the surface vibration acceleration of the 2-cylinder cylinder head.
  • the vibration signal waveforms of the cylinder head surface measured by the two cylinders should be the same, with only a fixed deviation in phase.
  • A1 and A2 have similar variation laws, and the vibration acceleration signal energy is the strongest near their respective combustion top dead centers.
  • the vibration acceleration signal has recovered to be stable, indicating that the combustion excitation of the adjacent cylinder does not affect the combustion process of the current cylinder.
  • FIG. 5 it shows the hammer excitation signal, the measured vibration acceleration signal, the measured vibration acceleration signal after 20 points of smoothing, and the vibration acceleration signal obtained by simulation.
  • the force hammer excitation was applied at point 2, and the hammer signal and the cylinder head vibration acceleration signal at measuring point 1 were measured; at the same time, the measured force hammer signal was applied to the model measuring point 2, and the vibration at measuring point 1 was simulated. acceleration signal.
  • the results in the figure show that the simulated vibration acceleration signal and the measured vibration acceleration signal have a similar trend in waveform, and the amplitude and phase are basically the same, which indicates that the established engine model can be used to analyze the characteristics and coupling relationship of each excitation response signal.
  • the invention adopts the compression wavelet transform technology to analyze the time-frequency characteristics of each excitation vibration response signal.
  • Figure 7 it is the time-frequency analysis result of the reciprocating inertial force excitation response signal.
  • Figures 7(a) and (b) are the time-frequency distributions of the vibration displacement and vibration velocity signals, respectively.
  • the time-frequency distributions of the vibration displacement signals and vibration velocity signals in the frequency domain are similar to the time-frequency distribution of the reciprocating inertial force signal.
  • the frequencies of the reciprocating inertial force and the second-order reciprocating inertial force radiate from the center to the surroundings;
  • Fig. 7(c) is the time-frequency distribution of the vibration acceleration signal, and
  • Fig. 7(d) is a partial enlarged view of Fig. 7(c).
  • the vibration acceleration signal also has a strong energy distribution around the frequencies of the first-order reciprocating inertial force and the second-order reciprocating inertial force. In addition, there are also some energy around 1000Hz and 2000Hz. It can be seen that the energy of the reciprocating inertial force excitation response signal is relatively concentrated in the frequency band, and runs through the entire working cycle of the internal combustion engine. In other words, both the reciprocating inertial force excitation response signal and the combustion excitation vibration response signal are coupled in the time-frequency domain.
  • the simulation After analyzing the characteristics of the vibration signal caused by the reciprocating inertial force excitation and the piston reversing excitation, it is believed that the main influence on the vibration response signal is the reciprocating inertial force excitation.
  • the simulation obtains the vibration velocity signal when the three excitations of piston reversing impact, reciprocating inertia force and cylinder pressure signal act together, and the simulation obtains the vibration velocity when the other two excitations work together without considering the piston reversing impact Signal, the vibration velocity signal when the other two excitations work together without considering the reciprocating inertial force.
  • the present invention takes the vibration velocity signal as an example, and further studies the correlation variation law between the reciprocating inertial force and the vibration signal at the key frequency.
  • the invention utilizes the partial coherence function to approximately determine the contribution of the reciprocating inertial force excitation response signal to the vibration signal.
  • Two inputs x 1 (t), x 2 (t), the partial coherence function of a single output signal y(t) can be defined as:
  • Figure 8(a) shows the contribution when the torque is 40N ⁇ m and the rotational speed is 800rpm, 1200rpm, 1650rpm, and 2200rpm, respectively.
  • Figure 8(b) shows that the rotational speed is 1650rpm, and the torque is 20N ⁇ m and 40N ⁇ m, respectively. , 60N ⁇ m, 80N ⁇ m, 100N ⁇ m contribution.
  • the contribution of the reciprocating inertial force is defined as: at a certain frequency, the ratio of the spectral amplitude of the vibration velocity signal when the reciprocating inertial force is applied alone to the amplitude of the vibration velocity signal when all excitations are applied. It can be seen from the figure that as the rotational speed increases, the contribution of the reciprocating inertial force to the vibration velocity signal increases significantly, and the contribution decreases slightly as the torque increases. This is because the relationship between the reciprocating inertial force and the rotational speed is a quadratic function.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)
  • Testing Of Engines (AREA)

Abstract

A method for solving the degree of contribution of a non-combustion excitation signal on the basis of a partial coherence function. The method comprises: step one, determining a main non-combustion excitation source; step two, determining non-combustion excitation having a primary effect on a vibration signal; and step three, according to the primary non-combustion excitation determined in step two, calculating the degree of contribution of an excitation signal to the vibration signal at a particular frequency by using a partial coherence function. The determination of the degree of contribution of a non-combustion excitation signal to a vibration signal is the basis for researching the elimination of interference in the vibration signal, is conducive to repairing partial combustion information in the vibration signal, and has an important theoretical significance and practical value in extracting combustion information using the vibration signal, monitoring the operating state of an internal combustion engine, closed-loop control over a combustion process, and fault diagnosis.

Description

基于偏相干函数求解非燃烧激励信号贡献度的方法A method for calculating the contribution of non-combustion excitation signal based on partial coherence function 技术领域technical field
本发明涉及一种贡献度的计算方法,特别涉及一种基于偏相干函数求解非燃烧激励信号贡献度的方法,属于内燃机技术领域。The invention relates to a method for calculating a contribution degree, in particular to a method for calculating the contribution degree of a non-combustion excitation signal based on a partial coherence function, and belongs to the technical field of internal combustion engines.
背景技术Background technique
内燃机振动信号是燃烧激励与往复惯性力、气门落座冲击力、活塞换向撞击力、活塞侧压力及曲轴主轴颈负荷等非燃烧激励共同作用的结果。且非燃烧激励在时域和频域上和燃烧激励存在耦合,导致非燃烧激励响应信号的剔除存在一定得困难,甚至会影响燃烧激励响应信号的完整度。The vibration signal of the internal combustion engine is the result of the combined action of combustion excitation and non-combustion excitation such as reciprocating inertial force, valve seating impact force, piston reversing impact force, piston side pressure and crankshaft main journal load. In addition, the non-combustion excitation is coupled with the combustion excitation in the time and frequency domains, which makes it difficult to eliminate the non-combustion excitation response signal, and even affects the integrity of the combustion excitation response signal.
山东大学唐娟在《内燃机学报》第29卷(2011)第5期中发表的《非燃烧激励对缸盖燃烧激励振动影响模拟分析》一文中模拟了活塞换向撞击、侧压力、曲轴主轴颈负荷单独作用及各激励和燃烧激励共同作用时引起的缸盖表面振动信号,并对各激励振动响应信号的时频特点进行了分析。结果表明,各激励引起的缸盖表面振动位移主要反映了2KHz以下的低频振动,振动速度和加速度能同时反映***的低频及高频振动,在时域和频域上各激励振动响应存在耦合。但活塞侧压力、曲轴主轴颈负荷对燃烧激励振动响应的影响可以忽略,活塞换向撞击对燃烧激励引起的振动位移和振动速度的影响可以忽略;但对燃烧激励引起的振动加速度的影响不能忽略。In the article "Simulation Analysis of Influence of Non-combustion Excitation on Cylinder Head Combustion Excitation Vibration" published by Tang Juan of Shandong University in Journal of Internal Combustion Engines, Vol. 29 (2011), Issue 5 The vibration signals of the cylinder head surface caused by the independent action and the combined action of each excitation and combustion excitation are analyzed, and the time-frequency characteristics of each excitation vibration response signal are analyzed. The results show that the vibration displacement of the cylinder head surface caused by each excitation mainly reflects the low frequency vibration below 2KHz, the vibration velocity and acceleration can reflect the low frequency and high frequency vibration of the system at the same time, and the vibration response of each excitation is coupled in the time and frequency domains. However, the influence of piston side pressure and crankshaft main journal load on the vibration response of combustion excitation can be ignored, and the influence of piston reversing impact on the vibration displacement and vibration velocity caused by combustion excitation can be ignored; but the influence on vibration acceleration caused by combustion excitation cannot be ignored. .
贡献度是指在某一频率处,单独施加主要激励时振动速度信号的频谱幅值与施加所有激励时振动速度信号幅值的比值。Contribution degree refers to the ratio of the spectral amplitude of the vibration velocity signal when the main excitation is applied alone to the amplitude of the vibration velocity signal when all excitations are applied at a certain frequency.
为此,对激励及其响应信号特性进行全面分析,研究各激励响应信号间耦合关系,探明振动信号中主要激励响应信号贡献度,将为非燃烧激励响应信号的剔除,实现燃烧信息的获取和燃烧过程控制提供理论依据。Therefore, the characteristics of the excitation and its response signals are comprehensively analyzed, the coupling relationship between the excitation and response signals is studied, and the contribution of the main excitation response signals in the vibration signal is verified, which will eliminate the non-combustion excitation response signals and realize the acquisition of combustion information. and combustion process control to provide a theoretical basis.
发明内容SUMMARY OF THE INVENTION
发明目的:针对现有技术中对各激励响应信号间耦合关系以及振动信号中主要激励响应信号贡献度研究不足,本发明提供了一种基于偏相干函数求解非燃烧激励信号贡献度的方法,根据非燃烧激励源特性和强度分析,确定主要非燃烧激励源;通过对振动信号特性进行分析,确定对振动信号起首要作用的非燃烧激励;利用偏相干函数计算该燃烧激励的贡献度及变化规律。Purpose of the invention: In view of the insufficient research on the coupling relationship between various excitation response signals and the contribution degree of the main excitation response signal in the vibration signal in the prior art, the present invention provides a method for solving the contribution degree of the non-combustion excitation signal based on the partial coherence function. Analyze the characteristics and intensity of non-combustion excitation sources to determine the main non-combustion excitation sources; determine the non-combustion excitation that plays a primary role in the vibration signal by analyzing the characteristics of the vibration signal; use the partial coherence function to calculate the contribution and variation of the combustion excitation .
技术方案:一种基于偏相干函数求解非燃烧激励信号贡献度的方法,包括以下步骤:Technical solution: a method for calculating the contribution of non-combustion excitation signals based on a partial coherence function, comprising the following steps:
步骤一、确定主要非燃烧激励源,分析气门开启关闭激励和相邻缸燃烧激励的时域特点, 求解往复惯性力和活塞换向撞击力的激励源强度,得到对振动信号影响敏感的激励源;Step 1: Determine the main non-combustion excitation source, analyze the time domain characteristics of valve opening and closing excitation and adjacent cylinder combustion excitation, solve the excitation source strength of reciprocating inertial force and piston commutation impact force, and obtain excitation sources sensitive to the influence of vibration signals ;
步骤二、确定对振动信号起首要作用的非燃烧激励,对实际内燃机结构进行简化,保留机体和缸盖,建立发动机仿真模型;在模型上选定测点施加力锤,通过对比仿真振动信号和实测振动信号验证模型有效性;Step 2: Determine the non-combustion excitation that plays a primary role in the vibration signal, simplify the actual internal combustion engine structure, retain the body and cylinder head, and establish an engine simulation model; select the measuring point on the model to apply a force hammer, and compare the simulation vibration signal and Measured vibration signal to verify the validity of the model;
通过压缩小波变换技术对各激励振动响应信号的时频特点进行了分析,得到对振动信号起首要作用的非燃烧激励;The time-frequency characteristics of each excitation vibration response signal are analyzed by compressive wavelet transform technology, and the non-combustion excitation which plays a primary role in the vibration signal is obtained;
步骤三、根据步骤二中确定的首要非燃烧激励,利用偏相干函数计算该激励信号在特定频率处对振动信号的贡献度。Step 3: According to the primary non-combustion excitation determined in Step 2, use the partial coherence function to calculate the contribution of the excitation signal to the vibration signal at a specific frequency.
所述步骤一中求解往复惯性力和活塞换向撞击力的激励源强度的方法为:根据测得的缸内压力信号、往复惯性力的激励曲线和活塞换向撞击力的激励曲线分别获取其最大值。The method for solving the excitation source strength of the reciprocating inertial force and the reversing impact force of the piston in the step 1 is as follows: according to the measured pressure signal in the cylinder, the excitation curve of the reciprocating inertial force and the excitation curve of the reversing impact force of the piston, respectively obtain their maximum value.
所述步骤二中建立发动机仿真模型时采用10节点四面体网格,机体通过4个支架固定在试验台架上,约束支架为刚性体不发生变形,对支架上各节点进行约束,限制其各个方向的自由度。In the second step, a 10-node tetrahedral mesh is used when establishing the engine simulation model. The body is fixed on the test bench through 4 brackets, and the brackets are restricted to be rigid bodies without deformation. degrees of freedom of orientation.
所述步骤二中仿真振动信号和实测振动信号模型有效性验证方法为:The method for verifying the validity of the simulated vibration signal and the measured vibration signal model in the second step is:
首先、在实际发动机上选择施力点施加力锤激励,测取力锤信号和测点处的振动加速度信号;First, select the force application point on the actual engine to apply the force hammer excitation, and measure the force hammer signal and the vibration acceleration signal at the measurement point;
然后,将测取的力锤信号施加到模型的相应的施力点上,仿真得到测点处的振动加速度信号;Then, apply the measured hammer signal to the corresponding force application point of the model, and simulate the vibration acceleration signal at the measurement point;
最后、得出力锤激励信号、实测振动加速度信号、20点光顺后的实测振动加速度信号和仿真得到的振动加速度信号图,得出验证结论。Finally, the excitation signal of the hammer, the measured vibration acceleration signal, the measured vibration acceleration signal after 20 points of smoothing, and the vibration acceleration signal graph obtained by simulation are obtained, and the verification conclusion is drawn.
所述步骤二中通过压缩小波变换技术对各激励振动响应信号的时频特点进行了分析的方法为:In the second step, the method of analyzing the time-frequency characteristics of each excitation vibration response signal through the compressive wavelet transform technology is as follows:
首先、测得往复惯性力信号和活塞换向撞击力的时频谱图,得出能量主要分布的频率范围;First, measure the time-frequency spectrum of the reciprocating inertial force signal and the reversing impact force of the piston, and obtain the frequency range of the main energy distribution;
然后、对激励响应信号时频分析,得到振动位移、振动速度和振动加速度信号的时频分布,得到单独施加往复惯性力激励和活塞换向激励引起的振动信号特性;Then, the time-frequency analysis of the excitation response signal is performed to obtain the time-frequency distribution of the vibration displacement, vibration velocity and vibration acceleration signals, and the vibration signal characteristics caused by the reciprocating inertial force excitation and the piston commutation excitation alone are obtained;
最后、模拟施加活塞换向撞击、往复惯性力和缸内压力信号三个激励共同作用时的振动速度信号,同时模拟不考虑活塞换向撞击时其他两个激励共同作用时的振动速度信号,不考虑往复惯性力时其他两个激励共同作用时的振动速度信号;得到在所有的非燃烧激励中起首要作用的非燃烧激励。Finally, simulate the vibration speed signal when the three excitations of piston reversing impact, reciprocating inertial force and cylinder pressure signal work together, and simulate the vibration speed signal when the other two excitations work together without considering the piston reversing impact. Considering the reciprocating inertial force, the vibration velocity signal when the other two excitations act together; obtain the non-combustion excitation that plays the primary role among all the non-combustion excitations.
利用偏相干函数计算该激励信号在特定频率处对振动信号的贡献度的方法为:The method of using the partial coherence function to calculate the contribution of the excitation signal to the vibration signal at a specific frequency is:
两个输入x 1(t)、x 2(t),单输出信号y(t)的偏相干函数可以定义为: Two inputs x 1 (t), x 2 (t), the partial coherence function of a single output signal y(t) can be defined as:
Figure PCTCN2021071201-appb-000001
Figure PCTCN2021071201-appb-000001
Figure PCTCN2021071201-appb-000002
称为x 2(t)与y(t)之间的偏相干函数,其中x 1(t)的影响已排除掉,S 2Y1(f)为x 2(t)与y(t) 的条件互功率谱,S 221(f)和S YY1(f)分别为去掉x 1(t)影响后x 2(t)与y(t)的条件自功率谱,且满足不等式:
Figure PCTCN2021071201-appb-000002
It is called the partial coherence function between x 2 (t) and y(t), where the influence of x 1 (t) has been excluded, and S 2Y1 (f) is the conditional relationship between x 2 (t) and y(t) power spectrum, S 221 (f) and S YY1 (f) are removed after x 1 (t) Effect x 2 (t) and y (t) from the power spectrum of conditions, and satisfies the inequality:
Figure PCTCN2021071201-appb-000003
Figure PCTCN2021071201-appb-000003
根据计算结果得出特定频率时贡献度的统计结果,最终得出结论。According to the calculation results, the statistical results of the contribution at a specific frequency are obtained, and finally a conclusion is drawn.
有益效果:本发明确定非燃烧激励信号对振动信号的贡献度,是研究剔除振动信号中干扰的基础,而且有助于修复振动信号中的部分燃烧信息,对实现利用振动信号提取燃烧信息、内燃机工作状态监测、燃烧过程闭环控制和故障诊断具有重要的理论意义和实用价值。Beneficial effects: The present invention determines the contribution of non-combustion excitation signals to the vibration signal, which is the basis for research to eliminate interference in the vibration signal, and helps to repair part of the combustion information in the vibration signal, which is helpful for realizing the extraction of combustion information by using the vibration signal, and the internal combustion engine. Working condition monitoring, closed-loop control of combustion process and fault diagnosis have important theoretical significance and practical value.
附图说明Description of drawings
图1为本发明内燃机气门开启、关闭工作相位图;Fig. 1 is the working phase diagram of valve opening and closing of an internal combustion engine of the present invention;
图2为本发明实测振动加速度信号图;Fig. 2 is the measured vibration acceleration signal diagram of the present invention;
图3为本发明缸内压力信号及计算得到的往复惯性力激励曲线;Fig. 3 is the pressure signal in the cylinder of the present invention and the reciprocating inertial force excitation curve obtained by calculation;
图4为本发明动机仿真模型图;Fig. 4 is the motive simulation model diagram of the present invention;
图5为本发明实测振动加速度信号和仿真得到的振动加速度信号的曲线图;Fig. 5 is the graph of the vibration acceleration signal obtained by the actual measurement vibration acceleration signal of the present invention and simulation;
图6为本发明往复惯性力信号的时频谱图;Fig. 6 is the time frequency spectrum diagram of the reciprocating inertial force signal of the present invention;
图7为本发明为本发明往复惯性力激励响应信号时频分布图;其中,图7(a)、(b)分别为振动位移和振动速度信号的时频分布图,图7(c)为振动加速度信号的时频分布图,图7(d)是图5(c)的局部放大图;Fig. 7 is the time-frequency distribution diagram of the reciprocating inertial force excitation response signal of the present invention; wherein, Fig. 7 (a), (b) are the time-frequency distribution diagrams of vibration displacement and vibration velocity signals respectively, Fig. 7 (c) is The time-frequency distribution diagram of the vibration acceleration signal, Fig. 7(d) is a partial enlarged view of Fig. 5(c);
图8为本发明二阶往复惯性力的贡献度统计图;其中,图8(a)是扭矩为40N·m,转速分别为800rpm、1200rpm、1650rpm、2200rpm时的贡献度统计图,图8(b)显示了转速为1650rpm,扭矩分别为20N·m、40N·m、60N·m、80N·m、100N·m时的贡献度统计图;Fig. 8 is a statistical diagram of the contribution degree of the second-order reciprocating inertial force of the present invention; wherein, Fig. 8(a) is a statistical diagram of the contribution degree when the torque is 40N m and the rotational speeds are respectively 800rpm, 1200rpm, 1650rpm and 2200rpm, Fig. 8( b) shows the contribution statistics when the rotational speed is 1650rpm and the torque is 20N·m, 40N·m, 60N·m, 80N·m, and 100N·m, respectively;
具体实施方式detailed description
下面结合附图以及具体实施例对本发明作进一步的说明,但本发明的保护范围并不限于此。The present invention will be further described below with reference to the accompanying drawings and specific embodiments, but the protection scope of the present invention is not limited thereto.
如图1所示,SD2100TA柴油机工作循环中两个缸气门开启、关闭时刻示意图。气门开启、关闭激励也是振动响应信号的重要激励源,对多缸机而言,相邻缸气门开启、关闭激励振动响应和燃烧激励振动响应还可能存在耦合。由图可见,两个缸在着火前都不存在其他非燃烧激励信号。因此,气门开启、关闭激励与燃烧激励在时域上不存在同步性,其响应信号的特性可以忽略。As shown in Figure 1, the schematic diagram of the opening and closing timing of the two cylinder valves in the working cycle of the SD2100TA diesel engine. Valve opening and closing excitations are also important excitation sources of vibration response signals. For a multi-cylinder engine, there may also be coupling between the valve opening and closing excitation vibration responses and combustion excitation vibration responses of adjacent cylinders. It can be seen from the figure that there are no other non-combustion excitation signals before the ignition of the two cylinders. Therefore, there is no synchronization between the valve opening and closing excitation and the combustion excitation in the time domain, and the characteristics of the response signal can be ignored.
如图2所示,显示了2200r/min,60N·m工况时实测的振动加速度信号,两缸机工作过程是两个缸交替进行,工作相位上,1缸发火时刻滞后2缸180℃A,2缸发火时刻滞后1缸540℃A。振动信号从激起到恢复需要一定的时间,转速越高,各激励出现时刻的间隔时间越 短,为验证相邻缸的燃烧激励是否对当前缸燃烧激励响应信号产生影响,在额定转速下测取了缸盖表面振动加速度信号并进行了分析。图中A1是1缸缸盖表面振动加速度,A2是2缸缸盖表面振动加速度。理论上,两个缸测得的缸盖表面振动信号波形应当相同,只会在相位上出现固定偏差。从图2中可以看出,A1、A2具有相似的变化规律,在各自的燃烧上止点附近振动加速度信号能量最强。从图中还可以看出,A1在270-300℃A转角范围内,振动加速度信号已经恢复平稳,表明相邻缸燃烧激励并未影响到当前缸的燃烧过程。As shown in Figure 2, it shows the measured vibration acceleration signal at 2200 r/min and 60 N m. The working process of the two-cylinder engine is that the two cylinders are alternately performed. In the working phase, the firing time of cylinder 1 lags behind the firing time of cylinder 2 by 180 °C A , the ignition time of cylinder 2 lags behind cylinder 1 by 540℃A. It takes a certain amount of time for the vibration signal to recover from excitation. The higher the rotational speed, the shorter the interval between the occurrences of each excitation. In order to verify whether the combustion excitation of the adjacent cylinder has an impact on the combustion excitation response signal of the current cylinder, the measurement is performed at the rated speed. The vibration acceleration signal of the cylinder head surface is taken and analyzed. In the figure, A1 is the surface vibration acceleration of the 1-cylinder cylinder head, and A2 is the surface vibration acceleration of the 2-cylinder cylinder head. Theoretically, the vibration signal waveforms of the cylinder head surface measured by the two cylinders should be the same, with only a fixed deviation in phase. It can be seen from Figure 2 that A1 and A2 have similar variation laws, and the vibration acceleration signal energy is the strongest near their respective combustion top dead centers. It can also be seen from the figure that within the range of 270-300°C A1, the vibration acceleration signal has recovered to be stable, indicating that the combustion excitation of the adjacent cylinder does not affect the combustion process of the current cylinder.
如图3所示,转速为1200r/min,扭矩为40N·m时缸内压力信号及计算得到的往复惯性力激励曲线,该工况下往复惯性力峰值与缸内压力峰值的比值为5.3%。由于随着转速的增加,往复惯性力增大,缸压也会发生变化。经计算,相同扭矩下,800r/min,1200r/min,1650r/min,2200r/min时,往复惯性力峰值与缸内压力峰值的比值为2.4%,5.3%,10.6%,20.7%,呈增大趋势,由此可见,往复惯性力激励对振动信号的影响,尤其是在中高转速时的影响难以忽略不计。类似的,计算了活塞换向激励的强度,通过与缸内压力信号对比,认为,二者对振动信号的影响不能忽略,仍需进一步分析。As shown in Figure 3, when the rotational speed is 1200r/min and the torque is 40N m, the pressure signal in the cylinder and the calculated reciprocating inertial force excitation curve, the ratio of the peak value of the reciprocating inertial force to the peak pressure in the cylinder under this working condition is 5.3% . Since the reciprocating inertial force increases with the increase of the rotational speed, the cylinder pressure will also change. After calculation, under the same torque, at 800r/min, 1200r/min, 1650r/min, 2200r/min, the ratio of the peak value of the reciprocating inertial force to the peak value of the pressure in the cylinder is 2.4%, 5.3%, 10.6%, 20.7%, increasing. It can be seen that the influence of the reciprocating inertial force excitation on the vibration signal, especially at medium and high speeds, is difficult to ignore. Similarly, the strength of the piston commutation excitation is calculated, and by comparing with the pressure signal in the cylinder, it is believed that the influence of the two on the vibration signal cannot be ignored, and further analysis is needed.
如图4所示,为了定量研究各激励响应信号特性和耦合关系,本文对实际柴油机结构进行简化,保留机体和缸盖,利用ABAQUS软件建立了发动机仿真模型,计算中采用10节点四面体网格,机体通过4个支架固定在试验台架上,假定约束支架为刚性体不发生变形,对支架上各节点进行约束,限制其各个方向的自由度。As shown in Figure 4, in order to quantitatively study the characteristics and coupling relationship of each excitation response signal, this paper simplifies the actual diesel engine structure, retains the body and cylinder head, and uses the ABAQUS software to establish the engine simulation model. 10-node tetrahedral mesh is used in the calculation. , the body is fixed on the test bench through 4 brackets, assuming that the constraint bracket is a rigid body without deformation, the nodes on the bracket are constrained to limit the degrees of freedom in all directions.
如图5所示,显示了力锤激励信号、实测振动加速度信号、20点光顺后的实测振动加速度信号和仿真得到的振动加速度信号,为了对所建立发动机模型进行验证,本文在实际发动机测点2施加力锤激励,测取了力锤信号和测点1处的缸盖振动加速度信号;同时,将测取的力锤信号施加到模型测点2处,仿真得到测点1处的振动加速度信号。图中结果表明,仿真振动加速度信号和实测振动加速度信号在波形上有近似的趋势,幅值和相位基本一致,这表明所建立的发动机模型可用于分析各激励响应信号的特点和耦合关系。As shown in Figure 5, it shows the hammer excitation signal, the measured vibration acceleration signal, the measured vibration acceleration signal after 20 points of smoothing, and the vibration acceleration signal obtained by simulation. The force hammer excitation was applied at point 2, and the hammer signal and the cylinder head vibration acceleration signal at measuring point 1 were measured; at the same time, the measured force hammer signal was applied to the model measuring point 2, and the vibration at measuring point 1 was simulated. acceleration signal. The results in the figure show that the simulated vibration acceleration signal and the measured vibration acceleration signal have a similar trend in waveform, and the amplitude and phase are basically the same, which indicates that the established engine model can be used to analyze the characteristics and coupling relationship of each excitation response signal.
为了对激励信号经机体后的激励响应信号特性进行分析,本发明采用压缩小波变换技术对各激励振动响应信号的时频特点进行了分析。In order to analyze the characteristics of the excitation response signal after the excitation signal passes through the body, the invention adopts the compression wavelet transform technology to analyze the time-frequency characteristics of each excitation vibration response signal.
如图6所示,为往复惯性力信号的时频谱图,往复惯性力激励的能量主要分布在27.5Hz和55Hz周围,表明一阶和二阶往复惯性力起主要作用。As shown in Figure 6, which is the time-frequency spectrum of the reciprocating inertial force signal, the energy excited by the reciprocating inertial force is mainly distributed around 27.5Hz and 55Hz, indicating that the first-order and second-order reciprocating inertial forces play a major role.
如图7所示,为往复惯性力激励响应信号时频分析结果。图7(a)、(b)分别为振动位移和振动速度信号的时频分布,频域上振动位移信号和振动速度信号的时频分布和往复惯性力信号的时频分布相似,以一阶往复惯性力和二阶往复惯性力频率为中心向周围发散;图7(c)为振动加速度信号的时频分布,图7(d)是图7(c)的局部放大图,从这两张图中可以看出,振动 加速度信号在一阶往复惯性力和二阶往复惯性力频率周围也有较强的能量分布,除此之外,在1000Hz和2000Hz附近也有部分能量。由此可见,往复惯性力激励响应信号的能量在频带上分布较为集中,且贯穿整个内燃机工作周期。换而言之,往复惯性力激励响应信号和燃烧激励振动响应信号在时频域上都存在耦合。As shown in Figure 7, it is the time-frequency analysis result of the reciprocating inertial force excitation response signal. Figures 7(a) and (b) are the time-frequency distributions of the vibration displacement and vibration velocity signals, respectively. The time-frequency distributions of the vibration displacement signals and vibration velocity signals in the frequency domain are similar to the time-frequency distribution of the reciprocating inertial force signal. The frequencies of the reciprocating inertial force and the second-order reciprocating inertial force radiate from the center to the surroundings; Fig. 7(c) is the time-frequency distribution of the vibration acceleration signal, and Fig. 7(d) is a partial enlarged view of Fig. 7(c). From these two figures It can be seen that the vibration acceleration signal also has a strong energy distribution around the frequencies of the first-order reciprocating inertial force and the second-order reciprocating inertial force. In addition, there are also some energy around 1000Hz and 2000Hz. It can be seen that the energy of the reciprocating inertial force excitation response signal is relatively concentrated in the frequency band, and runs through the entire working cycle of the internal combustion engine. In other words, both the reciprocating inertial force excitation response signal and the combustion excitation vibration response signal are coupled in the time-frequency domain.
对往复惯性力激励和活塞换向激励引起的振动信号特性分析后认为,对振动响应信号起主要影响的是往复惯性力激励,为进一步分析往复惯性力激励对燃烧激励振动响应信号输出结果的影响,模拟得到了施加活塞换向撞击、往复惯性力和缸内压力信号三个激励共同作用时的振动速度信号,同时模拟得到了不考虑活塞换向撞击时其他两个激励共同作用时的振动速度信号,不考虑往复惯性力时其他两个激励共同作用时的振动速度信号。结果表明,各激励共同作用和燃烧激励单独作用时振动信号输出结果的差异主要是由往复惯性力激励引起的,也就是说,在所有的非燃烧激励中,往复惯性力是其首要作用的。After analyzing the characteristics of the vibration signal caused by the reciprocating inertial force excitation and the piston reversing excitation, it is believed that the main influence on the vibration response signal is the reciprocating inertial force excitation. In order to further analyze the influence of the reciprocating inertial force excitation on the output result of the combustion excitation vibration response signal , the simulation obtains the vibration velocity signal when the three excitations of piston reversing impact, reciprocating inertia force and cylinder pressure signal act together, and the simulation obtains the vibration velocity when the other two excitations work together without considering the piston reversing impact Signal, the vibration velocity signal when the other two excitations work together without considering the reciprocating inertial force. The results show that the difference of the vibration signal output when the excitations work together and the combustion excitation alone is mainly caused by the reciprocating inertial force excitation, that is to say, in all non-combustion excitations, the reciprocating inertial force is the primary effect.
本发明以振动速度信号为例,进一步研究关键频率处往复惯性力与振动信号间的相关性变化规律。本发明利用偏相干函数近似确定往复惯性力激励响应信号对振动信号的贡献。The present invention takes the vibration velocity signal as an example, and further studies the correlation variation law between the reciprocating inertial force and the vibration signal at the key frequency. The invention utilizes the partial coherence function to approximately determine the contribution of the reciprocating inertial force excitation response signal to the vibration signal.
两个输入x 1(t)、x 2(t),单输出信号y(t)的偏相干函数可以定义为: Two inputs x 1 (t), x 2 (t), the partial coherence function of a single output signal y(t) can be defined as:
Figure PCTCN2021071201-appb-000004
Figure PCTCN2021071201-appb-000004
Figure PCTCN2021071201-appb-000005
称为x 2(t)与y(t)之间的偏相干函数,其中x 1(t)的影响已排除掉。S 2Y1(f)为x 2(t)与y(t)的条件互功率谱,S 221(f)和S YY1(f)分别为去掉x 1(t)影响后x 2(t)与y(t)的条件自功率谱。且满足不等式:
Figure PCTCN2021071201-appb-000005
is called the partial coherence function between x 2 (t) and y(t), where the influence of x 1 (t) has been excluded. S 2Y1 (f) of x 2 (t) with the condition y (t) is the cross power spectrum, S 221 (f) and S YY1 (f) are removed x 1 post (t) Effect x 2 (t) and y (t) Conditional autopower spectrum. and satisfy the inequality:
Figure PCTCN2021071201-appb-000006
Figure PCTCN2021071201-appb-000006
如图8所示,显示了在二阶往复惯性力频率处,往复惯性力的贡献度统计结果。其中,图8(a)是扭矩为40N·m,转速分别为800rpm、1200rpm、1650rpm、2200rpm时的贡献度,图8(b)显示了转速为1650rpm,扭矩分别为20N·m、40N·m、60N·m、80N·m、100N·m时的贡献度。往复惯性力贡献度定义为:在某一频率处,单独施加往复惯性力时振动速度信号的频谱幅值与施加所有激励时振动速度信号幅值的比值。从图中可以看出,随着转速升高,往复惯性力对振动速度信号贡献度有明显增加,随着扭矩增加贡献度略有减小。这是由于往复惯性力与转速之间是二次函数关系,因此,随着转速增加,往复惯性力权重增加较大;而转速固定时,随着扭矩增加,往复惯性力保持不变,其他激励能量稍微增加,导致往复惯性力贡献度略有减小。As shown in Figure 8, the statistics of the contribution of the reciprocating inertial force at the second-order reciprocating inertial force frequency are shown. Among them, Figure 8(a) shows the contribution when the torque is 40N·m and the rotational speed is 800rpm, 1200rpm, 1650rpm, and 2200rpm, respectively. Figure 8(b) shows that the rotational speed is 1650rpm, and the torque is 20N·m and 40N·m, respectively. , 60N·m, 80N·m, 100N·m contribution. The contribution of the reciprocating inertial force is defined as: at a certain frequency, the ratio of the spectral amplitude of the vibration velocity signal when the reciprocating inertial force is applied alone to the amplitude of the vibration velocity signal when all excitations are applied. It can be seen from the figure that as the rotational speed increases, the contribution of the reciprocating inertial force to the vibration velocity signal increases significantly, and the contribution decreases slightly as the torque increases. This is because the relationship between the reciprocating inertial force and the rotational speed is a quadratic function. Therefore, as the rotational speed increases, the weight of the reciprocating inertial force increases greatly; while when the rotational speed is fixed, as the torque increases, the reciprocating inertial force remains unchanged, and other excitations A slight increase in energy results in a slight decrease in the contribution of the reciprocating inertial force.
所述实施例为本发明的优选的实施方式,但本发明并不限于上述实施方式,在不背离本发明的实质内容的情况下,本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。The embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above-mentioned embodiments, and any obvious improvement, replacement or All modifications belong to the protection scope of the present invention.

Claims (6)

  1. 一种基于偏相干函数求解非燃烧激励信号贡献度的方法,其特征在于,包括以下步骤:A method for solving non-combustion excitation signal contribution based on partial coherence function, characterized in that it comprises the following steps:
    步骤一、确定主要非燃烧激励源,分析气门开启关闭激励和相邻缸燃烧激励的时域特点,求解往复惯性力和活塞换向撞击力的激励源强度,得到对振动信号影响敏感的激励源;Step 1: Determine the main non-combustion excitation source, analyze the time domain characteristics of the valve opening and closing excitation and the adjacent cylinder combustion excitation, solve the excitation source strength of the reciprocating inertial force and the piston commutation impact force, and obtain the excitation source sensitive to the influence of vibration signals. ;
    步骤二、确定对振动信号起首要作用的非燃烧激励,对实际内燃机结构进行简化,保留机体和缸盖,建立发动机仿真模型;在模型上选定测点施加力锤,通过对比仿真振动信号和实测振动信号验证模型有效性;Step 2: Determine the non-combustion excitation that plays a primary role in the vibration signal, simplify the actual internal combustion engine structure, retain the body and cylinder head, and establish an engine simulation model; select the measuring point on the model to apply a force hammer, and compare the simulation vibration signal and Measured vibration signal to verify the validity of the model;
    通过压缩小波变换技术对各激励振动响应信号的时频特点进行了分析,得到对振动信号起首要作用的非燃烧激励;The time-frequency characteristics of each excitation vibration response signal are analyzed by compressive wavelet transform technology, and the non-combustion excitation which plays a primary role in the vibration signal is obtained;
    步骤三、根据步骤二中确定的首要非燃烧激励,利用偏相干函数计算该激励信号在特定频率处对振动信号的贡献度。Step 3: According to the primary non-combustion excitation determined in Step 2, use the partial coherence function to calculate the contribution of the excitation signal to the vibration signal at a specific frequency.
  2. 根据权利要求1所述的基于偏相干函数求解非燃烧激励信号贡献度的方法,其特征在于,所述步骤一中求解往复惯性力和活塞换向撞击力的激励源强度的方法为:根据测得的缸内压力信号、往复惯性力的激励曲线和活塞换向撞击力的激励曲线分别获取其最大值。The method for solving the contribution degree of non-combustion excitation signal based on partial coherence function according to claim 1, characterized in that, in the step 1, the method for solving the excitation source strength of reciprocating inertial force and piston reversing impact force is: The obtained pressure signal in the cylinder, the excitation curve of the reciprocating inertial force and the excitation curve of the piston reversing impact force are obtained to obtain their maximum values respectively.
  3. 根据权利要求1所述的基于偏相干函数求解非燃烧激励信号贡献度的方法,其特征在于:所述步骤二中建立发动机仿真模型时采用10节点四面体网格,机体通过4个支架固定在试验台架上,约束支架为刚性体不发生变形,对支架上各节点进行约束,限制其各个方向的自由度。The method for calculating the contribution of non-combustion excitation signal based on partial coherence function according to claim 1, characterized in that: in the step 2, a 10-node tetrahedral mesh is used when establishing the engine simulation model, and the body is fixed on the On the test bench, the restraint bracket is a rigid body that does not deform, and each node on the bracket is constrained to limit its degrees of freedom in all directions.
  4. 根据权利要求1所述的基于偏相干函数求解非燃烧激励信号贡献度的方法,其特征在于,所述步骤二中仿真振动信号和实测振动信号模型有效性验证方法为:The method for solving non-combustion excitation signal contribution degree based on partial coherence function according to claim 1, it is characterized in that, in described step 2, the simulation vibration signal and actual measurement vibration signal model validity verification method is:
    首先、在实际发动机上选择施力点施加力锤激励,测取力锤信号和测点处的振动加速度信号;First, select the force application point on the actual engine to apply the force hammer excitation, and measure the force hammer signal and the vibration acceleration signal at the measurement point;
    然后,将测取的力锤信号施加到模型的相应的施力点上,仿真得到测点处的振动加速度信号;Then, apply the measured hammer signal to the corresponding force application point of the model, and simulate the vibration acceleration signal at the measurement point;
    最后、得出力锤激励信号、实测振动加速度信号、20点光顺后的实测振动加速度信号和仿真得到的振动加速度信号图,得出验证结论。Finally, the excitation signal of the hammer, the measured vibration acceleration signal, the measured vibration acceleration signal after 20 points of smoothing, and the vibration acceleration signal diagram obtained by simulation are obtained, and the verification conclusion is drawn.
  5. 根据权利要求1所述的基于偏相干函数求解非燃烧激励信号贡献度的方法,其特征在于:所述步骤二中通过压缩小波变换技术对各激励振动响应信号的时频特点进行了分析的方法为:The method for calculating the contribution degree of non-combustion excitation signals based on partial coherence function according to claim 1, characterized in that: in the second step, the time-frequency characteristics of each excitation vibration response signal are analyzed by compressive wavelet transform technology. for:
    首先、测得往复惯性力信号和活塞换向撞击力的时频谱图,得出能量主要分布的频率范围;First, measure the time-frequency spectrum of the reciprocating inertial force signal and the reversing impact force of the piston, and obtain the frequency range of the main energy distribution;
    然后、对激励响应信号时频分析,得到振动位移、振动速度和振动加速度信号的时频分布,得到单独施加往复惯性力激励和活塞换向激励引起的振动信号特性;Then, the time-frequency analysis of the excitation response signal is performed to obtain the time-frequency distribution of the vibration displacement, vibration velocity and vibration acceleration signals, and the vibration signal characteristics caused by the reciprocating inertial force excitation and the piston commutation excitation alone are obtained;
    最后、模拟施加活塞换向撞击、往复惯性力和缸内压力信号三个激励共同作用时的振动速度信号,同时模拟不考虑活塞换向撞击时其他两个激励共同作用时的振动速度信号,不考虑往复惯性力时其他两个激励共同作用时的振动速度信号;得到在所有的非燃烧激励中起首 要作用的非燃烧激励。Finally, simulate the vibration speed signal when the three excitations of piston reversing impact, reciprocating inertial force and cylinder pressure signal work together, and simulate the vibration speed signal when the other two excitations work together without considering the piston reversing impact. Considering the reciprocating inertial force, the vibration velocity signal when the other two excitations act together; obtain the non-combustion excitation that plays the primary role among all the non-combustion excitations.
  6. 根据权利要求1所述的基于偏相干函数求解非燃烧激励信号贡献度的方法,其特征在于,利用偏相干函数计算该激励信号在特定频率处对振动信号的贡献度的方法为:The method for calculating the contribution degree of a non-combustion excitation signal based on a partial coherence function according to claim 1, wherein the method for calculating the contribution degree of the excitation signal to the vibration signal at a specific frequency by using the partial coherence function is:
    两个输入x 1(t)、x 2(t),单输出信号y(t)的偏相干函数可以定义为: Two inputs x 1 (t), x 2 (t), the partial coherence function of a single output signal y(t) can be defined as:
    Figure PCTCN2021071201-appb-100001
    Figure PCTCN2021071201-appb-100001
    Figure PCTCN2021071201-appb-100002
    称为x 2(t)与y(t)之间的偏相干函数,其中x 1(t)的影响已排除掉,S 2Y1(f)为x 2(t)与y(t)的条件互功率谱,S 221(f)和S YY1(f)分别为去掉x 1(t)影响后x 2(t)与y(t)的条件自功率谱,且满足不等式:
    Figure PCTCN2021071201-appb-100002
    It is called the partial coherence function between x 2 (t) and y(t), where the influence of x 1 (t) has been excluded, and S 2Y1 (f) is the conditional relationship between x 2 (t) and y(t) power spectrum, S 221 (f) and S YY1 (f) are removed after x 1 (t) Effect x 2 (t) and y (t) from the power spectrum of conditions, and satisfies the inequality:
    Figure PCTCN2021071201-appb-100003
    Figure PCTCN2021071201-appb-100003
    根据计算结果得出特定频率时贡献度的统计结果,最终得出结论。According to the calculation results, the statistical results of the contribution at a specific frequency are obtained, and finally a conclusion is drawn.
PCT/CN2021/071201 2020-07-07 2021-01-12 Method for solving degree of contribution of non-combustion excitation signal on basis of partial coherence function WO2022007382A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202010647022.2 2020-07-07
CN202010647022.2A CN111964913B (en) 2020-07-07 2020-07-07 Method for solving contribution degree of non-combustion excitation signal based on partial coherence function

Publications (1)

Publication Number Publication Date
WO2022007382A1 true WO2022007382A1 (en) 2022-01-13

Family

ID=73361289

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2021/071201 WO2022007382A1 (en) 2020-07-07 2021-01-12 Method for solving degree of contribution of non-combustion excitation signal on basis of partial coherence function

Country Status (2)

Country Link
CN (1) CN111964913B (en)
WO (1) WO2022007382A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116306143A (en) * 2023-03-17 2023-06-23 中国航空发动机研究院 Double-rotor system unbalance fault simulation method and electronic equipment

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111964913B (en) * 2020-07-07 2021-09-10 江苏大学 Method for solving contribution degree of non-combustion excitation signal based on partial coherence function
CN113642204B (en) * 2021-07-06 2024-03-19 江苏大学 Method for correcting combustion starting point identification deviation based on combustion excitation contribution degree
CN113504143A (en) * 2021-07-08 2021-10-15 潍柴动力股份有限公司 Method and device for diagnosing abrasion of valve seat ring

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1194690A (en) * 1997-09-17 1999-04-09 Isuzu Motors Ltd Exciting device for simulated combustion
CN101839805A (en) * 2010-03-19 2010-09-22 同济大学 Method for quality inspection of active fault and diagnosis of intelligent fault of engine
CN103632017A (en) * 2013-12-24 2014-03-12 山东大学 Pattern recognition-based method for improving signal-to-noise ratio of vibration signals on surface of internal combustion engine
DE102014204154A1 (en) * 2014-03-06 2015-09-10 Bayerische Motoren Werke Aktiengesellschaft Test bench for testing a motor attachment
CN107101829A (en) * 2017-04-11 2017-08-29 西北工业大学 A kind of intelligent diagnosing method of aero-engine structure class failure
CN111964913A (en) * 2020-07-07 2020-11-20 江苏大学 Method for solving contribution degree of non-combustion excitation signal based on partial coherence function

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1194690A (en) * 1997-09-17 1999-04-09 Isuzu Motors Ltd Exciting device for simulated combustion
CN101839805A (en) * 2010-03-19 2010-09-22 同济大学 Method for quality inspection of active fault and diagnosis of intelligent fault of engine
CN103632017A (en) * 2013-12-24 2014-03-12 山东大学 Pattern recognition-based method for improving signal-to-noise ratio of vibration signals on surface of internal combustion engine
DE102014204154A1 (en) * 2014-03-06 2015-09-10 Bayerische Motoren Werke Aktiengesellschaft Test bench for testing a motor attachment
CN107101829A (en) * 2017-04-11 2017-08-29 西北工业大学 A kind of intelligent diagnosing method of aero-engine structure class failure
CN111964913A (en) * 2020-07-07 2020-11-20 江苏大学 Method for solving contribution degree of non-combustion excitation signal based on partial coherence function

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIANG XINGYU, SHU GEQUN,WANG GANGZHI, WANG YANGJUN: "A Study on 3D Vibration of the Crankshaft and Surface Vibration of the Crankcase based on Partial Coherence Function Analysis", AUTOMOTIVE ENGINEERING, vol. 28, no. 3, 25 March 2006 (2006-03-25), pages 271 - 275, XP055885550, DOI: 10.19562/j.chinasae.qcgc.2006.03.014 *
WANG YANG, CHENG YONG;JI SHAOBO;LAN XIN: "Coupling Relationship Analysis on the Vibration Signals of Cylinder-Head Surface", HUAZHONG KEJI DAXUE XUEBAO (ZIRAN KEXUE BAN)/JOURNAL OF HUAZHONG UNIVERSITY OF SCIENCE AND TECHNOLOGY (NATURAL SCIENCE EDITION), HUAZHONG KEJI DAXUE, CN, vol. 45, no. 4, 30 April 2017 (2017-04-30), CN , pages 78 - 83, XP055885544, ISSN: 1671-4512, DOI: 10.13245/j.hust.170415 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116306143A (en) * 2023-03-17 2023-06-23 中国航空发动机研究院 Double-rotor system unbalance fault simulation method and electronic equipment
CN116306143B (en) * 2023-03-17 2023-09-22 中国航空发动机研究院 Double-rotor system unbalance fault simulation method and electronic equipment

Also Published As

Publication number Publication date
CN111964913A (en) 2020-11-20
CN111964913B (en) 2021-09-10

Similar Documents

Publication Publication Date Title
WO2022007382A1 (en) Method for solving degree of contribution of non-combustion excitation signal on basis of partial coherence function
Yang et al. Fault detection in a diesel engine by analysing the instantaneous angular speed
Albarbar et al. Diesel engine fuel injection monitoring using acoustic measurements and independent component analysis
JP4871962B2 (en) Real-time estimation from vibration signals of engine combustion parameters
Vulli et al. Time-frequency analysis of single-point engine-block vibration measurements for multiple excitation-event identification
El-Ghamry et al. Indirect measurement of cylinder pressure from diesel engines using acoustic emission
Pruvost et al. Diesel engine combustion and mechanical noise separation using an improved spectrofilter
CN104763532A (en) Method and apparatus for recognizing knocking of an internal combustion engine, preferably of a gasoline engine
Hu et al. Multivariate statistical analysis strategy for multiple misfire detection in internal combustion engines
Han et al. Recursive engine in-cylinder pressure estimation using Kalman filter and structural vibration signal
Chen et al. Misfire detection in a dynamic skip fire engine
Al-Durra et al. A model-based estimator of engine cylinder pressure imbalance for combustion feedback control applications
Jafari et al. In-cylinder pressure reconstruction by engine acoustic emission
Mocanu et al. Estimation of main combustion parameters from the measured instantaneous crankshaft speed
Narayan Time-frequency analysis of Diesel engine noise
d'Ambrosio et al. Time-frequency analysis application to the evaluation of instantaneous combustion noise
Pan et al. Identification of start of combustion based on contribution level of principal component in vibration signal of cylinder head surface
Larsson et al. SI-engine cylinder pressure estimation using torque sensors
CN113642204B (en) Method for correcting combustion starting point identification deviation based on combustion excitation contribution degree
Anjum et al. Cyclic torque imbalance detection in gasoline engines using second order sliding mode
Narayan Correlation between cylinder pressure and noise emissions from diesel engines
Saxén et al. Power balancing of internal combustion engines–a time and frequency domain analysis
Wei et al. Prediction and analysis of diesel engine combustion noise using transfer function method
Grządziela et al. High and medium speed marine diesel engines injection installation vibration diagnostics
Grondin et al. Combustion parameters estimation and control using vibration signal: Application to the diesel HCCI engine

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21837055

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21837055

Country of ref document: EP

Kind code of ref document: A1