+高级检索
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于模态阻尼指标的空间杆系结构作动器位置优化  PDF

  • 万华平 1,2
  • 马强 1,2
  • 薛宇 1,2
  • 胡鹏华 1,2
  • 葛荟斌 1,2
  • 罗尧治 1,2
1. 浙江大学 建筑工程学院,浙江 杭州 310058; 2. 浙江省空间结构重点实验室,浙江 杭州 310058

中图分类号: TU311.3TU352V414.2

最近更新:2024-07-29

DOI: 10.16339/j.cnki.hdxbzkb.2024062

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

针对空间杆系结构作动器的位置优化问题,本文基于线性二次型最优控制理论(LQR),建立了空间杆系结构振动的主动控制系统,提出了目标控制模态的模态阻尼最大化的作动器位置优化准则. 采用四角锥网架结构和Levy型索穹顶结构来验证本文方法的有效性,并在控制效果和控制能量输入两方面与常用的模态控制力方法进行了对比分析. 结果表明:随着作动器数量的增加,本文方法的控制效果和控制能量输入比模态控制力方法提升更快;在作动器数量相同或控制能量输入相同时,本文方法对外界激励的控制效果优于模态控制力方法.

随着空间杆系结构设计日趋大型化和柔性化,低阻尼的杆系结构在地震作用、风荷载以及冲击荷载等动态激励作用下振动响应大,振动衰减慢,结构动力问题不容忽视,需要采取一定的措施加以控

1-2.主动振动控制技术改变了传统结构通过自身设计方案来被动抵抗外界荷载的设计思路,在轻量化的杆系结构中应用较为广3-4. 主动控制系统通过在结构中布置传感器监测结构变形状态,由控制器引导嵌入结构系统的作动器输入控制力,消耗结构系统在外界激励下的振动能量,可有效降低结构的振动水平,提高结构的抗振能力和防灾性5. 作动器是结构振动主动控制系统的重要组成部分,作动器的位置配置会直接影响振动控制效6. 因此,有必要对作动器的位置优化进行研究.

近年来,国内外许多学者从能量耗散角度对作动器位置优化开展研究. Li

7对线性二次型最优控制(LQR)的性能指标进行改进,得到了与初始状态无关的优化指标. 徐斌8构建了输入能量、溢出能量和有效阻尼响应时间的多目标评价体系,进行桁架结构作动器的位置优化. 王军9从状态空间方程出发,以可控性Gram矩阵特征值为基础,采用能量吸收率作为作动器位置优化的目标函数. 肖南10以LQR控制理论的性能泛函为基础,定义了作动器布设的敏感度,通过逐步删除结构中敏感度较小的作动器,实现了双层球面网壳结构的作动器位置优化. Angeletti11以作动器功耗的加权和为优化目标,确定了大型网格反射器的支撑结构速度反馈的最优增益以及作动器位置. 杨依领12对作动器位置矩阵进行奇异值分解,兼顾保留模态和截断模态对模态控制力的影响,提出一种综合模态控制力最大化的评价准则,并以挠性梁为例进行了压电作动器的位置优化.Li13提出作动器的最优位置是能产生最大的模态控制力的位置,并采用遗传算法进行了桁架结构的作动器位置优化,该优化准则计算方便、应用较广.

本文沿上述研究思路,从LQR控制理论的反馈增益出发,将模态控制力代入运动方程,推导了模态控制力对受控模态所产生的主动阻尼,提出了目标控制模态的模态阻尼最大化的作动器位置优化准则. 通过与文献[

13]中的优化准则(最大化目标控制模态的模态控制力)进行对比分析,说明本文方法布置作动器在控制效果和控制能量输入方面的优势.

1 基于LQR理论的结构主动控制系统建模

结构动力系统的运动方程为:

Mx¨+Cx˙+Kx=fd+Bsfa (1)

式中:MCK分别为结构质量矩阵、阻尼矩阵和刚度矩阵;x¨x˙x分别为结构加速度、速度和位移响应;fd为外界激励;Bs是作动器的位置矩阵;fa为控制力向量,用来模拟作动器的控制输入.

将运动方程转换为状态空间方程:

X˙=AX+Bfa+HfdY=EX (2)

式中:A=0I-M-1K-M-1C,其中I为单位对角矩阵;B=0M-1BsH=0M-1X=xx˙E为测量矩阵,对于全反馈控制的情况,E=I.

结构振动控制的二次型性能指标函数Z

14

Z=120+(XTQX+faTRfa)dt (3)

式中:Z1=XTQX项表征结构振动能量的大小,Q为半正定的权重矩阵; Z2=faTRfa项表征控制能量输入的大小,R为正定权重矩阵.

主动控制力fa和反馈增益矩阵G为:

fa=-GY (4)
G=R-1BP (5)

式中:P式(6)所示黎卡提代数方程的解.

PA+ATP+Q-PBR-1BTP=0 (6)

2 主动控制系统作动器位置优化

作动器的位置对结构主动振动控制系统的性能影响较大. 本文从LQR控制理论的反馈增益矩阵出发,使目标控制模态的模态阻尼最大化,以此作为作动器位置的优化准则.

2.1 优化准则

式(1)作模态变换,令x=Φq,可得到模态坐标系下的系统运动方程:

M¯q¨+C¯q˙+K¯q=ΦTfd+fma (7)

式中:Φ为模态矩阵;q为模态坐标向量;M¯=ΦTMΦ为模态质量矩阵;C¯=ΦTCΦ为模态阻尼矩阵;K¯=ΦTKΦ为模态刚度矩阵;fma=ΦTBsfa为模态控制力向量.

文献[

13]中作动器位置的优化准则为矩阵fmaTfma的特征值集中分布且最大,以此产生最大的模态控制力. 本文将式(4)中的反馈增益矩阵G分为位移反馈增益矩阵G1和速度反馈增益矩阵G2

fa=-GY=-[G1G2]xx˙ (8)

并将式(8)代入运动方程式(7),可将主动控制力转化为对模态阻尼和模态刚度的贡献:

M¯q¨+(C¯+C¯add)q˙+(K¯+K¯add)q=ΦTfd (9)

式中:C¯add=ΦTBsG2ΦT为主动模态阻尼矩阵,由速度反馈增益提供;K¯add=ΦTBsG1ΦT为主动模态刚度矩阵,由位移反馈增益提供. 由式(9)可知,基于LQR控制理论设计的控制力改变了结构的模态阻尼和模态刚度,以主动模态阻尼C¯add和主动模态刚度K¯add的形式施加在结构上,来降低结构振动水平. 文献[

14]指出LQR控制主要以改变结构阻尼为主,因此本文忽略主动模态刚度的影响.

作动器输入控制力所提供的主动模态阻尼C¯add越大,系统对外界激励的耗能能力就越强,对结构的振动控制就越有利. 空间杆系结构模态阶数较多,低阶模态在受迫振动时占主导地位,高阶模态不易被激发,可将全部模态Φ划分为目标控制模态Φc和剩余模态Φr. 记目标控制模态的模态阻尼矩阵为:

C¯call=C¯c+C¯cadd=ΦcTCΦcT+ΦcTBsG2ΦcT (10)

为避免矩阵C¯cadd奇异影响,本文取矩阵C¯call的特征值集中分布且最大时,作动器能提供最大的主动模态阻尼,对外界激励的耗能能力最强,因此作动器的最优位置集合p*为:

p*=argmaxp J(p) (11)
J=i=1ncλii=1nc(λi)nc (12)

式中:p为作动器的位置集合;J为本文作动器位置优化的目标函数;nc为目标控制模态Φc的数目;λi为矩阵C¯call的特征值向量λ的第i个元素;i=1ncλi项为矩阵C¯call所有特征值的算术平均值;i=1nc(λi)nc项为矩阵C¯call所有特征值的几何平均值,两项乘积最大可以保证矩阵C¯call所有特征值分布集中且最大.

2.2 序列前向选择优化方法

作动器位置优化为离散优化问题,该类优化问题求解已有许多成熟的方法. 本文采用序列前向选择方法,它是一种经典的贪心算法,通过寻找当前最优子集的方式来逐步获得最终寻优结

15.

对于作动器位置的选择,设L表示所选作动器位置的当前集合,Ω表示所有可供选择的位置集合,Λ表示当前未选位置的集合,na表示作动器的数量. 基于序列前向选择方法的作动器位置优化流程如 图1所示.

fig

图1  基于序列前向选择方法的作动器位置优化流程图

Fig.1  Flowchart of actuator placement optimization based on sequence forward selection method

3 数值验证

3.1 算例1:四角锥网架结构

四角锥网架结构(图2)用来验证本文作动器位置优化准则的有效性,该结构跨度为10 m,高度为0.6 m,由224根杆件组成,约束为周边固支. 四角锥网架结构具体参数如表1所示.

fig

(a)  模型图

fig

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

fig

(c)  侧视图(单位:m)

图2  四角锥网架

Fig.2  Square pyramid space grid

表1  四角锥网架结构参数
Tab.1  Structure parameters of square pyramid space grid
结构参数参数取值
总节点数 65
自由节点数 49
杆单元数 224
节点集中质量/kg 2
弹性模量/GPa 206
材料密度/(kg·m-3 7 850
杆件等效截面面积/m2 2×10-4
阻尼比 0.05

采用文献[

13]中的模态控制力方法与本文方法进行对比分析. 目标控制模态数nc=10,作动器数量na=4. 基于LQR理论的主动控制系统中QR的取值为:

Q=I147×1470147×1470147×1470147×147
R=3×10-14I4×4

除去结构周边被约束的16根杆件,考虑结构中间的208个单元位置为作动器可供选择的位置集合. 采用前述序列前向选择方法,当作动器数量na=4时,模态控制力方法的作动器位置优化结果如图3(a)所示,本文方法的作动器位置优化结果如图3(b)所示.

fig

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

fig

(b)  本文方法作动器位置优化结果

图3  作动器位置示意图(四角锥网架)

Fig.3  Schematic diagram of actuators placement

(square pyramid space grid)

为比较两种方法所选作动器位置对外界激励的振动控制效果,设计了两种荷载工况:

1)地震作用:结构竖向地震动输入为El-Centro地震波,峰值加速度取3g,作用时间取31.2 s.采用 式(13)将地震动转化为作用力施加到结构上.

fd=-MIgx¨g (13)

式中:Ig是影响向量;x¨g是地震动加速度.

2)风荷载:采用AR自回归模

16模拟风速时程,水平阵风采用Daveport谱,AR模型阶数取5,基本风压取0.45 kN/m2,模拟风速时长60 s. 采用

式(14)将风速时程转化为作用力施加到结构上弦节点.

fd=12CpρApV2 (14)

式中:Cp为风压分布系数,取-1;ρ为空气质量密度,取1.225 kg/m³;Ap为作用面积向量;V为风速时程向量.

为更加全面地展示模态控制力方法和本文方法对外界激励的振动控制效果,取控制性能指标Z的第一项Z1=XTQX来表征结构振动能量的大小. 图4为作动器数量na=4时,模态控制力方法和本文方法布置作动器与无控情况下结构振动能量Z1的对比图. 由图4可知,依据本文方法布置作动器,在地震作用和风荷载激励下均表现出良好的振动控制性能,结构振动能量Z1降低更为明显,优于模态控制力方法.

fig

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

fig

(b)  结构振动能量(风荷载)

图4  作动器数量na=4时结构振动能量对比图(四角锥网架)

Fig.4  Comparison of vibration energy with four actuators (square pyramid space grid)

为量化对比两种方法的控制效果,定义结构振动能量降低率η为:

η=Z1,无控-Z1,有控Z1,无控 (15)

式(15)可知,结构振动能量降低率η越大,则控制效果越好. 在地震作用和风荷载两种工况下,模态控制力方法的结构振动能量降低率η约为20%和13%,而本文方法约为46%和38%. 显然,在均采用 4个作动器时,本文方法的控制效果更优.

主动控制全程的能量消耗也是评价控制系统性能的重要指标,采用控制性能指标Z的第二项的积分Z2=faTRfa来表征整个控制过程中的能量输入. 采用序列前向选择方法,按照模态控制力方法和本文方法布置不同数量的作动器,对比控制效果与控制能量输入,结果如图5所示. 由图5可知,在两种荷载工况下,随着作动器数量增加,本文方法的控制效果提升较快,结构振动能量降低率η更快稳定在较高水平.当作动器数量na100时,本文方法的控制效果优于模态控制力方法;当作动器数量na>100时,本文方法的控制效果与模态控制力方法接近. 随着作动器数量的增加,本文方法的控制能量输入先增加后下降,最终随着控制效果的更快稳定,控制能量输入也更快达到稳定.当作动器数量na<20时,本文方法可用较高的控制能量输入维持更好的控制效果;在na=20附近,两种方法的控制能量输入相同,本文方法的控制效果优于模态控制力方法;随着作动器数量的增加,本文方法所需的能量输入更快开始下降,并小于模态控制力方法. 作动器的位置会影响主动控制力提供给结构的主动阻尼,主动阻尼越大,对外界激励的耗能能力就越强,因此在模态阻尼矩阵的特征值集中分布且最大时,控制效果与模态控制力方法相比更优.

fig

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

fig

(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型索穹顶结构(图6)用来进一步验证本文方法的有效性,该结构跨度为20 m,由96根拉索、 17根压杆,共113个单元组成,约束为周边固支,具体结构参数如表2所示. Levy型索穹顶结构具有良好的对称性,可分为11组杆件,采用二次奇异值分解方

17确定索穹顶结构整体可行预应力,预应力fp设置详见表3.

fig

(a)  模型图

fig

(b)  俯视图

fig

(c)  剖面图(单位:m)

图6  Levy型索穹顶

Fig.6  Levy-type cable dome

表2  Levy型索穹顶结构参数
Tab.2  Structure parameters of Levy-type cable dome
结构参数参数取值
总节点数 42
自由节点数 34
单元数 113
节点集中质量/kg 5
拉索弹性模量/GPa 181
拉索等效截面面积/m2 6×10-4
压杆弹性模量/GPa 206
压杆等效截面面积/m2 1×10-3
材料密度/(kg·m-3 7 850
阻尼比 0.05
表3  Levy型索穹顶预应力设置
Tab.3  Prestress setting of Levy-type cable dome
单元分组编号预应力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

目标控制模态数nc=8,作动器数量na=4. 基于LQR理论的结构主动控制系统中QR的取值为:

Q=I102×1020102×1020102×1020102×102
R=2×10-14I4×4

考虑96根拉索所在位置为作动器可供选择的位置集合. 采用序列前向选择方法,模态控制力方法和本文方法的作动器位置优化结果如图7所示.

fig

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

fig

(b)  本文方法作动器位置优化结果

图7  作动器位置示意图(Levy型索穹顶)

Fig.7  Schematic diagram of actuators placement (Levy-type cable dome)

荷载工况设置同算例1一致. 索穹顶结构由预应力提供刚度保持稳定,因此在振动控制过程中,作动器所在单元应避免预应力松弛,故预应力fp和控制力fa需满足如下约束:

fp+fa>0 (16)

式中:fpfa以受拉为正,受压为负.

图8为作动器数量na=4时,依据模态控制力方法和本文方法布置作动器和无控情况下结构振动能量的对比图. 由图8可知,在地震作用和风荷载两种工况下,依据本文方法在外圈环索布置作动器的振动控制效果更好. 具体地,在地震作用和风荷载激励下,模态控制力方法的结构振动能量降低率约为10%和8%,而本文方法约为33%和30%.

fig

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

fig

(b)  结构振动能量(风荷载)

图8  作动器数量na=4时结构振动能量对比图

Fig.8  Comparison of vibration energy with four actuators

(Levy型索穹顶)

(Levy-type cable dome)

图9为模态控制力方法和本文方法在布置不同数量作动器的情况下,控制效果与控制能量输入的对比图. 图9中曲线的变化趋势同算例1的图5一致,可知本文方法在作动器数量较少时,可用较高的控制能量输入维持更好的控制效果;在控制能量输入相同时,本文方法的控制效果优于模态控制力方法;随着作动器数量的继续增加,作动器提供给结构的主动阻尼增加了结构对外界激励的耗能能力,对控制效果的影响持续增强,因此本文方法的控制效果优于模态控制力方法.

fig

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

fig

(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)在作动器数量相同或控制能量输入相同时,本文方法的控制效果均优于模态控制力方法.

参考文献

1

QIN W HGAO HXI Zet alShaking table experimental investigations on dynamic characteristics of CFRP cable dome[J].Engineering Structures2023281115748 [百度学术] 

2

丁阳赵奕程大跨度空间钢管桁架结构的风振响应和风振控制研究[J].湖南大学学报(自然科学版)2006333):17-21 [百度学术] 

DING YZHAO Y CStudy of wind-induced dynamic response and wind vibration control for large-span spatial lattice structure[J].Journal of Hunan University (Natural Sciences)2006333):17-21(in Chinese) [百度学术] 

3

SENATORE GDUFFOUR PWINSLOW Pet alShape control and whole-life energy assessment of an ‘infinitely stiff’ prototype adaptive structure[J].Smart Materials and Structures2018271):015022 [百度学术] 

4

PENG H JLI FKAN Z YA novel distributed model predictive control method based on a substructuring technique for smart tensegrity structure vibrations[J].Journal of Sound Vibration2020471115171 [百度学术] 

5

蒋伟刘纲振动控制方法在大跨度空间结构中的应用与研究进展[J].智能建筑与智慧城市20206):9-14 [百度学术] 

JIANG WLIU GThe advance and application of vibration control method in long-span spatial structures[J].Intelligent Building & Smart City20206):9-14(in Chinese) [百度学术] 

6

潘兆东谭平周福霖基于一种新混合群算法(MOHO-SA)的结构主动控制多目标优化研究[J].湖南大学学报(自然科学版)2017445):20-26 [百度学术] 

PAN Z DTAN PZHOU F LMulti-objective optimization of structural active control system using a new hybrid swarm algorithm[J].Journal of Hunan University (Natural Sciences)2017445):20-26(in Chinese) [百度学术] 

7

LI W PHUANG HIntegrated optimization of actuator placement and vibration control for piezoelectric adaptive trusses[J].Journal of Sound and Vibration20133321):17-32 [百度学术] 

8

徐斌欧进萍姜节胜压电智能桁架作动器配置的优化设计[J].机械科学与技术20072611):1511-1516 [百度学术] 

XU BOU J PJIANG J SOptimal number and locations of actuators for piezoelectric intelligent trusses[J].Mechanical Science and Technology for Aerospace Engineering20072611):1511-1516(in Chinese) [百度学术] 

9

王军杨亚东张家应面向结构振动控制的压电作动器优化配置研究[J].航空学报2012333):494-500 [百度学术] 

WANG JYANG Y DZHANG J Yet alInvestigation of piezoelectric actuator optimal configuration for structural vibration control[J].Acta Aeronautica et Astronautica Sinica2012333):494-500(in Chinese) [百度学术] 

10

肖南容里董石麟双层球面网壳振动主动控制作动器位置优化[J].浙江大学学报(工学版)2010445):942-949 [百度学术] 

XIAO NRONG LDONG S LOptimal placements of actuators in vibration active control for double-layer spherical grid shell structures[J].Journal of Zhejiang University (Engineering Science)2010445):942-949(in Chinese) [百度学术] 

11

ANGELETTI FGASBARRI PSABATINI MOptimal design and robust analysis of a net of active devices for micro-vibration control of an on-orbit large space antenna[J].Acta Astronautica2019164241-253 [百度学术] 

12

杨依领娄军强魏燕定综合模态控制力下压电致动器的优化布局[J].浙江大学学报(工学版)2015495):841-847 [百度学术] 

YANG Y LLOU J QWEI Y Det alOptimal placement of piezoelectric actuators using synthetic modal control force[J].Journal of Zhejiang University (Engineering Science)2015495):841-847(in Chinese) [百度学术] 

13

LI D XLIU WJIANG J Pet alPlacement optimization of actuator and sensor and decentralized adaptive fuzzy vibration control for large space intelligent truss structure[J].Science China Technological Sciences2011544):853-861 [百度学术] 

14

欧进萍结构振动控制:主动、半主动和智能控制[M].北京科学出版社2003 [百度学术] 

OU J PStructural vibration control:active,semi-active and intelligent control[M].BeijingScience Press2003(in Chinese) [百度学术] 

15

XUE YLUO Y ZXU Xet alA robust method for pre-stress adjustment of cable-strut structures based on sparse regression[J].Engineering Structures2021246112987 [百度学术] 

16

POGGI PMUSELLI MNOTTON Get alForecasting and simulating wind speed in Corsica by using an autoregressive model[J].Energy Conversion and Management20034420):3177-3196 [百度学术] 

17

YUAN X FCHEN L MDONG S LPrestress design of cable domes with new forms[J].International Journal of Solids and Structures2007449):2773-2782 [百度学术] 

作者稿件一经被我刊录用,如无特别声明,即视作同意授予我刊论文整体的全部复制传播的权利,包括但不限于复制权、发行权、信息网络传播权、广播权、表演权、翻译权、汇编权、改编权等著作使用权转让给我刊,我刊有权根据工作需要,允许合作的数据库、新媒体平台及其他数字平台进行数字传播和国际传播等。特此声明。
关闭