一种钢轨服役安全性的计算方法及装置与流程

1.本发明涉及铁路安全运行技术领域，特别涉及一种钢轨服役安全性的计算方法及装置。


背景技术：

2.在铁路线路养护过程中，钢轨踏面低塌带来了许多线路病害，尤其是随着无缝线路的大量铺设，焊接接头由于硬度较低而出现低塌的现象普遍。踏面低塌使钢轨的连续性受到了影响，列车通过低塌处时钢轨会受到很大的冲击附加动力作用，根据对车轮簧下质量加速度的实测资料，这些附加力为正常荷载的2～3倍，严重时可达4～5倍。例如对于长度为1m、最大值为1mm的凹陷，轮重为7吨的列车以200km/h经过时轮重的最大值便超过了200kn，达到了静态轮重的3倍。随着接头低塌的发展，车轮对钢轨的附加冲击力不断增大，而巨大的冲击动力又加快了接头病害的发展，轨道发生强迫振动，使得钢轨变形的发展加快，轨道状态加速恶化。
3.此外，钢轨锈蚀是铁路线路上最常见的伤损，其中轨底锈蚀尤为严重。轨底存在较深的伤损坑会产生应力集中效应，因此由轨底锈蚀引发的钢轨失效案例在全部线路钢轨失效案例中占比最高，达到了约49％，对线路安全构成了很大的威胁。对于运煤线路钢轨，其长期处于煤粉及防冻液污染腐蚀环境中，钢轨轨底易于锈蚀而存在大量的伤损坑，此处成为应力集中点而极易萌生裂纹。
4.在钢轨失效案例中，发现包括重载、普速铁路在内的数起断裂失效钢轨存在一个共同特点：裂纹起源于轨底下表面的伤损坑处，而在伤损坑附近的闪光焊接头都存在着低塌现象。因而，踏面低塌对附近轨底伤损坑受力及钢轨服役状态的影响需引起重视。但是，现有技术中对于钢轨踏面低塌和附近轨底蚀坑两处伤损共同作用下的钢轨服役安全性的评估，目前尚无相关的研究和计算方法。
5.另外，现有技术中关于钢轨存在踏面凹陷时产生的钢轨附加冲击力的理论公式计算方法，没有考虑车辆自身重量的影响，更没有考虑对轨底伤损坑处受力的影响。同时，现有技术中关于钢轨存在踏面凹陷时产生的钢轨附加冲击力一般采用有限元仿真计算方法，建模及前处理操作过程复杂，计算耗时长，需要具有较强的仿真领域知识技能，不适用于铁路养护人员快速便捷的计算分析。并且，现有技术中关于钢轨不平顺状态下的钢轨受力计算方法，虽然考虑了车辆自身重量的影响，但其准静态计算方法并不能将钢轨存在明显踏面低塌伤损的情况包括在内，其基于经验值的动力增值系数也无法对低塌激励下的钢轨附加冲击力进行准确计算。
6.有鉴于此，本发明人根据多年从事本领域和相关领域的生产设计经验，经过反复实验设计出一种钢轨服役安全性的计算方法及装置，以期解决现有技术存在的问题。


技术实现要素：

7.本发明提出一种钢轨服役安全性的计算方法及装置，能够计算和评价钢轨同时存
在踏面低塌和附近轨底伤损坑两种伤损情况下的钢轨服役安全性。
8.为达到上述发明目的，本发明提出一种钢轨服役安全性的计算方法，其中，所述计算方法包括：
9.获取车辆载荷作用下钢轨伤损坑处的轨底最大动弯应力、踏面低塌产生的附加冲击作用在伤损坑处的附加动弯应力、以及钢轨伤损坑处的轨底所受温度应力；
10.基于所述轨底最大动弯应力、所述附加动弯应力和所述轨底温度应力，计算踏面低塌影响下钢轨伤损坑底部的最大拉应力；
11.根据所述最大拉应力及相应损伤下的疲劳极限评价钢轨服役安全性。
12.本发明还提出一种钢轨服役安全性的计算装置，其中，所述装置包括：
13.轨底最大动弯应力计算模块，通过所述钢轨截面静弯矩和所述轨底端面系数计算所述轨底最大动弯应力；
14.附加动弯应力计算模块，根据所述附加冲击力、所述钢轨基础参数以及所述水平距离系数计算所述附加冲击作用在伤损坑的附加动弯应力；
15.轨底所受温度应力，根据所述锁定轨温和所述测量轨温计算所述轨底温度应力；
16.最大拉应力计算模块，基于所述轨底最大动弯应力、所述附加动弯应力和所述轨底温度应力，计算踏面低塌影响下钢轨伤损坑底部的最大拉应力；
17.安全性评价模块，根据所述最大拉应力及相应损伤下的疲劳极限评价钢轨服役安全性。
18.本发明还提出一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如上任一所述方法。
19.本发明还提出一种计算机可读存储介质，其中，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上所述方法。
20.与现有技术相比，本发明提出的钢轨服役安全性计算方法具有如下特点和优点：
21.本发明提出的钢轨服役安全性计算方法考虑了钢轨存在焊接接头低塌及附近轨底伤损坑两种伤损的叠加影响，而非某一种伤损的单一影响；并且，通过建立的理论公式模型，评估了两种伤损叠加作用下的轨底伤损坑受力状态及钢轨服役安全性，能够有效解决铁路现场养护人员对此种状况钢轨服役安全状态的关切。
22.本发明提出的钢轨服役安全性计算方法对于轨底伤损坑受力状态及钢轨服役安全性的计算，相比于可以采用的有限元仿真方法，可以更为简单便捷的实现，不需要漫长的建模及运算时间及较高的有限元仿真知识技能门槛，更适用于现场铁路养护人员进行快速评估。
23.本发明提出的钢轨服役安全性计算方法通过计算不同数值的踏面低塌深度及轨底伤损坑应力集中系数，能够分析随着两处伤损的发展，钢轨伤损坑受力及疲劳极限的演变；并能够通过伤损坑受力与疲劳极限的对比，对现实中的类似钢轨失效案例的失效原因分析提供可靠的数据支撑。
附图说明
24.在此描述的附图仅用于解释目的，而不意图以任何方式来限制本发明公开的范
围。另外，图中的各部件的形状和比例尺寸等仅为示意性的，用于帮助对本发明的理解，并不是具体限定本发明各部件的形状和比例尺寸。本领域的技术人员在本发明的教导下，可以根据具体情况选择各种可能的形状和比例尺寸来实施本发明。
25.图1为钢轨轨底伤损坑受力组成示意图；
26.图2为钢轨受车辆载荷分布示意图；
27.图3a为钢轨踏面低塌形状示意图(一)；
28.图3b为钢轨踏面低塌形式示意图(二)；
29.图4为伤损坑示意图；
30.图5为本发明提出的钢轨服役安全性计算方法流程图；
31.图6为本发明实施例中钢轨受车辆载荷分布示意图；
32.图7为本发明实施例中失效钢轨轨底伤损坑示意图。
具体实施方式
33.结合附图和本发明具体实施方式的描述，能够更加清楚地了解本发明的细节。但是，在此描述的本发明的具体实施方式，仅用于解释本发明的目的，而不能以任何方式理解成是对本发明的限制。在本发明的教导下，技术人员可以构想基于本发明的任意可能的变形，这些都应被视为属于本发明的范围。
34.在本发明中，当元件被称为“固定于”或“设置于”另一个元件，它可以直接在另一个元件上或者可能存在居中元件。当一个元件被称为是“连接于”另一个元件，它可以是直接连接到另一个元件或者可能存在居中元件。
35.一种钢轨服役安全性的计算方法，如图1至图5所示，该计算方法包括：
36.获取车辆载荷作用下钢轨伤损坑处的轨底最大动弯应力、踏面低塌产生的附加冲击作用在伤损坑处的附加动弯应力、以及轨底所受温度应力；
37.基于轨底最大动弯应力、附加动弯应力和轨底温度应力，计算踏面低塌影响下钢轨伤损坑底部的最大拉应力；
38.根据所述最大拉应力及相应损伤下的疲劳极限评价钢轨服役安全性。
39.本发明提出的钢轨服役安全性计算方法考虑了钢轨存在焊接接头低塌及附近轨底伤损坑两种伤损的叠加影响，而非某一种伤损的单一影响；并且，通过建立的理论公式模型，评估了两种伤损叠加作用下的轨底伤损坑受力状态及钢轨服役安全性，能够有效解决铁路现场养护人员对此种状况钢轨服役安全状态的关切。
40.本发明提出的钢轨服役安全性计算方法对于轨底伤损坑受力状态及钢轨服役安全性的计算，相比于可以采用的有限元仿真方法，可以更为简单便捷的实现，不需要漫长的建模及运算时间及较高的有限元仿真知识技能门槛，更适用于现场铁路养护人员进行快速评估。
41.在本发明一个可选的实施方式中，踏面低塌包括焊接接头处低塌、钢轨踏面隐伤、擦伤、波磨等损伤引起的低塌、以及其他近似半波正弦形曲线形状的低塌。
42.在本发明一个可选的实施方式中，轨底伤损坑可以为锈蚀坑，也可以为钢轨在运输铺设过程中导致的轨底磕碰伤坑，以及钢轨生产过程中热状态下形成的轨底热刮伤、磨痕、折叠、氧化皮压入、轧痕等形成的伤损坑和冷状态下形成的划痕等形成的伤损坑。
43.在本发明一个可选的实施方式中，获取车辆载荷作用下钢轨伤损坑处的轨底最大动弯应力，包括；
44.获取钢轨伤损坑处的钢轨截面静弯矩和钢轨伤损坑处的轨底端面系数，通过钢轨截面静弯矩和轨底端面系数计算轨底最大动弯应力，则有，
[0045][0046]
式中，σ0为轨底最大动弯应力，m0为钢轨截面静弯矩，w1为轨底端面系数。
[0047]
在一个可选的例子中，w1可由标准tb 2034—1988附录中查得。
[0048]
需要说明的是，由于焊接接头不平顺激励的轮轨动态响应单独计算，因此，在上述公式(3)中车辆载荷对钢轨的作用不考虑车辆行驶中动力增值的影响。
[0049]
在该实施方式一个可选的例子中，获取所述钢轨伤损坑处的钢轨截面静弯矩，包括，
[0050]
确定车辆的最不利轮位、以及该最不利轮位的载荷；
[0051]
获取钢轨基础参数，以及钢轨伤损坑与所低塌的水平距离；
[0052]
根据最不利轮位的载荷、钢轨基础参数和所述水平距离，计算钢轨截面静弯矩。
[0053]
在一个可选的例子中，钢轨基础参数包括钢轨基础与钢轨刚比系数、钢轨基础弹性系数和钢轨竖向抗弯刚度，则有，
[0054][0055]
式中，m0为钢轨截面静弯矩，∑pμ为各车轮的当量载荷，p为最不利轮位的载荷，μ＝e-βx
(cosβx-sinβx)，x为伤损坑与低塌的水平距离，β为钢轨基础与钢轨刚比系数，sinβx)，x为伤损坑与低塌的水平距离，β为钢轨基础与钢轨刚比系数，钢轨基础弹性系数，ej
x
为钢轨竖向抗弯刚度。
[0056]
在本发明中，车辆载荷作用对伤损坑处钢轨截面静弯矩，采用轨道强度检算的基本力学原理，将钢轨视为连续弹性基础梁模型，且不考虑轨道自重，通过静态计算方法确定轨底的动弯应力。根据winkler假定进行求解，得到伤损坑处的钢轨截面静弯矩m0。
[0057]
在一个可选的例子中，最不利轮位，为车辆其前车的后转向架与后车的前转向架共四个轮位分别经过低塌处时，在伤损坑处产生的最大载荷的轮位。如图2所示，前车的后转向架与后车的前转向架共四个轮位(ⅰ～ⅳ)组成一单元，考虑到计算轮及其前后相临轮的影响，并根据力的独立作用原理，把各车轮作用在计算截面的力叠加起来，即分别求得各车轮当作计算轮时截面的合力，从而确定最不利轮位，即∑pμ最大的轮位。
[0058]
在另一个可选的例子中，最不利轮位，为六轴或八轴的车辆中，各车轮分别经过低塌处时，在伤损坑处产生的最大载荷的轮位。
[0059]
特别需要说明的是，对导致轨底损伤坑产生最大应力和最不利轮位的计算，定义为各分别经过低塌处时，而非经过锈蚀坑正上方时，在损伤坑处而非在低塌处产生的最大σpμ值，及产生此值时的计算轮位。
[0060]
在本发明一个可选的实施方式中，获取踏面低塌产生的附加冲击作用在轨底的附加动弯应力，包括，
[0061]
获取车辆的运行速度和踏面低塌的形状参数；
[0062]
基于踏面低塌的形状参数和车辆的运行速度，计算车辆经过踏面低塌时产生的附加冲击力；
[0063]
获取钢轨基础参数，以及钢轨伤损坑与所低塌的水平距离，并根据附加冲击力、钢轨基础参数以及水平距离系数计算附加冲击作用在伤损坑的附加动弯应力。
[0064]
在该实施方式一个可选的例子中，将踏面低塌形状看作半波正弦形曲线，如图3a、图3b所示，踏面低塌的形状参数包括最大低塌量、低塌长度和激励频率，则有，
[0065][0066]
式中：λ为最大低塌量；l为低塌长度，ω为激励频率，且ω＝2πv/l，v为行车速度，t为从起点算起的时间，t＝l/v；l为低塌距起点的水平距离，η为车轮经过低塌的某一时刻t时对应的轨面不平顺值。
[0067]
在一个可选的例子中，基于踏面低塌的形状参数和车辆的运行速度，计算踏面低塌产生的附加冲击力，则有，
[0068][0069]
式中：y为钢轨动挠度，表示车轮经过低塌的时间，表示簧下质量在轨道上自由振动的周期，为簧下质量的固有圆频率，k为轨下基础换算刚度，m为簧下质量，t为从起点算起的时间；
[0070]
需要说明的是，对于踏面低塌产生的动力响应，根据达朗贝尔原理建立的钢轨挠曲振动方程，以及边界条件(t＝0时，y＝0，dy/dt＝0)，可求解挠曲振动方程，即钢轨动挠度y，并且，为简化计算，不考虑轮轨接触变形，略去钢轨质量和摩擦力，则，钢轨动挠度y的计算公式便为公式(5)。
[0071]
因其中为钢轨基础弹性系数，lk为弹簧的有效换算长度，从而作用于钢轨上的附加动力pc为：
[0072]
pc＝ky＝8ej
x
β3y
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0073]
式中，k为轨下基础换算刚度，y为钢轨动挠度，ej
x
为钢轨竖向抗弯刚度，β为钢轨基础与钢轨刚比系数，且基础与钢轨刚比系数，且为钢轨基础弹性系数。
[0074]
进一步的，根据附加冲击力、钢轨基础参数以及水平距离系数计算附加冲击作用在伤损坑的附加动弯应力，则有，
[0075][0076]
式中，σ
cd
为附加动弯应力，pc为附加动力，μ＝e-βx
(cosβx-sinβx)，x为伤损坑与低塌的水平距离，β为钢轨基础与钢轨刚比系数，塌的水平距离，β为钢轨基础与钢轨刚比系数，钢轨基础弹性系数，ej
x
为钢轨竖向抗弯刚度，且β、和ej
x
属于钢轨基础参数，w1为轨底端面系数。
[0077]
需要说明的是，根据低塌产生的附加冲击力pc，由上述公式(2)和上述公式(3)便可获得上述公式(7)。
[0078]
在本发明一个可选的实施方式中，获取轨底所受温度应力，包括：
[0079]
获取钢轨的锁定轨温和测量轨温，根据锁定轨温和测量轨温计算所述轨底温度应力，则有，
[0080]
σ
t
＝2.45δt
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0081]
式中，σ
t
为轨底所受的温度应力，δt为所述锁定轨温与所述测量轨温的差值。
[0082]
在本发明一个可选的实施方式中，计算钢轨伤损坑底部的最大拉应力，包括：
[0083]
基于轨底最大动弯应力、附加动弯应力和轨底温度应力计算踏面低塌影响下钢轨锈蚀处的轨底最大应力；
[0084]
获取轨底伤损坑的形状参数，并根据该形状参数计算伤损坑的应力集中系数；
[0085]
根据钢轨锈蚀处的轨底最大应力和应力集中系数计算踏面低塌影响下钢轨伤损坑底部的最大拉应力。
[0086]
在该实施方式一个可选的例子中，基于轨底最大动弯应力、附加动弯应力和轨底温度应力计算踏面低塌影响下钢轨锈蚀处的轨底最大应力，则有，
[0087]
σd＝σ0+σ
cd
+σ
t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0088]
式中，σd为轨底的最大应力，σ0为轨底最大动弯应力，σ
cd
为附加动弯应力，σ
t
为轨底所受的温度应力。
[0089]
在本发明中，无缝线路钢轨的轨底受力由最大可能载荷(动载荷)作用下的轨底动弯应力及温度应力构成。其中最大可能载荷为考虑列车行驶中的动力作用使得轨道承受超出静载荷的动力增加值。对于失效轨存在较严重接头低塌的情况，为准确计算低塌附近轨底的受力状态，将轮轨系统在低塌激励下产生的垂向动态响应进行单独计算，因而车辆载荷作用在轨底的应力则通过静态计算获得。
[0090]
因此，对于低塌附近伤损坑处的轨底最大拉应力σd，如图1所示，应包括车辆静载荷作用在轨底的最大动弯应力σ0、车辆经过低塌时产生的垂向附加冲击作用在轨底的动弯应力σ
cd
以及温度应力σ
t
。
[0091]
进一步的，若分析钢轨的轨底锈蚀坑位于轨底中心线附近，或锈蚀坑离接头低塌处的焊缝距离较远，或不是接头低塌的其他踏面低塌情况，此时钢轨内部还存有残余应力，此情况下可在公式(1)的等号右侧再加上残余应力一项，也可同样采用本发明的计算方法对钢轨的服役安全性进行分析。
[0092]
在一个可选的例子中，轨底伤损坑的形状参数包括伤损坑的深度和伤损坑根部半径，则有，
[0093][0094]
式中，k
t
为应力集中系数，h为伤损坑的深度；ρ为伤损坑根部半径。
[0095]
在本发明中，对轨底损伤坑形状参数的确定，将轨底损伤坑的形状定义为在平行于钢轨长度方向截面上投影的最大轮廓。损伤坑深度h和底部半径ρ则在最大轮廓上测量得到。对轨底锈蚀坑应力集中系数的计算，将轨底伤损坑视为在平板边缘存在的一个缺口，根据断裂力学，对于一个无限长度的平板，如图4所示，沿长度方向存在一个远端拉应力，在平板边缘存在一个缺口，可将缺口根部形状等效为椭圆，从而可由公式(9)计算损伤坑的应力集中系数。
[0096]
同时，由于轨底伤损坑尺度相对于钢轨非常小，因此可将轨底伤损坑附近的应力视为沿钢轨纵向均匀分布的拉应力，拉应力方向为钢轨长度方向。测量目标伤损坑或断口上作为裂纹源的伤损坑由平行于钢轨长度方向上投影的最大轮廓取得，测量伤损坑最大轮廓的深度h和底部半径ρ，然后根据公式(9)便可计算伤损坑的应力集中系数k
t
。
[0097]
进一步的，已知钢轨锈蚀处的轨底最大应力和应力集中系数，则踏面低塌影响下钢轨伤损坑底部的最大拉应力为，
[0098]
σ
pit
＝k
t
·
σdꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0099]
式中，σ
pit
为伤损坑底部的最大拉应力，k
t
为应力集中系数，σd为伤损坑处轨底的最大应力。
[0100]
在本发明一个可选的实施方式中，根据所述最大拉应力及相应损伤下的疲劳极限评价钢轨服役安全性，包括，
[0101]
确认钢轨伤损坑底部的最大拉应力与钢轨在相应损伤下的疲劳极限是否符合判定标准公式，若符合所述判定标准公式，则判断钢轨处于安全服役状态；若不符合所述判定标准公式，则判断钢轨处于不安全服役状态。
[0102]
在该实施方式一个可选的例子中，判定标准公式为：
[0103]
σ
pit
≤σr·kt
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0104]
式中，σ
pit
为伤损坑底部的最大拉应力，k
t
为应力集中系数，σr为相应伤损下的疲劳极限。
[0105]
具体的，在钢轨服役过程中，钢轨在焊接接头踏面低塌及附近轨底伤损坑两处伤损的叠加影响下，会降低钢轨的疲劳寿命。当两处伤损达到一定程度时，伤损坑处受力较大而高于疲劳极限时，钢轨就会处于不安全的服役状态。因此，将轨底伤损坑处最大受力σ
pi
t与相应伤损下的疲劳极限σr进行比较，以对钢轨服役安全性进行判定。
[0106]
在该实施方式一个可选的例子中，相应损伤下的疲劳极限为包含踏面与轨底伤损坑影响的钢轨疲劳极限，则有，
[0107]
σr＝aσy+be+c
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)
[0108]
且，
[0109][0110]
c＝12.0k
t
r+28.2k
t-82.8r-191，
[0111]
式中：σy为钢轨的屈服强度；e为钢轨弹性模量；r为伤损坑处循环拉应力的应力比。
[0112]
对于服役钢轨，当车辆经过时每一个车轮行驶过都会对钢轨产生一次碾压，因而钢轨轨底会承受循环的动弯应力，在轨底伤损坑处则产生循环的拉应力作用，所以包含轨底伤损坑和踏面低塌的钢轨疲劳极限，通过多参数经验公式(11)进行评估。
[0113]
进一步的，伤损坑处循环拉应力的应力比r与踏面低塌相关，则有，
[0114][0115]
式中，σ
max
为车轮经过踏面低塌的时刻，伤损坑处的最大受力，σ
min
为轨底仅受温度应力而无低塌及车辆自重作用下的伤损坑处的受力，σ
t
为轨底所受的温度应力，σd为轨底的最大应力。
[0116]
即将σ
max
定义为车轮经过踏面时锈蚀坑处产生的最大应力σd；σ
min
定义为钢轨所受温度应力σ
t
，伤损坑处循环拉应力的应力比r为与踏面低塌λ相关的参数。
[0117]
本发明还提出一种钢轨服役安全性的计算装置，其中，该装置包括：
[0118]
轨底最大动弯应力计算模块，通过钢轨截面静弯矩和轨底端面系数计算轨底最大动弯应力；
[0119]
附加动弯应力计算模块，根据附加冲击力、钢轨基础参数以及水平距离系数计算附加冲击作用在伤损坑的附加动弯应力；
[0120]
轨底所受温度应力，根据锁定轨温和测量轨温计算轨底温度应力；
[0121]
最大拉应力计算模块，基于轨底最大动弯应力、附加动弯应力和轨底温度应力，计算踏面低塌影响下钢轨伤损坑底部的最大拉应力；
[0122]
安全性评价模块，根据最大拉应力评价钢轨服役安全性。
[0123]
本发明还提出一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，该处理器执行所述计算机程序时实现如上任一所述方法。
[0124]
本发明还提出一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现如上所述方法。
[0125]
现结合一实施例，详细说明本发明提出的钢轨服役安全性的计算方法及装置的具体实施过程：
[0126]
在本实施例中，失效钢轨为75kg/m钢轨，材质为某钢轨钢，热处理态，断裂位于线路的直线地段。如图5所示，钢轨上存在从轨底下表面萌生的与钢轨长度方向垂直的疲劳裂纹，距离裂纹60mm(x)处为一移动闪光焊接头，该接头踏面已发生低塌，最大低塌量λ为0.7mm，低塌范围l为90mm。轨底下表面的裂纹萌生点偏离轨腰中心线15mm。
[0127]
步骤1，计算车辆载荷作用下的伤损坑处轨底最大动弯应力σ0[0128]
该失效钢轨线路上运行的车辆是27t轴重c80e货车，结合前车的后转向架与后车的前转向架的载荷，寻找引起最大弯矩的最不利轮对，如图6所示。计算ⅰ～ⅳ车轮分别经过低塌时，蚀坑位置处的σpμ值，从而确定在蚀坑处引起最大弯矩的最不利轮对，结果见表1。
[0129]
表1各轮位σpμ计算
[0130]
[0131]
可看出当ⅲ轮经过低塌时，σpμ值最大，为100.8kn，因此其为最不利轮位。依据公式(2)～公式(3)进行计算，可得伤损坑处的轨底最大动弯应力σ0＝35.9mpa。
[0132]
相比之下，当接头踏面无低塌(λ＝0mm)时，车轮经过蚀坑正上方时轨底应力达到最大，按照表1的方式计算各轮位经过蚀坑上方时的σpμ值。结果表明当ⅱ(ⅲ)轮经过低塌时σpμ值最大，此时轨底最大动弯应力σ
0(λ＝0)
＝43.1mpa。
[0133]
步骤2，计算踏面低塌激励的附加冲击作用在轨底的动弯应力σ
cd
[0134]
依据公式(4)～公式(7)进行计算，可得低塌激励的附加冲击力pc为144.6kn下，伤损坑处轨底的附加动弯应力σ
cd
＝43.4mpa。
[0135]
步骤3，计算轨底温度应力σ
t
[0136]
失效钢轨线路的锁定轨温为28℃，发现钢轨产生裂纹时测得轨温为7℃。由公式(8)计算得到失效轨在当时温度下所受的温度应力σ
t
＝51.5mpa。
[0137]
步骤4，计算伤损坑处的轨底最大应力σd[0138]
由公式(1)及上面计算结果，可得车轮经过低塌时，附近蚀坑位置处的轨底应力最大，σd＝130.8mpa。
[0139]
当接头踏面无低塌(λ＝0mm)时，车轮经过蚀坑时的轨底应力位置处受力则只包含σ
0(λ＝0)
和σ
t
，此时轨底最大应力σ
d(λ＝0)
＝94.6mpa。
[0140]
步骤5，计算轨底伤损坑应力集中系数k
t
及最大应力σ
pit
[0141]
经测量，如图7所示，轨底裂纹萌生处的伤损坑深度h为0.43mm，底部半径ρ为0.20mm。
[0142]
依据公式(9)～公式(10)，可得伤损坑应力集中系数k
t
为3.92，轨底蚀坑底部所受的最大应力σ
pit
＝511.8mpa。而无低塌时，车轮经过蚀坑上方时蚀坑底部所受的σ
pit(λ＝0)
＝370.3mpa。
[0143]
步骤6，计算包含踏面低塌和轨底伤损坑影响的钢轨疲劳极限σr[0144]
依据公式(11)～公式(12)，计算含伤损坑钢轨的疲劳极限σr，结果如表2所示。
[0145]
表2含伤损坑钢轨轨底在不同踏面状态下的疲劳极限σr[0146][0147]
步骤7，钢轨服役安全性评估
[0148]
比较相同伤损状态下的σ
pit
及σr·kt
，评估钢轨服役安全性。由表2可知，在现有伤损坑及0.7mm踏面低塌情况下，σ
pit
大于σr·kt
，不满足公式(13)的判定标准，因此钢轨在服役中处于不安全状态，疲劳裂纹易于从轨底应力集中处即伤损坑底部萌生。
[0149]
相比之下，当踏面无低塌时，σ
pit
小于σr·kt
，满足公式(13)的判定标准，因此钢轨处于安全服役状态。通过对比可知，轨底伤损坑附近存在的较深的踏面低塌也会对伤损坑处受力产生显著影响，进而导致钢轨的服役安全性下降。
[0150]
本发明提出的钢轨服役安全性的计算方法，考虑了钢轨存在焊接接头低塌及附近轨底损伤坑(锈蚀坑)两种伤损的叠加影响，而非某一种伤损的单一影响。通过建立的理论
公式模型，评估了两种伤损叠加作用下的轨底锈蚀坑受力状态及钢轨服役安全性，能够有效解决铁路现场养护人员对此种状况钢轨服役安全状态的关切。
[0151]
本发明提出的钢轨服役安全性的计算方法，对于轨底锈蚀坑受力状态及钢轨服役安全性的计算，相比于可以采用的有限元仿真方法，本发明中的计算方法可以简单便捷的实现，不需要漫长的建模及运算时间及较高的有限元仿真知识技能门槛，因此更适用于现场铁路养护人员进行快速评估。
[0152]
本发明提出的钢轨服役安全性的计算方法，通过计算不同数值的踏面低塌深度及轨底锈蚀坑应力集中系数，分析随着两处伤损的发展，钢轨锈蚀坑受力及疲劳极限的演变。并通过锈蚀坑受力与疲劳极限的对比，能够对现实中的类似钢轨失效案例的失效原因分析提供可靠的数据支撑。
[0153]
针对上述各实施方式的详细解释，其目的仅在于对本发明进行解释，以便于能够更好地理解本发明，但是，这些描述不能以任何理由解释成是对本发明的限制，特别是，在不同的实施方式中描述的各个特征也可以相互任意组合，从而组成其他实施方式，除了有明确相反的描述，这些特征应被理解为能够应用于任何一个实施方式中，而并不仅局限于所描述的实施方式。

技术特征：
1.一种钢轨服役安全性的计算方法，其特征在于，所述计算方法包括：获取车辆载荷作用下钢轨伤损坑处的轨底最大动弯应力、踏面低塌产生的附加冲击作用在伤损坑处的附加动弯应力、以及轨底所受温度应力；基于所述轨底最大动弯应力、所述附加动弯应力和所述轨底温度应力，计算踏面低塌影响下钢轨伤损坑底部的最大拉应力；根据所述最大拉应力及相应损伤下的疲劳极限评价钢轨服役安全性。2.如权利要求1所述的钢轨服役安全性的计算方法，其特征在于，获取车辆载荷作用下钢轨伤损坑处的轨底最大动弯应力，包括；获取所述钢轨伤损坑处的钢轨截面静弯矩和所述钢轨伤损坑处的轨底端面系数，通过所述钢轨截面静弯矩和所述轨底端面系数计算所述轨底最大动弯应力，则有，式中，σ0为轨底最大动弯应力，m0为钢轨截面静弯矩，w1为轨底端面系数。3.如权利要求2所述的钢轨服役安全性的计算方法，其特征在于，获取所述钢轨伤损坑处的钢轨截面静弯矩，包括，确定车辆的最不利轮位、以及该最不利轮位的载荷；获取钢轨基础参数，以及所述钢轨伤损坑与所低塌的水平距离；根据所述最不利轮位的载荷、所述钢轨基础参数和所述水平距离，计算所述钢轨截面静弯矩。4.如权利要求3所述的钢轨服役安全性的计算方法，其特征在于，所述钢轨基础参数包括钢轨基础与钢轨刚比系数、钢轨基础弹性系数和钢轨竖向抗弯刚度，则有，式中，m0为钢轨截面静弯矩，∑p
μ
为各车轮的当量载荷，p为最不利轮位的载荷，μ＝e-βx
(cosβx-sinβx)，x为伤损坑与低塌的水平距离，β为钢轨基础与钢轨刚比系数，sinβx)，x为伤损坑与低塌的水平距离，β为钢轨基础与钢轨刚比系数，钢轨基础弹性系数，ej
x
为钢轨竖向抗弯刚度。5.如权利要求3所述的钢轨服役安全性的计算方法，其特征在于，所述最不利轮位，为所述车辆其前车的后转向架与后车的前转向架共四个轮位分别经过低塌处时，在所述伤损坑处产生的最大载荷的轮位。6.如权利要求1所述的钢轨服役安全性的计算方法，其特征在于，获取踏面低塌产生的附加冲击作用在轨底的附加动弯应力，包括，获取车辆的运行速度和所述踏面低塌的形状参数；基于所述踏面低塌的形状参数和所述车辆的运行速度，计算车辆经过所述踏面低塌时产生的附加冲击力；获取钢轨基础参数，以及所述钢轨伤损坑与所低塌的水平距离，并根据所述附加冲击力、所述钢轨基础参数以及所述水平距离系数计算所述附加冲击作用在伤损坑的附加动弯应力。
7.如权利要求5所述的钢轨服役安全性的计算方法，其特征在于，所述踏面低塌的形状参数包括最大低塌量、低塌长度和激励频率，则有，式中：λ为最大低塌量；l为低塌长度，ω为激励频率，且ω＝2πv/l，v为行车速度，t为从起点算起的时间，t＝l/v；l为低塌距起点的水平距离，η为车轮经过低塌的某一时刻t时对应的轨面不平顺值。8.如权利要求7所述的钢轨服役安全性的计算方法，其特征在于，基于所述踏面低塌的形状参数和所述车辆的运行速度，计算所述踏面低塌产生的附加冲击力，则有，式中：y为钢轨动挠度，表示车轮经过低塌的时间，表示簧下质量在轨道上自由振动的周期，为簧下质量的固有圆频率，k为轨下基础换算刚度，m为簧下质量，t为从起点算起的时间；因其中为钢轨基础弹性系数，l
k
为弹簧的有效换算长度，从而作用于钢轨上的附加动力p
c
为：p
c
＝ky＝8ej
x
β3y
ꢀꢀꢀꢀ
(6)式中，k为轨下基础换算刚度，y为钢轨动挠度，ej
x
为钢轨竖向抗弯刚度，β为钢轨基础与钢轨刚比系数，且钢轨刚比系数，且为钢轨基础弹性系数。9.如权利要求8所述的钢轨服役安全性的计算方法，其特征在于，根据所述附加冲击力、所述钢轨基础参数以及所述水平距离系数计算所述附加冲击作用在伤损坑的附加动弯应力，则有，式中，σ
cd
为附加动弯应力，p
c
为附加动力，μ＝e-βx
(cosβx-sinβx)，x为伤损坑与低塌的水平距离，β为钢轨基础与钢轨刚比系数，水平距离，β为钢轨基础与钢轨刚比系数，钢轨基础弹性系数，ej
x
为钢轨竖向抗弯刚度，且β、和ej
x
属于钢轨基础参数，w1为轨底端面系数。10.如权利要求1所述的钢轨服役安全性的计算方法，其特征在于，获取轨底所受温度应力，包括：获取钢轨的锁定轨温和测量轨温，根据所述锁定轨温和所述测量轨温计算所述轨底温度应力，则有，σ
t
＝2.45δt
ꢀꢀꢀꢀ
(8)式中，σ
t
为轨底所受的温度应力，δt为所述锁定轨温与所述测量轨温的差值。11.如权利要求1所述的钢轨服役安全性的计算方法，其特征在于，计算钢轨伤损坑底部的最大拉应力，包括：基于所述轨底最大动弯应力、所述附加动弯应力和所述轨底温度应力计算踏面低塌影
响下钢轨锈蚀处的轨底最大应力；获取轨底伤损坑的形状参数，并根据所述形状参数计算所述伤损坑的应力集中系数；根据所述钢轨锈蚀处的轨底最大应力和所述应力集中系数计算踏面低塌影响下所述钢轨伤损坑底部的最大拉应力。12.如权利要求11所述的钢轨服役安全性的计算方法，其特征在于，基于所述轨底最大动弯应力、所述附加动弯应力和所述轨底温度应力计算踏面低塌影响下钢轨锈蚀处的轨底最大应力，则有，σ
d
＝σ0+σ
cd
+σ
t
ꢀꢀꢀꢀ
(1)式中，σ
d
为轨底的最大应力，σ0为轨底最大动弯应力，σ
cd
为附加动弯应力，σ
t
为轨底所受的温度应力。13.如权利要求11所述的钢轨服役安全性的计算方法，其特征在于，所述轨底伤损坑的形状参数包括伤损坑的深度和伤损坑根部半径，则有，式中，k
t
为应力集中系数，h为伤损坑的深度；ρ为伤损坑根部半径。14.如权利要求13所述的钢轨服役安全性的计算方法，其特征在于，根据所述钢轨锈蚀处的轨底最大应力和所述应力集中系数计算踏面低塌影响下所述钢轨伤损坑底部的最大拉应力，则有，σ
pit
＝k
t
·
σ
d
ꢀꢀꢀꢀ
(10)式中，σ
pit
为伤损坑底部的最大拉应力，k
t
为应力集中系数，σ
d
为伤损坑处轨底的最大应力。15.如权利要求1所述的钢轨服役安全性的计算方法，其特征在于，根据所述最大拉应力及相应损伤下的疲劳极限评价钢轨服役安全性，包括，确认所述钢轨伤损坑底部的最大拉应力与钢轨在相应损伤下的疲劳极限是否符合判定标准公式，若符合所述判定标准公式，则判断钢轨处于安全服役状态；若不符合所述判定标准公式，则判断钢轨处于不安全服役状态。16.如权利要求15所述的钢轨服役安全性的计算方法，其特征在于，所述判定标准公式，为：σ
pit
≤σ
r
·
k
t
ꢀꢀꢀꢀ
(13)式中，σ
pit
为伤损坑底部的最大拉应力，k
t
为应力集中系数，σ
r
为相应伤损下的疲劳极限。17.如权利要求16所述的钢轨服役安全性的计算方法，其特征在于，所述相应损伤下的疲劳极限为包含踏面与轨底伤损坑影响的钢轨疲劳极限，则有，σ
r
＝aσ
y
+be+c
ꢀꢀꢀꢀ
(11)且，且，c＝12.0k
t
r+28.2k
t-82.8r-191，式中：σ
y
为钢轨的屈服强度；e为钢轨弹性模量；r为伤损坑处循环拉应力的应力比。18.如权利要求17所述的钢轨服役安全性的计算方法，其特征在于，伤损坑处循环拉应
力的应力比r与踏面低塌相关，则有，式中，σ
max
为车轮经过踏面低塌的时刻，伤损坑处的最大受力，σ
min
为轨底仅受温度应力而无低塌及车辆自重作用下的伤损坑处的受力，σ
t
为轨底所受的温度应力，σ
d
为轨底的最大应力。19.一种钢轨服役安全性的计算装置，其特征在于，所述装置包括：轨底最大动弯应力计算模块，通过所述钢轨截面静弯矩和所述轨底端面系数计算所述轨底最大动弯应力；附加动弯应力计算模块，根据所述附加冲击力、所述钢轨基础参数以及所述水平距离系数计算所述附加冲击作用在伤损坑的附加动弯应力；轨底所受温度应力，根据所述锁定轨温和所述测量轨温计算所述轨底温度应力；最大拉应力计算模块，基于所述轨底最大动弯应力、所述附加动弯应力和所述轨底温度应力，计算踏面低塌影响下钢轨伤损坑底部的最大拉应力；安全性评价模块，根据所述最大拉应力及相应损伤下的疲劳极限评价钢轨服役安全性。20.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至18任一所述方法。21.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1至18任一所述方法。

技术总结
本发明提出一种钢轨服役安全性的计算方法及装置，涉及铁路安全运行技术领域，该计算方法包括：获取车辆载荷作用下钢轨伤损坑处的轨底最大动弯应力、踏面低塌产生的附加冲击作用在伤损坑处的附加动弯应力、以及钢轨伤损坑处的轨底所受温度应力；基于轨底最大动弯应力、附加动弯应力和轨底温度应力，计算踏面低塌影响下钢轨伤损坑底部的最大拉应力；根据最大拉应力及相应损伤下的疲劳极限评价钢轨服役安全性。本发明提出的钢轨服役安全性的计算方法及装置能够计算和评价钢轨同时存在踏面低塌和附近轨底伤损坑两种伤损情况下的钢轨服役安全性。服役安全性。服役安全性。


技术研发人员：徐超 胡杰 杨其全 许鑫 王晨阳 张倩 吕晶 王玉婷 曹欣旺 孙旭文 冯博宇
受保护的技术使用者：中国铁道科学研究院集团有限公司金属及化学研究所 北京中铁科新材料技术有限公司 铁科金化科技有限公司 中国铁道科学研究院集团有限公司
技术研发日：2022.11.22
技术公布日：2023/4/5