Analysis of Influence Characteristics of Site Ground Motion in Ya'an Area, Sichuan
-
摘要: 为分析雅安地区场地结构引起的震动特性,选取该地区重大安评项目中的实际测试钻孔进行场地地震反应分析,采用等效线性化方法,并以NGA强震记录为输入,开展场地反应计算,在全面分析该地区场地地震特征的基础上,以峰值速度为参数指标,研究场地条件引起的峰值速度放大效应。研究结果表明,雅安地区场地卓越周期约为0.2 s,自振周期为0.2 s的建筑结构破坏相对严重;随着基岩峰值加速度的增大,峰值速度呈增大趋势,峰值速度与峰值加速度变化趋势并非完全一致,分析场地效应时应考虑多种地震动参数的变化规律;不同地震输入下峰值速度放大系数为1~1.8;峰值速度放大系数kv在不同强度基岩输入下,大致呈随等效剪切波速的增大而减小的趋势,软土场地更易引起场地峰值速度的放大,这种现象在芦山地震中较为明显。Abstract: In order to analyze the vibration characteristics caused by the site structure in the Ya'an area, the actual test boreholes in the major safety assessment projects in the area were selected for site seismic response analysis, The equivalent linearization method was used to carry out the site response calculation with the NGA strong earthquake records as the input. On the basis of a comprehensive analysis of the seismic characteristics of the site in this area, the peak velocity was used as the parameter index to study the magnification effect of the peak velocity caused by the site conditions. The main research results are as follows: The predominant period of the site in the Ya'an area is around the range of 0.2 s, and the building structure damage will be relatively severe if the natural vibration period is in this range; With the increase of the peak acceleration of the bedrock, the peak velocity shows an increasing trend, However, the changes in the peak velocity and the peak acceleration are not all consistent, so the variation laws of various ground motion parameters should be considered when analyzing the site effect.;The peak velocity magnification factor is between 1 and 1.5 under different seismic inputs; The site coefficient kv generally decreases with the increase of the equivalent shear wave velocity and soft soil filed is more likely to cause amplification of the peak velocity, especially in Lushan earthquake.
-
引言
21世纪初以来,伴随着城市化进程的不断加快,城市中可以使用的空间越来越少,通过建设和开发地下空间工程来提高城市空间利用率已经成为世界性的发展趋势。地下工程包括地下综合管廊、地下停车场、轨道交通、地下污水处理厂、地下商业综合体等(钱七虎,2017),其对于解决城市化进程中出现的土地紧张、环境污染、交通拥堵和能源浪费等问题都发挥着积极作用(陈晓强等,2010)。然而,地下结构的地震安全性并没有引起人们的重视,以往学者们普遍认为地下结构完全被土体所包围,地震时地下结构比地面结构安全(陈国兴等,2016),因而地下结构的抗震设计并没有受到充分的重视,这也导致地下结构抗震研究一直停滞不前。在1995年日本阪神大地震中,地下结构发生了严重破坏(Iida等,1996),科研人员才开始真正重视地下结构的抗震问题。
一直以来,中国缺乏专门针对地下结构进行抗震设计的统一规范,各种规范中所考虑的方法与参数各不相同。早期,大部分地下结构抗震设计沿用地面结构的静力设计方法,比如《铁路工程抗震设计规范(GB 50111—2006)》(中华人民共和国建设部,2006),就是用地面结构抗震的思想进行地下结构设计,但地下结构在地震过程中受地基土约束,变形也受土体控制(林皋, 1990a, 1990b),在抗震原理上与地面结构大不相同(Hashash等,2001)。
美国和日本对地下结构抗震开展研究较早,基于结构变形受土体变形控制这一核心思想,率先提出了很多实用的关于地下结构抗震简化设计方法(权登州等,2015),包括自由场变形法(Wang,1993)、反应位移法(川岛一彦,1994)和柔度系数法(Penzien,2000)。其中以反应位移法推广度最高,但在反应位移法的使用过程中,地基弹簧系数直接决定了计算结果是否合理,弹簧系数的取值需要通过6次有限元计算才能获得,计算过程比较繁琐且不准确。在规避弹簧系数的取值问题上,国内外诸多学者提出了改进方法,如片山幾夫等(1985)提出了反应加速度法;Tateishi(2005)提出了强制反应位移法和改进方法;刘如山等(2007)提出了有限元反应应力法,即将一维土层剪应力沿竖向微分,再将其作为水平体荷载离散到有限元节点上;刘晶波等(2008, 2013)提出了地下结构Pushover分析法和整体式反应位移法,前者借鉴地面结构Pushover方法,施加倒三角形荷载同时推动带结构的土体和自由场土体,通过自由场土体变形控制位移,后者直接在挖除结构的土-结构边界上施加位移,再把反力施加到带有结构的土体。这些简化分析方法为地下结构抗震研究提供了思路,但其适用性和精度还有待进一步检验。
本文以有限元软件ABAQUS为计算平台,以1个双跨箱形地下结构为算例,通过反应位移法、强制反应位移法、反应加速度法和仅在土体施加加速度4种方法进行对比计算,再以动力时程分析的结果为参考,分析4种简化分析方法所得弯矩、剪力和相对位移的计算精度,验证4种简化分析方法的适用性与准确性。
1. 常用地下结构抗震分析方法
1.1 动力时程分析法
动力时程分析方法以有限元软件为计算平台,建立土-结构相互作用模型。其基本模型是将土分割为二维平面应变有限元,采用梁单元模拟地下结构并与土体连接,仅考虑由场地底部往上传播的剪切波,侧边界采用固定竖向、释放水平向的方式,底边界则固定竖向,再从模型的底部输入地震动,进而分析得到水平地震作用下土-地下结构的地震反应规律,计算模型如图 1所示。动力时程分析法是采用逐步积分的方式求解地下结构在地震荷载作用下每一时刻的动力响应,能够较好地考虑材料的非线性、地震动频谱、持时和幅值情况。动力时程分析方法能够较准确地反应地下结构的地震反应,该方法常被作为其它简化分析方法的校核标准,然而其对土体动力本构关系要求较高,计算非线性问题耗时较长,求解难度较大,因而对使用人员的专业知识要求较高,同时还需要具备较好的计算平台,通常只被应用于重大工程项目的抗震分析中。国内地下建筑抗震设计规范《地下铁道建筑结构抗震设计规范(DG/TJ 08-2064-2009)》(同济大学,2009)和《城市轨道交通结构抗震设计规范(GB 50909—2014)》(中华人民共和国住房和城乡建设部,2014)都引入了动力时程分析法。
1.2 反应位移法
反应位移法采用梁-弹簧单元建模(不包含土体),其中梁单元模拟地下结构,弹簧单元(正向和切向)模拟结构和地基土之间的相互作用。首先,计算出结构顶、底板位置处所对应土层最大位移差,再把此时刻对应的位移差施加在弹簧远离结构的一端,同时在自由场计算中提取出该时刻对应的结构惯性力和地层剪应力,并直接施加在梁单元模型上,采用静力的方法计算结构响应,其模型见图 2。反应位移法的计算结果依赖于地基弹簧刚度的取值,通过有限元计算弹簧刚度的方法也比较繁琐,且其取值是否合理也有待商榷。同样,《地下铁道建筑结构抗震设计规范(DG/TJ 08-2064-2009)》和《城市轨道交通结构抗震设计规范(GB 50909—2014)》都借鉴了反应位移法。
1.3 强制反应位移法
土体边界施加强制变形是反应位移法的1种简化形式,进行简化主要是为了规避弹簧刚度的取值问题,计算过程和反应位移法相同。首先,需要通过自由场计算得到地震荷载下结构顶板和底板所处土层的最大位移差,此时结构变形最大,即处于极限响应状态,然后将该时刻自由场计算模型中全部土层的相对位移提取出来,并以强制位移的模式施加在有限元模型侧边界上,将动力问题简化为静力问题,从而得到结构响应,其计算模型如图 3所示,其中W为地下结构宽度,D为地下结构到有限元模型侧边的距离。Tateishi(2005)指出,采用强制反应位移法处理动力作用虽然能够满足自由场应变在静力模型边界及附近准确传递,但远离模型边界处会因为土体阻尼的影响而发生衰减现象,尤以结构周围土层最为明显,应变值可衰减20%—50%。由于地下结构刚度通常大于土体刚度,因而结构变形比自由场变形小,如果选取的侧边距过小,土体阻尼作用将无法得到发挥,则相当于把自由场变形直接施加到了地下结构上,从而导致地下结构响应过大。由此可见,找到合适的侧边距是强制反应位移法准确求解结构响应的关键。
1.4 反应加速度法
反应位移法也可以通过反应加速度法进行简化,与强制反应位移法相同,均是为了规避弹簧刚度取值失真产生较大的计算误差。首先,通过自由场计算给出地震荷载下结构顶板和底板所处土层的最大位移差,然后把该时刻模型中全部土层的加速度取出,最后将该加速度结果通过水平体力的方式赋予整个模型土体和地下结构,通过静力的方法给出结构响应,其计算模型如图 4所示。近年来,反应加速度法在国内得以大力推广,《城市轨道交通结构抗震设计规范(GB 50909—2014)》首次引入了反应加速度法。
1.5 仅在土地施加加速度的简化方法
Tateishi(2005)提出了强制反应位移法的改进方法,对整个土层强制施加水平体力从而让土层变形与自由场变形一致,即仅把水平加速度按体力的形式施加到土体上。仅在土体施加水平加速度也是对反应加速度法的进一步简化,计算求解的过程和反应加速度法一致,唯一的区别在于反应加速度法把最大变形状态时的土层加速度以体力的形式施加到土体和结构上,而该方法则进一步简化了这一处理过程,其模型见图 5。由于地下结构一般内部挖空,仅有周围一圈衬砌,纵向单位长度的隧道质量远远小于原土体的质量,加速度在结构上产生的体力较小,因此理论上对结构的地震响应影响甚微(杨智勇等,2012)。
2. 算例分析
2.1 模型建立
以某单一土层中埋深7m的双跨箱形地下结构为算例,计算其在地震荷载下的响应,该地下结构横截面尺寸为8m×4m,墙厚0.8m,中柱厚0.8m。依次采用动力时程法、反应位移法、强制反应位移法、反应加速度法和仅在土体施加加速度这5种方法计算该地下结构的内力和中柱顶底相对变形。本例中结构的横截面抗震分析可按平面应变问题处理,土-结构接触面采用ABAQUS软件中自带的绑定接触,土体与结构参数见表 1。使用实体单元对土体进行建模,地下结构则使用梁单元进行建模并与土体连接,土体选取摩尔-库伦本构模型,结构则取为线弹性本构模型。
表 1 土体与结构的物理力学参数Table 1. Physical and mechanical parameters of soil and structure类型 密度/g·cm-3 弹性模量/MPa 泊松比 内摩擦角/° 粘聚力/kPa 土体 1.9 30 0.3 30 30 结构 2.3 30000 0.167 采用动力时程分析法时,分别将结构宽度和高度的4倍取为侧边距和底边距,计算结果用作校核,模型见图 6、图 7。反应位移法不建立土体,改变模型尺寸进行强制反应位移法、反应加速度法和仅在土体施加加速度方法的试算。
采用动力时程分析法时,侧边界固定竖向、水平向释放,底边界为粘性边界,从模型底部输入峰值加速度为0.2g的El-Centro波,如图 8所示。然后,再采用峰值加速度为0.2g的Kobe波(图 9)验证分析所得的规律。自由场模型与土-结构动力时程分析模型采用相同的土体本构模型和土体参数,0.2g峰值加速度下土层位移差最大时刻的相对位移分布和加速度分布见图 10、图 11。
2.2 结构响应分析
2.2.1 El-Centro波下不同侧边距结果分析
首先,以计算出1的动力时程分析法结果作为后续简化分析方法的校核标准。根据《城市轨道交通结构抗震设计规范(GB 50909—2014)》要求,动力时程分析的侧边距应大于结构宽度的3倍及以上,本次计算模型侧边距取为结构宽度的4倍;其次,采用反应位移法单独建立梁-弹簧单元对结构响应进行分析;最后,改变模型土体尺寸,用3种需要建立土体的简化分析方法对结构响应进行计算。输入峰值加速度为0.2g的El-Centro波,则控制截面A、B、C、D的弯矩、剪力和顶、底相对位移计算结果见表 2、表 3和表 4,地下结构计算误差分析见图 12。
表 2 不同侧边距的弯矩计算结果Table 2. Bending moment of structure with different lateral marginsD/W 使用方法 弯矩/kN·m 截面A 截面B 截面C 截面D 1倍 强制反应位移法 292.9 -291.3 280.8 353.9 反应加速度法 319.4 -324.5 294.0 367.1 仅在土体施加加速度 296.5 -300.6 278.4 351.8 2倍 强制反应位移法 179.5 -181.3 200.7 282.1 反应加速度法 328.1 -33.3 301.0 372.2 仅在土体施加加速度 309.3 -313.6 287.9 360.0 3倍 强制反应位移法 100.4 -101.3 147.7 224.7 反应加速度法 331.7 -336.9 303.8 374.9 仅在土体施加加速度 314.5 -18.9 291.8 363.9 4倍 强制反应位移法 50.2 -50.4 113.8 187.4 反应加速度法 332.6 -37.8 304.3 375.5 仅在土体施加加速度 316.2 -20.6 292.8 365.0 反应位移法 286.5 -290.1 201.0 199.8 动力分析 282.5 -89.1 272.3 345.2 表 3 不同侧边距的剪力计算结果Table 3. Shear force of structure with different lateral marginsD/W 使用方法 剪力/kN 截面A 截面B 截面C 截面D 1倍 强制反应位移法 -182.6 -182.6 -213.5 -263.8 反应加速度法 -195.8 -204.0 -205.6 -271.7 仅在土体施加加速度 -186.6 -186.6 -196.8 -266.3 2倍 强制反应位移法 -112.8 -112.8 -169.0 -249.3 反应加速度法 -201.3 -209.4 -210.0 -274.3 仅在土体施加加速度 -194.7 -194.7 -202.0 -270.7 3倍 强制反应位移法 -63.0 -63.0 -144.4 -219.7 反应加速度法 -203.6 -211.7 -211.9 -276.2 仅在土体施加加速度 -198.0 -198.0 -204.2 -273.3 4倍 强制反应位移法 -31.4 -31.4 -128.2 -199.6 反应加速度法 -204.1 -212.3 -212.0 -276.4 仅在土体施加加速度 -199.0 -199.0 -204.5 -273.8 反应位移法 -174.8 -183.0 -106.2 -106.2 动力分析 -175.7 -180.6 -194.0 -266.7 表 4 不同侧边距的结构变形计算结果Table 4. Deformation of structure with different lateral marginsD/W 使用方法 中柱顶、底相对位移/mm 1倍 强制反应位移法 2.396 反应加速度法 5.011 仅在土体施加加速度 4.636 2倍 强制反应位移法 0.302 反应加速度法 5.340 仅在土体施加加速度 5.026 3倍 强制反应位移法 0.584 反应加速度法 5.354 仅在土体施加加速度 5.071 4倍 强制反应位移法 0.722 反应加速度法 5.346 仅在土体施加加速度 5.077 反应位移法 3.628 动力分析 5.596 在输入0.2g峰值加速度El-Centro波的情况下,反应位移法的计算精度在结构中柱部分较好。对比反应位移法和动力时程分析方法可知,弯矩误差在中柱上端为1.41%,下端最小仅为0.34%,剪力误差在中柱上端最小仅为0.51%,下端为1.32%。但是在结构的角部,计算误差较大,剪力和弯矩计算值都偏小,在C截面(右侧墙上端)弯矩误差值达到26.18%,剪力误差值达到45.25%,尤其D截面(左侧墙下端),弯矩误差值最大达到42.12%,剪力误差值最大达到60.17%。弯矩平均误差为13.8%,剪力平均误差23.3%。
强制反应位移法的计算精度依赖于静力模型土体侧边距的选取,随着土体侧边距取值的不断增加,地下结构的弯矩、剪力计算误差会相应的增大。当侧边距取为4倍结构宽度时,弯矩平均误差最大达到70.2%,剪力平均误差最大达到58.1%;然而,当侧边距取为1倍结构宽度时,弯矩平均误差仅为2.8%,剪力平均误差为2.5%。当侧边距为1倍结构宽度时,中柱顶、底相对变形误差为57.18%,侧边距为2倍结构宽度时误差最大为94.60%,但随着侧边距的增大,误差有所减小。
反应加速度法与仅在土体施加加速度的简化分析方法计算所得弯矩、剪力和中柱变形随着模型侧边距的增大,计算误差变化幅度很小。当侧边距取为1倍结构宽度时,反应加速度法的弯矩平均误差为10.1%,剪力平均误差为8.7%,中柱变形误差为10.5%,仅在土体施加加速度方法的弯矩平均误差为3.1%,剪力平均误差为2.4%,中柱变形误差为17.2%;当侧边距取为4倍结构宽度时,反应加速度法的弯矩平均误差为14.3%,剪力平均误差为12.7%,中柱变形误差为4.5%,仅在土体施加加速度方法的弯矩平均误差为9.2%,剪力平均误差为7.8%,中柱变形误差为9.3%。由此可见,以加速度模拟地震作用时土体侧边距对计算误差影响较小,2种方法变化趋势基本同步。由于仅在土体施加加速度的简化分析方法忽略了结构加速度,相同工况下该方法的计算数值略小于反应加速度法。施加加速度的2种方法对于结构变形的计算精度均较高,反应位移法次之,强制反应位移法的计算结果较不理想。
2.2.2 Kobe地震波验证分析
为了消除地震动不同带来的偶然性,改用0.2g峰值加速度的Kobe地震波验证上述分析规律。
从前文分析可知,反应位移法在应用过程中不需要建立土体;强制反应位移法在侧边距取为1倍结构宽度时计算误差较小;反应加速度法和仅在土体施加加速度的简化方法对于侧边界取值不敏感。因此,需要建立土体的3种静力模型统一采用1倍结构宽度作为侧边距。弯矩、剪力、中柱变形计算结果见表 5、表 6、表 7。
表 5 Kobe波下弯矩计算结果Table 5. Bending moment calculation results of structure under Kobe seismic wave使用方法 弯矩/kN·m 平均误差/% 截面A 截面B 截面C 截面D 动力分析 375.3 -383.7 332.1 410.6 反应位移法 400.2 -403.7 284.3 276.0 10.5 强制反应位移法 298.8 -310.4 275.9 364.3 18.0 反应加速度法 341.9 -349.2 308.4 384.7 8.0 仅在土体施加加速度 312.4 -317.7 288.0 364.5 15.0 表 6 Kobe波下剪力计算结果Table 6. Shear force calculation results of structure under Kobe seismic wave使用方法 剪力/kN 平均误差/% 截面A 截面B 截面C 截面D 动力分析 -230.8 -242.9 -218.3 -299.0 反应位移法 -242.6 -257.5 -159.9 -139.9 13.7 强制反应位移法 -190.4 -190.4 -180.2 -290.7 27.0 反应加速度法 -207.4 -222.2 -209.8 -281.2 4.4 仅在土体施加加速度 -196.9 -196.9 -198.1 -275.4 5.6 表 7 Kobe波下结构变形计算结果Table 7. Deformation calculation results of structure under Kobe seismic wave使用方法 中柱顶、底相对位移/mm 平均误差/% 动力分析 5.808 反应位移法 4.768 17.9 强制反应位移法 3.232 44.4 反应加速度法 5.371 7.5 仅在土体施加加速度 4.900 15.6 对于不同类型的地震动,4种简化分析方法的计算精度均比较稳定,弯矩误差在10%左右,反应加速度法平均弯矩误差最小,为8.0%,强制反应位移法平均弯矩误差最大,为18.0%;2种施加加速度简化分析方法的剪力误差极小,分别为4.4%和5.6%,强制反应位移法剪力误差最大,仅为27%;强制反应位移法的结构变形误差最大达到44.4%,其余3种方法精度较好,反应加速度法误差最小,仅为7.5%。
3. 4种简化分析方法的特点分析
4种简化分析方法的相同点在于脱离了动力分析,能够按静力的方式计算结构响应,使得计算效率大大提高,且在不同地震动输入条件下都能得到较好的精度;另外,4种简化分析方法均需要进行自由场动力时程分析,从而获取地层位移或者地层加速度作为静力计算的地震输入荷载,自由场变形是4种简化分析方法的基础,但是结构存在对自由场的影响却被忽略了。
4种简化分析方法的不同点在于,①采用反应位移法不需要建立模型土体,以弹簧单元模拟土-结构相互作用,但是地基弹簧刚度取值需要进行6次有限元计算,过程比较繁琐。并且,土体弹簧刚度值对结构内力计算有巨大影响,通过有限元计算弹簧刚度的方法无法反应土层自身的相互作用,计算过程中土体的等效弹簧之间是不相关的,造成的结果则是土层对结构特别是角部的约束作用减小,计算所得的角部内力值偏小。②强制反应位移法对静力模型侧边界宽度值的选取非常敏感,侧边距的选取对于结构内力计算精度有重大影响。从计算结果来看,采用1倍结构宽度作为侧边距时结构内力的计算误差最小。但是无论侧边距如何变化,对结构变形的计算误差均较大。在对结构内力进行校核时,采用该方法计算的结果较为可靠。③反应加速度法和仅在土层施加加速度的简化方法对于侧边界的选取不敏感,2种方法对于不同的侧边距得到的计算精度均比较稳定,计算精度随着侧边距的增加仅有轻微的变化,2种方法对结构内力、变形的计算误差均较小。
4. 结论
本文通过反应位移法、强制反应位移法、反应加速度法和仅在土体施加加速度4种简化分析方法,分别计算了1个双跨箱形地下结构在0.2g峰值加速度的EL-Centro波和Kobe波作用下的地震响应,再利用时程分析法进行校核,经过对比,得出以下结论:
(1)总体而言,4种简化分析方法计算得到的结构在地震作用下的内力、变形值接近,规律一致,不同类型的地震动作用下均能得到较好的计算精度,能合理地反映地下结构的动力特性。
(2)若要进一步推广使用反应位移法,则须使其基地弹簧系数的取值快捷、可靠。
(3)模拟地震等效荷载的自由场土层位移和土层加速度具有较好的计算精度,且便于施加,但是需要进一步考虑地下结构存在对自由场的影响。
(4)在使用强制反应位移法进行结构内力校核时,侧边距取为1倍结构宽度时效果最好,计算误差最小。而反应加速度法和仅在土体施加加速度的方法对侧边距的选取不敏感。
(5)强制反应位移法、反应加速度法和仅在土体施加加速度3种带土体模型的简化分析方法均省略了反应位移法中地基弹簧刚度系数的取值,同时更加真实地反应了土-结构间的协调作用,使用时比反应位移法更简便。
-
表 1 场地钻孔资料分布情况
Table 1. Site borehole data distribution map
参数分布 等效剪切波速 ${V}_{\mathrm{s}\mathrm{e} }$分布/m·s−1 覆盖层厚度 $ {d}_{\mathrm{s}\mathrm{e}} $分布/m 200~250 250~300 300~350 350~400 400~450 10~15 15~20 20~25 25~30 30~35 钻孔数/个 2 2 9 10 4 4 4 7 9 3 表 2 输入基岩时程信息(芦山地震强震数据)
Table 2. Input bedrock time-histories information(Data from Lushan earthquake)
地震名称 台站名称 记录编号(简称) 震级/级 震源距离/km PGA/Gal PGV/cm·s−1 芦山7.0级地震 51 TQL 051 TQL130420080202(TQL-0202) 7.0 13 280.7 12.62 51 PJD 051 PJD130420080201(PJD-0201) 7.0 13 153.6 9.87 51 LDG 051 LDG130420080202(LDG-0202) 7.0 13 86.8 5.73 51 HYT 051 HYT130420080201(HYT-0201) 7.0 13 70.7 4.39 51 MNC 051 MNC130420080202(MNC-0202) 7.0 13 25.3 1.39 51 PXZ 051 PXZ130420080202(PXZ-0202) 7.0 13 12.3 1.71 表 3 输入基岩时程信息(NGA强震数据)
Table 3. Input bedrock time-histories information(Data from NGA)
地震名称 台站名称 记录编号 震级/级 震源距离/km PGA/Gal PGV/cm·s−1 Northridge-01 LA - Griffith Park Observatory 0994 6.69 17.5 283.4 26.51 Northridge-01 Vasquez Rocks Park 1091 6.69 17.5 148.0 18.39 Northridge-01 Lake Hughes #4 - Camp Mend 1021 6.69 17.5 82.4 6.17 Loma Prieta Point Bonita 0789 6.93 17.48 71.0 12.86 Chi-Chi, Taiwan-04 TTN042 2929 6.20 18 25.1 1.63 Chi-Chi, Taiwan-06 ILA063 3390 6.30 16 12.2 0.96 -
薄景山, 李秀领, 刘德东等, 2003 a. 土层结构对反应谱特征周期的影响. 地震工程与工程振动, 23(5): 42—45Bo J. S. , Li X. L. , Liu D. D. , et al. , 2003 a. Effects of soil layer construction on characteristic periods of response spectra. Earthquake Engineering and Engineering Vibration, 23(5): 42—45. (in Chinese) 薄景山, 李秀领, 刘红帅, 2003 b. 土层结构对地表加速度峰值的影响. 地震工程与工程振动, 23(3): 35—40Bo J. S. , Li X. L. , Liu H. S. , 2003 b. Effects of soil layer construction on peak accelerations of ground motions. Earthquake Engineering and Engineering Vibration, 23(3): 35—40. (in Chinese) 陈波, 2013. 结构非线性动力分析中地震动记录的选择和调整方法研究. 北京: 中国地震局地球物理研究所.Chen B., 2013. Ground motion selection and modification methods for performing nonlinear dynamic analysis of buildings. Beijing: Institute of Geophysics, China Earthquake Administration. (in Chinese) 陈党民, 段蕊, 2013. 场地条件对地表峰值加速度的放大效应分析. 地震工程与工程振动, 33(1): 24—30Chen D. M. , Duan R. , 2013. Amplification effects of site conditions on ground peak accelerations. Journal of Earthquake Engineering and Engineering Vibration, 33(1): 24—30. (in Chinese) 冯希杰, 金学申, 2001. 场地土对基岩峰值加速度放大效应分析. 工程地质学报, 9(4): 385—388Feng J. X. , Jin X. S. , 2001. Analysis on amplification effect of bedrock peak acceleration caused by site soil. Journal of Engineering Geology, 9(4): 385—388. (in Chinese) 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局, 2010. GB 50011—2010 建筑抗震设计规范. 北京: 中国建筑工业出版社.Ministry of Housing and Urban-Rural Development of the People's Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, 2010. GB 50011—2010 Code for seismic design of buildings. Beijing: China Architecture & Building Press. (in Chinese) Akkar S. , Özen Ö. , 2005. Effect of peak ground velocity on deformation demands for SDOF systems. Earthquake Engineering & Structural Dynamics, 34(13): 1551—1571. O’Rourke T. D. , Stewart H. E. , Jeon S. S. , 2001. Geotechnical aspects of lifeline engineering. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering, 149(1): 13—26. doi: 10.1680/geng.2001.149.1.13 -