圆曲线与加权曲线的简解.docxVIP

圆曲线与加权曲线的简解 1 标的的解算与讨论 随着城市建设的快速发展，城市化以前所未有的速度发展起来。它的类型包括快速道路和交叉口道路。直线、圆形曲线和平滑曲线等形状。城市道路与绕城的高等级公路以及铁路的连接成网, 必然会经常涉及到线路及其平行线相交问题, 因此研究线路及其边线相交点坐标的解算就显得很有必要和有意义了。解算内容包括直线与直线、直线与缓和曲线、直线与圆曲线、圆曲线与缓和曲线等交点以及它们的平行线交点的解算, 这其中尤以圆曲线与缓和曲线及其平行线交点坐标的解算更复杂。现有文献已对缓和曲线与圆曲线的交点解算有研究, 但过程甚为麻烦, 且难于实现圆曲线与缓和曲线平行线交点的解算。以前相当长的一段时间内, 人们对缓和曲线平行线的性质认识不清, 直到九十年代中期, 文献首次基于严密正确弧长公式证明了缓和曲线平行线不再是缓和曲线的结论, 从而为其测设及曲线平行线交点的严密解算奠定了基础。本文所给出的圆曲线与缓和曲线平行线交点坐标的解算方法, 其前提就是基于上述结论而建立的缓和曲线平行线以其自身弧长为变量的参数方程。虽然, 实践中经常需求圆曲线平行线与缓和曲线平行线的交点, 但由于圆曲线的平行线与圆曲线本身为同心圆弧, 因此本文仅研究圆曲线与缓和曲线平行线交点坐标的解算问题。 2 缓和曲线方程的阶次 缓和曲线实质是数学上的科努螺旋线, 线路上用它来连接直线与圆曲线或不同半径的圆曲线。由于缓和曲线是高次曲线, 圆曲线是二次曲线, 因此求解两非线性曲线的交点并非容易。现有文献采用的都是联立曲线方程式的解法, 这样, 解算原理虽然简单, 但计算过程却十分复杂。文献介绍的解法是首先将圆曲线方程转换至缓和曲线的切线坐标系中, 然后将缓和曲线的参数方程代入圆曲线的二次方程, 由于缓和曲线方程本身即含有高次项, 因而就产生了至少10次以上的高次方程, 而且该方程的阶次还与缓和曲线方程的取项密切相关, 当对解算精度有较高要求时, 方程的阶次可达18阶以上。显然, 直接解算该高次方程并非易事。文献采用的是牛顿逼近法, 这不仅需多次迭代和计算十次方以上的数值, 而且在计算前还需人工给出两曲线交点的一个近似区间;解算后, 还需再将缓和曲线的切线坐标系中的坐标转换至所需的城市坐标系。文献给出的解法, 其原理与文献完全一样, 只不过它采用的是缓和曲线的直角坐标方程。 现有文献尚未讨论圆曲线与缓和曲线平行线交点的解算, 显然, 如果沿用上述所谓的“直接解法”, 计算过程将十分复杂;再考虑到缓和曲线平行线不是缓和曲线, 计算实际上将难以进行。因此, 本文首先寻找到一种求解圆曲线与缓和曲线交点的简便方法, 在此基础上利用缓和曲线平行线以其自身弧长为变量的参数方程, 从而给出圆曲线与缓和曲线平行线交点坐标的解算方法。 3 计算曲线和缓冲曲线之间的点坐标 3.1 外标系统1 文献给出的缓和曲线通用参数模型为: x=xΖΗ+(l-l540R2lS2)cosαΖΗ→JD±(l36RlS-l7336R3lS3)sinαΖΗ→JDy=yΖΗ+(l-l540R2lS2)sinαΖΗ→JD?(l36RlS-l7336R3lS3)cosαΖΗ→JD}(1)x=xZH+(l?l540R2lS2)cosαZH→JD±(l36RlS?l7336R3lS3)sinαZH→JDy=yZH+(l?l540R2lS2)sinαZH→JD?(l36RlS?l7336R3lS3)cosαZH→JD???(1) 式中:X、Y为缓和曲线上任意点i在地面坐标系中的坐标;αZH→JD为ZH至JD在地面坐标系中的方位角;R为与圆曲线相接圆的半径;lS为缓和曲线长;“±”表示曲线左偏时取正号, 右偏时取负号;“?”表示曲线左偏时取负号, 右偏时取正号;为至的弧长。 3.2 缓和曲线与圆曲线点的距离计算 如图1所示, 已知ZH点在地面线路坐标系 (或城市坐标系) 中的坐标 (XZH、YZH) ;ZH至JD的方位角αZH→JD;缓和曲线lS长和与缓和曲线相接圆的半径R;与缓和曲线相交圆的圆心O在地面坐标系中的坐标 (XO、YO) 及半径r, 现求缓和曲线段与圆曲线段相交点P在地面坐标系中的坐标 (XP、YP) 。 分析:要求得交点P, 只须求得ZH到P的弧长l, 即可由 (1) 式求得 (XP、YP) 。由缓和曲线方程求得的点必须同时满足圆曲线方程即图1中的点P到圆心O的距离等于圆半径r才是所要求的交点, 因此以某种方法取得弧长L的近似值, 以由弧长近似值求得近似点到圆心O的距离作为判断是否需要循环计算及最终结果的标准是十分巧妙的做法。这里采用r与O至近似点距离之差作为L的近似值, 如图1, 以ZH点为起点在缓和曲线上截取 (r-l1) 弧长得一点P1, 以 (r-l1) 代入缓和曲线方程可算得点P1的坐标