网刊加载中。。。

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

确定继续浏览么?

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

用于结构减震控制的拉索式惯容系统跨层布置优化研究

  • 薛松涛 1,2
  • 康建飞 1
  • 谢丽宇 1
1. 同济大学 结构防灾减灾工程系,上海 200092; 2. 日本东北工业大学 建筑系,仙台 982-8577

最近更新:2020-11-23

DOI:10.6052/1672-6553-2020-065

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

摘要

本文提出了一种拉索式惯容系统的跨层布置方法用于结构减震控制.将地震激励模拟为金井清谱随机激励,构建统一的状态空间方程,进行线性系统直接随机分析,获得系统响应统计量.利用布置效率确定惯容系统布置位置,考虑惯容系统整体出力影响惯容系统自身造价及结构柱额外受力,将惯容系统出力状况作为经济指标.考虑结构位移响应及加速度响应与结构性能的相关性,将结构的位移响应及加速度响应作为减震效果指标.利用设计参数的多目标优化,获取惯容系统最优参数的帕累托前沿,指导设计并比较不同布置方法.最后,利用一个10层的基准结构进行减震效果分析,验证跨层布置的惯容系统优于一般层间惯容系统.

引言

近年来,大量的振动控制装置被提

12,并应用于工程实际.在这些振动控制装置中,调谐质量阻尼器(TMD)设计简单,振动控制机构简洁,由质量元件、弹簧和阻尼器三个经典机械部件组成.TMD被证明在控制风荷载引起的结构振动方面十分有3,但是由于其调谐质量的大小受工程实际限制,TMD对地震激励下的结构振动控制方面效果较4.21世纪初,用于类比电路中电容的双端惯性单元概念的提5,为减小阻尼器的重量和尺寸提供了契机,惯性单元产生的惯性力与两个端子之间的相对加速度成正比,使得惯性单元能够利用楼层之间的加速度差进行减振.一种新型的民用结构惯容系统被广泛研6-10.在土木工程中,早在1999年,Arakaki等11就提出并独立实践了类似惯性器的概念,利用滚珠丝杠机构放大粘性阻尼器的有效输出力,用于抑制地震引起的振动.这是惯性阻尼器在土木工程中的首次应用.

针对惯容系统的设计及优化已经有许多学者进行了研究.对于单自由度结构,Ikago等

6基于不动点理论推导出了惯容系统优化设计的简单公式.Pan等912考虑单自由度(SDOF)结构固有的阻尼能力,提出了一种基于需求的惯容系统优化设计方法,以最小的控制成本满足性能要求.对于多自由度结构,Ikago等7研究了安装惯容系统的多自由度结构模态响应并给出了简化设计方法.Taflanidis等10提出了基于惯性放大装置的多目标优化方法.

本文研究了拉索式惯容系统的优化布置,提出拉索式惯容系统的跨层布置.第一节概述了拉索式惯容系统的基本组成.第二节建立安装有拉索式惯容系统的多自由框架结构分析模型.通过白噪声激励的拓展状态方程及观测方程,给出利用直接随机分析确定响应方差的方法.第三节提出惯容系统布置效率指标,确定惯容系统安装位置.考虑惯容系统的经济指标及减震效果指标,将惯容系统多目标优化转化为数学问题.第四节,利用一个十层的基准模型进行算例分析,研究层间布置、跨一层布置及跨两层布置的拉索式惯容系统的减震效率.分析不同减震指标对三种布置方案的位置及参数优化影响,并从阻尼需求、惯容质量需求、系统出力大小等多方面探讨拉索式惯容系统跨层布置优势.

1 拉索式惯容系统的基本组成

图1所示,该拉索式惯容系统由电涡流阻尼单元、旋转导体盘惯容单元及弹性拉索调谐单元组成.

图1 拉索式惯容系统示意图

Fig.1 Cable-bracing inertial system

在电涡流阻尼单元工作时,导体盘绕轴承转动,切割永磁铁静磁场的磁感线,在导体内产生涡电流,从而在与运动方向相反的方向发生磁性制动力,即电涡流阻尼力.如图2所示,旋转导体盘惯容单元是利用拉索平动位移向导体盘转动位移转换时的放大机制获取惯容效应的单元.弹性拉索调谐单元由刚性拉索及调谐弹簧串联组成,通过调整弹簧刚度实现惯容系统的调谐.

图2 电涡流阻尼单元示意图

Fig.2 The eddy current damping element

2 布置拉索式惯容系统的多层结构模型

考虑一个n自由度的多层框架结构,布置了拉索式惯容系统.如图3所示,第ib层和第id层分别为惯容阻尼单元安装层和拉索端部锚固层(id-ib1).

图3 布置拉索式惯容系统的多层结构示意图

Fig.3 Multi-story structure equipped with cable-bracing inertial system

拉索式惯容系统以倾斜角度θ直接安装于结构中时,其对上下层的水平方向控制力可通过下式计算:

fd(t)=kdy(t)cos2θ (1)
fb(t)=mdRcTx¨s(t)cos2θ+y¨(t)cos2θ+cdRcTx˙s(t)cos2θ+y˙(t)cos2θ (2)

其中,kd为弹簧刚度,md为惯容质量,cd为电涡流阻尼系数,y(t)为弹簧单元相对于id层的水平位移.xs(t)为结构各层位移矢量,RcT为第ib层和第id层相对位移关于各层位移矢量的转换矩阵.

根据(1)、(2)两式可知,拉索式惯容系统的水平控制力与其安装方式(如安装角度θ等)密切相关,故定义安装效率系数β.采用拉索斜向倾角为θ的直接安装方式时,β=cos2θ.引入等效设计参数,ke=βkdme=βmdce=βcd,水平方向控制力方程可写为:

fd(t)=key(t) (3)
fb(t)=meRcTx¨s(t)+y¨(t)+ceRcTx˙s(t)+y˙(t) (4)

布置拉索式惯容系统的多层框架结构动力方程可表达为:

(Ms+RcmeRcT)x¨s(t)+(Rcme)y¨(t)+(Cs+RcceRcT)x˙s(t)+Rccey˙(t)+Κsxs(t)=-MsRsx¨g(t) (5)
mey¨(t)+RcTx¨s(t)+cey˙(t)+RcTx˙s(t)+key(t)=0 (6)

其中,MsCsKs分别为结构的质量矩阵,阻尼矩阵和刚度矩阵,Rs为地震响应系数向量.x¨g(t)为地面运动加速度.

地震激励模拟为金井清

13

Sg(ω)=ωg4+4ζg2ω2ωg2ωg2-ω22+4ζg2ωg2ω2S0 (7)

其中,Sg,S0分别为地震激励的功率频谱密度和功率谱密度曲线的幅值系数,ωg,ζg分别为金井清谱中场地卓越周期圆频率和金井清谱中场地阻尼比.

地面激励的状态空间表达式:

x˙q(t)=Aqxq(t)+Eqw(t)x¨g(t)=Cqxq(t) (8)

其中,Aq=01-ωg2-2ζgωgEq=01

Cq=2πS0-ωg2-2ζgωg

地面激励模型(8)式与多层结构动力方程(5)、(6)式相结合,以x=xsTyx˙sTy˙xqTT为状态向量,可将布置拉索式惯容系统的多层结构分析模型表达为如下拓展的状态空间方程:

x˙(t)=Ax(t)+Ew(t)z(t)=Cx(t) (9)

其中,z(t)为输出向量,包含结构的相对位移响应,绝对加速度响应,惯容系统水平力的输出.矩阵ACE定义如下:

A=            0(n+1)×(n+1)            I(n+1)×(n+1)-Mt-1Ks0n×1Cs+RcceRcTRcce01×nkeceRcTceEsCq02×(2n+2)Aq
C=           Ts                           0n×(n+4)-In  0n×1Mt-1Ks0n×1Cs+RcceRcTRcce01×nkeceRcTce  0n×2               01×n              kd           01×(n+3)      E=0(2n+2)×1    EqEs=        0(n+1)×1-Mt-1MsRs  01×1
Mt=Ms+RcmeRcTRcmemeRcTme

Ts为层间位移传递矩阵,是一个维数为n的三对角方阵,主对角线上元素均为1,高对角线上元素均为-1,低对角线上元素均为0.

依据上述分析模型,可知系统输出z(t)是一个零均值的高斯随机向量,其协方差矩阵为:

Kzz=CPCT (10)

其中,P是状态协方差矩阵可通过求解如式(11)所示的李雅普诺夫方程得

14.

AP+PAT+EET=0 (11)

结构的相对位移响应方差,加速度响应方差及惯容系统水平力方差作为系统观测量均可由系统输出协方差矩阵Kzz直接求得.

3 跨层布置多目标优化

本文采用标准化的结构位移方差及加速度方差作为减震效果指标,将拉索式惯容系统的水平出力方差作为惯容系统的经济指标对惯容器设计参数kd,md,cd进行优化设计.

定义减震效果指标JH2及惯容系统经济指标Jfd

JH2=i=12nρi2σzi2σz0i2,    Jfd=σfd (12)

其中,σzi2,σz0i2分别为第i个受控系统输出的方差和未控系统输出的方差,包含结构相对位移和加速度,ρi为第i个系统输出的权重和,σfd为惯容系统的水平控制力标准差.

考虑仅位移响应为控制目标时权重因子ρi可表达为:

ρi=0.1i=1,2,,,100i=11,12,,,20 (13)

考虑仅加速度响应为控制目标时权重因子ρi可表达为:

ρi=0i=1,2,,,100.1i=11,12,,,20 (14)

同时考虑位移响应及加速度响应时权重因子ρi可表达为:

ρi=0.05i=1,2,,,20 (15)

针对设计参数kd,md,cd,定义无量纲频率比rd,阻尼比ζd,质量比μ

rd=ωdω1,ζd=ce2M1ω1,μ=meM1 (16)

其中,ωd,ω1,M1分别为惯容系统圆频率,结构基本圆频率及结构第一模态质量.

引入布置效率指标:

JE(ib)=1min rdminrdrdmaxζdminζdζdmax  μminμμmaxJH2 (17)

其中,rdmin,ζdmin,μmin分别为设计变量rd,ζd,μ的下界,rdmax,ζdmax,μmax分别为设计变量rd,ζd,μ的上界.由于减震效果JH2中权重因子ρi的不同定义,布置效率指标JE也存在相应的不同定义.

依据布置效率指标JE可初步确定惯容系统布置方式,进一步将拉索式惯容系统的优化问题表述为下面的数学优化问题:

minimizerd,ζd,μJH2, Jfdsubject to  rdminrdrdmaxζdminζdζdmaxμminμμmax (18)

4 十层框架结构算例分析

以日本隔震协会(JSSI)提供的十层基准结构进行算例分析.该基准结构的参数信息及无阻尼基本周期如表1-2所示.

表1 模型的参数信息
Table 1 Specifications of the analytical model
StoreyPrimary structure
mass[t]stiffness[kN/m]height[m]
10 875 158550 4
9 649 180110 4
8 656 220250 4
7 660 244790 4
6 667 291890 4
5 670 306160 4
4 676 328260 4
3 680 383020 4
2 682 383550 4
1 700 279960 6
表2 模型的无阻尼基本周期
Table 2 Undamped fundamental period of the model
ModePeriods/sAngular frequency
1 2.01 3.13
2 0.76 8.27
3 0.46 13.7

利用最不利的地震功率谱(地震谱达到峰值时的圆频率与结构第一模态频率一致)进行优化分析.地震动功率谱曲线如图4所示,金井清谱达到峰值的圆频率为3.13.

图4 地震动输入功率谱

Fig.4 Seismic input power spectra

充分利用拉索式惯容系统的安装便易性,拉索式惯容系统可采用的布置方式有:跨一层及跨两层安装(id-ib=2id-ib=3).在设计时,选取有效的变量区间,阻尼比ζd取值范围定义为ζd0,1,质量比μ取值范围定义为μ0,5,频率比rd取值范围定义为rd0,3.

不同的减震效果指标及惯容系统安装方式的布置效率如图5~7所示.

惯容系统三类布置效率分布较为相近,均为惯容单元安置于底层时最佳.综合指标及位移指标下,随着跨越层数的增加减震效率增加,在加速度指标下,跨两层布置最优减震效率略次于跨一层布置.

图5 惯容系统布置效率示意图(综合指标)

Fig.5 Layout efficiency of inerter system (comprehensive index)

图6 惯容系统布置效率示意图(位移指标)

Fig.6 Layout efficiency of inerter system (displacement index)

图7 惯容系统布置效率示意图(加速度指标)

Fig.7 Layout efficiency of inerter system (acceleration index)

依据布置效率指标,初步确定拉索式惯容系统下层安装位置均为底层,利用式(18),选取综合减震指标,进行多目标优化。由于减震效果指标与经济指标之间相互冲突,不存在同时满足两者最优的解,因此需要寻求一系列的最优解,使得任一指标在不损害另一指标的前提下达到最优,由这些最优设计点组成的曲线构成帕累托前沿,如图8所示。优化结果显示,层间布置惯容系统的最优减震效果JH2=0.462,跨一层布置惯容系统的最优减震效果JH2=0.421,跨两层布置惯容系统的最优减震效果JH2=0.411.跨一层及跨两层布置的拉索式惯容系统减震效果相近且均优于层间布置惯容系统.由帕累托前沿可知,在相同的减震效果下,层间布置水平控制力远大于跨一层布置及跨两层布置,即跨层布置相较于层间布置将能带来更多的经济效益.

图8 惯容系统的帕累托前沿

Fig.8 Pareto front for the inerter system

各减震指标下的设计参数如图9~11所示.由图9可知,惯容系统的频率比始终位于设计参数上界,说明在惯容系统设计时需考虑实际工程需要进行频率比选择.

图9 惯容系统帕累托前沿的频率比设计最优值

Fig.9 Design frequency ratio of Pareto optimal solutions

图10 惯容系统帕累托前沿的阻尼比设计最优值

Fig.10 Design damper ratio of Pareto optimal solutions

图11 惯容系统帕累托前沿的质量比设计最优值

Fig.11 Design mass ratio of Pareto optimal solutions

由图10-11所示阻尼比及质量比关系可知,对于同样的减震效果,跨一层布置及跨两层布置的惯容系统需要相对较小阻尼比和质量比,且随着质量比、阻尼比的增加减震效果显著提升.

依据平衡设计的原则,选定设计减震指标利用帕累托前沿便可获得较为经济的设计方案.以减震指标JH2=0.5为例,层间布置、跨一层布置及跨两层布置设计参数如表3所示.各层位移及加速度响应如图12-13所示.安装惯容系统后,结构响应显著降低,减震效果明显,且三种布置形式的惯容系统减震效果较为接近.

表3 惯容系统的设计参数(JH2=0.5
Table 3 Design parameters of the inerter system(JH2=0.5
LayoutFrequency ratioDamper ratioMass ratio
Inter layer 3.000 1.000 2.343
Cross 1 layer 3.000 0.506 0.634
Cross 2 layer 3.000 0.283 0.289

图12 指定减震指标下的结构位移响应

Fig.12 The structural displacement response against specified JH2

图13 指定减震指标下的结构加速度响应

Fig.13 The structural acceleration response against specified JH2

5 小结

本文研究了拉索式惯容系统的跨层布置,给出拉索式惯容系统在金井清谱随机激励下的参数优化方法,并验证了拉索式惯容系统通过便捷安装方式实现的跨层布置明显优于各类层间布置惯容系统.通过一个十层基准结构的算例分析,说明拉索式惯容系统的有效性.引入布置效率指标并分析确定惯容系统安装最佳位置为底层.通过多目标优化给出惯容系统的减震指标及经济指标的平衡设计的帕累托前沿.研究了惯容系统的优化参数分布.从最优减震效果、相同减震效果下的经济指标等多个方面论证了跨层布置的拉索式惯容系统优于层间布置惯容系统.

参考文献

1

王孝然申永军杨绍普.接地式三要素型动力吸振器的H优化.动力学与控制学报2016145):448~453 [百度学术

Wang X RShen Y JYang S PHoptimization of the grounded three-element type dynamic vibration absorber. Journal of Dynamics and Control2016145) : 448~453 ( in Chinese) [百度学术

2

郎君申永军杨绍普.半主动控制接地式动力吸振器参数优化及性能比较.动力学与控制学报2019172):168~177 [百度学术

Lang JShen Y JYang S P .Parameter optimization and performance comparison of semi-active ground-hook control DVAs. Journal of Dynamics and Control2019172):168~177 ( in Chinese) [百度学术

3

McNamara R J. Tuned mass dampers for buildings. ASCE Journal of Structural Engineering19771039):1785~1798 [百度学术

4

Kaynia A MVeneziano DBiggs J M.Seismic effectiveness of tuned mass dampers. ASCE Journal of the Structural Division19811078):1465~1484 [百度学术

5

Chen M Z QPapageorgiou CScheibe Fet al. The missing mechanical circuit element. IEEE Circuits and Systems Magazine200991):10~26 [百度学术

6

Ikago KSaito KInoue N.Seismic control of single-degree-of-freedom structure using tuned viscous mass damper. Earthquake Engineering & Structural Dynamics2012413):453~474 [百度学术

7

Ikago KSugimura YSaito Ket al.Modal response characteristics of a multiple-degree-of-freedom structure incorporated with tuned viscous mass dampers. Journal of Asian Architecture and Building Engineering2018112):375~382 [百度学术

8

Marian LGiaralis A.Optimal design of a novel tuned mass-damper-inerter (TMDI) passive vibration control configuration for stochastically support-excited structural systems.Probabilistic Engineering Mechanics201438156~164 [百度学术

9

Pan CZhang R F.Design of structure with inerter system based on stochastic response mitigation ratio. Structural Control and Health Monitoring2018256):e2169 [百度学术

10

Taflanidis A AGiaralis APatsialis D. Multi-objective optimal design of inerter-based vibration absorbers for earthquake protection of multi-storey building structures. Journal of the Franklin Institute201935614):7754~7784 [百度学术

11

Arakaki TKuroda HArima Fet al. Development of seismic devices applied to ball screw. Part 1Basic performance test of RD-series. AIJ Journal of Technology and Design199958):239~244 [百度学术

12

Pan CZhang R FLuo Het al. Demand-based optimal design of oscillator with parallel-layout viscous inerter damper. Structural Control and Health Monitoring2018251):e2051 [百度学术

13

Kanai K.Semi-empirical formula for the seismic characteristics of the ground. Bulletin of the Earthquake Research Institute195735309~325 [百度学术

14

Lutes L DSarkani S. Stochastic analysis of structural and mechanical vibrations. New JerseyPrentice Hall1997 [百度学术

微信公众号二维码

手机版网站二维码