工程科学学报，第41卷，第8期：1092-1102,2019年8月 Chinese Journal of Engineering,Vol.41,No.8:1092-1102,August 2019 DOI:10.13374/j.issn2095-9389.2019.08.015;http://journals.ustb.edu.cn 基于滚轴支座基础智能隔震结构的非光滑主动控制 刘彦辉)四，黄襄云)，黄小芳)，王清) 1)广州大学工程抗震研究中心，广州5104052)广州大学机械与电气工程学院，广州510006 ☒通信作者，E-mail:Liuyanhui2012@163.cm 摘要以滚轴支座基础隔震结构作为受控结构研究对象，在该隔震结构的隔震层施加主动控制装置，形成智能隔震体系， 以控制隔震层的位移，提高结构的安全性.在智能隔震结构中引入非光滑控制算法，基于隔震层位移和速度反馈，提出了智能 隔震结构的非光滑控制算法，进一步根据Lyapunov稳定理论，推导了在非光滑控制下智能控制闭环系统的全局有限时间稳定 性.结合一栋六层滚轴支座基础隔震结构，进行了非光滑主动控制算法下和LQG主动控制算法下的地震响应控制仿真分析. 结果表明，智能隔震结构可有效控制结构的位移，与被动隔震结构相比较上部结构的地震响应有一定程度的减小，同时提出 的非光滑控制算法与LQG控制算法相比较具有更好的控制效果，相比LQG控制算法通过较少的反馈量即可实现反馈控制， 且非光滑控制算法具有良好的稳定性. 关键词智能隔震；非光滑控制；LQG控制；稳定性；地震响应 分类号TU352.12 Nonsmooth active control method for base-smart isolated structures with roller bearings LIU Yan-hui,HUANG Xiang-yun,HUANG Xiao-fang,WANG Qing? 1)Earthquake Engineering Research Test Center,Guangzhou University,Guangzhou 510405,China 2)School of Mechanical and Electrical Engineering,Guangzhou University,Guangzhou 510006,China Corresponding author,E-mail:Liuyanhui2012@163.com ABSTRACT In applied structural control technology,base-isolated technology has become popular due to its advantage of simple shock absorption,stable performance,and reasonable cost.Currently,base isolation is extensively applied worldwide,and its role in mitigating the seismic response of structures continues to grow.Moreover,it has been proven effective in decreasing seismic response of structures under recent strong earthquakes.However,the displacement at the isolation layer is sometimes large under strong earth- quakes,which will decrease the safety of the structure and perhaps lead to the failure of the isolation layer.Therefore,in this study, the base-isolated structure with roller bearings is taken to investigate the seismic response control of structures,and the active control devices are added in the isolation layer of the isolated structure to decrease the seismic displacement at the isolation layer,so that a smart-isolated structure is formed.Nonsmooth control algorithm is introduced in the smart-isolated structure.Based on the feedback of the velocity and displacement of the isolation layer,nonsmooth control algorithm is proposed for designing the smart-isolated structure. Moreover,according to Lyapunov stable theory,the global finite time stability of intelligent control closed-loop system under nonsmooth control is deduced.A six-layer isolated structure with roller bearings is used as an example,and a simulation analysis of seismic re- sponse control is performed based on the nonsmooth active control algorithm and linear quadratic Gaussian (LQG)active control algo- rithm.The results show that the smart-isolated technology can effectively control seismic displacement at the isolation layer,and com- pared with the passive isolated technology,the superstructure seismic response is significantly decreased.Meanwhile,the results dem- onstrate that compared with the LQG control algorithm,the nonsmooth control algorithm has a better control effect and can implement 收稿日期：2018-07-27 基金项目：国家自然科学基金资助项目(51778163)：国家重点研发计划重点专项资助项目(2016Y下E0127600):广州市属高校科研项目重点资 助项目(1201610135)：广东省教育厅创新团队资助项目(2016 KCXTD016)
工程科学学报,第 41 卷,第 8 期:1092鄄鄄1102,2019 年 8 月 Chinese Journal of Engineering, Vol. 41, No. 8: 1092鄄鄄1102, August 2019 DOI: 10. 13374 / j. issn2095鄄鄄9389. 2019. 08. 015; http: / / journals. ustb. edu. cn 基于滚轴支座基础智能隔震结构的非光滑主动控制 刘彦辉1) 苣 , 黄襄云1) , 黄小芳1) , 王 清2) 1) 广州大学工程抗震研究中心, 广州 510405 2) 广州大学机械与电气工程学院, 广州 510006 苣通信作者,E鄄mail:Liuyanhui2012@ 163. com 摘 要 以滚轴支座基础隔震结构作为受控结构研究对象,在该隔震结构的隔震层施加主动控制装置,形成智能隔震体系, 以控制隔震层的位移,提高结构的安全性. 在智能隔震结构中引入非光滑控制算法,基于隔震层位移和速度反馈,提出了智能 隔震结构的非光滑控制算法,进一步根据 Lyapunov 稳定理论,推导了在非光滑控制下智能控制闭环系统的全局有限时间稳定 性. 结合一栋六层滚轴支座基础隔震结构,进行了非光滑主动控制算法下和 LQG 主动控制算法下的地震响应控制仿真分析. 结果表明,智能隔震结构可有效控制结构的位移,与被动隔震结构相比较上部结构的地震响应有一定程度的减小,同时提出 的非光滑控制算法与 LQG 控制算法相比较具有更好的控制效果,相比 LQG 控制算法通过较少的反馈量即可实现反馈控制, 且非光滑控制算法具有良好的稳定性. 关键词 智能隔震; 非光滑控制; LQG 控制; 稳定性; 地震响应 分类号 TU352郾 12 收稿日期: 2018鄄鄄07鄄鄄27 基金项目: 国家自然科学基金资助项目(51778163);国家重点研发计划重点专项资助项目(2016YFE0127600);广州市属高校科研项目重点资 助项目(1201610135);广东省教育厅创新团队资助项目(2016KCXTD016) Nonsmooth active control method for base鄄smart isolated structures with roller bearings LIU Yan鄄hui 1) 苣 , HUANG Xiang鄄yun 1) , HUANG Xiao鄄fang 1) , WANG Qing 2) 1) Earthquake Engineering Research & Test Center, Guangzhou University, Guangzhou 510405, China 2) School of Mechanical and Electrical Engineering, Guangzhou University, Guangzhou 510006, China 苣Corresponding author, E鄄mail: Liuyanhui2012@ 163. com ABSTRACT In applied structural control technology, base鄄isolated technology has become popular due to its advantage of simple shock absorption, stable performance, and reasonable cost. Currently, base isolation is extensively applied worldwide, and its role in mitigating the seismic response of structures continues to grow. Moreover, it has been proven effective in decreasing seismic response of structures under recent strong earthquakes. However, the displacement at the isolation layer is sometimes large under strong earth鄄 quakes, which will decrease the safety of the structure and perhaps lead to the failure of the isolation layer. Therefore, in this study, the base鄄isolated structure with roller bearings is taken to investigate the seismic response control of structures, and the active control devices are added in the isolation layer of the isolated structure to decrease the seismic displacement at the isolation layer, so that a smart鄄isolated structure is formed. Nonsmooth control algorithm is introduced in the smart鄄isolated structure. Based on the feedback of the velocity and displacement of the isolation layer, nonsmooth control algorithm is proposed for designing the smart鄄isolated structure. Moreover, according to Lyapunov stable theory, the global finite time stability of intelligent control closed鄄loop system under nonsmooth control is deduced. A six鄄layer isolated structure with roller bearings is used as an example, and a simulation analysis of seismic re鄄 sponse control is performed based on the nonsmooth active control algorithm and linear quadratic Gaussian (LQG) active control algo鄄 rithm. The results show that the smart鄄isolated technology can effectively control seismic displacement at the isolation layer, and com鄄 pared with the passive isolated technology, the superstructure seismic response is significantly decreased. Meanwhile, the results dem鄄 onstrate that compared with the LQG control algorithm, the nonsmooth control algorithm has a better control effect and can implement
刘彦辉等：基于滚轴支座基础智能隔震结构的非光滑主动控制 ·1093· feedback control for base-isolated structures by using fewer feedbacks.Furthermore,the nonsmooth control algorithm has great stability KEY WORDS smart isolation;nonsmooth control:linear quadratic Gaussian (LQG)control;stability;seismic response 随着建筑结构往高层化和复杂化的发展，在地 目前已经应用于航天姿态控制[]、高精度制导率设 震波激励下的建筑结构振动响应和安全话题已经引 计)、永磁同步机位置控制3]和轮式机器人4)等 起人们的关注.传统的抗震结构通过增加抗侧刚度 领域. 来减小地震的影响，但结构刚度越大，地震响应越 本文首先针对基础智能隔震结构，基于隔震层 大，这类结构体系非但不经济，且安全性无法保证. 位移和速度反馈，推导了智能隔震结构的非光滑控 对结构进行振动控制一直是土木工程领域重要的课 制算法，然后对该智能控制系统在地震波激励下进 题.为了解决这样的工程抗震问题，人们开始采用 行仿真分析，对比分析了基础隔震、LQG控制和非 结构振动控制方式来抑制结构振动的响应，从而提 光滑控制下的仿真结果，验证了非光滑控制算法的 高和改善结构的抗震性能.有学者提出了结构振动 有效性和稳定性 控制概念]，结构振动控制作为有效抑制动力荷载 1受控结构方程 作用下结构振动响应的方法.其中，基础隔震技术 在实际工程中的应用较为广泛，许多国家已将基础 1.1受控结构运动方程 隔震技术纳入抗震规范.基础隔震技术能有效减少 采用层剪切模型，n层的建筑结构简化为一维 上部结构地震响应，然而会使隔震层产生较大的位 具有n个自由度的结构模型，见图1，图中k,k1,k2, 移，这限制了隔震技术的广泛应用.因此采用主动 …,k。为层刚度；mb,m1,m2,…,mn为层质量.在单 控制、半主动控制方式来完善隔震结构的整体性能. 向水平地震波作用下，受控结构的运动方程为： 对于智能隔震控制系统而言，核心是控制算法，将控 MD+CD +F(D,D)=-MI,+EU (1) 制算法载入于计算机中，根据设置在结构上的传感 式中，D=[d。d,d2…dn]为结构位移向 器采集的状态信号，反馈于已经载入计算机中所设 量，其中d为隔震层位移，d,为结构第n层相对地 计的控制算法，从而获得控制力，然后通过伺服系统 将控制力转为驱动力施加给受控结构，以削弱结构 面的位移：D、D分别为速度矢量和加速度矢量；x。 的振动反应.具有代表性的控制算法有经典最优控 为地面地震加速度；M和C分别为结构质量和阻尼 制(linear quadratic regulator,LQR)[2】,线性二次型 矩阵；F(D,D)=KD+L为结构恢复力向量，K为 高斯控制(linar quadratic Guassican,LQG)[2),极点 刚度矩阵，其中L为塑性恢复力向量：I= 配置算法[s),瞬时最优控制(instantaneous optimal [11…1]x;E为位置矩阵；0= control,IOC)[等.然而，传统控制算法对参数摄动 [山14…u.]T为控制力输人，山n为第n个控 和外部干扰过于灵敏，发展和应用对系统参数影响 制力. 小且外干扰变化灵敏度更小、鲁棒性更强的结构振 1.2阻尼的选择 动控制算法就显得极为重要.滑膜控制[5-6]是一种 对于一维多自由度体系，上部结构为弹性，因此 非连续性的反馈控制，可以使得闭环系统在较短时 选用瑞雷阻尼(Rayleigh damping)来组合阻尼矩阵. 间内收敛，并具有较好鲁棒性的特性，但相对于系统 C=cM+qK (2) 状态变量而言是非连续的，控制过程中的高速切换 式中，c、9为比例系数，分别表示为： 容易产生抖振现象. C= 设计一种概念清晰、简单有效的反馈控制策略， 21c(w,5:-10) 0-0 (3) 使得闭环系统在有限时间内收敛且具备强鲁棒性的 20:0 特性，又不会产生由于非连续性控制引起的抖振现 q=w- 0: 象，将在工程应用中具有较好的发展前景，而非光滑 其中，专：和专，分别为结构第i和j阶振型阻尼比，0：、 控制正是这样一种算法.随着齐次性方法和Lya- 心，分别表示结构第ij阶频率.专：、专随着结构振幅 punov理论的研究和发展7-]，介于连续和非连续之 和材料的不同而发生变化，但受到频率的影响较小 间的非光滑控制有了一定的研究和应用.相比于传 在进行结构反应分析时，各阶的阻尼比只与结构材 统控制算法，非光滑控制具有3个优势[)，即快速 料的类型有关，因此可直接通过结构的材料阻尼比 的收敛性能、良好的抗干扰能力和广泛的应用范围. 来确定.阻尼比取值一般为0.05.若结构的材料相
刘彦辉等: 基于滚轴支座基础智能隔震结构的非光滑主动控制 feedback control for base鄄isolated structures by using fewer feedbacks. Furthermore, the nonsmooth control algorithm has great stability. KEY WORDS smart isolation; nonsmooth control; linear quadratic Gaussian (LQG) control; stability; seismic response 随着建筑结构往高层化和复杂化的发展,在地 震波激励下的建筑结构振动响应和安全话题已经引 起人们的关注. 传统的抗震结构通过增加抗侧刚度 来减小地震的影响,但结构刚度越大,地震响应越 大,这类结构体系非但不经济,且安全性无法保证. 对结构进行振动控制一直是土木工程领域重要的课 题. 为了解决这样的工程抗震问题,人们开始采用 结构振动控制方式来抑制结构振动的响应,从而提 高和改善结构的抗震性能. 有学者提出了结构振动 控制概念[1] ,结构振动控制作为有效抑制动力荷载 作用下结构振动响应的方法. 其中,基础隔震技术 在实际工程中的应用较为广泛,许多国家已将基础 隔震技术纳入抗震规范. 基础隔震技术能有效减少 上部结构地震响应,然而会使隔震层产生较大的位 移,这限制了隔震技术的广泛应用. 因此采用主动 控制、半主动控制方式来完善隔震结构的整体性能. 对于智能隔震控制系统而言,核心是控制算法,将控 制算法载入于计算机中,根据设置在结构上的传感 器采集的状态信号,反馈于已经载入计算机中所设 计的控制算法,从而获得控制力,然后通过伺服系统 将控制力转为驱动力施加给受控结构,以削弱结构 的振动反应. 具有代表性的控制算法有经典最优控 制( linear quadratic regulator, LQR) [2] ,线性二次型 高斯控制( linar quadratic Guassican, LQG) [2] ,极点 配置算法[3] ,瞬时最优控制( instantaneous optimal control, IOC) [4]等. 然而,传统控制算法对参数摄动 和外部干扰过于灵敏,发展和应用对系统参数影响 小且外干扰变化灵敏度更小、鲁棒性更强的结构振 动控制算法就显得极为重要. 滑膜控制[5鄄鄄6] 是一种 非连续性的反馈控制,可以使得闭环系统在较短时 间内收敛,并具有较好鲁棒性的特性,但相对于系统 状态变量而言是非连续的,控制过程中的高速切换 容易产生抖振现象. 设计一种概念清晰、简单有效的反馈控制策略, 使得闭环系统在有限时间内收敛且具备强鲁棒性的 特性,又不会产生由于非连续性控制引起的抖振现 象,将在工程应用中具有较好的发展前景,而非光滑 控制正是这样一种算法. 随着齐次性方法和 Lya鄄 punov 理论的研究和发展[7鄄鄄9] ,介于连续和非连续之 间的非光滑控制有了一定的研究和应用. 相比于传 统控制算法,非光滑控制具有 3 个优势[10] ,即快速 的收敛性能、良好的抗干扰能力和广泛的应用范围. 目前已经应用于航天姿态控制[11] 、高精度制导率设 计[12] 、永磁同步机位置控制[13] 和轮式机器人[14] 等 领域. 本文首先针对基础智能隔震结构,基于隔震层 位移和速度反馈,推导了智能隔震结构的非光滑控 制算法,然后对该智能控制系统在地震波激励下进 行仿真分析,对比分析了基础隔震、LQG 控制和非 光滑控制下的仿真结果,验证了非光滑控制算法的 有效性和稳定性. 1 受控结构方程 1郾 1 受控结构运动方程 采用层剪切模型,n 层的建筑结构简化为一维 具有 n 个自由度的结构模型,见图 1,图中 kb ,k1 ,k2 , …,kn 为层刚度;mb ,m1 ,m2 ,…,mn 为层质量. 在单 向水平地震波作用下,受控结构的运动方程为: M D ·· + C D · + Fs(D,D · ) = - MIx ·· g + EU (1) 式中, D = [db d1 d2 … dn ] T 为结构位移向 量,其中 db 为隔震层位移,dn 为结构第 n 层相对地 面的位移;D · 、D ·· 分别为速度矢量和加速度矢量; x ·· g 为地面地震加速度;M 和 C 分别为结构质量和阻尼 矩阵;Fs(D,D · ) = KD + L 为结构恢复力向量,K 为 刚 度 矩 阵, 其 中 L 为 塑 性 恢 复 力 向 量; I = [ 1 1 … 1 ] T 1 伊 n ; E 为 位 置 矩 阵; U = [u1 u2 … un ] T 为控制力输入,un 为第 n 个控 制力. 1郾 2 阻尼的选择 对于一维多自由度体系,上部结构为弹性,因此 选用瑞雷阻尼(Rayleigh damping)来组合阻尼矩阵. C = cM + qK (2) 式中,c、q 为比例系数,分别表示为: c = 2wiwj w 2 j - w 2 i (wj 孜i - wi 孜j) q = 2wiwj w 2 j - w 2 ( i 1 wi 孜j - 1 wj 孜 ) ì î í ï ï ï ï i (3) 其中,孜i 和 孜j 分别为结构第 i 和 j 阶振型阻尼比,wi、 wj 分别表示结构第 i、j 阶频率. 孜i、孜j 随着结构振幅 和材料的不同而发生变化,但受到频率的影响较小. 在进行结构反应分析时,各阶的阻尼比只与结构材 料的类型有关,因此可直接通过结构的材料阻尼比 来确定. 阻尼比取值一般为 0郾 05. 若结构的材料相 ·1093·
·1094· 工程科学学报，第41卷，第8期 式进行控制6).相比于传统的稳定闭环控制系统， 非光滑策略下的闭环控制系统能使受控对象具有较 强的抗干扰性能].此外，非光滑控制策略能实现 闭环控制系统的有限时间稳定. 2.1非光滑控制理论 针对一类二阶系统，可描述如下： (s1=s2 (6) 52=p 其中：S,和s2为系统状态：y为正奇数，当y=1时， 上述的系统为双积分系统.本文研究的受控结构系 统(5)可转化为如下表达形式： (X1=X2 (7) ☑m X2=M-(-MI +EU-CX2-F,) W 令F=M-1(-MⅢ元.+EU-CX2-F,),则系统 (7)正是y取值为1的系统，即双积分系统 田 针对受控系统，采用基于齐次性方法的非光滑 图1一维简化模型 控制方式，利用齐次性理论，对式(7)建立全局有限 Fig.1 Simplified one-dimensional model 时间闭环控制系统.该控制器为： 同，可取：=专=专，式(3)可简化为： =-k sign(s)Is,Ia-ksign(s2)Is2 I (8) 20" 其中，>0，石>0,0<a<,4,=y+1) 0:+0 y 1+1 (4) 2.2非光滑控制算法设计 9s、 10:+10 根据有限时间稳定理论18]，考虑到本建筑结构 该阻尼矩阵能够应用于大部分结构的动力分 的运动方程，将地震波加速度x。视为外部干扰项. 析.本文采用瑞雷阻尼组合运动方程的阻尼矩阵以 在结构的隔震层设置作动器，引入变量，则V= 进行后面的仿真分析. [心d…n]=[i1…in]T,选择如下 1.3受控结构状态方程 控制力： 确定结构的阻尼矩阵后，针对受控结构的运动 EU=MV+CX+F.+MI花. (9) 方程式(1)，定义状态变量X=[XX2],其中， 将式(8)代入式(5)，形成如下系统方程： X=[xbx]=D,X2=[x2 求=AX+BV (10) x2]T=D,则式(1)的状态空间方程为[5] 0 其中，A= TI 07.=1 因此，通过设计控制律 (X=AX+B(EU+P) 00 (5) Y=C,X+D,(EU+P) V实现结构振动的控制.系统(10)可解耦为如下n 01 0 个子系统： 式中，B=M小A= -M-K -M-C ,C= (x1h=x211=X21 (in=x2n (11) 01 「01 =x21=1 (i2n=n 0 I ,D= 0P=-MI -L 根据式(6)~(8)，对上述子系统采用如下非光 -M-K -M-C LM- 滑控制律： 2基于土木工程结构的非光滑控制算法设 i=-k1·sigm(x1b)lx1h1-k2·sign(x2h)1x2h1 (12) 计和系统稳定性分析 2.3非光滑控制系统稳定性分析 非光滑控制是对受控系统进行控制的一种非常 针对式(l0)和式(11)，选择子系统的Lyapunov 有效的方式，大部分非线性体系均可通过非光滑方 函数为[1]：
工程科学学报,第 41 卷,第 8 期 图 1 一维简化模型 Fig. 1 Simplified one鄄dimensional model 同,可取 孜i = 孜j = 孜,式(3)可简化为: c = 2wiwj wi + wj 孜 q = 2孜 wi + w ì î í ï ï ï ï j (4) 该阻尼矩阵能够应用于大部分结构的动力分 析. 本文采用瑞雷阻尼组合运动方程的阻尼矩阵以 进行后面的仿真分析. 1郾 3 受控结构状态方程 确定结构的阻尼矩阵后,针对受控结构的运动 方程式(1),定义状态变量 X = [X1 X2 ] T ,其中, X1 = [x1b x11 … x1n ] T = D,X2 = [ x2b x21 … x2n ] T = D · ,则式(1)的状态空间方程为[15] : X · = AX + B(EU + P) Y = CvX + Dv(EU + P { ) (5) 式中,B = 0 M é ë ê ê ù û ú - 1 ú ,A = 0 I - M - 1K - M - 1 é ë ê ê ù û ú ú C ,Cv = I 0 0 - M -1K I - M -1 é ë ê ê ê ù û ú ú ú C ,Dv = 0 0 M é ë ê ê ê ù û ú ú -1 ú ,P = - MIx ·· g - L. 2 基于土木工程结构的非光滑控制算法设 计和系统稳定性分析 非光滑控制是对受控系统进行控制的一种非常 有效的方式,大部分非线性体系均可通过非光滑方 式进行控制[16] . 相比于传统的稳定闭环控制系统, 非光滑策略下的闭环控制系统能使受控对象具有较 强的抗干扰性能[17] . 此外,非光滑控制策略能实现 闭环控制系统的有限时间稳定[11] . 2郾 1 非光滑控制理论 针对一类二阶系统,可描述如下: s · 1 = s 酌 2 s · 2 = { v (6) 其中:s1和 s2为系统状态;酌 为正奇数,当 酌 = 1 时, 上述的系统为双积分系统. 本文研究的受控结构系 统(5)可转化为如下表达形式: X · 1 = X2 X · 2 = M - 1 ( - MIx ·· g + EU - CX2 - Fs { ) (7) 令 v = M - 1 ( - MIx ·· g + EU - CX2 - Fs),则系统 (7)正是 酌 取值为 1 的系统,即双积分系统. 针对受控系统,采用基于齐次性方法的非光滑 控制方式,利用齐次性理论,对式(7)建立全局有限 时间闭环控制系统. 该控制器为: v = - k1 sign(s1 ) |s1 | 琢1 - k2 sign(s2 ) |s2 | 琢2 (8) 其中,k1 > 0,k2 > 0,0 < 琢1 < 1 酌 ,琢2 = (酌 + 1)琢1 1 + 琢1 . 2郾 2 非光滑控制算法设计 根据有限时间稳定理论[18] ,考虑到本建筑结构 的运动方程,将地震波加速度 x ·· g 视为外部干扰项. 在结构的隔震层设置作动器,引入变量 vb ,则 V = [vb d ·· 1 … d ·· n ] T = [ vb v1 … vn ] T ,选择如下 控制力: EU = M V + CX · + Fs + MIx ·· g (9) 将式(8)代入式(5),形成如下系统方程: X · = AX + B V (10) 其中,A = éI 0 ë ê ê ù û ú ú 0 0 ,B = é0 ë ê ê ù û ú ú I . 因此,通过设计控制律 V 实现结构振动的控制. 系统(10)可解耦为如下 n 个子系统: x · 1b = x2b x · 2b = v { b , x · 11 = x21 x · 21 = v { 1 ,…, x · 1n = x2n x · 2n = v { n (11) 根据式(6) ~ (8),对上述子系统采用如下非光 滑控制律: vb = - k1·sign(x1b ) | x1b | a1 - k2·sign(x2b ) | x2b | a2 (12) 2郾 3 非光滑控制系统稳定性分析 针对式(10)和式(11),选择子系统的 Lyapunov 函数为[19] : ·1094·
刘彦辉等：基于滚轴支座基础智能隔震结构的非光滑主动控制 ·1095· Fa)=7ln+宁 (13) L-sign()-sign() 其中，a,为非光滑控制参数.显然，F(xh,x)是正 [E+f(xbx2） 定的，当x1b≠0时，F(x,x2)是连续可微的，且 (16) E2+2(x1b,x2h） d1x+=(1+a)lxb1西sgn(xb.则式(13) dx 即 求导如下 p1+7=P2 (17) F()=ksign()+ p2+7=P1a1=P2a2 石2(-k1sign(x1b)lx1b1-k2sig(x2)lx21）= 其中e为大于0的常数取瓜=1=生时则 -k2x2sign(x2h)1xh1=-k21x211+(14) 式(17)即所建立的子系统(11)为负齐次系统，且有 显然，F(x,x)是半负定的，因此该子系统是 稳定的.定义如下集合： 负齐次度刀：<0，因此，所建立的子系统 H={(x,x2):F(x1b,x2)} (15) (11)为全局有限时间稳定，同理，系统(13)也为全 令x(t)=0,则可得x(t)=b(t)=0,代入 局有限时间稳定 闭环系统，可得xb(t)=0,因此集合H所包含的最 3滚轴支座智能隔震结构 大不变集为(xhx)=(00).由不变性原理可 得，系统(10)在平衡点(xx)=(00)为渐近 本文选取滚轴支座基础隔震结构作为智能隔震 稳定，也保证了整个智能隔震系统的全局渐近稳定 结构中的受控结构，在隔震结构的隔震层设置控制 根据齐次系统定义，令子系统f(xx如)= 器，形成智能隔震体系 3.1滚轴支座 [听]'=[x元]T,则有如下方程： 滚轴支座的单元见图2(a),由圆柱体滚轴、带 f(ex1b,ex2）= 有倾斜面的上滚座和下滚座三个部分组成.当地震 作用与圆柱体滚轴滚动方向一致时，通过上下滚座 -k1 sign(x1b)lx1b1a-石2sign(x2b)lx2b1西 错开运动和滚轴来回滚动来减小地震作用对上部结 0 构的影响.图2(b)为两个个滚轴支座单元组合后 E+1 的示意图，该支座能应对双向的地震作用. 上滚座 滚轴 下滚座 (a 问 图2滚轴支座单元和滚动示意图.()滚轴支座单元：(b)滚轴支座滚动 Fig.2 Unit and roll of roller bearing:(a)unit of roller bearing:(b)roll of roller bearing 图3中，k为滚轴支座的有效刚度，k,为滚轴支 式(18)中，x为滚轴支座位移，通过式(18)，可得到 座刚度，W为上层结构总重力，d、S分别为滚轴支 地震作用下结构隔震层恢复力 座的设计最大位移和接触弧的弧长，R为滚轴的半 3.2结构模型的选取 径，σ为滚轴支座的倾角.滚轴支座力学模型表达 以某实际工程为例，该结构总高为25.95m,上 式为： 部结构为6层，1至5层高为4.2m,第六层为3m, W F=2R'sgn () 隔震层高度为1.8m,结构形式为框架结构.结构设 lxbI≥S (18) 计使用年限：50a.建筑设防分类：重点设防类（乙 W-tgo·sg(x),lxI<S 类).建筑物抗震设防：8度，基本地震加速度为
刘彦辉等: 基于滚轴支座基础智能隔震结构的非光滑主动控制 F(x1b ,x2b ) = k1 1 + 琢1 | x1b | 1 + a1 + 1 2 x 2 2b (13) 其中,a1 为非光滑控制参数. 显然,F( x1b ,x2b )是正 定的,当 x1b 屹0 时,F ( x1b , x2b ) 是连续可微的,且 d dx1b | x1b | 1 + a1 = (1 + a1 ) | x1b | a1 sgn ( x1b ). 则式(13) 求导如下 F · (x1b ,x2b ) = k1 sign(x1b ) | x1b | a1 x2b + k2 ( - k1 sign(x1b ) | x1b | a1 - k2 sign(x2b ) | x2b | a2 ) = - k2 x2b sign(x2b ) | x2b | a2 = - k2 | x2b | 1 + a2 (14) 显然,F · (x1b ,x2b )是半负定的,因此该子系统是 稳定的. 定义如下集合: H = {(x1b ,x2b ):F(x1b ,x2b )} (15) 令 x2b (t) = 0,则可得 x2b ( t) = x · 1b ( t) = 0,代入 闭环系统,可得 x1b (t) = 0,因此集合 H 所包含的最 大不变集为( x1b x2b ) = (0 0). 由不变性原理可 得,系统(10)在平衡点( x1b x2b ) = (0 0)为渐近 稳定,也保证了整个智能隔震系统的全局渐近稳定. 根据 齐 次 系 统 定 义, 令 子 系 统 f ( x1b x2b ) = [f 1 f 2 ] T = [ x · 1b x · 2b ] T ,则有如下方程: f(着 籽1 x1b ,着 籽2 x2b ) = 着 籽2 x2b - k1着 籽1 a1 sign(x1b ) | x1b | a1 - k2 sign(x2b ) | x2b | a é ë ê ê ù û ú 2 ú = 着 籽1 + 1 0 0 着 籽2 é ë ê ê ù û ú + 1 ú· x2b - k1 sign(x1b ) | x1b | a1 - k2 sign(x2b ) | x2b | a é ë ê ê ù û ú 2 ú = 着 籽1 + 1 f 1 (x1b ,x2b ) 着 籽2 + 1 f 2 (x1b ,x2b é ë ê ê ù û ú ) ú (16) 即 籽1 + 浊 = 籽2 籽2 + 浊 = 籽1·a1 = 籽2·a { 2 (17) 其中,着 为大于 0 的常数,取 籽1 = 1,籽2 = 1 + a1 2 时,则 式(17)即所建立的子系统(11)为负齐次系统,且有 负齐次度 浊 = a2 - 1 2 < 0. 因此,所建立的子系统 (11)为全局有限时间稳定,同理,系统(13)也为全 局有限时间稳定. 3 滚轴支座智能隔震结构 本文选取滚轴支座基础隔震结构作为智能隔震 结构中的受控结构,在隔震结构的隔震层设置控制 器,形成智能隔震体系. 3郾 1 滚轴支座 滚轴支座的单元见图 2( a),由圆柱体滚轴、带 有倾斜面的上滚座和下滚座三个部分组成. 当地震 作用与圆柱体滚轴滚动方向一致时,通过上下滚座 错开运动和滚轴来回滚动来减小地震作用对上部结 构的影响. 图 2( b)为两个个滚轴支座单元组合后 的示意图,该支座能应对双向的地震作用. 图 2 滚轴支座单元和滚动示意图 郾 (a) 滚轴支座单元; (b) 滚轴支座滚动 Fig. 2 Unit and roll of roller bearing: (a) unit of roller bearing; (b) roll of roller bearing 图 3 中,ke为滚轴支座的有效刚度,kb为滚轴支 座刚度,W 为上层结构总重力,dmax、S 分别为滚轴支 座的设计最大位移和接触弧的弧长,R 为滚轴的半 径,滓 为滚轴支座的倾角. 滚轴支座力学模型表达 式为: Fb = W 2R xb·sgn (xb ), | xb |逸S W·tg滓·sgn (xb ), | xb | < { S (18) 式(18)中,xb为滚轴支座位移,通过式(18),可得到 地震作用下结构隔震层恢复力. 3郾 2 结构模型的选取 以某实际工程为例,该结构总高为 25郾 95 m,上 部结构为 6 层,1 至 5 层高为 4郾 2 m,第六层为 3 m, 隔震层高度为 1郾 8 m,结构形式为框架结构. 结构设 计使用年限:50 a. 建筑设防分类:重点设防类(乙 类). 建筑物抗震设防:8 度,基本地震加速度为 ·1095·
·1096· 工程科学学报，第41卷，第8期 力 W.tgo W2 位移 -W.Igo 图3滚轴支座力学模型 Fig.3 Roll diagram of roller bearing 0.3g.地震分组：第三组.建筑场地类别：Ⅱ类.建 作动器 筑物抗震等级：框架抗震等级二级.利用ETABS建 立结构有限元模型，见图4.该结构隔震层到第六层 的层质量分别取为809.2,705.0,705.0,696.9， 图6智能隔震结构简化模型 686.4,648.8和35.8t:刚度取值分别为3.14×10， Fig.6 Simplified model of smart-isolated structure 2.51×103,2.35×10,1.82×10°，1.59×10,3.71×10 表1 MATLAB模型与ETABS模型周期 kNm-1.在该结构隔震层设置滚轴支座，滚轴支座 Table 1 Period of MATLAB model and ETABS model 布置位置见图5，一共布置18个支座 振型号 MTALB周期/s ETABS周期/s 1.160 1.209 2 0.436 0.422 3 0.275 0.235 4 0.218 0.182 0.187 0.142 6 0.180 0.116 4 主动智能隔震结构的仿真分析 图4实际结构有限元模型 Fig.4 Finite element model of the actual structure 4.1地震波的选择 46.4m 本文选用三条地震波作为仿真控制系统的地震 波输入，分别为LK0676波，LK0714波和人工波 0 ◇ LK0676波时间间隔为0.005s,持续时间为48s. LK0714波时间间隔为0.003s,持续时间为35s.人 工波时间间隔为0.01s,持续时间为30s.三条地震 图5滚轴支座布置图 波时程曲线如图7~图9.三条地震波加速度反应 Fig.5 Diagram of roller bearing arrangement 谱与标准反应谱的对比见图10. 从图10中地震反应和规范反应谱对比可以看 将滚轴支座基础隔震结构简化为一维多自由度 出，三条地震波下基础隔震结构（本文选用的隔震 模型，见图6，在结构隔震层设置理想控制器，形成 结构模型)周期对应的加速度反应谱平均值与标准 智能隔震结构 反应谱相差不超过20%，满足《建筑抗震设计规范》 分别用MATLAB和ETABS建立非隔震结构模 要求. 型，计算得出的周期结果对比见表1，结果表明周期 4.2隔震层的响应分析 误差相差较小，可进一步深入分析隔震结构模型. 在结构基础和上部结构连接处设置隔震器，以
工程科学学报,第 41 卷,第 8 期 图 3 滚轴支座力学模型 Fig. 3 Roll diagram of roller bearing 0郾 3g. 地震分组:第三组. 建筑场地类别:II 类. 建 筑物抗震等级:框架抗震等级二级. 利用 ETABS 建 立结构有限元模型,见图 4. 该结构隔震层到第六层 的层 质 量 分 别 取 为 809郾 2, 705郾 0, 705郾 0, 696郾 9, 686郾 4,648郾 8 和 35郾 8 t;刚度取值分别为 3郾 14 伊 10 8 , 2郾 51 伊10 8 ,2郾 35 伊 10 8 ,1郾 82 伊 10 8 ,1郾 59 伊 10 8 ,3郾 71 伊 10 7 kN·m - 1 . 在该结构隔震层设置滚轴支座,滚轴支座 布置位置见图 5,一共布置 18 个支座. 图 4 实际结构有限元模型 Fig. 4 Finite element model of the actual structure 图 5 滚轴支座布置图 Fig. 5 Diagram of roller bearing arrangement 将滚轴支座基础隔震结构简化为一维多自由度 模型,见图 6,在结构隔震层设置理想控制器,形成 智能隔震结构. 分别用 MATLAB 和 ETABS 建立非隔震结构模 型,计算得出的周期结果对比见表 1,结果表明周期 误差相差较小,可进一步深入分析隔震结构模型. 图 6 智能隔震结构简化模型 Fig. 6 Simplified model of smart鄄isolated structure 表 1 MATLAB 模型与 ETABS 模型周期 Table 1 Period of MATLAB model and ETABS model 振型号 MTALB 周期/ s ETABS 周期/ s 1 1郾 160 1郾 209 2 0郾 436 0郾 422 3 0郾 275 0郾 235 4 0郾 218 0郾 182 5 0郾 187 0郾 142 6 0郾 180 0郾 116 4 主动智能隔震结构的仿真分析 4郾 1 地震波的选择 本文选用三条地震波作为仿真控制系统的地震 波输入,分别为 LK0676 波, LK0714 波和人工波. LK0676 波时间间隔为 0郾 005 s,持续时间为 48 s. LK0714 波时间间隔为 0郾 003 s,持续时间为 35 s. 人 工波时间间隔为 0郾 01 s,持续时间为 30 s. 三条地震 波时程曲线如图 7 ~ 图 9. 三条地震波加速度反应 谱与标准反应谱的对比见图 10. 从图 10 中地震反应和规范反应谱对比可以看 出,三条地震波下基础隔震结构(本文选用的隔震 结构模型)周期对应的加速度反应谱平均值与标准 反应谱相差不超过 20% ,满足《建筑抗震设计规范》 要求. 4郾 2 隔震层的响应分析 在结构基础和上部结构连接处设置隔震器,以 ·1096·
刘彦辉等：基于滚轴支座基础智能隔震结构的非光滑主动控制 ·1097· 10 20 20 30 40 50 时间/s 时间/s 图9人工波时程曲线 图7LK0676波时程曲线 Fig.9 Time history curve of artificial wave Fig.7 Time history curve of LK0676 wave 1.2 标准反应谱 1.0 LK0676 LK0714 08 人工波 平均值 06 0.4 0.2 非隔震 隔震 2 10 15.20253035 周期/s 时间s 图10反应谱对比 图8LK0714波时程曲线 Fig.8 Time history curve of LK0714 wave Fig.10 Comparison of response spectrum 阻止地震能量输入于上部结构中，从而确保上部结 将三条地震波依次施加于受控结构，经仿真分 构的安全.故在地震波激励下，隔震层的反应特征 析，得出基础隔震、非光滑主动控制和LQG主动控 能作为分析智能隔震或者基础隔震的结构振动抑制 制下的位移响应和加速度响应，见图11~13 效果的重要指标.由于隔震层设置了控制器，这里 从以上隔震层响应图中，在非光滑控制和LQG 对基于非光滑控制的智能隔震体系、基于LQG控制 控制下，隔震层的位移响应和加速度响应都比隔震 的智能隔震体系和基础隔震体系的隔震层响应进行 的有所下降.下面给出隔震、LQG控制和非光滑控 对比分析. 制下的隔震层响应峰值，见表2 120 …基础隔震 2.5m b) ,基础隔震 90 非光滑主动 非光滑主动 --LQG主动 -LQG主动 30 -30 60 10 20 256 10 20 30 40 50 时间s 时间s 图11LKO676波作用下结构隔震层响应曲线对比.(a)LKO676波作用下隔震层位移曲线对比：(b)LK0676波作用下隔震层加速度曲线 对比 Fig.11 Curve comparison of structural isolation layer response under LK0676 wave:(a)curve comparison of isolation layer displacement under LK0676 wave:(b)curve comparison of isolation layer acceleration under LK0676 wave
刘彦辉等: 基于滚轴支座基础智能隔震结构的非光滑主动控制 图 11 LK0676 波作用下结构隔震层响应曲线对比 郾 (a) LK0676 波作用下隔震层位移曲线对比; (b) LK0676 波作用下隔震层加速度曲线 对比 Fig. 11 Curve comparison of structural isolation layer response under LK0676 wave: ( a) curve comparison of isolation layer displacement under LK0676 wave; (b) curve comparison of isolation layer acceleration under LK0676 wave 图 7 LK0676 波时程曲线 Fig. 7 Time history curve of LK0676 wave 图 8 LK0714 波时程曲线 Fig. 8 Time history curve of LK0714 wave 阻止地震能量输入于上部结构中,从而确保上部结 构的安全. 故在地震波激励下,隔震层的反应特征 能作为分析智能隔震或者基础隔震的结构振动抑制 效果的重要指标. 由于隔震层设置了控制器,这里 对基于非光滑控制的智能隔震体系、基于 LQG 控制 的智能隔震体系和基础隔震体系的隔震层响应进行 对比分析. 图 9 人工波时程曲线 Fig. 9 Time history curve of artificial wave 图 10 反应谱对比 Fig. 10 Comparison of response spectrum 将三条地震波依次施加于受控结构,经仿真分 析,得出基础隔震、非光滑主动控制和 LQG 主动控 制下的位移响应和加速度响应,见图 11 ~ 13. 从以上隔震层响应图中,在非光滑控制和 LQG 控制下,隔震层的位移响应和加速度响应都比隔震 的有所下降. 下面给出隔震、LQG 控制和非光滑控 制下的隔震层响应峰值,见表 2. ·1097·
·1098· 工程科学学报，第41卷，第8期 120r 2.5 a 60 -60 ,基础隔震 一非光滑主动 基础隔震 非光滑主动 ---1QG主动 ---LQG主动 -1200 5 10 1520 25 30 35 -2.50 5 10 1520 253035 时间/s 时间 图12LK0714波作用下结构隔震层响应曲线对比.(a)LK0714波作用下隔震层位移曲线对比：(b)LK0714波作用下隔震层加速度曲线 对比 Fig.12 Curve comparison of structural isolation layer response under LK714 wave:(a)curve comparison of isolation laver displacement under LK714 wave:(b)curve comparison of isolation layer acceleration under LK714 wave I00 b) :基础隔震 …基础隔震 非光滑主动 一非光滑主动 50 ---L0G主动 ·=-LQG主动 50 -1006 10 20 30 10 20 30 时间s 时间s 图13人工波作用下结构隔震层响应曲线对比.(a)人工波作用下隔震层位移曲线对比：(b)人工波作用下隔震层加速度曲线 Fig.13 Curve comparison of structural isolation layer response under artificial wave:(a)curve comparison of isolation layer displacement under arti- ficial wave;(b)curve comparison of isolation layer acceleration under artificial wave 表2三条波分别作用下的隔震层响应峰值 Table 2 Peak values of isolation layer responses under three different waves LK0676波 LK0714 人工波 控制策略 位移/ 速度/ 加速度/ 位移/ 速度/ 加速度/ 位移/ 速度/ 加速度/ mm (ms1) (ms2) mm (m.s-1) (m52) mm (m.s-1) (ms2) 隔震 120.65 0.49 2.24 110.58 0.54 2.12 91.92 0.49 1.87 LQG 43.57 0.37 2.06 48.92 0.36 1.93 36.48 0.24 1.63 非光滑 40.10 0.32 1.48 39.56 0.27 1.59 29.44 0.21 1.25 从表2可以看出，三条地震波作用下，相比于隔 4.3上部结构响应峰值分析 震，LQG主动智能隔震结构隔震层位移分别下降了 图14~图19为上部结构各层响应峰值，从图 64%、56%和60%，非光滑主动智能隔震结构隔震 14~图19可以看出，相比于基础隔震，非光滑主动 层位移则分别下降了67%、64%和68%，其对位移 控制和LQG主动控制的层间位移峰值有大幅度的 峰值的控制效果优于LOG主动控制.两种控制算 下降，而且非光滑主动控制下层间位移峰值均比 法下受控结构的速度和加速度峰值均有所下降，尤 LQG主动控制的低，同时从结构加速度可以看出上 其是加速度峰值的控制，非光滑控制下的下降幅度 部结构基本为平动，但图19人工波作用下加速度峰 明显大于LQG控制.因此，对于隔震层响应峰值的 第3层加速度比其他层加速度略微偏小，分析主要 控制，非光滑主动控制效果明显优于LQG主动 是控制器控制力对上部结构动力响应相位有一定影 控制. 响，高阶振型有一定的参与.同样，两种主动控制方 式下结构的加速度峰值也有着明显的下降，且非光
工程科学学报,第 41 卷,第 8 期 图 12 LK0714 波作用下结构隔震层响应曲线对比 郾 (a) LK0714 波作用下隔震层位移曲线对比; (b) LK0714 波作用下隔震层加速度曲线 对比 Fig. 12 Curve comparison of structural isolation layer response under LK714 wave: ( a) curve comparison of isolation layer displacement under LK714 wave; (b) curve comparison of isolation layer acceleration under LK714 wave 图 13 人工波作用下结构隔震层响应曲线对比 郾 (a) 人工波作用下隔震层位移曲线对比; (b) 人工波作用下隔震层加速度曲线 Fig. 13 Curve comparison of structural isolation layer response under artificial wave: (a) curve comparison of isolation layer displacement under arti鄄 ficial wave; (b) curve comparison of isolation layer acceleration under artificial wave 表 2 三条波分别作用下的隔震层响应峰值 Table 2 Peak values of isolation layer responses under three different waves 控制策略 LK0676 波 LK0714 人工波 位移/ mm 速度/ (m·s - 1 ) 加速度/ (m·s - 2 ) 位移/ mm 速度/ (m·s - 1 ) 加速度/ (m·s - 2 ) 位移/ mm 速度/ (m·s - 1 ) 加速度/ (m·s - 2 ) 隔震 120郾 65 0郾 49 2郾 24 110郾 58 0郾 54 2郾 12 91郾 92 0郾 49 1郾 87 LQG 43郾 57 0郾 37 2郾 06 48郾 92 0郾 36 1郾 93 36郾 48 0郾 24 1郾 63 非光滑 40郾 10 0郾 32 1郾 48 39郾 56 0郾 27 1郾 59 29郾 44 0郾 21 1郾 25 从表 2 可以看出,三条地震波作用下,相比于隔 震,LQG 主动智能隔震结构隔震层位移分别下降了 64% 、56% 和 60% ,非光滑主动智能隔震结构隔震 层位移则分别下降了 67% 、64% 和 68% ,其对位移 峰值的控制效果优于 LQG 主动控制. 两种控制算 法下受控结构的速度和加速度峰值均有所下降,尤 其是加速度峰值的控制,非光滑控制下的下降幅度 明显大于 LQG 控制. 因此,对于隔震层响应峰值的 控制, 非光滑主动控制效果明显优于 LQG 主动 控制. 4郾 3 上部结构响应峰值分析 图 14 ~ 图 19 为上部结构各层响应峰值,从图 14 ~ 图 19 可以看出,相比于基础隔震,非光滑主动 控制和 LQG 主动控制的层间位移峰值有大幅度的 下降,而且非光滑主动控制下层间位移峰值均比 LQG 主动控制的低,同时从结构加速度可以看出上 部结构基本为平动,但图 19 人工波作用下加速度峰 第 3 层加速度比其他层加速度略微偏小,分析主要 是控制器控制力对上部结构动力响应相位有一定影 响,高阶振型有一定的参与. 同样,两种主动控制方 式下结构的加速度峰值也有着明显的下降,且非光 ·1098·
刘彦辉等：基于滚轴支座基础智能隔震结构的非光滑主动控制 ·1099· 滑主动控制下受控结构的加速度蜂值明显小于 LQG主动控制. 米 基础隔震 5 日一 非光滑主动 一[QG主动 3 3 米一基础隔震 一日一非光滑主动 0.5 1.0 1.5 20 25 △一LQG主动 加速度/m·8 图17LK0714波作用下加速度蜂值 0.005 0.0100.015 0.020 0.025 Fig.17 Peak acceleration under LK0714 层间位移lmm 图14LK0676波作用下上部结构层间位移峰值 Fig.14 Inter-story displacement peak of superstructure under 米一基础隔震 LK0676 日一非光滑主动 》 合一LQG主动 米一隔震 一日一非光滑主动 合一1QG主动 23 必 0.005 0.0100.015 0.020 0.025 层间位移mm 图18人工波作用下上部结构层间位移峰值 Fig.18 Inter-story displacement peak of superstructure under artifi- 1.0 15 2.0 2.5 cial wave 加速度/ms) 图15LK0676波作用下加速度峰值 Fig.15 Peak acceleration under LK0676 米一基础隔震 一日一非光滑主动 △一L0G主动 3 米 米一基础隔震 日一非光滑主动 0.5 1.0 1.5 2.0 △一LQG主动 加速度/(m·s 图19人工波作用下加速度峰值 0 0.005 0.0100.015 0.020 0.025 Fig.19 Peak acceleration under artificial wave 层间位移mm 图16LK0714波作用下上部结构层间位移峰值 看出，两种主动控制下的受控结构剪力小于基础 Fig.16 Inter-story displacement peak of superstructure under 隔震结构，且非光滑主动控制下结构隔震层和第 LK0714 一层的底部剪力均比LQG主动控制的有所下降. 4.4结构剪力分析 为了更加便于对比结构剪力峰值，将相关数据列 当智能隔震结构遭遇地震波激励时，主动控制 于表3，可以看出，两种主动控制算法下受控结构 智能隔震结构的隔震层底部和第一层底部所产生 剪力峰值小于基础隔震结构的剪力峰值，而非光 的剪力见图20(a)~(c).从图20(a)~(c)可以 滑控制算法对于结构剪力峰值的控制优于LQG控
刘彦辉等: 基于滚轴支座基础智能隔震结构的非光滑主动控制 滑主动控制下受控结构的加速度峰值明显小于 LQG 主动控制. 图 14 LK0676 波作用下上部结构层间位移峰值 Fig. 14 Inter鄄story displacement peak of superstructure under LK0676 图 15 LK0676 波作用下加速度峰值 Fig. 15 Peak acceleration under LK0676 图 16 LK0714 波作用下上部结构层间位移峰值 Fig. 16 Inter鄄story displacement peak of superstructure under LK0714 4郾 4 结构剪力分析 当智能隔震结构遭遇地震波激励时,主动控制 智能隔震结构的隔震层底部和第一层底部所产生 的剪力见图 20( a) ~ ( c) . 从图 20( a) ~ ( c) 可以 图 17 LK0714 波作用下加速度峰值 Fig. 17 Peak acceleration under LK0714 图 18 人工波作用下上部结构层间位移峰值 Fig. 18 Inter鄄story displacement peak of superstructure under artifi鄄 cial wave 图 19 人工波作用下加速度峰值 Fig. 19 Peak acceleration under artificial wave 看出,两种主动控制下的受控结构剪力小于基础 隔震结构,且非光滑主动控制下结构隔震层和第 一层的底部剪力均比 LQG 主动控制的有所下降. 为了更加便于对比结构剪力峰值,将相关数据列 于表 3,可以看出,两种主动控制算法下受控结构 剪力峰值小于基础隔震结构的剪力峰值,而非光 滑控制算法对于结构剪力峰值的控制优于 LQG 控 ·1099·
·1100· 工程科学学报，第41卷，第8期 制算法 从图24中可以看出，相比于基础隔震结构，主 表3隔震层底部剪力峰值 动控制下的受控结构能在很短时间收敛至平衡点， Table 3 Peak value of bottom shear force in the isolation layer 且非光滑控制下的受控结构响应收敛范围小于 剪力峰值/(103kN) 控制策略 LQG控制，说明非光滑控制算法的抗扰能力优于 LK0676波 LK0714波 人工波 LQG控制 隔震 9.602 9.098 8.029 LOG 8.910 8.317 7.107 5结论 非光滑 6.678 6.903 5.672 本文设计了基于滚轴支座基础智能隔震结构的 4.5主动控制力分析 非光滑控制算法，推导了非光滑控制算法下智能隔 前面对比了隔震、非光滑主动控制、LQG主动 震控制系统的稳定性，对比分析了基础隔震、LQG 控制下滚轴支座基础隔震结构的响应，可以看出智 控制和非光滑控制对结构响应的控制效果，得到以 能主动控制具有良好的控制效果，且非光滑主动控 下主要结论： 制在一定程度上比LQG主动控制对受控结构振动 (1)提出的针对轴支座基础智能隔震结构的非 的抑制效果要好.两种主动控制算法的主动控制力 光滑控制算法具有良好的减振性能，可有效控制结 曲线见图21~23 构地震响应，整体控制效果优于经典的QG控制 从图中可以看出，LK714波非光滑主动控制力 算法 在比LQG控制力输出稍大，但其他两条地震波下控 (2)提出的非光滑主动控制算法采用更少的反 制力峰值小于LQG主动控制 馈量—隔震层位移和速度，即可实现结构的振动 4.6控制器稳定性分析 控制，且控制效果良好，相比较经典LQG控制算法， 若结构的初始位移值为50mm,则基础隔震结 在实际工程中更容易实现. 构、LQG控制智能隔震结构、非光滑智能隔震结构 (3)提出的非光滑主动控制算法具有良好的稳 的复位情况见图24 定性，能以更快的速度使受控结构收敛于平衡点. 10 1.0r …基础隔震 (b) 一非光滑主动 05 ·=1.0G主动 15 15 …基础隔震 非光滑主动 -一-LQG主动 20 -100 5 10 15 20 时间/s 时间s .0 …基础隔震 非光滑主动 -LQG主动 0.5 10 20 时间s 图20三种波作用下结构隔震层底部剪力曲线.(a)LKO676波作用下结构隔震层底部剪力曲线：(b)LK0714波作用下结构隔震层底部 剪力曲线：（©）人工波作用下结构隔震层底部剪力曲线 Fig.20 Bottom shear force of structural isolation layer under three waves:(a)bottom shear force of structural isolation layer under;(b)bottom shear force of structural isolation layer under LK0714 wave;(c)bottom shear force of structural isolation layer under artificial wave
工程科学学报,第 41 卷,第 8 期 制算法. 表 3 隔震层底部剪力峰值 Table 3 Peak value of bottom shear force in the isolation layer 控制策略 剪力峰值/ (10 3 kN) LK0676 波 LK0714 波 人工波 隔震 9郾 602 9郾 098 8郾 029 LQG 8郾 910 8郾 317 7郾 107 非光滑 6郾 678 6郾 903 5郾 672 4郾 5 主动控制力分析 前面对比了隔震、非光滑主动控制、LQG 主动 图 20 三种波作用下结构隔震层底部剪力曲线 郾 (a) LK0676 波作用下结构隔震层底部剪力曲线; (b) LK0714 波作用下结构隔震层底部 剪力曲线; (c) 人工波作用下结构隔震层底部剪力曲线 Fig. 20 Bottom shear force of structural isolation layer under three waves: ( a) bottom shear force of structural isolation layer under; ( b) bottom shear force of structural isolation layer under LK0714 wave; (c) bottom shear force of structural isolation layer under artificial wave 控制下滚轴支座基础隔震结构的响应,可以看出智 能主动控制具有良好的控制效果,且非光滑主动控 制在一定程度上比 LQG 主动控制对受控结构振动 的抑制效果要好. 两种主动控制算法的主动控制力 曲线见图 21 ~ 23. 从图中可以看出,LK714 波非光滑主动控制力 在比 LQG 控制力输出稍大,但其他两条地震波下控 制力峰值小于 LQG 主动控制. 4郾 6 控制器稳定性分析 若结构的初始位移值为 50 mm,则基础隔震结 构、LQG 控制智能隔震结构、非光滑智能隔震结构 的复位情况见图 24. 从图 24 中可以看出,相比于基础隔震结构,主 动控制下的受控结构能在很短时间收敛至平衡点, 且非光滑控制下的受控结构响应收敛范围小于 LQG 控制,说明非光滑控制算法的抗扰能力优于 LQG 控制. 5 结论 本文设计了基于滚轴支座基础智能隔震结构的 非光滑控制算法,推导了非光滑控制算法下智能隔 震控制系统的稳定性,对比分析了基础隔震、LQG 控制和非光滑控制对结构响应的控制效果,得到以 下主要结论: (1)提出的针对轴支座基础智能隔震结构的非 光滑控制算法具有良好的减振性能,可有效控制结 构地震响应,整体控制效果优于经典的 LQG 控制 算法. (2)提出的非光滑主动控制算法采用更少的反 馈量———隔震层位移和速度,即可实现结构的振动 控制,且控制效果良好,相比较经典 LQG 控制算法, 在实际工程中更容易实现. (3)提出的非光滑主动控制算法具有良好的稳 定性,能以更快的速度使受控结构收敛于平衡点. ·1100·
刘彦辉等：基于滚轴支座基础智能隔震结构的非光滑主动控制 ·1101· 4000 参考文献 一非光滑主动 -LQG主动 [1]Zhou F L.Seismic Control in Engineering Structures.Beijing: 20D0 Seismological Press,1997 (周福霖.工程结构减震控制.北京：地震出版社，1997) [2]Ou J P.Structural Vibration Control-Actire,Semi-active and Intel- ligent control.Beijing:Science Press,2003 (欧进萍。结构振动控制：主动、半主动和智能控制.北京：科 -2000 学出版社，2003) [3]Wang J W,Yang Y F,Qian Y H.Design of controller for torsion 4000 20 30 40 50 vibration device based on pole assignment method.Exp Technol 时间/s Manage,2014,31(7):86 图21LK0676波用下控制力曲线 (王佳伟，杨亚非，钱玉恒.基于极点配置方法的扭转振动装 Fig.21 Curve of control force under LK0676 wave 置控制器设计.实验技术与管理，2014,31(7)：86) 4000 [4] Tong S W.Tang H P.Research on approximate discrete instanta- 一非光滑主动 =--1QG主动 neous optimal control algorithm and optimization of the actuator po- 2000 sitions.J Vib Eng,2016,29(5):920 (童少伟，唐怀平.近似离散瞬时最优控制及作动器位置优化 研究.振动工程学报，2016,29(5)：920) [5]Yang J N,Wu J C,Kawashima K,et al.Hybrid control of seis- mic-excited bridge structure.Earthquake Eng Struct Dyn,1995. -2000 24(11):1437 [6]Yang J N,Wu J C,Agrawal A K.Sliding mode control for seismi- 40006 35 cally excited linear structures.J Eng Mech,1995,121 (12): 10 15202530 时间/s 1386 图22LK0714波作用下控制力曲线 [7]Rosier L.Homogeneous Lyapunoy function for homogeneous con- Fig.22 Curve of control force under LK0714 wave tinuous vector field.Syst Control Lett,1992,19(6):467 [8]Kawski M.Stability and nilpotent approximations /Proceedings 一非光滑主动 of the 27th IEEE Conference on Decision and Control.Austin, --LQG主动 2000 1988:1244 [9]Hermes H.Asymptotically stabilizing feedback controls.I Differ- ential Equations 1,1991,92(1):76 [10]Wang S.Study on Non-smooth Control Algorithm for Smart Base Isolated Benchmark Building Dissertation].Guangrhou:Guan -1000 gzhou University,2016 (王思.基于Benchmark智能隔震结构的非光滑控制算法的 2500 研究与分析[学位论文].广州：广州大学，2016) 10 20 [11]Ma K M.Design of continuous non-smooth attitude control laws 3 时间s for spacecraft.J Astronautics,2012,33(6):713 图23人工波作用下控制力曲线 (马克茂.航天器连续非光滑姿态控制律设计.宇航学报， Fig.23 Curve of control force under artificial wave 2012.33(6):713) 50 [12]Ma K M.Non-smooth design and implementation of high-preci- …基础隔震 sion guidance laws.J Ball,2013,25(2):1 25月 一非光滑主动 (马克茂.高精度制导律的非光滑设计与实现.弹道学报。 ==LQG主动 2013,25(2):1) 03 [13]Wang J H,Wang Q,Ma K M.Non-smooth controller design for permanent magnet synchronous motors.Comput Simul,2016,33 25 (3):227 (王建晖，王清，马克茂.一种永磁同步电机非光滑控制器 506 设计.计算机仿真.2016,33(3)：227) 时间s [14]Li S H,Ding S H,Du H B,et al.Theory and Application of 图24基础隔震结构和智能隔震结构复位曲线 Non-smooth Control.Beijing:Science Press,2013 Fig.24 Reset curve of base-isolated structure and smart-isolated (李世华，丁世宏，都海波，等.非光滑控制理论与应用.北 structure 京：科学出版社，2013)
刘彦辉等: 基于滚轴支座基础智能隔震结构的非光滑主动控制 图 21 LK0676 波用下控制力曲线 Fig. 21 Curve of control force under LK0676 wave 图 22 LK0714 波作用下控制力曲线 Fig. 22 Curve of control force under LK0714 wave 图 23 人工波作用下控制力曲线 Fig. 23 Curve of control force under artificial wave 图 24 基础隔震结构和智能隔震结构复位曲线 Fig. 24 Reset curve of base鄄isolated structure and smart鄄isolated structure 参 考 文 献 [1] Zhou F L. Seismic Control in Engineering Structures. Beijing: Seismological Press, 1997 (周福霖. 工程结构减震控制. 北京: 地震出版社, 1997) [2] Ou J P. Structural Vibration Control鄄Active, Semi鄄active and Intel鄄 ligent control. Beijing: Science Press, 2003 (欧进萍. 结构振动控制: 主动、半主动和智能控制. 北京: 科 学出版社, 2003) [3] Wang J W, Yang Y F, Qian Y H. Design of controller for torsion vibration device based on pole assignment method. Exp Technol Manage, 2014, 31(7): 86 (王佳伟, 杨亚非, 钱玉恒. 基于极点配置方法的扭转振动装 置控制器设计. 实验技术与管理, 2014, 31(7): 86) [4] Tong S W, Tang H P. Research on approximate discrete instanta鄄 neous optimal control algorithm and optimization of the actuator po鄄 sitions. J Vib Eng, 2016, 29(5): 920 (童少伟, 唐怀平. 近似离散瞬时最优控制及作动器位置优化 研究. 振动工程学报, 2016, 29(5): 920) [5] Yang J N, Wu J C, Kawashima K, et al. Hybrid control of seis鄄 mic鄄excited bridge structure. Earthquake Eng Struct Dyn, 1995, 24(11): 1437 [6] Yang J N, Wu J C, Agrawal A K. Sliding mode control for seismi鄄 cally excited linear structures. J Eng Mech, 1995, 121 ( 12 ): 1386 [7] Rosier L. Homogeneous Lyapunov function for homogeneous con鄄 tinuous vector field. Syst Control Lett, 1992, 19(6): 467 [8] Kawski M. Stability and nilpotent approximations / / Proceedings of the 27th IEEE Conference on Decision and Control. Austin, 1988: 1244 [9] Hermes H. Asymptotically stabilizing feedback controls. J Differ鄄 ential Equations l, 1991, 92(1):76 [10] Wang S. Study on Non鄄smooth Control Algorithm for Smart Base Isolated Benchmark Building [Dissertation]. Guangzhou: Guan鄄 gzhou University, 2016 (王思. 基于 Benchmark 智能隔震结构的非光滑控制算法的 研究与分析[学位论文]. 广州: 广州大学, 2016) [11] Ma K M. Design of continuous non鄄smooth attitude control laws for spacecraft. J Astronautics, 2012, 33(6): 713 (马克茂. 航天器连续非光滑姿态控制律设计. 宇航学报, 2012, 33(6): 713) [12] Ma K M. Non鄄smooth design and implementation of high鄄preci鄄 sion guidance laws. J Ball, 2013, 25(2):1 (马克茂. 高精度制导律的非光滑设计与实现. 弹道学报, 2013, 25(2): 1) [13] Wang J H, Wang Q, Ma K M. Non鄄smooth controller design for permanent magnet synchronous motors. Comput Simul, 2016, 33 (3): 227 (王建晖, 王清, 马克茂. 一种永磁同步电机非光滑控制器 设计. 计算机仿真, 2016, 33(3): 227) [14] Li S H, Ding S H, Du H B, et al. Theory and Application of Non鄄smooth Control. Beijing: Science Press, 2013 (李世华, 丁世宏, 都海波, 等. 非光滑控制理论与应用. 北 京: 科学出版社, 2013) ·1101·