半主动控制凭借其兼备主动控制的自适应性和被动控制的稳定性等优势,在多层建筑振动控制领域中得到普遍应用[1].当前神经网络(NN)中的前馈控制主要采取反向传播的形式[2],然而这种形式为了调谐每个网络权重需要进行繁琐的学习,导致计算效率低下.对此本文在前馈控制中引入径向基函数网络(RBFN)[3],该网络只需一个隐层,隐层的传递函数是一个非线性半仿射函数,因此RBFN的学习效率要明显快于反向传播网络,且RBFN可以逼近任意非线性连续函数,并能消除局部最小化.
本文基于MGA提出了自适应神经网络控制(MGA_ANNC)策略,并设计了适合MRD的MGA_ANNC/MCO半主动控制策略.分别采用MGA_ANNC/MCO,MGA_ANNC和LQG对一座9层钢筋混凝土框架结构进行了振动控制,得出了一些结论,可为类似结构的减震设计提供参考.
1 控制装置与受控结构的力学模型 1.1 控制装置的力学模型本文采用调谐质量-磁流变阻尼器(TM-MRD)作为半主动控制装置,采用S-Bouc-Wen模型[4]模拟MRD的力学性能.
TM-MRD的控制力向量U由MRD的阻尼力向量Fb与弹簧的弹性恢复力向量Fc组成,即:
式中:kt为弹簧的恢复刚度;xt为TM-MRD与主结构的相对位移向量. 1.2 受控结构的力学模型受控多层剪切型框架结构的运动方程为
式中:M,C和K分别为质量、阻尼和刚度矩阵,,和x分别为加速度、速度和位移向量,且x=[x1,x2,…,xn]T;ω是影响系数矩阵;为地震加速度向量;U(t)为控制力向量.将式(2)两边都除以M可得
其中: 2 MGA_ANNC算法的设计 2.1 非线性动力系统的参考模型单输入/单输出的n(n≥1)阶系统可表示为
式中:x=[x1,x2,…,xn]T∈Rn为系统的状态向量;u∈R为控制信号;p(x)和g(x)均为光滑非线性函数;d(t)为有界的外部干扰.对输出进行n次微分,直到控制输入u出现,此时可获得式(5)的输入/输出形式为
假设1h(x)为越过紧集η⊂Rn且远离0的有界函数,即:
如果控制目标是让系统输出y去追踪参考轨迹yr,那么参考控制输入r可表示为
对于∀x∈S,如果p(x)和h(x)已知,假设1被满足的话,那么控制律可被定义为
将式(9)代入式(5),则线性化系统可变为
如果定义e=yr-y为追踪误差,那么参考控制输入(式(8))满足以下误差方程:
2.2 光滑控制输入在滑动模态控制系统[5]的开关控制输入中,易出现抖振.本文通过光滑与开关曲面邻近的薄边界层的不连续性,达到消除抖振的目的,其表达式为
式中μ为边界层的厚度. 2.3 适合非线性系统的MGA_ANNC算法由于标准遗传算法(GA)[6]的局部寻优能力较差,故作者提出的改进遗传算法(MGA)为
式中:m为隐节点的总和;为隐节点和输出节点间的可调网络权重,且满足||≤wmax;R(S)=[R1(S),R2(S) ,…,Rm(S)]T为径向基函数向量,它被定义为 式中,βk,Ak和λ分别为第k个隐单元的径向基函数、中心和宽度.假设2当x∈η⊂Rn,存在一个可调参数向量允许模糊系统(S,)= ×R(S)被近似为一个连续函数u,即:
在t时刻的估计采用表示,此时可定义控制输出的估计(S,)为
定义参数误差向量为we=-,则
定义建模误差为
进而得到:
将式(20)代入式(11)可获得误差动力方程为
定义增广误差为
选择式(22)中的λn-1,…,λ1,λ0和式(21)中的φn-1,…,φ1,φ0,使其满足:
如果定义es=[e,…,e(n-1)]T作为式(21)的状态向量,则式(17)又可表示为
式中:
当(l)为严格正实函数时,存在对称正定矩阵O和Q满足:
定义Lyapunov候选函数[7]为
式中:we=[we1,we2,…,wem]T;L11=h(x)·Im×m;τ是表示学习速率的正常数.如果eTsOes>μ2,那么存在:
则可推出,进而得到:
本文提出的改进自适应律为
式中,sat(Se/μ)为薄边界层函数,其表达式为将式(32)代入式(31),则式(31)可变为
当|Se|>μ,则
如果eTsOes>μ2,μ>0和σ>0,则有
由式(34)易推得≤-τ·‖es‖·σ,从而得出<0.故在所有的紧集η内,V将逐渐收敛至0.
3 半主动控制算法的设计为了解决限幅最优(CO)控制[8]的命令电压vi只能取0或的弊端,本文提出了MCO半主动控制算法,其施加于第i个MRD的vi可表示为
式中:H(·)为Heaviside阶跃函数;Fi为实测控制力;ζi为电压系数;为施加到电流驱动器的电压;Fci为理论控制力;Fmax为MRD的最大出力. 4 数值分析 4.1 系统参数以一座9层钢筋混凝土框架结构(以下简称结构)为算例,各楼层质量mj=843.6 t,各楼层刚度kj=1 352.9 MN·m-1.TM-MRD布置在结构顶层,质量块的质量mt=8.216×104 kg,kt=9.372×106 N·m-1.对于MRD,其标定的最大出力和最大冲程L分别为1 000 kN和±40 cm.主动控制装置采用压电陶瓷驱动器(PCA),其和L的取值与MRD相同.S=dTet=36e9-0.3et+ė9-0.6ėt,=[1,0.86,0.65,0.26,0.06,0,-0.08,-0.25,-0.38,-0.76,-0.96],μ=0.2,τ=10,ζi=1.12.
4.2 结果分析从表 1可以看出,无控制结构的弹塑性层间位移角峰值(J1)和弹塑性层间位移角均方根(J2)均超出了结构的弹性层间位移角限值(1/550 rad)[9],MGA_ANNC/MCO和MGA_ANNC对J1和J2的控制效果比LQG的都要好;MGA_ANNC/MCO对结构的加速度峰值(J3)、加速度均方根(J4)和基底剪力峰值(J5)的控制效果也要比LQG的显著;MGA_ANNC/MCO,MGA_ANNC和LQG控制下的结构曲率延性系数峰值(J6)、曲率延性系数均方根(J7)、构件耗能量(J8)和塑性铰个数(J9)等指标比无控制结构的均要有所降低.MGA_ANNC/MCO和MGA_ANNC的控制装置实际最大出力(J10)都接近其标定的最大出力;MGA_ANNC/MCO和MGA_ANNC的控制装置实际冲程(J11)均分别接近其标定的最大冲程.MGA_ANNC/MCO的性能接近MGA_ANNC;MGA_ANNC/MCO策略的所需外界能量(J12)比MGA_ANNC和LQG的都要少很多;MGA_ANNC/MCO,MGA_ANNC和LQG等策略所需的控制装置数目(J13)相等;MGA_ANNC/MCO和MGA_ANNC所需的传感器数目(J14)比LQG的要多;MGA_ANNC/MCO和MGA_ANNC所需的计算资源量(J15)均为零.
由图 1~图 3可以看出,MGA_ANNC/MCO对顶层加速度和弹塑性层间位移角峰值的控制效果要比LQG的显著;MGA_ANNC/MCO与MGA_ANNC的控制力时程很吻合且均无抖振出现.由表 2知,MGA_ANNC/MCO和MGA_ANNC控制下的结构响应值均明显小于相应的LQG的值,故前两者的鲁棒性要明显优于后者;MGA_ANNC的J12仍比MGA_ANNC/MCO的要大得多.
MGA_ANNC/MCO策略的减震效果要明显优于LQG策略;MGA_ANNC/MCO的实际最大出力和最大冲程均比LQG的要大得多,且前者无需任何计算资源;MGA_ANNC/MCO的鲁棒性要明显优于LQG;MGA_ANNC/MCO所需的外界能量要比MGA_ANNC的少得多,两者的控制力时程非常吻合且均无抖振发生.
[1] | Bitaraf M,Ozbulut O E,Hurlebaus S,et al.Application of semi-active control strategies for seismic protection of buildings with MR dampers[J].Engineering Structures,2010,32(10):3040-3047.(1) |
[2] | Bani-Hani K A,Sheban M A.Semi-active neuro-control for base-isolation system using magnetorheological (MR) dampers[J].Earthquake Engineering and Structural Dynamics,2006,35(9):1119-1144.(1) |
[3] | Qasem S N,Shamsuddin S M.Radial basis function network based on time variant multi-objective particle swarm optimization for medical diseases diagnosis[J].Applied Soft Computing,2011,11(1):1427-1438.(1) |
[4] | Zapateiro M,Karimi H R,Luo N,et al.Real-time hybrid testing of semiactive control strategies for vibration reduction in a structure with MR damper[J].Structural Control and Health Monitoring,2010,17(4):427-451.(1) |
[5] | Wu L,Zheng W X.Passivity-based sliding mode control of uncertain singular time-delay systems[J].Automatica,2009,45(9):2120-2127.(1) |
[6] | Peteghem V V,Vanhoucke M.A genetic algorithm for the preemptive and non-preemptive multi-mode resource-constrained project scheduling problem[J].European Journal of Operational Research,2010,201(2):409-418.(1) |
[7] | Diehl M,Amrit R,Rawlings J B.A Lyapunov function for economic optimizing model predictive control[J].IEEE Transactions on Automatic Control,2011,56(3):703-707.(1) |
[8] | Yoshida O,Dyke S J.Response control of full-scale irregular buildings using magnetorheological dampers[J].Journal of Structural Engineering,2005,131(5):734-742.(1) |
[9] | 中国建筑科学研究院.GB 50011—2010 建筑抗震设计规范[S].北京:中国建筑工业出版社,2010. (China Academy of Building Research.GB 50011-2010 Code for seismic design of buildings[S].Beijing:China Architecture and Building Press,2010.)(1) |