平曲线视距横净距的计算

平曲线视距横净距的计算发表人：王乃坤江树华单位：龙建路桥股份有限公司第二工程处日期：二OO四年十二月十五三十日平曲线视距横净距的计算王乃坤江树华（龙建路桥股份有限公司第二工程处)提要：本文介绍用计算机计算平曲线及相邻直线段上任一点的横净距的方法。

关键词：平曲线横净距计算机计算CalculatingabouttheCrossClearanceDistanceofHorizontalCurveStadiaWANGNai-kunJIANGShu-huaQUZhi-chengAbstract：Calculatingmethodofhorizontalcurveandcrossclearancedistancewithcomputerispresented.Keywords：HorizontalcurveCrossclearancedistanceComputer1前言如何准确计算平曲线及相邻直线段上任一点的横净距，是我们工程技术人员在实际工作中常遇到的问题。

近期我们成功地利用计算机程序解决了带缓和曲线的平曲线横净距计算，省时省力，起到了事半功倍的效果。

现介绍如下，仅供参考。

2横净距的计算方法2.1计算原理如图1所示，某交点转角为α，平曲线半径为R1，缓和曲线长为Ls1(我们将圆曲线作为Ls1＝0的特例处理)。

若行车道宽度为b，则计算横净距时的行车轨迹线（距未加宽时的行车道内侧边缘1.5m，图中虚曲线所示）与路中线的径向间距△R＝b／2－1.5。

M为平曲线和相邻直线段上的任一点，M所在断面的横净距可按下法计算：在M点的法线MN两侧的行车轨迹线上分别找一点A、B，使A、B两点间沿行车轨迹线的长度等于设计视距S，计算AB连线与MN的交点E到M点的距离值H；保持A、B两点间沿行车轨迹线的长度不变，使A、B两点沿行车轨迹线同步移动时，H 值也随之改变，最大的H值与△R之差即为M点的横净距。

2.2行车轨迹线参数的确定对于半径为R1、缓和曲线长度为L S1的平曲线，行车轨迹线上对应的曲线半径R2＝R1－△R，对应的缓和曲线长度L S２应是多少呢?我们分别用p、q、E、T表示圆曲线设置缓和曲线后的内移值、切垂距、外距、切线长，并分别用下标1、2对路中线、行车轨迹线上有关的量进行的区分。

如图2所示，因路中线与行车轨迹线的径向间距为△R，则应有：图二将有关各量代入并整理后可得:在R1R2L S1是定值的情况下，使△p、△q同时等于零是不可能的(L S1＝0时除外)。

也就是说，无论L S２取何值，由R2L S２确定的轨迹线(下文中称之为计算行车轨迹线)与理论行车轨迹线(距未加宽时的行车道内侧边缘1.5m)不可能完全重合，两轨迹线间的偏离程度可表示为：将有关各量代入并整理后可得：上式中，△E表示计算行车轨迹线与理论行车轨迹线的曲线中点间的距离，△T表示两轨迹线上的曲线起(终)点沿切线方向的距离。

既然计算行车轨迹线与理论行车轨迹线不能完全重合，那么我们能否选取适当的L S２值，使两轨迹线尽量接近，从而满足横净距计算的精度要求呢?由(1)、(2)式可以看出，△p对两轨迹线间偏离程度的影响远比△q的影响大，如果我们通过L S２的取值使△p等于零，即用牛顿迭代法按下式解算L S２此外，两迹线的偏离程度可表示为：△E＝0△T＝△q两轨迹线在圆曲线段重合，下面我们来看看两轨迹线在缓和曲线段的偏离情况。

? 如图3所示，A为路中线上缓和曲线段内的任一点，距缓和曲线起点ZH1的长度为L1，在xoy 坐标系中的坐标为(X A，Y A)，A点的法线AF与x轴的夹角为：沿法线AF由A点向曲线内侧移动△R得B点，则B点的坐标为：由路中线与理论行车轨迹线的关系知：B点在理论行车轨迹线上。

令C点为计算行车轨迹线上缓和曲线段内一点，距缓和曲线段起点ZH２的长度为L２，在x'o'y'坐标系中的坐标为(x'c，y'c)，经坐标轴平移后可得C点在xoy坐标系中的坐标：若X c＝X b，则(Y c－Y b)cos(γ－π/2)即可视为该位置处计算行车轨迹线与理论行车轨迹线间的径向偏离值。

按上述方法，我们用计算机对《公路曲线测设用表》(第二册)中“缓和曲线长度总表”列出的R1、Ls1的所有组合进行了有关计算。

计算结果表明：计算行车轨迹线与理论行车轨迹线在缓和曲线段的径向偏离值有以下特点：（1)两轨迹线在缓和曲线段中部的偏离值最大，向两端逐渐趋于零；（2)Ls1值不变时，R1值愈小，两轨迹线在缓和曲线段的偏离值就愈大；（3)R1值不变时，Ls1愈大，两轨迹线在缓和曲线段的偏离值就愈大；（4)表列各种组合中，两轨迹线在缓和曲线段的最大偏离值一般均小于0.01m，最大偏离值超过0.01m的组合只有以下4个：①R1＝65、Ls1＝60、b＝7.5时，最大偏离值为0.015m；②R1＝60、Ls1＝60、b＝7.5时，最大偏离值为0.018m；③R1＝35、Ls1＝35、b＝6时，最大偏离值为0.012m；④R1＝30、Ls1＝35、b＝6时，最大偏离值为0.016m。

通过对两轨迹线间偏离值的计算可知：由R２＝R1－△R和按(3)式解算的Ls２确定的计算行车轨迹线与理论行车轨迹线非常接近(Ls1＝0时两轨迹线重合)，若用计算行车轨迹线代替理论行车轨迹线来计算横净距，最大误差不会超过0.02m，完全能够满足横净距计算的精度要求。

2.3横净距的计算为方便计算，我们也用桩号表示行车轨迹线上各点间的相对位置关系，并假定行车轨迹线上曲线起点ZH2(ZY2)与路中线上曲线起点ZH1(ZY1)的桩号相同(两线上的桩号是相互独立的)。

这样对于任一给定的桩号，能且只能在行车轨迹线上找到一点与之对应。

如图1，欲求横净距的断面M的桩号为L0，我们不妨先假设行车轨迹线上视距起终点A、B的桩号分别为L0－S／2、L0＋S／2(在应用范围内，A、B不会落在法线MN的同一则)，此时沿行车轨迹线A、B间的长度等于设计视距S，计算AB与MN的交点到M点的距离值H；然后将A、B两点沿行车轨迹线向前移动0.01m，即A、B点的桩号同时增大0.01m，再计算MN与新直线AB的交点到M点的距离值H。

根据两次计算的距离值H的大小，使A、B向距离值H增大的方向同步等量逐渐移动(以保证沿行车轨迹线A、B间的长度恒等于设计视距S)直到距离值H开始减小为止，此时的距离值H与△R之差即为M点的横净距。

按上述方法，我们编写了计算横净距的程序。

3实例应用鹤大公路鸡牡段林口支线为山岭区三级公路，行车道宽6m，设计视距(二倍停车视距)S＝60m；JD2转角α＝76°32′，曲线半径100m，缓和曲线长Ls＝25m，曲线起点桩号为K2＋838.68，计算各有关断面的横净距。

程序运行时的操作及计算结果如下(带下画线者为输入的数据)：曲线起点桩号?838.68转角：?76°?32′?0〃平曲线半径?100缓和曲线长?25行车道宽度?6设计视距?60曲线主点桩号：ZH＝838.68HY＝863.68QZ＝917.97YH＝972.26HZ＝997.26断面桩号(＜0结束)?800横净距：0.00断面桩号(＜0结束)?820横净距：0.34断面桩号(＜0结束)?838.68横净距：1.41断面桩号(＜0结束)?840横净距：1.53断面桩号(＜0结束)?860横净距：3.52断面桩号(＜0结束)?863.68横净距：3.80断面桩号(＜0结束)?880横净距：4.45断面桩号(＜0结束)?900横净距：4.53断面桩号(＜0结束)?917.97横净距：4.53断面桩号(＜0结束)?920横净距：4.53断面桩号(＜0结束)?940横净距：4.53断面桩号(＜0结束)?960横净距：4.37断面桩号(＜0结束)?980横净距：3.14断面桩号(＜0结束)?997.26横净距：1.41断面桩号(＜0结束)?1000横净距：1.19断面桩号(＜0结束)?1020横净距：0.23断面桩号(＜0结束)?1040横净距：0.00断面桩号(＜0结束)?－1OK4源程序10CLEAR：DEFDBLA－H，20CLS：LOCATE11，1：INPUT〃曲线起点桩号：〃；ZH130PRINT〃转角：40LOCATE12，6：INPUTZJ1：LOCATE12，14：INPUTZJ2：LOCATE12，22：INPUTZJ3 50INPUT〃平曲线半径＝〃；R1：INPUT〃缓和曲线长＝〃；LS160INPUT〃平曲线半径＝〃；LMK：INPUT〃设计视距＝〃；S70ZJ＝(ZJ1＋ZJ2／60＋ZJ3／3600)*PI／180：DR＝LMK／2－1.5：R2＝R1－DR80IFLS1＞0THENGOSUB640：REM计算LS290R＝R1：LS＝LS1：GOSUB710：LY1＝LY：P1＝P：Q1＝Q：T1＝T100R＝R2：LS＝LS2：GOSUB710：LY2＝LY：P2＝P：Q2＝Q：T2＝T110DT＝Q1－Q2：ZH2＝ZH1：SS＝S／2／R2*R1120HY1＝ZH1＋LS1：YH1＝HY1＋LY1：HZ1＝YH1＋LS1：QZ1＝(ZH1＋HZ1)／2130HY2＝ZH2＋LS2：YH1＝HY2＋LY2：HZ2＝YH2＋LS2：QZ2＝(ZH2＋HZ2)／2140PRINT〃曲线主点桩号：〃：A＄＝〃######.##〃：IFLS1＝0THEN180150PRINT〃ZH＝〃；USINGA＄；ZH1：PRINT〃HY＝〃；USINGA＄；HY1160PRINT〃QZ＝〃；USINGA＄；QZ1：PRINT〃YH＝〃；USINGA＄；YH1170PRINT〃HZ＝〃；USINGA＄；HZ1：GOTO200180PRINT〃ZY＝〃；USINGA＄；ZH1：PRINT〃QZ＝〃；USINGA＄；QZ1190PRINT〃YZ＝〃；UXINGA＄；HZ1200INPUT〃断面桩号(＜0结束)〃；L0：IFL0＜0THEN450210IFL0＜＝ZH1－SORL0＞＝HZ1＋STHENH＝0：GOTO440220IFL0＞QZ1THENL0＝ZH1＋(HZ1－L0)：REM若欲求横净距断面在曲线后半部，则计算曲线前半部对称位置的横净距230IFL0－SS＞＝HY1ANDL0＋SS＜＝YH1THENH＝R2*(1－COS(S／2／R2))：GOTO440240R＝R1：LS＝LS1：P＝P1：Q＝Q1：X0＝T1＋T1*COS(ZJ)：Y0＝T1*SIN(ZJ)：IFABS(ZJ－PI)＜1E－10THENX0＝0Y0＝2*(R1＋P1)250IFL0＜＝ZH1THENA0＝PI／2：X＝L0－ZH1：Y＝0：GOTO300260IFL0＞＝HZ1THENA0＝ZJ＋PI／2：X＝HZ1－L0：Y＝0：GOSUB620：GOTO300270IFL0＞YH1THENL＝HZ1－L0：A0＝ZJ＋PI／2－L^2／2／R1／LS1：GOSUB600：GOSUB620：GOTO300280IFL0＞＝HY1THENL＝L0－HY1＋LS1／2：A0＝L／R1＋PI／2：GOSUB590：X＝X＋Q：Y＝Y＋P：GOTO300290L＝L0－ZH1：A0＝L^2／2／R1／LS1＋PI／2：GOSUB600300XM＝X－DT：YM＝Y－DR310IFA0＝PI／2THENA1＝1：B1＝0：C1＝－XMELSEK＝TAN(A0)：A1＝K：B1＝－1：C1＝YM－K*XM320R＝R2：LS＝LS2：P＝P2：Q＝Q2：X0＝T2＋T2*COS(ZJ)：Y0＝T2*SIN(ZJ)：IFABS(ZJ－PI)＜1E－10THENX0＝0：Y0＝2*(R2＋P2)330A＝L0－S／2：B＝L0＋S／2：GOSUB460：H1＝H：A＝A＋.01：B＝B＋.01：GOSUB460：A＝A－.01：B＝B－.01340REM分别以10m、1m、0.1m、0.01m的步距使A、B两点向H值增大的方向移动，以减少计算次数350D＝100：IFH＜H1THENF＝－1ELSEF＝1360IFH＝H1ANDB＜＝ZH2THENF＝1370IFH＝H1ANDA＞HZ2THENF＝－1380FORI＝1TO3：D＝D／10：D1＝D*F390FORJ＝1TO100：A＝A＋D1：B＝B＋D1：GOSUB460400IFH＜＝H1THENA＝A－D1：B＝B－D1：GOTO420410H1＝H：NEXTJ420NEXTI430H＝H1－DR440U＝CSRLIN：LOCATEU－1，48：PRINT〃横净距：〃；USING〃###.##〃；H：GOTO200450END 460REM计算视线AB与法线MN的交点E至M点的距离H值470G＝A：GOSUB520：XA＝X：YA＝Y480G＝B：BOSUB520：XB＝X：YB＝Y490A2＝YB－YA：B2＝XA－XB：C2＝(XB－XA)*YA－(YB－YA)*XA500AB＝A1*B2－A2*A2*B2：XE＝(B1*C2－B2*C1)／AB：YE＝－(A1*C2－A2*C1)／AB510H＝SQR((XE－XM)^2＋(YE－YM)^2)：RETURN520REM计算行车轨迹线已知桩号的任一点的坐标530IFG＜＝ZH2THENX＝G－ZH2：Y＝0：RETURN540IFG＞＝HZ2THENX＝HZ2－G：Y＝0：GOSUB620：RETURN550IFLS2＝0THENL＝G－ZH2：GOSUB590：RETURN560IFG＞YH2THENL＝HZ2－G：GOSUB600：GOSUB620：RETURN570IFG＞HY2THENL＝G－HY2＋LS2／2：GOSUB590：X＝X＋Q：Y＝Y＋P：RETURN580L＝G－ZH2：GOSUB600：RETURN590X＝R*SIN(L／R)：Y＝R*(1－COS(L／R))：RETURN600C＝R*LS：X＝L－L^5／40／C^2＋L^9／3456／C^4610Y＝L^3／6／C－L^7／336／C^3＋L^11／42240!／C^5：RETURN620XI＝－X*COS(ZJ)－Y*SIN(ZJ)：YI＝Y*COS(ZJ)－X*SIN(ZJ)630X＝XI＋X0：Y＝Y1＋Y0：TETURN640REM用牛顿迭代法计算LS2650LS2＝LS1660FORI＝1TO10670FX＝LS2^2／24／R2－LS2^4／2688／R2^3＋LS2^6／506880!／R2^5－LS1^2／24／R1＋LS1^4／2688／R1^3－LS1^6／506880!／R1^5680DX＝LS2／12／R2－LS2^3／672／R2^3＋LS2^5／84480!／R2^5690LS2＝LS2－FX／DX：PRINTLS2：NEXTI700RETURN710BT＝LS／2／R：LY＝R*ZJ－LS：IFLS＝0THEN730720L＝LS：GOSUB600：Q＝X－R*SIN(BT)：P＝Y－R*(1－COS(BT))730IFZJ＝PITHENRETURN740T＝(R＋P)*TAN(ZJ／2)＋Q：RETURN。