CN106564012B - 一种磨削加工颤振的检测方法 - Google Patents
一种磨削加工颤振的检测方法 Download PDFInfo
- Publication number
- CN106564012B CN106564012B CN201610936478.4A CN201610936478A CN106564012B CN 106564012 B CN106564012 B CN 106564012B CN 201610936478 A CN201610936478 A CN 201610936478A CN 106564012 B CN106564012 B CN 106564012B
- Authority
- CN
- China
- Prior art keywords
- signal
- flutter
- grinding
- autocorrelation sequence
- unbiased
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24B—MACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
- B24B49/00—Measuring or gauging equipment for controlling the feed movement of the grinding tool or work; Arrangements of indicating or measuring equipment, e.g. for indicating the start of the grinding operation
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/12—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring vibration
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24B—MACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
- B24B49/00—Measuring or gauging equipment for controlling the feed movement of the grinding tool or work; Arrangements of indicating or measuring equipment, e.g. for indicating the start of the grinding operation
- B24B49/006—Measuring or gauging equipment for controlling the feed movement of the grinding tool or work; Arrangements of indicating or measuring equipment, e.g. for indicating the start of the grinding operation taking regard of the speed
Abstract
一种磨削加工颤振检测方法,首先对测得的振动信号进行分段、加窗处理;然后,分别对每段信号进行最小熵解卷积滤波处理,消除传播路径对颤振特征信号的影响;接下来,通过无偏自相关序列计算,消除突发性干扰信号;最后,使用峭度指标对颤振特征信号进行量化,设置颤振阈值,通过颤振阈值跟量化后的颤振特征信号进行比较,确定磨床的颤振状态,本发明方法具有通用性和抗干扰性,颤振阈值不受加工设备,加工参数以及测试信号中突发性信号的影响,适合磨削颤振在线检测,可在机床企业推广使用,符合智能机床发展需求,具有广泛的应用前景。
Description
技术领域
本发明涉及机床加工过程颤振检测技术领域,特别涉及一种磨削加工颤振的检测方法。
背景技术
磨削加工作为常见的终加工工艺,对产品几何精度和表面质量具有决定性影响。磨削颤振是影响加工效率和表面质量最为严重的因素之一。磨削颤振的发生会加剧砂轮磨损,产生异常噪声,削减机床部件寿命,降低产品加工精度和表面质量。磨削颤振行为的演变具有阶段性,通常包括早期、上升期和稳定发展三个阶段。在磨削加工过程中,只有在颤振早期实现检测并采取抑制措施,方能避免颤振对产品加工精度和表面质量产生破坏。目前,已有的磨削颤振在线检测方法通用性差,主要体现在颤振阈值对加工设备和加工参数敏感。其原因在于:外接传感器测得的颤振信号是颤振特征信号与传播路径单元(轴承,轴,传感器等)卷积的结果,其中的颤振特征信号已被污染,信噪比降低。因此,直接使用低信噪比的测试信号难以实现早期颤振检测,消除传播路径对颤振特征信号的卷积影响,获得高信噪比的颤振特征信号才是解决该问题的关键。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种磨削加工颤振的检测方法,实现磨削颤振的早期检测,避免磨削颤振对产品加工质量的影响。
为了达到上述目的,本发明采取的技术方案为:
一种磨削加工颤振的检测方法,包括以下步骤:
步骤一,在磨床的砂轮架上安装一个加速度传感器;
步骤二,利用数据采集设备和加速度传感器,实时采集加速度传感器测得的振动信号x(t),其中,t代表时间;
步骤三,对振动信号x(t)进行分段、加窗处理,得到振动信号x(t1),x(t2),…,x(tm),其中,t1,t2,…,tm分别代表分段、加窗处理后信号所在的时间序列,m代表分段的段数;
步骤四,使用Ralph Wiggins的最小熵解卷积滤波方法对分段、加窗处理后的振动信号x(t1),x(t2),…,x(tm)分别进行解卷积处理,得到解卷积后信号y(t1),y(t2),…,y(tm),最小熵解卷积滤波方法的实现步骤如下:
①初始化一个输入信号的滤波器系数设置
②计算输入信号的Toeplitz自相关矩阵X0,计算公式为:
其中,N代表输入信号x向量的长度,L代表Toeplitz自相关矩阵X0的行数;
③计算输出信号计算公式为其中,代表自相关矩阵X0的转置矩阵,N代表y向量的长度;
④重新计算输入信号的滤波器系数计算公式为:
其中,n代表y向量中的第n个离散点,N代表y向量的长度;
⑤通过迭代方法,迭代计算②~④步骤,迭代结束的条件为:连续两次输出信号的峭度指标之差小于0.01,其中,峭度指标的计算公式为:
|·|代表绝对值运算,代表平均值运算;
⑥将第⑤步中的最后一组输出信号记为最小熵解卷积滤波方法的输出结果;
步骤五,使用无偏自相关序列的计算公式,分别计算解卷积后信号y(t1),y(t2),…,y(tm)的无偏自相关序列,得到R(τ1),R(τ2),…,R(τm),无偏自相关序列的计算公式为: 其中,n代表y向量中的第n个离散点,τ代表无偏自相关序列中的第τ个离散点,无偏自相关序列的长度为2N-1;
步骤六,分别截取无偏自相关序列R(τ1),R(τ2),…,R(τm)中长度为的一段信号,并分别计算每段信号的峭度指标,得到峭度指标序列K1,K2,…,Km,其中,无偏自相关序列R中τ的取值范围为
步骤七,磨床未发生颤振时,通过上述六个步骤得到的无偏自相关序列R(τ1),R(τ2),…,R(τm)都不存在明显冲击特征,其峭度指标K1,K2,…,Kn的值都靠近3;而磨床发生颤振时,通过上述六个步骤得到的无偏自相关序列R(τ1),R(τ2),…,R(τm)都存在明显冲击特征,其峭度指标K1,K2,…,Km的值都大于4;因此设置一个初始的颤振阈值T=4,如果计算得到的峭度指标序列K1,K2,…,Km中任何一个值大于设置的颤振阈值,那么判定该磨削过程发生颤振。
所述的颤振阈值T根据需求进行修订。
本发明的有益效果:
本发明根据磨削颤振信号的频谱具有在磨床工作台固有频率处存在转频为间隔的频率特点(颤振特征信号),考虑到传感器测得的颤振信号是颤振特征信号与传播路径卷积的结果,提出一种考虑传播路径影响的磨削加工颤振检测方法,即基于最小熵解卷积与无偏自相关的颤振检测方法,并通过齿轮磨削加工对该方法的有效性进行了验证。本发明方法可有效检测早期磨削颤振行为,相比于现有方法,该磨削颤振检测方法具有通用性和抗干扰性,颤振阈值不受加工设备,加工参数和测试信号中突变信号的影响,适合磨削颤振在线检测,可在机床企业推广使用,符合智能机床发展需求,具有广泛的应用前景。
附图说明
图1为本发明方法的流程图。
图2为实施例测得的磨床加工过程振动信号x(t)时域图。
图3为实施例磨床加工过程振动信号x(t)经本发明处理后得到的无偏自相关序列R(τ1),R(τ2),…,R(τm)时域图。
图4为实施例磨床加工过程振动信号x(t)经本发明处理后得到的峭度指标序列时域图。
具体实施方式
下面结合附图和实施案例对本发明进行详细描述。
某机床厂生产的齿轮磨床在齿轮加工过程中,出现异常的加工噪音,生产的齿轮在下线检测时发现加工质量不符合要求。经机床厂研究,齿轮下线检测不合格的主要原因是齿轮磨床在加工时出现了颤振行为;如果早期颤振行为可以被检测,那么调整加工参数就可以避免颤振行为对工件加工质量的影响,下面通过本发明解决早期颤振检测难题。
参照图1,一种磨削加工颤振的检测方法,包括以下步骤:
步骤一,在磨床的砂轮架上安装一个加速度传感器;
步骤二,利用数据采集设备和加速度传感器,实时采集加速度传感器测得的振动信号x(t),如图2所示,图2给出了加速度传感器测得的时域信号,从图2中可以观察到该磨削过程包括4个磨削阶段,并且测得的振动信号中存在多个突发性信号干扰,其中,t代表时间;
步骤三,对振动信号x(t)进行分段、加窗处理,得到振动信号x(t1),x(t2),…,x(tm),其中,t1,t2,…,tm分别代表分段、加窗处理后信号所在的时间序列,m代表分段的段数;
步骤四,使用RalphWiggins的最小熵解卷积滤波方法对分段、加窗处理后的振动信号x(t1),x(t2),…,x(tm)分别进行解卷积处理,消除传播路径对颤振特征信号的影响,得到解卷积后信号y(t1),y(t2),…,y(tm),最小熵解卷积滤波方法的实现步骤如下:
①初始化一个输入信号的滤波器系数设置
②计算输入信号的Toeplitz自相关矩阵X0,计算公式为:
其中,N代表输入信号x向量的长度,L代表Toeplitz自相关矩阵X0的行数;
③计算输出信号计算公式为其中,代表自相关矩阵X0的转置矩阵,N代表y向量的长度;
④重新计算输入信号的滤波器系数计算公式为:
其中,n代表y向量中的第n个离散点,N代表y向量的长度;
⑤通过迭代方法,迭代计算②~④步骤,迭代结束的条件为:连续两次输出信号的峭度指标之差小于0.01,其中,峭度指标的计算公式为:
|·|代表绝对值运算,代表平均值运算,n代表y向量中的第n个离散点,N代表y向量的长度;
⑥将第⑤步中的最后一组输出信号记为最小熵解卷积滤波方法的输出结果;
步骤五,使用无偏自相关序列的计算公式,分别计算解卷积后信号y(t1),y(t2),…,y(tm)的无偏自相关序列,消除突发性信号的干扰,得到R(τ1),R(τ2),…,R(τm),无偏自相关序列的计算公式为:其中,n代表y向量中的第n个离散点,τ代表无偏自相关序列中的第τ个离散点,无偏自相关序列的长度为2N-1;
步骤六,分别截取无偏自相关序列R(τ1),R(τ2),…,R(τm)中长度为的一段信号,并分别计算每段信号的峭度指标,得到峭度指标序列K1,K2,…,Km,如图3所示,图3为经过本步骤得到的无偏自相关序列,从图3中可以看到,磨削第1阶段不包含冲击信号,磨削第2,3,4阶段出现明显的周期性冲击信号,这说明磨削状态从第2阶段开始发生了变化,其中,无偏自相关序列R中τ的取值范围为
步骤七,磨床未发生颤振时,通过上述六个步骤得到的无偏自相关序列R(τ1),R(τ2),…,R(τm)都不存在明显冲击特征,其峭度指标K1,K2,…,Kn的值都靠近3;而磨床发生颤振时,通过上述六个步骤得到的无偏自相关序列R(τ1),R(τ2),…,R(τm)都存在明显冲击特征,其峭度指标K1,K2,…,Km的值都大于4;因此设置一个颤振阈值T=4,如图4所示,图4为经本步骤处理后得到的峭度指标序列时域图,从图4中可以观察到,当峭度指标序列在56秒时,其峭度指标大于阈值4,说明齿轮磨床在第2阶段中发生了早期磨削颤振,此时可对磨削颤振进行控制,如更换砂轮、改变加工参数,从而使磨削过程处于无颤振加工状态。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施方式仅限于此,对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单的推演或替换,都应当视为属于本发明由所提交的权利要求书确定专利的保护范围。
Claims (2)
1.一种磨削加工颤振的检测方法,其特征在于,包括以下步骤:
步骤一,在磨床的砂轮架上安装一个加速度传感器;
步骤二,利用数据采集设备和加速度传感器,实时采集加速度传感器测得的振动信号x(t),其中,t代表时间;
步骤三,对振动信号x(t)进行分段、加窗处理,得到振动信号x(t1),x(t2),…,x(tm),其中,t1,t2,…,tm分别代表分段、加窗处理后信号所在的时间序列,m代表分段的段数;
步骤四,使用Ralph Wiggins的最小熵解卷积滤波方法对分段、加窗处理后的振动信号x(t1),x(t2),…,x(tm)分别进行解卷积处理,得到解卷积后信号y(t1),y(t2),…,y(tm),最小熵解卷积滤波方法的实现步骤如下:
①初始化一个输入信号的滤波器系数设置
②计算输入信号的Toeplitz自相关矩阵X0,计算公式为:
其中,N代表输入信号x向量的长度,L代表Toeplitz自相关矩阵X0的行数;
③计算输出信号计算公式为其中,代表自相关矩阵X0的转置矩阵,N代表y向量的长度;
④重新计算输入信号的滤波器系数计算公式为:
其中,n代表y向量中的第n个离散点,N代表y向量的长度;
⑤通过迭代方法,迭代计算②~④步骤,迭代结束的条件为:连续两次输出信号的峭度指标之差小于0.01,其中,峭度指标的计算公式为:
|·|代表绝对值运算,代表平均值运算;
⑥将第⑤步中的最后一组输出信号记为最小熵解卷积滤波方法的输出结果;
步骤五,使用无偏自相关序列的计算公式,分别计算解卷积后信号y(t1),y(t2),…,y(tm)的无偏自相关序列,得到R(τ1),R(τ2),…,R(τm),无偏自相关序列的计算公式为: 其中,n代表y向量中的第n个离散点,τ代表无偏自相关序列中的第τ个离散点,无偏自相关序列的长度为2N-1;
步骤六,分别截取无偏自相关序列R(τ1),R(τ2),…,R(τm)中长度为的一段信号,并分别计算每段信号的峭度指标,得到峭度指标序列K1,K2,…,Km,其中,无偏自相关序列R中τ的取值范围为
步骤七,磨床未发生颤振时,通过上述六个步骤得到的无偏自相关序列R(τ1),R(τ2),…,R(τm)都不存在明显冲击特征,其峭度指标K1,K2,…,Kn的值都靠近3;而磨床发生颤振时,通过上述六个步骤得到的无偏自相关序列R(τ1),R(τ2),…,R(τm)都存在明显冲击特征,其峭度指标K1,K2,…,Km的值都大于4;因此设置一个初始的颤振阈值T=4,如果计算得到的峭度指标序列K1,K2,…,Km中任何一个值大于设置的颤振阈值,那么判定该磨削过程发生颤振。
2.根据权利要求1所述的一种磨削加工颤振的检测方法,其特征在于:所述的颤振阈值T根据需求进行修订。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610936478.4A CN106564012B (zh) | 2016-11-01 | 2016-11-01 | 一种磨削加工颤振的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610936478.4A CN106564012B (zh) | 2016-11-01 | 2016-11-01 | 一种磨削加工颤振的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106564012A CN106564012A (zh) | 2017-04-19 |
CN106564012B true CN106564012B (zh) | 2018-08-21 |
Family
ID=60414294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610936478.4A Active CN106564012B (zh) | 2016-11-01 | 2016-11-01 | 一种磨削加工颤振的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106564012B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107150260B (zh) * | 2017-05-17 | 2018-11-06 | 武汉工程大学 | 一种直角切削颤振解析建模方法 |
CN107560718B (zh) * | 2017-07-25 | 2020-01-10 | 苏州微著设备诊断技术有限公司 | 一种基于自相关峭度最大化的解卷积方法 |
CN111695582A (zh) * | 2019-04-10 | 2020-09-22 | 燕山大学 | 一种颤振纹理的检测方法及其装置 |
CN110091217B (zh) * | 2019-04-29 | 2021-04-27 | 长沙理工大学 | 一种基于机器视觉的车削颤振加速度采集系统及方法 |
CN111890126B (zh) * | 2020-07-03 | 2022-03-11 | 太原理工大学 | 一种基于声压能量峭度指标的早期车削颤振预警监测方法 |
CN113702042B (zh) * | 2021-07-14 | 2022-10-28 | 西南交通大学 | 基于一种增强最小熵解卷积的机械故障诊断方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103203670A (zh) * | 2013-04-03 | 2013-07-17 | 同济大学 | 基于最大信息熵与方向散度的磨削颤振预测方法 |
CN104390697A (zh) * | 2014-11-06 | 2015-03-04 | 西安交通大学 | 一种基于c0复杂度与相关系数的铣削颤振检测方法 |
CN105108584A (zh) * | 2015-07-21 | 2015-12-02 | 上海交通大学 | 一种车削颤振检测方法 |
CN105510032A (zh) * | 2015-12-11 | 2016-04-20 | 西安交通大学 | 基于谐躁比指导的解卷积方法 |
CN106021906A (zh) * | 2016-05-16 | 2016-10-12 | 常州信息职业技术学院 | 基于倒频谱分析的颤振在线监测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI322884B (en) * | 2007-03-27 | 2010-04-01 | Ind Tech Res Inst | Singal analysis method for vibratory interferometry |
-
2016
- 2016-11-01 CN CN201610936478.4A patent/CN106564012B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103203670A (zh) * | 2013-04-03 | 2013-07-17 | 同济大学 | 基于最大信息熵与方向散度的磨削颤振预测方法 |
CN104390697A (zh) * | 2014-11-06 | 2015-03-04 | 西安交通大学 | 一种基于c0复杂度与相关系数的铣削颤振检测方法 |
CN105108584A (zh) * | 2015-07-21 | 2015-12-02 | 上海交通大学 | 一种车削颤振检测方法 |
CN105510032A (zh) * | 2015-12-11 | 2016-04-20 | 西安交通大学 | 基于谐躁比指导的解卷积方法 |
CN106021906A (zh) * | 2016-05-16 | 2016-10-12 | 常州信息职业技术学院 | 基于倒频谱分析的颤振在线监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106564012A (zh) | 2017-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106564012B (zh) | 一种磨削加工颤振的检测方法 | |
Li et al. | Wavelet transform-based higher-order statistics for fault diagnosis in rolling element bearings | |
US8942939B2 (en) | Real-time detection system and the method thereof | |
CN106706118B (zh) | 基于盲源分离技术的运行工况传递路径分析方法 | |
Nair et al. | Permutation entropy based real-time chatter detection using audio signal in turning process | |
CN105058165A (zh) | 基于振动信号的刀具磨损量监测系统 | |
CN104990503A (zh) | 一种基于激光位移传感器的砂轮外圆跳动检测方法 | |
CN104657526B (zh) | 工具机的加工品质的预测方法 | |
CN101829951A (zh) | 外圆磨削工件表面质量可视化实时监测方法 | |
TW201332708A (zh) | 刀具顫振監控方法 | |
CN108205432B (zh) | 一种观测实验数据异常值的实时剔除方法 | |
CN104390697A (zh) | 一种基于c0复杂度与相关系数的铣削颤振检测方法 | |
CN114326594A (zh) | 一种基于机床数据的机床性能测评的系统及其测评方法 | |
Chen et al. | A pseudo wavelet system-based vibration signature extracting method for rotating machinery fault detection | |
Xu et al. | An enhanced multipoint optimal minimum entropy deconvolution approach for bearing fault detection of spur gearbox | |
Maamar et al. | Modal identification of a machine tool structure during machining operations | |
CN106525223A (zh) | 一种齿轮传动装置异响的下线检测方法 | |
Junior et al. | Feature extraction using frequency spectrum and time domain analysis of vibration signals to monitoring advanced ceramic in grinding process | |
CN112926698A (zh) | 一种大型旋转装备振动预测与装配评价方法 | |
Xu et al. | Peak-based mode decomposition for weak fault feature enhancement and detection of rolling element bearing | |
Liu et al. | Chatter detection in milling process based on time-frequency analysis | |
TW201521954A (zh) | 用於顫振監控之訊噪分離方法與裝置 | |
Wang et al. | Study of an efficient real-time monitoring and control system for BUE and cutter breakage for CNC machine tools | |
JP2018111171A (ja) | 異常予兆検知システム及び異常予兆検知方法 | |
CN107636583B (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 |