由于H ∞控制器中Riccati方程较难求解,本文通过利用李雅普诺夫稳定性理论及Matlab的LMI工具箱推导了一种LMI控制算法,并建立了一栋20层结构LMI最优主动控制的分析模型,并利用Matlab语言编制了求解软件,通过对El-Centro波下结构动力反应进行数值模拟,从楼层位移、速度、加速度,层间位移,AMD作动器控制力输出方面对LMI控制算法效果进行分析。计算结果表明:LMI控制能有效降低建筑结构楼层位移、速度、加速度、层间位移响应,但需要付出一定的控制成本。可以通过调整权重参数达到目标减振效果,权重越大,控制效果越好,控制成本越高。顶层控制力最大,当时减振效果较好,顶层控制力幅值为楼面重力荷载27.8%。 Because Riccati equation in H ∞ controller is more difficult to solve, in this paper, by using Lyapunov stability theory and LMI toolbox of Matlab, an LMI control algorithm was deduced, an analysis model of LMI optimal active control for a 20-story structure was established, and Matlab language was used to compile the solution software. The structural dynamic response under El-Centro wave was numerically simulated, and the effect of LMI control algorithm was analyzed from floor displacement, velocity, acceleration, inter-floor displacement and AMD actuator control output side. The calculation results show that the LMI control can effectively reduce the floor displacement, velocity, acceleration and inter-floor displacement response of the building structure, but it needs to pay a certain control cost. The target damping effect can be achieved by adjusting the weight parameters. The greater the weight is, the better the control effect is and the higher the control cost is. The top floor control force is the largest, and the damping effect is better when ; the amplitude of the top floor control force is 27.8% of the floor gravity load.
由于H∞控制器中Riccati方程较难求解,本文通过利用李雅普诺夫稳定性理论及Matlab的LMI工具箱推导了一种LMI控制算法,并建立了一栋20层结构LMI最优主动控制的分析模型,并利用Matlab语言编制了求解软件,通过对El-Centro波下结构动力反应进行数值模拟,从楼层位移、速度、加速度,层间位移,AMD作动器控制力输出方面对LMI控制算法效果进行分析。计算结果表明:LMI控制能有效降低建筑结构楼层位移、速度、加速度、层间位移响应,但需要付出一定的控制成本。可以通过调整权重参数 α 达到目标减振效果,权重 α 越大,控制效果越好,控制成本越高。顶层控制力最大,当 α = 1 e 1 2 时减振效果较好,顶层控制力幅值为楼面重力荷载27.8%。
建筑与土木工程,主动振动控制,龙格–库塔方法,状态变量,线性矩阵不等式,李雅普诺夫函数
Bingda Huang, Nan Ge, Yong Wang
Architectural Engineering Institute, North China University of Science and Technology, Tangshan Hebei
Received: Oct. 1st, 2020; accepted: Oct. 9th, 2020; published: Oct. 26th, 2020
Because Riccati equation in H∞ controller is more difficult to solve, in this paper, by using Lyapunov stability theory and LMI toolbox of Matlab, an LMI control algorithm was deduced, an analysis model of LMI optimal active control for a 20-story structure was established, and Matlab language was used to compile the solution software. The structural dynamic response under El-Centro wave was numerically simulated, and the effect of LMI control algorithm was analyzed from floor displacement, velocity, acceleration, inter-floor displacement and AMD actuator control output side. The calculation results show that the LMI control can effectively reduce the floor displacement, velocity, acceleration and inter-floor displacement response of the building structure, but it needs to pay a certain control cost. The target damping effect can be achieved by adjusting the weight parameters. The greater the weight is, the better the control effect is and the higher the control cost is. The top floor control force is the largest, and the damping effect is better when α = 1 e 1 2 ; the amplitude of the top floor control force is 27.8% of the floor gravity load.
Keywords:Architecture and Civil Engineering, Active Vibration Control, Runge-Kutta Method, State Variables, Linear Matrix Inequality, Lyapunov Function
Copyright © 2020 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
H∞控制是在H∞空间 [
系统的状态空间模型 [
x ˙ = A x + B u (1)
其中x是状态向量,u是控制输入,A和B分别是恰当维数的已知的常数矩阵。若采用状态反馈,则有: u = K x 。设计LMI控制器的目的就是使闭环系统 x ˙ = A x + B u = ( A + B K ) x 具有渐进稳定性,同时控制器具有鲁棒性。
考虑Lyapunov函数 [
( A + B K ) T P + P ( A + B K ) − P < 0 (2)
且 P = P T 。
对(2)式两端各乘以 P − 1 ,得:
P − 1 ( A + B K ) T + ( A + B K ) P − 1 − P − 1 < 0 (3)
令 X = P − 1 , Y = K P − 1 , K = Y P = X − 1 得:
( A X + B Y ) T + ( A X + B Y ) − X < 0 (4)
因此,根据Schusr补引理 [
[ ( A X + B Y ) T + ( A X + B Y ) 0 0 − X ] < 0 (5)
即上述的闭环系统是稳定的,由于(5)式是一个线性矩阵不等式 [
x ˙ = A x + B u + W x ¨ g (6)
就可求得LMI控制下的运动方程的地震动力反应,检验所设计的LMI控制器的效果。
因此,求取LMI控制算法可行解的过程如下:
setlmis([
X = lmivar(1,[2 1]); %定义矩阵变量
Y = lmivar(2,[1 2]); %定义矩阵变量
lmiterm([1 1 1 X],A,1,’s’); %描述线性矩阵不等式
lmiterm([1 1 1 Y],B,1,’s’); %描述线性矩阵不等式
lmiterm([1 2 2 X],-1,1); %描述线性矩阵不等式
lmis = getlmis; %完成写线性矩阵不等式
[tmin,xfeas] = feasp(lmis,[0,0,10,0,0],-1)%求取可行解
x = dec2mat(lmis,xfeas,X) %得到可行矩阵变量值
y = dec2mat(lmis,xfeas,Y)
K = y*inv(x) %求状态反馈控制器增益
对于被控结构的地震动力反应计算预测,则可采用的数值预测方法有Newmark法、经典四阶Runge-Kuta法等。
为评价LMI算法对建筑结构的减震效果,在这里为节约评价成本采用数值预测的方法 [
图1. 结构模型平面图
图2. 前三阶振型图
建立一栋20层建筑结构模型,其模型平面图如图1所示,为方便数值模拟将其简化为20个集中质量的多自由度体系。结构每个集中质量mi = 2933 t,每层层间侧移刚度ki = 28,950,000 kN/m,其一、二、三阶振型如图2所示。阻尼矩阵采用振型阻尼比为0.05时的Rayleigh阻尼的形式 [
从图3、图4可以看出,当 α = 1e10 - 1e12 时,LMI控制算法减振效果显著。通过调节系数 α 可以达到调节目标减震效果的效果,随着调节系数 α 相应的位移、层间位移、速度、加速度减震效果也随之增加,反之,随着调节系数 α 的减小,相应的控制效果也随之减弱,在实际工程设计中,可根据目标建筑结构抗震设防等级选取相应级别的地震波以通过不断调节参数试算以确定其参数 α 。
图3. 结构地震动力反应时程
图4. 调节系数 α 值与结构地震动力反应最大值
层间 (楼层) | α = 1e11 | α = 1e12 | 无控制 | |||||||
---|---|---|---|---|---|---|---|---|---|---|
层间位移/mm | 效率 /% | 楼层加速度/m·s−2 | 效率 /% | 层间位移/mm | 效率 /% | 楼层加速度/m·s−2 | 效率 /% | 层间位移/mm | 楼层加速度/m·s-2 | |
0~1(1) | 16.29 | 31.06 | 3.02 | −15.71 | 5.33 | 77.44 | 3.34 | −27.97 | 23.63 | 2.61 |
1~2(2) | 16.06 | 31.22 | 5.51 | −7.20 | 5.02 | 78.50 | 5.64 | −9.73 | 23.35 | 5.14 |
2~3(3) | 15.64 | 31.64 | 7.73 | −4.46 | 4.59 | 79.94 | 6.72 | 9.19 | 22.88 | 7.40 |
3~4(4) | 15.09 | 32.27 | 8.97 | 4.68 | 4.24 | 80.97 | 7.03 | 25.29 | 22.28 | 9.41 |
4~5(5) | 14.48 | 33.02 | 9.65 | 15.13 | 4.37 | 79.79 | 7.08 | 37.73 | 21.62 | 11.37 |
5~6(6) | 13.89 | 33.32 | 9.96 | 22.91 | 4.48 | 78.49 | 7.57 | 41.41 | 20.83 | 12.92 |
6~7(7) | 13.19 | 33.48 | 10.02 | 28.38 | 4.56 | 77.00 | 7.79 | 44.32 | 19.83 | 13.99 |
7~8(8) | 12.31 | 33.78 | 10.51 | 28.36 | 4.60 | 75.26 | 7.74 | 47.24 | 18.59 | 14.67 |
8~9(9) | 11.26 | 34.34 | 11.97 | 25.74 | 4.60 | 73.18 | 7.71 | 52.17 | 17.15 | 16.12 |
9~10(10) | 10.38 | 33.33 | 12.79 | 26.03 | 4.53 | 70.91 | 7.43 | 57.03 | 15.57 | 17.29 |
10~11(11) | 9.99 | 28.95 | 12.98 | 27.16 | 4.47 | 68.21 | 7.32 | 58.92 | 14.06 | 17.82 |
11~12(12) | 9.47 | 27.21 | 12.60 | 29.01 | 4.52 | 65.26 | 7.01 | 60.51 | 13.01 | 17.75 |
---|---|---|---|---|---|---|---|---|---|---|
12~13(13) | 8.89 | 27.07 | 11.81 | 31.18 | 4.31 | 64.64 | 6.68 | 61.07 | 12.19 | 17.16 |
13~14(14) | 8.31 | 26.53 | 10.63 | 34.30 | 3.90 | 65.52 | 6.58 | 59.33 | 11.31 | 16.18 |
14~15(15) | 7.65 | 25.44 | 10.47 | 31.92 | 3.65 | 64.42 | 6.96 | 54.75 | 10.26 | 15.38 |
15~16(16) | 6.83 | 24.20 | 10.57 | 32.42 | 3.31 | 63.26 | 6.59 | 57.86 | 9.01 | 15.64 |
16~17(17) | 5.80 | 23.08 | 12.16 | 25.26 | 2.88 | 61.80 | 6.65 | 59.13 | 7.54 | 16.27 |
17~18(18) | 4.56 | 22.05 | 13.82 | 22.88 | 2.33 | 60.17 | 6.70 | 62.61 | 5.85 | 17.92 |
18~19(19) | 3.14 | 21.50 | 14.98 | 22.26 | 1.66 | 58.50 | 7.38 | 61.70 | 4.00 | 19.27 |
19~20(20) | 1.60 | 21.18 | 15.56 | 22.32 | 0.92 | 54.68 | 9.40 | 53.07 | 2.03 | 20.03 |
表1. 控制方式与地震动力反应
注:某一控制方案的减震效率是指相对于无控制方案而言。
结合表1及图4可以看出,采用LMI控制后,随着调节系数 α 的增大,楼层加速度减震效果也随之增加,但底部结构的减震效果较不明显甚至出现增加建筑结构楼层加速度的情况,且在相同设计参数 α ,从而能够在一定程度上降低工程试算难度;在统一设计参数 α 下,作动器控制成本随着楼层的提高而提高,因此,最大作动器控制力在顶层,取 α = 1e12 时,顶层控制力的幅值达到8000 kN左右,为本层楼面重量(或等效重力荷载代表值)的27.8%。
从图5可以看出,随着设计参数的增大(控制效果越好),每一楼层的作动器控制力也随之增加即意味着控制成本的增加;与LQR、IOC等经典控制算法相比,LMI算法只有一个设计参数,从而能够在一定程度上降低工程试算难度;在统一设计参数下,作动器控制成本随着楼层的提高而提高,因此,最大作动器控制力在顶层,取时,顶层控制力的幅值达到8000 kN左右,为本层楼面重量(或等效重力荷载代表值)的27.8%。
图5. 调节系数 α 值与楼层控制力幅值
1) LMI控制有显著减振效果,调节系数 α 越大,控制效果越好,控制成本越高,调节系数 α 越小,控制效果越差,控制成本越低。
2) LMI控制,顶层控制力最大,当 α = 1e12 时减振效果较好时,顶层控制力幅值为楼面重力荷载27.8%。
黄炳达,葛 楠,王 永. 一种LMI建筑结构振动控制算法An LMI Vibration Control Algorithm for Building Structures[J]. 动力系统与控制, 2020, 09(04): 225-231. https://doi.org/10.12677/DSC.2020.94022