摘要
为全面科学地评估高速磁浮列车悬浮系统性能,基于模糊层次分析法设计了悬浮系统综合评估方法。针对悬浮系统运行规律和能力特点,提出了表征悬浮系统性能的评价指标体系;采用模糊层次分析法确定高速磁浮列车悬浮系统性能指标的权重,利用指数标度法和梯形分布法建立了悬浮系统多层模糊综合评估模型;开发了高速磁浮列车悬浮系统性能综合评估系统,并利用时速600 km高速磁浮样车在上海1.5 km高速磁浮试验线开展动态运行测试。结果表明基于模糊层次分析的悬浮系统综合评估方法能够对悬浮系统性能做出综合评估,实现了评估结果的可量化与可视化。
Abstract
In order to comprehensively and scientifically evaluate the performance of suspension system in high-speed maglev train, a comprehensive evaluation method was designed based on fuzzy analytic hierarchy process. According to the operation law and capability characteristics of suspension system, an evaluation index system was proposed to characterize the performance of suspension system. The fuzzy analytic hierarchy process was used to determine the weight of the performance index of the suspension system, and a multi-layer fuzzy comprehensive evaluation model of the suspension system was established based on exponential scaling method and trapezoidal distribution method. A comprehensive performance evaluation system of suspension system based on operation data was developed, and the dynamic operation test was carried out on the Shanghai 1.5 km test line by using a 600 km/h high-speed maglev train. Results show that the comprehensive evaluation method of suspension system based on fuzzy analytic hierarchy process can comprehensively evaluate the performance of suspension system and realize the quantification and visualization of the evaluation results.
高速磁浮技术是引发轨道交通产业变革的前瞻性、颠覆性技术,是世界科技争夺制高点,德国、日本、美国等国家均开展了相关技术研究,不断争夺和刷新速度纪录[1-3]。我国从“十五”开始支持高速磁浮技术发展,历经四个五年计划的科研攻关,实现了高速磁浮技术的国产化和自主创新。2021年7月20日,具有我国完全自主知识产权的时速600 km高速磁浮列车在青岛正式下线,标志着我国已掌握了常导电磁型高速磁浮的核心关键技术,已经具备自主化与产业化能力[4]。悬浮系统是时速600 km高速磁浮列车的“轮子”,列车利用悬浮系统产生的方向向上的电磁力来支撑列车车体,从而使列车无接触地悬浮于轨道上方。因此,悬浮系统性能是影响高速磁浮列车安全平稳运行的关键[5-6]。目前工程中主要采用基于预设阈值的性能评估方法,通过判断悬浮间隙是否超出预设阈值直接对悬浮系统性能做出评估。但是,该方法评估标准过于粗糙,只能对悬浮系统性能进行定性评估,无法量化。为实现磁浮列车悬浮系统性能的量化评估,倪菲等[7]基于云理论以及磁浮列车的运行数据,研究了数据驱动下的控制回路性能评估方法,利用多变量控制系统特征实现了评估性能的可视化;Xu等[8]提出了一种用于稳定性能监测的实时性能指标,并将监测到的稳定性结果进行性能分级;陈琛等[9-10]分析了随机不平顺下磁浮车辆的动力学特性,通过定义不同悬浮力条件下柔性轨道的振动固有频率和模态矩阵,提出了一种磁浮车辆垂向悬浮稳定性分析方法。上述方法虽然很好地解决了悬浮系统性能的量化评估问题,但是未形成全面完整的评估指标体系与评估标准,不能给出明确的性能评估等级,无法指导现场调试和后期维护,并且尚未实现工程应用。
模糊层次分析法能够综合考虑影响悬浮系统性能的关键因素,量化评估过程中的各种模糊信息,并合理给出表征悬浮系统性能的结果等级,全面、科学地评估悬浮系统性能。因此,针对悬浮系统运行规律和能力特点,本文提出了表征悬浮系统性能的评价指标,形成了全面完整的评估指标体系;采用指数标度法和梯形分布法建立了悬浮系统多层模糊综合评估模型,基于模糊层次分析法设计了悬浮系统综合评估方法;利用运行数据开发了高速磁浮列车悬浮系统性能综合评估系统,并利用时速600 km高速磁浮样车在上海1.5 km试验线开展了动态运行测试验证。
1 悬浮系统组成和工作原理
悬浮系统是高速磁浮列车的“轮子”,不同于传统轮轨列车依靠车轮支撑车体,高速磁浮列车通过悬浮系统来支撑整个车体。悬浮电磁铁通电后可以产生可控电磁场,悬浮电磁铁与轨道上的长定子铁芯相互吸引,方向垂直向上的电磁吸力将整个车体向上吸起从而使列车悬浮于轨道之上[11-12]。悬浮电磁铁与轨道之间的最终悬浮间隙被控制在10~12 mm。
高速磁浮列车悬浮系统由悬浮电磁铁、悬浮控制器以及悬浮传感器构成。悬浮传感器安装在悬浮电磁铁上,主要用于检测悬浮电磁铁与长定子轨道的间隙以及悬浮电磁铁的垂向加速度[13-14]。悬浮控制器利用悬浮传感器的量测信息,通过内置的悬浮控制算法计算得到控制输出量,并驱动悬浮电磁铁产生悬浮力使列车稳定悬浮。悬浮电磁铁是悬浮系统的执行机构,在控制电压的作用下产生动作,从而保证悬浮电磁铁始终工作在给定间隙[15-17]。高速磁浮列车悬浮系统控制框图如图1所示。

图1高速磁浮列车悬浮系统控制框图
Fig.1Control block diagram of suspension system of high-speed maglev train
2 悬浮系统的性能指标与性能评估
悬浮系统结构复杂,呈现非线性及强动态耦合特性,是一个多输入多输出、时变的不稳定系统。同时,悬浮系统运行过程中又面临轨道不平顺、线路参数时变、网络传输时延等多源复合扰动,电-磁-热-应力-流体等多场耦合效应复杂多变。悬浮系统为了能够快速跟踪线路轨道变化,始终处于动态调整状态。因此,为了保障和实现高速磁浮列车安全平稳地运行,悬浮系统必须具备良好的动态性能。具体为:
1)具备快速跟踪线路轨道变化,保证悬浮电磁铁始终工作在给定间隙附近的能力。
2)具备有效抑制轨道不平顺、线路参数时变、网络传输时延等多源复合扰动的不利影响,保证悬浮系统在多源复合扰动下能够快速恢复到稳定状态的能力。
3)具备减缓悬浮电磁铁的垂向振动加速度,保证悬浮系统运行总体平稳的能力。
综上所述,为了对悬浮系统性能做出全面、科学、准确的评价,全面掌握和反映悬浮系统的运行状态,必须针对悬浮系统以上三方面能力做出综合评估,提出适合高速磁浮列车悬浮系统的评价指标体系,利用不同维度的性能指标来表征悬浮系统性能的好坏,并作为评价悬浮系统性能优劣的依据。
2.1 确定性指标
确定性指标是用来表征悬浮系统跟踪线路轨道变化,保证悬浮电磁铁工作在给定间隙能力大小的指标。脉冲响应曲线指标和误差积分指标可以作为评价悬浮系统跟踪能力和控制精度好坏的指标[18]。脉冲响应曲线能够直观地反映出悬浮系统跟踪参考输入能力的好坏,而误差积分指标通过计算悬浮系统期望输出与实际输出之间误差的积分来反映系统稳态时的控制精度。
脉冲响应曲线指标主要包含峰值时间以及衰减速率。脉冲响应曲线的峰值时间越小,衰减速率越快,说明悬浮系统的跟踪能力越强。误差积分指标主要有绝对误差积分(integral of absolute error,IAE)、平方误差积分(integral of squared error,ISE)、时间乘绝对误差积分(integral of time multiplied by absolute error,ITAE)、时间乘平方误差积分(integral of time multiplied by squared error,ITSE)。误差积分指标值越小,说明悬浮系统的控制精度越高。
因此,通过绘制脉冲响应曲线获得峰值时间和衰减速率,利用误差积分准则计算得到误差积分指标,能够实现对悬浮系统的跟踪能力和控制精度的评估。
2.2 随机性指标
随机性指标是用来表征悬浮系统抑制轨道不平顺、线路参数时变、网络传输时延等多源复合扰动不利影响,保证悬浮系统在多源复合扰动下快速恢复到稳定状态能力大小的指标。最小方差指标可以作为评价悬浮系统抵抗外界干扰能力的指标[19]。
根据Harris对最小方差指标的定义,悬浮系统的最小方差指标如式(1)所示。
(1)
式中:为最小方差控制器输出误差的方差,可以通过计算最小方差不变项来获得;为实际悬浮控制器输出误差的方差。
由于大于,因此ηHarris的取值范围在0到1之间,ηHarris越接近1,悬浮系统的随机性指标越好,抵抗外界干扰的能力越强;ηHarris越接近0,悬浮系统的随机性指标越差,抵抗外界干扰的能力越弱。
为获得悬浮系统的最小方差指标,利用悬浮系统的输出数据建立自回归移动平均模型,采用长除法将其转化为无限阶的滑动平均模型(moving average model,MA),如式(2)所示。
(2)
式中,yt为悬浮系统输出,fi为常数项,at为平均值为零的白噪声。式(2)左右两边分别乘以at,at-1,···,at-d+1,分别取其数学期望,即可得到式(3)。
(3)
其中,为白噪声序列at的方差,γya(i)为yt与at-i的协方差。因此,输出方差的不变项如式(4)所示。
(4)
综上,悬浮系统的最小方差指标表达式如式(5)所示,利用输入输出数据可以计算得到悬浮系统的最小方差指标。
(5)
2.3 平稳性指标
平稳性指标是用来表征悬浮系统减缓悬浮电磁铁振动加速度,保证悬浮系统运行总体平稳能力大小的指标。振动加速度的欧式距离指标可以作为评价悬浮系统平稳性的指标[20]。
平稳性主要取决于悬浮电磁铁的振动加速度,振动加速度值越小,且振动加速度的变化率越小,振动加速度的欧式距离越小,悬浮系统的平稳性越高。
因此,通过对悬浮电磁铁的振动加速度进行特征提取,构造多维特征分布空间,对其特征进行标准化处理,根据超球体分布特点,获取振动加速度的欧式距离,可实现对悬浮系统平稳性的评估。
悬浮电磁铁振动加速度的数据特征主要有标准偏差、均方根、平均值、极差、方差,为了使针对平稳性的评估结果更可靠、准确,需要对上述特征进行标准化处理。但是,在标准化之前,有必要先明确特征分布,只有特征满足正态分布时,才能对其进行标准化。如果特征值不满足正态分布,利用Box-Cox将特征值的非正态分布转换为正态分布。特征数据全部转换为满足正态分布的数据后,利用z-score标准化方法进行标准化处理,如式(6)所示。
(6)
式中,X′k为第k个特征值标准化后的值,Xk为第k个特征值,为第k个特征平均值。
标准化后的数据特征分布在超球体中,根据超球体分布特点,采用到球心的欧式距离作为振动加速度的欧式距离(距离),如式(7)所示。
(7)
3 悬浮系统的性能指标与性能评估
影响高速磁浮列车悬浮系统性能的因素较多且相互耦合,导致性能评估结果比较模糊。为了对悬浮系统性能做出全面、科学准确的评价,全面掌握和反映悬浮系统的运行状态,本文采用模糊层次分析法确定高速磁浮列车悬浮系统性能指标的权重,基于指数标度法和梯形分布法建立悬浮系统多层模糊综合评估模型,通过计算各指标的隶属度,依据模糊理论最大隶属度法则确定高速磁浮列车悬浮系统性能等级。基于模糊层次分析法的综合评估流程图如图2所示。

图2基于模糊层次分析法的综合评估流程图
Fig.2Flow chart of comprehensive evaluation based on fuzzy analytic hierarchy process
3.1 确定评估指标与评估集
U=[U1 U2 U3]=[随机性指标 确定性指标 平稳性指标]
第二层因素集为:
U1=[u11]=[最小方差指标]
U2=[u21 u22 u23 u24 u25 u26]=[峰值时间 衰减速率 IAE ISE ITAE ITSE]
U3=[u31]=[振动加速度τ距离]
评估集为B={优秀,良好,一般,差}。
3.2 构建判断矩阵
由于确定性指标表征悬浮系统的跟踪能力和控制精度,随机性指标表征悬浮系统抵抗外界干扰的能力,平稳性指标表征悬浮系统的平稳性。因此,综合考虑U1,U2,U3对U的重要性,判断矩阵R1为:
将矩阵R1进行归一化处理:
通过求取R1矩阵的最大特征值,完成判断矩阵R1的一致性检验,确保构建的判断矩阵不存在逻辑问题。U1、U2、U3对U的权重为:

矩阵R2按照矩阵R1方式进行归一化处理:
通过求取R2矩阵的最大特征值,完成判断矩阵R2的一致性检验,u21、u22、u23、u24、u25、u26对U2的权重为:

3.3 建立隶属度函数
隶属度函数的建立是模糊综合评估法中至关重要的一步。确定隶属度函数的常用方法有:模糊分布法、德尔菲法、专家评分法等[27-30]。通过对悬浮系统性能指标的分析,隶属度函数的类型属于极小型。因此选择模糊分布法,并选择梯形分布来建立高速磁浮列车悬浮系统性能评估的隶属度函数,具有计算简便、分布合理的优点。
梯形分布隶属度函数示意如图3所示,x轴表示需要评判的参数值,y轴表示参数值的隶属度,隶属度在0到1之间。
悬浮系统性能“优秀”的隶属度函数如下所示:
(8)
悬浮系统性能“良好”“一般”的隶属度函数如下所示:
(9)
悬浮系统性能“差”的隶属度函数如下所示:
(10)
其中:y为隶属度值,x为各项指标值,a、b、c、d为各评估指标的参考点。
为了确定各评估指标的参考点,首先,通过分析悬浮系统可用数据集,特别是与模糊评价相关的悬浮间隙、加速度数据,观察数据的分布和特征,根据数据的范围和分布来选择拐点横坐标和斜率值。其次,结合项目组在磁浮列车悬浮系统设计、研制与调试的经验,以及磁浮运营线积累形成的专家经验,对得到参数组合进行调整和修正,最终得到如表1所示的各评估指标的参考点设置。

图3梯形分布隶属度函数示意图
Fig.3Schematic diagram of membership function of trapezoidal distribution
4 悬浮系统运行过程数据采集与处理
为实现基于模糊层次分析的悬浮系统综合评估,需要利用悬浮系统的输入输出数据。对于高速磁浮列车悬浮系统而言,通过车载数据采集和存储系统能够实时获取悬浮间隙、加速度、速度以及电流等运行过程数据。但是,获取的运行数据质量普遍不高,存在离群点污染严重等诸多问题,需要对数据做相应的处理,使数据满足时序分析要求。
表1隶属度函数的参考点
Tab.1 Reference point of membership function

4.1 数据采集
一节完整的时速600 km高速磁浮样车在上海同济大学嘉定校区1.5 km高速磁浮试验线进行动态运行测试。悬浮间隙、加速度、速度以及电流等运行过程数据可以通过周立功CAN USB-2E-U采集卡发送到车载数据存储系统。
采集的悬浮系统运行数据存在质量不高,离群点污染严重等诸多问题,因此,进行时间序列建模之前,需要针对数据进行预处理,具体包括异常值处理、平滑处理、中心化处理等。
4.2 数据处理
1)异常值处理。高速磁浮线路轨道为应对温升、沉降,相邻轨道梁间留有宽度约为86 mm的空隙。当悬浮传感器进入接缝区域,采集到的数据为输出最大饱和值20 mm。由于其不是真实测量值,悬浮控制系统内部也会将其作为异常值进行剔除,使用上个周期的悬浮间隙测量值来进行控制算法计算。如图4所示,对于分布在特定区域或限定范围之外的离群点数据,需要剔除异常值。本文采用拉依达准则判断量测数据是否含有粗大误差,并用量测数据的平均值替代离群点数据。如图5所示,悬浮间隙数据经异常值处理后,原始悬浮间隙数据中出现的不合理的“尖刺”现象得以消除。

图4磁浮轨道接缝图
Fig.4Maglev track joint diagram

图5悬浮间隙的数据处理效果图
Fig.5Data processing effect diagram of suspension gap
2)平滑处理。采集到的初始运行数据包含一定的噪声,直接影响拟合模型的准确性,采用移动平均法来对运行数据序列进行平滑处理,能够消除噪声的不利影响,也能减少瞬时随机波动的影响,从而更好地显示运行数据的规律特性。
3)中心化处理。当运行数据平均值不为零时,建立的时间序列模型为非中心化模型,需要进行中心化处理。通过对非中心化模型做一个简单的位移就可以转化为中心化模型。中心化模型将少一个参数且运算更方便。中心化运算不会影响序列值之间的相关关系,如图6所示,悬浮间隙数据经中心化处理后,其平均值由9.5变为0,并保持了原有变化规律。
5 悬浮系统性能的模糊综合评估
为实现评估结果的可量化与可视化,本文利用采集的悬浮间隙、加速度等实际运行数据,开发高速磁浮列车悬浮系统性能评估系统。该评估系统基于确定性指标、随机性指标以及平稳性指标,采用模糊层次分析法对悬浮系统性能进行综合评估,能够直接给出高速磁浮列车悬浮系统的性能指标与性能等级,非常便于指导现场调试和后期维护。高速磁浮列车悬浮系统模糊综合评估流程图如图6所示。

图6高速磁浮列车悬浮系统模糊综合评估流程图
Fig.6Flow chart of fuzzy comprehensive evaluation of suspension system of high-speed maglev train
采集时速600 km高速磁浮样车在上海同济大学嘉定校区1.5 km高速磁浮试验线在0~60 km/h速度区间全线路运行时的悬浮间隙、加速度数据,通过对运行数据进行处理和分析,基于运行数据计算得到高速磁浮列车悬浮系统全线路运行时的确定性指标、随机性指标以及平稳性指标,具体为:

随机性、确定性和平稳性的评判矩阵分别为:

对第二层进行模糊综合评判,即B2=ω2×A2,可得B2为:

再对第一层进行模糊综合评判,已知权重为ω1,且总评判矩阵为B=[A1 B2 A3]T,最终的模糊综合评估结果B为:

根据最大隶属度原则,高速磁浮列车悬浮系统全线路运行时的性能评估结果为“优秀”,并且“优秀”的隶属度为0.494 7。
基于相关APP Designer设计并开发的高速磁浮列车悬浮系统性能评估系统如图7所示,实现了评估结果的可量化与可视化。

图7全线路的评估结果
Fig.7Evaluation results of the whole line
6 结论
本文设计的基于模糊层次分析的高速磁浮列车悬浮系统综合评估方法能够综合考虑影响悬浮系统性能的关键因素,量化评估过程中的各种模糊信息,并合理给出表征悬浮系统性能的结果等级,全面、科学地评估悬浮系统性能。通过在高速磁浮试验线进行动态运行测试,基于本方法得到的性能评估结果能够有效反映悬浮系统性能和运行状态,验证了方法的有效性。
1)针对悬浮系统运行规律和能力特点,提出了适合高速磁浮列车悬浮系统的评价指标体系,采用确定性指标、随机性指标以及平稳性指标来表征悬浮系统性能的好坏,并作为评价悬浮系统性能优劣的依据。
2)采用模糊层次分析法确定了高速磁浮列车悬浮系统性能指标的权重,基于指数标度法和梯形分布法建立了悬浮系统多层模糊综合评估模型,通过计算各指标的隶属度,依据模糊理论最大隶属度法则确定了高速磁浮列车悬浮系统性能等级。
3)通过采集和处理悬浮间隙、加速度等运行数据,开发了基于数据驱动的高速磁浮列车悬浮系统性能综合评估系统,并利用时速600 km高速磁浮样车在上海1.5 km试验线开展动态运行测试,验证了方法的有效性和可行性。