Ranking the seismic input motion based on seismic response of soft soil tunnels
-
摘要: 为探讨“最不利地震动”概念在地下结构抗震设计中的适用性,以软土地铁区间隧道为对象建立相应的地层-结构动力分析模型。以直径变形率为分析指标,基于动力时程方法研究18条不同输入地震动作用下隧道结构动力响应的分布及差异性,得出基于隧道地震响应的输入地震动排序,并通过调幅手段对比分析了地面峰值加速度(PGA)和隧道埋深变化对隧道结构地震动响应排序的影响规律。最后,评价了不同输入地震动参数,包括峰值加速度、峰值速度、峰值位移、绝对累积速度(CAV)和阿里亚斯(Arias)强度(IA)与隧道地震响应之间的相关性。分析结果表明:① 随着PGA从0.5 m/s2增加到2 m/s2,地震动排序发生明显变化,并且不同输入地震动引起的隧道地震响应差异显著提高,最不利地震动引起的直径变形率与平均值的比值从1.1增加到1.9;② 隧道从浅埋到深埋的过程中,地震动排序结果基本保持不变;③ PGA为2 m/s2时,隧道地震响应与基岩面峰值速度(PBV)的相关性最好,相关系数达到0.94,其次是与基岩面峰值位移(PBD)和IA,相关系数分别为0.62和0.48,相关性最差的是基岩面峰值加速度(PBA)和CAV,相关系数仅为0.37和0.22。研究结论可为今后软土隧道的输入地震动选择提供科学依据。Abstract: The selection of input motion is of most importance in seismic design and analysis of underground structures. In recent years, the concept of “the most unfavorable ground motion” has been widely studied and applied in the field of seismic design for surface structures. However, there is no research on ranking the seismic input motion based on the seismic response of underground structures. To explore the applicability of the concept of “the most unfavorable ground motion” for the seismic design of underground structures, a circular tunnel in soft soil is taken as the research object, and the corresponding soil-structure dynamic analysis modelis established. Taking the ovaling deformation rate of the circular liner as the critical index of tunnel responses, the distribution as well as the difference of dynamic responses of the tunnel structure under excitations of 18 different input motions are studied with the dynamic time-history method, and thus the ranking of seismic input motions based on the tunnel response is obtained. By means of amplitude modulation, the effects of peak ground acceleration (PGA)and tunnel depth on the input motion ranking are further investigated. Finally , the correlation between different input motion parameters, including peak acceleration, peak velocity, peak displacement, absolute cumulative velocity (CAV) and Arias intensity (IA), and the tunnel seismic response is evaluated. Results show that, ① When PGA increases from 0.5 m/s2 to 2 m/s2, the ranking of ground motion changes apparently, and the diversity of tunnel responses caused by different input motions is significantly amplified, such that the ratio of the maximum diametrical deformation rate caused by the most unfavorable input motion to the mean value is increasing from 1.1 to 1.9; ② The input motion ranking results remain unchanged for different depths of the tunnel; ③ For the PGA of 2 m/s2, the correlation between the tunnel response and the peak bedrock velocity (PBV) matches well, i.e. the correlation coefficient is 0.94, then followed by the peak bedrock displacement (PBD) and the IA, i.e. the correlation coefficient is 0.62 and 0.48, respectively, and the worse is, the peak bedrock acceleration (PBA) and the CAV have a correlation coefficient of only 0.37 and 0.22, respectively. The conclusions provide a scientific basis for the selection of seismic input motions of soft soil tunnels.
-
引言
近年来地下结构震害频发,例如1999年土耳其迪兹杰(Düzce)地震导致博卢(Bolu)隧道区间发生倒塌(Kontoe et al,2008),2008年汶川地震中龙溪隧道部分区段发生完全坍塌(Yu et al,2016)。因此,地下结构抗震设计至关重要。
输入地震动的选取是地下结构抗震设计与分析面临的首要问题。目前,输入地震动的选取方法主要有以下几种(American Society of Civil Engineers,2013;中华人民共和国住房和城乡建设部,2014):根据场地信息和设防要求选取相似场地和地震危险性的实际地震记录;根据设计谱选用实际强震记录或人工合成地震动(刘帅等,2018);选择对结构最不利的地震动进行输入,即最不利地震动方法。
“最不利地震动”的概念最早由谢礼立和翟长海(2003)提出,是指在选择输入地震动进行设计、分析或试验研究时,在地面峰值加速度(peak ground acceleration,缩写为PGA)和场地类别均符合规范的前提下,选取能使结构趋于最不利状态的地震动。“最不利地震动”概念已经在地面建筑抗震设计领域中进行了广泛的研究,并得到了大量工程实例的验证,包括网壳结构、铁路梁桥和文物系统等等(范峰等,2003;樊圆,2018;胡进军等,2020,2021)。这些研究表明结构响应与PGA等地震动参数、绝对累积速度(cumulative absolute velocity,缩写为CAV)等能量参数和持时参数密切相关,并为相关结构的输入地震动选择提供了重要参考。
最不利地震动方法解决了目前地震动选取困难的问题,同时也为地震动对工程结构的破坏机理研究提供了新的途径。然而,目前最不利地震动方法研究主要针对地表结构,还未见关于地下结构地震动排序方面的研究。由于地下结构的地震响应机制和控制因素与地上结构存在本质区别,因此有必要分析与评价“最不利地震动”概念在地下结构抗震设计中的适用性。
本文拟以软土地区典型地铁区间隧道为例,采用动力时程方法研究基于隧道地震响应的最不利输入地震动排序,并对PGA和埋深变化对地震动排序的影响规律进行探讨。
1. 软土隧道模型
1.1 地铁区间隧道
假设隧道及地层满足平面应变条件,隧道结构地震响应以横断面剪切变形为主,从而建立二维有限元分析模型。选取典型区间隧道,横断面为直径6.0 m的圆形,管片厚度为0.3 m,采用C30混凝土,基本材料参数列于表1。隧道初始埋深设为6.0 m。
表 1 衬砌结构材料参数Table 1. Material parameters of lining structure密度/(kg·m−3) 弹性模量/(1010 N·m−2) 泊松比 混凝土C30 2 600 3 0.2 1.2 软土地层条件
设覆盖层厚度为70 m,土层假设为单层均质各向同性,最大剪切波速为200 m/s,重度为17.3 kN/m3,泊松比取为0.3。根据软土地区统计资料(张亚军等,2010)确定土层的动力特性参数及动剪切模量比G/Gmax和阻尼比ξ,其动力特性曲线如图1所示。
1.3 地层-结构动力分析模型
建立的地层-结构有限元模型如图2所示,其中土层两侧水平边界距隧道中心的距离取为100 m,覆盖层厚度为70 m。隧道结构采用梁单元模拟,可以考虑剪切、弯曲、拉压等作用;地层采用四结点双线性四边形单元;为简化分析,假设结构与周围土层之间完全黏结无滑移。为了考虑地震作用下土层非线性响应特征,采用等效线性化方法,基于自由场分析程序SHAKE91 (Idriss,Sun,1992)计算不同地震动输入下土层的动剪切模量G和阻尼比ξ。
计算过程分为两步:初始地应力分析和地震动力分析。在初始地应力分析步,模型侧向边界节点水平向和底部边界节点竖向均被约束,并施加重力;在地震动力分析步,模型侧向边界改为竖向约束,水平向设置等位移边界条件,且模型底部为黏弹性边界,地震动采用波动法进行输入。
1.4 输入地震动和分析指标
输入地震动基于PEER数据库中常用的强震记录,共18条地震动,如表2所示。为探讨地震动排序随PGA的变化规律,以PGA达到0.5,1和2 m/s2为目标,分别对应多遇、基本和罕遇三种地震动水准,采用SHAKE91通过反演方法实现模型底部的地震动输入。
表 2 地震动记录Table 2. Ground motion recordings地点 日期 台站 方向 Chi-Chi,Taiwan,China 1999-09-20 ILA004 N Chi-Chi,Taiwan,China 1999-09-20 ILA004 W Chi-Chi,Taiwan,China 1999-09-20 ILA044 W Niigata,Japan 2004-10-23 NIG014 EW Chuetsu-oki,Japan 2007-07-16 NIG014 NS Iwate,Japan 2008-06-13 MYG006 EW Tottori,Japan 2000-10-06 TTR008 EW El Mayor-Cucapah,Mexico 2010-04-04 El Centro Array #3 270° Christchurch,New Zealand 2011-02-21 Christchurch Resthaven S88°E Tottori,Japan 2000-10-06 SMN002 NS Iwate,Japan 2008-06-13 AKT015 EW El Mayor-Cucapah,Mexico 2010-04-04 El Centro Array #3 360° Chuetsu-oki,Japan 2007-07-16 NIG025 NS Iwate,Japan 2008-06-13 AKT016 EW Iwate,Japan 2008-06-13 AKT016 NS Imperial Valley-07 1979-10-15 El Centro Array #3 140° Iwate,Japan 2008-06-13 IWT020 EW Iwate,Japan 2008-06-13 IWT020 NS 注:表中方向数字代表地震动记录的方位角。 为分析隧道埋深变化对地震动排序的影响,计算条件保持不变,仅改变隧道的埋深(9,12,18和30 m),涵盖了隧道从浅埋(1倍隧道直径简称1D,下同)到深埋(5D)的变化过程。具体计算工况如表3所示。
表 3 计算工况Table 3. Computational cases编号 PGA/(m·s−2) 隧道埋深/m 隧道埋深直径比 工况1 0.5 6.0 1.0 工况2 1 6.0 1.0 工况3 2 6.0 1.0 工况4 2 9.0 1.5 工况5 2 12.0 2.0 工况6 2 18.0 3.0 工况7 2 30.0 5.0 本文以直径变形率作为评价隧道地震响应的分析指标,且计算结果仅分析由地震作用引起的隧道变形增量部分。因此,根据直径变形率即可实现对地震动响应的排序。
2. 计算结果分析
针对表3中的每个工况,计算可得到18条地震动输入下的隧道地震响应分布及差异性;通过对全部工况的对比分析,进而探究不同地震动强度和埋深变化对地震动排序的影响规律。
2.1 直径变形率累积分布函数
图3给出了隧道埋深为1D时,对应不同地震动强度的直径变形率累积分布函数。累积分布函数是概率密度函数的积分,对于所有实数x,累积分布函数定义为
$$ {F_X}(x) {\text{=}} P(X {\text{≤}} x) {\text{=}} \int_{ - \infty }^x {{f_X}(t){\rm{d}}t}{\text{,}} $$ (1) 式中,X为变量,P为概率,f为概率密度函数。
从图3中可知,随着地震动强度的增加,直径变形率的平均值逐渐增大,且不同输入地震动引起的隧道地震响应差异显著提高。当PGA为0.5 m/s2时,直径变形率平均值为0.21%,最大值仅为0.24%,最不利记录与最小记录的直径变形率几乎相同,仅相差0.05%;当PGA为2 m/s2时,直径变形率平均值为0.45%,最大值为0.87%,最不利记录的直径变形率达到了最小记录的3.5倍。随着PGA从0.5 m/s2增加到2 m/s2,最不利地震动引起的直径变形率与平均值的比值从1.1增加到1.9。
综上可知,输入地震动的选择对于小震工况(PGA=0.5 m/s2)影响不大,但大震工况(PGA=2 m/s2)必须面对输入地震动的选取问题,即选择对结构最不利的地震动进行输入。因此,最不利地震动方法成为解决强震工况下输入地震动选取难题的有效途径。
2.2 不同地震强度下地震动排序分布
为分析地震动强度变化对地震动排序的影响,采用颜色映射方法(胡进军等,2020)描述地震动排序变化情况。以工况1的排序结果为基准,红色对应直径变形率最大的地震动记录(简称最大响应记录),紫色对应直径变形率最小的地震动记录(简称最小响应记录),图4展示了对应不同地震动强度的颜色映射序列。
从图中可知,随着PGA的增加,相较于最小响应记录和最大响应记录,中间区段响应记录的排序变化更加显著。对于最小响应记录,当PGA为0.5 m/s2时,土层的地震响应较小,处于近似弹性状态或弱非线性段,因此排序没有随地震动强度发生明显变化;对于最大响应记录,当PGA为0.5 m/s2时,对应最不利地震记录的动剪切模量比G/Gmax已经达到0.69,且随着地震动强度的增加,地层响应的非线性特征更加显著,使得该记录下的隧道地震响应远远大于其它地震记录,因此最大响应记录的排序也没有发生明显变化。
对于中间区段响应记录,随着地震动强度的增加,排序的差异性变化更为显著。从而得出,即使对于同一隧道结构,不同地震水准下的最不利输入地震动排序仍是不同的,尤其是处于中间区段响应的地震动记录。
2.3 不同埋深条件下地震动排序分布
以PGA为2 m/s2为例,对应不同埋深的直径变形率累积分布函数如图5所示。从图中可知,随着埋深的增加,直径变形率的平均值逐渐增大,但不同输入地震动引起的隧道地震响应差异基本保持不变。
对应隧道埋深从1D,1.5D,2D,3D到5D的直径变形率平均值分别为0.47%,0.53%,0.60%,0.66%,0.68%。随着埋深的增加,直径变形率平均值的增长速率逐渐减小。隧道埋深从1D增加到3D时,即浅埋到深埋的转变,导致直径变形率的增长速率较大;当隧道埋深超过3D后,埋深的变化对直径变形率影响较小。
同样使用颜色映射的方法分析埋深变化对地震动排序的影响,以工况3的排序结果为基准,图6给出了对应不同埋深条件的不利地震动颜色映射序列。结果表明,隧道从浅埋到深埋过程中,地震动排序结果基本保持不变。
综合累积分布函数图和颜色映射图可知,当埋深变化时,不同输入地震动引起的直径变化率的改变几乎相同,因此直径变形率的相对大小关系基本保持不变,即在累积分布函数图上,表现为累积分布曲线的整体向右“平移”,而在颜色映射图上,则表现为基本一致的颜色映射序列。
3. 直径变形率与地震动参数相关性分析
为进一步分析地震动强度和埋深变化对地震动排序的影响规律,选取不同地震动参数包括峰值加速度、峰值速度、峰值位移、绝对累积速度(CAV)和阿里亚斯强度(IA)与直径变形率进行相关性分析。其中,后两个参数的计算公式分别为
$$ {\rm{CAV}} {\text{=}} \int_0^{{t_{\max }}} {|a(t)|{\rm{d}}t} {\text{,}}$$ (2) $$ I_{\rm{A}} {\text{=}} \frac{\pi }{{2g}}\int_0^t {{a^2}(t){\rm{d}}t} {\text{,}} $$ (3) 式中,tmax为地震动持时,a(t)为t时刻对应的地震动加速度,g为重力加速度。
由于不同输入地震动下的地层动剪切模量变化反映了地层的动力非线性特性,因此对动剪切模量与直径变形率进行相关性分析,以评价地层响应和隧道地震响应的相关性。
采用相关系数对上述参数和直径变形率进行相关性分析,对任意两个变量X和Y,相关系数的计算公式为
$$ r(X{\text{,}}Y) {\text{=}} \frac{{{\rm{Cov}}(X{\text{,}}Y)}}{{\sqrt {{\rm{Var}}[X]{\rm{Var}}[Y]} }} {\text{,}} $$ (4) 式中,Cov(X,Y)为X与Y的协方差,Var [ X ] 和Var [ Y ] 分别为X和Y的方差。
相关系数变化范围为1至−1,分别代表绝对正相关和绝对负相关。表4给出了不同工况下的相关系数,其中,工况1、工况2和工况3提供了埋深为1D时,三种地震动强度的对比关系;工况3、工况6和工况7提供了PGA为2 m/s2时,三种埋深条件下的对比分析。
表 4 直径变形率与地震动参数相关性Table 4. Correlation coefficients between diameter deformation rate and input motion parameters工况编号 PBA PBV PBD CAV IA 动剪切模量 工况1 −0.59 0.72 0.35 −0.27 −0.18 −0.99 工况2 −0.20 0.93 0.43 −0.01 0.22 −0.98 工况3 0.37 0.94 0.62 0.22 0.48 −0.94 工况6 0.39 0.94 0.61 0.23 0.50 −0.96 工况7 0.38 0.94 0.61 0.23 0.49 −0.96 结果表明,上述工况中直径变形率与动剪切模量的相关性均保持在较高的水平,且以绝对负相关为主,说明地层-结构动力相互作用是隧道地震响应的主要控制因素。对比工况3、工况6和工况7可见,地震动排序随埋深基本保持不变,这与2.3节分析结果基本一致,该现象可以解释为:在隧道埋深变化过程中,输入地震动保持不变,土层的动剪切模量也基本保持不变,由于动剪切模量与直径变形率的相关性较高,因此地震动排序结果基本保持不变。
此外,当PGA为0.5 m/s2和1 m/s2时,即表4中的工况1和工况2,直径变形率与输入地震动的峰值加速度,即基岩面峰值加速度(peak bedrock acceleration,缩写为PBA)均为负相关,其主要原因在于调幅处理降低了不同输入地震动加速度峰值的差异性,使得正相关性减小,以及隧道自身惯性力不是地震响应的主要控制因素,这正是地下结构与地面结构抗震的主要区别。
随着地震动强度的增加,隧道直径变形率与五种地震动参数之间的正相关性均逐渐提高。当PGA为2 m/s2时(表4中的工况3,6,7),隧道地震响应与基岩面峰值速度(peak bedrock velocity,缩写为PBV)相关性最好,相关系数达到0.94,其次是基岩面峰值位移(peak bedrock displacement,缩写为PBD)和IA,相关系数分别为0.62和0.48,相关性最差的是基岩面峰值加速度(peak bedrock acceleration,缩写为PBA)和CAV,相关系数仅为0.37和0.22。杜陆荣等(2021)计算了土坡在100条随机地震动记录作用下的响应,发现PGA,PGV和PGD与土坡变形位移具有较好的相关性,并且相关性从高到低依次为PGV,PGD和PGA,与本文计算结果相一致,说明了数值计算的可靠性和有效性。因此,对于强震工况,建议基于PBV参数的地震动排序作为输入地震动选择的依据。
4. 结论
本文旨在探讨“最不利地震动”概念在地下结构抗震设计中的适用性。以软土地铁区间隧道为对象建立地层-结构分析模型,采用动力时程方法,基于不同的输入地震动记录,以隧道直径变化率为响应分析指标,分析了隧道结构在不同地震动强度和埋深变化下的地震动排序规律,得出以下结论:
1) 随着PGA从0.5 m/s2增加到2 m/s2,地震动排序发生明显变化,并且不同的输入地震动引起的隧道响应差异显著提高,因此大震工况(PGA=2 m/s2)必须考虑输入地震动的选取问题,而最不利地震动方法为解决强震工况下输入地震动选取难题提供了有效途径;
2) 隧道从浅埋到深埋过程中,地震动排序结果基本保持不变,且不同输入地震动引起的隧道响应差异基本保持不变,因此隧道结构在选择输入地震动时,可不考虑埋深变化的影响;
3) PGA为2 m/s2时,隧道地震响应与PBV相关性最好,相关系数达到0.94,因此建议强震工况采用基于PBV参数的地震动排序作为输入地震动选择的依据。
本文研究结论可为今后软土隧道抗震设计的输入地震动选择提供科学依据,本方法还可以推广到其它类型地下结构和场地条件。
谢礼立院士曾多次到访同济大学讲授最不利地震动的概念,来庆辉博士为本研究提供了地震动记录和指导建议,作者在此一并表示感谢。
-
表 1 衬砌结构材料参数
Table 1 Material parameters of lining structure
密度/(kg·m−3) 弹性模量/(1010 N·m−2) 泊松比 混凝土C30 2 600 3 0.2 表 2 地震动记录
Table 2 Ground motion recordings
地点 日期 台站 方向 Chi-Chi,Taiwan,China 1999-09-20 ILA004 N Chi-Chi,Taiwan,China 1999-09-20 ILA004 W Chi-Chi,Taiwan,China 1999-09-20 ILA044 W Niigata,Japan 2004-10-23 NIG014 EW Chuetsu-oki,Japan 2007-07-16 NIG014 NS Iwate,Japan 2008-06-13 MYG006 EW Tottori,Japan 2000-10-06 TTR008 EW El Mayor-Cucapah,Mexico 2010-04-04 El Centro Array #3 270° Christchurch,New Zealand 2011-02-21 Christchurch Resthaven S88°E Tottori,Japan 2000-10-06 SMN002 NS Iwate,Japan 2008-06-13 AKT015 EW El Mayor-Cucapah,Mexico 2010-04-04 El Centro Array #3 360° Chuetsu-oki,Japan 2007-07-16 NIG025 NS Iwate,Japan 2008-06-13 AKT016 EW Iwate,Japan 2008-06-13 AKT016 NS Imperial Valley-07 1979-10-15 El Centro Array #3 140° Iwate,Japan 2008-06-13 IWT020 EW Iwate,Japan 2008-06-13 IWT020 NS 注:表中方向数字代表地震动记录的方位角。 表 3 计算工况
Table 3 Computational cases
编号 PGA/(m·s−2) 隧道埋深/m 隧道埋深直径比 工况1 0.5 6.0 1.0 工况2 1 6.0 1.0 工况3 2 6.0 1.0 工况4 2 9.0 1.5 工况5 2 12.0 2.0 工况6 2 18.0 3.0 工况7 2 30.0 5.0 表 4 直径变形率与地震动参数相关性
Table 4 Correlation coefficients between diameter deformation rate and input motion parameters
工况编号 PBA PBV PBD CAV IA 动剪切模量 工况1 −0.59 0.72 0.35 −0.27 −0.18 −0.99 工况2 −0.20 0.93 0.43 −0.01 0.22 −0.98 工况3 0.37 0.94 0.62 0.22 0.48 −0.94 工况6 0.39 0.94 0.61 0.23 0.50 −0.96 工况7 0.38 0.94 0.61 0.23 0.49 −0.96 -
杜陆荣,张江伟,迟明杰,陈苏. 2021. 地震动峰值特征参数对土坡地震响应的规律研究[J]. 地震学报,43(4):498–507. doi: 10.11939/jass.20200149 Du L R,Zhang J W,Chi M J,Chen S. 2021. Regularity research on the seismic response of characteristic parameters for ground motion peak to soil slope[J]. Acta Seismologica Sinica,43(4):498–507 (in Chinese).
范峰,钱宏亮,谢礼立. 2003. 最不利地震动在网壳结构抗震设计中的应用[J]. 世界地震工程,19(3):17–21. doi: 10.3969/j.issn.1007-6069.2003.03.004 Fan F,Qian H L,Xie L L. 2003. Applications of the severest ground motion to anti-seismic design for reticulated shells[J]. World Earthquake Engineering,19(3):17–21 (in Chinese).
樊圆. 2018. 基于破坏强度的典型桥梁结构地震动排序研究[D]. 哈尔滨: 中国地震局工程力学研究所: 51–70. Fan Y. 2018. Study on the Damage Potential Ordering of Design Ground Motion for Typical Bridge Structures[D]. Harbin: Institute of Engineering Mechanics: 51–70 (in Chinese).
胡进军,樊圆,来庆辉,谢礼立. 2020. 典型铁路梁桥抗震设计地震动破坏势排序[J]. 振动与冲击,39(11):14–21. Hu J J,Fan Y,Lai Q H,Xie L L. 2020. Ranking the damage potential of input ground motions for typical railway beam bridge[J]. Journal of Vibration and Shock,39(11):14–21 (in Chinese).
胡进军,来庆辉,梁琰,杨永强. 2021. 可移动文物隔震系统的地震动破坏强度排序[J]. 天津大学学报(自然科学与工程技术版),54(2):154–160. Hu J J,Lai Q H,Liang Y,Yang Y Q. 2021. Quantitative analysis of the damage potential of input ground motion based on movable cultural relic isolation systems[J]. Journal of Tianjin University:Science and Technology,54(2):154–160 (in Chinese).
刘帅,潘超,周志光. 2018. 对人造地震动反应谱求解及拟合的几个相关问题探讨[J]. 地震学报,40(4):519–530. Liu S,Pan C,Zhou Z G. 2018. Discussions on the response spectral solution and fitting of spectrum-compatible artificial seismic waves[J]. Acta Seismologica Sinica,40(4):519–530 (in Chinese).
谢礼立,翟长海. 2003. 最不利设计地震动研究[J]. 地震学报,25(3):250–261. doi: 10.3321/j.issn:0253-3782.2003.03.003 Xie L L,Zhai C H. 2003. Study on the most unfavorable design ground motion[J]. Acta Seismologica Sinica,25(3):250–261 (in Chinese).
张亚军,兰宏亮,崔永高. 2010. 上海地区土动剪切模量比和阻尼比的统计研究[J]. 世界地震工程,26(2):171–175. Zhang Y J,Lan H L,Cui Y G. 2010. Statistical studies on shear modulus ratios and damping ratios of soil in Shanghai area[J]. World Earthquake Engineering,26(2):171–175 (in Chinese).
中华人民共和国住房和城乡建设部. 2014. GB 50909—2014 城市轨道交通结构抗震设计规范[S]. 北京: 中国标准出版社: 24–27. Ministry of Housing and Urban-Rural Development of the People's Republic of China. 2014. GB 50909—2014 Code for Seismic Design of Urban Rail Transit Structures[S]. Beijing: Standards Press of China: 24–27 (in Chinese).
American Society of Civil Engineers. 2013. ASCE/SEI 7-10 Minimum Design Loads for Buildings and Other Structures[S]. Virginia: ASCE: 49–57.
Idriss I M, Sun J I. 1992. SHAKE91: A Computer Program for Conducting Equivalent Linear Seismic Response Analyses of Horizontally Layered Soil Deposits, User’S Guide[R]. California: University of California, Davis: 1–17.
Kontoe S,Zdravkovic L,Potts D M,Menkiti C O. 2008. Case study on seismic tunnel response[J]. Can Geotech J,45(12):1743–1764. doi: 10.1139/T08-087
Yu H T,Chen J T,Bobet A,Yuan Y. 2016. Damage observation and assessment of the Longxi tunnel during the Wenchuan earthquake[J]. Tunnell Underground Space Technol,54:102–116. doi: 10.1016/j.tust.2016.02.008
-
期刊类型引用(10)
1. 刘奕均,杨光亮,王嘉沛,谈洪波,周怀斌,申重阳. 基于粒子群算法的四川长宁地区印支构造界面反演. 地震地质. 2023(01): 172-189 . 百度学术
2. 罗嘉琦,帅向华,李继赓. 基于深度学习的倾斜摄影建筑物表面损毁信息提取. 中国地震. 2023(02): 271-281 . 百度学术
3. 龚丽文,张怀,陈石,David A.Yuen,陈丽娟,Brennan Brunsvik,尹光耀. 基于机器学习的长宁地震三维断层面几何特征建模. 地震学报. 2023(06): 1040-1054 . 本站查看
4. 扈桂让,张广伟,梁姗姗. 四川长宁地区地震活动空间分布特征研究. 中国地震. 2023(04): 883-892 . 百度学术
5. 李洪丽,刘财,田有,闫冬,李志强. 地震走时层析成像方法在研究强震孕震、发生机制中的应用. 地球与行星物理论评. 2022(02): 148-158 . 百度学术
6. 于书媛,方良好,宴金旭,倪红玉,陈靓,丁娟. 运用InSAR技术解算四川长宁M_S6.0地震三维形变场及时间序列形变分析. 地震工程学报. 2022(06): 1469-1477+1488 . 百度学术
7. 李昌珑,温增平. 2020年地震危险性分析热点回眸. 科技导报. 2021(01): 97-104 . 百度学术
8. 赵韬,王莹,马冀,邵若潼,徐一斐,胡景. 2021年青海玛多7.4级地震序列重定位和震源机制特征. 地震地质. 2021(04): 790-805 . 百度学术
9. 路珍,郭志,黄伟亮,张佩,李君,张媛媛. 2009年陕西高陵M_L4.8地震震源参数及构造特征分析. 大地测量与地球动力学. 2021(11): 1157-1162 . 百度学术
10. 王月,胡少乾,何骁慧,郭凯,解孟雨,邓世广,马亚伟. 2021年5月21日云南漾濞6.4级地震序列重定位及震源机制研究. 地球物理学报. 2021(12): 4510-4525 . 百度学术
其他类型引用(3)