• ISSN 1673-5722
  • CN 11-5429/P

地下结构地震反应的平面计算模型研究

殷琳 舒恩 楼梦麟 蔡海兵

殷琳,舒恩,楼梦麟,蔡海兵,2022. 地下结构地震反应的平面计算模型研究. 震灾防御技术,17(4):727−734. doi:10.11899/zzfy20220414. doi: 10.11899/zzfy20220414
引用本文: 殷琳,舒恩,楼梦麟,蔡海兵,2022. 地下结构地震反应的平面计算模型研究. 震灾防御技术,17(4):727−734. doi:10.11899/zzfy20220414. doi: 10.11899/zzfy20220414
Yin Lin, Shu En, Lou Menglin, Cai Haibing. Discuss on Plain Strain Model for Seismic Response of Underground Structure[J]. Technology for Earthquake Disaster Prevention, 2022, 17(4): 727-734. doi: 10.11899/zzfy20220414
Citation: Yin Lin, Shu En, Lou Menglin, Cai Haibing. Discuss on Plain Strain Model for Seismic Response of Underground Structure[J]. Technology for Earthquake Disaster Prevention, 2022, 17(4): 727-734. doi: 10.11899/zzfy20220414

地下结构地震反应的平面计算模型研究

doi: 10.11899/zzfy20220414
基金项目: 国家自然科学基金(91315301);安徽省2019年度高校学科(专业)拔尖人才学术资助项目
详细信息
    作者简介:

    殷琳,男,生于1982年。讲师。主要从事工程结构抗震防灾等研究。E-mail:315989441@qq.com

    通讯作者:

    楼梦麟,男,生于1947年。教授。主要从事工程结构抗震防灾等研究。E-mail: 93010@tongji.edu.cn

Discuss on Plain Strain Model for Seismic Response of Underground Structure

  • 摘要: 采用平面应变模型对地下结构进行地震反应分析时,其核心问题是中柱的二维等效简化。常用的简化方法是将中柱的材料性质(如弹性模量和密度)进行折减。在此基础上,进一步引入空间约束影响系数和三维还原系数,提出新的中柱二维等效简化方法。针对不同简化方法,分别建立对应的地下结构地震反应分析平面应变模型,计算各模型的地震反应。通过与三维模型计算结果进行对比分析,研究不同简化方法的合理性。计算结果表明,本研究建议的方法可有效提高地下结构平面应变模型的计算精度。
  • 有限元法是地下结构地震反应分析的常用方法(刘晶波等,2006杜修力等,2017陈国兴等,2019),采用有限元法对地下结构进行地震反应分析时,通常将三维模型简化为平面应变模型进行计算(王国波等,2015庄海洋等,2019)这是因为地震时,地下结构横断面变形以平面内的剪切变形为主(Hashash等,2001),平面应变模型适合描述此变形模式,且地下结构动力分析模型涉及一定范围的土体,三维数值模型自由度大,计算耗时长,对计算设备要求高。

    将地下结构地震反应分析的有限元模型简化为二维平面应变模型时,核心问题是中柱的模拟,即将中柱三维空间梁简化为二维平面梁时,相关力学参数合理等效问题。韩文星(2005)和田雪娟(2010)对该问题进行了研究,虽对简化模型相关参数取值问题进行了说明,但未分析简化后模型的计算误差水平,其他学者采用二维平面应变模型进行分析时也存在同样的问题。

    本文以大开地铁车站为例,分别建立二维和三维有限元模型,对地铁车站结构进行地震反应计算,并以三维模型的计算结果作为对比基准,分析不同二维简化模型计算误差,讨论不同简化方法的合理性。

    关于中柱等效二维平面应变的等效问题,目前常用以下方法:

    (1)方法1:调整弹性模量(庄海洋等,2009

    通过调整弹性模量,保证中柱等效前、后的抗弯刚度相同,弹性模量调整系数取为中柱纵向间距的倒数,弹性模量调整公式为:

    $$ {E^*} = {\alpha _1}E $$ (1)

    式中, ${E^*} $ 为等效后中柱弹性模量;α1为弹性模量调整系数,取中柱纵向间距的倒数;E为等效前中柱弹性模量。

    (2)方法2:调整弹性模量和密度(杜修力等,2018

    通过调整弹性模量和密度,保证中柱等效前、后的截面抗弯刚度和质量均相同,弹性模量和密度调整系数均取为中柱纵向间距的倒数,弹性模量调整公式如式(1)所示,密度调整公式为:

    $$ {\rho ^*} = {\alpha _2}\rho $$ (2)

    式中, ${\rho ^*} $ 为等效后中柱密度;α2为密度调整系数,取中柱纵向间距的倒数,数值上等于α1ρ为等效前中柱密度。

    为合理建立等效简化模型,应考虑以下因素:

    (1)对中柱横向刚度、密度和阻尼特性等动力特性进行等效。地下结构二维平面应变模型相当于平面模型沿纵向连续,且纵向尺寸为单位长度,因此将中柱弹性模量除以中柱纵向间距等效其横向刚度,这是现行方法1采用的等效方法。将中柱弹性模量除以中柱纵向间距等效其横向刚度,同时将中柱密度除以中柱纵向间距等效中柱密度,这是现行方法2采用的等效方法,本文采用该等效方法。由于考虑材料阻尼特性时,目前一般采用无量纲参数阻尼比,可不进行调整。

    (2)从有限单元建模角度来看,二维梁单元是在三维梁单元的基础上,每个结点增加了3个自由度的约束,2个结点的梁单元共增加了沿纵向的6个自由度的约束,从而增大了整体刚度。为降低刚度增大带来的影响,对材料弹性模量再乘以1个调整系数,本文称该系数为空间约束影响系数,通过多次计算,建议该系数取为0.8。

    (3)求解得到二维平面应变模型的地震反应后,需还原至三维模型。采用二维平面应变模型求解得到中柱端部截面内力(弯矩、剪力和轴力)后,应乘以相应系数得到三维模型中柱端部截面的真实内力。因二维平面应变模型相当于沿地下结构纵向取单位长度的范围,而实际上地震发生时,无论是水平还是垂直分量,每根中柱承担的地震作用范围应为中柱的纵向间距,所以真实的中柱截面内力应为二维模型所得结果再乘以数值上等于中柱纵向间距的系数,该系数称为三维还原调整系数,真实的中柱端部截面内力计算公式为:

    $$ {S ^*} = {\alpha _3}S $$ (3)

    式中, $ {S ^*} $ 为真实的中柱截面内力效应,包括柱端弯矩、剪力和轴力; $ {\alpha _3} $ 为三维还原系数,取中柱纵向间距; $ S $ 为二维等效模型所得的中柱截面内力效应。

    本文建议的中柱二维等效方法称为方法3。

    为分析不同影响因素在二维等效模型中的影响程度,将三维还原调整因素引入到方法1和方法2中,即将方法1和方法2的计算结果乘以三维还原系数 $ {\alpha _3} $ ,分别称为方法1a和方法2a,将方法3计算结果与方法1和方法2计算结果进行对比,总体上可判断本文建议方法的合理性。同时,根据方法2a与方法1a的对比结果可知质量密度调整系数的影响程度,根据方法3与方法2a的对比结果可知空间约束因素的影响程度。

    本文以大开地铁车站为例,说明3种不同中柱等效方法的合理性。大开地铁车站为单层双跨框架式结构,埋深4.8 m,横断面尺寸如图1所示,横断面宽17 m,高7.17 m,侧壁厚0.85 m,顶板厚0.80 m,底板厚0.85 m。车站周围土层物理性质如表1所示(曹炳政等,2002)。计算时忽略土体的非线性特性及车站与土体的接触非线性。

    图 1  大开地铁车站结构横断面(单位:毫米)
    Figure 1.  Cross section of Subway station(Unit: mm)
    表 1  场地土物理力学参数
    Table 1.  Physical parameters of site soil properties
    土质 深度/m 密度/t·m−3 剪切波速/m·s−1 最大剪切模量/MPa 泊松比
    人工填土 0~1.0 1.9 140 38.00 0.33
    全新世砂土 1.0~5.1 1.9 140 38.00 0.32
    全新世砂土 5.1~8.3 1.9 170 56.03 0.32
    更新世黏土 8.3~11.4 1.9 190 69.99 0.40
    更新世黏土 11.4~17.2 1.9 240 111.67 0.30
    更新世砂土 17.2~22.2 2.0 330 222.24 0.26
    下载: 导出CSV 
    | 显示表格

    为尽可能减小土层人工边界反射效应的影响,根据楼梦麟等(2016)的建议,建立二维有限元模型时,在地下结构两侧向外取5倍土层深度的有限土层,竖向人工边界设置为能够保证土层沿横向水平剪切位移的侧移边界(透射边界)。竖向范围取至第6层土底部至基岩面,该底部结点约束全部自由度,作为横向水平地震的输入面。

    为在振动分析中防止有限单元竖向尺寸的滤波效应对地震反应计算结果准确性的影响,土层有限元网格的竖向尺寸按最高有效截止频率为25 Hz的正弦波对应波长的1/16~1/8控制,单元横向尺寸不超过竖向的4倍。

    地下结构侧壁和中柱采用二维梁单元离散,地下结构周围土层采用四边形等参元离散,离散后二维模型共有6 356个四边形等参元、69个二维梁单元、6 592个结点,施加边界条件后共有12 860个自由度。二维模型的有限元网格示意如图2所示。

    图 2  有限元网格示意与监测点位置(单位:毫米)
    Figure 2.  The mesh of finite element and observation points(Unit: mm)

    为保证三维和二维计算模型一致,三维模型在竖平面中的计算范围和网格尺寸与二维模型保持一致。另外,在纵向取7跨中柱的地下结构和相应的有限土层,纵向有限元网格尺寸与横向一致。四边的竖向人工边界设置为能够保证土层沿横向水平剪切位移的侧移边界。

    三维模型中,地下结构侧壁和顶、底板采用三维板单元离散,中柱采用三维梁单元离散,土体采用八结点六面体等参元离散,离散后三维模型共有10 136个板单元、165个三维梁单元、444 920个八结点六面体等参元、467 192个结点,施加边界条件后共有1 381 252个自由度。有限元模型两侧结点约束垂向和纵向自由度,底部结点约束全部自由度,同时作为地震激励的输入面。

    本文选择了3条有代表性的基岩波作为激振地震波,分别为某大型工程场地安全性评价得到的基岩人工波、阪神地震中神户大学记录的Kobe波和汶川地震中四川郫县走石台记录的基岩波,分别称之为JY波、Kobe波和WC波,峰值加速度分别为1.00、0.27、1.20 m/s2,加速度时程、傅里叶幅值谱和反应谱如图3所示。

    图 3  地震波加速度时程、傅里叶幅值谱和反应谱
    Figure 3.  Time history of exciting and its Fourier spectrum and its response spectrum

    二、三维模型各部分采用合适的有限单元离散后的运动方程为:

    $$ {\boldsymbol{m}}\ddot {\boldsymbol{u}}\left( t \right) + {\boldsymbol{c}}\dot {\boldsymbol{u}}\left( t \right) + {\boldsymbol{ku}}\left( t \right) = - {\boldsymbol{mI}}{\ddot u_{g}}\left( t \right) $$ (4)

    式中, $ {\boldsymbol{u}}\left( t \right) $ $ \dot{\boldsymbol{ u}}\left( t \right) $ $ \ddot {\boldsymbol{u}}\left( t \right) $ 分别为体系各自由度的相对位移、相对速度和相对加速度反应列向量; $ {\ddot u_{g}}\left( t \right) $ 为基岩运动加速度向量; $ {\boldsymbol{m}} $ $ {\boldsymbol{c}} $ $ {\boldsymbol{k}} $ 分别为有限元离散后的质量矩阵、阻尼矩阵和刚度矩阵。

    采用逐步积分求解运动方程可得各自由度基本未知量,如位移、速度和加速度反应时程。由基本未知量可计算其他反应量,如梁单元的弯矩、剪力和轴力,实体单元的应力、应变等。

    三维模型模态分析结果如表2所示。根据各阶振型的横向(水平向)参与系数和计算所得振型可知,第1阶和第7阶是前2个横向(水平向)的剪切振型,其频率分别为2.66、2.89 Hz。

    表 2  三维模型的前7阶自振频率及横向(水平向)振型参与系数
    Table 2.  The first seven natural frequencies of three dimension model and modal participation factor of horizontal direction
    参数 阶序
    1 2 3 4 5 6 7
    自振频率/Hz 2.66 2.72 2.73 2.76 2.77 2.79 2.89
    参与系数/×104 1.00 0 0 0 0 0 0.48
    下载: 导出CSV 
    | 显示表格

    二维模型模态分析结果如表3所示,根据各阶振型的横向(水平向)参与系数和计算所得振型可知,第1阶和第3阶是前2个横向(水平向)的剪切振型,其频率分别为2.64、2.87 Hz。

    表 3  二维模型的前7阶自振频率及横向(水平向)振型参与系数
    Table 3.  The first seven natural frequencies of two dimension model and modal participation factor of horizontal direction
    参数 阶序
    1 2 3 4 5 6 7
    自振频率/Hz 2.64 2.79 2.87 3.24 3.45 3.95 4.20
    参与系数/×104 0.23 0 0.10 0 0.25 0 0.02
    下载: 导出CSV 
    | 显示表格

    时程分析时,选用Rayleigh比例阻尼作为阻尼模型,因本文算例的基频较高,Rayleigh阻尼系数的计算采用系统前2个水平剪切振动模态对应的圆频率(舒恩,2020),即计算阻尼系数时,三维模型采用2.66、2.89 Hz,二维模型采用2.64、2.87 Hz。

    三维和二维简化模型计算得到的地震反应分别如表4表5所示。表4中中柱地震反应未列出轴力,是因为单层双跨结构在仅有水平地震作用下可认为是反对称结构,中柱近似无轴力(许紫刚等,2019)。

    表 4  中柱地震反应峰值
    Table 4.  Peak seismic response of the inner column
    激励 考察点及反应量 三维模型 方法1 方法1a 方法2 方法2a 方法3
    JY波 柱顶
    (监测点P2
    弯矩Mz/kN·m 253.61 84.12
    (误差−66.83%)
    294.42
    (误差16.09%)
    92.21
    (误差−63.64%)
    322.72
    (误差27.25%)
    263.87
    (误差4.04%)
    剪力Fx/kN 60.40 13.71
    (误差−77.29%)
    48.00
    (误差−20.52%)
    22.54
    (误差−62.68%)
    78.89
    (误差30.63%)
    64.60
    (误差6.97%)
    柱底
    (监测点P3
    弯矩Mz/kN·m 246.16 95.92
    (误差−61.03%)
    335.74
    (误差36.39%)
    88.93
    (误差−63.87%)
    311.26
    (误差26.44%)
    254.02
    (误差3.19%)
    剪力Fx/kN 77.62 33.22
    (误差−57.20%)
    116.28
    (误差49.81%)
    27.61
    (误差−64.43%)
    96.62
    (误差24.48%)
    78.74
    (误差1.44%)
    Kobe波 柱顶
    (监测点P2
    弯矩Mz/kN·m 52.13 16.82
    (误差−67.72%)
    58.89
    (误差12.97%)
    18.51
    (误差−64.49%)
    64.78
    (误差24.27%)
    52.90
    (误差1.48%)
    剪力Fx/kN 12.57 2.90
    (误差−76.94%)
    10.15
    (误差−19.29%)
    4.55
    (误差−63.80%)
    15.93
    (误差26.71%)
    13.03
    (误差3.62%)
    柱底
    (监测点P3
    弯矩Mz/kN·m 51.07 19.46
    (误差−61.89%)
    68.12
    (误差33.38%)
    18.02
    (误差−64.72%)
    63.07
    (误差23.50%)
    51.40
    (误差0.65%)
    剪力Fx/kN 16.08 6.60
    (误差−58.95%)
    23.10
    (误差43.69%)
    5.61
    (误差−65.09%)
    19.64
    (误差22.17%)
    16.00
    (误差−0.46%)
    WC波 柱顶
    (监测点P2
    弯矩Mz/kN·m 307.13 96.51
    (误差−68.58%)
    337.79
    (误差9.98%)
    105.99
    (误差−65.49%)
    370.98
    (误差20.79%)
    303.11
    (误差−1.31%)
    剪力Fx/kN 73.78 16.79
    (误差−77.24%)
    58.78
    (误差−20.33%)
    26.11
    (误差−64.62%)
    91.37
    (误差23.84%)
    74.76
    (误差1.34%)
    柱底
    (监测点P3
    弯矩Mz/kN·m 304.37 112.25
    (误差−63.12%)
    392.87
    (误差29.08%)
    103.49
    (误差−66.00%)
    362.23
    (误差19.01%)
    295.34
    (误差−2.96%)
    剪力Fx/kN 96.71 38.91
    (误差−59.77%)
    136.17
    (误差40.80%)
    32.45
    (误差−66.45%)
    113.56
    (误差17.42%)
    92.43
    (误差−4.43%)
    下载: 导出CSV 
    | 显示表格
    表 5  关键点地震反应峰值
    Table 5.  Peak seismic response of observation points
    激励 考察点
    及反应量
    三维模型 方法1 方法2 方法3
    JY波 地表
    (监测点P1
    加速度a/m·s−2 4.22 4.49(误差6.32%) 4.48(误差6.21%) 4.49(误差6.32%)
    位移u/mm 14.53 15.58(误差7.21%) 15.47(误差6.43%) 15.50(误差6.69%)
    柱顶
    (监测点P2
    加速度a/m·s−2 4.53 4.99(误差10.22%) 4.96(误差9.49%) 4.98(误差9.81%)
    位移u/mm 13.69 15.38(误差12.29%) 15.23(误差11.23%) 15.30(误差11.7%)
    侧壁
    (监测点P4
    加速度a/m·s−2 1.64 1.62(误差−1.10%) 1.61(误差−2.13%) 1.61(误差−2.25%)
    位移u/mm 3.41 3.48(误差1.85%) 3.46(误差1.28%) 3.45(误差1.02%)
    Kobe波 地表
    (监测点P1
    加速度a/m·s−2 0.86 0.87(误差1.20%) 0.87(误差1.01%) 0.87(误差1.00%)
    位移u/mm 3.10 3.18(误差2.41%) 3.16(误差1.96%) 3.17(误差2.07%)
    柱顶
    (监测点P2
    加速度a/m·s−2 0.80 0.85(误差6.06%) 0.84(误差5.51%) 0.84(误差5.82%)
    位移u/mm 2.94 3.10(误差5.45%) 3.08(误差4.80%) 3.09(误差5.14%)
    侧壁
    (监测点P4
    加速度a/m·s−2 0.30 0.29(误差−2.89%) 0.29(误差−3.72%) 0.29(误差−3.96%)
    位移u/mm 0.78 0.77(误差−1.11%) 0.77(误差−1.50%) 0.77(误差−1.76%)
    WC波 地表
    (监测点P1
    加速度a/m·s−2 5.13 5.18(误差0.82%) 5.17(误差0.69%) 5.17(误差0.63%)
    位移u/mm 18.75 18.62(误差−0.72%) 18.55(误差−1.10%) 18.56(误差−1.03%)
    柱顶
    (监测点P2
    加速度a/m·s−2 4.80 4.87(误差1.54%) 4.86(误差1.27%) 4.87(误差1.44%)
    位移u/mm 17.81 18.33(误差2.94%) 18.23(误差2.36%) 18.28(误差2.66%)
    侧壁
    (监测点P4
    加速度a/m·s−2 1.39 1.40(误差1.12%) 1.39(误差0.62%) 1.39(误差0.47%)
    位移u/mm 4.45 4.33(误差−2.66%) 4.32(误差−2.96%) 4.30(误差−3.26%)
    下载: 导出CSV 
    | 显示表格

    表4可知,无论是何种地震波激励,方法3均优于方法1和方法2。需注意,方法1和方法2计算结果均为负误差,说明计算结果偏小,从工程安全性角度评价,方法1和方法2偏于不安全。方法1和方法2计算误差水平相当,方法1所得柱顶计算结果误差略大于方法2,但所得柱底计算结果误差略小于方法2,说明材料密度的等效对中柱地震反应的影响较小。方法1a的计算误差明显小于方法1,且方法2a的计算误差明显小于方法2,说明二维模型计算结果不能直接等效为三维真实情况,借用方法3的三维还原系数可改善相应的计算精度。方法3计算误差明显小于方法2a,说明空间约束因素对中柱计算结果有较大影响。

    表4表5可知,本文采用的3种二维等效模型对地下结构整体的运动反应影响较小,因采用平面应变问题求解地下结构的横向水平地震反应在力学基本原理上符合工程实际情况,但由于中柱平面应变等效是人为的主观假定,因此不同方法下计算得到的中柱截面内力相差较大,相比之下,方法3的近似计算精度远高于方法1和方法2。从本文算例结果来看,无论采用哪种方法,计算得到的加速度和位移误差水平相近,位移和加速度反应误差均在13%以内,当采用方法3时,柱弯矩和剪力误差≤7%,较方法1和方法2得到了明显改进。

    (1)对地下结构进行地震反应分析时,可将三维模型简化为平面应变二维模型进行计算。简化时应重点关注结构中柱的等效方法,如果简化方法选择不当,可能会给地下结构内力计算带来较大误差,应重视中柱的等效问题。

    (2)现行简化方法1和方法2会低估中柱地震反应,最大低估程度可达77.24%,从工程安全性角度评价,方法1和方法2偏于不安全。本文提出的方法3可明显减小中柱地震反应误差,使中柱弯矩和剪力误差≤7%。

    (3)方法1、方法2和方法3均可得到误差较小的加速度和位移反应,3种方法误差水平相近,方法3可使加速度和位移计算误差≤13%。

  • 图  1  大开地铁车站结构横断面(单位:毫米)

    Figure  1.  Cross section of Subway station(Unit: mm)

    图  2  有限元网格示意与监测点位置(单位:毫米)

    Figure  2.  The mesh of finite element and observation points(Unit: mm)

    图  3  地震波加速度时程、傅里叶幅值谱和反应谱

    Figure  3.  Time history of exciting and its Fourier spectrum and its response spectrum

    表  1  场地土物理力学参数

    Table  1.   Physical parameters of site soil properties

    土质 深度/m 密度/t·m−3 剪切波速/m·s−1 最大剪切模量/MPa 泊松比
    人工填土 0~1.0 1.9 140 38.00 0.33
    全新世砂土 1.0~5.1 1.9 140 38.00 0.32
    全新世砂土 5.1~8.3 1.9 170 56.03 0.32
    更新世黏土 8.3~11.4 1.9 190 69.99 0.40
    更新世黏土 11.4~17.2 1.9 240 111.67 0.30
    更新世砂土 17.2~22.2 2.0 330 222.24 0.26
    下载: 导出CSV

    表  2  三维模型的前7阶自振频率及横向(水平向)振型参与系数

    Table  2.   The first seven natural frequencies of three dimension model and modal participation factor of horizontal direction

    参数 阶序
    1 2 3 4 5 6 7
    自振频率/Hz 2.66 2.72 2.73 2.76 2.77 2.79 2.89
    参与系数/×104 1.00 0 0 0 0 0 0.48
    下载: 导出CSV

    表  3  二维模型的前7阶自振频率及横向(水平向)振型参与系数

    Table  3.   The first seven natural frequencies of two dimension model and modal participation factor of horizontal direction

    参数 阶序
    1 2 3 4 5 6 7
    自振频率/Hz 2.64 2.79 2.87 3.24 3.45 3.95 4.20
    参与系数/×104 0.23 0 0.10 0 0.25 0 0.02
    下载: 导出CSV

    表  4  中柱地震反应峰值

    Table  4.   Peak seismic response of the inner column

    激励 考察点及反应量 三维模型 方法1 方法1a 方法2 方法2a 方法3
    JY波 柱顶
    (监测点P2
    弯矩Mz/kN·m 253.61 84.12
    (误差−66.83%)
    294.42
    (误差16.09%)
    92.21
    (误差−63.64%)
    322.72
    (误差27.25%)
    263.87
    (误差4.04%)
    剪力Fx/kN 60.40 13.71
    (误差−77.29%)
    48.00
    (误差−20.52%)
    22.54
    (误差−62.68%)
    78.89
    (误差30.63%)
    64.60
    (误差6.97%)
    柱底
    (监测点P3
    弯矩Mz/kN·m 246.16 95.92
    (误差−61.03%)
    335.74
    (误差36.39%)
    88.93
    (误差−63.87%)
    311.26
    (误差26.44%)
    254.02
    (误差3.19%)
    剪力Fx/kN 77.62 33.22
    (误差−57.20%)
    116.28
    (误差49.81%)
    27.61
    (误差−64.43%)
    96.62
    (误差24.48%)
    78.74
    (误差1.44%)
    Kobe波 柱顶
    (监测点P2
    弯矩Mz/kN·m 52.13 16.82
    (误差−67.72%)
    58.89
    (误差12.97%)
    18.51
    (误差−64.49%)
    64.78
    (误差24.27%)
    52.90
    (误差1.48%)
    剪力Fx/kN 12.57 2.90
    (误差−76.94%)
    10.15
    (误差−19.29%)
    4.55
    (误差−63.80%)
    15.93
    (误差26.71%)
    13.03
    (误差3.62%)
    柱底
    (监测点P3
    弯矩Mz/kN·m 51.07 19.46
    (误差−61.89%)
    68.12
    (误差33.38%)
    18.02
    (误差−64.72%)
    63.07
    (误差23.50%)
    51.40
    (误差0.65%)
    剪力Fx/kN 16.08 6.60
    (误差−58.95%)
    23.10
    (误差43.69%)
    5.61
    (误差−65.09%)
    19.64
    (误差22.17%)
    16.00
    (误差−0.46%)
    WC波 柱顶
    (监测点P2
    弯矩Mz/kN·m 307.13 96.51
    (误差−68.58%)
    337.79
    (误差9.98%)
    105.99
    (误差−65.49%)
    370.98
    (误差20.79%)
    303.11
    (误差−1.31%)
    剪力Fx/kN 73.78 16.79
    (误差−77.24%)
    58.78
    (误差−20.33%)
    26.11
    (误差−64.62%)
    91.37
    (误差23.84%)
    74.76
    (误差1.34%)
    柱底
    (监测点P3
    弯矩Mz/kN·m 304.37 112.25
    (误差−63.12%)
    392.87
    (误差29.08%)
    103.49
    (误差−66.00%)
    362.23
    (误差19.01%)
    295.34
    (误差−2.96%)
    剪力Fx/kN 96.71 38.91
    (误差−59.77%)
    136.17
    (误差40.80%)
    32.45
    (误差−66.45%)
    113.56
    (误差17.42%)
    92.43
    (误差−4.43%)
    下载: 导出CSV

    表  5  关键点地震反应峰值

    Table  5.   Peak seismic response of observation points

    激励 考察点
    及反应量
    三维模型 方法1 方法2 方法3
    JY波 地表
    (监测点P1
    加速度a/m·s−2 4.22 4.49(误差6.32%) 4.48(误差6.21%) 4.49(误差6.32%)
    位移u/mm 14.53 15.58(误差7.21%) 15.47(误差6.43%) 15.50(误差6.69%)
    柱顶
    (监测点P2
    加速度a/m·s−2 4.53 4.99(误差10.22%) 4.96(误差9.49%) 4.98(误差9.81%)
    位移u/mm 13.69 15.38(误差12.29%) 15.23(误差11.23%) 15.30(误差11.7%)
    侧壁
    (监测点P4
    加速度a/m·s−2 1.64 1.62(误差−1.10%) 1.61(误差−2.13%) 1.61(误差−2.25%)
    位移u/mm 3.41 3.48(误差1.85%) 3.46(误差1.28%) 3.45(误差1.02%)
    Kobe波 地表
    (监测点P1
    加速度a/m·s−2 0.86 0.87(误差1.20%) 0.87(误差1.01%) 0.87(误差1.00%)
    位移u/mm 3.10 3.18(误差2.41%) 3.16(误差1.96%) 3.17(误差2.07%)
    柱顶
    (监测点P2
    加速度a/m·s−2 0.80 0.85(误差6.06%) 0.84(误差5.51%) 0.84(误差5.82%)
    位移u/mm 2.94 3.10(误差5.45%) 3.08(误差4.80%) 3.09(误差5.14%)
    侧壁
    (监测点P4
    加速度a/m·s−2 0.30 0.29(误差−2.89%) 0.29(误差−3.72%) 0.29(误差−3.96%)
    位移u/mm 0.78 0.77(误差−1.11%) 0.77(误差−1.50%) 0.77(误差−1.76%)
    WC波 地表
    (监测点P1
    加速度a/m·s−2 5.13 5.18(误差0.82%) 5.17(误差0.69%) 5.17(误差0.63%)
    位移u/mm 18.75 18.62(误差−0.72%) 18.55(误差−1.10%) 18.56(误差−1.03%)
    柱顶
    (监测点P2
    加速度a/m·s−2 4.80 4.87(误差1.54%) 4.86(误差1.27%) 4.87(误差1.44%)
    位移u/mm 17.81 18.33(误差2.94%) 18.23(误差2.36%) 18.28(误差2.66%)
    侧壁
    (监测点P4
    加速度a/m·s−2 1.39 1.40(误差1.12%) 1.39(误差0.62%) 1.39(误差0.47%)
    位移u/mm 4.45 4.33(误差−2.66%) 4.32(误差−2.96%) 4.30(误差−3.26%)
    下载: 导出CSV
  • 曹炳政, 罗奇峰, 马硕等, 2002. 神户大开地铁车站的地震反应分析. 地震工程与工程振动, 22(4): 102—107 doi: 10.3969/j.issn.1000-1301.2002.04.017

    Cao B. Z. , Luo Q. F. , Ma S. , et al. , 2002. Seismic response analysis of Dakai subway station in Hyogoken-Nanbu earthquake. Earthquake Engineering and Engineering Vibration, 22(4): 102—107. (in Chinese) doi: 10.3969/j.issn.1000-1301.2002.04.017
    陈国兴, 孙瑞瑞, 赵丁凤等, 2019. 海底盾构隧道纵向地震反应特征的子模型分析. 岩土工程学报, 41(11): 1983—1991

    Chen G. X. , Sun R. R. , Zhao D. F. , et al. , 2019. Longitudinal Seismic response characteristics of seabed shield tunnels using submodeling analysis. Chinese Journal of Geotechnical Engineering, 41(11): 1983—1991. (in Chinese)
    杜修力, 马超, 路德春等, 2017. 大开地铁车站地震破坏模拟与机理分析. 土木工程学报, 50(1): 53—62, 69 doi: 10.15951/j.tmgcxb.2017.01.007

    Du X. L. , Ma C. , Lu D. C. , et al. , 2017. Collapse simulation and failure mechanism analysis of the Daikai subway station under seismic loads. China Civil Engineering Journal, 50(1): 53—62, 69. (in Chinese) doi: 10.15951/j.tmgcxb.2017.01.007
    杜修力, 康凯丽, 许紫刚等, 2018. 地下结构地震反应的主要特征及规律. 土木工程学报, 51(7): 11—21 doi: 10.15951/j.tmgcxb.2018.07.002

    Du X. L. , Kang K. L. , Xu Z. G. , et al. , 2018. Main characteristics and rules of seismic response for underground structures. China Civil Engineering Journal, 51(7): 11—21. (in Chinese) doi: 10.15951/j.tmgcxb.2018.07.002
    韩文星, 2005. 软土地铁车站结构横向抗震设计方法研究. 上海: 同济大学.
    刘晶波, 李彬, 2006. 地铁地下结构抗震分析及设计中的几个关键问题. 土木工程学报, 39(6): 106—110 doi: 10.3321/j.issn:1000-131X.2006.06.019

    Liu J. B. , Li B. , 2006. Issues on the seismic analysis and design of subway structures. China Civil Engineering Journal, 39(6): 106—110. (in Chinese) doi: 10.3321/j.issn:1000-131X.2006.06.019
    楼梦麟, 董云, 张如林, 2016. 沉管隧道地震反应分析局部精细化建模中的几个问题. 岩土工程学报, 38(9): 1705—1712 doi: 10.11779/CJGE201609018

    Lou M. L. , Dong Y. , Zhang R. L. , 2016. Several problems in refined local modeling for seismic response analysis of immersed tunnel. Chinese Journal of Geotechnical Engineering, 38(9): 1705—1712. (in Chinese) doi: 10.11779/CJGE201609018
    舒恩, 2020. 增强地下结构中柱抗震能力的方法研究. 淮南: 安徽理工大学.

    Shu E., 2020. Study on the method of strengthening the seismic capacity of the middle column of underground structure. Huainan: Anhui University of Technology. (in Chinese)
    田雪娟, 2010. 地铁车站抗震分析. 北京: 北京交通大学.

    Tian X. J., 2010. Seismic analysis of subway station. Beijing: Beijing Jiaotong University. (in Chinese)
    王国波, 王亚西, 陈斌等, 2015. 隧道–土体–地表结构相互作用体系地震响应影响因素分析. 岩土力学与工程学报, 34(6): 1276—1287

    Wang G. B. , Wang Y. X. , Chen B. , et al. , 2015. Analysis of factors influencing seismic responses of tunnel-soil-ground structural system. Chinese Journal of Rock Mechanics and Engineering, 34(6): 1276—1287. (in Chinese)
    许紫刚, 杜修力, 许成顺等, 2019. 地下结构地震反应分析中场地瑞利阻尼构建方法比较研究. 岩土力学, 40(12): 4838—4847 doi: 10.16285/j.rsm.2018.1913

    Xu Z. G. , Du X. L. , Xu C. S. , et al. , 2019. Comparison of determination methods of site Rayleigh damping coefficients in seismic responses analysis of underground structures. Rock and Soil Mechanics, 40(12): 4838—4847. (in Chinese) doi: 10.16285/j.rsm.2018.1913
    庄海洋, 王修信, 陈国兴, 2009. 软土层埋深变化对地铁车站结构地震反应的影响规律研究. 岩土工程学报, 31(8): 1258—1266 doi: 10.3321/j.issn:1000-4548.2009.08.017

    Zhuang H. Y. , Wang X. X. , Chen G. X. , 2009. Earthquake responses of subway station with different depths of soft soil. Chinese Journal of Geotechnical Engineering, 31(8): 1258—1266. (in Chinese) doi: 10.3321/j.issn:1000-4548.2009.08.017
    庄海洋, 任佳伟, 王瑞等, 2019. 两层三跨框架式地铁地下车站结构弹塑性工作状态与抗震性能水平研究. 岩土工程学报, 41(1): 131—138

    Zhuang H. Y. , Ren J. W. , Wang R. , et al. , 2019. Elasto-plastic working states and seismic performance levels of frame-type subway underground station with two layers and three spans. Chinese Journal of Geotechnical Engineering, 41(1): 131—138. (in Chinese)
    Hashash Y. M. A. , Hook J. J. , Schmidt B. , et al. , 2001. Seismic design and analysis of underground structures. Tunnelling and Underground Space Technology, 16(4): 247—293. doi: 10.1016/S0886-7798(01)00051-7
  • 期刊类型引用(2)

    1. 许紫刚,李淳宇,徐长节,庄海洋,张季. 地下连续墙对软弱夹层场地地下结构地震反应影响研究. 世界地震工程. 2024(02): 78-87 . 百度学术
    2. 许紫刚,郑养宁,庄海洋,张季,韩润波. 基于不同分析方法的设置分体柱地铁车站结构抗震性能研究. 防灾减灾工程学报. 2024(06): 1359-1367 . 百度学术

    其他类型引用(2)

  • 加载中
图(3) / 表(5)
计量
  • 文章访问数:  143
  • HTML全文浏览量:  153
  • PDF下载量:  37
  • 被引次数: 4
出版历程
  • 收稿日期:  2021-11-11
  • 刊出日期:  2022-12-31

目录

/

返回文章
返回