文章编号:1004-3918(2008)05-0566-05
河南科学
HENANSCIENCEVol.26No.5May2008
遗传算法在结构振动控制中的应用
赵
摘
林1,张猛2
450046;2.郑州大学土木工程学院,郑州
450002)
(1.上海绿地集团郑州房地产事业部,郑州
在建筑结构系统辨识、振动控制算法、结构振动控制作动器优化布置和作动器参数优化要:对遗传算法(GA)
遗传算法是一种高效的搜索计算和优化计算方等研究方向的应用进行了较全面研究和计算机仿真分析.结果表明,
法,这一算法非常适合处理结构振动控制领域中的诸多复杂优化问题,能够精确地识别复杂结构系统的动态响应,实现结构振动控制系统的优化设置.
关键词:遗传算法;系统辨识;优化设计;作动器;控制算法中图分类号:TU311.3
文献标识码:A
遗传算法(GA)是计算智能技术的重要分支.它根据适者生存、优胜劣汰等自然进化规则来进行搜索计
算和问题求解.对许多用传统数学方法难以解决或明显失效的复杂问题,特别是优化问题,GA提供了一条行之有效的新途径.随着GA算法的不断成熟和计算机技术的迅速发展,近年来这种方法得到了土木工程界的重视和应用.在结构振动控制中,作动器的定位对控制效果是至关重要的,不恰当的布置会导致控制溢出和控制效果不良.本文对GA在建筑结构系统辨识、振动控制算法、结构振动控制作动器优化布置和作动器参数优化等研究方向的应用进行全面分析.
1应用GA的建筑结构系统参数辨识
建筑结构系统辨识技术在结构振动控制系统设计和结构破损诊断技术中都有重要意义.基于GA的结构系统参数辨识流程如图1所示.
图1
Fig.1
应用遗传算法的建筑结构系参数统辨识流程
Flowchartofbuildingstructuresystemidentificationwithgeneticalgorithms
设结构在地震作用下的运动方程为
!(t)\"(t)!(xMX+CX+KX(t)=-M{I}),gt
(1)
T
式中:M,C,K分别为结构质量、阻尼、刚度矩阵;X(t)=[x(),x(),…,x()]是结构各层相对于地面的位1t2tnt
!(x移;)为地震加速度.现应用遗传算法对结构各层刚度进行辨识.gt
T设!=[k1,k2,…,kn],将k1,k2,…,kn表示成二进制代码,似然准则取为
T
M=!(v(i)R-1v(i)+ln│R│).
i=1n
(2)
收稿日期:2008-01-16
基金项目:九五国家自然科学基金重大项目(59895410);国家自然科学基金项目(50078037)资助作者简介:赵
林(1978-),男,河南信阳人,工程师,工学博士,主要从事结构振动控制方面的研究.
2008年5月
L
赵林等:遗传算法在结构振动控制中的应用-567-
式中:R为协方差矩阵;v(i)=!(xj-xj)为残差信息序列;L是测点个数;n是外界激励采样点数;xj是设置在
\"j=1
j点的传感器测得的结构响应(位移、速度或加速度).为验证上述方法,本文选用某3层结构进行仿真分析.结构参数为:各层质量为1.0×105kg,第1层刚度k1为1.1×108N/m;第2层刚度k2为0.7×108N/m;第3层刚度k3为1.2×108N/m,各层阻尼为1.0×105N・s/m.要辨识的参数为结构各层刚度.GA的初始种群大小取为100个,交叉概率为0.9,变异概率为0.1,最大迭代次数为50次.仿真测量值为结构在EL-CENTRO(NS方向)
波作用下的结构位移响应.
图2地震输入与结构顶层位移响应的传递函数误差比较
经过辨识得到的结构第1层刚度为k1=1.1063×108N/m;第2层刚度为k2=0.6969×108N/m;第3层刚度为k3=1.1929×108N/m.可见辨识精度很高.地震输
Fig.2Comparisonoftransferfunctionerrorofseismicinputand
thestructuretoplayerdisplacement
x(!)
(!)3的辨识误差与理论真值如图2所示.入与结构顶层位移响应的传递函数H
!Xg(!)
2应用GA的结构振动控制算法
设线性时不变结构振动控制状态方程的差分形式为
!(!(xx{z(k)}=[A0]{(zk-1)}+[B0]{u(k-1)}+[B1]{u(k)}+[F0])+[F1]).gk-1gk
(3)
-1-2
=eA!t为2n×2n阶离散系统矩阵;[B0]=([A0][A]+1[A0](I-[A]))[B]为2n×p阶前一时刻的其中:[A0]
!t-1-2-1
控制矩阵;[B1]=(-[A0]+1[A0]([A]-I))[B]为2n×p阶当前时刻的控制矩阵;[F0]=([A0][A]+
!t
1[A]-2-1-2
[A]))[F]为2n×1阶前一时刻的地震矩阵;[F1]=(-[A0]+1[A0]([A]-I))[F]为2n×1阶0(I-
!t!t当前时刻的地震矩阵;n为结构自由度数;p为控制力向量维数.目标函数的确定可根据实际的要求进行,在此采用式(4)所示的目标函数
TT
J(k)={z(k)}Q{z(k)}+{u(k)}R{u(k)}.(4)
其中:Q为状态量的权矩阵;R为控制力的权矩阵.这样就可以利用遗传算法对结构进行在线控制.设[umini,
umaxi]为控制力向量中的元素ui的变化范围.用二进制字符串表示每个控制力分量,则控制力字符串表示值和实际控制力值之间有如下换算关系
(L)ui=umini+bin[umaxi-umini],(5)L2-1其中bin(L)为由L位字符串构成的二进制数.
在线控制的具体方法为:从第一采样时刻k=1开始,采集地震数据计算控制力.1)初始群体的产生.随机产生一组字符串,每一字符串表示一种控制力方案.每一字符串都由p个控制力编码级联而成;
2)将字符串译码成各个控制力向量,然后计算结构响应,计算目标函数值,并取适应度函数为1/((Jk)+1),得到每一组控制力向量的适应度值;3)应用选择、交叉、变异算子对控制力进行优化求解;
4)判断适应度值是否满足精度要求,或者迭代次数是否等于最大后代数,若不满足,则重复(2),(3)步,否则进行下一步;
5)选择群体中适应度最好的一组控制力施加于结构.-568-
河南科学
第26卷第5期
进行k=k+1时刻的计算,重复步骤(1)~(5).为了验证上述算法,以三层结构为例进行仿真分析.结构参数如表1所示.在结构顶层布置一套控制装置.地震输入为EL-CENTRO(NS方向)波.在遗传算法中,用12位长的二进制字符串对控制力进行标码,种群规模为30个,最大迭代次数取值为50,交叉概率为0.9,变异概率为0.35.当限定控制力小于1000时,结构控制效果如图3所示.
表1
Tab.1层号刚度/(N・m-1)质量/kg
某三层结构参数
1
23022200981
31333600981
Parametersofabuildingstructure
图3
Fig.3
应用遗传算法的结构位移控制效果
withgeneticalgorithms
2741700981
Comparisonofstructuredisplacementcontroleffect
3应用GA的结构振动控制作动器优化布置
在保证达到振动控制设计目标的前提下,如何优化布置作动器是一个重要问题.其目的是:①达到振动控制设计目标的要求;②使所布置的作动器充分发挥其控制能力.具有n个自由度的线形结构半主动控
#()的作用下,考虑半主动控制输入的运动方程为制系统,在一维水平地面加速度xgt
#+CX$+KX=-MIx%(MX)+BU.gt
(6)
式中:B为作动器分布矩阵;U为系统控制力列阵.那么确定作动器布置方案问题就是在某种准则下,对B
矩阵进行优化.本文提出应用改进遗传算法来解决上述问题,并应用此方法进行了高层结构振动半主动控
T
制作动器优化布置研究.对于(6)式描述的结构控制系统,令b=[b1,b2,…,ki,…,kn]为控制器位置向量,这样b就是由1和0构成的一个二进制编码串.与GA一般形式相比,此时不需要再计算十进制真值,根据这个二进制编码串就可以计算适应度函数了.设计两个适应度函数如下:
!0,goal<response,##
fit1=\"1#,goal>response,#$goal-response式中:goal为振动控制设计目标;response为振动控制设计目标要求控制的结构响应绝对最大值,如位移、速度或加速度.这一适应度函数保证种群向着振动控制设计目标进化.
X(t)QX(t)dt-%X(t)QX(t)dt%fit2=.
T
Tc
0
nc
nc
0
c
T
T
L式中:Xnc(t)为无控状态反应向量(也可以是无控时的位移或速度反应向量);X()为控制后的状态反应向量ct(也可以是控制后的位移或速度反应向量);Q为权矩阵,根据控制设计目标来确定;L为控制器布置总层数.这一目标函数保证了所布置的控制器都能充分发挥作用.
为加快收敛速度,同时保证种群的多样性和全局收敛,本文采用自适应交叉和变异概率计算方法[8].改进的遗传算法流程如图4所示.
图4改进遗传算法流程图
Fig.4Flowchartofupdatedgeneticalgorithms2008年5月赵林等:遗传算法在结构振动控制中的应用-569-
为验证上述方法,本文以一个11层框架结构的AVS/D半主动控制问题为例进行了计算机仿真分析.结构参数如表2所示.进行抗震设计后,在输入EL-CENTRO(NS方向)波的情况下,结构顶层最大位移反应为14.36cm,不能达到最大位移goal=10cm的设计要求.因此在进行抗震设计之后,应进行振动控制设
附加刚度为8.11×107N/m,附加阻尼为2.64×107N・s/m.计.所布置的AVS/D装置的参数是:
表2
Tab.2
层号刚度/(108N・m-1)重力荷载/kN
16.439633
28.118631
38.118631
某11层框架结构参数
48.11
58.118631
68.118413
76.518413
86.518413
96.518413
106.518413
116.516218
Parametersofaframeconstruction
8631
现应用本文提出的的改进遗传算法进行作动器优
化布置.要求在达到振动控制设计目标的前提下,得到作动器的最优布置方案.控制律采用推广的KOBORI开关控制律.计算得到的最优布置方案是在第1层、第2层、第8层和第9层布置作动器,此时顶层最大位移值为9.515cm.顶层位移控制效果如图5所示.
4应用GA的结构振动控制作动器参数寻优
对于各种结构振动控制系统都存在在达到目标的
前提下作动器参数最优的问题.对于TMD减振系统而言,其减振效果和TMD与主结构的频率比、TMD与主
图5顶层位移控制效果
结构的质量比、TMD系统的阻尼、位置和外激励的频谱
Fig.5Comparisonofstructuretoplayerdisplacement
组成有密切关系.因此寻找TMD最优参数是将之更有controleffectwithgeneticalgorithms效地应用于工程实践的关键之一.
本文应用遗传算法直接对影响TMD控制的参数进行编码,通过群体搜索策略和个体信息交换等进化操作,找出了影响TMD控制性能的最优或次优参数,解决了传统参数优化方法存在的易陷入局部极小、迭代时间长、难以找到最优或次优参数等弊病,为TMD在实际工程的应用提供了理论依据.
设在地震作用下主结构和TMD控制系统第一阶模态运动方程为
!\"\"\"\"\"\"\"#
1011
$’%%(%(%%(%%(&)
%q!v
*++++,
+
!\"\"\"\"\"\"\"#
2!1\"1-2#!d\"d
0
2!d\"d
$’%%(%(%%(%%(&)
%q!v
*++++,
+
!\"\"\"\"\"\"\"#
\"2-$\"21d0
2\"2
d
$’%%(%(%%(%%(&)
q
%*+Xg+=,+
%+vXg,
’(((()*++++,
(7)
假设地面运动加速度为频谱密度为S0的白噪声随机过程,则由结构随机地震反应理论,可以求得主结构相对地面位移反应的方差解析式为
22
(!S0/\"3)[(B0/A0)(A2A3-A1A4)+A(-2B0B2)+A1B2]3B012
,%=
(A(-A0A2)1A2A3-A1A4)32
y
(8)
(2+4’1’d(;A3=2’1+2’d((;A4=1;其中:A0=&2;A1=2’1(2+2’d(;A2=1+(1+$)(1+$)
B0=-(1+$)(1+$)(2;B1=-2’d(;B2=-1.
TMD与主结构质量比$一般在[0.5%,3%]之间变化,频率比(变化范围一般为[0.4,1.2],TMD系统
的阻尼比’d变化范围一般为[0,0.3].对质量比$选搜索2位有效数字,用来表示其的二进制编码长度为12位;频率比(选搜索4位有效数字,用来表示其的二进制编码长度为13位;TMD系统的阻尼比’d选搜索4位有效数字,用来表示其的二进制编码长度为13位;因此表征控制参数的一个染色体长度为37.种群个数取为100,进化代数取100代,染色体交叉概率取为0.9,变异概率为0.15.由于主结构相对地面位移反应的方差越小,表明TMD控制效果越好,所以取主结构相对地面位移反应的方差%y的倒数作为遗传算法的适应度函数.遗传算法求得的影响单TMD在地震荷载作用下振动控制的最优参数为:质量比$=0.03,阻尼比
频率比(=0.9592.’d=0.0856,2
-570-
河南科学
第26卷第5期
为了验证上述方法,本文以一栋30层高层结构为例进行仿真分析.结构每层质量为1.0×106kg;每层刚度为1.2×1010N/m;每层阻尼为2.0×105N・s/m.结构受El-Centro(NS方向)波作用,按照上述方法可得TMD参数为:TMD质量为9.0×105kg;刚度为2.635×107N/m;阻尼为4.1686×105N・s/m.结构在安装经优化过的TMD系统后与无控时的控制效果比较可见图6和图7.
图6结构顶层位移控制效果比较图
Fig.6
Comparisonofstructuretoplayerdisplacementcontroleffect
Fig.7
图7结构顶层加速度控制效果比较图
Comparisonofstructuretoplayeraccelerationcontroleffect
5结论
本文对遗传算法(GA)在建筑结构系统辨识、振动控制算法、结构振动控制作动器优化布置和作动器参数优化等研究方向的应用进行了较全面研究和计算机仿真分析.仿真分析表明,遗传算法是一种高效的搜索计算和优化计算方法,这一算法非常适合用来处理振动控制领域中的诸多复杂优化问题,它能够精确识别复杂结构系统的动态响应,能够实现结构振动控制系统作动器的优化布置和设计参数寻优.参考文献:
[1]HollandJH.Adaptationinnaturalandartificialsystems[M].AnnArbor:TheUniversityofMichiganPress,1975.[2]GoldbergDE.Geneticalgorithms[M].Addison-WesleyPublishingCompany,1989.[3]李人厚.智能控制理论和方法[M].西安:西安电子科技大学出版社,1999.[4]顾仲权,马扣根,陈卫东.振动主动控制[M].北京:国防工业出版社,1997.[5]何玉敖,郭[6]何玉敖,郭
婷.基于遗传算法的结构主动控制[J].振动工程学报,1999,12(2):182-187.婷.遗传算法-神经网络结构控制系统研究[J].振动工程学报,2001,14(2):192-195.
[7]何玉敖,何亚东.基于Lyapunov稳定性原理和遗传算法的结构半主动控制[J].土木工程学报,2000,33(6):88-93.[8]佟道林.基于遗传算法的结构系统辨识、半主动控制即控制优化研究[D].天津:天津大学硕士学位论文,2001.
ApplicationofGeneticAlgorithmsinStructureVibrationControl
ZhaoLin1,ZhangMeng2
(1.DepartmentofZhengzhouRealEstate,ShanghaiGreenlandGroup,Zhengzhou450046,China;
2.CollegeofCivilEngineering,ZhengzhouUniversity,Zhengzhou,450002,China)
Abstract:Theapplicationofgeneticalgorithmsinstructurevibrationcontrolisfullyinvestigated.Thefollowingfields,buildingstructuresystemidentification,vibrationcontrolalgorithm,optimalactuatorconfigurationandactuatorparameters,areinvolvedintheanalysis.Computingsimulationresultsindicatethatgeneticalgorithmsarehighlyeffective.Themethodadaptswelltodisposeofthecomplexproblemsinstructurevibrationcontrol.Itcandistinguishthedynamicresponseofcomplicatedstructuresystem’sprecisely,andbehelpfulintheoptimizationdesignofstructurevibrationcontrolsystem.
Keywords:geneticalgorithms;systemidentification;optimizationdesign;actuator;controlalgorithm
因篇幅问题不能全部显示,请点此查看更多更全内容