钢轨打磨过程中磨石参数对钢轨温度场影响研究

随着我国铁路运营密度的加大，钢轨疲劳伤损(钢轨剥离、波浪形磨损、斜裂纹、焊缝飞边等[1])成为钢轨的主要伤损形式。目前针对钢轨表面疲劳损伤最有效的修复方法就是进行钢轨打磨[2-4]。钢轨打磨属于复杂的表面材料的去除工艺，包括滑移、变形、剪切、材料去除、热量产生及传递等过程[5-7]。在钢轨打磨过程中，不合理的打磨参数设置将造成大量的磨削热堆积在钢轨轨头部位，这将会使钢轨表面的温度上升[8]。研究学者发现：打磨过程中轨头温度过高将会导致钢轨发蓝[9](图1)，甚至可能改变钢轨材料的微观组织[10-12]，形成打磨马氏体组织。国外学者开展了很多有关磨削温度场的研究:Jaeger[13]首次提出“移动热源法”用以计算剪切滑移区材料的温度场；Outwater等[14]将磨削区视为沿着工件表面移动的热源，使得温度场的计算更准确；Rowe等[15]建立了磨削传热模型并研究了不同磨削参数对温度的影响。国内学者通过仿真和计算研究了不同工艺参数对钢轨打磨温度变化的影响：聂蒙等[16]提出增加砂轮直径和提高作业速度均能降低打磨温度，且打磨头数量的增加有利于打磨温度的稳定；张青等[17]指出打磨车移速越快、磨石转速越低，相应打磨温度越低；Zhang等[18]认为打磨温度随打磨功率的增加而升高，可通过提高磨削区宽度和打磨车移动速度来限制打磨温度的升高。

虽然，国外学者对磨削过程的热影响做了大量研究，国内学者也对钢轨打磨过程的温度场做了探讨。然而，磨石作为打磨机构的重要组成部分，其参数对打磨温度场的影响规律尚不清楚。因此，本文在数值计算基础上，利用Abaqus有限元仿真方法对打磨过程中不同磨石粒度、不同进给深度以及不同磨石转速对打磨温度场的影响进行了分析。仿真过程中参考文献[18]中的移动热源法，将打磨磨石视为移动热源，以此来求解打磨过程中钢轨的温度场。研究结果有助于为现场钢轨打磨选取合适的磨石参数提供理论指导。

图1 钢轨发蓝[9]

1 钢轨打磨过程传热原理

与传统的外圆磨削温度模型不同，钢轨打磨温度模型属于端面磨削[17]。打磨过程中，钢轨瞬态温度场的微分形式为

(1)

式中：ρ为钢轨材料密度；c为比热容；λ为导热系数；τ为打磨时间(从打磨开始计时)；T为钢轨实时温度；div(gradT)为温度梯度/

2T的拉普拉斯算子。可以看出：式(1)为关于时间的一次微分方程，故需设置一定的初始条件。在钢轨打磨前钢轨的初始温度为

T(x,y,z,τ)=T0 τ=0

(2)

式(1)右边是与温度相关的二阶导数，需设置两个边界条件。在磨石与钢轨作用边界上，有热流的流入，其边界条件为

(3)

在其他边界上为

随着公职律师试点在全国的开展，各部门也探索开展了公职律师试点工作，但实际上成效不大。而国外公职律师普遍为专职岗位，不从事其他的行政工作，这是公职律师能够有效发挥作用的重要保证。据了解，现开展试点的税务局中没有一个税务局设立公职律师专职岗位，大多是在法规部门设立公职律师办公室，从其他岗位抽调人员参与公职律师工作，且没有固定工作机制，大多依照税务机关领导的指示来开展工作，因此很难对机关决策提供公正的法律意见，更别提开展监督了。且部分公职律师其行政工作已相当繁重，难以有精力和时间兼负公职律师的职责。所以，要更好地发挥公职律师的作用，必须使公职律师完全脱离其他事务性工作，使其专门从事公职律师工作。

(4)

式中：TE为环境温度；T为钢轨实时温度；n为钢轨界面向外的单位法向；h为对流系数，因对流方式不同，这个系数的取值也有很大的变化；ε为发射率，是表征辐射强度的物理量；σ为斯忒藩-玻尔兹曼常数；q为表面热通量，是打磨过程中磨石产生的打磨功率转化为的热量值。由式(3)和式(4)可以看出：接触区域有热流量的产生，而非接触区域钢轨仅通过辐射与对流与外界进行热量交换。

2 钢轨打磨有限元模型

2.1 打磨过程数值仿真模型

打磨磨石被视为持续发热的移动面热源，移动速度v与打磨车的行驶速度相同，典型的热源分布示意图见图2。由图2可知：热源模型主要有三种，分别是矩形分布、三角形分布和直角三角形分布。参照磨石与钢轨的接触过程，热源被假设成呈矩形分布。

图2 钢轨打磨热源分布模型[18]

钢轨打磨过程中打磨磨石由打磨电机驱动，根据能量守恒定律可知：磨石与钢轨接触区域由摩擦产生的热量是由打磨电机的功率转化而来。因此，磨石所产生的热流密度为

(5)

Q=Pm×τ×η

教学实践是建构多文本阅读课程体系的主战场。恰当、合理、科学的教学模式有利于教学工作的开展，有利于学生学习效率的提高。我们根据不同的文本寻找适合的教学模式，以期达到最好的教学效果。

(6)

s=w×l

(7)

式中：Q为打磨过程总发热量；Pm为钢轨打磨电机功率；s为接触区域面积，即打磨磨石与钢轨接触的长方形区域；w为接触区域宽度，取10 mm[17]；l为接触区域的长度，取100 mm[18](各参数参照图2)；η为总发热量进入钢轨的比例，这个值一般无法定量测量，根据文献[19]η值假设为75%。

2.2 磨石参数与打磨功率关系

钢轨打磨磨石由许多颗磨粒组成。磨粒的粒度由磨粒刚好可通过的筛网尺寸决定。随磨石粒度的增加，磨粒尺寸减小，磨石越细；相反，磨石粒度小，磨粒尺寸增加，对应磨石较粗。可假设磨粒为圆锥形进行磨削力的计算[20]，则单颗磨粒磨削时的磨削力可表示为

Fet=kθ10.84h2.46v-0.299

(8)

式中：k为常数，取值5.4×10-4；θ为磨粒半锥角；h为磨粒切削深度，即每旋转一周，磨石沿垂直打磨面方向上的进给深度；v为磨粒切削速度。

打磨过程中，整块磨石与钢轨的接触形式见图3，其中磨石外径R和内径r间的圆环面为磨石有效打磨区域。

图3 磨石钢轨接触示意[21]

在实际的钢轨打磨作业中，打磨车移动速度远小于磨石的旋转速度[21]，故以磨石的线速度(ωr)取代式(8)中的磨粒切削速度(v)。因此，沿半径方向在磨石的有效打磨区域取微元，微元的磨削力[22]为

儒学在东晋有所复兴，是相对魏晋南北朝尤其汉末、三国这段分裂动荡而又儒学极为衰微的时代而言的，其远未超过两汉这一经学的极盛时代。儒学在东晋的复兴是局部而又微弱的，即所谓儒学有所复兴而不强盛。这种局部而微弱的复兴，蕴育着其深层的自我否定因素与很大的外部挑战。这一两面性，突出表现在东晋的皇帝教育与官学教育上。

dFt=FetCslF(h)dr

(9)

(10)

(11)

式中：Cs为单位面积的磨粒个数；l为接触宽度为b时的接触弧长；F(h)为与进给深度相关的概率函数，当进给深度越深时，参与磨削的磨粒数就越多；μ为磨粒突出高度均值；σ为方差。

The scrap rate from optimized configuration is zero because the profile don′t need to cut. The cost saving based on material price is

（一）医疗服务定价原则与办法缺乏顶层设计。中国地域辽阔，经济发展水平差距大，财政补助政策差异更大，各地区医疗服务价格定价办法有多种，国家卫生、发改、财政等部门虽然出台了《医疗服务价格规范》，但没有统一的定价原则，定价办法。没有调价的程序、调价的依据、调价的范围、调价的幅度、调价的方法规制。从专家发表的研究资料中披露了多种定价方法，已经为政府主管部门制定价格定价原则和办法提供了思路。

利用磨削力式(9)，通过积分法可以得出打磨过程中的磨削力矩为

M=2

FetCslF(h)rdr

(12)

不同力矩下的打磨功率为

(13)

得到打磨功率后，将功率值带入式(5)～式(7)，便可得到相应的热流密度值。

2.3 有限元模型

利用Abaqus建立钢轨打磨有限元模型，选取材料为U71Mn的75 kg/m钢轨作为仿真钢轨模型，钢轨长度3 m，模型见图4。为了在减少计算量和提高求解精度之间找到平衡，通过多次仿真发现：仅建立轨头部分的有限元模型就能计算得到准确的温度场数值，且计算时间更少。使用热传导分析步，通过DFLUX子程序加载移动的表面热流密度，加载平面见图4，模拟打磨过程中磨石在钢轨表面的移动。设定模型的散热边界条件为对流及辐射，设定网格单元类型为计算热传导的DC3D20(20节点二次热传导单元)，并在热流加载区域使用更加精细化的网格。钢轨材料参数及热仿真参数见表1。

图4 钢轨打磨温度场计算有限云模型

表1 钢轨材料及仿真参数

参数数值参数数值密度/(kg·m-3)7 800对流系数·[W/(m2·℃)-1]10热传导系数/[W·(m·K)-1]51发射率0.8比热容/[W·(m·K)-1]500钢轨初始温度及环境温度/℃20

3 磨石参数对打磨钢轨温度场影响

3.1 进给深度

钢轨打磨过程中进给深度一般为8～14 μm。选择12#磨石,磨石转速为3 000 r/min，打磨平面宽度为10 mm，打磨列车速度为2.5 m/s，计算钢轨打磨中的温度变化，并绘制最高温度与进给深度的关系曲线，见图5。由图5可知，打磨钢轨的温度随磨石进给深度的增加快速升高。这是由于较高的打磨进给深度一方面导致单颗磨粒的磨削力增加，另一方面提高了参与磨削的磨粒个数。因此，随进给深度的增加，磨削力矩增加，打磨功率增加，引起更高的钢轨温升。

图5 进给深度对打磨温度影响

当进给深度为14 μm，磨石粒度为12#，磨石转速3 000 r/min，打磨平面宽度10 mm，列车速度5 m/s时选取钢轨打磨区域某一节点为温度监测点，考察该节点在打磨过程中的温度变化见图6。由图6可知，打磨开始后，节点温度先迅速升高到最高温度，然后慢慢下降，整个过程与打磨磨石相对该节点的位置变化相对应：当磨石靠近该节点时，打磨过程中的热效应导致该节点温度迅速上升；当磨石全部通过该节点时，打磨温度上升至535.7 ℃的最高值；当磨石远离该节点时，打磨温度逐渐下降。

图6 打磨过程中某节点的温度变化

打磨过程中整个钢轨的最高温度变化见图7。由于磨石转速较高，因此在打磨开始后，钢轨整体最大温度在0.04 s内迅速上升至535.7 ℃，并持续保持在稳定状态。由图7可知，钢轨温度迅速上升是磨石对钢轨持续的热流加载造成的，而当热流流入和热量损失达到平衡时，打磨热进入平稳状态，钢轨的最大温度保持稳定。

图7 钢轨整体温度随打磨时间变化

打磨过程中，钢轨整体温度仿真云图见图8。由图8(a)可知，钢轨温度呈椭圆形分布，在椭圆的中心部位最高，并向四周逐渐降低。这是由于打磨热产生于磨石与钢轨的接触区域，钢轨在越靠近接触区域的位置，其热传递效率越高，相应的温升越高、温度越大；而钢轨在远离接触区域的位置，温升主要来自于钢轨内部的热传导，由于传导过程中的热量损失，越远离接触区域的位置，对应温度越低。由图8(b)可知，温度云图呈碗状，温度在中心位置最高，并在远离中心位置的区域逐渐降低。

图8 打磨过程瞬时温度云图(单位：℃)

3.2 磨石转速

转速参数1 500～3 000 r/min、12#磨石,进给深度12 μm,打磨平面宽度10 mm,列车前进速度2.5 m/s下仿真得到打磨温度与磨削力矩随磨石转速变化的曲线，见图9。由图9可知，磨削力矩随磨石转速增加而降低，这是由于单颗磨粒的磨削力随磨石转速的升高而降低，导致磨石整体的磨削力矩下降。由图9可知，随打磨转速增加，打磨温度反而升高。由式(8)及式(13)中看出，这是由于转速升高对功率增加较大，而对磨削力下降的影响较小。

图9 磨石转速对打磨温度与磨削力矩影响

3.3 磨石粒度

明确磨石粒度与打磨温度之间的关系可以指导在现场打磨作业中选取合适粒度的打磨磨石。根据打磨磨石常用的粒度，选取12#、16#、24#、30#四种粒度参数进行仿真。其他仿真参数设置包括：进给深度为14 μm、磨石转速为3 000 r/min、打磨平面宽度为10 mm，列车前进速度为2.5 m/s。根据式(8)～式(13)和各粒度的特征参数[17]计算得不同磨石粒度对应的最大打磨温度，见图10。结果表明：随磨石粒度增加，钢轨温度呈现先上升后下降的变化趋势，四种磨石粒度中，16#粒度对应的最大温度最高。这一结果与文献[22]中“砂轮磨粒对磨削温度的影响”的有关描述相同，体现了仿真分析的正确性。这一规律与磨石粒度特征值的差异性有关：随着磨石粒度的增加，磨粒尺寸和半锥角减小，根据式(8)计算得对应的单颗磨粒磨削力减小；但随着磨石粒度的增加，磨粒数量也会增加，导致磨石整体的磨削力增加。因此，在12#～16#的粒度范围内，磨粒数量的增加相较单颗磨粒磨削力的减少对磨石整体磨削力的影响更大，而在16#～30#的粒度范围内则相反。

图10 磨石粒度与打磨温度关系曲线

4 结论

(1)基于单颗磨粒磨削力公式，利用数值积分方法，对打磨过程中不同粒度磨石所产生的打磨功率进行了计算。利用移动热源法将磨石视为在钢轨表面移动的热源，将功率等效为磨石打磨产生的热量加载在钢轨有限元模型上，对打磨过程中钢轨温度场进行了有限元仿真。

(2)当磨石进给深度从8 μm增加至14 μm，单颗磨粒的磨削力和参与磨削的磨粒个数均增加。导致打磨温度从164.5 ℃升高至535.7 ℃；由于转速对磨削力下降的贡献远小于直接对功率增加的贡献，因此，虽然在转速从1 500 r/min升高至3 000 r/min过程中磨削力矩从77.2 N·m下降至62.7 N·m，但温度却从242.5 ℃升高至352.9 ℃。

（1）全面做好前期性的水利施工预备。从前期开展的各项施工预备角度来讲，关键在于筹备当前水利施工必需的各类物资、人力资源以及施工资金等，确保做好妥善的前期施工预备。与此同时，技术人员有必要给出相应的水利设计图纸，并且拟定可行性较强的水利施工规划。针对各类必需的施工设备来讲，应当做好妥善性的设备养护以及设备检测。在此前提下，关于各类建材也要做好综合性的收管工作以及采购审核工作，确保能够录入详细性的产品资料与产品信息数据。

(3)随磨石粒度的增加，磨粒的半锥角减小，单颗磨粒的磨削减小，但磨石数量也会增加，导致磨石整体的磨削力矩增大，这种“竞争”的关系使得在磨石粒度从12#增加到16#过程中，打磨温度从535.7 ℃升高至656.7 ℃，当粒度从16#增加到30#过程，温度却逐渐下降到413.7 ℃。

摄像机参数及光靶结构尺寸给定后，影响解算模型的参数只有横向偏移量、光靶到摄像机的距离及光靶的偏转角等，为此，在Matlab软件中对解算模型进行仿真以了解各参数对模型的影响程度。

③疗效判定有效:自觉症状（疼痛、麻木、感觉异常）较前改善，体格检查（膝反射、跟腱反射）较前好转，但未达正常，神经传导速度较治疗前虽有所提升，但速度提升在5 m/s以下。显效:感觉无明显症状或症状有显著好转，膝反射及跟腱反射基本正常，神经传导速度接近或能够达到正常或加快速度≥5 m/s。无效:自觉症状无好转，膝反射及跟腱反射无好转，神经传导速度较治疗前无好转。显效+有效=总有效。

参考文献：

[1]刘月明，李建勇，蔡永林，等.钢轨打磨技术现状和发展趋势[J]．中国铁道科学，2014，35(4)：29-37.

LIU Yueming，LI Jianyong，CAI Yonglin，et al．Current State and Development Trend of Rail Grinding Technology[J].China Railway Science，2014，35(4)：29-37.

[2]GU K，LIN Q，WANG W，et al．Analysis on the Effects of Rotational Speed of Grinding Stone on Removal Behavior of Rail Material[J].Wear，2015，342：52-59.

[3]SATOH Y K，IWAFUCHI K．Effect of Rail Grinding on Rolling Contact Fatigue in Railway Rail Used in Conventional Line in Japan[J].Wear，2008，265(9/10)：1342-1348.

[4]KALOUSEK J，SROBA P，HEGELUND C．Analysis of Rail Grinding Tests and Implications for Corrective and Preventative Grinding[C]//The Fourth International Heavy Haul Railway Conference.Brisbane:IHHA，1989：193-204.

[5]GU R J，SHILLOR M，BARBER G C，et al．Thermal Analysis of the Grinding Process[J]．Mathematical & Computer Modelling，2004，39(9/10)：991-1003.

[6]HADAD M，SADEGHI B．Thermal Analysis of Minimum Quantity Lubrication-MQL Grinding Process[J].International Journal of Machine Tools & Manufacture, 2012，63：1-15.

[7]GARC

A E，MÉRESSE D，POMBO I，et al．Identification of Heat Partition in Grinding Related to Process Parameters，Using the Inverse Heat Flux Conduction Model[J].Applied Thermal Engineering，2014，66(1/2)：122-130.

[8]FANG C F，XU X P．Analysis of Temperature Distributions in Surface Grinding with Intermittent Wheels[J].The International Journal of Advanced Manufacturing Technology，2014，71(1)：23-31.

[9]RASMUSSEN C J，FAESTER S，DHAR S，et al.Surface Crack Formation on Rails at Grinding Induced Martensite White Etching Layers[J].Wear，2017，384/385：8-14.

[10]STEENBERGEN M．Rolling Contact Fatigue in Relation to Rail Grinding[J].Wear，2016，356/357：110-121.

[11]ZHANG J Z，TAN X M，LIU B，et al．Investigation for Convective Heat Transfer on Grinding Work-Piece Surface Subjected to a Mist/Air Impinging Jet [J]．Applied Mechanics and Materials，2012，249/250：434-442.

[12]TAHVILIAN A M，LIU Z H，CHAMPLIAUD H，et al．Experimental and Finite Element Analysis of Temperature and Energy Partition to the Workpiece While Grinding with a Flexible Robot[J].Journal of Materials Processing Technology，2013，213(12)：2292-2303.

[13]JAEGER J G．Moving Sources of Heat and the Temperature of Sliding Contacts[J].J. & Proc. Roy. Roc．New South Wales，1942，76：203-224.

[14]OUTWATER J O，SHAW M C．Surface Temperatures in Grinding[J].Transactions of the ASME，1952，74：73-78.

[15]ROWE W B，MORGAN M N，BLACK S C E，et al．A Simplified Approach to Control of Thermal Damage in Grinding[J].Cirp Annals，1996，45(1)：299-302.

[16]聂蒙，李建勇，沈海阔，等．钢轨打磨作业过程中磨削温度建模与仿真[J]．铁道学报，2013，35(10)：89-93.

NIE Meng，LI Jianyong，SHEN Haikuo，et al．Modeling and Simulation of Temperature in Process of Rail Grinding Operation[J].Journal of the China Railway Society，2013，35(10)：89-93.

[17]张青，崔晓璐，陈光雄，等．钢轨打磨热力耦合分析[J]．铁道学报，2015，37(5)：78-82.

ZHANG Qing，CUI Xiaolu，CHEN Guangxiong，et al．Thermal-mechanical Coupling Analysis in Rail Grinding[J].Journal of the China Railway Society，2015，37(5)：78-82.

[18]ZHANG Z J，SHANG W，DING H H，et al．Thermal Model and Temperature Field in Rail Grinding Process Based on a Moving Heat Source[J].Applied Thermal Engineering，2016，106：855-864.

[19]ROWE W B．Principles of Modern Grinding Technology[M].Oxford：William Andrew，2009.

[20]韩永超．基于钢轨打磨实验的磨削力影响因素仿真与实验研究[D]．成都：西南交通大学，2016.

[21]张科元．钢轨打磨小车动力学行为研究[D]．成都：西南交通大学，2016.

[22]任敬心，华定安．磨削原理[M]．北京：电子工业出版社，2011.