摘要
针对空间杆系结构作动器的位置优化问题,本文基于线性二次型最优控制理论(LQR),建立了空间杆系结构振动的主动控制系统,提出了目标控制模态的模态阻尼最大化的作动器位置优化准则. 采用四角锥网架结构和Levy型索穹顶结构来验证本文方法的有效性,并在控制效果和控制能量输入两方面与常用的模态控制力方法进行了对比分析. 结果表明:随着作动器数量的增加,本文方法的控制效果和控制能量输入比模态控制力方法提升更快;在作动器数量相同或控制能量输入相同时,本文方法对外界激励的控制效果优于模态控制力方法.
随着空间杆系结构设计日趋大型化和柔性化,低阻尼的杆系结构在地震作用、风荷载以及冲击荷载等动态激励作用下振动响应大,振动衰减慢,结构动力问题不容忽视,需要采取一定的措施加以控
近年来,国内外许多学者从能量耗散角度对作动器位置优化开展研究. Li
本文沿上述研究思路,从LQR控制理论的反馈增益出发,将模态控制力代入运动方程,推导了模态控制力对受控模态所产生的主动阻尼,提出了目标控制模态的模态阻尼最大化的作动器位置优化准则. 通过与文献[
1 基于LQR理论的结构主动控制系统建模
结构动力系统的运动方程为:
(1) |
式中:、和分别为结构质量矩阵、阻尼矩阵和刚度矩阵;、和分别为结构加速度、速度和位移响应;为外界激励;是作动器的位置矩阵;为控制力向量,用来模拟作动器的控制输入.
将运动方程转换为状态空间方程:
(2) |
式中:,其中为单位对角矩阵;;;;为测量矩阵,对于全反馈控制的情况,.
结构振动控制的二次型性能指标函数
(3) |
式中:项表征结构振动能量的大小,为半正定的权重矩阵; 项表征控制能量输入的大小,为正定权重矩阵.
主动控制力和反馈增益矩阵为:
(4) |
(5) |
式中:为
(6) |
2 主动控制系统作动器位置优化
作动器的位置对结构主动振动控制系统的性能影响较大. 本文从LQR控制理论的反馈增益矩阵出发,使目标控制模态的模态阻尼最大化,以此作为作动器位置的优化准则.
2.1 优化准则
对
(7) |
式中:为模态矩阵;为模态坐标向量;为模态质量矩阵;为模态阻尼矩阵;为模态刚度矩阵;为模态控制力向量.
文献[
(8) |
并将
(9) |
式中:为主动模态阻尼矩阵,由速度反馈增益提供;为主动模态刚度矩阵,由位移反馈增益提供. 由
作动器输入控制力所提供的主动模态阻尼越大,系统对外界激励的耗能能力就越强,对结构的振动控制就越有利. 空间杆系结构模态阶数较多,低阶模态在受迫振动时占主导地位,高阶模态不易被激发,可将全部模态划分为目标控制模态和剩余模态. 记目标控制模态的模态阻尼矩阵为:
(10) |
为避免矩阵奇异影响,本文取矩阵的特征值集中分布且最大时,作动器能提供最大的主动模态阻尼,对外界激励的耗能能力最强,因此作动器的最优位置集合为:
(11) |
(12) |
式中:为作动器的位置集合;为本文作动器位置优化的目标函数;为目标控制模态的数目;为矩阵的特征值向量的第个元素;项为矩阵所有特征值的算术平均值;项为矩阵所有特征值的几何平均值,两项乘积最大可以保证矩阵所有特征值分布集中且最大.
2.2 序列前向选择优化方法
作动器位置优化为离散优化问题,该类优化问题求解已有许多成熟的方法. 本文采用序列前向选择方法,它是一种经典的贪心算法,通过寻找当前最优子集的方式来逐步获得最终寻优结
对于作动器位置的选择,设表示所选作动器位置的当前集合,表示所有可供选择的位置集合,表示当前未选位置的集合,表示作动器的数量. 基于序列前向选择方法的作动器位置优化流程如

图1 基于序列前向选择方法的作动器位置优化流程图
Fig.1 Flowchart of actuator placement optimization based on sequence forward selection method
3 数值验证
3.1 算例1:四角锥网架结构
四角锥网架结构(

(a) 模型图

(b) 俯视图(单位:m)

(c) 侧视图(单位:m)
图2 四角锥网架
Fig.2 Square pyramid space grid
结构参数 | 参数取值 |
---|---|
总节点数 | 65 |
自由节点数 | 49 |
杆单元数 | 224 |
节点集中质量/kg | 2 |
弹性模量/GPa | 206 |
材料密度/(kg· | 7 850 |
杆件等效截面面积/ |
2×1 |
阻尼比 | 0.05 |
采用文献[
除去结构周边被约束的16根杆件,考虑结构中间的208个单元位置为作动器可供选择的位置集合. 采用前述序列前向选择方法,当作动器数量时,模态控制力方法的作动器位置优化结果如

(a) 模态控制力方法作动器位置优化结果

(b) 本文方法作动器位置优化结果
图3 作动器位置示意图(四角锥网架)
Fig.3 Schematic diagram of actuators placement
(square pyramid space grid)
为比较两种方法所选作动器位置对外界激励的振动控制效果,设计了两种荷载工况:
1)地震作用:结构竖向地震动输入为El-Centro地震波,峰值加速度取3g,作用时间取31.2 s.采用
(13) |
式中:是影响向量;是地震动加速度.
2)风荷载:采用AR自回归模
(14) |
式中:为风压分布系数,取-1;为空气质量密度,取1.225 kg/m³;为作用面积向量;为风速时程向量.
为更加全面地展示模态控制力方法和本文方法对外界激励的振动控制效果,取控制性能指标的第一项来表征结构振动能量的大小.

(a) 结构振动能量(地震作用)

(b) 结构振动能量(风荷载)
图4 作动器数量na=4时结构振动能量对比图(四角锥网架)
Fig.4 Comparison of vibration energy with four actuators (square pyramid space grid)
为量化对比两种方法的控制效果,定义结构振动能量降低率为:
(15) |
由
主动控制全程的能量消耗也是评价控制系统性能的重要指标,采用控制性能指标的第二项的积分来表征整个控制过程中的能量输入. 采用序列前向选择方法,按照模态控制力方法和本文方法布置不同数量的作动器,对比控制效果与控制能量输入,结果如

(a) 控制效果与控制能量输入对比(地震作用)

(b) 控制效果与控制能量输入对比(风荷载)
图5 不同数量作动器控制效果与控制能量输入对比图
Fig.5 Comparison of control effect and control energy input with different number of actuators (square pyramid space grid)
(四角锥网架)
3.2 算例2:Levy型索穹顶
Levy型索穹顶结构(

(a) 模型图

(b) 俯视图

(c) 剖面图(单位:m)
图6 Levy型索穹顶
Fig.6 Levy-type cable dome
结构参数 | 参数取值 |
---|---|
总节点数 | 42 |
自由节点数 | 34 |
单元数 | 113 |
节点集中质量/kg | 5 |
拉索弹性模量/GPa | 181 |
拉索等效截面面积/ |
6×1 |
压杆弹性模量/GPa | 206 |
压杆等效截面面积/ |
1×1 |
材料密度/(kg· | 7 850 |
阻尼比 | 0.05 |
单元分组编号 | 预应力fp /N |
---|---|
1 | 7 220 |
2 | 3 388 |
3 | 6 968 |
4 | 3 286 |
5 | 16 890 |
6 | 8 807 |
7 | 5 839 |
8 | 12 129 |
9 | -11 327 |
10 | -2 636 |
11 | -5 775 |
目标控制模态数,作动器数量. 基于LQR理论的结构主动控制系统中和的取值为:
考虑96根拉索所在位置为作动器可供选择的位置集合. 采用序列前向选择方法,模态控制力方法和本文方法的作动器位置优化结果如

(a) 模态控制力方法作动器位置优化结果

(b) 本文方法作动器位置优化结果
图7 作动器位置示意图(Levy型索穹顶)
Fig.7 Schematic diagram of actuators placement (Levy-type cable dome)
荷载工况设置同算例1一致. 索穹顶结构由预应力提供刚度保持稳定,因此在振动控制过程中,作动器所在单元应避免预应力松弛,故预应力和控制力需满足如下约束:
(16) |
式中:和以受拉为正,受压为负.

(a) 结构振动能量(地震作用)

(b) 结构振动能量(风荷载)
图8 作动器数量na=4时结构振动能量对比图
Fig.8 Comparison of vibration energy with four actuators
(Levy型索穹顶)
(Levy-type cable dome)

(a) 控制效果与控制能量输入对比(地震作用)

(b) 控制效果与控制能量输入对比(风荷载)
图9 不同数量作动器控制效果与控制能量输入对比图
Fig.9 Comparison of control effect and control energy input with different number of actuators (Levy-type cable dome)
(Levy型索穹顶)
4 结 论
本文研究了基于LQR理论的空间杆系结构振动主动控制系统,提出了目标控制模态的模态阻尼最大化的作动器位置优化准则. 采用四角锥网架结构和Levy型索穹顶结构验证了本文方法的有效性,并与常用的模态控制力方法进行对比,得到以下结论:
1)在作动器数量较少时,本文方法可用较高的控制能量输入维持更好的控制效果.
2)随着作动器数量的增加,本文方法的控制效果和控制能量输入比模态控制力方法提升更快.
3)在作动器数量相同或控制能量输入相同时,本文方法的控制效果均优于模态控制力方法.
参考文献
QIN W H,GAO H,XI Z,et al.Shaking table experimental investigations on dynamic characteristics of CFRP cable dome[J].Engineering Structures,2023,281:115748. [百度学术]
丁阳,赵奕程.大跨度空间钢管桁架结构的风振响应和风振控制研究[J].湖南大学学报(自然科学版),2006,33(3):17-21. [百度学术]
DING Y,ZHAO Y C.Study of wind-induced dynamic response and wind vibration control for large-span spatial lattice structure[J].Journal of Hunan University (Natural Sciences),2006, 33(3):17-21.(in Chinese) [百度学术]
SENATORE G,DUFFOUR P,WINSLOW P,et al.Shape control and whole-life energy assessment of an ‘infinitely stiff’ prototype adaptive structure[J].Smart Materials and Structures,2018, 27(1):015022. [百度学术]
PENG H J,LI F,KAN Z Y.A novel distributed model predictive control method based on a substructuring technique for smart tensegrity structure vibrations[J].Journal of Sound Vibration,2020,471:115171. [百度学术]
蒋伟,刘纲.振动控制方法在大跨度空间结构中的应用与研究进展[J].智能建筑与智慧城市,2020(6):9-14. [百度学术]
JIANG W,LIU G.The advance and application of vibration control method in long-span spatial structures[J].Intelligent Building & Smart City,2020(6):9-14.(in Chinese) [百度学术]
潘兆东,谭平,周福霖.基于一种新混合群算法(MOHO-SA)的结构主动控制多目标优化研究[J].湖南大学学报(自然科学版),2017,44(5):20-26. [百度学术]
PAN Z D,TAN P,ZHOU F L.Multi-objective optimization of structural active control system using a new hybrid swarm algorithm[J].Journal of Hunan University (Natural Sciences),2017,44(5):20-26.(in Chinese) [百度学术]
LI W P,HUANG H.Integrated optimization of actuator placement and vibration control for piezoelectric adaptive trusses[J].Journal of Sound and Vibration,2013,332(1):17-32. [百度学术]
徐斌,欧进萍,姜节胜.压电智能桁架作动器配置的优化设计[J].机械科学与技术,2007,26(11):1511-1516. [百度学术]
XU B,OU J P,JIANG J S.Optimal number and locations of actuators for piezoelectric intelligent trusses[J].Mechanical Science and Technology for Aerospace Engineering,2007, 26(11):1511-1516.(in Chinese) [百度学术]
王军,杨亚东,张家应,等.面向结构振动控制的压电作动器优化配置研究[J].航空学报,2012,33(3):494-500. [百度学术]
WANG J,YANG Y D,ZHANG J Y,et al.Investigation of piezoelectric actuator optimal configuration for structural vibration control[J].Acta Aeronautica et Astronautica Sinica,2012, 33(3):494-500.(in Chinese) [百度学术]
肖南,容里,董石麟.双层球面网壳振动主动控制作动器位置优化[J].浙江大学学报(工学版),2010,44(5):942-949. [百度学术]
XIAO N,RONG L,DONG S L.Optimal placements of actuators in vibration active control for double-layer spherical grid shell structures[J].Journal of Zhejiang University (Engineering Science),2010,44(5):942-949.(in Chinese) [百度学术]
ANGELETTI F,GASBARRI P,SABATINI M.Optimal design and robust analysis of a net of active devices for micro-vibration control of an on-orbit large space antenna[J].Acta Astronautica,2019,164:241-253. [百度学术]
杨依领,娄军强,魏燕定,等.综合模态控制力下压电致动器的优化布局[J].浙江大学学报(工学版),2015,49(5):841-847. [百度学术]
YANG Y L,LOU J Q,WEI Y D,et al.Optimal placement of piezoelectric actuators using synthetic modal control force[J].Journal of Zhejiang University (Engineering Science),2015, 49(5):841-847.(in Chinese) [百度学术]
LI D X,LIU W,JIANG J P,et al.Placement optimization of actuator and sensor and decentralized adaptive fuzzy vibration control for large space intelligent truss structure[J].Science China Technological Sciences,2011,54(4):853-861. [百度学术]
欧进萍.结构振动控制:主动、半主动和智能控制[M].北京:科学出版社,2003. [百度学术]
OU J P.Structural vibration control:active,semi-active and intelligent control[M].Beijing:Science Press,2003.(in Chinese) [百度学术]
XUE Y,LUO Y Z,XU X,et al.A robust method for pre-stress adjustment of cable-strut structures based on sparse regression[J].Engineering Structures,2021,246:112987. [百度学术]
POGGI P,MUSELLI M,NOTTON G,et al.Forecasting and simulating wind speed in Corsica by using an autoregressive model[J].Energy Conversion and Management,2003,44(20):3177-3196. [百度学术]
YUAN X F,CHEN L M,DONG S L.Prestress design of cable domes with new forms[J].International Journal of Solids and Structures,2007,44(9):2773-2782. [百度学术]