基于改进量子遗传算法的模型交互修正方法
Model interactive modification method based on improved quantum genetic algorithm
通讯作者:
收稿日期: 2021-06-12
基金资助: |
|
Received: 2021-06-12
Fund supported: | 国家自然科学基金资助项目(51878072);湖南省研究生科研创新资助项目(CX20190661);湖南省科技创新计划资助项目(2020RC4049) |
作者简介 About authors
向胜涛(1994—),男,博士生,从事钢-混组合结构大跨度桥梁施工控制的研究.orcid.org/0000-0002-3290-4050.E-mail:
针对传统有限元模型修正方法低效率、高成本且易陷入局部极值的缺点,提出基于改进量子遗传算法的静力多尺度有限元模型交互修正方法. 依据量子计算理论,对量子比特矢量态进行实数编码,以改进量子旋转门实现旋转角自适应更新,引入量子全局干扰交叉、变异、灾变等遗传操作,设计改进量子遗传算法. 以某钢-混组合梁桥为工程背景建立多尺度有限元模型,建立目标函数,对待修正区域进行分块处理. 利用最大互信息系数对待修正参数进行筛选,给出目标函数权重,通过Python语言实现了基于改进量子遗传算法的静力多尺度有限元模型交互修正. 结果表明,改进量子遗传算法相较于传统遗传算法、量子遗传算法具有更高的性能与精度,自动交互修正方法的效率较高,对材料弹性模量、厚度、车辆荷载等参数的修正与工程实际测试的情况基本吻合,目标函数修正结果相较于有限元计算的初始值,挠度误差降低至1.4%~14.3%,混凝土底板应力误差降低至2.6%~18.8%,钢梁应力误差降低至0%~11.1%.
关键词:
A static multi-scale finite element model interactive modification method based on improved quantum genetic algorithm was proposed aiming at the disadvantages of traditional finite element model modification methods, which are low efficiency, high cost and easy to fall into local extremums. The quantum bit vector states were encoded by real numbers according to the quantum computing theory, and the rotation angle was adaptatively updated by improving the quantum revolving gate. The improved quantum genetic algorithm was designed by introducing the quantum global interference crossover, mutation, catastrophe and other genetic operations. A multi-scale finite element model was established with a steel-concrete composite girder bridge as the engineering background. The objective function was established, the correction region was partitioned, and the maximum mutual information coefficient was used to screen the parameters and obtain the weight of the objective function. The interactive modification of static multi-scale finite element model based on improved quantum genetic algorithm was realized by Python language. Results show that the improved quantum genetic algorithm has higher performance and accuracy than the traditional genetic algorithm and quantum genetic algorithm, and the automatic interactive modification method is more efficient. The modification of material elastic modulus, thickness, vehicle load and other parameters accorded with the actual engineering test. The deflection error was reduced to 1.4%-14.3%, the stress error of concrete floor was reduced to 2.6%-18.8%, and the stress error of steel beam was reduced to 0%-11.1% compared with the initial finite element calculation.
Keywords:
本文引用格式
向胜涛, 王达.
XIANG Sheng-tao, WANG Da.
随着当前桥梁结构施工、运营、维护、健康监测等工作内容的精细化,对于结构计算的精度要求日益提高. 按照设计图纸建立的桥梁初始有限元(finite element, FE)模型相较于实际结构存在几何的简化与偏差、材料参数的不确定性、结构损伤与材料性能退化等问题,各类荷载下的结构响应往往难以与实际匹配,为此,对有限元模型进行修正[1].
模型修正本质上是优化问题,即根据实际结构对荷载的响应,修正FE模型的各参数,找到最优组合使计算结果与实际匹配. 按照目标数据的来源,可以将模型修正分为基于静力和动力的修正方法[2]. 相关学者在选取待优化参数、构造优化目标函数、改进优化算法等方面进行大量的研究. 秦世强等[3]基于菜园坝长江大桥荷载试验数据,结合Kriging代理模型及粒子群优化算法,对有限元模型进行修正,使得频率与位移的相对误差更小. Deng等[4]结合响应面和遗传算法对某简支梁桥及某既有桥梁的模型更新,利用响应面法进行参数更新的最佳试验设计,开展数值分析,从模拟结果中获得了结构响应和参数之间的显式关系. 方志等[5]基于参数灵敏度分析与ANSYS优化技术,利用静动力试验对某混凝土斜拉桥的静力位移和模态测试结果进行修正. Jaishi等[6]利用环境振动试验结果,考虑频率、模态振型、模态柔度,对某简支梁进行修正,开展敏感参数分析.
目前,模型修正方法均存在适用范围窄而不具普适性的缺陷[2]. 当使用代理模型时,修正结果的精度可能较差,需要通过有限元法(finite element method, FEM)进行验证,且待修正参数较多时可能难以得到最优解,部分情况下算法甚至不收敛. 当直接通过FEM进行迭代修正时,存在效率低下的问题[3]. 部分优化算法容易陷入局部极值,优化效率与精度不甚理想,且传统优化算法难以与FEM交互,对于修正参数与优化结果,往往须人工判断并手动输入FE模型,极大程度地降低了效率. 较多模型修正研究中对于结构不同部位的参数,采用统一的修正值进行修正[3-8],实际中不同部位的材料弹性模量、厚度可能均不同.
本文在相关学者研究的基础上,提出基于改进量子遗传算法(improved quantum genetic algorithm, IQGA)的静力多尺度FE模型交互修正方法. 对待修正参数较多的某钢-混组合梁桥进行模型修正,通过改进量子遗传算法,利用量子比特的矢量态进行实数遗传编码,利用改进量子旋转门实现旋转角自适应更新,引入量子全局干扰交叉、变异、灾变等操作,相较于其他传统算法具有更高的计算效率. 利用Python实现改进量子遗传算法与多尺度FE模型的交互,规避了对修正结果的人工判断与参数的手动输入,兼顾了计算效率与时间成本,提高了优化速度. 通过划分不同块建立单元集,在修正算法中分别考虑不同单元集的参数,采用最大互信息系数,验证了待修正参数与目标函数的相关性并进行筛选,采用加权平方误差函数作为适应度函数,提高了FE模型修正结果的准确性与计算效率.
1. 算法原理
1.1. 传统遗传算法
遗传算法(genetic algorithm,GA)是模拟生物进化和自然选择机理的优化算法,最早由Holland提出,通过随机选择、交叉、变异等遗传操作,使初始种群产生更适应环境的下一代个体,最后收敛于最优解[9]. 传统GA的主要流程如下.
1)确定编码方案并生成初始种群. 编码方案取决于问题性质与遗传算子设计,常用方案有二进制编码、实数编码、有序串编码等.
2)确定适应度函数. 通常以优化问题的目标函数作为个体适应度函数,但在某些情况下会将目标函数转化为标准适应度函数.
3)随机选择、交叉、变异. 随机选择是指从上一代群体中按一定概率选择个体到下一代群体中,被选中的概率与适应度函数值相关. 交叉是指2个个体的染色体进行交换组合,以产生可能更优秀的新个体. 变异是指对任意个体染色体中某一点或某些基因进行变异,以产生可能更优秀的新个体.
4)终止条件的确定. 当遗传算法运行至最大代数或适应度函数值达到可接受的范围内时,算法终止.
传统GA对于不能直接求解的复杂优化问题具有较好的适用性,且算法流程较为简单,但由于其需要计算大量个体的适应度函数值,单一遗传算法编码可能无法表示所有约束条件,对不可行解常采用阈值,算法效率通常相较于其他算法较低,容易出现早熟收敛与停滞的问题. 本文采用改进量子遗传算法,以进一步提高计算效率同时克服上述问题.
1.2. 改进量子遗传算法
1.2.1. 量子计算的基本概念
式中:
若干量子比特的有序集合可以构成量子寄存器,一般L位量子寄存器可以存储2L个不同的经典信息. 对量子比特进行幺正变换,称为量子门,Divincenzo[14]证明了对于任意复杂的幺正变换,均可以通过量子门的完备集合进行组合实现. 仅须满足单个量子比特在希尔伯特空间的任意转动和一对量子比特间的控制转动,即可实现量子态的幺正变换. 由于对一个叠加态的演化相当于对其中各个叠加成分的演化,即2L个态的演化,因而量子计算具有高度并行的特性.
1.2.2. 实数编码
基于
式中:
1.2.3. 量子旋转门选择策略
通过量子旋转门对个体基因(量子位)进行调整,对于量子门
式中:
令
考虑到目标函数在搜索点处的梯度,为了使梯度与转角步长呈负相关关系,参考何新贵等[16]考虑目标函数变化率所提出的新适应度函数,采用如下转角步长函数实现旋转角的自适应更新:
式中:
1.2.4. 量子全局干扰交叉
图 1
1.2.5. 量子变异和量子灾变
采用量子非门实现染色体的变异,通过指定变异概率选择染色体,随机选择一定数目的量子位基因,通过量子非门变换,使得指定量子位基因的2个概率幅互换. 通过量子变异,有利于算法控制早熟收敛,使得优化结果更有效.
在IQGA经过若干代进化后,群体可能已到达局部最优解趋向于早熟收敛,此时依靠交叉、变异操作可能无法使群体跳出局部最优陷阱. 该算法引入量子灾变机制,在保留当前最优解的前提下,重新初始化其他个体,使得算法迅速摆脱局部区域的搜索,从而极大程度地改善了IQGA的收敛性.
1.2.6. 基于MIC的待修正参数筛选方法
对于待修正参数,目前常根据经验进行选取,该选取方法具有一定的随意性,所选取的待修正参数在某些情况下可能对目标函数不敏感,导致待修正参数维度增加的负面效应. 该算法引入最大互信息系数(maximum mutual information coefficient, MIC)作为待修正参数筛选准则,在普适、公平的基础上有效衡量不同数据间的相关性[18]. 根据经验初步确定待修正参数后,通过FEM计算多组初始待修正参数与目标函数间的MIC值,舍去MIC较小(相关性较不显著)的待修正参数,可对待修正参数维度进行降维处理,减少优化算法的计算量. MIC法的基本步骤如下. 1)对随机变量X、Y进行排序并组成数据集,划分网格nx×ny并计算两变量互信息值,获取最大互信息值Im(见式(9)). 2)将最大互信息值归一化(见式(10)),选择不同尺度下的最大互信息值为MIC值(见式(11)).
式中:B(n)为网格分辨率,通常为数据量的0.6次幂.
1.2.7. 改进适应度函数
对于涉及的多目标优化问题,适应度函数通常可以选择均方误差函数(mean square error, MSE)来表征计算值与目标值的差异. 当计算至第
式中:
1.3. 算法性能对比
表 1 GA、QGA、IQGA算法的优化结果对比
Tab.1
函数 | 算法 | 最优搜索值 | 均值 | 最优理论值 |
| GA | 0.999 688 | 0.991 051 | 1 |
| QGA | 0.999 754 | 0.992 084 | 1 |
| IQGA | 0.999 991 | 0.999 112 | 1 |
| GA | 0.240 051 | 0.239 887 | 0.240 035 |
| QGA | 0.240 001 | 0.238 129 | 0.240 035 |
| IQGA | 0.240 035 | 0.239 911 | 0.240 035 |
图 2
图 2 GA、QGA、IQGA算法优化进程适应度函数值变化对比
Fig.2 Comparison of fitness function values of GA, QGA and IQGA optimization processes
式中:
1.4. 基于IQGA的参数化FE模型交互修正流程
利用ABAQUS自带宏管理功能,编辑Python脚本进行优化算法与FE模型计算的交互,可以实现结构参数化建模及自动修正,相较于传统的FEM计算与参数修正相对独立且须人工输入的方法,极大程度地提高了计算效率. 全流程算法如图3所示.
图 3
2. 实桥试验概况与FE模型
2.1. 工程背景和实桥试验
某高架桥全桥共四联,全长315.38 m,下部结构采用柱式花瓶墩,其中第2联为3×29.6 m钢-混组合梁,桥面全宽24.0 m. 采用9 cm沥青混凝土铺装,双向横坡坡度为1.5%,横向共布置8片工字钢主梁,间距为3.0 m,外侧桥面板挑臂长度为1.5 m,钢梁下翼缘横向平坡布置,标准横断面如图4所示.
图 4
竣工后对结构进行实桥试验,采用共8辆加载车进行加载,实测加载车的轴距Le、前轴重W1、后轴重W2、总重W如表2所示. 依据《公路桥梁荷载试验规程(JTG/T J21-01-2015)》可知,试验工况主要有中跨最大正弯矩(工况1),墩顶最大负弯矩(工况2),边跨最大正弯矩(工况3),主要测试内容包括跨中、边跨以及墩顶的边钢梁、中钢梁、混凝土板应变及挠度. 考虑到实桥的结构特点,其中挠度测点布置于桥面偏距为±10.5 m、±1.5 m处,钢梁应变测点布置于钢梁底板及顶板的下表面,混凝土应变测点布置于混凝土板底面.
表 2 加载车轴重
Tab.2
车辆编号 | Le/m | W1/kN | W2/kN | W/kN |
G1 | 4+1.4 | 80.5 | 322.1 | 402.6 |
G2 | 4+1.4 | 79.9 | 321.0 | 400.9 |
G3 | 4+1.4 | 78.6 | 319.0 | 397.6 |
G4 | 4+1.4 | 79.0 | 320.5 | 399.5 |
G5 | 4+1.4 | 78.9 | 321.0 | 399.9 |
G6 | 4+1.4 | 81.0 | 324.0 | 405.0 |
G7 | 4+1.4 | 40.5 | 160.1 | 200.6 |
G8 | 4+1.4 | 39.9 | 160.0 | 199.9 |
通过预加载以消除结构的非弹性变形,利用预加载检测数据采集设备是否处于正常工作状态,待测试结果回零表现良好后进行正式加载. 采用分级加载、卸载制度,对每一级加、卸载均进行测点数据采集,以确保加载过程安全、可靠.
2.2. 多尺度FE模型
在传统FE计算中,若对较大几何尺寸的结构建立单一尺度实体单元模型,则往往会造成FE模型单元数量庞大、计算效率低下的问题. 采用多尺度模型可以在保证精度可控的前提下大幅减少单元数量,对于替代实体单元的梁、杆、壳等单元,计算资源的消耗更低[8,19]. 依托设计图纸,建立ABAQUS多尺度FE模型,全桥包含实体、壳、杆3个尺度(见图5). 其中桥面板及铺装层采用C3D8R实体单元,钢梁采用S4R壳单元,钢筋采用T3D2桁架单元,对测试位置网格进行加密处理,以提高FEM的计算精度. 对于支座,通过建立弹簧,输入3个自由度方向的刚度实现对结构的相应约束,铺装层与桥面板、桥面板与工字钢采用绑定连接,普通钢筋采用内置区域的方式与桥面板单元耦合,初始材料参数容重
图 5
表 3 有限元模型初始材料参数
Tab.3
材料 | | E/MPa | | | 蠕变参数 | |||
A/10−11 | n | m | R2 | |||||
C50混凝土 | 24.2 | 34 500 | — | 0.2 | — | — | — | — |
Q345qd钢板 | 78.5 | 206 000 | — | 0.3 | — | — | — | — |
沥青混凝土 | 24.7 | — | 870 | 0.25 | 6.54 | 0.937 | −0.592 | 0.9326 |
3. 静力多尺度FE模型交互修正
3.1. 待修正参数的初步确定
根据施工过程中的实际测量结果,钢结构制造、安装精度相对较高,故不考虑将钢结构尺寸作为待修正参数. 混凝土桥面板由三部分组成:预制桥面板、湿接缝、现浇混凝土面板. 其中预制桥面板的尺寸精度相对较高,当预制板厚度出现误差时,在施工中通过混凝土标号相同的现浇面板加以找平,故将所有混凝土部分作为一个整体进行考虑. 混凝土护栏采用现浇方式进行施工,产生的抗力相较于主梁可以忽略,故不对混凝土弹性模量进行修正,仅考虑护栏高度带来的自重偏差. 铺装层较柔软,主要考虑实际厚度带来的自重偏差. 普通钢筋在实际定位中可能存在偏差,但由于钢筋定位偏差没有规律性,仅考虑弹性模量可能存在的偏差. 此外,加载车的轴重通过地磅确定,存在一定的测量误差. 综上所述,结合施工经验,选择钢材弹性模量X1、混凝土弹性模量X2、混凝土板总厚度X3、护栏高度X4、铺装层厚度X5、车辆荷载X6作为待修正参数.
其他传统修正方法中对相同材料或结构进行统一修正,未考虑材料特性及几何尺寸误差的离散性. 为了精确考虑不同部位的参数修正,对FE模型进行分块处理(见图6,铺装层为A1~A9,混凝土护栏为B1~B18,混凝土板为C1~C64,普通钢筋为D1~D64,钢纵梁为E1~E8,钢横梁为F1~F4). 对上述分块分别建立单元集,以便于后续优化算法中分别考虑不同块的材料、尺寸参数.
图 6
3.2. 目标函数的确定
基于实际的加载测试数据,以3种工况下跨中、边跨以及墩顶的边钢梁、中钢梁、混凝土板应力及挠度作为目标函数,编号如表4所示.
表 4 目标函数及权重
Tab.4
工况 | 位置 | 边梁目标函数 | 中梁目标函数 | 权重 |
1 | 中跨跨中 | 竖向挠度B1 | 竖向挠度Z1 | 0.075 |
混凝土底板应力B2 | 混凝土底板应力Z2 | 0.04 | ||
钢梁底板应力B3 | 钢梁底板应力Z3 | 0.05 | ||
钢梁顶板应力B4 | 钢梁顶板应力Z4 | 0.05 | ||
2 | 墩顶 | 混凝土底板应力B5 | 混凝土底板应力Z5 | 0.01 |
钢梁底板应力B6 | 钢梁底板应力Z6 | 0.03 | ||
钢梁顶板应力B7 | 钢梁顶板应力Z7 | 0.03 | ||
3 | 边跨跨中 | 竖向挠度B8 | 竖向挠度Z8 | 0.075 |
混凝土底板应力B9 | 混凝土底板应力Z9 | 0.04 | ||
钢梁底板应力B10 | 钢梁底板应力Z10 | 0.05 | ||
钢梁顶板应力B11 | 钢梁顶板应力Z11 | 0.05 |
3.3. 基于MIC的待修正参数筛选与目标函数权重确定
编制Python脚本,采用伪随机数方法实现上述各块初始待修正参数−X1、X2、X3、X4、X5、X6在修正范围内的随机修正,各个参数在随机修正下的组合即MIC法中的随机变量X(X = [X1, X2, X3, X4, X5, X6]),开展FEM计算,得到对应的目标函数组合Y(Y = [B1, B2, ···, B11, Z1, Z2, ···, Z11]). 采用式(9)~(11)计算网格nx×ny的MIC值,进行矩阵可视化,如图7所示. 图中,变量1~6依次对应X,变量7~28依次对应Y.
图 7
图 7 待修正参数与目标函数MIC矩阵热力图
Fig.7 MIC heat map of parameters to be modified and objective functions
从图7可以看出,X与Y的MIC值表现为强相关性,除了对角线自相关(MIC值为1.0)以外,最大为0.995,最小为0.141,故不舍去待修正参数,在后续IQGA-FEM交互修正中考虑所有块的参数修正.
为了给出目标函数的权重,由待修正参数与目标函数的MIC值矩阵热力图,分别计算目标函数组合Y各个元素对应随机变量X所有元素MIC值之和Xsum,进行归一化处理,结果见表5. 可知,Y各个元素对应随机变量X所有元素MIC值之和为[0.040 7, 0.051 4],考虑到实桥试验中挠度的测量精确度相较于应变更高,且中跨跨中加载重量最大,相应的系统误差、测量误差比例较小,相当于减小了误差所占的比重. 在混凝土底板应力测试中,由于现场条件的限制,应变片粘贴位置难以与理论达到同一精确度,因此通过试算对目标函数权重进行微调,并略微增大竖向挠度权重,适当减小混凝土底板应力的权重,对目标函数B1~B11、Z1~Z11采用的权重如表4所示.
表 5 目标函数各元素对应随机变量各元素MIC值之和
Tab.5
Y | Xsum | Y | Xsum | |
Z1 | 0.050 5 | B1 | 0.0461 | |
Z2 | 0.045 8 | B2 | 0.0451 | |
Z3 | 0.041 7 | B3 | 0.0451 | |
Z4 | 0.043 7 | B4 | 0.0514 | |
Z5 | 0.040 7 | B5 | 0.0457 | |
Z6 | 0.045 7 | B6 | 0.0429 | |
Z7 | 0.046 1 | B7 | 0.0445 | |
Z8 | 0.043 0 | B8 | 0.0432 | |
Z9 | 0.047 0 | B9 | 0.0477 | |
Z10 | 0.045 9 | B10 | 0.0467 | |
Z11 | 0.043 8 | B11 | 0.0477 |
3.4. IQGA-FEM交互修正结果
依照1.4节所示的流程实现模型的自动修正. 主要的初始参数如下:种群规模
表 6 参数修正结果
Tab.6
参数 | 块/编号 | 初始值 | 修正值 | 修正幅度/% | 参数 | 块/编号 | 初始值 | 修正值 | 修正幅度/% | |
X5/mm | A1 | 90 | 93 | 3.33 | X3/mm | C1 | 300 | 310 | 3.33 | |
A2 | 90 | 92 | 2.22 | C2 | 300 | 315 | 5.00 | |||
A3 | 90 | 95 | 5.56 | C3 | 300 | 316 | 5.33 | |||
X2/MPa | C1 | 34500 | 35900 | 4.06 | C4 | 300 | 316 | 5.33 | ||
C2 | 34500 | 36151 | 4.79 | C5 | 300 | 317 | 5.67 | |||
C3 | 34500 | 36151 | 4.79 | C6 | 300 | 315 | 5.00 | |||
C4 | 34500 | 36991 | 7.22 | C7 | 300 | 312 | 4.00 | |||
C5 | 34500 | 36991 | 7.22 | C8 | 300 | 316 | 5.33 | |||
C6 | 34500 | 37100 | 7.54 | C9 | 300 | 319 | 6.33 | |||
C7 | 34500 | 37890 | 9.83 | C10 | 300 | 315 | 5.00 | |||
C8 | 34500 | 37100 | 7.54 | C11 | 300 | 319 | 6.33 | |||
C9 | 34500 | 37201 | 7.83 | C12 | 300 | 313 | 4.33 | |||
C10 | 34500 | 37333 | 8.21 | C13 | 300 | 312 | 4.00 | |||
C11 | 34500 | 37332 | 8.21 | C14 | 300 | 319 | 6.33 | |||
C12 | 34500 | 37332 | 8.21 | C15 | 300 | 312 | 4.00 | |||
C13 | 34500 | 37333 | 8.21 | C16 | 300 | 319 | 6.33 | |||
C14 | 34500 | 37199 | 7.82 | C17 | 300 | 315 | 5.00 | |||
C15 | 34500 | 37101 | 7.54 | C18 | 300 | 304 | 1.33 | |||
C16 | 34500 | 36001 | 4.35 | C19 | 300 | 313 | 4.33 | |||
C17 | 34500 | 35100 | 1.74 | C20 | 300 | 303 | 1.00 | |||
C18 | 34500 | 37155 | 7.70 | C21 | 300 | 306 | 2.00 | |||
C19 | 34500 | 37161 | 7.71 | X4/mm | B1 | 1330 | 1320 | −0.75 | ||
C20 | 34500 | 36600 | 6.09 | B2 | 1330 | 1329 | −0.08 | |||
C21 | 34500 | 37091 | 7.51 | B3 | 1330 | 1328 | −0.15 | |||
X1/MPa | D1 | 206000 | 209970 | 1.93 | B10 | 1330 | 1330 | 0.00 | ||
D2 | 206000 | 209903 | 1.89 | B11 | 1330 | 1327 | −0.23 | |||
D3 | 206000 | 209918 | 1.90 | B12 | 1330 | 1329 | −0.08 | |||
D4 | 206000 | 209999 | 1.94 | X1/MPa | E1 | 206000 | 209800 | 1.84 | ||
D5 | 206000 | 212111 | 2.97 | E2 | 206000 | 209953 | 1.92 | |||
D6 | 206000 | 211058 | 2.46 | E3 | 206000 | 210471 | 2.17 | |||
D7 | 206000 | 210001 | 1.94 | E4 | 206000 | 210513 | 2.19 | |||
D8 | 206000 | 210000 | 1.94 | E5 | 206000 | 209699 | 1.80 | |||
D9 | 206000 | 210099 | 1.99 | E6 | 206000 | 210987 | 2.42 | |||
D10 | 206000 | 212000 | 2.91 | E7 | 206000 | 212100 | 2.96 | |||
D11 | 206000 | 211107 | 2.48 | E8 | 206000 | 210117 | 2.00 | |||
D12 | 206000 | 210011 | 1.95 | F1 | 206000 | 211098 | 2.47 | |||
D13 | 206000 | 212109 | 2.97 | F2 | 206000 | 209987 | 1.94 | |||
D14 | 206000 | 210000 | 1.94 | X6/kN | G1 | 402.6 | 384 | −4.62 | ||
D15 | 206000 | 210100 | 1.99 | G2 | 400.9 | 382.1 | −4.69 | |||
D16 | 206000 | 212010 | 2.92 | G3 | 397.6 | 380.1 | −4.40 | |||
D17 | 206000 | 212087 | 2.95 | G4 | 399.5 | 381.1 | −4.61 | |||
D18 | 206000 | 212033 | 2.93 | G5 | 399.9 | 391.5 | −2.10 | |||
D19 | 206000 | 212011 | 2.92 | G6 | 405 | 386.7 | −4.52 | |||
D20 | 206000 | 209410 | 1.66 | G7 | 200.6 | 191.4 | −4.59 | |||
D21 | 206000 | 211099 | 2.48 | G8 | 199.9 | 190.5 | −4.70 |
从表6可以看出,铺装层修正后厚度增加,说明了现场实际摊铺存在正误差. 修正后各材料弹性模量较初始值普遍增大,其中混凝土弹性模量最大增加了9.83%,钢材弹性模量最大增加了2.97%,这表明现场实际采用的混凝土弹性模量较标准值偏大,这与现场实验室对混凝土及钢材检测报告中弹性模量偏大的趋势保持一致. 模型未考虑钢梁的部分细节如拼接板、螺栓、局部肋板等构造,故算法对钢梁的弹性模量进行了小幅修正,这与实际情况较相符. 混凝土板修正后厚度普遍增加,说明在混凝土板预制及现场浇筑过程中存在一定的正误差,这与现场的预制板验收情况基本保持一致. 对护栏高度的修正存在离散性,此处未表现出明显的规律. 对车辆荷载的修正在合理范围内,普遍为负误差,说明对车辆称重的地磅存在系统性偏差.
表 7 目标函数计算结果与实测值的对比
Tab.7
函数值 | B1 | B2 | B3 | B4 | B5 | B6 | B7 | B8 | B9 | B10 | B11 |
实测值 | −7.0 | −3.2 | 44.3 | −4.1 | −1.6 | −28.7 | 5.3 | −8.0 | −3.8 | 43.3 | −5.4 |
初始值 | −8.8 | −4.0 | 48.8 | −4.5 | −2.1 | −29.1 | 6.0 | −10.0 | −4.2 | 46.7 | −5.5 |
修正值 | −8.0 | −3.5 | 46.8 | −4.3 | −1.9 | −28.8 | 5.6 | −9.1 | −3.9 | 44.7 | −5.3 |
函数值 | Z1 | Z2 | Z3 | Z4 | Z5 | Z6 | Z7 | Z8 | Z9 | Z10 | Z11 |
实测值 | −12.0 | −2.1 | 30.6 | −2.7 | −1.1 | −19.4 | 3.0 | −14.0 | −2.5 | 30.6 | −3.7 |
初始值 | −12.9 | −2.5 | 33.6 | −3.2 | −1.5 | −21.1 | 3.4 | −14.8 | −2.7 | 32.6 | −4.0 |
修正值 | −12.3 | −2.2 | 32.0 | −3.0 | −1.3 | −19.6 | 3.1 | −14.2 | −2.7 | 30.8 | −3.7 |
图 8
图 8 修正前、后目标函数误差的对比
Fig.8 Comparison of objective function error before and after modification
4. 结 论
(1)传统GA存在早熟收敛、停滞及计算效率的问题,结合量子计算与GA,采用QGA可以实现更优越的性能,但传统QGA受限于编码方式、量子旋转门策略等,计算效率有待提高.
(2)通过实数编码、改进量子旋转门策略、旋转角自适应更新、引入全局干扰交叉、量子变异、量子灾变等,实现了对传统QGA算法的改进. 提出IQGA算法,通过Schaffer函数与Bohachevsky函数对传统GA、QGA及本文IQGA算法的性能进行对比. 结果表明,IQGA在进化代数不超过50的情况下首先达到最优解,性能显著高于其余2种算法.
(3)通过Python语言,实现参数化建模. IQGA参数全自动交互修正,具有极高的效率. 通过引入加权平方误差函数作为适应度函数,可以分别考虑不同目标函数的权重. 通过对多尺度FE模型进行分块处理,规避了其他模型修正方法中对相同材料参数与结构尺寸统一修正的缺陷. 利用MIC法对待修正参数进行筛选,明确了待修正参数与目标函数之间的强相关性. 该模型修正方法考虑了实际结构材料性能、几何尺寸的离散性,结果更精确、可靠. 在最终的优化结果中,挠度误差降低至1.4%~14.3%,混凝土底板应力误差降低至2.6%~18.8%,钢梁应力误差降低至0%~11.1%,说明了本文方法的适用性及精度较高,可以为健康监测、施工监控、桥梁检测等工作提供参考.
参考文献
Operational modal testing and FE model tuning of a cable-stayed bridge
[J].DOI:10.1016/j.engstruct.2011.02.046 [本文引用: 1]
基于桥梁健康监测的有限元模型修正研究现状与发展趋势
[J].DOI:10.3969/j.issn.1671-8879.2014.04.009 [本文引用: 2]
Present research status and development trend of finite element model updating based on bridge health monitoring
[J].DOI:10.3969/j.issn.1671-8879.2014.04.009 [本文引用: 2]
基于试验数据的大跨度拱桥有限元模型修正
[J].
Finite element model updating of large-span arch bridge based on experimental data
[J].
Bridge model updating using response surface method and genetic algorithm
[J].DOI:10.1061/(ASCE)BE.1943-5592.0000092 [本文引用: 1]
基于多状态下静动态测试数据的斜拉桥模型修正
[J].DOI:10.3969/j.issn.1001-7372.2011.01.006 [本文引用: 1]
Cable-stayed bridge model updating based on static and dynamic test data of multi-state
[J].DOI:10.3969/j.issn.1001-7372.2011.01.006 [本文引用: 1]
Structural finite element model updating using ambient vibration test results
[J].DOI:10.1061/(ASCE)0733-9445(2005)131:4(617) [本文引用: 1]
Finite element model updating in structural dynamics by using the response surface method
[J].DOI:10.1016/j.engstruct.2010.04.019
基于响应面法的钢管混凝土组合桁梁桥多尺度有限元模型修正
[J].
Multi-scale finite element model updating of CFST composite truss bridge based on response surface method
[J].
An improved genetic algorithm for determining modified water-retention model for biochar-amended soil
[J].
Quantum-inspired evolutionary algorithms with a new termination criterion, H ε gate, and two-phase scheme
[J].DOI:10.1109/TEVC.2004.823467 [本文引用: 1]
Solving engineering optimization problems using an improved real-coded genetic algorithm (IRGA) with directional mutation and crossover
[J].DOI:10.1007/s00500-020-05545-9
一种解决组合优化问题的改进型量子遗传算法
[J].DOI:10.3321/j.issn:0372-2112.2007.10.034 [本文引用: 2]
A novel improved quantum genetic algorithm for combinatorial optimization problems
[J].DOI:10.3321/j.issn:0372-2112.2007.10.034 [本文引用: 2]
一种解决组合优化问题的量子遗传算法QGA
[J].DOI:10.3321/j.issn:0372-2112.2004.11.022 [本文引用: 1]
A quantum genetic algorithm to solve combinatorial optimization problem
[J].DOI:10.3321/j.issn:0372-2112.2004.11.022 [本文引用: 1]
Two-bit gates are universal for quantum computation
[J].DOI:10.1103/PhysRevA.51.1015 [本文引用: 1]
基于实数编码和目标函数梯度的量子遗传算法
[J].DOI:10.3321/j.issn:0367-6234.2006.08.002 [本文引用: 1]
Quantum genetic algorithm based on real encoding and gradient information of object function
[J].DOI:10.3321/j.issn:0367-6234.2006.08.002 [本文引用: 1]
利用目标函数梯度的遗传算法
[J].
Genetic algorithms using gradients of object functions
[J].
一种改进的量子遗传算法
[J].DOI:10.3969/j.issn.1671-1815.2012.12.014 [本文引用: 1]
An improved quantum genetic algorithm
[J].DOI:10.3969/j.issn.1671-1815.2012.12.014 [本文引用: 1]
Detecting novel associations in large data sets
[J].DOI:10.1126/science.1205438 [本文引用: 1]
Structural health rating (SHR)-oriented 3D multi-scale finite element modeling and analysis of Stonecutters Bridge
[J].DOI:10.12989/sss.2015.15.1.099 [本文引用: 1]
Updating multi-scale model of a long-span cable-stayed bridge
[J].
/
〈 |
|
〉 |
