DOI： 10.11817/j.issn.1672-7207.2019.11.029

基于改进粒子群算法的机车二系弹簧载荷分配优化

韩锟¹,杨广为¹,黄泽帆¹,鲍天哲²

（1. 中南大学 交通运输工程学院，湖南 长沙，410075；

2. 利兹大学 电子电气学院，英国 利兹，LS2 9JT）

摘 要：

簧载荷优化调整这一复杂的多变量优化问题，为进一步提高现有求解方法的优化效果和计算效率，将烟花算法融入粒子群算法，提出一种具有分层递阶结构的改进粒子群算法，算法为3层架构，其中，底层是基础层，为加入维变异算子的粒子群算法，是改进算法的基本框架；中间层是融合层，为引入烟花算法爆炸机制的粒子更新层，主要用于扩大算法搜索范围，提高全局搜索能力；顶层是扰动层，引入扰动因子，避免算法因陷入局部搜索而进行的大量无为冗余迭代，加快全局收敛速度。用典型测试函数对改进算法性能进行测试，并将其应用于机车二系弹簧载荷分配优化调整仿真实验。研究结果表明：改进算法与传统遗传算法、烟花算法和粒子群算法相比，全局搜索能力更强，鲁棒性更好，求解精度更高。

关键词:机车二系弹簧；载荷分配优化；改进粒子群算法；分层递阶结构；改进烟花算法

中图分类号:U260.72;TP301    文献标志码:A     文章编号:1672-7207（2019）11-2906-09

Optimization of locomotive secondary spring loads distribution based on improved particle swarm optimization algorithm

HAN　Kun¹, YANG　Guangwei¹, HUANG　Zefan¹, BAO　Tianzhe²

(1. School of Traffic & Transportation Engineering, Central South University, Changsha 410075, China；

2. School of Electronic and Electrical Engineering, University of Leeds, Leeds LS2 9JT, UK)

Abstract: The optimization of spring loads distribution can be summarized as a complicated nonlinear multidimensional optimization problem. In order to further improve the optimization effect and calculation efficiency of the existing solution methods, an improved particle swarm optimization algorithm(IPSO) with hierarchical structure was proposed by integrating the fireworks algorithm(FWA) into a particle swarm optimization algorithm(PSO). The IPSO is a three-layer architecture. The underlying layer was the basic framework of the IPSO and a classical PSO that adopted a dimensional mutation operator. The middle layer was a fusion layer used to expand the search range of IPSO. In essence, it was a particle update layer that introduced the explosion mechanism in FWA. The top layer was the disturbance layer, with a disturbance factor introduced to avoid redundant iteration due to local search and speed up the global convergence. The IPSO was tested on typical test function and applied to the simulation experiment of the distribution optimization and adjustment of the locomotive secondary spring loads. The results show that compared with traditional genetic algorithm(GA), FWA and PSO, the IPSO has stronger global search ability, better robustness and higher accuracy.

Key words: locomotive secondary springs; load distribution optimization; improved particle swarm algorithm; hierarchical structure; improved fireworks algorithm

机车轮(轴)重分配不均会使机车运行时轮轴间传递的最大切向力减小，造成机车牵引力下降，甚至会导致空转打滑、车轮踏面剥离、车轴断裂等事故，是影响行车安全的重要隐患^[1]。对具有两系悬挂结构的机车，通过加垫调节一、二系弹簧的变形量以改善轮(轴)重分配(以下简称调簧)，是减小甚至消除轮(轴)重分配偏差的有效方法^[2]。传统的人工调簧方法完全依赖经验判断，调整结果难以满足GB/T 3317—2016规定的轮(轴)重分配允许偏差要求。针对调簧这一复杂的多变量优化问题，在优化模型的表达方面，CURTIS等^[3]采用试验法，通过试验测试得到加垫高度与弹簧载荷分布间的关系，作为计算加垫高度的依据。该方法易于实现，但计算结果的准确性过度依赖试验数据，不具有通用性。潘迪夫等^[4]基于经典力学理论，得到了描述二系弹簧加垫高度与弹簧载荷分布间映射关系的调簧理论模型，具有一定的通用性。在此基础上，韩锟等^[5]提出了带误差补偿的混合模型，进一步完善了调簧理论模型，使模型包含了机车弹性悬挂结构中非线性因素对弹簧载荷分布的影响。上述成果为调簧提供了重要的理论依据。在优化模型的求解方面，近年来，各种智能优化方法如遗传算法、SUMT算法、人工鱼群算法、模拟退火算法等不断应用于调簧领域^[6-8]，这些算法理论上可以求解，但通常存在计算时间长、求解精度低等问题。为此，韩锟等^[9-10]提出了改进遗传算法，在一定程度避免了遗传算法早熟收敛的问题，求解精度有所提高。此外，潘迪夫等^[11-12]提出了多目标遗传算法以及免疫优化算法，在提高算法鲁棒性的同时，有效控制了加垫总高度，但算法复杂度也随之提高。针对上述问题，本文提出一种具有3层递阶结构的改进粒子群算法(improved particle swarm optimization，IPSO)，在充分利用粒子群算法参数设计简单、收敛速度快等优点的基础上对粒子的更新机制加以改进，通过引入维变异算子和烟花爆炸机制增加粒子的多样性，通过引入扰动因子避免算法因陷入局部搜索而进行大量冗余迭代，从而在提高算法全局搜索能力的同时加快算法的收敛速度，使调簧精度和效率共同提高。

1 改进的粒子群算法(IPSO)

1.1　粒子群算法改进的必要性以及改进粒子群算法总体思想

粒子群算法是20世纪末 KENNEDY等^[13-14]受鸟群合作觅食启发提出的一种智能进化算法。该算法通过模拟鸟群个体间对食物信息的共享，利用不断围绕距离食物最近的飞鸟展开搜索，逐渐缩小信息范围直到发现食物。粒子群算法原理简单，参数设置简单，收敛速度快。与遗传算法复杂的编码解码过程相比，粒子群算法更易实现，并具有良好的通用性和鲁棒性，已逐渐成为人们的研究热点。但传统粒子群算法容易发生早熟现象，其收敛精度还有待提高。为此，本文提出一种具有3层架构的分层递阶改进粒子群算法，如图1所示：底层是基础层，是改进算法的基本框架，其本质是加入了维变异算子的粒子群算法；中间层是融合层，在粒子群算法中引入烟花算法(fireworks algorithm，FWA)^[15-16]的爆炸机制，将粒子群算法中适应度较高的粒子作为烟花算法的爆炸烟花，以此增加粒子群算法的种群多样性，扩大搜索范围，提升算法的全局搜索能力，防止粒子群算法陷入局部最优，提高算法的收敛精度；顶层为扰动层，对粒子中的个体极值和全局极值施加扰动因子，避免算法因陷入局部搜索而进行大量冗余迭代，进一步加快收敛速度，提高算法效率。

图1　改进粒子群算法总体架构

Fig. 1　Framework of improved particle swarm optimization

1.2　IPSO算法底层的维变异算子

粒子群算法(particle swarm optimization，PSO)没有充分的精英机制，在淘汰迭代中产生劣解时，算法收敛速度会减慢。本文参考文献[17]，引入带维变异算子的PSO作为IPSO算法的基本框架，其变异策略为：在粒子群所在的D维搜索空间中，引入参数l_con(n)(n∈[ 1, D])表示各个粒子在第n维的当前位置(i=1, 2, , N，N为粒子群中的粒子总数)到粒子群在第n维上的质心D_n到的距离总和，可用下式表示：

图2　改进爆炸机制前后烟花爆炸结果示意图

Fig. 2　Diagram of fireworks explosion before and after improving explosion mechanism

(1)

l_con(n)描述了粒子群在第n维上的收敛程度，其值越小，表示粒子群在该维上的收敛程度越大。在粒子群算法的每一次迭代过程中，计算每个维度上的l_con(n)，找到这D个l_con(n)中最小值所在的维d，并对各个粒子在该维的当前位置以一定的概率按下式进行变异操作，使让所有粒子在该维上的位置重新均匀分布：

(2)

式中：α为(0,1)区间的随机数；x_max为所有粒子在d维的位置最大值。

1.3　IPSO算法中间层的烟花爆炸机制

为进一步提高全局和局部搜索能力，IPSO算法构建了融合烟花算法的中间层，对烟花算法的爆炸机制加以改进，用于更新粒子。在烟花算法中，每一个烟花个体都代表优化问题的1个候选解，算法在执行爆炸寻优操作时，以每个烟花作为爆炸中心，通过特定的爆炸策略生成子代火花，并在其中挑出优解作为下一代的候选解。IPSO算法参考文献[18-19]，改进的烟花爆炸机制如下。

在以烟花个体x_i为中心，半径为r的超立方体(D维空间)中随机产生M_g个火花的过程中，x_i先按式(3)进行处理：

(3)

式中：R为爆炸半径在坐标轴上的投影值；β为服从[0, 1]均匀分布的随机数；为方向向量。

从式(3)可以看出，改进后的爆炸机制使烟花不再仅沿坐标轴以固定半径爆炸，而且其爆炸方向和爆炸半径均随随机数β而变化，如图2所示。由图2可见：改进的爆炸机制丰富了由烟花爆炸得到的火花的多样性，使算法的搜索路径得到改善，不再仅沿坐标轴方向搜索，而且可以沿着任意角度搜索，拓展了算法在设计区间的搜索范围。

烟花算法的爆炸强度即产生的新火花的数目与烟花适应度紧密联系，通过控制爆炸强度使适应度高的烟花爆炸而生成更多新火花，可以提高寻优效率。减小爆炸幅度(爆炸半径)可以缩小烟花的爆炸范围，加快算法收敛。烟花算法是通过计算每个烟花个体的适应度函数值确定爆炸强度和爆炸幅度，计算方法如下：

(4)

(5)

式中：S_i为第i个烟花爆炸产生的火花数量；A_i为第i个烟花的爆炸半径；y_max和y_min分别为当前烟花种群适应度的最大值和最小值；f(x_i)为适应度函数，i=1, 2, , N；N为种群规模； 为最大爆炸半径；M为常数，用于控制N个烟花爆炸产生的火花总数。为避免在特殊情况下存在分母为0的可能，引入极小常数进行约束。

为了避免适应度高的烟花爆炸强度过大，产生压倒性的影响，通过下式对烟花爆炸产生的火花数量进行约束：

(6)

式中：为约束后第i个烟花爆炸产生的火花数量；a和b为常数；round为取整函数。

1.4　IPSO算法原理

如前所述，IPSO算法以带维变异算子的粒子群算法作为算法主体。但随着迭代次数增加，带维变异算子的粒子群算法会因粒子贫化而导致搜索能力下降，故IPSO算法在每次迭代过程中进行粒子更新操作时，引入1.3节所述的改变爆炸方式的烟花算法，将带维变异算子的粒子群算法更新得到的每个粒子个体极值作为烟花算法中的烟花个体，按式(3)给出的爆炸机制进行操作，以期增强IPSO算法在迭代过程中的种群多样性，提高算法全局搜索能力。爆炸操作中每个烟花个体产生的火花数量S_i和爆炸半径A_i分别按式(7)和式(8)计算：

(7)

(8)

式中：为第i个粒子第t次迭代时的个体极值的适应度；和分别为所有烟花个体中的适应度最大值和最小值。

最后，在IPSO算法中引入扰动因子，防止算法在寻优过程中因陷入局部搜索而出现扎堆现象，从而加快收敛速度。在第t次迭代过程中，对个体极值，全局极值以及迭代更新后粒子位置加入扰动因子，即

(9)

(10)

(11)

其中：，和为条件随机函数；和分别为个体极值和全局极值迭代停滞的次数；为对粒子位置进行扰动的次数；为更新后的粒子速度；，和分别为个体极值、全局极值和粒子位置扰动的限界值。

由式(12)～(16)计算得t+1次迭代更新后d维粒子的速度和位置。设为迭代次数t时的d维第i个粒子的个体极值，为迭代次数t时的d维第i个粒子的全局极值，t为迭代次数，T为最大迭代次数，为第t次迭代时的惯性权重，和分别为的最大值和最小值，和为学习因子(取值不大于2)，和分别为粒子的速度最大值和最小值，和分别为位置上、下限，有

                   (12)

(13)

(14)

(15)

(16)

IPSO算法的具体流程如下。

Step 1:设定种群规模、迭代次数、惯性系数、速度位置等算法参数。

Step 2:初始化粒子的速度和位置。

Step 3:计算粒子的适应度，再求出个体极值对应烟花算法的适应度。

图3　4种算法对基准测试函数的收敛曲线

Fig.3　Convergence curves of four algorithms running on benchmark function

Step 4:令t=t+1，开始迭代更新粒子的速度和位置。按照式(13)～(16)更新粒子的速度，再由速度更新位置。

Step 5:加入维变异算子，对粒子进行变异，计算适应度函数，依据种群粒子信息更新。

Step 6:将Step 5更新得到的作为烟花算法中的烟花个体，并对其进行爆炸操作，计算，并更新。

Step 7:根据与更新全局极值。

Step 8:引入扰动因子，计算适应度函数值，并再次更新。

Step 9:若t

2 算法性能测试

为了验证IPSO算法的性能，将IPSO算法和遗传算法(GA)、粒子群算法(PSO)及烟花算法(FWA)进行比较。选取3个最小值唯一且都为0的标准测试函数对算法性能进行测试。其中，为Rosenbrock函数，它是一个结构特征酷似峡谷的单峰函数，其全局最小值位于绵长且平坦的谷底，该函数可用于测试算法的全局搜索能力。

(17)

F₂(X)为Rastrigin函数，是1个局部极小值较多的多峰函数，可以测试算法跳出局部最优的能力：

(18)

F₃(X)为Ackleys函数，用来测试算法全局收敛速度：

(19)

为了比较算法的性能优势，IPSO，GA，PSO和FWA算法中均取最大迭代次数T=500，种群规模N=100。Rosenbrock函数对应维度为D=10，其他2种函数对应维度D=35。图3所示为各算法迭代过程中最佳适应度收敛曲线。

从图3可见：在这4种算法中，IPSO算法的搜索能力最强，收敛速度最快。这是因为IPSO算法充分利用了烟花算法的爆炸机制来增强种群多样性，降低了算法限入局部极值的概率；同时，引入的扰动因子避免了IPSO算法在局部极值区域的冗余迭代，加快了算法收敛速度。

为进一步验证算法的鲁棒性，这4种算法分别独立运行50次，得到的最优解、最差解、最优适应度均值、方差如表1所示。

表1　4种算法独立运行50次的测试结果

Table 1　Test results of four algorithms running independently 50 times

表1表明：在用4种算法对3个测试函数分别进行寻优运算的实验中，IPSO算法所得到的500次独立运行结果中的最优解、最差解，以及500次求解结果的均值和方差等均低于GA，PSO和FWA算法相应的求解结果，表明IPSO算法具有更强的全局搜索能力和鲁棒性。

3 IPSO算法在机车二系调簧中的应用

3.1　机车二系支承结构几何模型

本文实验平台为实验室自主研发设计的大功率机车车体称重调簧试验台^[4]，本实验所有数据均来源于现场采集，试验台为车体提供12个支承点，通过液压系统控制支承点抬升来模拟各二系弹簧处的加垫操作。试验台的位移控制精度达0.1 mm，足以保证试验的可靠性。

轴式为C₀-C₀的电传动机车二系支承结构的简化模型见图4。此模型只考虑二系弹簧垂向载荷分布，将车体视为一个整体，总重力为P，和分别为车体重心在X和Y方向的偏移量，a和b分别为沿X和Y方向距离车底中心O的距离，S₁和S₂为同一转向架同侧的二系弹簧横向间距，F₁～F₁₂分别为12个二系弹簧对车体的支承力。

图4　C₀-C₀轴式机车二系支承结构简化模型

Fig. 4　Simplified geometry model of secondary suspension system of C₀-C₀ locomotive

3.2　目标函数

文献[5-6]针对机车二系弹簧载荷分布建立了基于传统机理模型的优化模型，具体形式为

(20)

式中：目标函数J表示二系弹簧载荷分布均方差，J越小，表示载荷分布越均匀；和分别为加垫调整前后第j个二系弹簧的载荷；为所有二系弹簧的载荷平均值；n为机车二系弹簧数目，从图4可知n=12；为加垫引起载荷变化的增载系数；为第i个弹簧对应的加垫高度；为常数，为试验允许的最大加垫高度。

3.3　仿真试验结果

以大功率HXD1D型机车为试验对象，根据现场实测数据结合IPSO算法开展调簧仿真实验。为避免参数选择的盲目性，利用遗传算法对IPSO算法的参数进行全局优化，综合考虑计算复杂度和参数可靠度，得出一组最优参数组合：种群规模N=60，迭代次数G=200次，惯性权重最大值w_max=0.9，惯性权重最小值，学习因子，最大允许加垫高度=10 mm。4种算法仿真实验过程中的迭代收敛曲线如图5所示。

图5　4种算法的迭代收敛曲线

Fig. 5　Convergence curves of four algorithms

由图5可知：在迭代前期，IPSO算法的收敛速度明显比GA，PSO和FWA这3种算法的快；随着迭代次数增加，PSO和FWA算法趋于早熟收敛，而IPSO算法在迭代超过60次以后载荷方差还有明显下降，表明将粒子群算法与烟花算法相结合的IPSO算法很好地克服了传统算法早熟的问题；同时，在相同种群规模下，IPSO算法的收敛精度最高，具有较强的全局搜索能力。

HXD1D型机车共有12个二系支承点，IPSO算法优化前后，各支承点处二系弹簧的载荷分布情况和相应的加垫高度如表2所示。由表2可知：通过IPSO算法优化后，二系弹簧载荷分布方差由初始的1.300 降低到0.507 ，减小了61%；12个二系弹簧载荷的最大差值也由2.68 kN减小到2.03 kN，减小了24.3%。实验结果验证了IPSO算法优化二系弹簧载荷分布均匀性的有效性。

表2　HXD1D型机车初始二系载荷分布及优化结果

Table 2　Initial distribution and optimization results of secondary spring load of locomotive HXD1D

为了进一步验证IPSO算法在机车二系弹簧载荷分布优化中的有效性和鲁棒性，将其与遗传算法、粒子群算法、烟花算法进行比较，4种算法在种群规模为60，迭代次数为200次的条件下独立运行50次的结果见表3。从表3可以看出：IPSO算法的平均运行时间比GA略短，比PSO和FWA耗时要长；但IPSO算法50次求解得到的二系弹簧载荷分布方差的平均值为0.507 8，较PSO算法的0.744 9减小了31.8%，较FWA算法的0.633 5减小了19.8%，较GA算法的0.644 9减小了21.3%，在求解精度上具有显著优势；IPSO算法50次求解得到的二系弹簧载荷分布方差的标准差仅为3.113 5×10^-4，50次求解得到的二系弹簧载荷分布方差的最优解与最差解仅相差0.001 2，显示了良好的鲁棒性；与GA算法相比，IPSO算法在求解精度、求解效率、鲁棒性等方面均具有一定优势。

表3　4种算法独立运行50次的结果对比

Table 3　Results comparison of four algorithms independently running 50 times

PSO，FWA，GA和IPSO算法分别独立运行50次所得结果分布的箱式图见图6。图6中每个箱体的高度在一定程度上反映了50次结果数据的波动程度，箱体上方和下方的线段代表了该组数据的最大值和最小值，箱体内部的线代表该组数据中位数的值，箱子外部的点为数据中位于箱体外面且距离相应边界大于1.5倍箱体高度的“异常值”。从图6可以直观地看出：描述IPSO算法50次运行结果分布的箱体高度最小，表明该算法多次运行结果的波动远比其他3种算法的小，且50次运行结果中没有“异常值”，具有很好的鲁棒性；IPSO算法50次运行结果的中位数值在所有算法中最小，表明该算法求解精度最高。

图6　4种算法50次独立运行结果分布的箱式图

Fig. 6　Box-plot of four algorithms

4 结论

1) 提出了一种采用分层递阶结构的改进粒子群算法IPSO，算法底层为基础层，通过在粒子群算法中加入维变异算子，避免算法早熟收敛；中间层为融合层，融入烟花算法，采用改进的烟花爆炸机制进行粒子更新，充分利用烟花算法的优势平衡IPSO算法的全局和局部搜索能力；顶层为扰动层，在一定条件下对每次迭代过程中得到粒子的个体极值和全局极值进行加入扰动因子，以避免算法在局部极值区域进行冗余迭代，提高算法效率。

2) IPSO算法在全局搜索能力、求解精度、鲁棒性等方面的综合性能明显优于其他几种优化算法。

3) 经过IPSO算法优化后的载荷均方差相较于原始值减小了61%，改进效果明显；在算法独立运行50所得二系弹簧载荷分配方差的平均值与遗传算法、烟花算法和粒子群算法相比，分别减小21.3%，19.8%和31.8%，且50次运行结果的波动幅度远小于这些算法，充分证明了IPSO算法在求解机车二系弹簧载荷分配优化问题的有效性和可靠性。

参考文献：

[1] 潘迪夫, 陈军, 鲍天哲, 等. 多输出支持向量机混合模型在机车调簧中的应用[J]. 华南理工大学学报(自然科学版), 2018, 46(4): 22-28.

PAN Difu, CHEN Jun, BAO Tianzhe, et al. Application of multi-output support vector regression hybrid model in locomotive secondary spring loads adjustment[J]. Journal of South China University of Technology(Natural Science Edition), 2018, 46(4): 22-28.

[2] BAO Tianzhe, HAN Kun, PAN Difu. Method for secondary spring load equalization of railway vehicles with two-stage spring suspension: Modeling and optimal regulation[J]. Journal of Central South University, 2018, 25(4): 936-948.

[3] CURTIS D L, SKRZYPEZYK W G, THOMAS T J. Method of adjusting the distribution of locomotive axle loads: United States, 4793047[P]. 1988-12-27.

[4] 潘迪夫, 韩锟, 曾亚波, 等. 车体称重调簧试验装置及其应用[J]. 电力机车与城轨车辆, 2003, 26(5): 37-39.

PAN Difu, HAN Kun, ZENG Yabo, et al. Locomotive secondary spring load test device and its application[J]. Electric Locomotives & Mass Transit Vehicles, 2003, 26(5): 37-39.

[5] 韩锟, 潘迪夫. 机车二系弹簧载荷均匀性分配调整的混合建模方法[J]. 中南大学学报(自然科学版), 2012, 43(1): 378-383.

HAN Kun, PAN Difu. Hybrid modeling method for adjusting distribution of locomotive secondary spring loads[J]. Journal of Central South University(Science and Technology), 2012, 43(1): 378-383.

[6] 潘迪夫, 黎航, 韩锟. 基于遗传算法的机车二系支承载荷调整优化方法[J]. 中国铁道科学, 2005, 26(3): 83-87.

PAN Difu, LI Hang, HAN Kun. Optimization model of locomotive secondary spring load adjustment based on genetic algorithm[J]. China Railway Science, 2005, 26(3): 83-87.

[7] 杨本磊, 潘迪夫. 基于人工鱼群算法的机车二系支承载荷调整优化方法[J]. 计算机与现代化, 2011(1): 52-54, 120.

YANG Benlei, PAN Difu. Optimization model of locomotive secondary spring load adjustment based on artificial fish-swarm algorithm[J]. Computer and Modernization, 2011(1): 52-54, 120.

[8] 高久淳, 林建辉. 基于模拟退火算法的地铁车辆调簧优化算法设计[J]. 机械, 2014, 41(10): 29-32.

GAO Jiuchun, LIN Jianhui. Optimization method of spring adjustment model based on simulated annealing algorithm for metro vehicles[J]. Machinery, 2014, 41(10): 29-32.

[9] 韩锟, 潘迪夫, 韩洪飞. 机车车辆二系弹簧载荷分配优化的改进遗传算法[J]. 中南大学学报(自然科学版), 2016, 47(7): 2521-2527.

HAN Kun, PAN Difu, HAN Hongfei. An improved genetic algorithm for secondary spring load equalization of railway vehicles[J]. Journal of Central South University(Science and Technology), 2016, 47(7): 2521-2527.

[10] 李家才, 韩锟, 鲍天哲. 基于Bloch球面坐标的改进量子遗传算法及其应用[J]. 铁道科学与工程学报, 2016, 13(11): 2262-2269.

LI Jiacai, HAN Kun, BAO Tianzhe. An improved Bloch spherical quantum genetic algorithm and its application[J]. Journal of Railway Science and Engineering, 2016, 13(11): 2262-2269.

[11] 潘迪夫, 朱亚男. 基于多目标遗传算法的机车二系支承载荷调整优化方法[J]. 铁道科学与工程学报, 2011, 8(2): 76-80.

PAN Difu, ZHU Yanan. Optimization model for locomotive secondary spring load adjustment based on multi-objective genetic algorithm[J]. Journal of Railway Science and Engineering, 2011, 8(2): 76-80.

[12] PAN Difu, WANG Mengge, ZHU Yanan, et al. An optimization algorithm for locomotive secondary spring load adjustment based on artificial immune[J]. Journal of Central South University, 2013, 20(12): 3497-3503.

[13] KENNEDY J. Particle swarm optimization[M]//Encyclopedia of Machine Learning and Data Mining. Boston, MA: Springer, 2017: 967-972.

[14] KENNEDY J, EBERHART R. Particle swarm optimization[C]//Proceedings of ICNN'95-International Conference on Neural Networks. Perth, Australia: IEEE, 1995: 1942-1948.

[15] TAN Ying, ZHU Yuanchun. Fireworks algorithm for optimization[C]//Advances in Swarm Inteligence. Beijing, China: Springer, 2010: 355-364. DOI: 10.1007/978-3-642-13495-1_44.

[16] 谭营, 郑少秋. 烟花算法研究进展[J]. 智能系统学报, 2014, 9(5): 515-528.

TAN Ying, ZHENG Shaoqiu. Recent advances in fireworks algorithm[J]. CAAI Transactions on Intelligent Systems, 2014, 9(5): 515-528.

[17] 付国江, 王少梅, 刘舒燕, 等. 含维变异算子的粒子群算法[J]. 武汉大学学报(工学版), 2005, 38(4): 79-83.

FU Guojiang, WANG Shaomei, LIU Shuyan, et al. A PSO with dimension mutation operator[J]. Journal of Wuhan University of Hydraulic and Electric Engineering, 2005, 38(4): 79-83.

[18] LI Junzhi, ZHENG Shaoqiu, TAN Ying. Adaptive fireworks algorithm[C]//2014 IEEE Congress on Evolutionary Computation (CEC). Beijing, China: IEEE, 2014: 3214-3221.

[19] ZHENG Shaoqiu, JANECEK A, TAN Ying. Enhanced fireworks algorithm[C]//2013 IEEE Congress on Evolutionary Computation. Cancun, Mexico: IEEE, 2013: 2069-2077.

（编辑 陈灿华）

收稿日期： 2019 -03 -05; 修回日期： 2019 -05 -10

基金项目(Foundation item)：国家自然科学基金资助项目(51305467)(Project(51305467) supported by the National Natural Science Foundation of China)

通信作者：韩锟，博士，副教授，从事载运工具智能测控技术及性能优化研究；E-mail: hkun@csu.edu.cn

摘要:针对机车二系弹簧载荷优化调整这一复杂的多变量优化问题，为进一步提高现有求解方法的优化效果和计算效率，将烟花算法融入粒子群算法，提出一种具有分层递阶结构的改进粒子群算法，算法为3层架构，其中，底层是基础层，为加入维变异算子的粒子群算法，是改进算法的基本框架；中间层是融合层，为引入烟花算法爆炸机制的粒子更新层，主要用于扩大算法搜索范围，提高全局搜索能力；顶层是扰动层，引入扰动因子，避免算法因陷入局部搜索而进行的大量无为冗余迭代，加快全局收敛速度。用典型测试函数对改进算法性能进行测试，并将其应用于机车二系弹簧载荷分配优化调整仿真实验。研究结果表明：改进算法与传统遗传算法、烟花算法和粒子群算法相比，全局搜索能力更强，鲁棒性更好，求解精度更高。