CN102252636A - 一种多波长高铁轨道平顺性检测装置和检测方法 - Google Patents

一种多波长高铁轨道平顺性检测装置和检测方法 Download PDF

Info

Publication number
CN102252636A
CN102252636A CN 201110089802 CN201110089802A CN102252636A CN 102252636 A CN102252636 A CN 102252636A CN 201110089802 CN201110089802 CN 201110089802 CN 201110089802 A CN201110089802 A CN 201110089802A CN 102252636 A CN102252636 A CN 102252636A
Authority
CN
China
Prior art keywords
centerdot
inertial measurement
measurement cluster
computing machine
data
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
Application number
CN 201110089802
Other languages
English (en)
Other versions
CN102252636B (zh
Inventor
张金红
党进
陈欢
焦飞
祝运刚
刘小溪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
No 618 Research Institute of China Aviation Industry
Original Assignee
No 618 Research Institute of China Aviation Industry
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by No 618 Research Institute of China Aviation Industry filed Critical No 618 Research Institute of China Aviation Industry
Priority to CN201110089802A priority Critical patent/CN102252636B/zh
Publication of CN102252636A publication Critical patent/CN102252636A/zh
Application granted granted Critical
Publication of CN102252636B publication Critical patent/CN102252636B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Machines For Laying And Maintaining Railways (AREA)

Abstract

本发明属于铁路轨道检测领域,涉及一种多波长轨道不平顺检测装置和检测方法的改进。本发明的检测装置包括一个带有检测平台[15]的轨道检测车,在检测平台[15]的上表面安装有包括计算机,位移器,位移器数据采集卡[5]、OD数据采集卡[6]和电源的数据处理系统,其特征在于:所说的位移器是激光位移器,在检测梁[16]内部安装有两个激光位移器,在检测平台[15]的上表面安装有一套高低测量系统。本发明检测方法的步骤是:传感器数据采集存储的步骤;数据预处理的步骤;综合数据采集处理。本发明的检测装置和检测方法大大缩短了轨道上的工作时间,提高了检测效率;并满足了高铁对100米以上的长波不平顺的测量需求。

Description

一种多波长高铁轨道平顺性检测装置和检测方法
技术领域
本发明属于铁路轨道检测技术,涉及对铁路轨道高低、轨向等平顺性检测装置和方法的改进。
背景技术
目前,对铁路轨道平顺性的检测方法主要有弦测法、惯性基准法等。其中弦测法是指采用人工拉线的方法在轨道上真实搭建一条弦线,通过测量轨道轨顶面(轨距点)与该弦线的相对位移,评价轨道的高低平顺性。通过这种方法测得的轨道平顺性,波长与检测弦的长度密切相关,当需要分析多种波长的轨道平顺性时,需要更换检测装置重新测量,或者根据当前检测结果进行“以小推大”的转换,前者增加了检测人员的工作量,工作效率低,后者存在较大误差。惯性基准法是通过设计高阶滤波器实现对不同范围的轨道平顺性进行分析,目前技术较为成熟的可以分析70m以内的轨道不平顺。这种方法不能满足高铁对100米以上的长波不平顺测量需求。惯性基准法可参见《GJ-4型轨道检查车的原理与应用》,中国铁道出版社,2000,张未。
发明内容
本发明的目的是:提出一种能对轨道的高低、轨向进行多种波长连续检测的检测装置和检测方法,并能满足高铁对100米以上的长波不平顺测量需求。
本发明的技术方案是:一种高铁轨道平顺性检测装置,包括一个带有检测平台的轨道检测车,在检测平台的下表面固定有两个位移器,它们是第一位移器LPDT和第二位移器RPDT,第一位移器LPDT和第二位移器RPDT分别位于轮轴两端的正上方,第一位移器LPDT和第二位移器RPDT的轴线与轮轴的轴线正交,在轮轴的一端安装有里程仪OD,在检测平台的上表面安装有数据处理系统,该数据处理系统包括计算机、位移器数据采集控制盒、OD数据采集卡和电源;其特征在于:
(1)在检测平台下表面固定有一个截面为矩形的、中空的检测梁,检测梁与轮轴平行,检测梁横向的垂直平分面与轮轴横向的垂直平分面共面;检测梁的谐振频率不小于70Hz;在检测梁内部安装有两个位移器,分别是第三激光位移器GLPDT和第四激光位移器GRPDT,在检测梁的下表面、与第三激光位移器GLPDT和第四激光位移器GRPDT对应的位置各有一个通光孔,第三激光位移器GLPDT和第四激光位移器GRPDT发射的激光束分别穿过上述通光孔后照射到左右股钢轨的轨距点上;在检测梁内部还装有一个轨向惯性测量组件HIMU;在检测平台的上表面安装有一个轨向-里程预处理计算机,轨向惯性测量组件HIMU的输出端通过RS422总线与轨向-里程预处理计算机的串口连接;所说的第一位移器LPDT和第二位移器RPDT是激光位移器;
(2)在检测平台的上表面安装有一套高低测量系统,它由左惯性测量组件LIMU、右惯性测量组件RIMU和高低预处理计算机组成;左惯性测量组件LIMU和右惯性测量组件RIMU位于轮轴左右车轮的正上方,左惯性测量组件LIMU和右惯性测量组件RIMU的垂向轴线分别与所对应的左轨道、右轨道横截面的垂直中心线重合;左惯性测量组件LIMU的输出端通过RS422总线与高低预处理计算机的第一串口Com1连接;右惯性测量组件RIMU的输出端通过RS422总线与高低预处理计算机的第二串口Com2连接;
(3)第一激光位移器LPDT、第二激光位移器RPDT、第三激光位移器GLPDT和第四激光位移器GRPDT的输出端分别与位移器数据采集盒的相应输入端连接;在检测平台的上表面安装有一个位移器预处理计算机,位移器数据采集盒的输出端通过USB总线与位移器预处理计算机的USB口连接;
(4)高低预处理计算机、位移预处理计算机和轨向-里程预处理计算机通过局域网与计算机连接。
使用如上面所述的高铁轨道平顺性检测装置检测铁路轨道平顺性的方法,其特征在于,检测的步骤如下:
1、传感器数据采集存储的步骤:
1.1、高低数据采集存储:系统上电后,左惯性测量组件LIMU和右惯性测量组件开始工作,高低预处理计算机每隔m秒读取第一串口com1和第二串口com2的数据,将数据记为ωxL,ωyL,ωzL,AxL,AyL,AzL,ωxR,ωyR,ωzR,AxR,AyR,AzR,并进行保存;其中,ωxL表示左惯性测量组件LIMU的横滚角速率,ωyL表示左惯性测量组件LIMU的航向角速率,ωzL表示左惯性测量组件LIMU的俯仰角速率,AxL,AyL,AzL分别表示左惯性测量组件LIMU的横向、纵向、天向加速度;ωxR表示右惯性测量组件RIMU的横滚角速率,ωyR表示右惯性测量组件RIMU的航向角速率,ωzR表示右惯性测量组件RIMU的俯仰角速率,AxR,AyR,AzR分别表示右惯性测量组件RIMU的横向、纵向、天向加速度;m的取值范围是0.001秒~0.01秒;
1.2、轨向数据采集存储:系统上电后,轨向惯性测量组件HIMU开始工作,每隔m秒,轨向惯性测量组件HIMU向轨向-里程预处理计算机第一串口com1发送采样数据,记为ωxH,ωyH,ωzH,AxH,AyH,AzH;其中ωxH表示轨向惯性测量组件HIMU的横滚角速率,其中ωyH表示轨向惯性测量组件HIMU的航向角速率,其中ωzH表示轨向惯性测量组件HIMU的俯仰角速率,AxH,AyH,AzH分别表示轨向惯性测量组件HIMU的横向、纵向、天向加速度,并进行保存;
1.3、里程数据采集存储:系统上电后,里程测量组件OD开始工作,每隔m秒,里程测量组件OD向轨向-里程预处理计算机的第二串口com2发送一个采样数据,轨向-里程预处理计算机将采样数据记为Nod并进行保存;
1.4、激光位移器采集存储:系统上电后,激光位移器数据控制盒开始工作,每隔m秒,激光位移器数据控制盒向位移器预处理计算机发送一个采样数据包,每个采样数据包有4个采样数据,记为DWL、DWR、DGL、DGR,其中,DWL为第一激光位移器LPDT所测量位移的采样数据;DWR为第二激光位移器RPDT所测量位移的采样数据;DGL为第三激光位移器GLPDT所测量位移的采样数据;DGR为第四激光位移器GRPDT所测量位移的采样数据,并进行保存;
2、数据预处理的步骤:
2.1、高低预处理:高低预处理计算机对数据的预处理过程是:
2.1.1、左惯性测量组件LIMU数据预处理:
2.1.1.1、计算左惯性测量组件LIMU的当前姿态:保持轨道车静止m1秒,m1取值范围为100秒~200秒,在m1秒内的采样个数为n,n=m1/m取整数;
根据下式计算出左惯性测量组件LIMU的横滚角速率均值
ω ‾ xL = 1 n Σ i = 1 n ω xL ( i ) · · · [ 1 ]
式中,i为采样数据的序号;
计算航向角初值ψ0
Figure BSA00000470733400033
式中,
Figure BSA00000470733400034
表示当前位置的纬度;
2.1.1.2、地球自转和常值漂移补偿:
Figure BSA00000470733400036
式中,
Figure BSA00000470733400037
为补偿后的左惯性测量组件LIMU的横滚角速率,
Figure BSA00000470733400038
为补偿后的左惯性测量组件LIMU的航向角速率,δψ为航向角修正量,当i=1时,令δψ为0;
1.1.3、积分计算:
δψ ( i ) = δψ ( i - 1 ) + ω ~ zL ( i ) · m · π / 180 · · · [ 5 ]
θ L ( i ) = θ L ( i - 1 ) + ω ~ xL ( i ) · m · π / 180 · · · [ 6 ]
式中,θL(i)为左惯性测量组件LIMU的俯仰角;
2.1.2右惯性测量组件RIMU数据预处理:
2.1.2.1、计算右惯性测量组件RIMU的当前姿态:保持轨道车静止m1秒;根据下式计算出右惯性测量组件RIMU横滚角速率均值
ω ‾ xR = 1 n Σ i = 1 n ω xR ( i ) . . . [ 7 ]
2.1.2.2、地球自转和常值漂移补偿:
ω ~ xR ( i ) = ω xR ( i ) - 15 · cos L · cos ( ψ 0 + δψ ( i - 1 ) ) · · · [ 8 ]
式中,
Figure BSA00000470733400046
为补偿后的右惯性测量组件RIMU横滚角速率;
2.1.1.3、积分计算:
θ R ( i ) = θ R ( i - 1 ) + ω ~ xR ( i ) · m · π / 180 · · · [ 10 ]
式中,θR(i)为右惯性测量组件RIMU的俯仰角;
2.1.4、数据传输:高低预处理计算机每m秒通过局域网向计算机发送一次数据θR(i)和θR(i);
2.2、轨向-里程数据预处理:
2.2.1、计算轨向惯性测量组件HIMU的当前航向角:保持轨道车静止m1秒,计算航向角速率均值
Figure BSA00000470733400048
ω ‾ xH = 1 n Σ i = 1 n ω xH ( i ) · · · [ 11 ]
航向角初值ψ0H计算:
2.2.2、地球自转和常值漂移补偿:
式中,为补偿后的航向角速率,当i=1时,令ψH(i)=ψ0H为0;
2.2.3、积分计算
ψ H ( i ) = ψ H ( i - 1 ) + ω ~ zH ( i ) · m · π / 180 · · · [ 14 ]
式中,ψH(i)为轨向惯性测量组件HIMU的航向角;
2.2.4、数据传输:
轨向-里程预处理计算机每m秒通过局域网向计算机发送一次ψH(i)和Nod(i);
2.3、激光位移器数据预处理:
2.3.1、计算平均值:分别计算DWL、DWH、DGL、DGR四个采样数据中每个采样数据每连续5个数的平均值
Figure BSA00000470733400051
其中,
Figure BSA00000470733400052
为第一激光位移器LPDT所测量位移的采样数据平均值;
Figure BSA00000470733400053
为第二激光位移器RPDT所测量位移的采样数据平均值;
Figure BSA00000470733400054
为第三激光位移器GLPDT所测量位移的采样数据平均值;
Figure BSA00000470733400055
为第四激光位移器GRPDT所测量位移的采样数据平均值;
2.3.2、数据传输:位移预处理计算机每m秒通过局域网向计算机发送一次
3、综合数据采集处理:
3.1、综合数据采集:计算机通过局域网获取高低预处理计算机、位移预处理计算机和轨向-里程预处理计算机的数据并存储;
3.2、里程计算:根据下式计算当前里程:
里程仪的采样数据和事先标定好的刻度系数,
Lod(i)=Nod(i)·kod…………………………………………………[15]
式中,Lod(i)表示检测里程,kod是刻度系数,
Figure BSA00000470733400057
R是轨道车的车轮半径,其中p的取值范围是1024~4000;
3.3、时域-空域转换:
将步骤3.1计算机接收到的时域采样数据,根据步骤3.2得到的里程,转换成采样密度为0.125米的空域数据,记为ΘL、ΘR、ΨH、L、
Figure BSA00000470733400059
其中ΘL为空域的左惯性测量组件LIMU的俯仰角、ΘR为空域的右惯性测量组件RIMU的俯仰角;ΨH为空域的轨向惯性测量组件HIMU航向角,L为空域的里程;为空域的第一激光位移器LPDT所测量位移的采样数据平均值;为空域的第二激光位移器RPDT所测量位移的采样数据平均值;
Figure BSA000004707334000512
为空域的第三激光位移器GLPDT所测量位移的采样数据平均值;为空域的第四激光位移器GRPDT所测量位移的采样数据平均值,j为空域数据的序号;
3.4、当前波长平顺性解算:
3.4.1、波长输入:用户输入的波长λ1为数值10m、30m、42m、70m、300m其中之一;
3.4.2、平顺性解算:
3.4.2.1、高低平顺性解算:
h VL ( j ) = α [ dL 1 · tan ( dΘ 1 ) - 1 2 dL 2 · tan ( dΘ 2 ) ] + D ^ WL ( j ) · · · [ 16 ]
h VR ( j ) = α [ dL 1 · tan ( dΘ 3 ) - 1 2 dL 2 · tan ( dΘ 4 ) ] + D ^ WR ( j ) · · · [ 17 ]
式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
dL1=L(j)-L(j-d)……………………………………………………… [18]
dL2=L(j+d)-L(j-d)…………………………………………………… [19]
1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
d = λ 1 0.25 · · · [ 24 ]
3.4.2.2、轨向平顺性解算:
S L ( j ) = y ( j ) + D ^ GL ( j ) · · · [ 25 ]
S R ( j ) = y ( j ) + D ^ GR ( j ) · · · [ 26 ]
y ( j ) = β · [ dL 1 tan ( dψ 1 ) - 1 2 dL 2 tan ( dψ 2 ) ] · · · [ 27 ]
式中,y(j)为检测梁[16]的中心横向位移,β为幅值补偿系数,取0.03;
1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
至此,得到了轨道平顺性的检测结果。
本发明的优点是:使用本发明的检测装置和检测方法可以通过一次检测和多次解算,就可以得对10m、30m、42m、70m或300m多种波长的轨道不平顺进行检测,大大缩短了轨道上的工作时间,提高了检测效率;并满足了高铁对100米以上的长波不平顺的测量需求。
附图说明
图1是本发明检测装置的结构原理框图。
图2是本发明检测装置中传感器安装位置示意图。
具体实施方式
下面对本发明做进一步详细说明。规定方位如下,面对轨道检测车的前进方向,左手边为左方。参见图1、2。一种高铁轨道平顺性检测装置,包括一个带有检测平台15的轨道检测车,在检测平台15的下表面固定有两个位移器,它们是第一位移器LPDT9和第二位移器RPDT10,第一位移器LPDT9和第二位移器RPDT10分别位于轮轴17两端的正上方,第一位移器LPDT9和第二位移器RPDT10的轴线与轮轴17的轴线正交,在轮轴17的一端安装有里程仪OD14,在检测平台15的上表面安装有数据处理系统,该数据处理系统包括计算机1、位移器数据采集控制盒5、OD数据采集卡6和电源;其特征在于:
(1)在检测平台15下表面固定有一个截面为矩形的、中空的检测梁16,检测梁16与轮轴17平行,检测梁16横向的垂直平分面与轮轴横向的垂直平分面共面;检测梁16的谐振频率不小于70Hz;在检测梁16内部安装有两个位移器,分别是第三激光位移器GLPDT11和第四激光位移器GRPDT12,在检测梁16的下表面、与第三激光位移器GLPDT11和第四激光位移器GRPDT12对应的位置各有一个通光孔,第三激光位移器GLPDT11和第四激光位移器GRPDT12发射的激光束分别穿过上述通光孔后照射到左右股钢轨的轨距点上;在检测梁16内部还装有一个轨向惯性测量组件HIMU13;在检测平台15的上表面安装有一个轨向-里程预处理计算机4,轨向惯性测量组件HIMU13的输出端通过RS422总线与轨向-里程预处理计算机4的串口连接;所说的第一位移器LPDT9和第二位移器RPDT10是激光位移器;
(2)在检测平台15的上表面安装有一套高低测量系统,它由左惯性测量组件LIMU7、右惯性测量组件RIMU8和高低预处理计算机2组成;左惯性测量组件LIMU7和右惯性测量组件RIMU8位于轮轴17左右车轮的正上方,左惯性测量组件LIMU7和右惯性测量组件RIMU8的垂向轴线分别与所对应的左轨道、右轨道横截面的垂直中心线重合;左惯性测量组件LIMU7的输出端通过RS422总线与高低预处理计算机2的第一串口Com1连接;右惯性测量组件RIMU8的输出端通过RS422总线与高低预处理计算机2的第二串口Com2连接;
(3)第一激光位移器LPDT9、第二激光位移器RPDT10、第三激光位移器GLPDT11和第四激光位移器GRPDT12的输出端分别与位移器数据采集盒5的相应输入端连接;在检测平台15的上表面安装有一个位移器预处理计算机3,位移器数据采集盒5的输出端通过USB总线与位移器预处理计算机3的USB口连接;
(4)高低预处理计算机2、位移预处理计算机3和轨向-里程预处理计算机4通过局域网与计算机1连接。
本发明检测装置的检测原理是:利用3个光纤惯组分别测量轨道车的左、右端俯仰角和检测梁的航向角,结合里程和检测波长计算得到车体的高低和检测梁中心位移,结合第一激光位移器和第二激光位移器的测量结果,扣除车体相对于轮轴的振动量,计算得到左、右股轨道高低平顺性测量结果,结合第三激光位移器和第四激光位移器的测量结果,加上对应的左、右轨距点位移,计算得到左、右股轨道轨向平顺性测量结果。
使用如上面所述的高铁轨道平顺性检测装置检测铁路轨道平顺性的方法,其特征在于,检测的步骤如下:
1、传感器数据采集存储的步骤:
1.1、高低数据采集存储:系统上电后,左惯性测量组件LIMU7和右惯性测量组件8开始工作,高低预处理计算机2每隔m秒读取第一串口com1和第二串口com2的数据,将数据记为ωxL,ωyL,ωzL,AxL,AyL,AzL,ωxR,ωyR,ωzR,AxR,AyR,AzR,并进行保存;其中,ωxL表示左惯性测量组件LIMU7的横滚角速率,ωyL表示左惯性测量组件LIMU7的航向角速率,ωzL表示左惯性测量组件LIMU7的俯仰角速率,AxL,AyL,AzL分别表示左惯性测量组件LIMU7的横向、纵向、天向加速度;ωxR表示右惯性测量组件RIMU8的横滚角速率,ωyR表示右惯性测量组件RIMU8的航向角速率,ωzR表示右惯性测量组件RIMU8的俯仰角速率,AxR,AyR,AzR分别表示右惯性测量组件RIMU8的横向、纵向、天向加速度;m的取值范围是0.001秒~0.01秒;
1.2、轨向数据采集存储:系统上电后,轨向惯性测量组件HIMU13开始工作,每隔m秒,轨向惯性测量组件HIMU13向轨向-里程预处理计算机4第一串口com1发送采样数据,记为ωxH,ωyH,ωzH,AxH,AyH,AzH;其中ωxH表示轨向惯性测量组件HIMU13的横滚角速率,其中ωyH表示轨向惯性测量组件HIMU13的航向角速率,其中ωzH表示轨向惯性测量组件HIMU13的俯仰角速率,AxH,AyH,AzH分别表示轨向惯性测量组件HIMU13的横向、纵向、天向加速度,并进行保存;
1.3、里程数据采集存储:系统上电后,里程测量组件OD14开始工作,每隔m秒,里程测量组件OD14向轨向-里程预处理计算机4的第二串口com2发送一个采样数据,轨向-里程预处理计算机4将采样数据记为Nod并进行保存;
1.4、激光位移器采集存储:系统上电后,激光位移器数据控制盒5开始工作,每隔m秒,激光位移器数据控制盒5向位移器预处理计算机3发送一一个采样数据包,每个采样数据包有4个采样数据,记为DWL、DWR、DGL、DGR,其中,DWL为第一激光位移器LPDT9所测量位移的采样数据;DWR为第二激光位移器RPDT10所测量位移的采样数据;DGL为第三激光位移器GLPDT11所测量位移的采样数据;DGR为第四激光位移器GRPDT12所测量位移的采样数据,并进行保存;
2、数据预处理的步骤:
2.1、高低预处理:高低预处理计算机2对数据的预处理过程是:
2.1.1、左惯性测量组件LIMU7数据预处理:
2.1.1.1、计算左惯性测量组件LIMU7的当前姿态:保持轨道车静止m1秒,m1取值范围为100秒~200秒,在m1秒内的采样个数为n,n=m1/m取整数;
根据下式计算出左惯性测量组件LIMU7的横滚角速率均值
Figure BSA00000470733400091
ω ‾ xL = 1 n Σ i = 1 n ω xL ( i ) · · · [ 1 ]
式中,i为采样数据的序号;
计算航向角初值ψ0
Figure BSA00000470733400093
式中,
Figure BSA00000470733400094
表示当前位置的纬度;
2.1.1.2、地球自转和常值漂移补偿:
Figure BSA00000470733400095
Figure BSA00000470733400096
式中,
Figure BSA00000470733400097
为补偿后的左惯性测量组件LIMU7的横滚角速率,
Figure BSA00000470733400098
为补偿后的左惯性测量组件LIMU7的航向角速率,δψ为航向角修正量,当i=1时,令δψ为0;
1.1.3、积分计算:
δψ ( i ) = δψ ( i - 1 ) + ω ~ zL ( i ) · m · π / 180 · · · [ 5 ]
θ L ( i ) = θ L ( i - 1 ) + ω ~ xL ( i ) · m · π / 180 · · · [ 6 ]
式中,θL(i)为左惯性测量组件LIMU7的俯仰角;
2.1.2右惯性测量组件RIMU8数据预处理:
2.1.2.1、计算右惯性测量组件RIMU8的当前姿态:保持轨道车静止m1秒;根据下式计算出右惯性测量组件RIMU8横滚角速率均值
Figure BSA00000470733400101
ω ‾ xR = 1 n Σ i = 1 n ω xR ( i ) · · · [ 7 ]
2.1.2.2、地球自转和常值漂移补偿:
ω ~ xR ( i ) = ω xR ( i ) - 15 · cos L · cos ( ψ 0 + δψ ( i - 1 ) ) · · · [ 8 ]
式中,
Figure BSA00000470733400104
为补偿后的右惯性测量组件RIMU8横滚角速率;
2.1.1.3、积分计算:
θ R ( i ) = θ R ( i - 1 ) + ω ~ xR ( i ) · m · π / 180 · · · [ 10 ]
式中,θR(i)为右惯性测量组件RIMU8的俯仰角;
2.1.4、数据传输:高低预处理计算机2每m秒通过局域网向计算机1发送一次数据θL(i)和θR(i);
2.2、轨向-里程数据预处理:
2.2.1、计算轨向惯性测量组件HIMU13的当前航向角:保持轨道车静止m1秒,计算航向角速率均值
Figure BSA00000470733400106
ω ‾ xH = 1 n Σ i = 1 n ω xH ( i ) · · · [ 11 ]
航向角初值ψ0H计算:
2.2.2、地球自转和常值漂移补偿:
Figure BSA00000470733400109
式中,
Figure BSA000004707334001010
为补偿后的航向角速率,当i=1时,令ψH(i)=ψ0H为0;
2.2.3、积分计算
ψ H ( i ) = ψ H ( i - 1 ) + ω ~ zH ( i ) · m · π / 180 · · · [ 14 ]
式中,ψH(i)为轨向惯性测量组件HIMU13的航向角;
2.2.4、数据传输:
轨向-里程预处理计算机4每m秒通过局域网向计算机1发送一次ψH(i)和Nod(i);
2.3、激光位移器数据预处理:
2.3.1、计算平均值:分别计算DWL、DWR、DGL、DGR四个采样数据中每个采样数据每连续5个数的平均值
Figure BSA000004707334001012
其中,
Figure BSA000004707334001013
为第一激光位移器LPDT9所测量位移的采样数据平均值;
Figure BSA000004707334001014
为第二激光位移器RPDT10所测量位移的采样数据平均值;
Figure BSA000004707334001015
为第三激光位移器GLPDT11所测量位移的采样数据平均值;为第四激光位移器GRPDT12所测量位移的采样数据平均值;
2.3.2、数据传输:位移预处理计算机3每m秒通过局域网向计算机1发送一次
Figure BSA00000470733400112
3、综合数据采集处理:
3.1、综合数据采集:计算机1通过局域网获取高低预处理计算机2、位移预处理计算机3和轨向-里程预处理计算机4的数据并存储;
3.2、里程计算:根据下式计算当前里程:
里程仪的采样数据和事先标定好的刻度系数,
Lod(i)=Nod(i)·kod…………………………………………………[15]
式中,Lod(i)表示检测里程,kod是刻度系数,R是轨道车的车轮半径,其中p的取值范围是1024~4000;
3.3、时域-空域转换:
将步骤3.1计算机1接收到的时域采样数据,根据步骤3.2得到的里程,转换成采样密度为0.125米的空域数据,记为ΘL、ΘR、ΨH、L、
Figure BSA00000470733400114
Figure BSA00000470733400115
其中ΘL为空域的左惯性测量组件LIMU7的俯仰角、ΘR为空域的右惯性测量组件RIMU8的俯仰角;ΨH为空域的轨向惯性测量组件HIMU13航向角,L为空域的里程;
Figure BSA00000470733400116
为空域的第一激光位移器LPDT9所测量位移的采样数据平均值;
Figure BSA00000470733400117
为空域的第二激光位移器RPDT10所测量位移的采样数据平均值;
Figure BSA00000470733400118
为空域的第三激光位移器GLPDT11所测量位移的采样数据平均值;
Figure BSA00000470733400119
为空域的第四激光位移器GRPDT12所测量位移的采样数据平均值,j为空域数据的序号;
3.4、当前波长平顺性解算:
3.4.1、波长输入:用户输入的波长λ1为数值10m、30m、42m、70m、300m其中之一;
3.4.2、平顺性解算:
3.4.2.1、高低平顺性解算:
h VL ( j ) = α [ dL 1 · tan ( dΘ 1 ) - 1 2 dL 2 · tan ( dΘ 2 ) ] + D ^ WL ( j ) · · · [ 16 ]
h VR ( j ) = α [ dL 1 · tan ( dΘ 3 ) - 1 2 dL 2 · tan ( dΘ 4 ) ] + D ^ WR ( j ) · · · [ 17 ]
式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
dL1=L(j)-L(j-d)……………………………………………………… [18]
dL2=L(j+d)-L(j-d)…………………………………………………… [19]
1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
d = λ 1 0.25 · · · [ 24 ]
3.4.2.2、轨向平顺性解算:
S L ( j ) = y ( j ) + D ^ GL ( j ) · · · [ 25 ]
S R ( j ) = y ( j ) + D ^ GR ( j ) · · · [ 26 ]
y ( j ) = β · [ dL 1 tan ( dψ 1 ) - 1 2 dL 2 tan ( dψ 2 ) ] · · · [ 27 ]
式中,y(j)为检测梁[16]的中心横向位移,β为幅值补偿系数,取0.03;
1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
至此,得到了轨道平顺性的检测结果。
实施例1
在装有检测平台的轨道车上安装如上面所述的检测装置,其中传感器的技术指标分别为:左惯性测量组件LIMU7、右惯性测量组件RIMU8、轨向惯性测量组件HIMU13均选用光纤捷联惯组,第一位移器LPDT9、第二位移器RPDT10、第三激光位移器GLPDT11、第四激光位移器GRPDT 12均选用测量精度为0.05mm的LK-400G激光位移器,OD14选用3600脉冲/转的光电编码器。检测的步骤如下:
1、传感器数据采集存储的步骤:
1.1、高低数据采集存储:系统上电后,左惯性测量组件LIMU7和右惯性测量组件8开始工作,高低预处理计算机2每隔m秒读取第一串口com1和第二串口com2的数据,将数据记为ωxL,ωyL,ωzL,AxL,AyL,AzL,ωxR,ωyR,ωzR,AxR,AyR,AzR,并进行保存;其中,ωxL表示左惯性测量组件LIMU7的横滚角速率,ωyL表示左惯性测量组件LIMU7的航向角速率,ωzL表示左惯性测量组件LIMU7的俯仰角速率,AxL,AyL,AzL分别表示左惯性测量组件LIMU7的横向、纵向、天向加速度;ωxR表示右惯性测量组件RIMU8的横滚角速率,ωyR表示右惯性测量组件RIMU8的航向角速率,ωzR表示右惯性测量组件RIMU8的俯仰角速率,AxR,AyR,AzR分别表示右惯性测量组件RIMU8的横向、纵向、天向加速度;m的取值是0.005秒;
1.2、轨向数据采集存储:系统上电后,轨向惯性测量组件HIMU13开始工作,每隔m秒,轨向惯性测量组件HIMU13向轨向-里程预处理计算机4第一串口com1发送采样数据,记为ωxH,ωyH,ωzH,AxH,AyH,AzH;其中ωxH表示轨向惯性测量组件HIMU13的横滚角速率,其中ωyH表示轨向惯性测量组件HIMU13的航向角速率,其中ωzH表示轨向惯性测量组件HIMU13的俯仰角速率,AxH,AyH,AzH分别表示轨向惯性测量组件HIMU13的横向、纵向、天向加速度,并进行保存;
1.3、里程数据采集存储:系统上电后,里程测量组件OD14开始工作,每隔m秒,里程测量组件OD14向轨向-里程预处理计算机4的第二串口com2发送一个采样数据,轨向-里程预处理计算机4将采样数据记为Nod并进行保存;
1.4、激光位移器采集存储:系统上电后,激光位移器数据控制盒5开始工作,每隔m秒,激光位移器数据控制盒5向位移器预处理计算机3发送一次采样数据,每次有4个采样数据,记为DWL、DWR、DGL、DGR,其中,DWL为第一激光位移器LPDT9所测量位移的采样数据;DWR为第二激光位移器RPDT10所测量位移的采样数据;DGL为第三激光位移器GLPDT11所测量位移的采样数据;DGR为第四激光位移器GRPDT12所测量位移的采样数据,并进行保存;
2、数据预处理的步骤:
2.1、高低预处理:高低预处理计算机2对数据的预处理过程是:
2.1.1、左惯性测量组件LIMU7数据预处理:
2.1.1.1、计算左惯性测量组件LIMU7的当前姿态:保持轨道车静止m1秒,m1取值范围为150秒,在m1秒内的采样个数为n,n=m1/m取整数;
根据下式计算出左惯性测量组件LIMU7的横滚角速率均值
Figure BSA00000470733400131
ω ‾ xL = 1 n Σ i = 1 n ω xL ( i ) · · · [ 1 ]
式中,i为采样数据的序号;
计算航向角初值ψ0
Figure BSA00000470733400133
式中,
Figure BSA00000470733400134
表示当前位置的纬度;
2.1.1.2、地球自转和常值漂移补偿:
Figure BSA00000470733400141
Figure BSA00000470733400142
式中,
Figure BSA00000470733400143
为补偿后的左惯性测量组件LIMU7的横滚角速率,
Figure BSA00000470733400144
为补偿后的左惯性测量组件LIMU7的航向角速率,δψ为航向角修正量,当i=1时,令δψ为0;
1.1.3、积分计算:
δψ ( i ) = δψ ( i - 1 ) + ω ~ zL ( i ) · m · π / 180 · · · [ 5 ]
θ L ( i ) = θ L ( i - 1 ) + ω ~ xL ( i ) · m · π / 180 · · · [ 6 ]
式中,θL(i)为左惯性测量组件LIMU7的俯仰角;
2.1.2右惯性测量组件RIMU8数据预处理:
2.1.2.1、计算右惯性测量组件RIMU8的当前姿态:保持轨道车静止m1秒;根据下式计算出右惯性测量组件RIMU8横滚角速率均值
Figure BSA00000470733400147
ω ‾ xR = 1 n Σ i = 1 n ω xR ( i ) · · · [ 7 ]
2.1.2.2、地球自转和常值漂移补偿:
ω ~ xR ( i ) = ω xR ( i ) - 15 · cos L · cos ( ψ 0 + δψ ( i - 1 ) ) · · · [ 8 ]
式中,
Figure BSA000004707334001410
为补偿后的右惯性测量组件RIMU8横滚角速率;
2.1.1.3、积分计算:
θ R ( i ) = θ R ( i - 1 ) + ω ~ xR ( i ) · m · π / 180 · · · [ 10 ]
式中,θR(i)为右惯性测量组件RIMU8的俯仰角;
2.1.4、数据传输:高低预处理计算机2每m秒通过局域网向计算机1发送一次数据θL(i)和θR(i);
2.2、轨向-里程数据预处理:
2.2.1、计算轨向惯性测量组件HIMU13的当前航向角:保持轨道车静止m1秒,计算航向角速率均值
ω ‾ xH = 1 n Σ i = 1 n ω xH ( i ) · · · [ 11 ]
航向角初值ψ0H计算:
2.2.2、地球自转和常值漂移补偿:
Figure BSA000004707334001415
式中,
Figure BSA000004707334001416
为补偿后的航向角速率,当i=1时,令ψH(i)=ψ0H为0;
2.2.3、积分计算
ψ H ( i ) = ψ H ( i - 1 ) + ω ~ zH ( i ) · m · π / 180 · · · [ 14 ]
式中,ψH(i)为轨向惯性测量组件HIMU13的航向角;
2.2.4、数据传输:
轨向-里程预处理计算机4每m秒通过局域网向计算机1发送一次ψH(i)和Nod(i);
2.3、激光位移器数据预处理:
2.3.1、计算平均值:分别计算DWL、DWR、DGL、DGR四个采样数据中每个采样数据每连续5个数的平均值
Figure BSA00000470733400152
其中,
Figure BSA00000470733400153
为第一激光位移器LPDT9所测量位移的采样数据平均值;
Figure BSA00000470733400154
为第二激光位移器RPDT10所测量位移的采样数据平均值;为第三激光位移器GLPDT11所测量位移的采样数据平均值;为第四激光位移器GRPDT12所测量位移的采样数据平均值;
2.3.2、数据传输:位移预处理计算机3每m秒通过局域网向计算机1发送一次
Figure BSA00000470733400157
3、综合数据采集处理:
3.1、综合数据采集:计算机1通过局域网获取高低预处理计算机2、位移预处理计算机3和轨向-里程预处理计算机4的数据并存储;
3.2、里程计算:根据下式计算当前里程:
里程仪的采样数据和事先标定好的刻度系数,
Lod(i)=Nod(i)·kod…………………………………………………[15]
式中,Lod(i)表示检测里程,kod是刻度系数,
Figure BSA00000470733400158
R是轨道车的车轮半径,为0.915米,其中p的取值是3600;
3.3、时域-空域转换:
将步骤3.1计算机1接收到的时域采样数据,根据步骤3.2得到的里程,转换成采样密度为0.125米的空域数据,记为ΘL、ΘR、ΨH、L、
Figure BSA00000470733400159
Figure BSA000004707334001510
其中ΘL为空域的左惯性测量组件LIMU7的俯仰角、ΘR为空域的右惯性测量组件RIMU8的俯仰角;ΨH为空域的轨向惯性测量组件HIMU13航向角,L为空域的里程;
Figure BSA000004707334001511
为空域的第一激光位移器LPDT9所测量位移的采样数据平均值;为空域的第二激光位移器RPDT10所测量位移的采样数据平均值;
Figure BSA000004707334001513
为空域的第三激光位移器GLPDT11所测量位移的采样数据平均值;
Figure BSA000004707334001514
为空域的第四激光位移器GRPDT12所测量位移的采样数据平均值,j为空域数据的序号;
3.4、当前波长平顺性解算:
3.4.1、波长输入:用户输入的波长λ1为数值30m;
3.4.2、平顺性解算:
3.4.2.1、高低平顺性解算:
h VL ( j ) = α [ dL 1 · tan ( dΘ 1 ) - 1 2 dL 2 · tan ( dΘ 2 ) ] + D ^ WL ( j ) · · · [ 16 ]
h VR ( j ) = α [ dL 1 · tan ( dΘ 3 ) - 1 2 dL 2 · tan ( dΘ 4 ) ] + D ^ WR ( j ) · · · [ 17 ]
式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
dL1=L(j)-L(j-d)……………………………………………………… [18]
dL2=L(j+d)-L(j-d)…………………………………………………… [19]
1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
d = λ 1 0.25 · · · [ 24 ]
3.4.2.2、轨向平顺性解算:
S L ( j ) = y ( j ) + D ^ GL ( j ) · · · [ 25 ]
S R ( j ) = y ( j ) + D ^ GR ( j ) · · · [ 26 ]
y ( j ) = β · [ dL 1 tan ( dψ 1 ) - 1 2 dL 2 tan ( dψ 2 ) ] · · · [ 27 ]
式中,y(j)为检测梁[16]的中心横向位移,β为幅值补偿系数,取0.03;
1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
至此,得到了轨道平顺性的检测结果,在30米波长下的轨道左高低平顺性、右高低平顺性和轨向平顺性检测精度均为±1毫米,小于用户要求的理论精度±2毫米。
实施例2:
检测装置和使用的检测方法,从步骤1到步骤3.3与实施例1完全相同。
3.4、当前波长平顺性解算:
3.4.1、波长输入:用户输入的波长λ1为数值300m;
3.4.2、平顺性解算:
3.4.2.1、高低平顺性解算:
h VL ( j ) = α [ dL 1 · tan ( dΘ 1 ) - 1 2 dL 2 · tan ( dΘ 2 ) ] + D ^ WL ( j ) · · · [ 16 ]
h VR ( j ) = α [ dL 1 · tan ( dΘ 3 ) - 1 2 dL 2 · tan ( dΘ 4 ) ] + D ^ WR ( j ) · · · [ 17 ]
式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
dL1=L(j)-L(j-d)……………………………………………………… [18]
dL2=L(j+d)-L(j-d)…………………………………………………… [19]
1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
d = λ 1 0.25 · · · [ 24 ]
3.4.2.2、轨向平顺性解算:
S L ( j ) = y ( j ) + D ^ GL ( j ) · · · [ 25 ]
S R ( j ) = y ( j ) + D ^ GR ( j ) · · · [ 26 ]
y ( j ) = β · [ dL 1 tan ( dψ 1 ) - 1 2 dL 2 tan ( dψ 2 ) ] · · · [ 27 ]
式中,y(j)为检测梁[16]的中心横向位移,β为幅值补偿系数,取0.03;
1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
至此,得到了轨道平顺性的检测结果,在300米波长下的轨道左高低平顺性、右高低平顺性和轨向平顺性检测精度均为±8毫米,小于用户要求的理论精度±10毫米。

Claims (2)

1.一种高铁轨道平顺性检测装置,包括一个带有检测平台[15]的轨道检测车,在检测平台[15]的下表面固定有两个位移器,它们是第一位移器LPDT[9]和第二位移器RPDT[10],第一位移器LPDT[9]和第二位移器RPDT[10]分别位于轮轴[17]两端的正上方,第一位移器LPDT[9]和第二位移器RPDT[10]的轴线与轮轴[17]的轴线正交,在轮轴[17]的一端安装有里程仪OD[14],在检测平台[15]的上表面安装有数据处理系统,该数据处理系统包括计算机[1]、位移器数据采集控制盒[5]、OD数据采集卡[6]和电源;其特征在于:
(1)在检测平台[15]下表面固定有一个截面为矩形的、中空的检测梁[16],检测梁[16]与轮轴[17]平行,检测梁[16]横向的垂直平分面与轮轴横向的垂直平分面共面;检测梁[16]的谐振频率不小于70Hz;在检测梁[16]内部安装有两个位移器,分别是第三激光位移器GLPDT[11]和第四激光位移器GRPDT[12],在检测梁[16]的下表面、与第三激光位移器GLPDT[11]和第四激光位移器GRPDT[12]对应的位置各有一个通光孔,第三激光位移器GLPDT[11]和第四激光位移器GRPDT[12]发射的激光束分别穿过上述通光孔后照射到左右股钢轨的轨距点上;在检测梁[16]内部还装有一个轨向惯性测量组件HIMU[13];在检测平台[15]的上表面安装有一个轨向-里程预处理计算机[4],轨向惯性测量组件HIMU[13]的输出端通过RS422总线与轨向-里程预处理计算机[4]的串口连接;所说的第一位移器LPDT[9]和第二位移器RPDT[10]是激光位移器;
(2)在检测平台[15]的上表面安装有一套高低测量系统,它由左惯性测量组件LIMU[7]、右惯性测量组件RIMU[8]和高低预处理计算机[2]组成;左惯性测量组件LIMU[7]和右惯性测量组件RIMU[8]位于轮轴[17]左右车轮的正上方,左惯性测量组件LIMU[7]和右惯性测量组件RIMU[8]的垂向轴线分别与所对应的左轨道、右轨道横截面的垂直中心线重合;左惯性测量组件LIMU[7]的输出端通过RS422总线与高低预处理计算机[2]的第一串口Com1连接;右惯性测量组件RIMU[8]的输出端通过RS422总线与高低预处理计算机[2]的第二串口Com2连接;
(3)第一激光位移器LPDT[9]、第二激光位移器RPDT[10]、第三激光位移器GLPDT[11]和第四激光位移器GRPDT[12]的输出端分别与位移器数据采集盒[5]的相应输入端连接;在检测平台[15]的上表面安装有一个位移器预处理计算机[3],位移器数据采集盒[5]的输出端通过USB总线与位移器预处理计算机[3]的USB口连接;
(4)高低预处理计算机[2]、位移预处理计算机[3]和轨向-里程预处理计算机[4]通过局域网与计算机[1]连接。
2.使用如权利要求1所述的高铁轨道平顺性检测装置检测铁路轨道平顺性的方法,其特征在于,检测的步骤如下:
2.1、传感器数据采集存储的步骤:
2.1.1、高低数据采集存储:系统上电后,左惯性测量组件LIMU[7]和右惯性测量组件[8]开始工作,高低预处理计算机[2]每隔m秒读取第一串口com1和第二串口com2的数据,将数据记为ωxL,ωyL,ωzL,AxL,AyL,AzL,ωxR,ωyR,ωzR,AxR,AyR,AzR,并进行保存;其中,ωxL表示左惯性测量组件LIMU[7]的横滚角速率,ωyL表示左惯性测量组件LIMU[7]的航向角速率,ωzL表示左惯性测量组件LIMU[7]的俯仰角速率,AxL,AyL,AzL分别表示左惯性测量组件LIMU[7]的横向、纵向、天向加速度;ωxR表示右惯性测量组件RIMU[8]的横滚角速率,ωyR表示右惯性测量组件RIMU[8]的航向角速率,ωzR表示右惯性测量组件RIMU[8]的俯仰角速率,AxR,AyR,AzR分别表示右惯性测量组件RIMU[8]的横向、纵向、天向加速度;m的取值范围是0.001秒~0.01秒;
2.1.2、轨向数据采集存储:系统上电后,轨向惯性测量组件HIMU[13]开始工作,每隔m秒,轨向惯性测量组件HIMU[13]向轨向-里程预处理计算机[4]第一串口com1发送采样数据,记为ωxH,ωyH,ωzH,AxH,AyH,AzH;其中ωxH表示轨向惯性测量组件HIMU[13]的横滚角速率,其中ωyH表示轨向惯性测量组件HIMU[13]的航向角速率,其中ωzH表示轨向惯性测量组件HIMU[13]的俯仰角速率,AxH,AyH,AzH分别表示轨向惯性测量组件HIMU[13]的横向、纵向、天向加速度,并进行保存;
2.1.3、里程数据采集存储:系统上电后,里程测量组件OD[14]开始工作,每隔m秒,里程测量组件OD[14]向轨向-里程预处理计算机[4]的第二串口com2发送一个采样数据,轨向-里程预处理计算机[4]将采样数据记为Nod并进行保存;
2.1.4、激光位移器采集存储:系统上电后,激光位移器数据控制盒[5]开始工作,每隔m秒,激光位移器数据控制盒[5]向位移器预处理计算机[3]发送一个采样数据包,每个采样数据包有4个采样数据,记为DWL、DWR、DGL、DGR,其中,DWL为第一激光位移器LPDT[9]所测量位移的采样数据;DWR为第二激光位移器RPDT[10]所测量位移的采样数据;DGL为第三激光位移器GLPDT[11]所测量位移的采样数据;DGR为第四激光位移器GRPDT[12]所测量位移的采样数据,并进行保存;
2.2、数据预处理的步骤:
2.2.1、高低预处理:高低预处理计算机[2]对数据的预处理过程是:
2.2.1.1、左惯性测量组件LIMU[7]数据预处理:
2.2.1.1.1、计算左惯性测量组件LIMU[7]的当前姿态:保持轨道车静止m1秒,m1取值范围为100秒~200秒,在m1秒内的采样个数为n,n=m1/m取整数;
根据下式计算出左惯性测量组件LIMU[7]的横滚角速率均值
Figure FSA00000470733300031
ω ‾ xL = 1 n Σ i = 1 n ω xL ( i ) · · · [ 1 ]
式中,i为采样数据的序号;
计算航向角初值ψ0
式中,
Figure FSA00000470733300034
表示当前位置的纬度;
2.2.1.1.2、地球自转和常值漂移补偿:
Figure FSA00000470733300035
Figure FSA00000470733300036
式中,
Figure FSA00000470733300037
为补偿后的左惯性测量组件LIMU[7]的横滚角速率,
Figure FSA00000470733300038
为补偿后的左惯性测量组件LIMU[7]的航向角速率,δψ为航向角修正量,当i=1时,令δψ为0;
2.2.1.1.3、积分计算:
δψ ( i ) = δψ ( i - 1 ) + ω ~ zL ( i ) · m · π / 180 · · · [ 5 ]
θ L ( i ) = θ L ( i - 1 ) + ω ~ xL ( i ) · m · π / 180 · · · [ 6 ]
式中,θL(i)为左惯性测量组件LIMU[7]的俯仰角;
2.2.1.2右惯性测量组件RIMU[8]数据预处理:
2.2.1.2.1、计算右惯性测量组件RIMU[8]的当前姿态:保持轨道车静止m1秒;根据下式计算出右惯性测量组件RIMU[8]横滚角速率均值
Figure FSA000004707333000311
ω ‾ xR = 1 n Σ i = 1 n ω xR ( i ) · · · [ 7 ]
2.2.1.2.2、地球自转和常值漂移补偿:
ω ~ xR ( i ) = ω xR ( i ) - 15 · cos L · cos ( ψ 0 + δψ ( i - 1 ) ) · · · [ 8 ]
式中,为补偿后的右惯性测量组件RIMU[8]横滚角速率;
2.2.1.1.3、积分计算:
θ R ( i ) = θ R ( i - 1 ) + ω ~ xR ( i ) · m · π / 180 · · · [ 10 ]
式中,θR(i)为右惯性测量组件RIMU[8]的俯仰角;
2.2.1.4、数据传输:高低预处理计算机[2]每m秒通过局域网向计算机[1]发送一次数据θL(i)和θR(i);
2.2.2、轨向-里程数据预处理:
2.2.2.1、计算轨向惯性测量组件HIMU[13]的当前航向角:保持轨道车静止m1秒,计算航向角速率均值
Figure FSA00000470733300044
ω ‾ xH = 1 n Σ i = 1 n ω xH ( i ) · · · [ 11 ]
航向角初值ψ0H计算:
Figure FSA00000470733300046
2.2.2.2、地球自转和常值漂移补偿:
式中,
Figure FSA00000470733300048
为补偿后的航向角速率,当i=1时,令ψH(i)=ψ0H为0;
2.2.2.3、积分计算
ψ H ( i ) = ψ H ( i - 1 ) + ω ~ zH ( i ) · m · π / 180 · · · [ 14 ]
式中,ψH(i)为轨向惯性测量组件HIMU[13]的航向角;
2.2.2.4、数据传输:
轨向-里程预处理计算机[4]每m秒通过局域网向计算机[1]发送一次ψH(i)和Nod(i);
2.2.3、激光位移器数据预处理:
2.2.3.1、计算平均值:分别计算DWL、DWR、DGL、DGR四个采样数据中每个采样数据每连续5个数的平均值
Figure FSA000004707333000410
其中,
Figure FSA000004707333000411
为第一激光位移器LPDT[9]所测量位移的采样数据平均值;
Figure FSA000004707333000412
为第二激光位移器RPDT[10]所测量位移的采样数据平均值;
Figure FSA000004707333000413
为第三激光位移器GLPDT[11]所测量位移的采样数据平均值;
Figure FSA000004707333000414
为第四激光位移器GRPDT[12]所测量位移的采样数据平均值;
2.2.3.2、数据传输:位移预处理计算机[3]每m秒通过局域网向计算机[1]发送一次
Figure FSA00000470733300051
2.3、综合数据采集处理:
2.3.1、综合数据采集:计算机[1]通过局域网获取高低预处理计算机[2]、位移预处理计算机[3]和轨向-里程预处理计算机[4]的数据并存储;
2.3.2、里程计算:根据下式计算当前里程:
里程仪的采样数据和事先标定好的刻度系数,
Lod(i)=Nod(i)·kod…………………………………………………[15]
式中,Lod(i)表示检测里程,kod是刻度系数,
Figure FSA00000470733300052
R是轨道车的车轮半径,其中p的取值范围是1024~4000;
2.3.3、时域-空域转换:
将步骤2.3.1计算机[1]接收到的时域采样数据,根据步骤2.3.2得到的里程,转换成采样密度为0.125米的空域数据,记为ΘL、ΘR、ΨH、L、
Figure FSA00000470733300053
其中ΘL为空域的左惯性测量组件LIMU[7]的俯仰角、ΘR为空域的右惯性测量组件RIMU[8]的俯仰角;ΨH为空域的轨向惯性测量组件HIMU[13]航向角,L为空域的里程;
Figure FSA00000470733300055
为空域的第一激光位移器LPDT[9]所测量位移的采样数据平均值;
Figure FSA00000470733300056
为空域的第二激光位移器RPDT[10]所测量位移的采样数据平均值;
Figure FSA00000470733300057
为空域的第三激光位移器GLPDT[11]所测量位移的采样数据平均值;
Figure FSA00000470733300058
为空域的第四激光位移器GRPDT[12]所测量位移的采样数据平均值,j为空域数据的序号;
2.3.4、当前波长平顺性解算:
2.3.4.1、波长输入:用户输入的波长λ1为数值10m、30m、42m、70m、300m其中之一;
2.3.4.2、平顺性解算:
2.3.4.2.1、高低平顺性解算:
h VL ( j ) = α [ dL 1 · tan ( dΘ 1 ) - 1 2 dL 2 · tan ( dΘ 2 ) ] + D ^ WL ( j ) · · · [ 16 ]
h VR ( j ) = α [ dL 1 · tan ( dΘ 3 ) - 1 2 dL 2 · tan ( dΘ 4 ) ] + D ^ WR ( j ) · · · [ 17 ]
式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
dL1=L(j)-L(j-d)……………………………………………………… [18]
dL2=L(j+d)-L(j-d)…………………………………………………… [19]
1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
d = λ 1 0.25 · · · [ 24 ]
2.3.4.2.2、轨向平顺性解算:
S L ( j ) = y ( j ) + D ^ GL ( j ) · · · [ 25 ]
S R ( j ) = y ( j ) + D ^ GR ( j ) · · · [ 26 ]
y ( j ) = β · [ dL 1 tan ( dψ 1 ) - 1 2 dL 2 tan ( dψ 2 ) ] · · · [ 27 ]
式中,y(j)为检测梁[16]的中心横向位移,β为幅值补偿系数,取0.03;
1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
至此,得到了轨道平顺性的检测结果。
CN201110089802A 2011-04-12 2011-04-12 一种多波长高铁轨道平顺性检测装置和检测方法 Active CN102252636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110089802A CN102252636B (zh) 2011-04-12 2011-04-12 一种多波长高铁轨道平顺性检测装置和检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110089802A CN102252636B (zh) 2011-04-12 2011-04-12 一种多波长高铁轨道平顺性检测装置和检测方法

Publications (2)

Publication Number Publication Date
CN102252636A true CN102252636A (zh) 2011-11-23
CN102252636B CN102252636B (zh) 2012-09-05

Family

ID=44980057

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110089802A Active CN102252636B (zh) 2011-04-12 2011-04-12 一种多波长高铁轨道平顺性检测装置和检测方法

Country Status (1)

Country Link
CN (1) CN102252636B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103322948A (zh) * 2013-07-04 2013-09-25 上海宝锋工程技术有限公司 板坯表面凹凸点检测方法、检测装置与处理方法
CN104032629A (zh) * 2014-06-10 2014-09-10 南京理工大学 一种垂向轨道长波不平顺在线监测方法及系统
CN104358194A (zh) * 2014-10-27 2015-02-18 同济大学 基于激光跟踪仪的轨道静态平顺性测量与分析方法
CN104819851A (zh) * 2015-04-27 2015-08-05 上海海马汽车研发有限公司 一种汽车平顺性试验数据的处理方法及系统
CN105807084A (zh) * 2015-12-22 2016-07-27 成都安科泰丰科技有限公司 轨道几何参数检测同步及计算方法
CN106524924A (zh) * 2016-11-25 2017-03-22 闫东坤 一种基于光学惯性组合的轨检车位置姿态测量系统及方法
CN106643445A (zh) * 2016-12-30 2017-05-10 亿嘉和科技股份有限公司 一种轨道平整度测量方法
CN107200040A (zh) * 2016-03-17 2017-09-26 斯凯孚公司 用于确定轨道表面的垂直轮廓的方法和系统
CN108828068A (zh) * 2018-01-18 2018-11-16 武汉汉宁轨道交通技术有限公司 轨道表面波浪磨耗及伤损检测装置及方法
CN110274584A (zh) * 2019-08-08 2019-09-24 中建空列(北京)科技有限公司 空铁立柱测量装置
CN110530300A (zh) * 2019-09-05 2019-12-03 西南交通大学 一种轨道垂向不平顺测量方法
CN110954009A (zh) * 2019-12-20 2020-04-03 逸美德科技股份有限公司 轮毂端面变形检测方法及其装置
CN111380512A (zh) * 2018-12-28 2020-07-07 中国航空工业集团公司西安飞行自动控制研究所 基于惯性技术的高铁轨道长波测量方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07262717A (ja) * 1994-03-22 1995-10-13 Kenwood Corp 光ディスクプレーヤ
JP2595318B2 (ja) * 1988-07-11 1997-04-02 株式会社日立製作所 書替え可能形記録再生装置のデ−タ消去装置
CN101178303A (zh) * 2007-11-26 2008-05-14 刘洪云 利用激光检测铁路路轨的方法及装置
CN201575794U (zh) * 2009-11-27 2010-09-08 刘洪云 高速铁路轨道长波不平顺精密测量仪

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2595318B2 (ja) * 1988-07-11 1997-04-02 株式会社日立製作所 書替え可能形記録再生装置のデ−タ消去装置
JPH07262717A (ja) * 1994-03-22 1995-10-13 Kenwood Corp 光ディスクプレーヤ
CN101178303A (zh) * 2007-11-26 2008-05-14 刘洪云 利用激光检测铁路路轨的方法及装置
CN201575794U (zh) * 2009-11-27 2010-09-08 刘洪云 高速铁路轨道长波不平顺精密测量仪

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103322948A (zh) * 2013-07-04 2013-09-25 上海宝锋工程技术有限公司 板坯表面凹凸点检测方法、检测装置与处理方法
CN103322948B (zh) * 2013-07-04 2015-12-23 上海宝锋工程技术有限公司 板坯表面凹凸点检测方法、检测装置与处理方法
CN104032629A (zh) * 2014-06-10 2014-09-10 南京理工大学 一种垂向轨道长波不平顺在线监测方法及系统
CN104032629B (zh) * 2014-06-10 2016-03-02 南京理工大学 一种垂向轨道长波不平顺在线监测方法及系统
CN104358194A (zh) * 2014-10-27 2015-02-18 同济大学 基于激光跟踪仪的轨道静态平顺性测量与分析方法
CN104358194B (zh) * 2014-10-27 2016-08-24 同济大学 基于激光跟踪仪的轨道静态平顺性测量与分析方法
CN104819851A (zh) * 2015-04-27 2015-08-05 上海海马汽车研发有限公司 一种汽车平顺性试验数据的处理方法及系统
CN104819851B (zh) * 2015-04-27 2017-11-14 上海海马汽车研发有限公司 一种汽车平顺性试验数据的处理方法及系统
CN105807084A (zh) * 2015-12-22 2016-07-27 成都安科泰丰科技有限公司 轨道几何参数检测同步及计算方法
CN107200040A (zh) * 2016-03-17 2017-09-26 斯凯孚公司 用于确定轨道表面的垂直轮廓的方法和系统
CN107200040B (zh) * 2016-03-17 2022-10-28 斯凯孚公司 用于确定轨道表面的垂直轮廓的方法和系统
CN106524924A (zh) * 2016-11-25 2017-03-22 闫东坤 一种基于光学惯性组合的轨检车位置姿态测量系统及方法
CN106643445A (zh) * 2016-12-30 2017-05-10 亿嘉和科技股份有限公司 一种轨道平整度测量方法
CN106643445B (zh) * 2016-12-30 2019-04-16 亿嘉和科技股份有限公司 一种轨道平整度测量方法
CN108828068A (zh) * 2018-01-18 2018-11-16 武汉汉宁轨道交通技术有限公司 轨道表面波浪磨耗及伤损检测装置及方法
CN111380512A (zh) * 2018-12-28 2020-07-07 中国航空工业集团公司西安飞行自动控制研究所 基于惯性技术的高铁轨道长波测量方法
CN110274584A (zh) * 2019-08-08 2019-09-24 中建空列(北京)科技有限公司 空铁立柱测量装置
CN110274584B (zh) * 2019-08-08 2024-04-16 中建空列(北京)科技有限公司 空铁立柱测量装置
CN110530300A (zh) * 2019-09-05 2019-12-03 西南交通大学 一种轨道垂向不平顺测量方法
CN110954009A (zh) * 2019-12-20 2020-04-03 逸美德科技股份有限公司 轮毂端面变形检测方法及其装置
CN110954009B (zh) * 2019-12-20 2021-09-07 逸美德科技股份有限公司 轮毂端面变形检测方法及其装置

Also Published As

Publication number Publication date
CN102252636B (zh) 2012-09-05

Similar Documents

Publication Publication Date Title
CN102252636A (zh) 一种多波长高铁轨道平顺性检测装置和检测方法
CN102180187B (zh) 一种铁路轨道高低高精度检测装置和检测方法
CN102251451B (zh) 基于多源信息融合的轨道几何状态测量方法
CN103207403B (zh) 一种卫星导航与惯性测量组合轨道测量系统与方法
CN106274981B (zh) 一种轨道检测装置及检测方法
CN202644332U (zh) 轨道顶面不平顺检测装置
CN102252627B (zh) 一种高铁轨道轨距检测装置和检测方法
CN107244335B (zh) 一种波浪磨耗检测装置及方法
CN202368604U (zh) 一种动态检测铁路轨道左右钢轨高低的检测装置
CN101804822A (zh) 激光测量仪及测量方法
KR101358870B1 (ko) 궤도 레일의 종구배와 횡구배 및 레일 궤간 자동 측정 장치
CN105235713A (zh) 基于激光位移传感器的城轨车辆车轮直径在线检测方法
CN108755306A (zh) 一种轨道几何状态高精度测量装置
CN102337710A (zh) 一种gps轨道不平顺检测系统及其检测方法
CN107815935A (zh) 高速铁路轨道几何状态实时监测方法及系统
CN201746752U (zh) 轨道检测装置
CN204286406U (zh) 使用mems陀螺仪的便携式轨道几何量检测系统
CN101761004A (zh) 用于中低速磁悬浮f型轨的轨道检测仪
CN102175168B (zh) 中低速磁悬浮f型轨的轨道动态检测设备
CN104047212B (zh) 一种基于角度量测的轨道沉降自动测量装置及方法
CN205246128U (zh) 轨道几何参数在线检测装置
CN104567783A (zh) 使用mems陀螺仪的便携式轨道几何量检测系统及方法
Tanaka et al. Practical application of portable trolley for the continuous measurement of rail surface roughness for rail corrugation maintenance
CN113415308B (zh) 钢轨波磨检测方法及装置
CN105783777B (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