Analysis of Required Separation Distances of High-Pier Bridges Subjected to Non-stationary Ground Motions
-
摘要: 为了研究非平稳地震作用下高墩桥梁体防撞间隙需求,基于随机振动理论及虚拟激励法,对不同烈度下场地条件对非平稳间隙需求的影响进行了分析. 首先,建立了非平稳地震作用下相邻梁体相对位移需求与烈度间的数学关系;其次,基于理论计算的梁体间最大相对位移,确定碰撞间隙宽度需求以达到防止梁体间发生碰撞的目的;最后,以某大跨度连续刚构-连续梁体系为实际工程算例,研究了非平稳地震作用下桥梁结构在一致场地和非一致场地(实际场地)条件下的碰撞间隙需求量,且获得了不同烈度下非平稳碰撞间隙需求谱. 研究结果表明:非平稳地震作用下,硬土场地条件时,相对位移时变均方差的峰值最小,实际场地条件最大,约为硬土场地的4倍;实际场地条件的各烈度下非平稳碰撞间隙宽度需求均值比软土场地、中土场地和硬土场地分别大36%、69%和73%,均方差分别大45%、74%和78%;平稳地震激励比非平稳地震激励时碰撞间隙需求量大20%~30%.Abstract: To investigate the influence of non-stationary ground motions on required separation distances, the influence of site conditions on non-stationary gap under different intensities was analysed based on the random vibration theory and pseudo excitation method. First, the mathematical formulation clarifying the relationship between the required relative displacement and earthquake intensity was derived. Second, the minimum gap required between two adjacent girds to prevent seismic pounding was determined from the maximum relative displacement based on theoretical analysis. Finally, a real continuous rigid bridge was used to obtain the required gap length spectrum between two adjacent decks for considering non-stationary seismic excitations with different intensities under different local site conditions. Furthermore, the influence of actual site conditions and uniform site conditions on the required separation distance was investigated. Numerical results reveal that actual site conditions have the most obvious effect on the time-dependent mean square deviation of relative peak displacement, which is four times the corresponding value for firm site conditions. The mean values of the required separation distance for different intensities are larger than the corresponding values for uniform site conditions by 36%, 69%, and 73% for soft sites, median sites, and firm sites, respectively, whereas the mean square deviation is 45%, 74%, and 78%, respectively. The required separation distance under stationary seismic excitations is 20%–30% larger than that under non-stationary seismic excitations.
-
随着我国城市化进程不断加快,隧道数量和复杂性都在不断增加. 目前,我国已成为世界上隧道数量最多、条件最复杂、技术发展最快的国家[1-2]. 为解决交通便捷性需求,城市隧道结构逐渐丰富,各类分岔隧道数量也急剧增加.
火灾作为隧道中危害性最大的灾害,一旦发生,将造成严重的人员财产损失. 相较于普通地面建筑,隧道中火灾的烟气蔓延情况及控制手段都更为复杂,使得隧道发生火灾时人员的安全疏散及灭火救援难度高于普通地面建筑[3-4].
针对于单直隧道火灾,国内外进行了大量的研究. Alpert[5]首次提出了预测顶棚射流最高羽流温度的经验模型;Kurioka等[6]通过改变隧道的横截面以及火源的热释放速率进行实验,提出隧道顶棚下的最高温度模型;Li等[7]考虑连续火焰长度低于顶棚高度的情况下,通过理论分析及实验研究,探讨2种不同尺寸隧道中,不同热释放速率、纵向通风风速及隧道几何形状对顶棚下烟气最高温度的影响;王钟宽等[8]以美国Memorial隧道为原型,应用Froude准则建立1∶20缩尺隧道模型,探究不同坡度的隧道火灾自熄特性;姜学鹏等[9]提出侧部点式排烟模式隧道火灾临界风速的无量纲计算式.
随着地下空间的开发与应用,隧道结构逐渐丰富. 因此,针对于分岔隧道火灾的研究也快速展开. Yang等[10]采用理论分析和数值模拟相结合的方法,研究在具备机械排烟系统的倾斜分岔隧道中烟气流动的多重行为;Huang等[11]通过小尺寸分岔隧道火灾实验,研究隧道顶棚下气体最大超限温度,提出一个分为两部分的数学模型来预测分岔隧道的最大上限温度;Zhang等[12]研究T型分岔隧道结构火灾的顶棚最大温升,并提出在目前的模型结构规模和热释放率的研究范围内T型分岔隧道的最大顶烟温度随着无量纲支管宽度的增加而略有下降.
目前,空气幕也被广泛应用于隧道火灾排烟控制中,其工作原理是让高速喷出的空气流形成一道墙,以起到隔烟隔热的效果. 学者们研究了空气幕对火灾烟气蔓延的控制作用. Shu等[13]利用一系列典型的空气动力学性能曲线,建立一种在无风条件下计算空气幕效率系数的新模型,对分岔隧道火灾空气幕控烟有借鉴作用. Yu等[14]提出当动量比约为10∶1时,吹风角度垂直向下的空气幕封烟效果最佳,且其效果随缝宽和喷射角的增大而提升;当火源位于空气幕一侧时,空气幕的最佳喷射角为朝向火源倾斜30°. Li等[15]提出空气幕射流的均匀性受气幕单元的横截面影响,随着射流距离的增加,空气幕的射流速度先急剧下降,然后缓慢下降,直至为0.
隧道火灾中,单独使用纵向通风或空气幕对分岔隧道进行烟气的控制效果仍不尽如人意[16-17]. 而针对纵向通风与空气幕协同作用下的分岔隧道烟气控制及最高温度的研究还较少,故本文根据57组缩尺寸实验数据,对分岔隧道多向气流耦合下的烟气运动和最高温度变化进行研究.
1. 实验布置
如图1所示,在缩小比例(1∶10)分岔隧道模型中进行实验,隧道长10.0 m、宽1.0 m、高0.6 m,分岔隧道长5.0 m、宽1.0 m、高0.6 m,夹角45°. 以液化石油气为燃料,火源位于主隧道中间.
分岔隧道的岔道口处安装有空气幕. 该空气幕采用铁制栅格通过软性连接与风机相连,并通过底部支架固定于隧道内,从而保证气流的稳定性,且各个位置的风速分布均匀. 同时,可通过侧面螺丝及手柄调节空气幕的角度A在 −45°~45° (面向主隧道为正,分岔隧道口为负),空气幕喷射宽度(厚度T)为0~0.15 m. 通过变频器也可随意调节空气幕风速大小,使空气幕射流速度$ {v_{{\mathrm{air}}}} $在0~4.0 m/s. 空气幕的整体结构如图1(b)、(c)所示. 在分岔隧道的左侧安装变频风机,用于控制纵向通风. 为监测隧道中最高温度的分布,将127个K型热电偶沿着主隧道和分岔隧道的中心轴放置,热电偶设置细节如图2所示. 纪杰等[18]研究狭长空间内的烟气流动特性发现,在小尺寸隧道中进行火灾实验,雷诺数Re (式(1))和弗洛德数Fr (式(2))对流体运动具有决定性影响.
Re=vlμ, (1) Fr=glv2, (2) 式中:$ v $为纵向通风速度,m/s;$ l $为长度尺寸,m;$ \mu $为流体的运动黏度.
基于几何相似原理,热量的相似关系可表示为
Qm=Qf(LmLf)5/2, (3) 式中:$ {Q_{\mathrm{m}}} $为模型尺寸热释放速率,kW;$ {Q_{\mathrm{f}}} $为全尺寸热释放速率,kW;$ {L_{\mathrm{m}}} $为模型隧道的长度,m;$ {L_{\mathrm{f}}} $为全尺寸隧道的长度,m.
基于该相似关系,当全尺寸隧道中热释放速率 (HRR)为5、15、25 MW时,缩尺寸隧道模型中分别对应为15.9、47.9、77.7 kW.
本文总共进行了57组实验,变量包括热释放速率、纵向通风速度、空气幕角度、厚度和风速,拟探究纵向通风与空气幕之间的耦合效应. 详细实验工况见表1.
表 1 实验工况Table 1. Experimental conditions工况 HRR/kW 纵向通风/(m•s−1) 空气幕风速/(m•s−1) 角度/(°) 厚度/m 1~15 15.9,47.9,77.7 0.4,0.8,1.2,1.6,2.0 1.0 0 0.15 16~33 47.9 0.4,0.8,1.2 1.5,2.0,2.5,3.0,3.5,4.0 0 0.15 34~51 47.9 0.4,0.8,1.2 3.5 −45,−30,−15,15,30,45 0.15 52~57 47.9 0.4,0.8,1.2 3.5 30 0.05,0.10 2. 不同纵向通风下的隧道温度分布
2.1 主线隧道
当空气幕射流速度较小且恒定($ {v_{{\mathrm{air}}}} $=1.0 m/s)时,3种不同热释放速率下纵向通风大小对主隧道内温度分布影响如图3所示. 图中,横坐标x0为距离主隧道火源的距离. 由图可知,当纵向风速较小时($ v $=0.4 m/s),其对烟气回流层的阻隔作用较小,从而导致隧道上游出现温度波动的情况,而火源热释放速率越大,该现象越明显. 当纵向风速为0.8 m/s及以上时,对主线隧道内烟气的降温效果由纵向通风主导,能在隧道上游形成空气屏障,有效缩短烟气回流层长度. 烟气的蔓延过程为非定常的三维湍流运动[19],而空气幕的作用方向与烟气在顶棚形成的射流方向垂直,导致空气幕通风会对烟气的一维水平蔓延形成剪切作用,在此基础上再添加纵向通风,会对火灾烟气在主线隧道下游蔓延的湍流运动产生较大扰乱,且火源功率越小影响越大. 在3种热释放速率下,随着纵向风速的增加,不仅隧道内最高气温降低,而且烟气最高温度的位置也向隧道下游有一定偏移.
由图3(a)可知:当热释放速率为15.9 kW时,纵向风速从0.4 m/s增加到2.0 m/s,隧道顶棚最高温度从150.1 ℃降低到49.5 ℃;纵向风速达到2.0 m/s时,烟气温度已无明显峰值,火羽流在纵向通风的作用下向火源下游倾斜;当纵向通风增大时,火羽流的倾斜程度增加,对隧道顶棚的直接冲击作用减弱,因此,在隧道顶棚下方形成的高温区域逐渐减小,并且向火源下游偏移. 由图3(b)可知:当热释放速率为47.9 kW时,纵向通风从0.4 m/s增加到2.0 m/s,顶棚最高温度从381.4 ℃降低到83.6 ℃;当纵向风速为2.0 m/s时,隧道顶棚最高温度位置根据火源位置向下游偏移1.0 m. 由图3(c)可知:当热释放速率为77.7 kW时,纵向风速从0.4 m/s增加到2.0 m/s,最高温度从512.3 ℃降低到149.2 ℃,主线隧道下游沿程温度变化随纵向风速的增加而趋于稳定. 综上所述,当空气幕射流速度不变时,随着热释放速率的增大,主线隧道的纵向通风降温效果会达到阈值,因此,有必要探究纵向通风与空气幕协同作用时的最佳风速组合.
2.2 分岔隧道
当$ {v_{{\mathrm{air}}}} $=1.0 m/s时,不同热释放速率下纵向风速大小对分岔隧道温度分布影响如图4所示. 分岔隧道的最高温度随着热释放速率的增大而显著升高,因空气幕良好的隔热阻烟效果,温度在x0=0.8 m处呈断崖式下降. 当纵向风速不同时,分岔隧道中x0=−0.5 m处的温度差异较大. 纵向风速过小(如0.4、0.8 m/s)会导致高温烟气在火源附近滞留,并向分岔隧道蔓延. 然而,加大纵向风速虽能一定程度上避免烟气向分岔隧道蔓延,但会加剧主线隧道下游烟气的沉降效果.
综上,在城市隧道分岔段采用空气幕进行隔热阻烟具有一定的应用价值,其对火灾烟气的具体影响将在下文进行论述. 由于当纵向风速超过1.2 m/s时其降温效果不再显著变化,因此,后续实验选择在纵向风速为0.4、0.8、1.2 m/s条件下进行.
3. 空气幕对隧道最高温升的影响
3.1 空气幕射流速度
为探究不同纵向风速和空气幕射流速度条件下分岔隧道的最大温升,基于Chen等[20]的研究,建立式(4)~(7)所示关系.
ΔTmax=f(Q,Cp,g,v,ρ∞,T∞,H), (4) ΔTmaxT∞=γ(Q2/30Fr1/3)ε, (5) Q0=Qρ∞CPT∞g1/2H5/2, (6) ΔTmaxT∞∼Q0Fr, (7) 式中:$ \gamma 、\varepsilon $为常数,当$ \dfrac{{\mathop Q\nolimits_{\text{0}}^{{2 \mathord{\left/ {\vphantom {2 3}} \right. } 3}} }}{{F{r^{{1 \mathord{\left/ {\vphantom {1 3}} \right. } 3}}}}} < 1.35 $时,$ \gamma =1.77,\varepsilon = \text{1}\text{.2} $,当$ \dfrac{{\mathop Q\nolimits_{\text{0}}^{{2 \mathord{\left/ {\vphantom {2 3}} \right. } 3}} }}{{F{r^{{1 \mathord{\left/ {\vphantom {1 3}} \right. } 3}}}}} \geqslant 1.35 $时,$ \gamma =\text{2}\text{.54},\varepsilon =\text{0} $;$ \Delta {T_{\max }} $为隧道顶棚最高温升,K;$ Q $为热释放速率,kW;Cp为恒压比热容,kJ/(kg•K);$ {\rho _\infty } $为环境空气密度,kg/m3;$ {T_\infty } $为环境温度,K;H为燃烧器表面到隧道顶棚的距离,m;$ \mathop Q\nolimits_{\text{0}} $为无量纲热释放速率.
根据陈贻来[21]的研究,单个汽车燃烧时的热释放速率为3~5 MW,当隧道中单个汽车发生火灾并引燃两侧汽车时,火灾热释放速率最高可达15 MW. 而夏永旭等[22]研究表明,大巴车在隧道中平均热释放速率为8 MW. 因此,为更好地与实际相结合,探究隧道火灾极端情况下不同纵向风速与空气幕对隧道顶棚温度的影响,后续实验选取热释放速率为15 MW,即在缩尺寸隧道中热释放速率为47.9 kW时进行. 因此,在该情况下,$ {Q_0} $可以作为一个常数而忽略,故而式(7)可以写为
ΔTmaxT∞∼Fr∼v2gH. (8) 图5显示在不同空气幕射流速度下$ {{{v^2}} / {gH}} $与$ {{\Delta {T_{\max }}} / {{T_\infty }}} $之间的关系. 图中,x、y分别为拟合曲线横、纵坐标参量,余图同. 由此可以得出,在不同空气幕射流速度的影响下,分岔隧道的最大温升与纵向通风之间的关系如式(9)所示.
ΔTmaxT∞=m(v2gH)−n, (9) m=−0.2vair√gH+1.75, (10) n=−0.08vair√gH+0.18, (11) 式中:m、n为系数.
拟合的详细描述见图6.
将m和n代入式(9),便可以得到Q = 47.9 kW,空气幕角度为0° 时,纵向通风与空气幕协同作用下分岔隧道最大温升的理论模型,如式(12)所示.
ΔTmaxT∞=(−0.2vair√gH+1.75)(v2gH)0.08vair√gH−0.18. (12) 将不同空气幕射流速度下主隧道最大温升实验值与理论模型(式(12))中的预测值进行比较,如图7所示,其结果吻合度较高.
3.2 空气幕角度和厚度
根据图5可以看出,在3种纵向通风下,当空气幕射流速度达到3.5 m/s时,空气幕射流速度对主隧道最高温影响减弱,此时,空气幕和纵向通风产生的气流达到动态平衡. 因此,为探索空气幕的角度A和厚度T对隧道温度分布的影响,将空气幕的射流速度$ {v_{{\mathrm{air}}}} $设置为固定值3.5 m/s,以排除其影响.
不同空气幕角度下,$ {v^2}/(gH) $与$ {{\Delta {T_{\max }}} / {{T_\infty }}} $之间的关系如图8所示. 由此可以得出,在不同空气幕角度(−45° 除外)的影响下,分岔隧道最大温升与纵向通风的关系如图9所示,其拟合式如式(13)所示.
图 8 不同空气幕角度下$ {{{v^2}} \mathord{\left/ {\vphantom {{{v^2}} {(gH)}}} \right. } {(gH)}} $与$ {{\Delta {T_{\max }}} \mathord{\left/ {\vphantom {{\Delta {T_{\max }}} {{T_\infty }}}} \right. } {{T_\infty }}} $的关系Figure 8. Relationship between $ v^{2} / ({gH}) $and $ \Delta T_{\max } / T_{\infty} $at different angles of air curtainΔTmaxT∞=αv2gH+β, (13) α = −0.9sin2A+0.8sinA−0.3, (14) β = −0.61sinA+1.57, (15) 式中:α、β为拟合系数.
将系数$ \alpha $和$ \beta $代入式(13)中,得到热释放速率为47.9 kW,$ {v_{{\mathrm{air}}}} $=3.5 m/s时,纵向通风与空气幕协同作用下空气幕角度对分岔隧道最大温升影响的理论模型,如式(16)所示.
ΔTmaxT∞=(−0.9sin2A+0.8sinA−0.3)v2gH−0.61sinA+1.57. (16) 如图10所示,当火源功率47.9 kW,$ {v_{{\mathrm{air}}}} $=3.5 m/s 时,不同空气幕角度下,隧道最大温升拟合得到的理论公式较好与实验结果基本一致.
如图11所示,$ {v_{{\mathrm{air}}}} $=3.5 m/s,A=30° 条件下进行了3种不同空气幕厚度的实验,当T=0.10 m时分岔隧道降温效果最佳. 由图11可以得出,当vair=3.5 m/s,A=30° 时,在不同空气幕厚度的影响下,分叉隧道最大温升与纵向通风之间的关系如式(17)所示.
ΔTmaxT∞=av2gH+bv√gH+c, (17) 式中:a、b、c为拟合系数,拟合的详细描述见图12.
将系数a、b和c代入式(17)中,便可以得到热释放速率为47.9 kW时,纵向通风和空气幕协同作用下空气幕厚度对分岔隧道最大温升影响的理论模型,如式(18)所示.
ΔTmaxT∞=av2gH−bv√gH+24T2H2−13TH+2.7, (18) 式中:a = −26.88T2/H2−12T/H−0.44,b = −44.8T2/H2 +22T/H−1.3.
将不同空气幕喷射宽度下主隧道最大温升实验值与理论模型(式(18))中的预测值进行比较,如图13所示,具有较好的吻合性.
4. 结 论
1) 空气幕能有效协助纵向通风降低主线隧道温度,防止烟气进入分岔隧道.
2) 当空气幕射流速度较小时($ {v_{{\mathrm{air}}}} $=1.0 m/s),纵向风速的增加能有效防止烟气在分岔口积聚,提高空气幕对烟气的阻隔效率,降低隧道分岔点处的烟气温度.
3) 根据隧道火灾顶棚最大温升的无量纲经验相关公式,提出在固定热释放速率47.9 kW条件下,隧道最大温升与纵向通风和空气幕射流速度之间的关系,该模型的预测结果与实验结果吻合度较高.
4) 在固定热释放速率47.9 kW,空气幕射流速度3.5 m/s条件下,提出分岔隧道最大温升与纵向风速和空气幕角度之间的关系,其预测结果与实验结果吻合度较高.
5) 在固定热释放速率47.9 kW,空气幕射流速度3.5 m/s及空气幕角度30° 条件下,提出分岔隧道最大温升与纵向风速和空气幕厚度之间的关系,其预测结果与实验结果吻合度较高.
-
表 1 3类场地类型参数
Table 1. Type parameters for three sites
场地类型 ${\omega _{\rm{g}}}$ ${\xi _{\rm{g}}}$ ${\omega _{\rm{f}}}$ ${\xi _{\rm{f}}}$ 硬土 15.0 0.6 1.5 0.6 中土 10.0 0.4 1.0 0.6 软土 5.0 0.2 0.5 0.6 表 2 地震动强度比例因子
Table 2. Scale factor of ground motion intensity
场地类型 烈度/度 6 7 8 9 10 硬土 3.69 14.77 59.08 236.31 945.13 中土 3.97 15.89 63.53 254.07 1 016.17 软土 5.22 20.87 83.45 333.76 1 334.86 表 3 墩底工况分析
Table 3. Analysis of pier bottom condition
工况编号 烈度/度 1# 墩 2# 墩 3# 墩 4# 墩 5# 墩 6S 6 S S S S S 6M 6 M M M M M 6F 6 F F F F F 6R 6 M S F M S 7S 7 S S S S S 7M 7 M M M M M 7F 7 F F F F F 7R 7 M S F M S 8S 8 S S S S S 8M 8 M M M M M 8F 8 F F F F F 8R 8 M S F M S 9S 9 S S S S S 9M 9 M M M M M 9F 9 F F F F F 9R 9 M S F M S 10S 10 S S S S S 10M 10 M M M M M 10F 10 F F F F F 10R 10 M S F M S -
李乔, 赵世春. 汶川大地震工程震害分析[M]. 成都: 西南交通大学出版社, 2008: 18-27 KAWASHIMA K, UNJOH S, HOSHIKUMA J, et al. Damage of bridges due to the 2010 Maule,Chile,earthquake[J]. Journal of Earthquake Engineering, 2011, 15(7): 1036-1068 CHOUW N, HAO H. Pounding damage to buildings and bridges in the 22 February 2011 Christchurch earthquake[J]. International Journal of Protective Structures, 2012, 2(3): 123-140 HE L, SHRESTHA B, HAO H, et al. Experimental and three-dimensional finite element method studies on pounding responses of bridge structures subjected to spatially varying ground motions[J]. Advances in Structural Engineering, 2017, 20(1): 105-124 doi: 10.1177/1369433216646009 HAO H. A parametric study of the required seating length for bridge decks during earthquake[J]. Earthquake Engineering & Structural Dynamics, 2015, 27(1): 91-103 RUANGRASSAMEE A, KAWASHIMA K. Relative displacement response spectra with pounding effect[J]. Earthquake Engineering & Structural Dynamics, 2001, 30(10): 1511-1538 白凤龙. 空间变化地震动激励下大跨度结构的反应研究[D]. 大连: 大连理工大学, 2010 CHOUW N, HAO H. Significance of SSI and non-uniform near-fault ground motions in bridge response I:effect on response with conventional expansion joint[J]. Engineering Structures, 2008, 30(1): 141-153 CHOUW N, HAO H. Significance of SSI and non-uniform near-fault ground motions in bridge response II:effect on response with modular expansion joint[J]. Engineering Structures, 2008, 30(1): 154-162 李忠献,周莉,岳福青. 地震动空间效应与土-基础相互作用效应对隔震桥梁临界碰撞间隙的影响[J]. 土木工程学报,2010,43(7): 85-90LI Zhongxian, ZHOU Li, Yue Fuqiang. Effects of spatial variation of ground motion and siol-foundation interaction on critical pounding gap length of seismic isolated bridges[J]. Journal of Civil Engineering, 2010, 43(7): 85-90 李忠献,岳福青,周莉,等. 基于随机振动理论确定桥梁地震碰撞的临界间隙[J]. 地震工程与工程振动,2006,26(4): 156-161LI Zhongxian, YUE Fuqing, ZHOU Li, et al. Determination of critical gap length of seismic pounding for bridges based on random vibration theory[J]. Earthquake Engineering and Engineering Vibration, 2006, 26(4): 156-161 BI K, HAO H, CHOUW N. Influence of ground motion spatial variation,site condition and SSI on the required separation distances of bridge structures to avoid seismic pounding[J]. Earthquake Engineering & Structural Dynamics, 2011, 40(9): 1027-1043 贾宏宇,杜修力,李晰,等. 地震作用下高墩铁路桥梁梁体碰撞间隙宽度需求机理分析[J]. 工程力学,2017,34(2): 207-215JIA Hongyu, DU Xiuli, LI Xi, et al. Demand mechanism analysis on pounding separation distance of high pier railway bridges subjected to earthquake excitations[J]. Engineering Mechanics, 2017, 34(2): 207-215 JIA H Y, ZHANG D Y, ZHENG S X, et al. Local site effects on a high-pier railway bridge under tridirectional spatial excitations:Nonstationary stochastic analysis[J]. Soil Dynamics and Earthquake Engineering, 2013, 52(6): 55-69 doi: 10.1016/j.soildyn.2013.05.001 RUIZ P, PENZIEN J. Probabilistic study of the behavior of structures during earthquake[Z]. University of California, Berkeley: Earthquake Engineering. Research Center, 1969 ZHANG D Y, JIA H Y, ZHENG S X, et al. A highly efficient and accurate stochastic seismic analysis approach for structures under tridirectional nonstationary multiple excitations[J]. Computers & Structures, 2014, 145(C): 23-35 LIN J, ZHAO Y, ZHANG Y. Accurate and highly efficient algorithms for structural stationary/non-stationary random responses[J]. Computer Methods in Applied Mechanics and Engineering, 2001, 191(1): 103-111 doi: 10.1016/S0045-7825(01)00247-X DAVENPORT A G. Note on the distribution of the largest value of a random function with application to gust loading[J]. Process Institute Civil Engineering, 2015, 28(2): 187-196 江近仁,洪峰. 功率谱与反应谱的转换和人造地震波[J]. 地震工程与工程振动,1984,4(3): 1-11JIANG Jinren, HONG Feng. Conversion between power spectrum and response spectrum and artificial earthquakes[J]. Earthquake Engineering and Engineering Vibration, 1984, 4(3): 1-11 -