Optimization of JM3 Wheel Profile Considering Equivalent Conicity Dispersion
-
摘要:
轮轨匹配对轨道车辆动力学性能有着重要影响,针对现行标准JM3车轮型面与国内不同类型钢轨廓形匹配时等效锥度差异过大,以及在大轨底坡条件下匹配打磨钢轨时因等效锥度过低而引发机车晃车的问题,本文以减小车轮型面在不同轨底坡条件下与CN60型和CN60N型2种钢轨匹配对应的等效锥度离散度为优化目标,采用圆弧、直线组合的车轮型面描述方法,应用NSGA-Ⅱ遗传算法优化滚动圆附近两段圆弧圆心横向位置参数,对JM3型面进行优化,并对优化前后的车轮型面进行轮轨接触特性和机车动力学性能仿真对比分析. 结果表明:优化型面在与上述钢轨匹配时,轮对3 mm横移量处的名义等效锥度均在0.1左右,显著降低原JM3踏面等效锥度离散度,提高了车轮型面对不同钢轨廓形和线路条件的适应性;同时,优化后的型面对应机车蛇行稳定性、横向平稳性、曲线通过性能和磨耗性能指标较原型面均得到提升,消除了特定线路机车的低频晃车现象.
Abstract:Wheel-rail profile compatibility has an important influence on the dynamic performance of rail vehicles. The present standard JM3 wheel profile has a large difference in equivalent conicity when it is matched with different types of rail profiles in China. It has the problem of locomotive swaying caused by too low equivalent conicity when the profile is matched with the grinding rail with large rail cant. To address these issues, the optimization objective of reducing the equivalent conicity dispersion of the wheel profile matched with CN60 and CN60N rail profiles under different rail cants was set. The wheel profile was described by a method combining arcs and straight lines. The JM3 wheel profile was optimized by utilizing NSGA-II genetic algorithm to improve the lateral position parameters of the two arc centers near the rolling circle. The wheel-rail contact characteristics and the locomotive dynamic simulation of the wheel profile before and after optimization were compared. The results show that when the optimized wheel profile is matched with the rails above, the nominal equivalent conicities at the 3 mm transverse displacement of the wheelset are all about 0.1, which reduces the equivalent conicity dispersion of the original JM3 wheel profile and improves the adaptability of the wheel profile to different rail profiles and line conditions. At the same time, the locomotive hunting stability, lateral stability, curving performance, and wear performance index of the optimized profile are all improved compared with the original wheel profile. In addition, the phenomenon of low-frequency swaying of locomotives on specific lines is eliminated.
-
轮轨关系对轨道车辆的动力学性能有着重要影响,如蛇行运动稳定性、运行平稳性与舒适性以及运行安全性[1]. 近年来,随着我国铁路运输的快速发展,轨道车辆的运营速度不断提高,轮轨动态作用加剧,对于轮轨匹配提出了更高的要求. 车轮型面优化[2]和钢轨廓形优化[3]是改善轮轨匹配状态、提高车辆动力学性能的重要措施. 由于既有铁路线路大规模更换钢轨廓形成本高,车辆车轮型面种类繁多,各种类型的车辆存在跨线运行的情况,所以车轮型面优化相较于钢轨廓形优化有着更好的经济性与可行性.
对于车轮型面的生成方法有离散点曲线样条拟合方法[4]、车轮型面圆弧和直线特征参数表述法[5]以及旋转压缩微调法(rotary-scaling fine-tuning,RSFT)等[6]. 国内外学者在这些方法的基础上,为解决各种车辆动力学问题开展了大量车轮型面优化研究:Shen等[7]在忽略轮对侧滚角的前提下,基于轮轨接触角曲线反演车轮型面,并应用于独立轮对的开发;Jahed等[8]以滚动圆半径差(RRD)为目标,通过离散点和3次样条曲线表征优化区域,对车轮型面进行优化,并通过动力学软件仿真验证新踏面的优势;Polach[9]通过给定目标等效锥度以生成新的车轮型面,提高轮轨共形接触,减少了踏面的集中磨耗;Cui等[10]提出一种基于轮轨法向间隙的车轮型面正向优化方法,通过该方法得到的新踏面改善了接触点分布,降低接触应力,减少轮轨磨耗与滚动接触疲劳;Santamaria等[11]根据目标等效锥度曲线优化车轮型面,保证了车辆在不同的轨道配置下具有良好的动力学性能;干锋等[12]提出一种基于轮径差为目标反求车轮型面的方法,并通过S1002CN和LMA型面进行了验证;成棣等[13]以圆弧参数为设计变量,通过Isight软件对车轮型面进行多目标优化,改善了CRH3型动车组低锥度晃车问题;祁亚运等[14]采用旋转压缩微调法,以减小车轮磨耗和提高舒适度为目标对S1002CN型面进行多目标优化,通过动力学软件仿真得到良好的验证.
JM3车轮型面是基于中国标准60 kg/m钢轨(CN60型)匹配设计开发的,该型面与打磨钢轨廓形(CN60N型)匹配时,直线段因轮轨接触点集中于轨顶导致轮轨接触等效锥度较低,容易引发低频晃车;曲线段则因轮缘和钢轨内侧接触,导致曲线通过性能变差. 加之钢轨打磨偏差、实际线路轨底坡离散性[15]以及轨距离散性等因素,部分高速机车容易出现因轮轨匹配不良而导致车辆动力学性能恶化等问题,特别是在线路部分区段,钢轨轨底坡过大导致轮轨接触等效锥度过低,使机车出现晃车现象,降低了乘坐舒适性[16]. 针对传统车轮型面优化过程,很少考虑钢轨轨底坡的变化对轮轨接触等效锥度的影响,本文采用圆弧、直线等几何元素相组合的车轮型面描述方法[17],优化滚动圆附近两段圆弧圆心横向位置参数,针对现有JM3型面[18]进行优化,以减小在不同轨底坡条件下与CN60型和CN60N型钢轨廓形匹配时轮轨接触等效锥度离散度,降低钢轨廓形和轨底坡对等效锥度的影响,提高机车横向动力学性能,并对优化后的车轮型面进行车辆动力学性能仿真对比验证.
1. 车轮型面描述方法
常见的车轮型面如图1所示,图中:x、y分别为横、纵坐标,gˆg表示直线,A、B、…、N为节点. 该车轮型面由多个线段与圆弧组合而成,其中EI段为优化区域,由4段圆弧组成,分别是⌢EG、⌢GH、⌢HIA和⌢IAI. 在凹凸性方面,4段圆弧要满足从左到右依次为凹—凹—凹—凸. 令4段圆弧圆心及其坐标分别为O5(xO5,yO5)、O6(xO6,yO6)、O7(xO7,yO7)、O8(xO8,yO8),对应半径为R5、R6、R7、R8.
依据以下几何约束条件使用圆弧参数描述车轮型面[19]:
1) 车轮的轮缘角为70°;
2) EI段中相邻的两段圆弧相切,即⌢EG与⌢GH内切,⌢GH与⌢HIA内切,⌢HIA与⌢IAI外切;
3) 点E为⌢EG与直线gˆg的切点,点I为⌢IAI与直线IM的切点. 其中,点D与点M固定不动,点E与点I随着EI段中圆弧参数的变化而变化;
4) ⌢HIA经过踏面基点O(0, 0)以保证轮缘高度不变.
由⌢HIA经过基点O(0, 0)可得
yO7 = √R27−x2O7. (1) 由⌢GH与⌢HIA内切可得
yO6 = yO7−√(R7−R6)2−(xO7−xO6)2. (2) 由于⌢EG同时与直线gˆg和⌢GH相切,联立方程组可得
{xO5=−b1+√b21−4a1c12a1,yO5=yO6−√(R5−R6)2−(xO5−xO6)2, (3) 式中:a1=1+k2gˆg,kgˆg= −tan 70°;δ1=bgˆg+R5√1+k2gˆg−yO6,bgˆg=yD−kgˆgxD,xD和yD为点D的横、纵坐标;b1=2(kgˆgδ1−xR6);c1=x2O6+δ21−(R5−R6)2.
由于⌢IAI同时与直线IM和⌢HIA相切,联立方程组可得
{xO8=−b2+√b22−4a2c22a2,yO8=yO7−√(R8−R7)2−(xO8−xO7)2, (4) 式中:a2=1+k2IM,kIM=−0.1;δ2=bIM−R8√1+k2IM−yO7,bIM=yM−kIMxM,xM和yM为点M的横、纵坐标;b2=2(kIMδ2−xO7);c2=x2O7+δ22−(R7+R8)2.
可以看出用于描述EI段4段圆弧的12个圆弧参数(圆心位置和半径)可以缩减到用6个参数表示,即R5、R6、R7、R8和xO6、xO7.
2. 车轮型面多目标优化数学模型
由上述内容可知优化区域EI段可以由其中的4段圆弧半径R5、R6、R7、R8和中间2段圆弧的圆心横坐标xO6、xO7进行表述,在原有JM3型面基础上进行优化,为减小优化参数的数量,参考文献[14]的研究结果,同时保证优化前后轮轨接触应力不会出现较大变化,将中间4段圆弧的半径保持不变,即R5取14 mm,R6取100 mm,R7取450 mm,R8取250 mm,xO6、xO7为车轮型面优化设计的变量. 参考标准TB/T 449—2016[18]中各类车轮型面的尺寸,给定设计变量上下限,见表1.
表 1 设计变量上下限Table 1. Upper and lower limits of design variablesmm 圆弧参数 下限 上限 xO6 −9.0 −0.5 xO7 10.0 30.0 针对CN60型和CN60N型钢轨廓形,在现有JM3型面基础上进行车轮型面优化,降低轮轨匹配等效锥度受到不同钢轨廓形以及轨底坡随机性的影响. 本文针对CN60型和CN60N型钢轨的1/40和1/20 2种特定轨底坡工况,以车轮型面在第i种钢轨工况条件下匹配等效锥度λi与目标等效锥度λobj差值的2范数为目标函数,如式(5)所示.
fi(x)=‖λi−λobj‖2,i=1,2,3,4. (5) 式中:x为设计变量构成的向量.
计算等效锥度的匹配条件如表2所示.
表 2 计算等效锥度的匹配条件Table 2. Matching conditions for calculating equivalent conicity序号 钢轨 轨底坡 λi 1 CN60 1/40 λ1 2 1/20 λ2 3 CN60N 1/40 λ3 4 1/20 λ4 为保证优化后的车轮型面在2~3 mm的等效锥度内过渡平缓,目标等效锥度分别取轮对横移量2 mm和3 mm处的等效锥度. 由于标准JM3车轮型面与轨底坡1/40条件下CN60型钢轨匹配的车辆动力学性能优异,所以参照该轮轨匹配条件下的等效锥度,给定目标等效锥度,λobj=(0.09, 0.10)(λi=(λ2i, λ3i),λ2i、λ3i分别为轮对横移2 mm和3 mm处对应的等效锥度),xO6、xO7构成向量x=(x1, x2)=(xO6, xO7),建立车轮型面多目标优化模型,如式(6)所示.
minxmink⩽xk⩽xmaxkF(f1(x),f2(x),f3(x),f4(x)),k=1,2, (6) 式中:xmaxk和xmink分别为设计变量xk的上、下限,F(·)为总目标函数,fi(·)为目标函数.
由于NSGA-Ⅱ算法在保留不同侧重点精英解集的同时可以加快算法收敛速度,提高计算效率[20],本文采用该算法对4个目标函数同时进行寻优,得到图2所示的Pareto优化前沿,其中,水平坐标轴分别表示函数f1(x)和f2(x)的值,垂向坐标轴表示函数f3(x)的值,颜色深浅表示函数f4(x)的值. 根据计算结果,在解集中折中选取图2所示位置为最优解,此时对应的设计变量x =(−0.51, 26.85).
由此可得优化前后的车轮型面外形对比,如图3所示. 优化后的车轮型面记作JM3-opt,相比原标准JM3车轮型面,轮缘厚度保持不变,轮缘根部喉根圆曲线变化相对平缓,在横向位置 −10~10 mm内,曲线变化趋势相对较陡,10~40 mm曲线变化趋势相对缓和.
3. 优化结果分析
3.1 轮轨接触几何分析
图4和图5分别为原标准JM3型面和优化后的车轮型面JM3-opt在不同匹配条件下对应的轮轨接触点分布. 通过对比:在轮对横移量yw范围 −9~9 mm内,匹配CN60钢轨,当轨底坡为1/40时,优化后的车轮型面轮轨接触点分布与原型面相比,出现向钢轨内侧偏移,总体接触区域优化前后相差较小;当轨底坡为1/20时,优化后的车轮型面轮轨接触点分布与原型面相比范围更广,分布更加均匀. 匹配CN60N钢轨,当轨底坡为1/40时,优化前后的轮轨接触点分布差别较小;当轨底坡为1/20时,优化后的车轮型面轮轨接触点分布与原型面相比范围有所扩大,避免了接触点过于集中在轨顶区域.
图6为车轮型面优化前、后的等效锥度曲线对比;图7为优化前、后车轮型面在不同匹配条件下对应的名义等效锥度值. 通过图6(a)和图6(b)对比可以看出,优化后的车轮型面与原型面相比等效锥度整体有所增大,同时对钢轨类型和轨底坡的敏感度有所降低. 从图7可以看出:优化后的车轮型面JM3-opt与原标准JM3车轮型面相比,在轨底坡为1/40条件下,与上述2种钢轨廓形匹配对应的名义等效锥度相差都很小;在轨底坡为1/20条件下,与CN60钢轨匹配,名义等效锥度由0.04增大到0.13;与CN60N钢轨匹配,名义等效锥度由0.05增大到0.08. 优化后的车轮型面在不同匹配条件下,其名义等效锥度保持在0.1左右,减小了等效锥度的离散程度,同时有效解决了原JM3车轮型面在轨底坡为1/20时容易出现等效锥度过低的问题.
3.2 稳定性分析
为检验优化车轮型面JM3-opt的动力学性能,采用动力学仿真分析的方法对JM3车轮型面优化前后的动力学性能进行对比分析. 使用多体动力学商用软件SIMPACK建立某型动力集中动车组高速动力车动力学模型,整车模型包括轮对、轴箱、构架、空心轴、电机和牵引杆,车体共计25个刚体,90个自由度. 采用两系悬挂,同时考虑轮轨非线性接触关系,一、二系横向止挡的非线性以及各个减振器的非线性特性,其动力学模型如图8所示.
稳定性利用减速法[21]进行求解非线性临界速度,动力车初始速度为
1000 km/h,给轮对以初始激扰使其失稳,对机车施加恒定外力使其进行匀减速运动,考察运行时轮对横移的收敛情况来判断动力车的稳定性. 表3给出车轮型面优化前后在不同匹配条件下对应的非线性临界速度对比. 可以看出:优化后的车轮型面JM3-opt与原标准JM3车轮型面相比,在轨底坡为1/40条件下,分别与上述2种钢轨廓形匹配对应的非线性临界速度都很接近;在轨底坡为1/20条件下,与CN60钢轨匹配,非线性临界速度由215 km/h提升到345 km/h,如图9(a)所示;与CN60N钢轨匹配,非线性临界速度由200 km/h提升到530 km/h,如图9(b)所示.表 3 不同匹配条件的非线性临界速度Table 3. Nonlinear critical speeds under different matching conditionskm/h 钢轨 轨底坡 JM3 JM3-opt CN60 1/40 395 390 1/20 215 345 CN60N 1/40 660 670 1/20 200 530 由上述分析可以看出:优化后的车轮型面JM3-opt在与轨底坡为1/40的钢轨匹配时,非线性临界速度较原JM3车轮型面,差别不大;在与轨底坡为1/20的钢轨匹配时,JM3-opt车轮型面相较于原JM3车轮型面,非线性临界速度得到大幅提升,其中与CN60钢轨匹配时,非线性临界速度提升了60.5%;与CN60N钢轨匹配时,非线性临界速度提升了115.0%. 结果与前文的轮轨接触分析相契合:JM3-opt车轮型面在匹配轨底坡为1/40的钢轨时,名义等效锥度与原JM3车轮型面相差较小,故非线性临界速度差别也较小. 而与轨底坡为1/20的钢轨匹配时,原JM3车轮型面的名义等效锥度过小,导致其对应非线性临界速度过小. 而JM3-opt车轮型面对比原JM3车轮型面,名义等效锥度有所增加,非线性临界速度得到提升. 同时也可以看出名义等效锥度过大(或过小)都会使非线性临界速度降低.
3.3 平稳性分析
车辆平稳性计算采用武广线实测轨道谱作为轨道激励,速度区间设置为80~200 km/h. 车轮型面优化前后对应的横向平稳性指标对比如图10所示. 在该速度区间内,车轮型面优化前后对应的横向平稳性指标均在2.75以内,满足机车运行平稳性指标为优的标准.
由图10(a)可以看出:当轨底坡为1/40,与CN60钢轨匹配时,JM3-opt车轮型面对应的横向平稳性指标较优于原JM3车轮型面,在80~200 km/h的速度区间内,JM3-opt车轮型面的横向平稳性指标最大下降为2.8%,差异很小;与CN60N钢轨匹配时,JM3-opt车轮型面与原JM3车轮型面相比,横向平稳性指标基本一致,差异很小. 由图10(b)可以看出:当轨底坡为1/20,与CN60钢轨匹配时,当速度小于120 km/h,优化前后的车轮型面对应的横向平稳性指标差异很小;当速度大于120 km/h,JM3-opt车轮型面对应的横向平稳性指标要优于原JM3车轮型面. 对比原JM3车轮型面,在80~200 km/h的速度区间内,JM3-opt车轮型面的横向平稳性指标最大下降为7.9%;与CN60N钢轨匹配时,当速度小于160 km/h,优化前后的车轮型面对应的横向平稳性指标差异很小,当速度大于160 km/h,JM3-opt车轮型面对应的横向平稳性指标要优于原JM3车轮型面. 对比原JM3车轮型面,在80~200 km/h的速度区间里,JM3-opt车轮型面的横向平稳性指标最大下降为8.4%.
通过分析图10中JM3-opt车轮型面匹配不同类型钢轨廓形和轨底坡对应的横向平稳性指标差异很小,可以说明优化后的车轮型面对不同钢轨和轨底坡的适应程度有所改善.
当钢轨轨底坡条件为1/20时,车轮型面优化前后在速度200 km/h时对应的车体横向振动加速度对比,如图11所示. 与CN60钢轨匹配时,JM3-opt车轮型面对应的车体横向振动加速度对比原JM3车轮型面,振动主频由1.32 Hz增大到1.60 Hz,同时振动加速度幅值下降13.3%;与CN60N钢轨匹配时,JM3-opt车轮型面对应的车体横向振动加速度对比原JM3车轮型面,振动主频由1.18 Hz增大到1.39 Hz,同时振动加速度幅值下降35.2%,有效缓解了机车在特定线路发生低频晃车的现象.
3.4 曲线通过性能分析
车轮型面优化前、后对应的曲线通过性能指标如图12所示:优化前、后车轮型面对应的曲线通过性能指标均小于安全限值; 对于小半径曲线,JM3-opt车轮型面对应的曲线通过性能指标略优于原JM3车轮型面,但相差较小;随着曲线半径增大,JM3-opt车轮型面对应的曲线通过性能指标对比原JM3车轮型面,改善效果增大,其中,轮轴横向力和脱轨系数改善最明显. 根据图12,在R1 600曲线工况中,优化后的车型面相较于原JM3车轮型面. 与CN60钢轨匹配,当轨底坡为1/40时,轮轴横向力下降1.1 kN,脱轨系数下降0.01,优化前后差异很小;当轨底坡为1/20时,轮轴横向力下降5.9 kN,降低30.9%,脱轨系数下降0.08,降低33.3%. 与CN60N钢轨匹配,当轨底坡为1/40时,轮轴横向力下降3.7 kN,降低18.3%,脱轨系数下降0.06,降低19.4%;当轨底坡为1/20时,轮轴横向力下降5.4 kN,降低22.1%,脱轨系数下降0.09,降低24.3%.
3.5 磨耗性能分析
磨耗数W是评判车轮磨耗的重要指标[22],与车轮的实际磨耗量呈接近正比关系,如式(7)所示.
W=Txγx+Tyγy, (7) 式中:Tx、Ty分别为纵向和横向蠕滑力,γx、γy分别为纵向和横向蠕滑率.
优化前后车轮型面在不同曲线工况下磨耗数对比如图13所示,在与CN60钢轨匹配时,车轮型面优化前后对应的磨耗数大致相当;在与CN60N钢轨匹配时,JM3-opt车轮型面对应的磨耗数小于原JM3车轮型面,表明优化后的车轮型面抗磨耗性能优于原型面,随着曲线半径增大,抗磨耗性能得到进一步提升,且对于轨底坡为1/20时提升更为显著,在R1600曲线工况下磨耗数下降31.9%. 如图14所示,在R1600曲线工况,在轨底坡为1/20条件下与CN60N钢轨匹配,JM3-opt车轮型面在曲线段的磨耗数明显低于原JM3车轮型面,利用仿真手段证明了JM3-opt车轮型面相比原JM3车轮型面有更好的抗磨耗性能.
4. 结 论
1) 本文基于NSGA-Ⅱ遗传算法,利用圆弧、直线组合描述车轮型面的方法,以中间滚动圆附近的2段圆弧圆心横向位置参数xO6、xO7为设计变量,针对CN60和CN60N型钢轨廓形,以及1/20和1/402种轨底坡,通过给定目标等效锥度,对JM3型面进行多目标优化,以减小车轮型面在不同钢轨廓形和轨底坡条件下的等效锥度离散度. 该方法通过引入2个优化变量,在大大提高计算效率的同时也能满足各段曲线连接的光滑性以及凹凸性.
2) JM3-opt车轮型面在不同轨底坡条件下,与CN60和CN60N钢轨匹配时的名义等效锥度均在0.10左右,较原JM3车轮型面减小了等效锥度的离散程度,改善了车轮型面对不同钢轨和轨底坡的适应性,同时有效解决了原JM3车轮型面在轨底坡为1/20时容易出现等效锥度过低的问题.
3) 建立某型动力集中动车组动力车模型,对车轮型面优化前后的动力学性能进行仿真对比分析. 当轨底坡为1/40时,采用优化型面JM3-opt后动力车稳定性、平稳性、曲线通过性能较JM3型面差异很小;当轨底坡为1/20时,优化型面JM3-opt对应车辆动力学性能较原型面得到显著提升;优化型面JM3-opt在与CN60钢轨匹配时抗磨耗性能较JM3型面差异很小;与CN60N钢轨匹配时抗磨耗性能较JM3型面提升明显,曲线半径越大,改善效果越明显.
致谢:牵引动力国家重点实验室自主研究课题(2022TPL-T02)资助.
-
表 1 设计变量上下限
Table 1. Upper and lower limits of design variables
mm 圆弧参数 下限 上限 xO6 −9.0 −0.5 xO7 10.0 30.0 表 2 计算等效锥度的匹配条件
Table 2. Matching conditions for calculating equivalent conicity
序号 钢轨 轨底坡 λi 1 CN60 1/40 λ1 2 1/20 λ2 3 CN60N 1/40 λ3 4 1/20 λ4 表 3 不同匹配条件的非线性临界速度
Table 3. Nonlinear critical speeds under different matching conditions
km/h 钢轨 轨底坡 JM3 JM3-opt CN60 1/40 395 390 1/20 215 345 CN60N 1/40 660 670 1/20 200 530 -
[1] PERSSON I, IWNICKI S. Optimisation of railway wheel profiles using a genetic algorithm[J]. Vehicle System Dynamics, 2004, 41: 517-526. [2] 李立,崔大宾,金学松. 车轮型面优化的研究进展[J]. 西南交通大学学报,2009,44(1): 13-19.LI Li, CUI Dabin, JIN Xuesong. State of arts of research on railway wheel profile optimization[J]. Journal of Southwest Jiaotong University, 2009, 44(1): 13-19. [3] 吴磊,康彦兵,董勇,等. 考虑打磨量的重载钢轨打磨廓形优化设计[J]. 西南交通大学学报,2022,57(4): 805-812.WU Lei, KANG Yanbing, DONG Yong, et al. Optimal design of heavy-haul rail grinding profile considering grinding amount[J]. Journal of Southwest Jiaotong University, 2022, 57(4): 805-812. [4] SHEVTSOV I Y, MARKINE V L, ESVELD C. Optimal design of wheel profile for railway vehicles[J]. Wear, 2005, 258(7/8): 1022-1030. [5] 成棣,王成国,刘金朝,等. 以圆弧参数为设计变量的车轮型面优化数学模型研究[J]. 中国铁道科学,2011,32(6): 107-113.CHENG Di, WANG Chengguo, LIU Jinzhao, et al. Research on the mathematical model for wheel profile optimization with arc parameters as design variables[J]. China Railway Science, 2011, 32(6): 107-113. [6] YE Y G, QI Y Y, SHI D C, et al. Rotary-scaling fine-tuning (RSFT) method for optimizing railway wheel profiles and its application to a locomotive[J]. Railway Engineering Science, 2020, 28(2): 160-183. doi: 10.1007/s40534-020-00212-z [7] SHEN G, AYASSE J B, CHOLLET H, et al. A unique design method for wheel profiles by considering the contact angle function[J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2003, 217(1): 25-30. doi: 10.1243/095440903762727320 [8] JAHED H, FARSHI B, ESHRAGHI M A, et al. A numerical optimization technique for design of wheel profiles[J]. Wear, 2008, 264(1/2): 1-10. [9] POLACH O. Wheel profile design for target conicity and wide tread wear spreading[J]. Wear, 2011, 271(1/2): 195-202. [10] CUI D B, LI L, JIN X S, et al. Optimal design of wheel profiles based on weighed wheel/rail gap[J]. Wear, 2011, 271(1/2): 218-226. [11] SANTAMARIA J, HERREROS J, VADILLO E G, et al. Design of an optimised wheel profile for rail vehicles operating on two-track gauges[J]. Vehicle System Dynamics, 2013, 51(1): 54-73. doi: 10.1080/00423114.2012.711478 [12] 干锋,戴焕云,池茂儒,等. 铁道车辆车轮踏面反向优化设计方法[J]. 铁道学报,2015,37(9): 17-24.GAN Feng, DAI Huanyun, CHI Maoru, et al. A reverse optimal design method for wheel tread of railway vehicle[J]. Journal of the China Railway Society, 2015, 37(9): 17-24. [13] 成棣,孙琛,胡晓依,等. 面向低锥度晃车的CRH3平台动车组车轮型面优化[J]. 中国铁道科学,2020,41(6): 135-144.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. [14] 祁亚运,戴焕云,干锋,等. 基于车轮磨耗和舒适度的CRH3型动车组型面优化研究[J]. 振动与冲击,2021,40(18): 148-155.QI Yayun, DAI Huanyun, GAN Feng, et al. Wheel profile optimization of CRH3 type of EMU based on wheel wear and passenger comfort[J]. Journal of Vibration and Shock, 2021, 40(18): 148-155. [15] 杜星,陶功权,杨城,等. 轨底坡变化对高速车辆运行行为的影响[J]. 西南交通大学学报,2022,57(2): 286-294.DU Xing, TAO Gongquan, YANG Cheng, et al. Influence of different rail cants on dynamical characteristics of high-speed railway vehicles[J]. Journal of Southwest Jiaotong University, 2022, 57(2): 286-294. [16] SUN J F, CHI M R, JIN X S, et al. Experimental and numerical study on carbody hunting of electric locomotive induced by low wheel–rail contact conicity[J]. Vehicle System Dynamics, 2021, 59(2): 203-223. doi: 10.1080/00423114.2019.1674344 [17] 宋志坤,孙琛,成棣,等. 车轮型面圆弧参数及其对轮轨接触和车辆动力学影响研究[J]. 中国铁道科学,2019,40(6):104-113.SONG Zhikun, SUN Chen, CHENG Di, et al. Research on arc parameters of wheel profile and its influence on wheel-rail contact and vehicle dynamics[J]. China Railway Science, 2019, 40(6):104-113. [18] 国家铁路局. 机车车辆车轮轮缘踏面外形:TB/T 449—2016[S]. 北京:中国铁道出版社,2016. [19] 孙琛. CRH3型平台高速动车组车轮型面多目标优化研究[D]. 北京:北京交通大学,2020. [20] 张维煜,李凯,杨鑫. 基于参数优先级划分的飞轮电机多目标优化[J]. 西南交通大学学报,2023,58(4): 922-932.ZHANG Weiyu, LI Kai, YANG Xin. Multi-objective optimization for flywheel motors based on parameter priority division[J]. Journal of Southwest Jiaotong University, 2023, 58(4): 922-932. [21] 石怀龙,罗仁,曾京. 国内外高速列车动力学评价标准综述[J]. 交通运输工程学报,2021,21(1): 36-58.SHI Huailong, 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. [22] BRACCIALI A, MEGNA G. Contact mechanics issues of a vehicle equipped with partially independently rotating wheelsets[J]. Wear, 2016, 366/367: 233-240. doi: 10.1016/j.wear.2016.03.037 -