CN112284428B - 活塞发动机曲轴停机相位在线监测方法 - Google Patents

活塞发动机曲轴停机相位在线监测方法 Download PDF

Info

Publication number
CN112284428B
CN112284428B CN202010970888.7A CN202010970888A CN112284428B CN 112284428 B CN112284428 B CN 112284428B CN 202010970888 A CN202010970888 A CN 202010970888A CN 112284428 B CN112284428 B CN 112284428B
Authority
CN
China
Prior art keywords
signal
gear
crankshaft
rotation angle
formula
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
Application number
CN202010970888.7A
Other languages
English (en)
Other versions
CN112284428A (zh
Inventor
茆志伟
张进杰
江志农
王晨光
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Chemical Technology
Original Assignee
Beijing University of Chemical Technology
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 Beijing University of Chemical Technology filed Critical Beijing University of Chemical Technology
Priority to CN202010970888.7A priority Critical patent/CN112284428B/zh
Publication of CN112284428A publication Critical patent/CN112284428A/zh
Application granted granted Critical
Publication of CN112284428B publication Critical patent/CN112284428B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/12Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable using electric or magnetic means
    • G01D5/244Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable using electric or magnetic means influencing characteristics of pulses or pulse trains; generating pulses or pulse trains
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02BINTERNAL-COMBUSTION PISTON ENGINES; COMBUSTION ENGINES IN GENERAL
    • F02B77/00Component parts, details or accessories, not otherwise provided for
    • F02B77/08Safety, indicating, or supervising devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P3/00Measuring linear or angular speed; Measuring differences of linear or angular speeds
    • G01P3/42Devices characterised by the use of electric or magnetic means
    • G01P3/44Devices characterised by the use of electric or magnetic means for measuring angular speed
    • G01P3/48Devices characterised by the use of electric or magnetic means for measuring angular speed by measuring frequency of generated current or voltage
    • G01P3/481Devices characterised by the use of electric or magnetic means for measuring angular speed by measuring frequency of generated current or voltage of pulse signals
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Shafts, Cranks, Connecting Bars, And Related Bearings (AREA)

Abstract

活塞发动机曲轴停机相位在线监测方法属于设备故障监测与诊断技术领域。通过一个曲轴键相传感器和一个曲轴转速测量传感器分别检测键相脉冲信号和齿轮脉冲信号,所述键相脉冲信号为曲轴每转一圈过程中键相标记触发的脉冲信号,所述齿轮脉冲信号为齿轮旋转过程中每个轮齿触发的脉冲信号,所述键相脉冲信号和所述齿轮脉冲信号同步连续采集,通过对所述键相脉冲信号和所述齿轮脉冲信号进行耦合处理,自动计算得到准确的曲轴停机相位,无需停机后人工打开齿轮盘刻度视孔盖板进行齿轮刻度盘指针位置的检查,使得确定曲轴停机相位更加方便。

Description

活塞发动机曲轴停机相位在线监测方法
技术领域
本申请涉及设备状态监测技术领域,特别涉及一种活塞发动机曲轴停机相位在线监测方法。
背景技术
活塞式发动机是一种广泛应用的动力机械,对其进行在线监测,对提高活塞式发动机的运行安全性、可靠性、及可用性具有非常重要的意义和价值。基于热工参数、振动、声发射、缸内动态压力、瞬时转速、扭振以及在线油液监测在线监测技术技术一直是活塞发动机状态监测领域的研究热点,已形成了较为丰硕的监测技术成果。
活塞发动机曲轴是最为关键的运动部件之一,曲轴停机相位是指活塞发动机停机时,曲轴上某一标记位置相对于某一固定参考位置的角度差,用于表示曲轴停止运动后的具体方位。监测并记录每次停机后的曲轴停机相位,统计分析长期监测结果,有助于分析曲轴旋转运动的动平衡特性。理论上分析,曲轴每次停机相位应为0至360°之间的均匀分布,若实际运行中,若监测到停机相位概率分布严重偏离均匀分布,例如在某一较小区间内概率值远大于其它等宽度区域,则能够反应出曲轴的动平衡不良的问题。
但是活塞发动机相关监测技术中,尚缺乏针对上述曲轴停机相位的监测技术。
发明内容
本申请实施例提供了一种活塞发动机曲轴停机相位在线监测方法,有助于分析曲轴旋转运动的动平衡特性。所述技术方案如下:
获取曲轴键相传感器输出的键相脉冲信号和曲轴转速传感器输出的齿轮脉冲信号,其中所述键相脉冲信号为曲轴每转一圈过程中键相标记触发的脉冲信号,所述齿轮脉冲信号为齿轮旋转过程中每个轮齿触发的脉冲信号,所述键相脉冲信号和所述齿轮脉冲信号同步连续采集;
在所述键相脉冲信号中,寻找曲轴停止运动前的最后一个键相脉冲,并以所述最后一个键相脉冲触发点为参考起始点;在所述齿轮脉冲信号中,寻找曲轴停止运动前的最后一个齿轮脉冲,并以所述最后一个齿轮脉冲截止点为参考结束点;
在所述齿轮脉冲信号中,截取上述参考起始点与上述参考结束点之间的齿轮脉冲信号段,获得第一处理信号;
消除所述第一处理信号中的趋势分量,获得第二处理信号;
检测所述第二处理信号中的零点个数,以及其中第一个零点和最后一个零点的索引位置,得到第一零点位置和最后零点位置;
在所述第二处理信号中,以所述第一零点位置和所述最后零点位置为分割点,将所述第二处理信号分割成三段信号,依次为第一截取信号、第二截取信号和第三截取信号;
计算所述第一截取信号段对应的曲轴旋转角度,获得第一旋转角度;计算所述第二截取信号段对应的曲轴旋转角度,获得第二旋转角度;计算所述第三截取信号段对应的曲轴旋转角度,获得第三旋转角度;
将所述第一旋转角度、第二旋转角度和第三旋转角度进行加和,获得曲轴相对其键相位置的停机相位。
在一种可能的实现方式中,所述寻找曲轴停止运动前的最后一个键相脉冲,并以所述最后一个键相脉冲触发点为参考起始点,包括:求取所述键相脉冲信号的最大值和最小值;选择一个(0,1)之间的第一计算系数,通过公式一确定键相触发门限;以上升沿或下降沿触发的方式,确定所述键相脉冲信号最后一次穿越触发门限的键相脉冲即为最后一个键相脉冲;确定最后一个键相脉冲穿越触发门限时对应的数据点,获得所述参考起始点。
公式一:
Thrkey=α·keymax+(1-α)·keymin,α∈(0,1)
其中,Thrkey为所述键相触发门限值;keymax为所述键相脉冲信号最大值;keymin为所述键相脉冲信号最小值;α为所述第一计算系数,α∈(0,1)。该第一计算系数可在(0,1)之间任意取值,建议取值范围为(1/3,2/3)。
在一种可能的实现方式中,寻找曲轴停止运动前的最后一个齿轮脉冲,并以所述最后一个齿轮脉冲截止点为参考结束点,包括:曲轴停止旋转后,截取一段齿轮脉冲信号,获得第一齿轮脉冲停机信号,并求取所述第一齿轮脉冲停机信号的平均值、最大值和最小值;选择两个(0,1)之间的第二计算系数和第三计算系数,通过公式二确定齿轮脉冲信号停机幅值范围;确定停机过程中所述齿轮脉冲信号超出所述停机幅值范围的最后一个数据点,获得所述参考结束点;
公式二:
Thr1gear=(1+β)·gearmax-β·gearmean,β∈(0,1)
Thr2gear=(1+γ)·gearmin-γ·gearmean,γ∈(0,1)
其中,Thr1gear和Thr2gear分别为所述第一齿轮脉冲停机信号幅值上限和下限;gearmax、gearmin和gearmean分别表示所述第一齿轮脉冲停机信号最大值、最小值和平均值;β和γ分别为所述第二和第三计算系数,β∈(0,1),r∈(0,1)。
在一种可能的实现方式中,消除所述第一处理信号中的趋势分量,获得第二处理信号,包括:计算所述第一处理信号极大值上包络线;计算所述第一处理信号极小值下包络线;通过公式三计算所述第一处理信号趋势分量;将所述第一处理信号减去所述第一处理信号趋势分量,得到所述第二处理信号;
公式三:
Trend=(Lup+Llower)/2
其中,Trend为所述第一处理信号趋势分量;Lup和Llower分别为所述第一处理信号极大值上包络线和极小值下包络线。
在一种可能的实现方式中,计算所述第一截取信号段对应的曲轴旋转角度,获得第一旋转角度,包括:计算所述第一截取信号起始点幅值和所述第二处理信号的峰峰值;判断所述第一截取信号的单调性;当所述第一截取信号单调时,通过公式四计算得到所述第一旋转角度,当所述第一截取信号非单调时,通过公式五计算得到所述第一旋转角度;
公式四:
Figure BDA0002683996840000041
公式五:
Figure BDA0002683996840000042
在公式四和公式五中,Pha1为所述第一旋转角度;y1为所述第一截取信号起始点幅值;PP为所述第二处理信号的峰峰值;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
在一种可能的实现方式中,计算所述第二截取信号段对应的曲轴旋转角度,获得第二旋转角度,包括:检测所述第二处理信号中的零点个数,通过公式六计算得到所述第二旋转角度;
公式六:
Figure BDA0002683996840000043
其中,Pha2为所述第二旋转角度;k为所述第二处理信号中的零点个数;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
在一种可能的实现方式中,计算所述第三截取信号段对应的曲轴旋转角度,获得第三旋转角度,包括:计算所述第三截取信号终止点幅值和所述第二处理信号的峰峰值;判断所述第三截取信号的单调性;当所述第三截取信号单调时,通过公式七计算得到所述第三旋转角度,当所述第三截取信号非单调时,通过公式八计算得到所述第一旋转角度;
公式七:
Figure BDA0002683996840000051
公式八:
Figure BDA0002683996840000052
在公式七和公式八中,Pha3为所述第三旋转角度;y_为所述第三截取信号终止点幅值;PP为所述第二处理信号的峰峰值;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是根据本申请实施例提供的一种活塞发动机曲轴停机相位在线监测方法流程图;
图2是根据本申请实施例提供的一种活塞发动机曲轴停机相位在线监测方法流程示意图;
图3是根据本申请实施例提供的一种同步采集的活塞发动机键相脉冲信号和齿轮脉冲信号示意图;
图4是根据本申请实施例提供的一种曲轴停止旋转后截取的一段齿轮脉冲停机信号示意图;
图5是根据本申请实施例提供的一种活塞发动机齿轮脉冲信号趋势分量计算示意图;
图6是根据本申请实施例提供的一种活塞发动机去除趋势分量后的第二处理信号示意图;
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合附图对本申请实施方式作进一步地详细描述。
图1是本申请提供的一种活塞发动机曲轴停机相位在线监测方法流程图。参见图1,该方法包括:
101、同步连续采集获取曲轴键相传感器输出的键相脉冲信号和曲轴转速传感器输出的齿轮脉冲信号。
102、在所述键相脉冲信号中,寻找曲轴停止运动前的最后一个键相脉冲,并以所述最后一个键相脉冲触发点为参考起始点;在所述齿轮脉冲信号中,寻找曲轴停止运动前的最后一个齿轮脉冲,并以所述最后一个齿轮脉冲截止点为参考结束点。
在一种可能的实现方式中,所述寻找曲轴停止运动前的最后一个键相脉冲,并以所述最后一个键相脉冲触发点为参考起始点,包括:求取所述键相脉冲信号的最大值和最小值;选择一个(0,1)之间的第一计算系数,通过公式一确定键相触发门限;以上升沿或下降沿触发的方式,确定所述键相脉冲信号最后一次穿越触发门限的键相脉冲即为最后一个键相脉冲;确定最后一个键相脉冲穿越触发门限时对应的数据点,获得所述参考起始点。
公式一:
Thrkey=α·keymax+(1-α)·keymin,α∈(0,1)
其中,Thrkey为所述键相触发门限值;keymax为所述键相脉冲信号最大值;keymin为所述键相脉冲信号最小值;α为所述第一计算系数,α∈(0,1)。
在一种可能的实现方式中,寻找曲轴停止运动前的最后一个齿轮脉冲,并以所述最后一个齿轮脉冲截止点为参考结束点,包括:曲轴停止旋转后,截取一段齿轮脉冲信号,获得第一齿轮脉冲停机信号,并求取所述第一齿轮脉冲停机信号的平均值、最大值和最小值;选择两个(0,1)之间的第二计算系数和第三计算系数,通过公式二确定齿轮脉冲信号停机幅值范围;确定停机过程中所述齿轮脉冲信号超出所述停机幅值范围的最后一个采样点对应的数据点,获得所述参考结束点;
公式二:
Thr1gear=(1+β)·gearmax-β·gearmean,β∈(0,1)
Thr2gear=(1+γ)·gearmin-γ·gearmean,γ∈(0,1)
其中,Thr1gear和Thr2gear分别为所述第一齿轮脉冲停机信号幅值上限和下限;gearmax、gearmin和gearmean分别表示所述第一齿轮脉冲停机信号最大值、最小值和平均值;β和γ分别为所述第二和第三计算系数,β∈(0,1),r∈(0,1)。
103、在所述齿轮脉冲信号中,截取上述参考起始点与上述参考结束点之间的齿轮脉冲信号段,获得第一处理信号。
104、消除所述第一处理信号中的趋势分量,获得第二处理信号。
在一种可能的实现方式中,消除所述第一处理信号中的趋势分量,获得第二处理信号,包括:计算所述第一处理信号极大值上包络线;计算所述第一处理信号极小值下包络线;通过公式三计算所述第一处理信号趋势分量;将所述第一处理信号减去所述第一处理信号趋势分量,得到所述第二处理信号;
公式三:
Trend=(Lup+Llower)/2
其中,Trend为所述第一处理信号趋势分量;Lup和Llower分别为所述第一处理信号极大值上包络线和极小值下包络线。
105、检测所述第二处理信号中的零点个数,以及其中第一个零点和最后一个零点的索引位置,得到第一零点位置和最后零点位置;
106、在所述第二处理信号中,以所述第一零点位置和所述最后零点位置为分割点,将所述第二处理信号分割成三段信号,依次为第一截取信号、第二截取信号和第三截取信号;
107、计算所述第一截取信号段对应的曲轴旋转角度,获得第一旋转角度;计算所述第二截取信号段对应的曲轴旋转角度,获得第二旋转角度;计算所述第三截取信号段对应的曲轴旋转角度,获得第三旋转角度。
在一种可能的实现方式中,计算所述第一截取信号段对应的曲轴旋转角度,获得第一旋转角度,包括:计算所述第一截取信号起始点幅值和所述第二处理信号的峰峰值;判断所述第一截取信号的单调性;当所述第一截取信号单调时,通过公式四计算得到所述第一旋转角度,当所述第一截取信号非单调时,通过公式五计算得到所述第一旋转角度;
公式四:
Figure BDA0002683996840000081
公式五:
Figure BDA0002683996840000082
在公式四和公式五中,Pha1为所述第一旋转角度;y1为所述第一截取信号起始点幅值;PP为所述第二处理信号的峰峰值;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
在一种可能的实现方式中,计算所述第二截取信号段对应的曲轴旋转角度,获得第二旋转角度,包括:检测所述第二处理信号中的零点个数,通过公式六计算得到所述第二旋转角度;
公式六:
Figure BDA0002683996840000083
其中,Pha2为所述第二旋转角度;k为所述第二处理信号中的零点个数;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
在一种可能的实现方式中,计算所述第三截取信号段对应的曲轴旋转角度,获得第三旋转角度,包括:计算所述第三截取信号终止点幅值和所述第二处理信号的峰峰值;判断所述第三截取信号的单调性;当所述第三截取信号单调时,通过公式七计算得到所述第三旋转角度,当所述第三截取信号非单调时,通过公式八计算得到所述第一旋转角度;
公式七:
Figure BDA0002683996840000091
公式八:
Figure BDA0002683996840000092
在公式七和公式八中,Pha3为所述第三旋转角度;y_为所述第三截取信号终止点幅值;PP为所述第二处理信号的峰峰值;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
108、将所述第一旋转角度、第二旋转角度和第三旋转角度进行加和,获得曲轴相对其键相位置的停机相位。
下面以计算机设备通过本申请的活塞发动机曲轴停机相位在线监测方法确定一台型号为18PA6B活塞发动机某次停机过程的曲轴停机相位为例进行说明。
步骤(1):同步连续采集获取曲轴键相传感器输出的键相脉冲信号和曲轴转速传感器输出的齿轮脉冲信号。在本申请实施例中,同步采集获取的活塞发动机键相脉冲信号和齿轮脉冲信号如图3所示。
步骤(2):在所述键相脉冲信号中,寻找曲轴停止运动前的最后一个键相脉冲,并以所述最后一个键相脉冲触发点为参考起始点;在所述齿轮脉冲信号中,寻找曲轴停止运动前的最后一个齿轮脉冲,并以所述最后一个齿轮脉冲截止点为参考结束点。在本申请实施例中,求取图3所示键相脉冲信号的最大值keymax=1.8和最小值keymin=0;选择第一计算系数α=0.6,通过公式一确定键相触发门限Thrkey=1.08;可选的,采用上升沿触发方式,确定最后一个键相脉冲穿越触发门限时对应的数据点的索引地址为29670,获得所述参考起始点。
公式一:
Thrkey=α·keymax+(1-α)·keymin=1.08
曲轴停止旋转后,截取一段齿轮脉冲信号,获得第一齿轮脉冲停机信号,如图4所示,求取齿轮脉冲停机信号的平均值gearmean=1.4866、最大值gearmax=1.4972和最小值gearmin=1.4762;选择第二计算系数β=0.5和第三计算系数γ=0.5,通过公式二确定Thr1gear=1.5025和Thr2gear=1.4710,齿轮脉冲信号停机幅值范围为(1.4710,1.5025);确定停机过程中所述齿轮脉冲信号超出所述停机幅值范围的最后一个数据点的索引地址为32599,获得所述参考结束点;
公式二:
Thr1gear=(1+β)·gearmax-β·gearmean=(1+0.5)×1.4972-0.5×1.4866=1.5025
Thr2gear=(1+γ)·gearmin-γ·gearmean=(1+0.5)×1.4762-0.5×1.4866=1.4710
步骤(3):在所述齿轮脉冲信号中,截取上述参考起始点与上述参考结束点之间的齿轮脉冲信号段,获得第一处理信号。
步骤(4):消除所述第一处理信号中的趋势分量,获得第二处理信号。
计算所述第一处理信号极大值上包络线Lup;计算所述第一处理信号极小值下包络线Llower;通过公式三计算所述第一处理信号趋势分量Trend;将所述第一处理信号减去所述第一处理信号趋势分量,得到所述第二处理信号,处理过程示意图如图5所示;
公式三:
Trend=(Lup+Llower)/2
步骤(5):检测所述第二处理信号中的零点个数k=11,以及其中第一个零点的索引位置index1=9和最后一个零点的索引位置index_=2414,得到第一零点位置和最后零点位置,示意图如图5所示;
步骤(6):在所述第二处理信号中,以所述第一零点位置和所述最后零点位置为分割点,将所述第二处理信号分割成三段信号,依次为第一截取信号、第二截取信号和第三截取信号;
步骤(7):计算所述第一截取信号段对应的曲轴旋转角度,获得第一旋转角度;计算所述第二截取信号段对应的曲轴旋转角度,获得第二旋转角度;计算所述第三截取信号段对应的曲轴旋转角度,获得第三旋转角度。在本申请实施例中,首先,计算第一旋转角度,计算所述第一截取信号起始点幅值y1=-0.2048和所述第二处理信号的峰峰值PP=1.7830;确定所述第一截取信号是单调的,通过公式四计算得到所述第一旋转角度Pha1=0.0899°;
公式四:
Figure BDA0002683996840000111
然后,计算第二旋转角度,根据已经检测的所述第二处理信号中零点个数k=11,通过公式六计算得到所述第二旋转角度Pha2=7.8261°;
公式六:
Figure BDA0002683996840000112
最后,计算第三旋转角度,计算所述第三截取信号终止点幅值y_=0.5845和所述第二处理信号的峰峰值PP=1.7830;确定所述第三截取信号是非单调的,通过公式八计算得到所述第一旋转角度Pha3=0.5261°;
公式八:
Figure BDA0002683996840000113
步骤(8):将所述第一旋转角度、第二旋转角度和第三旋转角度进行加和,获得曲轴相对其键相位置的停机相位Pha=Pha1+Pha2+Pha3=0.0899°+7.8261°+0.5261°=8.4421°。
以上所述仅为本申请的可选实施例,并不用以限制本申请,凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (7)

1.一种活塞发动机曲轴停机相位在线监测方法,其特征在于,包括以下步骤:
同步连续采集曲轴键相传感器输出的键相脉冲信号和曲轴转速传感器输出的齿轮脉冲信号,其中所述键相脉冲信号为曲轴每转一圈过程中键相标记触发的脉冲信号,所述齿轮脉冲信号为齿轮旋转过程中每个轮齿触发的脉冲信号;
在所述键相脉冲信号中,寻找曲轴停止运动前的最后一个键相脉冲,并以所述最后一个键相脉冲触发点为参考起始点;在所述齿轮脉冲信号中,寻找曲轴停止运动前的最后一个齿轮脉冲,并以所述最后一个齿轮脉冲截止点为参考结束点;
在所述齿轮脉冲信号中,截取上述参考起始点与上述参考结束点之间的齿轮脉冲信号段,获得第一处理信号;
消除所述第一处理信号中的趋势分量,获得第二处理信号;
检测所述第二处理信号中的零点个数,以及其中第一个零点和最后一个零点的索引位置,得到第一零点位置和最后零点位置;
在所述第二处理信号中,以所述第一零点位置和所述最后零点位置为分割点,将所述第二处理信号分割成三段信号,依次为第一截取信号、第二截取信号和第三截取信号;
计算所述第一截取信号段对应的曲轴旋转角度,获得第一旋转角度;计算所述第二截取信号段对应的曲轴旋转角度,获得第二旋转角度;计算所述第三截取信号段对应的曲轴旋转角度,获得第三旋转角度;
将所述第一旋转角度、第二旋转角度和第三旋转角度进行加和,获得曲轴相对其键相位置的停机相位。
2.根据权利要求1所述的方法,其特征在于,所述寻找曲轴停止运动前的最后一个键相脉冲,并以所述最后一个键相脉冲触发点为参考起始点,包括:
求取所述键相脉冲信号的最大值和最小值;选择一个(0,1)之间的第一计算系数,通过公式一确定键相触发门限;以上升沿或下降沿触发的方式,确定所述键相脉冲信号最后一次穿越触发门限的键相脉冲即为最后一个键相脉冲;确定最后一个键相脉冲穿越触发门限时对应的数据点,获得所述参考起始点;
公式一:
Thrkey=α·keymax+(1-α)·keymin,α∈(0,1)
其中,Thrkey为所述键相触发门限值;keymax为所述键相脉冲信号最大值;keymin为所述键相脉冲信号最小值;α为所述第一计算系数,α∈(0,1)。
3.根据权利要求1所述的方法,其特征在于,所述寻找曲轴停止运动前的最后一个齿轮脉冲,并以所述最后一个齿轮脉冲截止点为参考结束点,包括:
曲轴停止旋转后,截取一段齿轮脉冲信号,获得第一齿轮脉冲停机信号,并求取所述第一齿轮脉冲停机信号的平均值、最大值和最小值;选择两个(0,1)之间的第二计算系数和第三计算系数,通过公式二确定齿轮脉冲信号停机幅值范围;确定停机过程中所述齿轮脉冲信号超出所述停机幅值范围的最后一个数据点,获得所述参考结束点;
公式二:
Thr1gear=(1+β)·gearmax-β·gearmean,β∈(0,1)
Thr2gear=(1+γ)·gearmin-γ·gearmean,γ∈(0,1)
其中,Thr1gear和Thr2gear分别为所述第一齿轮脉冲停机信号幅值上限和下限;gearmax、gearmin和gearmean分别表示所述第一齿轮脉冲停机信号最大值、最小值和平均值;β和γ分别为所述第二和第三计算系数,β∈(0,1),r∈(0,1)。
4.根据权利要求1所述的方法,其特征在于,所述消除所述第一处理信号中的趋势分量,获得第二处理信号,包括:
计算所述第一处理信号极大值上包络线;计算所述第一处理信号极小值下包络线;通过公式三计算所述第一处理信号趋势分量;将所述第一处理信号减去所述第一处理信号趋势分量,得到所述第二处理信号;
公式三:
Trend=(Lup+Llower)/2
其中,Trend为所述第一处理信号趋势分量;Lup和Llower分别为所述第一处理信号极大值上包络线和极小值下包络线。
5.根据权利要求1所述的方法,其特征在于,所述计算所述第一截取信号段对应的曲轴旋转角度,获得第一旋转角度,包括:
计算所述第一截取信号起始点幅值和所述第二处理信号的峰峰值;判断所述第一截取信号的单调性;当所述第一截取信号单调时,通过公式四计算得到所述第一旋转角度,当所述第一截取信号非单调时,通过公式五计算得到所述第一旋转角度;
公式四:
Figure FDA0003653652820000031
公式五:
Figure FDA0003653652820000032
在公式四和公式五中,Pha1为所述第一旋转角度;y1为所述第一截取信号起始点幅值;PP为所述第二处理信号的峰峰值;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
6.根据权利要求1所述的方法,其特征在于,所述计算所述第二截取信号段对应的曲轴旋转角度,获得第二旋转角度,包括:
检测所述第二处理信号中的零点个数,通过公式六计算得到所述第二旋转角度;
公式六:
Figure FDA0003653652820000041
其中,Pha2为所述第二旋转角度;k为所述第二处理信号中的零点个数;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
7.根据权利要求1所述的方法,其特征在于,所述计算所述第三截取信号段对应的曲轴旋转角度,获得第三旋转角度,包括:
计算所述第三截取信号终止点幅值和所述第二处理信号的峰峰值;判断所述第三截取信号的单调性;当所述第三截取信号单调时,通过公式七计算得到所述第三旋转角度,当所述第三截取信号非单调时,通过公式八计算得到所述第三旋转角度;
公式七:
Figure FDA0003653652820000042
公式八:
Figure FDA0003653652820000043
在公式七和公式八中,Pha3为所述第三旋转角度;y_为所述第三截取信号终止点幅值;PP为所述第二处理信号的峰峰值;N表示曲轴飞轮齿盘的齿数,为已知结构参数。
CN202010970888.7A 2020-09-15 2020-09-15 活塞发动机曲轴停机相位在线监测方法 Active CN112284428B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010970888.7A CN112284428B (zh) 2020-09-15 2020-09-15 活塞发动机曲轴停机相位在线监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010970888.7A CN112284428B (zh) 2020-09-15 2020-09-15 活塞发动机曲轴停机相位在线监测方法

Publications (2)

Publication Number Publication Date
CN112284428A CN112284428A (zh) 2021-01-29
CN112284428B true CN112284428B (zh) 2022-07-12

Family

ID=74420459

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010970888.7A Active CN112284428B (zh) 2020-09-15 2020-09-15 活塞发动机曲轴停机相位在线监测方法

Country Status (1)

Country Link
CN (1) CN112284428B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102063087B (zh) * 2010-11-18 2013-08-21 北京第二机床厂有限公司 曲轴工件回转测速及相位准停控制系统
DE102014221293A1 (de) * 2014-10-21 2016-04-21 Volkswagen Aktiengesellschaft Verfahren und Steuervorrichtung zum Erkennen eines Verbrennungsvorgangs einer Verbrennungskraftmaschine eines Hybridfahrzeugs
CN105510036A (zh) * 2015-12-02 2016-04-20 浙江吉利汽车研究院有限公司 一种检测发动机停机相位的传感器及方法
CN108593304A (zh) * 2018-03-23 2018-09-28 广东核电合营有限公司 核电应急柴油机多工况预警状态监测和诊断分析系统及方法
CN109900350A (zh) * 2019-03-14 2019-06-18 中国船舶重工集团公司第七0四研究所 转速与键相融合测试装置及其方法
CN111412978B (zh) * 2020-04-22 2021-06-08 北京化工大学 一种基于无故障振动信号的往复机械异常检测方法

Also Published As

Publication number Publication date
CN112284428A (zh) 2021-01-29

Similar Documents

Publication Publication Date Title
CN103743922B (zh) 一种基于dsp的发动机转速测量方法和装置
CN201615907U (zh) 一种旋转台低角速率检测装置
EP0184536B1 (en) Absolute compression test
CN206905879U (zh) 一种分离式叶片固有频率检测装置
CN102128674B (zh) 测量水轮机过流部件卡门涡振动的方法
CN103308707A (zh) 一种自适应转速测量方法
CN1955710A (zh) 发动机燃油系统故障的检测装置及方法
CN111947957A (zh) 一种振动信号采样方法、装置、存储介质及电子设备
CN112284428B (zh) 活塞发动机曲轴停机相位在线监测方法
CN111239432A (zh) 基于叶端定时的叶片到达时间测量方法及其采集系统
CN113431680A (zh) 一种基于缺齿齿盘轮齿脉冲的发动机监测信号整周期采集和瞬时转速计算方法
CN109344078A (zh) 应用fpga的嵌入式实时操作系统时间性能测试方法
CN109558041B (zh) 基于gpu加速的叶尖间隙信号采集、处理和传输方法
CN1112508C (zh) 一种柴油机喷射系统的故障诊断方法及其装置
CN105157914A (zh) 一种内燃机缸内压力信号时域转角度域的系统及方法
CN205333671U (zh) 一种基于双频计的风速采集器
US4520658A (en) Method of locating engine top dead center position
CN202948288U (zh) 一种基于金属带锯条颤振特性的锯切负载检测装置
CN110470265A (zh) 一种基于fpga和dsp的叶尖间隙与叶尖定时数据同步采集方法
CN109871661B (zh) 一种叶端定时信号丢失的识别补全方法
JPS6363850B2 (zh)
CN114091544A (zh) 基于振动信号与神经网络的tbm盘形滚刀磨损识别系统
CN101886980A (zh) 一种用于发动机实时燃烧分析的设备
CN201615827U (zh) 机械旋转系统的定转速触发非稳态信号采集装置
US4122706A (en) Apparatus for computing the mean indicated pressure of an engine

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