熊江陵 曾凡云 唐俊成 韓 揚(yáng)

(中國(guó)電建集團(tuán)中南勘測(cè)設(shè)計(jì)研究院有限公司 長(zhǎng)沙 410014)

在傳統(tǒng)的公路及城市道路設(shè)計(jì)中，通常將平面交叉口右轉(zhuǎn)彎車(chē)道邊線(xiàn)簡(jiǎn)化為單圓曲線(xiàn)形式，一方面是為了便于施工，另一方面原因是雙心圓及三心圓等復(fù)曲線(xiàn)無(wú)法通過(guò)固定的公式進(jìn)行直接計(jì)算。但是在高等級(jí)公路及渠化城市道路設(shè)計(jì)中，將交叉口右轉(zhuǎn)彎車(chē)道邊線(xiàn)設(shè)計(jì)為雙心圓或三心圓復(fù)曲線(xiàn)具有一定的必要性，因?yàn)樗苡行岣哂肄D(zhuǎn)彎通行能力和行車(chē)舒適性。蔡偉等[1]研究了不同類(lèi)型彎道路緣石對(duì)車(chē)輛在交叉口轉(zhuǎn)彎處行駛速度和軌跡的影響，結(jié)果表明在交叉口彎道處采用三心圓曲線(xiàn)進(jìn)行布設(shè)，不但節(jié)約用地資源，也有利于提高車(chē)輛在交叉口轉(zhuǎn)彎處行車(chē)的舒適性。JTG D20-2017《公路路線(xiàn)設(shè)計(jì)規(guī)范》[2]在10.4.3條規(guī)定：渠化平面交叉的右轉(zhuǎn)彎車(chē)道，在內(nèi)側(cè)路面邊緣應(yīng)采用三心圓復(fù)曲線(xiàn)。文獻(xiàn)[3]指出以鉸接列車(chē)控制設(shè)計(jì)時(shí)，相交路面的邊緣應(yīng)采用復(fù)曲線(xiàn)。在三心圓復(fù)曲線(xiàn)的計(jì)算方法方面，相關(guān)研究較少，朱家兵等[4]根據(jù)相交道路線(xiàn)形，分左直右直、左直右曲、左曲右直、左曲右曲等情況，提出一種以切線(xiàn)長(zhǎng)度為逼近目標(biāo)的三心圓復(fù)曲線(xiàn)計(jì)算方法，但該方法在曲線(xiàn)相交的交叉口中計(jì)算誤差較大(1 m)，并且文獻(xiàn)中沒(méi)有對(duì)進(jìn)、出口道均為曲線(xiàn)的形式進(jìn)行詳細(xì)論證。

為研究三心圓復(fù)曲線(xiàn)的通用計(jì)算方法，基于三心圓復(fù)曲線(xiàn)終點(diǎn)坐標(biāo)計(jì)算公式，提出以終點(diǎn)距離為逼近目標(biāo)的二分迭代法，并對(duì)該方法收斂性和迭代初始值區(qū)間選取進(jìn)行論證。通過(guò)在實(shí)際項(xiàng)目中的應(yīng)用結(jié)果表明：該算法能正確求解任意線(xiàn)形相交的交叉口右轉(zhuǎn)彎車(chē)道三心圓復(fù)曲線(xiàn)，具有更好的通用性、穩(wěn)定性和求解精度，同時(shí)也具有較高計(jì)算效率。

1 三心圓復(fù)曲線(xiàn)切線(xiàn)長(zhǎng)公式及應(yīng)用

根據(jù)文獻(xiàn)[5]中介紹的方法，進(jìn)、出口道為直線(xiàn)時(shí)，左、右側(cè)切線(xiàn)長(zhǎng)Tq、Th分別為

Tq=(R1-R2)sinα+[R3-(R3-R2) cosγ]/

sinδ-[R1-(R1-R2)cosα]/tanδ

(1)

Th=(R3-R2)sinγ+[R1-(R1-R2) cosα]/

sinδ-[R3-(R3-R2)cosγ]/tanδ

(2)

式中：R1、R2、R3分別為3段圓弧的半徑，R2小于R1、R3；α、β、γ分別為3段圓弧對(duì)應(yīng)的轉(zhuǎn)角；δ為三段圓弧旋轉(zhuǎn)角度之和，各參數(shù)示意圖見(jiàn)圖1。

圖1 三心圓復(fù)曲線(xiàn)切線(xiàn)長(zhǎng)公式參數(shù)示意圖

假定3段圓弧長(zhǎng)度分別為L(zhǎng)1、L2、L3，則有公式

α=(L1/R1)·(180°/π)

(3)

β=(L2/R2)·(180°/π)

(4)

γ=(L3/R3)·(180°/π)

(5)

δ=α+β+γ

(6)

通過(guò)確定R1、R2、R3，并且確定3段圓弧中任意2段的長(zhǎng)度，或選擇3段圓弧相等，利用式(1)～(6)可求出Tq和Th，進(jìn)而可以計(jì)算出A、D的位置和確定3段圓弧的準(zhǔn)確位置。

當(dāng)進(jìn)、出口道為曲線(xiàn)時(shí)，文獻(xiàn)[4]以切線(xiàn)長(zhǎng)T值為目標(biāo)進(jìn)行迭代逼近，目標(biāo)控制條件為T(mén)h-Th′<10-10，其中Th為根據(jù)公式(2)計(jì)算得到的切線(xiàn)長(zhǎng)，Th′為迭代點(diǎn)處實(shí)際的切線(xiàn)長(zhǎng)。經(jīng)測(cè)試，該方法結(jié)果存在較大誤差，原因是路線(xiàn)的圓弧半徑較大，δ的微小變化對(duì)Th結(jié)果影響較小，但對(duì)Th′影響較明顯。

2 基于終點(diǎn)距離逼近的二分迭代法

為解決上述三心圓復(fù)曲線(xiàn)計(jì)算的誤差問(wèn)題，并考慮相交道路線(xiàn)形的任意性，提出以終點(diǎn)距離為逼近目標(biāo)的二分迭代法。

由于相交道路幾何條件的不確定性，對(duì)于曲線(xiàn)相交的交叉口，三心圓復(fù)曲線(xiàn)僅在某一范圍存在可能解，可以視為局部收斂的非線(xiàn)性問(wèn)題[6]。局部收斂非線(xiàn)性問(wèn)題求解，需要解決如下3個(gè)問(wèn)題。

1) 選擇合適的迭代公式。

2) 選擇合適的迭代初始值。

3) 保證迭代的收斂性。

2.1 選擇迭代公式：基于終點(diǎn)距離逼近的迭代方程

當(dāng)相交道路線(xiàn)形為曲線(xiàn)時(shí)(見(jiàn)圖2)，三心圓復(fù)曲線(xiàn)終點(diǎn)切向量及坐標(biāo)計(jì)算公式如下。

圖2 三心圓復(fù)曲線(xiàn)終點(diǎn)切向量及坐標(biāo)計(jì)算參數(shù)示意圖

Ve=R(Vs,δ)

(7)

Pe=Ps+R1·(υ1+υ1′)+

R2·(υ2+υ2′)+R3·(υ3+υ3′)

(8)

式中：Ve為復(fù)曲線(xiàn)終點(diǎn)向量；Vs為起點(diǎn)向量；R(Vs,δ)為將向量Vs順時(shí)針旋轉(zhuǎn)δ角度；Pe為復(fù)曲線(xiàn)終點(diǎn)坐標(biāo)；Ps為復(fù)曲線(xiàn)起點(diǎn)坐標(biāo)；υ1、υ2、υ3分別為3段圓弧起點(diǎn)處指向圓心的單元向量，υ1′、υ2′、υ3′分別為3段圓弧終點(diǎn)處背向圓心的單位向量。

三心圓復(fù)曲線(xiàn)在終點(diǎn)需與邊線(xiàn)相切，即應(yīng)滿(mǎn)足如下2個(gè)條件。

1) 復(fù)曲線(xiàn)終點(diǎn)切向量與終點(diǎn)在邊線(xiàn)上投影處的切向量相等。

2) 復(fù)曲線(xiàn)終點(diǎn)與終點(diǎn)在邊線(xiàn)上投影點(diǎn)距離為0。

用公式表示為

D(Ve,Ve′)≤ε

(9)

D(Pe,Pe′)≤ε

(10)

以式(10)為逼近目標(biāo)建立二分法迭代方程，其中式(10)應(yīng)以滿(mǎn)足式(9)為前提，他們的因變量分別為起點(diǎn)切向量Vs和Ps起點(diǎn)坐標(biāo)。

2.2 選擇迭代初始值區(qū)間

對(duì)于迭代方程為y=f(x)的二分迭代法，需要確定迭代求解的初始值區(qū)間[X1,X2]，保證[X1,X2]對(duì)應(yīng)的目標(biāo)值y1，y2滿(mǎn)足y1×y2≤0。對(duì)于式(9)，假設(shè)復(fù)曲線(xiàn)終點(diǎn)向量目標(biāo)解為Ve*，則初始值區(qū)間[Vs1,Vs2]應(yīng)滿(mǎn)足復(fù)曲線(xiàn)終點(diǎn)向量Ve1、Ve2分別在Ve*的兩側(cè)。對(duì)于式(10)，初始值區(qū)間[Ps1,Ps2]應(yīng)滿(mǎn)足復(fù)曲線(xiàn)終點(diǎn)坐標(biāo)Pe1、Pe2分別在出口車(chē)道邊線(xiàn)的兩側(cè)。

由于路線(xiàn)幾何的不確定性，可能收斂區(qū)間[Ps1,Ps2]僅存在于路線(xiàn)中很短的一段距離，為了快速、準(zhǔn)確地找出初始值區(qū)間[Ps1,Ps2]，參考文獻(xiàn)[4]，可以先利用式(1)、(2)以切線(xiàn)長(zhǎng)為逼近目標(biāo)迭代求解到復(fù)曲線(xiàn)起點(diǎn)Ps，然后再在Ps附近找出符合上述條件的[Ps1,Ps2]，此處不再贅述。

2.3 迭代收斂性論證

迭代方程為y=f(x)在區(qū)間[X1,X2]能否收斂，除需要滿(mǎn)足對(duì)應(yīng)的目標(biāo)值y1×y2≤0外，還應(yīng)滿(mǎn)足y=f(x)在區(qū)間[X1,X2]連續(xù)。

根據(jù)式(3)～(8)，可以得到

Pe=Ps+f(Vs,δ)

(11)

當(dāng)三心圓復(fù)曲線(xiàn)起點(diǎn)Ps確定時(shí)，Vs能唯一確定，而δ的大小也與起點(diǎn)位置有關(guān)，也就是說(shuō)Pe可以表示為

Pe=F(Ps)

(12)

為保證行車(chē)的順暢性，路線(xiàn)在幾何上也必然是連續(xù)且光滑的，也就是說(shuō)Ps是連續(xù)變化的，所以基于終點(diǎn)距離逼近的迭代方程D(Pe,Pe′)≤i也具備連續(xù)性。

分別對(duì)相交道路為直-圓、圓-直、圓-圓等情況進(jìn)行測(cè)試，終點(diǎn)距離D(Pe,Pe′)值迭代收斂過(guò)程圖見(jiàn)圖3。

圖3 二分法迭代計(jì)算收斂過(guò)程圖

3 以終點(diǎn)距離逼近為目標(biāo)的二分法迭代計(jì)算流程圖

根據(jù)上述分析，以終點(diǎn)距離逼近為目標(biāo)的二分迭代法需要如下過(guò)程。

1) 根據(jù)切線(xiàn)長(zhǎng)預(yù)估復(fù)曲線(xiàn)起點(diǎn)位置。

2) 在預(yù)估的起點(diǎn)附近確定迭代區(qū)間初始值。

3) 以終點(diǎn)距離逼近為目標(biāo)進(jìn)行二分迭代計(jì)算。

算法流程圖設(shè)計(jì)見(jiàn)圖4。

圖4 三心圓復(fù)曲線(xiàn)終點(diǎn)逼近二分迭代法計(jì)算流程圖

其中三心圓復(fù)曲線(xiàn)終點(diǎn)切向量迭代逼近的算法包含在“根據(jù)Ps計(jì)算三心圓終點(diǎn)Pe”這一步驟中，迭代步驟與計(jì)算Pe相同，不再贅述。

4 項(xiàng)目應(yīng)用

金山大道位于廣西省河池市，設(shè)計(jì)長(zhǎng)度2 350 m，雙向八車(chē)道，設(shè)計(jì)速度60 km/h，為城市主干路。其中金山大道與開(kāi)元大道相交的交叉口位于K1+062.705處，相交處道路主線(xiàn)段為半徑625m的圓弧，支線(xiàn)段為直線(xiàn)，最小交角約為78°，該右轉(zhuǎn)區(qū)間行車(chē)道邊線(xiàn)選擇三心圓復(fù)曲線(xiàn)，R1、R2、R3分別為25,15,25 m，3段圓弧長(zhǎng)度相等。

基于上述三心圓復(fù)曲線(xiàn)算法，利用三維道路軟件OpenRoads Designer進(jìn)行二次開(kāi)發(fā)，完成“交叉口右轉(zhuǎn)邊線(xiàn)”設(shè)計(jì)功能。經(jīng)測(cè)試，該功能可以正確生成該交叉口右轉(zhuǎn)彎車(chē)道三心圓復(fù)曲線(xiàn)，邊線(xiàn)圖見(jiàn)圖5，三維實(shí)體模型圖見(jiàn)圖6。

圖5 交叉口右轉(zhuǎn)彎車(chē)道邊線(xiàn)圖

圖6 交叉口三維實(shí)體模型圖

5 結(jié)語(yǔ)

1) 利用以終點(diǎn)距離逼近為目標(biāo)的二分迭代法，計(jì)算三心圓復(fù)曲線(xiàn)，計(jì)算結(jié)果精度高，能夠?qū)?fù)曲線(xiàn)起、終點(diǎn)誤差控制在10-8m以?xún)?nèi)。

2) 該算法具有通用性，能夠應(yīng)用于進(jìn)、出口車(chē)道為直-直、直-曲、曲-直、曲-曲等不同情形。

3) 該算法在進(jìn)、出口車(chē)道足夠長(zhǎng)的前提下，能保證收斂性，并且迭代計(jì)算效率較高。

3) 僅通過(guò)修改迭代目標(biāo)公式，該算法框架可擴(kuò)展應(yīng)用于雙心圓復(fù)曲線(xiàn)的計(jì)算。

4) 該算法在逼近方法、初始值選取、收斂性論證等方面，對(duì)道路路線(xiàn)中復(fù)雜復(fù)曲線(xiàn)計(jì)算有一定參考價(jià)值。