CN109455197A - 轨道几何参数高精度检测方法和轨道几何参数高精度检测车 - Google Patents
轨道几何参数高精度检测方法和轨道几何参数高精度检测车 Download PDFInfo
- Publication number
- CN109455197A CN109455197A CN201811556962.XA CN201811556962A CN109455197A CN 109455197 A CN109455197 A CN 109455197A CN 201811556962 A CN201811556962 A CN 201811556962A CN 109455197 A CN109455197 A CN 109455197A
- Authority
- CN
- China
- Prior art keywords
- data
- coordinate data
- geometry parameter
- orbit geometry
- precision
- 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.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000001914 filtration Methods 0.000 claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 22
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 claims abstract description 14
- 229910052744 lithium Inorganic materials 0.000 claims abstract description 14
- 238000012937 correction Methods 0.000 claims abstract description 11
- 238000002203 pretreatment Methods 0.000 claims description 14
- 230000000903 blocking effect Effects 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 5
- 238000005299 abrasion Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 2
- 241000196324 Embryophyta Species 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 10
- 239000000758 substrate Substances 0.000 description 7
- 238000009434 installation Methods 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61K—AUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
- B61K9/00—Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
- B61K9/08—Measuring installations for surveying permanent way
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供一种轨道几何参数高精度检测方法和轨道几何参数高精度检测车,属于数据检测技术领域,轨道几何参数高精度检测车包括车体、轨道轮组、传感器组、倾角仪、数据采集卡、锂电池和工控机,车体包括底板、底部支架、检测梁、座椅、设备放置台和遮阳板,轨道轮组包括前轨道轮和后轨道轮,传感器包括两个激光传感器。轨道几何参数高精度检测方法为先对初始坐标数据依次进行滤波处理和误差校正处理后生成实际几何参数数据,滤波处理有效减小地面与杂草对初始坐标数据造成的干扰,从而保证最后生成的几何参数数据的准确性。
Description
技术领域
本发明涉及一种数据检测技术领域,具体地涉及一种基于有轨电车槽型轨的轨道几何参数高精度检测方法及轨道几何参数高精度检测车。
背景技术
近年来地铁、轻轨、现代有轨电车等城市轨道交通飞速发展,现代有轨电车以其在造价、能耗、使用寿命、建设周期、安全性等方面优于地铁和轻轨,能有效解决城市交通拥堵和交通污染等社会问题,适应现代城市低碳、快捷运输的需求,因而现代有轨电车在全国得到迅速推广和发展。目前全国已有十几个城市开通运营了有轨电车,并且超过100个城市提出了建设现代有轨电车线路的意向或规划。面对现代有轨电车迅猛的发展趋势,市场上却缺少对现代有轨电车槽型轨轨道几何参数高精度检测的设备,轨道作为现代有轨电车安全运行的基础,对其状态的检测是必不可少的。
现有一种轨道几何参数高精度检测方法采用激光位移传感器获取槽型轨的轮廓数据,并根据倾角仪获取的倾斜角数据、以及将获取的轮廓数据和预设的槽型轨模型数据进行比对而对槽型轨的轮廓数据进行校正,进而计算生成偏差更小的槽型轨几何参数数据。
现有的该种槽型轨几何参数检测方法存在的问题是,该方法在进行校正时,多采用轨顶点的数据作为运算基准,但实际上,轨顶点多被杂草或水泥地所覆盖,因而激光位移传感器获取的轮廓数据中轨顶点的坐标数据偏差较大,导致最后计算生成的轨道几何参数均存在较大偏差。
发明内容
本发明的第一目的在于提供一种有效减小检测偏差的轨道几何参数高精度检测方法。
本发明的第二目的在于提供一种有效减小检测偏差的轨道几何参数高精度检测车。
本发明第一目的提供的轨道几何参数高精度检测方法包括获取步骤,获取从激光传感器组输入的关于轨道点的初始坐标数据;预处理步骤,对初始坐标数据进行滤波处理后生成预处理坐标数据;校正步骤,根据预处理坐标数据对因车体倾斜角产生的误差进行校正后生成校正坐标数据。几何参数计算步骤,根据校正坐标数据生成实际几何参数数据。
由上述方案可见,处理器先对获取的初始坐标数据进行滤波处理,可对轨道上不同位置采用处理效果更好的滤波处理方式,从而获取偏离较小的预处理坐标数据,对经过滤波处理后的坐标数据进行校正以及参数计算,可有减小最终获得的轨道几何参数数据的计算偏差。
进一步的方案是,预处理步骤中,包括第一预处理步骤,对初始坐标数据整体进行限幅滤波处理生成第一预处理坐标数据。
由上可见,对初始坐标数据整体进行限幅滤波处理能有效去除划定坐标范围以外的干扰点,例如地面干扰点和杂草干扰点。
进一步的方案是,第一预处理步骤后包括第二预处理步骤,根据第一预设规则判断第一预处理坐标数据中的槽底坐标点区域,并对槽底坐标点区域进行高斯滤波处理。
由上可见,由于槽底坐标点区域的数据点呈围绕分布,且由于槽底区域磨损较小,采用高斯滤波能对槽底坐标点区域产生很好的滤波效果,随后计算以槽底区域的数据点作为参照能更有效较小计算偏差。
进一步的方案是,第一预处理步骤后包括第三预处理步骤,根据第一预设规则对第一预处理坐标数据中槽底坐标点区域以外的坐标点区域进行分段均值滤波处理。
由上可见,经测试,高斯滤波耗时较长,低通滤波、自适应滤波较均值滤波效果差,因此对槽底坐标点区域以外的坐标点区域采用均值滤波持恒滤波效果与处理效率。
进一步的方案是,第一预设规则为根据第一预处理坐标数据中各数据点的坐标数据对第一预处理坐标数据进行分段处理。
由上可见,例如以槽底点作为计算基准点,与槽底点坐标数据相差第一数值范围的其他数据点则位于轨腰左坐标点区域,与槽底点坐标数据相差第二数值范围的其他数据点则位于轨腰右坐标点区域。
进一步的方案是,几何参数计算步骤中包括,根据校正坐标数据中各数据点的坐标数据按第二预设规则计算生成关于轨道的轨宽数据。
由上可见,根据槽底点和轨腰两侧的数据点结合计算生成该槽型轨的轨宽值,偏差有效减小。
进一步的方案是,根据校正坐标数据按匹配算法计算生成关于轨道的磨损值数据。
由上可见,采用匹配算法将实测数据与标准数据拟合后再以左右两侧槽的圆心中点作为坐标原点重构坐标,再计算生成顶磨特征点和侧磨特征点,进而生成更准确的顶磨损值和侧磨损值。
本发明第二目的提供的轨道几何参数高精度检测车包括车体、轨道轮组、传感器组、倾角仪、数据采集卡和工控机,轨道轮组包括前轨道轮和后轨道轮,传感器包括两个激光传感器;轨道几何参数高精度检测车还包括锂电池,车体包括底板、底部支架、检测梁、座椅、设备放置台和遮阳板;座椅固定安装在底板上表面的中部,设备放置台固定连接在底板的上表面且位于座椅的朝向的前方,锂电池放置于底板的上表面上且位于座椅的后方,遮阳板和座椅的靠背共同从锂电池的上方对锂电池遮挡;底部支架包括框体,框体连接在底板的下表面上,且框体支撑于下表面的四周周边处;检测梁、前轨道轮和后轨道轮均安装于框体上,检测梁位于前轨道轮与后轨道轮之间;检测梁沿前轨道轮的转动轴线方向延伸;检测梁的延伸两端均安装有一个朝下设置的激光传感器,倾角仪安装在检测梁的延伸中部。
由上述方案可见,检测车上设置两个激光传感器、数据采集卡和工控机以实现上述的轨道几何参数高精度检测方法,对经过滤波处理后的坐标数据进行校正以及参数计算,可有减小最终获得的轨道几何参数数据的偏差;检测车上设置座椅以提供工作人员跟随检测车实时检测,发现异常可及时勘测,运算异常时及时中止纠错;遮阳板与座椅靠背共同作用对锂电池进行遮挡保护。
进一步的方案是,轨道几何参数高精度检测车还包括挡光件,检测梁的延伸两端均安装有一个挡光件,挡光件从水平方向对激光传感器遮挡。
更进一步的方案是,挡光件包括垂向设置的第一挡光壁、第二挡光壁和第三挡光壁,第一挡光壁、第二挡光壁和第三挡光壁相互连接成凵型,且第一挡光壁、第二挡光壁和第三挡光壁之间形成的开口朝向激光传感器。
由上可见,凵型的挡光件从检测车的横向外侧效遮挡光线,减少光线对激光传感器的影响,提高初始坐标数据的准确性。
附图说明
图1为本发明轨道几何参数高精度检测车实施例第一视角的结构图。
图2为本发明轨道几何参数高精度检测车实施例第二视角的结构图。
图3为本发明轨道几何参数高精度检测车实施中检测梁和传感器组的结构图。
图4为本发明轨道几何参数高精度检测车实施中底部支架的结构图。
图5为本发明轨道几何参数高精度检测方法实施例的流程图。
图6为本发明轨道几何参数高精度检测方法实施例中槽型轨的分段原理图。
以下结合附图及实施例对本发明作进一步说明。
具体实施方式
轨道几何参数高精度检测车实施例
参见图1和图2,图1为本发明轨道几何参数高精度检测车实施例第一视角的结构图,图2为本发明轨道几何参数高精度检测车实施例第二视角的结构图。轨道几何参数高精度检测车包括车体、轨道轮组、传感器组、倾角仪、数据采集卡、锂电池10和工控机71。车体包括底板1、底部支架8、检测梁3、座椅4、设备放置台7和遮阳板6。车体包括座椅41和座椅42,座椅41和座椅42在车体的横向上并排设置,座椅41具有后仰的靠背411,座椅42具有后仰的靠背421,座椅41和座椅42之间的置物盒43,置物盒43的可翻动上盖设置有内凹的圆形的放置位431。座椅41和座椅42固定安装在底板1上表面的中部,设备放置台7固定连接在底板1的上表面且位于座椅41的朝向的前方,工控机71放置在设备放置台7上;锂电池10放置于底板1的上表面上且位于座椅41和座椅42的后方。
遮阳板6包括朝向车体后方且朝上设置的第一板块61和朝向设置的第二板块62,第一板块61、第二板块62、靠背411和靠背421共同从锂电池10的上方对锂电池10遮挡。第一板块61和第二板块62可选择太阳能电池板支撑。
底部支架8连接在底板1的下表面上,且底部支架8支撑于底板1下表面的四周周边处;轨道轮组包括前轨道轮21和后轨道轮22,检测梁3、前轨道轮21和后轨道轮22均安装于底部支架上,检测梁3位于前轨道轮21与后轨道轮22之间;检测梁3沿前轨道轮21的转动轴线方向延伸。
结合图3,图3为本发明轨道几何参数高精度检测车实施中检测梁和传感器组的结构图。检测梁3包括一体成型的基板31、两块侧板32和多块肋板33;基板31和侧板32均沿前轨道轮21的转动轴线方向延伸,两块侧板32对称地连接在基板31的两侧,侧板32垂直于基板31,基板31和两块侧板32之间形成U型槽30,U型槽30沿前轨道轮21的转动轴线方向延伸,多块肋板33间隔地设置在U型槽30中,肋板33垂直于基板31且垂直于侧板32,肋板33连接在两块侧板32与基板31之间。基板31和侧板32在检测梁3的延伸方向的中部均具有较大面积,从而使检测梁3的延伸中部具有更大的截面积。
传感器组包括两个激光传感器51和两个双轴加速度传感器52,检测梁3的延伸两端均安装有一个朝下设置的激光传感器51和一个双轴加速度传感器52,倾角仪53安装在检测梁3的延伸中部。检测梁3的延伸两端均安装有一个挡光件9,挡光件9从检测车的横向外侧对激光传感器51遮挡。挡光件9包括垂向设置的第一挡光壁91、第二挡光壁92和第三挡光壁93,第一挡光壁91、第二挡光壁92和第三挡光壁93相互连接成凵型,且第一挡光壁91、第二挡光壁92和第三挡光壁93之间形成的开口朝向激光传感器51。
结合图4,图4为本发明轨道几何参数高精度检测车实施中底部支架的结构图。底部支架8一体成型,底部支架8包括一个框体801和设置在框体801中部的纵向梁802和横向梁803,框体801横截面呈工字型,框体801包括依次连接的上端面板81、垂向板82和下端面板83,垂向板82垂向设置,上端面板81和下端面板83均垂直于垂向板82设置。
上端面板81为表面平整的面板,用于与底板1固定连接;纵向梁802截面呈T型,纵向梁802的上面板与上端面板81位于同一平面。垂向壁82在检测车前进方向的前端和后端的部分上设置有一个弓形结构,弓形结构包括直线部821、弧形部822和形成于直线部821与弧形部822之间的开口823,开口823的内边上连接有一圈完整的弓形壁84,弧形部822朝地面方向从框体801的下端拱起。
框体801在垂直于检测车前进方向的横向的两侧均设置有一个安装框架85,安装框架85的壁体86垂直于垂向壁82,且安装框架85的壁体86从下端面板83朝地面方向下凸;横向梁803连接于两个安装框架85之间,横向梁803截面呈T型;检测梁3的延伸两端通过螺栓固定安装在两侧的两个安装框架85的壁体86上。
底部支架的工字型截面结构以及梁体的T型截面结构在减轻车体整体重量的同时保证检测车的力学强度;安装框架85的壁体86的下凸设置能提供激光传感器51更好的检测距离,同时保证连接处的力学强度。
轨道几何参数高精度检测方法实施例
参见图5,图5为本发明轨道几何参数高精度检测方法实施例的流程图。本实施例使用上轨道几何参数高精度检测车实施例中的检测***对轨道的几何参数数据进行检测和计算。工控机71内具有处理器,激光传感器51、双轴加速度传感器52和倾角仪53通过数据采集卡和工控机71连接。
首先由处理器执行步骤S1,获取来自激光传感器的初始坐标数据。轨道轮组的转轴上的光电编码器发出脉冲信号触发激光传感器51、双轴加速度传感器52和倾角仪53工作并获取对应的数据信息,并传输至处理器中,激光传感器51则获取关于左、右两个槽型轨的初始坐标数据,并将左、右两个槽型轨的初始坐标数据的数据点整合到同一坐标系中:
采集所得数据用坐标表示为左轨:(x10,y10),(x11,y11),…,(x1n,y1n),
右轨:(x20,y20),(x21,y21),…,(x2n,y2n);
矩阵表示为:左轨:X1,Y1,右轨:X2,Y2;
随后对两个2D激光传感器采集到的数据进行坐标转换整合到同一坐标系,将右轨坐标系转换至左轨坐标系,转换模型如下:
bx,by为两个激光传感器之间分别在X轴方向和Y轴方向的距离,m为尺度参数,α为两坐标系的角度差,已知z坐标不变,转换后右轨的检测点数据坐标为:(x′20,y′20),(x′21,y′21),…,(x′2n,y′2n),随后获得
原始组数据:左轨:X1,Y1右轨:X′2,Y2。
随后处理器执行步骤S2,对获取的初始坐标数据进行滤波处理。结合图6,图6为本发明轨道几何参数高精度检测方法实施例中槽型轨的分段原理图。首先对初始坐标数据进行限幅滤波处理并生成第一预处理坐标数据,对初始坐标数据整体进行限幅滤波处理能有效去除划定坐标范围以外的干扰点,例如地面干扰点和杂草干扰点。
然后处理器根据第一预处理坐标数据中各个数据点的坐标对第一预处理坐标数据进行分段处理,对左槽型轨的坐标数据点集和右槽型轨的坐标数据点集均划分为宽轨顶坐标点区域、轨腰坐标点区域、槽底坐标点区域和窄轨侧坐标点区域,图6所示为左轨11的轨面分段示意图,其包括nL个坐标数据点。首先获取左轨11槽底点110的坐标数据记录为CL,并将CL作为计算基准按以下公式组计算以对左轨数据点进行分段:
第一段(宽轨顶坐标点区域11a):
startL1=1
第二段(轨腰坐标点区域11b):
startL2=endL1+1
第三段(槽底坐标点区域11c):
startL3=endL2+1
第四段(窄轨侧坐标点区域11d):
startL4=endL3+1
随后处理器对第一预处理坐标数据进行第二次滤波处理,对槽底坐标点区域进行高斯滤波处理,对宽轨顶坐标点区域、轨腰坐标点区域和窄轨侧坐标点区域进行均值滤波处理。由于槽底坐标点区域的数据点呈围绕分布,且由于槽底区域磨损较小,采用高斯滤波能对槽底坐标点区域产生很好的滤波效果,随后计算以槽底区域的数据点作为参照能更有效较小计算偏差;经测试,高斯滤波耗时较长,低通滤波、自适应滤波较均值滤波效果差,因此对槽底坐标点区域以外的坐标点区域采用均值滤波持恒滤波效果与处理效率。滤波处理完成后,处理器生成预处理坐标数据。
随后处理器执行步骤S3,对预处理坐标数据进行校正并生成校正坐标数据。由于左轨道数据点的校正和左轨道数据点的校正规则相同,故仅以坐轨道点的校正进行具体说明。
先对预处理坐标数据中关于车体左右倾斜的误差进行校正:
A、首先计算出轨面顶部特征点。对于左侧轨道的轨道数据点,根据轨道数据点依次计算相邻轨面顶部特征点之差,可得:
Δy′1i=y′1i-y′1(i-1)i=1,2,…,n
B、选定相邻十个差值进行求和:
C、找出s1i值最小的一组数据对应的点作为校正所需的左轨轨面顶部特征点,即Δy′1s,Δy′1(s+1),…,Δy′1(s+9)。
D、同理可以得到右轨轨面顶部特征点,即
Δy′2s,Δy′2(s+1),…,Δy′2(s+9)。
E、根据提取到的两组轨面顶部特征点数据(以下简称轨顶面点),采用最小二乘法将所取的轨顶面点拟合成一条直线,可以得出其斜率k,故所得直线方程为:
y=kx+b
F、根据此斜率可以得出车体倾斜角度β:
β=tan-1k
G、然后根据此斜率可以校正轨道数据点误差。利用车体倾斜角度β对轨道数据点的横纵坐标进行误差校正,根据公式:
X*=Xcosβ+Ysinβ
Y*=Ycosβ-Xsinβ
经校正后得到的轨道数据点为:左轨:右轨:即第一校正坐标数据。
后对第一校正坐标数据中关于车体前后倾斜的误差进行校正:
A、提取轨宽特征点:
首先,找到左轨槽底最低点y1min;将直线y=y1min向上平移14mm后得到直线y=y1min+14,然后将轨腰两侧数据点进行直线拟合,得到直线y1a=k1ax+b1a,y1b=k1bx+b1b;最后将直线y=y1min+14与直线y1a=k1ax+b1a和y1b=k1bx+b1b分别联立,得到交点,即左轨的两个轨宽点A1(x1a,y1a),B1(x1b,y1b),可以得到右轨的轨宽点:A2(x2a,y2a),B2(x2b,y2b)。
B、根据计算所得轨宽点的横坐标,得到参考轨宽:d1=x1b-x1a、d2=x2b-x2a。
C、根据槽型轨的标准轨宽dz,计算得到偏角α:
D、根据偏角α校正车体前后倾斜即车体航向与轨向不一致产生的误差和轨道数据点的误差并生成校正坐标数据。
随后处理器执行步骤S4,根据校正坐标数据生成轨道几何参数数据。轨道几何参数数据包括轨道磨耗值。基于匹配算法,通过将无磨耗标准坐标数据和校正坐标数据(有磨耗实测坐标数据)进行拟合而实现对两个模型的匹配,再通过计算两者特征点差值而计算出槽型轨的垂直磨耗和侧面磨耗。其中,由于槽型轨凹槽部分磨耗较小,计算时选择凹槽数据点作为拟合基准。
假定,校正坐标数据为Xa={(xa1,ya1),(xa2,ya2),(xa2,ya2),(xa3,ya3)…(xan,yan)},首先获取轨道的槽底点x值xcd1:在Xa中找出y值最小的点(xmin,ymin),然后将直线y=ymin向上平移10mm生成直线y=ymin+10。在x值处于区间(xa1,xmin)和(xmin,xan)各找出最靠近直线y=ymin+10的两个点(xmin1,ymin1)和(xmin2,ymin2)并求出xmin1和xmin2的平均值作为槽底点x值:xcd1=(xmin1+xmin2)/2。
然后对校正坐标数据进行旋转校正处理。在Xa中找出x值处于区间(xcd1-17,xcd1-15)的点形成点集Xb,对Xb进行最小二乘法得出回归直线y=k11x+b11计算Xb的点与直线y=k11x+b11的距离并提取距离小于1mm形成点集Xc,对Xc进行最小二乘法得出回归直线y=k1x+b1。同理,在Xa中找出x值处于区间(xcd1+15,xcd1+17)的点形成点集Xd,经过两次直线回归得出回归直线y=k2x+b2根据两条回归直线方程求出倾斜角δδ=(tan-1k1+tan-1k2)/2,最后对Xa的点以坐标原点为旋转中心逆时针旋转δ,得出校正后的数据点集:
Xe={(xe1,ye1),(xe2,ye2),(xe2,ye2),(xe3,ye3)…(xen,yen)}
xei=xai*cosδ+yai*sinδ,i=1,2,3…n
yei=yai*cosδ-xai*sinδ,i=1,2,3…n
最后固定原点并重构坐标系。
利用求出Xe的槽底点x值xcd2在Xe中找出x值处于区间(xcd2-13,xcd2-10)的点形成点集
Xf={(xf1,yf1),(xf2,yf2),(xf2,yf2),(xf3,yf3)…(xfn,yfn)},对Xf进行圆拟合求出圆心坐标值(xc1,yc1)。本项目采取圆拟合算法如下,已知圆的标准方程为:
(x-xc1)2+(y-yc1)2=r2
其中,根据槽型轨的设计标准可知半径值r=13.8。
然后根据最小二乘法原理来求出圆心坐标值(xc1,yc1)的最优解,由于对于圆方程来说计算y值的残差的平方和计算过于复杂所以用数据点与圆心的距离与半径的差代替残差,即求的最小值。对xc1和yc1求偏导计算极小值,可得方程组:
上述方程组为二元非线性方程组,采取牛顿迭代法求解:
F(x)的Jacobian矩阵为:
根据牛顿迭代公式x(k+1)=x(k)-FJ(x(k))-1*F(x(k))可求出圆心坐标值(xc1,yc1)的近似值作为圆心坐标值(xc1,yc1)。
随后计算Xf中的点与圆心的距离并提取距离小于1mm形成点集Xg,用拟合圆算法再一次求Xg的圆心更新圆心坐标值(xc1,yc1)。在Xe中找出x值处于区间(xcd2+10,xcd2+13)的点形成点集Xh,对Xh用上述方法两次圆拟合求出圆心坐标值(xc2,yc2)。以两个圆心的中点作为坐标的原点,修正Xe得出点集
Xj={(xf1,yf1),(xf2,yf2),(xf2,yf2),(xf3,yf3)…(xfn,yfn)}即
上述所求的Xj为完成模型匹配的重构坐标轨道数据。
随后执行轨道磨耗值计算,假设重构坐标轨道数据为:
Xa={(xa1,ya1),(xa2,ya2),(xa2,ya2),(xa3,ya3)…(xan,yan)}
A.求出左轨标准顶磨特征点:
获取Xa槽底点x值xcd,在Xa中找出x值处于区间(xcd-37.785,xcd-35.785)的点形成点集Xb,对Xb进行最小二乘法得出回归直线y=k1x+b1,得左轨标准顶磨特征点ldm:
ldm=k1*(xcd-36.785)+b1
B.求出左轨标准侧磨特征点:
获取Xa中y值最大的点(xmax,ymax),后在Xa中找出x值处于区间(0,xcd)且y值处于区间(ymax-17,ymax-15)的点形成点集Xc,对Xc进行最小二乘法得出回归直线y=k2x+b2,则左轨标准顶磨特征点lcm:
lcm=(ymax-17-b2)/k2
C.同上A和B的算法求出左轨的实测顶磨特征点ldm′和实测侧磨特征点lcm′:
Ldm=l′dm-ldm
Lcm=l′cm-lcm
左轨道总磨耗Lmh则为Lmh=Ldm+Lcm/2。
本发明提供的轨道几何参数高精度检测方法中,处理器先对获取的初始坐标数据进行滤波处理,对轨道上不同位置采用处理效果更好的滤波处理方式,从而获取偏差较小的预处理坐标数据,对经过滤波处理后的坐标数据进行校正以及参数计算,可以减小最终获得的轨道几何参数数据的偏差。
最后需要强调的是,以上所述仅为本发明的优选实施例,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种变化和更改,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.轨道几何参数高精度检测方法,其特征在于,包括:
获取步骤,获取从激光传感器组输入的关于轨道点的初始坐标数据;
预处理步骤,对所述初始坐标数据进行滤波处理后生成预处理坐标数据;
校正步骤,根据所述预处理坐标数据对因车体倾斜角产生的误差进行校正后生成校正坐标数据;
几何参数计算步骤,根据所述校正坐标数据生成实际几何参数数据。
2.根据权利要求1所述的轨道几何参数高精度检测方法,其特征在于:
所述预处理步骤中,包括:
第一预处理步骤,对所述初始坐标数据整体进行限幅滤波处理生成第一预处理坐标数据。
3.根据权利要求2所述的轨道几何参数高精度检测方法,其特征在于:
执行第一预处理步骤后,还执行:
第二预处理步骤,根据第一预设规则判断所述第一预处理坐标数据中的槽底坐标点区域,并对所述槽底坐标点区域进行高斯滤波处理。
4.根据权利要求3所述的轨道几何参数高精度检测方法,其特征在于:
执行第一预处理步骤后,还执行:
第三预处理步骤,根据第一预设规则对所述第一预处理坐标数据中槽底坐标点区域以外的坐标点区域进行分段均值滤波处理。
5.根据权利要求3或4所述的轨道几何参数高精度检测方法,其特征在于:
所述第一预设规则为:
根据所述第一预处理坐标数据中各数据点的坐标数据对所述第一预处理坐标数据进行分段处理。
6.根据权利要求1至4任一项所述的轨道几何参数高精度检测方法,其特征在于:
所述几何参数计算步骤具体包括:
根据所述校正坐标数据中各数据点的坐标数据按第二预设规则计算生成关于轨道的轨宽数据。
7.根据权利要求1至4任一项所述的轨道几何参数高精度检测方法,其特征在于:
所述几何参数计算步骤中,包括
根据所述校正坐标数据按匹配算法计算后生成关于轨道的磨损值数据。
8.用于实现上述权利要求1至7任一项所述的轨道几何参数高精度检测方法的轨道几何参数高精度检测车,包括车体、轨道轮组、传感器组、倾角仪、数据采集卡和工控机,所述轨道轮组包括前轨道轮和后轨道轮,所述传感器组包括两个激光传感器;
其特征在于:
所述轨道几何参数高精度检测车还包括锂电池;
所述车体包括底板、底部支架、检测梁、座椅、设备放置台和遮阳板;
所述座椅固定安装在底板上表面的中部,所述设备放置台固定连接在所述底板的上表面且位于所述座椅的朝向的前方,所述锂电池放置于所述底板的上表面上且位于所述座椅的后方,所述遮阳板和所述座椅的靠背共同从所述锂电池的上方对所述锂电池遮挡;
所述底部支架连接在所述底板的下表面上,且所述底部支架支撑于所述下表面的四周周边处;
所述检测梁、所述前轨道轮和所述后轨道轮均安装于所述底部支架上,所述检测梁位于所述前轨道轮与所述后轨道轮之间;所述检测梁沿所述前轨道轮的转动轴线方向延伸;
所述检测梁的延伸两端均安装有一个朝下设置的所述激光传感器,所述倾角仪安装在所述检测梁的延伸中部。
9.根据权利要求8所述的轨道几何参数高精度检测车,其特征在于:
所述轨道几何参数高精度检测车还包括挡光件,所述检测梁的延伸两端均安装有一个所述挡光件,所述挡光件从所述检测车行驶方向的横向外侧对所述激光传感器遮挡。
10.根据权利要求9所述的轨道几何参数高精度检测车,其特征在于:
所述挡光件包括垂向设置的第一挡光壁、第二挡光壁和第三挡光壁,所述第一挡光壁、所述第二挡光壁和所述第三挡光壁相互连接成凵型,且所述第一挡光壁、所述第二挡光壁和所述第三挡光壁之间形成的开口朝向所述激光传感器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811556962.XA CN109455197A (zh) | 2018-12-19 | 2018-12-19 | 轨道几何参数高精度检测方法和轨道几何参数高精度检测车 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811556962.XA CN109455197A (zh) | 2018-12-19 | 2018-12-19 | 轨道几何参数高精度检测方法和轨道几何参数高精度检测车 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109455197A true CN109455197A (zh) | 2019-03-12 |
Family
ID=65613880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811556962.XA Pending CN109455197A (zh) | 2018-12-19 | 2018-12-19 | 轨道几何参数高精度检测方法和轨道几何参数高精度检测车 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109455197A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110046664A (zh) * | 2019-04-16 | 2019-07-23 | 中国石油大学(华东) | 一种油田安全突出问题检测方法 |
CN111307078A (zh) * | 2020-02-21 | 2020-06-19 | 暨南大学 | 一种基于四点弦测法的轨道不平顺检测方法及*** |
CN112964195A (zh) * | 2021-02-04 | 2021-06-15 | 暨南大学 | 基于激光三角法的供电轨几何参数综合检测方法及*** |
CN113104056A (zh) * | 2021-05-20 | 2021-07-13 | 成都盛锴科技有限公司 | 一种分体式工务巡检机器人平台 |
CN113779672A (zh) * | 2021-08-31 | 2021-12-10 | 北京铁科英迈技术有限公司 | 一种钢轨廓形磨耗计算方法、装置、设备及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090073428A1 (en) * | 2007-08-22 | 2009-03-19 | Steven Magnus | Rail measurement system |
JP2014017611A (ja) * | 2012-07-06 | 2014-01-30 | Toshiba Corp | 軌道点検装置 |
WO2014147563A1 (es) * | 2013-03-18 | 2014-09-25 | Universidad Eafit | Sistema y método para la inspección de los parámetros geométricos de ruedas de vehículos ferroviarios |
CN107380204A (zh) * | 2017-08-11 | 2017-11-24 | 暨南大学 | 轨道几何参数检测车及轨道几何参数检测方法 |
CN107554553A (zh) * | 2017-08-31 | 2018-01-09 | 常州路航轨道交通科技有限公司 | 基于二维激光位移传感器的轨道几何不平顺检测方法 |
CN207274714U (zh) * | 2017-08-11 | 2018-04-27 | 暨南大学 | 轨道数据检测车 |
CN209426782U (zh) * | 2018-12-19 | 2019-09-24 | 暨南大学 | 轨道几何参数高精度检测车 |
-
2018
- 2018-12-19 CN CN201811556962.XA patent/CN109455197A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090073428A1 (en) * | 2007-08-22 | 2009-03-19 | Steven Magnus | Rail measurement system |
JP2014017611A (ja) * | 2012-07-06 | 2014-01-30 | Toshiba Corp | 軌道点検装置 |
WO2014147563A1 (es) * | 2013-03-18 | 2014-09-25 | Universidad Eafit | Sistema y método para la inspección de los parámetros geométricos de ruedas de vehículos ferroviarios |
CN107380204A (zh) * | 2017-08-11 | 2017-11-24 | 暨南大学 | 轨道几何参数检测车及轨道几何参数检测方法 |
CN207274714U (zh) * | 2017-08-11 | 2018-04-27 | 暨南大学 | 轨道数据检测车 |
CN107554553A (zh) * | 2017-08-31 | 2018-01-09 | 常州路航轨道交通科技有限公司 | 基于二维激光位移传感器的轨道几何不平顺检测方法 |
CN209426782U (zh) * | 2018-12-19 | 2019-09-24 | 暨南大学 | 轨道几何参数高精度检测车 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110046664A (zh) * | 2019-04-16 | 2019-07-23 | 中国石油大学(华东) | 一种油田安全突出问题检测方法 |
CN111307078A (zh) * | 2020-02-21 | 2020-06-19 | 暨南大学 | 一种基于四点弦测法的轨道不平顺检测方法及*** |
CN112964195A (zh) * | 2021-02-04 | 2021-06-15 | 暨南大学 | 基于激光三角法的供电轨几何参数综合检测方法及*** |
CN113104056A (zh) * | 2021-05-20 | 2021-07-13 | 成都盛锴科技有限公司 | 一种分体式工务巡检机器人平台 |
CN113104056B (zh) * | 2021-05-20 | 2024-03-22 | 成都盛锴科技有限公司 | 一种分体式工务巡检机器人平台 |
CN113779672A (zh) * | 2021-08-31 | 2021-12-10 | 北京铁科英迈技术有限公司 | 一种钢轨廓形磨耗计算方法、装置、设备及存储介质 |
CN113779672B (zh) * | 2021-08-31 | 2024-04-12 | 北京铁科英迈技术有限公司 | 一种钢轨廓形磨耗计算方法、装置、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109455197A (zh) | 轨道几何参数高精度检测方法和轨道几何参数高精度检测车 | |
CN107389028B (zh) | 一种基于坐标投影的三维坐标转换方法及装置 | |
CN209426782U (zh) | 轨道几何参数高精度检测车 | |
CN107607044B (zh) | 一种基于激光位移传感器的车轮踏面磨耗检测方法 | |
CN100370222C (zh) | 确定轮缘空间位置的测量方法和装置,及车轮定位测量*** | |
CN103557788B (zh) | 一种高铁接触网接几何参数检测非接触式补偿及卡尔曼滤波修正方法 | |
CN100449259C (zh) | 车辆轮对直径在线检测方法及装置 | |
CN108917638A (zh) | 基于基准传递的地铁隧道三维变形监测的车载测量装置 | |
CN111532295A (zh) | 轨道交通移动智能运维检测*** | |
CN107380204A (zh) | 轨道几何参数检测车及轨道几何参数检测方法 | |
CN106225710B (zh) | 基于误差修正的列车轮踏面三维轮廓自动化测量方法和*** | |
CN101055230B (zh) | 用于无接触式3d车轮校准的设备、***和方法 | |
CN108362308A (zh) | 一种利用隧道环缝的移动激光测量***里程校正方法 | |
CN101144714A (zh) | 一种钢轨磨耗综合参数车载动态测量装置及方法 | |
CN207274714U (zh) | 轨道数据检测车 | |
CN105946898A (zh) | 一种基于激光测距的城轨列车车轮直径检测方法及*** | |
CN108457143A (zh) | 一种轨道线路坐标测量*** | |
CN109000568A (zh) | 一种钢轨超高的非接触式测量方法 | |
CN107685748A (zh) | 基于激光位移传感器的列车车轮尺寸参数在线检测方法 | |
CN108534765B (zh) | 一种轨道线路坐标测量方法 | |
CN106515784B (zh) | 一种轨道车辆车体尺寸测量方法 | |
CN103954458A (zh) | 一种非接触式四轮定位仪及其检测方法 | |
CN103344297B (zh) | 地下水封洞库的库容测量方法 | |
CN109334709A (zh) | 一种用于测量车轮轮缘参数的装置及方法 | |
CN114046789A (zh) | 轨检小车导航imu协同测量的铁轨检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |