卵形曲线计算

1、高等级公路卵形曲线的计算方法周烨摘要在高等级公路施工过程中，常遇到卵形曲线，而设计单位的出发点不同，中线的解算方法也大相径庭。本文着重从卵形中线几种计算方法入手，在此基础之上阐述了卵形曲线的测设。关键词卵形曲线复曲线匝道桥高等级公路卵形曲线是高等级公路、立交桥匝道常见的曲线形式，它由基本的三部分构成：第一圆曲线段、缓和曲线段和第二圆曲线段。中间段缓和曲线用来连接两个不同半径的圆曲线。其中线坐标解算方法有如下几种： 测绘信息网1补全缓和曲线我国公路上采用的缓和曲线为辐射螺旋线，夹在两圆曲线中间的缓和曲线为整个缓和曲线的一部分，缓和曲线上任一点半径与该点至该缓和曲线起点的距离乘积为一定值：r

2、15;la，假设r1r2，可由两圆半径及两圆间的缓和段长ls，求缓和曲线的总长l。llls(1)l就是夹在两圆曲线间缓和段省去的部分，由yh点补长l至o点，以o点为该缓和曲线起点，起点的切线方向为x轴，与之垂直的曲线内侧方向为y轴方向建立坐标系（图1）。缓和曲线公式（推导过程略）如下： 测绘信息网(2)(3)图1利用x、y值可以求得oyh弦与x轴的夹角：3。1为yh点的切线方位角，则ox的方位：1±。o点的坐标可由几何关系求得为（x0，y0）。缓和段上任一点统一坐标可求得：(4)y=yo+xsin±ycos(5)2曲率推算 测绘信息网缓和曲线段曲率半径由第一段圆曲线半径r1

3、变为第二段曲率半径r2（假设r1r2），则缓和曲线曲率半径变化为：(6)其中ls为中间段缓和曲线长，为求缓和曲线方程，现建立以缓和曲线起点为坐标原点，起点的切线方向为x轴，与之垂直的曲线内侧方向为y轴的坐标系（图2），设p点为缓和曲线上任一点，距原点的曲线长为l，该点附近的微分弧长为dl，缓和曲线偏角为，则有dx=dlcos(7)dy=dlsin(8)图2由于 测绘信息网将其代入上式并进行积分可得缓和曲线方程：(9)(10)中间缓和段统一坐标计算为：(11)yyyhxsin±ycos(12)为曲线yh点切线方位。3其它 测绘信息网连接两反曲线或在立交桥匝道上为使墩位美观，常采用缓和曲

4、线连接（图3）。其解算方法以yh点作为起点，以其切线为x轴建立坐标系，不考虑第一种情况中所讲的缓和曲线加长，而直接用式（2）、（3）进行计算，然后统一坐标。图3圆曲线坐标计算在此不再赘述。全站仪广泛地应用到路桥施工中，外业施工放样可在仪器匹配的支持下自动完成。复曲线内业计算、复核线路坐标则成为繁琐问题。4实例 测绘信息网已知r1为5000m，r2为90m，yh点切线方位1为3275659，yh点里程为3279，坐标为（61205.283，101834.119），hy点里程为416.28，坐标为：（61140.068，101892.317）。由已知数据可得：由式（2）、（3）计算填至附表第3、4

5、栏内，由式（4）、（5）公式计算填入第5、6栏内。里程延长至原点距离（m) x(m)y(m)xyk03600057957652398261185452101859215k03279090877761473461205283101834119作者简介：周烨，男，1996年毕业于辽宁工程技术大学测量工程专业，现任铁道部十九局二处助理工程师，先后参加过太峪隧道、南京长江二桥、宁台温高速公路等工程的施工建设。 作者单位：（铁道部十九局） 测绘信息网曲线任意里程中边桩坐标正反算(casio fx-4800p计算器)程序一、程序功能     &#

6、160;  本程序由一个主程序(tyqxjs)和两个子程正算子程序(sub1)、反算子程序(sub2)序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在casio fx-4500p计算器及 casio fx-4850p计算器上运行。二、源程序1.主程序(tyqxjs

7、)"1.sz => xy"："2.xy => sz"：n：u"x0"：v"y0"：o"s0"：g"f0"：h"ls"：p"r0"：r"rn"：q：c=1÷p：d=(p-r)÷(2hpr)：e=180÷：n=1=>goto 1：>goto 2lbl 1：sz：sz：w=abs(s-o)：prog&

8、#160;"sub1"：x"xs"=xy"ys"=ygoto 1lbl 2：xy：xy：i=x：j=y：prog "sub2"：s"s"=o+wz"z"=zgoto 22. 正算子程序(sub1)a=0.1739274226：b=0.3260725774：k=0.0694318442：l=0.3300094782：f=1-l：m=1-k：x=u+w(acos(g+qekw(c+kwd)+bcos(g+qelw(c+lwd)+b

9、cos(g+qefw(c+fwd)+acos(g+qemw(c+mwd)：y=v+w(asin(g+qekw(c+kwd)+bsin(g+qelw(c+lwd)+bsin(g+qefw(c+fwd)+asin(g+qemw(c+mwd)：f=g+qew(c+wd)+90：x=x+zcosf：y=y+zsinf3. 反算子程序(sub2)t=g-90：w=abs(y-v)cost-(x-u)sint)：z=0：lbl 0：prog "sub1"：l=t+qew(c+wd)：z=(j-y)cosl-(i-x)sinl：absz<1e-6=&g

10、t;goto1：>w=w+z：goto 0lbl 1：z=0：prog "sub1"：z=(j-y)÷sinf三、使用说明1、规定        (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，q=-1；当线元往右偏时，q=1；当线元为直线时，q=0。        (2) 当所求点位于中线时，z=0；当位于中线左铡时，z取负

11、值；当位于中线中线右侧时，z取正值。        (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。        (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。         (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以1

12、0的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。         (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明       输入部分：

13、0;     1. sz => xy      2. xy = > sz       n ?  选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算           

14、;   里程和边距。     x0 ？线元起点的x坐标      y0 ？线元起点的y坐标      s0 ？线元起点里程     f0 ？线元起点切线方位角     ls ？线元长度     r0&

15、#160;？线元起点曲率半径     rn ？线元止点曲率半径     q ？  线 元左右偏标志(左偏q=-1，右偏q=1，直线段q=0）      s ？ 正算时所求点的里程      z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)   &#

16、160;  x ？反算时所求点的x坐标      y ？反算时所求点的y坐标       显示部分：       xs=×××    正算时，计算得出的所求点的x坐标       ys=×××

17、;    正算时，计算得出的所求点的y坐标       s=×××    反算时，计算得出的所求点的里程       z=×××    反算时，计算得出的所求点的边距四、算例        某匝道的由

18、五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程s0、起点坐标x0 y0、起点切线方位角f0、线元长度ls、起点曲率半径r0、止点曲率半径rn、线 元左右偏标志q）如下：s0   x0       y0                   f0&#

19、160;           ls r0 rn     q500.000    19942.837    28343.561     125 16 31.00    269.256    &

20、#160;1e45 1e45      0769.256    19787.340    28563.378     125 16 31.00      37.492     1e45 221.75     -

21、1806.748    19766.566    28594.574     120 25 54.07     112.779    221.75 221.75     -1919.527    19736.072    

22、;28701.893       91 17 30.63       80.285    221.75 9579.228   -1999.812    19744.038    28781.659       

23、80 40 50.00     100.000    1e45 1e45       0         1、正算          （注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素） 

24、;           s=700  z=-5    计算得  xs=19831.41785   ys=28509.72590            s=700  z=0     计算得&#

25、160; xs=19827.33592   ys=28506.83837            s=700  z= 5    计算得  xs=19823.25398   ys=28503.95084          &

26、#160; s=780  z=-5    计算得  xs=19785.25749   ys=28575.02270            s=780  z=0     计算得  xs=19781.15561   ys=28572.16358

27、            s=780  z= 5    计算得  xs=19777.05373   ys=28569.30446                   

28、60;   s=870  z=-5    计算得  xs=19747.53609   ys=28654.13091            s=870  z=0     计算得  xs=19742.68648   ys=2

29、8652.91379            s=870  z= 5    计算得  xs=19737.83688   ys=28651.69668                  

30、; s=940  z=-5.123    计算得  xs=19741. 59118  ys=28722.05802            s=940  z=0            计算得  xs=1

31、9736.47687   ys=28722.35642            s=940  z= 3.009    计算得  xs=19733.47298   ys=28722.53168            

32、;       2、 反算             x=19831.418   y=28509.726   计算得  s=699.9999974  z= -5 .00018164       

33、0;     x=19827.336   y=28506.838   计算得  s=699.9996493  z= 0.000145136             x=19823.25398   y=28503.95084    计算得&#

34、160; s=699.9999985  z= 5.000003137           x=19785.25749   y=28575.02270  计算得  s=780.0000035  z= -5 .000001663         

35、0; x=19781.15561   y=28572.16358   计算得  s=780.0000025  z=- 0.000002979           x=19777.05373   y=28569.30446   计算得  s=780.0000016  z=

36、 4.99999578                     x=19747.536   y=28654.131   计算得  s=870.0001137  z= -4.99941049       

37、;    x=19742.686   y=28652.914   计算得  s=870.0003175  z=- 0.00041814           x=19737.837   y=28651.697   计算得  s=870.0002748 

38、 z= 4.999808656          x=19741.5912   y=28722.0580   计算得  s=939.9999786  z= -5.123024937          x=19736.4769   y=287

39、22.3564   计算得  s=939.9999862  z=- 0.000027710          x=19733.4730   y=28722.5317   计算得  s=940.0000238  z= 3.00898694另一个程序功能及原理1.功能说明：本程序由一个主程序(tyqxjs)和四个

40、子程正算子程序(sub1)、反算子程序(sub2)等构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序也可以在casio fx-4800p计算器运行。2计算原理：利用gauss-legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即dca×dcb<0=>该测点在其线元内）进行判断并利用该线

41、元要素反算中桩里程、支距，最后计算出放样数据。二、源程序1.主程序(tyqxjs)lbl 0：defm 50:"1.sz => xy，2.xy => sz"：n： n=1=>goto 1：>goto 2lbl 1：sz：szsz18=>j=1：prog “dat1”：goto 3sz26=>j=2：prog “dat1”：goto 3sz34=>j=3：prog “dat1”：goto

42、 3sz42=>j=4：prog “dat1”：goto 3             sz8（n+1）+2=>j=n：prog “dat1”：goto 3lbl 3：w=abs(s-o)：prog "sub1"："xs=":x=x"ys=":y=y”fwi=”: f=f-m :fdmsgoto 4

43、lbl 2：xy：xy：z4=x：z5 =y：n=0lbl  a：isz n：a=z8n+3-m：b=z8（n+1）+3-m：prog "zx1" ：z6×z7<0=>j=n：prog dat1" ：goto ba=z8n+3+m：b=z8（n+1）+3+m：prog "zx1" ：z6×z7<0=>j=n：prog dat1" ：goto b

44、：> goto albl  b：prog "sub2"："s=":s=o+w"z=":z=zgoto 2lbl 4:j=0:i=pol(x-z1，y-z2):f=j:f<0=>f=f+360  ”dist=”:  i”fw=”:fdms goto 12. 正算子程序(sub1)a=0.1184634425:b=0.2393143352:n=0.2844444444:k=0.046

45、9100770:l=0.2307653449:z3=0.5:x=u+w(acos(g+qekw(c+kwd)+bcos(g+qelw(c+lwd)+ncos(g+qez3w(c+z3wd)+bcos(g+qe(1-l)w(c+(1-l)wd)+acos(g+qe(1-k)w(c+(1-k)wd): y=v+w(asin(g+qekw(c+kwd)+bsin(g+qelw(c+lwd)+nsin(g+qez3w(c+z3wd)+bsin (g+qe(1-l)w(c+(1-l)wd)+asin (g+qe(1-k)w(c+(1-k)wd)：m”ang=”: 

46、f=g+qew(c+wd)+m：x=x+zcosf：y=y+zsinf 3. 反算子程序(sub2)m”ang=”:t=g-m：w=abs(y-v)cost-(x-u)sint)：z=0：lbl 0：prog "sub1"：l=t+qew(c+wd)：z=(z5-y)cosl-(z4-x)sinl：absz<1e-6=>goto1：>w=w+z：goto 0lbl 1：z=0：prog "sub1"：z=(z5-y)÷sinf4. 垂距计算子程序(zx1)z

47、6=(z5-z8n+1)cos a-(z4-z8n)sin az7=(z5-z8(n+1)+1)cos b-(z4-z8(n+1)sin b5.曲线元要素数据库：dat1lbl1:j=1=>u=z8:v=z9:o=z10:g=z11:h=z12:p=z13:r=z14:q=z15j=2=>u=z16:v=z17:o=z18:g=z19:h=z20:p=z21:r=z22:q=z23j=3=> u=z24:v=z25:o=z26:g=z27:h=z28:p=z29:r=z30:q=z31j=4=> u=z32:v

48、=z33:o=z34:g=z35:h=z36:p=z37:r=z38:q=z39                       j=n=> u=z8n:v=z8n+1:o=z8n+2:g=z8n+3:h=z8n+4:p=z8n+5:r=z8n+6:q=z8n+7(注：如有多个曲线元要素继续添加入数据库dat1中)lbl 2 ：u&qu

49、ot;x0"：v"y0"：o"s0"：g"f0"：h"ls"：p"r0"：r"rn"：q：c=1÷p：d=(p-r)÷(2hpr)：e=180÷三、使用说明1、规定        (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，q=-1；当线元往右偏时，q=1；当线元为直线时，q=0。   &

50、#160;    (2) 当所求点位于中线时，z=0；当位于中线左铡时，z取负值；当位于中线中线右侧时，z取正值。        (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。        (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。     &

51、#160;   (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。         (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线

52、的半径。(7)曲线元要素数据库（dat1）可根据线型不同分为各个线元段输入到dat1中，即分为直线段、缓和曲线、圆曲线等。(8)正算时可仅输入里程lp和边距dp及右交角ang全线计算，反算时通过输入测点的x、y坐标和右交角ang后计算器自动判断该点所属曲线元并利用该线元的曲线要素执行反算中桩里程lp及支距dp。2、输入与显示说明（一）、 输入部分：      1. sz => xy      2. xy =&#

53、160;> sz     1、  n ?  选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算              里程和边距。     2、x0 ？线元起点的x坐标      3、y0

54、 ？线元起点的y坐标      4、s0 ？线元起点里程     5、f0 ？线元起点切线方位角     6、ls ？线元长度     7、r0 ？线元起点曲率半径     8、rn ？线元止点曲率半径     9、q 

55、？  线 元左右偏标志(左偏q=-1，右偏q=1，直线段q=0)     10、s ？ 正算时所求点的里程     11、z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)   12、ang？正算边桩时左右边桩连线与线路中线的右交角   13、j？曲线元数据库曲线段判断系数（j=1、2.n）     14、x

56、 ？反算时所求点的x坐标     15、y ？反算时所求点的y坐标16、m?斜交右角17、z1-测站点x坐标18、z2- 测站点y坐标 19、a、b、n是gauss-legendre求积公式中的插值系数20、k 、l、z3 是gauss-legendre求积公式中的求积节点21、n=1时：z8、z9、z10、z11、z12、z13、z14、z15分别是各曲线元的x、y坐标、起始里程s0、初始方位角f0、线元长度ls、起点半径r0、终点半径rn、曲线偏向系数q   &

57、#160;      n=n时：z8n、z8n+1、z8n+2、z8n+3、z8n+4、z8n+5、z8n+6、z8n+7分别是各22、曲线元的x、y坐标、起始里程s0、初始方位角f0、线元长度ls、起点半径r0、终点半径rn、曲线偏向系数q23、正算中的z18、z26、z34、z42z8（n+1）+2分别为各线元中的终点里程24、主程序中的“defm 50”应视具体线元的数量，适当增加扩展变量。25、本次修改版不考虑casio fx-4500pa，因其容量太小，扩展变量实在太少而不能胜任。（二）、显示部分：

58、60;      xs=×××    正算时，计算得出的所求点的x坐标       ys=×××    正算时，计算得出的所求点的y坐标       s=×××    反算时，计算得出的所求

59、点的里程       z=×××    反算时，计算得出的所求点的边距 四、算例        某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程s0、起点坐标x0 y0、起点切线方位角f0、线元长度ls、起点曲率半径r0、止点曲率半径rn、线 元左右偏标志q）如下：   

60、60;s0             x0               y0                  f0  &

61、#160;             ls          r0      rn            q500.000    19942.837&#

62、160;   28343.561     125 16 31.00    269.256     1e45    1e45           0769.256    19787.340   

63、; 28563.378     125 16 31.00      37.492     1e45    221.75        -1806.748    19766.566    28594.574

64、60;    120 25 54.07     112.779    221.75  221.75       -1919.527    19736.072    28701.893      91 17

65、60;30.63      80.285     221.75  9579.228    -1999.812    19744.038    28781.659      80 40 50.00     100.000&

66、#160;    1e45     1e45           0         1、正算          （注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素） 

67、0;          s=700  z=-5    计算得  xs=19831.41785   ys=28509.72590            s=700  z=0     计算得 &

68、#160;xs=19827.33592   ys=28506.83837            s=700  z= 5    计算得  xs=19823.25398   ys=28503.95084           

69、 s=780  z=-5    计算得  xs=19785.25749   ys=28575.02270            s=780  z=0     计算得  xs=19781.15561   ys=28572.16358 

70、;           s=780  z= 5    计算得  xs=19777.05373   ys=28569.30446                    &#

71、160;  s=870  z=-5    计算得  xs=19747.53609   ys=28654.13091            s=870  z=0     计算得  xs=19742.68648   ys=28652.

72、91379            s=870  z= 5    计算得  xs=19737.83688   ys=28651.69668                  

73、0;s=940  z=-5.123    计算得  xs=19741. 59118  ys=28722.05802            s=940  z=0            计算得  xs=19736.47687   ys=28722.35642            s=940  z= 3.009    计算得  xs=19733.47298   ys=28722.53168