Dynamic properties of undisturbed sand in deep overburden of a dam foundation
-
摘要: 超深厚覆盖层地基是当前水利水电工程建设所面临的巨大挑战, 覆盖层土体的动力力学特性是进行地基土液化判别和坝体稳定分析的基础。本文以某拟建大型土石坝工程超深厚覆盖层地基中的砂土层为研究对象,对原状土样(通过深围井直接开挖取得的)和重塑土样开展了一系列的高围压条件下的不排水循环三轴试验。研究表明:原状砂土的动强度高于重塑砂土,且随着围压的增大二者之间的差别逐渐减小;原状土样和重塑土样的孔压、轴向应变发展规律类似;固结应力比对动强度指标及土样的破坏模式都有明显的影响,随着固结应力比的增大,动强度增大,土样破坏模式逐渐由液化破坏转变为塑性应变累积破坏。Abstract: The ultra-deep overburden foundation is a great challenge for the construction of water conservancy and hydropower projects. The dynamic mechanical properties of overburden soil are the basis for discriminating liquefaction of foundation soil and analyzing dam stability. The sand layer in the ultra-deep overburden foundation of a proposed large-scale earth rock dam project is taken as the research object. A series of undrained cyclic triaxial tests are conducted on the undisturbed soil samples (directly extracted from an ultra-deep in-situ test well) and remolded soil samples under high confining pressure conditions. The laboratory results demonstrate that the dynamic strength of the undisturbed sand is higher than that of the remolded sand, and the difference between them gradually decreases with the increase of the confining pressure. The evolution laws of pore pressure and axial strain of the undisturbed and remolded soil samples are similar. The consolidation stress ratio has obvious influences on both the dynamic strength index and the failure mode of the soil samples. With the increase of the consolidation stress ratio, the dynamic strength increases, and the failure mode of the soil samples gradually changes from liquefaction failure to plastic strain accumulation one.
-
0. 引言
随着城市化进程不断加速,城市地下空间的开发越来越受到人们的关注。大规模的新建地铁盾构隧道在掘进过程中,不可避免地出现穿越既有建筑桩基的工况[1]。隧道开挖引起的地层应力释放,造成既有桩基周围地层的应力场发生改变,从而导致桩基产生附加变形和附加内力,并造成桩基承载力发生变化,而过度的承载力损失将危及上部结构安全[2-3]。因此,准确评价隧道近接施工对既有摩擦桩基的影响是当前城市地下空间开发所面临的较为紧迫的问题。
当前,针对隧道施工引起的邻近桩基响应问题,学者们基于现场实测、模型试验[4-5]、数值模拟和理论解析等方法已取得了较为丰硕的成果。如李雪等[6]针对饱和砂土地区盾构隧道超近接高铁桥墩摩擦桩的工程问题,对既有高铁桥墩桩基变形进行了现场测试;Loganathan等[7]和孙庆等[8]采用离心机试验研究了隧道开挖对邻近桩基的瞬态和长期影响,指出隧道开挖影响下邻近桩基最大横向变形和附加弯矩均发生在隧道轴线附近。
隧道施工引起的邻近桩基响应问题,实质为隧道−土体−桩基之间相互作用。学者们已构建了十分丰富的解析计算模型。如ZHANG等[9]基于Pasternak地基模型及桩侧土体影响,研究了被动桩在盾构掘进下水平向变形响应;程康等[10]等进一步考虑桩基的横向剪切性能,将桩简化成Pasternak地基上的Timoshenko梁,进一步研究了盾构掘进引起的桩基水平变形;冯国辉等[11-13]分别将既有桩基视为Winkler地基和Pasternak地基中的Euler-Bernoulli梁,分析了盾构掘进下邻近桩基的横向附加响应。
本文在既有研究的基础上,结合隧道开挖前后地层应力场的变化,建立盾构近接施工下邻近摩擦桩基承载力解析计算模型,分析了隧道施工对既有摩擦桩基承载特性的影响机制,并进一步基于三维有限元分析论证了隧道洞周地层加固对既有桩基变形控制的有效性。
1. 盾构近接施工扰动下摩擦桩基承载特性分析
1.1 基本思路
盾构施工对邻近摩擦桩基承载特性的影响,实质在于隧道施工会引起周围土体应力场重分布,而导致摩擦桩侧法向正应力和切向剪应力发生变化,继而影响桩侧摩阻力。因此对于盾构施工扰动下邻近摩擦桩基承载力分析的基本思路:
(1)隧道施工前,地层内任一点的应力状态由围岩自重生成竖向应力σv0、水平向应力σh0,将该值作为桩侧地层应力状态初始值,可求得初始桩侧阻力Qs0。
(2)计算隧道施工后地层内任一点的竖向应力σv、水平应力σh和剪应力τxy,根据此时的应力状态计算隧道施工后的桩侧阻力Qs。
(3)结合步骤(1),(2),求得摩擦桩基承载力变化值ΔQs=Qs0−Qs。
1.2 解析模型的建立及基本假定
为获得盾构近接施工扰动下临近摩擦桩基的承载特性,建立如图 1所示的隧道-地层-桩基相互作用模型。图 1中x,y为弹性半空间任意一点的直角坐标,r,θ为相应的极坐标;h为隧道轴线埋深;σr,σθ,τrθ分别为任意一点土体径向应力、环向应力和剪应力;l为桩基长度;x0为桩基轴线与隧道轴线的水平间距;r0为隧道开挖半径。
在理论计算时作如下基本假定:
(1)为避免地层应力场计算时的奇异,隧道覆径比大于1.5。
(2)忽略土体与衬砌之间的摩擦力,即接触面处的剪应力满足:
τ|r=r0=0。 (1) (3)采用间隙参数w表征盾构开挖引起的等效地层损失,计算方法为
w=Vlπr022πr0=12Vlr0。 (2) 式中:Vl为地层损失率。
1.3 初始桩侧阻力Qs0
任在原岩应力状态下,桩侧任意一点的应力状态为
σv0=γ(h−y),σh0=kσv0,τxy0=0。} (3) 式中:γ为土体有效重度;k为侧压力系数。
根据摩擦桩基承载理论,此时桩基承载力为
Qs0=π d∫lfs0(y)dy=π d∫l(τxy0+ksσh0)dy。 (4) 式中:d为桩基直径;fs0为初始状态下的桩侧摩阻力分布函数;ks为桩侧界面摩擦系数,ks=tanδφ;δφ为界面摩擦角。
1.4 开挖扰动后桩侧阻力Qs
以干砂地层为例,BOBET基于应力函数法给出了此类地层中隧道开挖后任意一点的应力,即
σr=a0r2−12γh(1+k)+[14γr(k+3)−2c′1r−3+12c1r−1]sinθ+[12γh(1−k)−6a′2r−4−2b′2r−2]⋅cos2θ+[−14γr(1−k)−12c′3r−5−6d′3r−3]sin3θ, (5a) σθ=−a0r2−12γh(1+k)+[14γr(3k+1)+2c′1r−3−12c1r−1]sinθ+[−12γh(1−k)+6a′2r−4+2b′2r−2]⋅cos2θ+[14γr(1−k)+12c′3r−5+6d′3r−3]sin3θ, (5b) τrθ=[14γr(1−k)+2c′1r−3]cosθ−[12γh(1−k)+6a′2r−4+2b′2r−2]sin2θ−[14γr(1−k)−12c′3r−5−6d′3r−3]cos3θ。 (5c) 式中:a0,a′2,b′2,c1,c′1,c′3,d′1,d′3为常数,可根据问题的边界条件确定,其他符号意义同前。
进一步,可通过坐标变换将地层中一点应力状态(式(5))由极坐标系转换到直角坐标系,即
σx=σr+σθ2+σr−σθ2cos2θ−τrθsin2θ, (6a) τxy=σr−σθ2sin2θ+τrθcos2θ。 (6b) 将式(5)代入式(6),结合摩擦桩基承载理论,得到隧道开挖扰动后的桩基承载力,即
Qs=πd∫lfs(y)dy=πd∫l(τxy+ksσx)dy。 (7) 式中:fs为隧道开挖后桩侧摩阻力分布函数;σ′x,τ′xy分别表示直角坐标系下任意一点的有效水平应力和剪应力。
1.5 隧道开挖对摩擦桩基承载力的影响分析
以某地铁盾构隧道连续侧穿老旧建筑群施工为工程背景,基于上述隧道开挖前后摩擦桩基承载力理论分析,进行隧道开挖对邻近摩擦桩基承载特性的影响分析。选取典型工况:既有桩基及开挖隧道位于砂土地层中,地层重度约为18 k N/m3,侧压力系数约为0.32,根据地勘资料,计算时地层弹性模量和泊松比分别取30 MPa,0.3;盾构隧道开挖半径为3.14 m,轴线埋深为15 m,隧道衬砌采用C50钢筋混凝土管片,管片厚度为0.3 m,弹性模量和泊松比分别为34500 MPa,0.2;既有建筑物桩基直径为0.5 m。为反映隧道施工对摩擦桩基承载能力的影响,定义摩擦桩基承载力影响因子RQ:
RQ=QsQs0。 (8) RQ>1表示承载力提高,RQ<1表示承载力损失,当RQ=1时表示承载力无变化。
如图 2所示为桩隧横向间距为1.5r0,地层损失率为0.8%时,摩擦桩基承载力影响因子随桩长的变化曲线。由图可见:承载力影响因子变化曲线随桩长的增大而呈漏斗状分布,当l/h<1.0时,承载力影响因子随桩长的增大而迅速减小,而当l/h>1.0时,承载力影响因子随桩长的增大而又逐渐增大,并逐渐趋近于1.0;桩隧相对埋深越趋近1.0,承载力影响因子变化梯度越大,承载力损失越多;当桩基端部与隧道轴线相平时,即l/h=1.0时,其承载力损失最多,在本算例中此时承载力影响因子RQ=0.17,桩基承载力损失82.9%。
如图 3所示为当桩隧横向间距为3.0 r0,桩隧相对埋深l/h = 1.0时,摩擦桩基承载力影响因子随地层损失率的变化曲线。由图 3可见:摩擦桩基承载力损失因子与地层损失率呈线性相关关系,随着地层损失率的增大,摩擦桩基承载力影响因子随之线性减小;在本算例中,当地层损失率仅为0.4%时,承载力影响因子为0.98,而当地层损失率增大至1.8%时,摩擦桩基承载力影响因子减小至0.15,变化84.7%。
2. 盾构侧穿摩擦桩基建筑物控制技术
盾构隧道施工对周围岩土体造成扰动,应力平衡状态改变造成围岩应力释放,导致临近桩基产生沉降变形。综合前述计算分析,从地层加固措施入手减小盾构施工对邻近桩基的影响,通过注浆加固地层的方式减小盾构施工扰动在地层中的传播。
2.1 计算模型及计算参数
所建立的盾构-围岩-建筑物的三维数值模型如图 4所示,模型尺寸为长60 m (沿隧道开挖方向),宽120 m (垂直隧道开挖方向),高50 m (沿地层深度方向)。模型侧面约束水平位移,底部为固定边界,上表面为自由边界,共划分56972个单元,41673个节点。数值模型中盾构隧道管片外径6.0 m,内径5.4 m,厚度300 mm,宽度1.2 m,采用C50混凝土,考虑衬砌接头对管片刚度的削弱效应,管片刚度需折减15%。
模型中地层采用摩尔库仑模型,盾壳、隧道衬砌、同步注浆区(即后文所述的等代层)、袖阀管注浆加固区、建筑物梁板柱、桩基础均采用线弹性模型。地层、等代层、袖阀管注浆加固区采用实体单元模拟,盾壳、隧道衬砌、建筑楼板采用板单元模拟,建筑框架梁、立柱、建筑桩基采用梁单元模拟。参数取值见表 1。
表 1 结构参数取值Table 1. Values of structural parameters名称 密度
γ/(kN·m-3)弹性模量
E/GPa泊松比 盾壳 78.5 210 0.30 隧道衬砌 25.0 29.3 0.20 建筑框架
(梁、板、柱)25.0 20.4 0.20 桩基础 25.0 24.0 0.20 等代层 20.0 0.06 0.25 袖阀管注浆加固区 22.0 0.25 0.25 2.2 盾构隧道施工过程模拟
(1)同步注浆
在盾构机施工过程中,在盾尾管片脱出后,由于盾构机刀盘的超挖间隙、装配管片所需的操作间隙、盾构机壳的厚度3个因素将会在盾构管片和土体间产生盾尾间隙。为了减小施工扰动,控制地层变形,在盾构掘进过程中将会采取同步注浆对盾尾间隙进行填充。但是实际工程中,由于地层向盾尾间隙道的移动程度、注浆体注浆后的分布情况、隧道围岩壁面受扰动的范围和程度等情况相对复杂,因此,在建模过程中通常将注浆体简化为均质的、弹性的、等厚的等代层处理,等代层的厚度δ可为
δ=μΔ。 (9) 式中:Δ为盾尾间隙的理论计算值;μ为折减系数,取值范围为0.7~2.0,对于极软土层取上限值,对于硬土层取下限值。
根据盾构穿越地层条件及施工扰动范围及程度取等代层厚度为100 mm,弹性模量为60 MPa,泊松比为0.25,采用线弹性体模拟,并根据实际施工情况取注浆压力为300 kPa。
(2)土仓压力
对于土压平衡盾构机而言,掌子面的稳定一般依靠于密封土仓的土仓压力平衡掌子面前方的水土压力,并且考虑到盾构刀盘切削土体使土体剥落而导致土仓压力不稳定的情况,为了防止掌子面坍塌,根据施工经验,一般会预留10~20 kPa的预备压力,所以土仓压力一般按下式确定:
P0=Pc+Pw+p。 (10) 式中:P0为土仓压力;Pc为土压力;Pw为水压力;P为预备压力。其中Pc=∑K0γizi,Pc=γwh,K0为土体静止侧压力系数,γi为第i层土的重度,zi为第i层土的厚度,γw为水的重度,h为地下水位都盾构刀盘的高度。
通过计算取土仓压力为200 kPa均布施加于开挖面土体上。
(3)盾壳与地层摩擦力
盾构掘进过程,盾构机和地层发生相对运动而产生摩擦力,会使得盾构机周围土体由于剪切作用而向前运动,沿伸到地表则表现为地表隆起或者沉降。盾壳和地层间的摩擦力计算公式为
f=μs∑γizi。 (11) 式中:μs为盾构机与地层的摩擦系数,通过计算得到的摩擦力为153 kPa。
3. 计算结果分析
3.1 地层位移
在整个开挖过程中,注浆加固区附近的土体变形明显变小,当双线隧道开挖完成后,对于未进行注浆加固的工况来说,沉降发生区域在左线隧道和右线隧道间呈长条带状分布,而注浆加固后,沉降发生区域沿开挖方向呈现“两头大,中间小”分布,凸向隧道左线。选取隧道轴线纵向30 m处的横断面作为分析断面,对加固地层前后的地表沉降进行分析,计算结果如图 5所示。
由图 5可知,地表沉降曲线呈“U”字形,符合“Peck”曲线分布,并与实测数据有较好的拟合度。当右线隧道开挖完成时,地表沉降最大值出现在右线隧道正上方,当左线隧道开挖完成后,地层受二次开挖扰动,沉降值变大,沉降槽宽度加宽,最低点往左移动。对于距起始开挖面30 m处的横断面而言,处于加固区,受右线隧道开挖影响减小,沉降最大值出现在左线隧道正上方附近,注浆加固前后沉降明显变小,最大值分别为14.8,10.7 mm,降低27.8%。
3.2 桩基变形
选取距离起始开挖面最近的一排3根桩进行分析,距离隧道由近到远分别为桩基P1,P2,P3。由图 6可知,双线隧道开挖完后,由于土体的卸荷作用,桩基发生挠曲变形,桩基顶部向隧道方向移动,桩基底部远离隧道移动。离隧道越远,桩基水平位移越小,由于上部结构的约束作用,各单桩桩顶位移基本相同。在注浆加固后,桩基水平位移明显减小,桩顶最大水平位移由2.6 mm减少至1.4 mm,降低46%。
如图 7所示为加固地层和未加固地层工况下左右线贯通后桩基竖向位移沿桩身埋深分布曲线。由图 7可见:桩身沉降沿桩身基本相同,离隧道越远,沉降越小,桩基P1,P2沉降受隧道开挖影响较大,在注浆加固后,桩基沉降由5.7 mm减少至2.6 mm,降低54%。
4. 结论
针对盾构近接施工对邻近摩擦桩基的影响问题,分别建立了隧道开挖扰动下摩擦桩基承载力解析计算模型和三维数值分析模型,探明了隧道开挖对既有摩擦桩基承载力的影响机制,并论证了地层加固对桩基变形控制的有效性,得到3点结论。
(1)基于摩擦桩基承载力理论和隧道开挖前后地层应力场的分析,建立了隧道开挖影响下摩擦桩基承载力计算方法,并定义摩擦桩基承载力影响因子分析了不同工况下隧道施工对桩基承载力的影响。
(2)隧道施工对不同长度的摩擦桩基影响不同,当桩端埋深与隧道轴线埋深一致时,桩基承载力损失最为严重,而当桩长远大于隧道埋深时其承载力损失较小;隧道施工时的地层损失率越大,桩基承载损失越严重,二者呈正相关。
(3)建立了盾构侧穿浅层摩擦桩基建筑物施工的三维有限元模型,分析了盾构侧穿桩基施工时的桩基沉降变形规律。计算结果表明在采取洞周地层加固措施后,隧道围岩注浆加固后土体变形明显变小,可有效减少临近桩基水平位移和桩身的沉降。
-
图 1 深围井及现场取样图
Figure 1. Formation of in-situ test well[13] and photos of field sampling of soil materials
表 1 试验方案表
Table 1 Summary of test schemes
土样类型 固结应力比Kc 围压σ3/kPa 原状样 1.0,1.5,2.0 200,400,800,1600 重塑样 1.0 200,400,800,1600 说明:每个围压条件下4个试样,振动破坏周次分别在10,20,50,100左右。 -
[1] 余挺, 邵磊. 含软弱土层的深厚河床覆盖层坝基动力特性研究[J]. 岩土力学, 2020, 41(1): 267-277. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202001032.htm YU Ting, SHAO Lei. Study of dynamic characteristics of dam foundation on deep riverbed overburden with soft soil layer[J]. Rock and Soil Mechanics, 2020, 41(1): 267-277. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202001032.htm
[2] HAN B, ZDRAVKOVIC L, KONTOE S, et al. Numerical investigation of the response of the Yele rockfill dam during the 2008 Wenchuan earthquake[J]. Soil Dynamics and Earthquake Engineering, 2016, 88: 124-142. doi: 10.1016/j.soildyn.2016.06.002
[3] WANG W B, HÖEG K, ZHANG Y B. Design and performance of the Yele asphalt-core rockfill dam[J]. Canadian Geotechnical Journal, 2010, 47(12): 1365-1381. doi: 10.1139/T10-028
[4] 李鹏, 杨兴国, 薛新华, 等. 瀑布沟心墙堆石坝地震响应分析[J]. 世界地震工程, 2015, 31(1): 217-223. https://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201501032.htm LI Peng, YANG Xingguo, XUE Xinhua, et al. Earthquake response analysis of Pubugou corewall rockfill dam[J]. World Earthquake Engineering, 2015, 31(1): 217-223. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201501032.htm
[5] YU X, KONG X J, ZOU D G, et al. Linear elastic and plastic-damage analyses of a concrete cut-off wall constructed in deep overburden[J]. Computers and Geotechnics, 2015, 69: 462-473. doi: 10.1016/j.compgeo.2015.05.015
[6] 杨正权, 刘启旺, 刘小生, 等. 厄瓜多尔CCS水电站深厚覆盖层火山灰沉积土动强度特性试验研究[J]. 水利学报, 2014, 45(增刊2): 161-166. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB2014S2026.htm YANG Zhengquan, LIU Qiwang, LIU Xiaosheng, et al. Study on dynamic strength characteristics tests for deposited soils of volcanic ash[J]. Journal of Hydraulic Engineering, 2014, 45(S2): 161-166. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB2014S2026.htm
[7] 杨玉生, 刘小生, 李小泉, 等. 固结应力状态对超深厚覆盖层深埋砂土动强度参数的影响[J]. 水利学报, 2016, 47(4): 518-526. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201604007.htm YANG Yusheng, LIU Xiaosheng, LI Xiaoquan, et al. Effects of effective confining stresses on cyclic resistance ratio of deep buried sands in deep alluvial soils[J]. Journal of Hydraulic Engineering, 2016, 47(4): 518-526. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201604007.htm
[8] EL-SEKELLY W, DOBRY R, ABDOUN T. Assessment of state-of-practice use of field liquefaction charts at low and high overburden using centrifuge experiments[J]. Engineering Geology, 2023, 312: 106952. doi: 10.1016/j.enggeo.2022.106952
[9] 周燕国, 谭晓明, 陈捷, 等. 易液化深厚覆盖层地震动放大效应台阵观测与分析[J]. 岩土工程学报, 2017, 39(7): 1282-1291. doi: 10.11779/CJGE201707015 ZHOU Yanguo, TAN Xiaoming, CHEN Jie, et al. Observations and analyses of site amplification effects of deep liquefiable soil deposits by geotechnical downhole array[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(7): 1282-1291. (in Chinese) doi: 10.11779/CJGE201707015
[10] OKA L G, DEWOOLKAR M M, OLSON S M. Liquefaction assessment of cohesionless soils in the vicinity of large embankments[J]. Soil Dynamics and Earthquake Engineering, 2012, 43: 33-44. doi: 10.1016/j.soildyn.2012.06.012
[11] JANA A, STUEDLEIN A W. Dynamic in situ nonlinear inelastic response of a deep medium dense sand deposit[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2021, 147(6): 04021039. doi: 10.1061/(ASCE)GT.1943-5606.0002523
[12] WANG L, YANG Z Q, ZHAO J M, et al. Transmitting characteristics of seismic motion in super-deep overburden layer ground[J]. Shock and Vibration, 2021, 2021: 8898012.
[13] YU T, JIANG B N, HE S B, et al. Study on the seepage control of ultra-deep in situ test wells in deep overburden of riverbeds[J]. Arabian Journal of Geosciences, 2022, 15(9): 834. doi: 10.1007/s12517-022-10111-8
[14] WANG L, YANG Z, ZHAO J, et al. Seismic response analysis of earth-rock fill dam on deep overburden under viscoelastic boundary condition[J]. Revista Internacional De Métodos Numéricos Para Cálculo y Diseño En Ingeniería, 2021, 37(1): 1-8.
[15] 邵磊, 余挺, 徐强. 无限元在超深厚覆盖层坝基动力分析中的应用[J]. 地下空间与工程学报, 2017, 13(6): 1537-1543. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201706016.htm SHAO Lei, YU Ting, XU Qiang. Dynamic analysis on super-deep overburdened layer foundation based on infinite element method[J]. Chinese Journal of Underground Space and Engineering, 2017, 13(6): 1537-1543. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201706016.htm
[16] 罗启迅, 丁庆, 李建国, 等. 用于批量快速获取原状土样的施工方法: CN114739729A[P]. 2022-07-12. LUO Qixun, DING Qing, LI Jianguo, et al. Construction Method for Rapidly Obtaining Undisturbed Soil Samples in Batches: CN114739729A[P]. 2022-07-12. (in Chinese)
[17] 丁庆, 罗启迅, 田先忠, 等. 用于批量获取原状土样的砂浆材料及取样施工方法: CN114656229B[P]. 2023-04-28. DING Qing, LUO Qixun, TIAN Xianzhong, et al. Mortar Material for Batch Acquisition of Undisturbed Soil Samples and Sampling Construction Method: CN114656229B[P]. 2023-04-28. (in Chinese)
[18] 水电水利工程土工试验规程: DL/T 5355—2006[S]. 北京: 中国电力出版社, 2006. Code for Soil Tests for Hydropower and Water Conservancy Engineering: DL/T 5355—2006[S]. Beijing: China Electric Power Press, 2006. (in Chinese)
[19] SEED H B, IDRISS I M, LEE K L, et al. Dynamic analysis of the slide in the lower San Fernando Dam during the earthquake of February 9, 1971[J]. Journal of the Geotechnical Engineering Division, 1975, 101(9): 889-911. doi: 10.1061/AJGEB6.0000195