卵形曲线计算

高等级公路卵形曲线的计算方法周烨摘要在高等级公路施工过程中，常遇到卵形曲线，而设计单位的出发点不同，中线的解算方法也大相径庭。本文着重从卵形中线几种计算方法入手，在此基础之上阐述了卵形曲线的测设。关键词卵形曲线复曲线匝道桥高等级公路卵形曲线是高等级公路、立交桥匝道常见的曲线形式，它由基本的三部分构成：第一圆曲线段、缓和曲线段和第二圆曲线段。中间段缓和曲线用来连接两个不同半径的圆曲线。其中线坐标解算方法有如下几种：tp/w■mI■信息1补全缓和曲线我国公路上采用的缓和曲线为辐射螺旋线，夹在两圆曲线中间的缓和曲线为整个缓和曲线的一部分，缓和曲线上任一点半径与该点至该缓和曲线起点的距离乘积为一定值：RXL=A，假设R>R，可由两圆半径及两圆间的缓和段长ls，求缓和曲线的总长L。12Al=L-l（1）△l就是夹在两圆曲线间缓和段省去的部分，由YH点补长Al至。点，以。点为该缓和曲线起点，起点的切线方向为x轴，与之垂直的曲线内侧方向为y轴方向建立坐标系（图1）。缓和曲线公式（推导过程略）下：h.wwwil ■绘10R-E；- 3伯观1rr"'iiRL-s部"E十匝而iFDothermap・comothermap・com图1利用X、y值可以求得。一YH弦与x轴的夹角：B=36°a为YH点的切线方位角，则ox的方位：a=ai±3oo点的坐标可由几何关系求得为（X。，yo）。缓和段上任一点统一坐标可求得：y=yo+xsina土ycosa （5）2曲率推算缓和曲线段曲率半径由第一段圆曲线半径％变为第二段曲率半径r2（假设R1>R2），则缓和曲线曲率半径变化为：1 2其中l为中间段缓和曲线长，为求缓和曲线方程，现建立以缓和曲线起点为巫标原点，起点的切线方向为X轴，与之垂直的曲线内侧方向为y轴的坐标系（图2），设P点为缓和曲线上任一点，距原点的曲线长为1，该点附近的微分孤长为dl，缓和曲线偏角为B，则有dx二dlcosB （7）dy二dlsinB （8）图2由于/-将其代入上式并进行积分可得缓和曲线方程：=rm—R顶'一一E济如⑼1山皿・尾）'，I"（10）中间缓和段统一坐标计算为：（11）Y=yYJxsina+ycosa （12）a为曲线YH点切线方位。3其它连接两反曲线或在立交桥匝道上为使墩位美观，常采用缓和曲线连接（图3）。其解算方法以YH点作为起点，以其切线为x轴建立坐标系，不考虑第一种情况中所讲的缓和曲线加长，而直接用式（2）、（3）进行计算，然后统一坐标。圆曲线坐标计算在此不再赘述。全站仪广泛地应用到路桥施工中，外业施工放样可在仪器匹配的支持下自动完成。复曲线内业计算、复核线路坐标则成为繁琐问题。4实例已知R为5000m,R?为90m,YH点切线方位气为327—56—59,YH点里程为+1327.9,坐标为（61205.283,101834.119）,HY点里程为+416.28，坐标为:（61140.068,101892.317）。由已知数据可得：电0+*—心,少一"站、3七—卅一部二&_].N山式⑵哥）得；膈="我iy=0,000090-artL^II"Ax.ox矽仙'.’：'=方IJ=S若一苗一*由式（2）、（3）计算填至附表第3、4栏内，由式（4）、（5）公式计算填入第5、6栏内。里程延长至原点距离（m）△x(m)△y(m)xy、+360.0057.957.6523.98261185.452101859.215K0+327.909087.77614.73461205.283101834.119作者简介：周烨，男，1996年毕业于辽宁工程技术大学测量工程专业，现任铁道部十九局二处助理工程师，先后参加过太峪隧道、南京长江二桥、宁台温高速公路等工程的施工建设。作者单位：（铁道部十九局）曲线任意里程中边桩坐标正反算(CASIOfx-4800P计算器)程序一、 程序功能本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段一一直线、圆曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径)及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIOfx-4500P计算器及CASIOfx-4850P计算器上运行。二、 源程序主程序(TYQXjs)"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)H2HPR)：E=180；n：N=1=>Goto1：手〉Goto2A—」Lbl1：{SZ}：SZ：W=Abs(S-O)：Prog"SUB1"：X"XS"=X，Y"YS"=Y，Goto1」Lbl2：{XY}：XY：I=X：J=Y：Prog"SUB2"：S"S"=O+W，Z"Z"=Z，Goto2正算子程序(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))+Bcos(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：Lbl0：Prog"SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：手>W=W+Z：Goto0A^JLbl1：Z=0：Prog"SUB1"：Z=(J-Y)?sinF三、使用说明1、 规定以道路中线的前进方向(即里程增大的方向)区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、 输入与显示说明输入部分：SZ=>XYXY=>SZN?选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。X0？线元起点的X坐标Y0？线元起点的Y坐标S0？线元起点里程F0？线元起点切线方位角LS？线元长度R0？线元起点曲率半径RN？线元止点曲率半径Q？线元左右偏标志（左偏Q=-1，右偏Q=1，直线段Q=0）S？正算时所求点的里程Z？正算时所求点距中线的边距（左侧取负，值右侧取正值，在中线上取零）X？反算时所求点的X坐标Y？反算时所求点的Y坐标显示部分：XS=xxx正算时，计算得出的所求点的X坐标YS=xxx正算时，计算得出的所求点的Y坐标S=xxx反算时，计算得出的所求点的里程Z=xxx反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0Y0、起点切线方位角F0、线元长度

要素）/一、d.5（注意:略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元S=700Z=-5计算得XS=19831.41785YS=28509.72590S=700Z=0计算得XS=19827.33592YS=28506.83837S=700Z=5计算得XS=19823.25398YS=28503.95084S=780Z=-5计算得XS=19785.25749YS=28575.02270S=780Z=0计算得XS=19781.15561YS=28572.16358S=780Z=5计算得XS=19777.05373YS=28569.30446S=870Z=-5计算得XS=19747.53609YS=28654.13091S=870Z=0计算得XS=19742.68648YS=28652.91379S=870Z=5计算得XS=19737.83688YS=28651.69668S=940Z=-5.123计算得XS=19741.59118YS=28722.05802S=940Z=0计算得XS=19736.47687YS=28722.35642S=940Z=3.009计算得XS=19733.47298YS=28722.53168LS、起点曲率半径R0、止点曲率半径要素）/一、d.5（注意:略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元S=700Z=-5计算得XS=19831.41785YS=28509.72590S=700Z=0计算得XS=19827.33592YS=28506.83837S=700Z=5计算得XS=19823.25398YS=28503.95084S=780Z=-5计算得XS=19785.25749YS=28575.02270S=780Z=0计算得XS=19781.15561YS=28572.16358S=780Z=5计算得XS=19777.05373YS=28569.30446S=870Z=-5计算得XS=19747.53609YS=28654.13091S=870Z=0计算得XS=19742.68648YS=28652.91379S=870Z=5计算得XS=19737.83688YS=28651.69668S=940Z=-5.123计算得XS=19741.59118YS=28722.05802S=940Z=0计算得XS=19736.47687YS=28722.35642S=940Z=3.009计算得XS=19733.47298YS=28722.53168S0X0Y0F0LSR0RNQ500.00019942.83728343.5611251631.00269.2561E451E45 0769.25619787.34028563.3781251631.00 37.4921E45221.75-1806.74819766.56628594.5741202554.07112.779221.75221.75-1919.52719736.07228701.893911730.63 80.285221.759579.228-1999.81219744.03828781.659804050.00100.0001E451E45 0正算1、2、反算X=19831.418Y=28509.726计算得S=699.9999974Z=-5.00018164X=19827.336Y=28506.838计算得S=699.9996493Z=0.000145136X=19823.25398Y=28503.95084计算得S=699.9999985Z=5.000003137X=19785.25749Y=28575.02270计算得S=780.0000035Z=-5.000001663X=19781.15561Y=28572.16358计算得S=780.0000025Z=-0.000002979X=19777.05373Y=28569.30446计算得S=780.0000016Z=4.99999578X=19747.536Y=28654.131计算得S=870.0001137Z=-4.99941049X=19742.686Y=28652.914计算得S=870.0003175Z=-0.00041814X=19737.837Y=28651.697计算得S=870.0002748Z=4.999808656X=19741.5912Y=28722.0580计算得S=939.9999786Z=-5.123024937X=19736.4769Y=28722.3564计算得S=939.9999862Z=-0.000027710X=19733.4730Y=28722.5317计算得S=940.0000238Z=3.00898694另一个程序功能及原理功能说明：本程序由一个主程序(TYQXjs)和四个子程一一正算子程序(SUB1)、反算子程序(SUB2)等构成，可以根据曲线段——直线、圆曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径)及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序也可以在CASIOfx-4800P计算器运行。计算原理：利用Gauss-Legendre5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异(即DcaxDcbv0=＞该测点在其线元内)进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。二、源程序

1.主程序(TYQXjs)1.主程序(TYQXjs)Lbl0：Defm50:"1.SZ=>XY,Lbl1：{SZ}：SZS<Z[18]=>J=1：S<Z[26]=>J=2：S<Z[34]=>J=3：S<Z[42]=>J=4：2.XY=>SZ"：N：N=1=>Goto1：手〉Goto2A—Prog"DATl”：Prog"DAH”：Prog"DAH”：Prog"DATl”：Goto3AGoto3AGoto3AGoto3AS<Z[8(N+1)+2]=>J=n：Prog"DATl”：Goto3ALbl3：W=Abs(S-O)：Prog"SUB1"："XS=":X=X,"YS=":Y=Y,”FWI=”:F=F-M:F—DMS，Goto41Lbl2：{XY}：XY：Z[4]=X：Z[5]=Y：N=0LblA：IszN：A=Z[8N+3]-M：B=Z[8(N+1)+3]-M：Prog"ZX1":Z[6]xZ[7]v0=>J=N：ProgDAT1":GotoBAA=Z[8N+3]+M：B=Z[8(N+1)+3]+M：Prog"ZX1":Z[6]xZ[7]v0=>J=N：ProgDAT1":GotoB：手〉GotoAALblB：Prog"SUB2"："S=":S=O+W，"Z=":Z=Z』Goto2Lbl4:J=0:I=Pol(X-Z[1],Y-Z[2]):F=J:F<0=>F=F+360A”DIST=”：I，”FW='':F—DMS，Goto12.正算子程序(SUB1)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:Z[3]=0.5:X=U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(G+QEZ[3]W(C+Z[3]WD))+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+QEZ[3]W(C+Z[3]WD))+Bsin(G+QE(1-L)W(C+(1-L)WD))+Asin(G+QE(1-K)W(C+(1-K)WD)))：M”ANG=”:F=G+QEW(C+WD)+M：X=X+ZcosF：Y=Y+ZsinF3,反算子程序(SUB2)M”ANG=”:T=G-M：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl0：Prog"SUB1"：L=T+QEW(C+WD)：Z=(Z[5]-Y)cosL-(Z[4]-X)sinL：AbsZv1E-6=>Goto1：手>W=W+Z：Goto0A—」Lbl1：Z=0：Prog"SUB1"：Z=(Z[5]-Y)《sinF垂距计算子程序(ZX1)Z[6]=(Z[5]-Z[8N+1])COSA-(Z[4]-Z[8N])SINAZ[7]=(Z[5]-Z[8(N+1)+1])COSB-(Z[4]-Z[8(N+1)])SINB曲线元要素数据库：DAT1Lbl1:J=1=>U=Z[8]:V=Z[9]:O=Z[10]:G=Z[11]:H=Z[12]:P=Z[13]:R=Z[14]:Q=Z[15]AJ=2=>U=Z[16]:V=Z[17]:O=Z[18]:G=Z[19]:H=Z[20]:P=Z[21]:R=Z[22]:Q=Z[23]AJ=3=>U=Z[24]:V=Z[25]:O=Z[26]:G=Z[27]:H=Z[28]:P=Z[29]:R=Z[30]:Q=Z[31]AJ=4=>U=Z[32]:V=Z[33]:O=Z[34]:G=Z[35]:H=Z[36]:P=Z[37]:R=Z[38]:Q=Z[39]AJ=N=>U=Z[8N]:V=Z[8N+1]:O=Z[8N+2]:G=Z[8N+3]:H=Z[8N+4]:P=Z[8N+5]:R=Z[8N+6]:Q=Z[8N+7]A(注：如有多个曲线元要素继续添加入数据库DAT1中)Lbl2: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、规定以道路中线的前进方向(即里程增大的方向)区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。曲线元要素数据库(DAT1)可根据线型不同分为各个线元段输入到DAT1中，即分为直线段、缓和曲线、圆曲线等。正算时可仅输入里程Lp和边距Dp及右交角ANG全线计算，反算时通过输入测点的X、Y坐标和右交角ANG后计算器自动判断该点所属曲线元并利用该线元的曲线要素执行反算中桩里程Lp及支距Dp。2、输入与显示说明(一)、输入部分：SZ=>XYXY=>SZ1、 N?选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。2、 X0？线元起点的X坐标3、 Y0？线元起点的Y坐标4、 S0？线元起点里程5、 F0？线元起点切线方位角6、 LS？线元长度7、 R0？线元起点曲率半径8、 RN？线元止点曲率半径9、 Q？线元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0)10、 S？正算时所求点的里程11、Z?正算时所求点距中线的边距（左侧取负，值右侧取正值，在中线上取零）12、 ANG？正算边桩时左右边桩连线与线路中线的右交角13、 J？曲线元数据库曲线段判断系数（J=1、2.....n）14、 X？反算时所求点的X坐标15、 Y？反算时所求点的Y坐标16、 M?斜交右角17、 Z[1]-----测站点X坐标18、 Z⑵-----测站点Y坐标19、 A、B、N是Gauss-Legendre求积公式中的插值系数20、 K、L、Z[3]是Gauss-Legendre求积公式中的求积节点21、 N=1时：Z[8]、Z[9]、Z[10]、Z[11]、Z[12]、Z[13]、Z[14]、Z[15]分别是各曲线元的X、Y坐标、起始里程S0、初始方位角F0、线元长度Ls、起点半径R0、终点半径Rn、曲线偏向系数QN=n时：Z[8N]、Z[8N+1]、Z[8N+2]、Z[8N+3]、Z[8N+4]、Z[8N+5]、Z[8N+6]、Z[8N+7]分别是各22、曲线元的X、Y坐标、起始里程S0、初始方位角F0、线元长度Ls、起点半径R0、终点半径Rn、曲线偏向系数Q23、 正算中的Z[18]、Z[26]、Z[34]、Z[42]......Z[8（N+1）+2]分别为各线元中的终点里程24、 主程序中的“Defm50”应视具体线元的数量，适当增加扩展变量。25、 本次修改版不考虑CASIOfx-4500PA，因其容量太小，扩展变量实在太少而不能胜任。（二）、显示部分：XS=xxx正算时，计算得出的所求点的X坐标YS=xxx正算时，计算得出的所求点的Y坐标

S=xxx反算时，计算得出的所求点的里程Z=xxx反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：S0X0Y0F0LSR0RNQ500.00019942.83728343.5611251631.00269.2561E451E450769.25619787.34028563.3781251631.0037.4921E45221.75-1806.74819766.56628594.5741202554.07112.779221.75221.75-1919.52719736.07228701.893911730.6380.285221.759579.228-1999.8121、19744.038正算28781.659804050.00100.0001E451E450（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=19831.41785YS=28509.72590S=700Z=0计算得XS=19827.33592YS=28506.83837S=700Z=5计算得XS=19823.25398YS=28503.95084S=780Z=-5计算得XS