CN110781575A - 一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法 - Google Patents

一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法 Download PDF

Info

Publication number
CN110781575A
CN110781575A CN201910846936.9A CN201910846936A CN110781575A CN 110781575 A CN110781575 A CN 110781575A CN 201910846936 A CN201910846936 A CN 201910846936A CN 110781575 A CN110781575 A CN 110781575A
Authority
CN
China
Prior art keywords
length
sequence
chord
chords
determining
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
CN201910846936.9A
Other languages
English (en)
Other versions
CN110781575B (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.)
Jiangxi Agricultural University
Original Assignee
Jiangxi Agricultural University
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 Jiangxi Agricultural University filed Critical Jiangxi Agricultural University
Priority to CN201910846936.9A priority Critical patent/CN110781575B/zh
Publication of CN110781575A publication Critical patent/CN110781575A/zh
Application granted granted Critical
Publication of CN110781575B publication Critical patent/CN110781575B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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 With Unspecified Measuring Means (AREA)

Abstract

一种用于确定钢轨波磨多中点弦测量仪弦的数量及长度的方法，包括：根据空间采样率和逆滤波窗口长度确定有效弦长组合的方法；根据测量目标波段确定波长点子序列的方法；根据有效弦长组合及波长点子序列确定最优的弦的数量及长度的方法。本发明在现有方法的基础上，考虑了空间采样率和逆滤波窗口长度的影响，避免了测量时不响应盲区的出现，提高了计算效率，适合基于弦测法的钢轨波磨测量装置设计和制造过程中确定弦的数量及长度使用。

Description

一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法
技术领域
本发明涉及一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法，其直接应用于设计和生产基于多中点弦测量原理的钢轨波磨检测装置。
背景技术
目前，基于多中点弦的钢轨波磨检测装置已经能够实现快速、准确的进行波磨测量，其基本原理见下式：
Figure 152750DEST_PATH_IMAGE001
其中m为测量弦的数量，Σm为m个中点弦测量值的代数和，g m (x)表示第m根弦的中点弦测量值，Lm表示第m根测量弦的长度，f(x)表示轨道不平顺函数。由于采用了多个中点弦测量值的代数和作为最终测量结果，幅值增益-波长响应曲线得到有效改善。但是，目前基于多中点弦测量原理的钢轨波磨检测装置所用的弦长及组合方式的选择仅仅以幅值增益-波长响应曲线在给定待测波长内方差最小为目标，忽略了空间采样率和逆滤波窗口长度的影响，这种方法不但运算量极大，在测量弦的数目较多（大于2根）、待测目标波段较长（10mm-1000mm）时很难得到最优结果，而且当改变空间采样率和逆滤波窗口长度时会有可能出现对特殊波长波磨不响应的问题，无法得到正确的测量结果。鉴于钢轨波磨病害检测的重要性，优化基于多中点弦测量原理的钢轨波磨检测装置中弦的数量及长度的计算方法，既是检测理论上的创新，亦是实际工程的需要。
发明内容
本发明的目的是克服现有技术的不足，提出一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法，本发明是通过以下技术方案实现的。
（1）确定空间采样频率f s 以及逆滤波窗口的长度N。
根据下式得到序列λ i ：
Figure 941452DEST_PATH_IMAGE002
其中，f s =1/e，e为设定的采样间隔，N为逆滤波窗口包含的采样点数，n为1到N的整数序列，λ i 为n取不同值时对应的波长点。
（2）设定计算所需要的参数
取序列λ i 中大于10mm且小于1000mm的部分得到子序列λ j （j=1,2,3…），设定弦的数量m取值范围[Up，Uq]、每根弦L m 的长度取值范围[Wp，Wq]。
（3）确定有效的弦的数量m和每根弦的长度L m 的组合
将选定的弦数目m、每根弦的长度L m以及λ j 序列中的每个值代入到式(b)中计算，
Figure 821683DEST_PATH_IMAGE003
若对每个弦长组合，不论取λ j 序列中的何值均有至少一个tm≠0（m=1,2,3…），则将此时弦的数量m,每根弦的长度Lm标记为一个有效组合y，否则在取值范围[Up，Uq]内改变m的值或在取值范围[Wp，Wq]中改变每根弦Lm的长度得到一个新的弦长组合后再重新计算，由此可得到全部有效组合y k，k=1,2,3…。
（4）确定最优的组合Y
对全部有效组合y k 中每一个，按照式(c)分别计算λ取λ j序列中各值时的幅值增益H得到幅值增益序列H k ’ ，按照式(d)计算H k ’ 的方差S k 2 ：
Figure 555601DEST_PATH_IMAGE005
式(d)中S k 2 表示H k ’ 序列的方差值，x v表示H k ’ 序列中第v点的值，x表示H k ’ 序列中所有点的均值，z表示所得H k ’ 序列的总长度。
对得到的全部有效组合的S k 2 进行分析，取S k 2 最小时的组合y k 作为最优的组合Y，此时弦的数量以及各根弦的长度即为最优。
本发明所述算法是自动计算，或者是手动计算。
本发明在现有确定钢轨波磨多中点弦测量仪弦的数量及长度方法基础上添加了计算有效弦长及组合方式的算法，克服了现有方法计算最优弦长数量及长度时不考虑采样率及逆滤波窗口长度导致出现不响应盲区的缺陷，提高了计算效率，适合于基于弦测法的钢轨波磨测量仪设计和制造过程中确定弦的数量及长度使用。
附图说明
附图1为一种用于确定钢轨波磨多中点弦测量仪弦的数量及长度的方法的流程。
具体实施方式
本发明将通过以下实施例结合附图1作进一步说明。
根据实际需要，设定采样间隔为1mm，逆滤波窗口长度为8192个点，代入到式(a)中计算可得λi序列，保留λi序列中大于10mm且小于1000mm的部分得到子序列λ j ；设定弦的数目m=2，弦的数量变化范围为[2,3],弦长L变化范围为[10mm,1000mm], 将不同的弦的数量m和各弦的不同长度取值Lm进行组合，代入到式(b)中进行计算，得到全部有效组合y k ；根据式(c)可得到m=2时双中点弦幅值增益波长响应函数(e)和m=3时的三中点弦幅值增益波长响应函数(f)：
Figure 470248DEST_PATH_IMAGE006
Figure 255802DEST_PATH_IMAGE007
取上述计算得到的每一个有效组合y k 中的弦的数量m、每根弦的长度Lm，根据公式(e)和(f)分别计算λ取λ j序列各值时幅值增益H得到幅值增益序列H k ’ 后，按照式(d)计算H k ’ 的方差Sk 2。根据计算结果可知当选取m=3，L1=98mm，L2=236mm，L3=315mm的组合时方差最小，并且在采样间隔为1mm，逆滤波窗口长度为8192个点时，对10mm-1000mm波长范围内的波磨不存在无响应的盲区，此即为最优的弦长数量及长度。

Claims (2)

1.一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法，其特征是：
（1）确定空间采样频率f s 以及逆滤波窗口的长度N
根据下式得到序列λ i ：
其中，f s =1/e，e为设定的采样间隔，N为逆滤波窗口包含的采样点数，n为1到N的整数序列，λ i 为n取不同值时对应的波长点；
（2）设定计算所需要的参数
取序列λ i 中大于10mm且小于1000mm的部分得到子序列λ j （j=1,2,3…），设定弦的数量m取值范围[Up，Uq]、每根弦L m 的长度取值范围[Wp，Wq]；
（3）确定有效的弦的数量m和每根弦Lm长度的组合
将选定的弦数目m、每根弦的长度Lm以及λ j 序列中的每个值代入到式(b)中进行计算，若对每个弦长组合，不论取λ j 序列中的何值均有至少一个tm≠0（m=1,2,3…），则将此时弦的数量m、每个弦的长度Lm标记为一个有效组合y，否则在取值范围[Up，Uq]内改变m的值或在取值范围[Wp，Wq]中改变每根弦Lm的长度得到一个新的弦长组合后再重新计算，由此可得到全部有效组合y k，k=0,1,2…；
Figure 604972DEST_PATH_IMAGE002
（4）确定最优的组合Y
对全部有效组合y k 中每一个，按照式（c）分别计算λ取λ j序列各值时的幅值增益H得到幅值增益序列H k ’ ，按照式（d）计算H k ’ 的方差S k 2 ；
Figure 365118DEST_PATH_IMAGE003
Figure 569834DEST_PATH_IMAGE004
式（d）中S k 2 表示H k ’ 序列的方差值，x v表示H k ’ 序列中第v点的计算值，x表示H k ’ 序列中所有点的均值，z表示所得H k ’ 序列的总长度；
对得到的全部有效组合的S k 2 进行分析，取S k 2 最小时的组合y k 作为最优的组合Y, 此时弦的数量以及各根弦的长度即为最优。
2.根据权利要求1所述的一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法,其特征是：所述计算为自动计算或手动计算。
CN201910846936.9A 2019-09-09 2019-09-09 一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法 Active CN110781575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910846936.9A CN110781575B (zh) 2019-09-09 2019-09-09 一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910846936.9A CN110781575B (zh) 2019-09-09 2019-09-09 一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法

Publications (2)

Publication Number Publication Date
CN110781575A true CN110781575A (zh) 2020-02-11
CN110781575B CN110781575B (zh) 2023-10-24

Family

ID=69383405

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910846936.9A Active CN110781575B (zh) 2019-09-09 2019-09-09 一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法

Country Status (1)

Country Link
CN (1) CN110781575B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116452774A (zh) * 2023-06-09 2023-07-18 西南交通大学 基于改进cpd点云配准方法的钢轨波磨三维重建方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101144714A (zh) * 2007-10-26 2008-03-19 北京航空航天大学 一种钢轨磨耗综合参数车载动态测量装置及方法
CN105128888A (zh) * 2015-08-27 2015-12-09 江西日月明测控科技股份有限公司 基于多中点弦同步测量的轨道波浪形磨耗快速检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101144714A (zh) * 2007-10-26 2008-03-19 北京航空航天大学 一种钢轨磨耗综合参数车载动态测量装置及方法
CN105128888A (zh) * 2015-08-27 2015-12-09 江西日月明测控科技股份有限公司 基于多中点弦同步测量的轨道波浪形磨耗快速检测方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
殷华 等: "《基于ZFFT的轨道波磨高精度》", 《铁道科学与工程学报》 *
殷华 等: "《基于ZFFT的轨道波磨高精度》", 《铁道科学与工程学报》, vol. 14, no. 9, 30 September 2017 (2017-09-30) *
殷华 等: "《基于中点弦测模型的钢轨波磨量值估计》", 《振动.测试与诊断》 *
殷华 等: "《基于中点弦测模型的钢轨波磨量值估计》", 《振动.测试与诊断》, no. 05, 15 October 2016 (2016-10-15) *
殷华 等: "《轨道踏面损伤弦测响应研究》", 《铁道科学与工程学报》 *
殷华 等: "《轨道踏面损伤弦测响应研究》", 《铁道科学与工程学报》, vol. 15, no. 2, 28 February 2018 (2018-02-28) *
殷华;朱洪涛;魏晖;王志勇;谭卿杰;: "基于中点弦测模型的钢轨波磨量值估计", 振动.测试与诊断, no. 05 *
魏珲;刘宏立;马子骥;李艳福;: "基于组合弦测的钢轨波磨广域测量方法", 西北大学学报(自然科学版), no. 02 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116452774A (zh) * 2023-06-09 2023-07-18 西南交通大学 基于改进cpd点云配准方法的钢轨波磨三维重建方法
CN116452774B (zh) * 2023-06-09 2023-08-18 西南交通大学 基于改进cpd点云配准方法的钢轨波磨三维重建方法

Also Published As

Publication number Publication date
CN110781575B (zh) 2023-10-24

Similar Documents

Publication Publication Date Title
Huss et al. Homogenization of long-term mass-balance time series
CN101825567A (zh) 一种近红外光谱和拉曼光谱波长的筛选方法
CN106769937B (zh) 一种光谱数据处理方法
EP2351996B1 (de) Verfahren zur Ermittlung einer Kenngröße für die Korrektur von Messwerten eines Coriolis-Massedurchflussmessgeräts
CN110781575A (zh) 一种确定钢轨波磨多中点弦测量仪弦的数量及长度的方法
CN106295865A (zh) 一种水稻产量的预测方法
CN113420412B (zh) 基于成像光谱的土壤有机碳含量连续深度分布提取方法
CN110553631B (zh) 一种关于水位流量关系的水位测量系列误差分析方法
CN105005978A (zh) 一种基于Savitzky-Golay滤波器参数优化的光谱实时滤波方法
CN117349683A (zh) 基于光谱数据的汽车配件涂装色差异常检测系统
CN110334473B (zh) 一种基于模量衰变的沥青路面结构相对弯沉的计算方法
CN113806943A (zh) 一种基于关键水情变量与植被固碳速率关系的湿地植被固碳速率预测方法
CN113486295A (zh) 基于傅里叶级数的臭氧总量变化预测方法
CN116522085A (zh) 一种全自动拉索频率提取、定阶和索力识别方法与应用
CN101006481A (zh) 用于确定在道路网中的转向率的方法
CN112762885A (zh) 一种基于监测数据的桥梁实时挠度校验系数计算方法
CN111896505A (zh) 一种原位、非破坏、便携式的古代青砖主要元素检测方法
CN109059875B (zh) 一种驱动完美模型开展月尺度径流预报的方法
CN108170951B (zh) 基于采样数据时空匹配示踪试验的纵向离散系数确定方法
CN115271506A (zh) 一种区域降雨表征长系列构建分析一体化方法
CN114383646B (zh) 一种连续变化型被测量传感器分辨力的检测方法和设备
CN115270637A (zh) 一种基于gbrt的地下排水管道最大应力预测方法
CN115310370A (zh) 耦合深度学习和物理机制的区域性植被蒸腾预测方法
CN114577360A (zh) 一种拉曼测温雷达信号自动分析校准及反演方法
CN110895628B (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
GR01 Patent grant
GR01 Patent grant