基于剪切应变特征的饱和珊瑚砂超静孔压发展模型试验研究

吴琪; 王路阳; 刘启菲; 周正龙; 马维嘉; 陈国兴

doi:10.11779/CJGE20220956

基于剪切应变特征的饱和珊瑚砂超静孔压发展模型试验研究 English Version

吴琪^{1, 2,},
王路阳¹,
刘启菲¹,
周正龙¹,
马维嘉^3, ,,
陈国兴¹

1.
南京工业大学岩土工程研究所, 江苏南京 210009
2.
北京工业大学城市建设部, 北京 100124
3.
南京理工大学机械工程学院, 江苏南京 210094

基金项目:

国家自然科学基金项目 52008206

国家自然科学基金项目 52278503

国家自然科学基金项目 52208349

详细信息

作者简介:
吴琪(1991—)，男，博士，硕士生导师，副教授，主要从事珊瑚砂动力特性试验研究。E-mail：qw09061801@163.com

通讯作者:
马维嘉，E-mail: mwj9107@njust.edu.cn

中图分类号: TU432
计量
- 文章访问数: 279
- HTML全文浏览量: 46
- PDF下载量: 76
出版历程
- 收稿日期: 2022-08-03
- 网络出版日期: 2023-10-16
- 刊出日期: 2023-09-30

Experimental study on development model of excess pore pressure for saturated coral sand based on shear strain characteristics

1.
Institute of Geotechnical Engineering, Nanjing Tech University, Nanjing 210009, China
2.
Faculty of Architecture, Civil and Transportation Engineering, Beijing University of Technology, Beijing 100124, China
3.
School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China

摘要

摘要: 利用空心圆柱扭剪仪，针对南沙岛礁饱和珊瑚砂试样进行了一系列不排水循环扭剪试验，以探究细粒含量FC、相对密度D_r及循环应力比CSR对饱和珊瑚砂超静孔压比R_u增长特性的影响。试验结果表明：R_u-N曲线发展速率随FC及CSR的增大而增大，随D_r的增大而减小；饱和珊瑚砂的R_u-N曲线发展速率明显超出了陆域石英砂，传统的Seed孔压模型对珊瑚砂孔压发展规律的表征不完全适用。相同的FC和D_r下，不同CSR作用下的珊瑚砂R_u与剪应变幅值γ_a具有唯一性关系。建立了基于剪应变特征的孔压发展预测模型，分析表明，孔压模型参数A为土性相关的常数，密度修正的孔压模型参数的B/(D_r)^1.5与等效骨架孔隙比e_sk^*具有单一的负幂函数关系。
- 珊瑚砂 /
- 细粒含量 /
- 超静孔压 /
- 剪切应变 /
- 等效骨架孔隙比
Abstract: A series of undrained cyclic torsional shear tests are carried out on the saturated coral sand specimens from Nansha Islands by using the hollow cylinder torsion shear apparatus. The aim of the tests is to investigate the effects of fine contents(FC), relative density (D_r) and applied cyclic stress ratio CSR on the generated excess pore pressure ratio R_u of the specimens. The test results show that: the development rate of R_u-N curve increases with the increasing FC and CSR, but decreases with the increasing D_r. In addition, the development rate of R_u-N curve of the saturated coral sands significantly exceeds that of the terrestrial quartz sands, and the traditional Seed pore pressure model is not fully applicable to characterize the development pattern of the excess pore pressure of coral sands. For the same D_r and FC, R_u of the specimens under different CSR is uniquely related to the shear strain amplitude (γ_a). A pore pressure evaluation model based on the shear strain characteristics is established, and the analysis shows that the pore pressure model parameter A is a soil-specific constant, and the density-corrected pore pressure model parameter B/(D_r)^1.5 has a single negative power function relationship with the equivalent skeleton void ratio e_sk^*.
- coral sand /
- fines content /
- excess pore pressure /
- shear strain /
- equivalent skeleton void ratio

HTML全文

0. 引言

钙质砂广泛存在于世界低纬度地区的大陆架和沿海地区，是由贝壳和珊瑚等海洋生物的碎片或残骸经破碎沉积产生的，其主要成分为碳酸钙^[1]。钙质砂在形成过程中保留了原生生物骨骼的微观结构，其一般具有形状不规则、内孔隙发育、单颗粒强度低和易破碎的特点^[2]。钙质砂与结构物之间的接触在海洋岩土工程建设中较常见，如海洋风机和海洋钻井平台的桩基础^[3-4]。钙质砂作为一种力学性质不良的建筑材料，其与结构物之间的相互作用对构筑物的整体稳定性和安全运营十分重要。研究表明，土与结构物之间传递荷载时，结构物附近的薄层土体处最容易发生应变局部化，产生开裂、滑移和破坏，因此土-结构物界面的力学特性值得关注^[5-6]。

针对砂-结构物界面的剪切行为，目前已用学者对诸如砂的相对密实度^[7-9]、颗粒形状^[10-11]、颗粒尺寸^[12-13]、级配^[14]等影响因素开展研究。实际工程中，由于地震、交通和波浪荷载等对结构物的作用，砂-结构物界面常遭受循环荷载，通常需要计算界面的循环剪切刚度和阻尼比来反映其动力响应特征。已有研究表明，法向应力、剪切位移幅值、结构物材质、颗粒形状和粒径等都会对界面循环剪切刚度和阻尼比产生影响^[15-16]。

海洋平台的抗拔桩在贯入阶段需要尽可能的减小桩侧摩阻力，而在运营阶段则需要桩基提供更大的侧摩阻力以抵抗上拔力，因此需要一种能够根据使用工况调整侧摩阻力的桩基形式。近年来，蛇在运动时腹部鳞片与接触材质之间的摩擦各向异性（frictional anisotropy）逐渐受到人们的重视。Marvi等^[17]的研究表明蛇在运动过程中可通过肌肉调整腹部鳞片的张开角度来调整摩擦以保证不同的运动状态。此后，Martinez等^[18]根据仿生学原理，将多种蛇的腹部鳞片结构抽象化为蛇皮仿生界面，并开展了界面直剪试验，结果表明蛇皮仿生界面具有明显的摩擦各向异性。之后，Stutz等^[19]的研究进一步表明砂-蛇皮仿生界面的强度和剪胀性与仿生界面的表观形态有关，并认为鳞片的几何特征比能合理反应蛇皮仿生界面的力学特征。此外，O'Hara等^[20-21]对蛇皮仿生模型桩并开展的离心机试验表明桩侧摩阻力以及桩周土体扰动均与仿生桩的桩身鳞片排布方向有关。Zhong等^[22]通过离散元数值模拟方法同样证明了蛇皮仿生桩在竖向荷载作用下表现出了摩擦各向异性。上述研究的试验材料主要为石英砂，而钙质砂的力学特性显然与石英砂不同。Xiao等^[23]针对钙质砂-蛇皮仿生界面开展了静荷载直剪试验，表明除摩擦各向异性外，钙质砂-蛇皮仿生界面同样存在颗粒破碎各向异性（particle breakage anisotropy）。

综上所述，已有研究对钙质砂-蛇皮仿生界面的力学行为的研究有限，尤其是钙质砂-蛇皮仿生界面循环荷载作用下的动力响应问题一直未有学者进行系统性研究。因此，本文针对钙质砂-蛇皮仿生界面开展了一系列循环直剪试验，研究了法向应力、剪切位移幅值和蛇皮仿生界面表观形态对剪切刚度和阻尼比的影响，为蛇皮仿生界面在海洋平台抗拔桩的应用提供一定理论参考。

1. 试验介绍

1.1 试验设备

本试验采用南京泰克奥科技公司生产的液晶微控直剪仪TKA-DDS-10，界面剪切试验系统如图 1所示。

图 1 界面剪切试验设备

Figure 1. Test equipment for interface shear

下载: 全尺寸图片幻灯片

该直剪仪可通过电脑软件实现单调和循环加载并精确控制剪切速率和剪切位移，剪切速率变化范围为0.001~2.4 mm/min。法向荷载通过加载砝码和杠杆施加。法向位移通过法向位移数显表显示并直接传输给电脑软件记录，剪应力由水平荷载传感器测量，分辨率为1 N。上剪切盒内部尺寸为直径61.8 mm，高度为22 mm，如图 1（c）所示，下剪切盒尺寸为100 mm×100 mm×12.5 mm，如图 1（d）所示。采用下界面配置，即蛇皮仿生界面安装在下剪切盒中，钙质砂试样安装在上剪切盒中，剪切时上剪切盒保持不动，下剪切盒带动蛇皮仿生界面在导轨上水平移动。

1.2 钙质砂与蛇皮仿生界面

本研究所用钙质砂取自中国南海某岛礁，相对质量密度为2.78，通过X射线衍射分析（XRD）可知，该钙质砂的主要成分为碳酸钙，含量约为95.9%，其余杂质主要为碳酸镁，含量约为4%。试验前将钙质砂清洗以去除表面盐分并在105°烘箱中烘干24 h。通过标准筛对钙质砂进行筛分，选取粒径范围为0.4~0.6 mm的砂样开展试验，砂样的最大和最小孔隙比分别为1.196，0.817。图 2（a）为100倍率下扫描电子显微镜拍得的钙质砂照片，可以发现，钙质砂表面内孔隙发育，同时砂颗粒棱角较多，形状不规则。

图 2 试验材料与循环剪切加载路径

Figure 2. Test materials and cyclic shear loading path

下载: 全尺寸图片幻灯片

Martinez等^[18]的研究表明蛇皮仿生界面的鳞片结构为直角三角形时，界面表现出更强的摩擦各向异性，同时认为蛇腹部鳞片几何形状的非对称性带来的摩擦各向异性要比鳞片的柔性和鳞片张开角度变化更具有工程研究意义。基于上述原因，本文试验采用的蛇皮仿生界面的鳞片结构为直角三角形，如所示。同时，考虑到海洋工程中的桩基础一般采用钢材质，因此本试验中的蛇皮仿生界面采用不锈钢加工而成，其弹性模量为194 GPa，极限刚拉强度不小于515 MPa，洛氏硬度不大于92。蛇皮仿生钢板的尺寸为99.5 mm×99.5 mm×12.5 mm，其尺寸略小于下剪切盒内尺寸，以保证钢板可顺利安放在下剪切盒中。为保证循环剪切时蛇皮仿生钢板不在下剪切盒中移动，钢板安装完后用塑料薄片填充钢板与下剪切盒内壁之间的空隙。蛇皮仿生钢板的尺寸大于上剪切盒中砂样尺寸，可保证剪切过程中砂样与钢板之间的接触面积保持不变。从可以看出，蛇皮仿生界面的鳞片结构主要有鳞片高度 $H$ 和鳞片长度 $L$ 两个几何参数。为了探究蛇皮仿生界面表观几何形态对界面循环剪切动力响应的影响，本试验选择了6种蛇皮仿生界面，鳞片的高度分别为0.5，1.0 mm，对应的鳞片长度分别为2.5，5.0，10，5.0，10，20 mm。两种鳞片高度的蛇皮仿生界面的鳞片几何特征比 $L/H$ 分别为5，10，20。

1.3 试验方案与试验流程

本文研究的钙质砂-蛇皮仿生界面循环剪切试验采用带活塞的漏斗对干燥状态的钙质砂自然落砂装样，试样高度为20 mm，试样相对密实度为90%。装样时将所需的钙质砂按质量分为3等份，采用分层装样法以保证试样的均匀性。每层装样完毕后将试样表面轻轻抚平，垫上圆形木块，之后用橡胶锤轻轻敲击木垫块，将试样压实至设置的装样高度，尽量减小压实试样过程中产生的颗粒破碎。此外，用细铁丝在上一层的砂样表面划线，保证两层试样之间无明显分界。

采用砝码进行竖向加载，待固结稳定后进行界面循环剪切试验，所施加的加载路径如图 2（c）所示。第一个循环自平衡位置开始，到达设定的剪切位移幅值，之后反向剪切，如图 2（c）中的①→②→③→④，定义①和④分别为正方向的加载和卸载阶段，②和③分别为反方向的卸载和加载阶段。此外，由于蛇皮仿生界面上鳞片几何结构的非对称性，界面剪切试验分别为两种加载路径：正向剪切为尾向剪切（caudal shearing），反向剪切为颅向剪切（cranial shearing）；正向剪切为颅向剪切，反向剪切为尾向剪切。其中，尾向剪切的定义为砂试样顺着鳞片结构的斜边移动，模拟了蛇在向前运动时腹部鳞片闭合的情况；颅向剪切的定义为砂试样逆着鳞片结构的直角边移动，模拟了蛇在攀爬状态下腹部鳞片张开的情况^[18-19]。

钙质砂-蛇皮仿生界面循环剪切过程中剪切速率保持不变，设定为0.8 mm/min，为循环剪切加载的波形，其中 $T$ 为一个循环周期， ${A_{\text{w}}}$ 为剪切位移幅值。表 1为本文开展的所有循环剪切试验，共计24组。其中，试验组T1~T8主要分析了法向应力的影响，T9~T14主要分析了剪切位移幅值的影响，T15~T24主要分析了蛇皮仿生界面表观形态的影响。

表 1 界面循环剪切试验方案

Table 1. Cyclic test programs of interface

编号	表观形态	法向应力/kPa	位移幅值/mm	循环次数/次	加载路径
T1	H0.5L5.0	100	1.6	15	尾向→颅向
T2	H0.5L5.0	200	1.6	15
T3	H0.5L5.0	400	1.6	15
T4	H0.5L5.0	800	1.6	15
T5	H0.5L5.0	100	1.6	15	颅向→尾向
T6	H0.5L5.0	200	1.6	15
T7	H0.5L5.0	400	1.6	15
T8	H0.5L5.0	800	1.6	15
T9	H0.5L5.0	200	0.8	15	尾向→颅向
T10	H0.5L5.0	200	2.4	15
T11	H0.5L5.0	200	3.2	15
T12	H0.5L5.0	200	0.8	15	颅向→尾向
T13	H0.5L5.0	200	2.4	15
T14	H0.5L5.0	200	3.2	15
T15	H0.5L2.5	200	1.6	15	尾向→颅向
T16	H0.5L10	200	1.6	15
T17	H1.0L5.0	200	1.6	15
T18	H1.0L10	200	1.6	15
T19	H1.0L20	200	1.6	15
T20	H0.5L2.5	200	1.6	15	颅向→尾向
T21	H0.5L10	200	1.6	15
T22	H1.0L5.0	200	1.6	15
T23	H1.0L10	200	1.6	15
T24	H1.0L20	200	1.6	15
注：H0.5L5.0表示蛇皮仿生界面的表观形态为鳞片高度H = 0.5 mm，L = 5.0 mm，其他以此类推；尾向→颅向表示正向剪切为尾向剪切，反向剪切为颅向剪切的加载路径，颅向→尾向表示正向剪切为颅向剪切，反向剪切为尾向剪切的加载路径。

下载: 导出CSV

| 显示表格

2. 剪切刚度与阻尼比的定义

土体剪切刚度和阻尼比常被用来分析循环荷载下土体的动力响应^[16]。剪切刚度是指剪应力与所对应的剪切位移的比值，反映了土体抵抗剪切变形的能力。阻尼比是表征土体循环剪切过程中能量耗散快慢的参量，通常阻尼比越大循环剪切过程中能量耗散越快。在计算剪切刚度时本文采用Vieira等^[24]提出的割线刚度定义：

$K{\text{ = }}\frac{{{\tau _{\max }} - {\tau _{\min }}}}{{{u_{t, \max }} - {u_{t, \min }}}} 。$

(1)

式中： ${u_{t, \max }}$ 和 ${u_{t, \min }}$ 分别为循环剪应力滞回圈中最大和最小剪切位移，其绝对值一般为剪切位移幅值； ${\tau _{\max }}$ 和 ${\tau _{\min }}$ 分别为应力滞回圈中的最大和最小剪切位移幅值对应的剪应力。

按照式（1）计算得到的剪切刚度即为应力滞回圈中正、负剪应力峰值连线的斜率，如图 3所示。阻尼比的计算方法采用Nye等^[25]、刘飞禹等^{[15, 26]}和芮圣洁等^[16]提到的计算公式：

$D = \frac{1}{2}\left( {\frac{S}{{4{\text{π}}{S_1}}} + \frac{S}{{4{\text{π}}{S_2}}}} \right) 。$

(2)

图 3 剪切刚度和阻尼比示意图

Figure 3. Schematic diagram of shear stiffness and damping ratio

下载: 全尺寸图片幻灯片

式中： $S$ 为应力滞回圈的面积； ${S_1}$ 和 ${S_2}$ 分别为图 3中左右两个三角形阴影区域的面积。

3. 试验结果分析

3.1 钙质砂-蛇皮仿生界面典型的强度变化

以法向应力 ${\sigma _n}$ =200 kPa，蛇皮仿生界面表观形态为H0.5L5.0，剪切位移幅值 ${A_w}$ =1.6 mm，循环次数为15次的工况为例，本文试验在两种加载路径下界面循环剪切的强度变化如图 4所示。

图 4 两种加载路径下界面循环剪切强度

Figure 4. Change of cyclic shear strength under two loading paths

下载: 全尺寸图片幻灯片

从图 4（a），（b）可以看出，两种加载路径下钙质砂-蛇皮仿生界面的剪应力-剪切位移关系形成的滞回圈都呈现出明显的非对称性。尾向→颅向加载路径下正向剪切应力小于反向剪应力，而颅向→尾向加载路径下正向剪应力大于反向剪切应力。根据Xiao等^[23]的研究表明，钙质砂-蛇皮仿生界面尾向剪切时的强度小于颅向剪切，因此循环剪切时应力滞回圈表现出非对称性。此外，可以发现两种加载路径下的应力滞回圈形状基本互为反对称。

从图 4（c），（d）可以看出，两种循环加载路径下正向剪切时钙质砂发挥的界面摩擦角都呈现先减小后增大，在10次循环后逐渐稳定的趋势；反向剪切时第1个循环内发挥的界面摩擦角则都小于后续的其它循环，界面摩擦角呈现先随循环次数增加而增大，之后逐渐稳定的趋势。主要是由于首次循环正向剪切时钙质砂颗粒与仿生界面鳞片结构之间的咬合作用较大，在随后的剪切过程中咬合作用由于砂颗粒的运动有所减弱，随着循环次数的增加，试样逐渐被循环剪密，砂颗粒之间及砂颗粒与界面之间的咬合作用逐渐增加并趋于稳定。

3.2 法向应力的影响

研究表明，法向应力不仅影响界面的剪切强度，同时会影响界面的剪切刚度和阻尼比等动力响应特征^{[15-16, 26]}。图 5给出了剪切位移为1.6 mm，蛇皮仿生界面表观形态为H0.5L5.0，钙质砂-蛇皮仿生界面在两种循环加载路径下和不同法向应力下的剪切刚度和阻尼比。

图 5 不同法向应力下两种加载路径的剪切刚度和阻尼比

Figure 5. Shear stiffnesses and damping ratios of two loading paths under different normal stresses

下载: 全尺寸图片幻灯片

由图 5（a）可知，钙质砂-蛇皮仿生界面在两种加载路径和不同法向应力水平下的剪切刚度都随循环次数增大而增加，表明界面剪切过程中发生了剪切硬化现象，且前几次循环中剪切刚度增长较快，随着循环次数增加剪切刚度增长速率逐渐减小，剪切刚度趋于稳定，这与图 4（c），（d）中对界面摩擦角的分析结果相符合。此外，循环加载路径为尾向→颅向剪切时的剪切刚度与颅向→尾向剪切接近，说明加载路径对剪切刚度的影响较小。两种加载路径下界面剪切刚度都随法向应力的增加而增加，这主要是由于界面剪切本身是一种摩擦行为，界面剪切强度符合莫尔-库仑定理，因此法向应力的增加会导致界面剪切刚度的提升^[15-16]。

钙质砂-蛇皮仿生界面在两种加载路径和不同法向应力下的阻尼比随循环次数的变化规律如图 5（b）所示。从图 5（b）中可以看到，各工况下的界面阻尼比随循环次数的增加逐渐减小，其变化趋势与剪切刚度相反，表明随着循环次数增加界面剪切时的能量耗散速率逐渐减小，主要是由于钙质砂逐渐被循环剪密，砂颗粒的位置趋于稳定，循环剪切时输入的能量由于砂颗粒间的咬合作用增加而被存储。此外，由图 5（b）可知，法向应力越大，界面阻尼比越小，说明低法向应力下的界面循环剪切能量耗散更快。同时，可以发现较高法向应力下（400，800 kPa）两种循环加载路径下的阻尼比接近，较低法向应力下（100，200 kPa）两者之间有差异，但是随着循环次数的增加两种加载路径下的阻尼比逐渐趋同。

为了进一步探究钙质砂-蛇皮仿生界面循环剪切刚度和阻尼比与法向应力之间的关系，将剪切位移为1.6 mm，蛇皮仿生界面表观形态为H0.5L5.0时各法向应力下第1次、第5次、第10次和第15次循环时的剪切刚度和阻尼比绘制于。从，中可以看出，两种加载路径下4个循环次数时的剪切刚度和法向位移之间均呈现良好的线性关系，所有数据拟合的决定系数 ${R^2}$ 均接近于1.0，说明单次循环的剪切刚度随法向应力的增大等比例增大，这与芮圣洁等^[16]得到的结论一致。图 6（c），（d）的拟合结果显示4个循环次数时的阻尼比随法向应力的增加均呈对数下降。

图 6 剪切刚度和阻尼比与法向应力的关系

Figure 6. Relationship among shear stiffness, damping ratio and normal stress

下载: 全尺寸图片幻灯片

3.3 剪切位移幅值的影响

图 7为法向应力为200 kPa，蛇皮仿生界面表观形态为H0.5L5.0时两种循环加载路径下不同剪切位移幅值时的界面剪切刚度和阻尼比。从图 7（a）可以发现：无论是尾向→颅向循环剪切还是颅向→尾向循环剪切，剪切位移幅值越大，剪切刚度越小；剪切位移幅值较小时，剪切刚度在循环初期随循环次数增加而增加的趋势更加明显。图 7（b）表明，界面阻尼比同样受剪切位移幅值的影响，表现出随剪切位移幅值的增加而增加，表明剪切位移幅值的增大加快了循环剪切过程中的能量耗散，这主要是由于大剪切位移幅值下砂颗粒受剪切作用而发生的位置调整更加频繁，颗粒的翻滚更加剧烈，颗粒之间难以形成稳定的咬合作用，因此能量在颗粒的不断翻滚和位置调整中迅速耗散。此外，大剪切位移幅值下钙质砂颗粒破碎较多导致耗能量增加也可能是阻尼比增加的原因^[16]。

图 7 不同剪切位移幅值下两种加载路径的剪切刚度和阻尼比

Figure 7. Shear stiffnesses and damping ratios of two loading paths under different amplitudes of shear displacement

下载: 全尺寸图片幻灯片

图 8为法向应力为200 kPa，蛇皮仿生界面表观形态为H0.5L5.0时，各剪切位移幅值下循环次数为1，5，10，15时的界面剪切刚度和阻尼比关系。芮圣洁等^[16]比较了反比例函数和对数函数对剪切刚度和剪切位移幅值之间的拟合效果，发现反比例函数能更好地反映两者之间的关系。本文采用反比例函数对试验结果进行了最佳拟合，如图 8（a），8（b）所示，可以发现反比例函数同样适用于两种循环加载路径下的界面剪切刚度和剪切位移幅值之间的拟合关系。图 8（c），8（d）为阻尼比与剪切位移幅值之间的关系，可以发现两者之间的关系可以用对数函数实现良好的拟合。

图 8 剪切刚度和阻尼比与剪切位移幅值的关系

Figure 8. Relationship among shear stiffness, damping ratio and amplitude of shear displacement

下载: 全尺寸图片幻灯片

3.4 蛇皮仿生界面表观形态的影响

结构物表面的表观形态不同会导致不同的表面粗糙度，而界面剪切时表面粗糙度会对界面的力学特性产生重要影响，这主要是由于表面粗糙度的不同会改变界面的破坏模式：当粗糙度较低时，界面破坏模式为完全滑动破坏（a full sliding failure）；随着粗糙度的增加，由于土颗粒与粗糙面之间和咬合作用增加，界面破坏模式转变为土体破坏（soil failure）^[27-28]。砂-蛇皮仿生界面的静荷载试验结果表明，利用鳞片几何特征比 $L/H$ 表征蛇皮仿生界面的粗糙度可以合理地反映界面的力学行为，同时界面的强度和摩擦角随鳞片几何特征比的增加而减小，表明鳞片几何特征比的增加会降低蛇皮仿生界面的粗糙度^{[18-19, 23]}。图 9，10分别为法向应力为200 kPa，剪切位移为1.6 mm，两种鳞片高度下不同界面表观形态时的界面循环剪切刚度和阻尼比随循环次数的变化情况。从图中可以看出，在相同鳞片高度和相同加载路径下，剪切刚度随鳞片几何特征比的增加而降低，表明界面粗糙度减小可以降低界面的循环剪切强度。此外，阻尼比则随鳞片几何特征比的增加而增加，表明界面粗糙度较低时界面循环剪切时的能量耗散更快，主要是由于砂颗粒与低粗糙度的蛇皮仿生界面之间的咬合作用较弱，砂颗粒在剪切过程中更容易沿剪切方向平移运动，导致能量被迅速消耗。

图 9 H =0.5 mm时不同界面表观形态下两种加载路径的剪切刚度和阻尼比

Figure 9. Shear stiffnesses and damping ratios of two loading paths under different apparent shapes of interface at H=0.5 mm

下载: 全尺寸图片幻灯片

图 10 H =1.0 mm时不同界面表观形态下两种加载路径的剪切刚度和阻尼比

Figure 10. Shear stiffnesses and damping ratios of two loading paths under different apparent shapes of interface at H=1.0 mm

下载: 全尺寸图片幻灯片

为法向应力为200 kPa，剪切位移为1.6 mm，相同鳞片几何特征比（ $L/H{\text{ = }}5$ ）和不同鳞片高度（ $H$ 为0.5，1.0 mm）下两种循环加载路径的剪切刚度和阻尼比随循环次数的变化情况。从图 11中可以看出，相同鳞片几何特征比下界面剪切刚度随鳞片高度增加而增加，阻尼比随鳞片高度增加而降低，主要是由于蛇皮仿生界面上凸起鳞片的高度较大时可提供更大的粗糙度，钙质砂颗粒与蛇皮仿生界面之间的咬合作用更强，界面强度较大且循环剪切时的能量更多被存储下来。

图 11 相同

$L/H$ 不同鳞片高度下界面剪切刚度和阻尼比

Figure 11. Shear stiffnesses and damping ratios of interface with same

$L/H$ and different scale heights

下载: 全尺寸图片幻灯片

图 12，13分别为法向应力为200 kPa，剪切位移为1.6 mm，两种鳞片高度和两种循环加载路径下循环次数分别为1，5，10，15次时的界面剪切刚度-鳞片几何特征比关系和阻尼比-鳞片几何特征比关系。从图中可以发现，各剪切工况下的界面剪切刚度和鳞片几何特特征比之间呈幂函数下降趋势。除其中3组数据的拟合结果的决定系数小于0.9之外，其余的阻尼比和鳞片几何特征比之间的关系均可用幂函数实现良好的拟合。

图 12 剪切刚度与鳞片几何特征比的关系

Figure 12. Relationship between shear stiffness and scale geometric ratio

下载: 全尺寸图片幻灯片

图 13 阻尼比与鳞片几何特征比的关系

Figure 13. Relationship between damping ratio and scale geometric ratio

下载: 全尺寸图片幻灯片

4. 结论

本文基于界面直剪仪，对钙质砂-蛇皮仿生界面开展了一系列循环剪切试验，研究了法向应力水平、循环剪切位移幅值和蛇皮仿生界面表观形态对界面循环剪切动力响应的影响。得出4点结论。

（1）钙质砂-蛇皮仿生界面的应力滞回圈表现出非对称性，同时尾向→颅向加载路径和颅向→尾向加载路径下的应力滞回圈形状基本互为反对称。

（2）钙质砂-蛇皮仿生界面表现出循环硬化，剪切刚度随循环次数增加而增加，阻尼比随循环次数增加而降低。

（3）法向应力和剪切位移幅值对界面动力响应具有重要影响，表现为法向应力的增加会线性增加界面剪切刚度，同时阻尼比随法向应力增加呈对数下降趋势；剪切位移幅值的增加会导致剪切刚度呈反比例关系下降，而阻尼比则随剪切位移幅值的增加呈对数增加趋势。

（4）蛇皮仿生界面的表观形态对剪切刚度和阻尼比具有显著影响，表现为剪切刚度随鳞片几何特征比的增加呈幂函数增加，阻尼比则呈幂函数下降。同时，相同鳞片几何特征比下，鳞片高度越大，剪切刚度越大，阻尼比越小。

图 1 珊瑚砂电镜扫描图

Figure 1. SEM images of coral sand

下载: 全尺寸图片幻灯片

图 2 珊瑚砂级配曲线

Figure 2. Grain-size distribution curves of coral sand gradation curve

下载: 全尺寸图片幻灯片

图 3 D_r = 45%，CSR = 0.25条件下饱和珊瑚砂典型试验结果

Figure 3. Typical test results of saturated coral sand with D_r = 45%, CSR = 0.25

下载: 全尺寸图片幻灯片

图 4 饱和珊瑚砂R_u与N的关系曲线

Figure 4. Relationship between R_u and N for saturated coral sand

下载: 全尺寸图片幻灯片

图 5 饱和珊瑚砂超静孔压发展

Figure 5. Development patterns of excess pore water pressure for saturated coral sand

下载: 全尺寸图片幻灯片

图 6 剪应变幅值与超静孔压比的关系

Figure 6. Relationship between shear strain amplitude and excess pore pressure ratio

下载: 全尺寸图片幻灯片

图 7 参数A，B与FC关系曲线

Figure 7. Relationship between parameters A, B and FC

下载: 全尺寸图片幻灯片

图 8 密度修正的孔压模型参数B/(D_r)^1.5与e_sk^*关系曲线

Figure 8. Relationship between density-corrected pore pressure model parameter B/(D_r)^1.5 and e_sk^*

下载: 全尺寸图片幻灯片

表 1 纯砂粒、纯细粒及不同FC的珊瑚砂基本物理指标

Table 1 Basic physical indices of coral sand with pure sand, pure fine grain and different FC

物理指标	细粒含量FC/%
物理指标	0	6.41	10	20	30	100
e_max	1.79	1.72	1.70	1.65	1.62	1.52
e_min	1.12	0.99	0.91	0.77	0.69	0.89
d₅₀/mm	0.44	0.38	0.35	0.28	0.21	0.03
C_u	4.53	5.41	6.58	17.82	23.69	—
C_c	0.91	0.71	0.84	1.59	1.19	—
注：e_max为最大孔隙比；e_min为最小孔隙比；d₅₀为平均粒径；C_u为不均匀系数；C_c为曲率系数。

下载: 导出CSV

表 2 试验方案

Table 2 Test schemes

ID	D_r/%	FC/%	CSR	ID	D_r/%	FC/%	CSR	ID	D_r/%	FC/%	CSR	ID	D_r/%	FC/%	CSR	ID	D_r/%	FC/%	CSR
1	30	0	1.888	10	30	30	0.188	19	45	6.41	0.223	28	45	20	0.223	37	70	10	0.312
2	30	0	0.223	11	30	30	0.223	20	45	6.41	0.250	29	45	20	0.250	38	70	10	0.350
3	30	0	0.250	12	30	30	0.250	21	45	641	0.225	30	45	30	0.188	39	70	20	0.263
4	30	10	0.188	13	45	0	0.188	22	45	6.41	0.263	31	45	30	0.223	40	70	20	0.312
5	30	10	0.223	14	45	0	0.223	23	45	10	0.188	32	45	30	0.250	41	70	20	0.350
6	30	10	0.250	15	45	0	0.250	24	45	10	0.223	33	70	0	0.263	42	70	30	0.263
7	30	20	0.188	16	45	0	0.320	25	45	10	0.250	34	70	0	0.312	43	70	30	0.312
8	30	20	0.223	17	45	0	0.263	26	45	10	0.280	35	70	0	0.350	44	70	30	0.350
9	30	20	0.250	18	45	6.41	0.188	27	45	20	0.188	36	70	10	0.263

下载: 导出CSV

参考文献(30)

[1]	马维嘉, 陈国兴, 李磊, 等. 循环荷载下饱和南沙珊瑚砂的液化特性试验研究[J]. 岩土工程学报, 2019, 41(5): 981-988. doi: 10.11779/CJGE201905023 MA Weijia, CHEN Guoxing, LI Lei, et al. Experimental study on liquefaction characteristics of saturated coral sand in Nansha Islands under cyclic loading[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(5): 981-988. (in Chinese) doi: 10.11779/CJGE201905023
[2]	高冉, 叶剑红. 中国南海吹填岛礁钙质砂动力特性试验研究[J]. 岩土力学, 2019, 40(10): 3897-3908, 3919. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910024.htm GAO Ran, YE Jianhong. Experimental investigation on the dynamic characteristics of calcareous sand from the reclaimed coral reef islands in the South China Sea[J]. Rock and Soil Mechanics, 2019, 40(10): 3897-3908, 3919. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201910024.htm
[3]	王鸾, 汪云龙, 袁晓铭, 等. 人工场地吹填珊瑚土抗液化强度大粒径动三轴试验研究[J]. 岩土力学, 2021, 42(10): 2819-2829. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202110021.htm WANG Luan, WANG Yunlong, YUAN Xiaoming, et al. Experimental study on liquefaction resistance of hydraulic fill coralline soils at artificial sites based on large-scale dynamic triaxial apparatus[J]. Rock and Soil Mechanics, 2021, 42(10): 2819-2829. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202110021.htm
[4]	聂庆科, 白冰, 胡建敏, 等. 循环荷载作用下软土的孔压模式和强度特征[J]. 岩土力学, 2007, 28(增刊1): 724-729. https://cpfd.cnki.com.cn/Article/CPFDTOTAL-AGLU200710007148.htm NIE Qingke, BAI Bing, HU Jianmin, et al. The pore pressure model and undrained shear strength of soft clay under cyclic loading[J]. Rock and Soil Mechanics, 2007, 28(S1): 724-729. (in Chinese) https://cpfd.cnki.com.cn/Article/CPFDTOTAL-AGLU200710007148.htm
[5]	张建民. 砂土动力学若干基本理论探究[J]. 岩土工程学报, 2012, 34(1): 1-50. http://www.cgejournal.com/cn/article/id/14487 ZHANG Jianmin. New advances in basic theories of sand dynamics[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(1): 1-50. (in Chinese) http://www.cgejournal.com/cn/article/id/14487
[6]	孙海妹, 王兰民, 王平, 等. 饱和兰州黄土液化过程中孔压和应变发展的试验研究[J]. 岩土力学, 2010, 31(11): 3464-3468. doi: 10.3969/j.issn.1000-7598.2010.11.018 SUN Haimei, WANG Lanmin, WANG Ping, et al. Experimental study of development of strain and pore water pressure during liquefaction of saturated Lanzhou loess[J]. Rock and Soil Mechanics, 2010, 31(11): 3464-3468. (in Chinese) doi: 10.3969/j.issn.1000-7598.2010.11.018
[7]	徐斌, 孔宪京, 邹德高, 等. 饱和砂砾料振动孔压与轴向应变发展模式研究[J]. 岩土力学, 2006, 27(6): 925-928. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200606014.htm XU Bin, KONG Xianjing, ZOU Degao, et al. Study of dynamic pore water pressure and axial strain in saturated sand-gravel composites[J]. Rock and Soil Mechanics, 2006, 27(6): 925-928. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200606014.htm
[8]	CHEN G X, ZHAO D F, CHEN W Y, et al. Excess pore-water pressure generation in cyclic undrained testing[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2019, 145(7): 04019022. doi: 10.1061/(ASCE)GT.1943-5606.0002057
[9]	年廷凯, 焦厚滨, 范宁, 等. 南海北部陆坡软黏土动力应变-孔压特性试验[J]. 岩土力学, 2018, 39(5): 1564-1572, 1580. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201805003.htm NIAN Tingkai, JIAO Houbin, FAN Ning, et al. Experiment on dynamic strain-pore pressure of soft clay in the northern slope of South China Sea[J]. Rock and Soil Mechanics, 2018, 39(5): 1564-1572, 1580. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201805003.htm
[10]	HYODO M, HYDE A F L, ARAMAKI N. Liquefaction of crushable soils[J]. Géotechnique, 1998, 48(4): 527-543. doi: 10.1680/geot.1998.48.4.527
[11]	刘汉龙, 张宇, 郭伟, 等. 微生物加固钙质砂动孔压模型研究[J]. 岩石力学与工程学报, 2021, 40(4): 790-801. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202104012.htm LIU Hanlong, ZHANG Yu, GUO Wei, et al. A prediction model of dynamic pore water pressure for MICP-treated calcareous sand[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(4): 790-801. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202104012.htm
[12]	AIREY D W, FAHEY M. Cyclic response of calcareous soil from the North-West Shelf of Australia[J]. International Journal of Rock Mechanics & Mining Sciences & Geomechanics Abstracts 1991, 28: 101-121.
[13]	虞海珍, 汪稔. 钙质砂动强度试验研究[J]. 岩土力学, 1999, 20(4): 6-11. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX199904001.htm YU Haizhen, WANG Ren. The cyclic strength test research on calcareous sand[J]. Rock and Soil Mechanics, 1999, 20(4): 6-11. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX199904001.htm
[14]	虞海珍, 汪稔, 赵文光, 等. 波浪荷载下钙质砂孔压增长特性的试验研究[J]. 武汉理工大学学报, 2006, 28(11): 86-89. https://www.cnki.com.cn/Article/CJFDTOTAL-WHGY200611025.htm YU Haizhen, WANG Ren, ZHAO Wenguang, et al. Experimental research on development pattern of pore water pressure of carbonate sand under wave loads[J]. Journal of Wuhan University of Technology, 2006, 28(11): 86-89. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-WHGY200611025.htm
[15]	孙吉主, 黄明利, 汪稔. 内孔隙与各向异性对钙质砂液化特性的影响[J]. 岩土力学, 2002, 23(2): 166-169. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200202007.htm Sun Jizhu, Huang Mingli, Wang Ren. Influence of inner pore and anisotropy on liquefaction characteristics of calcareous sand[J]. Rock and Soil Mechanics, 2002, 23(2): 166-169. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200202007.htm
[16]	马维嘉, 陈国兴, 吴琪. 复杂加载条件下珊瑚砂抗液化强度试验研究[J]. 岩土力学, 2020, 41(2): 535-542, 551. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202002025.htm MA Weijia, CHEN Guoxing, WU Qi. Experimental study on liquefaction resistance of coral sand under complex loading conditions[J]. Rock and Soil Mechanics, 2020, 41(2): 535-542, 551. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202002025.htm
[17]	王刚, 查京京, 魏星. 循环三轴应力路径下钙质砂颗粒破碎演化规律[J]. 岩土工程学报, 2019, 41(4): 755-760. doi: 10.11779/CJGE201904020 WANG Gang, ZHA Jingjing, WEI Xing. Evolution of particle crushing of carbonate sands under cyclic triaxial stress path[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(4): 755-760. (in Chinese) doi: 10.11779/CJGE201904020
[18]	张家铭, 张凌, 蒋国盛, 等. 剪切作用下钙质砂颗粒破碎试验研究[J]. 岩土力学, 2008, 29(10): 2789-2793. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200810039.htm ZHANG Jiaming, ZHANG Ling, JIANG Guosheng, et al. Research on particle crushing of calcareous sands under triaxial shear[J]. Rock and Soil Mechanics, 2008, 29(10): 2789-2793. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200810039.htm
[19]	HALL E B, GORDON B B. Triaxial testing with large-scale high pressure equipment[J]. Laboratory Shear Testing of Soils, 1963, 361: 315-328.
[20]	WEI H Z, ZHAO T, HE J Q, et al. Evolution of particle breakage for calcareous sands during ring shear tests[J]. International Journal of Geomechanics, 2018, 18(2): 04017153.
[21]	王艳丽, 饶锡保, 潘家军, 等. 细粒含量对饱和砂土动孔压演化特性的影响[J]. 土木建筑与环境工程, 2011, 33(3): 52-56. https://www.cnki.com.cn/Article/CJFDTOTAL-JIAN201103010.htm WANG Yanli, RAO Xibao, PAN Jiajun, et al. Effects of fines content on evolutionary characteristics of dynamic pore water pressure of saturated sands[J]. Journal of Chongqing Jianzhu University, 2011, 33(3): 52-56. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JIAN201103010.htm
[22]	曾长女, 刘汉龙, 陈育民. 细粒含量对粉土动孔压发展模式影响的试验研究[J]. 岩土力学, 2008, 29(8): 2193-2198. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200808037.htm ZENG Changnü, LIU Hanlong, CHEN Yumin. Test study on influence of fine particle content on dynamic pore water pressure development mode of silt[J]. Rock and Soil Mechanics, 2008, 29(8): 2193-2198. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200808037.htm
[23]	吴杨, 崔杰, 李晨, 等. 细粒含量对岛礁吹填珊瑚砂最大动剪切模量影响的试验研究[J]. 岩石力学与工程学报, 2022, 41(1): 205-216. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202201017.htm WU Yang, CUI Jie, LI Chen, et al. Experimental study on the effect of fines on the maximum dynamic shear modulus of coral sand in a hydraulic fill island-reef[J]. Chinese Journal of Rock Mechanics and Engineering, 2022, 41(1): 205-216. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202201017.htm
[24]	HIGHT D W, GENS A, SYMES M J. The development of a new hollow cylinder apparatus for investigating the effects of principal stress rotation in soils[J]. Géotechnique, 1983, 33(4): 355-383.
[25]	Standard Test Methods for Minimum Index Density and Unit Weight of Soils and Calculation of Relative Density: ASTM D4254—14[S]. Annual Book of ASTM Standards, 2006.
[26]	Standard Test Methods for Maximum Index Density and Unit Weight of Soils Using A Vibratory Table: ASTM D4253—16 [S]. Annual Book of ASTM Standards, 2016.
[27]	SEED H B, LYSMER J, MARTIN P P. Pore-water pressure changes during soil liquefaction[J]. Journal of the Geotechnical Engineering, 1976, 102(4): 323-346.
[28]	LEE K L, ALBAISA A. Earthquake induced settlements in saturated sands[J]. Journal of the Geotechnical Engineering Division, 1974, 100(4): 387-406.
[29]	THEVANAYAGAM S, MARTIN G R. Liquefaction in silty soils-screening and remediation issues[J]. Soil Dynamics & Earthquake Engineering, 2002, 22(9): 1035-1042.
[30]	MOHAMMADI A, QADIMI A. A simple critical state approach to predicting the cyclic and monotonic response of sands with different fines contents using the equivalent intergranular void ratio[J]. Acta Geotechnica, 2015, 10(5): 587-606.