基于核函数-Wiener过程的轧辊非线性退化建模与剩余寿命预测.pdf

收稿日期:修回日期:录用日期:基金项目:国家重点研发计划项目(Y F B ),国家自然科学基金(,)作者简介:王汉禹(),硕士生,从事设备退化建模与寿命预测研究通信作者:潘尔顺,教授,博士生导师;E m a i l:p e s s j t u e d u c n 文章编号:()D O I:/j c n k i j s j t u 基于核函数 W i e n e r过程的轧辊非线性退化建模与剩余寿命预测王汉禹,陈震,周笛,陈兆祥,潘尔顺(上海交通大学 机械与动力工程学院,上海 )摘要:在轧钢生产过程中,由于磨损等原因,长时间复杂工况运行下的轧辊工作性能会呈现逐步衰退的特征考虑轧辊工作环境具有工况复杂、随机干扰强等特点,提出基于核函数 W i e n e r过程(KWP)退化模型,使用W i e n e r过程刻画轧辊退化趋势的强随机性特征,引入核函数捕捉轧辊的非线性退化路径,推导贝叶斯框架下参数估计的解析表达式,并构造轧辊可工作转动量的健康指标,进一步预测轧辊剩余寿命(RU L)结合某钢铁公司 热轧生产线现场数据,所构建模型拟合优度达到 ,剩余寿命预测误差低于 ,相较常见机器学习算法取得了更好效果,有助于提高设备运转效率与安全性并实现视情维护关键词:轧辊;性能退化;非线性建模;核函数;W i e n e r过程;剩余寿命预测中图分类号:TH 文献标志码:AN o n l i n e a rD e g r a d a t i o nM o d e l i n ga n dR e s i d u a l L i f eP r e d i c t i o n f o rR o l l e r sB a s e do nK e r n e l b a s e dW i e n e rP r o c e s sWANG H a n y u,CHENZ h e n,ZHO UD i,CHENZ h a o x i a n g,P ANE r s h u n(S c h o o l o fM e c h a n i c a lE n g i n e e r i n g,S h a n g h a i J i a oT o n gU n i v e r s i t y,S h a n g h a i ,C h i n a)A b s t r a c t:I nt h ep r o c e s so fs t e e l r o l l i n g,d u et ow e a ra n do t h e rr e a s o n s,t h ew o r k i n gp e r f o r m a n c eo f t h er o l lu n d e r l o n ga n dc o m p l e xw o r k i n gc o n d i t i o n sh a sag r a d u a ld e c l i n e C o n s i d e r i n gt h ec h a r a c t e r i s t i c so fc o m p l e xw o r k i n gc o n d i t i o n sa n ds t r o n gr a n d o mi n t e r f e r e n c eo ft h er o l lw o r k i n ge n v i r o n m e n t,t h i sp a p e rp r o p o s e dak e r n e l b a s e dW i e n e rp r o c e s s(KWP)d e g r a d a t i o nm o d e l t oc h a r a c t e r i z et h es t r o n gr a n d o m n e s so f t h er o l l d e g r a d a t i o n t r e n db yu s i n g t h eW i e n e rp r o c e s s,a n d t oc a p t u r e t h en o n l i n e a rd e g r a d a t i o np a t ho ft h er o l l b yu s i n g t h ek e r n a f u n c t i o n T h i sp a p e r d e r i v e s t h e a n a l y t i c a l e x p r e s s i o no f p a r a m e t e r e s t i m a t i o n i nt h eB a y e s i a nf r a m e w o r k,a n dc o n s t r u c t st h eh e a l t hi n d e xo ft h er o l lw o r k i n gr o t a t i o n,t h e np r e d i c t st h er e m a i n i n gu s e f u l l i f e(RU L)o f t h er o l l I nc o m b i n a t i o nw i t ht h e f i e l dd a t ao f h o t r o l l i n gp r o d u c t i o nl i n eo fa ni r o na n ds t e e lc o m p a n y,t h eg o o d n e s so ff i to ft h em o d e lb u i l t i s ,a n dt h er e s i d u a l l i f ep r e d i c t i o ne r r o ri sl e s st h a n C o m p a r e d w i t ht h ec o mm o n m a c h i n el e a r n i n ga l g o r i t h m,i th a sa c h i e v e db e t t e r r e s u l t s,w h i c hi sh e l p f u l t oi m p r o v et h eo p e r a t i n ge f f i c i e n c ya n ds a f e t yo fe q u i p m e n ta n da c h i e v em a i n t e n a n c ea sn e e d e d K e yw o r d s:r o l l e r;p e r f o r m a n c ed e g r a d a t i o n;n o n l i n e a rm o d e l i n g;k e r n e l f u n c t i o n;W i e n e rp r o c e s s;r e m a i n i n gu s e f u l l i f e(RU L)p r e d i c t i o n第 卷 第期 年月上 海 交 通 大 学 学 报J OUR NA LO FS HANGHA I J I A OT ON GUN I V E R S I T YV o l N o A u g 上海交通大学学报第 卷我国钢铁产量已多年位居世界第一,但作为典型的流程型制造行业,其生产效率和经济效益仍与世界先进水平存在较大差距轧辊是轧钢生产线上的关键部件和大型耗材,表面磨损是其主要失效形式,并对带钢表面质量及板形产生不利影响针对轧 辊 失 效 过 程 建 立 退 化 模 型 并 实 现 剩 余 寿 命(RU L)预测,对于提高轧辊使用寿命和产品质量尤为重要,也可为管理人员确定换辊周期与维修策略提供决策依据,最终达到节能减排、降低成本的目标,为推动绿色生产和碳中和做出贡献但由于轧辊工况恶劣,通常是在高压、高温、高速等条件下工作,且轧辊与金属表面间存在氧化铁皮、冷却水等多种影响因素,导致轧辊退化过程具有随机性、非线性、强耦合等特征,所以对轧辊进行准确的退化建模与RU L预测是当前研究难点目前,针对轧辊的退化建模研究,常用方法可分为两类:物理解析模型和基于数据驱动的方法物理解析模型通过研究轧辊材料属性、失效机理、负载条件等影响退化的因素,构建轧辊退化数学模型并进行RU L预测刘子英等和C a i等指出,轧辊退化过程涉及轧制力、带钢宽度、温度、速度、冷却水等众多影响因素,且这些因素具有耦合关系,并伴随机械磨损、热磨损、化学腐蚀等多种退化形式多数研究通常基于轧制力、轧制长度、轧辊材质、磨损距离等主要影响因素构建经典模型,也有研究在其基础上考虑轧辊表面温度、轧制速度、不均匀磨损等因素提出改进模型考虑轧辊磨损是一个极其复杂的过程,涉及物理、化学及材料学等众多学科,这无疑限制了物理解析模型在轧辊退化建模领域中的发展基于历史大数据,通过构建输入和输出间的映射关系,数据驱动方法可以避免物理模型完全还原复杂因果关系的缺陷目前,面向轧辊退化领域的数据驱动方法研究相对较少,仅少数学者将神经网络应用于轧辊 的 退 化 建 模 和RU L预 测如 吴 进等 采用反向传播(B P)神经网络与有限元分析软件,通过对摩擦因数进行预测得到了较高精度的轧制力预测结果 J i a o等 提出了一种新的基于递归神经网络的深度学习网络结构,从批量数据中提取健康指标预测轧辊的RU L神经网络方法在工程应用中需要大量轧辊历史数据来训练模型,否则难以较好地预测轧辊性能此外,神经网络模型难以捕捉退化数据的随机性,一般只能给出相关特征量的点估计,不能表征退化过程的不确定性鉴于轧辊退化的影响因素众多,难以建立有效的物理解析模型,而神经网络方法所需数据量大且对随机数据处理能力有限考虑随机过程能够较好表征设备退化量的随机性,可以给出RU L的概率分布,并能结合具有较强非线性映射能力的核函数,因此提出基于核函数 W i e n e r过程(KWP)退化模型,并在贝叶斯框架下推导参数估计的解析形式同时基于轧辊磨损机理,提出基于轧辊可工作转动量的健康指标,构建基于KWP退化模型的健康状态估计和RU L预 测 流 程所 提 方 法 在 某 钢 铁 公 司 热轧生产线实测数据上取得了较好的建模效果和准确的RU L预测,有助于提高生产线作业计划的稳定性和连续性 问题描述轧辊磨损是典型的逐步累积的非线性退化过程,其磨损速率随着轧制长度(k m)增加逐渐加快,主要原因包括三方面:轧辊表面经过特殊热处理形成耐磨层,耐磨层磨损后进入快速磨损期;随着轧辊直径减小,导致轧制相同距离的轧辊表面所受冲击频率更高;轧辊越细,单位轧辊表面所受压力更大,磨损更快受轧制压力、热冲击和热应力等多种因素耦合作用的影响,轧辊退化路径具有复杂的非线性和随机性特征考虑随机过程能够较好地处理设备退化过程中的随机性特征 ,并已在锂电池、数控机床、制动器、轴承等设备及部件上取得了较好的应用效果,它在轧辊退化建模领域具有巨大潜力,因此本文使用随机过程对轧辊进行可靠性建模并预测其RU L退化建模中常用的随机过程模型包括W i e n e r过程、G a mm a过程和I n v e r s eG a u s s i a n过程等,由于W i e n e r过程更适用于非单调退化路径的建模且具备良好的数学可追溯性,目前已成为可靠性技术研究和应用的重要方法之一 假设X(t)为t时刻的退化量,则满足以下性质的随机过程:X(t),t,定义为一元线性W i e n e r过程()任意时刻的增量服从正态分布,即:XX(tt)X(t)N(t,t)()式中:和分别表示漂移系数和扩散系数()任意两个不相交时间段内的增量相互独立相应的数学描述为:对于任意的ttn,则X(t)X(t),X(t)X(t),X(tn)X(tn)相互独立()X()线性W i e n e r过程难以描述非线性退化的动态第期王汉禹,等:基于核函数 W i e n e r过程的轧辊非线性退化建模与剩余寿命预测 特性,而轧辊历史数据表明其性能退化路径具有明显非线性特征,因此有必要进一步研究非线性W i e n e r过程模型典型非线性W i e n e r过程退化模型可由下式给出:X(t)(t)B(t)()式中:(t)为时间尺度转换函数;扩散系数可用于调整函数B(t)布朗运动的扩散程度非线性W i e n e r过程的研究大多是假设时间尺度转换函数满足指数形式,即(t)t然而轧辊的磨损是一个多因素相互耦合的复杂非线性过程,基于指数函数的W i e n e r过程只适用于小部分非线性退化问题,需研究形式更灵活的漂移函数以提高对复杂非线性问题的建模能力非线性建模可以通过核函数将低维数据映射到高维空间进行非线性变换,使其在高维空间中呈现线性特征 本文基于此思想在W i e n e r过程中引入核函数来解决轧辊退化过程中的非线性表征问题相较基于指数函数的W i e n e r过程,KWP具有如下优点:()针对轧辊退化趋势因工况复杂、随机干扰强造成的强随机性特征,将W i e n e r过程应用于轧辊的RU L预测中相较神经网络等人工智能方法,W i e n e r过程可以表征样本差异、测量误差、环境干扰等不确定信息,且推导得到的失效分布蕴含首达时的概念,为制定合理的预防维护方案提供决策支持表常用核函数及其公式T a b K e r n e l f u n c t i o n sa n df o r m u l a s核函数表达式参数线性核函数t ti无特殊参数高斯核函数e x p|tti|v()v为高斯核函数的宽度多项式核函数(t tir)d,d,nr为偏置系数;d为多项式阶数L a p l a c e核函数e x p|tti|o()oS i g m o i d核函数t a n h(t ti),()考虑现有随机过程模型对复杂非线性退化路径的刻画能力较弱,引入核函数构建基于KWP的退化模型该模型能够有效描述轧辊退化过程的随机性和非线性,有助于实现轧钢生产线从当前定期维护策略到基于状态维护策略的转变()相较神经网络等常见数据驱动方法,基于KWP可以获取RU L预测的概率密度函数,度量退化的不确定性,并为其他类似领域的RU L预测和实时状态监测问题提供参考 轧辊性能退化过程非线性建模 基于KWP的退化模型引入核函数代替式()中非线性W i e n e r过程的漂移函数项(t),构 建基于KWP的 退 化 模 型如下:X(t)(miwik(t,ti)w)B(t)()式中:k(t,ti)(t)(ti)为核函数;wi为权重系数;w为截距;ti为第i次测量的时间;m为测量总次数定义y(t|w)miwik(t,ti)w,wwwwmT式()中k(t,ti)为核函数的一般形式,常用核函数的表达式如表所示核函数性能取决于其形式和参数,也决定了模型对非线性问题的处理能力,从而影响KWP模型对轧辊性能预测的准确性与有效性为使核函数能够更好地刻画轧辊性能退化所特有的路径特征,在实际应用中可根据轧辊实测数据的预测效果,选取核函数嵌入到模型中的漂移函数项KW P模 型 的 独 立 增 量 为xxxxnT,其联合概率密度函数服从多元高斯分布:p(x|w,)()n|E|e x p(x(x)(xn)Tw)TE(x(x)(xn)Tw)()式 中:(xi)k(xi,x)k(xi,xn);E为 元 素 的 先 验 协 方 差 矩 阵,E(i,j)m i n(ti,tj)对于未知参数w和,其取值需根据轧辊退化历史数据进行估计 基于退化模型的贝叶斯参数估计针对模型参数估计,首先基于贝叶斯理论推导参数w和的后验分布,获得边际似然函数,再通过 上海交通大学学报第 卷最大化边缘似然函数,便可得到参数w和的最大似然估计值使用极大似然估计法估计参数易使模型出现过拟合情况,为避免模型出现过拟合现象,定义参数w的先验分布为均值等于的高斯分布,即:p(w|)mii e x piwi()式中:m为w的超参数向量根据多层贝叶斯先验分布,定义和的超先验分布服从G a mm a分布,即:p()miG a mm a(i|a,b)()p()G a mm a(|c,d)()式中:G a mm a(i|a,b)(a)baa eb a(a)ta etdt在选取无信息先验作为和的超先验分布时,参数a、b、c和d应足够小,如abcd 基于贝叶斯公式,未知参数w、和的后验分布可表示为p(w,|x)p(w|x,)p(,|x)()由贝叶斯定理,可得w计算式如下:p(w|x,)p(x|w,)p(w|)p(x|,)()(m)/|E|/e x p(w)TE(w)()参数和的估计值可以通过最大化边际似然函数p(,|x)得到,即使p(x|,)最大:p(x|,)p(x|w,)p(w|)dw()m A T|e x p xT(A T)x()式中:(x)(x)(xn)T;Ad i a g(a,am)式()的后验协方差和均值分别为C(Td i a g(,m)()(d i a g(,m)TE)TE x()设计迭代估计算法求得式()中超参数和的近似解按照M a c k a y 的方法令式()等于,整理得和的迭代更新值为n e wii/i()()n e wxmEi ii()式中:i为式()的第i个后验平均权重;定义iiEi i;Ei i为式()后验方差矩阵的第i个对角元素,可由和值计算得到迭代中若i很大,且wi被其先验概率严格约束时,则令Ei i i,i;若i很小,且wi适应数据时,令i针对单组轧辊退化向量样本x,为了得到相应的预测值x,基于概率论和贝叶斯定理可得:p(xx)p(xw,)p(w,x)dwdp(xw,)p(w,|x)dwdd()将式()代入式()中可得预测值x的条件概率为p(x|x)p(x|w,)p(w,|x)dwdd()因难以计算积分p(x)p(x|w,)p(w,)dwdd,无法直接得到式()中的后验概率p(w,|x),所以依据贝叶斯定理将后验概率分解为p(w,|x)p(w|x,)p(,|x)即p(x|x)p(x|w,)p(w|x,)p(,|x)dwdd()基于式()和()更新参数,同时计算式()和()更新后验统计C和,直到达到最大循环次数或者输出结果的梯度小于收敛条件,即 由于在迭代过程中易发生部分i趋于无穷的现象,根据核函数性质 删除迭代过程中所有im a x的i所对应的权重和核函数,并设m a x,通过此过程实现KWP退化模型的稀疏性,提升了工程分析效率以适用轧辊性能状态的在线监测记MP和MP为停止迭代后的最优值(MP),此时若给定一组新的轧辊退化数据x,其预测输出可由下式得到:p(x|x,MP,MP)p(x|w,MP)p(w|x,MP,MP)dw()由于等式右边的两个积分项均满足高斯分布,所以预测值x满足分布如下:p(x|x,MP,MP)N(x|y,)()式中:预测输出x的均值y和方差分别为第期王汉禹,等:基于核函数 W i e n e r过程的轧辊非线性退化建模与剩余寿命预测 yT(x)MP(x)TE(x)()由式()可知,本模型的预测误差由两部分组成:数据固有噪声和模型权重估计过程中的不确定性在一般的机器学习算法中,使用极大似然估计法估计参数易出现过拟合现象,KWP模型通过约束参数w和调整扩散参数 B(t),可以较好地规避过拟合的情况针对轧辊性能指标退化曲线具有随机性和非线性的特征,构建了KWP退化模型,并基于贝叶斯推理给出模型参数的估计值以及预测值的概率分布基于KWP的轧辊健康状态预测 健康指标及可解释性健康指标的构建是轧辊可靠性建模和寿命预测中不可缺少的一环,应满足以下条件:直观显示轧辊的磨损速度,判断轧辊是否进入快速磨损期;为轧辊下机周期的确定提供参考,以减少总停机时间并保证产品质量,进而提高钢铁生产线的经济效益;表征轧辊直径与工作距离间的耦合关系,为维修策略提供参考若采用累积距离或轧辊直径作为健康指标,并不能反映轧辊磨损速度随着轧辊直径减少而加剧的特征,对确定轧辊更换、维修决策的参考价值有限基于轧辊磨损机理,构建轧辊健康指标如下:轧辊的直径变化过程如图所示,F F 代表产线中存在的个精轧机组,R、R分别代表轧制前后的轧辊厚度,且轧辊的起始工作点Bn位于第n次工作周期开始前的表面Cn上,半径为rn,终止工作点Bn 位于轧辊工作周期结束后磨损的表面,半径为rn,rn为本次磨损半径弧BnBn 为轧辊实际工作里程,n为第n次工作周期实际转过的角度图轧辊磨损示意图F i g S c h e m a t i cd i a g r a mo f r o l l e rw e a r i n g考虑轧辊在单次工作周期前后直径变化十分微小,可近似为rnrn/rn,从而合理认为轧辊实际工作的路程BnBn 相当于CnBn,则第n次工作周期距离的计算式为ln(/)n(rnrn)()设轧辊在完整工作周期内的可工作总距离为l,则推导得轧辊在第n次工作周期下机后的剩余可工作距离lln(/)(n)(rrn),整理得:ly/ry(/)y()式中:ly为轧辊剩余可工作距离;ry为轧辊剩余可磨削直径;y为轧辊可工作转动量考虑ly/ryy,呈正相关关系,以轧辊可工作转动量作为健康指标,刻画轧机工作辊的可服役信息因此所构建指标既考虑了轧辊工作距离与轧辊直径的耦合关系,又具备工程可解释性与健康状态的集成表征性,为后续寿命预测和动态维护提供理论与技术支持 轧辊R U L预测将轧辊历史运行数据:n,n作为输入序列,基于KWP退化模型,通过循环迭代计算得到轧辊未来时刻的健康状态考虑轧辊的失效阈值,轧辊的可工作转动量退化轨迹由KWP退化模型描述,则轧辊寿命T为轧辊的性能退化量X(t)首次达到的时间(F P T):Ti n fT|X(t),t()式中:i n f 表示函数的下确界,则轧辊在任意时间tk的RU L被定义为Hki n fhk:X(tkhk)|X(tk)()W i e n e r过程的首达时间服从逆高斯分布,为估计轧辊F P T的概率密度,基于式()、(),推导得轧辊RU L的概率密度函数:f(t|w)te x p(y(t|w)tdy(t|w)dt()基于式()获得轧辊在任意时刻的RU L预测值,可为确定轧辊更换周期和磨削量提供技术支撑和决策依据,进而提高产品质量和轧辊使用寿命,降低停机损失 退化建模与健康预测流程基于KWP的模型构建和寿命预测流程如图所示输入轧辊历史数据并清洗、正则化,结合轧制距离与轧辊直径构建可工作转动量健康指标在模型训练过程中,首先初始化KWP退化模型的扩散参数和核函数参数根据式()()计算权重的后验协方差C和均值,并更新参数i和,直至 模型收敛得到参数最优值并剔除超参数中i 上海交通大学学报第 卷图KWP退化模型的轧辊健康状态估计方法F i g R o l l e rh e a l t hs t a t ee s t i m a t i o nb a s e do nKWPd e g r a d a t i o nm o d e lm a x的i所对应的权重和核函数,以实现KWP退化模型的稀疏性,即可获得完整预测模型最终结合轧制时间和失效阈值,给出RU L估计值及其概率分布函数,从而实现对轧辊可工作转动量的预测 基于某钢铁公司 热轧生产线的实验分析 实验数据轧辊实验数据来自某钢铁公司 热轧生产线精轧机组,该精轧机组的基本结构如图所示实地调研获得的部分数据如表所示基于轧辊可工作转动量健康指标,对数据进行整理、清洗,考虑数据完整性等因素后,以 热轧生产线F 精轧机图精轧机组结构示意图F i g F i n i s h i n gm i l l s t r u c t u r e表轧辊部分原始数据展示T a b R a wd a t ao f r o l l e r s轧辊编号上机时间下机时间生产量/t轧制距离/k m轧辊直径/mm I T :T :I T :T :I T:T:I T :T :I T :T :组C S 号轧辊的工作周期数据为例进行实验 实验结果与分析 评价指标使用均方根误差(RM S E)、平均绝对误差(MA E)、拟合优度(C D)评价算法的性能RM S E和MA E两者的值越小则表明算法预测精度越高,其公式分别为RM S En n t(x(t)x(t)()MA En n t|x(t)x(t)|()C D值越接近,说明模型预测性能越好,其公式为C Dn t(x(t)x)n t(x(t)x)()式中:n 表示时间序列总数;x(t)表示时刻t的健康状态实际值;x(t)为算法预测值;x为真实值均值第期王汉禹,等:基于核函数 W i e n e r过程的轧辊非线性退化建模与剩余寿命预测 由此,可通过计算绝对值和相对值全方位评估模型的建模效果与预测性能 核 函 数 选 择核 函 数 的 选 择 直 接 影 响KWP退化模型的预测精度,也同样影响模型稀疏性及训练时间基于C S 号轧辊的退化数据,根据核函数不同,实验共分组进行比较,计算各模型的RM S E、C D、MA E进行评估,如表所示表核函数选择T a b K e r n e l f u n c t i o ns e l e c t i o n核函数类型RM S EC DMA E高斯核函数 多项式核函数 线性核函数 L a p l a c e核函数 S i g m o i d核函数 根据表可知,高斯核函数的RM S E、C D和M S E均优其他核函数,说明原数据在通过高斯核函数映射到高维特征空间后线性特征更好,也证明了不同的核函数特性并不相同,不同的应用场景下性能差别很 大在迭代 过程中,基于 高斯核函数 的W i e n e r过 程 退 化 模 型 仅 使 用 了个 相 关 向 量(R V s),具有较好的稀疏性,从而测试时间更短,便于实时在线监测根据实验结果,轧辊性能的退化预测将高斯核函数嵌入W i e n e r过程中作为漂移函数为验证KWP退化模型相较其他常用算法的优越性,使用基于指数函数的W i e n e r过程算法和基于高斯核函数的相关向量机算法,对C S 号轧辊的非线性时间序列进行回归建模,种模型建模效果如图所示图KWP模型建模效果对比F i g M o d e l i n gr e s u l t so fKWP基于指数函数的W i e n e r过程见图,其预测值与真实值偏差较大,说明基于指数函数的W i e n e r过程算法在轧辊退化建模上的适应性较差其原因主要是指数函数构造过于简单,所能刻画的非线性数据较为有限,无法表征轧辊实际退化路径的复杂非线性特征基于高斯核函数的相关向量机算法预测效果见图高斯核函数相较指数函数能够捕捉更复杂非线性退化路径,获得更好的预测效果但是,基于高斯核函数的相关向量机算法在其退化建模上的鲁棒性较差考虑钢铁热轧生产线关键设备结构极端复杂、随机干扰强,轧辊性能退化过程需要有效地描述设备性能演化的随机特征因此,基于高斯核函数的相关向量机方法无法较好地处理退化路径中的随机性特征,难以描述真实的轧辊退化路径基于高斯核函数的W i e n e r退化过程的预测效果见图,可以看出模型可以较好地预测轧辊的退化趋势,使得预测值基本落在规定的误差范围内,预测效果大大优于基于指数函数的W i e n e r过程,并较基于高斯核函数的相关向量机算法有明显改善为了更好地说明本模型的优越性,统计个模型的RM S E、C D、MA E以进行量化分析,结果如表所示表模型预测误差T a b P r e d i c t i o ne r r o ro fm o d e l s模型RM S EC DMA E基 于 高 斯 核 函 数 的W i e n e r过程模型 指数W i e n e r过程模型 基于高斯核函数的相关向量机模型 由表可知,基于高斯核函数的W i e n e r过程模型,在拟合优度上较指数W i e n e r过程模型提高了 ,较 相 关 向 量 机 模 型 提 高 了,证 明 了KWP退化模型的优越性在计算过程中,相较相关向量机算法的 个R V s,基于高斯核函数的W i e n e r退化过程只需要个R V s,说明本模型的稀疏性远小于一般的相关向量机算法这证明KWP退化模型不仅对非线性和随机性强的数据具有较好的处理能力,并且有助于模型快速收敛,适用于轧辊性能状态的在线监测基于C S 号轧辊数据,采用KWP退化模型预测、批次的RU L,通过计算预测值与实际值之间的误差评价算法的性能,结果如表所示基于KWP退化模型的RU L预测值比较接近真实值,平均误差为 ,且随着时间推移误差不存在明显的扩散现象,所构建W i e n e r过程与核函数 上海交通大学学报第 卷表轧辊可工作转动量预测平均误差统计T a b A v e r a g e e r r o ro f r o l l e rw o r k i n gr o t a t i o np r e d i c t i o n时间序列RU L预测值RU L实际值误差/N o N o N o 相结合可以有效地预测轧辊的RU L以 批次时预测的RU L为例,由式()推导RU L的概率分布如图蓝色区间所示预计下线批次为 批次,包含了可能失效的概率分布范围,推导得RU L在、和 置信区间的取值分别为、和,可为制定轧辊维修和备件计划提供更具体的决策依据图基于 批次的R U L概率分布函数F i g E s t i m a t e dRU La t b a t c h e s基于轧辊可工作转动量健康指标,通过轧辊退化建模与寿命预测研究,可以有效区分轧辊磨损的平稳退化阶段和快速磨损阶段,表征退化过程慢变与快变过程,体现轧辊RU L的不同变化规律结合某钢铁公司 热轧线工程数据,基于KWP退化模型,捕捉到轧辊退化路径的非线性和随机性特征,刻画了轧辊退化的动态规律与不确定性范围,进一步预测轧辊服役过程的RU L值所构建方法通过数据驱动轧辊退化的特征参数,以机理模型刻画轧辊寿命的变化规律,融合两者表征精轧机组健康状态的退化过程和服役寿命特征结语针对轧辊退化路径中存在的随机性和非线性问题,提出了KWP退化模型一方面,W i e n e r过程充分刻画了轧辊工作中因工况复杂、随机干扰强而造成的退化趋势上的强随机性;另一方面,引入核函数将非线性数据映射到高维空间作时间变换,以解决复杂非线性问题并使模型稀疏化,两者结合使算法能够在更广泛的数据集上取得良好效果本文基于某钢铁公司 热轧生产线轧辊退化数据,相较非线性W i e n e r过程和相关向量机算法,所提方法在拟合优度上分别提高 和,RU L预测误差仅有 ,具有较优秀的建模能力与预测效果该结果充分表明所提模型具有较好适用性,可对轧辊运行状况进行分析、预测与决策,有助于轧钢生产线从当前的定期维护过渡到基于条件的维护,亦可为其他类似领域的RU L预测和实时状态监测问题提供参考由于轧钢生产线的产品种类具有多样性,KWP退化模型并未考虑轧钢排产、轧辊材质等具体信息,后续将融合轧辊生产批次信息构建基于多核函数的W i e n e r过程退化模型,开展预测性维护的相关研究参考文献:郑国栋,陈其慎,邢佳韵,等典型国家钢铁产业发展路径与启示J中国国土资源经济,():Z HE NGG u o d o n g,CHE N Q i s h e n,X I NGJ i a y u n,e ta lAp r e d i c t i o nm o d e l f o rw e a ro fw o r k i n gr o l l si nh o tr o l l i n gm i l l sJN a t u r a lR e s o u r c eE c o n o m i c so fC h i n a,():P R AKA S H GA B a y e s i a na p p r o a c ht od e g r a d a t i o nm o d e l i n ga n dr e l i a b i l i t ya s s e s s m e n to f r o l l i n ge l e m e n tb e a r i n gJC o mm u n i c a t i o n si nS t a t i s t i c s T h e o r ya n dM e t h o d s,():刘子英,孙彦广,宋向荣,等热连轧轧机工作辊磨损预报 模 型J热 加工 工艺,():L I U Z i y i n g,S UN Y a n g u a n g,S ONG X i a n g r o n g,e ta l P r e d i c t i o nm o d e l o fw o r kr o l lw e a r f o r h o t t a n d e mr o l l i n gm i l lJH o tW o r k i n gT e c h n o l o g y,():C A I JX,CHE N G X,Z HA O BJ,e ta l S t u d yo nt h ec o r r o s i o n m e c h a n i s m o ft h eo x i d es c a l eo nh o tr o l l e ds t e e l i na na t m o s p h e r i ce n v i r o n m e n tJA n t i C o r r o s i o nM e t h o d sa n dM a t e r i a l s,():L I UZ Y,GUAN Y P,WAN G F QM o d e ld e v e l o p m e n to fw o r kr o l lw e a r i nh o t s t r i pm i l lJM a t e r i a l sS c i e n c ea n dE n g i n e e r i n g,():杨阳,苏小平,赵春磊热轧辊表面疲劳寿命研究J热加工工艺,():YANG Y a n g,S U X i a o p i n g,Z HA O C h u n l e i R e s e a r c ho nf a t i g u el i f eo fh o tr o l ls u r f a c eJH o tW o r k i n gT e c h n o l o g y,():李长生,张晓明,刘相华,等轧制过程轧辊磨损数学模型试验研究J机械工程学报,():第期王汉禹,等:基于核函数 W i e n e r过程的轧辊非线性退化建模与剩余寿命预测 L IC h a n g s h e n g,Z HANGX i a o m i n g,L I U X i a n g h u a,e t a l E x p e r i m e n t a l i n v e s t i g a t i o no fm a t h e m a t i cm o d e l o nr o l l s w e a r i nr o l l i n gJ C h i n e s eJ o u r n a lo fM e c h a n i c a lE n g i n e e r i n g,():宋光义热轧带钢平整机工作辊磨损与工艺参数优化研究D北京:北京科技大学,S ON G G u a n g y i R e s e a r c h o n w o r kr o l l w e a ra n dp r o c e s sp a r a m e t e ro p t i m i z a t i o n f o rh o t s t r i ps k i n p a s sm i l lD B e i j i n g:U n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g yB e i j i n g,潘尔顺,陈震高可靠性产品退化建模研究综述J工业工程与管理,