CN107153728B - 一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法 - Google Patents
一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法 Download PDFInfo
- Publication number
- CN107153728B CN107153728B CN201710284246.XA CN201710284246A CN107153728B CN 107153728 B CN107153728 B CN 107153728B CN 201710284246 A CN201710284246 A CN 201710284246A CN 107153728 B CN107153728 B CN 107153728B
- Authority
- CN
- China
- Prior art keywords
- bispectrum
- grinding
- acoustic emission
- double
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
Abstract
本发明公开了一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法。该方法基于磨削加工声发射信号,通过计算得到其双谱,对双谱幅值空间进行划分,计算统计信号双谱落入每个幅值子空间的概率,根据子空间分布概率定义信号双谱幅值分布熵,最后通过双谱幅值分布熵有效提取不同砂轮磨损状态下的磨削声发射信号双谱特征。本发明提取的信号特征计算简易,可以定量描述不同砂轮磨损状态下的声发射信号双谱特点,为进一步确定砂轮磨损状态提供了评价指标。
Description
技术领域
本发明属于刀具状态监测技术领域,具体涉及一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法。
背景技术
磨削加工技术处于先进制造领域,是现代制造业中实现精密和超精密加工最有效、应用最广的制造技术。在磨削加工中,砂轮钝化则是影响加工效率与磨削质量的一个非常重要的因素。传统上,主要依靠经验,采用定时修整砂轮的方法来避免砂轮钝化的不利影响。这种方法有非常明显的局限性,并且严重阻碍了磨削加工装备向自动化、智能化方向发展。因此,开展砂轮磨损状态监测方法的研究对于提高磨削加工装备的自动化水平具有非常重要的意义。
金属磨削过程中,砂轮与工件的刮擦、磨粒崩碎、粘结剂破裂等都会产生声发射现象,因此声发射信号蕴含有丰富的磨削加工信息,被广泛应用于砂轮磨损状态监测领域。而且,随着砂轮的磨损,磨削过程从以刮擦、切削作用为主转变为以撞击、挤压作用为主,磨削声发射信号将表现出越来越强的非线性和非高斯性特征。而双谱是分析非线性、非高斯信号的有力工具,近年来在机械状态监测及故障诊断中受到了广泛的关注,并逐步应用于刀具状态监测领域。但是,由于计算得到的信号双谱往往是一个高维矩阵,数据量过于庞大,不适于后续的状态识别;另一方面,传统的双谱特征提取方法存在提取信息不全面,准确率低的问题,难以有效提取磨削声发射信号双谱中隐含的砂轮磨损状态信息。因此,迫切需要寻求一种更加有效的磨削声发射信号双谱特征,来反映不同砂轮状态下,信号双谱之间的特征差异。
发明内容
本发明的目的是针对以上问题,提供了一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法,该方法能够定量描述信号双谱幅值的分布特性,从而有效地提取不同状态磨削声发射信号双谱的差异特征,为进一步的磨损状态评价提供依据。
为达到上述目的,本发明采用如下的技术方案来实现:
一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法,包括以下步骤:
1)获取磨削加工过程中的声发射信号,并计算其双谱;
2)对磨削声发射信号的双谱进行幅值区间划分,确定需要划分的子空间数目以及各区间幅值范围;
3)计算信号双谱落入每个幅值子空间的概率,利用所得到的概率计算磨削声发射信号的双谱幅值分布熵特征。
本发明进一步的改进在于,步骤1)中,双谱的计算步骤为:
101)将采集到的长度为L的磨削声发射数据序列{x(n),n=1,2,…,L}分成M个小段,每段包含N个数据点,分段过程中,使相邻的段与段之间有50%的数据重叠率;
102)去除每段数据均值,并为便于FFT计算,对数据进行必要的补零;
103)依次完成每段数据的FFT计算,对于第i段数据{x(i)(p),p=1,2,…,N},得到:
104)根据各段数据FFT的计算结果,分别求得其分段双谱估计值:
105)对已求得的各段数据双谱估计进行统计平均,得声发射信号的双谱估计值:
本发明进一步的改进在于,步骤2)中,对磨削声发射信号双谱进行幅值区间划分的步骤为:
202)按以下公式确定划分的子空间数目m:
本发明进一步的改进在于,步骤3)中,双谱幅值分布熵的计算步骤为:
式中:count(·)为计数函数;n_sum为双谱矩阵总点数;Ak为所划分的幅值子空间;
通过双谱幅值分布熵来反映不同的砂轮磨损状态。
本发明具有以下的优点:
本发明提供的基于双谱幅值分布熵的砂轮磨损状态特征提取方法,采用双谱分析方法对磨削加工声发射信号进行处理,成功提取了砂轮磨损过程中表现出的非高斯性特征,并利用双谱幅值分布熵解决了特征信息的定量描述问题,有利于实现砂轮磨损状态的直观判别。
进一步,双谱计算过程中借助于快速傅里叶变换算法,减少了计算量,提高了计算效率。
进一步,双谱标准差能够反映双谱的波动情况,而利用双谱标准差对双谱幅值空间进行划分,则可以很好的匹配信号双谱的波动特征,合理反映双谱分布特性。
进一步,结合信息熵能够评价分布不均匀性的特点,利用双谱幅值分布熵这一量化指标有效提取不同砂轮磨损状态下磨削声发射信号双谱的特征差异,为后续进行状态识别提供可靠评价依据。
附图说明
图1为本发明的基本流程图;
图2为齿轮轴1不同磨齿阶段的声发射信号双谱等高线图;其中,图2(a)为磨削第2齿的声发射信号双谱等高线图,图2(b)为磨削第8齿的声发射信号双谱等高线图,图2(c)为磨削第14齿的声发射信号双谱等高线图;
图3为齿轮轴2不同磨齿阶段的声发射信号双谱等高线图;其中,图3(a)为磨削第2齿的声发射信号双谱等高线图,图3(b)为磨削第8齿的声发射信号双谱等高线图,图3(c)为磨削第14齿的声发射信号双谱等高线图;
图4为齿轮轴3不同磨齿阶段的声发射信号双谱等高线图;其中,图4(a)为磨削第2齿的声发射信号双谱等高线图,图4(b)为磨削第8齿的声发射信号双谱等高线图,图4(c)为磨削第14齿的声发射信号双谱等高线图;
图5为实验磨齿声发射信号双谱幅值分布熵均值变化曲线。
具体实施方式
本发明主要实现在磨削加工声发射信号双谱中,提取能够反映砂轮磨损状态变化的有效特征,图1为本发明的基本流程图。对于一个典型的磨齿加工过程,实验中,经过修整后的砂轮磨削加工到第14齿后将会发生钝化,失去有效加工能力,因此可以认为砂轮从磨削第1齿到第14齿经历了一个完整的砂轮工作寿命周期。采集加工每一齿时的磨削声发射信号,共得到43组砂轮完整工作寿命周期的磨削声发射数据,对其进行处理,并提取双谱幅值分布熵特征,具体实施步骤如下:
1)将采集到的长度为L的加工每一齿的磨削声发射数据序列{x(n),n=1,2,…,L}分成M个小段,每段包含N个数据点,分段过程中,使相邻的段与段之间有50%的数据重叠率;实验数据处理中选取L=102400,M=100,N=1024。
2)去除每段数据均值,并为便于FFT计算,对数据进行必要的补零;
3)依次完成每段数据的FFT计算,对于第i段数据{x(i)(p),p=1,2,…,N},得到:
4)根据各段数据FFT的计算结果,分别求得其分段双谱估计值:
其中:ω1,ω2满足关系0≤ω1≤ω2,ω1+ω2≤0,1,…,N/2;*表示取共轭;
5)然后,对已求得的各段数据双谱估计进行统计平均,得到声发射信号的双谱估计值:
其中:M为分段总数;
通过计算得到的不同砂轮磨损状态下的磨削声发射信号双谱如图2~4所示(限于篇幅,仅列举了3个工件典型加工阶段的信号双谱)。
从图中可以清晰的看出,当砂轮处于不同的磨损状态时(即磨削同一工件的不同齿序),磨削声发射信号表现出明显的变化;而加工不同的工件,在砂轮磨损状态相近时(即磨削不同工件的相同齿序),声发射信号双谱则表现出了较强的相似性。因此,磨削声发射信号双谱能够很好的反映砂轮的磨损情况。并且,从图中可以看出,随着砂轮磨损状态的改变,声发射双谱表现出的是幅值分布的变化。为了表征这种变化,下面进一步对磨削声发射双谱进行处理。
其中:max(·),min(·)分别为最大和最小值运算;
7)接着,按以下公式确定所划分子空间数目m:
计算得到的实验声发射信号双谱子空间划分数目m=50。
8)根据得到的子空间数目m,确定每个子空间的幅值范围为
式中:count(·)为计数函数;n_sum为双谱矩阵总点数;Ak为所划分的幅值子空间;
通过双谱幅值分布熵来反映不同的砂轮磨损状态。
图3为计算得到的43组实验声发射数据双谱幅值分布熵均值变化曲线。图中可以清晰的看出,随着加工的不断进行,双谱幅值分布熵呈现出明显上升的趋势。这说明由于砂轮的不断磨损,磨削声发射信号双谱中,落入高值空间的点数不断增多,幅值分布逐渐向高值扩展,信号非高斯特征越来越强,与实际情况相符。因此,双谱幅值分布熵能够有效量化表示磨削声发射信号双谱幅值分布的变化,进而反映出不同的砂轮磨损情况,可以作为表征砂轮磨损状态的有效特征指标。
以上对本发明的典型实例及原理进行了详细说明,对本领域的技术人员而言,依据本发明提供的思想,在具体的实施方式上会有所改变,但这些改变也应视为本发明的保护范围。
Claims (3)
1.一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法,其特征在于,包括以下步骤:
1)获取磨削加工过程中的声发射信号,并计算其双谱;双谱的计算步骤为:
101)将采集到的长度为L的磨削声发射数据序列{x(n),n=1,2,…,L}分成M个小段,每段包含N个数据点,分段过程中,使相邻的段与段之间有50%的数据重叠率;
102)去除每段数据均值,并为便于FFT计算,对数据进行必要的补零;
103)依次完成每段数据的FFT计算,对于第i段数据{x(i)(p),p=1,2,…,N},得到:
104)根据各段数据FFT的计算结果,分别求得其分段双谱估计值:
105)对已求得的各段数据双谱估计进行统计平均,得声发射信号的双谱估计值:
2)对磨削声发射信号的双谱进行幅值区间划分,确定需要划分的子空间数目以及各区间幅值范围;
3)计算信号双谱落入每个幅值子空间的概率,利用所得到的概率计算磨削声发射信号的双谱幅值分布熵特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710284246.XA CN107153728B (zh) | 2017-04-26 | 2017-04-26 | 一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710284246.XA CN107153728B (zh) | 2017-04-26 | 2017-04-26 | 一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107153728A CN107153728A (zh) | 2017-09-12 |
CN107153728B true CN107153728B (zh) | 2020-08-28 |
Family
ID=59792616
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710284246.XA Active CN107153728B (zh) | 2017-04-26 | 2017-04-26 | 一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107153728B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107979842B (zh) * | 2017-11-28 | 2019-11-19 | 电子科技大学 | 一种基于双谱特征和深度学习的手机个体识别方法 |
CN108460200B (zh) * | 2018-02-12 | 2019-07-16 | 青岛理工大学 | 基于性态双谱的多模态性能化抗震设计方法 |
CN113340997B (zh) * | 2021-05-11 | 2024-01-12 | 西安交通大学 | 基于声发射双通道极差的激光冲击强化缺陷在线检测方法 |
CN114280160A (zh) * | 2021-12-15 | 2022-04-05 | 上海机床厂有限公司 | 基于Shannon熵与声发射信号的CBN砂轮性能监测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105820A (zh) * | 2012-05-22 | 2013-05-15 | 华中科技大学 | 一种数控机床切削加工的刀具磨损状态识别方法 |
CN103646239A (zh) * | 2013-12-25 | 2014-03-19 | 武汉大学 | 基于极坐标傅里叶变换的旋转不变性图像特征提取方法 |
CN104723171A (zh) * | 2015-03-17 | 2015-06-24 | 洛阳理工学院 | 一种基于电流与声发射复合信号的刀具磨损监测方法 |
CN105158098A (zh) * | 2015-08-26 | 2015-12-16 | 南京航空航天大学 | 一种磨粒位姿可控的单颗磨粒磨削试验平台及其试验方法 |
CN105215852A (zh) * | 2015-05-29 | 2016-01-06 | 青岛理工大学 | 一种砂轮磨损及g比率的测量装置与方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170066104A9 (en) * | 2009-12-08 | 2017-03-09 | Allison Transmission Inc. | Method for Detecting And/Or Preventing Grind Burn |
-
2017
- 2017-04-26 CN CN201710284246.XA patent/CN107153728B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105820A (zh) * | 2012-05-22 | 2013-05-15 | 华中科技大学 | 一种数控机床切削加工的刀具磨损状态识别方法 |
CN103646239A (zh) * | 2013-12-25 | 2014-03-19 | 武汉大学 | 基于极坐标傅里叶变换的旋转不变性图像特征提取方法 |
CN104723171A (zh) * | 2015-03-17 | 2015-06-24 | 洛阳理工学院 | 一种基于电流与声发射复合信号的刀具磨损监测方法 |
CN105215852A (zh) * | 2015-05-29 | 2016-01-06 | 青岛理工大学 | 一种砂轮磨损及g比率的测量装置与方法 |
CN105158098A (zh) * | 2015-08-26 | 2015-12-16 | 南京航空航天大学 | 一种磨粒位姿可控的单颗磨粒磨削试验平台及其试验方法 |
Non-Patent Citations (1)
Title |
---|
"基于信息熵的马氏体相变声发射信号分析";戴光等;《金属热处理》;20120430;第37卷(第4期);第113-116页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107153728A (zh) | 2017-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107153728B (zh) | 一种基于双谱幅值分布熵的砂轮磨损状态特征提取方法 | |
CN109465676B (zh) | 一种刀具寿命预测方法 | |
CN109262368B (zh) | 一种刀具失效判定方法 | |
CN112032003B (zh) | 一种大型风电机组运行性能监测方法 | |
CN112171376B (zh) | 一种基于电流信号分割的机床工件实时统计方法 | |
CN111784093B (zh) | 一种基于电力大数据分析的企业复工辅助判断方法 | |
CN112380992B (zh) | 一种加工过程监控数据准确性评估与优化方法及装置 | |
CN116823826B (zh) | 一种数控机床刀具崩刃异常检测方法 | |
CN103203670B (zh) | 基于最大信息熵与方向散度的磨削颤振预测方法 | |
CN111958321B (zh) | 基于深度神经网络的数控机床刀具磨损程度识别方法 | |
CN116304751B (zh) | 用于动车组部件检修的运行数据处理方法 | |
CN110109431B (zh) | 一种压铸机oee信息智能获取系统 | |
CN112207631A (zh) | 刀具检测模型的生成方法、检测方法、系统、设备及介质 | |
CN110555243B (zh) | 一种铣削过程中机床主轴振动信息的二维图谱构建方法 | |
CN114523338A (zh) | 一种基于噪声分析的切削刀具磨损状态监测方法 | |
CN114326593A (zh) | 一种刀具寿命预测系统及方法 | |
CN101826070A (zh) | 一种基于关键点的数据序列线性拟合方法 | |
CN117116291A (zh) | 一种含沙水流冲击水轮机的声信号处理方法 | |
CN110744359A (zh) | 一种数控车床刀具磨损监测系统及方法 | |
Ferreira et al. | Inferential measurement of the dresser width for the grinding process automation | |
CN112903296B (zh) | 一种滚动轴承故障检测方法及系统 | |
CN115741235A (zh) | 基于五轴加工中心刀具的磨损预测与健康管理方法 | |
CN109426216B (zh) | 一种用于数控铣床的矩形平面加工深度识别方法 | |
CN107784015B (zh) | 一种基于电力系统在线历史数据的数据约简方法 | |
CN112434634A (zh) | 一种快速消除土木工程结构健康监测信号尖峰的方法及系统 |
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 |