Influence of permanent displacement of ground motion on seismic response of structures
-
摘要: 本文推导了基于位移激励计算单自由度体系拟速度谱公式,通过构造的脉冲位移时程对公式精度进行了验证;之后利用小波变换去除强震记录噪声而保留地震动永久位移,再基于去趋势项方法和滤波方法去除永久位移后,计算拟速度谱。算例结果表明:短周期段内,不保留永久位移的位移激励拟速度谱值与保留永久位移的位移激励拟速度谱值相差很小;中长周期段内,不保留永久位移的位移激励拟速度谱值总体上小于保留永久位移的位移激励拟速度谱值,且不保留永久位移时,滤波方法引起的拟速度谱降幅大于去趋势项方法所引起的拟速度谱降幅。因此,基于位移激励计算中长周期结构的地震反应时,应保留地震动永久位移,或基于去趋势项方法去除永久位移。Abstract: When calculating the seismic response of a structure based on the displacement excitation, the displacement excitation is obtained by integration of the velocity, which is obtained by integration of the acceleration records. At the same time, the acceleration noise must be removed. However, limited by the method for eliminating noise from acceleration records, the permanent displacement of the ground motion in the near-fault earthquake acceleration records is often weakened or even eliminated. This makes it impossible for researchers to judge the effect of permanent displacement of ground motion on the seismic response of structures excited by displacement. In this paper, the authors deduce the pseudo-velocity spectra formula of single degree of freedom system based on displacement excitation, verify the accuracy of the formula with the constructed pulse displacement time history, and then analyze the effect of the permanent displacement of the ground motion on the pseudo-velocity spectra. Firstly, the noise is executed on the acceleration records via wavelet transform, and gound motion with permanent displacement is remained. And then the permanent displacement is wiped out using reducing trend item method or filtering method from those records. Those excitations are used to calculate pseudo-velocity spectra. The results show that in short period, the pseudo-velocity spectra are similar whether ground motions contain the permanent displacement or not; in medium-long period, the pseudo-velocity spectra excited by non-permanent displacement excitation are lower than those by displacement excitation with permanent displacement, and this kind of difference are enlarged when the filtering method is used. Therefore, for seismic design of long-period structures, the ground motion with permanent ground displacement excitation should be used, or reducing trend item method is used to wipe out permanent displacement.
-
引言
国内外关于局部场地地震效应的研究始于20世纪初(胡聿贤等,1980;Trifunac,2016)。Wood (1908)在讨论1906年发生于美国旧金山MS7.8大地震的震害特征时,指出局部场地条件对震害有显著的影响,在这次地震中,位于软弱场地上的震害明显重于其它场地;Yamanaka和Yamada (2003)针对1923年日本关东MS8.1地震的震害研究发现,软土场地上木房屋的震害明显有别于刚性房屋;Seed等(1988)关于1985年墨西哥MS7.8地震的震害研究表明场地地震效应显著,墨西哥城下伏较厚的软弱沉积层是引起地震动显著放大的原因;Chávez-García和Bard (1994)对1985年墨西哥MS7.8地震作了进一步分析研究,大量的实验和理论结果表明,墨西哥城的严重震害是由盆地内深厚软弱的沉积层导致地震动显著放大引起。再如1960年的智利MS8.9地震(Doyel et al,1963)和1990年伊朗MS7.3大地震(Ishihara et al,1992)等都表现出了场地效应对地震动的影响。国内关于场地效应的研究始于20世纪70年代的云南通海地震。通海地震发生后,胡聿贤等对震区震害异常进行了宏观调查和统计,随后经过数年的地震震害、强震动记录和理论与实验分析研究了局部场地地震动效应(通海地震影响场调查组,1977),其关于场地效应的相关研究成果在我国抗震规范的制定中得到了体现。1975年辽宁海城MS7.3地震中,海城县和大石桥市的一般砖屋结构的震害远高于营口市和盘锦市,定性分析认为这一震害异常是由海城县和大石桥地区较厚的覆盖土层引起的(中国科学院工程力学研究所,1979)。此外,在河北唐山MS7.8地震(高振寰,胡碧茹,1987)和台湾集集MS7.6地震(Wang et al,2002)中均存在局部场地影响地震作用的现象。
基于地脉动的场地特征研究起始于20世纪30年代,随着基于同一测点地脉动水平垂直谱比法(H/V谱比法),即Nakamura方法的提出(Nakamura,1989),地脉动法在场地特征估计中得到了更广泛的应用(Yamazaki,Ansary,1997;Poovarodom et al,2017;Verdugo,2019;Benkaci et al,2021;Thornley et al,2021;Ramadan et al,2021)。温瑞智等(2013)为了进一步说明场地地形对四川省宝兴县防震减灾局强震台的影响,对台站附近进行了背景噪声测试,利用H/V谱比法计算了场地背景噪声的卓越周期和水平向放大系数,揭示了场地卓越频率随地势自山脚向山坡上的增高而逐渐增大,且放大系数也随之增大。单振东等(2014)采用H/V谱比法测试了思经乡希望小学附近场地及典型民居的特征频率,结合震后现场灾害考察,分析了场地特征与建筑物震害间的相关性,总结出思经乡异常区产生的原因为思经乡场地特征频率与现有民用两层建筑特征频率相近。
鉴于鲁甸龙头山集镇老行政机关办公区、老街区及骡马口村区在鲁甸MS6.5地震中所受震害明显重于老中学区(图1),除有房屋建筑结构本身的因素外,局部场地条件可能亦是主要影响因素。为此,本文拟对龙头山集镇老行政机关办公区、老街区及骡马口村区三个观测实验区进行地脉动观测,以期通过地脉动傅里叶幅值谱比了解各场地特征,即场地卓越频率和场地放大系数,确定各场地的地震动放大效应,并通过比较房屋建筑的动力特性,揭示房屋建筑震害加重机理。
1. 龙头山集镇地脉动测试
1.1 地脉动观测
老行政机关办公区、老街区和骡马口村区等三个观测实验区的分布位置如图2所示,各观测实验区地脉动观测点布局如图3所示,在每个观测实验区内将测点尽可能布设在震害较重的建筑场地或附近场地上。同时亦在龙头山集镇的两处基岩露头位置各布置一参考观测点,以作为传递函数法谱比参考点,其中老行政机关办公区、老街区观测实验区参考观测点均为基岩测点1 (图3a,3b),骡马口村区观测实验区参考观测点为基岩测点2 (图3c)。现场测试所采用的仪器为ETNA2加速度仪,观测时间为晚上11点后的5个小时,以尽量避免外界干扰。
图 3 老行政机关办公区(a)、老街区(b)和骡马口村区(c)测点布设示意图括号中数据为依据参考点谱比层得到的EW向和NS向场地放大系数Figure 3. Schematic layout of measuring points in the old administrative-office zone (a), old blocks (b) and Luomakou village (c)The data in parentheses are the amplification coefficients of the EW and NS directions obtained based on the spectral ratio of the reference point1.2 脉动数据处理
首先依据各测点记录起止时间确定出每个测试排列的有效记录时间,并据此从参考点连续记录中截取出与各测点相同的有效记录时间段内的记录数据。通过专用数据处理软件对每一测试排列中各测点记录及对应时段参考点记录进行处理(魏来,2018),处理内容如下:① 基线校正。由于温度、仪器布置不水平等导致实际记录存在一定的零位漂移,为此将通过趋势拟合得到每一测点及对应参考点记录的基线函数,并从每一测点及对应参考点实际记录中减去相应的基线函数值,实现每一测点及对应参考点记录的基线校正;② 数字滤波。为了减少噪音干扰,对经过基线校正的记录进行滤波处理,频带范围为0.1—25 Hz;③ 从滤波处理后的数据中连续截取40—43段数据,每个数据段包含8 192个采样点(时长40.96 s);④ 对截取的40—43段数据进行累加求和再平均,得到每一测点和对应参考点的平均记录,新的数据采样点为8 192个;⑤ 针对每一测点及对应参考点平均记录进行傅里叶变换,继而求其傅里叶振幅谱;⑥ 按照每一测试排列,对各测点及对应参考点平均记录的傅里叶振幅谱采用滑动平均法进行平滑。
1.3 测点谱比曲线特征
经脉动观测数据处理即可得到每一测试排列各测点及对应参考点三分量平滑傅里叶振幅谱,继而计算出各测点H/V谱比及各测点水平向傅里叶振幅谱与参考点水平向傅里叶振幅谱的谱比,即参考点谱比法谱比。参考点H/V谱比曲线如图4所示,从图4可知,基岩场点H/V谱比曲线呈一条幅值变化较小的曲线,无明显峰值,表现为基岩场地H/V谱比曲线特征,即单点H/V谱比曲线多为单峰值,而参考点谱比法谱比曲线多为双峰值。老行政机关办公区、老街区及骡马口村区典型观测点单点H/V谱比曲线与参考点谱比法谱比曲线如图5所示。
图 5 老行政机关办公区 (a)、老街区 (b) 和骡马口村区 (c) 测点及参考点谱比曲线图a,b参考点为基岩测点1,图c参考点为基岩测点2Figure 5. Spectral ratio curves of measuring points reference points in old administrative-office zone (a),old blocks (b) and Luomakou village (c)Reference points in Fig.(a) and Fig.(b) are base bedrock measuring point 1,while is 2 for Fig.(c)2. 场地特征与震害影响分析
经脉动测试数据处理分析得到了三个观测实验区各测点的地脉动谱比,其中包括参考点谱比法谱比及H/V谱比,据此可确定出三个观测实验区各测点对应场地的卓越频率和场地放大系数,结果列于表1,并在图3中标注了场地放大系数。
表 1 老行政机关办公区、老街区和骡马口村区各测点谱比计算结果Table 1. Statistics of the spectral ratio of the measured points in old administrative-office zone, old blocks and Luomakou village测点号 参考点谱比法卓越频率/Hz 参考点谱比法场地放大系数 H/V谱比法卓越频率/Hz EW NS EW NS EW NS X1 8.08 7.86 2.79 2.93 4.01 3.88 X2 4.15 4.35 2.09 2.67 4.00 3.91 X3 4.01 3.10 2.01 1.80 3.98 3.12 X4 4.12 4.66 2.53 2.43 3.91 4.57 X5 3.96 5.98 1.96 1.69 3.76 3.64 X6 4.64 4.66 1.39 1.47 4.76 3.96 X7 3.91 4.10 4.63 3.63 3.86 3.32 X8 3.71 3.81 2.82 3.08 3.47 3.50 X9 5.70 6.75 4.42 3.33 5.69 6.64 X10 5.64 6.23 4.20 4.00 3.96 4.48 X11 3.49 6.75 5.96 4.92 3.34 5.64 由表1和图3可以看出,老行政机关办公区各测点的参考点谱比法中南北方向上的场地放大系数多数与东西方向上的场地放大系数存在差异,且场地放大系数均大于1,表明老行政机关办公区的场地具有显著的场地放大效应,参考点谱比法卓越频率两方向有一定的差异,主要介于3—8 Hz频带范围内。各测点H/V谱比法卓越频率主要介于2—6 Hz频带范围内,基本与参考点谱比法优势频带相吻合。
老行政机关办公区多数房屋为2—3层砖混结构,其基频约为4—6 Hz (王广军,樊水荣,1989年;中华人民共和国建设部,1990),位于脉动测试分析得到的场地卓越频带范围内,因此,地震动作用下砖混结构地震反应的共振效应会加重其破坏程度。龙头山集镇老行政机关办公区位于洪积扇西缘侵蚀陡坡上方,多属残坡积堆积地貌,主要以台阶地形为主,这样的地貌特征使其对地震动具有显著的放大作用,这一点已被前述谱比分析所证实,老行政机关办公区场地的显著地震动放大作用导致了该区砖混结构建筑的严重破坏。
综上所述,老行政机关办公区房屋建筑所受的严重震害是由其场地显著的地震动放大效应及共振效应的共同作用引起的。
老街区各测点参考点谱比法南北方向上的场地放大系数与东西方向上的较接近,且场地放大系数均大于1 (表2)。此外,场地东西向与南北向的卓越频率存在一定的差异,主要位于3—8 Hz频带范围内。表明老街区场地具有明显的场地放大效应,且场地放大效应具有一定的方向性。位于幼儿园建筑附近的测点S3,S4,S5,S9,S10的场地放大系数为1—2,相对于其它测点(大于3)偏低,说明老街区内场地效应存在一定的差异,幼儿园附近场地的地震动放大效应比其它场地的小,这一特征可能是幼儿园附近房屋建筑的震害较周围房屋震害轻的原因。各测点H/V谱比法卓越频率与参考点谱比法的卓越频率差异较明显,主要介于3—6 Hz频带范围内,较参考点谱比法优势频带窄。
表 2 老街区各测点谱比计算结果Table 2. statistics of the spectral ratio of measured points in old blocks测点号 参考点谱比法卓越频率/Hz 参考点谱比法场地放大系数 H/V谱比法卓越频率/Hz EW NS EW NS EW NS S1 3.62 3.76 3.51 3.33 3.59 3.42 S2 5.35 4.35 4.22 4.31 4.08 3.88 S3 4.74 4.66 1.27 1.07 4.81 4.57 S4 4.66 7.64 1.65 1.66 2.34 4.93 S5 5.30 3.93 1.88 1.85 5.35 5.05 S6 4.27 3.98 3.18 3.95 4.74 3.47 S7 4.30 4.03 3.22 3.81 4.66 3.86 S8 8.55 3.39 2.71 3.64 5.18 3.86 S9 6.62 4.47 1.29 1.27 3.17 4.57 S10 6.98 5.54 1.68 1.28 4.03 3.08 老街区地震震害显现出一定的分区特征,这可能与近地表覆盖土层结构有关。科考期间进行的工程地质钻探结果显示,尽管两个钻孔(DZK4和DZK5)之间的距离不足30 m,均位于洪积扇中部,亦同属于中硬场地土,建筑场地类别为Ⅱ类,但松散覆盖层厚度明显不同,其中位于幼儿园场地上的DZK4钻孔揭露的松散覆盖层厚度为30.1 m,而位于老旧房屋处的DZK5钻孔揭露的松散覆盖层厚度仅为20.1 m,这使得地震动的放大效应不同,从而导致了地震震害的差异。
老街区砖混结构建筑多数为两层,其基频约为6 Hz (王广军,樊水荣,1989;中华人民共和国建设部,1990),位于测试区场地卓越频带范围内,因此,地震动作用下砖混结构地震反应的共振效应会加重其破坏。此外,老街区位于洪积扇东缘侵蚀陡坡上,多属残积和坡洪积堆积地貌,主要以坡地地形为主,复杂的场地条件使得场地动力特性表现出一定的差异。地脉动测试结果表明,幼儿园附近场地的地震动放大效应比其它场地小,因而幼儿园及其附近房屋建筑破坏程度较周边的轻,而幼儿园周边场地的显著地震动放大作用导致了砖混结构与土木结构的严重破坏甚至倒塌,证实了前文关于幼儿园建筑震害较周边房屋建筑轻的原因的推断。
土木结构多数为一层,其基频约为8 —10 Hz (王广军,樊水荣,1989;中华人民共和国建设部,1990),与场地脉动测试卓越频率范围不吻合,未能与结构地震反应产生共振效应,但显著的地震动放大效应导致了土木结构的严重破坏甚至倒塌。
综上所述,老街区场地效应存在一定的差异是导致震害差异的主要原因,幼儿园及其附近房屋建筑破坏程度较周边的轻;幼儿园周边场地的砖混建筑震害较重是由其场地的地震动放大效应及结构地震反应的共振效应共同作用引起的,土木建筑的破坏主要是由场地地震动放大效应引起的。
骡马口村区各测点参考点谱比法南北方向上的场地放大系数与东西方向上的场地放大系数差异明显,且谱比峰值均大于5 (表3),表明骡马口村区测试场地具有显著的场地地震动放大效应,且具有一定的方向性;参考点谱比法卓越频率两方向具有一定的差异,主要介于5—8 Hz频带范围内。各测点H/V谱比卓越频率亦主要介于5—8 Hz频带范围内,与参考点谱比法优势频带相吻合。
表 3 骡马口村区各测点谱比计算结果Table 3. Luomakou village measured spectral ratio statistics测点号 参考点谱比法卓越频率/Hz 参考点谱比法场地放大系数 H/V谱比法卓越频率/Hz EW NS EW NS EW NS L1 5.47 7.50 5.19 3.53 5.57 7.49 L2 5.53 5.36 5.81 6.25 5.44 5.04 L3 5.24 5.01 7.79 9.19 5.13 4.80 L4 5.47 5.30 5.95 6.25 5.45 5.04 由鲁甸县龙头山集镇震害调查可知,骡马口村区多数房屋建筑为一层土木与砖木结构,其基频约为8—10 Hz (王广军,樊水荣,1989;中华人民共和国建设部,1990),与场地脉动测试卓越频率范围不吻合,未有结构地震反应的共振效应,但显著的地震动放大效应可导致土木与砖木结构的严重破坏甚至倒塌。
由此可知,骡马口村区土木与砖木结构房屋(不包含木框架结构)地震震害较为严重,除其抗震性能较差外,也与显著的场地地震动放大效应有关。
3. 讨论与结论
基于鲁甸龙头山集镇老行政机关办公区、老街区及骡马口村区三个震害严重区域的地脉动测试记录的频谱分析,并结合各区内建筑结构特点、局部场地效应及震害进行综合分析,结果表明:
1) 老行政机关办公区场地具有显著的场地放大效应,且场地放大效应具有方向上的差异;参考点谱比法卓越频率东西与南北方向有一定的差异,主要介于3—8 Hz频带范围内;各测点H/V谱比法卓越频率主要位于2—6 Hz频带范围内,基本与参考点谱比法优势频带相吻合。老行政机关办公区房屋建筑震害较重是由其场地显著的放大效应及场地和结构体系的共振效应共同作用引起的。
2) 老街区场地具有明显的场地放大效应,且具有一定的方向性差异;老街区场地放大效应的差异性导致了震害的分区性;参考点谱比法场地两方向卓越频率存在一定的差异,主要介于3—8 Hz频带范围内;各测点H/V谱比法卓越频率与参考点谱比法的卓越频率差异较明显,主要介于3—6 Hz频带范围内,较参考点谱比法优势频带窄。老街区内砖混建筑震害较重是由其场地的放大效应及场地和结构体系的共振效应共同作用引起的,土木建筑主要是由场地放大效应引起的。
3) 骡马口村区测试场地具有显著的场地放大效应,且具有一定的方向性差异;参考点谱比法卓越频率两方向具有一定的差异,主要介于5—8 Hz频带范围内。各测点H/V谱比卓越频率亦主要介于5—8 Hz频带范围内,与参考点谱比法优势频带相吻合。骡马口村区测试场地上的房屋震害较重是由其场地显著的放大效应引起的。
-
图 4 加速度激励、无永久位移的位移激励及包含永久位移的位移激励
(a) 欧姆林公寓邦德角记录到的2010年4月4日美国卡莱克西科市地震;(b) 卡莱克西科市消防站记录到的2010年4月4日美国卡莱克西科市地震:(c) 林肯作物与食物研究中心记录到的2010年9月3日新西兰地震;(d) 里奇克莱斯特395号公路/布朗路桥记录到的2019年7月5日美国里奇克莱斯特地震
Figure 4. Acceleration excitations,displacement excitations containing permanent ground displacement and displacement excitations without permanent ground displacement
(a) Calexico earthquake in USA on 4 April 2010 recorded by the Bond’s Corner of Omlin Residence;(b) Calexico earthquake in USA on 4 April 2010 recorded by the Calexico Fire Station;(c) New Zealand earthquake on 3 September 2010 recorded by the Lincoln Crop and Food Research;(d) Ridgecrest earthquake in USA on 5 July 2019 recorded by the Ridgecrest-Highway 395/Brown Road Bridge
图 5 不同点位的拟速度谱对比图(单自由度体系阻尼比ζ=0.02)
(a) 欧姆林公寓邦德角;(b) 卡莱克西科市消防站;(c) 林肯作物与食物研究中心;(d) 里奇克莱斯特395号公路/布朗路桥
Figure 5. Comparing of pseudo-velocity spectra (damping ratio ζ=0.02)for different sites
(a) The Bond’s Corner,Omlin Residence;(b) The Calexico Fire Station;(c) The Lincoln Crop and Food Research; (d) The Ridgecrest-Highway 395/Brown Road Bridge
表 1 地震动信息
Table 1 Ground motion information
序号 地震 发震时间
年−月−日ML 观测位置 断层距/km 峰值加速度
/(cm·s−2)1 美国加州卡莱克西科市
(Calexico)地震2010−04−04 7.2 卡莱克西科市欧姆林(Omlin)
公寓邦德角(Bond’s Corner)35.7 120.32 2 美国加州卡莱克西科市
(Calexico)地震2010−04−04 7.2 卡莱克西科市消防站 23.1 266.45 3 新西兰地震 2010−09−03 7.0 新西兰林肯作物与食物研究中心 6.3 426.80 4 美国加州里奇克莱斯特
(Ridgecrest)地震2019−07−05 7.1 加州里奇克莱斯特395号
公路/布朗路桥22.3 237.74 -
陈勇,陈鲲,俞言祥. 2007. 用集集主震记录研究近断层强震记录的基线校正方法[J]. 地震工程与工程振动,27(4):1–7. Chen Y,Chen K,Yu Y X. 2007. Base line correction method for near-fault accelerograms using Chi-Chi main shock record[J]. Earthquake Engineering and Engineering Vibration,27(4):1–7 (in Chinese).
惠迎新,王克海. 2015. 基于多点激励位移输入模型的跨断层桥梁地震动输入方法[J]. 东南大学学报(自然科学版),45(3):557–562. Hui Y X,Wang K H. 2015. Earthquake motion input method for bridges crossing fault based on multi-support excitation displacement input model[J]. Journal of Southeast University (Natural Science Edition)
,45(3):557–562 (in Chinese). 李建波,陈健云,林皋. 2004. 相互作用分析中地震动输入长周期校正研究[J]. 大连理工大学学报,44(4):550–555. Li J B,Chen J Y,Lin G. 2004. Study of long-period correction of seismic accelerogram for dynamic interaction analysis[J]. Journal of Dalian University of Technology,44(4):550–555 (in Chinese).
李吉涛,杨庆山. 2010. 地震波基线漂移的处理方法[J]. 北京交通大学学报,34(1):95–99. Li J T,Yang Q S. 2010. A correction method for baseline drift of seismic wave[J]. Journal of Beijing Jiaotong University,34(1):95–99 (in Chinese).
荣棉水,李小军,卢滔. 2007. 局部地形对入射P波谱特性的影响分析[J]. 西北地震学报,29(4):297–302. Rong M S,Li X J,Lu T. 2007. Effect analysis of topography on the spectrum property of incident P waves[J]. Northwestern Seismological Journal,29(4):297–302 (in Chinese).
荣棉水,彭艳菊,喻畑,杨宇. 2014. 近断层强震观测记录基线校正的优化方法[J]. 土木工程学报,47(增刊2):300–306. Rong M S,Peng Y J,Yu T,Yang Y. 2014. Optimized baseline correction method for the near-fault observation strong motion records[J]. China Civil Engineering Journal,47(S2):300–306 (in Chinese).
王国权,周锡元. 2004. 921台湾集集地震近断层强震记录的基线校正[J]. 地震地质,26(1):1–14. Wang G Q,Zhou X Y. 2004. Baseline correction of near fault ground motion recordings of the 1999 Chi-Chi,Taiwan earthquake[J]. Seismology and Geology,26(1):1–14 (in Chinese).
王体强,王永志,袁晓铭,汤兆光,王海,段雪锋. 2019. 基于振动台试验的加速度积分位移方法可靠性研究[J]. 岩土力学,40(增刊1):565–573. Wang T Q,Wang Y Z,Yuan X M,Tang Z G,Wang H,Duan X F. 2019. Reliability analysis of acceleration integral displacement method based on shaking table tests[J]. Rock and Soil Mechanics,40(S1):565–573 (in Chinese).
谢俊举,温增平,高孟潭. 2013. 利用强震数据获取汶川地震近断层地面永久位移[J]. 地震学报,35(3):369–379. Xie J J,Wen Z P,Gao M T. 2013. Recovery of co-seismic deformation from strong motion records during the Wenchuan earthquake[J]. Acta Seismologica Sinica,35(3):369–379 (in Chinese).
张海明, 陈晓非. 2004. 根据数字强震加速度记录恢复位移波形的方法: 以1999年9月21日台湾集集地震为例[C]//中国大陆地震学与地球内部物理学研究进展: 庆贺曾融生院士八十寿辰. 北京: 地震出版社: 694–707. Zhang H M, Chen X F. 2004. The method of recovering displacement waveforms based on digital strong motion records: Take 1999·9·21 Chi-Chi earthquake for example[C]//Advances in Seismology and Physics of Earth’s Interior in China: In Honor of Academician Zeng Rongsheng’s 80th Birthday. Beijing: Seismological Press: 694–707 (in Chinese).
周宝峰,于海英,温瑞智,谢礼立. 2013. 一种识别永久位移的新方法[J]. 土木工程学报,46(增刊2):135–140. Zhou B F,Yu H Y,Wen R Z,Xie L L. 2013. A new way of permanent displacement identification[J]. China Civil Engineering Journal,46(S2):135–140 (in Chinese).
周雍年,章文波,于海英. 1997. 数字强震仪记录的长周期误差分析[J]. 地震工程与工程振动,17(2):1–9. Zhou Y N,Zhang W B,Yu H Y. 1997. Analysis of long-period for accelerograms recorded by digital seismographs[J]. Earthquake Engineering and Engineering Vibration,17(2):1–9 (in Chinese).
Baker J W. 2007. Quantitative classification of near-fault ground motions using wavelet analysis[J]. Bull Seismol Soc Am,97(5):1486–1501. doi: 10.1785/0120060255
Boore D M. 2001. Effect of baseline corrections on displacements and response spectra for several recordings of the 1999 Chi-Chi,Taiwan,earthquake[J]. Bull Seismol Soc Am,91(5):1199–1211.
Boore D M,Stephens C D,Joyner W B. 2002. Comments on baseline correction of digital strong-motion data:Examples from the 1999 Hector Mine,California,earthquake[J]. Bull Seismol Soc Am,92(4):1543–1560. doi: 10.1785/0120000926
Chanerley A A,Alexander N A. 2010. Obtaining estimates of the low-frequency ‘fling’,instrument tilts and displacement timeseries using wavelet decomposition[J]. Bull Earthq Eng,8(2):231–255. doi: 10.1007/s10518-009-9150-5
Chiu H C. 1997. Stable baseline correction of digital strong-motion data[J]. Bull Seismol Soc Am,87(4):932–944.
Hızal C,Turan G. 2017. Importance of static correction and damping in the analysis of a cable-stayed bridge subjected to displacement loading[J]. J Bridge Eng,22(6):04017009. doi: 10.1061/(ASCE)BE.1943-5592.0001041
Iwan W D,Moser M A,Peng C Y. 1985. Some observations on strong-motion earthquake measurement using a digital accelerograph[J]. Bull Seismol Soc Am,75(5):1225–1246.
Li Q,Wang R Q,Huang W F,Zheng G J. 2005. Method for morphological filtering in seismic data processing[J]. Petrol Sci,2(4):20–29.
Tsai H C. 1998. Modal superposition method for dynamic analysis of structures excited by prescribed support displacements[J]. Comput Struct,66(5):675–683. doi: 10.1016/S0045-7949(97)00108-9
Wang G Q,Boore D M,Igel H,Zhou X Y. 2003. Some observations on colocated and closely spaced strong ground-motion records of the 1999 Chi-Chi,Taiwan,earthquake[J]. Bull Seismol Soc Am,93(2):674–693. doi: 10.1785/0120020045
Xu G L,Zhang L X,Bai Y S,Sun H. 2020. Method for solving dynamic equilibrium equations based on displacement excitation[J]. Earthq Eng Eng Vib,19(2):423–433. doi: 10.1007/s11803-020-0571-0