CN112749506A - 基于ceemd和gwo-svr的振动信号预测方法 - Google Patents

基于ceemd和gwo-svr的振动信号预测方法 Download PDF

Info

Publication number
CN112749506A
CN112749506A CN202011517791.7A CN202011517791A CN112749506A CN 112749506 A CN112749506 A CN 112749506A CN 202011517791 A CN202011517791 A CN 202011517791A CN 112749506 A CN112749506 A CN 112749506A
Authority
CN
China
Prior art keywords
wolf
prediction
svr
gwo
prey
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
CN202011517791.7A
Other languages
English (en)
Other versions
CN112749506B (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.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN202011517791.7A priority Critical patent/CN112749506B/zh
Publication of CN112749506A publication Critical patent/CN112749506A/zh
Application granted granted Critical
Publication of CN112749506B publication Critical patent/CN112749506B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • G06N20/10Machine learning using kernel methods, e.g. support vector machines [SVM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Optimization (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Molecular Biology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

汽车覆盖件模具多采用镶块式模件拼接后整体加工,拼接区加工时易引发载荷突变产生冲击振动,影响拼接区的整体加工质量,为了提高拼接区的加工精度,对铣削过程的时域振动信号进行前瞻预测。首先基于互补式集合经验模态分解方法将铣削振动信号进行6层模态分解,得到各层本征模态函数及趋势序列;然后分别构建不同工况下的支持向量回归预测模型,采用灰狼优化算法对支持向量回归中的参数进行寻优分析;最后对时域振动信号进行重构和前瞻预测。试验结果表明,在淬硬钢拼接区铣削过程中,结合CEEMD和GWO‑SVR的铣削振动信号前瞻预测方法相较于其它传统方法具有更良好的预测效果,在预测时间为0.12秒时总体预测准确率达94%以上。

Description

基于CEEMD和GWO-SVR的振动信号预测方法
技术领域
本发明属于振动信号预测领域,特别是涉及淬硬钢拼接模具铣削过程的智能预测方法。
背景技术
汽车车身构件日益复杂,相应淬硬钢模具存在大量的凸凹、沟槽等结构,为克服模具加工制造和装配中存在的这种结构性困难,多采用镶块式拼接模件。在模具整体铣削加工过程中,由于在拼接处存在硬度差,铣削振动变化大,产生的冲击振动会在短期内导致刀具磨损加剧,模具表面加工质量下降等。如果能准确地预测模具拼接区的铣削振动趋势,通过快速改变铣削参数进行前瞻控制,有利于在拼接处减小铣削力的瞬时突变,减少冲击振动对模具加工精度和产品质量的影响。
在拼接模具整体铣削加工过程中,由于加工工况和切削环境复杂,测得的铣削振动信号中存在干扰信息和较强的非平稳特征,这些干扰信息的存在会严重影响铣削振动信号预测的精准度,因此需要对振动数据进行预处理以提高振动信号前瞻预测的准确性。
目前时频分析方法是分析非平稳和非线性信号的最常用的方法,如小波分析法、Hilbert-Huang变换等。小波变换应用过程具有局限性,其基函数和分解层数的选取缺乏自适应性。经验模态分解(EMD)方法是一种瞬时信号处理方法,避免了选取小波基函数的弊端,能够实现非平稳信号的分解并滤除噪声干扰,EMD方法以Hilbert变换为基础,将瞬时频域具有物理意义的信号定义为固有模态函数(IMF),在非平稳信号处理领域得到应用广泛。由于EMD分解方法存在模态混叠现象,因此Wu和Huang提出了集合经验模态分方法(EEMD),有效地抑制了信号分解过程中局部极值在短时间内的频繁跳动。
铣削振动信号预测常采用线性时间序列的预测方法,但在实际切削环境中,机床、刀具和工件使铣削呈现典型的非线性特征,因此铣削振动信号的非线性时间序列预测方法得到了迅速发展。对于非线性时间序列预测,人工神经网络(ANN)的特殊性在于其本身处理数据时具有自学习能力,而且对于非线性数据的预测能力强,然而ANN算法在操作中会出现一些问题,如收敛速度慢、难以确定网络结构和局部最优等。支持向量机是一种基于统计学习理论的学习算法,它的优势主要是在小样本、少信息、非线性及高维空间模式识别中,克服了人工神经网络极易陷于局部极小和过度拟合的问题,而且对于小样本数据分析来说,支持向量机拥有良好的学习能力和推广能力,被广泛应用于铣削颤振识别、粗糙度预测和铣削力预测中,已经成为研究者用于解决非线性时间序列预测的一种重要工具。
发明内容
本发明要解决的技术问题是提供基于CEEMD和GWO-SVR的振动信号预测模型,一方面通过加入正负白噪声的方式减少模态混叠对于振动信号EMD分解的影响,提高模态分量的平稳性;另一方面基于修正的GWO算法寻找支持向量回归机参数的最优解,总体上提高拼接处铣削振动信号预测的精准度。
为实现上述目的,本发明采用以下技术方案:
基于CEEMD和GWO-SVR的振动信号预测方法,所述方法的实现过程为:
步骤一、基于互补集合经验模态分解(CEEMD)将实时采集的时序振动信号进行6层分解,得到的各层本征模态函数及趋势序列;
步骤二、对得到的各层本征模态函数及趋势序列构建训练集和预测集,基于构建训练集和预测集对利用GWO优化支持向量回归的参数C和g进行寻优,建立GWO-SVR预测模型:
建立SVR预测模型:
为了使实际风险达到最小的效果,根据结构风险最小化原理,优化的结构风险目标函数为:
Figure BDA0002847988580000021
式中,||w||2为描述函数;C为经验风险和模型复杂度之间取一折中的常数;|yi-f(xi)|为ε不敏感损失函数,即:
Figure BDA0002847988580000031
SVR即为(4)式的优化问题的求解公式:
Figure BDA0002847988580000032
式中,ξi
Figure BDA0002847988580000033
为引入的松弛变量,目的是使(4)式有解。为了求解上述的优化问题,通过Lagrange函数,引入Lagrange乘子a、a*,将该问题转化为其对偶问题:
Figure BDA0002847988580000034
式中,当
Figure BDA0002847988580000035
非零时对应的训练样本为支持向量。求解此二次规划问题可求出a的值,同时求得w的值:
Figure BDA0002847988580000036
利用KKT(Karush-Kuhn-Tucker)条件计算出偏差b:
Figure BDA0002847988580000037
最后得到回归函数f(x)的表达式:
Figure BDA0002847988580000038
式中,K(xi,x)=Φ(xi)Φ(x)为一个满足Mercer条件的核函数。此函数可以越过具体形式从而达到非线性变换的非线性化操作,这是SVR的一个显著特点。目前SVR核函数中使用最多的是带有宽度为σ的径向基核函数(RBF核函数),使用径向基核函数的SVR预测模型的预测效果要优于使用其它核函数预测模型,径向基核函数为:
构建GWO的过程:
(1)追踪猎物
灰狼在狩猎过程中围绕猎物,将其行为做出以下定义:
Figure BDA0002847988580000041
Figure BDA0002847988580000042
式中,
Figure BDA0002847988580000043
代表灰狼经过第t次迭代后所得出的猎物当前所在位置,
Figure BDA0002847988580000044
代表灰狼个体经过第t次迭代后当前所在位置,这就是算法的当前局部最优解。
Figure BDA0002847988580000045
Figure BDA0002847988580000046
式中,
Figure BDA0002847988580000047
Figure BDA0002847988580000048
为包围猎物中的随机系数,
Figure BDA0002847988580000049
的区间为[0,2]且随着迭代次数的增加而降低,
Figure BDA00028479885800000410
Figure BDA00028479885800000411
是区间[0,2]上的随机数。
(2)包围猎物
灰狼有识别猎物的强能力,且α狼可以领导灰狼包围猎物,但在猎物(最优解)位置未知的情况下。为了模拟灰狼捕猎行为,假定α狼、β狼和δ狼了解猎物的潜在位置,保存三个最佳解决方案,并且要求其他灰狼(ω狼)作为搜索代理,根据最佳解决方案进行灰狼位置的迭代更新,公式为:
Figure BDA00028479885800000412
式中,
Figure BDA00028479885800000413
Figure BDA00028479885800000414
分别表示α狼、β狼和δ狼与ω狼之间的距离;
Figure BDA00028479885800000415
Figure BDA00028479885800000416
分别表示第t次迭代后α狼、β狼和δ狼的当前位置;
Figure BDA00028479885800000417
表示当前灰狼位置。
Figure BDA00028479885800000418
Figure BDA00028479885800000419
公式(17)代表ω狼朝三个潜在解α狼、β狼和δ狼的步长和方向;公式(14)代表ω狼的最终位置。
(3)攻击猎物
当猎物停止移动时狩猎完成,这一过程可以通过降低
Figure BDA0002847988580000051
来完成。当随机值
Figure BDA0002847988580000052
在[-1,1]时,搜索代理的下一位置可以在当前位置与猎物位置之间的位置进行更新。算法中,当
Figure BDA0002847988580000053
时强制狩猎。
(4)寻找猎物
灰狼常常分头搜寻猎物,同时根据α狼、β狼和δ狼的位置来攻击猎物,当
Figure BDA0002847988580000054
时,灰狼被迫离开猎物,
Figure BDA0002847988580000055
灰狼强制狼群攻击猎物,直到满足条件时, GWO算法停止。
建立GWO-SVR模型,其过程为:
(1)基于GWO对SVR中的参数C、g进行参数寻优过程中,首先进行初始化:种群规模N=20,最大迭代次数为200;确定SVR中参数C、g的范围,即C∈[0.01,100],g∈[0.01,100];
(2)将所述训练样本对输入到SVR模型中,对模型进行训练,将均方根误差RMSE作为适应度函数Fitness;
(3)找到适应度函数Fitness的最小值对应的个体,开始迭代寻优,并判断最小Fitness是否低于前一代最小Fitness;如低于,则保留最小Fitness值及其对应的坐标,并将其赋给初始坐标;如高于,则返回步骤(2);
(4)找到C、g的最佳值,建立GWO-SVR预测模型;
步骤三、采用滚动预测模型,铣削振动信号中的第1~100时序数据点预测第101点,用第2~101时序数据点预测第102点,输入记为X,输出记为 Y,一个X和Y记为一组样本,总样本数的60%构建为训练集,总样本数的40%构建为测试集;
步骤四、对经过CEEMD算法分解得到的加速度信号(IMFs信号)通过 GWO-SVR预测模型进行预测,并将预测曲线进行反归一化处理;然后将各个经过反归一化处理后的预测曲线进行预测信号结果重构,得到振动加速度总信号(IMFs信号)的预测曲线;
步骤五、选取两个训练集和测试集比例,一个是训练集为总体数据量的 60%,测试集为总体数据量的40%,另一个是训练集为总体数据量的20%,测试集为总体数据量的80%。当训练集时序数据减少,预测时序数据增加1倍时,观察所预测数据的精准度变化情况。
根据权利要求1所属基于CEEMD和GWO-SVR的振动信号预测方法,其特征在于:
在步骤一中,通过互补集合经验模态分解(CEEMD)将实时采集的时序振动信号进行6层分解,得到的各层本征模态函数及趋势序列:
步骤一一、获得拼接模具铣削振动信号;
步骤一二、基于EMD模型加入互补白噪声,建立CEEMD模型:
CEEMD算法具体步骤如下:
(1)将一组符号相反的噪声加入到目标信号中,每组噪声的幅值相同,然后产生一对新信号:
s+ i(t)=x(t)+ni +(t) (15)
s- i(t)=x(t)+ni -(t) (16)
式中,si(t)是新信号,ni +(t)是正噪声信号,ni -(t)是负噪声信号。
(2)分解新信号,得到两组含有白噪声的信号分量IMF1和IMF2,且这两组分量正负互补。
(3)将两组信号的IMF分量进行整合得到最终的IMF分量。
步骤一三、通过CEEMD分解振动信号获得IMFs。
本发明的有益效果是:
本发明提出了一种基于互补式集成经验模态分解结合改进灰狼群优化算法-支持向量机回归的铣削振动信号前瞻预测方法,展现了了GWO算法强大的计算力,发挥出GWO-SVR预测模型的预测性能。基于互补式集合经验模态分解方法将铣削振动信号进行6层模态分解,得到各层本征模态函数及趋势序列;然后分别构建不同工况下的支持向量回归预测模型,采用灰狼优化算法对支持向量回归中的参数进行寻优分析;最后对时域振动信号进行重构和前瞻预测。通过实际加工与理论对比验证精度,能在加工前进行参数优化,对加工提供指导
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所述的模具拼接区铣削振动信号示意图;
图2为本发明所述的CEEMD分解原理流程图;
图3为本发明所述的信号分解得到的IMFs图;
图4为本发明所述的GWO-SVR算法流程图。
图5为本发明所述的寻优算法参数迭代图;
图6为本发明所述的振动信号预测图;
图7为本发明所述的不同预测模型对比下的信号预测图;
图8为本发明所述的预测时长为0.12s的信号预测图;
图9为本发明所述的预测时长为0.24s的信号预测图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地阐述。
具体实施方式一:如图1-9所示,本实施方式针对所述的一种基于 CEEMD和GWO-SVR的铣削振动信号前瞻预测方法的实现过程如下:
1.模具拼接区铣削信号采集
实验用的机床为VDL-1000E型三轴立式加工中心,刀具为戴杰二刃整体硬质合金球头立铣刀(DV-OCSB2100-L140);工件材料为淬硬状态下模具钢 (Cr12MoV),工件尺寸为200mm×200mm×60mm,分成3个镶块,硬度分别通过不同的热处理工艺调至45HRC、50HRC、45HRC,通过2根长度为 200mm的不锈钢高强度铰制孔用螺栓连接为一个整体,拼接缝隙为0.2mm。曲面淬硬钢模具铣削采用小切深,小行距和小进给的切削参数,加工参数如表 1所示。
表1实验加工参数和刀具参数
Figure BDA0002847988580000071
Figure BDA0002847988580000081
铣削实验时,从硬度为45HRC的镶块向硬度为50HRC的镶块进行铣削加工,铣削采用顺铣切削,共计完成9组实验。每组实验时均通过加速度传感器、数据采集系统测试工件的X、Y、Z三个分量上的振动加速度信号,采集到的铣削振动信号为动态铣削力和过缝冲击力共同激励后产生;由于Z向系统刚度较大,铣削振动信号幅值较小;X、Y两个方向的铣削振动信号幅值较大,进给方向(机床Y方向)最大,因此本文分析铣削振动信号波形时只提取X、Y 两个方向的实验数据对其进行分析,过实验模具拼接缝后的铣削振动信号如图 1所示。如图1所示的铣削振动信号为主轴转速n=4000r/min、每齿进给量 fz=0.15mm/z、切深分别为ae=0.1mm、ae=0.2mm、ae=0.3mm和ae=0.4mm 时的四组铣削振动加速度实验数据,数据时序为铣削实验过程中的 39.43~39.73s
2.基于CEEMD的铣削振动信号分解
过拼接缝后的铣削振动信号中存在有干扰信息,同时在铣削加工过程中的振动信号具有较强的非平稳特征,这些干扰信息对准确预测拼接缝振动信号会产生较大影响,所以需要对数据进行预处理。
EMD算法是一种自适应式时频域信号分解方法,相对于传统信号分解模型, EMD展现出了适应强的特性,EMD将一组信号分解成了若干个固有模态分量和一个残余项,公式如下:
Figure BDA0002847988580000091
式中,x(t)是原始信号;m是IMF个数;ci(t)是第i个IMF分量;rm(t)是残余项。
但是EMD算法存在模态混叠问题,针对模态混叠问题,出现了EEMD算法来解决模态混叠问题,EEMD算法是在EMD分解的每一步骤中分别加入等幅的白噪声,但加入的白噪声会产生信息干扰,Ye等人提出了CEEMD算法,CEEMD 算法在信号中加入了N对相反噪声,然后分解加入噪声的新信号,最后得到不同的IMF分量,从而进行分解。CEEMD算法具体步骤如下:
(1)将一组符号相反的噪声加入到目标信号中,每组噪声的幅值相同,然后产生一对新信号:
s+ i(t)=x(t)+ni +(t) (18)
s- i(t)=x(t)+ni -(t) (19)
式中,si(t)是新信号,ni +(t)是正噪声信号,ni -(t)是负噪声信号。
(2)分解新信号,得到两组含有白噪声的信号分量IMF1和IMF2,且这两组分量正负互补。
(3)将两组信号的IMF分量进行整合得到最终的IMF分量。
CEEMD算法流程图如图2所示。
对加速度信号进行CEEMD分解之前,为了避免实验数据出现饱和现象,对实验数据进行归一化处理,使其样本数值处于[0,1]之间。
实验中,对铣削振动加速度信号进行CEEMD分解,得到5个频域稳定的固有模态分量IMF1-IMF5和一个剩余分量Res的波形图,不同工况下的CEEMD分解如图3所示。
3.支持向量机
为了使实际风险达到最小的效果,根据结构风险最小化原理,优化的结构风险目标函数为:
Figure BDA0002847988580000092
式中,||w||为描述函数;C为经验风险和模型复杂度之间取一折中的常数; |yi-f(xi)|为ε不敏感损失函数,即:
Figure BDA0002847988580000101
SVR即为(4)式的优化问题的求解公式:
Figure BDA0002847988580000102
式中,ξi
Figure BDA0002847988580000103
为引入的松弛变量,目的是使(4)式有解。为了求解上述的优化问题,通过Lagrange函数,引入Lagrange乘子a、a*,将该问题转化为其对偶问题:
Figure BDA0002847988580000104
式中,当
Figure BDA0002847988580000105
非零时对应的训练样本为支持向量。求解此二次规划问题可求出a的值,同时求得w的值:
Figure BDA0002847988580000106
利用KKT(Karush-Kuhn-Tucker)条件计算出偏差b:
Figure BDA0002847988580000107
最后得到回归函数f(x)的表达式:
Figure BDA0002847988580000108
式中,K(xi,x)=Φ(xi)Φ(x)为一个满足Mercer条件的核函数。此函数可以越过具体形式从而达到非线性变换的非线性化操作,这是SVR的一个显著特点。目前SVR核函数中使用最多的是带有宽度为σ的径向基核函数(RBF核函数),使用径向基核函数的SVR预测模型的预测效果要优于使用其它核函数预测模型[21],径向基核函数为:
Figure BDA0002847988580000111
4.改进的GWO-SVR预测模型
4.1灰狼优化算法
灰狼群狩猎过程分为以下几个阶段:
(1)追踪猎物
灰狼在狩猎过程中围绕猎物,将其行为做出以下定义:
Figure BDA0002847988580000112
Figure BDA0002847988580000113
式中,
Figure BDA0002847988580000114
代表灰狼经过第t次迭代后所得出的猎物当前所在位置,
Figure BDA0002847988580000115
代表灰狼个体经过第t次迭代后当前所在位置,这就是算法的当前局部最优解。
Figure BDA0002847988580000116
Figure BDA0002847988580000117
式中,
Figure BDA0002847988580000118
Figure BDA0002847988580000119
为包围猎物中的随机系数,
Figure BDA00028479885800001110
的区间为[0,2]且随着迭代次数的增加而降低,
Figure BDA00028479885800001111
Figure BDA00028479885800001112
是区间[0,2]上的随机数。
(2)包围猎物
灰狼有识别猎物的强能力,且α狼可以领导灰狼包围猎物,但在猎物(最优解)位置未知的情况下。为了模拟灰狼捕猎行为,假定α狼、β狼和δ狼了解猎物的潜在位置,保存三个最佳解决方案,并且要求其他灰狼(ω狼)作为搜索代理,根据最佳解决方案进行灰狼位置的迭代更新,公式为:
Figure BDA00028479885800001113
式中,
Figure BDA00028479885800001114
Figure BDA00028479885800001115
分别表示α狼、β狼和δ狼与ω狼之间的距离;
Figure BDA00028479885800001116
Figure BDA00028479885800001117
分别表示第t次迭代后α狼、β狼和δ狼的当前位置;
Figure BDA00028479885800001118
表示当前灰狼位置。
Figure BDA0002847988580000121
Figure BDA0002847988580000122
公式(17)代表ω狼朝三个潜在解α狼、β狼和δ狼的步长和方向;公式(18)代表ω狼的最终位置。
(3)攻击猎物
当猎物停止移动时狩猎完成,这一过程可以通过降低
Figure BDA0002847988580000123
来完成。当随机值
Figure BDA0002847988580000124
在[-1,1]时,搜索代理的下一位置可以在当前位置与猎物位置之间的位置进行更新。算法中,当
Figure BDA0002847988580000125
时强制狩猎。
(4)寻找猎物
灰狼常常分头搜寻猎物,同时根据α狼、β狼和δ狼的位置来攻击猎物,当
Figure BDA0002847988580000126
时,灰狼被迫离开猎物,
Figure BDA0002847988580000127
灰狼强制狼群攻击猎物,直到满足条件时, GWO算法停止。
4.2GWO-SVR模型
建立GWO-SVR模型,其过程为:
(1)基于GWO对SVR中的参数c、g进行参数寻优过程中,首先进行初始化:种群规模N=20,最大迭代次数为200;确定SVR中参数C、g的范围,即C∈[0.01,100],g∈[0.01,100];
(2)将所述训练样本对输入到SVR模型中,对模型进行训练,将均方根误差RMSE作为适应度函数Fitness;
(3)找到适应度函数Fitness的最小值对应的个体,开始迭代寻优,并判断最小Fitness是否低于前一代最小Fitness;如低于,则保留最小Fitness值及其对应的坐标,并将其赋给初始坐标;如高于,则返回步骤(2);
(4)找到C、g的最佳值,建立GWO-SVR预测模型.
GWO-SVR算法流程图如图4所示。
5.针对本发明提出的方法进行试验分析
对模具样件拼接缝中测得的铣削振动信号数据进行CEEMD分解,得到各个分量序列和残余项,然后基于GWO-SVR预测模型对每个分量进行训练和预测,训练和预测时基于GWO算法对SVR的参数c和g进行全局寻优,将各个分量的预测结果相加得到短期铣削振动总预测结果。
基于GWO对SVR中的参数c、g进行参数寻优过程中,首先进行初始化:种群规模N=20,最大迭代次数为200,c和g的范围从0.01到100。铣削振动加速度值被作为预测值输出,算法被用于训练模型时,GWO、PSO和GA三种寻优算法的适应度如图5所示。
三种算法对于适应度函数的计算不同,所以将会产生各自算法的全局最优参数,进而基于不同算法的支持向量机回归预测的精度也有所不同,不同寻优算法的最优参数寻优结果如表2所示。
表2不同方法参数寻优表
不同模型 c g
CEEMD-SVR 12 21
CEEMD-GA-SVR 28 36
CEEMD-PSO-SVR 45 53
CEEMD-GWO-SVR 76 88
图5中可以看出,PSO算法在迭代次数达到第78次左右适应度函数完全收敛得到最优性能,其收敛性较差;而GA算法和GWO算法寻优速度和收敛性均大于PSO算法,其中GA算法在迭代次数为第50次时适应度函数趋于收敛,GWO 算法则是在迭代次数为第5次时适应度函数完全收敛,可以看出GWO算法在寻优速度和收敛性方面性能优越,最终得到的全局最优惩罚因子c=76,全局最优核函数参数g=88。
建立SVR滚动预测模型,使SVR预测模型充分学习数据集内的数据规律,本文采用铣削振动信号中的第1~100时序数据点预测第101点,用第2~101时序数据点预测第102点,输入记为X,输出记为Y,一个X和Y记为一组样本,共计753组;前452组样本构建为训练集(占总样本数的60%),后301组样本构建为测试集(占总样本数的40%)。训练集是预测模型用来训练数据规律的数据集,当预测模型训练完成过后得到数据集的变化规律,然后预测模型根据预测集的步长预测步长相对的数据,得到一段模型预测后的数据,最后把预测数据和训练集的真实数据进行图形对比。
对经过CEEMD算法分解得到的加速度信号(IMFs信号)通过GWO-SVR预测模型进行预测,并将预测曲线进行反归一化处理;然后将各个经过反归一化处理后的预测曲线进行预测信号结果重构,得到振动加速度总信号(IMFs信号) 的预测曲线,预测曲线如图6所示。
从图6中可以看出,不同切深下的加速度信号预测结果图,可以明显地看出GWO-SVR模型对工件X方向的预测曲线拟合度高于工件Y方向的预测曲线拟合度,所以工件X方向的预测精度是高于工件Y方向的。预测准确度的计算采取的方法是先计算预测值与实际值的差值,然后与实际值进行比较,对比图如表3所示。
表3工件X、Y方向铣削振动加速度信号预测准确度
Figure BDA0002847988580000141
如表3所示,在工件X方向上,不同切深所带来的预测精度也是不同的,当切深ae=0.30mm时,加速度信号预测误差最小,当预测时间为0.24s时工件X方向准确度最高达到90.24%,当预测时间为0.12s时工件X方向准确度最高达到 96.96%,整体准确度(工件X方向和Y方向两者共同的准确度)为94.49%。分析中分别基于CEEMD-SVR模型、CEEMD-GA-SVR模型、CEEMD-PSO-SVR模型和 CEEMD-GWO-SVR模型对数据进行预测,得到了四条预测曲线,如图7所示。
从图7可以看出基于CEEMD-SVR模型的预测数据精准度要高于CEEMD- GA-SVR模型、CEEMD-PSO-SVR模型和CEEMD-GWO-SVR模型的预测精准度,且在数据集开始处和数据集中部差距比较明显,具体预测数据见表4。
本文选取两个训练集和测试集比例,一个是训练集为总体数据量的60%,测试集为总体数据量的40%,另一个是训练集为总体数据量的20%,测试集为总体数据量的80%。当训练集时序数据减少,预测时序数据增加1倍时,观察所预测数据的精准度变化情况如图8和图9所示。
当训练集的数据比例由60%下调到了20%,预测误差明显上升,预测精准度下降,这是由于训练集比例的降低会造成SVR所训练的样本个数的下降,训练样本数不足导致了SVR模型无法充分地学习数据样本的规律,从而导致预测精度的下降,为了进一步定量地评估预测结果,分别计算预测时间分别为为0.12s和0.24s的MSE(均方误差),计算结果如表4、表5所示。
表4预测时间为0.12S的不同优化算法的性能对比
Figure BDA0002847988580000151
表4预测时间为0.12S的不同优化算法的性能对比
Figure BDA0002847988580000152
Figure BDA0002847988580000161
从表4和表5的实验结果可以看出,预测时间为0.12s的预测步长模型的 MSE要优于预测时间为0.24s的预测步长模型的MSE,比如切深ae=0.30mm时,基于CEEMD-GWO-SVR模型的工件X方向0.12s预测步长的MSE为10.24要小于 0.24s预测步长MSE的15.31,所以0.12s预测步长模型的预测性能要优于0.24s预测步长模型的预测性能。
其中切深ae=0.30mm时同一预测模型的MSE要低于其他切深预测模型的 MSE,而且工件X方向MSE低于工件Y方向的MSE,如表4,当切深为ae=0.30mm 时工件X方向的MSE为10.24,当切深为ae=0.10mm时工件X方向的MSE为14.15,当切深为ae=0.20mm时工件X方向的MSE为17.38,切深为ae=0.40mm时工件X 方向的MSE为13.02;当切深ae=0.30mm时工件X方向的MSE为10.24,切深为 ae=0.30mm时工件Y方向MSE为29.60。
在对工件X方向加速度信号进行步长为0.12s的预测且切深ae=0.30mm时,CEEMD-SVR模型、CEEMD-GA-SVR模型、CEEMD-PSO-SVR模型和CEEMD- GWO-SVR模型各自的MSE为68.14、39.96、31.87和10.24,CEEMD-GWO-SVR 模型预测得到的MSE误差要低于CEEMD-SVR模型、CEEMD-GA-SVR模型、 CEEMD-PSO-SVR模型,所以相比于CEEMD-SVR模型、CEEMD-GA-SVR模型、 CEEMD-PSO-SVR模型,CEEMD-GWO-SVR模型具有更小的预测误差,证明了此方法更适用于预测铣削加速度信号。
通过对本发明的应用得出如下结论:
(1)提出了针对拼接模具铣削振动信号的CEEMD分解的方法,采取分解6层IMFs和残余项的方法,得到了稳定的子序列供后续的前瞻预测;
(2)采用了先进的灰狼优化算法对支持向量机回归(SVR)中的参数c和 g进行寻优,得到铣削振动信号规定时间段内的全局最优参数,然后将CEEMD 和GWO-SVR结合起来预测铣削拼接模具振动加速度信号。实验数据显示,在对工件X方向加速度信号进行步长为0.12s的预测且切深ae=0.30mm时, CEEMD-SVR模型、CEEMD-GA-SVR模型、CEEMD-PSO-SVR模型和CEEMD- GWO-SVR模型各自的MSE为68.14、39.96、31.87和10.24,CEEMD-GWO-SVR 模型的MSE低于其他三种模型的MSE,证明了本文所选方法相比其他方法具有更高的预测精度和更好的泛化能力,在预测步长为0.12s时总体预测准确率为 94.49%。
(3)在淬硬钢拼接区铣削过程中,当切深为0.10mm、0.20mm、0.30mm 和0.40mm时各自的MSE为14.15、17.38、10.24、13.02,得到结论,在切深 ae 0.30mm时,拼接区铣削振动信号预测误差最小;在振动加速度信号为 0.3s的时间段内,0.12s预测步长MSE为17.38,0.24s预测步长MSE为23.89。得到结论,在0.3s的数据集中,0.12s预测步长模型(预测集为全部数据集的 40%)的预测误差小于0.24s预测步长模型(预测集为全部数据集的80%)。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (2)

1.基于CEEMD和GWO-SVR的振动信号预测方法,所述方法的实现过程为:
步骤一、基于互补集合经验模态分解(CEEMD)将实时采集的时序振动信号进行6层分解,得到的各层本征模态函数及趋势序列;
步骤二、对得到的各层本征模态函数及趋势序列构建训练集和预测集,基于构建训练集和预测集对利用GWO优化支持向量回归的参数C和g进行寻优,建立GWO-SVR预测模型:
构建GWO的过程:
(1)追踪猎物
灰狼在狩猎过程中围绕猎物,将其行为做出以下定义:
Figure FDA0002847988570000011
Figure FDA0002847988570000012
式中,
Figure FDA0002847988570000013
代表灰狼经过第t次迭代后所得出的猎物当前所在位置,
Figure FDA0002847988570000014
代表灰狼个体经过第t次迭代后当前所在位置,这就是算法的当前局部最优解。
Figure FDA0002847988570000015
Figure FDA0002847988570000016
式中,
Figure FDA0002847988570000017
Figure FDA0002847988570000018
为包围猎物中的随机系数,
Figure FDA0002847988570000019
的区间为[0,2]且随着迭代次数的增加而降低,
Figure FDA00028479885700000110
Figure FDA00028479885700000111
是区间[0,2]上的随机数。
(2)包围猎物
灰狼有识别猎物的强能力,且α狼可以领导灰狼包围猎物,但在猎物(最优解)位置未知的情况下。为了模拟灰狼捕猎行为,假定α狼、β狼和δ狼了解猎物的潜在位置,保存三个最佳解决方案,并且要求其他灰狼(ω狼)作为搜索代理,根据最佳解决方案进行灰狼位置的迭代更新,公式为:
Figure FDA00028479885700000112
式中,
Figure FDA00028479885700000113
Figure FDA00028479885700000114
分别表示α狼、β狼和δ狼与ω狼之间的距离;
Figure FDA00028479885700000115
Figure FDA00028479885700000116
分别表示第t次迭代后α狼、β狼和δ狼的当前位置;
Figure FDA00028479885700000117
表示当前灰狼位置。
Figure FDA00028479885700000118
Figure FDA0002847988570000021
公式(17)代表ω狼朝三个潜在解α狼、β狼和δ狼的步长和方向;公式(7)代表ω狼的最终位置。
(3)攻击猎物
当猎物停止移动时狩猎完成,这一过程可以通过降低
Figure FDA0002847988570000022
来完成。当随机值
Figure FDA0002847988570000023
在[-1,1]时,搜索代理的下一位置可以在当前位置与猎物位置之间的位置进行更新。算法中,当
Figure FDA0002847988570000024
时强制狩猎。
(4)寻找猎物
灰狼常常分头搜寻猎物,同时根据α狼、β狼和δ狼的位置来攻击猎物,当
Figure FDA0002847988570000025
时,灰狼被迫离开猎物,
Figure FDA0002847988570000026
灰狼强制狼群攻击猎物,直到满足条件时,GWO算法停止。
建立GWO-SVR模型,其过程为:
(1)基于GWO对SVR中的参数C、g进行参数寻优过程中,首先进行初始化:种群规模N=20,最大迭代次数为200;确定SVR中参数C、g的范围,即C∈[0.01,100],g∈[0.01,100];
(2)将所述训练样本对输入到SVR模型中,对模型进行训练,将均方根误差RMSE作为适应度函数Fitness;
(3)找到适应度函数Fitness的最小值对应的个体,开始迭代寻优,并判断最小Fitness是否低于前一代最小Fitness;如低于,则保留最小Fitness值及其对应的坐标,并将其赋给初始坐标;如高于,则返回步骤(2);
(4)找到C、g的最佳值,建立GWO-SVR预测模型;
步骤三、采用滚动预测模型,铣削振动信号中的第1~100时序数据点预测第101点,用第2~101时序数据点预测第102点,输入记为X,输出记为Y,一个X和Y记为一组样本,总样本数的60%构建为训练集,总样本数的40%构建为测试集;
步骤四、对经过CEEMD算法分解得到的加速度信号(IMFs信号)通过GWO-SVR预测模型进行预测,并将预测曲线进行反归一化处理;然后将各个经过反归一化处理后的预测曲线进行预测信号结果重构,得到振动加速度总信号(IMFs信号)的预测曲线;
步骤五、选取两个训练集和测试集比例,一个是训练集为总体数据量的60%,测试集为总体数据量的40%,另一个是训练集为总体数据量的20%,测试集为总体数据量的80%。当训练集时序数据减少,预测时序数据增加1倍时,观察所预测数据的精准度变化情况。
2.根据权利要求1所属基于CEEMD和GWO-SVR的振动信号预测方法,其特征在于:
在步骤一中,通过互补集合经验模态分解(CEEMD)将实时采集的时序振动信号进行6层分解,得到的各层本征模态函数及趋势序列:
步骤一一、获得拼接模具铣削振动信号;
步骤一二、基于EMD模型加入互补白噪声,建立CEEMD模型:
CEEMD算法具体步骤如下:
(1)将一组符号相反的噪声加入到目标信号中,每组噪声的幅值相同,然后产生一对新信号:
s+ i(t)=x(t)+ni +(t) (8)
s- i(t)=x(t)+ni -(t) (9)
式中,si(t)是新信号,ni +(t)是正噪声信号,ni -(t)是负噪声信号。
(2)分解新信号,得到两组含有白噪声的信号分量IMF1和IMF2,且这两组分量正负互补。
(3)将两组信号的IMF分量进行整合得到最终的IMF分量。
步骤一三、通过CEEMD分解振动信号获得IMFs。
CN202011517791.7A 2020-12-21 2020-12-21 基于ceemd和gwo-svr的振动信号预测方法 Active CN112749506B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011517791.7A CN112749506B (zh) 2020-12-21 2020-12-21 基于ceemd和gwo-svr的振动信号预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011517791.7A CN112749506B (zh) 2020-12-21 2020-12-21 基于ceemd和gwo-svr的振动信号预测方法

Publications (2)

Publication Number Publication Date
CN112749506A true CN112749506A (zh) 2021-05-04
CN112749506B CN112749506B (zh) 2022-08-26

Family

ID=75648153

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011517791.7A Active CN112749506B (zh) 2020-12-21 2020-12-21 基于ceemd和gwo-svr的振动信号预测方法

Country Status (1)

Country Link
CN (1) CN112749506B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113222273A (zh) * 2021-05-26 2021-08-06 中铁十八局集团有限公司 一种基于gwo-fw-mkl-svr算法的tbm掘进速率预测方法
CN113300884A (zh) * 2021-05-19 2021-08-24 内蒙古农业大学 一种基于gwo-svr的分步网络流量预测方法
CN113567129A (zh) * 2021-07-28 2021-10-29 大连海天兴业科技有限公司 一种列车轴承振动信号基于ceemd的降噪方法
CN114398732A (zh) * 2022-01-11 2022-04-26 福州大学 一种cfrp铣削工艺参数的优化方法及系统
CN114676736A (zh) * 2022-05-13 2022-06-28 浙江工业大学 一种滚压工件表面质量预测方法
CN115327385A (zh) * 2022-07-29 2022-11-11 武汉理工大学 一种动力电池soc值估算方法及系统
US11879869B2 (en) 2022-05-13 2024-01-23 Zhejiang University Of Technology Method for predicting surface quality of burnishing workpiece

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107036817A (zh) * 2017-04-05 2017-08-11 哈尔滨理工大学 基于磷虾群算法的svr滚动轴承性能衰退预测方法
CN109242215A (zh) * 2018-10-26 2019-01-18 浙江工业大学之江学院 一种基于粒子群-支持向量机的旋转机械设备运行工况预测方法
US20190175041A1 (en) * 2017-12-11 2019-06-13 Adaptive, Intelligent And Dynamic Brain Corporation (Aidbrain) Method, module and system for analysis of physiological signal
CN111581829A (zh) * 2020-05-11 2020-08-25 辽宁大学 一种改进差分灰狼算法优化支持向量回归的带钢厚度预测方法
CN111859730A (zh) * 2020-06-05 2020-10-30 合肥通用机械研究院有限公司 一种燃料电池离心压缩机转子构形优化设计方法
CN111967187A (zh) * 2020-08-20 2020-11-20 新奥能源动力科技(上海)有限公司 一种预测燃气轮机启动过程输出数据的建模方法及设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107036817A (zh) * 2017-04-05 2017-08-11 哈尔滨理工大学 基于磷虾群算法的svr滚动轴承性能衰退预测方法
US20190175041A1 (en) * 2017-12-11 2019-06-13 Adaptive, Intelligent And Dynamic Brain Corporation (Aidbrain) Method, module and system for analysis of physiological signal
CN109242215A (zh) * 2018-10-26 2019-01-18 浙江工业大学之江学院 一种基于粒子群-支持向量机的旋转机械设备运行工况预测方法
CN111581829A (zh) * 2020-05-11 2020-08-25 辽宁大学 一种改进差分灰狼算法优化支持向量回归的带钢厚度预测方法
CN111859730A (zh) * 2020-06-05 2020-10-30 合肥通用机械研究院有限公司 一种燃料电池离心压缩机转子构形优化设计方法
CN111967187A (zh) * 2020-08-20 2020-11-20 新奥能源动力科技(上海)有限公司 一种预测燃气轮机启动过程输出数据的建模方法及设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SULING ZHU等: "PM2.5 forecasting using SVR with PSOGSA algorithm based on CEEMD,GRNN and GCA considering meteorological factors", 《ATMOSPHERIC ENVIRONMENT》, vol. 183, 30 June 2018 (2018-06-30), pages 20 - 32 *
XING LU等: "Prediction and analysis of cold rolling mill vibration based on a data-driven method", 《APPLIED SOFT COMPUTING JOURNAL》, vol. 96, 30 November 2020 (2020-11-30) *
宋祎: "基于CEEMD和特征融合的高速列车振动信号特征分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, 31 January 2017 (2017-01-31), pages 033 - 320 *
邓聪颖等: "基于SVR-GA算法的广义加工空间机床切削稳定性预测与优化研究", 《仪器仪表学报》, vol. 40, no. 10, 31 October 2019 (2019-10-31), pages 227 - 236 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113300884A (zh) * 2021-05-19 2021-08-24 内蒙古农业大学 一种基于gwo-svr的分步网络流量预测方法
CN113300884B (zh) * 2021-05-19 2022-07-08 内蒙古农业大学 一种基于gwo-svr的分步网络流量预测方法
CN113222273A (zh) * 2021-05-26 2021-08-06 中铁十八局集团有限公司 一种基于gwo-fw-mkl-svr算法的tbm掘进速率预测方法
CN113222273B (zh) * 2021-05-26 2022-06-07 中铁十八局集团有限公司 一种基于gwo-fw-mkl-svr算法的tbm掘进速率预测方法
CN113567129A (zh) * 2021-07-28 2021-10-29 大连海天兴业科技有限公司 一种列车轴承振动信号基于ceemd的降噪方法
CN114398732A (zh) * 2022-01-11 2022-04-26 福州大学 一种cfrp铣削工艺参数的优化方法及系统
CN114398732B (zh) * 2022-01-11 2024-09-06 福州大学 一种cfrp铣削工艺参数的优化方法及系统
CN114676736A (zh) * 2022-05-13 2022-06-28 浙江工业大学 一种滚压工件表面质量预测方法
CN114676736B (zh) * 2022-05-13 2022-08-26 浙江工业大学 一种滚压工件表面质量预测方法
US11879869B2 (en) 2022-05-13 2024-01-23 Zhejiang University Of Technology Method for predicting surface quality of burnishing workpiece
CN115327385A (zh) * 2022-07-29 2022-11-11 武汉理工大学 一种动力电池soc值估算方法及系统

Also Published As

Publication number Publication date
CN112749506B (zh) 2022-08-26

Similar Documents

Publication Publication Date Title
CN112749506B (zh) 基于ceemd和gwo-svr的振动信号预测方法
Çaydaş et al. An adaptive neuro-fuzzy inference system (ANFIS) model for wire-EDM
Chen et al. Predicting tool wear with multi-sensor data using deep belief networks
Markopoulos et al. Artificial neural network models for the prediction of surface roughness in electrical discharge machining
Jang et al. Modeling and parameter optimization for cutting energy reduction in MQL milling process
He et al. Milling tool wear prediction using multi-sensor feature fusion based on stacked sparse autoencoders
Kant et al. Predictive modelling and optimization of machining parameters to minimize surface roughness using artificial neural network coupled with genetic algorithm
Li et al. Multi-objective optimization of cutting parameters in sculptured parts machining based on neural network
Fei et al. Springback prediction for incremental sheet forming based on FEM-PSONN technology
Markopoulos et al. Modelling and optimization of machining with the use of statistical methods and soft computing
Rashid et al. Surface roughness prediction for CNC milling process using artificial neural network
Hamzaçebi et al. Taguchi method as a robust design tool
Akkuş et al. Predicting surface roughness of AISI 4140 steel in hard turning process through artificial neural network, fuzzy logic and regression models
Ming et al. Integrated ANN-LWPA for cutting parameter optimization in WEDM
CN114384861A (zh) 一种基于多层感知机的数控机床多工况切削参数优化方法
Lin et al. Application of the fuzzy-based Taguchi method for the structural design of drawing dies
Anitha et al. Multi-objective optimization of electrical discharge machining processes using artificial neural network.
Aich et al. Application of teaching learning based optimization procedure for the development of SVM learned EDM process and its pseudo Pareto optimization
Al Hazza et al. Investigation of the effect of cutting speed on the Surface Roughness parameters in CNC End Milling using Artificial Neural Network
Li et al. A data and knowledge-driven cutting parameter adaptive optimization method considering dynamic tool wear
CN116520772A (zh) 基于多任务迁移学习的薄板零件数控加工质量预测方法
Zhang et al. Tool life prediction model based on GA-BP neural network
Muthuram et al. Optimization of machining parameters using artificial Intelligence techniques
Zhang et al. Tool wear monitoring using a novel parallel BiLSTM model with multi-domain features for robotic milling Al7050-T7451 workpiece
Gao et al. A 6-sigma robust optimization method for stamping forming of automobile covering parts based on residual error and radial basis interpolation

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