王源，徐金輝，陳嶸，肖杰靈，王平

(西南交通大學高速鐵路線路工程教育部重點實驗室，四川成都610031)

基于中點弦測法的軌道不平順精確值數(shù)學模型研究

王源，徐金輝，陳嶸，肖杰靈，王平

(西南交通大學高速鐵路線路工程教育部重點實驗室，四川成都610031)

弦測法是測量軌道不平順的一種基本方法，原理簡單，使用方便，高效迅捷。傳統(tǒng)觀點是直接將弦測值作為軌道不平順的近似描述，這會不可避免地因基準線變動而產(chǎn)生較大誤差。針對該問題建立了一個描述中點弦測法本質(zhì)的數(shù)學模型，分析了軌道不平順與其弦測值之間的關(guān)系，構(gòu)造了一種計算軌道不平順精確值的迭代算法與快速算法，并采用數(shù)值仿真對弦測過程進行模擬。結(jié)果顯示:迭代算法總體誤差較小，傳遞函數(shù)較好，但由于迭代次數(shù)等原因會產(chǎn)生端點誤差;快速算法以犧牲計算內(nèi)存為代價能達到較高精度，絕對誤差在1 μm以內(nèi)，傳遞函數(shù)效果極好，從而證明了所建立的數(shù)學模型的正確性與計算結(jié)果的精確性。

中點弦測法 軌道不平順 傳遞函數(shù) 數(shù)學模型

從20世紀中葉起，各國鐵路大都采用弦測法測量軌道不平順，以弦線作為測量的基準線。由于該方法具有測量原理簡單、使用方便、裝備便宜等優(yōu)點，一度得到世界范圍的廣泛應(yīng)用［1］。傳統(tǒng)觀點認為，弦測法的傳遞函數(shù)是隨著弦長與軌道不平順波長的比值變化的，有較嚴重的“缺陷”，只有在部分情況下才能正確測量或近似反映軌道的平順狀態(tài)［2-3］?；诖?，文獻［4］分析了三點偏弦的軌面復(fù)原方法，文獻［5］研究了四點弦軌面復(fù)原方法，文獻［6］提出“以小推大”公式以改善弦測法的適用性。然而這些方法均未能從根本上改變弦測法的本質(zhì)缺陷。

本文認為，所謂“缺陷”是由于直接使用弦測值作為軌道不平順的近似描述產(chǎn)生的，這將不可避免地因基準線變動而產(chǎn)生較大誤差。針對此，本文建立了弦測法的數(shù)學模型，從理論上解釋了軌道不平順與其弦測值之間的數(shù)學關(guān)系，發(fā)現(xiàn)通過弦測值這一相對信息可以反算得出真實的軌道不平順，最后借助數(shù)值仿真驗證了該數(shù)學模型的正確性與計算結(jié)果的精確性。

1 弦測法的數(shù)學模型

弦測法測量原理如圖1所示，在直線段軌道上，設(shè)真實的軌道不平順y=f(x)，x為里程。圖中，AB為測量弦線，長度為L，C為弦線的中點，則C處的弦測正矢值為CD，由于軌道不平順幅值為毫米級，AB弦與水平向的夾角θ足夠小，可以認為CD值與CD'值近似相等。

圖1 弦測法示意

考慮到實際的軌道不平順只在有限里程范圍內(nèi)才有意義，可將f(x)定義在［xmin，xmax］，為了方便后面的數(shù)學推導，這里將軌道不平順推廣到(-∞，∞)，則f(x)可以表示為

里程x處的弦測值g(x)可以表示為［1］

這里，控制測量步長為L/2，可以在每個L/2處得到一個弦測值。式(1)稱為從不平順f(x)到弦測值g(x)的正變換過程。然而，問題的關(guān)鍵在于如何利用g(x)逆推得到f(x)，將式(1)變換如下

從式(2)可知，里程x處的不平順值與弦測值g(x)、不平順值f(x-L/2)和f(x+L/2)有關(guān)，通過弦測法只能測得g(x)，然而f(x-L/2)與f(x+L/2)均未知，故由式(2)無法直接得到真實的不平順f(x)，但是由式(2)可構(gòu)造一個遞推式

進一步遞推有

依次遞推下去，可以得到第n階公式

式中f0(x)是迭代的初值。

通過數(shù)學歸納法可證明式(5)即為式(3)的第n步遞推結(jié)果。分析式(5)可以發(fā)現(xiàn)，fn(x)由兩部分組成，其中第一部分只與g(x)有關(guān)，第二部分只與f0(x)有關(guān)，并且隨著n值的逐漸增大，由于初值f0(x)為有界收斂序列，可以證明

結(jié)合式(5)、式(6)、式(7)可以得到

從式(8)能夠發(fā)現(xiàn)f(x)可以僅用g(x)求出，通過式(8)可構(gòu)造一種迭代算法，以實現(xiàn)通過弦測值求出較準確的不平順值。

2 軌道不平順的逆推算法

2.1 軌道不平順的迭代算法

圖2為軌道不平順迭代算法的示意圖，圖中AC， BD，CE為測量弦線，O，N分別為AC，BD、CE中點，P為MN中點。BM，CO，DN分別為B，C，D三處的L長弦測值。

從圖2中能夠看出，C點的真實不平順|CC'|為

然后將CP值作為C點新的不平順近似值，再用該近似值修正B，D兩點的不平順值，根據(jù)式(10)的原理對每一個點的弦測值進行修正，進而得到不平順新的近似值，如此修正多次后，所得不平順近似值將收斂于真實不平順值，如式(11)所示。

該過程即構(gòu)成一套迭代算法，在迭代有限次之后即可發(fā)現(xiàn)不平順近似值趨于穩(wěn)定，此時對比原不平順與該迭代計算所得的近似不平順，即可發(fā)現(xiàn)兩者誤差很小。

2.2 逆推軌道不平順的快速算法

由于迭代過程計算時間較長，本文還針對此設(shè)計了一個快速算法，借助矩陣運算跳過如上迭代過程。如果把不平順f(x)作離散采樣，轉(zhuǎn)化為向量T，采樣間隔L/2，把測得的弦測值記作向量G，則式(2)可以用矩陣的形式寫為

式中T為轉(zhuǎn)換矩陣，如式(13)所示

把式(12)稱為正變換方程，也即弦測值測量過程。通過G向量逆推求出F向量稱為反變換方程，這里需要找到一個逆轉(zhuǎn)換矩陣S，使得

將式(14)帶入式(12)可得

也即

式中，I為單位矩陣，因而S為T的逆矩陣。

該方法本質(zhì)是求解線性方程組，在該過程中需要儲存一個N階三對角矩陣，N為采樣數(shù)目。在矩陣維數(shù)N較小時可以通過求轉(zhuǎn)換矩陣的逆矩陣的方式迅速求解，在維數(shù)N較大時對系統(tǒng)內(nèi)存需求較多，因而只能采取上述占用內(nèi)存較少的迭代算法，或者可以考慮使用其他線性方程組求解方法。從這點看來，上文所述迭代算法本質(zhì)上即為雅克比迭代格式的一種描述。

3 弦測法的數(shù)值仿真過程

為了驗證上面所述方法的可行性，本文采用數(shù)值仿真方法，通過MATLAB設(shè)計程序模擬弦測過程，并將測得的弦測值通過本文所述迭代算法與快速算法反算軌道不平順。仿真過程分為三部分:第一部分采用數(shù)值仿真模擬弦測過程;第二部分根據(jù)弦測數(shù)據(jù)逆推真實不平順;第三部分對比分析逆推不平順的效果、傳遞函數(shù)等。仿真流程如圖3所示。

圖3 仿真流程

第一部分模擬弦測過程，設(shè)置弦長為L，采樣間隔為L/2，算法流程如圖4所示。第二部分逆推過程的算法流程如圖5所示。第三部分分別從時域上的誤差與頻域上的傳遞函數(shù)兩個角度出發(fā)分析計算結(jié)果的正確性與準確性。

4 實例驗證

4.1 弦測過程數(shù)值仿真

本文弦測法仿真使用的軌道不平順數(shù)據(jù)通過美國六級軌向不平順譜反演獲得。采用文獻［7］提出的方法，首先根據(jù)軌道不平順功率譜密度與頻譜幅值的關(guān)系，得出不平順的頻域幅值，并給出隨機相位;然后根據(jù)實數(shù)離散傅里葉變換的共軛對稱性，將頻譜擴展完整;最后通過傅里葉逆變換得到軌道不平順的模擬時域樣本。不平順樣本截止波長為1～200 m，采樣間隔為0.25 m。圖6為美國六級軌向不平順譜的反演樣本，圖7為反演樣本的模擬譜與美國六級軌向不平順譜的比較，可以看出模擬譜與美國六級軌向不平順譜幾乎完全重合。

圖4 弦測法測量過程算法流程

圖5 逆推軌道不平順的迭代算法流程

圖6 美國六級軌向不平順譜的反演樣本

圖7 反演樣本的模擬譜與美國六級譜比較

反演不平順樣本長度取為1 km，用固定0.5 m弦線以0.25 m的步長測量得到弦測值，弦測結(jié)果如圖8所示。

圖80 .5 m弦測值

4.2 逆推軌道不平順

分別采用迭代算法與快速算法對弦測值進行反演。迭代算法的計算精度與迭代步數(shù)有關(guān)，對于不同的測量弦長，要達到同樣的精度其迭代步數(shù)是不一樣的。對于固定的軌道長度，測量弦線越短所需迭代步數(shù)越多。另外，由于迭代算法存在端點效應(yīng)，本文控制精度選用的是中間80%范圍內(nèi)的測量數(shù)據(jù)。對1 000 m軌道長度，采用0.5 m測量弦長，并設(shè)置迭代步數(shù)為10萬次，其計算效果如圖9(a)所示，圖9(b)為絕對誤差曲線，在100～900 m范圍內(nèi)數(shù)據(jù)吻合得很好，而兩邊出現(xiàn)邊緣效應(yīng)，總體而言在100～900 m范圍內(nèi)誤差＜0.5 mm，精度很高。

圖9 迭代算法逆推不平順與真實不平順的對比

通過快速算法計算得到的軌道不平順與真實不平順的結(jié)果如圖10(a)所示，絕對誤差曲線如圖10(b)所示。可見在整段范圍內(nèi)誤差在1 μm以內(nèi)。快速算法計算過程占用的內(nèi)存遠多于迭代算法，在工程中測量點過多時宜采用迭代算法。

圖10 快速算法逆推不平順與真實不平順對比

4.3 傳遞函數(shù)比較

圖11(a)為直接使用10 m弦線測得弦測值的傳遞函數(shù)，圖11(b)為先測得0.5 m弦測值，再通過“以小推大”公式計算所得10 m弦測值所對應(yīng)的傳遞函數(shù)。由于通過美國六級譜反演的不平順截止波長為1～200 m，因而傳遞函數(shù)的有效范圍為1～200 m。從圖中可見傳遞函數(shù)總體較差，僅個別波長處傳遞函數(shù)為1，說明使用10 m弦測值不能很好地描述軌道的不平順狀態(tài)。

圖11兩種方法的傳遞函數(shù)

圖12 為兩種逆推算法的傳遞函數(shù)。圖12(a)所示的迭代算法計算得到的不平順的傳遞函數(shù)在1～200 m范圍內(nèi)在1上下波動，大部分值為1，有個別波動是由于邊界效應(yīng)所致。圖12(b)所示的快速算法傳遞函數(shù)在1～200 m范圍內(nèi)穩(wěn)定在1左右，這也驗證了快速算法高精度的特點，當然這種高精度是以占用更多的計算資源獲得的。

圖12 兩種遞推算法的傳遞函數(shù)

5 結(jié)論

傳統(tǒng)觀點直接將弦測值作為不平順的近似，這實質(zhì)上是對弦測數(shù)據(jù)的一種極大浪費，本文研究了一種新的方法，借助測得的弦測值逆推真實的軌道不平順，并且能夠達到足夠的精度，傳遞函數(shù)效果極好。目前，本文涉及的數(shù)學模型為理論模型，僅適用于不計測量誤差時軌道不平順的推算，對于包含測量誤差的不平順反演則需要進一步研究。

［1］羅林，張格明，吳旺青，等.輪軌系統(tǒng)軌道不平順狀態(tài)控制［M］.北京:中國鐵道出版社，2006.

［2］羅林.軌道隨機干擾函數(shù)［J］.中國鐵道科學，1982，3(1): 74-100.

［3］高建敏.基于狀態(tài)轉(zhuǎn)移概率矩陣的軌道不平順發(fā)展預(yù)測研究［J］.鐵道建筑，2011(7):140-143.

［4］程櫻，許玉德，周宇，等.三點偏弦法復(fù)原軌面不平順波形的理論及研究［J］.華東交通大學學報，2011，28(1):42-46.

［5］毛曉君，許玉德，周宇.基于四點弦測法的軌面不平順檢測及復(fù)原方法［J］.華東交通大學學報，2013，30(5):13-17.

［6］朱洪濤，魏暉，熊瑞文，等.弦測法檢測軌向不平順的研究［J］.鐵道建筑，2005(10):63-64.

［7］陳果，翟婉明.鐵路軌道不平順隨機過程的數(shù)值模擬［J］.西南交通大學學報，1999，34(2):138-141.

Research on mathematical model of accurate value of track irregularity based on midpoint chord measurement method

WANG Yuan，XU Jinhui，CHEN Rong，XIAO Jieling，WANG Ping
(MOE Key Laboratory of High-speed Railway Engineering，Southwest Jiaotong University，Chengdu Sichuan 610031，China)

Chord measuring method is a basic method for measuring track irregularity，its theory is easily understood and it is convenient，quick and effective.T raditional method is considering the chord measuring value as approximate track irregularity，which inevitably leads to the big error because of the base line changing.In order to solve this problem，a mathematical model describing the nature of midpoint chord measuring method was established，the relationship between track irregularity and the chord measuring value was discussed，an iterative algorithm and a fast algorithm for calculating accurate track irregularity value were designed，the chord measuring process was simulated by numerical simulation.T he results showed that the iterative algorithm has less overall error，the transfer function is good and the endpoint error occurs because of such reasons as number of iterations，while the fast algorithm achieves a high accuracy at the expense of computing memory，the absolute error of which is less than 1 μm，the effect of transfer function is excellent，which proves the correctness of the mathematical model and the accuracy of the calculation results.

M idpoint chord measuring method;T rack irregularity;T ransfer function;M athematical model

U216

A

10.3969/j.issn.1003-1995.2015.05.35

1003-1995(2015)05-0139-05

(責任審編李付軍)

2014-09-28;

2014-12-30

國家自然科學基金——高鐵聯(lián)合基金項目(U1234201)

王源(1991—)，男，四川南充人，碩士研究生。