CN113815685B - 一种列车全程最优目标速度计算方法 - Google Patents
一种列车全程最优目标速度计算方法 Download PDFInfo
- Publication number
- CN113815685B CN113815685B CN202111287094.1A CN202111287094A CN113815685B CN 113815685 B CN113815685 B CN 113815685B CN 202111287094 A CN202111287094 A CN 202111287094A CN 113815685 B CN113815685 B CN 113815685B
- Authority
- CN
- China
- Prior art keywords
- train
- running
- time
- speed
- interval
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004891 communication Methods 0.000 claims abstract description 7
- 238000005265 energy consumption Methods 0.000 claims description 23
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 230000000295 complement effect Effects 0.000 claims description 6
- 238000005457 optimization Methods 0.000 description 7
- 239000010410 layer Substances 0.000 description 6
- 238000005192 partition Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 239000002356 single layer Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L15/00—Indicators provided on the vehicle or train for signalling purposes
- B61L15/0062—On-board target speed calculation or supervision
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L27/00—Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
- B61L27/20—Trackside control of safe travel of vehicle or train, e.g. braking curve calculation
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Train Traffic Observation, Control, And Security (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
本发明公开了一种列车全程最优目标速度计算方法,包括:S1:获取运行基础数据;S2:根据所述列车基础数据,计算全局最优目标速度和每个区间运行时间;S3:判断各区间运行时间是否满足通通时分约束,若是,输出所述全局最优目标速度;否则,进入步骤S4;S4:获取所有所述区间运行时间内与所述通通时分约束偏差最大的区间运行时间;S5:根据所述偏差最大的区间运行时间,调整该区间的局部运行速度;S6:根据所述局部运行速度,得到整个区段内列车的总运行时间;S7:判断总运行时间是否满足总时间约束,若是,返回步骤S3,否则,进入步骤S8;S8:调整所述全局目标速度并返回步骤S4。
Description
技术领域
本发明涉及列车运行控制技术领域,具体涉及一种列车全程最优目标速度计算方法。
背景技术
截至2020年,全国铁路营业里程达到14.63万公里,列车牵引能耗成为铁路运营的主要成本,列车牵引能耗是指由列车牵引过程所引起的能耗,受操纵策略影响,因此需要研究列车节能控制技术。
列车运行过程会受到线路条件约束和运行时刻表约束,其中运行时刻表约束又分为总运行时间约束和通通时分约束。列车在线路上的运行往往是长距离、跨多站的多区间运行方式,在遵循总运行时分的同时,为了减少线路占用,提高线路运输效率,区段总运行时间对于区段内各区间按运行速度不停车通过车站的运行时分做出的约束,称为通通时分约束。现有的列车节能控制研究主要分为两类:列车单区间节能优化和列车多区间节能优化。第一类优化问题仅考虑总运行时间约束进行建模求解,忽略了列车存在跨多站运行的实际情况;而列车多区间节能优化问题主要是对时间变量的建模和求解,这会大大增加问题求解的复杂度,目前已有算法主要包含双层优化和单层优化。双层优化通常采用“时间层-能量层”结合的方式,在计算过程中时间层需反复调用能量层,结构复杂且计算量较大;而单层优化则是将时间约束纳入列车节能优化模型中进行求解,导致模型更加复杂,且这些算法大多针对地铁列车、城际列车等短途运行且站站停列车节能控制,而非跨多站运行模式。
发明内容
本发明的目的在于提供一种列车全程最优目标速度计算方法,以能够提高获得全程最优目标速度的计算速率,并较低列车运行能耗。
本发明解决上述技术问题的技术方案如下:
本发明提供一种列车全程最优目标速度计算方法,包括:
S1:获取运行基础数据;
S2:根据所述列车基础数据,计算全局最优目标速度和每个区间运行时间;
S3:判断各区间运行时间是否满足通通时分约束,若是,输出所述全局最优目标速度;否则,进入步骤S4;
S4:获取所有所述区间运行时间内与所述通通时分约束偏差最大的区间运行时间;
S5:根据所述偏差最大的区间运行时间,调整该区间的局部运行速度;
S6:根据所述局部运行速度,得到整个区段内列车的总运行时间;
S7:判断总运行时间是否满足总时间约束,若是,返回步骤S3,否则,进入步骤S8;
S8:调整所述全局目标速度并返回步骤S4。
可选择地,所述步骤S1中,所述运行基础数据包括列车基本数据、线路基本数据、列车区段运行时间和通通时分约束,所述列车基本数据包括列车整备质量、基本运行阻力特性、列车牵引特性和制动特性;所述线路基本数据包括坡道信息、曲线信息、隧道信息。
可选择地,所述步骤S2包括:
S21:根据所述列车基础数据,建立列车运动学模型和列车运行能耗模型;
S22:根据所述列车运行学模型和所述列车运行能耗模型,构建哈密顿函数;
S23:根据所述哈密顿函数,利用庞特里亚金极大值原理,得到所述全局最优目标速度和每个区间运行时间。
可选择地,所述步骤S21中,所述列车运行学模型为:
其中,ut、ub分别是列车牵引力和制动力的控制系数,满足ut∈[0,1],ub∈[0,1]且ut·ub=0,F(v)为牵引力,B(v)为制动力,w(v)为列车基本运行阻力,g(x)线路附加阻力,v为列车运行速度,x为列车运行位移,t为列车运行时间;
可选择地,所述步骤S21中,所述列车运行能耗模型为:
subject to
v(h1)=0,v(hn+1)=0
其中,其中,为列车在区间(h1,hn+1)内任意位置最晚和最早到达时间,dM/dx和dN/dx为互补松弛因子,满足 F(v)为牵引力,v表示列车运行速度,x表示列车运行位移,t表示列车运行时间,h1表示列车运行区间(h1,hn+1)的起点位置,hn+1表示列车运行区间(h1,hn+1)的终点位置,v(h1)表示列车在h1点的速度,v(hn+1)表示列车在hn+1点的速度,J为能耗最低目标函数。
可选择地,所述步骤S22中,所述哈密顿函数为:
ut、ub分别是列车牵引力和制动力的控制系数,满足ut∈[0,1],ub∈[0,1]且ut·ub=0,F(v)为牵引力,B(v)为制动力,w(v)为列车基本运行阻力,g(x)线路附加阻力,为列车在区间(h1,hn+1)内任意位置最晚和最早到达时间,dM/dx和dN/dx为互补松弛因子,满足 λ1和λ2为伴随变量且v表示列车运行速度,x表示列车运行位移,t表示列车运行时间。
可选择地,所述步骤S23包括:
S231:利用庞特里亚金极大值原理,分别得到列车最优控制工况和恒速速度与伴随变量之间的关系;
S232:计算所述列车最大能力运行所需要的运行时间;
S233:判断所述列车最大能力运行所需要的运行时间是否等于预设运行总时间,若是,输出所述最大能力运行控制工况、所述恒速速度和各区间运行时间,否则,进入步骤S234;
S234:在所述整个区段内建立恒速区段,利用所述恒速区段,得到所述最优控制工况、所述全局最优速度和各区间运行时间。
可选择地,所述步骤S3中,通过如下判断条件判断所述各区间运行时间是否满足通通时分约束:
Ti≤ti,i=1,...,n
其中,Ti为列车在的区间i的运行时间,ti为通通时分约束,n为常数。
本发明具有以下有益效果:
本发明将通通时分约束纳入列车能耗模型中,一方面,能够得到能耗较低、运营成本较低的列车全程最优目标速度;另一方面,能够不扰乱列车运营秩序,确保列车的正常运行。
附图说明
图1为本发明所提供的列车通通时分约束示意图;
图2为本发明所提供的列车全程最优目标速度计算方法流程图;
图3为图2中步骤S2的分步骤流程图;
图4为图3中步骤S23的分步骤流程图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
实施例
首先对本发明中的部分名词作以解释:
全局最优目标速度:整个运行区段内,满足由庞特里亚金极大值原理(Pontryagin’s Maximum Principle,PMP)分析得到的列车能耗最小必要条件的部分牵引工况速度;
局部最优目标速度:某个运行区间内,满足通通时分约束且使得该区间能耗最小的恒速速度;
通通时分:列车从某一站点到另一站点的运行区段中,不停车运行所需的时间。
参考图1所示,Ti为列车在的区间i的运行时间,ti为通通时分约束,(h1,hn+1)代表一个列车运行区段,其中包含n个运行区间,在每个运行区间内需满足通通时分约束(Ti≤ti,i=1,...,n),整个区段内需满足总运行时间约束
本发明提供一种列车全程最优目标速度计算方法,参考图2所示,包括:
S1:获取运行基础数据;
这里,所述运行基础数据包括列车基本数据、线路基本数据、列车区段运行时间Td1和通通时分约束ti,所述列车基本数据包括列车整备质量M、基本运行阻力特性w(v)、列车牵引特性F(v)和制动特性B(v);所述线路基本数据包括坡道信息、曲线信息、隧道信息g(x)。
S2:根据所述列车基础数据,计算全局最优目标速度和每个区间运行时间;
可选择地,参考图3所示,所述步骤S2包括:
S21:根据所述列车基础数据,建立列车运动学模型和列车运行能耗模型;
这里,所述列车运行学模型为:
其中,ut、ub分别是列车牵引力和制动力的控制系数,满足ut∈[0,1],ub∈[0,1]且ut·ub=0,F(v)为牵引力,B(v)为制动力,w(v)为列车基本运行阻力,g(x)线路附加阻力,v为列车运行速度,x为列车运行位移,t为列车运行时间;
所述列车运行能耗模型为:
subject to
v(h1)=0,v(hn+1)=0
其中,其中,为列车在区间(h1,hn+1)内任意位置最晚和最早到达时间,dM/dx和dN/dx为互补松弛因子,满足 F(v)为牵引力,v表示列车运行速度,x表示列车运行位移,t表示列车运行时间,h1表示列车运行区间(h1,hn+1)的起点位置,hn+1表示列车运行区间(h1,hn+1)的终点位置,v(h1)表示列车在h1点的速度,v(hn+1)表示列车在hn+1点的速度,J为能耗最低目标函数。
S22:根据所述列车运行学模型和所述列车运行能耗模型,构建哈密顿函数;
所述哈密顿函数为:
ut、ub分别是列车牵引力和制动力的控制系数,满足ut∈[0,1],ub∈[0,1]且ut·ub=0,F(v)为牵引力,B(v)为制动力,w(v)为列车基本运行阻力,g(x)线路附加阻力,为列车在区间(h1,hn+1)内任意位置最晚和最早到达时间,dM/dx和dN/dx为互补松弛因子,满足 λ1和λ2为伴随变量且v表示列车运行速度,x表示列车运行位移,t表示列车运行时间。
S23:根据所述哈密顿函数,利用庞特里亚金极大值原理,得到所述全局最优目标速度和每个区间运行时间。
可选择地,参考图4所示,所述步骤S23包括:
S231:利用庞特里亚金极大值原理,分别得到列车最优控制工况和恒速速度与伴随变量之间的关系;
S232:计算所述列车最大能力运行所需要的运行时间;
S233:判断所述列车最大能力运行所需要的运行时间是否等于预设运行总时间,若是,输出所述最大能力运行控制工况、所述恒速速度和各区间运行时间,否则,进入步骤S234;
S234:在所述整个区段内建立恒速区段,利用所述恒速区段,得到所述最优控制工况、所述全局最优速度和各区间运行时间。
这里,预设运行总时间可以是任意一个给定的预设运行总时间,本领域技术人员可根据实际情况进行设定,本发明不做具体限制。
具体地,定义新的伴随变量θ(θ=λ2/v),因而哈密顿函数可以进一步转换为:
根据PMP,可得到5种列车最优控制工况,如表1所示:
表1列车运行模式和最优控制工况下的控制系数
在不考虑通通时分约束的情况下,即同时令dM/dx=0和dN/dx=0,得到伴随变量微分方程,及伴随变量与全局最优速度之间的关系,其中Vg为全局最优速度:
-Vg 2w'(Vg)=λ1
之后,求解列车全局最优速度和各区间运行时间
(1)计算列车最大能力运行所需要的总运行时间,与给定总运行时间进行比较;
(2)若两者相同则返回最大能力运行控制工况序列、恒速速度及各区间运行时间;
(3)若前者小于后者,初始化λ1与Vg,划分并建立坡道分区列表,在分区列表首尾分别插入区间起点恒速区和区间终点恒速区,遍历分区列表,实现满足总时间约束和限速约束的部分牵引工况恒速区之间连接,返回最优控制工况序列、Vg和对应的各区间运行时间Ti。
S3:判断各区间运行时间是否满足通通时分约束,若是,输出所述全局最优目标速度;否则,进入步骤S4;
通过如下判断条件判断所述各区间运行时间是否满足通通时分约束:
Ti≤ti,i=1,...,n
其中,Ti为列车在的区间i的运行时间,ti为通通时分约束,n为常数。
S4:获取所有所述区间运行时间内与所述通通时分约束偏差最大的区间运行时间;
主要通过搜索得到,即搜索时间偏差最大的区间i(maxΔT=|Ti-ti|),在不考虑总运行时间约束的情况下,计算满足该区间通通时分约束的局部最优速度Vi,其具体步骤如下:
建立第i个区间列车最小牵引能耗函数:
subject to
v(a)=v(d)=Vg,v(b)=v(c)=Vi
其中,
ψ(Vg)=Vg 2w'(Vg)
之后,分析全局、局部最优恒速速度与伴随变量之间的关系:
调整目标速度,计算整个区段的列车总运行时间。
(1)令Vs=V(1)=Vg,Vx=0,重置全局目标速度V(2)(满足Vx<V(2)<Vs,V(2)=1/2(Vs-Vx));
(2)局部最优速度调整方法为,若Ti>ti,增大Vi,否则减小,使得当列车速度为Vi时,满足第i个区间的通通时分约束;
S5:根据所述偏差最大的区间运行时间,调整该区间的局部运行速度;
S6:根据所述局部运行速度,得到整个区段内列车的总运行时间;
S7:判断总运行时间是否满足总时间约束,若是,返回步骤S3,否则,进入步骤S8;
S8:调整所述全局目标速度并返回步骤S4。
本发明具有以下有益效果:
本发明将通通时分约束纳入列车能耗模型中,一方面,能够得到能耗较低、运营成本较低的列车全程最优目标速度;另一方面,能够不扰乱列车运营秩序,确保列车的正常运行。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种列车全程最优目标速度计算方法,其特征在于,包括:
S1:获取运行基础数据;
S2:根据所述运行基础数据,计算全局最优目标速度和每个区间运行时间;
S3:判断各区间运行时间是否满足通通时分约束,若是,输出所述全局最优目标速度;否则,进入步骤S4;
S4:获取所有所述区间运行时间内与所述通通时分约束偏差最大的区间运行时间;
S5:根据所述偏差最大的区间运行时间,调整该区间的局部运行速度;
S6:根据所述局部运行速度,得到整个区段内列车的总运行时间;
S7:判断总运行时间是否满足总时间约束,若是,返回步骤S3,否则,进入步骤S8;
S8:调整所述全局最优目标速度并返回步骤S4;
所述步骤S2包括:
S21:根据列车基础数据,建立列车运动学模型和列车运行能耗模型;
S22:根据所述列车运行学模型和所述列车运行能耗模型,构建哈密顿函数;
S23:根据所述哈密顿函数,利用庞特里亚金极大值原理,得到所述全局最优目标速度和每个区间运行时间;
所述步骤S21中,所述列车运行能耗模型为:
subject to
v(h1)=0,v(hn+1)=0
2.根据权利要求1所述的列车全程最优目标速度计算方法,其特征在于,所述步骤S1中,所述运行基础数据包括列车基本数据、线路基本数据、列车区段运行时间和通通时分约束,所述列车基本数据包括列车整备质量、基本运行阻力特性、列车牵引特性和制动特性;所述线路基本数据包括坡道信息、曲线信息和隧道信息。
5.根据权利要求1所述的列车全程最优目标速度计算方法,其特征在于,所述步骤S23包括:
S231:利用庞特里亚金极大值原理,分别得到列车最优控制工况和恒速速度与伴随变量之间的关系;
S232:计算所述列车最大能力运行所需要的运行时间;
S233:判断所述列车最大能力运行所需要的运行时间是否等于预设运行总时间,若是,输出最大能力运行控制工况、所述恒速速度和各区间运行时间,否则,进入步骤S234;
S234:在所述整个区段内建立恒速区段,利用所述恒速区段,得到所述最优控制工况、所述全局最优目标速度和各区间运行时间。
6.根据权利要求1所述的列车全程最优目标速度计算方法,其特征在于,所述步骤S3中,通过如下判断条件判断所述各区间运行时间是否满足通通时分约束:
Ti≤ti,i=1,...,n
其中,Ti为列车在的区间i的运行时间,ti为通通时分约束,n为常数。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111287094.1A CN113815685B (zh) | 2021-11-02 | 2021-11-02 | 一种列车全程最优目标速度计算方法 |
LU502941A LU502941B1 (de) | 2021-11-02 | 2022-10-21 | Ein Verfahren zur Berechnung der optimalen Zielgeschwindigkeit eines Zuges während seiner gesamten Fahrt |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111287094.1A CN113815685B (zh) | 2021-11-02 | 2021-11-02 | 一种列车全程最优目标速度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113815685A CN113815685A (zh) | 2021-12-21 |
CN113815685B true CN113815685B (zh) | 2022-06-17 |
Family
ID=78919529
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111287094.1A Active CN113815685B (zh) | 2021-11-02 | 2021-11-02 | 一种列车全程最优目标速度计算方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113815685B (zh) |
LU (1) | LU502941B1 (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101875361A (zh) * | 2009-04-29 | 2010-11-03 | 华为技术有限公司 | 列车速度调整方法和装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06191408A (ja) * | 1992-12-25 | 1994-07-12 | Daido Signal Co Ltd | 列車運転時隔線図の作成方法及び該線図から運転時隔を算定する方法 |
CN108764571B (zh) * | 2018-05-28 | 2020-04-14 | 西南交通大学 | 一种重载列车运行的多目标优化方法 |
CN110738369A (zh) * | 2019-10-15 | 2020-01-31 | 西南交通大学 | 一种城市轨道交通列车的运行速度优化方法 |
CN111688760B (zh) * | 2020-06-23 | 2022-04-08 | 北京全路通信信号研究设计院集团有限公司 | 一种通过陡坡路段的列车快速节能优化方法及装置 |
CN112633598A (zh) * | 2020-12-30 | 2021-04-09 | 南京理工大学 | 城市轨道交通列车速度曲线与时刻表综合节能优化方法 |
-
2021
- 2021-11-02 CN CN202111287094.1A patent/CN113815685B/zh active Active
-
2022
- 2022-10-21 LU LU502941A patent/LU502941B1/de active IP Right Grant
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101875361A (zh) * | 2009-04-29 | 2010-11-03 | 华为技术有限公司 | 列车速度调整方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
LU502941B1 (de) | 2023-04-21 |
CN113815685A (zh) | 2021-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110782079B (zh) | 基于停站方案节能的高速铁路列车运行图的调整方法 | |
Li et al. | Energy minimization in dynamic train scheduling and control for metro rail operations | |
CN111311913B (zh) | 一种提高公路变窄路段通行效率的控制方法和系统 | |
CN108973998B (zh) | 一种基于mpc的异质车辆队列分布式节能控制方法 | |
Zhou et al. | Integrated optimization on train control and timetable to minimize net energy consumption of metro lines | |
CN111968377B (zh) | 面向节油和驾驶舒适的基于车辆网的车辆轨迹优化方法 | |
CN110182215A (zh) | 一种汽车经济性巡航控制方法及装置 | |
CN109703606B (zh) | 基于历史运行数据的高速列车智能驾驶控制方法 | |
CN113593275B (zh) | 一种基于公交信号优先的交叉口网联自动驾驶方法 | |
CN108510127B (zh) | 基于再生能利用的城轨列车运行过程优化方法 | |
CN106056238B (zh) | 列车区间运行轨迹的规划方法 | |
CN113299107A (zh) | 一种多目标融合的交叉口动态车辆网联速度引导控制方法 | |
CN112255918B (zh) | 汽车队列优化控制的方法及系统 | |
Xu et al. | An integrated micro-macro approach for high-speed railway energy-efficient timetabling problem | |
CN109398426B (zh) | 一种定时条件下基于离散蚁群算法的节能驾驶策略寻优方法 | |
CN104656452A (zh) | 一种基于矩阵离散算法的地铁列车优化控制方法及装置 | |
CN115743117A (zh) | 一种基于扰动观测的智能网联电动车队协同生态驾驶方法 | |
CN113815685B (zh) | 一种列车全程最优目标速度计算方法 | |
CN110364003A (zh) | 基于车联网的交叉路口双线路有轨电车信号优先控制方法 | |
CN108592930B (zh) | 基于大系统优化和车联网的新能源汽车智能辅助驾驶方法 | |
US20220185345A1 (en) | Control system for operating long vehicles | |
CN109101672B (zh) | 轨道交通系统列车运行能效研究方法 | |
CN113591301B (zh) | 一种城市轨道交通列车运行参数优化算法 | |
CN114187771A (zh) | 基于协同式自适应巡航控制的公交车驾驶控制方法及系统 | |
Liu et al. | Real-time Adjustment and optimization for Platoon Operation of High-speed Trains |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |