技术领域
[0001] 本
发明属于
铁路轨道检测技术,涉及对铁路轨道高低、轨向等平顺性检测装置和方法的改进。
背景技术
[0002] 目前,对铁路轨道平顺性的检测方法主要有弦测法、惯性基准法等。其中弦测法是指采用人工拉线的方法在轨道上真实搭建一条弦线,通过测量轨道轨顶面(轨距点)与该弦线的相对位移,评价轨道的高低平顺性。通过这种方法测得的轨道平顺性,波长与检测弦的长度密切相关,当需要分析多种波长的轨道平顺性时,需要更换检测装置重新测量,或者根据当前检测结果进行“以小推大”的转换,前者增加了检测人员的工作量,工作效率低,后者存在较大误差。惯性基准法是通过设计高阶
滤波器实现对不同范围的轨道平顺性进行分析,目前技术较为成熟的可以分析70m以内的轨道不平顺。这种方法不能满足高铁对100
米以上的长波不平顺测量需求。惯性基准法可参见《GJ-4型轨道检查车的原理与应用》,中国铁道出版社,2000,张未。
发明内容
[0003] 本发明的目的是:提出一种能对轨道的高低、轨向进行多种波长连续检测的检测装置和检测方法,并能满足高铁对100米以上的长波不平顺测量需求。
[0004] 本发明的技术方案是:一种高铁轨道平顺性检测装置,包括一个带有检测平台的轨道检测车,在检测平台的下表面固定有两个位移器,它们是第一位移器LPDT和第二位移器RPDT,第一位移器LPDT和第二位移器RPDT分别位于轮轴两端的正上方,第一位移器
LPDT和第二位移器RPDT的轴线与轮轴的轴线
正交,在轮轴的一端安装有里程仪OD,在检测平台的上表面安装有
数据处理系统,该
数据处理系统包括计算机、位移器
数据采集控制盒、OD数据采集卡和电源;其特征在于:
[0005] (1)在检测平台下表面固定有一个截面为矩形的、中空的检测梁,检测梁与轮轴平行,检测梁横向的垂直平分面与轮轴横向的垂直平分面共面;检测梁的谐振
频率不小于70Hz;在检测梁内部安装有两个位移器,分别是第三激光位移器GLPDT和第四激光位移器GRPDT,在检测梁的下表面、与第三激光位移器GLPDT和第四激光位移器GRPDT对应的
位置各有一个通光孔,第三激光位移器GLPDT和第四激光位移器GRPDT发射的
激光束分别穿过
上述通光孔后照射到左右股
钢轨的轨距点上;在检测梁内部还装有一个轨向惯性测量组件HIMU;在检测平台的上表面安装有一个轨向-里程预处理计算机,轨向惯性测量组件HIMU的输出端通过RS422总线与轨向-里程预处理计算机的串口连接;所说的第一位移器LPDT
和第二位移器RPDT是激光位移器;
[0006] (2)在检测平台的上表面安装有一套高低测量系统,它由左惯性测量组件LIMU、右惯性测量组件RIMU和高低预处理计算机组成;左惯性测量组件LIMU和右惯性测量组件RIMU位于轮轴左右
车轮的正上方,左惯性测量组件LIMU和右惯性测量组件RIMU的垂向轴
线分别与所对应的左轨道、右轨道横截面的垂直中心线重合;左惯性测量组件LIMU的输出端通过RS422总线与高低预处理计算机的第一串口Com1连接;右惯性测量组件RIMU的输
出端通过RS422总线与高低预处理计算机的第二串口Com2连接;
[0007] (3)第一激光位移器LPDT、第二激光位移器RPDT、第三激光位移器GLPDT和第四激光位移器GRPDT的输出端分别与位移器数据采集盒的相应输入端连接;在检测平台的上表面安装有一个位移器预处理计算机,位移器数据采集盒的输出端通过USB总线与位移器预处理计算机的USB口连接;
[0008] (4)高低预处理计算机、位移预处理计算机和轨向-里程预处理计算机通过局域网与计算机连接。
[0009] 使用如上面所述的高铁轨道平顺性检测装置检测铁路轨道平顺性的方法,其特征在于,检测的步骤如下:
[0011] 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秒;
[0012] 1.2、轨向数据采集存储:系统上电后,轨向惯性测量组件HIMU开始工作,每隔m秒,轨向惯性测量组件HIMU向轨向-里程预处理计算机第一串口com1发送
采样数据,记为ωxH,ωyH,ωzH,AxH,AyH,AzH;其中ωxH表示轨向惯性测量组件HIMU的横滚角速率,其中ωyH表示轨向惯性测量组件HIMU的航向角速率,其中ωzH表示轨向惯性测量组件HIMU的俯仰角速率,AxH,AyH,AzH分别表示轨向惯性测量组件HIMU的横向、纵向、天向加速度,并进行保存;
[0013] 1.3、里程数据采集存储:系统上电后,里程测量组件OD开始工作,每隔m秒,里程测量组件OD向轨向-里程预处理计算机的第二串口com2发送一个采样数据,轨向-里程预处理计算机将采样数据记为Nod并进行保存;
[0014] 1.4、激光位移器采集存储:系统上电后,激光位移器数据控制盒开始工作,每隔m秒,激光位移器数据控制盒向位移器预处理计算机发送一个采样数据包,每个采样数据包有4个采样数据,记为DWL、DWR、DGL、DGR,其中,DWL为第一激光位移器LPDT所测量位移的采样数据;DWR为第二激光位移器RPDT所测量位移的采样数据;DGL为第三激光位移器GLPDT所测量位移的采样数据;DGR为第四激光位移器GRPDT所测量位移的采样数据,并进行保存;
[0015] 2、数据预处理的步骤:
[0016] 2.1、高低预处理:高低预处理计算机对数据的预处理过程是:
[0017] 2.1.1、左惯性测量组件LIMU数据预处理:
[0018] 2.1.1.1、计算左惯性测量组件LIMU的当前
姿态:保持轨道车静止m1秒,m1取值范围为100秒~200秒,在m1秒内的采样个数为n,n=m1/m取整数;
[0019] 根据下式计算出左惯性测量组件LIMU的横滚角速率均值
[0020]
[0021] 式中,i为采样数据的序号;
[0022] 计算航向角初值ψ0:
[0023]
[0024] 式中,表示当前位置的纬度;
[0025] 2.1.1.2、地球自转和常值漂移补偿:
[0026]
[0027]
[0028] 式中, 为补偿后的左惯性测量组件LIMU的横滚角速率, 为补偿后的左惯性测量组件LIMU的航向角速率,δψ为航向角修正量,当i=1时,令δψ为0;
[0029] 1.1.3、积分计算:
[0030]
[0031]
[0032] 式中,θL(i)为左惯性测量组件LIMU的俯仰角;
[0033] 2.1.2右惯性测量组件RIMU数据预处理:
[0034] 2.1.2.1、计算右惯性测量组件RIMU的当前姿态:保持轨道车静止m1秒;根据下式计算出右惯性测量组件RIMU横滚角速率均值
[0035]
[0036] 2.1.2.2、地球自转和常值漂移补偿:
[0037]
[0038] 式中, 为补偿后的右惯性测量组件RIMU横滚角速率;
[0039] 2.1.1.3、积分计算:
[0040]
[0041] 式中,θR(i)为右惯性测量组件RIMU的俯仰角;
[0042] 2.1.4、数据传输:高低预处理计算机每m秒通过局域网向计算机发送一次数据θR(i)和θR(i);
[0043] 2.2、轨向-里程数据预处理:
[0044] 2.2.1、计算轨向惯性测量组件HIMU的当前航向角:保持轨道车静止m1秒,计算航向角速率均值
[0045]
[0046] 航向角初值ψ0H计算:
[0047]
[0048] 2.2.2、地球自转和常值漂移补偿:
[0049]
[0050] 式中, 为补偿后的航向角速率,当i=1时,令ψH(i)=ψ0H为0;
[0051] 2.2.3、积分计算
[0052]
[0053] 式中,ψH(i)为轨向惯性测量组件HIMU的航向角;
[0054] 2.2.4、数据传输:
[0055] 轨向-里程预处理计算机每m秒通过局域网向计算机发送一次ψH(i)和Nod(i);
[0056] 2.3、激光位移器数据预处理:
[0057] 2.3.1、计算平均值:分别计算DWL、DWH、DGL、DGR四个采样数据中每个采样数据每连续5个数的平均值 其中, 为第一激光位移器LPDT所测量位移的采样数据平均值; 为第二激光位移器RPDT所测量位移的采样数据平均
值; 为第三激光位移器GLPDT所测量位移的采样数据平均值; 为第四激光位移
器GRPDT所测量位移的采样数据平均值;
[0058] 2.3.2、数据传输:位移预处理计算机每m秒通过局域网向计算机发送一次
[0059] 3、综合数据采集处理:
[0060] 3.1、综合数据采集:计算机通过局域网获取高低预处理计算机、位移预处理计算机和轨向-里程预处理计算机的数据并存储;
[0061] 3.2、
里程计算:根据下式计算当前里程:
[0062] 里程仪的采样数据和事先标定好的刻度系数,
[0063] Lod(i)=Nod(i)·kod…………………………………………………[15]
[0064] 式中,Lod(i)表示检测里程,kod是刻度系数, R是轨道车的车轮半径,其中p的取值范围是1024~4000;
[0066] 将步骤3.1计算机接收到的时域采样数据,根据步骤3.2得到的里程,转换成采样
密度为0.125米的空域数据,记为ΘL、ΘR、ΨH、L、 其中ΘL为空域的左惯性测量组件LIMU的俯仰角、ΘR为空域的右惯性测量组件RIMU的俯仰角;ΨH为空域的轨向惯性测量组件HIMU航向角,L为空域的里程; 为空域的第一激光位移
器LPDT所测量位移的采样数据平均值; 为空域的第二激光位移器RPDT所测量位移
的采样数据平均值; 为空域的第三激光位移器GLPDT所测量位移的采样数据平均值;
为空域的第四激光位移器GRPDT所测量位移的采样数据平均值,j为空域数据的序
号;
[0067] 3.4、当前波长平顺性解算:
[0068] 3.4.1、波长输入:用户输入的波长λ1为数值10m、30m、42m、70m、300m其中之一;
[0069] 3.4.2、平顺性解算:
[0070] 3.4.2.1、高低平顺性解算:
[0071]
[0072]
[0073] 式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
[0074] dL1=L(j)-L(j-d)……………………………………………………… [18]
[0075] dL2=L(j+d)-L(j-d)…………………………………………………… [19]
[0076] dΘ1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
[0077] dΘ2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
[0078] dΘ3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
[0079] dΘ4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
[0080]
[0081] 3.4.2.2、轨向平顺性解算:
[0082]
[0083]
[0084]
[0085] 式中,y(j)为检测梁[16]的中心
横向位移,β为幅值补偿系数,取0.03;
[0086] dψ1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
[0087] dψ2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
[0088] 至此,得到了轨道平顺性的检测结果。
[0089] 本发明的优点是:使用本发明的检测装置和检测方法可以通过一次检测和多次解算,就可以得对10m、30m、42m、70m或300m多种波长的轨道不平顺进行检测,大大缩短了轨道上的工作时间,提高了检测效率;并满足了高铁对100米以上的长波不平顺的测量需求。
附图说明
[0090] 图1是本发明检测装置的结构原理
框图。
[0091] 图2是本发明检测装置中传感器安装位置示意图。
具体实施方式
[0092] 下面对本发明做进一步详细说明。规定方位如下,面对轨道检测车的前进方向,左手边为左方。参见图1、2。一种高铁轨道平顺性检测装置,包括一个带有检测平台15的轨道检测车,在检测平台15的下表面固定有两个位移器,它们是第一位移器LPDT9和第二位移器RPDT10,第一位移器LPDT9和第二位移器RPDT10分别位于轮轴17两端的正上方,第一位移器LPDT9和第二位移器RPDT10的轴线与轮轴17的轴线正交,在轮轴17的一端安装有里程仪OD14,在检测平台15的上表面安装有数据处理系统,该数据处理系统包括计算机1、位移器数据采集控制盒5、OD数据采集卡6和电源;其特征在于:
[0093] (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是激光位移
器;
[0094] (2)在检测平台15的上表面安装有一套高低测量系统,它由左惯性测量组件LIMU7、右惯性测量组件RIMU8和高低预处理计算机2组成;左惯性测量组件LIMU7和右惯
性测量组件RIMU8位于轮轴17左右车轮的正上方,左惯性测量组件LIMU7和右惯性测量组
件RIMU8的垂向轴线分别与所对应的左轨道、右轨道横截面的垂直中心线重合;左惯性测量组件LIMU7的输出端通过RS422总线与高低预处理计算机2的第一串口Com1连接;右惯
性测量组件RIMU8的输出端通过RS422总线与高低预处理计算机2的第二串口Com2连接;
[0095] (3)第一激光位移器LPDT9、第二激光位移器RPDT10、第三激光位移器GLPDT11和第四激光位移器GRPDT12的输出端分别与位移器数据采集盒5的相应输入端连接;在检测平台15的上表面安装有一个位移器预处理计算机3,位移器数据采集盒5的输出端通过
USB总线与位移器预处理计算机3的USB口连接;
[0096] (4)高低预处理计算机2、位移预处理计算机3和轨向-里程预处理计算机4通过局域网与计算机1连接。
[0097] 本发明检测装置的检测原理是:利用3个光纤惯组分别测量轨道车的左、右端俯仰角和检测梁的航向角,结合里程和检测波长计算得到车体的高低和检测梁中心位移,结合第一激光位移器和第二激光位移器的测量结果,扣除车体相对于轮轴的振动量,计算得到左、右股轨道高低平顺性测量结果,结合第三激光位移器和第四激光位移器的测量结果,加上对应的左、右轨距点位移,计算得到左、右股轨道轨向平顺性测量结果。
[0098] 使用如上面所述的高铁轨道平顺性检测装置检测铁路轨道平顺性的方法,其特征在于,检测的步骤如下:
[0100] 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秒;
[0101] 1.2、轨向数据采集存储:系统上电后,轨向惯性测量组件HIMU13开始工作,每隔m秒,轨向惯性测量组件HIMU13向轨向-里程预处理计算机4第一串口com1发送采样数据,记为ωxH,ωyH,ωzH,AxH,AyH,AzH;其中ωxH表示轨向惯性测量组件HIMU13的横滚角速率,其中ωyH表示轨向惯性测量组件HIMU13的航向角速率,其中ωzH表示轨向惯性测量组件
HIMU13的俯仰角速率,AxH,AyH,AzH分别表示轨向惯性测量组件HIMU13的横向、纵向、天向加速度,并进行保存;
[0102] 1.3、里程数据采集存储:系统上电后,里程测量组件OD14开始工作,每隔m秒,里程测量组件OD14向轨向-里程预处理计算机4的第二串口com2发送一个采样数据,轨向-里程预处理计算机4将采样数据记为Nod并进行保存;
[0103] 1.4、激光位移器采集存储:系统上电后,激光位移器数据控制盒5开始工作,每隔m秒,激光位移器数据控制盒5向位移器预处理计算机3发送一一个采样数据包,每个采样数据包有4个采样数据,记为DWL、DWR、DGL、DGR,其中,DWL为第一激光位移器LPDT9所测量位移的采样数据;DWR为第二激光位移器RPDT10所测量位移的采样数据;DGL为第三激光位移器GLPDT11所测量位移的采样数据;DGR为第四激光位移器GRPDT12所测量位移的采样数据,并进行保存;
[0104] 2、数据预处理的步骤:
[0105] 2.1、高低预处理:高低预处理计算机2对数据的预处理过程是:
[0106] 2.1.1、左惯性测量组件LIMU7数据预处理:
[0107] 2.1.1.1、计算左惯性测量组件LIMU7的当前姿态:保持轨道车静止m1秒,m1取值范围为100秒~200秒,在m1秒内的采样个数为n,n=m1/m取整数;
[0108] 根据下式计算出左惯性测量组件LIMU7的横滚角速率均值
[0109]
[0110] 式中,i为采样数据的序号;
[0111] 计算航向角初值ψ0:
[0112]
[0113] 式中,表示当前位置的纬度;
[0114] 2.1.1.2、地球自转和常值漂移补偿:
[0115]
[0116]
[0117] 式中, 为补偿后的左惯性测量组件LIMU7的横滚角速率, 为补偿后的左惯性测量组件LIMU7的航向角速率,δψ为航向角修正量,当i=1时,令δψ为0;
[0118] 1.1.3、积分计算:
[0119]
[0120]
[0121] 式中,θL(i)为左惯性测量组件LIMU7的俯仰角;
[0122] 2.1.2右惯性测量组件RIMU8数据预处理:
[0123] 2.1.2.1、计算右惯性测量组件RIMU8的当前姿态:保持轨道车静止m1秒;根据下式计算出右惯性测量组件RIMU8横滚角速率均值
[0124]
[0125] 2.1.2.2、地球自转和常值漂移补偿:
[0126]
[0127] 式中, 为补偿后的右惯性测量组件RIMU8横滚角速率;
[0128] 2.1.1.3、积分计算:
[0129]
[0130] 式中,θR(i)为右惯性测量组件RIMU8的俯仰角;
[0131] 2.1.4、数据传输:高低预处理计算机2每m秒通过局域网向计算机1发送一次数据θL(i)和θR(i);
[0132] 2.2、轨向-里程数据预处理:
[0133] 2.2.1、计算轨向惯性测量组件HIMU13的当前航向角:保持轨道车静止m1秒,计算航向角速率均值
[0134]
[0135] 航向角初值ψ0H计算:
[0136]
[0137] 2.2.2、地球自转和常值漂移补偿:
[0138]
[0139] 式中, 为补偿后的航向角速率,当i=1时,令ψH(i)=ψ0H为0;
[0140] 2.2.3、积分计算
[0141]
[0142] 式中,ψH(i)为轨向惯性测量组件HIMU13的航向角;
[0143] 2.2.4、数据传输:
[0144] 轨向-里程预处理计算机4每m秒通过局域网向计算机1发送一次ψH(i)和Nod(i);
[0145] 2.3、激光位移器数据预处理:
[0146] 2.3.1、计算平均值:分别计算DWL、DWR、DGL、DGR四个采样数据中每个采样数据每连续5个数的平均值 其中, 为第一激光位移器LPDT9所测量位移的采样数据平均值; 为第二激光位移器RPDT10所测量位移的采样数据平
均值; 为第三激光位移器GLPDT11所测量位移的采样数据平均值; 为第四激光
位移器GRPDT12所测量位移的采样数据平均值;
[0147] 2.3.2、数据传输:位移预处理计算机3每m秒通过局域网向计算机1发送一次
[0148] 3、综合数据采集处理:
[0149] 3.1、综合数据采集:计算机1通过局域网获取高低预处理计算机2、位移预处理计算机3和轨向-里程预处理计算机4的数据并存储;
[0150] 3.2、里程计算:根据下式计算当前里程:
[0151] 里程仪的采样数据和事先标定好的刻度系数,
[0152] Lod(i)=Nod(i)·kod…………………………………………………[15]
[0153] 式中,Lod(i)表示检测里程,kod是刻度系数, R是轨道车的车轮半径,其中p的取值范围是1024~4000;
[0154] 3.3、时域-空域转换:
[0155] 将步骤3.1计算机1接收到的时域采样数据,根据步骤3.2得到的里程,转换成采样密度为0.125米的空域数据,记为ΘL、ΘR、ΨH、L、 其中ΘL为空域的左惯性测量组件LIMU7的俯仰角、ΘR为空域的右惯性测量组件RIMU8的俯仰
角;ΨH为空域的轨向惯性测量组件HIMU13航向角,L为空域的里程; 为空域的第一激
光位移器LPDT9所测量位移的采样数据平均值; 为空域的第二激光位移器RPDT10所
测量位移的采样数据平均值; 为空域的第三激光位移器GLPDT11所测量位移的采样
数据平均值; 为空域的第四激光位移器GRPDT12所测量位移的采样数据平均值,j为
空域数据的序号;
[0156] 3.4、当前波长平顺性解算:
[0157] 3.4.1、波长输入:用户输入的波长λ1为数值10m、30m、42m、70m、300m其中之一;
[0158] 3.4.2、平顺性解算:
[0159] 3.4.2.1、高低平顺性解算:
[0160]
[0161]
[0162] 式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
[0163] dL1=L(j)-L(j-d)……………………………………………………… [18]
[0164] dL2=L(j+d)-L(j-d)…………………………………………………… [19]
[0165] dΘ1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
[0166] dΘ2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
[0167] dΘ3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
[0168] dΘ4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
[0169]
[0170] 3.4.2.2、轨向平顺性解算:
[0171]
[0172]
[0173]
[0174] 式中,y(j)为检测梁[16]的中心横向位移,β为幅值补偿系数,取0.03;
[0175] dψ1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
[0176] dψ2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
[0177] 至此,得到了轨道平顺性的检测结果。
[0179] 在装有检测平台的轨道车上安装如上面所述的检测装置,其中传感器的技术指标分别为:左惯性测量组件LIMU7、右惯性测量组件RIMU8、轨向惯性测量组件HIMU13均选用光纤捷联惯组,第一位移器LPDT9、第二位移器RPDT10、第三激光位移器GLPDT11、第四激光位移器GRPDT 12均选用测量
精度为0.05mm的LK-400G激光位移器,OD14选用3600脉冲/转的光电
编码器。检测的步骤如下:
[0180] 1、传感器数据采集存储的步骤:
[0181] 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秒;
[0182] 1.2、轨向数据采集存储:系统上电后,轨向惯性测量组件HIMU13开始工作,每隔m秒,轨向惯性测量组件HIMU13向轨向-里程预处理计算机4第一串口com1发送采样数据,记为ωxH,ωyH,ωzH,AxH,AyH,AzH;其中ωxH表示轨向惯性测量组件HIMU13的横滚角速率,其中ωyH表示轨向惯性测量组件HIMU13的航向角速率,其中ωzH表示轨向惯性测量组件
HIMU13的俯仰角速率,AxH,AyH,AzH分别表示轨向惯性测量组件HIMU13的横向、纵向、天向加速度,并进行保存;
[0183] 1.3、里程数据采集存储:系统上电后,里程测量组件OD14开始工作,每隔m秒,里程测量组件OD14向轨向-里程预处理计算机4的第二串口com2发送一个采样数据,轨向-里程预处理计算机4将采样数据记为Nod并进行保存;
[0184] 1.4、激光位移器采集存储:系统上电后,激光位移器数据控制盒5开始工作,每隔m秒,激光位移器数据控制盒5向位移器预处理计算机3发送一次采样数据,每次有4个采样数据,记为DWL、DWR、DGL、DGR,其中,DWL为第一激光位移器LPDT9所测量位移的采样数据;DWR为第二激光位移器RPDT10所测量位移的采样数据;DGL为第三激光位移器GLPDT11所测量
位移的采样数据;DGR为第四激光位移器GRPDT12所测量位移的采样数据,并进行保存;
[0185] 2、数据预处理的步骤:
[0186] 2.1、高低预处理:高低预处理计算机2对数据的预处理过程是:
[0187] 2.1.1、左惯性测量组件LIMU7数据预处理:
[0188] 2.1.1.1、计算左惯性测量组件LIMU7的当前姿态:保持轨道车静止m1秒,m1取值范围为150秒,在m1秒内的采样个数为n,n=m1/m取整数;
[0189] 根据下式计算出左惯性测量组件LIMU7的横滚角速率均值
[0190]
[0191] 式中,i为采样数据的序号;
[0192] 计算航向角初值ψ0:
[0193]
[0194] 式中,表示当前位置的纬度;
[0195] 2.1.1.2、地球自转和常值漂移补偿:
[0196]
[0197]
[0198] 式中, 为补偿后的左惯性测量组件LIMU7的横滚角速率, 为补偿后的左惯性测量组件LIMU7的航向角速率,δψ为航向角修正量,当i=1时,令δψ为0;
[0199] 1.1.3、积分计算:
[0200]
[0201]
[0202] 式中,θL(i)为左惯性测量组件LIMU7的俯仰角;
[0203] 2.1.2右惯性测量组件RIMU8数据预处理:
[0204] 2.1.2.1、计算右惯性测量组件RIMU8的当前姿态:保持轨道车静止m1秒;根据下式计算出右惯性测量组件RIMU8横滚角速率均值
[0205]
[0206] 2.1.2.2、地球自转和常值漂移补偿:
[0207]
[0208] 式中, 为补偿后的右惯性测量组件RIMU8横滚角速率;
[0209] 2.1.1.3、积分计算:
[0210]
[0211] 式中,θR(i)为右惯性测量组件RIMU8的俯仰角;
[0212] 2.1.4、数据传输:高低预处理计算机2每m秒通过局域网向计算机1发送一次数据θL(i)和θR(i);
[0213] 2.2、轨向-里程数据预处理:
[0214] 2.2.1、计算轨向惯性测量组件HIMU13的当前航向角:保持轨道车静止m1秒,计算航向角速率均值
[0215]
[0216] 航向角初值ψ0H计算:
[0217]
[0218] 2.2.2、地球自转和常值漂移补偿:
[0219]
[0220] 式中, 为补偿后的航向角速率,当i=1时,令ψH(i)=ψ0H为0;
[0221] 2.2.3、积分计算
[0222]
[0223] 式中,ψH(i)为轨向惯性测量组件HIMU13的航向角;
[0224] 2.2.4、数据传输:
[0225] 轨向-里程预处理计算机4每m秒通过局域网向计算机1发送一次ψH(i)和Nod(i);
[0226] 2.3、激光位移器数据预处理:
[0227] 2.3.1、计算平均值:分别计算DWL、DWR、DGL、DGR四个采样数据中每个采样数据每连续5个数的平均值 其中, 为第一激光位移器LPDT9所测量位移的采样数据平均值; 为第二激光位移器RPDT10所测量位移的采样数据平
均值; 为第三激光位移器GLPDT11所测量位移的采样数据平均值; 为第四激光
位移器GRPDT12所测量位移的采样数据平均值;
[0228] 2.3.2、数据传输:位移预处理计算机3每m秒通过局域网向计算机1发送一次
[0229] 3、综合数据采集处理:
[0230] 3.1、综合数据采集:计算机1通过局域网获取高低预处理计算机2、位移预处理计算机3和轨向-里程预处理计算机4的数据并存储;
[0231] 3.2、里程计算:根据下式计算当前里程:
[0232] 里程仪的采样数据和事先标定好的刻度系数,
[0233] Lod(i)=Nod(i)·kod…………………………………………………[15]
[0234] 式中,Lod(i)表示检测里程,kod是刻度系数, R是轨道车的车轮半径,为0.915米,其中p的取值是3600;
[0235] 3.3、时域-空域转换:
[0236] 将步骤3.1计算机1接收到的时域采样数据,根据步骤3.2得到的里程,转换成采样密度为0.125米的空域数据,记为ΘL、ΘR、ΨH、L、 其中ΘL为空域的左惯性测量组件LIMU7的俯仰角、ΘR为空域的右惯性测量组件RIMU8的俯仰
角;ΨH为空域的轨向惯性测量组件HIMU13航向角,L为空域的里程; 为空域的第一激
光位移器LPDT9所测量位移的采样数据平均值; 为空域的第二激光位移器RPDT10所
测量位移的采样数据平均值; 为空域的第三激光位移器GLPDT11所测量位移的采样
数据平均值; 为空域的第四激光位移器GRPDT12所测量位移的采样数据平均值,j为
空域数据的序号;
[0237] 3.4、当前波长平顺性解算:
[0238] 3.4.1、波长输入:用户输入的波长λ1为数值30m;
[0239] 3.4.2、平顺性解算:
[0240] 3.4.2.1、高低平顺性解算:
[0241]
[0242]
[0243] 式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
[0244] dL1=L(j)-L(j-d)……………………………………………………… [18]
[0245] dL2=L(j+d)-L(j-d)…………………………………………………… [19]
[0246] dΘ1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
[0247] dΘ2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
[0248] dΘ3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
[0249] dΘ4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
[0250]
[0251] 3.4.2.2、轨向平顺性解算:
[0252]
[0253]
[0254]
[0255] 式中,y(j)为检测梁[16]的中心横向位移,β为幅值补偿系数,取0.03;
[0256] dψ1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
[0257] dψ2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
[0258] 至此,得到了轨道平顺性的检测结果,在30米波长下的轨道左高低平顺性、右高低平顺性和轨向平顺性检测精度均为±1毫米,小于用户要求的理论精度±2毫米。
[0259] 实施例2:
[0260] 检测装置和使用的检测方法,从步骤1到步骤3.3与实施例1完全相同。
[0261] 3.4、当前波长平顺性解算:
[0262] 3.4.1、波长输入:用户输入的波长λ1为数值300m;
[0263] 3.4.2、平顺性解算:
[0264] 3.4.2.1、高低平顺性解算:
[0265]
[0266]
[0267] 式中,hVL为轨道的左高低,hVR为轨道的右高低,α为幅值补偿系数,取0.25;
[0268] dL1=L(j)-L(j-d)……………………………………………………… [18]
[0269] dL2=L(j+d)-L(j-d)…………………………………………………… [19]
[0270] dΘ1=ΘL(j)-ΘL(j-d)…………………………………………………[20]
[0271] dΘ2=ΘL(j+d)-ΘL(j-d)………………………………………………[21]
[0272] dΘ3=ΘR(j)-ΘR(j-d)…………………………………………………[22]
[0273] dΘ4=ΘR(j+d)-ΘR(j-d)………………………………………………[23]
[0274]
[0275] 3.4.2.2、轨向平顺性解算:
[0276]
[0277]
[0278]
[0279] 式中,y(j)为检测梁[16]的中心横向位移,β为幅值补偿系数,取0.03;
[0280] dψ1=ΨH(j)-ΨH(j-d)…………………………………………………[28]
[0281] dψ2=ΨH(j+d)-ΨH(j-d)………………………………………………[29]
[0282] 至此,得到了轨道平顺性的检测结果,在300米波长下的轨道左高低平顺性、右高低平顺性和轨向平顺性检测精度均为±8毫米,小于用户要求的理论精度±10毫米。