高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于位置敏感探测器的轨距在线监测系统研究

刘子英 张靖 邓芳明

引用本文:
Citation:

基于位置敏感探测器的轨距在线监测系统研究

    作者简介: 刘子英(1964-),女,教授,研究方向为检测与监测控制。E-mail:151696128@qq.com.
  • 基金项目:

    国家自然科学基金资助项目 51767006

    江西省应用研究培育计划资助项目 20181BBE58015

    江西省重点研发计划资助项目 20181BBE50019

    江西省教育厅科学技术资助项目 180310

  • 中图分类号: TN247

Research of on-line monitoring of rail distance system based on position sensitive detector

  • CLC number: TN247

  • 摘要: 为了实时了解铁路轨道的轨距变化情况,研制了一种实时在线监测轨距变化的测量系统。采用位置敏感探测器和激光发射器作为测量元件实时采集轨距的变化信息,建立数学模型计算轨距的偏移量,通过组合算法来获取两根钢轨移动状态后的偏移量,并通过二次插值算法对位置敏感探测器进行非线性补偿以提高位置敏感探测器的精度和准确度;采用控制器局域网络总线和通用分组无线服务技术将一定区域内的轨距变化信息发送到监控中心进行显示;通过和现有二级铁路轨距尺分别对同一模拟轨道进行了测量实验。结果表明,该系统的测量精度为0.32mm,达到了目前轨道检测行业所规定的1mm精度要求。
  • Figure 1.  Structure diagram of gauge monitoring system

    Figure 2.  Gauge measurement schematic

    Figure 3.  Calculated value of gauge offset Δx

    Figure 4.  Calculation model of Δx1

    Figure 5.  2-D PSD detection schematic

    Figure 6.  Compensation grid for quadratic interpolation algorithm

    Figure 7.  Test scenario of rail gauge monitoring system

    Figure 8.  Curved curve of rail gauge measurement system and the measured data

    Table 1.  Comparison between the ruler of railway gauge and the system

    test location measurement method average value/mm difference/mm standard deviation/mm
    1 ruler of railway gauge 1434.14 -0.35 0.09
    system 1434.49 0.30
    2 ruler of railway gauge 1435.18 -0.29 0.09
    system 1435.47 0.29
    3 ruler of railway gauge 1433.17 -0.37 0.10
    system 1433.54 0.29
    4 ruler of railway gauge 1436.17 0.61 0.12
    system 1435.46 0.27
    5 ruler of railway gauge 1433.17 -0.31 0.11
    system 1433.48 0.27
    6 ruler of railway gauge 1435.19 -0.34 0.10
    system 1435.53 0.32
    7 ruler of railway gauge 1434.19 -0.28 0.10
    system 1434.47 0.31
    8 ruler of railway gauge 1434.18 0.22 0.09
    system 1434.40 0.28
    下载: 导出CSV
  • [1]

    WANG Y M, LU M. Interpretation of JJG 219-2015'track gage for standard gauge rail way'[J]. Railway Quality Control, 2016, 44(8):13-17(in Chinese).
    [2]

    CHEN D Sh, TIAN X Y. Development of China's high-speed railway track detectiontechnology[J]. Railway Engineering, 2008(12):82-86(in Chinese).
    [3]

    XIONG Sh Y, CHEN Ch J, WANG F, et al. Research on a new method of dynamic gauge detection[J]. Journal of Railway Science and Engineering, 2018, 15(7):1825-1831(in Chinese).
    [4]

    ZHENG Sh B, CHAI X D, HANG G G, et al. Research on gauge detection method based on machine vision[J]. Urban Mass Transit, 2010, 13(9):73-76(in Chinese).
    [5]

    ZHENG Sh B, CHAI X D, AN X C, et al. Railway track gauge inspection method based on computer vision[C]//2012 IEEE International Conference on Mechatronics and Automation. New York, USA: IEEE, 2012: 1292-1296.
    [6]

    MIN Y Zh, WANG H X, KANG F, et al. Research on railway gauge detection system based on image sensor[J]. Laser Technology, 2015, 39(3):344-348(in Chinese).
    [7]

    WANG Y. Synchronous trigger design of laser camera based on GJ-6 track detection system[J]. China Railway Science, 2018, 39(4):139-144(in Chinese).
    [8]

    SHI H M, ZHANG J K. Research on gauge detection system based on laser triangulation principle[J]. Chinese Journal of Scientific Instrument, 2013, 34(9):1934-1940(in Chinese). 
    [9]

    TANG W, LONG Y, XIE Y, et al. A channel rail gauge detection system of modern tram[C]//2018 IEEE 15th International Confe-rence on Networking, Sensing and Control (ICNSC). New York, USA: IEEE, 2018: 1-5.
    [10]

    KARTTUNEN K, KABO E, EKBERG A. Estimation of gauge corner and flange root degradation from rail, wheel and track geometries[J]. Wear, 2016, 366/367:294-302. doi: 10.1016/j.wear.2016.03.030
    [11]

    BAE J, KIM J. A rail detection algorithm based on pair particles filtering[C]//2017 Asia-Pacific Signal and Information Processing Association Annual Summit and Conference (APSIPA ASC). Hong Kong, China: APSIPA, 2017: 1440-1443.
    [12]

    ASHAR A Z, GANESH N, NARAYAN K S. Hybrid perovskite-based position-sensitive detectors[J]. Advanced Electronic Materials, 2018, 2(4):1-4.
    [13]

    KONG Q. Research and implementation of phase laser ranging technology[D]. Chengdu: Sichuan Normal University, 2018: 10-11(in Chinese).
    [14]

    SONG H, ZHANG Y F, WU Ch P, et al. Underwater phase laser ranging calibration method[J].Infrared and Laser Technology, 2019, 48(4):0406008(in Chinese). doi: 10.3788/IRLA201948.0406008
    [15]

    HAI Y, LI Y, MA X H, et al. Optical system design for small size laser ranging[C]//2015 International Conference on Optoelectronics and Microelectronics (ICOM). New York, USA: IEEE, 2015: 115-118.
    [16]

    XU Y, ZHU X, CHEN S H, et al. Optical design of dynamic focusing system for laser processing[C]//2017 Conference on Lasers and Electro-Optics Pacific Rim (CLEO-PR). New York, USA: IEEE, 2017: 1-2.
    [17]

    ZHANG M, LIANG Y B. Experimental study on two-dimensional quadrilateral position sensitive detector[J].Infrared and Laser Engineering, 2013, 42(2):459-464(in Chinese).
    [18]

    ZHANG H, CHEN Y S, GEN T W, et al. Study on main factors affecting position detection accuracy of four-quadrant detector[J]. Chinese Journal of Lasers, 2015, 42(12): 1217002 (in Chinese). doi: 10.3788/CJL201542.1217002
    [19]

    JIANG X D, YU J Y, ZHU L K. Research of combined navigation technology based on position sensitive detectors[J]. Laser Techno-logy, 2019, 43(3): 335-340 (in Chinese). 
    [20]

    HUANG P, CHEN B, LI Y. Displacement detection system based on phase sensitive detection and position detector[J]. Laser Technology, 2017, 41(2):235-239(in Chinese). 
    [21]

    XU C, TIAN S Y, YANG J Ch. Two-dimensional PSD nonlinear correction algorithm based on fuzzy neural network[J]. Electronic Measurement Technology, 2017, 40(2):50-54(in Chinese).
    [22]

    LIN Q S, YANG X J, WANG J X, et al. Two-dimensional PSD nonlinear correction of improved BP neural network[J]. Laser Technology, 2012, 36(1):124-126(in Chinese).
    [23]

    LI Q, LI T Z, XIE A Zh. Two-dimensional PSD nonlinear correction based on newton interpolation algorithm[J]. Journal of Shandong University of Technology(Science and Technology Edition), 2015, 29(3):69-72(in Chinese).
  • [1] 蒋晓东于纪言朱立坤 . 基于位置敏感探测器的组合导航技术研究. 激光技术, 2019, 43(3): 335-340. doi: 10.7510/jgjs.issn.1001-3806.2019.03.009
    [2] 王红敏罗自赢李蕾彭俊玲 . 平行度误差激光准直法测量技术研究. 激光技术, 2022, 46(5): 674-679. doi: 10.7510/jgjs.issn.1001-3806.2022.05.016
    [3] 范志刚张凤生左保军何瑾 . 位置敏感探测器受杂光干扰特性的研究. 激光技术, 2004, 28(4): 442-444.
    [4] 吴斌王晓蕾 . 一种可变长度光电基准尺. 激光技术, 2016, 40(5): 734-737. doi: 10.7510/jgjs.issn.1001-3806.2016.05.024
    [5] 曹玲杨玲珍郭雄英周希坚 . 非线性光纤环形镜掺铒光纤激光器的实验研究. 激光技术, 2007, 31(6): 613-615,619.
    [6] 闵永智王红霞康飞党建武 . 基于图像式传感器的铁路轨距检测系统研究. 激光技术, 2015, 39(3): 344-348. doi: 10.7510/jgjs.issn.1001-3806.2015.03.014
    [7] 杨振刚陈海清王磊任温馨 . 利用内腔变形反射镜补偿激光器畸变的研究. 激光技术, 2007, 31(5): 449-451.
    [8] 陈蔚曹益平龚赤坤曾钦勇甘春泉 . 对光纤阵列激光告警器定向算法的改进. 激光技术, 2008, 32(3): 252-254,258.
    [9] 齐艺超陈伟穆春元祝宁华 . 基于粒子群自整定PID算法的激光器温度控制系统. 激光技术, 2019, 43(5): 650-654. doi: 10.7510/jgjs.issn.1001-3806.2019.05.012
    [10] 张小光李乙钢冯鸣吕可诚 . 一种级联喇曼光纤激光器数值模拟的新算法. 激光技术, 2010, 34(2): 250-252. doi: 10.3969/j.issn.1001-3806.2010.02.029
    [11] 梁巍巍黄振宇张文攀殷瑞光刘艳芳 . 激光导引头四象限探测器偏差信号特性研究. 激光技术, 2014, 38(4): 569-573. doi: 10.7510/jgjs.issn.1001-3806.2014.04.027
    [12] 董吉辉胡企铨 . 空间碎片的激光雷达探测. 激光技术, 2007, 31(2): 185-187,191.
    [13] 张宇赵远吴晓敏陈钟贤孙秀冬 . 多光谱探测与激光多光谱探测技术的进展. 激光技术, 2007, 31(2): 188-191.
    [14] 何建平周智吴源华欧进萍 . 光纤布里渊与布喇格光栅共线技术的温度互补偿. 激光技术, 2010, 34(1): 13-16. doi: 10.3969/j.issn.1001-3806.2010.01.004
    [15] 刘兵陶炜柯尊贵冯力天袁菲李晓峰 . 相干激光雷达平衡式相干探测技术研究. 激光技术, 2015, 39(1): 46-49. doi: 10.7510/jgjs.issn.1001-3806.2015.01.009
    [16] 薛芳谭勇邬志强王玉诏 . 基于激光遥感的酒驾遥感探测技术研究. 激光技术, 2019, 43(1): 93-97. doi: 10.7510/jgjs.issn.1001-3806.2019.01.019
    [17] 马泳冀航梁琨林宏 . 载波调制激光雷达技术在海洋探测中的应用. 激光技术, 2008, 32(4): 346-349.
    [18] 马泳林宏冀航董天临 . 基于边缘探测技术的激光单稳频指标分析. 激光技术, 2007, 31(5): 469-472.
    [19] 黄沛陈博李岩 . 基于相敏检波和位置探测器的位移检测系统. 激光技术, 2017, 41(2): 235-239. doi: 10.7510/jgjs.issn.1001-3806.2017.02.018
    [20] 马聪叶志生张贵忠任国权撒昱 . 医用激光器闪光灯的准谐振触发预燃技术. 激光技术, 2009, 33(2): 191-194.
  • 加载中
图(8) / 表(1)
计量
  • 文章访问数:  3942
  • HTML全文浏览量:  2528
  • PDF下载量:  14
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-11-16
  • 录用日期:  2019-05-16
  • 刊出日期:  2020-03-25

基于位置敏感探测器的轨距在线监测系统研究

    作者简介: 刘子英(1964-),女,教授,研究方向为检测与监测控制。E-mail:151696128@qq.com
  • 华东交通大学 电气与自动化工程学院,南昌 330013
基金项目:  国家自然科学基金资助项目 51767006江西省应用研究培育计划资助项目 20181BBE58015江西省重点研发计划资助项目 20181BBE50019江西省教育厅科学技术资助项目 180310

摘要: 为了实时了解铁路轨道的轨距变化情况,研制了一种实时在线监测轨距变化的测量系统。采用位置敏感探测器和激光发射器作为测量元件实时采集轨距的变化信息,建立数学模型计算轨距的偏移量,通过组合算法来获取两根钢轨移动状态后的偏移量,并通过二次插值算法对位置敏感探测器进行非线性补偿以提高位置敏感探测器的精度和准确度;采用控制器局域网络总线和通用分组无线服务技术将一定区域内的轨距变化信息发送到监控中心进行显示;通过和现有二级铁路轨距尺分别对同一模拟轨道进行了测量实验。结果表明,该系统的测量精度为0.32mm,达到了目前轨道检测行业所规定的1mm精度要求。

English Abstract

    • 铁路轨道基础设施是列车运行的基础,轨道的状态检测有利于保障列车顺利通行。铁路轨道作为列车运行的基础,经常会受到列车轮对的横向冲击,使得钢轨横移,增加了两钢轨的间距,导致轮轨之间的关系发生变化,严重时会造成列车发生脱轨或爬轨事故。

      虽然JJG 219-2015《标准轨距铁路轨距尺》[1]中对轨尺做了一个标准限定,但是采用轨尺检测轨距的方法需要人工现场检测,劳动强度大,且不具有实时性。目前主要采用非接触式测量方式对轨距进行检测。国外铁路发达国家在轨检车上安装光发射器向钢轨上发射光束,并通过光接收器或摄像头接收钢轨上反射回来的光信息计算轨距[2]。但是这种方法需要轨检车在轨道上运行才能对轨距进行检测,不具有实时性。我国借鉴了国外检测轨距的先进技术,并研发了适用于国内轨道的轨距检测系统。参考文献[3]中提出了一种钢轨轮廓特征曲线来检测轨距的方法,采用两个激光摄像式传感器,并建立了空间姿态非线性标定模型[3]。但是该方法对激光摄像式传感器的安装要求严格,车体振动或摆动都会增加该系统的测量误差。参考文献[4]和参考文献[5]中提出了一种基于图像处理技术的轨距检测方法,采用4台CCD摄像机采集钢轨半截面轮廓的图像信息,并通过图像分析、立体视觉算法来计算轨距值[4-5]。但是该方法容易受到车体振动的影响,增加了图像处理的难度。参考文献[6]和参考文献[7]中采用激光技术与相机结合成一个视觉轨距检测系统,通过摄像机接收激光扫描后的钢轨轮廓来测量轨距[6-7],但是该方法采用的阈值分割法处理的图像会丢失大量的信息,造成测量误差偏高。参考文献[8]和参考文献[9]中提出了一种将2维激光传感器安装在铁轨两栖综合检测车上对轨距进行测量的系统[8-9],采用激光三角原理对采集到的数据进行处理,实现了非接触式动态实时检测,但是该方法需要移动检测车才能对轨距进行测量。参考文献[10]和参考文献[11]中提出了一种基于元模型的轨距估算方法,将所测得的轨道、车轮轮廓的数据创建测量情景对轨距进行测量,但是该方案是针对一种车型设计出来的,对于其它车型的模拟情景需要另设模型[10-11]

      本文中提出了一种基于位置敏感探测器(position sensitive detector, PSD)的轨距在线监测系统,该系统通过确定照射到PSD上的激光束位置来计算轨距的变化量,并将检测到的轨距实时变化信息传输到监控中心进行显示。该装置的探测器铺设在轨道和接触网支柱上,不受车体振动的影响,且安装简单、实时性强、检测精度高。

    • 利用照射到PSD上的激光束会改变其光电流分布的特性[12],将PSD安装在接触网支柱上,激光发射器安装在钢轨上,当激光发射器随着钢轨偏移时,照射到PSD上的光斑位置也发生了变化,再根据光斑引起的光电流输出值来计算光斑移动位置,进而计算出钢轨偏移量。铁路沿线的监测节点之间采用控制器局域网络(controller area network, CAN)总线连接起来,每220个监测节点中间设置中央处理器,对各个节点采集到的轨距变化信息进行处理分类,并采用通用分组无线服务技术(general packet radio service, GPRS)将处理后的信息发送至监控中心,提醒工作人员轨距变化超限的位置及轨距变化量,设置在监测节点上的指示灯也会跟随着轨距变化是否超出阈值进行指示。当轨距超出预设值时,控制器控制指示灯亮起,表明该处轨距变化超出了阈值,需要修整,铁路职工到达现场后,根据指示灯的状态即可知道需要修整的轨距位置。图 1为轨距监测系统结构示意图。

      Figure 1.  Structure diagram of gauge monitoring system

    • 安装在轨底上的激光发射器向安装在接触网支柱上的PSD发射激光束,当激光发射器随着钢轨移动时,打到PSD上的光斑位置也发生了变化。PSD上光斑的移动位置反映在了钢轨移动位置上,通过计算两股钢轨的偏移量既可计算出轨距的变化量,图 2为轨距测量原理图。

      Figure 2.  Gauge measurement schematic

      激光发射器1、激光发射器2分别发射出的激光束到达PSD1和PSD2之间的距离AB, CD分别为S1, S2。采用相位式激光测距原理来计算出S1S2[13]

      设定钢轨校准时激光发射器1和激光发射器2到接触网支柱之间的水平距离AO, CO分别为x1x2,激光发射器1的发射光与水平方向上的夹角α之间的关系为:

      $ \cos \alpha = \frac{{{x_1}}}{{{S_1}}} = \frac{{{x_1}}}{{{\lambda _0}\left( {N + \Delta N} \right)}} $

      (1)

      式中,x1为激光发射器1到接触网支柱的距离,S1为激光发射器1发出的激光束的距离,λ0为调制光波波长,N为相移的整周期数,ΔN为相移不足整周期数的余数[14-16]

      同理可以得出激光发射器2的发射光与水平方向上的夹角β之间的关系。由于钢轨在水平移动和竖直方向上移动时,水平夹角不发生改变,所以确定了x1x2的值即确定了αβ的值。

      设定钢轨校准时激光发射器1和激光发射器2打到PSD1和PSD2上的点O1O2为初始点。由于1号钢轨和2号钢轨会发生同样的移动,所以先对1号钢轨移动状态进行研究,2号钢轨同理。1号钢轨向右移动或向上凸起时,激光发射器1打到PSD1上的光斑位置高于O1点,PSD1上的光斑移动距离为:

      $ {f_1}\left( y \right) = \int_B^{{B_1}} k {\rm{d}}y $

      (2)

      式中,k为光斑从y1y1′的单位量,本设计中取PSD1的测量精度,BO1的纵坐标,B1为激光发射器1打到PSD1上的光斑向上移动后的纵坐标。

      1号钢轨上凸的移动量即为光斑在PSD1上的移动距离,向右移动量为:

      $ {f_1}\left( x \right) = \frac{{{f_1}\left( y \right)}}{{\tan \alpha }} = \frac{{\int_B^{{B_1}} k {\rm{d}}y}}{{\tan \left\{ {\arccos \left[ {\frac{{{x_1}}}{{{\lambda _0}\left( {N + \Delta N} \right)}}} \right]} \right\}}} $

      (3)

      1号钢轨向左移动或沉降时,激光发射器1打到PSD1上的光斑位置低于O1点,PSD1上的光斑移动距离为:

      $ {f_2}\left( y \right) = \left| {\int_B^{{B_2}} k {\rm{d}}y} \right| $

      (4)

      式中,B2为激光发射器1打到PSD1上的光斑向下移动后的纵坐标。

      1号钢轨沉降的移动量即为光斑在PSD1上的移动距离,向左移动量为:

      $ {f_2}\left( x \right) = \frac{{{f_2}\left( y \right)}}{{\tan \alpha }} = \frac{{ - \int_B^{{B_2}} k {\rm{d}}y}}{{\tan \left\{ {\arccos \left[ {\frac{{{x_1}}}{{{\lambda _0}\left( {N + \Delta N} \right)}}} \right]} \right\}}} $

      (5)

      同理可以得出2号钢轨上凸的移动量f3(y),向右移动量f3(y),向左移动量f4(y),沉降的移动量f4(y)。所以轨距的偏移量Δx根据1号钢轨和2号钢轨的移动状态来进行组合计算。根据轨距定义:轨距是铁轨(钢轨)顶面16mm范围内两股钢轨作用之间的最小距离可以得出,所以Δx为1号、2号钢轨变化后的位移偏移量。采用组合算法来计算1号钢轨和2号钢轨移动状态后的偏移量,图 3为轨距的偏移量Δx计算值。

      Figure 3.  Calculated value of gauge offset Δx

      图 3可知,计算其中一个组合的轨道偏移量,其它的同理可得。本设计中选取1号钢轨的右移上凸和2号钢轨的右移上凸的组合Δx1进行计算。图 4为Δx1的计算模型。

      $ \begin{array}{*{20}{c}} {{L_1}(x) = \sqrt {{{\left[ {{f_1}(y) - {f_3}(y)} \right]}^2} + {{\left[ {\left( {{x_1} - {x_2}} \right) - {f_3}(x) + {f_1}(x)} \right]}^2}} = }\\ {\sqrt {{{\left[ {\int_{{y_1}}^{y_1^\prime } k {\rm{d}}y - \int_{{y_2}}^{y_2^\prime } k {\rm{d}}y} \right)}^2} + \left[ {\left( {{x_1} - {x_2}} \right) - \frac{{\int_{{y_2}}^{{{y'}_2}} k {\rm{d}}y}}{{\tan \left\{ {\arccos \left[ {\frac{{{x_2}}}{{{\lambda _0}(N + \Delta N)}}} \right]} \right\}}} + \frac{{\int_{{y_1}}^{{{y'}_1}} k {\rm{d}}y}}{{\tan \left\{ {\arccos \left[ {\frac{{{x_1}}}{{{\lambda _0}(N + \Delta N)}}} \right]} \right\}}}} \right]} } \end{array} $

      (6)

      $ \Delta {x_1} = {L_1}\left( x \right) - L $

      (7)

      Figure 4.  Calculation model of Δx1

      式中,L1(x)为1号、2号钢轨变化后的轨距,L为两钢轨校准后的轨距。

    • 位置敏感探测器(PSD)是一种对入射光敏感器件,照射到器件上的光产生光电流,通过计算光电流的分布来确定器件敏感面的光的位置,进而计算光移动的距离[17]X1, Y1, X2, Y2为PSD的4个电极管脚,具有一定能量的光子入射到PSD的光敏面上的任意一点(X, Y)时,光电流在光敏面上以面电流的形式扩散,最终流向PSD的4个电极管脚,根据4个电极管脚上的电流来计算光敏面上光斑的移动位置[18-19],即:

      $ \left\{ {\begin{array}{*{20}{l}} {X = {L_x}\frac{{\left( {{I_{{X_1}}} + {I_{{Y_2}}}} \right) - \left( {{I_{{X_2}}} + {I_{{Y_1}}}} \right)}}{{2\left( {{I_{{X_1}}} + {I_{{Y_1}}} + {I_{{X_2}}} + {I_{{Y_2}}}} \right)}}}\\ {Y = {L_y}\frac{{\left( {{I_{{X_1}}} + {I_{{Y_1}}}} \right) - \left( {{I_{{X_2}}} + {I_{{Y_2}}}} \right)}}{{2\left( {{I_{{X_1}}} + {I_{{Y_1}}} + {I_{{X_2}}} + {I_{{Y_2}}}} \right)}}} \end{array}} \right. $

      (8)

      式中,LxLy为PSD的边长,XY为入射光点在PSD光敏面上的坐标。IX1, IY1, IX2, IY2分别为X1, Y1, X2, Y2上的输出电流。由于激光发射器1和激光发射器2的光信号调制频率为f0,且PSD的4个电极管脚上采集到的信号均为交流信号,所以采用相敏检波技术[20]对频率为f0的激光在PSD的光敏面上产生的交流信号进行解调,即在PSD的处理电路上添加相敏检波电路,只对激光发射器的光信号进行解调,进而得出光斑的位置。该技术只解调预设好频率为f0的激光信号,可以消除其它光源的干扰,使其信噪比和抗干扰能力得以提高。图 5为2维PSD探测原理图。

      Figure 5.  2-D PSD detection schematic

      钢轨校准时PSD1和PSD2上的初始点O1, O2的光斑坐标为(X1, Y1),(X2, Y2),那么PSD1上的入射光点从(X1, Y1)向上移动到(X1′, Y1′)、向下移动到(X1″, Y1″)的距离即分别为激光发射器1打到PSD1上的光斑移动的距离f1(y)和f2(y),同理(X2, Y2)向上移动到(X2′, Y2′)、向下移动到(X2″, Y2″)的距离分别为f3(y)和f4(y)。考虑到钢轨在竖直方向上发生的位移很微小,且对轨距产生的影响非常小,所以PSD上的X方向的变化忽略不计。从而可得:

      $ \left\{ {\begin{array}{*{20}{l}} {{f_1}(y) = \int_{{Y_1}}^{Y_1^\prime } k {\rm{d}}y}\\ {{f_2}(y) = \left| {\int_{{Y_1}}^{Y_1^{\prime \prime }} k {\rm{d}}y} \right|}\\ {{f_3}(y) = \int_{{Y_2}}^{Y_2^\prime } k {\rm{d}}y}\\ {{f_4}(y) = \left| {\int_{{Y_2}}^{Y_2^{\prime \prime }} k {\rm{d}}y} \right|} \end{array}} \right. $

      (9)
    • 由于PSD容易受除目标物以外的背景光和电路暗电流的影响, 导致PSD器件输出信号的精度和可靠性降低[21-22]。此外受其制造工艺和感光面积等因素的影响,导致它的光敏面上阻性不均匀,导致PSD输出存在非线性,这一缺陷增加了PSD的测量误差。所以本设计中采用二次插值算法对PSD进行非线性补偿,提高PSD的测量精度和数据的置信度。

      二次插值算法是一种用于一元函数f(α)在确定的初始区间内搜索极小点的一种方法。利用直线逼近代替初始区间内的实际曲线[23],对PSD的线性误差进行校正。PSD表面是一种分布许多离散点形成的矩形网络形状。任意给定一个测量点M(x, y),取最靠近M点的9个定点的方阵作为补偿网格,如图 6所示。

      Figure 6.  Compensation grid for quadratic interpolation algorithm

      取任意一个矩形区域A1A2A3A4这4个顶点对应的误差函数值分别为Ex(A1),Ex(A2),Ex(A3),Ex(A4), 其中A1点的插值函数为:

      $ {\varphi _1}(x,y) = \frac{{\left( {{x_{i + 1}} - x} \right)\left( {{y_{j + 1}} - y} \right)}}{{\left( {{x_{i + 1}} - {x_i}} \right)\left( {{y_{j + 1}} - {y_j}} \right)}} $

      (10)

      式中, xiyj为网格点的实际测量值。同理可以得出其它3点的插值函数。令该矩形内任意一点(x, y)的误差函数为Ex(x, y)。那么同理可以通过A1, A2, A3, A4, A5, A6, A7, A8, A9这9个顶点的误差函数的线性组合来表示M(x, y)点的位置误差函数Ex(x, y)和Ey(x, y),即:

      $ \begin{array}{*{20}{c}} {{E_x}(x,y) = \sum\limits_{q = j}^{j + 2} {\sum\limits_{p = i}^{i + 2} {\prod\limits_{k = 1,k \ne p}^{i + 2} {\left( {\frac{{x - {x_k}}}{{{x_p} - {x_k}}}} \right)} } } \times }\\ {\prod\limits_{l = j,l \ne q}^{j + 2} {\left( {\frac{{y - {y_l}}}{{{y_q} - {y_l}}}} \right)} {E_x}\left( {{x_p},{y_q}} \right)} \end{array} $

      (11)

      $ \begin{array}{*{20}{c}} {{E_y}(x,y) = \sum\limits_{q = j}^{j + 2} {\sum\limits_{p = i}^{i + 2} {\prod\limits_{k = 1,k \ne p}^{i + 2} {\left( {\frac{{x - {x_k}}}{{{x_p} - {x_k}}}} \right)} } } \times }\\ {\prod\limits_{l = j,l \ne q}^{j + 2} {\left( {\frac{{y - {y_l}}}{{{y_q} - {y_l}}}} \right)} {E_y}\left( {{x_p},{y_q}} \right)} \end{array} $

      (12)

      式中, p, q表示A1~A4这9个顶点最接近M(x, y)的点所在的位置,即顶点的坐标点的角标; k, l表示与位置p, q最接近的位置。其中,把φ(p, q, k, l)=$\prod\limits_{k=1, k \neq p}^{j+2}\left(\frac{x-x_{k}}{x_{q}-x_{k}}\right) \prod\limits_{l=j, l \neq q}^{j+2}\left(\frac{y-y_{l}}{y_{q}-y_{l}}\right)$称为系数函数。

      当补偿网格是均匀的时候,补偿网格的9个顶点的坐标可以用两个1维矩阵来表示[23],即:X=[X1, X2, X3], Y=[Y1, Y2, Y3]。由X1X2X3Y1Y2Y3这6个坐标值来确定系数函数φ(p, q, k, l)。可根据下式取点阵来作为二次插值算法的插值:

      $ \left\{ {\begin{array}{*{20}{l}} {\left| {{x_{i + 1}} - x} \right| < \Delta x}\\ {\left| {{y_{j + 1}} - y} \right| < \Delta y} \end{array}} \right. $

      (13)

      $ \left\{ {\begin{array}{*{20}{l}} {{x_{i + 2}} = {x_{i + 1}} + \Delta x\sin x}\\ {{y_{j + 2}} = {y_{j + 1}} + \Delta y\sin y}\\ {{x_i} = {x_{i + 1}} - \Delta x\sin x}\\ {{y_j} = {y_{j + 1}} - \Delta y\sin y} \end{array}} \right. $

      (14)

      式中,Δx,Δy为PSD光敏面xy方向上的步距。且二次插值算法的插值余项为:

      $ \begin{array}{*{20}{c}} {{R_n}(x,y) = }\\ {\frac{1}{{(n + 1)!}}\left( {\frac{{{\partial ^{n + 1}}E}}{{\partial {x^{n + 1}}}} + \frac{{{\partial ^{n + 1}}E}}{{\partial {y^{n + 1}}}}} \right){\phi _{n + 1}}(x,y)} \end{array} $

      (15)

      式中,ϕn+1(x, y)为构建出来的插值函数,因为A2, A3, A4都在过A4的两条直线上面,所以可以构造ϕn+1(x, y)如下:

      $ {\phi _1}(x,y) = {C_1}\left( {{x_{i + 1}} - x} \right)\left( {{y_{j + 1}} - y} \right) $

      (16)

      式中,$C_{1}=\left[\left(x_{i+1}-x_{i}\right)\left(y_{j+1}-y_{j}\right)\right]^{-1}$。

      由于ExEy的每一阶导数可近似小于1%,取样间隔为1mm,那么误差估计值为:

      $ \begin{array}{*{20}{c}} {{R_2}(x,y) = \frac{1}{{(2 + 1)!}} \times (1\% + 1\% ) \times }\\ {1 \times 1 = 3.33 \times {{10}^{ - 4}}{\rm{mm}}} \end{array} $

      (17)

      因此采用二次插值算法对PSD进行非线性补偿,大大提高了PSD的精度, 有利于检测轨距细微的变化。

    • 为了验证系统的检测轨距的精度和实时性,采用本系统和现有标准轨尺在学校实验基地的轨道上进行测量实验分析。由于基地钢轨不能移动,没法测量轨距变化后的偏移量,所以采用滑轨来模拟钢轨移动,记录每次移动滑块后的距离与标准轨尺检测两滑块之间的距离做对比,进而计算出轨距监测系统的精度。图 7为轨距监测系统测试场景图。

      Figure 7.  Test scenario of rail gauge monitoring system

      轨尺采用二级铁路轨距尺,《标准轨距铁路轨距尺》TB/T 1924-2008,规格为GJC-JG0,轨距尺的轨距测量范围为1427mm ~1470mm,轨距测量精度为0.20mm。任意选取长度为200m的两股钢轨作为检测对象,并选取8个不同位置分别进行多次实验,且轨尺和系统测试的位置、次数相同,得到本系统和轨尺测得的数据如表 1所示。

      Table 1.  Comparison between the ruler of railway gauge and the system

      test location measurement method average value/mm difference/mm standard deviation/mm
      1 ruler of railway gauge 1434.14 -0.35 0.09
      system 1434.49 0.30
      2 ruler of railway gauge 1435.18 -0.29 0.09
      system 1435.47 0.29
      3 ruler of railway gauge 1433.17 -0.37 0.10
      system 1433.54 0.29
      4 ruler of railway gauge 1436.17 0.61 0.12
      system 1435.46 0.27
      5 ruler of railway gauge 1433.17 -0.31 0.11
      system 1433.48 0.27
      6 ruler of railway gauge 1435.19 -0.34 0.10
      system 1435.53 0.32
      7 ruler of railway gauge 1434.19 -0.28 0.10
      system 1434.47 0.31
      8 ruler of railway gauge 1434.18 0.22 0.09
      system 1434.40 0.28

      通过对比分析了本系统与轨尺分别测量的8组实验数据, 可以得出轨距检测系统与轨尺对同一位置的轨距测量平均差值最大可达到0.61mm,轨距检测系统的标准差最大为0.32mm。

    • 通过MATLAB对第一处位置测得的数据进行曲线拟合,以此来对比分析轨距检测系统与轨尺在测量过程中产生的误差。图 8为轨距检测系统与轨尺所测数据拟合曲线。从图中可以看出,轨尺检测的数据点比较集中,拟合曲线趋于平缓,其中包含了读数误差和由于轨尺与钢轨发生机械摩擦造成的误差以及读数误差。而轨距检测系统检测到的数据点相对比较发散,使得它的拟合曲线相对变化比较快,这是由于轨距检测系统在做实验过程中采用了滑轨来模拟钢轨移动,使得测量过程中出现了人为误差和机械误差。由于本系统采用相敏检波技术消除了背景光的影响,所以背景光对轨距检测系统造成的影响是很小的,可以忽略不计。再从表 1中可以得出,轨距检测系统的标准差最大为0.32mm,而轨尺的标准差最大达到了0.12mm,明显轨距检测系统的误差接近于标准轨尺的测量误差,表明轨距检测系统具有实时性强、检测精度高等特点。

      Figure 8.  Curved curve of rail gauge measurement system and the measured data

    • 本文中设计了一套采用PSD和激光发射器相结合的轨距在线监测系统,在每根接触网支柱上安装两个相同的PSD,并在每根接触网支柱下方对应的两股钢轨轨底上安装相同的激光发射器,通过激光发射器在PSD上移动的位置来计算轨距的变化量,并采用CAN总线技术将一定区域的轨距变化量信息汇集起来,经GPRS技术发送至监控中心显示。利用二次插值算法对PSD的非线性误差进行补偿,将PSD的测量精度提高到3.33×10-4mm。最后通过对比分析轨距检测系统与轨尺检测同一轨道的数据得出,轨距检测系统的误差较低,精度翻了1倍,标准差最大为0.32mm,达到了目前轨道检测行业所规定的1.00mm精度要求。

参考文献 (23)

目录

    /

    返回文章
    返回