CN105834835B - 一种基于多尺度主元分析的刀具磨损在线监测方法 - Google Patents
一种基于多尺度主元分析的刀具磨损在线监测方法 Download PDFInfo
- Publication number
- CN105834835B CN105834835B CN201610269068.9A CN201610269068A CN105834835B CN 105834835 B CN105834835 B CN 105834835B CN 201610269068 A CN201610269068 A CN 201610269068A CN 105834835 B CN105834835 B CN 105834835B
- Authority
- CN
- China
- Prior art keywords
- scale
- data
- principal component
- training sample
- sample set
- 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
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 33
- 238000012544 monitoring process Methods 0.000 title claims abstract description 32
- 238000012549 training Methods 0.000 claims abstract description 39
- 238000012360 testing method Methods 0.000 claims abstract description 30
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 19
- 238000001228 spectrum Methods 0.000 claims description 39
- 239000011159 matrix material Substances 0.000 claims description 25
- 238000012545 processing Methods 0.000 claims description 21
- 238000001914 filtration Methods 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 11
- 238000003754 machining Methods 0.000 claims description 10
- 230000002159 abnormal effect Effects 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 5
- 239000000463 material Substances 0.000 abstract description 9
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000005299 abrasion Methods 0.000 description 6
- 238000000605 extraction Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012567 pattern recognition method Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
Classifications
-
- 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/09—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool
- B23Q17/0952—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool during machining
- B23Q17/0957—Detection of tool breakage
-
- 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/09—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool
- B23Q17/0952—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool during machining
- B23Q17/0971—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool during machining by measuring mechanical vibrations of parts of the machine
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开一种基于多尺度主元分析的刀具磨损在线监测方法,将原始数据集分为训练样本和测试样本,通过小波分解将训练样本分解成多个尺度,利用主元分析对各个尺度及整体进行降维处理构建模型;将测试样本依据训练样本的层数进行小波分解,将各尺度下的数据加载到相应尺度的主元模型上计算出各尺度下数据的SPE和T2统计量;对显著尺度上的信号利用小波重构方法组成一个新的测试样本,加载到整体主元模型上,通过计算SPE和T2统计量并判断是否超出控制限,如果超出则说明过程刀具磨损异常,否则为刀具磨损正常,从而监测处刀具的磨损状态。本发明提高了加工过程中刀具磨损的在线识别的准确率,特别是对难加工材料存在复杂非线性信号的具有强的适应性。
Description
技术领域
本发明涉及一种机床刀具的状态监测与辨识技术领域,特别是涉及。
背景技术
目前,刀具磨损状态监测方法主要是基于传感器信号建立信号特征与刀具磨损状态之间关系模型,并利用模式识别技术间接监测刀具磨损状态的方法。但是这种方法的缺点是信号中噪声较多、信息量大、计算困难以及磨损模型难建立等。虽然支持向量机(SVM)、人工神经网络(ANN)以及隐马尔科夫模型(HMM)等模式识别方法在刀具磨损状态监测中取得了很好的应用效果,但是还存在着很多缺点,尤其是在难加工材料加工过程的刀具磨损状态监测中。首先,这些方法需要有大量包含多种磨损状态的实验数据来训练模型,这就意味着一种切削参数下的实验要进行很多次。对于普通材料来说,这样的实验是可行的;但是,在难加工材料加工过程中,难加工材料不但材料昂贵,而且加工费时费力,刀具极易发生磨损,我们很难获得所有刀具磨损状态下的数据。更重要的是,刀具磨损状态是多种多样的,由于训练样本以及模型的复杂性,建模过程是非常耗时的。这一点显然是无法满足实际加工的在线监测需要。因此提供精确可靠的刀具状态监测系统是解决此类加工难题行之有效的方法之一。
发明内容
针对目前刀具状态监测的不足,本发明提出了一种基于多尺度主元分析的刀具磨损状态监测方法,以主元分析(PCA)为基础结合小波分析的多尺度主元分析(MSPCA)实现刀具状态的在线监测。
本发明提出了一种基于多尺度主元分析(MSPCA)的刀具磨损状态监测方法基于多尺度主元分析的刀具磨损在线监测方法,该方法包括以下步骤:
步骤一,获取在线监测时间内,刀具加工运行的原始信号数据;包括加工过程中刀具对工件两个侧面和顶面的三个方向的切削力信号数据和振动信号数据,将所述原始信号数据构成样本集;
步骤二,通过对上述采集的原始信号数据进行滤波降噪、提取时域和频域12个特征参数,将样本集分为训练样本集和测试样本集;
步骤三,通过小波分解将训练样本集分解成多个尺度,小波分解后,得到一个近似矩阵AL和L个细节矩阵,AL为目标层的尺度函数系数组成的向量,L个细节矩阵为每一层产生的小波函数系数组成的向量;
步骤四,利用主元分析对多个尺度下的所述训练样本集的数据分别进行建模,构建整体主元模型以及计算相应的SPE值及控制限和T2统计量及控制限;
步骤五,步骤五,将测试样本依据训练样本的层数进行小波分解,将各分解尺度下的数据加载到相应尺度的主元模型上计算出各尺度下数据的SPE和T2统计量:
在i时刻,SPE统计量表示为:
其中,Xi为i时刻测得的数据,I为单位矩阵,Pk为主元模型的负荷矩阵;
在i时刻,T2统计量以及控制限表示为:
其中,O为前k个主元的特征值组成的对角矩阵,Xi为i时刻测得的数据,Pk为主元模型的负荷矩阵;
判断显著尺度,舍弃非显著尺度;
步骤六,将显著尺度上的信号利用小波重构方法组成一个新的测试样本;
步骤七,所述新的测试样本将加载到整体主元模型上,计算SPE控制限和T2统计量,并判断是否超出相应控制图的控制限;
SPE控制限表示为:
其中,λ为前k个主元的特征值组成的对角矩阵,k是主元个数,m是变量个数,Cα是正态分布在假设检验水平为α下的临界值;
T2控制限表示为:
其中,k为主元个数,m为变量个数,Fk,m-k,α为对应于检验水平为α、自由度为k、n-1个样本条件下的F分布临界值;
如果超出,则说明过程刀具磨损异常;否则,说明过程刀具磨损正常。
所述步骤二中的12个特征参数包括6个时域特征指标即包括均值、均方根、峰峰值、方差、峭度、偏度和6个频域特征指标即功率谱和、功率谱均值、功率谱峰值、功率谱方差、功率谱峭度、功率谱偏度。
所述步骤四的利用主元分析对多个尺度下的所述训练样本集的数据分别进行建模的步骤,具体包括以下处理:
获取正常刀具加工过中的信号数据;
对训练样本集X′进行标准化处理;
将训练样本集中数据X进行协方差矩阵分解,数据X表示为:
X=SPT
其中:S=[s1,s2,s3,…,sm]称为得分矩阵,顺序按向量的模的大小排列,提取的是前后数据间的相互关系。P=[p1,p2,p3,…,pm]称为负荷矩阵,反映了变量间的相互关系;
计算各主元贡献率,依序将主元的贡献率累加,达到阈值后,确定主元个数。
与现有技术相比,本发明充分考虑到刀具磨损多尺度、非线性信号,依据难加工材料刀具磨损的特点,通过对预处理和提取特征之后的数据进行小波分解,分解到多个尺度上,克服了主元分析只适合于单尺度主元分析的缺点;再利用主元分析降低各尺度数据的维数,有效减小变量间的相关性,减少模型计算量。通过主元模型及相应的SPE值及控制限和T2统计量及控制限,可以准确的识别刀具磨损处于正常或异常的状态;提高了加工过程中刀具磨损的在线识别的准确率,特别是对难加工材料存在复杂非线性信号的具有强的适应性,对刀具磨损识别具有重要意义,在提高加工表面质量也具有重要意义。
附图说明
图1为本发明的一种基于多尺度主元分析的刀具磨损在线监测方法流程图;
图2为本发明的一种基于多尺度主元分析的刀具磨损在线监测方法的多尺度主元分析示意图;
图3为本发明一种基于多尺度主元分析的刀具磨损在线监测方法的多离散小波分解示意图。
具体实施方式
本发明与其他模式识别方法不同,它不需要训练样本包含多种刀具磨损状态下的数据,只是利用刀具正常切削时的信号进行建模。其次,本发明利用离散小波分解从多个尺度上监测了刀具磨损状态是否发生,十分适合处理加工过程中多尺度、非线性、非平稳信号。更重要的是,该方法基于主元分析(PCA)方法克服了多传感器融合技术带来的信息冗余、变量相关度高等缺点,降低了数据维度,简化模型运算过程,对于刀具磨损状态在线监测具有重要意义。
以下通过待定的具体实例并结合附图说明本发明的实施方式,本领域技术人员可以由本说明书所揭示的内容轻易地了解本发明的其他优点和功效。本发明亦可通过其他不同的具体实例加以实行或应用,本说明书中的各项细节亦可基于不同观点与应用,在不背离本发明的精神下进行各种修饰与变更。
如图1所示,本发明一种刀具磨损在线监测方法,包括如下步骤:
步骤101,获取一定时间内的刀具加工运行的原始信号数据;在该步骤中,所述原始信号数据包括:1)通过一个测力仪采集加工过程中刀具对工件的力的作用,由于刀具磨损程度的不同,对工件施加的力是变化的,通过采集的力信号的数据,做数据分析,可以得到刀具磨损的程度;2)选用三个加速度传感器,分别安装在工件两个侧面和顶面,测量三个方向的振动信号,用于采集加工过程中工件的振动数据,同样,由于刀具磨损程度不同,产生的振动程度也不同;在本发明中,通过结合这两种传感器采集的数据信息,对刀具的磨损程度做分析判断;
步骤102,对上述采集的数据信息进行滤波降噪和特征提取,得到样本集,将样本集分为训练样本集和测试样本集;
步骤103,选择适当的小波函数对训练样本进行小波分解,分解到多个尺度下;
步骤104,对多个尺度下的训练样本数据分别进行主元分析(PCA)建模并且构建整体主元模型以及相应的SPE值和T2控制限;
步骤105,将测试样本依据训练样本层数分解,将各尺度下的数据加载到相应尺度的主元模型上计算出各尺度下数据的SPE控制限和T2统计量,判断显著尺度,舍弃非显著尺度;步骤106,将显著尺度上的信号利用小波重构方法组成一个新的测试样本;
步骤107,新的测试样本将加载到整体主元模型上,计算SPE控制限和T2统计量并判断是否超出相应控制图的控制限,如果超出则说明过程刀具磨损异常,否则为正常;
所述步骤102具体包括以下处理:
步骤2.1,对力信号进行低通滤波处理,振动信号进行带通滤波处理;
步骤2.2,提取时域特征参数和频域特征参数;
步骤2.3,计算6个时域特征指标,该时域特征指标包括均值、均方根、峰峰值、方差、峭度、偏度;
步骤2.4,计算6个频域特征指标,该频域特征指标包括功率谱和、功率谱均值、功率谱峰值、功率谱方差、功率谱峭度、功率谱偏度;
步骤2.5,训练样本的信息为正常磨损状态下的数据;
步骤2.6,测试样本的信息为一半正常磨损状态下的数据,另一半为异常磨损状态下的数据;
步骤2.7,训练样本和测试样本的数量相同。
在时域范围内,即在时间轴上做特征提取。具体是对每一列信号以一定的长度N提取信号,例如一列信号的长度是m*N长度,我们以长度N来提取信号,那么这一列就可以得到m个特征数,如果我们有6列数,这提取完特征就的到的是m×6的一个矩阵:
(1)均值:计算一段信号的幅值平均值。针对信号X={x1,x2,…,xN}其均值表达式为:
N为信号长度,xi每个信号点对应的值;
(2)均方根:计算一段信号的有效值,其表达式为:
N为信号长度,xi每个信号点对应的值;
(3)峰峰值:计算一段信号中最大值与最小值的差值,其表达式为:
Pk=max(xi)-min(xi)
Xi为每个信号点对应的值;
(4)方差:计算一段信号的方差值,其表达式为:
N为信号长度,Xi为每个信号点对应的值,X为均值;
(5)峭度:计算信号波形陡峭程度的物理量。峭度值的改变反映了系统发生了变化。在机械加工中,峭度值越大,说明加工信号越不正常,发生故障的可能性也就越大。其表达式为:
其中,N为信号长度,Xi为每个信号点对应的值,X为均值,σ为标准差;
(6)偏度:描述信号对称性分布的物理量,反映了信号分布倾斜的方向和程度大小。偏度越大,信号对称性越差。其表达式为:
其中,N为信号长度,Xi为每个信号点对应的值,X为均值,σ为标准差。
信号的频域特征主要反映的是系统在幅值和波动性上的变化,但是却无法获得信号频域上的信息。因此,本发明采用功率谱来提取相关频域特征。功率谱的计算方法主要包括参数估计法和非参数估计法。本文采用的是非参数估计法中的Welch法来计算信号的功率谱。首先进行快速傅立叶变换(FFT数据长度为1024)得到信号的频谱,本文采用加矩形窗的方法,从而有效防止频谱发生泄漏,然后通过计算得到信号的功率谱密度Sx(f)。其中,功率谱密度横轴代表频率,纵轴代表对应频率的功率值。
假设信号X={x1,x2,…,xN},信号的自功率谱Sx(f)和自相关函数Rx(i)之间有下面的表达式:
Rx(i)自相关函数,Ts为信号采样间隔,f为频率;
其中,Ts为信号采样间隔。下面介绍本文所选取的频率特征:
(1)功率谱和:计算信号在一个频段内每个频率所对应的功率值之和。其表达式为:
(2)功率谱均值:计算信号在一定频段内的功率谱均值,其表达式为:
(3)功率谱峰值:计算一定频段内信号功率谱的最大值,其表达式为:
Pp=max(S(f))
(4)功率谱方差:计算一定频段内功率谱的方差,其表达式为:
(5)功率谱峭度:计算一定频段内功率谱的峭度值,其表达式为:
(6)功率谱偏度:计算一定频段内功率谱的偏度值,其表达式为:
所述步骤103包括以下具体处理:
步骤3.1,分解尺度L的选择可以根据L=log2n-5(其中n为样本个数)来确定;
步骤3.2,小波分解后,通过Mallat算法得到一个近似矩阵AL和L个细节矩阵(从D1到DL),Mallat算法本质上相当于用低通滤波器和带通滤波器表示尺度函数和小波函数,分解过程示意如图3所示。其表达式为:
Ai=HAi-1,Di=GAi-1
其中,Di为在尺度i上小波函数系数组成的向量,Ai为在尺度i上尺度函数系数组成的向量,G表示带通滤波器,H表示低通滤波器。图3中以一个三层Mallat算法为例,X为原始数据,G表示带通滤波器,H表示低通滤波器。1/2为二抽取,抽取后的数据采样率和数据点数均变为原来的二分之一。首先数据通过第一次的带通滤波和低通滤波处理,经二抽取后得到第一层的高频系数和低频系数。接着对第一层的低频系数再进行带通滤波和低通滤波处理,从而得到第二层的高频系数和低频系数。依次类推,分解到目标层数为止。AL则为目标层的尺度函数系数组成的向量,L个细节矩阵为每一层产生的小波函数系数组成的向量;
所述步骤104包括以下具体处理:
步骤4.1,获取正常刀具加工过中的信号数据;
步骤4.2,数据样本X′进行标准化处理;
步骤4.3,将数据X进行协方差矩阵分解,数据X用下式表示:X=SPT;
步骤4.4,计算各主元贡献率(CPV),将前几个主元的贡献率累加,达到一定阈值后即可确定主元个数,其阈值一般设定在80%;
步骤4.5,计算SPE统计量以及控制限;
步骤4.6,计算T2统计量以及控制限。
以下将通过具体实施例说明本发明的步骤:
步骤1,获取原始信号数据。通过测力仪和加速度传感器分别测量三个方向切削力信号和三个方向振动信号,要求刀具正常磨损的信号数据是异常磨损得信号数据的3倍;
步骤2,对数据进行滤波降噪和特征提取,得到样本集,将样本集分为训练样本集和测试样本集。对切削力信号做低通滤波处理,振动信号做带通滤波处理。对完成滤波的信号提取时域和频域特征。提取6个时域特征指标包括均值、均方根、峰峰值、方差、峭度、偏度。提取6个频域特征指标包括功率谱和、功率谱均值、功率谱峰值、功率谱方差、功率谱峭度、功率谱偏度。取正常磨损数据的三分之二为训练样本,正常磨损数据剩下的三分之一和异常磨损信号组成测试样本。此时训练样本中的数据全部为刀具正常磨损信号数据,训练样本的数据前一半为正常磨损信号数据,后一半为异常磨损信号数据,且两组数据量相同;
步骤3,选择适当的小波函数对训练样本进行小波分解,分解到多个尺度下。通过多分辨率小波分析(Mallat算法)分解训练样本,分解层数L=log2n-5(其中n为样本个数)。分解后得到一个近似矩阵AL和L个细节矩阵(从D1到DL)(L已经表示过分解层数,因此,分解L层则得到L个细节矩阵(如图3所示,分解三层则得到三个细节矩阵,D1、D2、D3));
步骤4,对多个尺度下的训练样本数据分别进行主元分析(PCA)建模并且构建整体主元模型以及相应的SPE和T2控制限。利用主元分析对训练样本计算各个主元的贡献率(CPV),累加前几个主元的贡献率达到80%,则认为这些主元信息包括大部分的训练样本的信息,从而大幅度降低计算难度。经降维处理后,各尺度所选择的主元向量被用来构建主元模型。主元模型主要包括负载矩阵和得分矩阵。负载矩阵主要用来将新的测量数据投影到各主元方向,得分矩阵则是用来计算SPE值和T2统计量及各尺度下的控制限。同时,为了实现最后的结果监测,仅仅依靠各尺度下的控制限是不够的。我们还需要针对刀具正常磨损状态下的整体数据进行主元分析,构建整体主元模型以及相应的SPE值和T2控制限;
步骤5,将测试样本依据训练样本层数分解,将各尺度下的数据加载到相应尺度的主元模型上计算出各尺度下数据的SPE和T2统计量,判断显著尺度,舍弃非显著尺度。训练样本模型构建过程,得到了L个尺度下的主元模型和一个整体主元模型以及相对应的控制限,首先将测试数据进行小波分解,这里分解层数与建模过程相同。然后,将各尺度下的数据加载到相应尺度的主元模型上计算出各尺度下数据的SPE和T2统计量,判断哪些尺度为非显著尺度,在小波重构时舍弃非显著尺度;
步骤6,通过步骤5,选择出测试样本中的显著尺度。将显著尺度上的信号利用小波重构方法组成一个新的测试样本;
步骤7,新的测试样本将加载到步骤4的整体主元模型上,计算SPE值和T2统计量并判断是否超出相应控制图的控制限,如果超出则说明过程刀具磨损异常,否则为刀具磨损正常,从而监测处刀具的磨损状态。
综上所述,本发明一种刀具磨损在线监测方法是根据加工过程刀具磨损的特点,即存在多尺度、非线性、非平稳信号,特别是难加工材料加工过程中,刀具极易发生磨损,我们很难获得所有刀具磨损状态下的数据,多尺度主元分析(MSPCA)结合了主元分析(PCA)和小波分析,前者具有强大的数据降维以及提取变量间线性关系的能力,可以降低变量间的线性相关性、减小模型计算量,后者可以提取变量多尺度特征和变量自相关性,克服了主元分析只适用于单尺度的局限,提高了刀具磨损在线监测的准确率,对于刀具磨损状态在线监测具有重要意义。
上述实例仅示例性说明本发明的原理及其功效,而非用于限制本发明。任何本发明领域技术人员均可在不违背本发明的精神及范畴下,对上述实例进行修饰与改变。因此,本发明的权利保护范围,应如权利要求书所列。
Claims (3)
1.一种基于多尺度主元分析的刀具磨损在线监测方法,其特征在于,该方法包括以下步骤:
步骤一,获取在线监测时间内,刀具加工运行的原始信号数据;包括加工过程中刀具对工件两个侧面和顶面的三个方向的切削力信号数据和振动信号数据,将所述原始信号数据构成样本集;
步骤二,通过对上述采集的原始信号数据进行滤波降噪、提取时域和频域12个特征参数,将样本集分为训练样本集和测试样本集;
步骤三,通过小波分解将训练样本集分解成多个尺度,小波分解后,得到一个近似矩阵AL和L个细节矩阵,AL为目标层的尺度函数系数组成的向量,L个细节矩阵为每一层产生的小波函数系数组成的向量;
步骤四,利用主元分析对多个尺度下的所述训练样本集的数据分别进行建模,构建整体主元模型以及计算相应的SPE值及控制限和T2统计量及控制限;
步骤五,将测试样本依据训练样本的层数进行小波分解,将各分解尺度下的数据加载到相应尺度的主元模型上计算出各尺度下数据的SPE控制限和T2统计量:
在i时刻,SPE控制限表示为:
其中,Xi为i时刻测得的数据,I为单位矩阵,Pk为主元模型的负荷矩阵;
在i时刻,T2统计量以及控制限表示为:
其中,O为前k个主元的特征值组成的对角矩阵,Xi为i时刻测得的数据,Pk为主元模型的负荷矩阵;
判断显著尺度,舍弃非显著尺度;
步骤六,将显著尺度上的信号利用小波重构方法组成一个新的测试样本;
步骤七,所述新的测试样本将加载到整体主元模型上,计算SPE控制限和T2统计量,并判断是否超出相应控制图的控制限
SPE控制限表示为:
其中,λ为前k个主元的特征值组成的对角矩阵,k是主元个数,m是变量个数,Cα是正态分布在假设检验水平为α下的临界值;
T2控制限表示为:
其中,k为主元个数,m为变量个数,Fk,m-k,α为对应于检验水平为α、自由度为k、n-1个样本条件下的F分布临界值;
如果超出,则说明过程刀具磨损异常;否则,说明过程刀具磨损正常。
2.如权利要求1所述的一种基于多尺度主元分析的刀具磨损在线监测方法,其特征在于,所述步骤二中的12个特征参数包括6个时域特征指标即包括均值、均方根、峰峰值、方差、峭度、偏度和6个频域特征指标即功率谱和、功率谱均值、功率谱峰值、功率谱方差、功率谱峭度、功率谱偏度。
3.如权利要求1所述的一种基于多尺度主元分析的刀具磨损在线监测方法,其特征在于,所述步骤四的利用主元分析对多个尺度下的所述训练样本集的数据分别进行建模的步骤,具体包括以下处理:
获取正常刀具加工过中的信号数据;
对训练样本集X′进行标准化处理;
将训练样本集中数据X进行协方差矩阵分解,数据X用下式表示:
X=SPT
其中:S=[s1,s2,s3,…,sm]称为得分矩阵,顺序按向量的模的大小排列,提取的是前后数据间的相互关系;P=[p1,p2,p3,…,pm]称为负荷矩阵,反映了变量间的相互关系;
计算各主元贡献率,依序将主元的贡献率累加,达到阈值后,确定主元个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610269068.9A CN105834835B (zh) | 2016-04-26 | 2016-04-26 | 一种基于多尺度主元分析的刀具磨损在线监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610269068.9A CN105834835B (zh) | 2016-04-26 | 2016-04-26 | 一种基于多尺度主元分析的刀具磨损在线监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105834835A CN105834835A (zh) | 2016-08-10 |
CN105834835B true CN105834835B (zh) | 2018-06-19 |
Family
ID=56590232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610269068.9A Active CN105834835B (zh) | 2016-04-26 | 2016-04-26 | 一种基于多尺度主元分析的刀具磨损在线监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105834835B (zh) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107263211B (zh) * | 2017-05-26 | 2018-10-30 | 温州大学 | 一种基于多传感融合的刀具状态监测方法 |
CN107253094B (zh) * | 2017-05-26 | 2018-12-28 | 温州大学 | 一种基于稳态子空间分析多维信号的刀具状态监测方法 |
CN107588947B (zh) * | 2017-08-23 | 2020-04-28 | 电子科技大学 | 考虑运行工况和信息简化的风电机组传动链状态监测方法 |
JP2019072806A (ja) * | 2017-10-17 | 2019-05-16 | オムロン株式会社 | 切削加工装置 |
TWI650625B (zh) | 2017-11-16 | 2019-02-11 | 財團法人工業技術研究院 | 刀具磨耗檢測裝置、其檢測方法及刀具磨耗補償方法 |
CN108520125A (zh) * | 2018-03-29 | 2018-09-11 | 上海理工大学 | 一种预测刀具磨损状态的方法及系统 |
CN109333159B (zh) * | 2018-09-11 | 2021-04-13 | 温州大学苍南研究院 | 刀具磨损状态在线监测的深度核极限学习机方法及系统 |
TWI665051B (zh) * | 2018-12-10 | 2019-07-11 | National Chin-Yi University Of Technology | 用於工具機之刀具磨損之檢測方法 |
CN109605127A (zh) * | 2019-01-21 | 2019-04-12 | 南京航空航天大学 | 一种刀具磨损状态识别方法及系统 |
CN110153799B (zh) * | 2019-05-14 | 2021-02-09 | 华中科技大学 | 一种基于永磁扰动探头的铣刀破损检测方法、装置及应用 |
CN110288046B (zh) * | 2019-07-02 | 2022-11-18 | 南京恩瑞特实业有限公司 | 一种基于小波神经网络与隐马尔科夫模型的故障预测方法 |
CN110610484B (zh) * | 2019-08-21 | 2022-04-12 | 西安理工大学 | 一种基于旋转投影变换的印刷网点质量检测方法 |
CN110653661A (zh) * | 2019-09-30 | 2020-01-07 | 山东大学 | 基于信号融合和多重分形谱算法的刀具状态监测识别方法 |
CN110856437B (zh) * | 2019-10-09 | 2021-06-15 | 浙江工业大学 | 一种smt生产过程控制图模式识别方法 |
CN110850230B (zh) * | 2019-11-07 | 2022-04-29 | 天津大学 | 一种基于多小波理论的直流配电网极间短路故障特征量提取方法 |
CN111008356B (zh) * | 2019-11-13 | 2023-06-16 | 成都理工大学 | 一种基于WTSVD算法扣除背景的γ能谱集分析方法 |
CN112008495B (zh) * | 2020-07-28 | 2021-09-07 | 成都飞机工业(集团)有限责任公司 | 一种基于振动监测的刀具破损识别方法 |
CN113344395B (zh) * | 2021-06-14 | 2022-06-21 | 西北工业大学 | 一种基于动态pca-svm的加工质量监测方法 |
CN113458873B (zh) * | 2021-07-01 | 2022-03-11 | 太原科技大学 | 一种刀具磨损量和剩余寿命预测的方法 |
CN113601261B (zh) * | 2021-08-10 | 2022-06-14 | 中国科学院合肥物质科学研究院 | 一种面向刀具的在线快速优化模型的监测方法 |
CN114273978B (zh) * | 2021-12-27 | 2023-03-21 | 深圳市玄羽科技有限公司 | 刀具的磨损检测方法、电子设备及计算机可读存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102689230A (zh) * | 2012-05-09 | 2012-09-26 | 天津大学 | 基于条件随机场模型的刀具磨损状态监测方法 |
CN103456107A (zh) * | 2013-08-19 | 2013-12-18 | 福建升腾资讯有限公司 | 一种快速检测pos支付终端的方法 |
CN103962888A (zh) * | 2014-05-12 | 2014-08-06 | 西北工业大学 | 一种基于小波去噪和希尔伯特-黄变换的刀具磨损监测方法 |
CN104708497A (zh) * | 2015-03-17 | 2015-06-17 | 洛阳理工学院 | 一种基于电流与声发射复合信号的刀具磨损监测系统 |
CN105196114A (zh) * | 2015-11-05 | 2015-12-30 | 西安科技大学 | 基于小波分析和神经网络的刀具磨损实时在线监测方法 |
CN105397566A (zh) * | 2015-12-22 | 2016-03-16 | 大连民族大学 | 一种加工中心刀具磨损在线检测与智能补偿系统及方法 |
CN105479269A (zh) * | 2015-12-28 | 2016-04-13 | 东莞市巨冈机械工业有限公司 | 一种数控机床的刀具磨损监控方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20050030925A (ko) * | 2005-03-03 | 2005-03-31 | 김용연 | 가공공정에서 표면결함 및 공구상태에 대한 음향신호 분석시스템 |
-
2016
- 2016-04-26 CN CN201610269068.9A patent/CN105834835B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102689230A (zh) * | 2012-05-09 | 2012-09-26 | 天津大学 | 基于条件随机场模型的刀具磨损状态监测方法 |
CN103456107A (zh) * | 2013-08-19 | 2013-12-18 | 福建升腾资讯有限公司 | 一种快速检测pos支付终端的方法 |
CN103962888A (zh) * | 2014-05-12 | 2014-08-06 | 西北工业大学 | 一种基于小波去噪和希尔伯特-黄变换的刀具磨损监测方法 |
CN104708497A (zh) * | 2015-03-17 | 2015-06-17 | 洛阳理工学院 | 一种基于电流与声发射复合信号的刀具磨损监测系统 |
CN105196114A (zh) * | 2015-11-05 | 2015-12-30 | 西安科技大学 | 基于小波分析和神经网络的刀具磨损实时在线监测方法 |
CN105397566A (zh) * | 2015-12-22 | 2016-03-16 | 大连民族大学 | 一种加工中心刀具磨损在线检测与智能补偿系统及方法 |
CN105479269A (zh) * | 2015-12-28 | 2016-04-13 | 东莞市巨冈机械工业有限公司 | 一种数控机床的刀具磨损监控方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105834835A (zh) | 2016-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105834835B (zh) | 一种基于多尺度主元分析的刀具磨损在线监测方法 | |
EP3136255A1 (en) | Method and system for automatically detecting defects on a rotating shaft | |
US20140149325A1 (en) | System monitor and method of system monitoring | |
KR101891175B1 (ko) | 동특성 변화를 이용한 건물 구조물의 안전도 평가방법 및 시스템 | |
CN103310113B (zh) | 一种基于频带分离和数据建模的通用血糖预测方法 | |
JP6792746B2 (ja) | 周波数領域のセグメント特徴解析による状態識別法 | |
CN111354366A (zh) | 一种异常声音检测方法及异常声音检测装置 | |
CN116237817B (zh) | 基于物联网的五轴联动数控机床智能监测系统 | |
CN109238455B (zh) | 一种基于图论的旋转机械振动信号监测方法及系统 | |
CN111170103B (zh) | 设备故障识别方法 | |
CN111678699B (zh) | 一种面向滚动轴承早期故障监测与诊断方法及系统 | |
CN107748209B (zh) | 一种结构损伤检测方法 | |
CN112393906B (zh) | 一种地铁车辆转向架轴承微弱信号故障诊断、分类和健康评估的方法 | |
CN112720071B (zh) | 多能域信号智能化融合的刀具实时状态监测指标构造方法 | |
CN115741235A (zh) | 基于五轴加工中心刀具的磨损预测与健康管理方法 | |
CN115127806A (zh) | 一种基于多传感器振动信号的齿轮箱故障诊断方法及装置 | |
Medina et al. | Deep learning-based gear pitting severity assessment using acoustic emission, vibration and currents signals | |
CN114492527A (zh) | 基于模糊神经网络与主成分分析表面粗糙度在线预测方法 | |
CN110108474A (zh) | 一种旋转机械运行稳定性在线监测与评估方法及系统 | |
CN108615018A (zh) | 基于时域直方图特征提取的对象物状态识别方法 | |
Al-Bugharbee et al. | Autoregressive modelling for rolling element bearing fault diagnosis | |
CN108647604A (zh) | 基于频域直方图特征提取的对象物状态识别方法 | |
KR101724151B1 (ko) | 앰비규어티 패턴인식을 이용한 금속 이물질의 충격 위치와 질량 동시추정 방법 | |
Yanez-Borjas et al. | Methodology based on statistical features and linear discriminant analysis for damage detection in a truss-type bridge | |
Deuschle et al. | Robust sensor spike detection method based on dynamic time warping |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP02 | Change in the address of a patent holder |
Address after: 300350 Haijing garden, Haihe Education Park, Jinnan, Tianjin, 135, Tianjin University. Patentee after: Tianjin University Address before: 300072 Tianjin City, Nankai District Wei Jin Road No. 92 Patentee before: Tianjin University |
|
CP02 | Change in the address of a patent holder |