引用本文
王笃国, 贾其军, 尤红兵, 翟奋楼, 薛齐齐, 喻烟. 2020. 内马铁路一期工程局部不规则二维场地地震反应分析. 震灾防御技术, 15(4): 708-717, DOI:10.11899/zzfy20200405.
权限
内马铁路一期工程局部不规则二维场地地震反应分析
王笃国2) 贾其军1) 尤红兵2) 翟奋楼3) 薛齐齐3) 喻烟4)
1) 中国路桥工程有限责任公司, 北京 100011
2) 中国地震灾害防御中心, 北京 100029
3) 北京市计算中心, 北京 100101
4) 中国地震局地壳应力研究所, 北京 100085
[基金项目]:地震科技星火计划项目(XH18060);北京市科委重大科技研发攻关项目(Z181100003918005);北京市科学技术研究院创新团队计划(IG201506C2)
[收稿日期]:2020-05-20
[作者简介]:王笃国, 男, 生于1979年。博士, 正高级工程师。主要从事岩土工程与地震工程科研工作。E-mail: wangduguo@163.com
摘要

基于内马铁路一期工程地质资料,利用ABAQUS有限元软件建立某局部不规则二维场地有限元计算模型,利用Python语言进行二次开发,编制二维等效线性化计算程序。开展50年超越概率63%、10%和2%水准下的土层地震反应分析计算,对不规则地形不同位置处的地震动峰值加速度和频谱特性进行对比分析。研究结果表明,不规则地形对地震动特性的影响显著。本文研究结果对内马铁路一期工程地震动参数确定具有指导意义,同时可为跨越不规则地形工程结构抗震设计提供参考。

关键词: 铁路工程  二维等效线性化  不规则地形  地震动参数  土层反应  


引言

肯尼亚内罗毕至马拉巴新建标轨铁路一期工程(以下简称“内马铁路一期工程”)位于非洲东部,全线位于东非高原区,线路横跨埃塞俄比亚-肯尼亚裂谷肯尼亚段,是世界最大的断层陷落带,沿线地势高低起伏较大,第四系覆盖层较厚。本工程所在区域地震地质条件复杂,研究基础薄弱,铁路跨越多个不规则局部地形,需开展专题研究确定场地地震动参数。

局部不规则地形、土层分布不均匀性是影响地表地震动的重要因素,对该问题的研究一直是地震工程领域备受关注的课题。工程上通常采用一维等效线性化方法确定水平成层场地地震动参数,但对于地形起伏较大的场地,仅进行一维地震反应分析是远远不够的,需采用二维或三维模型(Kamiyama等,2002车伟等,2008)进行研究。战吉艳等(2013)采用ABAQUS有限元软件,利用修正Martin-Seed-Davidenkov黏弹性动力本构模型,研究了横向不均匀深厚软土场地非线性地震效应;赖强林等(2014)研究了复杂V型谷地一维模型与二维模型地震动参数的差异;王永幸等(2015)杨笑梅等(2017)分别给出了二维土层地震反应分析频域、时域等效线性化解法。上述研究均认为,对于复杂不规则或横向不均匀场地,采用二维非线性模型能得到更合理的地震动参数。基于此,针对内马铁路一期工程某典型局部不规则场地,利用ABAQUS有限元软件建立二维场地计算模型,利用Python语言进行二次开发,采用等效线性化方法进行土层地震反应分析计算,得到场地地震响应,并探讨二维不规则地形地震动分布规律,为内马铁路一期工程抗震设防提供科学依据。

1 工程概况与地震动输入
1.1 工程概况

拟分析的局部不规则地形工程场地位于东非大裂谷西翼山区,第四系覆盖层较厚,达80—100m,基岩以第四系玄武岩、凝灰岩为主。根据已有工程地质资料,确定工程地质剖面图如图 1所示,水平向长90m,竖向深100m。


图 1 工程地质剖面图 Fig. 1 Engineering geological profile

将工程场地各钻孔土层资料进行整理,可知该场地土体类别主要为黏土和细砂,钻孔ZK1与ZK2土体类型和各项参数(陆鸣等,2016)分别见表 1表 2,其中土层沿浅部至深部依次进行编号,动三轴序号为10的土体为计算基底,其剪切模量比通常取1.0,阻尼比取0.05。

表 1 ZK1钻孔土层分布资料 Table 1 Soil layer distribution data of ZK1 borehole
表 2 ZK2钻孔土层分布资料 Table 2 Soil layer distribution data of ZK2 borehole

场地9种不同层位土体剪切模量和阻尼比随剪应变变化曲线(陆鸣等,2016)如图 2所示,图中土样编号1—9分别对应表 1表 2中动三轴序号1—9,为不同土类不同深度土样动三轴试验结果。


图 2 9种土层剪切模量和阻尼比随剪应变变化曲线 Fig. 2 G/Gmax~γ and D~γ curves for 9 representative soil specimens
1.2 地震动输入

假定地震动为垂直入射,取值参考《肯尼亚内罗毕-马拉巴铁路CK0-CK120段地震危险性分析报告》(陆鸣等,2016),工程场地50年超越概率63%、10%和2%的水平向基岩加速度反应谱如图 3所示,每个概率水准下拟合3条人工地震动时程样本,拟合谱和目标谱周期控制点上的误差不超过5%,50年超越概率63%水准下的3条人工地震动时程样本曲线如图 4所示。


图 3 50年不同超越概率下水平向基岩加速度反应谱(阻尼比为5%) Fig. 3 Acceleration response spectrum for rock site with different exceedance probabilities in 50 years (damping ratio 5%)

图 4 50年超越概率63%水准下的3条人工地震动时程样本曲线 Fig. 4 Acceleration response history curves for rock site with 63% probability of exceedance in 50 years
2 二维等效线性化模型
2.1 二维等效线性化方法

本文采用等效线性化模型描述场地非线性特性,二维模型等效剪应变取值采用王笃国等(2016)给出的方法。二维计算模型示意如图 5所示,计算区域包括水平成层介质区和不规则地形区,计算区域两侧设置黏弹性人工边界,基岩基底为水平向。采用时域等效线性化方法进行计算,阻尼矩阵的选取采用Park等(2004)给出的完整Rayleigh阻尼公式,即[C]=α[M]+β[K],目标频率f的选取参考马俊玲等(2013)的研究结果,采用第1阶和第3阶自振频率,阻尼系数αβ取值公式为:

$ \alpha =2\zeta \left(\frac{{\omega }_{1}{\omega }_{3}}{{\omega }_{1}+{\omega }_{3}}\right),\beta =2\zeta \left(\frac{1}{{\omega }_{1}+{\omega }_{3}}\right)$ (1)

图 5 二维模型计算示意图 Fig. 5 Schematic diagram for two dimensional simulation model

式中,ζ为土体阻尼比;ω1=2πf1ω3=2πf3f1f3可由下式求得:

${f_n} = \frac{{{v_{{\rm{seq}}}}}}{{{\rm{4}}H}}(2n - 1), \;{v_{{\rm{seq}}}} = \frac{H}{{\sum\limits_{i = 1}^N {({h_i}/{v_{{\rm{s}}i}})} }}$ (2)

式中,n分别取值1和3;hi为第i层土体厚度;vsi为第i层土体剪切波速;$H = \sum\limits_{i = 1}^N {{h_i}} $

计算时对每个单元剪切模量和阻尼比与上次迭代计算结果进行对比,所有单元误差均小于5%时停止迭代,进行二维线弹性场地地震反应分析计算,得到场地地震响应。

2.2 人工边界地震动输入

计算区域两侧采用黏弹性人工边界,地震动输入时需将波动输入转换为等效荷载施加在人工边界上,并保证人工边界上的位移和应力与原自由场相同。本文采用刘晶波等(1998)给出的边界节点等效荷载计算公式:

$p(t) = \sigma (x, y, z, t) + c\dot u(x, y, z, t) + ku(x, y, z, t)$ (3)

式中,$p(t)$为边界点处等效荷载;$\sigma (x, y, z, t)$$\dot u(x, y, z, t)$$u(x, y, z, t)$分别为边界点处自由场求解得到的应力、速度和位移;ck分别为黏弹性人工边界法向和切向弹簧阻尼系数,求解方法参见刘晶波等(1998)

自由场求解采用一维等效线性化方法,一维计算模型可由计算区域两侧水平成层介质区得到。求得计算区域两侧自由场地震响应后,根据式(3)可得到等效荷载,将等效荷载分别施加于两侧人工边界,实现黏弹性人工边界地震动输入。

2.3 场地有限元模型建立

利用ABAQUS有限元软件建立模型,采用Springs/Dashpots单元对黏弹性边界进行模拟,Springs/Dashpots单元为弹簧阻尼器单元,切向刚度KbT=0.5G/R,法向刚度KbN=G/R,切向阻尼系数CbT=ρvs,法向阻尼系数CbN=ρvp,其中G为自由场不同土层等效线性化计算得到的等效剪切模量,R近似取不规则地形中心点到人工边界点的距离,vsvp分别为自由场不同土层等效线性化计算得到的剪切波波速和压缩波波速。土体单元大部分采用四边形平面应变单元CPE4R进行模拟,极少部分采用三角形平面应变单元CPE3进行模拟,本次计算剖面单元总数为33096个,节点总数为36369,长度、宽度方向单元尺寸基本小于1m。利用Python语言进行二次开发,实现了二维等效线性化方法的迭代算法,本次计算采用隐式H.H.T法求解动力方程,有限元模型如图 6所示,图中L1、L2、L3、R1、R2、R3为不规则地形不同位置处计算点,模型尺寸为375m×88m(长×宽)。


图 6 二维场地有限元模型 Fig. 6 Two dimensional FEM site model
3 计算结果与分析

在50年超越概率63%、10%和2%水准下,取每个概率水准3条时程样本开展场地土层地震反应计算(共9种工况),分析峰值加速度和频谱特性,探讨局部不规则地形地震动分布规律。

3.1 峰值加速度

6个计算点50年超越概率63%、10%和2%水准下地震动峰值加速度计算结果见表 3

表 3 6个计算点50年超越概率63%、10%和2%水准下地震动峰值加速度(单位:gal) Table 3 Peak ground motions for 6 points with 63%、10% and 2% probabilities of exceedance in 50 years (Unit: gal)

表 3可知,同一工况下不同位置处地震动峰值加速度存在一定变化,底部峰值加速度较小,顶部峰值加速度较大,总体呈越靠近顶部峰值加速度越大的趋势;以底部峰值加速度为基准点,50年超越概率63%、10%和2%水准下的顶部峰值加速度放大系数分别为1.24—1.30、1.14—1.19和1.11—1.15,表明随着地震动输入强度的增加,放大系数呈略减小的趋势。

3.2 频谱特性

50年超越概率63%、10%和2%水准下不同位置处地震动加速度放大系数谱如图 7所示,限于篇幅问题,本文仅给出不规则地形左侧L1、L2和L3计算点每个水准下1条样本时程计算结果。由图 7可知:(1)3个工况下的地震动加速度放大系数谱均呈双峰或多峰特性;(2)同一概率水准下不同位置处的地震动加速度放大系数谱谱形存在一定差异,总体上看,顶部计算点长周期段地震动加速度放大系数较大。


图 7 不规则地形左侧不同位置处地震动加速度放大系数谱(阻尼比为5%)(一) Fig. 7 Acceleration amplification facor spectrum for different positions in left side of irregular topography (damping ratio 5%)
4 结论

基于内马铁路一期工程某局部不规则地形工程地质资料,利用ABAQUS有限元软件建立二维场地有限元模型,利用Python语言进行二次开发,实现二维等效线性化方法的迭代算法,共进行3个概率水准、9种工况下的土层反应计算,得出以下结论:

(1)不规则地形不同位置处地震动峰值加速度存在一定差异,通常情况下顶部峰值加速度高于底部,且随着地震动输入强度的增大,加速度放大系数呈略减小的趋势,50年超越概率63%、10%和2%水准下的顶部峰值加速度放大系数分别为1.24—1.30、1.14—1.19和1.11—1.15。

(2)不规则地形不同位置处地震动放大系数谱呈双峰或多峰特性,但谱形有所不同,主要体现在长周期段放大效应不同,顶部计算点长周期段地震动加速度放大系数较大。

本文给出的结果为内马铁路一期工程确定抗震设防要求提供了科学依据,同时可为类似跨越局部不规则地形工程结构地震动参数的确定提供参考。

参考文献
车伟, 罗奇峰, 2008. 复杂地形条件下地震波的传播研究[J]. 岩土工程学报, 30(9): 1333-1337. DOI:10.3321/j.issn:1000-4548.2008.09.012
赖强林, 杨笑梅, 2014. 复杂V型谷地条件下地震动数值模拟的比较[J]. 地震工程与工程振动, 34(S1): 41-46.
刘晶波, 吕彦东, 1998. 结构-地基动力相互作用问题分析的一种直接方法[J]. 土木工程学报, 31(3): 55-64.
陆鸣, 李小军, 周正华, 等, 2016. 肯尼亚内罗毕至马拉巴铁路CK0-CK120段地震危险性分析[M]. 北京: 中国地震局地壳应力研究所, 96-104.
马俊玲, 丁海平, 2013. 土层地震反应分析中不同阻尼取值的影响比较[J]. 防灾减灾工程学报, 33(5): 517-523, 547.
王笃国, 赵成刚, 2016. 地震波斜入射时二维成层介质自由场求解的等效线性化方法[J]. 岩土工程学报, 38(3): 554-561.
王永幸, 杨仕升, 高鹏飞, 2015. 基于等效线性化的二维地震反应分析研究[J]. 广西大学学报(自然科学版), 40(3): 692-698.
杨笑梅, 赖强林, 2017. 二维土层地震反应分析的时域等效线性化解法[J]. 岩土力学, 38(3): 847-856.
战吉艳, 陈国兴, 刘建达, 等, 2013. 远场大地震作用下大尺度深软场地的非线性地震效应分析[J]. 岩土力学, 34(11): 3229-3238.
Kamiyama M., Satoh T., 2002. Seismic response analysis of laterally inhomogeneous ground with emphasis on strains[J]. Soil Dynamics and Earthquake Engineering, 22(9-12): 877-884. DOI:10.1016/S0267-7261(02)00110-0
Park D., Hashash Y. M. A., 2004. Soil damping formulation in nonlinear time domain site response analysis[J]. Journal of Earthquake Engineering, 8(2): 249-274.