CN102180187A - 一种铁路轨道高低高精度检测装置和检测方法 - Google Patents
一种铁路轨道高低高精度检测装置和检测方法 Download PDFInfo
- Publication number
- CN102180187A CN102180187A CN2011100898054A CN201110089805A CN102180187A CN 102180187 A CN102180187 A CN 102180187A CN 2011100898054 A CN2011100898054 A CN 2011100898054A CN 201110089805 A CN201110089805 A CN 201110089805A CN 102180187 A CN102180187 A CN 102180187A
- Authority
- CN
- China
- Prior art keywords
- data
- shifter
- group
- centerdot
- time domain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Machines For Laying And Maintaining Railways (AREA)
Abstract
本发明属于铁路轨道检测技术,涉及对铁路轨道高低连续检测装置和方法的改进。本发明检测装置包括一个带有检测平台[7]的轨道检测车,在检测平台[7]的上表面安装有包括计算机[1]、位移器数据采集卡[2]、OD数据采集卡[6]和电源的数据处理***,其特征在于:所说的第一位移器PDT1和第二位移器PDT2是激光位移器,有一个由三个光纤陀螺以及FOG数据采集卡[3]组成的检测平台高低测量***。本发明检测方法的步骤是:陀螺数据预处理;里程数据预处理;位移器数据预处理;时域空域转换;检测平台[7]高低计算;轨道高低计算。本发明不受检测速度影响,能对左、右轨的高低进行分别测量,测量精度高,适用范围广。
Description
技术领域
本发明属于铁路轨道检测技术,涉及对铁路轨道高低连续检测装置和方法的改进。
背景技术
目前,对铁路轨道高低的检测方法主要有弦测法、惯性基准法、激光准直法、三维测量以及视像法等。
1、弦测法是指采用人工拉线的方法在轨道上真实搭建一条弦线,通过测量轨道轨顶面与该弦线的相对位移,评价轨道高低平顺性。这种方法原理简明,但是这种方法的传递函数不恒为1,当需要将测量结果转换成其他弦长的不平顺值时,存在误差;另外。测量精度与工作人员的熟练程度密切相关,且工作效率偏低,不能适应密集的铁路建设、维修。
2、惯性基准法是指通过加速度计二次积分测量得到轨道基准值(载体上下振动频率大大高于自身频率时,质量块不能追随而保持静止的位置。这个静止位置即为惯性基准),与位移器测量到的瞬时载体与基准值的相对位移组合得到轨道平顺测量值。这种方法能够比较如实地反映0.1~50m波长范围内的轨道不平顺。当检测速度过低时,加速度信号比较微弱,信噪比低,并且要作积分运算,低频信号容易引起积分饱和,因此惯性基准法不适合做低速测量。
3、激光准值法是指在两个测量基点之间通过光学视准线建立一条基准线,测量各测点偏离该基准线的偏离值。这种方法对环境要求较为苛刻,不适合强光、雨雪、沙尘等天气情况。
4、三维测量法借助于精密全站仪实现,光学跟踪测量法,指采用全站仪辅助测量的轨道平顺性检测方法。这种方法的基本思路是通过全站仪给出轨道测量点的精确三维坐标,将测量得到的三维坐标与设计值比较,得到轨道平顺性检测结果。这种测量方法的代表产品有瑞士“安伯格”和德国“GEDO”,是目前高铁建设、维护过程中普遍采用的方法,具备测量精度高的优点,但存在以下几点不足:
(1)、由于光在不同介质中的传播速度不同,对空气能见度有要求;
(2)受全站仪跟踪距离限制,每隔固定距离需要重新建站,工作效率较低;
(3)测量结果为离散点,不能反映轨道状态全貌。
5、光学摄像法是大型轨检车T17、GJ-5以及GJ-4改中采用测量方法。光学摄像传感器安装在检测车的检测梁上,在检测过程中,采用4台CCD摄像机对左右两轨高速拍照,之后通过Hub回传到车体内部的PC机中,通过图像重构的方法重现轨道几何外形,分析轨道的高低、轨向、水平等各项参数。该方法具备检测精度高、不受速度影响等优点,是目前最为先进的检测理论。但是它极易受外界光干扰,一旦受干扰整个***均不能工作。
发明内容
本发明的目的是:提出一种不受检测速度影响、能对左、右轨的高低进行分别测量、测量精度高、适用范围广的铁路轨道高低连续检测装置和方法。
本发明的技术方案是:一种铁路轨道高低高精度检测装置,包括一个带有检测平台的轨道检测车,在检测平台的下表面固定有两个位移器,它们是第一位移器PDT1和第二位移器PDT2,第一位移器PDT1和第二位移器PDT2分别位于轮轴两端的正上方,第一位移器PDT1和第二位移器PDT2的轴线与轮轴的轴线正交,在轮轴的一端安装有里程仪OD,在检测平台的上表面安装有数据处理***,该数据处理***包括计算机、位移器数据采集卡、OD数据采集卡和电源,位移器数据采集卡通过USB通讯线与计算机连接,里程仪OD的输出端与OD数据采集卡的输入端连接,OD数据采集卡通过RS422串口通讯线与计算机的第一串口COM1连接;其特征在于:
(1)所说的第一位移器PDT1和第二位移器PDT2是激光位移器,在数据处理***中有两个激光位移器控制盒,它们是PDT1控制盒和PDT2控制盒,第一位移器PDT1的输出端与PDT1控制盒的输入端连接,PDT1控制盒的输出端与位移器数据采集卡的第一输入端连接,第二位移器PDT2的输出端与PDT2控制盒的输入端连接,PDT2控制盒的输出端与位移器数据采集卡的第二输入端连接;
(2)有一个由第一光纤陀螺FOG1、第二光纤陀螺FOG2和第三光纤陀螺FOG3三个光纤陀螺以及FOG数据采集卡组成的检测平台高低测量***,三个光纤陀螺均安装在检测平台的上表面,第一光纤陀螺FOG1与第一位移器PDT1同轴,第三光纤陀螺FOG3与第二位移器PDT2同轴,第二光纤陀螺FOG2位于第一光纤陀螺FOG1和第三光纤陀螺FOG 3的中心,三个光纤陀螺的轴线相互平行并共面;第一光纤陀螺FOG1、第二光纤陀螺FOG2和第三光纤陀螺FOG 3的输出端分别与FOG数据采集卡的第一至第三输入端连接,FOG数据采集卡通过RS422串口通讯线与计算机的第二串口COM2连接。
使用如上面所述的检测装置检测铁路轨道高低的方法,其特征在于,检测的步骤如下:
1、陀螺数据预处理:
1.1、获取陀螺数据初值:保持检测平台静止m秒,m取值范围为100s至200s,计算机每5ms读取一次第二串口COM2,采集通过FOG数据采集卡传送来的第一光纤陀螺FOG1至第三光纤陀螺FOG 3的静止陀螺数据 分别计算m秒内的均值作为陀螺数据初值,记为ω10、ω20、ω30;
1.2、陀螺数据数据采集:启动轨道检测车在轨道上行驶,计算机每5ms读取一次第二串口COM2,采集通过FOG数据采集卡传送来的第一光纤陀螺FOG1至第三光纤陀螺FOG 3的实时陀螺数据ω1、ω2、ω3;
1.3、时域滤波:
将采集到的第一光纤陀螺FOG1至第三光纤陀螺FOG3的实时陀螺数据ω1、ω2、ω3分别做n点滑动平均,n取值范围为50~150的整数,得到时域滤波后的陀螺数据ω11、ω21、ω31;滑动平均的计算公式为:
式中,当k<n时,有ω11(k)=ω1(k),ω21(k)=ω2(k),α31(k)=ω3(k),k为自然数,是采样的序列号;
1.4、陀螺数据补偿:对经过n点滑动平均的陀螺数据ω11、ω21、ω31按照下式进行补偿,得到补偿后的陀螺数据ω12、ω22、ω32;
ω12(k)=ω11(k)-ω10-f(ψ(k-1))………………………………[4]
ω32(k)=ω31(k)-ω30-f(ψ(k-1))………………………………[5]
ω22(k)=ω21(k)-ω20…………………………………………… [6]
式中,ω12(1)=ω11(1),ω22(1)=ω21(1),ω32(1)=ω31(1),
f(ψ(k-1))为相对航向角函数,ψ(1)=0;
1.5、相对角度计算:对补偿后的陀螺数据ω12、ω22、ω32进行积分运算,得到左相对俯仰角θL、右相对俯仰角θR、相对航向角ψ,计算公式如下:
θL(k)=θL(k-1)+ω12·Ts1………………………………………[7]
θR(k)=θR(k-1)+ω32·Ts1………………………………………[8]
ψ(k)=ψ(k-1)+ω22·Ts1…………………………………………[9]
式中,采样间隔Ts1=5ms;
2、里程数据预处理:在轨道检测车运行过程中,计算机每5ms读取一次第一串口COM1,采集通过OD数据采集卡传送来的里程仪累积脉冲数Nod,根据计算机内事先标定好的里程仪刻度系数kod计算里程lod;
lod=Nod×kod………………………………………………………[10]
3、位移器数据预处理:
3.1、获取位移器数据初值:保持检测平台静止m秒,计算机每5ms读取一次USB通信端口,采集通过位移器数据采集卡传送来的第一位移器PDT1和第二位移器PDT2的测量值和分别计算m秒内和的平均值作为初值d10和d20;
3.2、位移器数据采集:在轨道检测车运行过程中,计算机每5ms读取一次USB通信端口,采集通过位移器数据采集卡传送来的第一位移器PDT1和第二位移器PDT2的测量值d1和d2;
3.3、位移器数据补偿:补偿后得到时域内位移器数据的预处理数值dL和dR,补偿公式为:
dL(k)=d1(k)-d10………………………………………………… [11]
dR(k)=dR(k)-dR0………………………………………………… [12]
4、时域空域转换:根据步骤1.5计算得到的相对角度θL、θR和步骤2计算得到的里程lod以及步骤3计算得到的位移器数据的预处理数值dL和dR进行时域空域转换,转换的步骤如下:
4.1、将根据步骤1.5计算得到的相对角度θL、θR和步骤2计算得到的里程lod以及步骤3计算得到的位移器数据的预处理数值dL和dR组成时域数据组[lod,θL,θR,dL,dR];将对应的空域数据记为[L、ΘL、ΘR、DL、DR],将K作为时域数据组的序列号,将j作为空域数组的序列号;
4.2、将第一组时域数据组K1作为第一组空域数据组j1;
4.3、读取下一组时域数据组K2,即k=2时的时域数据;
4.4、计算第二组时域数据K2中的时域里程lod(2)与第一组空域数据的空域里程L(1)之差dL;
4.5、判断下列条件是否成立:
|dL-Lmd|<0.05m……………………………………………[13]
式中,Lmd表示检测密度,Lmd为0.125m;若成立,则将第二组时域数据作为第二组空域数据;否则,按照以下方法进行计算:
4.6、读取下一组时域数据组K3,即k=3时的时域数据;
4.7、计算第三组时域数据K3中的时域里程lod(3)与第一组空域数据的空域里程L(1)之差dL;
4.8、判断判定条件是否成立,进行时域/空域数据转换:
|dL-Lmd|<0.05m……………………………………………[13]
若成立,则将第三组时域数据作为第二组空域数据;否则,继续读取下一组时域数据组并进行计算和判断,直到上述判定条件成立为止,设满足判定条件的时域数据为第Kη组时域数据,η为大于3的整数,则将作为第二组空域数据组j2;
4.9、读取时域数据组的下一组时域数据组,计算该组时域数据Kη+1中的时域里程lod(η+1)与第二组空域数据的空域里程L(2)之差dL;根据步骤4.8的方法进行时域/空域数据转换,直到检测结束为止;
5、检测平台高低计算:计算公式为:
式中,hVL为检测平台的左高低,hVR为检测平台的右高低,α为幅值补偿系数,取0.25;λ为检测弦长,为10m或30m;
dL1=L(j)-L(j-d);
dL2=L(j+d)-L(j-d);
dΘ1=ΘL(j)-ΘL(j-d);
dΘ2=ΘL(j+d)-ΘL(j-d);
dΘ3=ΘR(j)-ΘR(j-d);
dΘ4=ΘR(j+d)-ΘR(j-d);
6、轨道高低计算:计算公式为:
hL(j)=hVL(j)+DL(j)……………………………………………[16]
hR(j)=hVR(j)+DR(j)……………………………………………[17]
式中,hL为轨道左高低,hR为轨道右高低。
本发明的优点是:不受检测速度影响,能对左、右轨的高低进行分别测量,测量精度高,适用范围广。
附图说明
图1是本发明测量装置的结构原理框图。
图2是本发明测量装置中传感器安装位置示意图。
具体实施方式
下面对本发明做进一步详细说明。规定方位如下,面对轨道检测车的前进方向,左手边为左方。参见图1、2,一种铁路轨道高低高精度检测装置,包括一个带有检测平台7的轨道检测车,在检测平台7的下表面固定有两个位移器,它们是第一位移器PDT1和第二位移器PDT2,第一位移器PDT1和第二位移器PDT2分别位于轮轴8两端的正上方,第一位移器PDT1和第二位移器PDT2的轴线与轮轴8的轴线正交,在轮轴8的一端安装有里程仪OD,在检测平台7的上表面安装有数据处理***,该数据处理***包括计算机1、位移器数据采集卡2、OD数据采集卡6和电源,位移器数据采集卡2通过USB通讯线与计算机1连接,里程仪OD的输出端与OD数据采集卡6的输入端连接,OD数据采集卡6通过RS422串口通讯线与计算机1的第一串口COM1连接;其特征在于:
(1)所说的第一位移器PDT1和第二位移器PDT2是激光位移器,在数据处理***中有两个激光位移器控制盒,它们是PDT1控制盒4和PDT2控制盒5,第一位移器PDT1的输出端与PDT1控制盒4的输入端连接,PDT1控制盒4的输出端与位移器数据采集卡2的第一输入端连接,第二位移器PDT2的输出端与PDT2控制盒5的输入端连接,PDT2控制盒5的输出端与位移器数据采集卡2的第二输入端连接;
(2)有一个由第一光纤陀螺FOG1、第二光纤陀螺FOG2和第三光纤陀螺FOG 3三个光纤陀螺以及FOG数据采集卡3组成的检测平台高低测量***,三个光纤陀螺均安装在检测平台7的上表面,第一光纤陀螺FOG1与第一位移器PDT1同轴,第三光纤陀螺FOG3与第二位移器PDT2同轴,第二光纤陀螺FOG2位于第一光纤陀螺FOG1和第三光纤陀螺FOG 3的中心,三个光纤陀螺的轴线相互平行并共面;第一光纤陀螺FOG1、第二光纤陀螺FOG2和第三光纤陀螺FOG3的输出端分别与FOG数据采集卡3的第一至第三输入端连接,FOG数据采集卡3通过RS422串口通讯线与计算机1的第二串口COM2连接。
本发明检测装置的检测原理是:利用陀螺分别测量轨道车的左右端相对俯仰角,结合里程和检测波长计算得到车体的高低,结合第一激光位移器和第二激光位移器的测量结果,扣除车体相对于轮轴的振动量,计算得到左右股轨道的高低。
使用如上面所述的检测装置检测铁路轨道高低的方法,其特征在于,检测的步骤如下:
1、陀螺数据预处理:
1.1、获取陀螺数据初值:保持检测平台7静止m秒,m取值范围为100s至200s,计算机1每5ms读取一次第二串口COM2,采集通过FOG数据采集卡3传送来的第一光纤陀螺FOG1至第三光纤陀螺FOG 3的静止陀螺数据分别计算m秒内的均值作为陀螺数据初值,记为ω10、ω20、ω30;
1.2、陀螺数据数据采集:启动轨道检测车在轨道上行驶,计算机1每5ms读取一次第二串口COM2,采集通过FOG数据采集卡3传送来的第一光纤陀螺FOG1至第三光纤陀螺FOG3的实时陀螺数据ω1、ω2、ω3;
1.3、时域滤波:
将采集到的第一光纤陀螺FOG1至第三光纤陀螺FOG3的实时陀螺数据ω1、ω2、ω3分别做n点滑动平均,n取值范围为50~150的整数,得到时域滤波后的陀螺数据ω11、ω21、ω31;滑动平均的计算公式为:
式中,当k<n时,有ω11(k)=ω1(k),ω21(k)=ω2(k),ω31(k)=ω3(k),k为自然数,是采样的序列号;
1.4、陀螺数据补偿:对经过n点滑动平均的陀螺数据ω11、ω21、ω31按照下式进行补偿,得到补偿后的陀螺数据ω12、ω22、ω32;
ω12(k)=ω11(k)-ω10-f(ψ(k-1))………………………………[4]
ω32(k)=ω31(k)-ω30-f(ψ(k-1))………………………………[5]
ω22(k)=ω21(k)-ω20…………………………………………… [6]
式中,ω12(1)=ω11(1),ω22(1)=ω21(1),ω32(1)=ω31(1),
f(ψ(k-1))为相对航向角函数,ψ(1)=0;
1.5、相对角度计算:对补偿后的陀螺数据ω12、ω22、ω32进行积分运算,得到左相对俯仰角θL、右相对俯仰角θR、相对航向角ψ,计算公式如下:
θL(k)=θL(k-1)+ω12·Ts1………………………………………[7]
θR(k)=θR(k-1)+ω32·Ts1………………………………………[8]
ψ(k)=ψ(k-1)+ω22·Ts1…………………………………………[9]
式中,采样间隔Ts1=5ms;
2、里程数据预处理:在轨道检测车运行过程中,计算机1每5ms读取一次第一串口COM1,采集通过OD数据采集卡6传送来的里程仪累积脉冲数Nod,根据计算机内事先标定好的里程仪刻度系数kod计算里程lod;
lod=Nod×kod………………………………………………………[10]
3、位移器数据预处理:
3.1、获取位移器数据初值:保持检测平台7静止m秒,计算机1每5ms读取一次USB通信端口,采集通过位移器数据采集卡2传送来的第一位移器PDT1和第二位移器PDT2的测量值和分别计算m秒内和的平均值作为初值d10和d20;
3.2、位移器数据采集:在轨道检测车运行过程中,计算机1每5ms读取一次USB通信端口,采集通过位移器数据采集卡2传送来的第一位移器PDT1和第二位移器PDT2的测量值d1和d2;
3.3、位移器数据补偿:补偿后得到时域内位移器数据的预处理数值dL和dR,补偿公式为:
dL(k)=d1(k)-d10…………………………………………………[11]
dR(k)=dR(k)-dR0…………………………………………………[12]
4、时域空域转换:根据步骤1.5计算得到的相对角度θL、θR和步骤2计算得到的里程lod以及步骤3计算得到的位移器数据的预处理数值dL和dR进行时域空域转换,转换的步骤如下:
4.1、将根据步骤1.5计算得到的相对角度θL、θR和步骤2计算得到的里程lod以及步骤3计算得到的位移器数据的预处理数值dL和dR组成时域数据组[lod,θL,θR,dL,dR];将对应的空域数据记为[L、ΘL、ΘR、DL、DR],将K作为时域数据组的序列号,将j作为空域数组的序列号;
4.2、将第一组时域数据组K1作为第一组空域数据组j1;
4.3、读取下一组时域数据组K2,即k=2时的时域数据;
4.4、计算第二组时域数据K2中的时域里程lod(2)与第一组空域数据的空域里程L(1)之差dL;
4.5、判断下列条件是否成立:
|dL-Lmd|<0.05m……………………………………………[13]
式中,Lmd表示检测密度,Lmd为0.125m;若成立,则将第二组时域数据作为第二组空域数据;否则,按照以下方法进行计算:
4.6、读取下一组时域数据组K3,即k=3时的时域数据;
4.7、计算第三组时域数据K3中的时域里程lod(3)与第一组空域数据的空域里程L(1)之差dL;
4.8、判断判定条件是否成立,进行时域/空域数据转换:
|dL-Lmd|<0.05m……………………………………………[13]
若成立,则将第三组时域数据作为第二组空域数据;否则,继续读取下一组时域数据组并进行计算和判断,直到上述判定条件成立为止,设满足判定条件的时域数据为第Kη组时域数据,η为大于3的整数,则将作为第二组空域数据组j2;
4.9、读取时域数据组的下一组时域数据组,计算该组时域数据Kη+1中的时域里程lod(η+1)与第二组空域数据的空域里程L(2)之差dL;根据步骤4.8的方法进行时域/空域数据转换,直到检测结束为止;
5、检测平台7高低计算:计算公式为:
式中,hVL为检测平台7的左高低,hVR为检测平台7的右高低,α为幅值补偿系数,取0.25;λ为检测弦长,为10m或30m;
dL1=L(j)-L(j-d);
dL2=L(j+d)-L(j-d);
dΘ1=ΘL(j)-ΘL(j-d);
dΘ2=ΘL(j+d)-ΘL(j-d);
dΘ3=ΘR(j)-ΘR(j-d);
dΘ4=ΘR(j+d)-ΘR(j-d);
6、轨道高低计算:计算公式为:
hL(j)=hVL(j)+DL(j)……………………………………………[16]
hR(j)=hVR(j)+DR(j)……………………………………………[17]
式中,hL为轨道左高低,hR为轨道右高低。
实施例
在装有检测平台的轨道车上安装如权利要求书1中所述检测装置,其中关键传感器的技术指标分别为:FOG1、FOG2、FOG 3均选用0.05°/h的陀螺,PDT1和PDT2均选用测量精度为0.05mm的LK-400G激光位移器,OD选用3600脉冲/转的光电编码器;对采集来的数据作如下处理:
1、陀螺数据预处理:
1.1、获取陀螺数据初值:保持检测平台[7]静止m秒,m取值范围为100s至200s,计算机[1]每5ms读取一次第二串口COM2,采集通过FOG数据采集卡[3]传送来的第一光纤陀螺FOG1至第三光纤陀螺FOG 3的静止陀螺数据分别计算m秒内的均值作为陀螺数据初值,记为ω10、ω20、ω30;
1.2、陀螺数据数据采集:启动轨道检测车在轨道上行驶,计算机[1]每5ms读取一次第二串口COM2,采集通过FOG数据采集卡[3]传送来的第一光纤陀螺FOG1至第三光纤陀螺FOG 3的实时陀螺数据ω1、ω2、ω3;
1.3、时域滤波:
将采集到的第一光纤陀螺FOG1至第三光纤陀螺FOG3的实时陀螺数据ω1、ω2、ω3分别做n点滑动平均,n取值范围为50~150的整数,得到时域滤波后的陀螺数据ω11、ω21、ω31;滑动平均的计算公式为:
式中,当k<n时,有ω11(k)=ω1(k),ω21(k)=ω2(k),ω31(k)=ω3(k),k为自然数,是采样的序列号;
1.4、陀螺数据补偿:对经过n点滑动平均的陀螺数据ω11、ω21、ω31按照下式进行补偿,得到补偿后的陀螺数据ω12、ω22、ω32;
ω12(k)=ω11(k)-ω10-f(ψ(k-1))……………………………………[4]
ω32(k)=ω31(k)-ω30-f(ψ(k-1))……………………………………[5]
ω22(k)=ω21(k)-ω20………………………………………………… [6]
式中,ω12(1)=ω11(1),ω22(1)=ω21(1),ω32(1)=ω31(1),
f(ψ(k-1))为相对航向角函数,ψ(1)=0;
1.5、相对角度计算:对补偿后的陀螺数据ω12、ω22、ω32进行积分运算,得到左相对俯仰角θL、右相对俯仰角θR、相对航向角ψ,计算公式如下:
θL(k)=θL(k-1)+ω12·Ts1……………………………………………[7]
θR(k)=θR(k-1)+ω32·Ts1……………………………………………[8]
ψ(k)=ψ(k-1)+ω22·Ts1………………………………………………[9]
式中,采样间隔Ts1=5ms;
2、里程数据预处理:在轨道检测车运行过程中,计算机[1]每5ms读取一次第一串口COM1,采集通过OD数据采集卡[6]传送来的里程仪累积脉冲数Nod,根据计算机内事先标定好的里程仪刻度系数kod计算里程lod;
lod=Nod×kod……………………………………………………………[10]
3、位移器数据预处理:
3.1、获取位移器数据初值:保持检测平台[7]静止m秒,计算机[1]每5ms读取一次USB通信端口,采集通过位移器数据采集卡[2]传送来的第一位移器PDT1和第二位移器PDT2的测量值和分别计算m秒内和的平均值作为初值d10和d20;
3.2、位移器数据采集:在轨道检测车运行过程中,计算机[1]每5ms读取一次USB通信端口,采集通过位移器数据采集卡[2]传送来的第一位移器PDT1和第二位移器PDT2的测量值d1和d2;
3.3、位移器数据补偿:补偿后得到时域内位移器数据的预处理数值dL和dR,补偿公式为:
dL(k)=d1(k)-d10…………………………………………………[11]
dR(k)=dR(k)-dR0…………………………………………………[12]
4、时域空域转换:根据步骤1.5计算得到的相对角度θL、θR和步骤2计算得到的里程lod以及步骤3计算得到的位移器数据的预处理数值dL和dR进行时域空域转换,转换的步骤如下:
4.1、将根据步骤1.5计算得到的相对角度θL、θR和步骤2计算得到的里程lod以及步骤2.3计算得到的位移器数据的预处理数值dL和dR组成时域数据组[lod,θL,θR,dL,dR];将对应的空域数据记为[L、ΘL、ΘR、DL、DR],将K作为时域数据组的序列号,将j作为空域数组的序列号;
4.2、将第一组时域数据组K1作为第一组空域数据组j1;
4.3、读取下一组时域数据组K2,即k=2时的时域数据;
4.4、计算第二组时域数据K2中的时域里程lod(2)与第一组空域数据的空域里程L(1)之差dL;
4.5、判断下列条件是否成立:
|dL-Lmd|<0.05m……………………………………………[13]
式中,Lmd表示检测密度,Lmd为0.125m;若成立,则将第二组时域数据作为第二组空域数据;否则,按照以下方法进行计算:
4.6、读取下一组时域数据组K3,即k=3时的时域数据;
4.7、计算第三组时域数据K3中的时域里程lod(3)与第一组空域数据的空域里程L(1)之差dL;
4.8、判断判定条件是否成立,进行时域/空域数据转换:
|dL-Lmd|<0.05m……………………………………………[13]
若成立,则将第三组时域数据作为第二组空域数据;否则,继续读取下一组时域数据组并进行计算和判断,直到上述判定条件成立为止,设满足判定条件的时域数据为第Kη组时域数据,η为大于3的整数,则将作为第二组空域数据组j2;
4.9、读取时域数据组的下一组时域数据组,计算该组时域数据Kη+1中的时域里程lod(η+1)与第二组空域数据的空域里程L(2)之差dL;根据步骤2.4.8的方法进行时域/空域数据转换,直到检测结束为止;
5、检测平台[7]高低计算:计算公式为:
式中,hVL为检测平台[7]的左高低,hVR为检测平台[7]的右高低,α为幅值补偿系数,取0.25;λ为检测弦长,为10m或30m;
dL1=L(j)-L(j-d);
dL2=L(j+d)-L(j-d);
dΘ1=ΘL(j)-ΘL(j-d);
dΘ2=ΘL(j+d)-ΘL(j-d);
dΘ3=ΘR(j)-ΘR(j-d);
dΘ4=ΘR(j+d)-ΘR(j-d);
6、轨道高低计算:计算公式为:
hL(j)=hVL(j)+DL(j)……………………………………………[16]
hR(j)=hVR(j)+DR(j)……………………………………………[17]
式中,hL为轨道左高低,hR为轨道右高低。
处理得到的结果,30m波长的轨道左高低测量精度为±1mm,右高低测量精度为±1mm。
Claims (2)
1.一种铁路轨道高低高精度检测装置,包括一个带有检测平台[7]的轨道检测车,在检测平台[7]的下表面固定有两个位移器,它们是第一位移器PDT1和第二位移器PDT2,第一位移器PDT1和第二位移器PDT2分别位于轮轴[8]两端的正上方,第一位移器PDT1和第二位移器PDT2的轴线与轮轴[8]的轴线正交,在轮轴[8]的一端安装有里程仪OD,在检测平台[7]的上表面安装有数据处理***,该数据处理***包括计算机[1]、位移器数据采集卡[2]、OD数据采集卡[6]和电源,位移器数据采集卡[2]通过USB通讯线与计算机[1]连接,里程仪OD的输出端与OD数据采集卡[6]的输入端连接,OD数据采集卡[6]通过RS422串口通讯线与计算机[1]的第一串口COM1连接;其特征在于:
(1)所说的第一位移器PDT1和第二位移器PDT2是激光位移器,在数据处理***中有两个激光位移器控制盒,它们是PDT1控制盒[4]和PDT2控制盒[5],第一位移器PDT1的输出端与PDT1控制盒[4]的输入端连接,PDT1控制盒[4]的输出端与位移器数据采集卡[2]的第一输入端连接,第二位移器PDT2的输出端与PDT2控制盒[5]的输入端连接,PDT2控制盒[5]的输出端与位移器数据采集卡[2]的第二输入端连接;
(2)有一个由第一光纤陀螺FOG1、第二光纤陀螺FOG2和第三光纤陀螺FOG 3三个光纤陀螺以及FOG数据采集卡[3]组成的检测平台高低测量***,三个光纤陀螺均安装在检测平台[7]的上表面,第一光纤陀螺FOG1与第一位移器PDT1同轴,第三光纤陀螺FOG3与第二位移器PDT2同轴,第二光纤陀螺FOG2位于第一光纤陀螺FOG1和第三光纤陀螺FOG3的中心,三个光纤陀螺的轴线相互平行并共面;第一光纤陀螺FOG1、第二光纤陀螺FOG2和第三光纤陀螺FOG3的输出端分别与FOG数据采集卡[3]的第一至第三输入端连接,FOG数据采集卡[3]通过RS422串口通讯线与计算机[1]的第二串口COM2连接。
2.使用如权利要求1所述的检测装置检测铁路轨道高低的方法,其特征在于,检测的步骤如下:
2.1、陀螺数据预处理:
2.1.1、获取陀螺数据初值:保持检测平台[7]静止m秒,m取值范围为100s至200s,计算机[1]每5ms读取一次第二串口COM2,采集通过FOG数据采集卡[3]传送来的第一光纤陀螺FOG1至第三光纤陀螺FOG3的静止陀螺数据分别计算m秒内的均值作为陀螺数据初值,记为ω10、ω20、ω30;
2.1.2、陀螺数据数据采集:启动轨道检测车在轨道上行驶,计算机[1]每5ms读取一次第二串口COM2,采集通过FOG数据采集卡[3]传送来的第一光纤陀螺FOG1至第三光纤陀螺FOG3的实时陀螺数据ω1、ω2、ω3;
2.1.3、时域滤波:
将采集到的第一光纤陀螺FOG1至第三光纤陀螺FOG 3的实时陀螺数据ω1、ω2、ω3分别做n点滑动平均,n取值范围为50~150的整数,得到时域滤波后的陀螺数据ω11、ω21、ω31;滑动平均的计算公式为:
式中,当k<n时,有ω11(k)=ω1(k),ω21(k)=ω2(k),ω31(k)=ω3(k),k为自然数,是采样的序列号;
2.1.4、陀螺数据补偿:对经过n点滑动平均的陀螺数据ω11、ω21、ω31按照下式进行补偿,得到补偿后的陀螺数据ω12、ω22、ω32;
ω12(k)=ω11(k)-ω10-f(ψ(k-1))……………………………………[4]
ω32(k)=ω31(k)-ω30-f(ψ(k-1))……………………………………[5]
ω22(k)=ω21(k)-ω20………………………………………………… [6]
式中,ω12(1)=ω11(1),ω22(1)=ω21(1),ω32(1)=ω31(1),
f(ψ(k-1))为相对航向角函数,ψ(1)=0;
2.1.5、相对角度计算:对补偿后的陀螺数据ω12、ω22、ω32进行积分运算,得到左相对俯仰角θL、右相对俯仰角θR、相对航向角ψ,计算公式如下:
θL(k)=θL(k-1)+ω12·Ts1……………………………………………[7]
θR(k)=θR(k-1)+ω32·Ts1……………………………………………[8]
ψ(k)=ψ(k-1)+ω22·Ts1………………………………………………[9]
式中,采样间隔Ts1=5ms;
2.2、里程数据预处理:在轨道检测车运行过程中,计算机[1]每5ms读取一次第一串口COM1,采集通过OD数据采集卡[6]传送来的里程仪累积脉冲数Nod,根据计算机内事先标定好的里程仪刻度系数kod计算里程lod;
lod=Nod×kod……………………………………………………………[10]
2.3、位移器数据预处理:
2.3.1、获取位移器数据初值:保持检测平台[7]静止m秒,计算机[1]每5ms读取一次USB通信端口,采集通过位移器数据采集卡[2]传送来的第一位移器PDT1和第二位移器PDT2的测量值和分别计算m秒内和的平均值作为初值d10和d20;
2.3.2、位移器数据采集:在轨道检测车运行过程中,计算机[1]每5ms读取一次USB通信端口,采集通过位移器数据采集卡[2]传送来的第一位移器PDT1和第二位移器PDT2的测量值d1和d2;
2.3.3、位移器数据补偿:补偿后得到时域内位移器数据的预处理数值dL和dR,补偿公式为:
dL(k)=d1(k)-d10…………………………………………………[11]
dR(k)=dR(k)-dR0…………………………………………………[12]
2.4、时域空域转换:根据步骤2.1.5计算得到的相对角度θL、θR和步骤2.2计算得到的里程lod以及步骤2.3计算得到的位移器数据的预处理数值dL和dR进行时域空域转换,转换的步骤如下:
2.4.1、将根据步骤2.1.5计算得到的相对角度θL、θR和步骤2.2计算得到的里程lod以及步骤2.3计算得到的位移器数据的预处理数值dL和dR组成时域数据组[lod,θL,θR,dL,dR];将对应的空域数据记为[L、ΘL、ΘR、DL、DR],将K作为时域数据组的序列号,将j作为空域数组的序列号;
2.4.2、将第一组时域数据组K1作为第一组空域数据组j1;
2.4.3、读取下一组时域数据组K2,即k=2时的时域数据;
2.4.4、计算第二组时域数据K2中的时域里程lod(2)与第一组空域数据的空域里程L(1)之差dL;
2.4.5、判断下列条件是否成立:
|dL-Lmd|<0.05m……………………………………………[13]
式中,Lmd表示检测密度,Lmd为0.125m;若成立,则将第二组时域数据作为第二组空域数据;否则,按照以下方法进行计算:
2.4.6、读取下一组时域数据组K3,即k=3时的时域数据;
2.4.7、计算第三组时域数据K3中的时域里程lod(3)与第一组空域数据的空域里程L(1)之差dL;
2.4.8、判断判定条件是否成立,进行时域/空域数据转换:
|dL-Lmd|<0.05m……………………………………………[13]
若成立,则将第三组时域数据作为第二组空域数据;否则,继续读取下一组时域数据组并进行计算和判断,直到上述判定条件成立为止,设满足判定条件的时域数据为第Kη组时域数据,η为大于3的整数,则将作为第二组空域数据组j2;
2.4.9、读取时域数据组的下一组时域数据组,计算该组时域数据Kη+1中的时域里程lod(η+1)与第二组空域数据的空域里程L(2)之差dL;根据步骤2.4.8的方法进行时域/空域数据转换,直到检测结束为止;
2.5、检测平台[7]高低计算:计算公式为:
式中,hVL为检测平台[7]的左高低,hVR为检测平台[7]的右高低,α为幅值补偿系数,取0.25;λ为检测弦长,为10m或30m;
dL1=L(j)-L(j-d);
dL2=L(j+d)-L(j-d);
dΘ1=ΘL(j)-ΘL(j-d);
dΘ2=ΘL(j+d)-ΘL(j-d);
dΘ3=ΘR(j)-ΘR(j-d);
dΘ4=ΘR(j+d)-ΘR(j-d);
2.6、轨道高低计算:计算公式为:
hL(j)=hVL(j)+DL(j)……………………………………………[16]
hR(j)=hVR(j)+DR(j)……………………………………………[17]
式中,hL为轨道左高低,hR为轨道右高低。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110089805 CN102180187B (zh) | 2011-04-12 | 2011-04-12 | 一种铁路轨道高低高精度检测装置和检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110089805 CN102180187B (zh) | 2011-04-12 | 2011-04-12 | 一种铁路轨道高低高精度检测装置和检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102180187A true CN102180187A (zh) | 2011-09-14 |
CN102180187B CN102180187B (zh) | 2013-04-10 |
Family
ID=44566525
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110089805 Active CN102180187B (zh) | 2011-04-12 | 2011-04-12 | 一种铁路轨道高低高精度检测装置和检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102180187B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103174072A (zh) * | 2011-12-26 | 2013-06-26 | 同济大学 | 基于三点偏弦法的轨面检测方法 |
CN103552571A (zh) * | 2013-11-18 | 2014-02-05 | 南车洛阳机车有限公司 | 综合检测车的轨道检测*** |
CN104260751A (zh) * | 2014-09-28 | 2015-01-07 | 北京理工大学 | 一种多传感器融合的高铁轨道中心线检测***和检测方法 |
CN104536462A (zh) * | 2015-01-09 | 2015-04-22 | 西安应用光学研究所 | 基于光纤陀螺积分手段的位置控制方法 |
CN104567783A (zh) * | 2014-12-01 | 2015-04-29 | 华南理工大学 | 使用mems陀螺仪的便携式轨道几何量检测***及方法 |
CN105424057A (zh) * | 2015-10-30 | 2016-03-23 | 东南大学 | 一种基于改进lmd的光纤陀螺振动信号分析方法 |
CN109477712A (zh) * | 2016-07-27 | 2019-03-15 | 中国涂料株式会社 | 三维表面粗糙度评价装置、三维表面粗糙度评价方法、三维表面粗糙度数据获取装置和三维表面粗糙度数据获取方法 |
CN109606417A (zh) * | 2018-12-04 | 2019-04-12 | 武汉理工大学 | 一种铁路轨道超高检测装置 |
CN113548068A (zh) * | 2021-07-23 | 2021-10-26 | 中车长春轨道客车股份有限公司 | 轨面不平顺检测装置及检测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4421741A1 (de) * | 1994-06-22 | 1996-01-04 | Ruhrkohle Ag | Vorrichtung zum Messen einer Höhendifferenz an Gleisanlagen |
JPH1030919A (ja) * | 1996-07-16 | 1998-02-03 | Pasuko Doro Gijutsu Center:Kk | 道路縦断プロファイル測定装置及び方法 |
CN2473089Y (zh) * | 2001-04-28 | 2002-01-23 | 上海理工大学附属二厂 | 轨道综合测量仪 |
CN201158385Y (zh) * | 2008-02-25 | 2008-12-03 | 长沙悦诚机电科技有限公司 | 高速铁路轨道精测仪 |
CN201161592Y (zh) * | 2008-03-06 | 2008-12-10 | 北京光电技术研究所 | 轨道平顺度检测车 |
KR20100068595A (ko) * | 2008-12-15 | 2010-06-24 | 한국철도기술연구원 | 관성센서를 이용한 궤도의 수평 틀림 측정 시스템 및 그 방법 |
CN201575794U (zh) * | 2009-11-27 | 2010-09-08 | 刘洪云 | 高速铁路轨道长波不平顺精密测量仪 |
-
2011
- 2011-04-12 CN CN 201110089805 patent/CN102180187B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4421741A1 (de) * | 1994-06-22 | 1996-01-04 | Ruhrkohle Ag | Vorrichtung zum Messen einer Höhendifferenz an Gleisanlagen |
JPH1030919A (ja) * | 1996-07-16 | 1998-02-03 | Pasuko Doro Gijutsu Center:Kk | 道路縦断プロファイル測定装置及び方法 |
CN2473089Y (zh) * | 2001-04-28 | 2002-01-23 | 上海理工大学附属二厂 | 轨道综合测量仪 |
CN201158385Y (zh) * | 2008-02-25 | 2008-12-03 | 长沙悦诚机电科技有限公司 | 高速铁路轨道精测仪 |
CN201161592Y (zh) * | 2008-03-06 | 2008-12-10 | 北京光电技术研究所 | 轨道平顺度检测车 |
KR20100068595A (ko) * | 2008-12-15 | 2010-06-24 | 한국철도기술연구원 | 관성센서를 이용한 궤도의 수평 틀림 측정 시스템 및 그 방법 |
CN201575794U (zh) * | 2009-11-27 | 2010-09-08 | 刘洪云 | 高速铁路轨道长波不平顺精密测量仪 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103174072B (zh) * | 2011-12-26 | 2016-04-20 | 同济大学 | 基于三点偏弦法的轨面检测方法 |
CN103174072A (zh) * | 2011-12-26 | 2013-06-26 | 同济大学 | 基于三点偏弦法的轨面检测方法 |
CN103552571A (zh) * | 2013-11-18 | 2014-02-05 | 南车洛阳机车有限公司 | 综合检测车的轨道检测*** |
CN104260751A (zh) * | 2014-09-28 | 2015-01-07 | 北京理工大学 | 一种多传感器融合的高铁轨道中心线检测***和检测方法 |
CN104260751B (zh) * | 2014-09-28 | 2016-06-15 | 北京理工大学 | 一种多传感器融合的高铁轨道中心线检测***和检测方法 |
CN104567783A (zh) * | 2014-12-01 | 2015-04-29 | 华南理工大学 | 使用mems陀螺仪的便携式轨道几何量检测***及方法 |
CN104536462A (zh) * | 2015-01-09 | 2015-04-22 | 西安应用光学研究所 | 基于光纤陀螺积分手段的位置控制方法 |
CN105424057A (zh) * | 2015-10-30 | 2016-03-23 | 东南大学 | 一种基于改进lmd的光纤陀螺振动信号分析方法 |
CN105424057B (zh) * | 2015-10-30 | 2019-04-09 | 东南大学 | 一种基于改进lmd的光纤陀螺振动信号分析方法 |
CN109477712A (zh) * | 2016-07-27 | 2019-03-15 | 中国涂料株式会社 | 三维表面粗糙度评价装置、三维表面粗糙度评价方法、三维表面粗糙度数据获取装置和三维表面粗糙度数据获取方法 |
US11162786B2 (en) | 2016-07-27 | 2021-11-02 | Chugoku Marine Paints, Ltd. | Three-dimensional surface roughness evaluating device, three-dimensional surface roughness evaluating method, three-dimensional surface roughness data acquiring device, and three-dimensional surface roughness data acquiring method |
CN109606417A (zh) * | 2018-12-04 | 2019-04-12 | 武汉理工大学 | 一种铁路轨道超高检测装置 |
CN113548068A (zh) * | 2021-07-23 | 2021-10-26 | 中车长春轨道客车股份有限公司 | 轨面不平顺检测装置及检测方法 |
CN113548068B (zh) * | 2021-07-23 | 2023-09-12 | 中车长春轨道客车股份有限公司 | 轨面不平顺检测装置及检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102180187B (zh) | 2013-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102180187B (zh) | 一种铁路轨道高低高精度检测装置和检测方法 | |
CN107554553B (zh) | 基于二维激光位移传感器的轨道几何不平顺检测方法 | |
CN102252636B (zh) | 一种多波长高铁轨道平顺性检测装置和检测方法 | |
CN103207403B (zh) | 一种卫星导航与惯性测量组合轨道测量***与方法 | |
CN102251451B (zh) | 基于多源信息融合的轨道几何状态测量方法 | |
CN102337710B (zh) | 一种gps轨道不平顺检测***及其检测方法 | |
CN102252627B (zh) | 一种高铁轨道轨距检测装置和检测方法 | |
CN202644332U (zh) | 轨道顶面不平顺检测装置 | |
CN101758835A (zh) | 构架式轨距和轨向测量装置 | |
CN202368604U (zh) | 一种动态检测铁路轨道左右钢轨高低的检测装置 | |
CN101804822A (zh) | 激光测量仪及测量方法 | |
CN115597535B (zh) | 基于惯性导航的高速磁悬浮轨道不平顺检测***及方法 | |
CN110657788B (zh) | 一种起重机轨道平顺性动态检测方法 | |
CN107815935A (zh) | 高速铁路轨道几何状态实时监测方法及*** | |
CN203020332U (zh) | 一种卫星导航与惯性测量组合轨道测量*** | |
CN110174089A (zh) | 高速铁路轨面局部沉降检测方法及*** | |
CN102445176A (zh) | 高速列车运行姿态参数测量*** | |
CN103644888A (zh) | 一种用于检测桥梁形变的惯性基准测量方法 | |
CN104047212B (zh) | 一种基于角度量测的轨道沉降自动测量装置及方法 | |
CN110106755A (zh) | 利用姿态重构铁轨几何形态的高铁轨道不平顺性检测方法 | |
CN102175168B (zh) | 中低速磁悬浮f型轨的轨道动态检测设备 | |
CN113415308B (zh) | 钢轨波磨检测方法及装置 | |
CN110658543B (zh) | 一种基于非接触式测量的高速铁路轨道几何参数检测方法 | |
CN111854742B (zh) | 基于多源信息融合的移动物体的测速定位方法和*** | |
CN111778791A (zh) | 低速轨道检测***与方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |