平曲线视距横净距的计算

平曲线视距横净距的计算蒋玲玲王平安张泽桂(河南省交通科研所)(河南省交通学校)(河南省交通工程局)提要本文介绍用计算机计算平曲线及相邻直线段上任一点的横净距的方法。关键词平曲线横净距计算机计算(HenanCommunicationsResearchInstitute)AbstractCalculatingmethodofhorizontalcurveandcrossclearanceKeywordsHorizontalcurveCrossclearancedistanceComputer1前言用计算机计算圆曲线及相邻直线段上任一点的横净距的方法，我们曾在本刊1997年第4期发表的《圆曲线视距横净距的计算》一文中介绍过。近期我们又成功地解决了带缓和曲线的平曲线的横净距计算。现介绍如下，供参考。2横净距的计算方法2.1计算原理如图1所示，某交点转角为α，平曲线半径为，缓和曲线长为我们将圆曲线作为＝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、缓和曲线长度为LS1的平曲线，行车轨迹线上对应的曲线半径R2＝R1－△R，对应的缓和曲线长度LS２应是多少呢我们分别用p、q、E、T表示圆曲线设置缓和曲线后的内移值、切垂距、外距、切线长，并分别用下标1、2对路中线、行车轨迹线上有关的量进行的区分。如图2所示，因路中线与行车轨迹线的径向间距为△R，则应有：图二将有关各量代入并整理后可得:在R1R2LS1是定值的情况下，使△p、△q同时等于零是不可能的(LS1＝0时除外)。也就是说，无论LS２取何值，由R2LS２确定的轨迹线(下文中称之为计算行车轨迹线)与理论行车轨迹线(距未加宽时的行车道内侧边缘1.5m)不可能完全重合，两轨迹线间的偏离程度可表示为：将有关各量代入并整理后可得：上式中，△E表示计算行车轨迹线与理论行车轨迹线的曲线中点间的距离，△T表示两轨迹线上的曲线起(终)点沿切线方向的距离。既然计算行车轨迹线与理论行车轨迹线不能完全重合，那么我们能否选取适当的LS２值，使两轨迹线尽量接近，从而满足横净距计算的精度要求呢?由(1)、(2)式可以看出，△p对两轨迹线间偏离程度的影响远比△q的影响大，如果我们通过LS２的取值使△p等于零，即用牛顿迭代法按下式解算LS２此外，两迹线的偏离程度可表示为：△E＝△T＝两轨迹线在圆曲线段重合，下面我们来看看两轨迹线在缓和曲线段的偏离情况。如图3所示，A为路中线上缓和曲线段内的任一点，距缓和曲线起点ZH1的长度为L1，在xoy坐标系中的坐标为(XA，YA)，A点的法线AF与x轴的夹角为：沿法线AF由A点向曲线内侧移动△R得B点，则B点的坐标为：由路中线与理论行车轨迹线的关系知：B点在理论行车轨迹线上。令C点为计算行车轨迹线上缓和曲线段内一点，距缓和曲线段起点ZH２的长度为L２，在x'o'y'坐标系中的坐标为(x'c，y'c)，经坐标轴平移后可得C点在xoy坐标系中的坐标：若Xc＝Xb，则(Yc－Yb)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横净距的计算为方便计算，我们也用桩号表示行车轨迹线上各点间的相对位置关系，并假定行车轨迹线上曲线起点与路中线上曲线起点的桩号相同(两线上的桩号是相互独立的)。这样对于任一给定的桩号，能且只能在行车轨迹线上找到一点与之对应。如图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；转角α＝76°32′，曲线半径100m，缓和曲线长Ls＝25m，曲线起点桩号为k2＋838.68，计算各有关断面的横净距。程序运行时的操作及计算结果如下(带下画线者为输入的数据)：曲线起点桩号转角：〃平曲线半径缓和曲线长行车道宽度设计视距曲线主点桩号：ZH＝HY＝QZ＝YH＝HZ＝断面桩号(＜0结束横净距：断面桩号(＜0结束)?820横净距：断面桩号(＜0结束横净距：断面桩号(＜0结束)?840横净距：断面桩号(＜0结束)?860横净距：断面桩号(＜0结束横净距：断面桩号(＜0结束)?880横净距：断面桩号(＜0结束)?900横净距：断面桩号(＜0结束横净距：断面桩号(＜0结束)?920横净距：断面桩号(＜0结束横净距：断面桩号(＜0结束)?960横净距：断面桩号(＜0结束)?980横净距：断面桩号(＜0结束横净距：断面桩号(＜0结束)?1000横净距：断面桩号(＜0结束)?1020横净距：断面桩号(＜0结束)?1040横净距：断面桩号(＜0结束)?－4源程序10CLEAR：DEFDBLA－H，K－Z：PI＝20CLS：LOCATE11，1：INPUT〃曲线起点桩号：〃；30PRINT〃转角：40LOCATE12，6：INPUTZJ1：LOCATE12，14：INPUTZJ2：LOCATE12，22：INPUT50INPUT〃平曲线半径＝〃；R1：INPUT〃缓和曲线长＝〃；60INPUT〃平曲线半径＝〃；LMK：INPUT〃设计视距＝〃；70ZJ＝(ZJ1＋ZJ2／60＋ZJ3／3600)*PI／180：DR＝LMK／2－1.5：R2＝R1－80IFLS1＞0THENGOSUB640：REM计算90R＝R1：LS＝LS1：GOSUB710：LY1＝LY：P1＝P：Q1＝Q：T1＝100R＝R2：LS＝LS2：GOSUB710：LY2＝LY：P2＝P：Q2＝Q：T2＝110DT＝Q1－Q2：ZH2＝ZH1：SS＝S／2／120HY1＝ZH1＋LS1：YH1＝HY1＋LY1：HZ1＝YH1＋LS1：QZ1＝(ZH1＋HZ1)／130HY2＝ZH2＋LS2：YH1＝HY2＋LY2：HZ2＝YH2＋LS2：QZ2＝(ZH2＋HZ2)／140PRINT〃曲线主点桩号：〃：A＄＝〃######.##〃：IFLS1＝150PRINT〃ZH＝〃；USINGA＄；ZH1：PRINT〃HY＝〃；USINGA＄；160PRINT〃QZ＝〃；USINGA＄；QZ1：PRINT〃YH＝〃；USINGA＄；170PRINT〃HZ＝〃；USINGA＄；HZ1：180PRINT〃ZY＝〃；USINGA＄；ZH1：PRINT〃QZ＝〃；USINGA＄；190PRINT〃YZ＝〃；UXINGA＄；200INPUT〃断面桩号(＜0结束)〃；L0：IFL0＜210IFL0＜＝ZH1－SORL0＞＝HZ1＋STHENH＝0：G220IFL0＞QZ1THENL0＝ZH1＋(HZ1－L0)：REM若欲求横净距断面在曲线后半部，则计算曲线前半部对称位置的横净距230IFL0－SS＞＝HY1ANDL0＋SS＜＝YH1THENH＝R2*(1－COS(S／2／R2))：GOTO240R＝R1：LS＝LS1：P＝P1：Q＝Q1：X0＝T1＋T1*COS(ZJ)：Y0＝T1*SIN(ZJ)：IFABS(ZJ－PI)＜1E－10THENX0＝0Y0＝2*(R1＋250IFL0＜＝ZH1THENA0＝PI／2：X＝L0－ZH1：Y＝0：260IFL0＞＝HZ1THENA0＝ZJ＋PI／2：X＝HZ1－L0：Y＝0：GOSUB620：270IFL0＞YH1THENL＝HZ1－L0：A0＝ZJ＋PI／2－L^2／2／R1／LS1：GOSUB600：GOSUB620：280IFL0＞＝HY1THENL＝L0－HY1＋LS1／2：A0＝L／R1＋PI／2：GOSUB590：X＝X＋Q：Y＝Y＋P：290L＝L0－ZH1：A0＝L^2／2／R1／LS1＋PI／2：300XM＝X－DT：YM＝Y－310IFA0＝PI／2THENA1＝1：B1＝0：C1＝－XMELSEK＝TAN(A0)：A1＝K：B1＝－1：C1＝YM－320R＝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＋P2330A＝L0－S／2：B＝L0＋S／2：GOSUB460：H1＝H：A＝A＋.01：B＝B＋.01：GOSUB460：A＝A－.01：B＝B－340REM分别以10m、1m、0.1m、0.01m的步距使A、B两点向H值增大的方向移动，以减少计算次数350D＝100：IFH＜H1THENF＝－1ELSEF＝360IFH＝H1ANDB＜＝ZH2THENF＝370IFH＝H1ANDA＞HZ2THENF＝－380FORI＝1TO3：D＝D／10：D1＝390FORJ＝1TO100：A＝A＋D1：B＝B＋D1：400IFH＜＝H1THENA＝A－D1：B＝B－D1：410H1＝H：430H＝H1－440U＝CSRLIN：LOCATEU－1，48：PRINT〃横净距：〃；USING〃###.##〃；H：GOTO460REM计算视线AB与法线MN的交点E至M点的距离H值470G＝A：GOSUB520：XA＝X：YA＝480G＝B：BOSUB520：XB＝X：YB＝490A2＝YB－YA：B2＝XA－XB：C2＝(XB－XA)*YA－(YB－500AB＝A1*B2－A2*A2*B2：XE＝(B1*C2－B2*C1)／AB：YE＝－(A1*C2－A2*C1)／510H＝SQR((XE－XM)^2＋(YE－YM)^2)：520REM计算行车轨迹线已知桩号的任一点的坐标530IFG＜＝ZH2THENX＝G－ZH2：Y＝0：540IFG＞＝HZ2THENX＝HZ2－G：Y＝0：GOSUB620：550IF