改进粒子群算法在系杆拱桥成桥索力优化中的应用

1、改进粒子群算法在系杆拱桥成桥索力优化中的应用系杆拱桥是一种基于拱肋、吊杆、系杆，外部静定内部超静定的三元结构，因其整体刚度较大，力学性能较好、施工便捷、造价低等特点被国内外广泛应用。系杆拱桥以合理成桥状态为设计目标，即在恒载作用下系杆、拱肋的受力合理、主桥线形平顺，这就需要吊杆具有合适的成桥索力。因此，系杆拱桥的优化问题便可转化成吊杆索力的优化问题。目前系杆拱桥吊杆的索力求解方法主要沿用一些斜拉桥上的索力求解方法，Leonhardt率先提出了刚性支撑连续梁法并最早应用于美国P-K桥的分析中；肖汝诚等1，2结合影响矩阵推导了弯曲最小能量法的公式；梁鹏等3提出通过调整结构刚度来近似求解最小弯曲能量

2、法。传统的调索方法因受限于各自的局限性而不具有普遍适用性，如刚性支撑连续梁法不能考虑多条件的约束问题，因而用其计算的索力值往往不合理；弯曲最小能量法在结合影响矩阵情况下，想要精确求解计算会十分繁琐，而采用近似求解方法也需要不断调节结构的刚度才能使其达到预期结果。为了使索力求解的结果更加准确，很多学者将数学方法与桥梁模型相结合4，5，然而通过数学方法所得的目标函数往往十分复杂而不易求解。为了避免传统调索方法的缺点，且能够快速有效地求解出索力，有学者提出将计算机智能优化算法与索力结合的方法，陈志军等6、吴霄等7分别采用传统的粒子群和遗传算法对斜拉桥索力进行了优化，但传统的优化算法因其自身的局限性导

3、致结果往往不够理想。为了使结果更加准确，本文以实际工程项目为依托，通过对粒子群优化算法加以改进并与影响矩阵相结合，在兼顾计算效率和多目标约束的条件下求解吊杆索力，并与其他多种方法进行分析比较，从而验证该方法的可行性。1 基于改进粒子群算法的自动调索1.1理论基础粒子群算法最早是由Kennedy和Eberhart8提出的一种基于群鸟觅食的仿生智能算法，其原理是群鸟在觅食过程中，每只鸟的初始位置和飞行方向都处于随机的状态，而且不知道最佳的觅食点在何处。鸟群通过相互学习、信息共享，并在每次觅食过程中结合自身经验和种群之间的信息传递不断调整速度和位置，最终步步逼近食物。粒子寻优示意图如图1所示。图1粒

4、子寻优示意图Fig.1Schematic diagram of the particle swarm optimization如果将该方法运用到实际应用中，则是把鸟群抽象为m个没有质量和体积的粒子，再将其延伸至N维空间，粒子在N维空间中的位置表示为Xi=(xi1,xi2,xij)，速度表示为Vi=(vi1,vi2,vij)，且其位置和速度需限制在Xmin,Xmax和Vmin,Vmax范围内。每个粒子都有一个由目标函数决定的适应度值fitness，并且知道自己迄今为止发现的最好位置pb=(pi1,pi2,pij)和群体发现的最好位置gb=(g1,g2,gj)。粒子通过追踪pb和gb来更新自己的速

5、度和位置。速度和位置更新公式为vk+1ij=wvkij+c1r1pkijxkij+c2r2gkixkij（1）xk+1ij=xkij+vk+1ij（2）式中：i表示粒子的个数，i=1，2，m；j表示粒子的维度，j=1，2，N；k是迭代步数； w是惯性因子；c1、c2是学习因子；r1、r2是0，1范围内的均匀随机数；vmax为用户自定义常数，一般取vmax=aXmax，vmin=aXmin，通常a为0.10.2。1.2算法优化由于传统粒子群算法的程序实现过程十分简单，且需要调整的参数较少，因此在随机优化算法中具有明显的优势，并被广泛应用于工程领域。但同时它也存在收敛速度慢、局部搜索能力差等缺点。

6、为了提高算法的适用范围，使其结果更加准确可靠，需要对传统的算法进行改进9-11。研究结果表明12-16，w较大时算法的全局搜素能力较强，w较小时算法的局部搜索能力较强，因此算法的优化问题便可转化成一个w值的选取问题，即当粒子目标值趋于局部最优时，需要增大惯性因子，当粒子目标值比较分散时，则需减小惯性因子。为了合理选取w，本文将权重进行自适应调整，通过判断粒子当前的目标函数值U的好坏来自动调整w：当UUavg时，粒子的目标函数值要比平均目标值差，为了使其向较好的搜索区靠拢，对应的惯性因子要较大。当UavgUUavg时，粒子的全局寻优能力和局部寻优能力都较好，故不改变惯性因子。当UUavg时，粒子

7、的目标函数值要优于平均目标值，因此选用较小的惯性因子来保护该粒子，本文改进后w的表达式可表示为w=w(wwmin)UUavgUminUavg,UUavgw,UavgUUavgwmax,UavgU（3）式中：wmax、wmin分别表示w的最大值和最小值；Umin、Uavg分别表示当前所有粒子的最小适应度值和平均适应度值。由于本文是以弯曲能量的最小值为目标，因此Uavg是所有小于Uavg的适应度值的平均值。1.3目标函数目标函数下每一个粒子都有一个适应度值，它作为每次迭代的临时储存变量，并不具有记忆性，优化算法可通过粒子每次迭代的适应度值来评价结果的好坏，进而不断寻优，最终靠近目标。实践证明，通过

8、限制拱桥结构的弯曲能量所求解的吊杆索力能使桥梁结构处于一个较好的受力状态，因此本文以吊杆拱桥的弯曲能量为目标，并将结构离散化，离散后的弯曲应变能公式为U=i=1mLi2EiIiM2i（4）式中：U为结构的弯曲应变能；m为离散单元的数量；Li、Ei、Ii、Mi分别为单元的长度、材料的弹性模量、截面惯性矩和弯矩。为了能够建立起设计变量（拱桥的吊杆索力）与优化目标间的函数关系，可利用影响矩阵原理作为连接二者的桥梁，进而达到索力优化的目的。影响矩阵公式具体表示为AY=D（5）式中：Y为施调向量，在本文中是指吊杆索力组成的列向量；D为受调向量，是指结构中关心截面上若干独立元素所组成的列向量，这些元素一般

9、指关心截面的内力、位移，通过将这些元素调整到期望状态来求解施调向量Y；A为影响矩阵，是指当施调向量的某一向量发生单位变化时引起的受调向量D的变化量。从式（4）可以看出，当结构确定时，其单元及材料的具体特性已经确定，即Li、Ei、Ii已经确定，因而只能通过改变截面弯矩Mi来调整应变能的大小。将截面的弯矩改变量作为受调向量，将结构调索前离散单元的弯矩组成的向量作为M0，将结构单元需要达到的弯矩作为目标弯矩M，由每根吊杆索力S组成的向量X作为施调向量，即X=S1,S2,SnT，则调索后的弯矩为M=M0+CX（6）将式（6）中M0移至等号左边，MM0便对应弯矩的改变量，对应式（5）中的D。由每根吊杆单

10、位索力引起的弯矩变化值组成影响矩阵C，X则对应式（5）中的施调向量项。再将式（6）代入式（4）后可得结构的目标函数为U=M0TBM0+M0TBCX+XTCTBM0+XTCTBCX（7）式中：B为系数矩阵，是对角阵，表示单元柔度对单元弯矩的加权矩阵，其元素bii=li/(2EiIi)。1.4条件约束对于系杆拱桥，其弯矩主要由系梁承受，因此必须对系梁弯矩加以约束以使其落在允许的范围内。弯矩的约束条件可表示为MminMMmax。根据影响矩阵理论，在对弯矩进行限制的同时也相当于对索力进行了约束，从而使调索后的最优解具有了合理性，即当MminMMmax时，可得索力约束条件为C1(M1minM01)S1C

11、1(M1maxM01)C1(M2minM02)S2C1(M2maxM02)C1(MjminM0j)SiC1(MjmaxM0j)（8）式中：Si表示第i根吊杆的索力；Mjmax、Mjmin分别代表系梁上第j个单元所允许的最大、最小弯矩；M0j表示系梁第j个单元的初始弯矩。同时为了满足索力均匀原则，需将相邻索力的差值限定在可控范围内，另对索力做以下约束：i=119(Si+1SiSi+1)2（9）式中：为索力均匀度的评价阈值。1.5优化调索路线本文将索力调整和计算机智能优化算法结合，并考虑了惯性权重的实时变化，使索力调整过程自动化，从而实现索力的高效寻优。其具体步骤如下：1）利用有限元软件建立拱桥模

12、型，并提取相关结构数据及影响矩阵作为基础数据库。2）建立目标函数，利用Matlab软件编程粒子群算法并搜索在约束范围内的最优解，如果满足停止条件，则输出当前解，否则更新权重继续搜索。3）将搜索值代回有限元模型，校核结果的有效性。优化后的调索路线具体计算流程如图2所示。图2计算流程图Fig.2Calculation flow chart2 工程应用2.1工程概况本文以某钢管混凝土拱桥为实际工程依托，主桥总长150 m，计算跨径146.28 m，桥梁总宽 12.9 m，拱轴线为二次抛物线，矢跨比为1/5。拱肋采用哑铃型钢管砼，采用Q345q钢材，内充C50微膨胀混凝土。系梁采用箱形断面，系梁高为2

13、40 cm，宽为120 cm，系梁和横梁为预应力混凝土结构，桥面2%横坡通过横梁高度的变化进行调整，吊杆间距为7.30 m，每片拱肋设吊杆19根，吊杆编号由小里程到大里程依次为119号，如图3所示。图3吊杆编号示意图Fig.3Schematic diagram of boom numbering2.2参数设置本文系杆拱桥采用有限元软件ANSYS建模，如图4所示。除吊杆采用桁架单元外，其余均采用梁单元。恒载作用下主梁和拱肋的弯矩值M0通过软件的数据提取功能来获取。通过分别改变每根吊杆的初拉力，使其发生单位1的变化，得到关键节点的弯矩变化值从而组成影响矩阵。模型的主要材料特性见表1。图4系杆拱有限

14、元模型Fig.4Finite element model of tied arch表1材料特性Table 1Material properties材料弹性模量/GPa泊松比线膨胀系数/(10-6/)容重/(kNm-3)C5034.50.21025.0Q345q20.60.31278.5Strand186019.50.31278.5利用Matlab工具编写改进前和改进后的粒子群算法程序并进行迭代求解，改进前后除惯性权重外，其他系数取值相同。目前w较典型的取值范围为0.20.917，在此范围内取值时算法局部和全局收敛能力较好。为使算法能在一开始便找到较好的解，需要取范围内较大值，因此本文对改进前的

15、惯性权重w恒取0.9，对改进后惯性权重系数按式（3）计算。学习因子c1、c2取1.49，由文献18可知，学习因子取1.49时算法的自我学习能力和社会学习能力均较好。种群规模m取50，索力范围取500，1 500，速度范围取索力范围的15%，具体数据见表2。由图5、6可以看出，在迭代次数都为200次的条件下，改进后的算法收敛速度得到明显提高，在收敛精度方面也有一定的改善。表2改进后粒子群算法参数Table 2Improved PSO parametersmaxmaxminminc1c1c2c2mmSminSminSmaxSmaxVminVminVmaxVmax0.90.61.491.495050

16、01 50075225图5算法改进前后效果对比Fig.5Comparison of the effect before and after the improvement of algorithm图6算法改进前后索力对比Fig.6Comparison of cable forces before and after the improvement of algorithm2.3优化结果分析控制拱桥的内力和线形使其在安全的范围内是施工监控过程的重要环节。为了能清楚地展示粒子群优化算法在实际工程中的效果，本节还另使用近似最小弯曲能量法、刚性支撑连续梁法及未知荷载系数法对索力进行求解，4种方法求得的

17、索力见图7。由图可知，利用改进后的粒子群优化算法求得的索力值较其他3种方法更加均匀适中。图7索力对比Fig.7Comparison of cable forces图8系梁弯矩Fig.8Bending moment of tie beam将4种方法求得的索力分别代入有限元模型中，得到各索力下系梁弯矩和竖向位移，如图8、9所示。从图8可以看出，用不同方法约束的系梁弯矩在数值和均匀度上都有较大的差异，其中未知荷载系数法约束的弯矩上下峰值差距较大，在边吊杆处弯矩较大，而刚性支撑连续梁法约束的弯矩在中吊杆处出现了骤增现象，这是因为在吊杆数量有限的情况下，这2种方法不能同时满足多目标约束，可能会导致结构其

18、余部位出现状态异常。近似最小弯曲能量法和改进粒子群优化法优化后的索力均没有出现索力突变现象。而本文的改进粒子群优化法优化后的系梁弯矩较近似最小弯曲能量法更加均匀，且索力值相差不大，适合系杆拱桥成桥索力布置。从图9可以看出，刚性支撑连续梁法和近似最小弯曲能量法优化后的竖向位移较大，这是因为这2种方法只是分别从结构的受力和弯曲能量进行了单一约束，而没有考虑结构的线形。未知荷载系数法是通过限制结构位移的方法来求解索力，因此位移较小，而改进粒子群优化法由于在考虑结构内力的同时对位移也加以约束，因此优化后的竖向位移较好，但本文的改进粒子群优化法较未知荷载系数法优化后的竖向位移更加均匀，峰值更小。4种方法对结构具体影响见表3。图9竖向位移Fig.9Vertical displacement表34种优化方法对结构的具体影响Table 3The specific influence of four optimization methods on the structure求解方法最大正弯矩/(kNm)最大负弯矩/(kNm)位移下限/mm位移上限/mm未知荷载系数法2 319.04-875.64-1.286.45改进粒子群优化法220.76-441.23-1.815.57刚性支撑连续梁法824.7