CN103399849B - 一种基于改进型粒子群算法的道路三维线形自动优化方法 - Google Patents
一种基于改进型粒子群算法的道路三维线形自动优化方法 Download PDFInfo
- Publication number
- CN103399849B CN103399849B CN201310252969.3A CN201310252969A CN103399849B CN 103399849 B CN103399849 B CN 103399849B CN 201310252969 A CN201310252969 A CN 201310252969A CN 103399849 B CN103399849 B CN 103399849B
- Authority
- CN
- China
- Prior art keywords
- point
- length
- dimensional linear
- plane
- road
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于改进型粒子群算法的道路三维线形自动优化方法,首先建立道路三维线形优化目标函数,依据设定线路平面及纵断面的约束条件、道路建设单价、线路取点间隔及线形优化目标区域,采用改进型粒子群算法对利用ArcGIS软件获得导向线方案的三维线形的平面交点及纵断面变坡点参数同时进行优化,获得自动优化后的道路三维线形;自动优化后的三维线形能较好的绕避障碍物和适应地形,作为工程设计人员提供前期设计参考,降低设计人员劳动强度。
Description
技术领域
本发明涉及公路或铁路的选线领域,特别是一种基于改进型粒子群算法的道路三维线形自动优化方法。
背景技术
依据铁(公)路的基本的走向及主要技术标准,自动选线目的在于同时实现铁路或公路的平面和纵断面的自动设计。自动选线过程中需要从满足约束条件的多个线路方案中找到一条最优的方案(通常是经济费用最省)。由于其目标函数要通过数字高程模型内插才能获得,并且包含了各种费用,通常该函数是不可微的。由于起点终点之间满足约束条件的线路方案有无穷多个,在优化的过程中,需要考虑项目建设中众多复杂的因素,而各因素间相互作用、互相制约,因此,自动选线是一个十分复杂的问题。现有技术中对该问题的优化通常是对于线路的平面和纵断面的选线分别进行优化,容易陷入局部最优解,得到的选线结果只能作为局部优化且这种模型的计算效率不高。
随着地理信息***(GIS)技术的发展,应用GIS技术进行道路自动选线。然而现有的道路路选线设计方法由于不能够直接对地理信息进行识别、分析和利用,选线结果不尽人意,难以实现实际应用。GIS现有的地理信息快速的获取方法、科学的表达模式以及强大的空间检索分析能力为自动选线投入实用提供了技术支持。如文献[1]:Min Wook Kang,Paul Schonfeld,NingYang.Prescreening and Repairing in a Genetic Algorithm forHighway Alignment Optimization[J].Computer-Aided Civil and InfrastructureEngineering.24(2009)109–119。然而,现有技术中并未提供一种利用GIS同时进行线路的平面和纵断面的选线优化。
发明内容
本发明提供了一种基于改进型粒子群算法的道路三维线形自动优化方法,其目的在于,克服现有技术的不足,同时优化线路的平面和纵断面的选线,避免在优化过程中出现的局部最优问题。
一种基于改进型粒子群算法的道路三维线形自动优化方法,首先建立道路三维线形优化目标函数,依据设定线路平面及纵断面的约束条件、道路建设单价、线路取点间隔及线形优化目标区域,采用改进型粒子群算法对利用ArcGIS软件获得导向线方案的三维线形的平面交点及纵断面变坡点参数同时进行优化,获得自动优化后的道路三维线形;
所述三维线形的平面交点及纵断面变坡点参数是指平面交点的平面坐标(x,y)、半径R及缓长l,纵断面变坡点里程mil及高程z;
三维线形上所有平面交点及纵断面变坡点参数作为改进型粒子群算法中的粒子k;
所述改进型粒子群算法模型为:
其中,vkd t为粒子k第d维在第t次迭代的速度分量;ω为权重系数设定值;xkd t为粒子k第d维当前的位置分量;t为迭代的次数;int为向下取整函数,原型为int(number),取小于或等于number的最大的整数;ceil为向上取整函数,原型为ceil(number),取大于或者等于number的最小整数;
为整数区间内均匀分布的随机数,当时, 当时,
为整数区间内均匀分布的随机数,当时, 当时,
c1为控制粒子在飞向自身最好位置方向的步长;c2为控制粒子在飞向全局最好位置方向的步长;rand1、rand2是两个介于[0,1]范围的随机数;
pbestkd t表示粒子k在第t次迭代之前的历史最好位置;gbestkd t表示种群在第t次迭代之前的全局历史最好位置;
c1取2.0,c2取2.0。
所述的基于改进型粒子群算法的道路三维线形自动优化方法,具体包括以下步骤:
步骤1:建立道路三维线形目标函数F:
式中,i为横断面序号;b为桥梁的序号;t为隧道的序号;j为土地类型序号;IK为线形设计区间中横断面的个数;TK为线形设计区间中隧道座数;BK为设计区间中桥梁座数;YK为用地类型总个数;C为土石方工程单价;Cb为第b座桥梁的施工单价;Ct为第t座隧道的施工单价;Cj为第j种类型土地的用地单价;Si为第i个设计横断面的面积;Δw为两相邻横断面间的距离;Lb为第b座桥梁的长度;Lt为第t座隧道的长度;Aj为第j种土地类型的用地面积;
沿三维线形平面的法向量方向对线路进行剖分获得线路的横断面;
步骤2:设定三维线形平面及纵断面的约束条件、道路建设单价、三维线形平面取点间隔及三维线形优化目标区域xmin≤x≤xmax,ymin≤y≤ymax;
(1)平面约束条件
(a)最小曲线半径约束,平面交点Pi平曲线半径为r1i,不小于设定的最小曲线半径rmin,rmin≤r1i(i=1,2,...,n-1);
(b)最短圆曲线长度约束,两相邻缓和曲线间的圆曲线长度不小于设定的最短圆曲线长度
式中:为圆曲线长度;αi为圆曲线夹角;
(c)最短夹直线约束,平面夹直线最小长度为平面夹直线长度Li,i+1-Ti-Ti+1必须大于设定的平面夹直线最小长度为
式中:Li,i+1为平面交点i至交点i+l间的距离;Ti为平面交点i处的切线长;Ti+1为平面交点i+1处的切线长;
(2)纵断面约束条件
(a)坡长约束,纵断面设计坡长须满足最小坡长的限制;两相邻竖曲线的切线长之和小于该坡段长度,即:
其中,lcj-1,lcj分别为纵断面变坡点j-1及j的里程;分别为纵断面变坡点j-1及j处的切线长;
(b)坡度约束,纵断面第j段的设计坡度ij应该满足规范规定的允许坡度,即:
式中:hj,hj-1分别为纵断面变坡点j与j-1的设计高程;为纵断面坡段j的坡段长;imin,imax分别为规定的最小、最大坡度;
(c)相邻坡段坡度差约束,相邻坡段的坡度差的绝对值不超过允许的最大值Δimax,即
(d)控制点高程约束
线路在指定点处的设计高程hj应当满足最大最小高程(hmax,hmin)要求,即:
hmin≤hj≤hmax
(3)线形优化目标区域设定为xmin≤x≤xmax,ymin≤y≤ymax,同时满足以下条件:
1)禁止通过设定区域,即禁区;
设置禁区的地价费用或地面高程为无穷大,即从目标区域中排除禁区;
2)必须通过设定区域;
设平面必经区域中心点为O,计算O点距离线路中心线的距离do,设定较小值ε,并满足do≤ε,则线路可以从该区域通过,即该区域中必须包含有所选线路;
步骤3:从三维线形优化目标区域中,采用改进型粒子群算法寻找最优线形,适应度函数为F;
步骤a)采用由ArcGIS软件获得导向线方案的三维线形作为初始线形,初始线形上所有平面交点及纵断面变坡点参数作为粒子,设置粒子群种群大小、迭代次数及初始种群,利用改进型粒子群算法搜索更新初始线形上的所有交点的位置,输出最好解x1,即粒子位置构成的最佳三维线形;
步骤b)若步骤a)中搜索得到的三维线形中包含非三维线形优化目标区域中的位置,则采用静态搜索更新步骤a)得到的最好解x1在三维线形优化目标区域中的速度,利用速度更新最好解x1,获得第二最好解x2;
步骤c)以步骤b)获得的第二最好解x2作为初始探测点进行轴搜索,获得第三最好解x3;
步骤d)将步骤c)获得的第三最好解x3输入ArcGIS软件,自动获得道路最优线形。
所述步骤3中步骤b)的静态搜索的搜索过程如下:
对于落入不可行区域的粒子,速度更新方程中,设定惯性权重系数ω=0,使得粒子新的速度与原飞行无关,只取决于粒子的当前位置、历史最好位置和粒子群的历史最好位置,即速度更新方程为:
vkd t+1=c1·rand1·(pbestkd t-xkd t)+c2·rand2·(gbestkd t-xkd t) (16)
其中:vkd t为粒子k第d维在第t次迭代的速度分量;ω为权重系数设定值;xkd t为粒子k第d维当前的位置分量;t为迭代的次数;c1控制粒子在飞向自身最好位置方向的步长;c2控制粒子在飞向全局最好位置方向的步长;rand1、rand2是两个介于[0,1]范围的随机数;pbestkd t表示粒子k在第t次迭代之前的历史最好位置;gbestkd t表示种群在第t次迭代之前的全局历史最好位置;
c1取2.0,c2取2.0。
所述步骤3中步骤c)的轴搜索的搜索过程如下:
以第二最好解x2作为第一轮探测即初始探测点为x0,依次沿d维坐标轴的方向进行以下探测过程搜索得到多个探测点,从多个探测点中选择使得目标函数值F最小的点作为轴搜索得到的最好点;
设放大因子α>1,缩减因子β∈(-1,0),沿各轴方向的步长为δ1,δ2,...,δd;
每轮探测的起点和终点分别为y(1)和ye(即yn+1),令y(1)=x(k),x(k)为第k次转轴后的搜索起点,d(i)表示d维坐标轴中第i个坐标轴的方向,按如下规则进行探测;
y(i+1)=y(i)+δid(i)
从以x(k)为基点沿各个方向探测所得的点y(j)对应的目标函数值F中,选择一个目标函数值最小的对应点作为好解将好解与x(k)点对应的目标函数值F比较进行比较,按照下式更新步长:
判断是否满足以下任一条件,若满足,即完成探测,获得当前得到的好解并将作为第三最好解x3,否则进入下一轮探测:
(1)沿各轴方向不能找到比初始探测点x0的目标函数值更小的点y(j);
(2)各个方向的步长δj均达到设定的最小值ε,ε的取值范围是0.01-0.1。
有益效果
本发明提供了一种基于改进型粒子群算法的道路三维线形自动优化方法,通过对同一目标区域进行铁路和公路自动选线优化测试结果分析可知:应用本发明提出的理论和方法能快速高效地同时进行线路的平面和纵断面的选线优化,优化后的三维线形能很好的绕避障碍物和适应地形,可以作为工程设计人员提供前期设计参考,降低设计人员劳动强度。
①方案多样性好、效率高。最终生成并非单一方案,而是一个包含绕避方案和架设桥隧方案等不同特点的最优方案群。利用Rosenbrock方法与PSO算法,保证了优化后的线路方案的多样性,为线路设计人员提供更多有参考价值的比选方案,避免了由于设计人员主观因素而遗漏可行方案。
②全局最优化搜索方法,避免陷入局部最优。改进的PSO算法从理论上保证了线形优化在三维搜索空间中的全局最优,更容易搜索到全局最优方案。
③适用于大范围选线。ArcGIS地理信息***拥有强大的空间分析能力和海量数据处理能力,能方便快捷地处理大范围的地形地貌,为大范围选线提供强有力的支撑体系。改进的PSO算法能依据线路必经地区将大范围选线问题分成多个相互关联的段落,分别进行线路搜索,并在程序设计时利用多线程等技术同时搜索,提高搜索效率。
附图说明
图1为RPSO整数算法流程图;
图2为纵断面上的变坡点在平面上的投影示意图;
图3为本发明的应用示意图;
图4为本发明的线路寻优过程示意图;
图5为实例的地理环境信息示意图;
图6为自动选线参数设置示意图;
图7为线形优化方案群示意图;
图8参数对比分析图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明:
本发明一种基于改进型粒子群算法的道路三维线形自动优化方法,首先建立道路三维线形优化目标函数,依据设定线路平面及纵断面的约束条件、道路建设单价、线路取点间隔及线形优化目标区域,采用改进型粒子群算法对利用ArcGIS软件获得导向线方案的三维线形的平面交点及纵断面变坡点参数同时进行优化,获得自动优化后的道路三维线形;
所述三维线形的平面交点及纵断面变坡点参数是指平面交点的平面坐标(x,y)、半径R及缓长l,纵断面变坡点里程mil及高程z;
三维线形上所有平面交点及纵断面变坡点参数作为改进型粒子群算法中的粒子K;
所述改进型粒子群算法模型为:
其中,vkd t为粒子k第d维在第t次迭代的速度分量;ω为权重系数设定值;xkd t为粒子k第d维当前的位置分量;t为迭代的次数;int为向下取整函数,原型为int(number),取小于或等于number的最大的整数;ceil为向上取整函数,原型为ceil(number),取大于或者等于number的最小整数;
为整数区间内均匀分布的随机数:
当时,
当时,
为整数区间内均匀分布的随机数:
当时,
当时,
c1为控制粒子在飞向自身最好位置方向的步长;c2为控制粒子在飞向全局最好位置方向的步长;rand1、rand2是两个介于[0,1]范围的随机数;
pbestkd t表示粒子k在第t次迭代之前的历史最好位置;gbestkd t表示种群在第t次迭代之前的全局历史最好位置;
c1取2.0,c2取2.0。
改进型粒子群算法流程图如图1所示;
本实例的原始数据是内蒙古自治区乌兰察布市某地区5公里区段内的地形以及房屋、河流、道路、用地等地物信息资料。资料收集中部分数据可直接应用既有的AutoCAD图形文件及数据,并基于ArcGIS9.3软件环境完成资料的收集和处理。地形及地物的相关属性信息则依据实地收集资料输入***,建立了选线区域的信息数据库——数字高程模型、房屋、河流、道路以及用地等信息数据,从而得到了区域内智能选线的原始资料,如图5所示,图中所示***基于ArcGIS9.3软件环境开发及应用。
本次选线是在选线区域内应用本发明设计一条长度约5km的三级公路。由图5可看出本选线区域为山岭重丘区,地形丘陵起伏、沟壑纵横,无不良地质状况。根据具体情况设置参数:设计时速为60km/h;最小曲线半径为200m;曲线间最小直线长度为360m;平曲线最小长度300m;最大纵坡差为10%;土石方单价为20元/m3;桥梁单价为30000元/m2;隧道单价为25000元/m2,具体参数见图6。
如图3和图4所示为本实例应用本发明的示意图。
步骤1:建立道路三维线形目标函数F:
式中,i为横断面序号;b为桥梁的序号;t为隧道的序号;j为土地类型序号;IK为线形设计区间中横断面的个数;TK为线形设计区间中隧道座数;BK为线形设计区间中桥梁座数;YK为用地类型总个数;C为土石方工程单价;Cb为第b座桥梁的施工单价;Ct为第t座隧道的施工单价;Cj为第j种类型土地的用地单价;Si为第i个设计横断面的面积;Δw为两相邻横断面间的距离;Lb为第b座桥梁的长度;Lt为第t座隧道的长度;Aj为第j种土地类型的用地面积;
沿三维线形平面的法向量方向对线路进行剖分获得线路的横断面;
步骤2:设定三维线形平面及纵断面的约束条件、道路建设单价、三维线形平面取点间隔及三维线形优化目标区域xmin≤x≤xmax,ymin≤y≤ymax;
(1)平面约束条件
(a)最小曲线半径约束,平面交点Pi平曲线半径为r1i,不小于设定的最小曲线半径rmin,rmin≤r1i(i=1,2,...,n-1);
(b)最短圆曲线长度约束,两相邻缓和曲线间的圆曲线长度不小于设定的最短圆曲线长度
式中:为圆曲线长度;αi为圆曲线夹角;
(c)最短夹直线约束,平面夹直线最小长度为平面夹直线长度Li,i+1-Ti-Ti+1必须大于设定的平面夹直线最小长度为
式中:Li,i+1为平面交点i至交点i+l间的距离;Ti为平面交点i处的切线长;Ti+1为平面交点i+1处的切线长;
(2)纵断面约束条件
(a)坡长约束,纵断面设计坡长须满足最小坡长的限制;两相邻竖曲线的切线长之和小于该坡段长度,即:
其中,lcj-1,lcj分别为纵断面变坡点j-1及j的里程;分别为纵断面变坡点j-1及j处的切线长;
(b)坡度约束,纵断面第j段的设计坡度ij应该满足规范规定的允许坡度,即:
式中:hj,hj-1分别为纵断面变坡点j与j-1的设计高程;为纵断面坡段j的坡段长;imin,imax分别为规定的最小、最大坡度;
(c)相邻坡段坡度差约束,相邻坡段的坡度差的绝对值不超过允许的最大值Δimax,即
(d)控制点高程约束
线路在指定点处的设计高程hj应当满足最大最小高程(hmax,hmin)要求,即:
hmin≤hj≤hmax
(3)线形优化目标区域设定为xmin≤x≤xmax,ymin≤y≤ymax,同时满足以下条件:
1)禁止通过设定区域,即禁区;
设置禁区的地价费用或地面高程为无穷大,即从目标区域中排除禁区;
2)必须通过设定区域;
设平面必经区域中心点为O,计算O点距离线路中心线的距离do,设定较小值ε,并满足do≤ε,则线路可以从该区域通过,即该区域中必须包含有所选线路;
步骤3:从三维线形优化目标区域中,采用改进型粒子群算法寻找最优线形,适应度函数为F;
步骤a)采用由ArcGIS软件获得导向线方案的三维线形作为初始线形,初始线形上所有平面交点及纵断面变坡点参数作为粒子,设置粒子群种群大小、迭代次数及初始种群,利用改进型粒子群算法搜索更新初始线形上的所有交点的位置,输出最好解x1,即粒子位置构成的最佳三维线形;
步骤b)若步骤a)中搜索得到的三维线形中包含非三维线形优化目标区域中的位置,则采用静态搜索更新步骤a)得到的最好解x1在三维线形优化目标区域中的速度,利用速度更新最好解x1,获得第二最好解x2;
步骤c)以步骤b)获得的第二最好解x2作为初始探测点进行轴搜索,获得第三最好解x3;
步骤d)将步骤c)获得的第三最好解x3,输入ArcGIS软件,自动获得道路最优线形。
所述步骤3中步骤b)的静态搜索的搜索过程如下:
对于落入不可行区域的粒子,速度更新方程中,设定惯性权重系数ω=0,使得粒子新的速度与原飞行无关,只取决于粒子的当前位置、历史最好位置和粒子群的历史最好位置,即速度更新方程为:
vkd t+1=c1·rand1·(pbestkd t-xkd t)+c2·rand2·(gbestkd t-xkd t)
其中:vkd t为粒子k第d维在第t次迭代的速度分量;ω为权重系数设定值;xkd t为粒子k第d维当前的位置分量;t为迭代的次数;c1控制粒子在飞向自身最好位置方向的步长;c2控制粒子在飞向全局最好位置方向的步长;rand1、rand2是两个介于[0,1]范围的随机数;pbestkd t表示粒子k在第t次迭代之前的历史最好位置;gbestkd t表示种群在第t次迭代之前的全局历史最好位置;
c1取2.0,c2取2.0。
所述步骤3中步骤c)的轴搜索的搜索过程如下:
以第二最好解x2作为第一轮探测即初始探测点为x0,依次沿d维坐标轴的方向进行以下探测过程搜索得到多个探测点,从多个探测点中选择使得目标函数值F最小的点作为轴搜索得到的最好点;
设放大因子α>1,缩减因子β∈(-1,0),沿各轴方向的步长为δ1,δ2,...,δd;
每轮探测的起点和终点分别为y(1)和ye(即yn+1),令y(1)=x(k),x(k)为第k次转轴后的搜索起点,按如下规则进行探测;
y(i+1)=y(i)+δid(i)
从以x(k)为基点沿各个方向探测所得的点y(j)对应的目标函数值F中,选择一个目标函数值最小的对应点作为好解将好解与x(k)点对应的目标函数值F比较进行比较,按照下式更新步长:
判断是否满足以下任一条件,若满足,即完成探测,获得当前得到的好解并将作为第三最好解x3,否则进入下一轮探测:
(1)沿各轴方向不能找到比初始探测点x0的目标函数值更小的点y(j);
(2)各个方向的步长δj均达到设定的最小值ε,ε的取值范围是0.01-0.1。
图2为纵断面上的变坡点在平面上的投影示意图;
选线基本参数设置完成后,可以通过点取或录入的方式设置设计线路的起终点。然后,通过现有选线***,生成线路的导向线方案,如图7所示,其中展示了方案2、方案6、方案8和方案9的导向线走向。从导向线方案中选取方案8进行线形优化。
对方案8进行线形优化,配置合理的线形参数,经过采用本发明提供的方法,生成的线路的三维线形优化方案,如图8所示。
经三维线形优化后可以得到线路的平纵面数据信息。线路的平面交点数据如表2所示,纵断面变坡点数据如表3所示。
表2平面交点数据
表3纵断面变坡点数据
为了对***智能选线的运行效率进行测试,将本实例与文献[1]中的实例进行了对比分析,具体的运行结果数据对比见表4。
表4运行效率对比图
从上表可以看出,与文献[1]的实例规模相比,本实例的选线区域面积为其2倍,线路长度为其4倍。但本实例的运行时间比文献[1]缩短了近一倍,而且迭代次数较少,可以迅速收敛。所以与相对文献[1]比较,本文方法的运行效率较高。
Claims (4)
1.一种基于改进型粒子群算法的道路三维线形自动优化方法,其特征在于,首先建立道路三维线形优化目标函数,依据设定线路平面及纵断面的约束条件、道路建设单价、线路取点间隔及线形优化目标区域,采用改进型粒子群算法对利用ArcGIS软件获得导向线方案的三维线形的平面交点及纵断面变坡点参数同时进行优化,获得自动优化后的道路三维线形;
所述三维线形的平面交点及纵断面变坡点参数是指平面交点的平面坐标(x,y)、半径R及缓长l,纵断面变坡点里程mil及高程z;
三维线形上所有平面交点及纵断面变坡点参数作为改进型粒子群算法中的粒子k;
所述改进型粒子群算法模型为:
其中,vkd t为粒子k第d维在第t次迭代的速度分量;ω为权重系数设定值;xkd t为粒子k第d维当前的位置分量;t为迭代的次数;int为向下取整函数,原型为int(number),取小于或等于number的最大的整数;ceil为向上取整函数,原型为ceil(number),取大于或者等于number的最小整数;
为整数区间内均匀分布的随机数,当时, 当时,
为整数区间内均匀分布的随机数,当时, 当时,
c1为控制粒子在飞向自身最好位置方向的步长;c2为控制粒子在飞向全局最好位置方向的步长;
pbestkd t表示粒子k在第t次迭代之前的历史最好位置;gbestkd t表示种群在第t次迭代之前的全局历史最好位置;
c1取2.0,c2取2.0。
2.根据权利要求1所述的基于改进型粒子群算法的道路三维线形自动优化方法,其特征在于,包括以下步骤:
步骤1:建立道路三维线形目标函数F:
式中,i为横断面序号;b为桥梁的序号;t为隧道的序号;j为土地类型序号;IK为线形设计区间中横断面的个数;TK为线形设计区间中隧道座数;BK为设计区间中桥梁座数;YK为用地类型总个数;C为土石方工程单价;Cb为第b座桥梁的施工单价;Ct为第t座隧道的施工单价;Cj为第j种类型土地的用地单价;Si为第i个设计横断面的面积;Δw为两相邻横断面间的距离;Lb为第b座桥梁的长度;Lt为第t座隧道的长度;Aj为第j种土地类型的用地面积;
沿三维线形平面的法向量方向对线路进行剖分获得线路的横断面;
步骤2:设定三维线形平面及纵断面的约束条件、道路建设单价、三维线形平面取点间隔及三维线形优化目标区域xmin≤x≤xmax,ymin≤y≤ymax;
(1)平面约束条件
(a)最小曲线半径约束,平面交点Pi平曲线半径为r1i,不小于设定的最小曲线半径rmin,rmin≤r1i(i=1,2,...,n-1);
(b)最短圆曲线长度约束,两相邻缓和曲线间的圆曲线长度不小于设定的最短圆曲线长度
式中:为圆曲线长度;αi为圆曲线夹角;
(c)最短夹直线约束,平面夹直线最小长度为平面夹直线长度Li,i+1-Ti-Ti+1必须大于设定的平面夹直线最小长度为
式中:Li,i+1为平面交点i至交点i+l间的距离;Ti为平面交点i处的切线长;Ti+1为平面交点i+1处的切线长;
(2)纵断面约束条件
(a)坡长约束,纵断面设计坡段长度须满足最小坡长的限制;两相邻竖曲线的切线长之和小于纵断面设计坡段长度,即:
其中,lcj-1,lcj分别为纵断面变坡点j-1及j的里程;分别为纵断面变坡点j-1及j处的切线长;
(b)坡度约束,纵断面第j段的设计坡度ij应该满足规范规定的允许坡度,即:
式中:hj,hj-1分别为纵断面变坡点j与j-1的设计高程;为纵断面坡段j的坡段长;imin,imax分别为规定的最小、最大坡度;
(c)相邻坡段坡度差约束,相邻坡段的坡度差的绝对值不超过允许的最大值Δimax,即
(d)控制点高程约束
线路在指定点处的设计高程hj应当满足最大最小高程(hmax,hmin)要求,即:
hmin≤hj≤hmax
(3)线形优化目标区域设定为xmin≤x≤xmax,ymin≤y≤ymax,同时满足以下条件:
1)禁止通过设定区域,即禁区;
设置禁区的地价费用或地面高程为无穷大,即从目标区域中排除禁区;
2)必须通过设定区域;
设平面必经区域中心点为O,计算O点距离线路中心线的距离do,设定较小值ε,并满足do≤ε,则线路可以从该区域通过,即该区域中必须包含有所选线路;
步骤3:从三维线形优化目标区域中,采用改进型粒子群算法寻找最优线形,适应度函数为F;
步骤a)采用由ArcGIS软件获得导向线方案的三维线形作为初始线形,初始线形上所有平面交点及纵断面变坡点参数作为粒子,设置粒子群种群大小、迭代次数及初始种群,利用改进型粒子群算法搜索更新初始线形上的所有交点的位置,输出最好解x1,即粒子位置构成的最佳三维线形;
步骤b)若步骤a)中搜索得到的三维线形中包含非三维线形优化目标区域中的位置,则采用静态搜索更新步骤a)得到的最好解x1在三维线形优化目标区域中的速度,利用速度更新最好解x1,获得第二最好解x2;
步骤c)以步骤b)获得的第二最好解x2作为初始探测点进行轴搜索,获得第三最好解x3;
步骤d)将步骤c)获得的第三最好解x3输入ArcGIS软件,自动获得道路最优线形。
3.根据权利要求2所述基于改进型粒子群算法的道路三维线形自动优化方法,其特征在于,所述步骤3中步骤b)的静态搜索的搜索过程如下:
对于落入不可行区域的粒子,速度更新方程中,设定惯性权重系数ω=0,使得粒子新的速度与原飞行无关,只取决于粒子的当前位置、历史最好位置和粒子群的历史最好位置,即速度更新方程为:
vkd t+1=c1·rand1·(pbestkd t-xkd t)+c2·rand2·(gbestkd t-xkd t) (16)
其中:vkd t为粒子k第d维在第t次迭代的速度分量;xkd t为粒子k第d维当前的位置分量;t为迭代的次数;c1控制粒子在飞向自身最好位置方向的步长;c2控制粒子在飞向全局最好位置方向的步长;rand1、rand2是两个介于[0,1]范围的随机数;pbestkd t表示粒子k在第t次迭代之前的历史最好位置;gbestkd t表示种群在第t次迭代之前的全局历史最好位置;
c1取2.0,c2取2.0。
4.根据权利要求2所述基于改进型粒子群算法的道路三维线形自动优化方法,其特征在于,所述步骤3中步骤c)的轴搜索的搜索过程如下:
以第二最好解x2作为第一轮探测即初始探测点为x0,依次沿d维坐标轴的方向进行以下探测过程搜索得到多个探测点,从多个探测点中选择使得目标函数值F最小的点作为轴搜索得到的最好点;
设放大因子α>1,缩减因子β∈(-1,0),沿各轴方向的步长为δ1,δ2,...,δd;
每轮探测的起点和终点分别为y(1)和ye,令y(1)=x(k),x(k)为第k次转轴后的搜索起点,d(i)表示d维坐标轴中第i个坐标轴的方向,按如下规则进行探测;
y(i+1)=y(i)+δid(i)
从以x(k)为基点沿各个方向探测所得的点y(j)对应的目标函数值F中,选择一个目标函数值最小的对应点作为好解将好解与x(k)点对应的目标函数值F比较进行比较,按照下式更新步长:
判断是否满足以下任一条件,若满足,即完成探测,获得当前得到的好解并将作为第三最好解x3,否则进入下一轮探测:
(1)沿各轴方向不能找到比初始探测点x0的目标函数值更小的点y(j);
(2)各个方向的步长δj均达到设定的最小值ε,ε的取值范围是0.01-0.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310252969.3A CN103399849B (zh) | 2013-06-24 | 2013-06-24 | 一种基于改进型粒子群算法的道路三维线形自动优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310252969.3A CN103399849B (zh) | 2013-06-24 | 2013-06-24 | 一种基于改进型粒子群算法的道路三维线形自动优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103399849A CN103399849A (zh) | 2013-11-20 |
CN103399849B true CN103399849B (zh) | 2017-04-12 |
Family
ID=49563480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310252969.3A Expired - Fee Related CN103399849B (zh) | 2013-06-24 | 2013-06-24 | 一种基于改进型粒子群算法的道路三维线形自动优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103399849B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106126834B (zh) * | 2016-06-28 | 2019-07-19 | 西安交通大学 | 一种基于惯导角度量测的轨道平面最佳线形确定方法 |
CN107248006B (zh) * | 2017-05-27 | 2020-06-16 | 北方工业大学 | 基于分层递阶的地铁线路客流协调控制方法 |
CN107632616B (zh) * | 2017-09-05 | 2019-08-09 | 华中科技大学 | 一种基于三维空间曲线的无人机协同路径规划方法 |
CN108734760A (zh) * | 2017-09-12 | 2018-11-02 | 武汉天际航信息科技股份有限公司 | 三维网格线条拉直***及三维网格线条拉直方法 |
CN107704678B (zh) * | 2017-09-29 | 2022-03-15 | 天津师范大学 | 一种脉冲成形滤波器的滚降系数的最优选择方法 |
CN107743104B (zh) * | 2017-09-29 | 2021-08-31 | 天津师范大学 | 基于扫描和缩放的广义频分复用技术中脉冲成形滤波器的滚降系数选择方法 |
CN107729663A (zh) * | 2017-10-23 | 2018-02-23 | 陈国� | 一种基于多层约束体系的公路路线三维优化方法 |
CA3093667C (en) | 2018-05-08 | 2023-02-21 | Landmark Graphics Corporation | System and method for application of elastic property constraints to petro-elastic subsurface reservoir modeling |
CN109165427B (zh) * | 2018-08-06 | 2022-11-04 | 中国铁路设计集团有限公司 | 运营高铁纵断面线形调整方案的优化计算方法 |
CN109447437B (zh) * | 2018-10-17 | 2021-10-29 | 中南大学 | 一种包含土石方调配的公(铁)路纵断面自动构建方法 |
CN113204856B (zh) * | 2020-05-14 | 2022-03-18 | 中国铁路设计集团有限公司 | 运营高铁轨道空间线形优化设计方法 |
CN111709079B (zh) * | 2020-06-16 | 2022-02-11 | 中交第四航务工程勘察设计院有限公司 | 一种悬浮隧道参数化横断面几何优化形状及参数寻优算法 |
CN112883474B (zh) * | 2021-03-02 | 2021-09-24 | 中南大学 | 一种既有轨道线路智能改建的布设方法、***、终端及可读存储介质 |
CN113554467B (zh) * | 2021-07-26 | 2022-07-22 | 中南大学 | 一种基于协同进化的铁路三维线形智能设计方法 |
CN114154382B (zh) * | 2021-11-18 | 2022-11-15 | 中国公路工程咨询集团有限公司 | 高分遥感三维可视化铁路选线辅助决策***及方法 |
CN115113390B (zh) * | 2022-07-22 | 2024-05-17 | 中国科学院光电技术研究所 | 一种基于改进粒子群算法的光束整形方法 |
CN116305678B (zh) * | 2022-12-08 | 2023-12-15 | 西南交通大学 | 基于低碳排放的高密城区铁路线路精细优化方法 |
CN116433860A (zh) * | 2023-03-27 | 2023-07-14 | 腾讯科技(深圳)有限公司 | 路口面生成方法、装置、设备、存储介质和程序产品 |
CN117708961B (zh) * | 2024-02-05 | 2024-04-30 | 中南大学 | 一种既有铁路三维空间线位整体智能重构方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102752761A (zh) * | 2012-06-19 | 2012-10-24 | 江苏科技大学 | 基于粒子群的无线传感器网络移动节点覆盖优化方法 |
-
2013
- 2013-06-24 CN CN201310252969.3A patent/CN103399849B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102752761A (zh) * | 2012-06-19 | 2012-10-24 | 江苏科技大学 | 基于粒子群的无线传感器网络移动节点覆盖优化方法 |
Non-Patent Citations (3)
Title |
---|
公路平面线形交互设计***的研究;繆鹍;《工程设计CAD与智能建筑》;20011231(第12期);全文 * |
基于改进粒子群优化算法的三维空间路径规划;陈家照 等;《计算机应用与软件》;20101130;第27卷(第11期);全文 * |
基于线元的公路平面线形交互设计方法研究;繆鹍等;《中国公路学报》;20010731;第14卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103399849A (zh) | 2013-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103399849B (zh) | 一种基于改进型粒子群算法的道路三维线形自动优化方法 | |
CN103996089B (zh) | 基于gis的输电线路最优路径生成方法 | |
Zhao et al. | Truck traffic speed prediction under non-recurrent congestion: Based on optimized deep learning algorithms and GPS data | |
Jong et al. | An evolutionary model for simultaneously optimizing three-dimensional highway alignments | |
Pu et al. | A three‐dimensional distance transform for optimizing constrained mountain railway alignments | |
CN103971184B (zh) | 基于gis空间地理信息***的输电线路路径生成方法 | |
CN106649821A (zh) | 空间目标索引构建、碰撞预警、区域和最近邻查询方法 | |
CN112734176B (zh) | 一种充电站建站方法、装置、终端设备和可读存储介质 | |
Maji et al. | Multi‐objective highway alignment optimization using a genetic algorithm | |
CN106197460A (zh) | 一种应用gps出行数据进行出行目的地预测的方法 | |
CN106127855B (zh) | 一种海量路网实时精细化三维模型的渲染方法 | |
CN105678034A (zh) | 行人过街时间模型及信号交叉口人行横道宽度优化方法 | |
CN104406590B (zh) | 一种基于道路等级的最短路径规划方法 | |
CN105374204A (zh) | 一种城市道路交通检测器布点的方法 | |
CN105067004A (zh) | 一种基于地形的路径搜索方法 | |
An et al. | Research on the construction of carbon emission evaluation system of low-carbon-oriented urban planning scheme: taking the West New District of Jinan city as example | |
CN104282142A (zh) | 一种基于出租车gps数据的公交站台设置方法 | |
Megarestya et al. | The Horizontal Curved Geometric Design with Autocad® Civil 3D on Jalan Muara Wahau, East Kalimantan | |
CN108122185B (zh) | 考虑布控区域边界点的移动污染源排放遥测站点选址方法 | |
CN104331424B (zh) | 一种用户移动轨迹中敏感轨迹模式的净化方法 | |
CN103870631A (zh) | 一种基于3s技术的智能输电网络布设模型构建方法 | |
CN105887593B (zh) | 多年冻土地区公路三维导向线选定方法 | |
CN109800903A (zh) | 一种基于出租车轨迹数据的盈利路线规划方法 | |
CN104463442A (zh) | 一种城乡建设集聚性的探测方法 | |
Wahdan et al. | Design an optimum highway route using remote sensing data and GIS-based least cost path model, case of Minya-Ras Ghareb and Minya-Wahat-Bawiti highway routes, Egypt |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170412 Termination date: 20180624 |
|
CF01 | Termination of patent right due to non-payment of annual fee |