一种基于多尺度主元分析的刀具磨损在线监测方法

文档序号:10480255阅读:1048来源:国知局
一种基于多尺度主元分析的刀具磨损在线监测方法
【专利摘要】本发明公开一种基于多尺度主元分析的刀具磨损在线监测方法,将原始数据集分为训练样本和测试样本,通过小波分解将训练样本分解成多个尺度,利用主元分析对各个尺度及整体进行降维处理构建模型;将测试样本依据训练样本的层数进行小波分解,将各尺度下的数据加载到相应尺度的主元模型上计算出各尺度下数据的SPE和T2统计量;对显著尺度上的信号利用小波重构方法组成一个新的测试样本,加载到整体主元模型上,通过计算SPE和T2统计量并判断是否超出控制限,如果超出则说明过程刀具磨损异常,否则为刀具磨损正常,从而监测处刀具的磨损状态。本发明提高了加工过程中刀具磨损的在线识别的准确率,特别是对难加工材料存在复杂非线性信号的具有强的适应性。
【专利说明】
一种基于多尺度主元分析的刀具磨损在线监测方法
技术领域
[0001] 本发明涉及一种机床刀具的状态监测与辨识技术领域,特别是涉及。
【背景技术】
[0002] 目前,刀具磨损状态监测方法主要是基于传感器信号建立信号特征与刀具磨损状 态之间关系模型,并利用模式识别技术间接监测刀具磨损状态的方法。但是这种方法的缺 点是信号中噪声较多、信息量大、计算困难以及磨损模型难建立等。虽然支持向量机(SVM)、 人工神经网络(ANN)以及隐马尔科夫模型(HMM)等模式识别方法在刀具磨损状态监测中取 得了很好的应用效果,但是还存在着很多缺点,尤其是在难加工材料加工过程的刀具磨损 状态监测中。首先,这些方法需要有大量包含多种磨损状态的实验数据来训练模型,这就意 味着一种切削参数下的实验要进行很多次。对于普通材料来说,这样的实验是可行的;但 是,在难加工材料加工过程中,难加工材料不但材料昂贵,而且加工费时费力,刀具极易发 生磨损,我们很难获得所有刀具磨损状态下的数据。更重要的是,刀具磨损状态是多种多样 的,由于训练样本以及模型的复杂性,建模过程是非常耗时的。这一点显然是无法满足实际 加工的在线监测需要。因此提供精确可靠的刀具状态监测系统是解决此类加工难题行之有 效的方法之一。

【发明内容】

[0003] 针对目前刀具状态监测的不足,本发明提出了一种基于多尺度主元分析的刀具磨 损状态监测方法,以主元分析(PCA)为基础结合小波分析的多尺度主元分析(MSPCA)实现刀 具状态的在线监测。
[0004] 本发明提出了一种基于多尺度主元分析(MSPCA)的刀具磨损状态监测方法基于多 尺度主元分析的刀具磨损在线监测方法,该方法包括以下步骤:
[0005] 步骤一,获取在线监测时间内,刀具加工运行的原始信号数据;包括加工过程中刀 具对工件两个侧面和顶面的三个方向的切削力信号数据和振动信号数据,将所述原始信号 数据构成样本集;
[0006] 步骤二,通过对上述采集的原始信号数据进行滤波降噪、提取时域和频域12个特 征参数,将样本集分为训练样本集和测试样本集;
[0007] 步骤三,通过小波分解将训练样本集分解成多个尺度,小波分解后,得到一个近似 矩阵AL和L个细节矩阵,AL为目标层的尺度函数系数组成的向量,L个细节矩阵为每一层产 生的小波函数系数组成的向量;
[0008] 步骤四,利用主元分析对多个尺度下的所述训练样本集的数据分别进行建模,构 建整体主元模型以及计算相应的SPE值及控制限和T 2统计量及控制限;
[0009] 步骤五,步骤五,将测试样本依据训练样本的层数进行小波分解,将各分解尺度下 的数据加载到相应尺度的主元模型上计算出各尺度下数据的SPE和T 2统计量:
[0010] 在i时刻,SPE统计量表示为:
[0011] Smi =x,a-pkp!}xf
[0012]其中,xdi时亥鋼得的数据,I为单位矩阵,Pk为主元模型的负荷矩阵;
[0013 ]在i时刻,T2统计量以及控制限表示为:
[0014] Ti- = XiPlO sPlX]
[0015] 其中,0为前k个主元的特征值组成的对角矩阵,X1Si时刻测得的数据,Pk为主元模 型的负荷矩阵;
[0016] 判断显著尺度,舍弃非显著尺度;
[0017] 步骤六,将显著尺度上的信号利用小波重构方法组成一个新的测试样本;
[0018] 步骤七,所述新的测试样本将加载到整体主元模型上,计算SPE控制限和T 2统计 量,并判断是否超出相应控制图的控制限
[0019] SPE控制限表示为:
角矩阵,k是主元个数,m是变量个数,Ca是正态分布在假设检验水平为α下的临界值;[0022] T2控制限表示为:
[0020
[0021 ,λ为前k个主元的特征值组成的对
[0023]
[0024]其中,k为主元个数,m为变量个数,Fk,m-k,a为对应于检验水平为α、自由度为k、n-l 条件下的F分布临界值;
[0025]如果超出,则说明过程刀具磨损异常;否则,说明过程刀具磨损正常。
[0026] 所述步骤二中的12个特征参数包括6个时域特征指标即包括均值、均方根、峰峰 值、方差、峭度、偏度和6个频域特征指标即功率谱和、功率谱均值、功率谱峰值、功率谱方 差、功率谱峭度、功率谱偏度。
[0027] 所述步骤四的利用主元分析对多个尺度下的所述训练样本集的数据分别进行建 模的步骤,具体包括以下处理:
[0028]获取正常刀具加工过中的信号数据;
[0029]对训练样本集f进行标准化处理;
[0030] 将训练样本集中数据X进行协方差矩阵分解,数据X表示为:
[0031] X = SPt
[0032] 其中:S=[S1,S2,S3,~,Sm]称为得分矩阵,顺序按向量的模的大小排列,提取的是 前后数据间的相互关系。P =[Pl,P2,P3,…,Pm]称为负荷矩阵,反映了变量间的相互关系;
[0033] 计算各主元贡献率,依序将主元的贡献率累加,达到阈值后,确定主元个数。
[0034] 与现有技术相比,本发明充分考虑到刀具磨损多尺度、非线性信号,依据难加工材 料刀具磨损的特点,通过对预处理和提取特征之后的数据进行小波分解,分解到多个尺度 上,克服了主元分析只适合于单尺度主元分析的缺点;再利用主元分析降低各尺度数据的 维数,有效减小变量间的相关性,减少模型计算量。通过主元模型及相应的SPE值及控制限 和T2统计量及控制限,可以准确的识别刀具磨损处于正常或异常的状态;提高了加工过程 中刀具磨损的在线识别的准确率,特别是对难加工材料存在复杂非线性信号的具有强的适 应性,对刀具磨损识别具有重要意义,在提高加工表面质量也具有重要意义。
【附图说明】
[0035] 图1为本发明的一种基于多尺度主元分析的刀具磨损在线监测方法流程图;
[0036] 图2为本发明的一种基于多尺度主元分析的刀具磨损在线监测方法的多尺度主元 分析示意图;
[0037] 图3为本发明一种基于多尺度主元分析的刀具磨损在线监测方法的多离散小波分 解示意图。
【具体实施方式】
[0038] 本发明与其他模式识别方法不同,它不需要训练样本包含多种刀具磨损状态下的 数据,只是利用刀具正常切削时的信号进行建模。其次,本发明利用离散小波分解从多个尺 度上监测了刀具磨损状态是否发生,十分适合处理加工过程中多尺度、非线性、非平稳信 号。更重要的是,该方法基于主元分析(PCA)方法克服了多传感器融合技术带来的信息冗 余、变量相关度高等缺点,降低了数据维度,简化模型运算过程,对于刀具磨损状态在线监 测具有重要意义。
[0039]以下通过待定的具体实例并结合【附图说明】本发明的实施方式,本领域技术人员可 以由本说明书所揭示的内容轻易地了解本发明的其他优点和功效。本发明亦可通过其他不 同的具体实例加以实行或应用,本说明书中的各项细节亦可基于不同观点与应用,在不背 离本发明的精神下进行各种修饰与变更。
[0040] 如图1所示,本发明一种刀具磨损在线监测方法,包括如下步骤:
[0041] 步骤101,获取一定时间内的刀具加工运行的原始信号数据;在该步骤中,所述原 始信号数据包括:1)通过一个测力仪采集加工过程中刀具对工件的力的作用,由于刀具磨 损程度的不同,对工件施加的力是变化的,通过采集的力信号的数据,做数据分析,可以得 到刀具磨损的程度;2)选用三个加速度传感器,分别安装在工件两个侧面和顶面,测量三个 方向的振动信号,用于采集加工过程中工件的振动数据,同样,由于刀具磨损程度不同,产 生的振动程度也不同;在本发明中,通过结合这两种传感器采集的数据信息,对刀具的磨损 程度做分析判断;
[0042]步骤102,对上述采集的数据信息进行滤波降噪和特征提取,得到样本集,将样本 集分为训练样本集和测试样本集;
[0043] 步骤103,选择适当的小波函数对训练样本进行小波分解,分解到多个尺度下;
[0044] 步骤104,对多个尺度下的训练样本数据分别进行主元分析(PCA)建模并且构建整 体主元模型以及相应的SPE值和T2控制限;
[0045] 步骤105,将测试样本依据训练样本层数分解,将各尺度下的数据加载到相应尺度 的主元模型上计算出各尺度下数据的SPE控制限和T 2统计量,判断显著尺度,舍弃非显著尺 度;步骤106,将显著尺度上的信号利用小波重构方法组成一个新的测试样本;
[0046]步骤107,新的测试样本将加载到整体主元模型上,计算SPE控制限和T2统计量并 判断是否超出相应控制图的控制限,如果超出则说明过程刀具磨损异常,否则为正常; [0047] 所述步骤102具体包括一下处理:
[0048]步骤2.1,对力信号进行低通滤波处理,振动信号进行带通滤波处理;
[0049]步骤2.2,提取时域特征参数和频域特征参数;
[0050]步骤2.3,计算6个时域特征指标,该时域特征指标包括均值、均方根、峰峰值、方 差、峭度、偏度;
[0051 ]步骤2.4,计算6个频域特征指标,该频域特征指标包括功率谱和、功率谱均值、功 率谱峰值、功率谱方差、功率谱峭度、功率谱偏度;
[0052]步骤2.5,训练样本的信息为正常磨损状态下的数据;
[0053]步骤2.6,测试样本的信息为一半正常磨损状态下的数据,另一半为异常磨损状态 下的数据;
[0054] 步骤2.7,训练样本和测试样本的数量相同。
[0055] 在时域范围内,即在时间轴上做特征提取。具体是对每一列信号以一定的长度N提 取信号,例如一列信号的长度是m*N长度,我们以长度N来提取信号,那么这一列就可以得到 m个特征数,如果我们有6列数,这提取完特征就的到的是mX6的一个矩阵:
[0056] (1)均值:计算一段信号的幅值平均值。针对信号乂={^1,12,一,幼}其均值表达式 为:
[005:
[0058] N为信号长度,X1每个信号点对应的值)
[0059] (2)均方根:计算一段信号的有效值,其表达式为:
[0060]
[0061] N为信号长度,^每个信号点对应的值;
[0062] (3)峰峰值:计算一段信号中最大值与最小值的差值,其表达式为:
[0063] Pk=max(xi)-min(xi)
[0064] Xi为每个信号点对应的值;
[0065] (4)方差:计算一段信号的方差值,其表达式为:
[0066]
[0067] N为信号长度,乂1每个信号点对应的值,X为均值;
[0068] (5)峭度:计算信号波形陡峭程度的物理量。峭度值的改变反映了系统发生了变 化。在十-说明加工信号越不正常,发生故障的可能性也就越大。其表 达式为
[0069]
[0070] 其中,N为信号长度,^每个信号点对应的值,X为均值,σ标准差;
[0071 ] (6)偏度:描述信号对称性分布的物理量,反映了信号分布倾斜的方向和程度大 小。偏度越大,信号对称性越差。其表达式为:
[0072]
[0073] 其中,N为信号长度,^每个信号点对应的值,X均值,σ标准差。
[0074]信号的频域特征主要反映的是系统在幅值和波动性上的变化,但是却无法获得信 号频域上的信息。因此,本发明采用功率谱来提取相关频域特征。功率谱的计算方法主要包 括参数估计法和非参数估计法。本文采用的是非参数估计法中的Welch法来计算信号的功 率谱。首先进行快速傅立叶变换(FFT数据长度为1024)得到信号的频谱,本文采用加矩形窗 的方法,从而有效防止频谱发生泄漏,然后通过计算得到信号的功率谱密度S x(f)。其中,功 率谱密度横轴代表频率,纵轴代表对应频率的功率值。
[0075] 假设信号父=&1^,~,幼},信号的自功率谱5办)和自相关函数以1)之间有下面 的表达式·
[0076]
[0077] Rx(i)自相关函数,Ts为信号采样间隔,f为频率;
[0078] 其中,%为信号采样间隔。下面介绍本文所选取的频率特征:
[0079] (1)功率谱和:计算信号在一个频段内每个频率所对应的功率值之和。其表达式 为:
[0080]
[00811 &算信号在一定频段内的功率谱均值,其表达式为:
[0082:
[0083] (3)功率谱峰值:计算一定频段内信号功率谱的最大值,其表达式为:
[0084] Pp=max(S(f))
[0085] (4)功率谱方差:计算一定频段内功率谱的方差,其表达式为:
[0091] 所述步骤103包括以下具体处理:[0092]步骤3.1,分解尺度L的选择可以根据L= log2n-5(其中η为样本个数)来确定;
[0086]
[0087]
[0088]
[0089]
[0090]
[0093] 步骤3.2,小波分解后,通过Mallat算法得到一个近似矩阵AL和L个细节矩阵(从D1 到Dl),Mallat算法本质上相当于用低通滤波器和带通滤波器表示尺度函数和小波函数,分 解过程示意如图3所示。其表达式为:
[0094] Ai = HAi-i,Di = GAi-i
[0095] 其中,D1S在尺度i上小波函数系数组成的向量,仏为在尺度i上尺度函数系数组成 的向量,G表示带通滤波器,H表示低通滤波器。图3中以一个三层Mallat算法为例,X为原始 数据,G表示带通滤波器,H表示低通滤波器。1/2为二抽取,抽取后的数据采样率和数据点数 均变为原来的二分之一。首先数据通过第一次的带通滤波和低通滤波处理,经二抽取后得 到第一层的高频系数和低频系数。接着对第一层的低频系数再进行带通滤波和低通滤波处 理,从而得到第二层的高频系数和低频系数。依次类推,分解到目标层数为止。AL则为目标 层的尺度函数系数组成的向量,L个细节矩阵为每一层产生的小波函数系数组成的向量; [0096] 所述步骤104包括以下具体处理:
[0097]步骤4.1,获取正常刀具加工过中的信号数据;
[0098]步骤4.2,数据样本X'进行标准化处理;
[0099]步骤4.3,将数据X进行协方差矩阵分解,数据X可用下式表示:X = SPT
[0100] 步骤4.4,计算各主元贡献率(CPV),将前几个主元的贡献率累加,达到一定阈值后 即可确定主元个数,其阈值一般设定在80% ;
[0101] 步骤4.5,计算SPE统计量以及控制限;
[0102] 步骤4.6,计算T2统计量以及控制限;
[0103] 以下将通过具体实施例说明本发明的步骤:
[0104] 步骤1,获取原始信号数据。通过测力仪和加速度传感器分别测量三个方向切削力 信号和三个方向振动信号,要求刀具正常磨损的信号数据是异常磨损得信号数据的3倍。
[0105] 步骤2,对数据进行滤波降噪和特征提取,得到样本集,将样本集分为训练样本集 和测试样本集。对切削力信号做低通滤波处理,振动信号做带通滤波处理。对完成滤波的信 号提取时域和频域特征。提取6个时域特征指标包括均值、均方根、峰峰值、方差、峭度、偏 度。提取6个频域特征指标包括功率谱和、功率谱均值、功率谱峰值、功率谱方差、功率谱峭 度、功率谱偏度。取正常磨损数据的三分之二为训练样本,正常磨损数据剩下的三分之一和 异常磨损信号组成测试样本。此时训练样本中的数据全部为刀具正常磨损信号数据,训练 样本的数据前一半为正常磨损信号数据,后一半为异常磨损信号数据,且两组数据量相同。
[0106] 步骤3,选择适当的小波函数对训练样本进行小波分解,分解到多个尺度下。通过 多分辨率小波分析(Mallat算法)分解训练样本,分解层数L=log 2n-5(其中η为样本个数)。 分解后得到一个近似矩阵AL和L个细节矩阵(从Dl到DL)(L已经表示过分解层数,因此,分解 L层则得到L个细节矩阵(如图3所示,分解三层则得到三个细节矩阵,Dl、D2、D3))。
[0107] 步骤4,对多个尺度下的训练样本数据分别进行主元分析(PCA)建模并且构建整体 主元模型以及相应的SPE和T2控制限。利用主元分析对训练样本计算各个主元的贡献率 (CPV),累加前几个主元的贡献率达到80%,则认为这些主元信息包括大部分的训练样本的 信息,从而大幅度降低计算难度。经降维处理后,各尺度所选择的主元向量被用来构建主元 模型。主元模型主要包括负载矩阵和得分矩阵。负载矩阵主要用来将新的测量数据投影到 各主元方向,得分矩阵则是用来计算SPE值和T 2统计量及各尺度下的控制限。同时,为了实 现最后的结果监测,仅仅依靠各尺度下的控制限是不够的。我们还需要针对刀具正常磨损 状态下的整体数据进行主元分析,构建整体主元模型以及相应的SPE值和T2控制限。
[0108] 步骤5,将测试样本依据训练样本层数分解,将各尺度下的数据加载到相应尺度的 主元模型上计算出各尺度下数据的SPE和T2统计量,判断显著尺度,舍弃非显著尺度。训练 样本模型构建过程,得到了 L个尺度下的主元模型和一个整体主元模型以及相对应的控制 限,首先将测试数据进行小波分解,这里分解层数与建模过程相同。然后,将各尺度下的数 据加载到相应尺度的主元模型上计算出各尺度下数据的SPE和T 2统计量,判断哪些尺度为 非显著尺度,在小波重构时舍弃非显著尺度。
[0109] 步骤6,通过步骤5,选择出测试样本中的显著尺度。将显著尺度上的信号利用小波 重构方法组成一个新的测试样本。
[0110] 步骤7,新的测试样本将加载到步骤4的整体主元模型上,计算SPE值和T2统计量并 判断是否超出相应控制图的控制限,如果超出则说明过程刀具磨损异常,否则为刀具磨损 正常,从而监测处刀具的磨损状态。
[0111] 综上所述,本发明一种刀具磨损在线监测方法是根据加工过程刀具磨损的特点, 即存在多尺度、非线性、非平稳信号,特别是难加工材料加工过程中,刀具极易发生磨损,我 们很难获得所有刀具磨损状态下的数据,多尺度主元分析(MSPCA)结合了主元分析(PCA)和 小波分析,前者具有强大的数据降维以及提取变量间线性关系的能力,可以降低变量间的 线性相关性、减小模型计算量,后者可以提取变量多尺度特征和变量自相关性,克服了主元 分析只适用于单尺度的局限,提高了刀具磨损在线监测的准确率,对于刀具磨损状态在线 监测具有重要意义。
[0112]上述实例仅示例性说明本发明的原理及其功效,而非用于限制本发明。任何本发 明领域技术人员均可在不违背本发明的精神及范畴下,对上述实例进行修饰与改变。因此, 本发明的权利保护范围,应如权利要求书所列。
【主权项】
1. 一种基于多尺度主元分析的刀具磨损在线监测方法,其特征在于,该方法包括以下 步骤: 步骤一,获取在线监测时间内,刀具加工运行的原始信号数据;包括加工过程中刀具对 工件两个侧面和顶面的三个方向的切削力信号数据和振动信号数据,将所述原始信号数据 构成样本集; 步骤二,通过对上述采集的原始信号数据进行滤波降噪、提取时域和频域12个特征参 数,将样本集分为训练样本集和测试样本集; 步骤三,通过小波分解将训练样本集分解成多个尺度,小波分解后,得到一个近似矩阵 AL和L个细节矩阵,AL为目标层的尺度函数系数组成的向量,L个细节矩阵为每一层产生的 小波函数系数组成的向量; 步骤四,利用主元分析对多个尺度下的所述训练样本集的数据分别进行建模,构建整 体主元模型以及计算相应的SPE值及控制限和T2统计量及控制限; 步骤五,将测试样本依据训练样本的层数进行小波分解,将各分解尺度下的数据加载 到相应尺度的主元模型上计算出各尺度下数据的SPE控制限和T2统计量: 在i时刻,SPE控制限表示为:其中,X1Si时刻测得的数据,I为单位矩阵,Pk为主元模型的负荷矩阵; 在i时刻,T2统计量以及控制限弃·其中,〇为前k个主元的特征值组成的对角矩阵,X1Si时刻测得的数据,Pk为主元模型的 负荷矩阵; 判断显著尺度,舍弃非显著尺度; 步骤六,将显著尺度上的信号利用小波重构方法组成一个新的测试样本; 步骤七,所述新的测试样本将加载到整体主元模型上,计算SPE控制限和T2统计量,并 判断是否超出相应控制图的控制限 SPE控制限表示为: 阵其中,k为主元个数,m为变量个数,Fk,m-k,a为对应于检验水平为α、自由度为k、n-l条件下 的F分布临界值; 如果超出,则说明过程刀具磨损异常;否则,说明过程刀具磨损正常。2. 如权利要求1所述的一种基于多尺度主元分析的刀具磨损在线监测方法,其特征在 于,所述步骤二中的12个特征参数包括6个时域特征指标即包括均值、均方根、峰峰值、方 差、峭度、偏度和6个频域特征指标即功率谱和、功率谱均值、功率谱峰值、功率谱方差、功率 谱峭度、功率谱偏度。3. 如权利要求1所述的一种基于多尺度主元分析的刀具磨损在线监测方法,其特征在 于,所述步骤四的利用主元分析对多个尺度下的所述训练样本集的数据分别进行建模的步 骤,具体包括以下处理: 获取正常刀具加工过中的信号数据; 对训练样本集f进行标准化处理; 将训练样本集中数据X进行协方差矩阵分解,数据X可用下式表示: X = SPt 其中:S=[S1,S2,S3,~,Sm]称为得分矩阵,顺序按向量的模的大小排列,提取的是前后 数据间的相互关系。P =[Pl,P2,P3,…,Pm]称为负荷矩阵,反映了变量间的相互关系; 计算各主元贡献率,依序将主元的贡献率累加,达到阈值后,确定主元个数。
【文档编号】B23Q17/09GK105834835SQ201610269068
【公开日】2016年8月10日
【申请日】2016年4月26日
【发明人】王国锋, 吴丽蕊, 杨星焕, 宋庆月
【申请人】天津大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1