CN110686892A - 在线柴油机弹性传动齿轮状态检测方法 - Google Patents
在线柴油机弹性传动齿轮状态检测方法 Download PDFInfo
- Publication number
- CN110686892A CN110686892A CN201911010125.1A CN201911010125A CN110686892A CN 110686892 A CN110686892 A CN 110686892A CN 201911010125 A CN201911010125 A CN 201911010125A CN 110686892 A CN110686892 A CN 110686892A
- Authority
- CN
- China
- Prior art keywords
- diesel engine
- rotating speed
- vibration
- instantaneous
- value
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/04—Testing internal-combustion engines
- G01M15/05—Testing internal-combustion engines by combined monitoring of two or more different engine parameters
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明涉及一种在线柴油机弹性传动齿轮状态检测方法,属于内燃机及海洋设备技术领域。
背景技术
柴油机弹性齿轮检测实际上技术核心是特征提取及特征分析,是智能柴油机发展的主要方向。
现代大型柴油机是一个包含材料学,机械学,电子控制等多个多学科为一体的复杂的非线性系统,它包含多个子系统,且个各子系统相互作用,柴油机在运行期间是由机械结构和电控单元协调个部件按一定时序工作,然而随着机械结构的磨损或电控单元的不准确或不稳定,就导致各系统不能及时有效的工作,系统间又是相互耦合的,一个系统的不正常,会影响到其他系统,表现出其它的故障现象。导致故障的起因和故障表现不能一一对应,出现一因多果、和一果多因的现象,所以很多情况都是柴油机某些零部件彻底损坏后才会被发现或定位,此时已经造成比较大的损失。柴油机弹性传动齿轮,与普通齿轮相比,其内部有供缓冲的弹簧,加之柴油机运行时缸内爆炸振动大,现有针对齿轮的振动故障诊断此时是失效的。所以柴油机弹性传动齿轮状态的管理从“事后维修”需要向“事前维护”方向转变。
发明内容
本发明的目的在于提供一种在线柴油机弹性传动齿轮状态检测方法,实现对弹性齿轮微小故障的检测,从而提醒用户进行提前检修更换从而减少更大的损失。
本发明的目的通过以下技术方案予以实现:
一种在线柴油机弹性传动齿轮状态检测方法,包括:
1、通过传感器对柴油机的瞬时转速、齿轮轴振动、功率负荷进行监测;
2、步骤1采集的信号转化成标准电信号;
3、由数据采集卡对标准电信号进行采集及数模转换;
4、通过对柴油机瞬时转速及气缸振动的耦合,将其中振动在时域内转化到与曲轴转角对应的角度域内进行阶比分析;
6、带入公式将柴油机功率P kW、转速n RMP变成裕度系数
h为常数,根据各机型而定;
7、测取任意状态下阶比数据矩阵B,与ηA对比,当B>ηA时,弹性传动齿轮存在异常。
本发明的目的还可以通过以下技术措施来进一步实现:
前述在线柴油机弹性传动齿轮状态检测方法,其中步骤4通过以下完成:
1)瞬时转速异常值剔除
柴油机瞬时转速测量时,测量的是由信号调理板转化后的脉冲频率,然而在实际测量中,由于传感器的振动或者调理板受到电磁干扰,会出现脉冲丢失或者脉冲增多现象,表现在数据上面,就是在两个上止点信号之间的瞬时转速值或周期值不是齿轮的齿数;然而增加的脉冲周期值和柴油机转速无关,通常表现为特别大的异常值,有这些原因的存在,使得直接使用测量数据进行阶比跟踪会对计算精度产生很大影响;剔除明显异常值后,使用差分统计寻找其它异常值时,首先将瞬时转速后向差分
X(i)=x(i)-x(i+1) (1)
x(i)为第i个瞬时转速值,X(i)为第i与第i-1齿之间的转速变化量,因为转速的变化不可能无限制大的,正常情况下符合正态分布,剔除异常值属于单正态总体方差置信区间问题,设总体X~N(μ,σ2),其中μ,σ2未知,X1,X2,…,Xn是取自总体X的一个样本.求方差σ2的置信度为1-α的置信区间,σ2的无偏估计为S2,有
对给定的置信水平1-α,由
于是方差σ2的1-α置信区间为
而方差σ的1-α置信区间
当找出异常点后,要将异常点所占用的时间合并给其后一个正常测量值,并将所有值前移一个,设x(i)为异常值,还原后x(i)为
2)基于重采样的计算阶比
将振动信号从时域转化到角度域内是为了准确定位和获取故障段内信号,是特征提取的第一步,它是从等时间间隔采样的数据出发进行数字重采样获得到带有角度标识的等角度间隔数据的过程,即计算阶比跟踪;使用的方法是重采样,其流程为:
(1)利用上至点信号截取出720度内的瞬时转速数据N和振动信号X(t);
(2)由瞬时转速Ni求出每个脉冲间隔Ti,其中i=1,2,3…269.,
式中z-----齿轮盘齿数;
(3)求出每个齿角度段内的瞬时周期值Γi,即每个齿段内柴油机转过单位度数所需时间
(4)确定每个齿的中间曲轴转角度数Ri;
(6)根据振动信号X(t)的采样频率求出采样间隔时间
t=1/f (10)
式中f-----振动采样频率;
k=T0/t,m=T0%t (12)
与现有技术相比,本发明的有益效果是:本发明的在线柴油机弹性传动齿轮状态检测方法,将振动信号从时域转化到角度域内是为了准确定位和获取故障段内信号,实现对弹性齿轮微小故障的检测,从而提醒用户进行提前检修更换从而减少更大的损失。
附图说明
图1是本发明数据采集及传递示意图;
图2是柴油机一工作循环瞬时转速测量值图;
图3是本发明重采样计算流程图;
图4是本发明系统数据处理流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
本发明在线柴油机弹性传动齿轮状态检测方法的数据采集及传递如图1所示,该方法包括:
1、通过传感器对柴油机的瞬时转速、齿轮轴振动、功率负荷进行监测。
2、其中瞬时转速监测,通过MP-981磁电式感器检测飞轮齿圈,其输出频率与转速相关的正弦波信号,通过史密斯触发器转化成脉冲信号,随后输入数据采集卡。齿轮轴振动监测,通过安装在齿轮轴端部的GST压电式加速度传感器型号:CA-YD-107CM,输出电荷信号,直接接入数据采集卡。功率负荷监测通过安装于调速器输出轴上的位移传感器,将油门的位移通过0-40mA送入数据采集卡。
3、由数据采集卡对标准电信号进行采集及数模转换;采用市面上较成熟的产品NI公司的PXI系类8106,主机是PXI系类1050,模拟输入数据采集卡采用6251用于温度压力的采集,脉冲定时器采用6602用于瞬时转速采集,4472用于振动加速度采集;
4、通过对柴油机瞬时转速及气缸振动的耦合,将其中振动在时域内转化到与曲轴转角对应的角度域内进行阶比分析;通过以下步骤完成:
4.1瞬时转速异常值剔除
柴油机瞬时转速测量时,测量的是由信号调理板转化后的脉冲频率,然而在实际测量中,由于传感器的振动或者调理板受到电磁干扰,会出现脉冲丢失或者脉冲增多现象,表现在数据上面,就是在两个上止点信号之间的瞬时转速值或周期值不是齿轮的齿数;然而增加的脉冲周期值和柴油机转速无关,通常表现为特别大的异常值,有这些原因的存在,使得直接使用测量数据进行阶比跟踪会对计算精度产生很大影响;如图2所示是柴油机一个工作循环内的瞬时转速测量值,实验中测量齿轮盘为148个齿应该得到296个值,实际测量中共得到301个测量值,多出了5个异常值,由于柴油机转速不会突变,所以可以很明显看出三个异常值,除了这三个外,还要剔除两个异常值,在寻找其它异常值时,使用差分统计;首先将瞬时转速后向差分
X(i)=x(i)-x(i+1) (4-1)
x(i)为第i个瞬时转速值,X(i)为第i与第i-1齿之间的转速变化量,因为转速的变化不可能无限制大的,正常情况下符合正态分布,剔除异常值属于单正态总体方差置信区间问题,设总体X~N(μ,σ2),其中μ,σ2未知,X1,X2,…,Xn是取自总体X的一个样本.求方差σ2的置信度为1-α的置信区间,σ2的无偏估计为S2,有
对给定的置信水平1-α,由
于是方差σ2的1-α置信区间为
而方差σ的1-α置信区间
对于置信区间1-α的选择,可以根据异常值个数确定,对于上文中除去两个异常值时,可选1-α=0.99,此时不在置信区间的两个即有99%的概率是异常点。
当找出异常点后,要将异常点所占用的时间合并给其后一个正常测量值,并将所有值前移一个,设x(i)为异常值,还原后x(i)为
4.2基于重采样的计算阶比
将振动信号从时域转化到角度域内是为了准确定位和获取故障段内信号,是特征提取的第一步,它是从等时间间隔采样的数据出发进行数字重采样获得到带有角度标识的等角度间隔数据的过程,即计算阶比跟踪;使用的方法就是重采样,其对运算量的要求较大,如图3为计算过程:
下面为根据具体情况设计的算法将插值和抽取合并唯一,其流程为:
1)利用上至点信号截取出720度内的瞬时转速数据N和振动信号X(t);
2)由瞬时转速Ni求出每个脉冲间隔Ti,其中i=1,2,3…269.,
式中z-----齿轮盘齿数;
3)求出每个齿角度段内的瞬时周期值Γi,即每个齿段内柴油机转过单位度数所需时间
4)确定每个齿的中间曲轴转角度数Ri;
对横坐标为R纵坐标为Γ做三次样条插值,得出新的坐标数组和
6)根据振动信号X(t)的采样频率求出采样间隔时间
t=1/f (4-10)
式中f-----振动采样频率;
k=T0/t,m=T0%t (4-12)
6、带入公式将柴油机功率P kW、转速n RMP变成裕度系数
h为常数,根据各机型而定。
7、测取任意状态下阶比数据矩阵B,与ηA对比,当B>ηA时,弹性传动齿轮存在异常。
本发明整个系统数据处理流程如图4所示。
除上述实施例外,本发明还可以有其他实施方式,凡采用等同替换或等效变换形成的技术方案,均落在本发明要求的保护范围内。
Claims (2)
2.如权利要求1所述的在线柴油机弹性传动齿轮状态检测方法,其特征在于,所述步骤4通过以下完成:
1)瞬时转速异常值剔除
柴油机瞬时转速测量时,测量的是由信号调理板转化后的脉冲频率,然而在实际测量中,由于传感器的振动或者调理板受到电磁干扰,会出现脉冲丢失或者脉冲增多现象,表现在数据上面,就是在两个上止点信号之间的瞬时转速值或周期值不是齿轮的齿数;然而增加的脉冲周期值和柴油机转速无关,通常表现为特别大的异常值,有这些原因的存在,使得直接使用测量数据进行阶比跟踪会对计算精度产生很大影响;剔除明显异常值后,使用差分统计寻找其它异常值时,首先将瞬时转速后向差分
X(i)=x(i)-x(i+1) (1)
x(i)为第i个瞬时转速值,X(i)为第i与第i-1齿之间的转速变化量,因为转速的变化不可能无限制大的,正常情况下符合正态分布,剔除异常值属于单正态总体方差置信区间问题,设总体X~N(μ,σ2),其中μ,σ2未知,X1,X2,…,Xn是取自总体X的一个样本.求方差σ2的置信度为1-α的置信区间,σ2的无偏估计为S2,有
对给定的置信水平1-α,由
于是方差σ2的1-α置信区间为
而方差σ的1-α置信区间
当找出异常点后,要将异常点所占用的时间合并给其后一个正常测量值,并将所有值前移一个,设x(i)为异常值,还原后x(i)为
2)基于重采样的计算阶比
将振动信号从时域转化到角度域内是为了准确定位和获取故障段内信号,是特征提取的第一步,它是从等时间间隔采样的数据出发进行数字重采样获得到带有角度标识的等角度间隔数据的过程,即计算阶比跟踪;使用的方法是重采样,其流程为:
(1)利用上至点信号截取出720度内的瞬时转速数据N和振动信号X(t);
(2)由瞬时转速Ni求出每个脉冲间隔Ti,其中i=1,2,3…269.,
式中z-----齿轮盘齿数;
(3)求出每个齿角度段内的瞬时周期值Γi,即每个齿段内柴油机转过单位度数所需时间
(4)确定每个齿的中间曲轴转角度数Ri;
(6)根据振动信号X(t)的采样频率求出采样间隔时间
t=1/f (10)
式中f-----振动采样频率;
k=T0/t,m=T0%t (12)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911010125.1A CN110686892B (zh) | 2019-10-23 | 2019-10-23 | 在线柴油机弹性传动齿轮状态检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911010125.1A CN110686892B (zh) | 2019-10-23 | 2019-10-23 | 在线柴油机弹性传动齿轮状态检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110686892A true CN110686892A (zh) | 2020-01-14 |
CN110686892B CN110686892B (zh) | 2021-07-13 |
Family
ID=69113883
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911010125.1A Active CN110686892B (zh) | 2019-10-23 | 2019-10-23 | 在线柴油机弹性传动齿轮状态检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110686892B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114459766A (zh) * | 2022-01-13 | 2022-05-10 | 武汉理工大学 | 一种海洋平台原油发电机组油头工作状态监测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175439A (zh) * | 2011-01-25 | 2011-09-07 | 上海交通大学 | 针对旋转机械的阶次分析实现方法 |
CN102998110A (zh) * | 2012-11-29 | 2013-03-27 | 西安交通大学 | 一种基于阶比-全息谱原理的旋转机械故障特征提取方法 |
CN105157914A (zh) * | 2015-09-06 | 2015-12-16 | 山东大学 | 一种内燃机缸内压力信号时域转角度域的系统及方法 |
CN105510023A (zh) * | 2015-11-24 | 2016-04-20 | 国网内蒙古东部电力有限公司电力科学研究院 | 基于散度指标的变工况风电行星齿轮箱故障诊断方法 |
CN107941513A (zh) * | 2017-12-25 | 2018-04-20 | 北京建筑大学 | 一种列车走行部轴承非平稳运维的时频阶比跟踪方法 |
JP6417824B2 (ja) * | 2014-09-26 | 2018-11-07 | 沖電気工業株式会社 | 偏波依存性損失エミュレータ及び偏波依存性損失エミュレート方法 |
CN109682597A (zh) * | 2019-02-01 | 2019-04-26 | 朗坤智慧科技股份有限公司 | 一种齿轮箱振动信号处理和分析方法 |
-
2019
- 2019-10-23 CN CN201911010125.1A patent/CN110686892B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175439A (zh) * | 2011-01-25 | 2011-09-07 | 上海交通大学 | 针对旋转机械的阶次分析实现方法 |
CN102998110A (zh) * | 2012-11-29 | 2013-03-27 | 西安交通大学 | 一种基于阶比-全息谱原理的旋转机械故障特征提取方法 |
JP6417824B2 (ja) * | 2014-09-26 | 2018-11-07 | 沖電気工業株式会社 | 偏波依存性損失エミュレータ及び偏波依存性損失エミュレート方法 |
CN105157914A (zh) * | 2015-09-06 | 2015-12-16 | 山东大学 | 一种内燃机缸内压力信号时域转角度域的系统及方法 |
CN105157914B (zh) * | 2015-09-06 | 2017-06-30 | 山东大学 | 一种内燃机缸内压力信号时域转角度域的系统及方法 |
CN105510023A (zh) * | 2015-11-24 | 2016-04-20 | 国网内蒙古东部电力有限公司电力科学研究院 | 基于散度指标的变工况风电行星齿轮箱故障诊断方法 |
CN107941513A (zh) * | 2017-12-25 | 2018-04-20 | 北京建筑大学 | 一种列车走行部轴承非平稳运维的时频阶比跟踪方法 |
CN109682597A (zh) * | 2019-02-01 | 2019-04-26 | 朗坤智慧科技股份有限公司 | 一种齿轮箱振动信号处理和分析方法 |
Non-Patent Citations (3)
Title |
---|
卫耿源: "70B变速箱阶比谱知识库特征分析与研究", 《中国优秀硕士学位论文全文数据库》 * |
蒋大伟: "基于多传感器信息融合的柴油机配气机构故障诊断技术研究", 《中国优秀硕士学位论文全文数据库》 * |
费红姿 等: "柴油机气门故障特征提取方法研究", 《内燃机工程》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114459766A (zh) * | 2022-01-13 | 2022-05-10 | 武汉理工大学 | 一种海洋平台原油发电机组油头工作状态监测方法 |
CN114459766B (zh) * | 2022-01-13 | 2024-06-11 | 武汉理工大学 | 一种海洋平台原油发电机组油头工作状态监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110686892B (zh) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US4215404A (en) | Automatic device for diagnostic checkup of vehicles | |
CN110686890B (zh) | 在线柴油机气阀状态检测方法 | |
CN104819841B (zh) | 基于内置编码信息的单传感柔性角度域平均方法 | |
CN108015797B (zh) | 一种rv减速机传动误差在线监测方法 | |
CN110285045B (zh) | 一种核电厂rcp主泵振动相位监测系统和方法 | |
CN112539866B (zh) | 一种基于视觉深度学习的螺栓轴力监测系统及监测方法 | |
CN108692711B (zh) | 基于低空探空火箭的海洋数据处理的实现方法 | |
CN110686892B (zh) | 在线柴油机弹性传动齿轮状态检测方法 | |
CN110470382A (zh) | 一种无转速同步传感器下的叶片自动定位方法 | |
CN110686830B (zh) | 在线柴油机活塞环状态检测方法 | |
US5311123A (en) | Method of measuring the instantaneous shaft velocity of a rotary machine | |
CN112381860A (zh) | 一种旋转叶片动频测量的无标记计算机视觉方法 | |
CN109883692B (zh) | 基于内置编码器信息的广义差分滤波方法 | |
US6082187A (en) | Method for detecting a power loss condition of a reciprocating internal combustion engine | |
CN117972600A (zh) | 一种基于多维故障特征学习的风电机组关键部件异常探测方法 | |
CN103760376A (zh) | 一种基于振动原理的发动机转速测量仪及其测试方法 | |
CN110686879B (zh) | 在线柴油机气缸套状态检测方法 | |
CN105090084B (zh) | 风机在线监测系统及其监测方法 | |
CN115859210A (zh) | 一种带电机设备停机尾音智能检测方法及系统 | |
CN112857806B (zh) | 一种基于移动窗口时域特征提取的轴承故障检测方法 | |
CN112393907B (zh) | 基于扫频分析技术的风电机组轴承典型故障自动诊断方法 | |
CN114166540A (zh) | 一种基于虚拟仪器的滚装船动力装置综合监测系统 | |
CN111272104A (zh) | 适用于极端环境下自动屏蔽位置传感器故障读数头的方法 | |
CN118626752B (zh) | 一种变计数时钟的转速采样计算方法 | |
RU2078324C1 (ru) | Способ определения технического состояния двигателей внутреннего сгорания и экспертная система для его осуществления |
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 | ||
CP03 | Change of name, title or address |
Address after: 212211 No. 402, Changjiang Road, Zhenjiang City, Jiangsu Province Patentee after: China Shipbuilding Power Zhenjiang Co.,Ltd. Address before: 402 Changjiang Road, Runzhou District, Zhenjiang City, Jiangsu Province Patentee before: CSSC MARINE POWER Co.,Ltd. |
|
CP03 | Change of name, title or address |