一种精确预测钢轨轨温的方法

1.本发明属于轨道交通技术领域，具体涉及一种精确预测钢轨轨温的方法。


背景技术：

2.随着我国科学技术和生产力不断发展，高速铁路、重载铁路技术不断更新，列车通过钢轨接头时会产生很大的轮轨冲击力，对轨道结构产生很大的破坏作用，造成轨道部件破损，同时加剧了机车车辆的振动、车辆部件的破损，增加能耗和降低旅客舒适度。为了改善钢轨接头的工作状态，无缝线路问世，为大量减少钢轨接头创造了条件，因而具有行车平稳，机车车辆及轨道维修费用降低，设备使用寿命延长，适合高速行车等优点，是铁路轨道现代化的一项重要技术措施，也是当前高速、重载铁路必需的条件。但是，无缝线路钢轨在变温环境下不能沿纵向伸缩，因此造成钢轨内部温度应力增大，存在轨道结构高温胀轨和低温拉断的危险。所以定期对无缝线路纵向温度力进行检测具有尤为重要的意义。
3.目前的无缝线路温度检测主要分为有损检测和半有损检测方法，这两种方法会对原有钢轨的稳定性造成破坏，并且增加了检测过程风险。第一类半有损检测方法只能检测钢轨内部存在的拉应力，并未得到广泛应用。第二类半有损检测方法能保证轨道结构完整和稳定，但精度较低，耗费人力物力较大。近年来，一些无损检测法也应用到无缝线路上，如超声导波法、x射线法等，但此类方法并不能反映整个钢轨截面上的温度应力分布情况，同时受钢轨自身状态影响，其检测准确性有待深究。所以，亟需研究一种更为精确的无损钢轨温度检测方法。


技术实现要素：

4.本发明的目的在于解决现有方法对无缝线路钢轨轨温预测不佳的现状，提一种预测钢轨轨温的方法，其特征在于：所述方法包括如下步骤：
5.利用激励源在安装加速度传感器附近钢轨轨顶处施加振动，并获取钢轨振动加速度信号；
6.计算振动的频响函数曲线，得出不同纵向力f和自振频率之间的关系，并计算得到纵向力f；
7.将纵向力f代入f＝kδt*a，计算得到温度变化量δt，其中k为常数，a为钢轨横截面积；
8.根据锁定轨温与温度变化量δt得到钢轨实际轨温。
9.在一个实施方案中，所述方法还包括在钢轨轨顶上安装若干振动加速度传感器，每两组相邻轨枕等间距安装一个加速度传感器。
10.在一个实施方案中，所述不同纵向力和自振频率之间的关系为二次关系式。
11.在一个实施方案中，所述二次关系式为f＝af2+bf+c，其中，a为二次项系数，b为一次项系数，c为常数，f为自振频率。
12.在一个实施方案中，所述计算振动的频响函数曲线，包括，先计算激励f(t)的自相
关函数r
ff
(τ)，再计算f(t)与响应x(t)的互相关函数r
fx
(τ)，分别对r
ff
(τ)及r
fx
(τ)作傅里叶变换，得到f(t)的自功率谱s
ff
(f)及f(t)与x(t)的互功率谱s
fx
，频响函数满足
13.在一个实施方案中，每个加速度传感器通过数据线传输到一个或者多个采集仪，用于进行集中数据分析。
14.本发明的优点在于：
15.(1)针对目前方法解决现有方法对无缝线路钢轨温度预测不佳的现状，提出一种精确预测钢轨轨温的方法，可以简便无损的检测钢轨纵向力。
16.(2)目前已有的文献检测方法是假设钢轨纵向力和自振频率为线性关系，本发明基于试验和理论推导，提出一种纵向力
‑
自振频率二次关系拟合式，计算所得钢轨轨温结果更加符合试验值，提高了检测精度。
附图说明
17.图1是预测钢轨轨温方法流程图；
18.图2是轨道结构简化示意图；
19.图3是钢轨纵向力和自振频率的试验关系结果图；
20.图4是钢轨自振频率的平方与纵向力关系拟合图；
21.图5是钢轨纵向力和自振频率的二次关系拟合曲线图。
具体实施方式
22.为了使本发明的技术方案和优点更加清楚，下面结合附图和具体实施例对本发明进行详细描述。
23.本发明预测钢轨温度的方法包括：
24.利用激励源4在安装加速度传感器3附近钢轨1轨顶处施加振动，并获取钢轨1振动加速度信号；
25.计算振动的频响函数曲线，得出不同纵向力f和自振频率之间的关系，并计算得到纵向力f；
26.将纵向力f代入f＝kδt*a，计算得到温度变化量δt，其中k为常数，a为钢轨横截面积；
27.根据锁定轨温与温度变化量δt得到钢轨实际轨温。
28.其中，无缝线路锁定轨温是指无缝线路的零应力轨温，是指无缝线路钢轨被完全锁定时的轨温，此时钢轨内部的温度力为零。温度变化量δt有正负，正表示温度升高，负表示温度降低。
29.在一个方案中，所述方法还包括在钢轨1轨顶上安装若干振动加速度传感器3，每两组相邻轨枕等间距安装一个加速度传感器3，优选至少布置5个传感器。
30.在一个方案中，所述激励源为力锤4，通过在安装加速度传感器3附近钢轨1轨顶处进行敲击，使钢轨1产生振动加速度信号。
31.采用本发明的方案，通过振动激励可以测得温度变化量，进而可以测得钢轨的实
际温度，其为完全的无损的探测方式，相比于现有技术更加简单、便利且节约成本。而且本发明的方案非常巧妙的利用了温度与激振力的关系，现有技术中都是根据温度变化来计算激振力，本发明反其道而行之。
32.在一个方案中，所述不同纵向力和自振频率之间的关系为二次关系式。
33.在一个方案中，所述二次关系式为f＝af2+bf+c，其中，a为二次项系数，b为一次项系数，c为常数，f为自振频率。
34.经过研究和实验返现，激振力与自振频率的二次关系式更加精确，所以采用本发明的方案，能够更加准确计算激振力，进而可以更加精确计算温度。
35.在一个方案中，每个加速度传感器通过数据线传输到一个或者多个采集仪，用于进行集中数据分析。
36.如图2，示出了实施本发明方法的装置结构图，其中包括钢轨1、扣件2、振动加速度传感器3、力锤4。在钢轨1轨顶上安装若干振动加速度传感器3，每两组相邻轨枕等间距安装一个传感器3，用于测试钢轨1的振动加速度信号。
37.在一个方案中，计算振动的频响函数曲线，包括，先计算激励f(t)的自相关函数r
ff
(τ)，再计算f(t)与响应x(t)的互相关函数r
fx
(τ)，分别对r
ff
(τ)及r
fx
(τ)作傅里叶变换，得到f(t)的自功率谱s
ff
(f)及f(t)与x(t)的互功率谱s
fx
，频响函数满足
38.在一个方案中，钢轨单位长度质量为60kg/m时，经过计算，所述拟合曲线中二次项系数a和一次项系数b和常数项c的值按下表选取：
[0039][0040]
实施例1:
[0041]
下面结合附图说明书本发明预测钢轨轨温的方法的一个具体实施1，所述钢轨单位长度质量为60kg/m，轨枕间距为0.57m；所述系统包括钢轨1、扣件2、振动加速度传感器3、力锤4。在钢轨1轨顶上安装若干振动加速度传感器3，每两组相邻轨枕等间距安装一个传感器3。每个加速度传感器3通过数据线传输到一个或者多个采集仪(图中未实)，进行集中数据分析。
[0042]
随后计算激励f(t)的自相关函数r
ff
(τ)，再计算f(t)与响应x(t)的互相关函数r
fx
(τ)，分别对r
ff
(τ)及r
fx
(τ)作傅里叶变换，得到f(t)的自功率谱s
ff
(f)及f(t)与x(t)的互功
率谱s
fx
，频响函数满足
[0043]
通过实验分析计算得到不同钢轨纵向力和自振频率的关系，频率随纵向力的变化呈现出一个曲线变化的趋势，对其按照f＝af2+bf+c进行二次关系拟合，f表示钢轨弯曲振动第一阶pinned
‑
pinned共振频率，f表示钢轨纵向力。
[0044]
钢轨单位长度质量为60kg/m时，所述拟合曲线中二次项系数a和一次项系数b和常数项c的值按下表选取：
[0045][0046][0047]
试验轨枕间距为0.57m，根据上表查询，确定其二次关系拟合式为f＝
‑
0.00015196x2
‑
0.01798x+174.94603。
[0048]
通过计算进一步得到钢轨弯曲振动第一阶pinned
‑
pinned共振频率实际值f＝1013hz，将其代入拟合曲线f＝
‑
0.00015196x2
‑
0.01798x+174.94603，得到纵向力为0.79569mn。
[0049]
将纵向力的值代入f＝kδt*a，从而获得钢轨温度变化量δt＝41.43℃，即钢轨的温度变化量超过锁定轨温41.43℃，钢轨内部出现压应力。
[0050]
根据现场线路条件查得锁定轨温为22℃，因此可得钢轨实际温度为22℃+41.43℃＝63.43℃。
[0051]
下面通过公式推导验证钢轨自振频率和纵向力的二次关系实验结果的合理性，
[0052]
一段钢轨可以认为是一段梁，承受轴向力的钢轨的自由运动方程为
[0053][0054]
进行变量分离，得到
[0055][0056]
将式(2)展开为两个独立的常微分方程，即
[0057][0058]
eiφ
″″
(x)+nφ
″
(x)
‑
mω2φ(x)＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0059]
将式(4)处理为
[0060]
φ
″″
(x)+g2φ
″
(x)
‑
a4φ(x)＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0061]
其中：m表示质量，u表示垂向位移。n表示钢轨纵向力，x表示x方向一维梁，e表示弹性模量，i表示惯性矩，w表示圆频率，φ(x)为人为所定义函数。a和g为自定义参数
[0062]
式中各个参数满足
[0063][0064]
根据常微分方程通解形式，可以写出方程(4)的通解
[0065]
φ(x)＝asinδx+bcosδx+csinhεx+dcoshεx
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0066]
式中，a～d表示四个常数，由梁边界条件确定，且各项参数满足
[0067][0068]
根据简支梁两端铰接的边界条件，可以解得考虑轴向力时简支梁得自振频率为
[0069][0070]
由式(9)可以看出，当纵向力为正时，钢轨自振频率会减小，相当于降低了钢轨的刚度，且压力越大，频率降低得越多；当纵向力为负时，钢轨的自振频率会增大，相当于提高了梁的刚度，从中可以得出钢轨的纵向力和自振频率并非线性关系，而更靠近二次关系。同时，从图4可以看出，钢轨的自振频率的平方与纵向力呈线性关系，线性拟合相关系数高达0.99983(如下表)，也就是说，实质上钢轨的自振频率和力呈二次关系，证明了本发明方法的可靠性。
[0071][0072]
相关系数用公式进行计算，r越接近1，表示拟合越好，两者相关性越强，将上式中数据代入公式
[0073][0074]
可得钢轨的自振频率的平方与纵向力线性相关系数为0.99983
[0075]
上文所列出的一系列的详细说明仅仅是针对本发明的优选实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。