高速动车组线性涡流制动系统特性仿真研究_杨利强.pdf

第 1 期学术专论高速动车组线性涡流制动系统特性仿真研究杨利强1徐凯歌2刘赛赛3(1徐州地铁集团有限公司，221018，徐州;2徐州地铁运营有限公司，221018，徐州;3南京中车浦镇海泰制动设备有限公司，211899，南京第一作者，正高级工程师)摘要基于涡流制动原理建立涡流制动力的数学模型，并利用 ANSYS Maxwell 软件建立 LECB(线性涡流制动)三维仿真模型。根据控制变量法研究列车速度、气隙、励磁电流等因素对涡流制动特性的影响，并分析了常用制动和紧急制动工况下的电磁特性。研究结果表明:线性涡流制动力受速度的影响明显，低速时制动力快速上升并达到幅值，然后随着速度的增加，制动力下降并趋于平稳;励磁电流、励磁线圈匝数与线性涡流制动力成正相关，气隙、钢轨材料电导率与线性涡流制动力成负相关;相同条件下，励磁线圈材料为铝时，线性涡流制动系统产生的制动力大小优于励磁线圈材料为铜时产生的制动力。关键词高速动车组;涡流制动;有限元分析;控制变量法中图分类号U270357U266DOI:1016037/j1007869x202301001Characteristic Simulation of Linear Eddy CurrentBraking System of High-speed EMUYANG Liqiang，XU Kaige，LIU SaisaiAbstractA mathematical model of eddy current braking forceis established based on eddy current braking principle，and theLECB(linear eddy current brake)3D simulation model is es-tablishedusingANSYSMaxwellsoftwareAccordingtocontrol variable method，the influence of factors such as trainspeed，air gap，excitation current on eddy current braking char-acteristics is studied，and the electromagnetic characteristicsunder common braking and emergency braking conditions areanalyzed esearch results show that the linear eddy currentbraking force is obviously affected by the speed At low speed，the braking force rises rapidly and reaches the amplitude，andthen decreases and becomes stable as speed increases The exci-tation current and the number of excitation coil turns are posi-tively correlated with the linear eddy current braking force，while the air gap and the conductivity of rail material are nega-tively correlated with the linear eddy current braking force Un-der the same conditions，the braking force generated by the lin-ear eddy current braking system with excitation coil material ofaluminum is larger than with material of copperKey wordshigh-speed EMU;eddy current braking;finite el-ement analysis;control variable methodFirst-authors addressXuzhou Metro Group Co，Ltd，221018，Xuzhou，China线性涡流制动是一种不依靠车轮与钢轨之间的摩擦进行制动的方式，其在高速下可提供稳定、可靠的制动力，已逐渐成为高速列车上一种新型的制动方式1-2。但线性涡流制动系统受列车速度、气隙等因素影响，想要获得安全可靠的制动力就要为线性涡流制动系统选取合适的参数范围。关于线性涡流制动的研究，国内外学者已进行过大量的试验与分析。法国研发的 AGV(高速动车组)列车运用的是线性涡流制动系统进行制动，当列车运行速度在200 km/h 以上时，每列列车的瞬时制动力可达到 20 kN3-4。德国从 1985 年开始在ICE-V(试验型城际列车特快)试验车上进行试验。结合 ICE-V 试验车的问题及法国 AGV 列车的经验，德国铁路公司在 ICE-3(德国高速动车组 ICE 的第三代)试验车上对改良后的线性涡流制动系统进行了试验并取得了良好的效果。其研究结果表明:ICE3 试验车使用线性涡流制动系统进行制动时，若动车速度升高，其制动力也能保持较稳定的水平，而制动功率则跟随试验车速度的变化而变化(近似的可以看成线性关系)，但线性涡流制动热负荷却比摩擦制动产生的热负荷低很多。这为列车线性涡流制动技术的推广提供了应用基础5-6。1964 年，日本建成的新干线是世界上第一条高速铁路，长度约为 550 km，最快速度可达 210 km/h。大约 6 年后，日本开始在新干线上进行一系列线性涡流制动试验。1990 年，日本又开始研究电磁涡流制动技术和永磁涡流制动技术。日本在其高速列车上运用电制动与空气制动组成的联合制动系统进行制动，列车的拖车采用旋转型涡流制动系统作为12023 年电制动系统。旋转型涡流制动在高速列车上的成功应用表明，在拖车上使用旋转型涡流制动与空气制动相结合的制动方式是行之有效的7。文献 8以“迎流的”有限元法为基础，推导了涡流制动力公式，并得到其制动力曲线。由该曲线可以看出，随着列车速度的增加，其制动力呈现出先变大后减小的趋势，并在中间某速度处达到最大值。文献 9对永磁涡流制动器进行了分析，求解出其制动力公式，并分析了不同参数变化对制动力的影响。文献 10 利用德国 T07 型磁浮列车的线性涡流制动装置，推导出其制动力的关系式，得到了制动力随列车速度变化的曲线。文献 11采用子域法计算了线性涡流制动力与速度、励磁电流、气隙长度、板厚和电导率等参数之间的解析公式，并通过有限元仿真验证了理论计算的正确性。由上述研究可知，线性涡流制动系统产生的制动力与所选取的不同参数有关，因此研究不同参数对线性涡流制动系统的影响至关重要。本文通过 ANSYS Maxwell 仿真软件对线性涡流制动系统进行建模仿真，分析了线性涡流制动系统在不同车速、气隙、励磁电流及励磁线圈匝数时产生的制动力大小，为线性涡流制动系统的设计提供借鉴与参考。1线性涡流制动系统的数学模型高速列车线性涡流制动系统是由电磁系统和机械系统组成的。电磁系统包括牵引变流器、励磁电磁铁等，其中励磁电磁铁由磁轭、电枢线圈和铁心组成。机械系统一般由升降风缸、套筒等组成。励磁电磁铁与钢轨平行，其包括 8 个磁极，线圈的电流通入方向使磁极按照“N-S-N-S”交替排列。当高速列车利用线性涡流制动系统进行制动时，钢轨切割电磁铁产生的磁场使钢轨中感应出涡流并形成涡流磁场，该新生磁场将阻碍原本电磁铁产生的磁场变化而使气隙中的磁场畸变产生 2 个力的分量，其中一个横向电磁力会阻碍列车运动，从而达到列车制动的目的。由于线性涡流制动力的计算涉及到很多复杂因素，因此，本文根据线性涡流制动系统的工作原理及其在不同工况下的工作过程，采用磁路分析法对涡流磁场进行相应简化，进而建立线性涡流制动力的数学模型。由图 1 的线性涡流制动磁场示意图可知，磁极与磁极之间构成闭合回路，由于每两个磁极之间形成的回路相同且对称，因此分析电磁铁的电磁机理时，可将多个磁极等效成一对磁极进行分析。图 1线性涡流制动磁场示意图Fig1Diagram of linear eddy current braking magnetic field根据等效原理，可将磁极与钢轨之间产生的涡流等效成一个半径为 r、集肤深度为 的涡流计算微元。涡流区域磁通是按照余弦规律变化的，因此在钢轨表面产生的感应电动势 E 为:E=ddt=B2r2vsinvt()(1)式中:线性涡流区域磁通;B 气隙处的磁感应强度;v 列车速度;相邻电磁铁之间的极距;t 变化时间。感应涡流在钢轨表面的集肤深度为:=2 v(2)式中:钢轨的电阻率，/m;钢轨磁导率，H/m。由式(1)和式(2)可得，涡流有效值 为:=2BvD232(3)式中:D 涡流区域直径，mm。按照等效磁路原理，对磁路模型各部分磁阻进行等效，得到气隙处的磁感应强度为:B=32 0NI032 l0+2ke0vD2(4)式中:0 空气磁导率，H/m;I0 励磁电流，A;ke 折算系数;N 线圈安匝数，匝;l0 气隙长度，mm。2第 1 期学术专论从能量守恒的角度出发，可以得到涡流区域内产生的功率有效值 P 为:P=B23v2D4256 2(5)根据线性涡流制动力 F=Pv，可得:F=25/203/21/2NI0()2D4v16 l0+ke0D2v()2(6)从式(6)中可以看出，线性涡流制动力与列车速度、气隙长度、励磁电流等诸多因素有关。此外，涡流密度衰减率 可以表示为:=exp 2z()(7)式中:z 深度，mm;涡流频率，Hz;钢轨电导率，S/m。由式(7)可知，和 对电流衰减起到了重要作用。当 不变时，或 越大，则在同一 z 处衰减得越大12，因此 和 的取值也会影响制动力的大小。综上可知，列车速度、气隙、励磁电流、励磁线圈匝数、钢轨材料电导率及励磁线圈材料等能够影响线性涡流制动特性和电磁特性。2线性涡流制动系统仿真模型21模型建立利用 ANSYS Maxwell 软件建立线性涡流制动系统仿真模型(以下简称“仿真模型”)，如图 2 所示。仿真模型尺寸如表 1 所示。图 2线性涡流制动系统仿真模型Fig2Simulation model of linear eddy current braking system由于钢轨是铁磁性物质，该仿真模型中取其相对磁导率为 200 H/m，电导率为 5106S/m;磁轭与磁心的相对磁导率取值为 1 360 H/m，电导率为0153 8 S/m;励磁线圈选铜材料。表 1线性涡流制动系统仿真模型尺寸Tab 1Simulation model size of linear eddy currentbraking system名称尺寸铁心115 mm75 mm110 mm励磁线圈170 mm130 mm102 mm磁轭1 460 mm130 mm40 mm磁极两端极板35 mm130 mm110 mm磁矩(中心距)182 mm为保证仿真模型计算的准确率，要对模型进行剖分;由于线形涡流制动系统会在钢轨表面产生涡流，所以针对钢轨要按照集肤深度进行剖分。线性涡流制动系统网格剖分图如图 3 所示。图 3线性涡流制动系统网格剖分Fig3Grid division of linear eddy current braking system22制动特性仿真常用制动时，线性涡流制动系统励磁电流 I078 A;紧急制动时 I0为95 A。气隙长度 l0变化范围为 68 mm，制动特性试验参数如表 2 所示。表 2制动特性试验参数Tab2Experimental parameters of braking characteristics仿真试验条件参数仿真试验 1:改变气隙长度l0为 6 mm、7 mm、8 mm;I0为 78 A;匝数为 280 匝仿真试验 2:改变励磁电流l0为 7 mm;I0为 60 A、78 A、95 A;匝数为280 匝仿真试验 3:变匝数l0为 7 mm;I0为 78 A;匝数为 280 匝、380匝、480 匝、580 匝仿真试验 4:变钢轨电导率l0为 7 mm;钢轨电导率为 2、4;I0为78A;匝数为 280 匝仿真试验 5:变励磁线圈材料(铜换铝)l0为 7 mm;I0为 78 A;匝数为 280 匝221励磁电流对制动力的影响设定不同励磁电流 I0研究其制动力与速度变化之间的关系，如图 4 所示。由图 4 可知，不同励磁电流下的曲线变化趋势基本一致。当气隙长度保持不变时，增加励磁电流值，制动力数值也随之增大，说明励磁电流与制动力成正相关，这也符合制动力数学模型。在励磁电流分别为 60 A、78 A、95A 时，其相应的最大制动力分别约为 10538 0 kN、15100 9 kN、19887 8 kN。由此可见，当励磁电流约增加 30%时，其制动力增加了 44%。32023 年图 4不同励磁电流下制动力随速度变化曲线Fig 4Curve of braking force changing with speed fordifferent excitation currents222气隙长度对制动力的影响保持励磁电流不变，通过改变气隙长度 l0来观察制动力的变化，如图 5 所示。图 5不同气隙长度下制动力随速度变化曲线Fig 5Curve of braking force changing with speed fordifferent air gap lengths由图 5 可知，随着气隙长度的增加，其制动力曲线峰值不断下降。气隙长度越小，制动力曲线的变化趋势越大，且制动力呈明显增大的趋势。这是因为磁场强度与磁路中磁阻大小有关，气隙长度增大造成了磁阻增大、磁场强度减小，使得涡流密度降低，因此制动力有所下降。当气隙长度分别为 6mm、7 mm、8 mm 时，其相应的最大制动力分别为16412 6 kN、15100 9 kN、13964 5 kN。由此可见，当气隙长度减小约 33%时，其制动力减小了 172%。此外，综合图4 和图5 可知，涡流制动力受速度影响较为明显，速度越大，制动力曲线上升速度越快，当速度为10 km/h 左右时就可以达到涡流制动力的最大值，此时动车速度称为临界速度。当速度持续增大到高速区时，涡流制动力呈下降趋势，最后趋于某个稳定值，总体来说在高速状态时呈良好的制动状态。223励磁线圈匝数对制动力的影响为研究励磁线圈匝数变化对制动力的影响，选取列车速度在 50 km/h、250 km/h、500 km/h 时进行分析，研究不同匝数对制动力变化的影响。在相同的条件下，取励磁线圈为 280 匝、380 匝、480 匝、580 匝分别进行计算，得到 3 个速度下不同励磁线圈匝数的制动力变化曲线，如图 6 所示。图 6不同列车速度下励磁线圈匝数对制动力变化的影响曲线Fig6Curve of the number of excitation coil turns influencingbraking force change at different train speeds由图 6 可知，当速度为 50 km/h 时，制动力增加速度最快;当速度为 500 km/h 时，制动力增加速度最缓慢。由此可见，当励磁线圈匝数增加时，低速运行的列车制动力增加更快，符合制动力的变化规律。在列车运行速度为 50 km/h 下，励磁线圈匝数为 580 匝时的制动力为 33500 0 kN;励磁线圈匝数为 280 匝时的制动力为 10700 0 kN。由此可见，励磁线圈匝数增大 107%，其制动力增加了 214%，可见制动力增加的倍数远大于匝数增加的倍数，即增加线圈匝数可以有效提高制动力。224钢轨电导率对制动力的影响保持励磁电流为 78 A、气隙长度为 7 mm、励磁线圈匝数为 280 匝不变，将钢轨电导率 分别变为原来的 2 倍、4 倍，并对模型进行计算，得到不同磁导率下制动力随速度变化的曲线，如图 7 所示。图 7不同钢轨电导率下制动力随速度变化曲线Fig 7Curve of braking force changing with speed fordifferent rail electrical conductivity由图 7 可知，当钢轨的电导率越大，涡流制动力反而越小。由此可知，钢轨的材料特性对涡流制动力有较大的影响。选择钢轨材料时，在符合工程条4第 1 期学术专论件的前提下，应尽可能选择电导率较小的材料，才能有更好的制动性能。225励磁线圈材料对制动力的影响在保持励磁电流为 78 A、气隙长度为 7 mm、励磁线圈匝数为 280 匝时，将励磁线圈材料由铜换成铝，将相同结构的铜线圈与铝线圈的涡流制动系统产生的制动力进行对比，得到不同速度下的制动力曲线，如图 8 所示。图 8采用铜线圈和铝线圈时的制动力随速度的变化曲线Fig8Curve of braking force changing with speed whenadopting copper coil and aluminum coil从图 8 中可以看出，在匝数相同的情况下，不同速度下铜线圈产生的涡流制动力与铝线圈的涡流制动力变化趋势一致，但铝线圈的涡流制动力比铜线圈的涡流制动力略大，说明在线圈匝数、气隙长度、线圈励磁电流相同的条件下，选用铝线圈获得的涡流制动力效果更好。23电磁特性仿真本节主要分析列车速度为400 km/h 时，常用制动和紧急制动下的磁感应强度分布情况。通过上述分析，线性涡流制动系统参数选择励磁电流 I0为78 A、气隙长度 l0为 7 mm、励磁线圈材料为铝线圈。其分析结果如图 9图 12 所示。图 9励磁电流为 78 A 时涡流制动系统磁场分布云图(常用制动工况)Fig9Nephogram of eddy current braking system magneticfield distribution with excitation current of 78 A(nor-mal braking condition)图 10励磁电流为 78 A 时气隙磁场分布(常用制动工况)Fig10Air gap magnetic field distribution with excitation cur-rent of 78 A(normal braking condition)图 11励磁电流为 95 A 时涡流制动系统磁场分布云图(紧急制动工况)Fig11Nephogram of eddy current braking system magneticfield distribution with excitation current of 95 A(ur-gent braking condition)图 12励磁电流为 95 A 时气隙磁场分布(紧急制动工况)Fig12Air gap magnetic field distribution with excitation cur-rent of 95 A(urgent braking condition)由图 9 可知，在磁极附近磁场强度最强，随着与磁极距离越远，磁场强度呈衰减趋势。由于相邻磁极之间形成闭合回路，在单个磁极上方磁轭处两侧的磁感应强度大于中间的磁感应强度。此外，在相邻磁场的叠加作用下，钢轨表面产生了较强的磁感应强度。由图 11 可知，紧急制动工况下，磁感应强度的分布情况与常用制动工况下基本一致。同时，由于紧急制动工况下励磁电流增大，使得涡流制动系统的整体磁感应强度有所增强。52023 年由图 10 和图 12 可知，两种制动工况下的磁感应强度的变化趋势基本一致，但紧急制动工况的磁感应强度值要大于常用制动工况。气隙长度处的磁场强度在相邻磁极之间磁场较强，在涡流制动装置两侧存在磁场泄露的问题，在短距离下会迅速衰减为 0。在常用制动情况下，气隙长度处的磁场强度最高可达到 877640 kA/m，紧急制动下可以达到1 030595 kA/m。通过上述分析说明，紧急制动工况下，涡流制动系统可以产生更强的磁感应强度和更高的制动力，以保证列车的运行安全。3结论1)本仿真研究得到的制动力曲线变化趋势符合涡流制动力数学模型中各参数对制动力的影响规律，从而验证了仿真模型的正确性。2)涡流制动力随着速度的增加先增大后减小，在某个中间速度值达到临界点。气隙长度增加，涡流制动力增大;励磁电流增加，涡流制动力增大。3)当励磁线圈匝数增大后，涡流制动力增加的倍数大于励磁线圈匝数增加的倍数，说明增加线圈匝数可以有效提高涡流制动力。4)当励磁线圈材料为铝时，与铜线圈相比，其制动力有所增加。在符合工程的条件下，可以考虑将励磁线圈材料选为铝线圈。5)当增加钢轨材料电导率时，发现涡流制动力反而下降，说明钢轨材料电导率与涡流制动力成负相关。因此在符合工程的条件下，选择电导率较小的钢轨材料可以获得更好的涡流制动效果。6)紧急制动工况下，涡流制动系统可以提供更强的磁场、产生更大的涡流制动力，以保障列车的运行安全。参考文献1张建柏，彭辉水，倪大成，等高速列车制动技术综述J机车电传动，2011(4):1ZHANG Jianbai，PENG Huishui，NI Dacheng，et al Overvie-wing braking technology of the high-speed trainsJ ElectricDrive for Locomotives，2011(4):12李培署，王风洲我国高速动车组制动技术现状及未来技术发展探讨J 铁道车辆，2018，56(9):1LI Peishu，WANG Fengzhou Present conditions of brakingtechnology for high speed multiple units in China and discussionof technical development in futureJ olling Stock，2018，56(9):13应之丁涡流制动技术在高速列车上的应用J电力机车与城轨车辆，2004，27(5):19YING ZhidingApplication of eddy current brake technology onhighspeed trainsJElectric Locomotives Mass Transit Ve-hicles，2004，27(5):194TEIB Ludwig，谢小海线性涡流制动装置的运用经验J国外机车车辆工艺，2013(3):38TEIB Ludwig，XIE Xiaohai Application experience of lineareddy current braking deviceJForeign Locomotive ollingStock Technology，2013(3):385KUNZ M，王渤洪，黄济荣ICE3 的线性涡流制动:技术观点和运行经验J变流技术与电力牵引，2005(2):4KUNZ M，WANG Bohong，HUANG Jirong Integration of theICE 3s linear eddy-current brake in the infrastructure-technicalaspects and operational experienceJConverter Technology Electric Traction，2005(2):46GABE Johannes，张彦儒 ICE3 列车直线涡流制动装置J国外机车车辆工艺，2003(5):1GABE Johannes，ZHANG Yanru The linear eddy-currentbrake on the ICE3operational concept and first experienceJForeign Locomotive olling Stock Technology，2003(5):17邓长青城轨车辆涡流制动技术研究D成都:西南交通大学，2011DENG Changqingesearch on eddy current braking technologyof urban rail vehiclesDChengdu:Southwest Jiaotong Univer-sity，20118朱仙福，张秀荣高速列车轨道涡流制动的制动力分析与计算J上海铁道大学学报，1996(4):1ZHU Xianfu，ZHANG Xiurong Analysis and calculation of bra-king force on rail eddy current braking of high speed trainsJJournal of Shanghai Tiedao University(Natural Science Edi-tion)，1996(4):19张圣楠永磁涡流制动的电磁分析与设计J内蒙古科技与经济，2005(13):118ZHANG Shengnan Electromagnetic analysis and design of per-manent magnet eddy current brakingJInner Mongolia ScienceTechnology and Economy，2005(13):118 10龙鑫林，张丙飞，佘龙华磁悬浮列车的涡流制动J兵工自动化，2007，26(9):58LONG Xinlin，ZHANG Bingfei，SHE Longhua Eddy brake ofmaglev trainJ Ordnance Industry Automation，2007，26(9):58 11CHEN C，XU J，WU X Analytical calculation of braking forceof super-high speed maglev eddy current braking systemC2019 22nd International Conference on Electrical Machines andSystems(ICEMS)Harbin:IEEE，2019:1 12HEALD M AMagnetic braking:improved theoryJAmericanJournal of Physics，1988，56(6):521(收稿日期:20220824)6