李書亮 劉成龍 鄭子天

(1.港珠澳大橋管理局，廣東珠海 519015;2.西南交通大學(xué)地球科學(xué)與環(huán)境工程學(xué)院，四川成都 610031;3.中鐵二院工程集團有限責(zé)任公司，四川成都 610031)

軌道幾何形位［1］是指軌道各部分的幾何形狀、相對位置和基本尺寸，包括軌距、水平、軌向以及外軌超高等。軌道的幾何形位與機車車輛輪對的幾何尺寸必須密切配合，因而軌道幾何形位的控制對于保證列車運行安全至關(guān)重要。隨著鐵路列車提速及高速鐵路技術(shù)的廣泛應(yīng)用，為保持高速列車運行的平穩(wěn)性和舒適性，必須對軌道的幾何形位實行嚴(yán)格控制。

實際線路設(shè)計過程中，軌道幾何形位與設(shè)計里程相對應(yīng)。在軌道幾何形位的檢測過程中，實測軌道上某點的幾何形位，如果能夠得到該點對應(yīng)的中線里程，便可得到該點設(shè)計的幾何形位，之后計算實測幾何形位與設(shè)計幾何形位之間的差值，這些差值是軌道靜態(tài)幾何形位檢測的主要指標(biāo)。因此，精確計算軌道測點對應(yīng)的線路中線里程，是實現(xiàn)對軌道靜態(tài)幾何形位進行檢測的前提條件。

目前，國內(nèi)只有少數(shù)文獻介紹軌道點實測坐標(biāo)與其線路里程之間函數(shù)關(guān)系模型的建立方法［2-3］，文獻資料［3］利用二次函數(shù)的極大極小值定理推導(dǎo)出軌道點坐標(biāo)與其線路里程之間嚴(yán)密的數(shù)學(xué)模型，并可用程序自動計算。本文則利用平面線形函數(shù)上任一點處切線斜率與其法線斜率之積為常數(shù)的特性，推導(dǎo)出了軌道點坐標(biāo)與其線路里程之間嚴(yán)密的數(shù)學(xué)模型，并利用某高速鐵路平面線形的設(shè)計參數(shù)和軌道點坐標(biāo)，編寫程序自動計算其線路里程。

1 軌道點里程計算模型

軌道的平面線形一般包括直線、緩和曲線和圓曲線［1］，為了實現(xiàn)計算軌道測點對應(yīng)的中線里程，首先要判斷軌道測點所在的線形。

1.1 判斷軌道點所在位置

軌道點里程是指軌道點所對應(yīng)的到設(shè)計線路中線上最近點的里程［3］。如圖1所示，點P為線路五大樁點(或者直線點)，→PX為線路前進方向，點O為在軌道控制網(wǎng)CPⅢ下實測的任意軌道點，點B為設(shè)計線路中線上距離點O最近的點，α為直線間的夾角。根據(jù)里程的定義，點B的里程即為軌道點O的里程。

通過五大樁點P(或直線點P，其坐標(biāo)可通過設(shè)計文件得到)與點O(通過實測得到其坐標(biāo))之間的坐標(biāo)關(guān)系，可以計算PO的距離和α，根據(jù)最短距離原則，可定位距離點O最近的五大樁點P(或直線點P)。再根據(jù)α的取值可以判斷點O是在點P的里程增加方向還是減小方向，即可判斷軌道點所在的線形。具體判斷方法如下:

圖1 軌道點與其在中線上投影點間關(guān)系

①若α＞90°，點B位于點P里程減小方向，并標(biāo)記定位因子fq=－1;

②若α=90°，點B與點P重合，點P的里程即為軌道點O的里程;

③若α＜90°，點B位于點P里程增加方向，并標(biāo)記fq=1。

根據(jù)上述計算和點P屬性(直緩點、緩圓點、圓緩點等)，可以判斷點B所處位置的線路線形。例如，當(dāng)點P的屬性為緩圓點并且α＞90°，那么點B就位于第一緩和曲線段上。

1.2 軌道點里程計算原理

如圖1所示，假設(shè)過點P(xp，yp)和B(xB，yB)直線的線形函數(shù)為 f(x，y)，當(dāng)點 O(xO，yO)到函數(shù) f(x，y)有最近點B時，一定滿足以下條件

式(1)中，kPB=

即過點B的切線斜率與其法線斜率的乘積為－1。

由于軌道點O和五大樁點P坐標(biāo)已知，通過式(1)可以得到點B的坐標(biāo)，進一步可以計算得到點B與點P之間的距離，再根據(jù)點P的已知里程，可以得到點B的里程，也就是軌道點O的里程。下面將按此思路計算軌道點里程的方法稱為“斜率法”。

1.3 軌道點里程計算方法

線路平面線形主要由直線段、圓曲線段和緩和曲線段組成，以下分別就這三種線形介紹斜率法計算軌道點里程的詳細方法。

(1)點B位于直線段上

若點B所在的線段為直線段，則函數(shù)f(x，y)可以表示為

把式(2)代入式(1)，有

因為點B(xB，yB)是直線方程f(x，y)上的點，因此有:

聯(lián)立式(3)與式(4)求解

求出B點的坐標(biāo)后，可以按照下式計算BP的距離dBP

則直線段上點B的里程LB為點P的里程LP與fq×dBP之和，即

(2)點B位于圓曲線段上

若點B所在的線段是圓曲線段，則首先建立以緩圓點為坐標(biāo)原點和過緩圓點切線為X軸的圓曲線獨立坐標(biāo)系，圓曲線函數(shù)f(x，y)的參數(shù)方程可以表示為

把點O(xO，yO)的坐標(biāo)通過坐標(biāo)轉(zhuǎn)換的方法旋轉(zhuǎn)到圓曲線坐標(biāo)系，然后將式(8)代入式(1)，有

則

把式(10)代入式(8)可求得點B的坐標(biāo)，同時有

則圓曲線段上點B的里程LB為點P的里程LP與fq× l之和［4］，即

(3)點B位于緩和曲線段上

若點B所在的線段是緩和曲線段，則首先建立以緩直點為坐標(biāo)原點和過緩直點切線為X軸的緩和曲線獨立坐標(biāo)系，緩和曲線函數(shù)f(x，y)的高次項參數(shù)方程［5－6］可以表示為

式(13)中，l為曲線長，l0為緩和曲線長度。

同樣的道理，把點O(xO，yO)的坐標(biāo)通過坐標(biāo)轉(zhuǎn)換的方法旋轉(zhuǎn)到緩和曲線坐標(biāo)系，然后將式(13)代入式(1)有

將式(14)化簡并忽略l的高次項后有

式(15)中，(xO，yO)為軌道點在緩和曲線獨立坐標(biāo)系中的坐標(biāo)，可通過坐標(biāo)轉(zhuǎn)換得到［7］。

式(15)可以采用簡單迭代、Newton法與弦截法等方法求解［8-9］。當(dāng)|li+1－li|＜ ξ，即第 i+1 與 i的差值小于給定的限差時，停止迭代;把得到的l代入式(16)即可以得到緩和曲線段上點B的里程LB為點P的里程LP與fq×l之和。

2 實際應(yīng)用

依據(jù)推導(dǎo)的數(shù)學(xué)模型與算法，編制了相應(yīng)的軌道實測點對應(yīng)的中線里程計算程序［10］。程序在讀取設(shè)計文件(如表1所示)后，通過輸入線路附近任意軌道點的坐標(biāo)，就能夠計算出線路中線上距離其最近的點的設(shè)計坐標(biāo)及其里程。

表1 某高速鐵路部分曲線設(shè)計要素

將表1的數(shù)據(jù)還原成線路圖，在線路圖上任取一點并過該點作線路的垂線，即可獲得垂點的設(shè)計坐標(biāo)，將垂點的理論坐標(biāo)與計算出的坐標(biāo)進行對比(如表2所示)。

表2 垂點計算設(shè)計坐標(biāo)與計算坐標(biāo)對比分析

表2的結(jié)果說明，垂點計算值與垂點設(shè)計值基本一致，也驗證了數(shù)學(xué)模型的正確性。有了垂點準(zhǔn)確的坐標(biāo)，就可以進一步計算垂點的里程。

3 結(jié)論

(1)利用平面線形函數(shù)上任一點處切線斜率與其法線斜率積為常數(shù)求取里程的方法在數(shù)學(xué)理論上是嚴(yán)密的，且在直線段和圓曲線段上尋找的最近點直接就是理論解，無需迭代計算;在緩和曲線段上，僅需要少量迭代就可以尋找到最優(yōu)數(shù)值解。

(2)利用點與點之間的距離、點與點之間的方位和線路走向的角度等幾何關(guān)系判斷目標(biāo)點在線路上的位置，避免了人為判斷與迭代求解的多值性問題。

(3)本文提出的里程計算方法簡單有效，有利于編制程序?qū)崿F(xiàn)，給高速鐵路軌道測量生產(chǎn)帶來較大便利。

［1］郝瀛．鐵道工程［M］．北京:中國鐵道出版社，2002

［2］李全信，張予東．直線與線路曲線交點坐標(biāo)及里程的通用解算方法［J］．鐵路航測，2002(1):124

［3］鄭子天，劉成龍，曹體濤．軌道點線路里程的一種計算方法［J］．測繪科學(xué)技術(shù)學(xué)報，2013(5):448451

［4］李全信．確定地面點與線路中線相對關(guān)系的統(tǒng)一數(shù)學(xué)模型．［J］．測繪通報，2002(8):3437

［5］王兆祥．鐵道工程測量［M］．北京:中國鐵道出版社，2008．

［6］高淑照，史東晏．緩和曲線非線性方程的快速解算［J］．測繪通報，2006(3):2830

［7］鄭子天．客運專線軌檢儀研制過程中的關(guān)鍵算法研究［D］．成都:西南交通大學(xué)，2009:2233．

［8］王景海．用簡化牛頓法進行緩和曲線段坐標(biāo)反算解算［J］．鐵道勘察，2009(1):1517

［9］許曦，劉慶元，余加勇．基于牛頓法的緩和曲線加樁計算［J］．測繪通報，2004(4):4345

［10］張秋隼．曲線中點和邊線點的坐標(biāo)計算程序［J］．鐵道勘察，2004(2):5961