基于粒子群算法的改建道路纵断面自动设计研究_董华珍.pdf

摘要：为改变现有改建道路纵断面设计软件拉坡设计人工耗时长、随意性大、主观性强、设计方案参差不齐的状况，首次提出基于粒子群算法的改建道路纵断面自动设计，以提高设计效率、减少工程投资。在满足相关规范、路面改造方案前提下，以改造纵断面填挖方高度总和最小为目标函数，以纵坡长度、纵向坡度、竖曲线半径为约束条件，以变坡点的里程、高程为变量，利用粒子群算法原理，构建了基于粒子群算法的改建道路纵断面自动设计模型。研究结果表明：该设计模型达到了改建道路纵断面自动适应贴合现状道路高程、路面改造工程量最小的目的，当粒子群算法迭代次数达到 400 多次时，粒子适应度函数趋于收敛，当迭代次数设为 500 次时，可以获得满意的成果。实例验证表明：与传统的人工设计相比，算法自动设计的效率更高，自动设计过程仅需 23 s，拟合原地面线的效果更好，能节省近一半的路面加铺工程量，优化设计效果显著。基于粒子群算法的改建道路纵断面自动设计，可为道路改造设计及设计咨询提供理论基础与技术支持。关键词：道路工程；纵断面自动设计；粒子群算法；改建道路中图分类号：U412文献标志码：A文章编号：1009-7716（2023）01-0205-05基于粒子群算法的改建道路纵断面自动设计研究收稿日期：2022-03-01作者简介：董华珍（1981），女，硕士，高级工程师，从事市政路桥设计与咨询工作。董华珍1，聂 涔2（1.广州地铁设计院施工图咨询有限公司，广东 广州 510010；2.广州地铁设计研究院股份有限公司，广东 广州 510010）DOI:10.16799/ki.csdqyfh.2023.01.0510引言随着我国国民经济的持续快速发展、城市化进程的加快，机动化程度得到很快的发展，私家车己逐渐进入普通市民的家庭当中。根据公安部发布的统计数据，截至 2021 年 6 月，我国机动车保有量达到2.92 亿辆，其中千人机动车保有量已从 2002 年的62 辆1增至 2021 年的 200 辆，增幅超过 倍。机动车保有量的激增使得原有的城市道路满足不了人们对城市道路的需求。近年来，城市道路交通已经从迅猛发展的高潮期逐步过渡到改造维护的存量时代。随着城市土地资源的越发紧张，旧路改造已经成为改善现状道路交通条件的关键手段。在道路改造时，既有道路的平面线形基本定格，受限于道路红线的约束，道路改造一般情况下会保持原有的平面线形不变。受常年累月的交通荷载、自然环境、路面结构病害等因素的影响，道路纵断面往往发生较大的沉降变化，导致竖向行车舒适性大大降低。因此，纵断面设计是道路改造较关键的环节。从 20 世纪 60 年代开始，计算机技术开始应用于公路设计中。经过多年的发展，形成了诸多成熟的市政道路设计软件，如 HY-SZDL、EICAD、RDCADS、纬地道路、OpenRoads Designer 等。市场上应用的这些道路设计软件，仅是提供人机交互操作平台，软件的纵断面设计方法由设计师在软件界面上手动拉坡、反复调整拉坡方案，需要耗费大量的设计精力处理复杂的设计条件。同时还受限于设计人员的技术水平和经验的主观性，往往难以取得理想的设计方案。旧路改造工程要求纵断面线形与现状道路的地面线尽量贴合，保证纵断面填挖方工程量尽量小，传统的人机交互拉坡设计方法难以满足设计的需要。将复杂的设计条件抽象成数学模型，利用人工智能算法进行优化求解，已经成为新的设计潮流。目前关于纵断面自动设计的研究，多集中在铁路与公路方面2-4。铁路与公路一般的纵断面受限于地面的起伏情况，一般情况下纵断面的起伏较大，属于比较粗犷的纵断面设计类型。而市政道路改造纵断面波动较小，属于精细化方面的设计，然而目前关于这方面的研究很少，而且基本上是传统的数值解法，不够智能化、自动化5-6。笔者通过分析改建道路纵断面设计要点和现有道路设计软件纵断面设计功能局限性，首次将粒子群算法用于城市道路改建当中，以道路路面改造工程量最小为目标，建立了基于粒子群算法道路改造纵断面自动设计模型。该模型可大大提升设计效率，提高设计质量，减少工程投资。1数学模型道路纵断面是由一系列的变坡点组成的，抽象成数学参数为变坡点里程、高程和竖曲线半径。其中，变URBAN ROADS BRIDGES&FLOOD CONTROL第 1 期（总第 285 期）2023 年 1 月科技研究205坡点里程和高程是确定纵断面位置和坡度的关键因素。在道路改造的纵断面设计中，竖曲线半径起到平顺过渡坡度的作用，其取值一般受限于保证竖曲线半径、竖曲线长度满足规范最低要求即可，整体上起到一个缓和与辅助的作用。因此，本文不将竖曲线半径作为优化参数，后续在满足规划约束的条件下，通过变坡点的相邻坡度差来适应和匹配竖曲线半径。因此，将每个变坡点的里程 mi1、高程 hi1变量编码成粒子群向量，如下：Xi=（mi1，hi1，mi2，hi2，min，hin）1.1目标函数城市道路改造工程的纵断面设计思路，主要是根据路面结构改造的具体方案，纵坡的设计尽量贴合拟改造道路的现状高程，尽量减少填挖方的工程量，以节省工程投资5。改造的纵断面设计抽象成数学模型，即求解多变量函数的条件极值问题。也就是在满足道路改造所遵循的规范与标准等前提下，得到工程改造量最小的设计方案。实际工程中，计算工程总费用较为复杂，需要考虑土石方、桥涵、隧道、横断面形式等因素。但在道路改造工程中，主要的工程量都是集中在道路纵断面中线上的填挖方总和。因此，为了简化数学模型，将目标函数定义为：纵断面中线上的填挖总和。其表达式定义如下：S=x0g（x）（x）-h（x）dx式中：f（x）为纵断面设计标高；h（x）为纵断面地面标高；g（x）为权函数，针对不同的横断面时，可以调整不同的断面影响值。1.2约束条件1.2.1坡长约束城市道路路线设计规范对最小坡长与最大坡长进行了规定，旧路的一般路段设计坡长应能满足规范要求的最小坡长限制，且应考虑路线尽端道路起（讫）点一端可不受最小坡长限制。一般情况下，设计速度不变，旧路的原设计最大坡长能满足规范要求，故暂不考虑最大坡长约束。因此，本模型坡长约束考虑最小坡长满足规范要求，即 LLmin，且起终点及与旧路交叉口路段不受最小坡长限制。1.2.2坡度约束为满足竖向排水的要求，规范规定道路最小纵坡为 0.3%。但在实际改造的情况中，受到施工精度、路基沉降等影响，道路实际坡度难以避免会有小于0.3%的情况。对于最小纵坡小于 0.3%的情形，按规范要求设置相应排水措施。因此，最小纵坡不做要求。道路的最大纵坡应满足道路设计规范的要求，即：iimax。1.2.3竖曲线半径为保证行车视距合理，竖向离心加速度舒适可行，竖曲线半径应满足规范规定，即：RRmin。一般竖曲线长度需保证车辆 3 s 行程的要求，有 IImin。1.2.4平纵组合城市道路路线设计规范 对设计速度不小于60 km/h 的道路强调线形组合设计，保证指标均衡、视觉良好、安全舒适；对设计速度小于 60 km/h 的道路的线形组合则没有硬性要求。故本研究暂不考虑对设计速度小于 60 km/h 的道路的线形组合的约束。1.3粒子群算法求解粒子群算法（particle swarm optimization,PSO）是一种源于鸟群觅食行为的算法7。该算法具有计算速度快、并行性高、收敛速度快、参数少、稳定性强等特点，在工程领域内多变量求优化解上有广泛的应用8-9。算法中将每个粒子组编成一个 n 维向量，每个粒子根据自己当前的位置和速度，参考群体中的经验，每时每刻都进行动态调整，通过不断的进化迭代，达到群体目标最优的状态。具体的变量编码形式如下：Xi=（Xi1，Xi2，Xi3，Xin）表示粒子 i 当前的位置；Vi=（Vi1，Vi2，Vi3，Vin）表示粒子 i 当前的速度；Pi=（Pi1，Pi2，Pi3，Pin）表示是粒子 i 的个体最好位置；Pg=（Pg1，Pg2，Pg3，Pgn）表示整个群体中的最好位置。粒子进化的方程如下：Vij（t+1）=wVij（t）+c1r1Pij（t）-Xij（t）+c2r2Pgj（t）-Xij（t）Xij（t+1）=Xij（t）+Vij（t）式中：Vij（t+1）为粒子 i 中第 j 个变量在 t+1 时刻的速度；Xij（t）为粒子 i 中第 j 个变量在 t 时刻的位置；Pij（t）为粒子 i 中第 j 个变量在 t 时刻的最优解；Pgj（t）为整个种群中第 j 个变量在 t 时刻最优解；w 为速度变化权重；c1，c2为学习因子；r1，r2为（0，1）区间内均匀分布的随机数。同时，为了保证粒子不超越搜索空间，通常给粒子的速度设定一个变化范围。Vij（-Vmax，Vmax）董华珍，等：基于粒子群算法的改建道路纵断面自动设计研究2023 年第 1 期206基于粒子群算法求解的具体步骤如下。Step1：种群初始化以原始的设计方案为基础，将每个变坡点的里程 mi1、高程 hi1变量编码成粒子群向量，如下：Xi=（mi1，hi1，mi2，hi2，min，hin）竖曲线半径取值在满足规范的要求下，根据相邻坡段的坡差进行确定，一般情况下取整千的数值。具体确定方式如下：R=max（1000ceil（lmin/（1000i），Rmin）式中：lmin表示最小竖曲线长度；i 表示变坡点 i 的相邻坡差的绝对值；Rmin表示最小竖曲线半径；ceil（）表示向上取整的函数。对初始线路表示的纵断面向量分别进行随机扰动，形成多条新的纵断面方案。Step2：初始化适应度计算分别计算 step1 中生成的每个粒子的适应度。为了进一步简化目标函数，采用数学中微元求和的思想，在纵断面上每隔 1 m 取一个点，进行纵断面中线上填挖方总量的累积求和。具体适应度函数的计算方法如下：minS=ni=1g（i）f（i）-h（i）-M式中：f（i）为里程 i 处的纵断面设计高程；h（i）为里程i 处的地面线高程；g（i）为里程 i 处权函数，若道路横断面形式一致，则取 g（i）=1；M 为约束惩罚项，若粒子不满足相应的约束条件，则 M 取一个很大的值，若粒子满足规范约束，则 M=0。将初始的粒子群中的各个粒子的评价值作为各个粒子的个体历史最优解 j，并寻找各子群中的最优解 PL和总群体中的最优解 Pg。Step3：迭代优化。不断迭代以下步骤，直到满足收敛条件或达到最大迭代次数。（1）对每一个粒子，按照粒子进化方程更新。（2）计算每一个粒子的适应度。（3）若粒子的当前适应度优于其历史最优适应度，则记该粒子的当前适应度为该粒子的历史个体最优适应度。同时记当前位置为该粒子历史个体最优位置。（4）寻找当前各子群中的最优解和总群体中的最优解，若分别优于各自的历史最优解，则更新 PL、Pg。算法的终止条件：当算法进行迭代 N 次后，适应度函数值趋于平缓，最终收敛于一个较优值时，即 N代之后的个体最优解即本次优化的纵断面设计最优方案。经多次试验验证，当迭代次数达到 400 多次时，适应度函数趋于收敛；设定迭代次数为 500 次时，可以获得满意的成果。2实例验证及分析为验证粒子群算法在道路改造纵断面优化设计中的可行性和有效性，本文进行了计算机编程求解验证。笔者选取两个经施工图审查合格的主干道改造项目（设计速度均为 60 km/h）进行优化测试，并将优化方案（即基于粒子群算法的纵断面自动设计方案）与审查合格方案（以下称为“原始方案”）进行对比。根据 城市道路路线设计规范 要求，60 km/h 设计速度的道路最小纵坡长度 150 m，凹形竖曲线最小半径 1 500 m，凸形竖曲线最小半径 1 800 m，竖曲线最小长度 120m（现状路接顺段不受此限）。粒子群的参数设置：粒子数 n=20，c1=c2=1.5，最大迭代次数 N=500，对应变坡点里程的粒子最大速度为 5，对应变坡点高程的粒子最大速度为 0.1。2.1实例验证实例 1：某城市某旧路改造项目，道路等级为城市主干路，设计车速 60 km/h，道路北起工业大道，南至金斗大桥，全长 2 500 m。道路红线宽 50 m，车行道宽 26 m