1、第 卷第 期 年 月兵工学报 ：基于智能优化相关向量机的综合传动变工况直驶传递效率预测张金豹，盖江涛，安媛媛，桂林，朱炳先，邹天刚（中国北方车辆研究所，北京）摘要：传递效率是表征履带车辆综合传动装置性能的重要指标之一，但综合传动装置内部复杂的动力传递路径和外部工作环境导致其传递效率预测模型建模困难，无法保证传递效率的准确计算。基于上述分析，提出一种基于综合传动装置直驶模式下多种工况传递效率数据驱动的传递效率预测方法。对 台综合传动装置不同工况下的传递效率数据进行统计分析并得到统计特征；通过统计特征选择智能优化相关向量机进行传递效率预测建模，进而实现待测综合传动装置在变工况直驶模式下的传递效率预
2、测。结果表明，变工况下测试样本传递效率预测结果的平均绝对误差和均方根误差分别控制在.和.以内，验证了该方法的可行性和正确性。关键词：综合传动；直驶；传递效率预测；变工况；相关向量机 中图分类号：.文献标志码：文章编号：（）收稿日期：基金项目：武器装备预先研究项目（）；基础产品创新科研项目（），（，）：，（），（）（）.，：；第 期基于智能优化相关向量机的综合传动变工况直驶传递效率预测；引言履带车辆综合传动装置内部包括大量齿轮、轴和轴承等关键零部件，通过不同的功率传递路径及配合来实现履带车辆的直驶、转向和制动等动作。为满足实际工况如爬坡、跨越壕沟和高速行驶等需求，需要设置多种工作模式及对应多个挡
3、位，由此引起的功率损失包括旋转零部件搅油损失、摩擦生热损失以及润滑系统的功率损失等。通过对综合传动装置传递效率进行预测评估，可以找到功率损耗的主要影响因素并加以改进，达到综合传动装置质量持续提升的目的。由于综合传动装置的个体差异性，并受外部复杂工作环境如路面、气压等和内部因素如油液高度、温度等的影响，其传递效率为一随机变量，并且不同工况下的传递效率并不一定服从同一概率分布，且传递效率统计特征呈现时变性、异方差性和“尖峰厚尾”等特性。目前工程应用中主要使用的概率分布有正态分布、对数正态分布以及威布尔分布等。通过正态分布的 图，可以判断传递效率是否服从正态分布，以及其是否具有“尖峰厚尾”的特性。与
4、指数分布相比，厚尾分布的尾部更缓慢的趋向于，能够对传递效率测试中出现的极值进行统计，因此可以在预测建模过程中更全面的考虑传递效率数据的统计特性，不仅能保证预测的精度，还能为传递效率的提升提供技术支持。由于传统概率统计预测模型基于正态分布假设，因此传统概率统计预测模型将忽略传递效率尾部特征，反映的是综合传动装置一致性的范围，无法体现综合传动装置传递效率的差异性。除概率统计预测外，基于数据驱动的机器学习方法如人工神经网络、支持向量机（）等能够通过非线性映射实现预测，其中 通过引入移位损失函数实现与厚尾分布的良好结合。相关向量机作为 的一种改进形式，在保证 良好性能的基础上还能够实现区间预测，且已在
5、疲劳裂纹和轴承寿命预测等领域进行应用。为提升 的预测性能，途径之一就是核函数的自适应选择，通过多种核函数的组合及智能优化算法，实现自适应核函数的 参数估计。基于上述分析，本文首先对综合传动装置各挡位传递效率的统计特性进行分析，进而提出一种基于智能优化 的综合传动装置变工况直驶传递效率预测方法，为复杂工况下综合传动装置的效率优化和提升奠定基础。综合传动传递效率特性分析.综合传动直驶模式传递效率测试综合传动装置是机电液一体化的复杂系统，如图 所示，由前传动总成、液力变矩器、行星变速机构、左 右汇流排、液压操纵系统等 多种部件组成。其工作原理如下：动力由前传动输入，一路功率流（直驶）经液力变矩器传至
6、行星变速机构，再经主轴和液力减速器传至两侧汇流排。另一路功率流（转向）由前传动带动联体泵马达，再由马达传至两侧汇流排太阳轮。直驶时，零轴（汇流排太阳轮）制动，动力经汇流排输出至两侧主动轮。图 综合传动装置原理图 .变工况传递效率数据分析根据国家军用标准装甲车辆综合传动装置台架试验方法（），共对 台综合传动装置进行变工况传递效率测试。传递效率测试试验台如图 所示，综合传动装置输入端通过弹性联轴器与电机相连，两侧输出端通过弹性联轴器分别与负载连接，输入端和输出端均安装转速和扭矩传感器。测试试验装置中，电机用于模拟发动机输入，以下全功率加载指的是实际工作过程中发动机所能提供的最大功率输入。根据项目试
7、验测试大纲，全功率 加载测试兵 工 学 报第 卷图 综合传动装置传递效率测试试验台 时直驶前进挡位 个，分别为 ；倒挡挡位 个，分别为 和；每个挡位下均设有 个测试转速点，分别为、和 ，测试方案如表 所示；较全功率 加载测试，全功率测试时不包括 挡和 挡在转速 和 下的测试，以及 和 挡位。表 多挡位试验方案设计 速度（）挡位 在本文中，定义每一个挡位速度负载的组合称为一个工况点。对于每一台测试样机，按照以下顺序进行全功率 加载和全功率加载下的传递效率测试：）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）
8、（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）。综合传动装置传递效率计算公式为 （）式中：、和 分别为综合传动装置的输入扭矩和两侧输出扭矩；、和 分别为综合传动装置的输入转速和两侧输出转速。各工况点的传递效率测量结果如图 所示，在、和 挡位的传递效率随着转速的升高而升高，而在 挡位则随着转速的升高而降低，这是因为在、和 挡位车辆传递装置处于液力机械工作模式，在 挡位则处于纯机械工作模式。绘制不同挡位全功率加载下 时传递效率的正态分布 图，如图 所示。从图 中可以看出，在每个 图中，在中间数据部分线性拟合良好，但两端数据均呈现不同程度的偏移，这说明传递效率的概率统计为“
9、尖峰厚尾”。如图 所示，计算 个工况点下传递效率的峭度值，有 个工况点的峭度值大于，其中一个略小于，进一步验证了传递效率“尖峰厚尾”的统计特征。在基于机器学习的预测过程中，输入变量中挡位、转速为固定离散变量，输入扭矩与传递效率为随机变量，受到综合传动装置设计质量和运行状态的影响。如图 所示，输入扭矩与传递效率的概率分布并不一致，虽然在趋势上呈现粗糙的线性关系，但离散点依旧很多，无法通过简单的线性模型进行传递效率预测。因此本文关注的重点就是如何解决这些在概率分布“厚尾”处的传递效率预测。基于 的综合传动变工况直驶传递效率预测.相关向量机对于给定的训练样本集，其均值输出 第 期基于智能优化相关向量
10、机的综合传动变工况直驶传递效率预测图 各工况下传递效率统计规律 定义为（；）（）（）（）式中：（，）为加权值；（），（），（），（）为核函数矩阵，其中（）（，）为核函数。进而求得目标预测函数为（；），（，）（）式中：为样本噪声，进而有（，）（，），），则训练样本集的似然函数为（，）（）（）式中：（，）为 真 实 目 标 值；（），（），（）为 （）的预设矩阵，其中（），（，），（，），（，）。如果直接通过最大化似然函数来估计权重向量 会导致模型的过学习。因此在 模型参数估计中，使用贝叶斯框架求解极大似然函数，得到超参数 和 的估计值 和，其中 为 方差的倒数。对于一个新的测试样本，就可得到相应
11、输出的概率分布为（，）（，）（，）（）该测试样本的预测均值和方差分别为（）（）（）（）（）其中权重向量 的后验均值和协方差分别为 （）（）（）式中：（，）。基于贝叶斯框架的相关向量机学习方法可以提供概率式的预测，能够较好地获得预测值及其概率分布。由于单个核函数的模型表达能力有限，为获得更为准确的预测结果，本文采用由高斯核函数和 核函数线性组合的核函数，其中高斯核函数表达式为（，）()（）核函数表达式为（，）（）（）得到的组合核函数为（，）（，）（，）（）式中：为核函数带宽；和 为表达式系数；和 为组合核函数的加权系数。由于智能优化算法具有自组织自适应等优点，能够实现对非线性方程多个参数的同时求
12、解，因此本文采用粒子群优化（）算法对模型的多个参数进行估计。.预测流程基于 模型的综合传动装置变工况直驶传递效率预测具体步骤如下：）按照测试顺序进行各工况点下的综合传动装置直驶传递效率采集。将综合传动装置的挡位、转速、工况和负载作为 的输入，将传递效率作兵 工 学 报第 卷图 各挡位在 ，全功率加载下的传递效率 图 为 的输出，代入 进行训练，如表 所示。其中挡位、转速和工况均为指定的离散变量，而负载输入扭矩则为与综合传动装置本身运行状态相关的随机变量，并与传动效率一一对应。本次测试共有 台样机，随机选取 台样机作为测试样本，其他作为训练样本。表 智能优化 的输入和输出 输入测试变量设定值输出
13、挡位 ；转速、传递工况液力（、和 ）、机械（）效率负载（输入扭矩）全功率、全功率 ）明确 预测模型中需评估的参数，包括高斯核函数宽度，表达式系数 和，以及组合核函数的加权系数。在本文进行的多次 参数估计中，求解得到的目标函数值接近，但估计的参数波动性较大，因此本文后续的计算结果不再给出估计参数的具体值。）初始化参数，包括粒子群规模，最大迭代次数，随机初始化粒子位置，即 预测模型中所有优化参数对应一个初始粒子位置。将所有初始化粒子代入 预测模型，通过贝叶斯框架求解极大似然函数得到其他参数。）采用 折交叉验证的最小均方根误差作为优化目标函数 第 期基于智能优化相关向量机的综合传动变工况直驶传递效率
14、预测图 各工况下传递效率的峭度值 （）（）图 各挡位在 ，全功率加载下的传递效率和输入扭矩的概率关系图 式中：为试验值；为拟合值。本文中 。通过跟踪当前最优的粒子来不断搜索和迭代，最终得到所求的参数值。综合传动装置变工况传递效率预测结果与分析 图 为 个训练样本各工况点及相应传递效率代入 进行训练的目标函数收敛过程。文中采用均方根误差和平均绝对误差对智能优化 预测模型的拟合精度进行评价，如图 所示，其中平均绝对误差表达式如下：（）式中：为试验值；为拟合值。从图 中可以看出，训练样本的 和 分别控制在.和.以内。对随机选取的 个测试样本进行变工况传递效率预测，结果如图 所示。从图 可以看出，各个
15、工兵 工 学 报第 卷图 粒子群优化收敛过程 图 训练样本的拟合误差和测试样本的预测误差 况点的传递效率预测情况良好，其 和 均小于训练样本，如图 所示。另外，对应预测区间，的间距小，说明了预测结果稳定。进一步对一挡和六挡在转速区间 依次升高，负载依次为、和 的传递效率进行了预测，如图 所示。可以看出，预测趋势与试验测试结果基本符合，说明本图 测试样本变工况直驶传递效率预测结果 第 期基于智能优化相关向量机的综合传动变工况直驶传递效率预测文研究的有益性。图 不同转速及对应输入扭矩下的传递效率趋势预测 结论本文通过智能优化 对综合传动变工况直驶传递效率预测进行建模和预测。结果表明：）从不同负载不
16、同转速不同挡位下传递效率的统计特征出发，得到传递效率服从厚尾分布，如果采用传统概率预测模型将无法对传递效率进行准确的预测。）基于粒子群优化的双核，能够将变工况直驶模式下测试样本传递效率预测值的 和 分别控制在.和.以内，并给出预测区间 ，。）通过粒子群优化双核 提取不同转速及对应输入扭矩下传递效率的发展趋势，能够对没有测量的工况点进行定量分析。参考文献（）徐保荣，张金豹，姚李刚，等 液力机械综合传动装置低温阻力矩特性研究与验证 兵工学报：，：（）李义军，毕小平，马宁，等 车辆传动装置产热建模与计算分析 中国工程机械学报，（）：，（）：（）王斌，宁斌，陈辛波，等 齿轮传动搅油功率损失的研究进展
17、机械工程学报，（）：，（）：（），（）：，（）：，（）：，（）：张金豹，邹天刚，高秀才，等 基于中位秩加权的车辆传动装置直驶 模 式 传 递 效 率 性 能 评 估 兵 工 学 报，（）：，（）：（）：，（）：，：，（）：杜刚，赵冬梅，刘鑫，等 基于重尾分布的风电功率波动特性概率分布 电力自动化设备，（）：，（）：，（）郑建锋，王应明 基于 神经网络的效率置信区间预测模型研究 计算机工程与应用，（）：，（）：（），（）：，兵 工 学 报第 卷 ，（）：，（）：，：范强飞，廖爱华，丁亚琦 基于 和 的滚动轴承剩余寿命预测 上海工程技术大学学报，（）：，（）：（）刘月峰，赵光权，彭喜元 多核相关向量机优化模型的锂电池剩余寿命预测方法 电子学报，（）：，（）：（）雷亚国，陈吴，李乃鹏，等 自适应多核组合相关向量机预测方法及其在机械设备剩余寿命预测中的应用 机械工程学报，（）：，（）：（），：张金豹 基于全寿命数据的滚动轴承运行状态评估和剩余寿命预测 哈尔滨：哈尔滨工业大学，：，（）作者简介：张金豹（），男，博士。：盖江涛（），通信作者，男，研究员。：安媛媛（），女，博士，工程师。：桂林（），男，研究员。：朱炳先（），男，工程师。：邹天刚（），男，研究员。：
