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

确定继续浏览么?

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

非平稳地震激励下惯容耗能结构动力响应解析法  PDF

  • 李创第 1
  • 王瑞勃 1
  • 葛新广 2
  • 江丽富 1
1. 广西科技大学土木建筑工程学院,广西 柳州 545006; 2. 柳州工学院土木建筑工程学院,广西 柳州 545616

中图分类号: TU311.3O324

最近更新:2024-12-03

DOI:10.16385/j.cnki.issn.1004-4523.2024.11.007

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

摘要

针对惯容减震系统在非平稳地震激励下动力响应研究不足的问题,提出设有混联Ⅰ型惯容减震系统的多自由度耗能结构动力响应时变方差的解析解法。根据惯容减震系统的本构关系、结构的动力方程及非平稳地震激励,综合利用复模态法和虚拟激励法,将惯容耗能结构解耦为一阶系统,以便获得结构位移、速度、层间剪力等响应量的统一解。采用二次式分解法将统一解的时变功率谱密度函数转化为惯容耗能结构复模态特征值、模态系数、时变模态强度系数和含有圆频率平方项的二次式乘积的线性组合。在此基础上,利用非平稳模态谱矩在无限积分区间有解析解的特征,推导出非平稳地震激励下耗能结构响应时变方差的解析解。通过采用突加型白噪声激励对结构动力响应进行分析,验证了所提动力响应功率谱和时变方差的正确性。同时,基于突加型Kanai‑Tajimi模型的框架结构动力响应研究,分析了惯容系统参数对减震效果的影响。所提方法可适用于线性结构在其他非平稳调制函数的地震动下的响应分析。

引 言

惯容器作为实用有效的被动控制装置,在结构中可采用与传统消能装置相同的安装方式。在外部激励作用下,惯容器通过两端点之间的加速度差值产生作用力,从而降低结构的动力响

1⁃3。IKAGO4提出了调谐黏性质量阻尼器(Tuned Viscous Mass Damper,TVMD),并验证了单自由度结构中TVMD在谐波地震激励下抗震的有效性。LAZAR5将惯容⁃弹簧⁃阻尼3个元件相互连接构成调谐惯容阻尼器(Tuned Inerter Damper,TID),经研究表明TID在地震激励作用下能够达到极好的减震效果。LI6探索了调谐黏性质量阻尼器(TVMD)、调谐惯容阻尼器(TID)和调谐质量阻尼器(Tuned Mass Damper,TMD)3种控制系统在白噪声激励下对结构响应的影响,同时基于不同控制目标对3种控制系统进行了优化设计。DE DOMENICO7以双过滤白噪声激励作为地震加速度,研究了基础隔震结构中调谐惯容质量阻尼器(Tuned Mass Damper Inerter,TMDI)的抗震性能,研究表明TMDI的抗震效果比TMD更好,同时能够避免长周期地震作用下可能发生的共振行为。为进一步加强惯容器的振动控制效果,国内外学者将惯容元件、耗能元件和弹簧元件连接协同工作构成惯容减震系统,其中经典的惯容减震系统有:串联型惯容减震系8、混联Ⅰ型惯容减震系统(SPIS⁃Ⅰ9和混联Ⅱ型惯容减震系统(SPIS⁃Ⅱ10。文献[11⁃12]对3种惯容减震系统在平稳随机地震激励作用下的耗能增效现象进行了研究,阐明了耗能增效理论并推导了与之相关的解析表达式。为促进惯容减震系统在工程中的应用,文献[13⁃14]对其振动控制性能和参数设15也进行了研究,研究结果表明,惯容减震系统可以显著降低结构的动力响应,在结构的振动控制中具有良好的应用前景。

上述研究基于平稳地震激励(如白噪声激励、Kanai⁃Tajimi激励、双过滤白噪声激励等),虽简化了分析过程,但平稳地震激励与实际地震过程存在一定差异。实际地震过程具有明显的非平稳特性,故拟合地震的随机模型应呈非平稳状态。目前已经提出了多种拟合非平稳随机过程的功率谱模

16⁃18,其中,PRIESTLEY18提出的演变功率谱模型,用“频率”和“能量”的概念解释了演变随机过程的时变均方能量分布,该模型在地震工程的动力响应分析中尤为重要。演变功率谱模型可分为均匀调制非平19⁃21和完全非平22两类。其中,均匀调制非平稳的谱密度函数等于与时间t相关的调制函数和与圆频率ω相关的平稳地震激励谱密度函数的乘积。现有的调制函数主要包括阶跃型调制函21、余弦型调制函21⁃23、Shinozuka⁃Sato调制函22⁃23及Iyengar型调制函23等。时域法和频域法是求解非平稳随机结构动力响应的主要方法。方同2022利用复模态法和杜哈梅积分,求解了响应及激励的协方差,分析了结构基于非平稳激励的动力响应。由于现存的地震激励模型以功率谱形式存在,采用时域法分析结构随机响应时,需对地震激励的功率谱密度函数进行转换,增加了计算量。文献[24⁃26]从频域角度提出虚拟激励法,通过构造虚拟激励,将求解非平稳响应功率谱密度函数的问题转化为确定性外荷载作用下的瞬态响应分析,简化了计算过程。响应的功率谱密度函数和时变方差是非平稳随机响应分析的重要参数。虚拟激励法计算响应时变方差是在已知响应功率谱的前提下,对其在整个频率范围内进行积分,为获得准确的结果需取较大的积分上限值和非常小的积分步长,导致计算效率较低。葛新广27在研究平稳地震过程中提出功率谱二次式分解法,该方法能避开平稳激励下响应谱矩和方差的复杂计算,但并未研究非平稳激励下耗能结构的动力特性。

本文为研究设有惯容减震系统的耗能结构在非平稳地震作用下的动力特性,同时为提高时变方差的计算效率和计算精度,以设有混联Ⅰ型惯容减震系统的耗能结构为例,对均匀调制非平稳随机地震激励下的结构响应做完整的推导,并获得功率谱及时变方差的解析解。应用多自由度数值算例及框架结构算例,证明了本文所求动力响应解析解的有效性,对参数及减震效果也进行了分析。

1 动力方程

惯容减震系

8⁃9(简称惯容系统)由阻尼、弹簧、惯容3种力学元件连接组成,将其设置于多自由度结构中构成惯容耗能结构,该结构在地震加速度x¨g作用下的动力方程为:

Mx¨+Cx˙+Kx+I1Fs=-MΙx¨g (1)

式中 MKC分别表示结构的质量矩阵、抗侧刚度矩阵及阻尼矩阵;M=diagmi,其中mi为结构第i层的质量(i=1,2,,n)Kn×n=K1+K2-K200-K2K2+K30000-KnKn,其中Kn为结构第n层抗侧刚度;C依据Rayleigh阻尼模型计算,C=aM+bK,其中ab为Rayleigh系数;xx˙x¨分别表示结构的位移向量、速度向量及加速度向量;I1FsΙ分别表示惯容系统的影响向量、惯容力向量及常数向量。

混联Ⅰ型惯容系统设置于结构第i层的构造简图如图1所示。

图1  混联Ⅰ型惯容系统第i层的构造简图

Fig.1  Schematic layout of SPIS-Ⅰ in the ith storey

式(1)中的惯容力为:

Fs,1=kd1x˜1,1+cd1x˜˙1,1=md1x˜¨1,2,  x˜1,1+x˜1,2=x1 (2)
Fs,i=kdix˜i,1+cdix˜˙i,1=mdix˜¨i,2,  x˜i,1+x˜i,2=xi-xi-1 (3)

式中 Fs,1Fs,i分别为首层和第i层的惯容力;x˜i,1x˜i,2分别为惯容系统中弹簧元件和惯容元件两端的相对位移,二者之和等于结构的层间位移。kdicdimdi分别表示结构第i层惯容系统的刚度系数、阻尼系数及惯容系数。

式(2)和(3)用矩阵形式表示为:

Fs=Kdx˜1+Cdx˜˙1=Mdx˜¨2 (4)
x˜1+x˜2=I1Tx (5)

式中 CdKdMd分别表示惯容系统的阻尼系数矩阵、刚度系数矩阵及惯容系数矩阵;x˜1为弹簧元件的位移向量,x˜1=x˜1,1  ...  x˜i,1  ...  x˜n,1Tx˜2为惯容元件的位移向量,x˜2=x˜1,2  ...  x˜i,2  ...   x˜n,2T

根据式(4)和(5),惯容系统的各参数矩阵与结构加速度矩阵x¨之间的关系为:

Mdx˜¨1+Kdx˜1+Cdx˜˙1=MdI1Tx¨ (6)

式(6)代入式(1),动力方程为:

Mx¨+Cx˙+Kx+I1(Kdx˜1+Cdx˜˙1)=-MΙx¨g (7)

2 非平稳结构响应统一解

2.1 复模态分析

混联Ⅰ型惯容减震系统被设置于多自由度结构中,构成非对称、非经典阻尼系统,需运用复模态

202228式(7)解耦,现引入状态变量y

y=xx˙x˜1x˜˙1T (8)

联立式(6),(7)和 (8),状态方程为:

M0y˙+K0y=αx¨g (9)

式中 M0K0为4n×4n阶矩阵,其值分别为:

M0=CMI1Cdo1o1-MdI1To1MdMo1o1o1o1o1Mdo1
K0=Ko1I1Kdo1o1o1KdCdo1-Mo1o1o1o1o1-Md

α=-MI  o2  o2  o2TM0K0α中包含的o1n×n阶零矩阵,o2n×1阶零矩阵。

式(9)中存在右模态矩阵U、左模态矩阵V及特征值矩阵p,三者具有以下关系:

p=-VTM0U-1VTK0U (10)

复模态变换:

y=Uz (11)

式中 z为模态变量。

联立式(6),(10)和(11),并在方程两边乘以VT可得:

z˙-pz=ηx¨g (12)

式中 η为模态系数矩阵,η=VTα/(VTM0U)

式(12)的分量形式为:

z˙k-pkzk=ηkx¨g;  k=1,2,,4n (13)

式中 pk为振动特征值,其值为复数,实部为负。

2.2 均匀调制非平稳地震激励下的结构响应

工程上具有多种非平稳地震激励,PRIESTLEY

18提出的演变功率谱密度模型用“频率”与“能量”的概念解释了演变随机过程的时变均方能量分布,其谱密度函数Φx¨g(ω,t)可表示为:

Φx¨g(ω,t)=a(ω,t)a*(ω,t)Sx¨g(ω)=a(ω,t)2Sx¨g(ω) (14)

式中 Sx¨g(ω)为平稳随机过程的谱密度函数;a(ω,t)为调制函数,a*(ω,t)a(ω,t)的共轭项。

式(14)中调制函数a(t)取单位阶跃函数:

a(t)=1,  t00,  t<0 (15)

根据虚拟激励

24⁃26x¨g(ω,t)可表示为:

x¨g(ω,t)=Φx¨g(ω,t)eiωt (16)

式(16)代入式(13),求解非齐次线性方程,获得模态变量zk

zk(ω,t)=ηkSx¨g(ω)0tepk(t-τ)a(τ)eiωτdτ (17)

化简式(17)得:

zk(ω,t)=Ik(ω,t)Sx¨g(ω) (18)

由式(15),(17)及(18)得Ik(ω,t)为:

Ik(ω,t)=ηk0tepk(t-τ)a(τ)eiωτdτ=ηkeiωt-epktiω-pk (19)

联立式(8),(9)及(18),获得结构各层相对地面的位移及速度分别为:

xi(ω,t)=k=14nui,kzk(ω,t) ; i=1,2,,n (20a)
x˙i(ω,t)=k=14nun+i,kzk(ω,t) ; i=1,2,,n (20b)

式中 ui,k表示Ui行第k列的元素。

层间位移及层间速度分别为:

Δxi(ω,t)=k=14nui,k-ui-1,kzk(ω,t) ; i=2,3,,n (21a)
Δx˙i(ω,t)=k=14n(un+i,k-un+i-1,k)zk(ω,t) ; i=2,3,,n (21b)

式中 Δx1=x1Δx˙1=x˙1,其中x1x˙1分别为结构第1层的层间位移和层间速度。

层间剪力为:

V1(ω,t)=K1k=14nu1,kzk(ω,t) ;  i=1 (22a)
Vi(ω,t)=Kik=14nui,k-ui-1,kzk(ω,t) ; i=2,3,,n (22b)

层间位移角为:

θ1(ω,t)=1h1k=14nu1,kzk(ω,t) ;  i=1 (23a)
θi(ω,t)=1hik=14nui,k-ui-1,kzk(ω,t) ; i=2,3,,n (23b)

式中 h表示层高。

由式(20)~(23)可知,在非平稳随机激励下耗能结构的动力响应(位移、速度、层间位移、层间速度、层间剪力及层间位移角)表达式相似,可统一表示为:

X¯l(ω,t)=k=14nκl,kzk(ω,t) ; l=1,2,,6 (24)

式中 X¯l(ω,t)为响应量;κl,k为各响应分量对应的模态系数:κ1,k为位移模态系数,κ2,k为速度模态系数,κ3,k为层间位移模态系数,κ4,k为层间速度模态系数,κ5,k为层间剪力模态系数,κ6,k为层间位移角模态系数。

3 结构动力响应功率谱

由随机振动理

22可知,结构动力响应的自功率谱密度函数为:

SX¯lX¯l(ω,t)=k=14nj=14nκl,kκl,jzk(ω,t)zj*(ω,t) (25)

式中 zj*(ω,t)zk(ω,t)互为共轭。

式(18)代入式(25)

SX¯lX¯l(ω,t)=k=14nj=14nκl,kκl,jAjk(ω,t)Sx¨g(ω) (26)

其中:

Ajk(ω,t)=Ik(ω,t)Ij*(ω,t) (27)

存在恒等式:

l=1Nj=1Ndlbj=l=1Ndlbl+l=1N-1j=l+1N(dlbj+bldj) (28)

依据式(28)式(26)分解为:

SX¯lX¯l(ω,t)= Sx¨g(ω)×j=14nκj2Ajj(ω,t)+j=14n-1k=j+14nκjκkAjk(ω,t)+Akj(ω,t) (29)

依据式(27)和(19)对参数Ajj(ω,t)Ajk(ω,t)+Akj(ω,t)进行求解,将其转化为包含复模态特征值和圆频率平方项的二次式,具体过程见附录A,结果为:

Ajj(ω,t)=ηj2a1j(t)+a2j(t)cos(ωt)pj2+ω2 (30)
Ajk(ω,t)+Akj(ω,t)=2ηjηkpk+pja5jkωsin(ωt)pk2+ω2-ωsin(ωt)pj2+ω2+a3jk+a4jkcos(ωt)pkpk2+ω2+pjpj2+ω2 (31)

式中 a1j(t)=1+e2pjta2j(t)=-2epjta3jk=1+e(pk+pj)ta4jk=-epkt-epjta5jk=epkt-epjt

参数a1j(t)a2j(t)a3jk(t)a4jk(t)a5jk(t)仅与复振动特征值pj和时间t有关,与平稳地震激励参数无关,故将以上参数称之为时变模态强度系数。将时变模态强度系数、复振动特征值pj及模态系数ηj等代入式(29),即可得到动力响应功率谱密度函数的二次式分解。

4 结构动力响应时变方差

由随机振动理

22可知,动力响应时变方差MX¯l,0与响应0阶谱矩相同,动力响应变化率的时变方差MX¯˙l,0与2阶谱矩结果一致,耗能结构动力响应的时变方差表达式为:

MX¯l,0(t)=-SX¯X¯(ω,t)dω (32a)
MX¯˙l,0(t)=-SX¯˙X¯˙ω,tdω (32b)

式中 MX¯l,0为动力响应的时变方差(0阶谱矩);MX¯˙l,0为动力响应变化率的时变方差(2阶谱矩)。

式(29)代入式(32a)得:

MX¯l,0=2j=14nκj2ηj20Sx¨g(ω)Ajj(ω,t)dω+2j=14n-1k=j+14nκjκk0Sx¨g(ω)Ajk(ω,t)+Akj(ω,t)dω (33)

式(33)进行简

29

MX¯l,0=2j=14nκj2ηj2(a1jβ1j+a2jβ2j)+4j=14n-1k=j+14nκjκkηjηkpk+pja3jkpjβ1j+pkβ1k+a4jkpjβ2j+pkβ2k+a5jkβ3k-β3j (34)

其中:

β1j(t)=0Sx¨g(ω)1pj2+ω2dω (35a)
β2j(t)=0Sx¨g(ω)cos(ωt)pj2+ω2dω (35b)
β3j(t)=0Sx¨g(ω)ωsin(ωt)pj2+ω2dω (35c)

由式(35a)~(35c)可知,β1j(t)β2j(t)β3j(t)与时间t、复振动特征值pj和平稳随机过程的谱密度函数Sx¨g(ω)有关,本文将其称为非平稳模态谱矩。

5 算例验证

5.1 突加白噪声激励下惯容耗能结构动力分析

某6层结构中设有混联Ⅰ型惯容系统,如图2所示。主体结构的质量m1 = m2 = 4.25×105 kg,3~6层质量m3 ~ m6 = 3.00×105 kg;抗侧刚度k1 = k2 = 2.58×107 N/m,3~6层刚度k3 ~ k6 = 1.45×107 N/m;阻尼依据Rayleigh阻尼计算即C=aM+bK,其中a= 0.13885,b=0.014298,表示Rayleigh阻尼常数。混联Ⅰ型惯容系统中的参数分别有惯容系数md=μdm1、刚度系数kd=κdK1和阻尼系数cd=2ξdmdkd/mdm1K1为结构第一层的质量和刚度;μdκdξd分别表示惯质比、刚度比和名义阻尼比)。

图2  混联Ⅰ型惯容耗能结构简图

Fig.2  Schematic diagram of the energy-dissipation structure with SPIS-Ⅰ

假定μd=0.04κd=0.1927ξd=0.043以验证上述所求响应功率及时变方差的正确性,即混联Ⅰ型惯容系统的惯容系数 md = 1.7×104 kg,刚度系数kd=4.971×106 N/m,阻尼系数cd=2.523×104 N·s/m。均匀调制非平稳地震激励模型的谱密度函数等于有关时间t的调制函数与有关圆频率的平稳地震激励谱密度函数的乘积,如式(14)所示。选取调制函数为单位阶跃型,平稳地震激励为白噪声模型进行分析。白噪声模型是最早模拟地面运动的随机地震激励,其谱密度函数在整个频域范围内为常值S0=1.2×10-4 m2·s-3,能够极大地简化计算过程,即谱密度函数为:

Sx¨g(ω)=S0 (36)

式(36)代入式(29),即可得到耗能结构基于突加型白噪声激励的位移功率谱Sx(ω,t)、层间位移功率谱SΔx(ω,t)和速度功率谱Sx˙(ω,t)。各响应功率谱密度函数与结构的自振频率ω、振动复特征值pj、各响应对应的模态系数κl,k及参数Ajj(ω,t)等有关。为验证本文所求上述响应功率谱的正确性,将本文方法与传统方法(虚拟激励法)依据式(25)所得结果进行比较,结构顶层的各响应功率谱如图3所示。由图3可知,两种方法所求位移、层间位移和速度的功率谱是关于频率ω及时间t的函数,两种结果在频率ω和时间t范围内的变化趋势完全相同,在时间t > 20 s后逐渐减小随后趋于稳定,符合突加型非平稳地震激励的特性,由此表明本文方法所求响应功率谱正确有效。

图3  基于突加型白噪声激励的响应功率谱

Fig.3  Response power spectrum based on sudden-additive white noise excitation

式(34)所示,惯容耗能结构在突加型白噪声激励下的响应时变方差及速度时变方差与非平稳模态谱矩β1,jβ2,jβ3,j密切相关,且β1,jβ2,jβ3,j在积分运算时原函数存在,所得结果如下:

β1,j(t)=S001pj2+ω2dω=πS02pj (37a)
β2,j(t)=S00cos(ωt)pj2+ω2dω=πe-pjt2pjS0 (37b)
β3,j(t)=0ωsin(ωt)pj2+ω2dω=πe-pjt2S0 (37c)

将式(37a)~(37c)代入式(34),即可获得在突加型白噪声激励下结构动力响应时变方差的解析解。为验证所求结果的有效性,借助MATLAB工具将本文方法与传统分析方法(虚拟激励法)得到的结果进行对比。传统方法依据式(32a)~(32b)对响应时变方差进行求解,即对响应功率谱在频率域ω中取(0,∞)进行积分,为保证计算结果的准确性,本文取积分上限为500 rad/s,划分区间的积分步长Δω依次取2,1,0.5和0.1 rad/s。对比本文方法和传统方法(4种不同工况)所求结构顶层的各响应时变方差的吻合程度,如图4~9所示。

图4  位移时变方差

Fig.4  Time-varying variance of displacement

图5  层间位移时变方差

Fig.5  Time-varying variance of interlayer displacement

图6  速度时变方差

Fig.6  Time-varying variance of velocity

图7  层间速度时变方差

Fig.7  Time-varying variance of interlayer velocity

图8  层间剪力时变方差

Fig.8  Time-varying variance of interlayer shear force

图9  层间位移角时变方差

Fig.9  Time-varying variance of interlayer deflection angle

由图4~9可知,在固定积分上限的情形下,积分步长的大小影响了传统方法计算结果的准确性和精度。传统方法在Δω=2 rad/s工况下,获得的结构顶层的位移、层间位移、速度、层间速度、层间剪力和层间位移角的时变方差准确度不高,且两种方法存在较大偏差。传统方法在Δω=1 rad/s工况下,上述响应时变方差的准确性及计算精度有一定提升,两种方法之间的偏差相对减小。在Δω=0.5 rad/s工况下,传统方法与本文方法结果较为接近,且随时间t的增大偏差逐渐增大,两种结果不完全相同。在Δω=0.1 rad/s工况下,频率域ω被5000等分,时间域t被80等分,传统方法与本文方法结果完全一致。经分析可知,传统方法所求结果的计算精度与积分步长Δω的取值密切相关,Δω取值越小,越逼近本文结果,由此反映了本文方法所求各响应的时变方差正确、有效且精度较高。由于传统方法及本文方法都是在响应功率谱的基础上对时变方差进行求解,且两种方法所得结果最终相同,再次证明了本文推导的响应功率谱表达式合理。

5.2 框架结构的动力响应分析

为进一步研究本文方法的实用性,采用与文献[

30]相似的10层框架结构进行分析,如图10所示。结构中柱子尺寸为600 mm×600 mm,梁横截面尺寸为300 mm×700 mm,楼板厚10 mm。结构跨度6 m,层高为3 m,梁、板、柱均采用C30混凝土,材料弹性模量E= 3.0×1010 N/m2,假设楼板为刚性,结构各层间的抗侧刚度等于各柱的抗侧刚度之和,单个柱的抗侧刚度为12EI/h3。结构竖向荷载主要由结构自重、楼板的铺装自重(1.0 kN/m2)、楼板自重(2.5 kN/m2)及竖向活荷载(2.0 kN/m2)构成。各节点质量依据集中质量法进行计算,荷载从属宽度为6.0 m,活荷载系数为0.5,恒荷载系数为1.0,阻尼比为0.05。将混联Ⅰ型惯容系统设置于多层框架结构中,各层设置惯容系统的型号相同。现对混联Ⅰ型惯容耗能框架结构的动力响应进行研究。

图10  耗能框架结构

Fig.10  Energy-dissipation frame structure

为分析框架结构基于均匀调制非平稳函数的动力特性,取调制函数为阶跃性调制函数,平稳地震激励选用考虑了基岩土层的过滤作用且具有明确频谱特性的Kanai⁃Tajimi平稳随机激励。依据文献[

31]中提出的Kanai⁃Tajimi模型,参数取值分别为S0= 5.75 cm2/s4/Hz,ωg2 = 242.0 s-2ξg= 0.63988,Kanai⁃Tajimi模型的谱密度函数为:

Sx¨g(ω)=S0ωg4+4ξg2ωg2ω2(ωg2-ω2)2+4ξg2ωg2ω2 (38)

获得突加型Kanai⁃Tajimi激励下响应时变方差的前提是已知响应功率谱,响应功率谱与激励功率谱之间存在定关系,2.2节对调制函数进行了处理,现对式(38)进行简化,结果如下式所示:

Sx¨g(ω)=S0r=12cg,rω2+ωg,r2 (39)

式中 ωg,12=-1-2ξg2ωg2+2iωg2ξgξg2-1ωg,22=-1-2ξg2ωg2-2iωg2ξgξg2-1cg,1=ωg4-4ξg2ωg2ωg,12ωg,22-ωg,12cg,2=ωg4-4ξg2ωg2ωg,22ωg,12-ωg,22

5.2.1 响应功率谱正确性分析

响应功率谱可显示耗能结构响应随频率及时间的振动特性,现对框架结构的动力响应功率谱进行分析。采用与5.1节相同的惯质比、刚度比和阻尼比验证本节框架结构响应功率谱的正确性。将式(39)代入式(29),即可获得混联Ⅰ型耗能框架结构基于突加型Kanai⁃Tajimi激励的层间剪力功率谱SV(ω,t)、层间位移角功率谱Sθ(ω,t)、位移功率谱Sx(ω,t)、层间位移功率谱SΔx(ω,t)及速度功率谱Sx˙(ω,t)图11为两种方法所求结构顶层的各响应功率谱的对比。由图11可知,本文方法与传统方法所求Sx(ω,t)SΔx(ω,t)Sx˙(ω,t)SV(ω,t)Sθ(ω,t)同样受时间t与频率ω的共同影响,且两种方法所求各响应的功率谱峰值与变化趋势完全一致,由此反映了所求结果的正确性,同时也表明本文方法同样适用于突加型Kanai⁃Tajimi激励的动力响应分析。其次,由图11可知,响应功率谱Sx(ω,t)SΔx(ω,t)Sx˙(ω,t)SV(ω,t)Sθ(ω,t)在频率ω = 3 rad/s及时间t = 20 s左右出现最大值,随后受到突加型地震激励的影响趋于平稳。

图11  基于突加型Kanai-Tajimi激励的响应功率谱

Fig.11  Response power spectrum based on sudden-additive Kanai-Tajimi excitation

5.2.2 响应时变方差分析

在已知响应功率谱的前提下,依据式(35a)~(35c)求解了基于突加型Kanai⁃Tajimi激励的β1,jβ2,jβ3,j,计算过程见附录B,结果为:

β1,j(t)=πS02r=12cg,rωg,r2-pj21pj-1ωg,r (40a)
β2,j(t)=πS02r=12cg,rωg,r2-pj2e-pjtpj-e-ωg,rtωg,r (40b)
β3,j(t)=πS02r=12cg,rωg,r2-pj2e-pjt-e-ωg,rt (40c)

将式(40a)~(40c)代入式(34),即可获得基于突加型Kanai⁃Tajimi地震激励的响应时变方差解析解。图12~17分别为结构顶层各响应时变方差基于本文方法与虚拟激励法(取不同积分步长Δω)的对比图。由图12~17可知,随积分步长Δω的减小,传统方法(虚拟激励法)所获得的位移、层间位移、速度、层间速度、层间剪力和层间位移角的时变方差越逼近本文结果。传统方法在Δω ≤ 0.5 rad/s工况下所求结果与本文所求结果完全相同,由此表明了本文所求各响应时变方差的正确性。此外,借助MATLAB工具,传统方法依据式(32a)~(32b)计算响应时变方差,在工况Δω=2 rad/s,Δω=1 rad/s,Δω=0.5 rad/s和Δω=0.1 rad/s下分别耗时150.23,287.845,497.137,537.92 s,而本文方法耗时54.01 s,由此反映了本文方法在保证计算精度的同时也提高了计算效率。

图12  位移时变方差

Fig.12  Time-varying variance of displacement

图13  层间位移时变方差

Fig.13  Time-varying variance of interlayer displacement

图14  速度时变方差

Fig.14  Time-varying variance of velocity

图15  层间速度时变方差

Fig.15  Time-varying variance of interlayer velocity

图16  层间剪力时变方差

Fig.16  Time-varying variance of interlayer shear force

图17  层间位移角时变方差

Fig.17  Time-varying variance of interlayer displacement angle

5.2.3 参数及减震效果分析

上文验证了本文方法所求框架结构在突加型Kanai⁃Tajimi激励下的响应方差正确有效,为研究惯容系统主要参数对减震效果的影响,基于式(40a)~(40c)及(34),通过改变惯质比μd、刚度比κd及阻尼比ξd分析位移方差及速度方差的变化趋势。μd的取值范围为0.001~0.4,κd取值范围为0.001~0.06,ξd的取值范围为0.01~1。

图18展示了惯质比、刚度比和阻尼比在t=20 s时对结构顶层的位移方差和速度方差的影响。由图18(a)和(b)可知,阻尼比ξd=0.2时,位移方差和速度方差随惯质比的增大而逐渐减小,位移方差和速度方差随刚度比的增大先减小后增大,且在κd=0.02左右出现极小值。由图18(c)和(d)可知,惯质比μd=0.4时,位移方差和速度方差随阻尼比的增大逐渐增大,随刚度比的增大先减小后增大。分析表明,位移方差和速度方差在其中两个参数的共同作用下出现极小值,此时减震效果较好。

图18  对结构响应方差的影响

Fig.18  Influence of the parameters on the structural response variance

图18改变参数使位移方差和速度方差同时出现极小值,此时惯质比μd、刚度比κd及阻尼比ξd分别为0.4,0.023及0.2,依据式(40a)~(40c)和(34),对比该情形下耗能结构与无控结构的响应量(位移、速度、层间位移和层间剪力)如图19所示。由图19可知,在结构中设置混联Ⅰ型惯容系统可使结构顶层的位移、速度、层间位移及层间剪力显著减小,各响应最大值分别减小26.71%,31.15%,26.09%及31.16%。

图19  减震效果分析

Fig.19  Analysis of damping effects

6 结 论

本文研究了混联Ⅰ型惯容耗能结构在均匀调制非平稳随机地震激励作用下的动力响应,为进一步了解结构的振动特性,推导了结构响应时变功率谱及时变方差的解析解。基于解析解的表达式,借助与动力学相关的数值算例及框架结构的算例,验证了本文结果的正确性,分析了参数对减震效果的影响。所得结论如下:

(1) 将非平稳地震激励下的响应时变功率谱密度函数二次式分解转化为由振动特征值,模态系数及时变模态强度系数等参数构成的线性组合,给定模态响应系数可得到与之对应的响应时变功率谱,具有明确的物理意义。

(2) 经本文方法的运算,将求解响应时变方差的积分问题转化为求解非平稳模态谱矩βlj问题。非平稳模态谱矩βlj的积分运算存在原函数,进而可获得响应时变方差的解析解,简化了计算过程。响应的时变方差可应用于分析结构在地震激励下的可靠度、安全性及舒适性。

(3) 通过动力学数值算例与框架结构算例分析表明,本文与虚拟激励法所求响应时变功率谱完全一致,验证了本文方法的准确性。本文所求响应时变方差与虚拟激励法最精确的值相同,且计算效率较高,即在提高精度的同时也提高了计算效率。

(4) 通过参数对位移时变方差、速度时变方差和减震效果的影响分析表明,改变惯质比、刚度比及阻尼比,动力响应方差出现极小值,此时混联Ⅰ型惯容系统具有良好的减震性能。

附录A 时变模态强度系数及Akj(ω,t)Ajj(ω,t)的推导

依据式(27)及(19)推导Akj(ω,t)Ajj(ω,t),其中Ajj(ω,t)的表达式为:

Ajj(ω,t)=1-epjt+iωt-epjt-iωt+e2pjtηj2pj2+ω2 (A1)

欧拉公式:

e-iωt=cos(ωt)-isin(ωt),    eiωt=cos(ωt)+isin(ωt) (A2)

运用欧拉公式可将式(A1)简化为:

Ajj(ω,t)=ηj2a1j(t)+a2j(t)cos(ωt)pj2+ω2 (A3)

同理Akj(ω,t)为:

Akj(ω,t)=ηkηj1-epkt+iωt-epjt-iωt+e(pk+pj)t1iω-pj1-iω-pk=-ηkηjpk+pj1pj2+ω21-e(pj-iω)t-e(pk+iω)t+e(pk+pj)t-pj-iω-ηjηkpk+pj1pk2+ω2[1-e(pj-iω)t-e(pk+iω)t+e(pk+pj)t]-pk+iω (A4)

Ajk(ω,t)Akj(ω,t)推导过程相同,故运用欧拉公式推导Ajk(ω,t)+Akj(ω,t)

Ajk(ω,t)+Akj(ω,t)=-ηjηkpk+pj1pk2+ω2[pk(e(pk-iω)t+e(pk+iω)t+e(pj+iω)t+e(pj-iω)t-2-2e(pk+pj)t)+iω(e(pk-iω)t-e(pk+iω)t+e(pj+iω)t-e(pj-iω)t)]-ηjηkpk+pj1pj2+ω2[pj(e(pk-iω)t+e(pk+iω)t+e(pj+iω)t+e(pj-iω)t-2-2e(pk+pj)t)+iω(e(pk+iω)t-e(pk-iω)t+e(pj-iω)t-e(pj+iω)t)]=ηjηkpk+pj2pk2+ω2[pk(1+e(pk+pj)t)-pk(epkt+epjt)cos(ωt)-epkt-epjtωsin(ωt)]+ηjηkpk+pj2pj2+ω2[pj(1+e(pk+pj)t)-pj(epkt+epjt)cos(ωt)-epjt-epktωsin(ωt)]=2ηjηkpk+pj×a3jk+a4jkcos(ωt)pkpk2+ω2+pjpj2+ω2+a5jkωsin(ωt)pk2+ω2-ωsin(ωt)pj2+ω2 (A5)
β1,j(t)=S0r=120cg,rω2+ωg,r21pj2+ω2dω=S0r=12cg,rωg,r2-pj201pj2+ω2-1ω2+ωg,r2dω=πS02r=12cg,rωg,r2-pj21pj-1ωg,r (B1)
β2,j(t)=S0r=120cg,rω2+ωg,r2cos(ωt)pj2+ω2dω=S0r=12cg,rωg,r2-pj20cos(ωt)pj2+ω2-cos(ωt)ω2+ωg,r2dω=πS02r=12cg,rωg,r2-pj2e-pjtpj-e-ωg,rtωg,r (B2)
β3,j(t)=S0r=120cg,rω2+ωg,r2sin(ωt)pj2+ω2dω=S0r=12cg,rωg,r2-pj20sin(ωt)pj2+ω2-sin(ωt)ω2+ωg,r2dω=πS02r=12cg,rωg,r2-pj2e-pjt-e-ωg,rt (B3)

参考文献

1

李壮壮申永军杨绍普. 基于惯容-弹簧-阻尼的结构减振研究[J]. 振动工程学报2018316): 1061-1067. [百度学术] 

LI ZhuangzhuangSHEN YongjunYANG Shaopuet al. Study on vibration mitigation based on inerter-spring-damping structure[J]. Journal of Vibration Engineering2018316): 1061-1067. [百度学术] 

2

路畅李春祥曹黎媛. 基于拓扑布置的超高层建筑TTMDI风致振动控制[J]. 振动与冲击2022419): 244-252. [百度学术] 

LU ChangLI ChunxiangCAO Liyuan. Wind-induced vibration control of super high-rise building TTMDI system based on topological layout[J]. Journal of Vibration and Shock2022419): 244-252. [百度学术] 

3

LI Z CXU KMA R Set al. A novel lever-arm tuned mass damper inerter (LTMDI) for vibration control of long-span bridges[J]. Engineering Structures2023293116662. [百度学术] 

4

IKAGO KSAITO KINOUE N. Seismic control of single-degree-of-freedom structure using tuned viscous mass damper[J]. Earthquake Engineering and Structural Dynamics2012413): 453-474. [百度学术] 

5

LAZAR I FNEILD S AWAGG D J. Using an inerter‐based device for structural vibration suppression[J]. Earthquake Engineering and Structural Dynamics2014438): 1129-1147. [百度学术] 

6

LI Y FLI S YCHEN Z Q. Optimal design and effectiveness evaluation for inerter-based devices on mitigating seismic responses of base isolated structures[J]. Earthquake Engineering and Engineering Vibration2021204): 1021-1032. [百度学术] 

7

DE DOMENICO DRICCIARDI G. Optimal design and seismic performance of tuned mass damper inerter (TMDI) for structures with nonlinear base isolation systems[J]. Earthquake Engineering and Structural Dynamics20184712): 2539-2560. [百度学术] 

8

张瑞甫曹嫣如潘超. 惯容减震(振)系统及其研究进展[J]. 工程力学20193610): 8-27. [百度学术] 

ZHANG RuifuCAO YanruPAN Chao. Inerter system and its state-of-the-art[J]. Engineering Mechanics20193610): 8-27. [百度学术] 

9

KRENK SHØGSBERG J. Tuned resonant mass or inerter-based absorbers: unified calibration with quasi-dynamic flexibility and inertia correction[J]. Proceedings of the Royal Society A20164722185): 20150718. [百度学术] 

10

李创第江丽富王瑞勃. 单自由度混联Ⅱ型惯容系统随机地震动响应分析[J]. 应用数学和力学2023443): 260-271. [百度学术] 

LI ChuangdiJIANG LifuWANG Ruiboet al. Responses of SDOF structures with SPIS-Ⅱ dampers under random seismic excitation[J]. Applied Mathematics and Mechanics2023443): 260-271. [百度学术] 

11

ZHANG R FZHAO Z PPAN Cet al. Damping enhancement principle of inerter system[J]. Structural Control and Health Monitoring2020275): e2523. [百度学术] 

12

潘超韩笑张瑞甫. 基于最大耗能增效原则的惯容减震系统解析设计公式[J]. 工程力学2023404): 91-101. [百度学术] 

PAN ChaoHAN XiaoZHANG Ruifuet al. Closed-form design formula for inerter system based on the principle of maximum damping enhancement[J]. Engineering mechanics2023404): 91-101. [百度学术] 

13

张瑞甫张璐琦潘超. 考虑正常使用功能的非线性黏滞阻尼惯容系统多指标减震控制[J]. 工程科学与技术2023555): 14-22. [百度学术] 

ZHANG RuifuZHANG LuqiPAN Chaoet al. Multi-objective seismic control effect of inerter system with nonlinear viscous damping considering functionality of buildings[J]. Advanced Engineering Sciences2023555): 14-22. [百度学术] 

14

潘超刘媛张瑞甫. 惯容减震系统性能成本控制解析设计方法[J]. 建筑结构学报20224311): 107-116. [百度学术] 

PAN ChaoLIU YuanZHANG Ruifuet al. Performance-cost design method of inerter system based on closed-form formulae[J]. Journal of Building Structures20224311): 107-116. [百度学术] 

15

潘超张瑞甫王超. 单自由度混联Ⅱ型惯容减震体系的随机地震响应与参数设计[J]. 工程力学2019361): 129-137. [百度学术] 

PAN ChaoZHANG RuifuWANG Chaoet al. Stochastic seismic response and design of structural system with series-parallel-Ⅱ inerter system[J]. Engineering Mechanics2019361): 129-137. [百度学术] 

16

WANG DFAN Z LHAO S Wet al. An evolutionary power spectrum model of fully nonstationary seismic ground motion[J]. Soil Dynamics and Earthquake Engineering20181051-10. [百度学术] 

17

何浩祥丁佳伟闫维明. 地震动全非平稳时频谱统计模型及其在结构分析中的应用[J]. 振动工程学报2021346): 1257-1267. [百度学术] 

HE HaoxiangDING JiaweiYAN Weiming. Statistical model of time-frequency spectrum of earthquake and application in structural vibration analysis[J]. Journal of Vibration Engineering2021346): 1257-1267. [百度学术] 

18

PRIESTLEY M B. Evolutionary spectra and non-stationary processes[J]. Journal of the Royal Statistical Society: Series B1965272): 204-229. [百度学术] 

19

霍慧陈国海王文培. 平稳/非平稳激励下中厚圆柱壳随机振动响应的基准解[J]. 力学学报2022543): 762-776. [百度学术] 

HUO HuiCHEN GuohaiWANG Wenpeiet al. Benchmark solutions of random vibration responses for moderately thick cylindrical shells under stationary/nonstationary excitations[J]. Chinese Journal of Theoretical and Applied Mechanics2022543): 762-776. [百度学术] 

20

方同孙木楠张天舒. 相关演变随机激励下响应演变谱矩阵的表达式[J]. 应用力学学报2008253): 351-354. [百度学术] 

FANG TongSUN MunanZHANG Tianshu. General expression of response evolutionary spectrum matrix due to correlated evolutionary random excitations[J]. Chinese Journal of Applied Mechanics2008253): 351-354. [百度学术] 

21

李创第王博文昌明静. 广义Maxwell阻尼耗能系统非均匀响应精细积分精确格式[J]. 桂林理工大学学报2021414): 783-790. [百度学术] 

LI ChuangdiWANG BowenCHANG Mingjing. Precise integration scheme for nonuniform response of energy dissipation system with generalized Maxwell damping[J]. Journal of Guilin University of Technology2021414): 783-790. [百度学术] 

22

方同. 工程随机振动[M]. 北京国防工业出版社1995. [百度学术] 

23

李创第王博文昌明静. 一般线性黏弹性阻尼器保护系统非均匀与完全非平稳地震响应解析分析[J]. 振动工程学报2022355): 1084-1100. [百度学术] 

LI ChuangdiWANG BowenCHANG Mingjing. Analytical analysis of non-uniform and completely nonstationary seismic response of a general linear viscoelastic damper protection system[J]. Journal of Vibration Engineering2022355): 1084-1100. [百度学术] 

24

林家浩张亚辉. 随机振动的虚拟激励法[M]. 北京科学出版社2004. [百度学术] 

25

XING ChenxiWANG HaoXU Zidonget al. Stochastic analysis of a large-span continuous girder high-speed railway bridge under fully non-stationary earthquake[J]. Applied Sciences20221224): 12684. [百度学术] 

26

张步云戴涛TAN Chin-An. 悬架系统振动特性的非平稳虚拟激励法研究[J]. 振动、测试与诊断2022422): 227-234. [百度学术] 

ZHANG BuyunDAI TaoTAN Chin-Anet al. Non‑stationary pseudo excitation method for analyzing vibration characteristics of suspension system[J]. Journal of Vibration, Measurement and Diagnosis2022422): 227-234. [百度学术] 

27

葛新广龚景海李创第. 功率谱二次正交化法在随机地震动响应的应用[J].振动工程学报2022353): 616-624. [百度学术] 

GE XinguangGONG JinghaiLI Chuangdiet al. Application of quadratic orthogonalization method of response power spectrum to random ground motion response[J]. Journal of Vibration Engineering2022353): 616-624. [百度学术] 

28

邹万杰刘美华李创第. 基于胡聿贤谱的带支撑广义Maxwell阻尼隔震结构随机响应分析[J]. 力学学报2022549): 2601-2615. [百度学术] 

ZOU WanjieLIU MeihuaLI Chuangdiet al. Seismic response analysis of generalized Maxwell damping isolated structure with braces under Hu Yuxian spectrum excitation[J]. Chinese Journal of Theoretical and Applied Mechanics2022549): 2601-2615. [百度学术] 

29

葛新广. 分数导数型阻尼耗能结构地震动响应的分析方法研究[D]. 上海上海交通大学2022. [百度学术] 

GE Xinguang. Research on analysis methods for seismic response of energy dissipative structure with dampers modeled by fractional derivative[D]. ShanghaiShanghai Jiao Tong University2022. [百度学术] 

30

CHEN MLIANG X JYANG Z Wet al. Analytical Study on the random seismic responses of an asymmetrical suspension structure[J]. Buildings2023136): 1435. [百度学术] 

31

HOUSNER G W. Characteristics of strong-motion earthquakes[J]. Bulletin of the Seismological Society of America1947371): 19-31. [百度学术]