CN111114562B - 一种机车车辆及其加权参数粘着控制方法 - Google Patents

Abstract

本发明公开了一种机车车辆及其加权参数粘着控制方法，根据机车运行状态获取牵引电机转矩给定值和加权粘着控制参考值，生成牵引电机转矩给定控制值；计算加权粘着控制反馈值；将加权粘着控制参考值与加权粘着控制反馈值在加权粘着控制PI闭环控制器中计算，得出加权粘着控制值；依据加权粘着控制值和牵引电机转矩给定控制值二者中较小者控制牵引电机转矩；从T0时刻开始，按m*G(S)对牵引电机转矩给定值进行滤波，直到牵引电机转矩给定控制值恢复到牵引电机转矩最大值的a％时，按G(S)/n对牵引电机转矩给定值进行滤波处理。本发明最大限度发挥机车车辆的粘着力，有效防止牵引空转或制动滑行，使机车车辆在当前的轨面状态下发挥最大的牵引力或制动力。

Description

一种机车车辆及其加权参数粘着控制方法

技术领域

本发明涉及铁路机车车辆电力牵引与传动控制技术领域，尤其涉及一种机车车辆及其加权参数粘着控制方法。

背景技术

高速铁路、重载运输是铁路现代化的重要标志，为了实现高速、重载的目标，最大限度发挥机车车辆牵引与制动性能至关重要，而轮轨粘着利用状况直接影响机车车辆牵引与制动性能的发挥，当轮对产生的轮周牵引力或制动力大于轮轨间的粘着力时车轮就会发生空转或打滑，使牵引力或制动力急剧减小，还会出现轮轨发热、轮轨擦伤，严重时还会影响机车车辆的安全运行，危害极大。轮轨之间的粘着是一个具有不确定性的复杂时变系统，需要采用特定的控制方法，有效防止牵引空转或制动滑行，并使机车车辆在当前的轨面状态下发挥最大的牵引力或制动力。

在机车车辆防空转滑行及粘着控制领域，常用的方案有差动继电器法、阈值法、多参数联合控制、增稳控制等方法。然而这些控制方法都采用侧重于对空转与滑行的识别和主动减载策略上，难以选取合适的减载时刻、减载率百分比及减载持续时间，均难以最大限度发挥粘着力。

具体地，现有的一种控制方法中，通过实时检测机车轮对加速度，当加速度超过保护阈值时，开始卸载力矩，并在卸载过程中不断搜索加速度峰值，直至检测到加速度的峰值，立即停止卸载力矩，使机车轮对粘着重新恢复。该方法以实时检测的机车运行速度为评判依据，以机车轮对加速度超过某一阈值为切入点，实施机车减载措施，在机车速度传感器正常情况下，该方法很适用，但当其传感器故障时，仅凭借单轴速度、加速度就判断其发生空转滑行，就会造成误判，由于误判减少机车牵引力，影响机车运行。

基于此，现有技术仍然有待改进。

发明内容

为解决上述技术问题，本发明实施例提出一种机车车辆及其加权参数粘着控制方法。以解决现有技术粘着力无法最大限度发挥的技术问题，最大化地利用轮轨粘着力，并且有效防止牵引空转或制动滑行。

本发明实施例所公开的一种加权参数粘着控制方法，其包括：

步骤一 根据机车运行状态，通过查询牵引/电制转矩特性曲线获取牵引电机转矩给定值Tqref，经过传递函数G(S)的一阶低通滤波器生成牵引电机转矩给定控制值Tqrefout，并且，限制Tqmin≤Tqrefout≤Tqmax；

步骤二 根据机车运行状态，通过查询牵引/电制加权粘着控制参考值曲线，获取加权粘着控制参考值TqAdref；

步骤三 基于加权粘着控制反馈因子Adiv1、Adiv2及Adiv3，计算加权粘着控制反馈值TqAdfdb＝Adiv1+Adiv2+Adiv3；

步骤四 将所述加权粘着控制参考值TqAdref与加权粘着控制反馈值TqAdfdb在加权粘着控制PI闭环控制器TqAd中计算，得出加权粘着控制值TqAdout，且限制Tqmin≤TqAdout≤Tqmax；

步骤五 依据所述加权粘着控制值TqAdout和牵引电机转矩给定控制值Tqrefout二者中较小者控制牵引电机转矩Tqout；

步骤六 当TqAdout<Tqrefout时，令Tqrefout＝TqAdout，并记录该时刻T0的前Ts秒内Tqout的最大值Tqoutmax；

从T0时刻开始，按传递函数为m*G(S)一阶低通滤波器，对牵引电机转矩给定值Tqref进行滤波，直到Tqrefout恢复到Tqoutmax的a％时，再按传递函数为G(S)/n一阶低通滤波器，对牵引电机转矩给定值Tqref进行滤波处理，直到Tqrefout等于Tqref；

其中，1＜m＜10；1＜n＜10。

在车辆运行过程中，上述步骤周期性的循环执行。

进一步地，所述机车运行状态包括牵引/电制指令、手柄级位、牵引电机转速中的至少一种。

具体地，

进一步地，所述传递函数G(S)为：

其中，k为增益，ω_c为截止角频率。

进一步地，所述步骤二还包括：

实时监测牵引电机转速传感器状态，

检测并校正当前有效牵引电机转速传感器的转速信号RPM[i](i＝1,2…n)，其中n为机车车辆车轴数，

计算整车牵引电机平均转速值RPMav、最大牵引电机转速值RPMmax、最小牵引电机转速值RPMmin，

计算各转向架的牵引电机平均转速值RPMTav，最大牵引电机转速值RPMTmax、最小牵引电机转速值RPMTmin，

计算车体对地折算转速估算值RPMg。

具体地，包括：

1)检测各个轴牵引电机实际转速RPMR[i](i＝1,2…n),其中n为机车车辆车轴数，计算各个轴牵引电机实际平均速度RPMRav；

2)检测各个轴牵引电机实际电流AMP[i](i＝1,2…n),其中n为机车车辆车轴数，计算各个轴牵引电机平均电流AMPav；

3)当满足

RPMRav>＝RPMRLow，

且RPMR[i]<＝k1*RPMRav，

且AMP[i]<k2*AMPav时,

延时t1秒后判断第i牵引电机转速传感器故障，该故障在控制系统重新上电后自动清除，其中参数RPMRLow一般选取机车车辆速度为3km/h对应的牵引电机转速，k1一般选取0.6～0.9之间的实数，k2一般选取1.1～1.25之间的实数；

4)将无转速传感器故障轴位的转速信号RPMR[i]经一阶低通滤波器滤波后，经过轮径校核，得到有效牵引电机转速传感器的转速信号RPM[i]，将RPM[i]累加除以无速传感器故障轴数得出牵引电机平均转速值RPMav,比较RPM[i]中最大值为RPMmax，最小值为RPMmin；分别比较各转向架RPM[i]中最大值为RPMTmax，最小值为RPMTmin，同时计算各转向架牵引电机平均转速值RPMTav；

5)对RPMav进行定斜率滤波处理得到车体对地折算转速估算值RPMg：在牵引状态下，如果RPMg+ΔN<RPMav则RPMg＝RPMg+ΔN，否则RPMg＝RPMav；在电制状态下，如果RPMg>RPMav+ΔN则RPMg＝RPMg-ΔN，否则RPMg＝RPMav；

其中，ΔN为机车车辆在一个执行周期内以最大加速度2Km/h/s运行时，对应的轮对转速，单位为转每秒。

如果机车车辆装备雷达测速装置，则RPMg的值直接从雷达测速信号获取；如果机车车辆装备非动力轴测速装置，则RPMg的值直接从非动力轴测速信号获取。

进一步地，所述车体对地折算转速估算值RPMg通过雷达测速装置或非动力轴测速装置获取。

进一步地，所述粘着控制反馈因子

Adiv1＝δ1*AdF1，

Adiv2＝δ2*AdF2，

Adiv3＝δ3*AdF3，

其中，δ1、δ2、δ3为大于1且小于5加权系数，

AdF1为瞬态加/减速度反馈参数，AdF2为速度差反馈参数，AdF3为加速度限制反馈参数。

进一步地，步骤四中，所述加权粘着控制PI闭环控制器TqAd采用的传递函数为

TqAdout＝TqAdout+Kp[e(k)-e(k-1)]+Ki*e(k),

其中，Kp、Ki为PI闭环控制器TqAd的比例参数和积分参数，

e(k)为本次TqAdref与TqAdfdb的差值，

e(k-1)为上一次TqAdref与TqAdfdb的差值。

即从T0时刻开始，对牵引电机转矩给定值Tqref进行滤波，直到Tqrefout恢复到Tqoutmax的a％时，所述Kp、Ki分别增大到原Kp、Ki参数的a倍。

进一步地，步骤四中，当所述加权粘着控制反馈值TqAdfdb与加权粘着控制参考值TqAdref之间的差值大于β*TqAdref(0.5<β<1)时,控制系统执行相应轴位撒砂指令。

进一步地，所述AdF1，AdF2和AdF3按车控、架控和轴控的方式进行计算。

具体地，当电机控制采用车控方式时，所述瞬态加(减)速度反馈参数AdF1为本次采样点的整车牵引电机平均转速值RPMav(tk)与上次采样点的整车牵引电机平均转速值RPMav(tk-1)的差值(牵引工况)，或为上次采样点的整车牵引电机平均转速值RPMav(tk-1)与本次采样点的整车牵引电机平均转速值RPMav(tk)的差值(电制工况)，如差值小于零则取零；速度差反馈参数AdF2为整车最大牵引电机转速RPMmax与整车平均牵引电机转速RPMav的差值(牵引工况)或整车平均转速RPMav与整车最小转速RPMmin的差值(电制工况)，如差值小于零则取零；加速度限制反馈参数AdF3为整车平均转速RPMav与车体对地折算转速估算值RPMg的差值(牵引工况)或车体对地折算转速估算值RPMg与整车平均转速RPMav的差值(电制工况)，如差值小于零则取零。

当电机控制采用架控方式时，分别对前、后转向架进行计算。当所计算的转向架为前架时，所述瞬态加(减)速度反馈参数AdF1为本次采样点的前转向架牵引电机平均转速值RPMFav(tk)与上次采样点的前转向架牵引电机平均转速值RPMFav(tk-1)的差值(牵引工况)，或上次采样点的前转向架牵引电机平均转速值RPMFav(tk-1)与前转向架牵引电机平均转速值RPMFav(tk)的差值(电制工况)，如差值小于零则取零；速度差反馈参数AdF2为前转向架最大牵引电机转速RPMFmax与整车牵引电机平均转速RPMav的差值(牵引工况)或整车牵引电机平均转速RPMav与前转向架最小转速RPMFmin的差值(电制工况)，如差值小于零则取零；加速度限制反馈参数AdF3为前转向架平均转速RPMFav与车体对地折算转速估算值RPMg的差值(牵引工况)，或车体对地折算转速估算值RPMg与前转向架平均转速RPMFav的差值(电制工况)，如差值小于零则取零。相同的，当所计算的转向架为后架时，所述瞬态加(减)速度反馈参数AdF1为本次采样点的后转向架牵引电机平均转速值RPMBav(tk)与上次采样点的后转向架牵引电机平均转速值RPMBav(tk-1)的差值，如差值小于零则取零；速度差反馈参数AdF2为后转向架最大牵引电机转速RPMBmax与整车牵引电机平均转速RPMav的差值(牵引工况)或整车牵引电机平均转速RPMav与后转向架最小转速RPMBmin的差值(电制工况)，如差值小于零则取零；加速度限制反馈参数AdF3为后转向架平均转速RPMBav与车体对地折算转速估算值RPMg的差值(牵引工况)，或车体对地折算转速估算值RPMg与后转向架平均转速RPMBav的差值(电制工况)，如差值小于零则取零。

当电机控制采用轴控方式时，分别对各轴进行计算。所述瞬态加(减)速度反馈参数AdF1为本次采样点的当前轴牵引电机转速值RPM[i](tk)与上次采样点当前轴牵引电机转速值RPM[i](tk-1)的差值的绝对值，如差值小于零则取零；速度差反馈参数AdF2为当前轴牵引电机转速RPM[i]与整车牵引电机转速平均值RPMav的差值(牵引工况)或整车牵引电机转速平均值RPMav与当前轴牵引电机转速RPM[i]的差值(电制工况)，如差值小于零则取零；加速度限制反馈参数AdF3为当前轴牵引电机转速RPM[i]与车体对地折算转速估算值RPMg的差值(牵引工况)，或车体对地折算转速估算值RPMg与当前轴牵引电机转速RPM[i]的差值(电制工况)，如差值小于零则取零。

本发明还公开了一种机车或车辆，其采用上述的加权参数粘着控制方法。

采用上述技术方案，本发明至少具有如下有益效果：

本发明通过对机车多运行参数的综合判断，确定最佳减载时刻，减载率百分比及减载持续时间，最大限度发挥机车车辆的粘着力，有效防止牵引空转或制动滑行，使机车车辆在当前的轨面状态下发挥最大的牵引力或制动力。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例的主程序流程图。

图2为本发明实施例的计算转速相关值程序流程图。

图3为本发明实施例的计算加权粘着控制反馈值TqAdfdb子程序流程图。

图4为本发明实施例的架控模式下AdF1、AdF2、AdF3子程序流程图。

图5为本发明实施例的轴控模式下AdF1、AdF2、AdF3子程序流程图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明实施例进一步详细说明。

需要说明的是，本发明实施例中所有使用“第一”和“第二”的表述均是为了区分两个相同名称非相同的实体或者非相同的参量，可见“第一”“第二”仅为了表述的方便，不应理解为对本发明实施例的限定，后续实施例对此不再一一说明。

如图1所示，本发明实施例公开了一种加权参数粘着控制方法的主程序流程图，周期性(一般选取10—20ms)调用此流程图对应软件首先，根据机车车辆运行状态，通过查询牵引/电制转矩特性曲线获取牵引电机转矩给定值Tqref，经过传递函数为G(S)的一阶低通滤波器，生成所述牵引电机转矩给定控制值Tqrefout，且限制Tqmin≤Tqrefout≤Tqmax，作为机车车辆牵引/电制工况电机扭矩给定的一个参考值；其次，根据机车车辆运行状态获取加权粘着控制参考值TqAdref，加权粘着控制反馈值TqAdfdb，将加权粘着控制参考值TqAdref与加权粘着控制反馈值TqAdfdb送入加权粘着控制PI闭环控制器TqAd，得出加权粘着控制值TqAdout,且限制Tqmin≤TqAdout≤Tqmax，作为机车车辆牵引/电制工况电机扭矩给定的另一个参考值；然后，对Tqrefout和TqAdout进行比较，由两者中较小者控制牵引电机转矩Tqout。当TqAdout<Tqrefout时，令Tqrefout＝TqAdout，并记录该时刻T0的前Ts秒内Tqout的最大值Tqoutmax。从T0时刻开始，按传递函数为m*G(S)一阶低通滤波器，对牵引电机转矩给定值Tqref进行滤波，直到Tqrefout恢复到Tqoutmax的a％时，再按传递函数为G(S)/n一阶低通滤波器，对牵引电机转矩给定值Tqref进行滤波处理，直到Tqrefout等于Tqref。

在本发明中，为了实现机车车辆粘着的闭环控制，构造加权粘着控制PI闭环控制器TqAd是其中的关键。因此在软件中设计了加权粘着控制参考值曲线，并根据牵引/电制指令、手柄级位、牵引电机转速等指令和状态参数，实时查询牵引/电制加权粘着控制参考值曲线获取加权粘着控制参考值TqAdref，TqAdref是多参数通过加权相加得到的用作PI闭环控制器TqAd的给定值。在计算加权粘着控制反馈值TqAdfdb的过程中，首先是检测并校正当前有效牵引电机转速RPM[i]，并计算其他相关转速值,程序流程图如图2所示，由检测到的当前有效牵引电机转速RPM[i]，计算出整车平均速度RPMRav。并检测各个轴牵引电机实际电流AMP[i]，计算牵引电机平均电流AMPav。逐个轴判断牵引电机转速传感器的有效性，判断条件如下：当满足RPMRav>＝RPMRLow(RPMRLow一般选取机车车辆速度为3km/h对应的牵引电机转速)且RPMR[i]<＝k1*RPMRav且AMP[i]<k2*AMPav时,延时t1秒后判断第i牵引电机转速传感器故障，其中：k1一般选取0.6～0.9之间的实数，k2一般选取1.1～1.25之间的实数。标记下速度传感器故障轴数，并计算得到整车无速度传感器故障的轴数Sn，前架无速度传感器轴速SFn，后架无速度传感器轴速SBn。将无转速传感器故障轴位的转速信号RPMR[i]经一阶低通滤波器滤波后，经过轮径校核，得到有效牵引电机转速传感器的转速信号RPM[i]。将RPM[i]累加除以无速传感器故障轴数Sn得出牵引电机平均转速值RPMav,比较RPM[i]中最大值为RPMmax，最小值为RPMmin；将前转向架下牵引电机RPM[i]累加除以前转向架无速传感器故障轴数SFn得出前架牵引电机平均转速值RPMFav,比较前转向架下牵引电机RPM[i]中最大值为RPMFmax，最小值为RPMFmin；相同原理得到后架牵引电机平均转速值RPMBav，后架最大牵引电机转速值为RPMBmax，后架最小牵引电机转速值为RPMBmin。

对RPMav进行定斜率滤波处理得到车体对地折算转速估算值RPMg：在牵引状态下，如果RPMg+ΔN<RPMav则RPMg＝RPMg+ΔN，否则RPMg＝RPMav；在电制状态下，如果RPMg>RPMav+ΔN则RPMg＝RPMg-ΔN，否则RPMg＝RPMav。如果机车车辆装备雷达测速装置，则RPMg的值直接从雷达测速信号获取；如果机车车辆装备非动力轴测速装置，则RPMg的值直接从非动力轴测速信号获取。

计算加权粘着控制反馈值TqAdfdb程序流程图如图3所示。该流程需考虑机车车辆的牵引系统控制方式：车控、架控、轴控。若采用车控方式，首先判断机车车辆处于牵引还是制动工况。当机车车辆处于牵引工况时，机车车辆当前牵引电机平均转速值RPMav大于上一周期牵引电机平均转速值RPMavLast时，AdF1＝RPMav-RPMavLast，否则AdF1为零；最大牵引电机转速值RPMmax大于牵引电机平均转速值RPMav时：AdF2＝RPMmax-RPMav，否则AdF2为零；牵引电机平均转速值RPMav大于车体对地折算转速估算值RPMg时：AdF3＝RPMav-RPMg，否则AdF3为零。当机车车辆处于制动工况时，机车车辆当前牵引电机平均转速值RPMav小于上一周期牵引电机平均转速值RPMavLast时，若AdF1＝RPMavLast-RPMav，否则AdF1为零；最小牵引电机转速值RPMmin小于牵引电机平均转速值RPMav时：AdF2＝RPMav-RPMmin，否则AdF2为零；牵引电机平均转速值RPMav小于车体对地折算转速估算值RPMg时：AdF3＝RPMg-RPMav，否则AdF3为零。若采用架控方式，AdF1、AdF2、AdF3子程序流程图如图4所示，其中：jmax为机车车辆转向架数量，需循环判断，当机车车辆处于牵引工况时，判断当前转向架牵引电机平均转速值RPMTav[j](j＝1,…jmax)是否大于该转向架上一周期牵引电机平均转速值RPMTavLast[j](j＝1,…jmax)，若RPMTav[j]>RPMJavLast[j]，则AdF1＝RPMTav[j]-RPMJavLast[j]，否则AdF1为零；判断该转向架牵引电机最大转速值RPMTmax[j](j＝1,…jmax)是否大于整车牵引电机平均转速值RPMav，若RPMTmax[j]>RPMav，则AdF2＝RPMTmax[j]-RPMav，否则AdF2为零；判断该转向架牵引电机平均转速值RPMTav[j](j＝1,…jmax)是否大于车体对地折算转速估算值RPMg，若RPMTav[j]>RPMg，则AdF3＝RPMTav[j]-RPMg，否则AdF3为零。当机车车辆处于制动模式时，判断当前转向架牵引电机平均转速值RPMTav[j](j＝1,…jmax)是否小于该转向架上一周期牵引电机平均转速值RPMTavLast[j](j＝1,…jmax)，若RPMTav[j]<RPMJavLast[j]，则AdF1＝RPMJavLast[j]-RPMTav[j]，否则AdF1为零；判断该转向架牵引电机最小转速值RPMTmin[j](j＝1,…jmax)是否小于整车牵引电机平均转速值RPMav，若RPMTmin[j]<RPMav，则AdF2＝RPMav-RPMTmin[j]，否则AdF2为零；判断该转向架牵引电机平均转速值RPMTav[j](j＝1,…jmax)是否小于车体对地折算转速估算值RPMg，若RPMTav[j]<RPMg，则AdF3＝RPMg-RPMTav[j]，否则AdF3为零。若采用轴控方式，AdF1、AdF2、AdF3子程序流程图如图5所示，其中：imax为机车车辆轴数，需循环判断，当机车车辆处于牵引工况时，判断当前有效牵引电机转速传感器的转速信号RPM[i](i＝1,2…imax)与是否大于其上一周期该轴有效牵引电机转速传感器的转速信号RPMLast[i](i＝1,2…imax)，若RPM[i]>RPMLast[i]，则AdF1＝RPM[i]-RPMLast[i]，否则AdF1为零；判断当前有效牵引电机转速传感器的转速信号RPM[i](i＝1,2…imax)是否大于整车牵引电机平均转速值RPMav，若RPM[i]>RPMav，则AdF2＝RPM[i]-RPMav，否则AdF2为零；判断当前有效牵引电机转速传感器的转速信号RPM[i](i＝1,2…imax)是否大于车体对地折算转速估算值RPMg，若RPM[i]>RPMg，则AdF3＝RPM[i]-RPMg，否则AdF3为零。若为制动模式时，判断当前有效牵引电机转速传感器的转速信号RPM[i](i＝1,2…imax)与是否小于上一周期该轴有效牵引电机转速传感器的转速信号RPMLast[i](i＝1,2…imax)，若RPM[i]<RPMLast[i]，则AdF1＝RPMLast[i]-RPM[i]，否则AdF1为零；判断当前有效牵引电机转速传感器的转速信号RPM[i](i＝1,2…imax)是否小于整车牵引电机平均转速值RPMav，若RPM[i]<RPMav，则AdF2＝RPMav-RPM[i]，否则AdF2为零；判断当前有效牵引电机转速传感器的转速信号RPM[i](i＝1,2…imax)是否小于车体对地折算转速估算值RPMg，若RPM[i]<RPMg，则AdF3＝RPMg-RPM[i]，否则AdF3为零。

由上述计算得到的AdF1、AdF2、AdF3参数值，乘以各自的加权系数δ1、δ2、δ3(大于1且小于5加权系数)，从而得到加权粘着控制反馈值TqAdfdb＝δ1*AdF1+δ2*AdF2+δ3*AdF3。

需要特别指出的是，上述各个实施例中的各个组件或步骤均可以相互交叉、替换、增加、删减，因此，这些合理的排列组合变换形成的组合也应当属于本发明的保护范围，并且不应将本发明的保护范围局限在所述实施例之上。

以上是本发明公开的示例性实施例，上述本发明实施例公开的顺序仅仅为了描述，不代表实施例的优劣。但是应当注意，以上任何实施例的讨论仅为示例性的，并非旨在暗示本发明实施例公开的范围(包括权利要求)被限于这些例子，在不背离权利要求限定的范围的前提下，可以进行多种改变和修改。根据这里描述的公开实施例的方法权利要求的功能、步骤和/或动作不需以任何特定顺序执行。此外，尽管本发明实施例公开的元素可以以个体形式描述或要求，但除非明确限制为单数，也可以理解为多个。

所属领域的普通技术人员应当理解：以上任何实施例的讨论仅为示例性的，并非旨在暗示本发明实施例公开的范围(包括权利要求)被限于这些例子；在本发明实施例的思路下，以上实施例或者不同实施例中的技术特征之间也可以进行组合，并存在如上所述的本发明实施例的不同方面的许多其它变化，为了简明它们没有在细节中提供。因此，凡在本发明实施例的精神和原则之内，所做的任何省略、修改、等同替换、改进等，均应包括在本发明实施例的保护范围之内。

Claims (12)

1.一种加权参数粘着控制方法，其特征在于，包括：

步骤一 根据机车运行状态获取牵引电机转矩给定值Tqref，经过传递函数G(S)生成牵引电机转矩给定控制值Tqrefout；

步骤二 根据机车运行状态获取加权粘着控制参考值TqAdref；

步骤三 基于加权粘着控制反馈因子Adiv1、Adiv2及Adiv3，计算加权粘着控制反馈值TqAdfdb＝Adiv1+Adiv2+Adiv3；

步骤四 将所述加权粘着控制参考值TqAdref与加权粘着控制反馈值TqAdfdb在加权粘着控制PI闭环控制器TqAd中计算，得出加权粘着控制值TqAdout；

步骤五 依据所述加权粘着控制值TqAdout和牵引电机转矩给定控制值Tqrefout二者中较小者控制牵引电机转矩Tqout；

步骤六 当TqAdout<Tqrefout时，令Tqrefout＝TqAdout，并记录该时刻T0的前Ts秒内Tqout的最大值Tqoutmax；

从T0时刻开始，按传递函数为m*G(S)一阶低通滤波器，对牵引电机转矩给定值Tqref进行滤波，直到Tqrefout恢复到Tqoutmax的a％时，再按传递函数为G(S)/n一阶低通滤波器，对牵引电机转矩给定值Tqref进行滤波处理，直到Tqrefout等于Tqref；

其中，1＜m＜10；1＜n＜10；

所述步骤二还包括：

实时监测牵引电机转速传感器状态，

检测并校正当前有效牵引电机转速传感器的转速信号RPM[i](i＝1,2…n)，其中n为机车车辆车轴数，

计算整车牵引电机平均转速值RPMav、最大牵引电机转速值RPMmax、最小牵引电机转速值RPMmin，

计算各转向架的牵引电机平均转速值RPMTav，最大牵引电机转速值RPMTmax、最小牵引电机转速值RPMTmin，

计算车体对地折算转速估算值RPMg；

所述粘着控制反馈因子

Adiv1＝δ1*AdF1，

Adiv2＝δ2*AdF2，

Adiv3＝δ3*AdF3，

其中，δ1、δ2、δ3为大于1且小于5加权系数，

AdF1为瞬态加/减速度反馈参数，AdF2为速度差反馈参数，AdF3为加速度限制反馈参数；

所述AdF1，AdF2和AdF3按车控、架控和轴控的方式进行计算。

2.根据权利要求1所述的方法，其特征在于，所述机车运行状态包括牵引/电制指令、手柄级位、牵引电机转速中的至少一种。

3.根据权利要求1所述的方法，其特征在于，所述传递函数G(S)为：

其中，k为增益，ω_c为截止角频率。

4.根据权利要求1所述的方法，其特征在于，所述车体对地折算转速估算值RPMg通过雷达测速装置或非动力轴测速装置获取。

5.根据权利要求1所述的方法，其特征在于，步骤四中，所述加权粘着控制PI闭环控制器TqAd采用的传递函数为

TqAdout＝TqAdout+Kp[e(k)-e(k-1)]+Ki*e(k),

其中，Kp、Ki为PI闭环控制器TqAd的比例参数和积分参数，

e(k)为本次TqAdref与TqAdfdb的差值，

e(k-1)为上一次TqAdref与TqAdfdb的差值。

6.根据权利要求1所述的方法，其特征在于，步骤四中，当所述加权粘着控制反馈值TqAdfdb与加权粘着控制参考值TqAdref之间的差值大于β*TqAdref时,控制系统执行相应轴位撒砂指令，其中，0.5<β<1。

7.根据权利要求1所述的方法，其特征在于，当电机控制采用车控方式时，

牵引工况下，所述瞬态加/减速度反馈参数为本次采样点的整车牵引电机平均转速值与上次采样点的整车牵引电机平均转速值的差值，如差值小于零则取零；

电制工况下，所述瞬态加/减速度反馈参数为上次采样点的整车牵引电机平均转速值与本次采样点的整车牵引电机平均转速值的差值，如差值小于零则取零。

8.根据权利要求1所述的方法，其特征在于，当电机控制采用车控方式时，

牵引工况下，所述速度差反馈参数为整车最大牵引电机转速与整车平均牵引电机转速的差值，如差值小于零则取零；

电制工况下，所述速度差反馈参数为整车平均转速与整车最小转速的差值，如差值小于零则取零。

9.根据权利要求1所述的方法，其特征在于，当电机控制采用车控方式时，

牵引工况下，加速度限制反馈参数为整车平均转速与车体对地折算转速估算值的差值，如差值小于零则取零；

电制工况下，加速度限制反馈参数为车体对地折算转速估算值与整车平均转速的差值，如差值小于零则取零。

10.根据权利要求1所述的方法，其特征在于，当电机控制采用架控方式时，分别对前、后转向架进行计算。

11.根据权利要求1所述的方法，其特征在于，当电机控制采用轴控方式时，分别对各轴进行计算。

12.一种车辆，其特征在于，采用权利要求1-11任意一项所述的加权参数粘着控制方法。

Images