重要机械零件的剩余寿命预测，特别是对航空航天、大型生产线、船舶、铁路公路运输等重要设备尤为重要.精准的剩余寿命预测为制定适当的预防性维修和更换策略提供保证，既能保证机械设备在高可靠性条件下运行，又能避免由于过频维护造成的浪费.因此机械零件(如齿轮、轴承)的剩余寿命预测越来越得到学者关注.

随着材料制造工艺的进步，重要机械零件往往寿命较长且昂贵，大批量做全寿命试验收集大量状态数据并不现实.在少量试验或使用数据条件下建立适当模型准确预测在役零件剩余寿命是当前研究重点.针对各类设备的剩余寿命预测模型[1-2]大致可分为：专家系统[3]、物理模型、统计和随机过程模型[4-5]、机器学习算法模型[6-7]等.以SVM为代表的机器学习算法模型[8-10]对研究对象退化轨迹适应性强，具有很好的鲁棒性，得到广泛应用.

SVM在小子样预测等问题上表现突出，在时间序列预测、结构可靠性分析、回归分析、经济金融等领域有广泛的应用，机械零件剩余寿命预测研究中也被较频繁使用.Emmanuel等[6]通过分析轴承工作过程中振动信号提取退化信息，建立基于SVM的轴承剩余寿命点预测模型；Loutas等[7]，Caesarendraa等[8]以非线性退化零件或系统为研究对象，将贝叶斯方法应用到SVM建模过程中.

以上文献主要针对剩余寿命点估计取得了很好的效果.在基于可靠性的维护策略制定中，研究对象寿命概率分布情况或其置信区间对维护策略制定也至关重要.

本文在假定机械零件退化指标及对应阀值已知情况下，分别将零件退化指标和与其对应的剩余寿命作为SVM输入变量和输出变量，使用零件全寿命退化数据建立基于SVM机械零件剩余寿命预测模型，并根据SVM概率估计理论推导剩余寿命一定置信水平下的置信区间.

机械零件从开始工作到失效大致要经历3个过程：稳定运行，性能持续退化，最终失效.图 1为某双列滚动轴承整个寿命期间振动均方根值(root mean square,RMS)变化情况，在稳定运行阶段其RMS没有明显变化，初始损伤后其退化现象随着运行时间较明显加剧，当退化量达到某一数值时最终失效.机械零件在初始损伤前零件退化不明显，较难建立数学模型进行分析，因此主要研究初始损伤之后其剩余寿命分布情况.

设某机械零件退化程度的指标为x(t)，t为累计运行时间.不失一般性假设，随着工作时间t的增长，x(t)呈现上升“趋势”，这是由于退化特征提取准确性、工作环境、传感器精度、其他零件等因素影响x(t)不一定表现严格单调.当x(t)等于或大于预设阀值x0时认为该零件已无法满足正常工作要求，应立即更换或进行必要维修.则t时刻的剩余寿命可定义为

其中，inf{}表示下确界.

支持向量机是建立在统计学习理论基础上以结构风险最小化为原则的机器学习算法，在处理小样本问题中表现出良好性质.样本集为

式中：xn为输入变量，可以为标量或l(l>1)维向量；yn为与xn对应的预期输出值，一般为标量，回归函数为

式中：w为与x维数相同的权值向量；b∈R；w·x表示w与x的点积.w和b需通过求解如下优化问题得到：

C和ε分别为惩罚因子和不敏感因子，ξn和ξ*n为松弛因子.根据Karush-Kuhn-Tucker条件，以上问题可转化为凸二次优化问题.

当输入变量x与y之间是非线性关系时，通过引入核函数将x映入高维空间，将低维空间非线性问题转化为高维空间的线性问题.此时二次优化问题只需将式(2)中的xn·x换成K(xn,x)，回归函数为

式中，K(·,·)是满足Mercer条件的核函数，核函数有多种形式，常见的核函数有线性核函数、多项式核函数、Gauss径向基核函数、sigmoid核函数，本文采用Gauss径向基核函数，其形式为

式中，γ是核函数参数.

根据SVM理论，输入x与输出y之间可以表示为

其中，ζ为随机变量，为了得到输出y的置信区间，假设ζ与输入x相互独立.理论上这种假设与实际情况不相符，对于固定的f̂ (x)f^(x) 输出y的残差ζ必然与x存在很大程度的相关性，当假设ζ与x相互独立时，所计算得到的y的置信区间在实际中并不会产生很大误差[10]，能够满足工程需要.基于独立性假设的理论推导复杂程度及计算复杂度有极大简化.

设 通常情况ζ服从零均值的高斯分布或拉普拉斯分布[10]，更多情况下ζ更好地服从拉普拉斯分布.参数为σ的拉普拉斯分布密度函数为

以ζn(n=1,…,N)为样本采用极大似然法估计σ，极大似然函数为

令，得

令，则

所以SVM输出变量y的置信度为p0的置信区间为

若tn(n=1,2,…,N)为传感器采集数据时间点，各时间点对应的机械零件的退化指标为x(tn),剩余寿命为y(tn)，若将x(tn)作为输入变量，将y(tn)作为输出变量，即可训练SVM得到x与y的对应关系模型，进而使用该模型实现机械零件剩余寿命区间估计.

核函数参数γ和惩罚因子C是影响SVM训练效率和推广能力的核心参数，因此γ和C的选择直接影响剩余寿命预测精度.SVM模型推广能力和预测精度是通过其在验证集上的预测表现来衡量的.

设验证集T为

SVM模型推广能力及预测精度可由其在验证集T上的预测均方误差(mean squared error，MSE)度量，

交叉验证法(cross validation，CV)可用来验证训练所得SVM模型推广能力及预测精度.K-fold CV的基本思想是将训练集S随机分成K组(一般采用近似均分方式)，对于同一组γ和C值，K个子集分别做一次验证集，其余的K-1个子集的并集做训练集，得到K个模型及对应的MSEk，取K个MSE值的均值作为与γ和C对应SVM模型的推广能力及预测精度度量.

目前还没有系统准确地确定最佳γ和C的方法，被比较广泛使用的算法有蚁群算法、粒子群算法、网格搜索算法等智能或近似算法.由于蚁群和粒子群算法伴随往往以较大的计算量为代价才能确定较好的(γ,C)，且不能保证算法每次搜索都收敛于同一点，网格搜索算法在处理二维优化问题时稳定性较高，计算量处于可接受范围，本文采用变步长网格搜索算法确定最佳(γ,C)，保证在不明显增加计算量情况下搜索到较精确(γ̂ ,Ĉ )(γ^,C^).

变步长网格搜索步骤：

1) 初步确定γ和C的取值范围.为保证最后搜索到最优(γ,C)，初始范围应适当大些，如γ和C范围都是[10-15,1015].

2) 设定搜索步长.在双对数直角坐标系中设定搜索步长，如0.1.这样就在双对数直角坐标系中生成了一个二维网格，每个交叉点都对应一组(γ,C)取值，对每一组(γ,C)采用交叉验证法计算其对应的MSE值.在双对数直角坐标系中绘制等高线.

3) 缩小γ和C范围.由步骤2)中等高线缩小搜索范围，在双对数直角坐标系中缩小步长，重复步骤2)直至γ和C的范围满足预先设定的精度要求.

本文采用IMS数据库中滚动轴承全寿命试验数据验证所提出基于SVM剩余寿命区间估计模型.轴承的全寿命试验分3组进行，每组4个轴承.试验转速为2000r/min，加载在轴承径向载荷为27kN，采用强制润滑.轴承型号为美国莱克斯诺工业集团ZA-2115双列滚动轴承，加速度传感器安装在轴承座上.每10min采集一次轴承振动信号，每次采集时长为1s，收集20480个加速度值.

针对轴承振动信号退化特征提取方法主要有时域分析、频域分析和时频域分析，结合PAC及ISOMAP等线性或非线性降维方法得到最终的能够反映轴承工作状态的特征量.本文借鉴文献[9]选用相对均方根值(RRMS)作为轴承退化特征.经过平滑处理[9]后第二组试验2，3，4号轴承RRMS如图 2所示.定义RRMS值1.05和1.5分别为轴承衰退期起始阀值和最终失效阀值，当RRMS处于1.05和1.5之间时认为轴承处于衰退期.2，3，4号轴承的起始衰退时间分别为100，103和118h.3个轴承起始衰退后RRMS与剩余寿命对应集分别为S2,S3,S4.

设ε=0.5，初始γ和C范围都是[10-8,108].将S2∪S3作为训练集，采用3.2中改进步长网格搜索确定(γ,C)，得到最佳(γ,C)为

将(γ̂ 1,Ĉ 1)(γ^1,C^1)及S2∪S3代入SVM模型中进行训练得到回归模型1，并对σ做极大似然估计：

用S4作为验证集预测4号轴承剩余寿命，并得到其置信度p0=0.9的置信区间，如图 3a所示.分别以S2∪S4和S3∪S4作为训练集计算(γ̂ 2,Ĉ 2)(γ^2,C^2)和(γ̂ 3,Ĉ 3)(γ^3,C^3)，进而得到回归模型2和回归模型3并估计σ2,σ3，再分别以S3和S2为验证集得到置信度p0=0.9的剩余寿命区间估计结果，如图 3b和3c所示.表 1为3个回归模型参数估计结果.

由图 3可知，基于SVM得到的轴承剩余寿命预测均值与准确值可能存在较大偏差(见图 3b)，导致剩余寿命真实值处在置信区间上下限附近甚至不在置信区间内.导致预测均值误差较大的主要原因是：①轴承的制造加工误差.滚动轴承特别是双列滚动轴承对制造加工误差极敏感，导致其退化轨迹存在较大差异(如图 2所示)；②由于目前还没有准确计算SVM模型参数方法，本文采用应用较广泛的交叉验证法计算SVM模型均方误差作为确定模型精度指标.该方法具有一定随机性，因而不能在理论上保证所确定的参数为全局最佳.

基于支持向量机的剩余寿命点估计能够较准确地预测零件剩余寿命，但当全寿命状态数据较少且零件对制造加工误差敏感时预测结果并不稳定，有时预测值与真实剩余寿命有较大偏差.预测精度随真实剩余寿命变化有明显波动，剩余寿命估计曲线不光滑.本文提出的剩余寿命区间估计方法能够准确估计零件剩余寿命所在置信区间，有效弥补了单纯点估计的不稳定性的不足，能够为制定合理的维修或更换计划提供参考.