摘要
为改善传统导线法存在导线与平曲线偏离较大等不足,提出“两点法”基本型平曲线形设计方法。仅需利用导线法中的直缓点与缓直点两个端点即可获得由缓和曲线及圆曲线共3段线元组成的基本型平曲线形,其平曲线总偏转角α、总切线长T为已知值。先采用迭代法确定缓和曲线长度LS (对称基本型)或采用夹逼法获得圆曲线半径R(非对称基本型)的理论最小值、最大值,并结合平面线形技术要求确定LS或R的设计取值范围;然后基于LS或R进行遍历循环,通过迭代法获得缓和曲线偏转角β,进而获得各线元其余参数。研究结果表明,“两点法”设计可获得系列“平曲线族”,线形分布规律,互不交叉重叠,且全部曲线位于两端细窄、中间宽厚的“特定设计区域”;LS、R、圆曲线长LC及平曲线总长LH等参数呈规律性变化,其中LS与R、LC的增减趋势相反,LH平均值近似固定,非对称基本型的总切线长度之比须在(0.5,2.0)范围内方可有解。“两点法”突显坐标位置点的控制作用,线形布设高效,设计成果为特定区域内的系列平曲线族而非单条曲线,有益于结合实际地形及设计意图进行线形比选、优选和路线智能化设计。
近几十年来我国公路及铁路建设迅猛发展,已形成较为完善的路线平面线形设计方法,其中导线
在复杂地形及立交匝道线形设计实践中,路线设计方法进一步发展,如针对导线法改进的附合导线
导线法与曲线法两种方法相结合是较好的创新思路,如曲直
综上,目前平面线形设计方法面临的主要问题有:①两种基本方法均存在固有不足;②设计成果通常为单一平曲线,不利于进行多种适宜线形的比选。两点线元法对此有所改进,但仍局限于曲线法优势。
近年来,随着高新技术及人工智能的发展,线形优
如

图1 导线法基本型平曲线组成及参数
Fig. 1 Composition and parameters for basic horizontal curves using traverse method

图2 对称基本型线形及参数示意图
Fig. 2 Alignment and parameter diagram of symmetric basic horizontal curve
由
(8) |
(9) |
由对称性有LS==、β=β1=β2、q=q1=q2及R==,其中β=LS/(2R);由
(10) |
(11) |
由缓和曲线基本特性可知,当LS增大时,β增加而R将减小;当LS达到理论最大值LSmax0时,相应的β亦达到理论最大值βmax0=0.5α,而R达到理论最小值Rmin0,由此可得LSmax0=2βmax0Rmin0=αRmin0,其中Rmin0由
(12) |
当已知α、TS及R时,由式(
(14) |
(15) |
依据式(
(1)已知α、TS、LSmin及LSmax;设定LS的遍历步长Δ (如Δ=20 m),设定β计算精度εβ(如εβ=1.0×1
(2)由当前LS计算常数c= 0.5LS/TS;设定β初始值β0=0.1α,令βp=β0。
(3)将β=βp带入迭代
(4)计算偏转角差值Δβ=|βq-βp|,若Δβ>εβ,则令βp=βq,返回步骤(3),进行下一次迭代计算。
(5)若Δβ≤εβ,迭代结束,获得满足精度要求的偏转角β=βq。
获得偏转角β后,由当前LS值可得R=0.5LS/β、圆心角αC=α-2β及圆曲线长LC=RαC等其余参数并确定出“LS+LC+LS”单条对称型平曲线,LS遍历结束即可获得系列平曲线族。

图3 非对称基本型线形及参数示意图
Fig. 3 Diagram of alignment and parameters of asymmetric basic horizontal curve
依据式(
(16) |
(17) |
(18) |
为探寻函数f的变化规律,在α∈[1°,179°]范围内以及β1∈[0,1.4]、β2∈[0,1.4]条件下,以步长 Δα=1°、Δβ=0.01 rad,由
(19) |

图4 α、f关系图
Fig. 4 Diagram of correlation between α and f
在已知α、T、T条件下,
求解缓和曲线偏转角β1、β2的计算过程均基于针对R的遍历循环,参数α、T、、R均已知。
(1)设定β1、β2的计算精度ε(如ε=1.0×1
(2)由
(22) |
(23) |
(3)依据
(24) |
(4)由
(25) |
(26) |
(5)依据
(27) |
(6)计算偏转角差值Δ=|β1q-β1p|;若Δ≤ε,则获得满足精度要求的β1=β1q及β2=β2q,计算结束;否则令β1p=β1q,返回步骤(2)进入下一轮迭代计算,直至满足偏转角差值计算精度要求。
在上述每轮计算过程中,需对β1、β2进行两次迭代求解,故称为“双迭代法”。
前述“双迭代法”的迭代式(
(1) 与“双迭代法”相同,设定β1、β2的计算精度
ε及初始值。
(2) 由
(3) 由
gy(β2)= c1,计算获得常数c3,整理得式(
(28) |
(29) |
(4) 依据
(30) |
(5) 由
(31) |
(32) |
(6)由
(33) |
(34) |
(7)由
(35) |
(8)计算偏转角差值Δ=|β1q-β1p|;若Δ≤ε,则获得满足精度要求的β1=β1q及β2=β2q,计算结束;否则令β1p=β1q,返回步骤(2)进入下一轮迭代计算,直至满足偏转角差值计算精度要求。
上述“改进双迭代法”的迭代计算式(

图5 对称基本型平曲线“两点法”设计计算示例
Fig. 5 Example of the “two-point method” design and calculation for symmetric basic horizontal curve

图6 非对称基本型平曲线“两点法”设计计算示例
Fig.6 Example of the “two-point method” design and calculation for asymmetric basic horizontal curve
由式(

图7 对称基本型平曲线“两点法”设计主要参数及平曲线族
Fig. 7 Main parameters and horizontal curve family of the “two-point method” design for symmetric basic horizontal curve
LS | β | R | LC | LS/LC | LH | LH/LHA | ND |
---|---|---|---|---|---|---|---|
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 |
为验证
ND | βp | fx(β) | fy(β) | βq | Δβ=|βq-βp | |
---|---|---|---|---|---|
1 | 0.100 000 000 | 0.999 000 463 | 0.333 0953 14 | 0.097 079 331 |
2.920 67×1 |
2 | 0.097 079 331 | 0.999 057 971 | 0.333 109 010 | 0.096 340 856 |
7.384 75×1 |
3 | 0.096 340 856 | 0.999 072243 | 0.333 112 409 | 0.096 154 162 |
1.866 94×1 |
7 | 0.096 092 019 | 0.999 077 027 | 0.333 113548 | 0.096 091 256 |
7.624 83×1 |
8 | 0.096 091 256 | 0.999 077 042 | 0.333 113 552 | 0.096 091 064 |
1.927 54×1 |
9 | 0.096 091 064 | 0.999 077 045 | 0.333 1135 53 | 0.096 091 015 |
4.872 78×1 |
ND | βp | fx(β) | fy(β) | βq | Δβ=|βq-βp | |
---|---|---|---|---|---|
1 | 4.000 000 000 | 0.230 830 968 | 0.110 482 307 | 0.859 944 449 |
3.140 06×1 |
2 | 0.85 9944 449 | 0.928 538 573 | 0.316 135 071 | 0.290 930 962 |
5.690 13×1 |
3 | 0.290 9309 62 | 0.991 569 020 | 0.331 323 495 | 0.145 645 470 |
1.452 85×1 |
12 | 0.096 091 827 | 0.999 077 031 | 0.333 113 549 | 0.096 091 208 |
6.190 96×1 |
13 | 0.096 091208 | 0.999 077 043 | 0.333 1135 52 | 0.096 091 051 |
1.565 06×1 |
14 | 0.096 091 051 | 0.999 077 046 | 0.333 113 553 | 0.096 091 012 |
3.956 44×1 |
依据式(
ND | β1min | β1max | β1 | β2 | R1 | GL | R2 | Δ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×1 | β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×1 | β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×1 | β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×1 | β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×1 | β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×1 | (结束) |
采用“改进双迭代法”的计算结果如
R | β1 | β2 | αC | LS1 | LS2 | LC | LH | LH/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 |
ND | β1p、c1、c3 | ND1 | β2p | β2q | Δ | β2q、c4、c5 | ND2 | β1m | β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,Δ=|Δ-Δ|=5.953×1 | ||||||||||
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,Δ=|Δ-Δ|=4.239×1 | ||||||||||
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.002×1 | ||||||||||
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,Δ=|Δ-Δ|=7.000×1 |
注: 限于篇幅,每个“双迭代”轮回中未列出第4步至倒数第2步的中间计算结果。
平曲线族及主要参数如

图8 非对称基本型平曲线“两点法”设计平曲线族及主要参数
Fig. 8 Main parameters and horizontal curve family of the “two-point method” design for asymmetric basic horizontal curve
(1) 存在规律性分布的含3段线元的非对称型平曲线族,各曲线不重叠、不交叉。
(2) 仅当比值=/在(0.5,2.0)范围内才可设置基本型平曲线,且圆曲线半径R及L、L均存在理论最小值与理论最大值;受、及α限制,L、L与R不能同时取整。
(3) 与对称型平曲线相似,非对称型平曲线族亦仅分布在中间宽、两边窄的特定区域范围内,圆曲线族位于内侧、缓和曲线族位于外侧,圆曲线族占据绝大部分区域;但平曲线族分布形状不对称,中间较宽区域偏向于总切线长较短的一侧。
(4) 当R增加时,平曲线参数均呈现规律性变化,其中LC呈增加趋势,L、L、LH呈减少趋势,LH相对变化在0.5%以内,即其近似为恒定值。
(5) L、L呈现同步变化趋势,即:若L (或L)增大,则L(或L)亦增加。
(6) 总切线长相对较大的一侧,其缓和曲线长亦较长;反之,总切线长相对较小的一侧,其缓和曲线长相对较短。

图9 某二级公路平面线形设计对比分析
Fig. 9 Comparison and analysis of horizontal alignment drawing of second-class highway
由图可知:①由ZH(直缓)、HZ(缓直)两点可构建“缓和曲线+圆曲线+缓和曲线”3段线元组成的系列平曲线族;②平曲线族各曲线参数及几何分布规律,构成两端狭窄、中间宽厚的“特定区域”,线形布设灵活多变、便于比选优选。
(1)仅需起点(直缓点)与终点(缓直点)两点即可确定由“第1缓和曲线+圆曲线+第2缓和曲线”3段线元组成的基本型平曲线;在顺次进行平面线形布设时,仅需拟定当前设计段的终点即可。
(2)仅当平曲线总切线长的比值在(0.5,2.0)范围内时,圆曲线半径、缓和曲线长度才有理论解;对称基本型的总切线长比值始终为1,故总有理论解。
(3)线形设计成果为3段线元组成的系列“平曲线族”而非单条曲线,各曲线互不交叉重叠,且分布在两侧狭窄、中间宽厚的“特定区域”内。
(4)圆曲线半径、缓和曲线长度、圆曲线长度及平曲线总长,均存在理论最小值、最大值;平曲线总长的变化值在0.5%以内,其平均值近似为固定值。
(5)当圆曲线半径增加时,圆曲线长呈增加趋势,而缓和曲线长呈减少趋势,平曲线总长亦呈减少趋势但变化很小;第1、第2缓和曲线长度具有同步增减趋势,总切线长较大一侧对应的缓和曲线长亦较长。
(6)“两点法”突显坐标点位的控制作用,呈规律性分布的“平曲线族”和“特定区域”有益于路线优选、线形拟合重构和智能化路线设计。
作者贡献声明
李玉华:提出两点法设计基本思想、计算原理,设计计算软件编制,论文撰写。
王奥凯:工程应用示例,软件调试。
刘佳音:计算原理验证,设计方法检验。
马洋洋:软件调试,文字校对。
吴树铭:软件调试。
参考文献
许金良, 程建川, 栗志海. 道路勘测设计[M].5版. 北京:人民交通出版社, 2018. [百度学术]
XU Jinliang, CHENG Jianchuan, LI Zhihai. Road survey and design[M] 5th ed. Beijing: China Communications Press, 2018. [百度学术]
吴国雄.公路平面线形曲线型设计方法综述[J].公路,1996(9): 21. [百度学术]
WU Guoxiong. A summary of curve design method for highway horizontal alignment[J]. Highway, 1996 (9): 21. [百度学术]
宋文. “附合导线法”设计匝道平面线形[J]. 公路,1991(9):11 [百度学术]
SONG Wen. Design of horizontal alignment for interchange ramps by “connecting traverse method”[J]. Highway, 1991(9):11. [百度学术]
丁建明,李方.公路平面线形设计的五单元导线法[J].东南大学学报,1998,28(2): 155. [百度学术]
DING Jianming, LI Fang. Five-element traverse method for highway horizontal alignment design[J]. Journal of Southeast University, 1998, 28(2): 155. [百度学术]
吴国雄. 一种新的公路平面线形曲线型设计方法—综合法[J].重庆交通学院学报,1996,15(增刊):59. [百度学术]
WU Guoxiong. A new curve method for designing highway horizontal alignment - comprehensive method[J]. Journal of Chongqing Jiaotong Institute, 1996, 15(S):59. [百度学术]
许金良,杨少伟,徐其福.互通式立交平面线形控制线元法[J].西安公路交通大学学报,1997,17(4):5. [百度学术]
XU Jinliang, YANG Shaowei, XU Qifu. The controlled segmentary element design method of interchange[J]. Journal of Xi’an Highway University, 1997, 17(4): 5. [百度学术]
李玉华,孙依人,刘佳音,等. 路线平面线形设计计算的两点线元法原理[J]. 同济大学学报(自然科学版),2021,49(4):507. [百度学术]
LI Yuhua, SUN Yiren, LIU Jiayin, et al. Principle of two-point segmentary element method for horizontal alignment design and calculation[J]. Journal of Tongji University (Natural Science), 2021, 49(4): 507. [百度学术]
LI W, PU H. A method for automatically recreating the horizontal alignment geometry of existing railways[J]. Computer-Aided Civil and Infrastructure Engineering, 2019, 34 (1):71. [百度学术]
CAMACHO-TORREGROSA F J, PEREZ-ZURIAGA A M, CAMPOY-UNGRIA J M, et al. Use of heading direction for recreating the horizontal alignment of an existing road[J]. Computer-Aided Civil and Infrastructure Engineering, 2015, 30(4): 282. [百度学术]
缪鹍,詹振炎.基于直线约束的道路线形设计通用方法[J].中国公路学报,2002,15(3): 15. [百度学术]
MIU Kun , ZHAN Zhenyan. Universal method for highway’s horizontal alignment based on restrictions of straight line[J]. China Journal of Highway and Transport, 2002, 15(3): 15. [百度学术]
吴国雄, 陈欣斗. 道路立交平面线形圆心连线闭合导线设计法[J]. 重庆交通学院学报,2003,22(4):20. [百度学术]
WU Guoxiong, CHEN Xindou. Closed-traverse method of center-center lines of circles in horizontal alignment design of road interchange[J]. Journal of Chongqing Jiaotong University, 2003, 22(4): 20. [百度学术]
张航,韦金君,张肖磊.道路平面线形拟合方法比较研究[J].武汉理工大学学报(交通科学与工程版),2018,42(4):594. [百度学术]
ZHANG Hang, WEI Jinjun, ZHANG Xiaolei. Comparative study on road plane alignment fitting methods[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2018, 42(4): 594. [百度学术]
BOSURGI G, D’ANDREA A. A polynomial parametric curve (PPC-CURVE) for the design of horizontal geometry of highways[J]. Computer-Aided Civil and Infrastructure Engineering, 2012, 27(4): 303. [百度学术]
徐进, 罗庆, 毛嘉川,等. 轨迹‒速度协同控制的山区公路平面线形设计方法[J]. 中国公路学报, 2013, 26(4):43. [百度学术]
XU Jin, LUO Qing, MAO Jiachuan, et al. Method for horizontal geometry design of mountainous roads based on trajectory-speed cooperative control[J]. China Journal of Highway and Transport, 2013, 26(4): 43. [百度学术]
于玲,秦翔,包龙生. 基于遗传算法的山区高速公路平面线形优化设计方法[J]. 公路交通科技,2018,165(9): 175. [百度学术]
YU Ling, QIN Xiang, BAO Longsheng. Optimization design method for horizontal alignment of mountain expressway based on genetic algorithm[J]. Journal of Highway and Transportation Research and Development, 2018, 165(9): 175. [百度学术]
韩元利,靖仕元,陈燕平. 铁路智能选线过程中平曲线最优配置模型研究[J]. 铁道工程学报,2021,275(8):1. [百度学术]
HAN Yuanli, JING Shiyuan, CHEN Yanping. Rsearch on the optimal parametric configuration of horizontal curve on railway intelligent design[J]. Journal of Railway Engineering Society, 2021, 275(8): 1. [百度学术]