技术领域
[0001] 本
发明涉及铁路设计方法，具体涉及一种既有铁路纵面线位整体智能重构方法。
背景技术
[0002] 铁路是国家重要的
基础设施，是国民经济的大动脉，在我国经济和社会发展中起着重要的作用。既有铁路增改建是我国未来铁路建设领域面临的重点任务之一，既有线路纵面重构是既有线增改建设计的重要环节；且长期的运营维护导致纵断面线形偏离原设计，造成轨道不平顺，影响行车安全及旅客舒适度，因此在常规铁路工作中，也面临着大量的铁路整正工作。开展既有线纵断面重构方法研究对于既有铁路整正及增改建设计工作具有重大意义。
[0003] 纵断面线形的重构，一些学者运用拟合方法进行既有轨面的平顺，但是由于实测点数目较多，拟合出来的线形并不能满足实际工作的需要。大部分既有线重构首先要进行直曲线元的识别，传统识别方法主要基于实测点正矢、坡度差、近似
曲率及方位
角变化进行直曲线元的识别，但是纵断面线形中，竖曲线半径大而长度小，对于实测点的偏差极为敏感，上述方法
精度均不能满足实际工程的需求；关于纵断面线形拟合，主要基于最小二乘原理进行拟合，但是大多没有考虑约束条件的限制，导致拟合的线形无法满足规范设计要求；且大多研究限于单个曲线的重构。
[0004] 因此，迫切需要一种能够实现线路纵断面整体智能重构的方法，满足实测点到重构线形抬落道量最小，降低线路整正与改建的工程成本，提高纵面重构的效率与精度。
发明内容
[0005] 针对现有铁路纵面线形拟合中存在的问题，本发明提出了一种既有铁路纵面线位整体智能重构方法，能满足实测点到重构线形抬落道量最小，降低线路整正与改建的工程成本，提高纵面重构的效率与精度。
[0006] 本发明的技术方案为：
[0007] 一种铁路纵面线形整体重构设计方法，包括以下步骤：
[0008] S1：以整条线路变坡点数量、各变坡点
位置(包括里程与高程)、竖曲线半径为优化变量，以全线抬落道量平方和最小为目标函数，并考虑各种约束条件(包括最大纵坡限制、最大及最小竖曲线半径、两相邻坡度的最大及最小坡度代数差、最小坡长、竖曲线与缓和曲线不重叠、竖曲线不与
道岔重叠、竖曲线不设置在明桥上这些约束条件)，建立纵面线位整体重构优化模型；
[0009] S2：将实测点初步区分为一系列直线与竖曲线点群，初步求解变坡点数量；
[0010] S3：对初步区分的直线与竖曲线点群进行拟合，然后根据拟合结果重新调整实测点所属的点群，基于点线一致思想进行震荡
迭代，即每次拟合之后重新调整实测点归属的点群，不断迭代，直至点群中实测点的归属不再发生变化，进一步区分直线与竖曲线，求解变坡点数量精确值及其余优化变量的值，并初步拟合出铁路纵面整体线位(直线-竖曲线-直线-竖曲线-直线…)；
[0011] S4：根据初步拟合出的铁路纵面整体线位(步骤S3得到的变坡点数量精确值及其余优化变量的值)确定优化模型的初始解，采用
粒子群优化算法求解优化模型的最优解，实现既有线纵面整体智能重构。
[0012] 进一步地，所述步骤S1具体包括：
[0013] S1-2：选取实测点抬落道量dj平方和最小为目标函数。
[0014] 根据公式(1)计算实测点到重构的直坡段的抬落道量：
[0015] dj＝kxj+b-yj (1)
[0016] 式中：j＝1，...，n'，dj为直坡段上第j个实测点处的抬落道量；n'表示直坡段上的实测点个数；
[0017] k为重构直坡段的斜率；
[0018] b为重构直坡段的切距；
[0019] xj为直坡段第j个实测点的里程；
[0020] yj为直坡段第j个实测点的高程。
[0021] 根据公式(2)计算实测点到重构的竖曲线的抬落道量：
[0022]
[0023] 式中：j＝1，...，n″，dj为竖曲线上第j个实测点处的抬落道量；n'表示竖曲线上的实测点个数；
[0024] X0为重构竖曲线圆心里程；
[0025] Y0为重构竖曲线圆心高程；
[0026] R为竖曲线半径；
[0027] S1-3：约束条件包括：
[0028] (1)重构纵面线位应高于计算轨面线，满足规范要求的最小道床厚度要求；
[0029] (2)重构的竖曲线应该与相邻的直坡段相切，满足线路平顺性要求；
[0030] (3)最小坡长、最大坡度、两相邻坡度的最小及最大代数差、最大竖曲线半径、最小竖曲线半径满足铁路线路设计规范的要求；
[0031] (4)竖曲线与缓和曲线不重叠、竖曲线与道岔不重叠以及竖曲线不设置在明桥上；
[0032] S1-4：构建纵面线位整体优化模型：
[0033] f[(X1,Y1，R1)；(X2,Y2,R2)；…(Xi,Yi,Ri)；(XN,YN,RN)]＝min∑dj2[0034]
[0035] 式中：Xi和Yi为第i个变坡点的里程和高程，Ri为第i个变坡点对应的竖曲线的半径，yj″为第j个实测点处的重构纵面线位高程，yj'为第j个实测点处的计算轨面线高程；Li为第i个坡段长度，Lmin为铁路线路设计规范允许的最小坡段长度，Li根据公式(4)计算得到：
[0036]
[0037] ii为第i个坡段的坡度，imax为铁路线路设计规范允许的最大坡度；ii根据公式(5)计算得到：
[0038]
[0039] Δii为第i个变坡点处坡度代数差，Δimin和Δimax分别为铁路线路设计规范允许的两相邻坡度的最小代数差和最大代数差。根据公式(6)计算Δii：
[0040] Δii＝ii+1-ii (6)
[0041] Rmin、Rmax分别为铁路线路设计规范允许的最小及最大竖曲线半径。
[0042] 竖曲线与直坡段相切的约束条件通过公式(37-42)具体所述原理满足。
[0043] 进一步地，所述步骤S2具体包括：
[0044] S2-1：设既有轨面线上实测点Pj坐标为(xj,yj)，其中，j＝2，3，...，l，xj为实测点Pj的里程，yj为实测点Pj的高程，根据公式(7)将高程放大e倍，得到坐标
[0045]
[0046] 其中，e根据经验设置；
[0047] S2-2：对于实测点Pj,j＝2，3，...，l-1，自P2开始，分别计算当前实测点Pj与前一个实测点Pj-1之间的距离Δxj-1，与后一个实测点Pj+1之间的距离Δxj，并判断距离是否满足(8)的约束要求；
[0048] Δxj≥dT且Δxj-1≥dT (8)
[0049] 式中：Δxj-1＝xj-xj-1,Δxj＝xj+1-xj，dT为距离
阈值，取经验值；
[0050] S2-3：若不满足式(8)，则剔除点Pj，令j＝j+1，重新进行步骤S2-2判断，直到找到满足约束要求的实测点，遍历所有的实测点，剔除所有不满足约束要求的实测点，记满足要求的实测点数为n；
[0051] S2-4：自P2至Pn-1，依次根据公式(9)计算各个实测点处的坡度值，实测点Pj处的坡度值Gj'计算公式为：
[0052]
[0053] 式中：
[0054] 式中：Δxj-1＝xj-xj-1,Δxj＝xj+1-xj；
[0055] 设置第1个实测点(起点)P1处的坡度值G1'＝G2'，最后1个实测点(终点)Pn处的坡度值Gn'＝Gn-1'；此处认为起点处坡度值与第二个实测点的坡度值相同，终点处坡度值与其前一个实测点处坡度值相同。因为坡度值的计算仅仅是为了初步区分直线与竖曲线上的实测点，一般条件下，直线上的点必须大于等于2个，曲线上的点必须大于等于3个，所以第1、2个实测点，第n-1、n个实测点应该属于同一个线元，坡度值设置一样可以保证其一定属于同一个线元；
[0056] S2-5：自P2至Pn-1，依次根据公式(10)计算各个实测点处的坡度变化率，实测点Pj处的坡度变化率Gj″计算公式为：
[0057]
[0058] 设置最后一个实测点Pn处的坡度变化率Gn″＝Gn-1″；
[0059] S2-6：根据式(11)，将实测点Pj(j＝1，...，n)处坡度变化率绝对值的平均值作为区分直线实测点或竖曲线实测点的阈值f；
[0060]
[0061] S2-7：根据各实测点的坡度变化率和阈值f粗略区分直线点群与竖曲线点群。若满足两个以上连续的实测点的坡度变化率绝对值|Pj″|均小于阈值f，则把这些连续的实测点归于一个直线点群；若不满足，则把这些连续的实测点归于一个竖曲线点群，由此建立多个直线点群Qz和竖曲线点群QY，初步将整个纵面线位划分为M个点群，其中，直线点群有个，竖曲线点群 个。
[0062] 进一步地，所述步骤S3具体包括：
[0063] S3-1：从第一个点群开始线位重构，如果第i个点群Qi为竖曲线点群QY，则进行步骤S3-2，基于拟合的原理对线位进行重构；如果为直线点群，则搜索下一个点群，其中i＝1,2,…,M；
[0064] S3-2：将竖曲线点群Qi拟合成竖曲线Ci，并基于
正交最小二乘原理，拟合与该竖曲线点群相邻的两个直线点群；
[0065] (1)如该竖曲线点群Qi为第一个竖曲线点群，则将与其相邻的两个直线点群Qi-1和Qi+1拟合成两个直坡段Li-1和Li+1；
[0066] (2)如该竖曲线点群Qi非第一个竖曲线点群，则仅将与其相邻的后一个直线点群Qi+1拟合成直坡段Li+1，与其相邻的前一个直线点群Qi-1已拟合完成的直坡段Li-1则固定不动；
[0067] 根据重构纵面线位高程应大于或等于计算轨面线高程这一约束条件，对所拟直坡段进行调整，并计算与竖曲线Ci相邻的两个直坡段Li-1和Li+1的交点坐标，即变坡点坐标(Xi,Yi)；
[0068] S3-3：前后直坡段固定，基于正交最小二乘原理，考虑竖曲线应与两端直坡段相切的约束条件，对竖曲线点群Qi进行拟合，并根据重构纵面线位高程应大于或等于计算轨面线高程这一约束条件，对所拟竖曲线进行调整；
[0069] S3-4：计算调整后竖曲线范围(竖曲线切线长)T＝|Rtanβ|，式中，R为拟合得到的竖曲线半径，β为竖曲线转角的一半；
[0070] S3-5：过变坡点分别作直线Li-1与直线Li+1的垂线VLi-1、VLi+1，斜率分别为1/ki-1和1/ki+1，垂线VLi-1、VLi+1的方程：
[0071] y＝(-1/ki-1)x+bi-1 (12)
[0072] y＝(-1/ki+1)x+bi+1 (13)
[0073] 将变坡点(Xi,Yi)分别带入式(12)、(13)，求得bi-1，bi+1值，建立垂线VLi-1、VLi+1方程：
[0074] VLi-1:y＝(-1/ki-1)x+[Yi-(1/ki-1)Xi] (14)
[0075] VLi+1:y＝(-1/ki+1)x+[Yi-(1/ki+1)Xi] (15)
[0076] S3-6：判断竖曲线点群Qi和与其相邻的直线点群Qi-1内的实测点(x,y)，是否满足约束要求(16)：
[0077]
[0078] 如果满足要求，则将其分到竖曲线点群Qi内，否则，将其划分到直线点群Qi-1内。
[0079] S3-7：判断竖曲线点群Qi和与其相邻的直线点群Qi+1内的实测点(x,y)，是否满足约束要求(17)：
[0080]
[0081] 如果满足要求，则将其分到竖曲线点群Qi内，否则，将其划分到直线点群Qi+1内。
[0082] S3-8：根据步骤S3-6、S3-7结果对竖曲线点群Qi及其前后相邻的直线点群Qi-1和Qi+1内的实测点重新划分点群，对重新分类后的实测点重复进行步骤S3-2到S3-7，直至这些实测点所属的点群均不再发生变化。
[0083] S3-9：遍历所有点群，直至所有的点群内实测点不再发生变化，实现整个线路范围内的直线点群与竖曲线点群的重新划分，将解出的变坡点数量作为精确值。
[0084] S3-10：连接整个线路范围内的直线点群与竖曲线点群拟合形成的线位，生成整体线位初步方案(直线-竖曲线-直线-竖曲线-直线…)。
[0085] 进一步地，所述步骤S4粒子群
优化算法基本原理为：
[0086] 粒子群优化算法将优化问题中的可行解抽象为多维搜索空间中的粒子，对每个粒子赋予位置、速度、适应度、学习与记忆能
力等特征。所有的粒子跟随当前群体中的最优粒子在解空间内“飞行”。在飞行的过程中，粒子间按照公式(18)所述规则对粒子的位置与速度进行更新，从而形成新的粒子群。
[0087] 粒子群中第m个粒子第t次迭代过程中的速度向量与位置向量更新公式为：
[0088]
[0089] 式中： 和 分别为第m个粒子第t次迭代后的速度向量和位置向量，第l个粒子的速度向量的初始值 迭代次数1≤t≤1000， 为第m个粒子的初始位置向量(初始粒子群中第m个粒子的位置向量)；w为速度权重系数，c1、c2为
加速常数，其中，c1表示粒子向自身最优位置运动的步长，c2为粒子向整个粒子群最优位置运动的步长， 为[0,1]范围的随机数， 为第m个粒子的个体最优解(历史最优位置)，且 gt为粒子群的全局最优解(历史最优位置)，g0为初始粒子群中的最优位置；
[0090] 第一次更新粒子位置向量与速度向量的公式(20)：
[0091]
[0092] 式中： 为粒子l跟随全局最优粒子位置“飞行”1次后速度向量， 为[0,1]范围的随机数， 为第m个粒子跟随全局最优粒子位置“飞行”1次后的位移向量。
[0093] 进一步地，所述步骤S4具体实现为：
[0094] S4-1：初始化一个包含B个粒子的粒子群，每个粒子的位置向量代表一个纵面线位整体重构方案；第m个粒子的位置向量记为sm＝(X1m,Y1m,R1m,X2m,Y2m,R2m,…,Xim,Yim,Rim,…,XNm,YNm,RNm)T，第m个粒子的位置向量记为vm，sm和vm的维度均为4N；其中，sm中的元素Xim、Yim和Rim分别表示第m个纵面线位整体重构方案中第i个变坡点的里程、高程和对应的竖曲线半径，N为变坡点总个数，N个变坡点包括步骤3计算得到的所有变坡点，以及第1个实测点P1和最后1个实测点Pn，其中P1和Pn对应的竖曲线半径设置为零；以N个变坡点中第i个变坡点的里程Xi、高程Yi和对应的竖曲线半径Ri为基准，在各自允许的变化范围(各变量允许的变化范围根据约束条件确定)内随机变化得到Xim、Yim和Rim的初始值，由此得到第m个粒子的位置向量sm的初始值 第m个粒子的速度向量vm的初始值 设置各个粒子的个体最优解为其初始位置向量；
[0095] S4-2：遍历粒子群中所有粒子，依次根据公式(20)计算粒子群中各个粒子的位置向量代表的纵面线位整体重构方案下所有实实测点抬落道量平方和，即粒子群中各个粒子的适应度，进而评估解(各粒子位置向量)的优劣；
[0096]
[0097] 若计算出的粒子适应度优于该粒子的个体最优解pbestj，则将该粒子的个体最优解pbestj更新为该粒子当前的位置向量；各粒子的个体最优解pbestj的初始值为其初始位置向量；
[0098] 若计算出的粒子的适应度优于全局最优解Gbest，则将粒子群的全局最优解Gbest更新为该粒子当前的位置向量；粒子群中全局最优解Gbest的初始值为初始粒子群中的全局最优解；
[0099] S4-3：根据公式(18)迭代更新粒子群中的各个粒子的速度向量与位置向量，并重复步骤S4-2计算各个粒子的适应度，更新粒子群中全局最优解与粒子个体最优解；不断迭代，每次迭代形成新的粒子群位置，即优化模型产生新的解集，直至全局最优解前后两次迭代不发生变化或达到最大迭代次数为止，得到最优的纵面线位整体重构方案。
[0100] 所述步骤S3-2中基于最小二乘原理，进行直坡段线路拟合方法具体如下：
[0101] 设直线点群Qi-1拟合形成的直坡段方程为：
[0102] y＝ki-1x+bi-1 (21)
[0103] 式中两个待拟合参数，ki-1代表斜率，bi-1代表截距。对于直线点群Qi-1中包含的任意实测点Pj(xj,yj)到直坡段的抬落道量(拨量)为：
[0104] dj＝ki-1xj+bi-1-yj(j＝1，2，...，ni-1)
[0105] (22)[0106] 式中：ni-1表示直线点群Qi-1内实测点个数；
[0107] 基于最小二乘原理的直线拟合需使式(23)的值最小，即抬落道量平方和作为直坡段重构效果的评价函数，则：
[0108]
[0109] 当目标函数取得最小时，直坡段完成重构，则对ki-1和bi-1求偏导数，令偏导数为零。
[0110]
[0111]
[0112] 求得最优重构直坡段ki-1和bi-1最佳估计值为：
[0113]
[0114]
[0115] 所述步骤S3-2中，计算轨面线具体计算原理如下：
[0116] 根据外业实测的既有轨面实测点数据、既有道床厚度及既有轨道高度测点数据以及线路最小道床厚度、轨道高度计算出实测点的计算轨面线数据，作为线路纵面线位拟合的约束条件。
[0117] (1)建立既有轨面线实测点数据集{Pj(xj,yj)|j＝1，2，...，n}、既有道床厚度测点数据集{Bi(xi,ti)|i＝1，2，...，k}和既有轨道高度测点数据集{Ci(xi,ci)|i＝1，2，...，k}，其中，(xj,yj)表示第j个既有轨面线实测点Pj处的里程和高程，Bi(xi,ti)表示第i个既有道床厚度测点Bi处的里程和道床厚度，(xi,ci)表示第i个既有道床厚度测点Ci处的里程和轨道高度；既有道床厚度测点与既有轨道高度测点数目相同，且测点里程一一对应；
[0118] 从既有轨面线实测点P1处开始，遍历所有的既有轨面线实测点，对于Pj，判断是否存在与其里程相等的既有道床厚度测点数据和既有轨道高度测点数据，否存在，则可以直接得到该实测点出的既有道床厚度tj和既有轨道高度cj；否则，根据Pj前后与其里程最接近的两个既有道床厚度测点Bi-1和Bi，以及Pj前后与其里程最接近的两个既有轨道高度测点Ci-1和Ci，利用线性内插公式(28)计算出实实测点Pj处的既有道床厚度tj和既有轨道高度cj：
[0119]
[0120] (2)根据公式(29)计算出实测点Pj处的既有道床底面高度Yj(j＝1，...，n)；
[0121] Yj＝yj-tj-cj(29)
[0122] (3)根据公式(30)计算出实测点Pj处的计算轨面线高程yj'(j＝1，...，n)；
[0123] yj′＝Yj+t′j+c′j (30)
[0124] 其中，t'j和c'j为铁路线路设计规范允许的最小道床厚度与轨道高度；
[0125] 所述步骤S3-2中，根据重构纵面线位高程应大于或等于计算轨面线高程这一约束条件，对所拟直坡段进行调整，具体原理如下：
[0126] 首先对属于前直坡段上的实测点进行检测。
[0127] 根据式(31)计算前直坡段上的实测点Pj(xj,yj)处的重构纵面线位高程y'j'：
[0128]
[0129] 遍历前直坡段上所有的实测点，根据公式(32)判断实测点Pj处的重构纵面线位高程y'j'是否小于计算轨面线高程yj'：
[0130] Δyj＝y″j-yj'＞0 (32)
[0131] 设满足公式(32)的实测点有a个，则遍历a个实测点，寻找使Δyj达到最大值的实测点Pj0(xj0,yj0)，记此时Δyj＝Δyj0，将直坡段整体向上平移Vyj0，得到直坡段的方程为：
[0132] y＝ki-1'x+bi-1'+Δyj0 (33)
[0133] 同理对后直坡段拟合直线进行调整。
[0134] 所述步骤S3-2中，计算变坡点坐标(Xi,Yi)具体原理如下：
[0135] 设两直线点集Qi-1和Qi+1拟合的直线Li-1和Li+1方程分别为：
[0136] y＝ki-1x+bi-1 (34)
[0137] y＝ki+1x+bi+1 (35)
[0138] 则两直线的交点坐标(Xi,Yi)为：
[0139]
[0140] 所述步骤S3-3中，基于正交最小二乘原理，并考虑竖曲线边界与直线相切这一约束条件，对竖曲线点群进行拟合，具体原理如下：
[0141] 设变坡点处角平分线所在的直线的表达式如下：
[0142] kx-y+b＝0 (37)
[0143] 为其方向单位向量(单位向量，且模等于1)；
[0144] 根据公式(38)求得经过变坡点处角平分线所在直线的斜率k。
[0145]
[0146] 根据公式(39)计算竖曲线圆心的坐标：
[0147]
[0148] 以拨道平方和最小为约束条件，建立目标函数如下：
[0149]
[0150] 根据式(41)化简目标函数：
[0151]
[0152] 令：
[0153] a＝ni(A2+B2-1)
[0154]
[0155]
[0156] 根据式(42)求解目标函数最小值，得到竖曲线半径最佳估计值：
[0157]
[0158] 根据式(43)求解竖曲线范围T：
[0159] T＝|R tanβ| (43)
[0160] 所述步骤S5-3中，根据重构纵面线位应大于或等于计算轨面线高程这一约束条件，进行竖曲线的调整，具体原理如下：
[0161] 根据公式(31)、(32)计算求出竖曲线范围内使Δyj达到最大值的实测点Pj1(xj1,yj1)，记该点处计算轨面线比重构纵面线位高出Δyj1。设该点处对应的计算轨面线的数据为yj1'，设调整后的竖曲线半径为R'，圆心坐标为(X0，Y0')，为满足抬落道量平方和最小的目标，设Pj1'(xj1,yj1')在调整后的竖曲线上。
[0162] 此时竖曲线满足：
[0163]
[0164] 式中仅有Y0'为未知数，求得：
[0165]
[0166] 式中：
[0167] 则调整后的竖曲线方程为：
[0168] (x-X0)2+(y-Y0')2＝R'2 (46)
[0169] 有益效果：
[0170] 本发明提供的一种既有铁路纵面线位整体智能重构方法，能满足实测点到重构线形抬落道量最小，降低线路整正与改建的工程成本，提高纵面重构的效率与精度。
[0171] 由于变坡点数量、各变坡点位置(包括里程与高程)及竖曲线半径可以将纵面线形唯一确定，纵面重构的目的在于使重构的线路尽可能接近既有轨面线，即既有轨面线上实测点距重构线路的距离和尽可能小，所以本发明以变坡点个数、各变坡点位置、竖曲线半径为优化变量，以全线抬落道量平方和最小为目标函数，考虑最大纵坡限制、最大及最小竖曲线半径、两相邻坡度的最大及最小坡度代数差、最小坡长、竖曲线与缓和曲线不重叠、竖曲线不与道岔重叠、竖曲线不设置在明桥上这些约束条件，构建纵面线位整体重构模型。因自变量的个数不确定，会给模型求解带来很大的困难，且初始解的选取会影响到优化模型的收敛性。因此，模型分两阶段进行求解，第一阶段精确求解变坡点个数，并求解出其它优化变量的初始值，给出初步线路重构方案；第二阶段将变坡点个数作为已知量，根据第一阶段求解的优化变量的初始值产生多个可行解，将每个可行解分别作为一个粒子的位置向量，引入粒子群优化算法对各个粒子的位置向量进行优化，得到最优位置向量，作为模型最优解，实现纵面线位整体重构。分两阶段逐步求解的方法，避免了在很大的空间内搜索可行解带来的复杂性，且第一阶段模型求解中，首先基于高程放大了一定倍数的实测点数据，计算既有轨面线实测点处的坡度变化率，进行直线与竖曲线点群的初步区分，初步确定变坡点个数，实验证明，将实测点数据中的高程放大一定倍数来计算坡度变化率，减少了纵面坡度代数差较小引起的坡度变化率对测量误差的敏感性，提高了区分直线与竖曲线点群的精度。并提出基于点线一致思想的震荡迭代算法进一步区分直线与竖曲线点群，求解变坡点个数精确值；每次迭代基于最小二乘拟合原理进行直线与竖曲线拟合，且保证实测点所属的点群与实测点参与拟合的线形相一致，求得模型初始解，完成线位初步重构；优化变量初始值的选取较为合理，提高了第二阶段模型求解的精度与效率。
附图说明
[0173] 图2是实测点加密示意图。
[0174] 图3是考虑竖曲线与前后直坡段相切为约束条件的竖曲线拟合示意图。
[0175] 图4是考虑实测点处重构纵面线位应该大于或等于计算轨面线的约束条件的直坡段调整示意图。
[0176] 图5是考虑实测点处重构纵面线位应该大于或等于计算轨面线的约束条件的竖曲线调整示意图。
[0177] 图6是基于点线一致思想震荡迭代区分直线与竖曲线点群示意图；图6(a)、图6(b)和图6(c)分别为初步拟合得到的点群、第1次迭代、第N次迭代得到的点群；
[0178] 图7是粒子群算法全线优化步骤流程图。
具体实施方式
[0179] 本
实施例以石
门-长沙既有铁路纵面线位重构设计为例对本发明的技术方案做进一步说明。如图1所示，本发明所涉及的既有铁路纵面线位重构设计方法包括以下步骤：
[0180] S1：以变坡点数量N、各变坡点位置(包括里程X与高程Y)、各竖曲线半径R为优化变量，以全线整体抬落道量平方和∑dj2最小为目标函数，并考虑各种约束条件，建立纵面线位整体重构优化模型；其中约束条件包括：
[0181] (1)重构的纵面线位应高于计算轨面线，满足规范要求的最小道床厚度要求。
[0182] (2)重构的竖曲线应该与相邻的直坡段相切，满足线路平顺性要求。
[0183] (3)石门-长沙既有铁路为Ⅰ级客货共线铁路，设计速度120km/h，全线位于湖南省境内，沿线大部分
地貌属于丘陵区。查规范得到：最小坡长选取400m，最大坡度取12‰，最小竖曲线半径取10000m，最大竖曲线半径取40000m，最大坡度代数差取8‰，最小坡度代数差取2.8‰。
[0184] (4)竖曲线与缓和曲线不重叠、竖曲线与道岔不重叠以及竖曲线不设置在明桥上，这些约束条件通过约束检测进行调整。
[0185] 本实施例构建的纵面线位整体重构优化模型为：
[0186] f[(X1,Y1，R1)；(X2,Y2,R2)；…(Xi,Yi,Ri)；(XN,YN,RN)]＝min∑dj2[0187]
[0188] S2：将实测点高程放大e倍，基于高程放大e倍后的实测点数据求实测点处坡度变化率，将实测点粗略区分为一系列直线与竖曲线点群，初步求解变坡点数量；本实施例中设定e＝10；
[0189] S2-1：根据公式(7)计算各个实测点处的高程放大10倍的坐标(xj,yje)。本实施例中共有1500个既有轨面数据实实测点数据；
[0190] S2-2：对于实测点Pj,j＝2，3，...，1499，自P2点开始，分别计算当前实测点Pj与前一个实测点Pj-1之间的距离Δxj-1，与后一个实测点Pj+1之间的距离Δxj，并判断距离是否满足(8)的约束要求，本实施例中设定dT＝15m；
[0191] S2-3：若不满足式(8)，则剔除点Pj，令j＝j+1，重新S2-2进行判断，直到找到满足约束要求的实测点，遍历所有的实测点，剔除所有不满足约束要求的实测点。满足约束条件的点数n＝1356；
[0192] S2-4：自P2至Pn-1，依次依据式(9)计算各个实测点处的坡度值Gj'；设置第1个实测点P1处的坡度值G1'＝G2'，最后1个实测点Pn处的坡度值Gn'＝Gn-1'；
[0193] S2-5：自P1至Pn-1，依次根据公式(10)计算各个实测点处的坡度的变化率Gj″；设置最后一个实测点Pn处的坡度变化率Gn″＝Gn-1″；
[0194] S2-6：根据式(11)，将实测点Pj(j＝1，2，...，1356)处坡度变化率绝对值的平均值作为区分直线实测点或竖曲线实测点的阈值f，计算得到阈值为4×10-6；
[0195] S2-7：根据各实测点的坡度变化率和阈值f粗略区分直线点群与竖曲线点群。若满足两个以上连续的实测点的坡度变化率绝对值|Gj″|均小于阈值f，则把这些连续的实测点归于一个直线点群；若不满足，则把这些连续的实测点归于一个竖曲线点群，由此建立多个直线点群和竖曲线点群。本实施例中按上述方法将整个纵面线位划分为直线点群120个，竖曲线点群119个，共有点群239个，解出变坡点数量初始值为119个。
[0196] S3：对初步区分的直线与竖曲线点群进行拟合，然后根据拟合结果重新调整实测点所属的点群，基于点线一致思想进行震荡迭代，使多次拟合后，点群中实测点的归属不再发生变化，进一步区分直线与竖曲线，并初步拟合出纵面整体线位(直线-竖曲线-直线-竖曲线-直线…)，求解出的变坡点数量作为精确值。
[0197] 其中直线与竖曲线的拟合考虑实测点处重构纵面线位数据应大于计算轨面线数据这一约束条件，首先对直坡段进行拟合与调整，然后在直坡段固定的基础上进行竖曲线拟合与调整。计算轨面线计算过程如下：
[0198] (1)根据实测既有轨面数据Pj(xj,yj)(j＝1，...，1356)、既有道床厚度数据Bi(xi,bi)(i＝1，...，1200)和既有轨道高度数据Ci(xi,ci)(i＝1，...，1200)，既有道床厚度与既有轨道高度实测点数目相同，且实测点里程一一对应，对实测点数据进行加密处理。
[0199] 对于各既有轨面实测点Pj(xj,yj)(j＝1，...，1356)，判断是否存在与其里程相等的既有道床厚度实测点数据和既有轨道高度实测点数据，否存在，则可以直接得到该实测点出的既有道床厚度tj和既有轨道高度cj；否则根据Pj前后与其里程最接近的两个既有道床厚度实测点Bi-1和Bi，以及Pj前后与其里程最接近的两个既有轨道高度实测点Ci-1和Ci，利用线性内插公式(28)计算出实实测点Pj处的既有道床厚度tj和既有轨道高度cj。
[0200] 遍历所有实测点，使每个既有轨面实实测点Pj处均有相应的既有道床厚度数据和既有轨道高度数据。
[0201] (2)根据公式(29)计算实测点既有道床底面高度Yj(j＝1，...，1356)；
[0202] (3)按照铁路线路设计规范允许的最小道床厚度与轨道高度，根据公式(30)计算出实实测点Pj处的计算轨面高程数据yj'(j＝1，...，1356)。
[0203] S4：以初步拟合的整体线位作为初值，基于粒子群算法对优化模型进行精确求解，完成纵面线位整体智能重构。本实施例中初值为设置粒
子群大小为50，最大迭代次数100。最终得到优化模型最优解为:(97700,48.84,0,
980025.3,49.30,30085,…,98625.12,50.45,18500.30,159500,54.09,0)。
