• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

颗粒表面粗糙度对材料小应变动力特性的影响

张涛, 吴健, 魏骁, 杨仲轩

张涛, 吴健, 魏骁, 杨仲轩. 颗粒表面粗糙度对材料小应变动力特性的影响[J]. 岩土工程学报, 2024, 46(8): 1783-1790. DOI: 10.11779/CJGE20230390
引用本文: 张涛, 吴健, 魏骁, 杨仲轩. 颗粒表面粗糙度对材料小应变动力特性的影响[J]. 岩土工程学报, 2024, 46(8): 1783-1790. DOI: 10.11779/CJGE20230390
ZHANG Tao, WU Jian, WEI Xiao, YANG Zhongxuan. Effects of surface roughness of particles on small-strain dynamic properties of granular materials[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(8): 1783-1790. DOI: 10.11779/CJGE20230390
Citation: ZHANG Tao, WU Jian, WEI Xiao, YANG Zhongxuan. Effects of surface roughness of particles on small-strain dynamic properties of granular materials[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(8): 1783-1790. DOI: 10.11779/CJGE20230390

颗粒表面粗糙度对材料小应变动力特性的影响  English Version

基金项目: 

浙江省“尖兵”“领雁”研发攻关计划项目 2022C03180

国家自然科学基金青年科学基金项目 52108351

浙江省建设厅项目 2022K073

浙江交工协同创新联合研究中心项目 ZDJG2021001

详细信息
    作者简介:

    张涛(1998—),男,硕士研究生,主要从事土动力学方面的研究工作。E-mail: zhangtaot@zju.edu.cn

    通讯作者:

    魏骁, E-mail: weixiaos@zju.edu.cn

  • 中图分类号: TU435

Effects of surface roughness of particles on small-strain dynamic properties of granular materials

  • 摘要: 土体的小应变剪切模量和阻尼比是表征土体动力学特性的重要参数,不仅受到土体密实度和应力状态的影响,还受到土体颗粒级配、形状等颗粒特征的影响。颗粒表面粗糙度是重要的土体颗粒特征之一,然而关于颗粒表面粗糙度对土体小应变动力特性影响的研究较为匮乏。利用能量注入式虚拟质量共振柱设备,系统地测试了颗粒表面粗糙度不同的玻璃珠所成试样的小应变剪切模量和阻尼比;采用三维干涉显微镜测量了玻璃珠的表面粗糙度,并量化表征粗糙度对试样小应变剪切模量和阻尼比的影响。试验结果表明,在相同孔隙比和有效应力条件下,试样的小应变剪切模量随颗粒表面粗糙度增大而减小,而小应变阻尼比受颗粒表面粗糙度影响的规律不明显。研究结果表明,当材料的颗粒形状、级配等因素相近时,颗粒表面粗糙度对材料小应变剪切模量的影响不应被忽略。
    Abstract: The small-strain shear modulus and damping ratio are the important dynamic properties of granular soils, which are influenced by the packing density and stress states of the soil element as well as the particle characteristics including gradation, shape of particles, etc. As one of the important particle characteristics, however, the surface roughness of particles is often ignored by the previous investigations regarding the small-strain dynamic properties of granular soils. In this study, the small-strain dynamic properties of specimens of glass beads with different surface roughnesses are measured using the energy injecting virtual mass resonant column system. The surface roughness of the glass beads is quantified with the assistance of the three-dimensional interference microscope. The experimental results show that the small-strain shear modulus decreases with the increasing surface roughness of particles, while there is no obvious effects of surface roughness on the small-strain damping ratio. Clearly, the effects of surface roughness on the small-strain dynamic properties cannot be ignored particularly when other influencing factors, such as particle shape and gradation, are similar.
  • 土壤液化是国际上具有普遍意义的典型震害现象,是土木工程防灾减灾中的复杂难题之一。真实场地的液化响应规律及特征,仍是当今液化研究中没有解决好的突出问题之一,尚需深入研究[1]

    试验是液化研究中不可或缺的手段。然而,迄今尚没有一种试验,能够完全复现真实场地的液化响应过程及行为特征。室内液化试验,具有大量且突出的优点。可是,无论动三轴等土单元试验,还是振动台或动力离心机等模型试验,目前都不能同时、全面地对可液化土所处的现场条件和影响因素进行模拟,与真实液化场地存在差距。作为室内试验的一种补充,直接开展现场液化试验或现场液化监测,逐渐成为当今土壤液化研究的必需手段[2]。采用室内和现场液化试验两种不同的研究手段,通过不断地“相互对比、验证、补充”,可以进一步完善、修正当前的液化分析方法。

    相比室内液化试验方面的丰富成果,人工激振下现场液化试验的研究成果尚少,目前国际上总体处于探索阶段[2]。中国虽有现场液化试验方面的研究成果并取得一些开创性的进展[3-7],但是总体仍处于起步阶段,亟待进一步发展。

    值得注意的是,目前国内外已有的人工激振下现场液化试验,无论是基于原状土还是重塑土,可液化土的上覆压力基本不超过20 kPa,绝大多数在10 kPa以内,大致对应1~2 m的土层埋深。根据对1962—2010年中国大陆地区地震液化实例的统计,液化土层的埋深集中在0~6 m,约占总数的80%,其中2~4 m最为集中。而其地下水位埋深集中在0~3 m,约占总数的90%,其中1~2 m最为集中[8]。以此粗略推算,真实地震中液化土层所对应的上覆压力范围应集中在0~120 kPa,最为集中的均值约50 kPa左右。众所周知,上覆压力(或埋深)是控制土壤液化响应的主要因素之一。显然,当今开展的现场液化试验中可液化土的上覆压力较低,与真实液化场地存在较大差距。而且,由于试验操作的扰动,有效围压50 kPa以下的室内土单元试验成果相对很少。这也导致了已有现场液化试验成果难以和室内液化试验,尤其是动三轴等土单元试验,在相近的动、静应力水平下进行对比分析。可以说,过低的上覆压力,使得目前通过现场液化试验研究土壤液化问题的价值大打折扣。

    现有的振动台或动力离心机液化试验中,模型土的上覆压力和真实液化场地也会存在差距,往往通过堆载重物的方式来改善[9-10]。动力离心机试验,当应力水平无法满足时,在ng重力环境下,采用比土的密度大几倍的材料(如钢材)作为堆载重物,通过较小的堆载高度就会达到理想效果。然而,在振动台或现场试验等1g重力环境下,当所需上覆压力较大时,通过堆载重物的方式并不理想,甚至不可行。毕竟,试验中堆载重物的高度不可能无限制增加,而且实际加压效果、可操作性等都会遇到问题。若采用传统的“千斤顶+反力架”的方式,似乎非常容易实现对模型土施加上覆压力,且能完全满足液化试验需求的压力上限值。然而,对于液化试验而言,这种传统的加压方式并不可行。因为液化发生前,土体一定会率先出现沉降变形,微小的沉降都会导致通过千斤顶施加的荷载无法传递,造成上覆压力“瞬间消失”。因此,需专门发展试验技术,才能开展考虑上覆压力影响的现场液化试验。

    鉴于此,笔者在已有现场液化试验研究基础上[6-7],围绕能够开展“多工况、可重复、易操作、低成本”试验的基本设想,专门研制施加上覆压力的设备,发展考虑上覆压力影响的现场液化试验技术,设计完成2次不同上覆压力(25,50 kPa)的饱和砂、干砂对比试验,分析其动力响应规律及特征,并与已有认识对比,验证试验技术的可靠性。本文研究成果,可为当今人工激振下现场液化试验的发展提供技术借鉴,也为在1g重力环境下,以更接近真实场地的模拟条件开展土壤液化问题研究,提供一种简单实用、可操作性强的试验新方法。

    为评价考虑上覆压力影响的现场液化试验是否成功实施,本文首先提出试验孔压比的概念,并阐述其必要性和合理依据。

    由于现场液化试验是在1g重力环境下进行,上覆压力必须通过特定技术手段施加给试验土体。在动荷载开始前,预先施加给土体的上覆压力能够较容易地保持平稳和恒定状态。但是,在动荷载作用下,人为施加的上覆压力会发生很大变化,尤其是采用机械设备以“接触力传递”方式施加的上覆压力。理论上,液化的土体会软化和变形,通过机械设备一直维持初始施加的上覆压力的平稳,几乎不可能,因为液化的土体无法提供良好的反力支撑。但是,在液化发展阶段,尤其是超静孔压急剧上升时,如果预先施加的上覆压力下降(损失)过快、过大,则必然会对现场液化试验结果产生较大影响,甚至失去考虑和控制上覆压力的意义。因此,在现场液化试验中,施加给土体的上覆压力应保持大致的稳定。只有预先施加的上覆压力损失有限,以此得到的现场液化试验结果才有价值。而且,需要对现场液化试验中上覆压力的损失程度,进行定量评价和判断。

    一方面,根据已有试验结果和液化认识,当孔压比达到0.4~0.6左右,是液化变形及对地震动产生影响等动力响应特征的阈值[6, 11-12]。另一方面,从液化机理上,当土壤液化的条件均能充分满足时,超静残余孔压的发展是不可能大于其上覆压力的。也就是说,如果设备预先施加的上覆压力下降(损失)过快、过大,超静孔压一定不会发展过高。

    因此,在土壤液化的条件均能满足的情况下,只要现场液化试验中的超静孔压与预先施加的初始上覆压力的最大比值能够达到0.4~0.6及以上,就可获取该上覆压力作用下土壤液化的主要响应特征。本文将该最大比值,定义为试验孔压比,以区别通常描述液化特征的孔压比。因此,试验孔压比,可初步判断预先施加的上覆压力是否大致稳定,从而评价可考虑上覆压力影响的现场液化试验是否成功实施(有价值)。

    当然,当土体抗液化强度较高或动荷载较小时,液化发展程度本身有限,就不会发生较大的土体软化或变形,预先施加的上覆压力自然也会保持较平稳状态。此时,则无需采用试验孔压比进行评价。

    本文现场液化试验,所用的砂土全部采用福建标准砂。通过室内土工试验实测,其物理指标如表 1所示,其颗粒级配曲线如图 1所示。

    表  1  试验砂物理指标
    Table  1.  Physical indexes of sands used in tests
    最大干密度ρdmax/(g·cm-3) 最小干密度ρdmin/(g·cm-3) 有效粒径d10/mm 特征粒径d30/mm 控制粒径d60/mm 不均匀系数Cu
    1.953 1.538 0.309 0.475 0.889 2.877
    下载: 导出CSV 
    | 显示表格
    图  1  试验用砂颗粒级配曲线
    Figure  1.  Grain-size distribution curves of sands used in the tests

    本文现场液化试验所用的动力加载设备,与笔者以往开展的现场液化试验相同[6-7]。该设备是自主研发,主要有由激振装置、刚体基础、变频控制器、电动机等组成,能在试验场地中激励产生水平向的等幅正弦波地震动。加速度传感器型号为LC0113M,量程±2g,灵敏度2500 mv/g左右。孔压传感器型号为HC-25,量程0~50 kPa,灵敏度0.2%FS;土压力传感器型号HC-16,量程0~150 kPa,灵敏度0.2%FS。位移传感器型号为HCW-20,量程0~15 cm,灵敏度0.1%FS。动态数据采集设备,采用日本东京测器研究所生产的TMR多功能记录仪,最高数据采样率可达100 kHz。部分试验设备及传感器如图 2所示。由图可知,所用的加速度、孔压、土压力传感器体积小,能够满足土体“点测量”的基本要求。

    图  2  现场试验用的部分试验设备及传感器
    Figure  2.  Photos of experimental apparatus and sensors used in in-situ tests

    施加上覆压力的设备,是本文现场液化试验能否成功实施的关键硬件。笔者自主研发了适用于现场液化试验的施加上覆压力的设备。该设备主要由钢筋混凝土反力架、液压千斤顶、加载钢板及持压装置组成,其示意图及现场实景如图 3所示。该设备,是在传统的“反力架+千斤顶”模式的基础上,增加了自主研发的持压装置,解决了因土体沉降导致上覆压力“瞬间消失”的问题。持压装置的工作原理,是通过合理刚度的弹簧的回弹来补偿沉降,从而实现“伺服”控制。通过该设备,现场液化试验所需的上覆压力,不但可以被定量施加,而且能够在液化发展阶段保持大致稳定的状态。该设备是无源工作,操作简单,尤其适合野外试验环境。

    图  3  施加上覆压力设备的示意图及现场实景
    Figure  3.  Sketch and photo of apparatus for applying overburden pressure in in-situ liquefaction tests

    本文试验的主要目的,是在已有重塑土现场液化试验基础上,围绕“多工况、可重复、易操作、低成本”的基本设想,发展考虑上覆压力影响的现场液化试验技术并验证其可靠性。考虑到本文试验的探索性,并为后续针对性地开展具体液化问题的研究奠定认识基础,本次试验方案设计遵循以下基本原则:①设置同工况的干砂试验,作为对比参照;②仅把上覆压力作为试验变量,其他影响液化因素暂不考虑。如前文所述,考虑当前上覆压力基本现状,优先设计上覆压力分别为25,50 kPa的两次工况的现场液化试验。

    基于上述设想,在试验场地直接开挖2个相同尺寸的试坑,然后在坑内分别制备干砂、饱和砂。人工激励产生的水平向震动,通过刚体基础周围的土层传递,以“地震动”的形式直接对试坑内的砂土施加动荷载,因此不存在振动台试验的动力边界效应问题。由于试坑尺寸的大小不会对试验技术可行性产生实质影响,经综合考虑,试坑的尺寸确定为40 cm×40 cm×50 cm。试坑和刚体基础的空间相对位置关系如图 4所示。

    图  4  现场液化试验的试坑布置图
    Figure  4.  Sketch of test pits used in in-situ liquefaction tests

    本文设计试坑内砂土的初始相对密实度Dr为40%。砂土的制备技术及操作方法,与笔者之前进行的现场液化试验相同[6-7]。饱和砂,仍采用水沉法制备,其内部的传感器布置方案如图 5所示。干砂内部,除不放置孔压传感器外,其他布置均与饱和砂相同。考虑到试验的探索性,砂土内部的传感器均采用“备份”原则布置。即在砂土中心和左边的不同深度测点,逐一对应,布置相同类型的传感器。在干砂试坑一侧的地表,与砂土表面中心对应(均距离基础边缘70 cm),单独埋设1个加速度传感器,其实测加速度作为试验的输入地震动。位移传感器固定于反力架上端,通过加载钢板来测量砂土的地表沉降。

    图  5  饱和砂土的传感器布置图
    Figure  5.  Sketch of sensor distribution for saturated sands

    为表述方便,本文上覆压力为25 kPa工况的试验命名为Test-L试验,上覆压力为50 kPa工况的试验命名为Test-H试验。实际中,施加至设计上覆压力后,砂土均有一定程度的沉降。根据沉降值估算,在Test-L试验中,干砂的密实度变为49%,饱和砂密实度变为46%。在Test-H试验中,干砂的密实度为54%,饱和砂变为50%。砂土密实度的变化程度,均在可接受的范围。动荷载频率,与笔者之前的现场液化试验相同,仍为15 Hz。输入地震动的加速度幅值为0.11g左右,数据采样频率为100 Hz。两个工况现场试验分别加载。

    根据笔者观察,加载期间,无论Test-L试验,还是Test-H试验,饱和砂土均出现明显的喷砂冒水现象。图 6为Test-H试验后的饱和砂土表面实景,可以看出大量的砂堆积在加载钢板的一角,是加载期间喷出。相比笔者以前的现场液化试验,本文试验中饱和砂土的喷水高度、喷砂量均有大幅增加。这应该是,较大的上覆压力,导致砂土液化时超静孔压明显升高的结果。试验宏观现象初步表明,两次现场液化试验均取得了成功。加载期间,除去Test-L试验干砂中的土压力时程数据因电源线意外故障导致采集失败外,其他测点均成功采集到了数据。

    图  6  现场液化试验后的饱和砂实景(Test-H工况)
    Figure  6.  Photo of saturated sands in in-situ liquefaction tests after dynamic loadings (case of Test-H)

    图 7分别给出了两个工况现场试验中试坑外侧的地表加速度时程。经对比,试坑内砂土底部的加速度与其一致,限于篇幅不再给出,并将该数据视为试验的地震动输入加速度。可以看出,两次试验的输入加速度的幅值几乎一致,均在0.11g~0.12g左右。

    图  7  试验的输入加速度时程
    Figure  7.  Time histories of input acceleration in tests

    超静孔压是液化发展程度最直接的反映和标志。由于试坑内饱和砂土中的孔压传感器距离真实地表仅15 cm,如前文所述,倘若没有额外施加的上覆压力,超静的残余孔压是不可能超过1.5 kPa的(即大致15 cm高度饱和砂的有效自重)。图 8给出了两次现场液化试验饱和砂中心测点的超静孔压时程。结合图 7的地表加速度时程可知,在地震动开始的最初3 s内(约45个荷载循环),超静孔压迅速发展并达到最大值,之后孔压开始逐渐消散,这与以往的现场液化试验孔压响应规律和认识一致[6-7]。所不同的是,超静孔压到达最大值后,便表现为消散趋势且速率较快,这应该与实际排水路径过短(此处仅15 cm)有关。从图 8中可以看出,Test-L试验中超静孔压的最大值约为17 kPa,与其初始上覆压力23 kPa的比值,即试验孔压比,为0.75。Test-H试验中超静孔压的最大值约为25 kPa,与初始上覆压力50 kPa的比值为0.5。超静孔压的响应表明,现场液化试验中,设备施加给饱和砂的上覆压力能够保持大致稳定。这也说明两次现场液化试验均是成功的,可以结合土压力、地表沉降等响应做进一步分析。

    图  8  现场液化试验的超静孔压时程(中心测点)
    Figure  8.  Time histories of excess pore pressure in in-situ liquefaction tests(at central measuring point)

    另外,从理论角度,上覆压力是影响液化响应的重要因素。相比Test-L试验,在相同的加速度幅值(0.11g左右)和相对密度(50%左右)下,Test-H试验中的试验孔压比相对较小,应该还与饱和砂本身的上覆压力过大(50 kPa)有较大关系。后文,结合上覆压力的响应,这点也被证实。为后续分析需要,图 9分别给出两次现场液化试验中饱和砂中心和左边测点的超静孔压时程。为图示清晰,仅保留残余孔压。由图 9可知,同一工况试验中,两测点的孔压时程几乎完全一致。

    图  9  现场液化试验饱和砂中心和左边测点超静孔压时程
    Figure  9.  Time histories of excess pore pressure at central and left measuring points in the in-situ liquefaction tests

    砂土中的上覆压力通过土压传感器实测得到。上覆压力,是本文试验非常关心的问题,直接关系试验技术的可行性,进行重点分析。

    图 10给出了两次现场试验的砂土中心和左边测点的土压力时程,也即上覆压力时程。低上覆压力工况Test-L试验中,因未能成功采集到加载期间的干砂的土压力时程数据,故无法给出。根据试验记录,加载前,Test-L试验干砂中心和左边测点实测上覆压力分别为23,20 kPa。从图 10可以看出,饱和砂中心和左边测点的上覆压力为25,21 kPa,两者基本一致。在高上覆压力工况Test-H试验中,干砂中心和左边测点的上覆压力分别为60,40 kPa,对应的饱和砂上覆压力约为61,42 kPa,中心和左边测点的上覆压力值差别较大。

    图  10  现场液化试验的上覆压力时程
    Figure  10.  Time histories of overburden pressure in in-situ liquefaction tests

    砂土的上覆压力,是由液压千斤顶通过加载钢板中心位置施加给砂土。笔者分析认为,中心和边侧压力值差别较大,主要原因应是加载板的刚度不够大导致了应力集中。只是荷载较小的情况下,加载板的刚度影响尚不明显。从理论角度,施加的上覆压力越大时,加载板的受力就越接近于点荷载模式。根据弹性力学的Boussinesq解答,砂土中应力分布的规律就应该是中间大,两侧小。定性上,实测结果与理论解答相符。由于土压力、孔压传感器体积很小,能够实现“点”测量,上覆压力和孔压均可实测得到,因此考虑上覆压力对液化响应的影响时,可以做到精确分析。基于简化考虑,本文施加上覆压力时,将中心和左边测点的上覆压力平均值视为最终的上覆压力。实际中,Test-L试验干砂、饱和砂的最终上覆压力分别为22,23 kPa,Test-H试验干砂、饱和砂的最终上覆压力分别为50,51 kPa,达到试验设计预期。

    为更好地分析上覆压力响应,略去实测数值的动态波动,并同时给出孔压比时程,如图 11所示,并将孔压急剧上升阶段前后的5~6 s时程详细给出。图 11给出的孔压比是每一时刻的残余孔压与对应的平均上覆压力的比值,这与通常液化分析中的孔压比的概念完全一致。

    图  11  现场液化试验的上覆压力与孔压比时程
    Figure  11.  Time histories of overburden pressure and excess pore pressure ratio in in-situ liquefaction tests

    图 11可以看出,地震动开始时(两次试验的时程上分别为20.5,13.0 s处),上覆压力均有短暂的小幅度增加,推测仍与加载钢板刚度不够大引起的弹性回弹有关。在动荷载作用下,干砂只会发生沉降,不存在土体软化问题。从其上覆压力响应可以看出,在整个加载过程中,Test-H试验中干砂的平均上覆压力缓慢下降,最终从50 kPa下降至40 kPa左右,而且中心和左边测点处的上覆压力下降趋势和程度基本一致。根据试验记录,Test-L试验中干砂的平均上覆压力,最终从22 kPa下降至14 kPa左右。干砂的上覆压力响应表明,施加上覆压力设备的工作性能是较为理想的。相比之下,饱和砂的上覆压力,随着超静孔压的急剧上升和土体软化,在短时间内出现了较为明显的下降,尤其是饱和砂的中心测点处。对于液化试验,重点关注孔压急剧上升阶段。当超静孔压发展到最大值时,Test-L试验中的平均上覆压力从23 kPa跌至18 kPa左右,孔压比最大值约为0.85;Test-H试验中的平均上覆压力从51 kPa跌至35 kPa左右,孔压比最大值约为0.68。由于Test-L和Test-H工况试验中,仅存在上覆压力上的差别,即使考虑加载期间上覆压力的下降,从图 11中可知,当孔压比发展到最大时,Test-H试验中的上覆压力最低值也大于Test-L试验中的初始上覆压力。因此,这充分说明,Test-H试验中孔压发展程度较低(即孔压比相对较小),是上覆压力较高导致。也就是说,Test-H试验中的试验孔压比虽然为0.5,若在更大幅值的加速度作用下,还会增大的。这也表明,其他条件相同时,上覆压力越大,饱和砂的抗液化强度越高,这与目前对液化的理论认识和实际震害调查结果一致。

    图 11还可以看出,同一工况试验中,饱和砂中心和左边测点的上覆压力响应不同,在初始压力值、下降程度两方面均有差异。尤其是,在孔压上升阶段,中心测点处上覆压力的下降程度明显大于左边测点,这应该与边侧的饱和砂受到非液化土的约束有关。为进一步分析上覆压力对液化响应的影响,将两测点处的孔压比分别计算,即为该测点处每一时刻的残余孔压与对应的上覆压力的比值,如图 12所示。为了图示清晰,仅给出孔压急剧上升阶段前后的5 s时程。由图 12可知,Test-L试验饱和砂中心和左边测点的最大孔压比分别为1和0.7,Test-H试验饱和砂中心和左边测点的最大孔压比分别为0.75和0.6。Test-L试验中,砂土中心测点处的孔压比为1,这说明此处的砂土已经完全液化。

    图  12  现场液化试验饱和砂中心和左边测点的孔压比时程
    Figure  12.  Time histories of excess pore pressure ratio at central and left measuring points in in-situ liquefaction tests

    值得注意的是,Test-L和Test-H工况试验中,两者上覆压力的差别本身并不大,仅有25 kPa,但是在地震动加速度幅值0.11g下,饱和砂的液化程度(孔压比)却有较大差异。而且,对于同一工况试验而言,饱和砂中心和两侧的液化发展程度也不一样。从图 12中也可看到,Test-L和Test-H试验中,中心和左边测点的孔压比时程曲线均存在交叉现象。由此可知,液化发展初期(边侧上覆压力较小),边侧的孔压比是大于中心处的,随着发展(边侧上覆压力变大),边侧的孔压比是小于中心处的。这说明,在此工况条件下,上覆压力对液化的影响较为明显。尤其值得注意的是,Test-L工况试验中,中心和左侧处的上覆压力差值仅有4 kPa(相当于40 cm厚的饱和土层),但是时程曲线上仍表现出交叉现象,定性上的响应规律与Test-H工况是一致的。这更表明,在上覆压力水平较低(至少50 kPa以下)时,上覆压力这一因素对液化响应的影响是敏感的。

    上覆压力50 kPa,大致对应真实场地中埋深2.5~5 m土层的埋藏条件。可以说,真实地震中,普遍且大量存在着上覆压力为50 kPa及以下的液化场地。但是,目前液化研究中,上覆压力(有效围压)50 kPa以下的液化试验成果相对很少,很多问题尚缺乏清晰认识。而且,当前一些实用的液化分析方法,考虑埋深影响的修正,也多是以上覆压力(有效围压)100 kPa作为基准参照。因此,低上覆压力水平(0~50 kPa)下的液化响应问题,应需重视并进行深入探究。

    本文现场液化试验的反力架是稳定可靠的固定端,因此能够方便、准确地获取整个加载期间的砂土地表沉降时程。

    图 13给出了两次现场试验的砂土地表沉降位移和孔压比时程,并将孔压急剧上升阶段前后的5~6 s时程详细给出。从图 13中可以看出,无论Test-L试验还是Test-H试验中,孔压迅速发展阶段,饱和砂发生的地表沉降了约为3.5~4.5 mm,将该沉降值换算成体应变约为1%(砂土高度300 mm)。孔压刚开始发展(触发)时的沉降太小,实测的位移难以准确甄别,但是可以确定当地表沉降在0.3~0.5 mm时,即体应变约为0.1%时,孔压比发展到约0.1左右。由此可知,当体应变在0.1%以下时,孔压就已经开始发展。当体应变在0.1%~1%时,孔压快速地大幅度发展。本文得到的孔压与应变关系的定量结果,与笔者之前的现场液化结果基本一致[6],也与目前国际上普遍认同的Dobry等人有关应变与孔压关系的定量认识基本一致[13]。需说明的是,Dobry等试验结果的剪应变γ,是根据土样轴向应变ε,通过γ=(1+ν)ε的理论关系换算得到,其泊松比ν取为0.5。由于前面的系数仅为1.5,因此剪应变和侧限条件下体应变的量级不会发生变化。因此,定量认识上两者是一致的。这表明,本文试验中获取的地表沉降响应是可靠的,能够反映液化的主要变形特征,并可做进一步分析。

    图 13中还可以发现,地震动开始阶段,同等工况下的干砂、饱和砂的沉降量和沉降速率几乎完全一致。但是,当饱和砂的孔压比发展到0.6左右时,两者的沉降响应开始出现明显差异。此后一段时间,饱和砂的沉降量和沉降速率要明显小于同等工况下的干砂,这可能与饱和砂的超静孔隙水无法及时排出有较大关系,而干砂内的气体却能随时排出,从而表现出这种沉降差异现象。随着孔压的逐渐消散,饱和砂的沉降量和沉降速率明显增大,并最终大于干砂。由此可知,饱和砂液化后的沉降主要发生在孔压消散阶段。这一沉降规律,与他人的动力离心机液化试验一致,也与真实地震现场观察到的现象一致[14-15]。这也说明,本文现场液化试验,能够较理想地再现真实液化场地的地表沉降响应特征。

    图  13  现场液化试验的地表沉降与孔压比时程
    Figure  13.  Time histories of soil surface settlement and excess pore pressure ratio in in-situ liquefaction tests

    通过对宏观试验现象以及超静孔压、土压力(上覆压力)、地表沉降等响应的分析,并与已有的液化认识对比,本文设计完成的两个不同上覆压力工况的现场液化试验,均可得到液化的主要响应特征和规律,并能够以此获取新认识。这说明,本文提出的考虑上覆压力影响的现场液化试验技术,总体上是可行和可靠的,能够为在1g重力环境下,以更接近真实场地的模拟条件开展液化相关问题研究,提供一种可能的新手段。

    本文在已有现场液化试验工作基础上,围绕“多工况、可重复、易操作、低成本”的基本设想,提出了一套考虑上覆压力影响的现场液化试验技术,设计完成了2次不同上覆压力(25,50 kPa)的饱和砂、干砂对比现场试验,分析了其动力响应规律及特征,并与已有认识对比,验证了试验技术的可靠性。主要得到以下3点结论。

    (1)新提出的试验孔压比这一概念,可作为评价考虑上覆压力影响的现场液化试验是否成功的基本指标。两个工况的现场试验表明,新提出的考虑上覆压力影响的现场液化试验技术可行,能够得到土壤液化的主要响应特征。这为1g重力环境下,以更接近真实场地的模拟条件开展液化相关问题研究,提供了一种可能的新手段。

    (2)在低上覆压力水平(0~50 kPa)下,大致对应真实场地0~5 m的埋深,上覆压力是影响土壤液化响应的敏感因素。实际地震中,普遍且大量存在上覆压力50 kPa及以下条件的液化场地,而目前与该静应力水平一致的液化研究成果并不丰富,相关认识也不清晰。低上覆压力水平(0~50 kPa)下的土壤液化响应问题,应需重视并进行深入探究。

    (3)当饱和砂的体应变在0.1%以下时,孔压已经开始发展;当饱和砂体应变在0.1%~1%时,孔压迅速且大幅度发展。当孔压比发展到0.6左右时,同等工况下的干砂、饱和砂的地表沉降响应出现差异。并且,液化砂土的沉降主要发生在孔压消散阶段。

  • 图  1   不同打磨时长的玻璃珠在光学显微镜与扫描电镜下图像

    Figure  1.   Images of glass beads subjected to different durations of abrasion obtained by optical and scanning electron microscopes

    图  2   3种玻璃珠表面纹理

    Figure  2.   Surface textures of three types of glass beads

    图  3   3种玻璃珠的Sq测试结果

    Figure  3.   Measured Sq for three types of glass beads

    图  4   能量注入式虚拟质量共振柱系统示意图

    Figure  4.   Schematic diagram of energy injecting virtual mass resonant column system

    图  5   孔隙比与有效应力对Gmax影响

    Figure  5.   Effects of void ratio and effective stress on Gmax

    图  6   不同有效应力下粗糙度对Gmax影响

    Figure  6.   Effects of roughness on Gmax under different effective confining pressures

    图  7   有效应力对3种表面粗糙度不同的颗粒材料Gmax的影响

    Figure  7.   Effects of effective stress on Gmax of granular materials with different particle surface roughnesses

    图  8   Sq对拟合参数An的影响

    Figure  8.   Effects of Sq on fitting parameters A and n

    图  9   孔隙比对Dmin影响

    Figure  9.   Effects of void ratio on Dmin

    图  10   有效应力对Dmin影响

    Figure  10.   Effect of effective stress on Dmin

    图  11   表面粗糙度对颗粒材料Dmin的影响

    Figure  11.   Effects of surface roughness on Dmin of granular materials

    图  12   未打磨的不同粒径玻璃珠试样的Gmax

    Figure  12.   Gmax of original glass beads with different sizes

    图  13   Gmax/f(e)随Sq/R变化规律

    Figure  13.   Variation of Gmax/f(e) with Sq/R

    图  14   天然砂粒与玻璃珠Gmax

    Figure  14.   Gmax of natural sand particles and glass beads

    表  1   试验方案

    Table  1   Test schemes

    编号 初始孔隙比e0 有效应力/kPa
    GB-2-A 0.620, 0.622, 0.637, 0.652 50, 100, 200, 400
    GB-2-B 0.628, 0.631, 0.632, 0.633, 0.658, 0.662 50, 100, 200, 400
    GB-2-C 0.637, 0.654, 0.660, 0.669 50, 100, 200, 400
    下载: 导出CSV

    表  2   Gmax拟合参数

    Table  2   Fitting parameters of Gmax

    颗粒类别 d A n R2
    GB-2-A 609.28 0.370 0.971
    GB-2-B 1.935 557.38 0.380 0.978
    GB-2-C 420.88 0.486 0.995
    下载: 导出CSV

    表  3   Dmin拟合参数表

    Table  3   Fitting parameters of Dmin

    颗粒类别 D1 κ R2
    GB-2-A 0.0062 0.451 0.846
    GB-2-B 0.0069 0.580 0.910
    GB-2-C 0.0060 0.281 0.867
    下载: 导出CSV

    表  4   颗粒材料相关信息[7, 16-19]

    Table  4   Information of granular materials [7, 16-19]

    颗粒名称 Sq/nm OR Cu d50/mm
    Leighton Buzzard砂 300 0.87 1.45 0.84
    Hostun砂 1972 0.84 1.48 0.34
    丰浦砂 1847 0.85 1.39 0.216
    2 mm玻璃珠 64, 602, 1176 0.960 1.02 2
    1 mm玻璃珠 170 0.955 1.06 1
    下载: 导出CSV
  • [1] 杨文保, 吴琪, 陈国兴. 长江入海口原状土动剪切模量预测方法探究[J]. 岩土力学, 2019, 40(10): 3889-3896. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910023.htm

    YANG Wenbao, WU Qi, CHEN Guoxing. Dynamic shear modulus prediction method of undisturbed soil in the estuary of the Yangtze River[J]. Rock and Soil Mechanics, 2019, 40(10): 3889-3896. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910023.htm

    [2]

    WICHTMANN T, TRIANTAFYLLIDIS T. Influence of the grain-size distribution curve of quartz sand on the small strain shear modulus Gmax[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2009, 135(10): 1404-1418. doi: 10.1061/(ASCE)GT.1943-5606.0000096

    [3]

    LIU X, YANG J. Shear wave velocity in sand: effect of grain shape[J]. Géotechnique, 2018, 68(8): 742-748. doi: 10.1680/jgeot.17.T.011

    [4]

    SARKAR D, GOUDARZY M, WICHTMANN T. Influence of particle size on the small-strain stiffness in granular soils: experimental observations and micromechanical interpretation[J]. Géotechnique, 2022, 19(6): 545-554.

    [5] 黄志鹏, 魏骁, 潘坤, 等. 颗粒特性对材料小应变动力特性的影响[J]. 岩土工程学报, 2022, 44(5): 889-897. doi: 10.11779/CJGE202205012

    HUANG Zhipeng, WEI Xiao, PAN Kun, et al. Effects of particle characteristics on small-strain dynamic properties of granular materials[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(5): 889-897. (in Chinese) doi: 10.11779/CJGE202205012

    [6]

    MITCHELL J K, SOGA K. Fundamentals of Soil Behavior[M]. New York: John Wiley & Sons, 2005.

    [7]

    CAVARRETTA I. The Influence of Particle Characteristics on the Engineering Behaviour of Granular Materials[D]. London : Imperial College London, 2009.

    [8]

    SANTAMARINA C, CASCANTE G. Effect of surface roughness on wave propagation parameters[J]. Géotechnique, 1998, 48(1): 129-136. doi: 10.1680/geot.1998.48.1.129

    [9]

    SHARIFIPOUR M, DANO C. Effect of grains roughness on waves velocities in granular packings[C]//1st Euro- Mediterranean Symposium in Advances on Geomaterials and Structures. Hammamet, 2006.

    [10]

    OTSUBO M, O'SULLIVAN C. Experimental and DEM assessment of the stress-dependency of surface roughness effects on shear modulus[J]. Soils and Foundations, 2018, 58(3): 602-614. doi: 10.1016/j.sandf.2018.02.020

    [11]

    YANG Z X, WEN Y X, PAN K. Previbration signature on dynamic properties of dry sand[J]. Journal of Testing and Evaluation, 2019, 47(3): 2167-2192. doi: 10.1520/JTE20180037

    [12] 蔡正银, 李相崧. 材料状态对干砂小应变特性的影响[J]. 岩土力学, 2004, 25(1): 10-14. doi: 10.3969/j.issn.1000-7598.2004.01.003

    CAI Zhengyin, LI Xiangsong. Effects of material state on the small strain behavior of dry sand[J]. Rock and Soil Mechanics, 2004, 25(1): 10-14. (in Chinese) doi: 10.3969/j.issn.1000-7598.2004.01.003

    [13]

    MENQ F Y, STOKOE K H. Linear dynamic properties of sandy and gravelly soils from large-scale resonant tests[C]// Deformation Characteristics of Geomaterials. IS Lyon, 2003.

    [14]

    SHARIFIPOUR M, DANO C, HICHER P. Wave velocities in assemblies of glass beads using bender-extender elements[C]// 17th ASCE Engineering Mechanics Conference. Newark, 2004.

    [15]

    PATEL A, BARTAKE P P, SINGH D N. An empirical relationship for determining shear wave velocity in granular materials accounting for grain morphology[J]. Geotechnical Testing Journal, 2008, 32(1): 1-10.

    [16]

    ALSHIBLI, KHALID A, et al. Quantifying morphology of sands using 3D imaging. Journal of Materials in Civil Engineering[J]. Journal of Materials in Civil Engineering, 2015, 27(10): 04014275. doi: 10.1061/(ASCE)MT.1943-5533.0001246

    [17]

    ALTUHAFI F, O'SULLIVAN C, CAVARRETTA I. Analysis of an image-based method to quantify the size and shape of sand particles[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2013, 139(8): 1290-1307. doi: 10.1061/(ASCE)GT.1943-5606.0000855

    [18]

    GU X Q, YANG J, HUANG M S. Laboratory measurements of small strain properties of dry sands by bender element[J]. Soils and Foundations, 2013, 53(5): 735-745. doi: 10.1016/j.sandf.2013.08.011

    [19]

    GOUDARZY M, KÖNIG D, SCHANZ T. Interpretation of small and intermediate strain characteristics of Hostun sand for various stress states[J]. Soils and Foundations, 2018, 58(6): 1526-1537. doi: 10.1016/j.sandf.2018.09.002

  • 期刊类型引用(1)

    1. 梁越,李越,许彬,杨牛虎,张斌,潘剑,尹朝恒. 考虑围压影响的间断级配散粒土潜蚀规律试验研究. 河海大学学报(自然科学版). 2025(01): 80-86 . 百度学术

    其他类型引用(1)

图(14)  /  表(4)
计量
  • 文章访问数:  483
  • HTML全文浏览量:  50
  • PDF下载量:  107
  • 被引次数: 2
出版历程
  • 收稿日期:  2023-05-09
  • 网络出版日期:  2024-03-24
  • 刊出日期:  2024-07-31

目录

/

返回文章
返回