脱轨美国一运有化学品列车脱轨起火

第17卷第11期铁道科学与工程学报Volume17Number112020年11月JournalofRailwayScienceandEngineeringNovember2020DOI:10.
19713/j.
cnki.
431423/u.
T20191174无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算龚凯1,向俊2,刘林芽1,周钦悦1,杨海明2,吴宝游3(1.
华东交通大学铁路环境振动与噪声教育部工程研究中心,江西南昌330013;2.
中南大学土木工程学院,湖南长沙410075;3.
中铁电气化局集团有限公司,北京100036)摘要:为研究无缝线路钢轨鼓胀引起的脱轨规律,以"单波形"钢轨鼓胀为例,建立无缝线路钢轨鼓胀状态下货物列车-轨道系统空间振动计算模型;按照列车脱轨能量随机分析方法,提出无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算方法,计算脱轨全过程,分析鼓胀状态对轮轨几何接触状态和相对位置的影响.
研究结果表明:鼓胀波长一定时脱轨系数和轮重减载率随鼓胀幅值增大而增大,鼓胀幅值一定时脱轨系数和轮重减载率随鼓胀波长减小而增大,但脱轨系数和轮重减载率对脱轨缺乏控制作用;转向架与钢轨横向相对位移随鼓胀波长的减小或幅值的增大而增大,得到考虑工况下的最大值为72.
6mm.
上述结果可为研发具有实时判别脱轨功能的检测装置提供基础数据.
关键词:铁道工程;无缝线路钢轨鼓胀;脱轨全过程;车轮脱轨检测装置中图分类号:U213.
2文献标志码:A文章编号:16727029(2020)11276707CalculationoffreighttrainderailmentprocesscausedbyrailbuckingofCWRGONGKai1,XIANGJun2,LIULinya1,ZHOUQinyue1,YANGHaiming2,WUBaoyou3(1.
EngineeringResearchCenterofRailwayEnvironmentalVibrationandNoise,MinistryofEducation,EastChinaJiaotongUniversity,Nanchang330013,China;2.
SchoolofCivilEngineering,CentralSouthUniversity,Changsha410075,China;3.
ChinaRailwayElectrificationBureauGroupCo.
,Ltd.
,Beijing100036,China)Abstract:InordertostudyfreighttrainderailmentlawcausedbytherailbuckingofCWR,thespatialvibrationcalculationmodelforthefreighttrain-tracksystemundertherailbuckingofCWRwasestablishedbytakingthe"singlewaveform"railbuckingasanexample.
Accordingtotherandomenergyanalysismethodoftrainderailment,acalculationmethodforthewholederailmentprocessoffreighttraincausedbytherailbuckingofCWRwasproposed.
Thetrainderailmentprocesseswerecalculatedandtheeffectsoftherailbuckingconditiononthewheel-railgeometriccontactstateandrelativepositionwereanalyzed.
Resultsshowthatthederailmentcoefficientandwheelloadreductionratioincreasewiththeincreaseofbuckingamplitudewhenthewavelengthofrailbuckingisfixed.
Meanwhile,thederailmentcoefficientandwheelloadreductionratioincreasewiththedecreaseofbuckingwavelengthwhenthebuckingamplitudeisfixed.
However,thederailmentcoefficientand收稿日期:20191226基金项目:国家自然科学基金委员会与神华集团有限公司联合资助项目(U1261113);国家自然科学基金资助项目(51578238,52068028);江西省优秀科技创新团队建设计划项目(20152BCB24007);华东交通大学科研启动基金资助项目(10012003416034);赣教高字[2018]江西省教育厅科学技术项目(GJJ170392)通信作者:向俊(1968),男,湖南溆浦人,教授,博士,从事列车脱轨控制、列车轨道(桥梁)系统振动及铁路轨道结构等研究;Email:jxiang@csu.
edu.
cn铁道科学与工程学报2020年11月2768wheelloadreductionratiolackcontrolfunctionforthederailmentcausedbytherailbucking.
Thelateralrelativedisplacementbetweenbogieandrailincreasewiththedecreaseofbuckingwavelengthortheincreaseofbuckingamplitudeoftherail,andthemaximumrelativedisplacementofthestudycaseis72.
6mm.
Theaboveresultscanprovidethebasicdataforthedevelopmentofadetectiondevicewithareal-timederailmentdetectionfunction.
Keywords:railwayengineering;railbuckingofCWR;wholeprocessofderailment;wheelderailmentdetectiondevice无缝线路因钢轨连接接头少、钢轨平顺性好,轨道结构维修费用低等特点,在高速、重载线路上大量铺设.
众所周知,轨道结构常年暴露在自然环境中,受温度影响显著,尤其是在烈日高照下,钢轨内部产生的温度应力使得在钢轨薄弱部分出现钢轨鼓胀现象,严重时甚至发生胀轨跑道.
国内外因钢轨鼓胀引起的列车脱轨事故仍时有报道,如2002年美国由弗罗里达州开出的一趟列车因钢轨鼓胀而发生脱轨[1]、2016年由台东开往彰化的651次"莒光号"列车因高温引起的钢轨变形而发生脱轨[2].
在列车运行安全事故中,列车脱轨给人们带来了严重的损失,为预防钢轨鼓胀引起的列车脱轨,国内外学者针对这一问题开展了一些研究[312].
Masayuki[3]调查了轨道横向不平顺对两轴列车行车安全的影响,提出对轨道横向不平顺应加强检测与修复;Cheli等[4]采用数值仿真及结合试验研究了钢轨缺陷对城轨电车行车安全的影响;Steenbergen等[5]探讨了钢轨焊接接头不平顺对列车行车中的轮轨相互作用力的影响;Ishida等[6]研究了钢轨垂向周期性不平顺对车辆脱轨安全性的影响,提出了一种新的脱轨评价准则;Mohammadzadeh等[79]系统研究了轨道不平顺对车辆动态响应的影响;XIAO等[1011]研究了钢轨鼓胀对高速列车行车安全的影响,并基于车轮抬升量、脱轨系数和轮重减载率等安全性指标建立了脱轨安全域准则.
实际上,要预防无缝线路钢轨鼓胀引起的列车脱轨,最理想的预防措施是保证无缝线路的稳定性不使其发生钢轨鼓胀,可是因复杂环境、人为或管理不善等因素使得因钢轨鼓胀引起的列车脱轨事故仍偶有发生,难以完全杜绝.
并且,在实际行车中,由于货车车辆数较多,在钢轨鼓胀引起的部分车辆脱轨后,驾驶员未能及时察觉,导致事故扩大,造成无法估量的损失.
因此,为最大程度减小无缝线路钢轨鼓胀引起的列车脱轨事故带来的损失,确保列车在脱轨第一时间判别车轮是否脱轨,并采取措施及时停车十分必要.
为此,基于列车轨道系统空间振动计算模型[12]、列车脱轨能量随机分析方法[13],以"单波形"钢轨鼓胀为例,建立无缝线路钢轨鼓胀状态下货物列车轨道系统空间振动计算模型,提出无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算方法,计算列车脱轨全过程,分析无缝线路钢轨鼓胀对车轮与钢轨间接触状态的影响,得到车轮脱轨瞬间轮轨相对位置和几何尺寸关系,为研发具有实时判别脱轨功能的检测装置提供基础数据.
1无缝线路钢轨鼓胀状态下货物列车轨道系统空间振动计算模型实践表明,无缝线路钢轨鼓胀后的线形多为"单波形"[14].
钢轨鼓胀前、后的无缝线路模型如图1所示.
图1钢轨鼓胀前、后的无缝线路模型Fig.
1ModelofCWRbeforeandaftertherailbucking在图1中,黑色虚线为钢轨鼓胀前的线形,红色实线为钢轨鼓胀部分线形,钢轨鼓胀路段设置在轨道长度的中间位置.
x为钢轨在X方向的位置;Hb为钢轨在Y方向的鼓胀幅值;Lb为钢轨鼓胀波长;L0为钢轨鼓胀轨道计算起点到钢轨鼓胀起点间的距离.
这里,以正弦函数为例模拟钢轨鼓胀线形,第11期龚凯,等:无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算2769如式(1)所示.
00b00bb0b0π()()sin0xLxLyxHLxLLLxLL≤≤(1)式中:y(x)为钢轨鼓胀的横向位移.
根据式(1)可计算轨道上任意位置处钢轨鼓胀引起的钢轨横向位移,将该位移与轨道初始不平顺叠加来模拟无缝线路钢轨鼓胀后的轨道线路.
基于上述轨道线路,设t时刻,一列由1辆机车+M辆货车组成的货物列车运行在长度为L的轨道上.
假定将长度为L轨道离散为N个轨段单元,针对每个单元建立具有34个自由度的空间振动模型,在模型中,将钢轨和轨枕分别视为弹性点支承Euler梁和放置于道床上的弹性变形体,扣件及道床顶面的弹性均采用线性弹簧与黏滞阻尼器模拟,并且不考虑道床振动.
进而可导出第j个轨段单元的空间振动势能ΠTj,并将每一个单元势能进行叠加得到长度L范围内轨道空间振动总势能ΠT[12],如式(2)所示.
TT1NjjΠΠ(2)同时,针对车体、前转向架和后转向架均考虑3个平动和3个转动自由度,每个轮对考虑横摆、沉浮2个自由度,形成具有26个自由度机车或车辆模型,进而导出第i辆机车或车辆的空间振动势能ΠVi,并将每一辆机车或车辆的势能叠加得到货物列车空间振动总势能ΠV[12],如式(3)所示.
VV1MiiΠΠ(3)由式(2)和式(3)可得无缝线路钢轨鼓胀状态下货物列车轨道系统空间振动总势能Π,如式(4)所示:VTΠΠΠ(4)为实现列车脱轨全过程计算,本文取消传统的轮轨"密贴"假定,采用轮轨位移衔接条件[13]作为列车与轨道系统间的纽带并满足轮轨"游间"要求,轮轨位移衔接条件如式(5)~式(6)所示:Yw=Yt+Yior+ΔYwt(5)Zw=Zt+Zior+ΔZwt(6)式中:Yw,Zw分别为车轮横向和竖向位移;Yt,Zt分别为钢轨横向和竖向位移;Yior,Zior分别为钢轨横向、竖向不平顺;ΔYwt,ΔZwt分别为车轮与钢轨间的横向和竖向相对位移.
同时,采用弹性系统动力学总势能不变值原理[15]和"对号入座"法则[13],形成无缝线路钢轨鼓胀状态下系统的矩阵方程如式(7)所示:KCMP(7)式中:[K],[C],[M]和{P}分别为刚度、阻尼、质量矩阵和荷载列阵;{},{}和{}分别为系统自由度的位移、速度和加速度列阵.
针对式(7)采用Wilson-θ逐步积分法求解,并编制计算程序.
为模拟车轮轮缘逐渐爬上钢轨顶部中点的全过程,将钢轨鼓胀路段设置在计算线路的中间位置.
2无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算方法研究表明[11,1617],列车低速脱轨表现为爬轨脱轨,而高速脱轨表现为跳轨脱轨.
本文研究重点是爬轨脱轨,仍基于文献[13]提出的脱轨机理开展无缝线路钢轨鼓胀引起的货物列车脱轨全过程研究.
列车脱轨机理表明[13],随着输入列车轨道系统能量不断增大,系统的横向振动逐渐增大,车轮轮缘逐渐爬上钢轨顶部,横向振动状态失去稳定,即列车脱轨.
文献[13]表明,采用构架实测(或人工)蛇行波及其标准差σp分别作为列车轨道系统横向振动激振源和输入能量.
可是,这里的标准差σp是正常行车时测得的,因列车脱轨试验不易实现,脱轨状态下构架蛇行波难以测得,相应的标准差也无法统计.
但由压杆失稳临界荷载计算方法可知,通过增大荷载,可试算得到压杆失稳临界状态时临界荷载,这一荷载是使压杆发生失稳的输入能量,比压杆处于失稳临界状态前的任一输入能量都大.
同理,使列车发生脱轨的构架蛇行波标准差σp1大于列车未脱轨时的构架蛇行波标准差σp,比σp1更大的系统输入能量不存在,因为列车已经脱轨.
为此,提出无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算方法如下:铁道科学与工程学报2020年11月27701)基于第1节中无缝线路钢轨鼓胀状态下货物列车轨道系统空间振动计算模型,将某一钢轨鼓胀时"单波形"线形的波长和幅值输入模型.
2)以实测、统计得到的标准差σp为基础,假定一个比σp更大的标准差σp1,采用Monte-Carlo法随机模拟出一条构架蛇行波输入步骤(1)中的模型,每计算一个时刻的系统振动响应,就按照车轮脱轨几何准则[13]判断该时刻车轮悬浮量是否达到25mm,若达到,则判定列车脱轨,计算停止;若未达到,则继续计算,直到满足上述要求为止.
此时,得到的系统空间振动响应反映了无缝线路钢轨鼓胀引起的货物列车脱轨全过程.
3)若步骤2)仍未能判定列车脱轨,则通过减小"单波形"线形波长或增大幅值,重复步骤1)~2),直到实现无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算为止.
3计算实例与分析研究及实践表明[13],我国自提速以来空载货车发生脱轨概率较大.
本节以1辆DF4型机车+16辆C62型空载敞车为例,计算车速V=80km/h时列车在300m有砟轨道上的脱轨全过程.
有砟轨道结构由60kg/m钢轨、II型混凝土轨枕及碎石道砟组成,轨枕间距为0.
5435m.
根据前述模型及计算方法,参考文献[11]中采用的钢轨鼓胀波长和幅值,分别计算钢轨鼓胀波长分别为20,30,40和50m,幅值分别为40,60,80和100mm时货物列车脱轨全过程,得到脱轨车轮在脱轨瞬间的悬浮量Δz,脱轨车辆Mi,脱轨车轮在脱轨瞬间的脱轨系数P及轮重减载率D,脱轨车轮对应的转向架与钢轨横向相对位移ΔH等动力响应,计算结果见表1所示.
由表1可知,因无缝线路钢轨鼓胀引起的脱轨车辆Mi为第1车(机车)或第2车,均位于列车编组前部;并且,在钢轨鼓胀波长Lb一定时,随着鼓胀幅值Hb的增大,钢轨不平顺逐渐急剧,脱轨车辆Mi由第2车变为第1车(机车).
可见,钢轨鼓胀对列车行车安全具有直接影响.
在表1的计算工况中,当波长Lb一定时,随着幅值Hb的增大,车轮脱轨瞬间脱轨系数增大;当幅值Hb一定时,随着波长Lb的增大,脱轨系数减小.
脱轨系数P最大为1.
14,最小为0.
18,均未超过规范GB5599—85[18]中要求的脱轨系数限值1.
2,但车轮发生脱轨.
因此,钢轨鼓胀波长及幅值对脱轨系数具有一定的影响,但脱轨系数难以判别车轮是否脱轨.
表1不同鼓胀波长及幅值条件下货物列车脱轨全过程计算结果Table1CalculationresultsofthewholederailmentprocessoffreighttrainatdifferentrailbuckingwavelengthandamplitudeLb/mHb/mmΔz/mmMiPDΔH/mm2040252车0.
510.
4264.
560252车1.
140.
7572.
680251车0.
410.
2954.
3100251车0.
480.
3659.
13040252车0.
420.
3860.
560252车0.
530.
6268.
280251车0.
350.
2452.
8100251车0.
390.
2858.
74040252车0.
410.
3556.
560252车0.
470.
3958.
880252车0.
520.
4363.
6100251车0.
220.
1848.
55040252车0.
350.
2852.
660252车0.
420.
3654.
980252车0.
470.
4157.
2100251车0.
180.
1546.
4此外,当鼓胀波长Lb一定时,鼓胀幅值Hb的增大,车轮在脱轨瞬间的轮重减载率增大;当鼓胀幅值Hb一定时,鼓胀波长Lb增大,相应的轮重减载率减小.
在表1各工况中,当车轮脱轨时,轮重减载率D最大为0.
75,超过规范GB5599—85[18]要求的限值0.
65,而D最小为0.
15,未超过限值要求.
可见,钢轨鼓胀波长及幅值对轮重减载率也具有一定的影响,但轮重减载率难以在其未超限时判别车轮是否脱轨.
为给研发具有实时判别脱轨功能的检测装置提供基础数据,本文计算了无缝线路钢轨鼓胀状态下车轮脱轨瞬间转向架与钢轨横向相对位移ΔH,计算结果如表1所示.
由表1可知,当波长Lb一定第11期龚凯,等:无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算2771时,机车、车辆在脱轨瞬间脱轨车轮对应的转向架与钢轨横向相对位移ΔH随幅值Hb的增大而增大;而当幅值Hb一定时,机车、车辆在脱轨瞬间脱轨车轮对应的转向架与钢轨横向相对位移ΔH随波长Lb的减小而增大.
可见,在钢轨鼓胀后,随着幅值Hb的增大及波长Lb的减小,钢轨逐渐形成较为急剧的不平顺,此时,机车或车辆振动剧烈.
此外,因机车质量大于车辆,则机车横摆幅度较小,机车脱轨车轮对应转向架与钢轨横向相对位移ΔH也比车辆的相应值小.
综上所述,在本文考虑的工况中,当钢轨鼓胀波长Lb为20m、幅值Hb为60mm时,转向架与钢轨横向相对位移ΔH最大,为72.
6mm.
为直观分析无缝线路钢轨鼓胀引起的货物列车脱轨全过程,钢轨鼓胀波长Lb为20m、幅值Hb为60mm时脱轨车轮的车轮悬浮量、脱轨系数、轮重减载率和转向架与钢轨横向相对位移时程曲线如图2~5所示.
在该工况下,根据钢轨鼓胀波长及轨道计算长度可得钢轨鼓胀路段设置在距轨道起点140~160m处.
图2第2车第1轴右轮悬浮量时程曲线Fig.
2Timehistoryofrightwheelliftvalueofthefirstaxleonthesecondwagon由图2表明,在列车运行至第7s前,第2车第1轴右轮悬浮量较小,而在第7s后第2车第1轴右轮悬浮量急剧增大至25mm,达到车轮脱轨几何准则[13]判定车轮脱轨要求.
按照列车运行速度,可得列车在第7s时已运行至80/3.
6*7.
0=155.
6m处,此时,列车已进入钢轨鼓胀路段.
可见,钢轨鼓胀对车轮悬浮量影响明显,对行车安全具有直接影响.
此外,由图3~4表明,车轮脱轨前脱轨系数及轮重减载率最大值分别为3.
63和0.
79(图中的负号表示与坐标方向相反),均超过规范限值[18]要求,但车轮未脱轨.
当车轮脱轨时脱轨系数为1.
14,轮重减载率为0.
75,脱轨系数虽未超过限值,但接近限值1.
2,而轮重减载率超过限值要求.
可见,钢轨鼓胀对脱轨系数和轮重减载率均具有影响,但脱轨系数和轮重减载率在其超过规范限值时难以判别车轮是否脱轨.
最后,由图5表明在脱轨瞬间,与第2车第1轴右轮相应的转向架与钢轨横向相对位移增大明显.
图3第2车第1轴右轮脱轨系数时程曲线Fig.
3Timehistoryofderailmentcoefficientofthefirstaxleonthesecondwagon图4第2车第1轴右轮轮重减载率时程曲线Fig.
4Timehistoryofwheelloadreductionrateofthefirstaxleonthesecondwagon铁道科学与工程学报2020年11月2772图5第2车前转向架与右侧钢轨横向相对位移时程曲线Fig.
5Timehistoryoflateralrelativedisplacementbetweenfrontbogieandrightrailonthesecondwagon4结论1)钢轨鼓胀波长及幅值对脱轨系数、轮重减载率均有影响,但脱轨系数、轮重减载率在其超过规范限值时难以判别钢轨鼓胀状态下车轮是否脱轨.
2)转向架与钢轨横向相对位移随鼓胀幅值的增大及波长的减小而增大明显.
3)通过计算得到了波长20m,幅值60mm时转向架与钢轨横向相对位移为72.
6mm,此为本文分析结果中的最大值.
该值反映了轮轨在脱轨瞬间的相对位置,并且考虑到检测装置不易安装在车轮表面,故将其安装在转向架上,以转向架的动态行为作为阀值,既克服检测装置不易安装,也起到了第一时间判别车轮脱轨的作用,进而为研发具有实时判别脱轨功能的检测装置提供参考.
4)本文研究的是列车低速时钢轨鼓胀引起的车轮爬轨脱轨,对于列车高速时钢轨鼓胀引起的车轮跳轨脱轨研究将在后续文章中报道.
参考文献:[1]BrabieD.
Ontheinfluenceofrailvehicleparametersonthederailmentprocessanditsconsequences[D].
Stockholm:KTH,Sweden,2005.
[2]俞菁.
台湾发生火车脱轨事件疑因高温致铁轨变形[EB/OL].
20160606.
http://news.
youth.
cn/gn/201606/t20160606_8087735.
htm.
YUJing.
TrainderailmentinTaiwanissuspectedtocausedbyrailwaytrackdeformationduetohightemperatures[EB/OL].
20160606.
http://news.
youth.
cn/gn/201606/t20160606_8087735.
htm.
[3]MasayukiI.
Astudyoftrackmaintenanceforthederailmentduetotheinteractionbetweentrackandvehicle[J].
JapaneseNationalRailwayTechnicalResearchInstituteQuarterlyReports,1978,19(1):16.
[4]CheliF,CorradiR,DianaG,etal.
Effectoftrackgeometricaldefectsonrunningsafetyoftramcarvehicles[J].
VehicleSystemDynamics,2006,44(Suppl1):302312.
[5]SteenbergenM,EsveldC.
Relationbetweenthegeometryofrailweldsandthedynamicwheel-railresponse:numericalsimulationsformeasuredwelds[J].
ProceedingoftheInstitutionofMechanicalEngineers,PartF:JournalofRailandRapidTransit,2006,220(4):409423.
[6]IshidaH,MatsuoM,FujiokaT.
Safetyassessmentmethodofrailwayvehicleunderoscillatorywheelloadfluctuation[J].
JournalofEnvironmentandEngineering,2007,2(2):407418.
[7]MohammadzadehS,SangtarashhaM,MolatefiH.
Anovelmethodtoestimatederailmentprobabilityduetotrackgeometricirregularitiesusingreliabilitytechniquesandadvancedsimulationmethods[J].
ArchiveofAppliedMechanics,2011,81(11):16211637.
[8]SuarezB,FelezJ,AntonioLozanoJ,etal.
Influenceofthetrackqualityandthepropertiesofthewheel-railrollingcontactonvehicledynamics[J].
VehicleSystemDynamics,2013,51(2):301320.
[9]ChoiIUmJH,LeeJS,etal.
Theinfluenceoftrackirregularitiesontherunningbehaviorofhigh-speedtrains[J].
ProceedingoftheInstitutionofMechanicalEngineers,PartF:JournalofRailandRapidTransit,2013,227(1):94102.
[10]XIAOXinbiao,JINXuesong,WENZefeng.
Effectoftangenttrackbuckleonvehiclederailment[J].
MultibodySystemDynamics,2011,25(1):141.
[11]肖新标.
复杂环境状态下高速列车脱轨机理研究[D].
成都:西南交通大学,2013.
XIAOXinbiao.
Studyonhigh-speedtrainderailmentmechanisminsevereenvironment[D].
Chengdu:SouthwestJiaotongUniversity,2013.
[12]XIANGJun,LIDejian,ZENGQingyuan.
Simulationofspatiallycouplingdynamicresponseoftraintracktime第11期龚凯,等:无缝线路钢轨鼓胀引起的货物列车脱轨全过程计算2773variantsystem[J].
JournalofCentralSouthUniversityofTechnology,2003,10(3):226230.
[13]曾庆元,向俊,周志辉,等.
列车脱轨分析理论与应用[M].
长沙:中南大学出版社,2006:204276.
ZENGQingyuan,XIANGJun,ZHOUZhihui,etal.
Analysistheoryandapplicationoftrainderailment[M].
Changsha:CentralSouthUniversityPress,2006:204276.
[14]陈秀方.
轨道工程[M].
北京:中国建筑工业出版社,2017:244246.
CHENXiufang.
Trackengineering[M].
Beijing:ChinaArchitecture&BuildingPress,2017:244246.
[15]ZENGQingyuan,LOUPing,XIANGJun.
Theprincipleoftotalpotentialenergywithstationaryvalueinelasticsystemdynamicsanditsapplicationtotheanalysisofvibrationanddynamicstability[J].
JournalofHuangzhongUniversityofScience&Technology(UrbanScienceEdition),2002,19(1):714.
[16]翟婉明.
车辆轨道耦合动力学[M].
4版.
北京:科学出版社,2015:427434.
ZHAIWanming.
Vehicle-trackcoupleddynamics[M].
4thed.
Beijing:SciencePress,2015:427434.
[17]陈淮,杜晓伟,曾庆元.
车辆脱轨问题的探讨[J].
郑州工业大学学报,2010,21(18):21922195.
CHENHuai,DUXiaowei,ZENGQingyuan.
Discussionaboutderailmentproblemofvehicles[J].
JournalofZhengzhouUniversityofTechnology,2010,21(18):21922195.
[18]GB5599—85,铁道车辆动力学性能评定和试验鉴定规范[S].
GB5599—85,Railwayvehicle-sepcificationforevaluationthedynamicperformanceandaccreditationtest[S].
(编辑涂鹏)