缓和曲线参数再推导
〇、前言一、缓和曲线特性二、缓和曲线参数计算
〇、前言
在20世纪90年代之前,我国的高速公路还没有大规模建设,公路建设基本以二、三级甚至更低等级的公路为主,因此测量放样采用以前教科书上的切线支距法、偏角法等已经可以满足测量精度和工期的需要。随着20世纪80年代末90年代初高速公路建设事业的兴起,公路平面线型设计多样化,线型组合复杂,加之建设进度要求加快,相应精度要求也比普通公路要高,切线支距法、偏角法等已经不能满足测量精度和工期的需求了。因此必须采用坐标法精确计算平曲线上任意点的坐标,然后再采用坐标法进行放样。在计算坐标的过程当中,直线和圆曲线的计算方法简单,难点在于计算缓和曲线上任意点的坐标。在我国的公路、铁路工程平曲线设计中,普遍采用辐射螺旋线(又称回旋曲线)作为缓和曲线的线型计算模型(因此以下的缓和曲线均指回旋曲线)。为了计算缓和曲线上任意点的坐标,我们有必要从缓和曲线特性来推导出其代表参数的计算公式,进而为坐标的精确计算提供基础依据。
一、缓和曲线特性
对于完整缓和曲线,设起点为ZH点,其半径r=∞,设终点为HY点,其半径r=R,起终点间缓和曲线全长为Ls;设p为缓和曲线上任意一点,曲率半径为r,该点至起点的曲线长度为l,根据回旋线特性:回旋线是半径与曲线长度成反比的曲线,则rl=RLs=C。其中C为常数,称为回旋线半径变化率,为今后计算方便,引入回旋曲线参数A,令A2=C。即当缓和曲线终点圆曲线半径(回旋线曲率半径)R及回旋线长度Ls已知时,C及A即可唯一确定。 ![缓和曲线-回旋线要素示意图](https://img-blog.csdnimg.cn/20190710160425147.png)
二、缓和曲线参数计算
回旋曲线参数A的计算:根据前述缓和曲线特性,已经得出回旋曲线参数A的计算公式即:A2=C=rl=RLs。回旋线中心角(缓和曲线角)的计算:根据几何关系可知,回旋线上任意点p与缓和曲线起点之间的曲线长度l所对应的曲线中心角β即切线角,与p点切线与起点切线之间的夹角相等。在该p点取一微分弧段dl,所对应的中心角为dβ,于是有: ~~
d
β
=
d
l
r
=
l
d
l
C
=
l
d
l
R
L
s
⇒
β
=
∫
l
d
l
R
L
s
=
l
2
2
R
L
s
d\beta =\frac{dl}{r}=\frac{ldl}{C}=\frac{ldl}{RLs}\Rightarrow \beta =\int \frac{ldl}{RLs}=\frac{l^2}{2RLs}
dβ=rdl=Cldl=RLsldl⇒β=∫RLsldl=2RLsl2 ~~ 则当l=Ls即缓和曲线终点处对应的缓和曲线圆心角为: β0=
L
s
2
2
R
L
s
=
L
s
2
R
=
A
2
2
R
2
=
L
s
2
2
A
2
\frac{Ls^2}{2RLs}=\frac{Ls}{2R}=\frac{A^2}{2R^2}=\frac{Ls^2}{2A^2}
2RLsLs2=2RLs=2R2A2=2A2Ls2局部坐标参数计算 : 建立以ZH点为原点,过该点的切线为x轴,法线为y轴的坐标系,则回旋线上任意点p的坐标为x,y,在该p点取一微分弧段dl,其在坐标轴上的投影分别为dx,dy,则有:
d
x
=
d
l
⋅
c
o
s
β
=
r
d
β
⋅
c
o
s
β
=
A
2
β
⋅
c
o
s
β
⇒
x
=
A
2
⋅
∫
c
o
s
β
β
d
β
dx=dl·cosβ=rdβ·cosβ=\frac{A}{\sqrt{2β}}·cosβ\Rightarrow x=\frac{A}{\sqrt{2}}·∫\frac{cosβ}{\sqrt{β}}dβ
dx=dl⋅cosβ=rdβ⋅cosβ=2β
A⋅cosβ⇒x=2
A⋅∫β
cosβdβ
d
y
=
d
l
⋅
s
i
n
β
=
r
d
β
⋅
s
i
n
β
=
A
2
β
⋅
s
i
n
β
⇒
y
=
A
2
⋅
∫
s
i
n
β
β
d
β
dy=dl·sinβ=rdβ·sinβ=\frac{A}{\sqrt{2β}}·sinβ\Rightarrow y=\frac{A}{\sqrt{2}}·∫\frac{sinβ}{\sqrt{β}}dβ
dy=dl⋅sinβ=rdβ⋅sinβ=2β
A⋅sinβ⇒y=2
A⋅∫β
sinβdβ 而cos(β)和sin(β)的级数展开公式分别为:
c
o
s
(
β
)
=
1
−
β
2
2
!
+
β
4
4
!
+
…
…
+
(
−
1
)
i
⋅
β
2
i
(
2
i
)
!
cos(β)=1-\frac{β^2}{2!}+\frac{β^4}{4!}+……+(-1)^i·\frac{β^{2i}}{(2i)!}
cos(β)=1−2!β2+4!β4+……+(−1)i⋅(2i)!β2i(注:i=0从第一项1开始)
s
i
n
(
β
)
=
β
−
β
3
3
!
+
β
5
5
!
+
…
…
+
(
−
1
)
i
⋅
β
2
i
+
1
(
2
i
+
1
)
!
sin(β)=\beta-\frac{β^3}{3!}+\frac{β^5}{5!}+……+(-1)^i·\frac{β^{2i+1}}{(2i+1)!}
sin(β)=β−3!β3+5!β5+……+(−1)i⋅(2i+1)!β2i+1(注:i=0从第一项β开始) 则有: 将其入前述x,y的积分公式即可求出x,y的数值公式: 将
r
=
R
L
s
l
r=\frac{RLs}{l}
r=lRLs代入得: 由前式即可在已知R,Ls以及缓和曲线任意点p至起点缓和曲线长l的情况下计算p点的坐标,然后再通过坐标平移转换即可求出该点在全局坐标系下的坐标X,Y;由
β
=
l
2
2
R
L
s
\beta=\frac{l^2}{2RLs}
β=2RLsl2即可求出p点的缓和曲线角,切线方位角也可由缓和曲线起点切线方位角+β的方式计算出来。至此,缓和曲线上任意点的坐标及方位角计算公式已推导完毕。 从该式我们可以看出,如果在x值计算时取i=1,则
x
=
l
−
l
5
40
R
2
L
s
2
x=l-\frac{l^5}{40R^2Ls^2}
x=l−40R2Ls2l5,在y值计算时取i=0,则
y
=
l
3
6
R
L
s
y=\frac{l^3}{6RLs}
y=6RLsl3,此两公式即为我们在教科书以及各种参考书上普遍见到的x,y值计算公式,也就是舍弃了高次项的近似公式表达式。 为了描述方便,我们假定两个参数: 测量施工人员可以根据设计精度要求(在实际施工测量中,若坐标精度达到1mm则完全满足规范和相关控制需求),通过判断前述两个参数xi、yi的值是否小于精度要求值来确定x、y参数计算所需取的项数。曲线内移距p和切线增长值q的计算 : 由几何特性可知,p=y-r(1-cosβ),q=x-rsinβ,将x、y、cosβ、sinβ的级数展开式代入并合并同类项后(推导过程略)得到: 为了防止在计算过程中 数值溢出,也可将前述公式修改为: 当i取值0时,
p
=
L
s
2
24
R
,
q
=
L
s
2
p=\frac{Ls^2}{24R},q=\frac{Ls}{2}
p=24RLs2,q=2Ls也是我们在教科书以及各种参考书上普遍见到的p,q值计算公式,也就是舍弃了高次项的近似公式表达式。 同理,为了描述方便,我们假定pi,qi两个参数,测量施工人员可以根据设计精度要求(在实际施工测量中,若坐标精度达到1mm则完全满足规范和相关控制需求),通过判断下述两个参数pi、qi的值是否小于精度要求值来确定p、q参数计算所需取的项数。
|