CN110348128B - 一种基于k型曲线的公路平曲线设计方法 - Google Patents

一种基于k型曲线的公路平曲线设计方法 Download PDF

Info

Publication number
CN110348128B
CN110348128B CN201910630637.1A CN201910630637A CN110348128B CN 110348128 B CN110348128 B CN 110348128B CN 201910630637 A CN201910630637 A CN 201910630637A CN 110348128 B CN110348128 B CN 110348128B
Authority
CN
China
Prior art keywords
curve
tangent
shaped
point
calculated
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910630637.1A
Other languages
English (en)
Other versions
CN110348128A (zh
Inventor
潘兵宏
王开明
赵胜林
温长鹏
陈林圻
胡炜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changan University
Original Assignee
Changan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changan University filed Critical Changan University
Priority to CN201910630637.1A priority Critical patent/CN110348128B/zh
Publication of CN110348128A publication Critical patent/CN110348128A/zh
Application granted granted Critical
Publication of CN110348128B publication Critical patent/CN110348128B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Image Generation (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开一种基于K型曲线的公路平曲线设计方法,其内容包括:(1)K型曲线及k值的定义;(2)K型曲线的特点,如K型曲线在局部坐标系下的参数方程,K型曲线任意点曲率半径以及K型曲线长度。(3)K型曲线在公路平曲线交点法设计中的应用,主要对对称单曲线、非对称单曲线以及卵形曲线三种常用线型的平曲线设计要素进行计算。K型曲线与回旋线相比,可根据地形合理选用k值,以达到顺应地形、减小填挖及避免不良线型组合的目的,同时车辆运行轨迹的曲率变化率是连续的,可使道路平面线形与汽车的重心轮迹线完全重合,更加符合车辆的运行轨迹特性与驾驶员的操作特性。

Description

一种基于K型曲线的公路平曲线设计方法
技术领域
本发明属于公路路线设计领域,具体为一种利用K型曲线替代回旋线作为缓和曲线的公路平曲线设计方法。
背景技术
缓和曲线是设置在直线与圆曲线之间、半径不同的同向圆曲线之间的曲率连续变化的曲线。现行《公路路线设计规范》(JTG D20-2017)仍规定在公路平曲线设计中使用回旋线作为缓和曲线。
回旋线基本公式如下:
A2=RL
其中:
A表示回旋线参数;
R表示终点曲率半径,m;
L表示起终点回旋线长度,m。
由回旋线基本公式可以看出,回旋线任意点的曲率随弧长呈线性变化,说明回旋线的曲率是连续的,但是曲率变化率并不连续,因此回旋线的平面线形与汽车重心轮迹线并不重合,不符合车辆的运行轨迹曲率和曲率变化率均为连续的特性,不符合驾驶员的操作特性,同时,回旋线形式较为固定,在工程实际应用中灵活性较差。
发明内容
针对上述问题,本发明的目的在于提出K型曲线作为新线种,替代回旋线作为缓和曲线应用于公路上的平曲线设计,使道路线形曲率变化率连续,更加符合车辆运行轨迹和驾驶员操作特性,且设计应用上更加灵活。
本发明是通过以下技术方案来实现:
一种基于K型曲线的公路平曲线设计方法,包括:
(1)定义K型曲线:
曲线起点O,曲线终点D,该曲线上任一点P具有如下特征:
P点的弦切角θ和反向弦切角
Figure BDA0002128597740000024
满足如下关系:
Figure BDA0002128597740000021
其中:
θ为曲线上P点和曲线起点O之间的弦线与曲线起点O切线的夹角;
t为过曲线上P点的直线GH与该曲线起点O切线的夹角;
Figure BDA0002128597740000025
为曲线上P点和曲线终点D之间的弦线与曲线终点D切线的夹角;
t'为过曲线上P点的直线GH与曲线终点D切线的夹角(称为反向控制角);
k1、k2为常数,表征曲线扁平程度,采用式(2)计算:
Figure BDA0002128597740000022
其中:
a为曲线起点O的切线与曲线起点O和曲线终点D之间弦线的弦切角;
b为曲线终点D的切线与曲线起点O和曲线终点D之间弦线的反向弦切角;
βF为曲线终点D切线与曲线起点O切线的航向角;
将k1记为k,0<k<0.5,则有式(3)
Figure BDA0002128597740000023
将具有上述特征的曲线称为K型曲线;
(2)以K型曲线作为缓和曲线设计公路平曲线。
优选的,K型曲线在局部坐标系下的x、y坐标参数方程为:
Figure BDA0002128597740000031
其中:
n=1-2k;
m表示K型曲线起点和终点之间的弦长。
进一步的,K型曲线任意点曲率半径为:
Figure BDA0002128597740000032
进一步的,K型曲线作为缓和曲线时,K型曲线的长度应按式(13)计算:
Figure BDA0002128597740000033
再进一步的,步骤(2)中,当公路平曲线为对称单曲线时,K型曲线起点的曲率半径为无穷大,K型曲线终点的曲率半径为所接的圆曲线半径R;已知圆曲线半径R、平曲线转角Φ,设计方法为:
当K型曲线起点和终点之间的弦长m值已知时,根据K型曲线起点和终点的曲率半径计算公式有:
Figure BDA0002128597740000034
式中:R1为K型曲线起点曲率半径,取10000m;
R2为K型曲线终点曲率半径,取圆曲线半径R;
因为:
Figure BDA0002128597740000041
所以代入(14),得:
Figure BDA0002128597740000042
上式中,令:
Figure BDA0002128597740000043
则:
Figure BDA0002128597740000044
因此,式(15)变成:
Figure BDA0002128597740000045
因此有:
Figure BDA0002128597740000046
利用牛顿下山法解(17)得U值,进一步算得k值;将k值代入式(15)得到a值;
圆曲线的曲线内移值p与切线增长值q采用式(18)、(19)计算:
Figure BDA0002128597740000047
Figure BDA0002128597740000048
平曲线切线长T采用式(20)计算:
Figure BDA0002128597740000049
式中:Φ为平曲线转角;
平曲线长度L采用式(21)计算:
Figure BDA00021285977400000410
外距E采用式(22)计算:
Figure BDA00021285977400000411
切曲差J采用式(23)计算:
J=2T-L (23)
再进一步的,步骤(2)中,当公路平曲线为非对称单曲线时,K型曲线起点的曲率半径为无穷大,终点的曲率半径为所接的圆曲线半径R,将圆曲线两侧的K型曲线分别记为第一K型曲线和第二K型曲线,已知圆曲线半径R,两侧K型曲线参数分别为kY1和kY2,平曲线转角Φ,设计方法为:
根据式(15),计算两侧K型曲线的参数a1、a2、m1和m2
Figure BDA0002128597740000051
则平曲线两侧K型曲线内移值p1和p2由式(24)和(25)分别计算得到:
Figure BDA0002128597740000052
Figure BDA0002128597740000053
两侧K型曲线切线增长值q1和q2由式(26)和(27)分别计算:
Figure BDA0002128597740000054
Figure BDA0002128597740000055
对于平曲线两侧切线长T1和T2的计算如下:
Figure BDA0002128597740000056
Figure BDA0002128597740000057
平曲线长度L按式(32)进行计算:
Figure BDA0002128597740000058
式中:a1为第一K型曲线的起点弦切角;
a2为第二K型曲线的终点弦切角;
kY1为第一K型曲线参数;
kY2为第二K型曲线参数;
LS1为第一K型曲线长度;
LS2为第二K型曲线长度;
外距E按下式(33)计算:
Figure BDA0002128597740000061
E=(R+p1)Sin(δ1)-R (33)
切曲差J可采用式(34)计算:
J=T1+T2-L (34)
再进一步的,步骤(2)中,当公路平曲线为卵形曲线时,K型曲线起点的曲率半径为R1,终点的曲率半径为所接的圆曲线半径R2,已知两圆曲线半径分别为R1和R2,两圆曲线所对应的曲线转角分别为Φ1和Φ2,两侧K型曲线的参数kR1和kR2值、中间K型曲线的kf值,切基线长度lCD,设计方法为:
根据式(15)计算得到两侧K型曲线的参数a1、a2、m1和m2,以及中间K型曲线的参数af和mf
Figure BDA0002128597740000062
两侧圆曲线内移值p1和p2以及切线增长值q1和q2分别由式(24)~式(27)计算;
Figure BDA0002128597740000071
Figure BDA0002128597740000072
Figure BDA0002128597740000073
Figure BDA0002128597740000074
对于第一圆曲线侧切线T1,被第一圆曲线终点分割的切基线长T2和T3,以及第二圆曲线侧切线T4,其计算过程分别如下:
Figure BDA0002128597740000075
T3=lCD-T2 (37)
Figure BDA0002128597740000076
Figure BDA0002128597740000077
Figure BDA0002128597740000078
Figure BDA0002128597740000079
Figure BDA00021285977400000710
其中,MN为两个圆曲线圆心之间的距离,N′N为圆心N到圆心M与第一圆曲线终点连线的垂直距离,H′D′为第二圆曲线发生内移前的切线长。
令:R1-R2-p2=RP;T3-T′3=ΔT,
整理式(38)~(47)得:
[1+Ctan22)]μ2+2(ΔT·Ctan(Ф2)-RP)μ+RP2-MN2+ΔT2=0 (48)
解方程(48)得μ,代入式(38),计算得到T4
平曲线长度L按式(51)进行计算:
Figure BDA00021285977400000711
切曲差J采用式(52)计算:
Figure BDA0002128597740000081
与现有技术相比,本发明具有以下有益的技术效果:
与传统的回旋线相比,K型曲线的弦切角与其对应的航向角的比值为一参数k(0<k<0.5)。即该缓和曲线航向角随弧长的比值可根据需要进行变化,在路线平面线型设计中可根据地形合理灵活地选用k值,以达到顺应地形、减小填挖及避免不良线型组合的目的。且K型曲线的参数方程二阶可导,即将K型曲线作为公路缓和曲线时,车辆运行轨迹的曲率变化率是连续的,可使道路平面线形与汽车的重心轮迹线完全重合。与传统回旋线相比更加符合车辆的运行轨迹特性与驾驶员的操作特性。
附图说明
图1为K型曲线示意图;
图2为Csc[a+nt]3(nSin[a-t]+Sin[a+nt])函数图;
图3为对称单曲线计算示意图;
图4为非对称单曲线计算示意图;
图5为卵形曲线计算示意图。
图6为K型曲线应用实例
具体实施方式
下面结合具体的实施例对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
K型曲线定义如下:
图1中,曲线起点O的曲率半径为R1(为无穷大时,表示直线),曲线终点D的曲率半径为R2,曲线的长度为LS,该曲线上任一点P具有如下特征:
P点与弦切角和反向弦切角(图1)之间满足如下关系:
Figure BDA0002128597740000091
其中:
θ为P点的弦切角,即曲线上P点和曲线起点O之间的弦线与曲线起点O切线的夹角,取值范围[0,a];
t为过曲线上P点的直线GH与该曲线起点O切线的夹角(称为控制角);
Figure BDA0002128597740000092
为P点的反向弦切角,即曲线上P点和曲线终点D之间的弦线与曲线终点D切线的夹角,取值范围[b,0];
t′为过曲线上P点的直线GH与曲线终点D切线的夹角(称为反向控制角);
k1、k2为常数,表征曲线扁平程度,采用式(2)计算:
Figure BDA0002128597740000093
其中:
a为曲线起点O的切线与起终点OD之间的弦切角(弧度);
b为曲线终点D的切线与起终点OD之间的反向弦切角(弧度);
βF为曲线的总航向角,即曲线终点D切线与曲线起点O切线的航向角(弧度)。
ΔODE中外角βF=a+b,所以k1+k2=1,将k1记为k,则
Figure BDA0002128597740000094
将具有上述特征的曲线称为K型曲线。
K型曲线可替代回旋线作为另一种形式的缓和曲线应用于公路平面线形设计,其中k表示缓和曲线上任意点P的弦切角θ(或
Figure BDA0002128597740000101
)与航向角t(或t′)的比值参数,从而被定义为K型曲线,即
Figure BDA0002128597740000102
K型曲线其特征为:
①K型曲线在局部坐标系下的x、y坐标参数方程推导如下:
如图1在ΔPOD中,根据正弦定理有:
Figure BDA0002128597740000103
式中:
LOP为曲线上O点与任意点P之间的矢径(距离);
LOD为O点和D点之间的距离,即曲线的弦长(m)。
于图1建立局部直角坐标系统xoy。
因为在ΔPOD中,根据正弦定理有:
Figure BDA0002128597740000104
式中:
LOP为曲线上O点与任意点P之间的矢径(距离);
LOD为O点和D点之间的距离,即曲线的弦长(m)。
因此,任一点P的坐标用下述方程描述:
Figure BDA0002128597740000105
又因为
Figure BDA0002128597740000106
因此得:
Figure BDA0002128597740000111
令:
m=LOD
d=1-k
n=1-2k
则:k=d-n
所以:
Figure BDA0002128597740000112
整理后有:
Figure BDA0002128597740000113
而2d-n=2-2k-1+2k=1,所以:
Figure BDA0002128597740000114
上式就是K型曲线在局部坐标系下的参数方程。
②K型曲线任意点的曲率半径随弧长的变化关系为:
K型曲线上任意一点的曲率半径r为:
Figure BDA0002128597740000115
式中:
x′、x″分别为x的一阶、二阶导数(对角度t求导);
y′、y″分别为y的一阶、二阶导数(对角度t求导)。
所以
Figure BDA0002128597740000121
对上式的分母进行分析:
(1)为满足曲线半径从大圆渐变到小圆的要求,n应大于零,即要求:
0<k<0.5
而:n=1-2k,
因此:0<n<1
所以:(n2-1)<0
又因为:Csc[a+nt]3(nSin[a-t]+Sin[a+nt])
上式中:t∈[0,π],
Figure BDA0002128597740000124
在该区间内,函数的图形如图2所示;
所以:Csc[a+nt]3(nSin[a-t]+Sin[a+nt])>0
所以:m2(n2-1)Csc[a+nt]3(nSin[a-t]+Sin[a+nt])<0
因此K型曲线上任意点的曲率半径采用下式计算:
Figure BDA0002128597740000122
由式(11)可知,K型曲线任意点曲率半径与弧长具有三角函数关系,且二阶可导,因此K型曲线的曲率变化率连续,与汽车的重心轮迹线完全重合。与传统回旋线相比更加符合车辆的运行轨迹特性与驾驶员的操作特性。
③根据K型曲线的参数方程,k型曲线为连续的光滑曲线,一阶导数存在,因此,根据弧长积分原理可以求弧长。
根据微分弧长公式:
Figure BDA0002128597740000123
所以曲线弧长公式为:
Figure BDA0002128597740000131
式中:
x′为x的一阶导数;
y′为y的一阶导数。
代入后,得到:
Figure BDA0002128597740000132
上式无法积分,因此采用不定积分的方法不可行,可借助软件采用数值积分的方法进行求解。
在实际应用中,K型曲线的曲线长度与弦长的差异不太大,为控制曲线长度,可以用弦长来控制曲线长度,即在设计中,可先给出弦长m。
基于上述K型曲线各参数的推导,现给出三种以K型曲线代替回旋线作为缓和曲线的平曲线常用设计方法。
一、对称单曲线
图3中,在对称单曲线设计中缓和曲线采用K型曲线时,K型曲线起点的曲率半径为无穷大(∞),在公路中,认为半径大于等于10000m时,与直线没有区别,因此,取K型曲线起点的曲率半径R1=10000m;K型曲线终点的曲率半径为所接的圆曲线半径R。
当K型曲线的参数m值已知时,根据K型曲线起点和终点的曲率半径计算公式有:
Figure BDA0002128597740000133
式中:R1为K型曲线起点曲率半径,取10000m;
R2为K型曲线终点曲率半径,R。
因为:
Figure BDA0002128597740000141
所以代入(14),得:
Figure BDA0002128597740000142
上式中,令:
Figure BDA0002128597740000143
则:
Figure BDA0002128597740000144
因此,上式可以变成:
Figure BDA0002128597740000145
因此有:
Figure BDA0002128597740000146
利用牛顿下山法解(17)可得U值,进一步可算的k值;将k值代入式(15)得到a值。
根据图1,圆曲线曲线内移值p与切线增长值q可采用式(18)、(19)计算:
Figure BDA0002128597740000147
Figure BDA0002128597740000148
根据图3,平曲线切线长T可采用式(20)计算:
Figure BDA0002128597740000149
式中:Φ为交点的偏角(Rad),即平曲线转角。
平曲线长度L可采用式(21)计算:
Figure BDA00021285977400001410
式中:β0为K型曲线终点的航向角,采用下式计算:
Figure BDA0002128597740000151
外距E可采用式(22)计算:
Figure BDA0002128597740000152
切曲差J可采用式(23)计算:
J=2T-L (23)
二、非对称单曲线
如图4,非对称单曲线为单圆曲线两侧设置了非对称的缓和曲线,记为第一K型曲线和第二K型曲线,两侧K型曲线的参数k不一样。K型曲线起点的曲率半径为无穷大(∞),终点的曲率半径为所接的圆曲线半径R。
当两侧K型曲线的参数kY1和kY2值已知时,根据式(15)和(16),可以计算两侧K型曲线的参数a1、a2和m1、m2。则可以计算非对称单曲线的曲线要素。
则公路平曲线两侧缓和曲线内移值p1、p2可由式(24)、(25)分别计算:
Figure BDA0002128597740000153
Figure BDA0002128597740000154
两侧缓和曲线切线增长值q1、q2可由式(26)、(27)分别计算:
Figure BDA0002128597740000155
Figure BDA0002128597740000156
对于平曲线两侧切线长T1和T2的计算。由图4可知,OA、OB是过圆心O,垂直于交点线的垂线。通过圆心O分别作交点线的平行线,与交点线的交点分别为D、F,即:OD∥CF、OF∥CD。则有:
∠ODA=∠OFB=Φ (28)
由图4可知:
T1=CD-DA+q1 (29)
式中:
Figure BDA0002128597740000161
Figure BDA0002128597740000162
所以:
Figure BDA0002128597740000163
同理可以推导出切线长度T2
Figure BDA0002128597740000164
而平曲线长L按下式进行计算:
L=(Φ-β0102)R+LS1+LS2
式中:β01、β02为K型曲线终点的航向角,采用下式计算:
Figure BDA0002128597740000171
所以平曲线长度L:
Figure BDA0002128597740000172
式中:a1为非对称单曲线第一K型曲线的起点(ZH)弦切角(rad);
a2为非对称单曲线第二K型曲线的终点(HZ)弦切角(rad);
kY1为非对称单曲线第一K型曲线参数;
kY2为非对称单曲线第二K型曲线参数;
LS1为非对称单曲线第一K型曲线长度(m);
LS2为非对称单曲线第二K型曲线长度(m);
对于外距E,因非对称单曲线两边切线不等长,曲线中点可取圆曲线中点或全曲线中点。为计算和测设方便,可取交点和圆心的连线与圆曲线的交点作为曲线中点(QZ)(图4),其要素按下式(33)计算:
Figure BDA0002128597740000173
E=(R+p1)Sin(δ1)-R (33)
对于切曲差J,可采用式(34)计算:
J=T1+T2-L (34)
三、卵形曲线
用一段缓和曲线连接两个同向圆曲线的组合形式称为卵形曲线。缓和曲线采用K型曲线时,其起点的曲率半径为R1,终点的曲率半径为所接的圆曲线半径R2
如图5,当两侧和中间缓和曲线采用K型曲线,且两侧K型曲线的参数kR1和kR2值、中间K型曲线的kf值已知时,根据式(15)和式(16),可以计算两侧K型曲线的参数a1、a2和m1、m2,中间K型曲线的参数af和mf。则可计算卵形曲线的曲线要素。
两侧圆曲线内移值p1、p2和切线增长值q1、q2可分别由式(24)~式(27)计算。
对于第一圆曲线侧切线T1、被第一圆曲线终点分割的切基线长T2+T3,其计算过程分别如下,如图5示,MG′、NH′是过圆心M、N垂直基线CD的垂线,其长度分别为R1+p1、R2+p2,且MG′//NH′。则可知内移前的圆弧
Figure BDA0002128597740000181
分别与BC、C′G′和H′D′、D′E相切,则有:
Figure BDA0002128597740000182
而:
Figure BDA0002128597740000183
G′G=p1
Figure BDA0002128597740000184
所以:
Figure BDA0002128597740000185
已知基线CD的长度为lCD,由图5可知:
T3=lCD-T2 (37)
而对于T4的计算,则由图5可知,
Figure BDA0002128597740000191
又由图5可知:
μ=HH′=GG″=MG-MN′-N′G″′-G″G″′ (39)
而:
MG=R1
N′G″′=R2
G″G″′=H′H″=p2
所以:
μ=R1-MN′-R2-p2 (40)
又因为:
Figure BDA0002128597740000192
所以:
Figure BDA0002128597740000193
而:
Figure BDA0002128597740000194
Figure BDA0002128597740000195
因此:
Figure BDA0002128597740000196
而两个圆心之间的距离MN为:
MN2=(xm-xn)2+(ym-yn)2 (43)式中两个圆心的坐标采用中间缓和曲线上曲率半径分别为R1、R2的两个圆的圆心坐标:
Figure BDA0002128597740000201
Figure BDA0002128597740000202
式中:
kf为中间K型曲线设计的参数;
mf、af为中间K型曲线设计的参数,根据下列方程组求解:
Figure BDA0002128597740000203
所以:
Figure BDA0002128597740000204
令:R1-R2-P2=RP;T3-T′3=ΔT
整理上式得:
[1+Ctan22)]μ2+2(ΔT·Ctan(Φ2)-RP)μ+RP2-MN2+ΔT2=0 (48)
解方程(48)得μ。代入式(38),即可计算T4
平曲线长L可按式(51)进行计算:
从图5可知:
Figure BDA0002128597740000205
式中:β1、β2为两侧K型曲线终点的航向角,采用下式计算:
Figure BDA0002128597740000211
所以卵形曲线的长度L为:
Figure BDA0002128597740000212
对于切曲差J,可采用式(52)计算:
Figure BDA0002128597740000213
对本发明中利用K型曲线进行对称单曲线、非对称单曲线和卵形曲线设计方法进行总结如下。
对于对称单曲线的设计,其设计步骤如下:
步骤1:已知圆曲线半径R、起终点弦长m、平曲线转角Φ;
步骤2:利用式(17)求解k,并带入式(15)求解a;
步骤3:利用式(18)和式(19)求解曲线内移值p和切线增长值q;
步骤4:利用式(20)~(23)分别求解出切线长T,平曲线长L,外距E和切曲差J。
对于非对称单曲线的设计,其计算步骤如下:
步骤1:已知圆曲线半径R,两侧K型曲线参数分别为kY1和kY2,平曲线转角Φ;
步骤2:利用式(15)分别求解两侧K型曲线的参数a1、a2和m1、m2
步骤3:利用式(24)~(27)分别求解两侧K型曲线内移值p1和p2,切线增长值q1和q2
步骤4:利用式(30)~(31)分别求解切线长T1、T2
步骤5:利用式(32)~(34)分别求解平曲线长L、外距E和切曲差J。
对于卵形曲线设计,其计算步骤如下:
步骤1:已知两端圆曲线半径R1、R2,两圆曲线所对应曲线转角Φ1和Φ2,两侧K型曲线的参数kR1和kR2值、中间K型曲线的kf值,切基线长度lCD
步骤2:根据式(15)和式(16),可以计算两侧K型曲线的参数a1、a2和m1、m2,中间K型曲线的参数af和mf
步骤3:利用式(24)~(27)分别求解两侧K型曲线内移值p1和p2,切线增长值q1和q2
步骤4:根据式(36)~(37)分别求解切线长T1、T2和T3
步骤5:利用式(44)~(45)求解两圆心坐标xm、ym和xn、yn
步骤6:利用式(43)求解圆心距离MN;
步骤7:利用式(48)求解μ;
步骤8:利用式(38)求解切线长T4
步骤9:利用式(51)和(52)分别求解平曲线长L和切曲差J。
计算实例
本发明以某三级公路为例,对K型曲线在公路平曲线中的应用特点加以说明。如图6所示,该三级路为山区道路,路线为克服高差,需顺应山区地形进行自然展线,并利用图6中山脊地形进行大角度转向,图中点划线为导向线,路线中线沿着导向线前进可有效控制填挖,并在该处设置交点JD1,实线为利用回旋线作为缓和曲线下的设计线,相应设计指标有:交点转角Φ=114°44’10.0”,圆曲线半径R=40m,回旋线长度LH=40m。
因此设置回旋线为缓和曲线下平曲线测设元素计算结果如下:
圆曲线内移值pH和切线增长值qH为:
Figure BDA0002128597740000231
Figure BDA0002128597740000232
缓和曲线角βH
Figure BDA0002128597740000233
切线长TH
Figure BDA0002128597740000234
平曲线长LPH
Figure BDA0002128597740000235
外距EH
Figure BDA0002128597740000236
切曲差JH
JH=2TH-LPH=49.691m
图6中虚线为以K型曲线作为缓和曲线下的平曲线,为与回旋线做对比,与回旋线对应平曲线相同的圆曲线半径R=40m,交点转角Φ=114°44’10.0”,并采用回旋线长度LH对应弦长m=39.557m,由此计算的K型曲线作为缓和曲线下平曲线测设元素为:
利用式(17)求解得K型曲线k值为:k=0.137
并根据式(15)求得a=0°42’59.6”
由式(18)和(19)分别解得圆曲线内移值pK和切线增长值qK
pK=0.3275m,qK=35.919
根据式(20)~(23)分别求得,
切线长TK=98.900m
平曲线长LK=151.948m
外距EK=34.786m
切曲差JK=45.852m
结合图6可以看出,与回旋线相比,K型曲线对应平曲线可更快地从直线过渡至圆曲线,且外距EK小于回旋线对应平曲线外距EH,即K型曲线对应平曲线更靠近导向线,在相同的平面要素指标下能减少更多的挖方,因此K型曲线在该工程条件下更为适用。

Claims (4)

1.一种基于K型曲线的公路平曲线设计方法,其特征在于,包括以下步骤:
(1)定义K型曲线:
曲线起点O,曲线终点D,该曲线OD上任一点P具有如下特征:
P点的弦切角θ和反向弦切角
Figure FDA0003830591140000011
满足如下关系:
Figure FDA0003830591140000012
其中:
θ为曲线上P点和曲线起点O之间的弦线与曲线起点O切线的夹角;
t为过曲线上P点的直线GH与该曲线起点O切线的夹角;
Figure FDA0003830591140000013
为曲线上P点和曲线终点D之间的弦线与曲线终点D切线的夹角;
t'为过曲线上P点的直线GH与曲线终点D切线的夹角;
k1、k2为常数,表征曲线扁平程度,采用式(2)计算:
Figure FDA0003830591140000014
其中:
a为曲线起点O的切线与曲线起点O和曲线终点D之间弦线的弦切角;
b为曲线终点D的切线与曲线起点O和曲线终点D之间弦线的反向弦切角;
βF为曲线终点D切线与曲线起点O切线的航向角;
将k1记为k,0<k<0.5,则有式(3)
Figure FDA0003830591140000015
将满足上述条件的曲线OD称为K型曲线;
(2)以K型曲线作为缓和曲线设计公路平曲线;
K型曲线在局部坐标系下的x、y坐标参数方程为:
Figure FDA0003830591140000021
其中:
n=1-2k;
m表示K型曲线起点和终点之间的弦长;
K型曲线作为缓和曲线时,K型曲线的长度应按式(13)计算:
Figure FDA0003830591140000022
所述步骤(2)中,当公路平曲线为对称单曲线时,K型曲线起点的曲率半径为无穷大,K型曲线终点的曲率半径为所接的圆曲线半径R;已知圆曲线半径R、平曲线转角Φ,设计方法为:
当K型曲线起点和终点之间的弦长m值已知时,根据K型曲线起点和终点的曲率半径计算公式有:
Figure FDA0003830591140000023
式中:R1为K型曲线起点曲率半径,取10000m;
R2为K型曲线终点曲率半径,取圆曲线半径R;
因为:
Figure FDA0003830591140000024
所以代入(14),得:
Figure FDA0003830591140000025
上式中,令:
Figure FDA0003830591140000026
则:
Figure FDA0003830591140000027
因此,式(15)变成:
Figure FDA0003830591140000031
因此有:
Figure FDA0003830591140000032
利用牛顿下山法解(17)得U值,进一步算得k值;将k值代入式(15)得到a值;
圆曲线的曲线内移值p与切线增长值q采用式(18)、(19)计算:
Figure FDA0003830591140000033
Figure FDA0003830591140000034
平曲线切线长T采用式(20)计算:
Figure FDA0003830591140000035
式中:Φ为平曲线转角;
平曲线长度L采用式(21)计算:
Figure FDA0003830591140000036
外距E采用式(22)计算:
Figure FDA0003830591140000037
切曲差J采用式(23)计算:
J=2T-L (23)。
2.根据权利要求1所述的基于K型曲线的公路平曲线设计方法,其特征在于,K型曲线任意点曲率半径为:
Figure FDA0003830591140000038
3.根据权利要求1所述的基于K型曲线的公路平曲线设计方法,其特征在于,所述步骤(2)中,当公路平曲线为非对称单曲线时,K型曲线起点的曲率半径为无穷大,终点的曲率半径为所接的圆曲线半径R,将圆曲线两侧的K型曲线分别记为第一K型曲线和第二K型曲线,已知圆曲线半径R,两侧K型曲线参数分别为kY1和kY2,平曲线转角Φ,设计方法为:
根据式(15),计算两侧K型曲线的参数a1、a2、m1和m2
Figure FDA0003830591140000041
则平曲线两侧K型曲线内移值p1和p2由式(24)和(25)分别计算得到:
Figure FDA0003830591140000042
Figure FDA0003830591140000043
两侧K型曲线切线增长值q1和q2由式(26)和(27)分别计算:
Figure FDA0003830591140000044
Figure FDA0003830591140000045
对于平曲线两侧切线长T1和T2的计算如下:
Figure FDA0003830591140000046
Figure FDA0003830591140000047
平曲线长度L按式(32)进行计算:
Figure FDA0003830591140000048
式中:a1为第一K型曲线的起点弦切角;
a2为第二K型曲线的终点弦切角;
kY1为第一K型曲线参数;
kY2为第二K型曲线参数;
LS1为第一K型曲线长度;
LS2为第二K型曲线长度;
外距E按下式(33)计算:
Figure FDA0003830591140000051
E=(R+p1)Sin(δ1)-R (33)
切曲差J可采用式(34)计算:
J=T1+T2-L (34)。
4.根据权利要求1所述的基于K型曲线的公路平曲线设计方法,其特征在于,步骤(2)中,当公路平曲线为卵形曲线时,K型曲线起点的曲率半径为R1,终点的曲率半径为所接的圆曲线半径R2,已知两圆曲线半径分别为R1和R2,两圆曲线所对应的曲线转角分别为Φ1和Φ2,两侧K型曲线的参数kR1和kR2值、中间K型曲线的kf值,切基线长度lCD,设计方法为:
根据式(15)计算得到两侧K型曲线的参数a1、a2、m1和m2,以及中间K型曲线的参数af和mf
Figure FDA0003830591140000052
两侧圆曲线内移值p1和p2以及切线增长值q1和q2分别由式(24)~式(27)计算;
Figure FDA0003830591140000053
Figure FDA0003830591140000054
Figure FDA0003830591140000055
Figure FDA0003830591140000061
对于第一圆曲线侧切线T1,被第一圆曲线终点分割的切基线长T2和T3,以及第二圆曲线侧切线T4,其计算过程分别如下:
Figure FDA0003830591140000062
T3=lCD-T2 (37)
Figure FDA0003830591140000063
Figure FDA0003830591140000064
Figure FDA0003830591140000065
Figure FDA0003830591140000066
Figure FDA0003830591140000067
其中,MN为两个圆曲线圆心之间的距离,N′N为圆心N到圆心M与第一圆曲线终点连线的垂直距离,H′D′为第二圆曲线发生内移前的切线长;
令:R1-R2-p2=RP;T3-T′3=ΔT,
整理式(38)~(47)得:
[1+Ctan22)]μ2+2(ΔT·Ctan(Φ2)-RP)μ+RP2-MN2+ΔT2=0 (48)
解方程(48)得μ,代入式(38),计算得到T4
平曲线长度L按式(51)进行计算:
Figure FDA0003830591140000068
切曲差J采用式(52)计算:
Figure FDA0003830591140000069
CN201910630637.1A 2019-07-12 2019-07-12 一种基于k型曲线的公路平曲线设计方法 Active CN110348128B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910630637.1A CN110348128B (zh) 2019-07-12 2019-07-12 一种基于k型曲线的公路平曲线设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910630637.1A CN110348128B (zh) 2019-07-12 2019-07-12 一种基于k型曲线的公路平曲线设计方法

Publications (2)

Publication Number Publication Date
CN110348128A CN110348128A (zh) 2019-10-18
CN110348128B true CN110348128B (zh) 2022-10-14

Family

ID=68175112

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910630637.1A Active CN110348128B (zh) 2019-07-12 2019-07-12 一种基于k型曲线的公路平曲线设计方法

Country Status (1)

Country Link
CN (1) CN110348128B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111765875A (zh) * 2020-07-13 2020-10-13 江苏东交智控科技集团股份有限公司 圆曲线测量放线方法
CN113643407A (zh) * 2021-08-18 2021-11-12 北京赛目科技有限公司 两点间道路曲线确定方法、装置、电子设备及存储介质
CN114818071B (zh) * 2022-04-29 2024-04-26 中国铁路设计集团有限公司 一种铁路线路平面曲线的计算及绘制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106245477A (zh) * 2016-08-25 2016-12-21 山东省交通规划设计院 基于极轴计算法的公路匝道接线设计方法
WO2017113416A1 (zh) * 2015-12-31 2017-07-06 深圳配天智能技术研究院有限公司 加工轨迹平滑转接的方法及加工装置
WO2017113195A1 (zh) * 2015-12-30 2017-07-06 深圳配天智能技术研究院有限公司 一种加工路径规划方法、加工路径规划装置及数控机床
CN108520120A (zh) * 2018-03-19 2018-09-11 大连理工大学 一种道路路线平面线形设计的“两点”法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017113195A1 (zh) * 2015-12-30 2017-07-06 深圳配天智能技术研究院有限公司 一种加工路径规划方法、加工路径规划装置及数控机床
WO2017113416A1 (zh) * 2015-12-31 2017-07-06 深圳配天智能技术研究院有限公司 加工轨迹平滑转接的方法及加工装置
CN106245477A (zh) * 2016-08-25 2016-12-21 山东省交通规划设计院 基于极轴计算法的公路匝道接线设计方法
CN108520120A (zh) * 2018-03-19 2018-09-11 大连理工大学 一种道路路线平面线形设计的“两点”法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
公路立交匝道卵型平面曲线线形设计改进方法;景天虎;《山西建筑》;20070510(第14期);全文 *

Also Published As

Publication number Publication date
CN110348128A (zh) 2019-10-18

Similar Documents

Publication Publication Date Title
CN110348128B (zh) 一种基于k型曲线的公路平曲线设计方法
CN102918312B (zh) 一种改进的管间密封
CN108520120B (zh) 一种道路路线平面线形设计的“两点”法
EP2108911B1 (en) Heat exchanger
CN113177256B (zh) 一种道路曲线的曲线要素通用获取方法
US8425212B2 (en) Positive displacement flowmeter and helical gear
CN109297501B (zh) 一种gis单线路网自动生成双线路网的方法
CN104540997B (zh) 具有缓和曲线的道路弯道及其形成方法
CN110222395A (zh) 一种纤维缠绕复合材料壳体封头厚度计算方法
US20070179698A1 (en) Traction control method for a tracked vehicle
CN108716532B (zh) 多段耦合型曲线齿轮齿形及其设计方法
CN103950544B (zh) 通用亚声速扩压器设计方法
CN102900672B (zh) 无油螺杆压缩机的转子
CN105127521A (zh) 一种用于加工双渐开线齿轮的滚刀及加工方法
CN108757447B (zh) 一种分段圆弧型罗茨转子及其型线设计方法
CN110802280A (zh) 一种渐开线螺旋锥齿轮齿面设计方法
CN103438180A (zh) 一种定传动比点接触曲线齿锥齿轮的修形方法
CN107273616B (zh) 一种基于“交点法”的高速铁路相邻线路平面设计方法
USRE29115E (en) Contact gearing
CN208749866U (zh) 一种形状记忆抱闸环
CN207147305U (zh) 用于换热器管箱的密封导向垫及换热器管箱
WO2001018406A1 (en) Improved tube for heat exchangers
Zhang et al. Controllable Clothoid Path Generation for Autonomous Vehicles
CN110954123A (zh) 一种基于阿克曼约束的路径规划方法
CN103089616A (zh) 一种内啮合齿廓副

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant