高速列车动力学性能参数优化迭代设计

姜杰; 丁国富; 邹益胜; 张海柱; 黄海于; 黎荣; 张剑

doi:10.3969/j.issn.0258-2724.20220071

高速列车动力学性能参数优化迭代设计

doi: 10.3969/j.issn.0258-2724.20220071

西南交通大学机械工程学院，四川成都 610031

基金项目: 国家重点研发计划（2017YFB1201201）

详细信息

作者简介:
姜杰（1988—），男，博士研究生，研究方向为复杂产品概念设计与多学科优化，E-mail：jj13688158635@126.com

通讯作者:
丁国富（1972—），男，教授，博士，研究方向为数字化设计与制造，E-mail：dingguofu@swjtu.edu.cn

中图分类号: U238
计量
- 文章访问数: 39
- HTML全文浏览量: 22
- PDF下载量: 17
- 被引次数: 0
出版历程
- 收稿日期: 2022-01-22
- 修回日期: 2022-05-17
- 网络出版日期: 2025-04-11

Iterative Optimization Design for Dynamic Performance Parameters of High-Speed Trains

School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China

摘要

摘要:
针对高速列车动力学性能参数优化求解问题，首先，搭建高速列车动力学性能优化迭代设计框架，提取动力学属性设计参数并建立动力学性能分析模型；其次，基于设计参数重要度分析和自组织映射缩减设计空间维度，通过多学科领域耦合仿真计算生成性能参数试验样本集；最后，构建多工况下的多目标优化模型，并在此基础上建立基于贝叶斯优化随机森林的高速列车多工况代理模型，通过改进NSGA-Ⅱ算法找出满意的设计参数集. 以某工况为例，实验结果表明：优化后的横向平稳性、垂向平稳性、轮轨垂向力、轮轴横向力、脱轨系数、轮重减载率和倾覆系数性能分别提升1.14%、3.19%、2.86%、2.30%、8.33%、2.77%、8.11%，验证了所提迭代设计方法有效可行，对复杂装备正向创新设计具有一定参考价值.
- 高速列车 /
- 动力学性能参数 /
- 多领域耦合仿真 /
- 优化迭代
Abstract:
To optimize the dynamic performance parameters of high-speed trains, an iterative design framework for dynamic performance optimization of high-speed trains was established. First, dynamic design parameters were extracted, and a dynamic performance analysis model was constructed. Next, the design space was reduced through importance analysis of the design parameters combined with self-organizing mapping, and an experimental dataset of performance parameters was generated through multidisciplinary coupled simulation. Finally, a multi-objective optimization model was established under multiple working conditions. Based on Bayesian optimization and random forest, a surrogate model for multiple working conditions was developed. An improved non-dominated sorting genetic algorithm Ⅱ (NSGA-Ⅱ) was employed to identify a set of optimal design parameters. Following optimization, the experiment conducted under a representative working condition demonstrates performance improvements of 1.14%, 3.19%, 2.86%, 2.30%, 8.33%, 2.77%, and 8.11% in lateral ride comfort, vertical ride comfort, vertical wheel-rail force, lateral wheelset force, derailment coefficient, wheel load reduction ratio, and overturning coefficient, respectively. These findings validate the feasibility and effectiveness of the proposed iterative design method, providing references for the forward innovative design of complex equipment.
- high-speed trains /
- dynamic performance parameters /
- multi-domain coupled simulation /
- iterative optimization

HTML全文

现行电气化铁路采用换相（异相）单边供电^[1-2]，以减少负序对三相电网的影响^[3]，进而降低电能质量问题中的三相电压不平衡^[4]. 为此，需要在接触网的各供电分区之间增设电分相以进行隔离. 然而，电分相^[5]的设置导致了无电区，给铁路运行带来诸多问题^[6-10]. 同相供电可以同时解决牵引变电所的负序问题及其出口处的电分相问题. 在文献[11]提出同相供电方案之后，研究者们基于无功型对称补偿技术^[12-17]进一步开发了牵引变电所有功型同相供电方法^[18-19]. 该方法能以最小补偿容量使三相电压不平衡度达标^[20]，同时消除变电所出口处的电分相. 为进一步消除分区所处和全线的电分相和无电区，不仅需要实现两个牵引变电所间的双边供电和更大范围的贯通供电^[18,21-22]，还需应对铁路贯通供电可能对电网带来的新问题.

1. 贯通供电

贯通供电是两个及以上牵引变电所同相供电和分区所双边供电的结合，能够确保现代交直交列车的牵引和（电）制动都需要的牵引网不间断供电. 最简单的贯通供电是采用单相牵引变压器的牵引变电所同相供电和分区所双边供电的结合，如图1所示. 图中，分区所纵向断路器K需要合闸才能实现贯通供电，否则，仅为牵引变电所的同相供电，属于单边供电. 电网中，将分区所断路器K合闸的操作称为合环，合环需满足电网规定的条件，否则无法进行贯通供电.

图 1 同相单边供电和贯通供电空载示意

Figure 1. Co-phase single-end power supply and interconnected power supply under no-load condition

下载: 全尺寸图片幻灯片

虽然俄罗斯、乌克兰等前苏联国家^[2]及韩国的电气化铁路一直采用双边供电，但其牵引变电所仍然采用异相供电，即换相接入电网，并不属于本文定义的贯通供电.

图1（a）所示的同相单边供电牵引网通过单个牵引变电所与电网的单个公共连接点连接后向列车供电. 而图1（b）所示的贯通供电牵引网则通过多个牵引变电所与电网的多个公共连接点连接后向列车供电，此时，牵引网成为电网中的连接支路，牵引供电系统和电网的拓扑结构发生变化. 以牵引网空载情况为例，双边供电的牵引变电所和牵引网串联后在电网不同公共连接点之间又增生了一条新的支路，形成新的电流分量，称为均衡电流或环路电流，或者说增生了一个新的功率分量，称为穿越功率. 与贯通供电相比，由于同相单边供电分区所（SP）的纵向断路器K分闸，牵引网不贯通，空载工况下不会产生穿越功率.

图2为同相单边和贯通供电的列车牵引用电和再生发电工况示意. 从图中可见，同相单边供电时，不同牵引变电所的列车牵引用电功率和列车再生发电功率在贯通供电时得以叠加，虽然列车的用电和发电总量不变，但贯通供电后，从电网用电或向电网发电的最大功率会减小，对应的负序最大值也会减小. 从这个意义上说，贯通供电优于单边供电.

图 2 同相单边和贯通供电运行工况示意

Figure 2. Operation under co-phase single-end power supply and interconnected power supply

下载: 全尺寸图片幻灯片

但与同相单边供电相比，贯通供电给电网带来新的问题，即能否合环以及合环后的穿越功率问题. 本文讨论电网的正常状态，关于故障状态和对应的继电保护技术，可参见文献[18]、[23-24].

2. 贯通供电合环

以图1（b）中选择两个牵引变电所构成的双边供电合环情形加以讨论，如图3所示. 假设实施双边供电的两个变电所SS₁、SS₂的电源分别来自有（等效）电气连接的电网公共连接点PCC₁、PCC₂，当牵引变电所SS₁、SS₂之间分区所的纵向断路器K合闸，就构成双边供电.

图 3 双边供电等效结构图

Figure 3. Equivalent structure of bilateral power supply

下载: 全尺寸图片幻灯片

实现双边供电必须在同一电网内进行. 电气化铁路的双边供电与电网的电磁合环有所不同：电网电磁合环为三相合环，而双边供电是单相合环；电网电磁合环的开关是常开的，而双边供电的纵向开关是常闭的. 尽管如此，由于双边供电和电网的电磁合环都将改变电网的拓扑结构和运行方式，因此，必须满足电网规定的合环条件.

《青海省电力系统调度规程》^[25]规定了电网的解合环操作，其中合环操作：“必须相位相同，电压差一般允许在20%之内，相角差一般不超过20°，以确保合环时不因环路电流过大引起潮流的变化而超过继电保护、系统稳定和设备容量等方面的限额. 合环后应当及时通知有关调度及厂站. 对于较复杂环网的操作，应当先进行模拟计算. ”

《中国南方电网电力调度管理规程》^[26]解合环操作规定：1）必须确保相序相位正确才能合环；2）解、合环操作必须确保解、合环后潮流不超过稳定极限、设备不过负荷、电压在正常范围内，不引起继电保护和安全自动装置误动；3）具备条件时，合环操作应使用同期装置；4） 500 kV合环时电压差一般不应超过额定电压10%，220 kV不应超过额定电压20%；500 kV系统合环一般应检同期合环，有困难时应启用合环开关的同期装置检查相角差；500 kV合环时相角差一般不应超过20°，220 kV一般不应超过25°.

图3中PCC₁和PCC₂之间的电压差和相角差可以于K合闸之前在分区所测得，也可以在已知PCC₁和PCC₂之间的电气连接关系时通过潮流计算得到^[3].

实际电网中，图3所示PCC₁和PCC₂之间的电气连接方式多种多样，作为举例，这里考虑220 kV输电线路对铁路两个牵引变电所的专用线供电方式，则三相等值电路可由图4描述. 图中，Z_A、Z_B和Z_C为输电线的三相阻抗，L为牵引变电所之间的距离.

图 4 三相等值电路示意

Figure 4. Three-phase equivalent circuit

下载: 全尺寸图片幻灯片

首先，分析分区所纵向断路器K合闸之前PCC₁和PCC₂之间的电压差和相角差是否满足合环规程的规定.

当末端负荷电流 $\dot I$ 通过电阻为R、电抗为X，即阻抗为 $Z = R + {\text{j}}X$ 的线路时，首端电压 ${\dot U_{\text{f}}}$ 与末端电压 ${\dot U_{\text{e}}}$ 的电压降 $\Delta \dot U$ 、电压差 $\Delta U$ 和相位差θ的相量图如所示^[1,15]. $\Delta U$ 和θ分别为^[1,27]

图 5 电压差与相角差示意

Figure 5. Voltage difference and phase angle difference

下载: 全尺寸图片幻灯片

$\begin{split} &\Delta U = \left| {\dot U}_{\mathrm{f}} \right| - \left| {\dot U}_{\mathrm{e}} \right|= \\ & \quad \sqrt {{{\left( {U_{\mathrm{e}} + XI\sin \;\varphi + RI\cos\; \varphi } \right)}^2} + {{\left( {XI\cos \;\varphi - RI\sin \;\varphi } \right)}^2}} - \\& \quad U_{\mathrm{e}} \approx XI\sin \;\varphi + RI\cos\; \varphi, \end{split}$

(1)

$\theta = \arctan \frac{{XI\cos \varphi - RI\sin \varphi }}{{XI\sin \varphi + RI\cos \varphi }},$

(2)

式中： $\varphi$ 为电流 $\dot I$ 的功率因数角，以滞后为正.

220 kV输电线路最大输送功率为500 MV•A^[3]. 计算中，输电线采用2 × LGJQ-500二分裂导线，取单位长阻抗z₀=0.0312 + j0.3052 Ω/km，PCC₂处的输送功率分别取200、300、500 MV•A. 考虑相邻牵引变电所的间距L与输电线匹配，L分别取50、60、70、80 km，输送功率的功率因数分别取0.95和0.90，ΔU和θ的计算结果分别列于表1和表2. 其中，ΔU的百分值以额定电压220 kV为基准.

表 1 功率因数为0.95（滞后）时的电压差ΔU和相角差θ

Table 1. Voltage difference (ΔU) and phase angle difference (θ) at a power factor of 0.95 (lagging)

L/km	传输功率 200 MV•A			传输功率 300 MV•A			传输功率 500 MV•A
L/km	ΔU/kV	ΔU/%	θ/（^o）	ΔU/kV	ΔU/%	θ/（^o）	ΔU/kV	ΔU/%	θ/（^o）
50	3.49	2.74	3.23	5.38	4.23	4.78	9.44	7.43	7.74
60	4.23	3.33	3.85	6.56	5.16	5.69	11.60	9.13	9.16
70	4.99	3.93	4.47	7.77	6.12	6.58	13.85	10.90	10.53
80	5.77	4.54	5.08	9.02	7.10	7.45	16.17	12.73	11.85

下载: 导出CSV

| 显示表格

表 2 功率因数为0.90（滞后）时的电压差ΔU和相角差θ

Table 2. Voltage difference (ΔU) and phase angle difference (θ) at a power factor of 0.90 (lagging)

L/km	传输功率 200 MV•A			传输功率 300 MV•A			传输功率 500 MV•A
L/km	ΔU/kV	ΔU/%	θ/（^o）	ΔU/kV	ΔU/%	θ/（^o）	ΔU/kV	ΔU/%	θ/（^o）
50	4.41	3.47	2.99	6.74	5.30	4.41	11.63	9.16	7.10
60	5.33	4.20	3.56	8.17	6.44	5.23	14.19	11.17	8.37
70	6.26	4.93	4.13	9.64	7.59	6.04	16.81	13.24	9.60
80	7.21	5.68	4.68	11.13	8.76	6.84	19.50	15.35	10.78

下载: 导出CSV

| 显示表格

从图5和表1、2可见，输送（视在）功率的功率因数越高，电压差越小，但相角差越大. 总的看：在牵引变电所间距不超过80 km下（通常如此），双边供电时分区所的合环电压差不大于16.00%，相角差不大于12.00°，低于合环规程的规定值，符合要求.

实践中，双边供电除了仿真计算之外，在合环前应选择牵引网处于空载状态，并检测分区所两侧接触网的电压差和相角差，在允许的电压差和相角差范围内，才能令纵向断路器K合闸.

双边供电一旦合环就可长期运行.

3. 穿越功率分析与计算

牵引供电系统的正常状态可分为牵引、再生和空载3种工况. 其中，空载工况下，双边供电与单边供电才有本质区别，双边供电的牵引供电系统与电网（等效）输电线并联，造成均衡电流或者穿越功率. 从图6可见，空载工况下，双边供电的穿越功率在一个变电所（如SS₁）表现为用电状态，在另一个变电所（如SS₂）则表现为发电状态，更多牵引变电所贯通供电时，穿越功率可能穿越得更远.

图 6 空载工况取流示意

Figure 6. Current distribution under no-load condition

下载: 全尺寸图片幻灯片

3.1 等效电路与均衡电流计算

对应图1和图4，将归算到电网侧的三相等效电路示于图7. 图中，Z_d为三相输电线的相阻抗，即Z_A=Z_B=Z_C=Z_d； ${Z_{{\text{J1}}}}$ 和 ${Z_{{\text{J}}2}}$ 为牵引变电所SS₁和SS₂的进线相阻抗； $Z_{{\text{T1,e}}}$ 和 $Z_{{\text{T}}2,{\mathrm{e}}}$ 为归算到电网侧的单相牵引变压器阻抗； $Z_{\text{q,e}}$ 为归算到电网侧的牵引网阻抗；LC为电力机车.

图 7 双边供电的三相等效电路

Figure 7. Three-phase equivalent circuit for bilateral power supply

下载: 全尺寸图片幻灯片

考虑牵引网空载工况，选择A相或B相可得计算均衡电流的单相等效电路，如所示. 图中，用 ${Z_{\text{d}}} 、{Z_{\text{q}}}$ 分别表示电网输电线与牵引供电系统（包括进线相阻抗）的归算阻抗. 设总电流为 $\dot I$ ，通过电网输电线的电流为 ${\dot I_{\text{d}}}$ ，通过牵引供电系统的电流为 ${\dot I_{\text{q}}}$ ，则有支路电压方程如式（3）所示.

图 8 均衡电流示意

Figure 8. Equalizing current

下载: 全尺寸图片幻灯片

${\dot I_{\text{q}}}{{\textit{Z}}_{\text{q}}} = {{\textit{Z}}_{\text{d}}}{\dot I_{\text{d}}}.$

(3)

设Z_J1=Z_J2=Z_J，则求得均衡电流I_q为

${\dot I_{\text{q}}} = \frac{{{{\textit{Z}}_{\text{d}}}}}{{{{\textit{Z}}_{\text{q}}}}}{\dot I_{\text{d}}} = \frac{{{{\textit{Z}}_{\text{d}}}}}{{2{{\textit{Z}}_{\text{J}}} + ({{\textit{Z}}_{\text{T}}} + \dfrac{1}{2}{{Z}_{\text{q,e}}})k_{\text{T}}^2}}{\dot I_{\text{d}}},$

(4)

式中：k_T为牵引变压器变比，即电网进线的线电压与牵引母线额定电压的比值.

定义分流比为η，即均衡电流与电网输电线中的电流之比，则由式（4）得

$\eta = \left| {\frac{{{{\dot I}_{\text{q}}}}}{{{{\dot I}_{\text{d}}}}}} \right| = \left| {\frac{{{{\textit{Z}}_{\text{d}}}}}{{{{\textit{Z}}_{\text{q}}}}}} \right| = \left| {\frac{{{{\textit{Z}}_{\text{d}}}}}{{2{{\textit{Z}}_{\text{J}}} + ({{\textit{Z}}_{\text{T}}} + \dfrac{1}{2}{{Z}_{\text{q,e}}})k_{\text{T}}^2}}} \right|.$

(5)

精确计算需要掌握具体数据. 假设电网PCC₁和PCC₂之间220 kV输电线的长度为50 km，PCC₁（PCC₂）到牵引变电所SS₁（SS₂）的进线长度为10 km，输电线单位长阻抗同前，牵引变压器均为单相变压器，额定容量均为31.5 MV•A，短路电压百分比为10.5%，则求得归算到牵引侧的漏抗Z_T= 0.2134 + j2.52 Ω，相邻牵引变电所间距50 km，直供方式牵引网空载电压为27.5 kV，变比k_T=8，单线铁路、单链型悬挂牵引网单位长阻抗z=0.232 + j0.515 Ω/km，代入式（5）可得η=1.44%. 工程上，可以用

$\eta \approx \frac{1}{{k_{\text{T}}^2}},$

(6)

进行估算，计算得 $\eta$ =1/64≈1.56%，接近且略大于精确计算值.

3.2 穿越功率计算与测量

设电网输电线的额定电压为U_N，由式（3）得到穿越功率S_q为

${S _{\text{q}}} = {U_{\text{N}}}{I_{\text{q}}} = \frac{{{{\textit{Z}}_{\text{d}}}}}{{{{\textit{Z}}_{\text{q}}}}}{U_{\text{N}}}{I_{\text{d}}},$

(7)

或者已知电网输电线三相功率S时，由式（8）估算得牵引网（单相）穿越功率.

${S _{\text{q}}} \approx \frac{S}{{\sqrt 3 k_{\text{T}}^{\text{2}}}}.$

(8)

一般110 kV输电线最大输送功率为50 MV•A，则直供方式时，变比k_T=4，由式（8）估算得双边供电最大穿越功率为1.8 MV•A，占比3.6%；220 kV输电线最大输送功率为500 MV•A，则直供方式时，k_T=8，由式（8）估算得双边供电穿最大越功率为4.6 MV•A，占比1.8%.

上述定义和讨论的均衡电流和穿越功率均沿线路传输，故称为纵向分量. 穿越功率可进一步分为有功分量和无功分量（一般为感性），都属于纵向分量. 因此，如果穿越功率存在，则在牵引变电所进线、牵引馈线以及牵引网的任何方便的部位都能测量得到. 还有一个像负荷一样的横向分量，这就是线路的容性无功功率，也叫充电功率. 穿越功率取决于电网对牵引供电系统的供电方式和电网潮流，充电功率则沿输电线和牵引网分布，取决于其长度和电压等级. 线路充电功率在分区所纵向断路器合闸前后的总量不变，但合闸后将在系统中重新分配，此时也存在一个分流点. 当分流点与分区所不重合时，分区所会测到这个分量不是0，而当分流点与分区所重合时，分区所的这个分量就是0.

要判断分区所合环后是否存在穿越功率，则应在牵引网空载时测量牵引网的有功分量，才能准确反映穿越功率情况：若检测到的有功分量为0，则穿越功率中的（纵向）无功分量亦为0，穿越功率即为0；而此时若检测到的无功分量不为0，只能说明存在充电功率（横向分量），不是穿越功率. 下面会看到，当电网对铁路牵引变电所进行树形供电时，理论上穿越功率为0.

4. 穿越功率解决方案

正常运行状态下，双边供电与单边供电在潮流分布上的唯一区别就是穿越功率问题. 由上分析可知，穿越功率占比很小，不会对电网潮流造成大的影响. 但是，为不产生穿越功率或不向电网发电（分量），仅达到与单边供电相同或更好的效果，必须关注并解决穿越功率的问题.

1）树形供电

电网的变电站用同一（分段）母线公共连接点（PCC）为多个铁路变电所供电，这种供电方式称为树形供电^[19,22,28]，如图9所示. 此时，对应图3，可认为PCC₁与PCC₂合并，电网输电线长度为0，即令式（6）中Z_d=0，得η=0. 因此，铁路实施树形双边供电不会产生均衡电流和穿越功率，而负序等电能质量影响与铁路单边供电的处理方式相同.

图 9 树形供电方式示意

Figure 9. Tree-structured power supply

下载: 全尺寸图片幻灯片

我国西北地区地广人稀，电网对铁路牵引变电所多进行树形供电，为开展和实施树形双边供电提供了得天独厚的条件. 国铁集团于2020年在科技研究开发计划中立项开展“电气化铁路牵引变电所双边供电试验”研究，通过综合比选，考虑电网电源条件，试验段选取在海拔约3000 m、线路坡度16‰的格库铁路（青海段）东柴山至花土沟车站区段，长度约75 km. 格库铁路为国铁Ⅰ级客货共线的单线双向电气化铁路，设计速度120 km/h，使用HXD1C系列电力机车. 2021年7月试验取得圆满成功，仿真和试验证明在树形供电方式下，双边供电不产生穿越功率，具备投运条件.

交通运输部国家“交通强国”试点中的“巴准铁路贯通式同相供电工程化”项目也是在电网树形供电下的铁路贯通供电应用实例^[26]. 该铁路位于内蒙古鄂尔多斯市，为国铁Ⅰ级复线电气化重载运煤铁路，东起点岱沟站，向西经海勒斯壕南站与包神铁路在巴图塔站接轨，全长128 km. 改造前设有4座牵引变电所，贯通供电改造后由电网的川掌变电站树形供电至保留的四道柳和纳林川两个牵引变电所，牵引变电所均采用组合式同相供电方案^[11]：主变更换为不等边Scott接线，其M座额定容量为40 MV•A，T座同相供电装置采用5 MV•A运行 + 5 MV•A备用方式，2014年6月25日8时58分四道柳和纳林川变电所之间的分区所合环成功，2024年9月实现全线贯通，这是世界首条取消全部电分相和无电区的贯通供电线路.

2）合建所

将牵引变电所与电力（动照）变电所合建在一起，牵引变压器和电力（动照）变压器接入同一进线母线，如图10的SS₂所示. 假设穿越功率由变电所SS₁流向SS₂，只要电力变压器的电力负荷大于穿越功率，则穿越功率就被“淹没”，合建所总体对电网表现为用电状态，而不会向电网发电. 由于穿越功率通常较小，只需电力负荷大于该值，即可解决穿越功率向电网发电的问题.

图 10 牵引电力合建所示意

Figure 10. Co-built power and traction substation

下载: 全尺寸图片幻灯片

3）穿越功率利用

首先识别穿越功率，例如通过光纤实时传递双边供电两端牵引馈线的有功功率，设流向牵引网为正，当该功率非0且两者之和为0时，则判定该功率为穿越功率，再在穿越功率末端的牵引变电所牵引母线上将穿越功率通过交/直/交（AC/DC/AC）装置将穿越功率转移到三相电力用电母线上，供电力负荷利用，或者通过交/直（AC/DC）装置与储能装置（ES）连接进行储能利用，为了更好利用和控制储能装置容量，一般串接直/直（DC/DC）变换装置，如图11所示. 图中，TT为牵引变压器，PTa为电压互感器，CTa为电流互感器.

图 11 穿越功率利用示意

Figure 11. Through power utilization

下载: 全尺寸图片幻灯片

多数情况下，电网对铁路的供电方式是非树形的，这时就有穿越功率产生^[12-14]，例如，广州地铁18号和22号线采用同相供电，在陇枕变电所和陈头岗变电所进行了贯通供电试验，试验测试显示，穿越功率最大为2.2 MV•A.

5. 多功能的智能牵引变电所

与单边供电相比，贯通供电仍会产生负序功率和再生发电功率，但其最大值通常较小. 此外，贯通供电可能还会产生穿越功率. 由于再生发电功率和穿越功率的发电分量可以在牵引变电所进行统一处理，因此应提出一种多功能的牵引变电所供电方案，如图12所示. 图中：TSB为牵引母线；部分功率直流变换器（TPC）和三端口变换器（TPC）均为单相变流器，两者串联构成同相供电装置；DPC为三相交直变流器构成的能量转换装置，与TPC配合向铁路三相用电系统送电；ES为储能装置；RES为新能源发电装置；HPF为高通滤波器；PT为电压互感器；CT1和CT2为电流互感器；A、B、C为牵引变电所进行三相母线；a、b、c则为铁路三相用电系统的母线.

图 12 智能牵引变电所方案

Figure 12. Scheme of intelligent traction substation

下载: 全尺寸图片幻灯片

该方案借助单相组合式同相供电方案，仍然以Scott牵引变压器的M座牵引绕组为供电主通道，通过共用直流母线，增加新的潮流调节通道，最大程度复用变流器，减低变流器容量，增加性价比，使全寿命周期成本降至最低. 由于这种多功能供电方案可以根据供电需求进行灵活配置，同时实现牵引供电系统状态数据离线和在线智慧预测与监测，故亦可称为智能牵引变电所.

智能牵引变电所按照需要可以配置以下几个部分：1） M座牵引绕组构成单相牵引变压器，为供电主通道，简单、经济、可靠；2）同相供电装置，由TPC、PPC串联构成，按需配置，使负序达标；3）发电功率及穿越功率发电分量利用通道由TPC和DPC或ES构成，联通铁路三相用电或融入储能系统，能源高效融合利用，按需配置；4）新能源发电通道，按需配置，接入地点可在变电所或车站，实现铁路沿线新能源消纳；5）高通滤波通道，按需配置，消除高次谐波和谐振.

一般同相供电装置的容量仅占M座牵引绕组（主变）容量1/4甚至更小，穿越和再生功率利用通道容量也在5 MV•A及以下，可有效减少设备重复与投资.

智能牵引变电所具有故障导向安全特性，即使图12中的②～⑤通道退出运行，也能够保持牵引供电.

此外，借助单相组合式同相供电方案，还可以构造共用交流牵引母线的智能牵引变电所方案.

智能牵引变电所具有以下多种功能：1）可以补偿负序，使三相电压不平衡达标，不对电网造成影响；2）可以补偿无功功率，提高功率因数，提高和稳定网压；3）可以消除高次谐波和谐振；4）可以调控系统潮流，将剩余再生发电功率和穿越功率发电分量转变为用电功率，不向电网发电，不影响电网发电潮流调度；5）可以促进铁路沿线新能源消纳等.

以上功能可以根据工程实际需要进行灵活组合.

6. 结　论

1）在一般牵引变电所间距（不超过80 km）下，双边供电时分区所的合环电压差不大于16.00%，相角差不大于12.00°，符合电网合环规定.

2）与单边供电相比，双边供电和贯通供电可能会产生穿越功率. 该穿越功率可以通过计算和测量获得. 以110 kV输电线为例，其最大输送功率为50 MV•A，采用直供方式时，牵引网的最大穿越功率为1.8 MV•A，占电网输送功率的3.6%；而220 kV输电线的最大输送功率为500 MV•A，直供方式下牵引网的最大穿越功率为4.6 MV•A，占比为1.8%.

3）消除穿越功率向电网发电的影响，可以采取电网树形供电、合建变电所和穿越功率利用等方案加以解决.

4）可以进一步发展多功能的智能牵引变电所，根据实际需求进行灵活配置，综合解决负序、再生发电和穿越功率等对电网的影响问题，从而实现电气化铁路的更大范围贯通供电，推动电网和铁路的双赢发展.

致谢：黄小红、王辉、张戬、王帅博士参加了相关研究工作，在此表示感谢.

图 1 高速列车动力学性能参数闭环设计方法框架

Figure 1. Framework of closed-loop design method for dynamic performance parameters of high-speed train

下载: 全尺寸图片幻灯片

图 2 高速列车三维设计模型与分析模型的映射与参数提取

Figure 2. Mapping and parameter extraction of 3D design model and analysis model of high-speed train

下载: 全尺寸图片幻灯片

图 3 SOM神经网络结构

Figure 3. SOM neural network structure

下载: 全尺寸图片幻灯片

图 4 高速列车多工况下的多目标优化迭代思路

Figure 4. Iterative approach to multi-objective optimization under multiple working conditions for high-speed train

下载: 全尺寸图片幻灯片

表 1 高速列车典型服役工况设计列表

Table 1. Typical service conditions designed for high-speed train

工况	运行速度/ （km•h⁻¹）	曲线半径/ m	超高/ mm	轨道不平顺
工况一	300	∞		京津谱
工况二	250	5000	80	秦沈谱
工况三	300	7000	100	京津谱
工况四	350	9000	140	武广谱

下载: 导出CSV

表 2 高速列车设计空间缩减结果表

Table 2. Results of design space reduction of high-speed train

序号	参数符号	参数名称	初始取值范围	缩减后的范围
1	${x_1}$	车轮直径/mm	790～920	795～870
2	${x_2}$	轮对质量/kg	1200～2200	1800～2160
3	${x_3}$	轮对侧滚转动惯量/（kg•m²）	500～750	535～712
4	${x_4}$	轮对摇头转动惯量/（kg•m²）	500～800	522～764
5	${x_5}$	一系弹簧纵向刚度/（kN•m⁻¹）	800～1150	860～1080
6	${x_6}$	一系弹簧垂向刚度/（kN•m⁻¹）	900～1300	950～1230
7	${x_7}$	一系垂向阻尼/（kN•s•m⁻¹）	10～30	11.3～25.5
8	${x_8}$	轴箱转臂节点纵向刚度/（MN•m⁻¹）	5～15	6.5～14.2
9	${x_9}$	轴箱转臂节点横向刚度/（MN•m⁻¹）	4～10	4.2～9.5
10	${x_{10}}$	抗蛇行减振器横向跨距/mm	2400～2800	2450～2750
11	${x_{11}}$	空气弹簧垂向刚度/（kN•m⁻¹）	120～450	150～350
12	${x_{12}}$	空气弹簧横向刚度/（kN•m⁻¹）	150～400	160～350
13	${x_{13}}$	二系横向减振器节点刚度/（MN•m⁻¹）	10～50	12.0～45.6
14	${x_{14}}$	二系横向阻尼/（kN•s•m⁻¹）	30～65	34～60
15	${x_{15}}$	抗蛇行减振器节点刚度/（MN•m⁻¹）	5～13	7.2～12.5
16	${x_{16}}$	抗蛇行减振器阻尼（取非线性曲线的三点平均值）/（kN•s•m⁻¹）	50～3500	75～3200

下载: 导出CSV

表 3 优化后的随机森林模型的超参数值

Table 3. Hyperparametric values of optimized random forest model

	树的棵数/棵	拆分内部节点最小样本数/个	叶子节点最小样本数/个	最大深度
工况一	127	2	11	9
工况二	132	3	15	9
工况三	326	2	10	9
工况四	271	4	7	12

下载: 导出CSV

表 4 各工况贝叶斯优化超参数的随机森林代理模型综合平均值

Table 4. Comprehensive average value of random surrogate forest model hyperparameters optimized by Bayesian method under each working condition

	工况一	工况二	工况三	工况四
MAE	0.2324	0.2864	0.3533	0.3924
MSE	0.3527	0.4292	0.3976	0.4935
拟合系数	0.9981	0.9968	0.9945	0.9883

下载: 导出CSV

表 5 4种工况性能优化对比

Table 5. Comparison of performance optimization under four working conditions

性能指标	工况一		工况二		工况三		工况四
性能指标	优化前	优化后	优化前	优化后	优化前	优化后	优化前	优化后
横向平稳性	1.76	1.74	1.95	1.93	1.83	1.80	1.79	1.75
垂向平稳性	1.88	1.82	1.98	1.95	1.86	1.82	1.92	1.90
轮轨垂向力/kN	77.56	75.34	83.46	81.82	79.67	77.72	78.82	77.94
轮轴横向力/kN	6.52	6.37	12.65	11.45	10.53	9.83	7.75	7.23
脱轨系数	0.12	0.11	0.23	0.21	0.18	0.17	0.15	0.14
轮重减载率	0.36	0.35	0.41	0.4	0.38	0.37	0.37	0.36
倾覆系数	0.37	0.34	0.40	0.38	0.38	0.36	0.38	0.36

下载: 导出CSV

参考文献(26)

[1]	石怀龙,罗仁,曾京. 国内外高速列车动力学评价标准综述[J]. 交通运输工程学报,2021,21(1): 36-58 SHI Huai-long, LUO Ren, ZENG Jing. Review on domestic and foreign dynamics evaluation criteria of high-speed train[J]. Journal of Traffic and Transportation Engineering, 2021, 21(1): 36-58.
[2]	安东,邹益胜,赵春发,等. 高速磁浮列车动力学性能参数多目标优化方法研究[J]. 机械科学与技术,2022,41(03): 466-472 AN Dong, ZOU Yishen, ZHAO Chun Fa, et al. Study on multi-objective optimization method of dynamic performance parameters of high-speed maglev train[J]. Mechanical Science and Technology for Aerospace Engineering, 2022, 41(03): 466-472.
[3]	姚远,陈相旺,李广,等. 高速列车抗蛇行减振器参数的多目标优化研究[J]. 西南交通大学学报,2021,56(06): 1298-1304 YAO Yuan, CHEN Xiangwang, LI Guang, et al. Multi-objective optimization of yaw damper parameters for high-speed train[J]. Journal of Southwest Jiaotong University, 2021, 56(06): 1298-1304.
[4]	CHEN X W, YAO Y, SHEN L J, et al. Multi-objective optimization of high-speed train suspension parameters for improving hunting stability[J]. International Journal of Rail Transportation, 2022, 10(2): 159-176. doi: 10.1080/23248378.2021.1904444
[5]	Zhu Y, Lyu Y, Olofsson U. Mapping the friction between railway wheels and rails focusing on environmental conditions[J]. Wear, 2015, 324-325: 122-128. doi: 10.1016/j.wear.2014.12.028
[6]	宫金良,胡光学,张彦斐. 以刚度为目标的微位移放大模块闭环设计方法[J]. 机械工程学报,2012,48(15): 58-64 GONG Jinliang, HU Guangxue, ZHANG Yanfei. Closed loop design method of micro-driving displacement amplifier module targeting for stiffness[J]. Journal of Mechanical Engineering, 2012, 48(15): 58-64.
[7]	罗钜,郭福永. 民用涡扇发动机总体多学科设计优化研究[J]. 中国机械工程,2019,30(15): 1813-1820 doi: 10.3969/j.issn.1004-132X.2019.15.007 LUO Ju, GUO Fuyong. Study on multi-disciplinary design optimization on civil turbofan engines[J]. China Mechanical Engineering, 2019, 30(15): 1813-1820. doi: 10.3969/j.issn.1004-132X.2019.15.007
[8]	ZHANG W H, Shen Z Y, Zeng J. Study on dynamics of coupled systems in high-speed trains[J]. Vehicle System Dynamics, 2013, 51(7): 966-1016. doi: 10.1080/00423114.2013.798421
[9]	Angelo G, Marzolla M. New trends in parallel and distributed simulation: From many-cores to Cloud Computing[J]. Simulation Modelling Practice and Theory, 2014, 49: 320-335. doi: 10.1016/j.simpat.2014.06.007
[10]	姚远,李广,梁树林,等. 基于健壮稳定性的高速列车悬挂参数优化匹配方法[J]. 铁道学报,2021,43(8): 35-44 doi: 10.3969/j.issn.1001-8360.2021.08.005 YAO Yuan, LI Guang, LIANG Shulin, et al. Optimal matching method for suspension parameters of high-speed train bogie based on robust hunting stability[J]. Journal of the China Railway Society, 2021, 43(8): 35-44. doi: 10.3969/j.issn.1001-8360.2021.08.005
[11]	姚远,任铖铭,陈相旺,等. 基于频域平稳性的高速机车悬挂参数优化匹配[J]. 西南交通大学学报,2022,57(6): 1259-1267 YAO Yuan, REN Cheng Ming, CHEN Xiang Wang, et al. Suspension parameters optimum matching of high-speed locomotive based on frequency domain stationarity[J]. Journal of Southwest Jiaotong University, 2022, 57(6): 1259-1267.
[12]	YAO Y, CHEN X W, LI H, et al. Suspension parameters design for robust and adaptive lateral stability of high-speed train[J]. Vehicle System Dynamics, 2023, 61(4): 943-967. doi: 10.1080/00423114.2022.2062012
[13]	姚远,陈相旺,李广,等. 高速列车抗蛇行减振器参数的多目标优化研究[J]. 西南交通大学学报,2021,56(6): 1298-1304 YAO Yuan, CHEN Xiangwang, LI Guang, et al. Multi-objective optimization of yaw damper parameters for high-speed train[J]. Journal of Southwest Jiaotong University, 2021, 56(6): 1298-1304.
[14]	于日伟,周长城,赵雷雷. 高速列车转向架-车体-座椅垂向耦合振动机理及悬挂参数联合优化[J]. 机械工程学报,2018,54(8): 57-67 doi: 10.3901/JME.2018.08.057 YU Riwei, ZHOU Changcheng, ZHAO Leilei. Vertical coupled vibration mechanism of bogie-body-seat system and joint optimization of suspension parameters for high-speed train[J]. Journal of Mechanical Engineering, 2018, 54(8): 57-67. doi: 10.3901/JME.2018.08.057
[15]	解欢,杨岳,童林军,等. 基于混合代理模型的高速轨道车辆悬挂参数多目标优化[J]. 铁道科学与工程学报,2016,13(10): 2056-2063 doi: 10.3969/j.issn.1672-7029.2016.10.025 XIE Huan, YANG Yue, TONG Linjun, et al. Multi-objective optimization of the suspension parameters for high speed rail vehicle based ona hybrid surrogate model[J]. Journal of Railway Science and Engineering, 2016, 13(10): 2056-2063. doi: 10.3969/j.issn.1672-7029.2016.10.025
[16]	姜杰,杨旭锋,丁国富. 基于多峰优化Kriging模型与距离相关系数的高速列车动力学参数多输出灵敏度分析[J]. 工程科学与技术,2024,56(4): 250-260 JIANG Jie, YANG Xufeng, DING Guofu. Sensitivity analysis for dynamic parameters of high-speed train based on multimodal-optimization improved kriging model and distance correlation[J]. Advanced Engineering Sciences, 2024, 56(4): 250-260.
[17]	许平,杨丽婷,姚曙光,等. 城轨列车方锥式防爬吸能结构碰撞力学参数设计及多目标优化[J]. 中南大学学报(自然科学版),2022,53(5):1689-1699. XU Ping,YANG Liting,YAO Shuguang,et al. Collision mechanics parameter design and multi-objective optimization of square cone anti-climbing energy-absorbing structure for urban rail trains[J]. Journal of Central South University(Science and Technology), 2022, 53(5):1689-1699.
[18]	张剑. 基于代理模型技术的高速列车性能参数设计及优化[D]. 成都:西南交通大学博士学位论文,2015 ZHANG Jian. The high-speed train performance parameter design and optimization based on surrogate model technology[D]. Chengdu: Southwest Jiaotong University. College of Mechanical Engineering, 2015.
[19]	何庆,利璐,李晨钟,等. 基于Kriging模型的在役高速列车悬挂参数近似贝叶斯估计[J]. 机械工程学报,2023,59(12): 139-148 doi: 10.3901/JME.2023.12.139 HE Qing, LI Lu, LI Chenzhong et al. Approximate bayesian estimation of suspension parameters of in-service high-speed trains based on kriging surrogate model[J]. Journal of Mechanical Engineering, 2023, 59(12): 139-148. doi: 10.3901/JME.2023.12.139
[20]	肖乾,罗超,欧阳志许,等. 基于RBF神经网络代理模型的车辆/轨道参数多目标优化[J]. 机械强度,2021,43(2): 319-326 XIAO Qian, LUO Chao, OUYANG Zhixu, et al. Multi-objective optimization of vehicle/track parameters based on RBF neural network surrogate model[J]. Journal of Mechanical Strength, 2021, 43(2): 319-326.
[21]	饶坝,马术文,张海柱,等. 高速列车转向架动力学设计指标分解冲突博弈协调方法[J]. 铁道科学与工程学报,2024,21(4): 1334-1344 RAO Ba, MA Shuwen, ZHANG Haizhu, et al. Game coordination method for decomposition conflict of dynamic design indicators of high-speed train bogie[J]. Journal of Railway Science and Engineering, 2024, 21(4): 1334-1344.
[22]	Johnsson A, Berbyuk V, Enelund M. Pareto optimization of railway bogie suspension damping to enhance safety and comfort[J]. Vehicle System Dynamics, 2012, 50(9):1379-1407.
[23]	成棣,孙琛,胡晓依,等. 面向低锥度晃车的CRH3平台动车组车轮型面优化[J]. 中国铁道科学,2020,41(6): 135-144 doi: 10.3969/j.issn.1001-4632.2020.06.15 CHENG Di, SUN Chen, HU Xiaoyi, et al. Wheel profile optimization of CRH3 EMU oriented to carbody shaking caused by low equivalent conicity[J]. China Railway Science, 2020, 41(6): 135-144. doi: 10.3969/j.issn.1001-4632.2020.06.15
[24]	Zou Yi sheng, Ding Guo fu, Zhang Wei Hua, . et al. Collaborative simulation method with spatiotemporal synchronization process control[J], Chinese Journal of Mechanical Engineering, 2016, 29 (6): 1074-1082
[25]	王伟. 轨道不平顺对地铁车辆动力学性能的影响分析[D]. 成都:西南交通大学硕士学位论文,2014 WANG Wei. Analysis of metro vehicle dynamics performance on track irregularity[D]. Chengdu: Southwest Jiaotong University, 2014.
[26]	Liaw A, Wiener M. Classification and regression by random Forest[J]. R news, 2002, 2(3): 8-22.

施引文献

附加材料(0)

访问统计

点击查看大图

图(4) / 表(5)

计量

文章访问数: 39
HTML全文浏览量: 22
PDF下载量: 17
被引次数: 0

1. 贯通供电
2. 贯通供电合环
3. 穿越功率分析与计算
3.1 等效电路与均衡电流计算
3.2 穿越功率计算与测量
4. 穿越功率解决方案
5. 多功能的智能牵引变电所
6. 结　论

高速列车动力学性能参数优化迭代设计

doi: 10.3969/j.issn.0258-2724.20220071

作者简介: 姜杰（1988—），男，博士研究生，研究方向为复杂产品概念设计与多学科优化，E-mail：jj13688158635@126.com

通讯作者: 丁国富（1972—），男，教授，博士，研究方向为数字化设计与制造，E-mail：dingguofu@swjtu.edu.cn

计量

出版历程