测力轮对状态估计计算方法研究

测力轮对状态估计计算方法研究

任愈;陈建政

【摘要】Accurate measurement of wheel/rail force is the premise and basis of monitoring the train operation status.The discontinuous and continuous measurement methods of instrumented wheelset were analysed,the variable characteristics of wheel /rail contact state were studied,the space model for wheel /rail contact force state estimation was established and the wheel /rail force calculation method based on state estimation was proposed.The algorithm was verified by simulation calculations.Moreover,a continuous measurement was achieved in actual track test by using a simple discontinuous instrumented wheelset.%对轮轨间相互作用力的准确测量是监测列车运行安全状况的前提和基础。分析了测力轮对间断测量和连续测量方法，通过研究轮轨接触状态的变化特征，建立了轮轨接触力估计状态空间模型，提出基于状态估计的测力轮对轮轨力计算方法；在此基础上进行了仿真计算验证了算法的有效性，并在线路试验中以简单的间断测量测力轮对实现了轮轨力的连续测量。

【期刊名称】《振动与冲击》

【年(卷),期】2015(000)009

【总页数】5页(P111-115)

【关键词】测力轮对;状态估计;间断测量;连续测量

【作者】任愈;陈建政

【作者单位】西南交通大学牵引动力国家重点实验室，成都 600031;西南交通大

学牵引动力国家重点实验室，成都 600031

【正文语种】中文

【中图分类】U270.7;U216.3;TH721

脱轨是列车运行中最危险的安全事故，也是列车提速中需要研究的关键问题之一。由于轮轨系统的实际运行行为复杂，影响脱轨的因素众多，至今对脱轨机理尚未完全掌握。但可以肯定的是，轮轨力变化是引起脱轨的直接原因。高精度的轮轨力检测，可以为脱轨机理研究、车辆动力学研究和列车安全检测提供重要的基础信息，对轮轨间相互作用力的准确测量是监测列车运行安全状况的前提和基础。

目前，最直接、最准确的轮轨力测量技术就是测力轮对方法，它以车轮作为轮轨力的检测传感器，在车轮辐板有限点处位置贴片组桥，通过测量应变信号实现轮轨接触力的检测。列车运行过程中，由于车轮和轨道间时变、非平稳的轮轨作用力和轮对转动对辐板应变信号的调制作用，使得辐板应变信号较为复杂。不同的组桥方式对应的计算方法也有所不同，从测量方式来看，主要分为间断式测量和连续式测量两种。

间断测量实现方式和计算方法都相对简单，并且可以最大程度的避免横向力、垂向力和轮轨接触点位置三者之间耦合对应变信号的影响。目前在国内测力轮对实测试验中使用较多的是间断式的测力轮对。

间断式测量通过特殊的轮对制作方法，使特定电桥输出只受横向力或垂向力影响。其中，垂向桥输出如图1所示。

垂向桥输出波形中波峰和波谷交替出现，每个波峰和波谷之间轮对相位差为180°；同一辐板上两组垂向桥之间相位差为90°。通过搜索轮对旋转过程中数据的峰值，可以得到间隔角度为90°间断测量的间断点时间坐标，从间断点时间坐标取得间断

点的横向桥和垂向桥输出，代入参数方程进行计算，得到间断点的轮轨横向力和垂向力。

间断测量方法实际上是一种近似的空间采样方法，采样频率为车轮转速的1/4，决定了其采样频率较低，其主要作用是用于评估列车运行中的一些关键参数(脱轨系数、轮重减载率等)，对于轮轨间的高频冲击很难进行有效测量；另外，如果频率较高的冲击发生在靠近实际间断点位置，很容易导致间断点位置判断不准确，从而影响间断测量计算结果精度。

连续测量方法可以有效克服间断测量的弊端，是测力轮对技术的发展趋势。早期连续测力轮对方法是通过仿真和试验选择合适的组桥半径，分别建立横向桥、垂向桥和作用点位置桥，近似的认为轮轨接触力对车轮辐板应变的影响不耦合，其存在的问题是针对不同的轮对，需要设计不同的组桥方案，并且很难找到特别合适的贴片半径。

目前常用的方法是建立特定的应变桥，使电桥输出波形近似正余弦波形或三角波，通过特定方程进行计算[1-2]。如果不能得到良好的输出波形，一般利用三角函数关系式的性质将复杂的非线性方程组简化为线性方程组，通过计算得到轮轨作用力[3-4]。这种方法的问题是应变组桥方法繁琐，并且计算过程相对复杂，对检测精度和可靠性存在一定影响。

一般情况下，采用解方程组方式的测力轮对轮轨接触状态计算方法都是假设采集得到的应变桥输出信号即为应变信号真值。但是实际上，采集到的信号都叠加了相对高频的测试噪声，而噪声的来源包括多个因素，很难通过改进某一具体环节对其进行消除。

常规消除噪声干扰的方法是对测试信号进行选频滤波，以滤波后的数据进行计算。实际上，对于测力轮对这一特殊的应用环境，辐板电桥信号不仅受到轮轨力影响，也被轮对转动所调制，见图2。

假定车轮辐板应变信号经过组桥后输出载波信号为较简单形式，即仅含有一次谐波z(t)=acos(ω0t)，同时将轮轨力描述为傅里叶级数形式：

根据式(1)，任意一个轮轨力信号分量在经过轮对转动调制后都变为含有θn+ω0

和θn-ω0这两个频率分量的信号。可以看出，无论滤波频率选择在什么范围，都会导致有更高频率的轮轨力与车轮转动载波信号的|θn-ω0|分量混合在其中，而更高频率信号中的分量已经被滤波消除了，则滤波后的信号输出为：

从信号测试的角度分析，轮轨接触力的变化是由车辆在行驶过程中受到随机激励引起的，而随机激励在任意不同时刻的取值之间都是无关的。而从微观的角度分析，在车辆行驶阶段的任意间隔较短时间内，轮轨接触状态在时间上存在一定的延续性。也就是说，在正常行驶状态下，只要保证采样频率相对于需要检测的轮轨力的频率足够高，则可以认为当前状态是在前一时刻状态基础上有一些规律性的变化。

状态空间法(State-space Techniques)是现代控制理论中建立在状态变量描述基础上对控制系统分析和综合的方法[5-9]。状态变量是能完全描述系统运动的一组变量。如果系统的外输入为已知，那么由这组变量的现时值就能完全确定系统在未来各时刻的运动状态。状态空间模型是动态时域模型，以隐含着的时间为自变量。状态空间模型包括两个模型：一是状态转移模型，反映动态系统在输入变量作用下在某时刻所转移到的状态，如式(3)；二是量测模型，描述输出变量与状态变量和输

入变量间的变换关系，如式(4)。两者构成了状态空间模型。

通过分析辐板应变输出与轮轨接触状态的关系，可以认为导致轮对应变桥输出改变的作用因素包括横向力Q、垂向力P、作用点位置X和轮对旋转角度θ，即可设轮轨接触状态向量为x=[Q,P,X,θ]T。设uk为状态的激励变化量，wk为状态转移噪声。在常规条件下，可以近似地假设轮轨接触状态变化为：

确定状态向量后，需确定各成员变量的递推规律进行统计和推导。

图3为在西南交通大学滚振试验台上进行的测力轮对滚动试验中的轮轨力输出，