CN107144234A - 一种城轨车辆轮对踏面轮廓拟合方法 - Google Patents
一种城轨车辆轮对踏面轮廓拟合方法 Download PDFInfo
- Publication number
- CN107144234A CN107144234A CN201710266743.7A CN201710266743A CN107144234A CN 107144234 A CN107144234 A CN 107144234A CN 201710266743 A CN201710266743 A CN 201710266743A CN 107144234 A CN107144234 A CN 107144234A
- Authority
- CN
- China
- Prior art keywords
- mtd
- mtr
- mrow
- msubsup
- mtable
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Tires In General (AREA)
Abstract
本发明公开了一种城轨车辆轮对踏面轮廓曲线拟合方法。该方法为:步骤1,获取轮对踏面数据:利用激光位移传感器获取轮对踏面轮廓离散数据;步骤2,踏面轮廓数据预处理:将轮对踏面轮廓离散数据进行预处理,融合成轮对踏面轮廓;步骤3,确定待拟合区域:根据轮对踏面轮廓特征点的位置选取拟合区域;步骤4,插值处理:对拟合区域进行插值处理;步骤5,曲线拟合:对插值处理结果使用最小平方法进行曲线拟合;步骤6,获取轮对参数:根据轮对尺寸参数定义准则获取轮对参数。本发明方法增加了轮对踏面轮廓上的有效点数,提高了测量精度。
Description
技术领域
本发明属于交通安全工程技术领域,特别是一种城轨车辆轮对踏面轮廓曲线拟合方法。
背景技术
轮对在城轨车辆运行过程中的起着重要的作用,承载着整个列车的负荷。轮对是与钢轨耦合的关键部位,在列车行驶过程中受到撞击、循环应力、高低温等作用的影响,因此轮对故障时而发生,需要定期检查城轨车辆轮对的运行状况,保障列车的安全运营。为了实时掌握轮对尺寸变化情况,获得准确的轮对尺寸参数,国内外已开发出不同的基于激光位移传感器的轮对尺寸检测系统。由于激光位移传感器采集的轮对踏面数据是以离散点坐标呈现的,需要借助离散点处理方法构造出平滑的轮对踏面曲线,精确地复现轮对踏面关键信息。采用不同的曲线拟合方法对最终的轮对尺寸检测结果有着重要的影响,曲线过拟合或欠拟合都会带来较大的检测误差。
针对轮对尺寸检测系统探测轮对踏面曲线拟合方法对检测系统检测结果影响的研究,已经取得了一些研究成果。哈尔滨工业大学的何平等人提出了一种基于最小二乘拟合的轮对轮缘在线检测系统,但该方法存在曲线过拟合、计算结果精度低的问题;南京理工大学的王晓浩等人提出一种基于拉格朗日乘数法的轮对外形轮廓线拟合方法,并利用实测数据进行了轮对踏面轮廓线拟合,但是该方法计算复杂,效率低,不适用于现场在线检测。
发明内容
本发明的目的在于提供一种精确、高效的城轨车辆轮对踏面轮廓曲线拟合方法,从而提高轮对尺寸检测系统的检测精度。
实现本发明目的的技术解决方案为:一种城轨车辆轮对踏面轮廓曲线拟合方法,包括以下步骤:
步骤1,获取轮对踏面数据:利用激光位移传感器获取轮对踏面轮廓离散数据;
步骤2,踏面轮廓数据预处理:将轮对踏面轮廓离散数据进行预处理,融合成轮对踏面轮廓;
步骤3,确定待拟合区域:根据轮对踏面轮廓特征点的位置选取拟合区域;
步骤4,插值处理:对拟合区域进行插值处理;
步骤5,曲线拟合:对插值处理结果使用最小平方法进行曲线拟合;
步骤6,获取轮对参数:根据轮对尺寸参数定义准则获取轮对参数。
本发明与现有技术相比,其显著优点是:(1)通过三次样条插值处理,增加轮对踏面轮廓上的有效点数,改善了处理数据的方法;(2)最小平方法寻找最优曲线作为踏面轮廓线,提高了测量精度。
附图说明
图1为本发明城轨车辆轮对踏面轮廓曲线拟合方法流程图。
图2为轮对轮缘区域划分图。
图3为插值处理轮对踏面曲线图。
图4为最小平方法拟合后轮对踏面曲线图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
结合图1,本发明城轨车辆轮对踏面轮廓曲线拟合方法,包括以下步骤:
步骤1,获取轮对踏面数据:利用激光位移传感器获取轮对踏面轮廓离散数据,具体为:在城轨车辆行进的轨道内外两侧低于钢轨位置分别安装第一激光位移传感器U1、第二激光位移传感器U2,该两个激光位移传感器与铅垂线均呈α角,第一激光位移传感器U1和第二激光位移传感器U2的激光发射面相交于同一平面,且该平面与水平面呈θ角,α、θ∈[42°,48°];当城轨车辆通过激光位移传感器时,第一激光位移传感器U1和第二激光位移传感器U2分别探测到轮对踏面的不同位置,获取的离散点坐标分别为其中i、j分别为U1、第二激光位移传感器U2测得的离散点数编号,且i=1,2,...,m,j=1,2,...,n,m、n分别是第一激光位移传感器U1、第二激光位移传感器U2测得的离散点个数。
步骤2,踏面轮廓数据预处理:将轮对踏面轮廓离散数据进行预处理,融合成轮对踏面轮廓,具体过程如下:
(2.1)由于激光位移传感器是有仰角安装,在传感器自身坐标系下采集的轮对踏面曲线也呈现倾斜状态,先将踏面离散数据进行坐标旋转,使传感器自身坐标系的纵轴平行于轮对内端面,根据式(1)、(2)进行离散数据坐标旋转:
此时两激光位移传感器探测数据的坐标分别为
(2.2)由于城轨车辆轮对有轴箱、刹车片和其他机械结构的干扰,需要对原始的踏面数据进行消噪处理,消除干扰物的影响;由于两个激光位移传感器采集的踏面数据互相独立,隔离在两个坐标系下,需要将两个激光传感器的探测数据融合在同一坐标系中,调节两坐标系中端面及踏面距离将第一激光位移传感器U1的数据融合到第二激光位移传感器U2的坐标系中,由公式(3)平移处理:
构成完整的踏面轮廓,各离散点坐标为(x(k),y(k)),其中k=m+n,d、l分别为U1所在坐标系距离U2所在坐标系沿x,y轴的距离大小分量。
步骤3,确定待拟合区域:由于在轮对尺寸计算过程中有效特征点仅在轮廓线某些区段上,不需要完全拟合整个轮对踏面,根据轮对踏面轮廓特征点的位置选取拟合区域,具体过程如下:
如图2所示,提取轮对内端面,以内端面曲线为基准线lb,第一轮缘待拟合区域为[l-d2,l-d1],其中d1∈(2,7),d2∈(17,25);第二轮缘拟合区域为[l-d3,l-d4],其中d3、d4根据实际轮缘情况取值,确定拟合区域;第一踏面拟合区域为[l-70-d5,l-70+d6],其中d5、d6∈(8,25)。
步骤4,插值处理:对拟合区域进行插值处理,具体如下:
(4.1)踏面轮廓曲线上区间[a,b]内各离散点均在函数上,在该区间内给定有关划分a=x0<x1<…<xp=b,且ξ(x)为该区间内每个子区间都满足三次多项式的条件函数;和ξ(x)满足以下条件:
1)
2)
3)
4)p为该区间的划分数目;
(4.2)ξ(x)有p分段,在计算过程中有4p个未知参数,插值条件和样条条件有4p-2个条件,借助两边界条件:
1)
2)δ0,δp,ρ0,ρp分别是离散点所在函数在x0,xp的一阶、二阶导数值;
(4.3)如图3所示,提提取满足插值函数ξ(x)的踏面轮廓插入点(x(i),y(i)),获取平滑的离散踏面曲线。
步骤5,曲线拟合:对插值处理结果使用最小平方法进行曲线拟合,具体如下:
由步骤4知轮对踏面轮廓的插值函数为ξ(x),待拟合区域的各离散点为(x(i),y(i)),设拟合函数为Pm(x)=a0+a1x+…+amxm。a0,a1,...,am为多项式拟合系数;根据最小平方法即寻找拟合曲线的多项式方程Pm(x),如图4所示,p为区间划分数目。
步骤6,获取轮对参数:根据轮对尺寸参数定义准则获取轮对参数;在步骤5获得轮对踏面拟合曲线的基础上,根据轮对踏面不同区段的特点分别提取第一、二轮缘基准点和第一踏面基准点的坐标,根据轮对尺寸参数的定义获取轮对尺寸参数。
下面结合具体实施例对本发明作进一步详细说明。
实施例1
以某地铁公司现场安装的一套轮对尺寸在线检测系统采集的轮对踏面数据为研究对象,获取多组踏面检测数据,以其中一组数据进行方法说明。轨道内、外侧的激光位移传感器采集的数据分别为如式(1)所示:
式中i=61,j=92,表明U1、U2这一帧分别采集到61、92个踏面数据点。将原始的踏面数据点进行坐标旋转,旋转角度为传感器安装倾角45°,获得的踏面点分别为如式(2)、(3)进行计算:
对旋转后的踏面点做滤除干扰处理后,进行两传感器坐标融合处理,如式(4)所示:
两传感器坐标融合处理后,重构出轮对踏面处切面曲线,获得在同一坐标系下的踏面轮廓坐标。然后根据基准点的位置,确定待插值处理的踏面区域。提取轮对内端面,以内端面曲线为基准线l=-268.3,第一轮缘待拟合区域为[l-22,l-5];第二轮缘拟合区域为[l-38,l-26];第一踏面拟合区域为[l-70-15,l-70+15],划分结果如图2所示。
在三个基准点位置区域进行插值处理,部分第一、二轮缘待拟合区域和第一踏面区域的插值函数系数矩阵如表1、2、3所示:
表1第一轮缘待拟合区域插值函数系数矩阵
表2第二轮缘待拟合区域插值函数系数矩阵
表3第一踏面待拟合区域插值函数系数矩阵
在对应的插值函数中提取插入轮对踏面的插值点,将插值点与原采集的踏面点共同组成一个踏面整体的离散点,利用最小平方法对拟合区域的离散点进行拟合处理,获得待拟合区域的拟合曲线,第一轮缘区域对应的拟合函数为:P3(x)=0.0x3+0.2x2+66.6x+8059.4,第二轮缘区域对应的拟合函数为:P3(x)=0.0x3-20x2-4890x-493340,第一踏面区域对应的拟合函数为:P3(x)=0.0x3-0.1x2-20.5x-2033.2,插值结果图像如图3所示。
对轮对踏面曲线进行拟合处理后,提取各基准点位置的坐标值,第一轮缘基准点为(-282.7132,303.4258),第二轮缘基准点为(-297.4342,322.2367),第一踏面处基准点为(-338.2721,332.2367)。根据轮对尺寸参数计算关系式获得轮对尺寸的各个参数。由该组踏面数据计算得出的轮缘高度为28.8110mm,轮缘厚度29.1621mm。
Claims (7)
1.一种城轨车辆轮对踏面轮廓曲线拟合方法,其特征在于,包括以下步骤:
步骤1,获取轮对踏面数据:利用激光位移传感器获取轮对踏面轮廓离散数据;
步骤2,踏面轮廓数据预处理:将轮对踏面轮廓离散数据进行预处理,融合成轮对踏面轮廓;
步骤3,确定待拟合区域:根据轮对踏面轮廓特征点的位置选取拟合区域;
步骤4,插值处理:对拟合区域进行插值处理;
步骤5,曲线拟合:对插值处理结果使用最小平方法进行曲线拟合;
步骤6,获取轮对参数:根据轮对尺寸参数定义准则获取轮对参数。
2.根据权利要求1所述的城轨车辆轮对踏面轮廓曲线拟合方法,其特征在于,步骤1所述的利用激光位移传感器获取轮对踏面轮廓离散数据,具体为:在城轨车辆行进的轨道内外两侧低于钢轨位置分别安装第一激光位移传感器U1、第二激光位移传感器U2,该两个激光位移传感器与铅垂线均呈α角,第一激光位移传感器U1和第二激光位移传感器U2的激光发射面相交于同一平面,且该平面与水平面呈θ角,α、θ∈[42°,48°];当城轨车辆通过激光位移传感器时,第一激光位移传感器U1和第二激光位移传感器U2分别探测到轮对踏面的不同位置,获取的离散点坐标分别为其中i、j分别为U1、第二激光位移传感器U2测得的离散点数编号,且i=1,2,...,m,j=1,2,...,n,m、n分别是第一激光位移传感器U1、第二激光位移传感器U2测得的离散点个数。
3.根据权利要求1所述的城轨车辆轮对踏面轮廓曲线拟合方法,其特征在于,步骤2所述的将轮对踏面轮廓离散数据进行预处理,融合成轮对踏面轮廓,具体如下:
(2.1)先将踏面离散数据进行坐标旋转,使传感器自身坐标系的纵轴平行于轮对内端面,根据式(1)、(2)进行离散数据坐标旋转:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>x</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>x</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>x</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>x</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
此时两激光位移传感器探测数据的坐标分别为
(2.2)对原始的踏面数据进行消噪处理,调节两坐标系中端面及踏面距离将第一激光位移传感器U1的数据融合到第二激光位移传感器U2的坐标系中,由公式(3)平移处理:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>x</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>x</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>x</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>x</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mi>d</mi>
<mo>,</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>y</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mi>l</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
构成完整的踏面轮廓,各离散点坐标为(x(k),y(k)),其中k=m+n,d、l分别为U1所在坐标系距离U2所在坐标系沿x,y轴的距离大小分量。
4.根据权利要求1所述的城轨车辆轮对踏面轮廓曲线拟合方法,其特征在于,步骤3所述的根据轮对踏面轮廓特征点的位置选取拟合区域,具体如下:
提取轮对内端面,以内端面曲线为基准线lb,第一轮缘待拟合区域为[l-d2,l-d1],其中d1∈(2,7),d2∈(17,25);第二轮缘拟合区域为[l-d3,l-d4],其中d3、d4根据实际轮缘情况取值,确定拟合区域;第一踏面拟合区域为[l-70-d5,l-70+d6],其中d5、d6∈(8,25)。
5.根据权利要求1所述的城轨车辆轮对踏面轮廓曲线拟合方法,其特征在于,步骤4所述的对拟合区域进行插值处理,具体如下:
(4.1)踏面轮廓曲线上区间[a,b]内各离散点均在函数上,在该区间内给定有关划分a=x0<x1<…<xp=b,且ξ(x)为该区间内每个子区间都满足三次多项式的条件函数;和ξ(x)满足以下条件:
1)
2)
3)
4)p为该区间的划分数目;
(4.2)ξ(x)有p分段,在计算过程中有4p个未知参数,插值条件和样条条件有4p-2个条件,借助两边界条件:
1)
2)δ0,δp,ρ0,ρp分别是离散点所在函数在x0,xp的一阶、二阶导数值;
(4.3)提取满足插值函数ξ(x)的踏面轮廓插入点(x(i),y(i)),获取平滑的离散踏面曲线。
6.根据权利要求1所述的城轨车辆轮对踏面轮廓曲线拟合方法,其特征在于,步骤5所述的对插值处理结果使用最小平方法进行曲线拟合,具体如下:
轮对踏面轮廓的插值函数为ξ(x),待拟合区域的各离散点为(x(i),y(i)),设拟合函数为Pm(x)=a0+a1x+…+amxm,a0,a1,...,am为多项式拟合系数;根据最小平方法即获取拟合曲线多项式方程Pm(x),p为区间划分数目。
7.根据权利要求1所述的城轨车辆轮对踏面轮廓曲线拟合方法,其特征在于,步骤6所述的根据轮对尺寸参数定义准则获取轮对参数,具体如下:
在步骤5获得轮对踏面拟合曲线的基础上,根据轮对踏面不同区段的特点分别提取第一、二轮缘基准点和第一踏面基准点的坐标,根据轮对尺寸参数的定义获取轮对尺寸参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710266743.7A CN107144234A (zh) | 2017-04-21 | 2017-04-21 | 一种城轨车辆轮对踏面轮廓拟合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710266743.7A CN107144234A (zh) | 2017-04-21 | 2017-04-21 | 一种城轨车辆轮对踏面轮廓拟合方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107144234A true CN107144234A (zh) | 2017-09-08 |
Family
ID=59774567
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710266743.7A Pending CN107144234A (zh) | 2017-04-21 | 2017-04-21 | 一种城轨车辆轮对踏面轮廓拟合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107144234A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109242776A (zh) * | 2018-09-11 | 2019-01-18 | 江苏君英天达人工智能研究院有限公司 | 一种基于视觉系统的双车道线检测方法 |
CN109885913A (zh) * | 2019-01-18 | 2019-06-14 | 北京工业大学 | 拟合滚刀磨损轮廓的方法、装置、设备及存储介质 |
CN110689250A (zh) * | 2019-09-20 | 2020-01-14 | 深圳供电局有限公司 | 一种用户用电数据的处理方法、系统及计算机可读介质 |
CN113028966A (zh) * | 2021-03-10 | 2021-06-25 | 哈尔滨市科佳通用机电股份有限公司 | 车轮踏面擦伤及不圆度在线动态检测方法 |
CN114396897A (zh) * | 2022-01-19 | 2022-04-26 | 成都理工大学 | 一种岩石结构面粗糙度高精度测绘装置及其测绘方法 |
CN114426039A (zh) * | 2022-03-14 | 2022-05-03 | 浙江师范大学 | 一种基于城轨列车车轮踏面的点云数据处理方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1183451A (ja) * | 1997-09-09 | 1999-03-26 | Sumitomo Metal Ind Ltd | 鉄道用車輪踏面の形状測定装置 |
CN101619965A (zh) * | 2009-08-14 | 2010-01-06 | 北京领邦仪器技术有限公司 | 轮对自动测量装置 |
CN104163188A (zh) * | 2014-08-28 | 2014-11-26 | 南京理工大学 | 一种城轨列车轮对踏面轮廓线的获取方法及装置 |
CN105043248A (zh) * | 2015-04-29 | 2015-11-11 | 南京理工大学 | 基于激光位移传感器的城轨车轮参数测量装置及方法 |
CN105203552A (zh) * | 2015-09-18 | 2015-12-30 | 南京拓控信息科技有限公司 | 一种360°踏面图像检测系统及其检测方法 |
CN106017354A (zh) * | 2016-07-26 | 2016-10-12 | 广州地铁集团有限公司 | 列车轮对踏面三维轮廓自动化测量方法和系统 |
-
2017
- 2017-04-21 CN CN201710266743.7A patent/CN107144234A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1183451A (ja) * | 1997-09-09 | 1999-03-26 | Sumitomo Metal Ind Ltd | 鉄道用車輪踏面の形状測定装置 |
CN101619965A (zh) * | 2009-08-14 | 2010-01-06 | 北京领邦仪器技术有限公司 | 轮对自动测量装置 |
CN104163188A (zh) * | 2014-08-28 | 2014-11-26 | 南京理工大学 | 一种城轨列车轮对踏面轮廓线的获取方法及装置 |
CN105043248A (zh) * | 2015-04-29 | 2015-11-11 | 南京理工大学 | 基于激光位移传感器的城轨车轮参数测量装置及方法 |
CN105203552A (zh) * | 2015-09-18 | 2015-12-30 | 南京拓控信息科技有限公司 | 一种360°踏面图像检测系统及其检测方法 |
CN106017354A (zh) * | 2016-07-26 | 2016-10-12 | 广州地铁集团有限公司 | 列车轮对踏面三维轮廓自动化测量方法和系统 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109242776A (zh) * | 2018-09-11 | 2019-01-18 | 江苏君英天达人工智能研究院有限公司 | 一种基于视觉系统的双车道线检测方法 |
CN109242776B (zh) * | 2018-09-11 | 2023-04-07 | 江苏君英天达人工智能研究院有限公司 | 一种基于视觉系统的双车道线检测方法 |
CN109885913A (zh) * | 2019-01-18 | 2019-06-14 | 北京工业大学 | 拟合滚刀磨损轮廓的方法、装置、设备及存储介质 |
CN110689250A (zh) * | 2019-09-20 | 2020-01-14 | 深圳供电局有限公司 | 一种用户用电数据的处理方法、系统及计算机可读介质 |
CN110689250B (zh) * | 2019-09-20 | 2023-04-14 | 深圳供电局有限公司 | 一种用户用电数据的处理方法、系统及计算机可读介质 |
CN113028966A (zh) * | 2021-03-10 | 2021-06-25 | 哈尔滨市科佳通用机电股份有限公司 | 车轮踏面擦伤及不圆度在线动态检测方法 |
CN113028966B (zh) * | 2021-03-10 | 2021-11-12 | 哈尔滨市科佳通用机电股份有限公司 | 车轮踏面擦伤及不圆度在线动态检测方法 |
CN114396897A (zh) * | 2022-01-19 | 2022-04-26 | 成都理工大学 | 一种岩石结构面粗糙度高精度测绘装置及其测绘方法 |
CN114426039A (zh) * | 2022-03-14 | 2022-05-03 | 浙江师范大学 | 一种基于城轨列车车轮踏面的点云数据处理方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107144234A (zh) | 一种城轨车辆轮对踏面轮廓拟合方法 | |
CN104915963B (zh) | 一种用于plcc元件的检测与定位方法 | |
CN107607044B (zh) | 一种基于激光位移传感器的车轮踏面磨耗检测方法 | |
CN106969749B (zh) | 一种地铁隧道断面形变的检测方法 | |
CN105946898B (zh) | 一种基于激光测距的城轨列车车轮直径检测方法及系统 | |
CN104163188B (zh) | 一种城轨列车轮对踏面轮廓线的获取方法及装置 | |
CN104228875B (zh) | 一种城轨列车轮对尺寸在线检测方法及装置 | |
CN102628682A (zh) | 基于工业ct扫描技术的体积测量方法 | |
CN107248156A (zh) | 分段点自动提取的踏面轮廓拟合方法 | |
CN103295232B (zh) | 基于直线和区域的sar图像配准方法 | |
CN105783779B (zh) | 基于三层匹配的钢轨轮廓实时形态识别与失真校准方法 | |
CN105292180A (zh) | 基于多种传感器的非接触式轮对尺寸在线检测方法和装置 | |
CN106091951A (zh) | 一种城轨列车轮缘参数在线检测系统及方法 | |
CN107101594B (zh) | 一种提取轮轨车轮空间轮缘最低点的方法 | |
CN102749061A (zh) | 基于动态模板的钢轨磨耗测量方法 | |
AU2022349296A1 (en) | Road surface technical condition detection method and device based on three-dimensional contour | |
CN103839262A (zh) | 一种基于直线和fft的sar图像配准方法 | |
CN104501719A (zh) | 一种确定列车轮对轴线的方法 | |
CN105835902B (zh) | 一种基于激光位移传感器的车轮直径的检测方法 | |
CN104008257A (zh) | 一种针对具有复杂曲面的构件的逆向设计方法 | |
CN110232682A (zh) | 一种基于图像的轨道异物探测方法 | |
CN106813569A (zh) | 一种基于线结构光的汽车轮胎三维定位方法 | |
CN107202543A (zh) | 城轨列车车轮直径的检测方法 | |
CN114240916A (zh) | 钢轨外观状态的多偏振光点云数据融合方法及装置 | |
CN112414318B (zh) | 基于结构光钢轨磨耗快速测量算法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170908 |