CN107449600A - 一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统 - Google Patents

一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统 Download PDF

Info

Publication number
CN107449600A
CN107449600A CN201710859314.0A CN201710859314A CN107449600A CN 107449600 A CN107449600 A CN 107449600A CN 201710859314 A CN201710859314 A CN 201710859314A CN 107449600 A CN107449600 A CN 107449600A
Authority
CN
China
Prior art keywords
mrow
mfrac
msub
band saw
band
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
Application number
CN201710859314.0A
Other languages
English (en)
Other versions
CN107449600B (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.)
Beihua University
Original Assignee
Beihua University
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 Beihua University filed Critical Beihua University
Priority to CN201710859314.0A priority Critical patent/CN107449600B/zh
Publication of CN107449600A publication Critical patent/CN107449600A/zh
Application granted granted Critical
Publication of CN107449600B publication Critical patent/CN107449600B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING 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
    • G01D21/00—Measuring or testing not otherwise provided for
    • G01D21/02—Measuring two or more variables by means not covered by a single other subclass

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统，包括：获取待检测木工带锯条的主轴转速、进料速度和锯条张紧力、带锯条横向振动位移；根据获取主轴转速、进料速度和锯条张紧力、带锯条横向振动位移，得到基于杜芬方程和改进的QR分解方法的Lyapunov指数；根据Lyapunov指数与Lyapunov指数范围对照表进行比较，判断木工带锯条故障状况。通过监测木工带锯条数，进行Lyapunov指数检测，根据Lyapunov指数，进行木工带锯条损伤判断，并且Lyapunov指数为基于杜芬方程最大Lyapunov指数，能进行微弱信号检测，将其运用到木工带锯条掉齿故障诊断中，实现木工带锯条故障实时在线的预报，避免木工带锯条故障发现不及时而导致安全故障的产生及木材加工质量下降、木材极大浪费。

Description

一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统
技术领域
本发明涉及故障检测技术领域，尤其涉及一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统。
背景技术
木工带锯机是在木材加工企业中广泛使用的机床。它的质量优劣直接影响企业的经济效益。作为旋转机械木工带锯机设备监测诊断已成为林业企业经济发展的迫切需要。从目前木工机械在线故障诊断方法来看，解决其诊断问题的关键之一是如何对监测诊断中获得的机械振动动态信号的非平稳性非线性进行有效地分析。
木工机械设备在运行过程中产生的多发故障，如剥落、摩擦、松动、爬行、冲击、裂纹、断裂、掉齿、喘振等，都将导致其输出的动态信号非平稳性的出现。而且木工设备运行中的驱动力、阻尼力、弹性力的非线性以及系统本身的结构反映在动态信号上也具有非线性。即使木工机械稳态运行，当出现摩擦、冲击等故障时，其转动件的阻尼、刚度、弹性力等会发生变化，呈现出非线性，使输出的振动信号变得非平稳。以上情况表明，从实际工程中获得的动态信号，它们的平稳性是相对的、局部的，而非平稳性非线性是绝对的、广泛的。木工带锯机状态监测与故障诊断的非线性方法具有特殊性，分析系统的动态特性，找出系统的非线性特征，达到对木工设备进行诊断、监测的目的是非常重要的。
现有的木工带锯条产生裂纹和掉齿缺陷只能事后诊断，在加工木材时由于其采用的是闭式切削，极不容易看见带锯条状态，只能根据声音诊断，由于经验问题很难精确、及时；从而造成木工带锯条在锯切中“放炮”-断裂，从而影响带锯机的安全工作，影响操作者的人身安全。以往的研究方法采用木工带锯条横向振动位移和带锯条横向振动频率数值范围进行判断，其存在范围交叉问题，很难做到精确诊断。
发明内容
针对上述缺陷或不足，本发明的目的在于提供一种在线木工带锯条裂纹和掉齿故障检测诊断方法。
本发明的技术方案为：
一种在线木工带锯条裂纹和掉齿故障检测诊断方法,包括：
1)、实时获取待检测木工带锯条的主轴转速、进料速度和带锯条张紧力、带锯条横向振动位移；
2)、根据所获取的主轴转速、进料速度和锯条张紧力、带锯条横向振动位移，得到基于杜芬方程和改进的QR分解方法的Lyapunov指数；
3)、根据Lyapunov指数与Lyapunov指数范围对照表进行比较，判断木工带锯条的故障状况。
所述步骤2为：
2.1定义Lyapunov指数为：
式中F(x)-n维动力系统在带锯条横向振动位移，x＝xi为迭代点处的锯条张紧力；λ-为原动力系统的Lyapunov指数；
2.2、基于杜芬方程的Lyapunov指数计算，受迫Duffing方程形式为：
式中x-所测木工带锯条横向振动位移；
a-阻尼系数矩阵s-1；C-系统库伦阻尼系数矩阵N*sm-1；m-系统质量系数矩阵；b-圆频率矩阵Hz2，k-系统弹簧刚度矩阵；m-系统弹簧刚度矩阵；
c-二维非自治系统木工带锯条横向振动位移方程的常系数；
在系数b矩阵中的b0＝2πn，n-木工带锯轮回转速度，也称木工带锯条横向振动基频；其他为以b0为整数倍的系数；
-受迫振动圆频率，V-木工带锯机进料系统速度；R-进料跑车车轮的直径；
2.3设公式(2-2)解为λ1(t)和λ2(t)，则用改进的QR分解方法得到
式中
求解方程得到得到
式中λ1和λ2是受迫Duffing方程的两个Lyapunov指数；
λ1和λ2之间关系为：
对于杜芬方程的Jacobian矩阵，其对角元素J11≡0，J22是常数，令J22＝c,c为常数,积分式(2-4)得到：
λ1(t)+λ2(t)＝ct+d，c、d为任意常数。
所述步骤1)中，所述锯条张紧力通过数字式力传感器获得；主轴转速通过数字式非接触转速测定仪获取。
所述Lyapunov指数范围对照表为木工带锯条的正常情况下以及故障时的Lyapunov指数范围。
一种在线木工带锯条裂纹和掉齿故障检测诊断系统，包括：监测系统、处理单元、以及判断单元；
所述监测系统，用于实时获取待检测木工带锯条的主轴转速、进料速度和锯条张紧力、木工带锯条横向振动位移；
所述处理单元，用于根据所获取的主轴转速、进料速度和锯条张紧力、木工带锯条横向振动位移，得到基于杜芬方程和改进的QR分解方法的Lyapunov指数；
所述判断单元，用于根据Lyapunov指数与Lyapunov指数范围对照表进行比较，判断木工带锯条的故障状况。
所述处理单元具体用于：
2.1定义Lyapunov指数为：
式中F(x)-n维动力系统在带锯条横向振动位移，x＝xi为迭代点处的锯条张紧力；λ-为原动力系统的Lyapunov指数；
2.2、基于杜芬方程的Lyapunov指数计算，受迫Duffing方程形式为：
式中x-所测木工带锯条横向振动位移；
a-阻尼系数矩阵s-1；C-系统库伦阻尼系数矩阵N*sm-1；m-系统质量系数矩阵；
b-圆频率矩阵Hz2，k-系统弹簧刚度矩阵；m-系统弹簧刚度矩阵；
c-二维非自治系统木工带锯条横向振动位移方程的常系数；
在系数b矩阵中的b0＝2πn，n-木工带锯轮回转速度，也称木工带锯条横向振动基频；其他为以b0为整数倍的系数；
-受迫振动圆频率，V-木工带锯机进料系统速度；R-进料跑车车轮的直径；
2.3设公式(2-2)解为λ1(t)和λ2(t)，则用改进的QR分解方法得到
式中
求解方程得到得到
式中λ1和λ2是受迫Duffing方程的两个Lyapunov指数；
λ1和λ2之间关系为：
对于杜芬方程的Jacobian矩阵，其对角元素J11≡0，J22是常数，令J22＝c,c为常数,积分式(2-4)得到：
λ1(t)+λ2(t)＝ct+d，c、d为任意常数。
所述监测系统包括：木工带锯机踞轮主轴转速监测装置、木工带锯机进料速度监测装置、木工带锯条张紧力监测装置、以及木工带锯条横向振动位移参数监测装置；
所述木工带锯机踞轮主轴转速监测装置通过数字式非接触转速测定仪获取主轴转速；
所述木工带锯机进料速度监测装置通过数字式非速度传感器获取进料速度；
所述木工带锯条张紧力监测装置通过数字式力传感器获取带锯条张紧力；
所述木工带锯条横向振动位移参数监测装置通过电涡流传感器和振动测试数据采集仪获取带锯条横向振动位移。
与现有技术比较，本发明的有益效果为：
本发明提供了一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统，通过实时监测木工带锯条的参数，实时进行Lyapunov指数的检测，根据Lyapunov指数，进行损伤判断，并且该Lyapunov指数为基于杜芬方程和改进的QR分解方法的最大Lyapunov指数，能够进行微弱信号检测，由于在木工带锯条故障幅度较弱时由于故障特征信号淹没在高频振动和噪声中较难分辨，难以检测出信噪比很低的故障特征信号。即功率谱分析或倒频谱分析对高频振动和噪声无免疫力，并且对微弱的故障特征信号不敏感，影响了诊断结果的准确性。在木工带锯条故障监测中，被监测信号伴随着大量噪声，使早期故障特征信号信噪比很低，传统的滤波方法难以实现这种非平稳随机信号的信噪分离。本发明用Duffing方程的非平衡相变对小信号的敏感性及对噪声的免疫力，计算其最大Lyapunov指数，能检测出故障信号，并将其运用到木工带锯条掉齿故障的诊断中，实现故障的实时在线的预报，避免了故障发现不及时而导致判断产生极大误差的现象，提高了故障判断的准确度，并且避免安全故障的产生及木材加工质量的下降、木材的极大浪费。
附图说明
图1是本发明在线木工带锯条裂纹和掉齿故障检测诊断方法流程框图
图2是本发明在线木工带锯条裂纹和掉齿故障检测诊断系统结构图；
图3是本发明在线木工带锯条裂纹和掉齿故障检测诊断监测系统结构图。
具体实施方式
下面将结合附图对本发明做详细描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。
如图1所示，本发明提供了一种在线木工带锯条裂纹和掉齿故障检测诊断方法，包括：
1)、实时获取待检测木工带锯机的主轴转速、进料速度和木工带锯条张紧力、木工带锯条横向振动位移；
本发明通过非接触式电涡流传感器(或激光数字传感器)、数字式速度传感器、数字式力传感器的探头检测，采集木工带锯机的主轴转速、进料速度和木工带锯条张紧力、木工带锯条横向振动位移数据。
2)、根据所获取的主轴转速、进料速度和锯条张紧力、木工带锯条横向振动位移数据，实时得到基于杜芬方程和改进的QR分解方法的Lyapunov指数；
Lyapunov指数不仅仅用来判别混沌现象，还应用到混沌相关的其它领域。Lyapunov指数的值表明了系统混沌的程度，为系统的预测和决策提供了重要信息。所以，Lyapunov指数不仅是判别混沌存在与否的重要指标，也可以用来求取系统从混沌态跃变到周期态的阈值，为系统状态判断提供有效途径。
具体描述用Lyapunov指数确定系统阈值的方法，并基于杜芬方程计算最大Lyapunov指数，进行微弱信号检测。
所谓混沌，最一般意义上的定义是确定性系统的内在随机性，其本质特征是对初值的敏感依赖性。从数学上讲，对于各模型参数都已知的确定性系统，由确定的初始值可以推知该系统的长期行为，甚至可以追溯其以往的状态。然而大量的实例表明，有很多系统对其初值有着很强的依赖性，初值上的“差之毫厘”往往会导致系统长期行为的“失之千里”。因而在实际误差的影响下，这种系统将表现出极大的随机性。然而，这种随机性与由于系统本身具有随机项或随机系数而产生的随机现象是完全不同的。对于一个真正的随机系统，从某一特定时刻的量无法知道以后任何时候量的确定值，即系统行为在短期内也是不可预测的。而对于确定性系统，它的短期行为是完全确定的，但由于对初值依赖的敏感性，使得确切运动在长期内不可预测。这种由确定性系统的内在故有的随机性所引起的随机现象，我们称之为混沌现象。
具体地，所述步骤2)为：
2.1定义Lyapunov指数为：
式中F(x)-n维动力系统在带锯条横向振动位移，x＝xi为迭代点处的锯条张紧力；λ-为原动力系统的Lyapunov指数；
2.2、基于杜芬方程的Lyapunov指数计算，受迫Duffing方程形式为：
式中x-所测木工带锯条横向振动位移；
a-阻尼系数矩阵s-1；C-系统库伦阻尼系数矩阵N*sm-1；m-系统质量系数矩阵；
b-圆频率矩阵Hz2，k-系统弹簧刚度矩阵；m-系统弹簧刚度矩阵；
c-二维非自治系统木工带锯条横向振动位移方程的常系数；
在系数b矩阵中的b0＝2πn，n-木工带锯轮回转速度，也称木工带锯条横向振动基频；其他为以b0为整数倍的系数；
-受迫振动圆频率，V-木工带锯机进料系统速度(所测切削木材时进料跑车的速度m*s-1)；R-进料跑车车轮的直径；
公式(2-2)这是一个二维非自治系统，记为：g(t)＝F cos(ωt)；
2.3设公式(2-2)解为λ1(t)和λ2(t)，则用改进的QR分解方法得到
式中
求解方程得到得到
式中λ1和λ2是受迫Duffing方程的两个Lyapunov指数；
λ1和λ2之间关系为：
对于杜芬方程的Jacobian矩阵，其对角元素J11≡0，J22是常数，令J22＝c,c为常数,积分式(2-4)得到：
λ1(t)+λ2(t)＝ct+d，c、d为任意常数。
3)、根据Lyapunov指数与Lyapunov指数范围对照表进行比较，判断木工带锯条的故障状况。
本发明还提供了一种在线木工带锯条裂纹和掉齿故障检测诊断系统，如图2所示，包括：监测系统1、处理单元2、以及判断单元3；
所述监测系统1，用于实时获取待检测木工带锯条的主轴转速、进料速度和锯条张紧力；
如图3所示，本发明中所述监测系统1包括：木工带锯机踞轮主轴转速监测装置11、木工带锯机进料速度监测装置12、木工带锯条张紧力监测装置13、以及木工带锯条横向振动位移参数监测装置14；
所述木工带锯机踞轮主轴转速监测装置11通过数字式非接触转速测定仪获取主轴转速；
所述木工带锯机进料速度监测装置12通过数字式非接触速度传感器获取进料速度；
所述木工带锯条张紧力监测装置13通过数字式力传感器获取带锯条张紧力
所述木工带锯条横向振动位移参数监测装置14通过电涡流传感器(或数字激光位移传感器)和振动测试数据采集仪获取木工带锯条横向振动位移。
监测系统1将所获取的数据通过计算机、大数据振动数据分析软件及计算的处理单元进行分析，最后通过打印机进行打印。
所述处理单元2，用于根据所获取的主轴转速、进料速度和锯条张紧力，得到基于杜芬方程的Lyapunov指数；所述处理单元(2)具体用于：
2.1定义Lyapunov指数为：
式中F(x)-n维动力系统在带锯条横向振动位移，x＝xi为迭代点处的锯条张紧力；λ-为原动力系统的Lyapunov指数；
2.2、基于杜芬方程的Lyapunov指数计算，受迫Duffing方程形式为：
式中x-所测木工带锯条横向振动位移；
a-阻尼系数矩阵s-1；C-系统库伦阻尼系数矩阵N*sm-1；m-系统质量系数矩阵；
b-圆频率矩阵Hz2，k-系统弹簧刚度矩阵；m-系统弹簧刚度矩阵；
c-二维非自治系统木工带锯条横向振动位移方程的常系数；
在系数b矩阵中的b0＝2πn，n-木工带锯轮回转速度，也称木工带锯条横向振动基频；其他为以b0为整数倍的系数；
-受迫振动圆频率，V-木工带锯机进料系统速度；R-进料跑车车轮的直径；
2.3设公式(2-2)解为λ1(t)和λ2(t)，则用改进的QR分解方法得到
式中
求解方程得到得到
式中λ1和λ2是受迫Duffing方程的两个Lyapunov指数；
λ1和λ2之间关系为：
对于杜芬方程的Jacobian矩阵，其对角元素J11≡0，J22是常数，令J22＝c,c为常数,积分式(2-4)得到：
λ1(t)+λ2(t)＝ct+d，c、d为任意常数。
所述判断单元3，用于根据Lyapunov指数与Lyapunov指数范围对照表进行比较，判断木工带锯条的故障状况。
所述监测系统1包括：木工带锯机踞轮主轴转速监测装置11、木工带锯机进料速度监测装置12、木工带锯条张紧力监测装置13、以及木工带锯条横向振动位移参数监测装置14；
所述木工带锯机踞轮主轴转速监测装置11通过数字式非接触转速测定仪获取主轴转速；
所述木工带锯机进料速度监测装置12通过数字式非速度传感器获取进料速度；
所述木工带锯条张紧力监测装置13通过数字式力传感器获取带锯条张紧力；
所述木工带锯条横向振动位移参数监测装置14通过电涡流传感器(或数字激光位移传感器)和振动测试数据采集仪获取带锯条横向振动位移。
对于本领域技术人员而言，显然能了解到上述具体事实例只是本发明的优选方案，因此本领域的技术人员对本发明中的某些部分所可能作出的改进、变动，体现的仍是本发明的原理，实现的仍是本发明的目的，均属于本发明所保护的范围。

Claims (7)

1.一种在线木工带锯条裂纹和掉齿故障检测诊断方法，其特征在于，包括：
1)、实时获取待检测木工带锯条的主轴转速、进料速度和带锯条张紧力、带锯条横向振动位移；
2)、根据所获取的主轴转速、进料速度和锯条张紧力、带锯条横向振动位移，得到基于杜芬方程和改进的QR分解方法的Lyapunov指数；
3)、根据Lyapunov指数与Lyapunov指数范围对照表进行比较，判断木工带锯条的故障状况。
2.根据权利要求1所述的在线木工带锯条裂纹和掉齿故障检测诊断方法，其特征在于，所述步骤2为：
2.1定义Lyapunov指数为：
<mrow> <mi>&amp;lambda;</mi> <mo>=</mo> <munder> <mrow> <mi>l</mi> <mi>i</mi> <mi>m</mi> </mrow> <mrow> <mi>n</mi> <mo>&amp;RightArrow;</mo> <mi>&amp;infin;</mi> </mrow> </munder> <mfrac> <mn>1</mn> <mi>n</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>l</mi> <mi>n</mi> <msub> <mrow> <mo>|</mo> <mfrac> <mrow> <mi>d</mi> <mi>F</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>x</mi> </mrow> </mfrac> <mo>|</mo> </mrow> <mrow> <mi>x</mi> <mo>=</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式中F(x)-n维动力系统在带锯条横向振动位移，x＝xi为迭代点处的锯条张紧力；λ-为原动力系统的Lyapunov指数；
2.2、基于杜芬方程的Lyapunov指数计算，受迫Duffing方程形式为：
<mrow> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>a</mi> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>b</mi> <mi>x</mi> <mo>+</mo> <msup> <mi>cx</mi> <mn>3</mn> </msup> <mo>=</mo> <mi>F</mi> <mi> </mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;omega;</mi> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式中x-所测木工带锯条横向振动位移；
a-阻尼系数矩阵s-1；C-系统库伦阻尼系数矩阵N*sm-1；m-系统质量系数矩阵；
b-圆频率矩阵Hz2，k-系统弹簧刚度矩阵；m-系统弹簧刚度矩阵；
c-二维非自治系统木工带锯条横向振动位移方程的常系数；
在系数b矩阵中的b0＝2πn，n-木工带锯轮回转速度，也称木工带锯条横向振动基频；其它为以b0为整数倍的系数；
-受迫振动圆频率，V-木工带锯机进料系统速度；R-进料跑车车轮的直径；
2.3设公式(2-2)解为λ1(t)和λ2(t)，则用改进的QR分解方法得到
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>v</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>v</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mn>2</mn> <mi>&amp;theta;</mi> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mn>2</mn> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <mi>d</mi> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mi>sin</mi> <mn>2</mn> <mi>&amp;theta;</mi> <mo>+</mo> <msup> <mi>sin</mi> <mn>2</mn> </msup> <mi>&amp;theta;</mi> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> 1
式中y1＝xy3＝t
求解方程得到得到
式中λ1和λ2是受迫Duffing方程的两个Lyapunov指数；
λ1和λ2之间关系为：
<mrow> <mfrac> <mrow> <msub> <mi>d&amp;lambda;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>d&amp;lambda;</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>J</mi> <mn>11</mn> </msub> <mo>+</mo> <msub> <mi>J</mi> <mn>22</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
对于杜芬方程的Jacobian矩阵，其对角元素J11≡0，J22是常数，令J22＝c,c为常数,积分式(2-4)得到：
λ1(t)+λ2(t)＝ct+d，c、d为任意常数。
3.根据权利要求1所述的在线木工带锯条裂纹和掉齿故障检测诊断方法，其特征在于，所述步骤1)中，所述木工带锯条张紧力通过数字式力传感器获得；主轴转速通过数字式非接触转速测定仪获取。
4.根据权利要求1所述的在线木工带锯条裂纹和掉齿故障检测诊断方法，其特征在于，所述Lyapunov指数范围对照表为木工带锯条的正常情况下以及故障时的Lyapunov指数范围。
5.一种在线木工带锯条裂纹和掉齿故障检测诊断系统，其特征在于，包括：监测系统(1)、处理单元(2)、以及判断单元(3)；
所述监测系统(1)，用于实时获取待检测木工带锯条的主轴转速、进料速度和锯条张紧力、木工带锯条横向振动位移；
所述处理单元(2)，用于根据所获取的主轴转速、进料速度和锯条张紧力、木工带锯条横向振动位移，得到基于杜芬方程和改进的QR分解方法的Lyapunov指数；
所述判断单元(3)，用于根据Lyapunov指数与Lyapunov指数范围对照表进行比较，判断木工带锯条的故障状况。
6.根据权利要求5所述的在线木工带锯条裂纹和掉齿故障检测诊断系统，其特征在于，所述处理单元(2)具体用于：
2.1定义Lyapunov指数为：
<mrow> <mi>&amp;lambda;</mi> <mo>=</mo> <munder> <mrow> <mi>l</mi> <mi>i</mi> <mi>m</mi> </mrow> <mrow> <mi>n</mi> <mo>&amp;RightArrow;</mo> <mi>&amp;infin;</mi> </mrow> </munder> <mfrac> <mn>1</mn> <mi>n</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>l</mi> <mi>n</mi> <msub> <mrow> <mo>|</mo> <mfrac> <mrow> <mi>d</mi> <mi>F</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>x</mi> </mrow> </mfrac> <mo>|</mo> </mrow> <mrow> <mi>x</mi> <mo>=</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式中F(x)-n维动力系统在带锯条横向振动位移，x＝xi为迭代点处的锯条张紧力；λ-为原动力系统的Lyapunov指数；
2.2、基于杜芬方程的Lyapunov指数计算，受迫Duffing方程形式为：
<mrow> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>a</mi> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>b</mi> <mi>x</mi> <mo>+</mo> <msup> <mi>cx</mi> <mn>3</mn> </msup> <mo>=</mo> <mi>F</mi> <mi> </mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;omega;</mi> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> 2
式中x-所测木工带锯条横向振动位移；
a-阻尼系数矩阵s-1；C-系统库伦阻尼系数矩阵N*sm-1；m-系统质量系数矩阵；
b-圆频率矩阵Hz，k-系统弹簧刚度矩阵；m-系统弹簧刚度矩阵；
c-二维非自治系统木工带锯条横向振动位移方程的常系数；
在系数b矩阵中的，b0＝2πn n-木工带锯轮回转速度，也称木工带锯条横向振动基频；其他为以b0为整数倍的系数；
-受迫振动圆频率，V-木工带锯机进料系统速度；R-进料跑车车轮的直径；
2.3设公式(2-2)解为λ1(t)和λ2(t)，则用改进的QR分解方法得到:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>v</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>v</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mn>2</mn> <mi>&amp;theta;</mi> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mn>2</mn> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <mi>d</mi> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mi>sin</mi> <mn>2</mn> <mi>&amp;theta;</mi> <mo>+</mo> <msup> <mi>sin</mi> <mn>2</mn> </msup> <mi>&amp;theta;</mi> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced>
式中y1＝xy3＝t
求解方程得到得到
式中λ1和λ2是受迫Duffing方程的两个Lyapunov指数；
λ1和λ2之间关系为：
<mrow> <mfrac> <mrow> <msub> <mi>d&amp;lambda;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>d&amp;lambda;</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>J</mi> <mn>11</mn> </msub> <mo>+</mo> <msub> <mi>J</mi> <mn>22</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
对于杜芬方程的Jacobian矩阵，其对角元素J11≡0，J22是常数，令J22＝c,c为常数,积分式(2-4)得到：
λ1(t)+λ2(t)＝ct+d，c、d为任意常数。
7.根据权利要求5所述的在线木工带锯条裂纹和掉齿故障检测诊断系统，其特征在于，所述监测系统(1)包括：木工带锯机踞轮主轴转速监测装置(11)、木工带锯机进料速度监测装置(12)、木工带锯条张紧力监测装置(13)、以及木工带锯条横向振动位移参数监测装置(14)；
所述木工带锯机踞轮主轴转速监测装置(11)通过数字式非接触转速测定仪获取主轴转速；
所述带锯机进料速度监测装置(12)通过数字式非接触速度传感器获取进料速度；
所述木工带锯条张紧力监测装置(13)通过数字式力传感器获取带锯条张紧力；
所述木工带锯条横向振动位移参数监测装置(14)通过电涡流传感器和振动测试数据采集仪获取带锯条横向振动位移。
CN201710859314.0A 2017-09-21 2017-09-21 一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统 Expired - Fee Related CN107449600B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710859314.0A CN107449600B (zh) 2017-09-21 2017-09-21 一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710859314.0A CN107449600B (zh) 2017-09-21 2017-09-21 一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统

Publications (2)

Publication Number Publication Date
CN107449600A true CN107449600A (zh) 2017-12-08
CN107449600B CN107449600B (zh) 2019-10-01

Family

ID=60497934

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710859314.0A Expired - Fee Related CN107449600B (zh) 2017-09-21 2017-09-21 一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统

Country Status (1)

Country Link
CN (1) CN107449600B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021170429A1 (de) 2020-02-27 2021-09-02 Gebrüder Linck Maschinenfabrik "Gatterlinck" GmbH & Co. KG Vorrichtung zur vorhersage eines blattabrisses eines bandsägeblattes einer bandsäge

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2557273Y (zh) * 2002-07-03 2003-06-25 陈柏梁 卧式带锯床用工件检测装置
CN2932345Y (zh) * 2006-05-25 2007-08-08 宝山钢铁股份有限公司 一种带锯条张紧力、垂直度测量装置
CN102156873A (zh) * 2010-12-31 2011-08-17 北京航空航天大学 一种基于混沌的机械零部件早期单点故障检测与分类方法
CN203502787U (zh) * 2013-10-15 2014-03-26 丽水学院 一种金属带锯床远程无线故障诊断装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2557273Y (zh) * 2002-07-03 2003-06-25 陈柏梁 卧式带锯床用工件检测装置
CN2932345Y (zh) * 2006-05-25 2007-08-08 宝山钢铁股份有限公司 一种带锯条张紧力、垂直度测量装置
CN102156873A (zh) * 2010-12-31 2011-08-17 北京航空航天大学 一种基于混沌的机械零部件早期单点故障检测与分类方法
CN102156873B (zh) * 2010-12-31 2013-01-30 北京航空航天大学 一种基于混沌的机械零部件早期单点故障检测与分类方法
CN203502787U (zh) * 2013-10-15 2014-03-26 丽水学院 一种金属带锯床远程无线故障诊断装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张焰明: "基于Lyapunov指数分析的液压泵故障诊断方法及实验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021170429A1 (de) 2020-02-27 2021-09-02 Gebrüder Linck Maschinenfabrik "Gatterlinck" GmbH & Co. KG Vorrichtung zur vorhersage eines blattabrisses eines bandsägeblattes einer bandsäge
DE102020105223A1 (de) 2020-02-27 2021-09-02 Gebrüder Linck, Maschinenfabrik "Gatterlinck" GmbH & Co. KG Vorrichtung zur Vorhersage eines Blattabrisses eines Bandsägeblattes einer Bandsäge
DE102020105223B4 (de) 2020-02-27 2021-11-04 Gebrüder Linck, Maschinenfabrik "Gatterlinck" GmbH & Co. KG Vorrichtung zur Vorhersage eines Blattabrisses eines Bandsägeblattes einer Bandsäge

Also Published As

Publication number Publication date
CN107449600B (zh) 2019-10-01

Similar Documents

Publication Publication Date Title
CN103575523B (zh) 基于FastICA-谱峭度-包络谱分析的旋转机械故障诊断方法
CN104764869B (zh) 一种基于多维特征量的变压器气体故障诊断和报警方法
TWI449883B (zh) 結構體安全性之分析方法
CN102175768A (zh) 一种基于振动信号的高铁钢轨伤损探测方法及装置
Dube et al. Vibration based condition assessment of rolling element bearings with localized defects
CN112781820B (zh) 一种滚刀性能退化趋势评估方法
CN104458173B (zh) 钢框架结构突变损伤的识别方法及系统
CN105137277A (zh) 一种变压器突发短路时绕组机械状态的诊断方法
Jena et al. Radial ball bearing inner race defect width measurement using analytical wavelet transform of acoustic and vibration signal
US6873918B2 (en) Control embedded machine condition monitor
He et al. Bearing condition evaluation based on the shock pulse method and principal resonance analysis
De Almeida et al. New technique for evaluation of global vibration levels in rolling bearings
CN107449600B (zh) 一种在线木工带锯条裂纹和掉齿故障检测诊断方法及系统
Lybeck et al. Validating prognostic algorithms: a case study using comprehensive bearing fault data
Zamorano et al. Analysis in the time-frequency domain of different depths of a crack located in a change of section of a shaft
Aye Statistical approach for tapered bearing fault detection using different methods
CN112857798B (zh) 一种基于频谱分析的多轴系机械传动系统质量评估方法及装置
Klinchaeam et al. Fault detection of a spur gear using vibration signal with multivariable statistical parameters.
Ogundare et al. Review of fault detection techniques for health monitoring of helicopter gearbox
CN103076174A (zh) 一种基于锁相放大算法的轴承故障及故障程度诊断方法
Sanderson et al. The applicability of online monitoring for the assurance of fixed jacket structures’ structural integrity
CN115791040A (zh) 一种基于多指标融合的轧辊磨床颤振故障诊断方法
Vispute et al. Detection of Gear Tooth Defects by Using MATLAB & FEA Technique-Review
Singh et al. Time domain analysis in condition monitoring-a
Munde et al. A review on bearing fault detection by vibration signature analysis

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20191001

Termination date: 20210921