76
新疆有色金属增刊1
三次抛物线缓和曲线的计算
顾大军
(新疆国地测绘工程有限责任公司
摘
要
乌鲁木齐830063)
根据三次抛物线的数学方程,给H{了在高速公路中利用三次抛物线作为缓和曲线的计算公式,包括曲线要素、中桩坐标计算,
并以一算例说明。
关键词
缓和曲线
三次抛物线
曲线要素
坐标
1
引言
入三次抛物线,以抛物线的起点(半径无穷大处)一ZH点为坐标原点,该点的切线为x轴,垂直切线方向为y轴,建立起如图所示的切线坐标系,三次抛物线的方程为:
车辆在曲线路段上行驶会产生离心力,所以在曲线上要用超高的办法来克服离心力。离心力的大小与曲线半径等因素有关,半径愈小,离心力愈大,超高也就愈大,故一定半径的曲线上应有一定量的超高。因此,在直线与圆曲线连接处,就会有如何设置超高的问题。为了解决这一问题,在直线和圆曲线之问采用一段过渡的曲线。在直线连接处,它的半径等于∞,随着距离的增加,半径逐渐减小,到与圆曲线连接处,它的半径等于圆曲线的半径R。这样随着半径的减小,使相应的超高逐渐增大,起到了过渡的作用。这种曲线称缓和曲线。
缓和曲线的线型有很多种,目前我国公路上广泛采用的缓和曲线是回旋曲线,曲线方程为A一 ̄/R・l。考虑到高速公路的特点,根据研究成果。缓和曲线也有多种形式,如三次抛物线形、三次抛物线余弦改善型、七次四项式型、半波正弦型等,考虑到三次抛物线线形简单、设计方便、平立面有效长度长、现场铺设和养护容易的特点,目前在高速公路建设中大部分采用三次抛物线形。本文主要介绍三次抛物线缓和曲线的计算。
y一恚
相似,与抛物线的长度L0,圆曲线的半径R有关。
(1)
式中:C为三次抛物线参数,和回旋曲线中的A
设P为抛物线线上任意一点,则P点处的’转向角或切线偏角(P点处的切线与ZH点的切线的夹角)B为:
.
tanp=y'-菱
P点处的曲率半径为:
p—arctan(差)
(2)
。一盟拦:!!二/夏X4::3
下面给出由弧长s计算(x,y)的方法。
(3)
在实际计算过程中,往往不知道x,只知道曲线上P点里程,即曲线弧长s=P点里程一ZH点里程。
设各点至ZH(或HZ)的弧长为si,在曲线上取
2三次抛物线型缓和fHl线方程
微小段ds,根据高等数学微积分,则有:
ds一√五_可一√l+(塞)2dx一 ̄/rF严dx
ds一√1+砑X4
dx
(4)
叔靠‰一o√i。+1嘉如
上式积分不可积,并且是个复式积分,比较难计算,我们反过来求,由于上式积分不可积,用幂级数展成多项式并舍去不影响精度项:
图l
三次抛物线型缓和曲线
则s=c,+丢一志+志一熹
如图l所示,在直线和半径为R的圆曲线间插
11苎!.一!!!!!!!!mm.m
2009年111411.2C8。
新疆有色金属
一立≥+…
X5
两端积分得.s;x+丽X5一嘉+丽X13
现在要解决的问题是知道s求x,写成方程形式:
一————I————————————————————————。!!!!
77
一
圆曲线ly=缩R短j所(a栅--213b)角.n变/为180一p’圆曲线长为:
圆曲线的切线长为:
(6’…
一s=0
v—卜一‘ro’o——o—‘’一
Xg
X13
X17
t:R.tan(与霉’
圆曲线的外矢距为:
5。40c2
1152C'。13312C6111411.2C8
(7)
e:R.sec(≤霉)一R。
由HY点的切线坐标yo可计算:
s1一y0/tanp;s2一Yo/sinp
则AB长度为:AB=t+s:,在三角形ABJD中,由正弦定理得:
B—JD=AB.sin(.180+2a-2尽>/sin(.180--a)
此式为非线性多项式方程,有多个解,根据缓和曲线使用的抛物线中一段,一个s只对应唯一的x,并且0
≤《s,所以解算此方程,可以用迭代法求x在[o,s)区
间的符合精度的唯一解,如牛顿法或埃得金迭代法・利用程序很容易实现,得到x后可根据(1)求n
当s—k时,可求得缓圆HY点的切线坐标(xo,yo)。上面的讨论是在知道C值的情况下,在实际问题中.往往给定k和R,那么如何求C呢?
实际上C是由(3)式决定的:
A—JD:AB.sinfl/sim(毕)
C・(1+蔫)亏
口
v
4
这样,插入缓和抛物线的切线长为:T一】‘D+B—JD一两,外矢距为:E—e+A—JD,曲线全长为:L=ZLo+k,
(8)
R一————』坠,_一
xo
如知道JD点的里程,各主点的里程很容易求得。
而在计算】‘o时又需要C.这实际上需要一个迭代计算,先给定初值Co,由Lo计算】‘o,由(3)计算三次抛物线k处的曲率半径p,如p≠R,需对Co进行修正,p随C的增大而增大,当p>R时,减小C;当p
<R时,增大C,直到p=R为止。
通过研究发现,C与回旋曲线中的A较接近和相似.可以设定C=Lo・K,K是一个与R比较接近的数值(当缓和曲线比较短时,K与R接近,如Lo一
490,R一9000时,K=8989.318,而Lo一150,R=9000时,K一8
999.0;再如Lo一490,R=900时,K=
777.411;而k=150,R一900时,K=889.892),所以在程序设计时,令K。一R,Co—Lo・R,计算L。处的曲率半径po,对K。修正为K。。K。=K。+(R--po),依次类推进行迭代:K;一KH+(R—pi-1),直到p—R为止。通过多个数据试验,这个迭代是收敛的,说明
此法可行。
图2
从图2看出,也可按照回旋曲线方法计算。圆曲线内移值:AR=yo+R・COSB--R切线增长值:q=xo—R・sinlB
带有三次抛物线缓和曲线的圆曲线切线长为:
3带有缓和抛物线的圆曲线要素计算
下面来计算曲线的其它要素,这里只讨论对称插入缓和抛物线。
如图2所示,两直线交点JD,路线偏角为a,在直线和圆曲线之间插入一段缓和抛物线,圆曲线必须往内侧内移AR,HY点处的切偏角B为:
T一(R+△R)・tg罢+q
厶
外矢距:E=(R+AR)sec弓“---R
厶
从后面的算例看,两者结果一致。
从上面可以看出,计算曲线要素,关键是由抛物
p=arctan(瓦X02)
(9)
线长Lo求HY点的切线坐标(Xo,Yo)。
78
顾大军:i次抛物线缓和曲线的计算增刊1
4
曲线上任意点测量坐标系中坐标计算
目前施工测量点位放样全部采用全站仪极坐标
5
算
例
在和库高速公路跨河大桥建设中,交点Jth。到Jn。的方位角为119。35731”38,J聩。到D晓的方位角为133。38’24”34,JD36的坐标(39990.5462,9675.9215),路线在
JD3。处转向角a=14。02’52,’96,里程为K7+728.5087,设
法或GPSRTK,它们的关键是求出路线在控制测量
坐标系XOY中的坐标。设计时交点在测量坐标系中坐标是已知,路线每条直线的方位角是可以得到:
fxzH=Xm—T・cosFzH—JD【YzH=YJD——T・sinFzH~JD
设任意点P的里程为LP,当ZH<LP<YH,P点位于抛物线上,计算P点到ZH点的曲线长:s=LP—ZH,由s用(7)和(1)式计算P点在切线坐标系中坐标(X。,Yp),再由坐标平移旋转公式计算P点在测量坐标系中坐标(X,,Y,):
fXP—XzH+xp・cosFzn~JD+yP・sinFz.一JDlYP=YzH+Xp・sinFzH~JD±yP・cosFzH—JD
计取三次抛物线缓和曲线长度L。为490m,而圆曲线半径为9
000
rfl。
经过计算得:
C一490.0×8989.318426T一1353.8285。E=69.1752LY=1716.2739。L一2696.2739
ZH—K6+374.6802.HZ=K9+070.9541
缓圆点HY和圆缓点YH处切偏角Bwv—pYH一
1。33739”4
P点处的切线方位角(指向路线前进方向):
FP—FzH
圆曲线内移1.1108,切线增长244.799
K6十700(40797.3152,8780.9205)K7+800(39888.7219,9696.4146)
9
JD士flP,flP—arCtan(杀)。
当P点为HY点时。S=L0,利用上式可计算出HY点的测量坐标(X。,,YH。)和FH,。
当YH<LP<HZ时,计算方法与上面相似。当HY<IJ<YH时,P点位于圆曲线上,我们以HY点为原点,HY的切线为x轴,垂直切线方向为Y轴,建立缓圆点切线坐标系。P在此坐标系中坐标为:
6结语
从上面可以看出,三次抛物线作为缓和曲线应用在高速公路中计算过程与回旋曲线相似,区别主要在根据弧线长求切线坐标方面,需采用迭代法求出。另外从算例中看出,同样长度的缓和曲线和同样的圆曲线半径,三次抛物线的切线、圆曲线长度要比回旋曲
线短。
o一吾.型,其中s—LP—HY
fxP=R・sinO
【YP=R(1一cos0)
再Fh坐标平移旋转公式计算P点在测量坐标系中坐标(XP,YP):
fXP—XHY+xP・cosFHY+yp・sinFHY
参考文献
C13李青岳,陈永奇.工程测量学.测绘出版社,1995.C23孙志忠,袁尉平,阐震初.数值分析.南京:东南大学
出版社,2002・
IWp
YHY+xP.sinFHY±y。.cosFHY
C3]李远富・线路勘测设计・北京:高等教育出版社,
2004・
上面公式中,当曲线右偏时取“+”,左偏时取“一”。
收稿:2009—06—08
(上接75页)段是含有少量的问断性完整地层的风化层,使用无同相PHP+KHm泥浆护壁,依然用西75mm绳索取心钻具钻进到403.35ITI,顺利终孔。
12线ZK34孔,倾角80。,方位角25。,设计孔深
500
5
结语
总之,后期的护壁技术、钻进工艺更加完善,钻进
比较顺利,一般孑L深400rfl左右,含有长孔段强风化
m。该孑L深241~272m,274~326.49m段见强
层的钻孔一个半月终孔,生产效率大提高,经济效益显著。我们通过不断摸索,总结经验,针对丰山矿区各孔风化程度的不同,结合设备能力、经济因素,采用不同方法护壁,总结出前述4种护壁措施,基本上解决了丰山矿区长孑L段强风化层的护壁难题。
‘
风化层。由于某些原因。致使该孑L非常复杂,钻孑L严重超径,经常发生钻杆折断事故。考虑到YL一600型钻机的能力,在使用低固相泥浆护壁、钻穿风化层
后于326.79m处下人中73mm反丝套管,用中59
mm绳索取心钻具钻至566.10m终孔。
收稿:2009—06—01
三次抛物线缓和曲线的计算
作者:作者单位:刊名:英文刊名:年,卷(期):
顾大军
新疆国地测绘工程有限责任公司,乌鲁木齐,830063新疆有色金属
XINJIANG YOUSE JINSHU2009,32(z1)
参考文献(3条)
1. 李远富 线路勘测设计 2004
2. 孙志忠;袁尉平;闻震初 数值分析 20023. 李青岳;陈永奇 工程测量学 1995
本文链接:http://d.g.wanfangdata.com.cn/Periodical_xjysjs2009z1030.aspx
76
新疆有色金属增刊1
三次抛物线缓和曲线的计算
顾大军
(新疆国地测绘工程有限责任公司
摘
要
乌鲁木齐830063)
根据三次抛物线的数学方程,给H{了在高速公路中利用三次抛物线作为缓和曲线的计算公式,包括曲线要素、中桩坐标计算,
并以一算例说明。
关键词
缓和曲线
三次抛物线
曲线要素
坐标
1
引言
入三次抛物线,以抛物线的起点(半径无穷大处)一ZH点为坐标原点,该点的切线为x轴,垂直切线方向为y轴,建立起如图所示的切线坐标系,三次抛物线的方程为:
车辆在曲线路段上行驶会产生离心力,所以在曲线上要用超高的办法来克服离心力。离心力的大小与曲线半径等因素有关,半径愈小,离心力愈大,超高也就愈大,故一定半径的曲线上应有一定量的超高。因此,在直线与圆曲线连接处,就会有如何设置超高的问题。为了解决这一问题,在直线和圆曲线之问采用一段过渡的曲线。在直线连接处,它的半径等于∞,随着距离的增加,半径逐渐减小,到与圆曲线连接处,它的半径等于圆曲线的半径R。这样随着半径的减小,使相应的超高逐渐增大,起到了过渡的作用。这种曲线称缓和曲线。
缓和曲线的线型有很多种,目前我国公路上广泛采用的缓和曲线是回旋曲线,曲线方程为A一 ̄/R・l。考虑到高速公路的特点,根据研究成果。缓和曲线也有多种形式,如三次抛物线形、三次抛物线余弦改善型、七次四项式型、半波正弦型等,考虑到三次抛物线线形简单、设计方便、平立面有效长度长、现场铺设和养护容易的特点,目前在高速公路建设中大部分采用三次抛物线形。本文主要介绍三次抛物线缓和曲线的计算。
y一恚
相似,与抛物线的长度L0,圆曲线的半径R有关。
(1)
式中:C为三次抛物线参数,和回旋曲线中的A
设P为抛物线线上任意一点,则P点处的’转向角或切线偏角(P点处的切线与ZH点的切线的夹角)B为:
.
tanp=y'-菱
P点处的曲率半径为:
p—arctan(差)
(2)
。一盟拦:!!二/夏X4::3
下面给出由弧长s计算(x,y)的方法。
(3)
在实际计算过程中,往往不知道x,只知道曲线上P点里程,即曲线弧长s=P点里程一ZH点里程。
设各点至ZH(或HZ)的弧长为si,在曲线上取
2三次抛物线型缓和fHl线方程
微小段ds,根据高等数学微积分,则有:
ds一√五_可一√l+(塞)2dx一 ̄/rF严dx
ds一√1+砑X4
dx
(4)
叔靠‰一o√i。+1嘉如
上式积分不可积,并且是个复式积分,比较难计算,我们反过来求,由于上式积分不可积,用幂级数展成多项式并舍去不影响精度项:
图l
三次抛物线型缓和曲线
则s=c,+丢一志+志一熹
如图l所示,在直线和半径为R的圆曲线间插
11苎!.一!!!!!!!!mm.m
2009年111411.2C8。
新疆有色金属
一立≥+…
X5
两端积分得.s;x+丽X5一嘉+丽X13
现在要解决的问题是知道s求x,写成方程形式:
一————I————————————————————————。!!!!
77
一
圆曲线ly=缩R短j所(a栅--213b)角.n变/为180一p’圆曲线长为:
圆曲线的切线长为:
(6’…
一s=0
v—卜一‘ro’o——o—‘’一
Xg
X13
X17
t:R.tan(与霉’
圆曲线的外矢距为:
5。40c2
1152C'。13312C6111411.2C8
(7)
e:R.sec(≤霉)一R。
由HY点的切线坐标yo可计算:
s1一y0/tanp;s2一Yo/sinp
则AB长度为:AB=t+s:,在三角形ABJD中,由正弦定理得:
B—JD=AB.sin(.180+2a-2尽>/sin(.180--a)
此式为非线性多项式方程,有多个解,根据缓和曲线使用的抛物线中一段,一个s只对应唯一的x,并且0
≤《s,所以解算此方程,可以用迭代法求x在[o,s)区
间的符合精度的唯一解,如牛顿法或埃得金迭代法・利用程序很容易实现,得到x后可根据(1)求n
当s—k时,可求得缓圆HY点的切线坐标(xo,yo)。上面的讨论是在知道C值的情况下,在实际问题中.往往给定k和R,那么如何求C呢?
实际上C是由(3)式决定的:
A—JD:AB.sinfl/sim(毕)
C・(1+蔫)亏
口
v
4
这样,插入缓和抛物线的切线长为:T一】‘D+B—JD一两,外矢距为:E—e+A—JD,曲线全长为:L=ZLo+k,
(8)
R一————』坠,_一
xo
如知道JD点的里程,各主点的里程很容易求得。
而在计算】‘o时又需要C.这实际上需要一个迭代计算,先给定初值Co,由Lo计算】‘o,由(3)计算三次抛物线k处的曲率半径p,如p≠R,需对Co进行修正,p随C的增大而增大,当p>R时,减小C;当p
<R时,增大C,直到p=R为止。
通过研究发现,C与回旋曲线中的A较接近和相似.可以设定C=Lo・K,K是一个与R比较接近的数值(当缓和曲线比较短时,K与R接近,如Lo一
490,R一9000时,K=8989.318,而Lo一150,R=9000时,K一8
999.0;再如Lo一490,R=900时,K=
777.411;而k=150,R一900时,K=889.892),所以在程序设计时,令K。一R,Co—Lo・R,计算L。处的曲率半径po,对K。修正为K。。K。=K。+(R--po),依次类推进行迭代:K;一KH+(R—pi-1),直到p—R为止。通过多个数据试验,这个迭代是收敛的,说明
此法可行。
图2
从图2看出,也可按照回旋曲线方法计算。圆曲线内移值:AR=yo+R・COSB--R切线增长值:q=xo—R・sinlB
带有三次抛物线缓和曲线的圆曲线切线长为:
3带有缓和抛物线的圆曲线要素计算
下面来计算曲线的其它要素,这里只讨论对称插入缓和抛物线。
如图2所示,两直线交点JD,路线偏角为a,在直线和圆曲线之间插入一段缓和抛物线,圆曲线必须往内侧内移AR,HY点处的切偏角B为:
T一(R+△R)・tg罢+q
厶
外矢距:E=(R+AR)sec弓“---R
厶
从后面的算例看,两者结果一致。
从上面可以看出,计算曲线要素,关键是由抛物
p=arctan(瓦X02)
(9)
线长Lo求HY点的切线坐标(Xo,Yo)。
78
顾大军:i次抛物线缓和曲线的计算增刊1
4
曲线上任意点测量坐标系中坐标计算
目前施工测量点位放样全部采用全站仪极坐标
5
算
例
在和库高速公路跨河大桥建设中,交点Jth。到Jn。的方位角为119。35731”38,J聩。到D晓的方位角为133。38’24”34,JD36的坐标(39990.5462,9675.9215),路线在
JD3。处转向角a=14。02’52,’96,里程为K7+728.5087,设
法或GPSRTK,它们的关键是求出路线在控制测量
坐标系XOY中的坐标。设计时交点在测量坐标系中坐标是已知,路线每条直线的方位角是可以得到:
fxzH=Xm—T・cosFzH—JD【YzH=YJD——T・sinFzH~JD
设任意点P的里程为LP,当ZH<LP<YH,P点位于抛物线上,计算P点到ZH点的曲线长:s=LP—ZH,由s用(7)和(1)式计算P点在切线坐标系中坐标(X。,Yp),再由坐标平移旋转公式计算P点在测量坐标系中坐标(X,,Y,):
fXP—XzH+xp・cosFzn~JD+yP・sinFz.一JDlYP=YzH+Xp・sinFzH~JD±yP・cosFzH—JD
计取三次抛物线缓和曲线长度L。为490m,而圆曲线半径为9
000
rfl。
经过计算得:
C一490.0×8989.318426T一1353.8285。E=69.1752LY=1716.2739。L一2696.2739
ZH—K6+374.6802.HZ=K9+070.9541
缓圆点HY和圆缓点YH处切偏角Bwv—pYH一
1。33739”4
P点处的切线方位角(指向路线前进方向):
FP—FzH
圆曲线内移1.1108,切线增长244.799
K6十700(40797.3152,8780.9205)K7+800(39888.7219,9696.4146)
9
JD士flP,flP—arCtan(杀)。
当P点为HY点时。S=L0,利用上式可计算出HY点的测量坐标(X。,,YH。)和FH,。
当YH<LP<HZ时,计算方法与上面相似。当HY<IJ<YH时,P点位于圆曲线上,我们以HY点为原点,HY的切线为x轴,垂直切线方向为Y轴,建立缓圆点切线坐标系。P在此坐标系中坐标为:
6结语
从上面可以看出,三次抛物线作为缓和曲线应用在高速公路中计算过程与回旋曲线相似,区别主要在根据弧线长求切线坐标方面,需采用迭代法求出。另外从算例中看出,同样长度的缓和曲线和同样的圆曲线半径,三次抛物线的切线、圆曲线长度要比回旋曲
线短。
o一吾.型,其中s—LP—HY
fxP=R・sinO
【YP=R(1一cos0)
再Fh坐标平移旋转公式计算P点在测量坐标系中坐标(XP,YP):
fXP—XHY+xP・cosFHY+yp・sinFHY
参考文献
C13李青岳,陈永奇.工程测量学.测绘出版社,1995.C23孙志忠,袁尉平,阐震初.数值分析.南京:东南大学
出版社,2002・
IWp
YHY+xP.sinFHY±y。.cosFHY
C3]李远富・线路勘测设计・北京:高等教育出版社,
2004・
上面公式中,当曲线右偏时取“+”,左偏时取“一”。
收稿:2009—06—08
(上接75页)段是含有少量的问断性完整地层的风化层,使用无同相PHP+KHm泥浆护壁,依然用西75mm绳索取心钻具钻进到403.35ITI,顺利终孔。
12线ZK34孔,倾角80。,方位角25。,设计孔深
500
5
结语
总之,后期的护壁技术、钻进工艺更加完善,钻进
比较顺利,一般孑L深400rfl左右,含有长孔段强风化
m。该孑L深241~272m,274~326.49m段见强
层的钻孔一个半月终孔,生产效率大提高,经济效益显著。我们通过不断摸索,总结经验,针对丰山矿区各孔风化程度的不同,结合设备能力、经济因素,采用不同方法护壁,总结出前述4种护壁措施,基本上解决了丰山矿区长孑L段强风化层的护壁难题。
‘
风化层。由于某些原因。致使该孑L非常复杂,钻孑L严重超径,经常发生钻杆折断事故。考虑到YL一600型钻机的能力,在使用低固相泥浆护壁、钻穿风化层
后于326.79m处下人中73mm反丝套管,用中59
mm绳索取心钻具钻至566.10m终孔。
收稿:2009—06—01
三次抛物线缓和曲线的计算
作者:作者单位:刊名:英文刊名:年,卷(期):
顾大军
新疆国地测绘工程有限责任公司,乌鲁木齐,830063新疆有色金属
XINJIANG YOUSE JINSHU2009,32(z1)
参考文献(3条)
1. 李远富 线路勘测设计 2004
2. 孙志忠;袁尉平;闻震初 数值分析 20023. 李青岳;陈永奇 工程测量学 1995
本文链接:http://d.g.wanfangdata.com.cn/Periodical_xjysjs2009z1030.aspx