0引言在桥梁健康监测中,传感器布设是首先要解决关键技术之一.
		传感器布设即通过尽可能少的传感器来获取最全面可靠的桥梁健康
		状况信息.在这方面已经有许多学者做了大量的研究工作,常用的方
		法有有效独立法、模态动能法、奇异值分解法、Guyan 模型缩减法
		、遗传算法等等[13].然而,现有的优化布设都是在传感器数量
		确定的基础上进行位置优化,对传感器的数量如何优化研究得非常
		少,而传感器的数量决定了桥梁健康监测系统的造价及最终监测结
		果的精度和可靠度,因此,传感器布设非常有必要从定量和定位两方
		面同时进行研究.通过研究发现,传感器的数量与所选的模态阶数是
		对应的,模态阶数选取的不同传感器的数量就会不同.因此,确定传
		感器的数量之前要进行模态阶数的选取.由此,本文基于振型贡献率
		和遗传算法,提出三步法对传感器的优化布设进行研究,具体步骤为
		:a. 确定振型模态阶数;b. 确定传感器的数量;c. 确定传感器的位
		置.本文的研究方法使传感器的布设更加全面系统,在满足桥梁健康
		监测需要的同时,可以最大程度地节约造价,具有广阔的工程应用前
		景.1模态阶数确定1.1振型参与系数概念模态阶数的选取是确定传
		感器数量进而确定传感器位置的前提.如果选取的模态过多,不仅会
		耗费大量的计算时间和空间,而且难以对计算结果的正确性作出准
		确的界定.那么如何选取模态阶数呢?为了解决这个问题,一些学者
		提出了振型参与系数的概念,即为该振型的影响占所有振型的影响
		的比值来计算[4].τm=∑Miφim∑Miφ2im(1)
式(1)中:τm为振型参与系数;m为振型阶数;Mi为i节点位置质量;
		φim为i节点位置的第m阶振型向量.1.2振型贡献率概念振型参与系
		数虽在一定程度上反映了各阶振型反应在总反应中所占的比例,但
		该参数有正负号不便比较.为解决此问题, 提出了振型贡献率的概
		念,即振型的有效质量占总质量的百分比,其本质是基于有效质量法
		来确定振型的选取阶数.这个概念已成功应用于抗震设计中[46
		],但在桥梁传感器的布置中还无人应用.本文将引入此概念来确定
		桥梁结构传感器布置所需的模态阶数.有效质量由式(2)计算:Mm=
		[∑φimMi]2∑2φimMi(2)
式(2)中:Mm为各振型有效质量.定义有效质量与总质量的比率为振
		型贡献率或模态贡献率:γm=MmM(3)γm为第j 阶模态模态贡献率.
		实际计算中,若取式中k 阶模态参与计算,k<n ,则此k阶模态的总贡
		献率为:∑ki=1γi=R(4)一般抗震设计规范中会规定各振型的R值
		要取到90%以上.这是为了保证计算地震作用时包含足够的主要振型
		.在桥梁结构传感器的优化布置中建议取R≥0.7,从而确定所需的模
		态阶数.2遗传算法设计2.1遗传算法简介遗传算法是模拟生物在自
		然环境中的遗传和进化过程而形成的一种自适应概率性全局收敛的
		优化算法.本文遗传算法设计中,编码方法采用二进制编码,选择算
		子采用比例选择,交叉算子采用单点交叉,变异采用倒位遗传算子.
		为了增加种群的多样性,避免早熟提高收敛速度,本文还引入了最优
		保留策略使每一代适应度值最优个体不参与交叉和变异直接遗传到
		下一代,引入自适应交叉及变异概率对简单的遗传算法进行了改进
		[7].遗传算法具体操作步骤如下:①进化代数初始化;②采用二
		进制编码,形成初始种群;③计算适应度值;④采用比例选择算子进
		行选择操作,之后附带最优保留策略;⑤采用单点交叉算子及自适应
		交叉概率进行交叉操作;⑥采用倒位算子及自适应变异概率进行变
		异操作;⑦终止条件判断,不满足返回步骤③.第12期金中凡,等:基
		于振型贡献率和遗传算法的传感器布设研究
武汉工程大学学报第32卷
2.2适应度函数设计遗传算法在操作进化过程中以适应度函数来评
		价染色体或种群的优劣,本文采用模态确认准则MAC矩阵作为适应度
		函数MAC矩阵,通过表达式(5)定义[8].MACij=(φTiφj)2
		(φTiφi)(φTjφj)(5)
式(5)中:φi和φj分别为第i 阶和第j 阶振型模态向量.MAC矩
		阵的非对角线元素反映了各振型向量间的正交性.这就要求应使得
		MAC矩阵的非对角线元素越小越好,而遗传算法中,个体的适应度是
		个最大化问题.因此,将MAC矩阵进行转化得适应度函数(优化目标函
		数)如式(6).F(t)=1-f(t)if f(t)<1
0if f(t)>1(6)
式(6)中,f(t)=max{MACij},(i≠j).2.3改进交叉设计在确定传
		感器位置时,遗传算法中的传感器的数量是固定的.由于本文采用的
		编码方法是二进制编码,普通的交叉操作会改变二进制编码中1的数
		量即传感器的数量,因此,本文采用单点交叉来解决这一问题.具体
		做法如下:设父代的两个个体分别为A和B,随机产生一个交叉点i,
		新产生的二进制编码串需满足式(7)条件才进行交叉操作.∑
		nk=iAk=∑nk=iBk(7)
式(7)中:Ak及Bk为父代个A和B中的第k个二进制编码;n为个体的
		总长度.3传感器数量及位置确定 由于遗传算法为概率搜索算法,因
		此每次遗传操作所得的结果不尽相同,为此本文在确定传感器数量
		时,在遗传算法中加入循环策略,即多次循环遗传算法,利用MATLAB
		软件编制了遗传算法优化传感器数量程序,程序框图如图1所示.图1
		优化传感器数量程序框图
Fig.1Program diagram on optimizing sensor number通过此程
		序可提取随着传感器数量增加n次循环遗传算法操作的适应度函数
		平均值,并以传感器数量N为横坐标,以适应度函数平均值f为纵坐标
		,绘制N-f曲线.通过N-f曲线的突变点可确定出最经济的传感器数量
		.传感器的数量确定以后,即可采用本文改进的遗传算法进行传感器
		位置的优化.4算例分析本文以珲乌高速公路辅道洮儿河大桥一联4
		×25 m连续梁桥为研究对象,利用大型计算软件ANSYS10.0建立了有
		限元模型,如图2所示.图2桥梁的有限元模型
Fig.2Finite element model of bridge因为遗传算法计算适应度
		函数时需要模态振型数据,所以在此分析该座桥梁的动力特性.本桥
		监测方案重点是桥梁的面内竖向振动,安放的传感器为竖向加速度
		传感器.洮儿河大桥计算模型共划分了104个单元,105个节点,在所
		有节点自由度上布设传感器是不可能的,因此本文沿桥纵向每隔
		3.125 m (L/8) 做为预布传感器测点,共提取了33个节点的振型数
		据.基于振型贡献率和遗传算法,本文采用三步法对该桥传感的布设
		进行了研究,具体步骤如下:第一步:确定振型模态阶数.根据振型
		参与因子的概念,采用ANSYS对洮儿河大桥进行模态分析及谱分析,
		在整体坐标系竖向施加单位位移谱激励,获得前十二阶振型参与因
		子,由式(2)~(4)计算该桥振型贡献率及其累积值R,计算结果列于表
		1中.表1振型贡献率及累积值
Table 1Modal mass ratio and cumulative percentage
模态阶数振型参
与因子振型贡献率累积值10.000.000.002177.3211.330.113 
		330.000.000.113 34401.8058.170.695 050.000.000.695 
		0654.441.070.705 670.000.000.705 68-5.970.010.705 
		890.000.000.705 81063.101.430.720 1110.000.000.720 
		112177.3711.340.833 5对于梁式桥而言,低阶模态较易获得,本文
		建议R取至0.7.从表1结果可以看出,前四阶振型模态总贡献率已达
		0.695,已经非常接近模态贡献率累积值R≥0.7的要求,因此可选取
		前四阶振型数据作为遗传算法中计算适应度函数时需要模态振型数
		据.1~4阶振型动力特性如图3所示.图3桥梁动力特性计算结果
Fig.3Calculation result of dynamic character of bridge第
		二步:确定传感器的数量.选取前四阶振型按图1所示程序进行传感
		器数量的计算,传感器数量初值N0=1;遗传算法循环次数n=20;遗传
		算法中群体规模取50;进化代数T=100.将每次循环遗传算法所得适
		应度函数平均值f用不同符号表示,以传感器数量N为横坐标,以f为
		纵坐标,绘制N-f点据分布图,如图4所示.图4传感器数量与适应度函
		数关系曲线
Fig.4Relationship between the number of sensors and 
		fitness由图4可见,拟合曲线的割线斜率随着传感器数量的增加逐
		渐减小,当N≥4时趋于水平,因此,最经济传感器数量可确定为4.第
		三步:确定传感器的位置.在确定传感器的数量为4后,利用本文设
		计的遗传算法对归流河大桥传感器的位置进行优化,遗传种群规模
		取100,最大进化代数取100,优化所得历代最佳适应度进化曲线见图
		5,传感器具体位置如图6中圆点所示.图5历代最佳适应度进化曲线
Fig.5Evolution curve of maximum fitness in every 
		iteration图6传感器位置布置图
Fig.6Placement map of sensor location当传感器数量N=4时,
		适应度函数值为1,即MAC矩阵非对角元最大值为0,说明改进遗传算
		法大大改善了测取目标模态间的振型可分性,优化结果合理.另外,
		从图6也可以看出,4个传感器分别布置在桥梁每跨的跨中位置,这与
		静力计算结果也吻合.5结语大多数研究者把传感器布设的重点集中
		于传感器位置的确定,而对传感器数量的研究却非常少.本文基于振
		型贡献率和遗传算法从定量和定位两个方面对桥梁结构传感器的布
		设进行了研究,得到如下结论:a. 将振型贡献率用于模态阶数的选
		取从而确定传感器数量是切实可行的,对于梁式桥而言振型贡献率
		累积值取至0.7可满足桥梁结构传感器布设要求.b. 利用遗传算法,
		通过传感器数量与适应度函数关系曲线可以确定出最经济的传感器
		数量,大降低桥梁健康监测造价.c. 从桥梁动力计算出发,以MAC矩
		阵为适应度函数的遗传算法可布设出合理经济的传感器位置,动力
		计算结果与静力计算相吻合.通过应用于工程实例的结果表明,本文
		所提出来的三步法为实现桥梁结构传感器的布设提供了一种更加全
		面和有效的途径.