网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基本型平曲线“两点法”设计计算原理  PDF

  • 李玉华 1
  • 王奥凯 1
  • 刘佳音 2
  • 马洋洋 1
  • 吴树铭 1
1. 大连理工大学 建设工程学院,辽宁 大连 116023; 2. 大连海洋大学 海洋与土木工程学院,辽宁 大连 116023

中图分类号: U412.33U412.31

最近更新:2025-02-10

DOI:10.11908/j.issn.0253-374x.23223

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

为改善传统导线法存在导线与平曲线偏离较大等不足,提出“两点法”基本型平曲线形设计方法。仅需利用导线法中的直缓点与缓直点两个端点即可获得由缓和曲线及圆曲线共3段线元组成的基本型平曲线形,其平曲线总偏转角α、总切线长T为已知值。先采用迭代法确定缓和曲线长度LS (对称基本型)或采用夹逼法获得圆曲线半径R(非对称基本型)的理论最小值、最大值,并结合平面线形技术要求确定LSR的设计取值范围;然后基于LSR进行遍历循环,通过迭代法获得缓和曲线偏转角β,进而获得各线元其余参数。研究结果表明,“两点法”设计可获得系列“平曲线族”,线形分布规律,互不交叉重叠,且全部曲线位于两端细窄、中间宽厚的“特定设计区域”;LSR、圆曲线长LC及平曲线总长LH等参数呈规律性变化,其中LSRLC的增减趋势相反,LH平均值近似固定,非对称基本型的总切线长度之比须在(0.5,2.0)范围内方可有解。“两点法”突显坐标位置点的控制作用,线形布设高效,设计成果为特定区域内的系列平曲线族而非单条曲线,有益于结合实际地形及设计意图进行线形比选、优选和路线智能化设计。

近几十年来我国公路及铁路建设迅猛发展,已形成较为完善的路线平面线形设计方法,其中导线

1和曲线2应用最为广泛。导线法一般设计过程是先拟定直线导线,然后在相邻导线间设置由 “第1缓和曲线+圆曲线+第2缓和曲线”共3段线元构成的基本型平曲线形。导线法敷设便捷,简单易行,但存在导线与实际平曲线偏离较大等问题,较适合于平原微丘区或立交桥主线等简单线形设计。曲线法的一般设计方法是先布设圆曲线,然后采用直线、缓和曲线或圆曲线3种基本线元相互连接,由此构建复杂的平曲线形。曲线法依靠线元灵活布线、设计过程“所见即所得”,无赘余线形,在山岭重丘区、立交匝道等复杂线形设计中极具优势;但在调整线元参数时,其两端的切线方向易改变,由此将产生“牵连效应”。由于导线法依靠直线控制基本走向,所含3段平曲线线元参数的调整不会改变两端的切线方向,这正好弥补曲线法的不足。

在复杂地形及立交匝道线形设计实践中,路线设计方法进一步发展,如针对导线法改进的附合导线

3、五单元导线4等,针对曲线法改进的综合5、控制线元6等。这些方法尽管有利于线形优化,但增加了选线的复杂性。两点线元7基于已确定的终点坐标、切线方向角、曲线半径3个参数,仅需再拟定1个端点便可确定6类线元,同时生成10余条线元曲线,有益于线形比选及线形拟合重8-9

导线法与曲线法两种方法相结合是较好的创新思路,如曲直

2、直线约束通用方10、圆心连线闭合导线11等,但这些方法偏重于完善曲线法的优势;采用样条函12、多项式缓和曲线13等方法拓展了传统平曲线形类型,提高了线形曲率的连续14,但增加了设计计算难度。

综上,目前平面线形设计方法面临的主要问题有:①两种基本方法均存在固有不足;②设计成果通常为单一平曲线,不利于进行多种适宜线形的比选。两点线元法对此有所改进,但仍局限于曲线法优势。

近年来,随着高新技术及人工智能的发展,线形优

15和智能选线16备受关注。本文提出的“两点法17-18综合了导线法的便捷性与曲线法的灵活性,可获得规律性分布的“平曲线族”及“特定区域”的线形,便于设计比选,有助于路线优化和智能选线设计。

1 导线法基本型平曲线参数及计算

1.1 平曲线组成及参数

图1所示,导线法基本型平曲线参数包括:平曲线偏转角α,圆曲线半径R、圆心角αC及圆曲线长LC,第1、第2缓和曲线长LS1LS2及相应的缓和曲线偏转角β1β2,第1、第2总切线长TS1TS2,因设置缓和曲线而引起的切线长增加值q1q2及相应的圆曲线内移值ΔR1ΔR2,其中αRLS1LS2通常为已知参数;点D为导线交点,点C为圆曲线圆心,点A、点B分别为过圆心C向第1、第2导线作垂线的垂足;T1T2为不设缓和曲线(仅含圆曲线)时的名义切线长,其相应的圆曲线半径分别为|CA¯|=R+ΔR1、|CB¯|=R+ΔR2

图1  导线法基本型平曲线组成及参数

Fig. 1  Composition and parameters for basic horizontal curves using traverse method

1.2 平曲线参数计算

1.2.1 基本参数

按式(1)~ (3)计算基本参数β1β2q1q2ΔR1ΔR2,其中关于缓和曲线偏转角β的多项式函数 gxβ)、gyβ)、fxβ)、fyβ)如式(4)、(5)所示。

β1=LS1/(2R)β2=LS2/(2R) (1)
q1=Rgx(β1)q2=Rgx(β2) (2)
ΔR1=Rgy(β1)-RΔR2=Rgy(β2)-R (3)
gx(β)=2βfx(β)-sin βgy(β)=2β2fy(β)+cos β (4)
fx(β)=i=1(-1)i-1β2i-2(2i-2)!(4i-3)fy(β)=i=1(-1)i-1β2i-2(2i-1)!(4i-1) (5)

1.2.2 切线长

总切线长TS1TS2的计算公式如式(6)所示。

TS1=q1+T1=q1+R+ΔR2sin α-R+ΔR1tan αTS2=q2+T2=q2+R+ΔR1sin α-R+ΔR2tanα (6)

1.2.3 其他参数

圆心角αC式(7)计算,圆曲线长LC=C

αC=α-(β1+β2) (7)

2 对称基本型“两点法”线形及参数计算

2.1 线形设计与计算概述

图2为对称基本型“两点法”线形设计及主要参数示意图,点B为已设计的上一段平面线形终点,亦为当前平曲线设计段起点,其切线方向角αB已知。“两点法”线形设计时,仅需再拟定当前平曲线设计段的终点E即可获得“LS+LC+LS”对称基本型系列平曲线族;利用对称平曲线的几何特点,先计算获得α及总切线长TS,然后依据路线平面技术指标要求确定LSR的可行值范围,再以LSR为遍历参数,通过TSαLSRβ等参数之间的特定关系建立关于β的显示迭代公式,最后计算得到β值及其余参数。

图2  对称基本型线形及参数示意图

Fig. 2  Alignment and parameter diagram of symmetric basic horizontal curve

2.2 总切线长及参数关系

图2可知,在拟定当前平曲线设计段的终点E后,对称基本型“两点法”平曲线参数αTS即可通过式(8)、(9)计算确定,其中αBE为矢量BE的方向角(弧度),α0为由起点B切线方向旋转至矢量BE方向的偏转角。

α0=αBE-αBα=2α0 (8)
TS=0.5|BE|/cos(α/2) (9)

由对称性有LS=LS1=LS2β=β1=β2q=q1=q2ΔR=ΔR1=ΔR2,其中β=LS/(2R);由式(6)可得TS关于RLSα的关系式(10)、(11)。

TS=q+(R+ΔR)tan(α/2)=R f(α,β) (10)
f(α,β)=gx(β)+gy(β)tan(α/2) (11)

2.3 理论缓和曲线长度及圆曲线半径

2.3.1 缓和曲线理论最大长度及圆曲线理论最小半径

由缓和曲线基本特性可知,当LS增大时,β增加而R将减小;当LS达到理论最大值LSmax0时,相应的β亦达到理论最大值βmax0=0.5α,而R达到理论最小值Rmin0,由此可得LSmax0=2βmax0Rmin0=αRmin0,其中Rmin0式(12)计算;需注意应满足Rmin0Rmin,其中Rmin为路线平面技术指标要求的圆曲线最小半径。

Rmin0=TS/f(α,βmax0=0.5α) (12)

2.3.2 缓和曲线理论最小长度及圆曲线理论最大半径

依据路线平面线形技术要求,当R大于不设超高的圆曲线最小半径RminC时,可不设置缓和曲线,即LS理论最小长度LSmin0=0,易知此时R达到理论最大值Rmax0β达到理论最小值βmin0=0;将β=βmin0=0带入式(10),可获得Rmax0的计算式(13)。需注意应检验是否满足Rmax0RmaxC的条件,否则LSmin0应采用缓和曲线最小长度LSmin,即LSmin0=LSmin>0,此时Rmax0需依据式(10)计算确定。

Rmax0=TSf(α,βmin0=0)=TStan(α/2) (13)

2.4 线形参数求解过程

当已知αTSR时,由式(10)、(11)可得到式(14),该式仅含未知参数β。为便于求解β,将R=LS/(2β)带入式(14)得到式(15);此时若LS已知,则变量c为常数,式(15)为关于β的显示迭代计算公式。

f(α,β)=gx(β)+gy(β)tan(α/2)=TS /R=c0 (14)
β=c[gx(β)+gy(β)tan(α/2)]c=0.5LS/TS (15)

依据式(12)、(13)可获得LS的理论值范围[LSmin0LSmax0],再结合路线平面设计技术指标要求可确定LS的可行值范围[LSminLSmax];若以一定步长ΔLSLS进行遍历,则可基于常数c=0.5LS/TS,由式(15)迭代计算得到β。主要步骤简述如下:

(1)已知αTSLSminLSmax;设定LS的遍历步长ΔLS (如ΔLS=20 m),设定β计算精度εβ(如εβ=1.0×10-7);开始针对LS进行循环遍历。

(2)由当前LS计算常数c= 0.5LS/TS;设定β初始值β0=0.1α,令βp=β0

(3)将β=βp带入迭代公式(15),计算获得新的偏转角βq=c f αβp)。

(4)计算偏转角差值Δβ=|βq-βp|,若Δβ>εβ,则令βpq,返回步骤(3),进行下一次迭代计算。

(5)若Δβεβ,迭代结束,获得满足精度要求的偏转角β=βq

获得偏转角β后,由当前LS值可得R=0.5LS/β、圆心角αC=α-2β及圆曲线长LC=C等其余参数并确定出“LS+LC+LS”单条对称型平曲线,LS遍历结束即可获得系列平曲线族。

3 非对称基本型平曲线参数及计算

3.1 线形设计与计算概述

图3为非对称基本型“两点法”线形设计及主要参数示意图,与对称基本型“两点法”线形设计有以下不同:①LS1LS2,因此β1β2q1q2ΔR1ΔR2TS1TS2;②除需拟定终点E之外,还需指定终点E的切线方向αE。先计算平曲线偏转角αTS1TS2,并获得圆曲线半径R的理论值范围[Rmin0 Rmax0],再结合路线平面技术指标要求确定R的可行值区间[RminRmax],然后针对R进行遍历,采用“双迭代法”计算获得β1β2,最后确定LS1LS2LCαC等其余参数。

图3  非对称基本型线形及参数示意图

Fig. 3  Diagram of alignment and parameters of asymmetric basic horizontal curve

3.2 总切线长计算分析

依据式(2)、(3)、(6)可获得TS1TS2αRβ1β2的关系式(16),其中Rβ1β2应同时满足如式(17)所示的几何限定条件及路线平面技术指标要求;由式(16)可得函数fTS=TS1/TS2αβ1β2的关系式(18)

TS1R=gx(β1)+gy(β2)sinα-gy(β1)tanαTS2R=gx(β2)+gy(β1)sinα-gy(β2)tanα (16)
0<RminRRmax0<βminβ1βmax<α0<βminβ2βmax<αβ1+β2α (17)
fTS=TS1TS2=gx(β1)+gy(β2)sinα-gy(β1)tanαgx(β2)+gy(β1)sinα-gy(β2)tanα (18)

为探寻函数fTS的变化规律,在α∈[1°,179°]范围内以及β1∈[0,1.4]、β2∈[0,1.4]条件下,以步长 Δα=1°、Δβ=0.01 rad,由式(18)计算[β1 β2]域内的fTS,并绘制αfTS关系如图4所示。由图4可知,对任意偏转角α,比值fTS=TS1/TS2限定在一定范围[fTS minfTS max],且fTS minfTS max互为倒数关系;当α较小、接近于0°时,fTS min=0.5、fTSmax=2.0;当α较大、接近于180°时,fTS min=fTS max=1.0。如果总切线长TS1TS1相差过大,当比值fTS小于0.5或大于2.0时,将不存在满足要求的非对称型平曲线。可获得平曲线偏转角αfTS max的二次多项式函数拟合关系如式(19)所示,拟合相关系数R2=0.999 425 4。

fTSmax(α)=2+0.022 774 6α-0.104 851 2α2 (19)

图4  αfTS关系图

Fig. 4  Diagram of correlation between α and fTS

3.3 圆曲线半径理论范围

3.3.1 理论圆曲线最小半径

对于完整型缓和曲线,起点处LS =0、β=0、R=∞;随LS增加,β亦增大,与缓和曲线终点相接的圆曲线半径R及圆曲线长LC均将减小。由此可知,当第1、第2缓和曲线偏转角之和达到最大值即β1+ β2=α时,相应的RLC将分别达到理论最小值Rmin0LCmin0,其中LCmin0=0。引入上述条件,由式(16)整理得到式(20),可采用“夹逼法”求解获得Rmin0及相应的β1β2

R1=TS1sinαgx(β1)sinα-gy(β1)cosα+gy(β2)R2=TS2sinαgx(β2)sinα-gy(β2)cosα+gy(β1)β1+β2=αR1=R2=Rmin0 (20)

3.3.2 理论圆曲线最大半径

由前述分析可知,当TS1>TS2时,若设有缓和曲线LS1LS2,则必有LS1> LS2,故当圆曲线半径R达到极大值Rmax0时,必有LS2=0、β2=0及gxβ2)=0、 gyβ2)=1,带入式(16)整理可得到式(21);若TS1< TS2,则定有LS1=0、β1=0及gxβ1)=0、gyβ1)=1,带入式(16)整理可获得与式(21)类似的关系式。式(21)亦可采用“夹逼法”求解获得Rmax0β1β2,与前述Rmin0求解方法类似。

R1=TS1sinαgx(β1)sinα-gy(β1)cosα+1R2=TS2sinαgy(β1)-cosα=R1=Rmax0 (21)

3.4 参数求解方法

3.4.1 概述

在已知αTS1TS2条件下,式(16)含2个方程、3个互相独立的未知数Rβ1β2,故不存在唯一解。由上文分析可知,依据几何限定条件及式(20)~ (21)可计算获得圆曲线半径的理论最小值Rmin0、理论最大值Rmax0,再考虑路线平面技术指标要求的圆曲线最小半径Rmin1、最大半径Rmax1,则可最终确定R的可行值范围[RminRmax],并同时获得缓和曲线长度及偏转角范围[LSminLSmax]、[βminβmax];由此针对R进行遍历循环,即可在已知αTS1TS2R条件下求解式(16),获得β1β2,继而得到其余参数。

3.4.2 缓和曲线偏转角求解的双迭代法

求解缓和曲线偏转角β1β2的计算过程均基于针对R的遍历循环,参数αTS1TS2R均已知。

(1)设定β1β2的计算精度εΔβ(如εΔβ=1.0×10-7 rad),设置β1β2迭代计算初始值β1pβ2p (如β1p=0.25αβ1p=βmin)。

(2)由式(16)的第1个方程,带入β1p计算获得常数c1,整理后得式(22)、(23)。

c1=sinα[TS1/R-gx(β1p)]+                         cosα gy(β1p)a2+b2 (22)
gy(β2)=c1=2β22fy(β2)+cosβ2 (23)

(3)依据式(23)得到β2的显示迭代计算式(24),迭代求解获得满足精度要求的β2,记β2=β2q

β2=c1-cos β22β2fy(β2) (24)

(4)由式(16)的第2个方程,带入β2=β2q计算获得常数c2并整理得式(25)、(26)。

c2=sinαTS2R-gx(β2q)+gy(β2q)tanα (25)
gy(β1)=c2=2β12fy(β1)+cosβ1 (26)

(5)依据式(26)得到β1的显示迭代计算式(27),迭代求解获得满足精度要求的β1新值,记β1=β1q

β1=c2-cos β12β1fy(β1) (27)

(6)计算偏转角差值Δβ1=|β1q-β1p|;若Δβ1εΔβ,则获得满足精度要求的β1=β1qβ2=β2q,计算结束;否则令β1p=β1q,返回步骤(2)进入下一轮迭代计算,直至满足偏转角差值计算精度要求。

在上述每轮计算过程中,需对β1β2进行两次迭代求解,故称为“双迭代法”。

3.4.3 缓和曲线偏转角求解的改进双迭代法

前述“双迭代法”的迭代式(24)、(27)均为偏转角β2β1的近似二次式,计算表明在求解过程中对β2β1的迭代初始值限定较为苛刻,每个迭代计算轮回中β1收敛性尚可,但β2收敛性欠佳。采用以下改进方法,可明显提高迭代求解的收敛性。

(1) 与“双迭代法”相同,设定β1β2的计算精度

εΔβ及初始值。

(2) 由式(16)的第1个方程,带入β1p计算获得常数c1,有gyβ2)= c1,参见式(22)、(23)。

(3) 由式(16)的第2个方程,带入β1p并替换

gyβ2)= c1,计算获得常数c3,整理得式(28)、(29)。

c3=TS2R-gy(β1p)sin α+c1tan α (28)
gx(β2)=c3=2β2fx(β2)-sin β2 (29)

(4) 依据式(29)得到β2的显示迭代计算式(30),迭代求解获得满足精度要求的β2,记β2=β2q

β2=c2+sin β22fx(β2) (30)

(5) 由式(16)的第2个方程,带入β2=β2q计算获得常数c4并整理得式(31)、(32)。

c4=(sin α )TS2R-gx(β2q)+gy(β2q)tan α (31)
gy(β1)=c4 (32)

(6)由式(16)的第1个方程,带入β2=β2q并替换gyβ1)= c4,计算获得常数c5,整理得式(33)、(34)。

c5=TS1R-gy(β2q)sin α+c4tan α (33)
gx(β1)=c5=2β1fx(β1)-sin β1 (34)

(7)由式(34)得到β1的显示迭代计算式(35),迭代求解获得满足精度要求的β1新值,记β1=β1q

β1=c5+sin β12fx(β1) (35)

(8)计算偏转角差值Δβ1=|β1q-β1p|;若Δβ1εΔβ,则获得满足精度要求的β1=β1qβ2=β2q,计算结束;否则令β1p=β1q,返回步骤(2)进入下一轮迭代计算,直至满足偏转角差值计算精度要求。

上述“改进双迭代法”的迭代计算式(30)、(35)为偏转角β2β1的近似一次式,计算分析表明β2β1可取的初始值范围较大,迭代算法的收敛性、可靠性得以提高。另外还可采用“夹逼+迭代”的综合方法,限于篇幅不再赘述。

4 设计计算示例

4.1 基本资料

图5为对称基本型平曲线“两点法”计算示意图,其中,起点B(100.0,500.0)、终点E(844.0,685.5)、起点方向角αB=30°。计算得TS=398.838 6 m、α0=16°、α=2α0=32°(0.558 505 361 rad)。

图5  对称基本型平曲线“两点法”设计计算示例

Fig. 5  Example of the “two-point method” design and calculation for symmetric basic horizontal curve

图6为非对称基本型平曲线“两点法”计算示意图,其中,起点B(100.0,500.0)、终点E(844.0,685.5)、起点方向角αB=30°,拟定终点E的切线方向角αE=-5°。计算得α0=16°、α=35°(0.610 865 238 rad)、 TS1=435.230 6 m、TS2=368.481 5 m;可知fTS =TS1/TS2=1.181 1,带入式(19)fTSmax=1.880 6,显然fTS< fTSmax,故式(16)有解;注意到TS1> TS2,则在所有平曲线族中,将有LS1> LS2

图6  非对称基本型平曲线“两点法”设计计算示例

Fig.6  Example of the “two-point method” design and calculation for asymmetric basic horizontal curve

4.2 对称基本型计算结果

4.2.1 计算结果

由式(12)、(13)分别得到圆曲线的理论最小半径Rmin 0=700.959 m、理论最大半径Rmax0= 1 390.915 m,缓和曲线理论最大偏转角βmax0= 0.5α=0.279 252 680 rad,缓和曲线理论最大长度LS max0=2βmax0Rmin0=391.489 m。取缓和曲线长LS为20 m倍数,针对LS进行遍历循环,由式(15)迭代计算获得缓和曲线偏转角β,继而得到其余参数。平曲线族及主要参数如图7所示,LS遍历及相应参数计算结果如表1所示,其中LHA为平曲线总长的平均值,ND为缓和曲线偏转角迭代计算次数。

图7  对称基本型平曲线“两点法”设计主要参数及平曲线族

Fig. 7  Main parameters and horizontal curve family of the “two-point method” design for symmetric basic horizontal curve

表1  对称基本型平曲线“两点法”主要参数计算结果
Tab. 1  Results of main parameters of the “two-point method” for symmetrical basic horizontal curve
LSβRLCLS/LCLHLH/LHAND
0 0 1390.915 776.834 0.00 776.834 0.995 0
20.000 0.007 374 5 1356.029 737.350 0.03 777.350 0.995 6 6
40.000 0.015 138 7 1321.117 697.851 0.06 777.851 0.996 3 6
360.000 0.237 704 0 757.244 62.925 5.72 782.925 1.002 8 21
380.000 0.263 334 1 721.517 22.971 16.54 782.971 1.002 8 22
391.489 0.279 252 5 700.959 0 782.978 1.002 8 29

为验证式(15)迭代计算的收敛性,以LS=200 m、c=0.5LS/TS=0.250 727 99为例,设置β初始值为βp=0.1、4.0 rad两种极端条件,迭代计算过程如表2表3所示,由表可知迭代计算收敛性很好。

表2  对称基本型平曲线“两点法”LS=200 m计算过程(βp=0.1)
Tab. 2  Calculation process of the “two-point method” for symmetrical basic horizontal curve at LS=200 m ( βp=0.1 )
NDβpfxβfyββqΔβ=|βqp |
1 0.100 000 000 0.999 000 463 0.333 0953 14 0.097 079 331 2.920 67×10-3
2 0.097 079 331 0.999 057 971 0.333 109 010 0.096 340 856 7.384 75×10-4
3 0.096 340 856 0.999 072243 0.333 112 409 0.096 154 162 1.866 94×10-4
7 0.096 092 019 0.999 077 027 0.333 113548 0.096 091 256 7.624 83×10-7
8 0.096 091 256 0.999 077 042 0.333 113 552 0.096 091 064 1.927 54×10-7
9 0.096 091 064 0.999 077 045 0.333 1135 53 0.096 091 015 4.872 78×10-8
表3  对称基本型平曲线“两点法”LS =200 m计算过程(βp =4.0)
Tab. 3  Calculation process of the “two-point method” for symmetrical basic horizontal curve at LS =200 m (βp=4.0)
NDβpfxβfyββqΔβ=|βqp |
1 4.000 000 000 0.230 830 968 0.110 482 307 0.859 944 449 3.140 06×100
2 0.85 9944 449 0.928 538 573 0.316 135 071 0.290 930 962 5.690 13×10-1
3 0.290 9309 62 0.991 569 020 0.331 323 495 0.145 645 470 1.452 85×10-1
12 0.096 091 827 0.999 077 031 0.333 113 549 0.096 091 208 6.190 96×10-7
13 0.096 091208 0.999 077 043 0.333 1135 52 0.096 091 051 1.565 06×10-7
14 0.096 091 051 0.999 077 046 0.333 113 553 0.096 091 012 3.956 44×10-8

4.2.2 线形参数特点

(1)存在规律性分布的含3段线元的对称型平曲线族,各曲线不重叠、不交叉。

(2)圆曲线半径R、缓和曲线长LS均存在理论最小值及最大值;受αTS几何条件限制,LSR不能同时取整。

(3)平曲线族仅分布于中间宽、两边窄的特定区域范围,其中圆曲线族位于内侧、缓和曲线族位于外侧,圆曲线族占据绝大部分区域。

(4)LS增加时,平曲线参数均呈规律性变化,其中RLC呈减少趋势,平曲线总长LH呈缓慢增加趋势,迭代计算次数呈增加趋势;LH相对变化值在0.5%以内,即其值近似恒定。

4.3 非对称基本型计算结果

4.3.1 圆曲线半径理论值

依据式(20)~(21),采用“夹逼法”确定圆曲线理论最小半径Rmin0 (见表4)、理论最大半径Rmax0β1初始边界为β1min=0、β1max=α=0.610 865 238,每次计算均依据R1R2的大小关系(见表4 “GL”列)对β1minβ1max进行替换,逐渐缩减“夹逼”范围;本示例分别经过23次、26次逼近计算,获得Rmin0=643.178 m、Rmax0=1164.917 m, 计算精度ΔR<1.0×10-5

表4  非对称基本型平曲线“两点法”圆曲线理论最小半径Rmin0计算过程
Tab. 4  Calculation process of the theoretical minimum radius Rmin0 of the “two-point method”" circular curve of asymmetric basic horizontal curve
NDβ1minβ1maxβ1β2R1GLR2ΔR=| R1 -R2|β1替换
1 0 0.610 865 238 0.305 432 619 0.305 432 619 696.7366 4 > 589.881 71 1.068 549×10-2 β1min
2 0.305 432 619 0.610 8652 38 0.458 148 929 0.152 716 310 598.301 58 < 704.706 65 1.064 051×102 β1max
3 0.305 432 619 0.458 148 929 0.381 790 774 0.229 074 464 643.883 63 > 642.354 75 1.528 881×100 β1min
21 0.382 892 407 0.382 892 990 0.382 892 698 0.227 972 540 643.177 94 < 643.178 23 2.963 525×10-4 β1max
22 0.382 892 407 0.382 892 698 0.382 892 553 0.227 972 685 643.178 03 < 643.178 12 9.420 483×10-5 β1max
23 0.382 892 407 0.382 892 553 0.382 892 480 0.227 972 758 643.1780 8 > 643.178 07 6.868 999×10-6 (结束)

4.3.2 其他参数计算结果

采用“改进双迭代法”的计算结果如表5所示,表中LHA为平曲线总长的平均值。以R=900 m为例,表6列出了“改进双迭代法”的部分计算过程,其中ND为“双迭代”轮回数,ND1ND2分别为每个迭代轮回的两次迭代次数。

表5  非对称基本型平曲线“两点法”主要参数计算结果(改进双迭代法)
Tab. 5  Calculation results of main parameters of the “two-point method” for asymmetric basic horizontal curve (improved double iteration method)
Rβ1β2αCLS1LS2LCLHLH/LHA
Rmin0 0.382 892 4 0.227 972 8 0 492.54 293.25 0 785.790 1.002 4
660 0.363 467 7 0.216 248 6 0.031 148 9 479.78 285.45 20.558 785.784 1.002 4
680 0.341 834 2 0.202 829 0 0.066 202 0 464.89 275.85 45.017 785.759 1.002 4
900 0.173 403 1 0.087 230 4 0.350 2318 312.13 157.01 315.209 784.349 1.000 6
1 140 0.067 566 4 0.006 614 7 0.536 684 1 154.05 15.08 611.820 780.953 0.996 3
1 160 0.0607 86 3 0.001 2845 0.548 794 4 141.02 2.98 636.602 780.602 0.995 8
Rmax 0.059 155 8 0 0.551 709 4 137.82 0 642.696 780.519 0.995 7
表6  非对称基本型平曲线“两点法”R=900 m计算过程(改进双迭代法)
Tab. 6  Calculation process of the “two-point method” for asymmetric basic horizontal curve at R=900 m (improved double iteration method)
NDβ1p、c1、c3ND1β2pβ2qΔβ2β2q、c4、c5ND2β1mβ1qΔβ1q
1

(令β1p=0.25α

β1p=0.152 716 3

c1=1.012 182 6

c3=0.104 752 4

1 0.152 716 3 0.128 737 8 0.023 978 5

β2q=0.104 790 9

c4=0.995 403 0

c5= 0.158 536 0

1 0.152 716 3 0.155 692 4 0.002 976 1
2 0.128 7378 0.116 760 8 0.011 9770 2 0.155 692 4 0.157 180 6 0.001 488 2
3 0.116 760 8 0.110 775 0 0.005 985 8 3 0.157 180 6 0.157 924 8 0.000 744 2
19 0.104 790 9 0.104 790 8 0.000 000 1 16 0.158 668 8 0.158 668 9 0.000 000 1
第1次双迭代计算结果:β1p= 0.152 716 31、β1q=0.158 668 95,Δβ1=|Δβ1q-Δβ1p|=5.953×10-3,需继续迭代计算
2

β1p=β1q

β1p=0.158 668 9

c1=1.009 029 2

c3=0.099 711 1

1 0.158 668 9 0.129 182 4 0.029 486 4

β2q=0.099 744 3

c4=0.998 153 7

c5= 0.162 7641

1 0.158 668 9 0.160 788 4 0.002 119 4
2 0.129 182 4 0.114 458 2 0.014 724 2 2 0.160 788 4 0.161 848 2 0.001 059 8
3 0.114 458 2 0.107 100 0 0.007 358 1 3 0.161 848 2 0.162 378 1 0.000 529 9
19 0.099 744 4 0.099 744 3 0.000 000 1 16 0.162 907 8 0.162 907 9 0.000 000 1
第2次双迭代计算结果:β1p=0.158 668 90、β1q=0.162 908 09,Δβ1=|Δβ1q-Δβ1p|=4.239×10-3,需继续迭代计算
33

β1p=β1q

β1p=0.173 403 0

c1=1.001 268 0

c3=0.087 208 2

1 0.173 403 0 0.130 262 9 0.043 140 1

β2q=0.087 230 4

c4=1.005 005 9

c5=0.173 229 4

1 0.173 403 0 0.173 403 0 0.000 000 1
2 0.130 262 9 0.108 735 9 0.021 527 0 2 0.173 403 0 0.173 403 1 0
3 0.108 735 9 0.097 980 8 0.010 755 1
19 0.087 230 5 0.087 230 4 0.000 000 1
第33次双迭代计算结果:β1p=0.173 402 97、β1q=0.173 403 07,Δβ1=|Δβ1q-Δβ1p|=1.002×10-7,需继续迭代计算
34

β1p=β1q

β1p=0.173 403 1

c1=1.001 268 0

c3=0.087 208 1

1 0.173 403 1 0.130 262 9 0.043 140 2

β2q=0.087 230 4

c4=1.005 005 9

c5=0.173 229 5

1 0.173 403 1 0.173 403 1 0
2 0.130 262 9 0.108 735 9 0.021 527 1
3 0.108 735 9 0.097 980 7 0.010 755 1
19 0.087 230 4 0.087 230 4 0
第34次双迭代计算结果:β1p=0.173 403 07、β1q=0.173 403 14,Δβ1=|Δβ1q-Δβ1p|=7.000×10-8,满足精度要求,计算结束

注:  限于篇幅,每个“双迭代”轮回中未列出第4步至倒数第2步的中间计算结果。

平曲线族及主要参数如图8所示。

图8  非对称基本型平曲线“两点法”设计平曲线族及主要参数

Fig. 8  Main parameters and horizontal curve family of the “two-point method” design for asymmetric basic horizontal curve

4.3.3 线形参数特点

(1) 存在规律性分布的含3段线元的非对称型平曲线族,各曲线不重叠、不交叉。

(2) 仅当比值fTS=TS1/TS2在(0.5,2.0)范围内才可设置基本型平曲线,且圆曲线半径RLS1LS2均存在理论最小值与理论最大值;受TS1TS2α限制,LS1LS2R不能同时取整。

(3) 与对称型平曲线相似,非对称型平曲线族亦仅分布在中间宽、两边窄的特定区域范围内,圆曲线族位于内侧、缓和曲线族位于外侧,圆曲线族占据绝大部分区域;但平曲线族分布形状不对称,中间较宽区域偏向于总切线长较短的一侧。

(4) 当R增加时,平曲线参数均呈现规律性变化,其中LC呈增加趋势,LS1LS2LH呈减少趋势,LH相对变化在0.5%以内,即其近似为恒定值。

(5) LS1LS2呈现同步变化趋势,即:若LS1 (或LS2)增大,则LS2(或LS1)亦增加。

(6) 总切线长相对较大的一侧,其缓和曲线长亦较长;反之,总切线长相对较小的一侧,其缓和曲线长相对较短。

5 工程应用对比

图9a为某二级公路K6+959.85~K7+744.434段采用导线法的平面线形设计图,其中JD13为对称基本型、JD14为非对称基本型,两者构成S型平曲线,点ABC分别为平曲线ZH点(直缓点)和HZ点(缓直点),点B为S型曲线拐点。作为对比,图9b为采用“两点法”进行相应平面线形布设的过程示意图。

图9  某二级公路平面线形设计对比分析

Fig. 9  Comparison and analysis of horizontal alignment drawing of second-class highway

由图可知:①由ZH(直缓)、HZ(缓直)两点可构建“缓和曲线+圆曲线+缓和曲线”3段线元组成的系列平曲线族;②平曲线族各曲线参数及几何分布规律,构成两端狭窄、中间宽厚的“特定区域”,线形布设灵活多变、便于比选优选。

6 结论

(1)仅需起点(直缓点)与终点(缓直点)两点即可确定由“第1缓和曲线+圆曲线+第2缓和曲线”3段线元组成的基本型平曲线;在顺次进行平面线形布设时,仅需拟定当前设计段的终点即可。

(2)仅当平曲线总切线长的比值在(0.5,2.0)范围内时,圆曲线半径、缓和曲线长度才有理论解;对称基本型的总切线长比值始终为1,故总有理论解。

(3)线形设计成果为3段线元组成的系列“平曲线族”而非单条曲线,各曲线互不交叉重叠,且分布在两侧狭窄、中间宽厚的“特定区域”内。

(4)圆曲线半径、缓和曲线长度、圆曲线长度及平曲线总长,均存在理论最小值、最大值;平曲线总长的变化值在0.5%以内,其平均值近似为固定值。

(5)当圆曲线半径增加时,圆曲线长呈增加趋势,而缓和曲线长呈减少趋势,平曲线总长亦呈减少趋势但变化很小;第1、第2缓和曲线长度具有同步增减趋势,总切线长较大一侧对应的缓和曲线长亦较长。

(6)“两点法”突显坐标点位的控制作用,呈规律性分布的“平曲线族”和“特定区域”有益于路线优选、线形拟合重构和智能化路线设计。

作者贡献声明

李玉华:提出两点法设计基本思想、计算原理,设计计算软件编制,论文撰写。

王奥凯:工程应用示例,软件调试。

刘佳音:计算原理验证,设计方法检验。

马洋洋:软件调试,文字校对。

吴树铭:软件调试。

参考文献

1

许金良程建川栗志海. 道路勘测设计[M].5版. 北京人民交通出版社2018. [百度学术] 

XU JinliangCHENG JianchuanLI Zhihai. Road survey and design[M] 5th ed. BeijingChina Communications Press2018. [百度学术] 

2

吴国雄公路平面线形曲线型设计方法综述[J].公路19969): 21. [百度学术] 

WU Guoxiong. A summary of curve design method for highway horizontal alignment[J]. Highway19969): 21. [百度学术] 

3

宋文. “附合导线法”设计匝道平面线形[J]. 公路19919):11 [百度学术] 

SONG Wen. Design of horizontal alignment for interchange ramps by “connecting traverse method”[J]. Highway19919):11. [百度学术] 

4

丁建明李方公路平面线形设计的五单元导线法[J].东南大学学报1998282): 155. [百度学术] 

DING JianmingLI Fang. Five-element traverse method for highway horizontal alignment design[J]. Journal of Southeast University1998282): 155. [百度学术] 

5

吴国雄. 一种新的公路平面线形曲线型设计方法—综合法[J].重庆交通学院学报199615增刊):59. [百度学术] 

WU Guoxiong. A new curve method for designing highway horizontal alignment - comprehensive method[J]. Journal of Chongqing Jiaotong Institute199615S):59. [百度学术] 

6

许金良杨少伟徐其福互通式立交平面线形控制线元法[J].西安公路交通大学学报1997174):5. [百度学术] 

XU JinliangYANG ShaoweiXU Qifu. The controlled segmentary element design method of interchange[J]. Journal of Xi’an Highway University1997174): 5. [百度学术] 

7

李玉华孙依人刘佳音. 路线平面线形设计计算的两点线元法原理[J]. 同济大学学报(自然科学版)2021494):507. [百度学术] 

LI YuhuaSUN YirenLIU Jiayinet al. Principle of two-point segmentary element method for horizontal alignment design and calculation[J]. Journal of Tongji University (Natural Science)2021494): 507. [百度学术] 

8

LI WPU H. A method for automatically recreating the horizontal alignment geometry of existing railways[J]. Computer-Aided Civil and Infrastructure Engineering2019341):71. [百度学术] 

9

CAMACHO-TORREGROSA F JPEREZ-ZURIAGA A MCAMPOY-UNGRIA J Met al. Use of heading direction for recreating the horizontal alignment of an existing road[J]. Computer-Aided Civil and Infrastructure Engineering2015304): 282. [百度学术] 

10

缪鹍詹振炎基于直线约束的道路线形设计通用方法[J].中国公路学报2002153): 15. [百度学术] 

MIU KunZHAN Zhenyan. Universal method for highway’s horizontal alignment based on restrictions of straight line[J]. China Journal of Highway and Transport2002153): 15. [百度学术] 

11

吴国雄陈欣斗. 道路立交平面线形圆心连线闭合导线设计法[J]. 重庆交通学院学报2003224):20. [百度学术] 

WU GuoxiongCHEN Xindou. Closed-traverse method of center-center lines of circles in horizontal alignment design of road interchange[J]. Journal of Chongqing Jiaotong University2003224): 20. [百度学术] 

12

张航韦金君张肖磊道路平面线形拟合方法比较研究[J].武汉理工大学学报(交通科学与工程版)2018424):594. [百度学术] 

ZHANG HangWEI JinjunZHANG Xiaolei. Comparative study on road plane alignment fitting methods[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering)2018424): 594. [百度学术] 

13

BOSURGI GD’ANDREA A. A polynomial parametric curve (PPC-CURVE) for the design of horizontal geometry of highways[J]. Computer-Aided Civil and Infrastructure Engineering2012274): 303. [百度学术] 

14

徐进罗庆毛嘉川. 轨迹‒速度协同控制的山区公路平面线形设计方法[J]. 中国公路学报2013264):43. [百度学术] 

XU JinLUO QingMAO Jiachuanet al. Method for horizontal geometry design of mountainous roads based on trajectory-speed cooperative control[J]. China Journal of Highway and Transport2013264): 43. [百度学术] 

15

于玲秦翔包龙生. 基于遗传算法的山区高速公路平面线形优化设计方法[J]. 公路交通科技20181659): 175. [百度学术] 

YU LingQIN XiangBAO Longsheng. Optimization design method for horizontal alignment of mountain expressway based on genetic algorithm[J]. Journal of Highway and Transportation Research and Development20181659): 175. [百度学术] 

16

韩元利靖仕元陈燕平. 铁路智能选线过程中平曲线最优配置模型研究[J]. 铁道工程学报20212758):1. [百度学术] 

HAN YuanliJING ShiyuanCHEN Yanping. Rsearch on the optimal parametric configuration of horizontal curve on railway intelligent design[J]. Journal of Railway Engineering Society20212758): 1. [百度学术] 

17

李玉华赵延庆周长俊. 一种道路路线平面线形设计的“两点”法ZL201810255466.4 [P].2022-04-08. [百度学术] 

LI YuhuaZHAO YanqingZHOU Changjunet al. A“two-point method” for designing horizontal alignment of road routesZL201810255466.4 [P]. 2022-04-08. [百度学术] 

18

刘也嘉. 道路平面线形设计计算的两点法[D]. 大连大连理工大学2018. [百度学术] 

LIU Yejia. Two-point method for design and calculation of road horizontal alignment[D]. DalianDalian University of Technology2018. [百度学术]