饱和砂-黏及砂-粉-黏混合土动孔压特性

李涛; 唐小微; 曾铃; 杨钢

doi:10.11779/CJGE2023S20044

饱和砂-黏及砂-粉-黏混合土动孔压特性 English Version

李涛^{1, 2,},
唐小微²,
曾铃¹,
杨钢²

1.
长沙理工大学土木工程学院, 湖南长沙 410114
2.
大连理工大学海岸和近海工程国家重点实验室, 辽宁大连 116024

基金项目:

国家自然科学基金项目 52308339

海岸和近海工程国家重点实验室开放基金项目 LP2110

桥梁工程安全控制教育部重点实验室开放基金项目 23ZDXK10

湖南省教育厅一般项目 22C0143

详细信息

作者简介:
李涛(1989—)，男，博士，讲师，主要从事混合砂土液化试验及理论等方面的研究工作。E-mail: litao@csust.edu.cn

中图分类号: TU435
计量
- 文章访问数: 137
- HTML全文浏览量: 28
- PDF下载量: 28
出版历程
- 收稿日期: 2023-11-29
- 网络出版日期: 2024-04-19
- 刊出日期: 2023-11-30

Dynamic pore water pressure characteristics of saturated sand-clay and sand-silt-clay mixtures

1.
School of Civil Engineering, Changsha University of Science & Technology, Changsha 410114, China
2.
State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China

摘要

摘要: 含细粒砂土在吹填土中广泛分布，其在长期服役过程中可能经受波浪及车辆荷载的循环作用，掌握含细粒砂土动孔隙水压力的变化特性对于评价吹填土在循环荷载作用下的安全性及稳定性尤为重要。通过对饱和重塑砂-黏混合土及砂-粉-黏混合土开展室内固结不排水循环三轴试验，研究了不同黏粒含量、细粒含量、黏粉比和孔隙比下砂-黏混合土及砂-粉-黏混合土动孔压的发展特征。结合不同混合土内部颗粒赋存状态光学显微镜图像及颗粒赋存状态理论，揭示了黏粒含量、细粒含量、黏粉比和孔隙比对混合土动孔压发展的影响规律及影响机理。试验结果表明：对于砂骨架孔隙比与孔隙比恒定的砂-黏混合土来说，随着黏粒含量的增加，动孔压在振动后期的波动性增长趋势分别增强、减弱；对于砂-粉-黏混合土来说，随着细粒中粉粒含量增多（黏粒含量减少），FC=5%试样动孔压在振动后期的波动性增长趋势呈现非单调变化规律，不同孔隙比下FC=10%试样动孔压在振动后期的波动性增长趋势单调增强；提出了适用于砂-黏混合土以及砂-粉-黏混合土动孔压发展的Mod-Baziar模型，该模型不仅扩展了动孔压的预测范围，而且具有更高的预测精度，其可以更好反映动孔压的整体发展趋势以及孔压在振动后期的细节变化；当采用Mod-Baziar模型分别对不同孔隙比下的纯砂、砂-黏混合土以及砂-粉-黏混合土动孔压进行预测时，模型中各个参数的取值范围依次增大。初步分析认为上述试验结果是由不同颗粒间复杂的黏结、摩擦及填充作用所致。
- 混合土 /
- 黏粒含量 /
- 细粒含量 /
- 黏粉比 /
- 动三轴试验 /
- 动孔压模型
Abstract: The sand-fines mixtures are widely distributed in hydraulic fill soil and may be subjected to the cyclic action of wave and vehicle loads during their long-term service. It is very important to know the dynamic pore water pressure changing characteristics of the sand-fines mixtures for evaluating the safety and stability of hydraulic fill soil under cyclic loading. The consolidated undrained cyclic triaxial tests are carried out on the saturated reconstituted sand-clay and sand-silt-clay mixtures to investigate the development features of dynamic pore water pressure of mixed soil with different clay contents, fines contents, clay silt ratios and void ratios. By combining with the optical microscope images of soil particles and the theory of soil particles occurrence modes, the influence laws and mechanisms of clay content, fines content, clay silt ratio and void ratio on the dynamic pore water pressure of mixed soils are revealed. It is find that for the sand-clay mixtures (constant e and constant e_s), with the increasing clay contents, the fluctuation growth trend of dynamic pore water pressure at the later stage of vibration increases and weakens respectively. For the sand-silt-clay mixtures, with the increasing silt contents (decreasing clay contents), the fluctuation growth trend of the dynamic pore water pressure at the later stage of vibration of specimens with FC of 5% changes nonmonotonously. While the fluctuation growth trend of the dynamic pore water pressure at the later stage of vibration of specimens with FC of 10% increases monotonically under different void ratios. The Mod-Baziar model applied to the development of dynamic pore water pressure of sand-clay mixtures and sand-silt-clay mixtures is proposed, which extends the prediction range of dynamic pore water pressure and has higher prediction accuracy. The Mod-Baziar model can better reflect the overall development trend and detail change at the later stage of dynamic pore water pressure. By comparing the values of each parameter in the Mod-Baziar model, it is found that when the Mod-Baziar model is used to predict the dynamic pore water pressure of sand, sand-clay mixtures and sand-silt-clay mixtures with different void ratios respectively, the value range of each parameter in the model increases successively. The preliminary analysis suggests that the above results are attributed to the complex bonding, friction and filling between different particles.
- mixed soil /
- clay content /
- fines content /
- clay silt ratio /
- dynamic triaxial test /
- dynamic pore water pressure model

HTML全文

0. 引言

场地土性质的改变会引起地震波的放大或者衰减，因此在进行地震响应的研究时要选择符合实际的场地模型。然而，以往对场地地震响应问题的研究，是将场地土简化为单相弹性地基或饱和土地基，且没有考虑热效应的影响。但是，在自然界中，非饱和土是土体更为常见的一种状态，且随着温室效应的加重，温度的变化越来越显著。此外，随着各国工业化和现代化进程的加快，以及热力学工程，化工工程和核废料管理等学科领域的大力发展，温度变化等热效应对多孔介质中波的传播特性的影响越来越受到众多学者的关注。由于热和力之间的相互耦合效应影响着波的形成和传播，并且多相介质之间的温度和应力-应变的耦合比经典情况下更复杂，这些因素叠加起来使得波的传播也更为复杂。因此，考虑热效应影响，建立一个能够反映实际场地土特性的非饱和土场地模型，并研究其地震地面运动具有重要意义。

自Biot^[1]首次提出流体饱和多孔弹性介质的动力学模型以来，针对等温条件下地震波在饱和多孔介质中的传播特性，已有大量的研究结果^[2-5]。伴随着混合物理论与多孔介质理论的发展，学者们^[6-7]结合Biot理论各自建立了非饱和多孔介质的动力方程，确定了多孔介质中存在3种压缩波和一种剪切波。另外，已有学者利用数值方法，对非饱和多孔介质的动力方程进行求解，并研究弹性波在非饱和土层分界面上的传播特性。Chen等^[8-9]研究了面波和P波在单相介质和非饱和多孔介质交界面上的波动特性，Tomar等^[10]研究了地震波在单相介质和非饱和多孔介质交界面上的透反射问题。

但是随着全球变暖趋势在逐步增强，在温度变化等热效应影响下，对于岩土介质中波的传播特性的研究逐渐进入到人们的视野中。早期的热弹性理论并不能分析有限速度热波的传播。Lord等^[11]基于前人的研究成果，建立了一种新型的热弹性理论，修正了这一缺陷。Pecker等^[12]研究了在温度变化条件下，地震波在饱和多孔弹性介质中的传播特性。Youssef^[13]结合前人的研究成果，进一步建立了广义的热弹性理论。柳鸿博等^[14]根据广义的热弹性模型，在考虑热效应和流体迂曲度的影响下，建立了非饱和多孔热弹性介质的波动方程。随后又有多位学者在前人研究的基础上，对波在多孔热弹性介质分界面上的透反射问题进行研究。郑荣跃等^[15]采用第一类Biot的理论模型，研究了平面SV波在饱和多孔热弹性介质分界面上的反射问题。柳鸿博等^[16]在热-流-固耦合动力响应模型的基础上，研究平面S波在饱和多孔热弹性介质边界上的反射特性。Zhou等^[17]结合多孔介质理论和广义热弹性模型，研究了平面SV波在非饱和多孔热弹性介质自由边界上的反射问题。Liu等^[18]采用了考虑孔隙流体与孔隙材料固相相互作用的非等温动力学模型，利用混合物理论分析了均匀平面P₁波在非饱和多孔热弹性介质边界处的反射特性。Wei等^[19]研究了各向同性、均匀的热弹性介质与多孔热弹性介质在平面界面上，由斜入射纵波所引起的反射和透射现象。Hou等^[20]研究了在任意入射角下，两种不相似的饱和多孔热弹性介质界面上，非均匀平面波的反射和透射特性。Liu等^[21]在热流-固体耦合的非饱和土波动控制方程的基础上，进一步研究了热弹性固体介质与非饱和多孔热弹性介质之间的平面界面上，由均匀平面P波入射所产生的反射和透射现象。

虽然目前国内外学者对多孔介质中波的传播，以及在多孔介质分界面上的波动响应问题做了一定的研究，但是对于地震波入射下场地的地震响应问题的研究较少。李伟华等^[22]首次将三相介质的动力模型引入到自由场地地震地面运动的分析上，研究了饱和度对地震运动的影响规律。随后，Li等^[23]基于非饱和多孔介质的波动特性，研究了平面P波和SV波对非饱和土层-基岩系统地面运动的影响。但是前人对场地地震响应的研究一直是基于等温的条件，鲜有文献对在考虑热效应作用下，平面P波入射到非饱和土自由场地对地震响应的影响问题进行研究，然而热效应对波的传播特性有着显著影响，且热效应对弹性波传播的影响关系到许多地震学和天体物理学问题，故而对考虑热效应作用下非饱和土自由场地地震响应分析的研究具有更为深远的意义。

本文对非饱和多孔热弹性介质波动方程以及单相热弹性介质波动方程进行求解，并建立非饱和土自由场地模型，采用亥姆霍兹矢量分解原理对场地波场进行分析，结合波的势函数和边界条件，求得土层中各波幅系数的解，进而确定场地中各点的位移和应力，最后研究各物理参数对场地地震响应所产生的影响。

1. 场地的波动方程

1.1 单相热弹性介质

基岩层采用单相热弹性介质来模拟。根据广义热弹性模型，单相热弹性介质的波动方程为^[21]

$\mu_{\mathrm{e}} \nabla^2 \boldsymbol{u}^{\mathrm{s}}+\left(\lambda_{\mathrm{e}}+\mu_{\mathrm{e}}\right) \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{s}}\right)-3 K_{\mathrm{b}} \beta_{\mathrm{s}} \nabla T=\rho^{\mathrm{e}} \ddot{\boldsymbol{u}}^{\mathrm{s}},$

(1a)

$\begin{gathered} 3 K_{\mathrm{b}} \beta_{\mathrm{s}} T_0 \nabla \cdot\left(\dot{\boldsymbol{u}}^s+\tau_{\mathrm{qe}} \ddot{\boldsymbol{u}}^{\mathrm{s}}\right)+\rho^{\mathrm{e}} c_{\mathrm{se}}\left(\dot{T}+\tau_{\mathrm{qe}} \ddot{T}\right) \\ =K_{\mathrm{e}} \nabla^2\left(T+\tau_{\theta \mathrm{e}} \dot{T}\right)。 \end{gathered}$

(1b)

式中：K_b= ${\lambda _{\text{e}}}$ +2μ_e/3为体积模量； ${\lambda _{\text{e}}}$ 和 ${\mu _{\text{e}}}$ 为单相介质的拉梅常量； ${\beta _{\text{s}}}$ 为热膨胀系数； ${\rho ^{\text{e}}}$ 为密度；c_se为固相比热； ${\tau _{{\text{qe}}}}$ 和 ${\tau _{\theta {\text{e}}}}$ 分别为热通量相位延迟时间和温度梯度相位延迟时间；T和 ${T_0}$ 分别为开尔文介质温度和初始温度；K_e为单相介质的热传导系数；u^s为固相位移矢量。

根据Helmholtz矢量的分解原理，固相位移矢量可分解为

$\boldsymbol{u}^{\mathrm{s}}=\nabla \psi_{\mathrm{s}}+\nabla \times \boldsymbol{H}_{\mathrm{s}}。$

(2)

式中：Ψ_s和H_s分别为固相的标量和矢量势函数。

固相介质以及热波的势函数可假设为

${\psi _{\text{s}}} = {A^{\text{s}}}\exp \left[ {{\text{i}}({k_{\text{p}}}x - \omega t)} \right] \text{，}$

(3a)

${\boldsymbol{H}_{\text{s}}} = {\boldsymbol{B}^{\text{s}}}\exp \left[ {{\text{i}}({k_{\text{s}}}x - \omega t)} \right] \text{，}$

(3b)

$T = {A^{\text{T}}}\exp \left[ {{\text{i}}({k_{\text{p}}}\boldsymbol{x} - \omega t)} \right] 。$

(3c)

式中：A，B分别为P波（包含T波）和S波的波幅系数；k_p和k_s分别为P波（包括T波）和S波的波数；ω表示频率。

将式（2）代入式（1a），（1b），并结合式（3a）~（3c），可得热弹性介质中体波的特征方程：

$\left| {\begin{array}{*{20}{c}} {{k_{11}}}&{{k_{12}}} \\ {{k_{21}}}&{{k_{22}}} \end{array}} \right| = 0 \text{，}$

(4a)

$\left| {{\rho ^e}{\omega ^2} - {\mu _{\text{e}}}k_{\text{s}}^2} \right| = 0 。$

(4b)

其中，式（4a），（4b）分别为热弹性介质中压缩波（包含热波）和剪切波的特征方程，求解可得单相热弹性介质中各热弹性波的波数。方程中的元素详见文献[21]。

1.2 非饱和多孔热弹性介质

根据广义热弹性理论，非饱和多孔热弹性介质的波动方程为^[14]

$\begin{aligned} & \mu \nabla^2 \boldsymbol{u}^{\mathrm{s}}+(\bar{\lambda}+\mu) \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{s}}\right)+D_1 \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{f}}\right)+ \\ & D_2 \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{a}}\right)+D_3 \nabla T=\rho \ddot{\boldsymbol{u}}^{\mathrm{s}}+\rho^1 \ddot{\boldsymbol{u}}^{\mathrm{f}}+\rho^{\mathrm{g}} \ddot{\boldsymbol{u}}^{\mathrm{a}}, \end{aligned}$

(5a)

$\begin{aligned} & B_1 \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{s}}\right)+B_2 \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{f}}\right)+ \\ & B_3 \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{a}}\right)+B_4 \nabla T=\rho^1 \ddot{\boldsymbol{u}}^{\mathrm{s}}+\vartheta^{\mathrm{i}} \ddot{\boldsymbol{u}}^{\mathrm{f}}+v^{\mathrm{l}} \dot{\boldsymbol{u}}^{\mathrm{f}}, \end{aligned}$

(5b)

$\begin{aligned} & B_5 \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{s}}\right)+B_6 \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{f}}\right)+ \\ & B_7 \nabla\left(\nabla \cdot \boldsymbol{u}^{\mathrm{a}}\right)+B_4 \nabla T=\rho^{\mathrm{g}} \ddot{\boldsymbol{u}}^{\mathrm{s}}+\vartheta^{\mathrm{g}} \ddot{\boldsymbol{u}}^{\mathrm{a}}+v^{\mathrm{g}} \dot{\boldsymbol{u}}^{\mathrm{a}}\text{，} \end{aligned}$

(5c)

$\begin{aligned} & C_1 \nabla \cdot\left(\dot{\boldsymbol{u}}^{\mathrm{s}}+\tau_{\mathrm{q}} \ddot{\boldsymbol{u}}^{\mathrm{s}}\right)+C_2 \nabla \cdot\left(\dot{\boldsymbol{u}}^{\mathrm{f}}+\tau_{\mathrm{q}} \ddot{\boldsymbol{u}}^{\mathrm{f}}\right)+ \\ & \quad C_3 \nabla \cdot\left(\dot{\boldsymbol{u}}^{\mathrm{a}}+\tau_{\mathrm{q}} \ddot{\boldsymbol{u}}^{\mathrm{a}}\right)+C_4\left(\dot{T}+\tau_{\mathrm{q}} \ddot{T}\right)=K \nabla^2\left(T+\tau_\theta \dot{T}\right) 。 \end{aligned}$

(5d)

式中：λ=λ+D₀，λ和μ表示Lame常数；u^f=nS^l(u^l-u^s)为土体孔隙中液相相对于固体骨架的相对位移；u^a=nS^g(u^g-u^s)为孔隙中气相相对于固体骨架的相对位移，其中u^s，u^l和u^g为固相、液相和气相的绝对位移向量，n为土体的孔隙率；S^l为孔隙内液态水的饱和度、S^g为孔隙内气体的饱和度；ρ为非饱和多孔介质的总密度，并且满足 $\rho = {\bar \rho ^{\text{s}}} + {\bar \rho ^{\text{l}}} + {\bar \rho ^{\text{g}}}$ ，其中 ${\bar \rho ^{\text{s}}}$ ， ${\bar \rho ^{\text{l}}}$ 和 ${\bar \rho ^{\text{g}}}$ 表示固-液-气三相的表观密度，且3个表观密度可以分别表示为 ${\bar \rho ^{\text{s}}}$ =(1-n)ρ^s， ${\bar \rho ^l}$ =nS^lρ^l， ${\bar \rho ^{\text{g}}}$ = nS^gρ^g，其中ρ^s，ρ^l，ρ^g分别为固相、液相和气相的密度； ${\tau _{\text{q}}}$ 为热通量的相位延迟时间；K为热传导系数； ${\tau _\theta }$ 为温度梯度的相位延迟时间。此外，其他符号的表达式见文献[14]。

根据Helmholtz势函数分解原理，三相位移矢量可以分解为

$\boldsymbol{u}^\alpha=\nabla \psi_\alpha+\nabla \times \boldsymbol{H}_\alpha 。$

(6)

式中： ${\psi _\alpha }$ 和 ${\boldsymbol{H}_\alpha }$ 为各相介质的势函数（α=s，f，a），s表示固相、f表示液相、a表示气相。

假设平面波的势函数解如下：

${\psi _\alpha } = {A^\alpha }\exp \left[ {{\text{i}}({k_{\text{p}}}\boldsymbol{x} - \omega t)} \right] \text{，}$

(7a)

${\boldsymbol{H}_\alpha } = {\boldsymbol{B}^\alpha }\exp \left[ {{\text{i}}({k_{\text{s}}}\boldsymbol{x} - \omega t)} \right] \text{，}$

(7b)

$T = {A^{\text{T}}}\exp \left[ {{\text{i}}({k_{\text{p}}}\boldsymbol{x} - \omega t)} \right] 。$

(7c)

式中： ${A^\alpha }$ 和 ${\boldsymbol{B}^\alpha }$ （α=s，f，a）分别为各相中P波（包括T波）和S波对应的势函数幅值；k_p和k_s分别为压缩波（包括热波）和剪切波的复波数；ω为角频率。

将式（6）代入式（5a）~（5d），并结合式（7a），（7b），（7c），可得到热弹性体波的特征方程如下：

$\left| {\begin{array}{*{20}{c}} {{b_{11}}}&{{b_{12}}}&{{b_{13}}}&{{b_{14}}} \\ {{b_{21}}}&{{b_{22}}}&{{b_{23}}}&{{b_{24}}} \\ {{b_{31}}}&{{b_{32}}}&{{b_{33}}}&{{b_{34}}} \\ {{b_{41}}}&{{b_{42}}}&{{b_{43}}}&{{b_{44}}} \end{array}} \right| = 0 \text{，}$

(8a)

$\left| {\begin{array}{*{20}{c}} {{c_{11}}}&{{c_{12}}}&{{c_{13}}} \\ {{c_{21}}}&{{c_{22}}}&{{c_{23}}} \\ {{c_{31}}}&{{c_{32}}}&{{c_{33}}} \end{array}} \right| = 0 。$

(8b)

式（8a），（8b）分别为压缩波（包括热波）和剪切波的特征方程，方程中的元素详见文献[14]。通过求解方程，可得非饱和土层中各热弹性波的波数。

2. 波场分析

假设频率为ω的平面P波，从基岩层以任意角度θ_ip入射至单相热弹性介质和非饱和多孔热弹性介质的分界面上时，在基岩层中产生如图 1所示的3种反射波（反射P波、反射S波、反射T波）和非饱和土层中5种透射波（透射P₁波、透射P₂波、透射P₃波、透射T波、透射S波），随着地震波在非饱和土层中继续传播，在非饱和土层和地表交界面处产生反射P₁波、反射P₂波、反射P₃波、反射T波、反射S波。

图 1 非饱和土自由场地简化模型

Figure 1. Simplified model for free site of unsaturated soils

下载: 全尺寸图片幻灯片

2.1 波的势函数

基岩层和非饱和土层中的波场函数如下所示：

（1）单相热弹性介质中体波的势函数

$\begin{gathered} \mathrm{P} \text { 波: } \psi_{\mathrm{e}}=A_{\mathrm{ip1}}^{\mathrm{s}} \exp \left[\mathrm{i} k_{\mathrm{ip1}}\left(l_{\mathrm{ip1}} x-n_{\mathrm{ip1}} z-c_{\mathrm{ip1}} t\right)\right]+ \\ \sum\limits_{q=1}^2 A_{\mathrm{1rpq} q}^{\mathrm{s}} \exp \left[\mathrm{i} k_{\mathrm{1rpq}}\left(l_{\mathrm{1rpq}} x+n_{\mathrm{1rpq} q} z-c_{\mathrm{1rpq}} t\right)\right], \end{gathered}$

(9a)

$\begin{gathered} &T 波：T_{\mathrm{e}}=A_{\mathrm{ip1}}^{\mathrm{s}} \delta_{\mathrm{Tp} 1}^{\mathrm{e}} \exp \left[\mathrm{i} k_{\text {ip1 }}\left(l_{\mathrm{ip1}} x-n_{\mathrm{ipl}} z-c_{\mathrm{ip1}} t\right)\right]+\\ &\sum\limits_{q=1}^2 A_{ \mathrm{1rpq}}^{\mathrm{s}} \delta_{\mathrm{Tpq}}^{\mathrm{e}} \exp \left[\mathrm{i} k_{\mathrm{1rpq}}\left(l_{\mathrm{1rpq}} x+n_{\mathrm{1rpq}} z-c_{\mathrm{1rpq}} t\right)\right], \end{gathered}$

(9b)

$S波： {H_{\text{e}}} = {B_{{\text{1rs}}}}\exp \left[ {{\text{i}}{k_{1rs}}({l_{{\text{1rs}}}}x + {n_{{\text{1rs}}}}z - {c_{{\text{1rs}}}}t)} \right] 。$

(9c)

（2）非饱和土层中体波的势函数

$\begin{aligned} & \mathrm{P} \text { 波: } \psi_{\mathrm{u}}^\alpha=\sum\limits_{n=1}^4 A_{\mathrm{tp} n}^\alpha \exp \left[\mathrm{i} k_{\mathrm{tp} n}\left(l_{\mathrm{tp} n} x-n_{\mathrm{tp} n} z-c_{\mathrm{tp} n} t\right)\right]+ \\ & \sum\limits_{n=1}^4 A_{2 \mathrm{rp} n}^\alpha \exp \left[\mathrm{i} k_{2 \mathrm{rp} n}\left(l_{2 \mathrm{rp} n} x+n_{2 \mathrm{rp} n} z-c_{2 \mathrm{rp} n} t\right)\right], \end{aligned}$

(10a)

$\begin{aligned} & \mathrm{T} \text { 波: } T_{\mathrm{u}}=\sum\limits_{n=1}^4 A_{\mathrm{tp} n}^{\mathrm{s}} \delta_{\mathrm{Tp} n} \exp \left[\mathrm{i} k_{\mathrm{t p} n}\left(l_{\mathrm{tp} n} x-n_{\mathrm{tp} n} z-c_{\mathrm{tp} n} t\right)\right]+ \\ & \sum\limits_{n=1}^4 A_{2 \mathrm{rp} n}^{\mathrm{s}} \delta_{\mathrm{Tp} n} \exp \left[\mathrm{i} k_{2 \mathrm{rp} n}\left(l_{2 \mathrm{rp} n} x+n_{2 \mathrm{rp} n} z-c_{2 \mathrm{rp} n} t\right)\right], \end{aligned}$

(10b)

$\begin{gathered} \mathrm{S} \text { 波: } H_{\mathrm{u}}^\alpha=B_{2 \mathrm{ts}}^\alpha \exp \left[\mathrm{i}k_{2 \mathrm{ts}}\left(l_{2 \mathrm{ts}} x-n_{2 \mathrm{ts}} z-c_{2 \mathrm{ts}} t\right)\right]+ \\ B_{2 \mathrm{rs}}^\alpha \exp \left[\mathrm{i} k_{2 \mathrm{rs}}\left(l_{2 \mathrm{rs}} x+n_{2 \mathrm{rs}} z-c_{2 \mathrm{rs}} t\right)\right] \quad 。 \end{gathered}$

(10c)

式中： $\alpha$ 取s，f，a，表示固、液、气三相；下标i，r，t分别表示入射波、反射波和透射波；下标e和1代表基岩层，u和2代表上覆非饱和土层；q表示单相热弹性介质中不同的反射和透射P波(包括T波)，其中q=1，2，分别代表单相热弹性介质中的反射P波、反射T波；n表示的是非饱和土层中不同的反射和透射P波（包括T波），其中n=1，2，3，4，分别代表反射（透射）P₁波、反射（透射）P₂波、反射（透射）P₃波和反射（透射）T波。k_ip1，c_ip1分别表示入射P波的波数和波速；k_1rpq，k_1rs，c_1rpβ，c_1rs分别表示半无限基岩层中2种反射压缩波（P波、T波）和1种反射剪切波（S波）的波数以及对应的波速；k_tpn，k_2ts分别表示非饱和土层中4种透射波（包括P波和T波）和1种透射S波的波数；c_tpn，c_2ts分别表示非饱和土层中4种透射纵波（包括热波）和1种透射S波的波速；k_2rpn、k_2rs分别表示非饱和土层中4种反射波（包括P波和T波）和1种反射S波的波数；c_2rpn，c_2rs分别表示非饱和土层中4种反射波（包括P波和T波）和1种反射S波的波速；l，n分别表示x，z方向的方向矢量； ${\text{i}} = \sqrt { - 1}$ 。

由Snell定理可得

$\begin{aligned} l_{\mathrm{ip1} } k_{\mathrm{ip} 1} & =l_{1 \mathrm{rp1} } k_{1 \mathrm{p} 1}=l_{\mathrm{1rp} 2} k_{1 \mathrm{rp} 2}=l_{\mathrm{1rs}} k_{1 \mathrm{rs}}=l_{\mathrm{tp1} } k_{\mathrm{tp} 1} \\ & =l_{\mathrm{tp} 2} k_{\mathrm{tp} 2}=l_{\mathrm{tp} 3} k_{\mathrm{tp} 3}=l_{\mathrm{tp} 4} k_{\mathrm{tp} 4}=l_{2 \mathrm{ts}} k_{2 \mathrm{ts}} 。 \end{aligned}$

(11)

由式（4a），（4b）经过推导可得单相热弹性介质中P波和T波的势函数幅值比例为

${\delta }_{\text{Tp}n}^{\text{e}}=\frac{{A}_{\text{T}}}{{A}_{\text{s}}}=-\frac{{k}_{11}}{{k}_{12}}(n=1,2\text{) }。$

(12)

由式（8a），（8b）经过推导可得非饱和多孔热弹性介质中势函数的幅值比例关系为

${\delta }_{\text{fp}n}=\frac{{A}_{}^{\text{f}}}{{A}_{}^{\text{s}}}=\frac{{d}_{11}^{}{d}_{15}^{}-{d}_{12}^{}{d}_{16}^{}}{{d}_{13}^{}{d}_{16}^{}-{d}_{14}^{}{d}_{15}^{}}\text{ }(n=1,2,3,4)\text{ }\text{，}$

(13a)

${\delta }_{\text{ap}n}=\frac{{A}_{}^{\text{a}}}{{A}_{}^{\text{s}}}=\frac{{d}_{21}^{}{d}_{25}^{}-{d}_{22}^{}{d}_{26}^{}}{{d}_{23}^{}{d}_{26}^{}-{d}_{24}^{}{d}_{25}^{}}\text{ }(n=1,2,3,4)\text{ }\text{，}$

(13b)

${\delta }_{\text{Tp}n}=\frac{{A}_{}^{\text{T}}}{{A}_{}^{\text{s}}}=\frac{{d}_{31}^{}{d}_{35}^{}-{d}_{32}^{}{d}_{36}^{}}{{d}_{33}^{}{d}_{36}^{}-{d}_{34}^{}{d}_{35}^{}}\text{ }(n=1,2,3,4)\text{ }\text{，}$

(13c)

${\delta }_{\text{fs}}=\frac{{B}_{}^{\text{f}}}{{B}_{}^{\text{s}}}=-\frac{{c}_{12}}{{c}_{22}},{\delta }_{\text{as}}=\frac{{B}_{}^{a}}{{B}_{}^{\text{s}}}=\frac{{c}_{21}{c}_{12}-{c}_{11}{c}_{22}}{{c}_{13}{c}_{22}}\text{ }。$

(13d)

式中：元素d₁₁~d₃₆的表达式详见文献[21]。

2.2 边界条件

（1）为保证非饱和土层与基岩层的边界处（z=H）紧密接触，两种不同介质边界处的应力和位移是连续的。此外，该边界在两种不同介质下的等温条件通过温度和热通量的连续性来保证。由于不可渗透的热弹性介质的存在，非饱和多孔热弹性介质中的孔隙流体不可渗透通过该分界面，故适用于该界面的边界条件可以表示如下^[21]：

交界面上应力连续：

${\sigma }_{zz}={\sigma }_{zz}^{\text{e}},{\sigma }_{xz}={\sigma }_{xz}^{\text{e}}\text{ }。$

(14a)

交界面上位移连续：

${u}_{z}^{\text{s}}={u}_{z}^{\text{e}},{u}_{x}^{\text{s}}={u}_{x}^{\mathrm{e}}\text{ }。$

(14b)

交界面上温度以及温度变化率连续：

$T={T}^{\text{e}},K\frac{\partial T}{\partial z}={K}^{\text{e}}\frac{\partial {T}^{\text{e}}}{\partial z}\text{ }。$

(14c)

交界面上液相、气相相对于固体骨架的法向位移为0：

${u}_{z}^{\text{f}}={u}_{z}^{\text{a}}=0\text{ }。$

(14d)

（2）考虑自由地表为自由透水、透气和绝热的边界。由于地球表面弹性波的反射，一般认为该边界表面没有外力，孔隙内的液体和气体可以在边界表面自由流动，介质内外的温度在边界表面绝缘^[17-18]，因此适用于该界面的边界条件如下所示：

交界面上正应力、剪应力、液相压力、气相压力以及温度变化率为0

${\sigma _{zz}} = {\sigma _{xz}} = {p_{\text{w}}} = {p_{\text{a}}} = \frac{{\partial T}}{{\partial z}} = 0。$

(15)

单相介质中的应力位移表达式为

$u_x^{\text{e}} = \frac{{\partial {\psi _{\text{e}}}}}{{\partial x}} - \frac{{\partial {H_{\text{e}}}}}{{\partial z}} \text{，}$

(16a)

$u_z^{\text{e}} = \frac{{\partial {\psi _{\text{e}}}}}{{\partial z}} + \frac{{\partial {H_{\text{e}}}}}{{\partial x}} \text{，}$

(16b)

$\sigma _{xz}^{\text{e}} = {\mu _{\text{e}}}\left( {2\frac{{{\partial ^2}{\psi _{\text{e}}}}}{{\partial x\partial z}} + \frac{{{\partial ^2}{H_{\text{e}}}}}{{\partial {x^2}}} - \frac{{{\partial ^2}{H_e}}}{{\partial {z^2}}}} \right) \text{，}$

(16c)

${\sigma }_{zz}^{\text{e}}={\lambda }_{\text{e}}\left(\frac{{\partial }^{2}{\psi }_{\text{e}}}{\partial {x}^{2}}+\frac{{\partial }^{2}{\psi }_{\text{e}}}{\partial {z}^{2}}\right)+2{\mu }^{\text{e}}\left(\frac{{\partial }^{2}{\psi }_{\text{e}}}{\partial {z}^{2}}+\frac{{\partial }^{2}{H}_{\text{e}}}{\partial x\partial z}\right)-3{K}_{\text{b}}{\beta }_{\text{T}}T\text{ }。$

(16d)

非饱和多孔介质中的应力位移表达式为

$u_x^\alpha = \frac{{\partial \psi _{\text{u}}^\alpha }}{{\partial x}} - \frac{{\partial H_{\text{u}}^\alpha }}{{\partial z}} \text{，}$

(17a)

$u_z^\alpha = \frac{{\partial \psi _\rm u^a}}{{\partial z}} + \frac{{\partial H_{\text{u}}^\alpha }}{{\partial x}} \text{，}$

(17b)

$\begin{gathered} \sigma_{z z}=\bar{\lambda}\left(\frac{\partial^2 \psi_\rm u^{\mathrm{s}}}{\partial x^2}+\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{s}}}{\partial z^2}\right)+2 \mu\left(\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{s}}}{\partial z^2}+\frac{\partial^2 H_{\mathrm{u}}^{\mathrm{s}}}{\partial x \partial z}\right)+ \\ D_1\left(\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{f}}}{\partial x^2}+\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{f}}}{\partial z^2}\right)+D_2\left(\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{a}}}{\partial x^2}+\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{a}}}{\partial z^2}\right)+D_3 T_{\mathrm{u}}, \end{gathered}$

(17c)

${\sigma _{xz}} = \mu \left( {2\frac{{{\partial ^2}\psi _{\text{u}}^{\text{s}}}}{{\partial x\partial z}} + \frac{{{\partial ^2}H_{\text{u}}^{\text{s}}}}{{\partial {x^2}}} - \frac{{{\partial ^2}H_{\text{u}}^{\text{s}}}}{{\partial {z^2}}}} \right) \text{，}$

(17d)

$\begin{aligned} p_g= & -B_5\left(\frac{\partial^2 \psi_\rm u^{\mathrm{s}}}{\partial x^2}+\frac{\partial^2 \psi_\rm u^{\mathrm{s}}}{\partial z^2}\right)-B_6\left(\frac{\partial^2 \psi_\rm u^{\mathrm{f}}}{\partial x^2}+\frac{\partial^2 \psi_\rm u^{\mathrm{f}}}{\partial z^2}\right)- \\ & B_7\left(\frac{\partial^2 \psi_\rm u^{\mathrm{a}}}{\partial x^2}+\frac{\partial^2 \psi_\rm u^{\mathrm{a}}}{\partial z^2}\right)-B_8 T_\rm u\text{，} \end{aligned}$

(17e)

$\begin{aligned} p_{\mathrm{w}}= & -B_1\left(\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{s}}}{\partial x^2}+\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{s}}}{\partial z^2}\right)-B_2\left(\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{f}}}{\partial x^2}+\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{f}}}{\partial z^2}\right)- \\ & B_3\left(\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{a}}}{\partial x^2}+\frac{\partial^2 \psi_{\mathrm{u}}^{\mathrm{a}}}{\partial z^2}\right)-B_4 T_{\mathrm{u}}。 \end{aligned}$

(17f)

将波场函数（9a）~（9c）以及（10a）~（10c）代入边界条件（14a）~（14d）以及（15）可求得如下线性方程组：

$\boldsymbol{F N}=A_{\mathrm{ip}} \boldsymbol{G}。$

(18)

式中：

$N={\left[{A}_{\text{1rp1}}^{\text{s}},{A}_{\text{1rp2}}^{\text{s}},{B}_{\text{1rs}}^{\text{s}},{A}_{\text{tp1}}^{\text{s}},{A}_{\text{tp2}}^{\text{s}},{A}_{\text{tp3}}^{\text{s}},{A}_{\text{tp4}}^{\text{s}},{B}_{\text{2ts}}^{\text{s}},{A}_{\text{2rp1}}^{\text{s}},{A}_{\text{2rp2}}^{\text{s}},{A}_{\text{2rp3}}^{\text{s}},{A}_{\text{2rp4}}^{\text{s}},{B}_{\text{2rs}}^{\text{s}}\right]}^{T}，$ 矩阵F和矩阵G中的元素在附录A中给出。

3. 自由地表的位移

波幅系数确定后，场地各点的应力和位移均可确定。将式（10a）~（10c）代入到水平方向位移u_x和竖直方向位移u_z的表达式中可得

$\begin{aligned} u_x= & A_{\mathrm{tp1}}^{\mathrm{s}} \mathrm{i} k_{\mathrm{tp} 1} l_{\mathrm{tp} 1}+A_{\mathrm{tp} 2}^{\mathrm{s}} \mathrm{i} k_{\mathrm{tp} 2} l_{\mathrm{tp} 2}+A_{\mathrm{tp} 3}^{\mathrm{s}} k_{\mathrm{tp} 3} l_{\mathrm{tp} 3}+A_{\mathrm{tp} 4}^{\mathrm{s}} \mathrm{i} k_{\mathrm{tp} 4} l_{\mathrm{tp} 4}+ \\ & A_{2 \mathrm{rp} 1}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rp} 1} l_{2 \mathrm{rp} 1}+A_{2 \mathrm{rp} 2}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rp} 2} l_{2 \mathrm{rp} 2}+A_{2 \mathrm{r} 3}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rp} 3} l_{2 \mathrm{rp} 3}+ \\ & A_{2 \mathrm{rp} 4}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rp} 4} l_{2 \mathrm{rp} 4}+B_{2 \mathrm{ts}}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{ts}} n_{2 \mathrm{ts}}-B_{2 \mathrm{rs}}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rs}} n_{2 \mathrm{rs}}, \end{aligned}$

(19a)

$\begin{aligned} u_z= & -A_{\mathrm{tp} 1}^{\mathrm{s}} \mathrm{i} k_{\mathrm{tp} 1} n_{\mathrm{tp} 1}-A_{\mathrm{tp} 2}^{\mathrm{s}} \mathrm{i} k_{\mathrm{tp} 2} n_{\mathrm{tp} 2}-A_{\mathrm{tp} 3}^{\mathrm{s}} \mathrm{i} k_{\mathrm{tp} 3} n_{\mathrm{tp} 3}-A_{\mathrm{tp} 4}^{\mathrm{s}} k_{\mathrm{tp} 4} n_{\mathrm{tp} 4}+ \\ & A_{2 \mathrm{rp} 1}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rp} 1} n_{2 \mathrm{rp} 1}+A_{2 \mathrm{rp} 2}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rp} 2} n_{2 \mathrm{rp} 2}+A_{2 \mathrm{rp} 3}^{\mathrm{s}} k_{2 \mathrm{rp} 3} n_{2 \mathrm{rp} 3}+ \\ & A_{2 \mathrm{rp} 4}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rp} 4} n_{2 \mathrm{rp} 4}+B_{2 \mathrm{ts}}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{ts}} l_{2 \mathrm{ts}}+B_{2 \mathrm{rs}}^{\mathrm{s}} \mathrm{i} k_{2 \mathrm{rs}} l_{2 \mathrm{rs}}。 \end{aligned}$

(19b)

本文采用自由场表面的水平和竖向位移放大系数(u_x/u₀和u_z/u₀)来表征非饱和土自由场地的表面位移，其中位移放大系数用对应方向的位移幅值与入射波的位移幅值u₀之比来表示。

4. 数值算例

4.1 验证

本文的模型是基岩上覆厚度为H的非饱和土自由场地简化模型。在厚度为H的非饱和土层中，当平面P波从基岩入射，各热弹性波将在基岩和非饱和土层分界面处，以及非饱和土层和地表交界面处产生多种反射波和透射波。由P波入射激发的自由场中的平面波如图 1所示。为了验证本文求解过程的正确性，令本文模型中的非饱和土层厚度H→∞，此时本文的模型退化为与文献[21]一致的模型，即P波在热弹性介质与非饱和多孔热介质交界面上的透反射模型，非饱和土层处于半空间中，同时本文的解可以退化为当P波入射到热弹性介质与非饱和多孔热弹性介质交界面时，各反射波和透射波振幅反射率和透射率的解。在验证计算中取与上述文献[21]相同的物理参数，通过与上述文献的对比，由图 2可以看出，本文的计算结果与文献解答相比，透射P₁波和透射S波的振幅透射率的计算结果基本一致，进而说明了本文方法的正确性。

图 2 本文数值解与文献解的比较

Figure 2. Comparison between numerical solutions in this study and literature solutions

下载: 全尺寸图片幻灯片

随后通过数值计算，研究热效应作用下P波入射非饱和土自由场地的地震响应问题，具体分析了非饱和土层中的介质温度等热物性参数对场地地震响应所产生的影响，在数值计算中所取的热物性参数^[21]如表 1，2所示。本文分析中采用无量纲频率ω₀，其中ω₀=ω/ω₁，ω₀=1，ω₁为上覆土层完全饱和时的固有频率, 土层厚度H=20 m。

表 1 单相热弹性介质的物理力学参数

Table 1. Physical-mechanical parameters of single-phase thermoelastic media

参数	λ_e/Pa	μ_e/Pa	ρ^e/(kg·m^-3)	β_s/K^-1	K_e/(J·s^-1·m^-1·K^-1)
量值	12.0×10⁹	8.0×10⁹	2700	4.0×10^-4	2.5

参数	T_e/K	τ_qe/s	τ_θe/s	T₀/K	c_se/(J·kg^-1·K^-1)
量值	293.2	2.0×10^-7	1.5×10^-7	300	1000

下载: 导出CSV

| 显示表格

表 2 非饱和热弹性介质的物理力学参数

Table 2. Physical-mechanical parameters of unsaturated thermoelastic media

参数	λ/Pa	μ/Pa	ρ^s/(kg·m^-3)	β_T/K^-1	K/(J·s^-1·m^-1·K^-1)
量值	4.4×10⁹	2.8×10⁹	2650	1.0×10^-4	0.6

参数	T/K	τ_q/s	τ_θ/s	T/K	c_s/(J·kg^-1·K^-1)
量值	293.2	2.0×10^-7	1.5×10^-7	300	1000

参数	c_l/(J·kg^-1·K^-1)	c_g/(J·kg^-1·K^-1)	β_ST/K^-1	β_WT/K^-1	β_WP/(Pa^-1)
量值	4180	1900	7.8×10^-6	2.1×10^-4	4.58×10^-10

参数	S^l	S_sat	S_res	n	k/m²
量值	0.6	1.0	0.05	0.4	1.0×10^-10

参数	ρ^l/(kg·m^-3)	K_s/kPa	P_g/kPa	m	d
量值	1000	3.5×10⁷	101.3	0.5	2

下载: 导出CSV

| 显示表格

4.2 热传导系数的影响

热传导系数表征了岩石的导热性，不仅是研究地壳热状态的重要参数，也被广泛应用于地震勘探和采矿等领域。为了研究不同热传导系数下和等温模型下地表位移放大系数随入射角的变化曲线，忽略热效应的影响，令波动控制方程中温度的变化量以及温度随时间的导数为0，此时本文的热弹性模型退化成等温模型，即本文的解退化成等温模型下的解^[24]。图 3给出了不同热传导系数下地表位移放大系数随入射角的变化曲线。首先从图 3中可以看出在考虑热效应和不考虑热效应两种理论模型下所得到的位移放大系数有着明显差异，在热效应条件下，水平和竖向位移放大系数的变化趋势与无热效应时的趋势相似，但热效应作用下的水平和竖向位移放大系数较小。这是因为介质中的温度与应力应变之间的耦合效应相比于等温情况下更为复杂。在等温模型中，温度变化量和温度梯度及温度随时间的导数均为0，且忽略了热传导方程，而热波（T波）又是随着热的传导而产生的。由于热波的出现，热效应条件下波的波速、透反射的能量有别于等温时的情况，使得热效应和无热效应两种理论模型下波的传播特性有明显的不同，然而波的传播特性又影响着场地的地震响应。因此，这些因素相互叠加导致了热效应条件下位移放大系数与等温条件下有显著不同。其次，从图 3中可以发现，地表水平和竖向位移放大系数均随着入射角的增大，先增大后减小，水平和竖向位移放大系数分别在入射角θ=45°和θ=5°时达到峰值。

图 3 不同热传导系数下位移放大系数随入射角度变化的曲线

Figure 3. Variation curves of displacement magnification factor with incident angle under different thermal conductivities

下载: 全尺寸图片幻灯片

随着热传导系数的增大，地表水平位移放大系数逐渐增大，而竖向位移放大系数几乎不发生变化。图 3所表现出的位移放大系数在热效应和无热效应条件下变化趋势的一致性以及两者在数值上的差异，既验证了本文结果的有效性，又说明在非等温条件下研究自由场地地震响应具有重要意义。

4.3 饱和度的影响

为了分析热效应下饱和度变化对地表位移放大系数的影响规律，图 4给出了不同饱和度下地表位移放大系数随入射角的变化曲线。首先从图 4可知，随着入射角的增大，水平位移放大系数和竖向位移放大系数均呈现先增加后减小的趋势，分别在θ=45°和θ=5°时达到峰值，另外当入射角θ=0°时，地表水平位移放大系数为0，而当入射角θ=90°时，地表水平位移和竖向位移放大系数均为0。其原因是：平面P波的能量主要沿着传播的方向传递，当平面P波沿着基岩和非饱和土层的分界面水平入射时（θ=90°），由于没有垂直于分界面的分量，非饱和土层中的反射和透射波消失^[21]，故非饱和土表面不会产生位移。其次，从图 4可以发现随着饱和度的增大，地表水平位移放大系数逐渐减小，竖向位移放大系数逐渐增大。其原因是：P₂波是由液相的作用产生的，P₃波是由孔隙液体和气体的压力差产生的，当土体接近完全饱和时，土体孔隙中液相的作用力达到最大，P₂波的波速达到最大值，而随着气体的不断消散，P₃波的波速减小到最小值。此外，入射P波的能量主要分配给P₁波和S波，而多孔介质孔隙中存在的流体相会阻碍P₁波和S波的传播，液相的阻碍能力强于气相^[24]。因此这些因素的共同作用使得不同饱和度下的土体，由于气体含量的不同对场地地震地面运动的影响也有所不同。此外，本文所得的这一变化趋势与文献[25]所得到的等温条件下饱和度变化对地表位移放大系数的影响规律一致。由此可见，图 4所表现出的地表位移放大系数随饱和度变化的规律，再一次验证了本文计算结果的有效性，又说明非饱和土中气相的存在对地表位移放大系数的影响不可忽视。把场地土视为饱和土地基从而没有考虑土体饱和度对场地地震地面运动的影响，这与实际现象是有出入的。

图 4 不同饱和度下位移放大系数随入射角度变化的曲线

Figure 4. Variation curves of displacement amplification coefficient with incident angle under different saturation

下载: 全尺寸图片幻灯片

4.4 热膨胀系数的影响

热膨胀系数主要涉及描述固相应力应变关系的物理方程和描述热传导的热传导方程，为了分析热膨胀系数变化对地表位移放大系数的影响规律，图 5给出了不同热膨胀系数下位移放大系数随入射角的变化曲线。首先从图 5可知，地表水平位移放大系数和竖向位移放大系数随着入射角度的增大，均呈现先增加后减小的趋势。其次，从图 5可以看出，随着热膨胀系数的增大，水平位移放大系数逐渐增大，竖向位移放大系数也逐渐减小，且变化的幅值较大。Liu等^[21]研究发现反射P波、反射S波以及透射P₁波和透射S波所占据的能量较多，且当热膨胀系数β_T增加1×10^-4 K⁻¹时，透射P₁波的振幅比和能量比分别增加约5%和3%，透射S波的振幅比和能量比分别减少约0.5%和0.07%，而反射P波的振幅比和能量最多分别下降了约3%和4%，但反射S波的振幅比和能量比几乎不受影响。这些因素的共同作用产生了图 5所示的不同热膨胀系数下地表位移放大系数随入射角度变化的曲线。因此，热膨胀系数对地表位移放大系数有明显的影响，在工程选址，大型小区建设等实际工程中，应关注热膨胀系数的影响。

图 5 不同热膨胀系数下位移放大系数随入射角度变化的曲线

Figure 5. Variation curve of amplification factor of displacement with incident angle for different coefficients of thermal expansion

下载: 全尺寸图片幻灯片

4.5 介质温度的影响

为了分析介质温度变化对地表位移放大系数的影响规律，图 6给出了不同介质温度下地表位移放大系数随入射角的变化曲线。首先由图 6可知，地表水平和竖向位移放大系数均随着入射角的增大先增大后减小，当温度较低时水平和竖向位移放大系数分别在入射角θ=45°和θ=5°达到峰值。当介质温度达到313.2 K时，水平位移放大系数在θ=55°时达到峰值。其次，从图 6可以发现，随着介质温度的不断增大，地表水平位移放大系数逐渐减小，竖向位移放大系数逐渐增大，其中水平位移放大系数减小的幅值较大而竖向位移放大系数增大的幅值较小。其原因是：介质温度的变化，影响着各热弹性波的传播特性，从而影响着场地的地震地面运动。介质温度的变化对非饱和土层中的透射P₂波、透射P₃波、透射T波和反射T波的能量比和振幅比均产生显著影响，即随着温度的增大，其振幅比和能量比显著增加。但是，随着介质温度的增大，透射P₁波和S波以及反射的P波和S波的振幅比和能量比基本保持不变^[21]。这些因素的共同作用使得介质温度增大时，水平位移放大系数减小的幅值较大，而竖向位移放大系数增大的幅值较小。因此，在实际工程中，应关注土层的温度变化情况。

图 6 不同温度下位移放大系数随入射角度变化的曲线

Figure 6. Variation curve of displacement amplification coefficient with incident angle at different temperatures

下载: 全尺寸图片幻灯片

4.6 热通量相位延迟时间的影响

相对于等温条件下的多孔介质理论，在多孔热弹性介质理论下将会产生T波。热通量相位延迟时间作用于决定热波（T波）的波动方程式（5d），对T波的波速产生影响。为了分析热通量相位延迟时间变化对地表位移放大系数的影响规律，保持其他参数不变，图 7给出了不同热通量相位延迟时间下地表位移放大系数随入射角的变化曲线。

图 7 不同热通量相位延迟时间下位移放大系数随入射角度变化的曲线

Figure 7. Variation curves of displacement amplification coefficient with incident angle under different phase lags of heat flux

下载: 全尺寸图片幻灯片

首先从图 7可知，地表水平位移放大系数和竖向位移放大系数随着入射角度的增大，均呈现先增加后减小的趋势。其次，从图 7中可以发现，随着热通量相位延迟时间的增大，水平位移放大系数和竖向位移放大系数仅发生微小变化。其原因是：热通量相位延迟时间仅作用于决定T波的波动控制方程式（5d），对P波和S波的传播特性几乎不产生影响，且T波占据的能量和振幅比P波和S波的小好几个数量级。因此，热通量相位延迟时间的变化对位移放大系数的影响较小。

4.7 土层深度变化对土层应力的影响

为了分析土层深度变化对土层应力的影响，图 8给出了土层竖向应力随土层深度变化的曲线。图 8考虑了4种土层深度，分别为H=5，10，15，20 m。首先，应力的幅值随入射角的增大而增大，在入射角θ=58°时达到峰值。其次，随着土层深度H的增大，竖向应力逐渐增大。值得注意的是，当H=5 m时，相比于其他土层深度的情况，此时的应力较小。由此可见，在热效应条件下，土层的深度变化对应力变化的影响显著，当土层深度较小时，应力较小，当土层深度较大时，应力较大。因此，在实际工程抗震设计中，要关注土层深度的影响。

图 8 土层竖向应力随深度变化的曲线

Figure 8. Variation curves of vertical stress of soil layers with depth

下载: 全尺寸图片幻灯片

5. 结论

本文基于非饱和多孔热弹性介质和单相热弹性介质中波的传播理论，研究平面P波入射非饱和土自由场地的地震响应，分析了热传导系数等热物性参数变化对场地地震响应所产生的影响，主要得到以下3点结论。

（1）有热效应和无热效应两种理论模型下所得到的位移放大系数有显著的差距。热膨胀系数的增大将引起水平和竖向位移放大系数分别增大和减小。

（2）随着介质温度的增大，地表水平位移放大系数逐渐减小，竖向位移放大系数逐渐增大。热通量相位延迟时间对地表位移放大系数的影响较小。此外，随着土层深度的增大，土层中的竖向应力逐渐增大。

（3）随着饱和度的增大，地表水平位移放大系数逐渐减小，而竖向位移放大系数逐渐增大，说明非饱和土中气相对场地地震响应的影响不可忽视。

本文的研究结果可以为干旱、半干旱、土坝、铁路、港口平台以及海相沉积土分布等地区提供抗震理论指导，以便给地震小区划分和工程选址等工作提供理论依据，从而防范温室效应作用下场地的工程地震灾变风险，是保障工程抗震设计可靠性的基础科学手段，对工程应用具有重要的指导价值和现实意义。

附录A

$\begin{aligned} & \left.f_{11} \sim f_{12}=\left[-\left(\lambda^{\mathrm{e}}+2 \mu^{\mathrm{e}} n_{1 r p n}^2\right) k_{1 r p n}^2-3 k_b^{\mathrm{e}} \beta_T^{\mathrm{e}} \delta_{T p q}^{\mathrm{e}}\right)\right] \exp \left(\mathrm{i} k_{1 r p q} n_{1 r p q} H\right), \\ & f_{13}=-2 \mu^{\mathrm{e}} l_{1 r s} n_{1 r s} k_{1 r s}^2 \exp \left(\mathrm{i} k_{1 r s} n_{1 r s} H\right), \end{aligned}$

$\begin{aligned} & f_{14} \sim f_{17}=\left[\left(\bar{\lambda}+2 \mu n_{t p n}^2+D_1 \delta_{f p n}+D_2 \delta_{a p n}\right) k_{t p n}^2-D_3 \delta_{T p n}\right] . \\ & \exp \left(-\mathrm{i} k_{t p n} n_{t p n} H\right), \\ & f_{18}=-2 \mu l_{2 t s} n_{2 t s} k_{2 t s}^2 \exp \left(-\mathrm{i} k_{2 t s} n_{2 t s} H\right), \end{aligned}$

$\begin{aligned} & f_{19} \sim f_{1(12)}=\left[\left(\bar{\lambda}+2 \mu n_{2 r p n}^2+D_1 \delta_{2 r p n}+D_2 \delta_{2 r p n}\right) k_{2 r p n}^2-D_3 \delta_{Tp n}\right] . \\ & \exp \left(\mathrm{i} k_{2 r p n} n_{2 r pn} H\right), \\ & f_{1(13)}=2 \mu l_{2 r s} n_{2 r s} k_{2 r s}^2 \exp \left(\mathrm{i} k_{2 r s} n_{2 r s} H\right) \text {, } \\ & f_{21} \sim f_{22}=2 \mu^{\mathrm{e}} l_{1 rp q} n_{1 r p q} k_{1 r p q}^2 \exp \left(\mathrm{i} k_{1 r p q} n_{1 r p q} H\right), \\ & f_{23}=\mu^{\mathrm{e}}\left(l_{1 r s}^2-n_{1 r s}^2\right) k_{1 r s}^2 \exp \left(\mathrm{i}_{1 r s} n_{1 r s} H\right), \\ & f_{24} \sim f_{27}=2 \mu l_{ {tpn }} n_{ {tpn }} k_{ {tpn }}^2 \exp \left(-\mathrm{i} k_{ {tpn }} n_{ {tpn }} H\right), \\ & f_{28}=\mu\left(n_{2 t s}^2-l_{2 ts}^2\right) k_{2 ts}^2 \exp \left(-\mathrm{i} k_{2 ts} n_{2 ts} H\right) \text {, } \\ & f_{29} \sim f_{2(12)}=-2 \mu l_{2 r p n} n_{2 r pn} k_{2 { rpn }}^2 \exp \left(\mathrm{i} k_{2 r pn } n_{2 r p n} H\right) \text {, } \\ & f_{2(13)}=\mu\left(n_{2 r s}^2-l_{2 r s}^2\right) k_{2 r s}^2 \exp \left(\mathrm{i} k_{2 r s} n_{2 rs} H\right) \text {, } \\ & f_{31} \sim f_{32}=n_{1 r p q} k_{1 rp n} \exp \left(\mathrm{i} k_{1 rpq} n_{1 r p q} H\right) \text {, } \\ & f_{33}=l_{1 r s} k_{1 r s} \exp \left(\mathrm{i} k_{1 r s} n_{1 r s} H\right) \text {, } \\ & f_{34} \sim f_{37}=n_{ {tpn }} k_{t p n} \exp \left(-\mathrm{i} k_{t p n} n_{t p n} H\right), \\ & f_{43}=n_{1 r s} k_{1 r s} \exp \left(\mathrm{i} k_{1 r s} n_{1 r s} H\right) \text {, } \\ & f_{38}=-l_{2 t s} k_{2 t s} \exp \left(-\mathrm{i} k_{2 t s} n_{2 t s} H\right) \text {, } \\ & f_{39} \sim f_{3(12)}=-n_{2 r p n} k_{2 rpn} \exp \left(\mathrm{i}_{2 rpn} n_{2 r p n} H\right), \\ & f_{3(13)}=-l_{2 rs} k_{2 rs} \exp \left(\mathrm{i}_{2 rs} n_{2 r s} H\right) \text {, } \\ & f_{41} \sim f_{42}=-l_{1 rpq} k_{1 rpq} \exp \left(\mathrm{i} k_{1 rpq} n_{1 rpq} H\right), \\ & f_{44} \sim f_{47}=l_{t p n} k_{t p n} \exp \left(-\mathrm{i} k_{t p n} n_{t p n} H\right), \\ & f_{48}=n_{2 t s} k_{2 t s} \exp \left(-\mathrm{i} k_{2 ts} n_{2 t s} H\right) \text {, } \\ & f_{49} \sim f_{4(12)}=l_{2 r p n} k_{2 r p n} \exp \left(i k_{2 r p n} n_{2 r p n} H\right), \\ & f_{4(13)}=-n_{2 rs} k_{2 r s} \exp \left(\mathrm{i} k_{2 r s} n_{2 r s} H\right) \text {, } \\ & f_{51} \sim f_{52}=-\delta_{Tpq}^{\mathrm{e}} \exp \left(\mathrm{i} k_{1 r p q} n_{1 r p q} H\right), \\ & f_{54} \sim f_{57}=\delta_{T p n} \exp \left(-\mathrm{i} k_{t p n} n_{t p n} H\right), \\ & f_{59} \sim f_{5(12)}=\delta_{Tp n} \exp \left(\mathrm{i} k_{2 r p n} n_{2 r p n} H\right), \\ & f_{53}=0, \quad f_{58}=0, \quad f_{5(13)}=0 \text {, } \\ & f_{61} \sim f_{62}=K^{\mathrm{e}} n_{1 r p q} k_{1 r p q} \delta_{T p q}^{\mathrm{e}} \exp \left(\mathrm{i} k_{1 rp q} n_{1 r p q} H\right), \quad f_{63}=0, \\ & f_{64} \sim f_{67}=K n_{t p n} k_{t p n} \delta_{T p n} \exp \left(-\mathrm{i} k_{t p n} n_{t p n} H\right), \\ & f_{68}=0, f_{6(13)}=0 \text {, } \\ & f_{69} \sim f_{6(12)}=-K n_{2 r p n} k_{2 r p n} \delta_{T p n} \exp \left(i k_{2 r p n} n_{2 r p n} H\right), \\ & f_{71}=f_{72}=f_{73}=0 \text {, } \\ & f_{74} \sim f_{77}=n_{t p n} k_{t p n} \delta_{f p n} \exp \left(-\mathrm{i} k_{\mathrm{tpn}} n_{t p n} H\right), \end{aligned}$

$\begin{aligned} &\begin{aligned} & f_{7 8}=-l_{2 t s} k_{2 t s} \delta_{f s} \exp \left(-\mathrm{i} k_{2 ts} n_{2 t s} H\right), \\ & f_{79} \sim f_{7(12)}=-n_{2 r p n} k_{2 r p n} \delta_{f p n} \exp \left(\mathrm{i} k_{2 r p n} n_{2 r pn} H\right) \text {, } \\ & f_{7(13)}=-l_{2 r s} k_{2 r s} \delta_{f s} \exp \left(i k_{2 r s} n_{2 r s} H\right), \quad f_{81} \sim f_{83}=0 \text {, } \\ & f_{84} \sim f_{87}=n_{t p n} k_{t pn} \delta_{a p n} \exp \left(-\mathrm{i} k_{t p n} n_{t p n} H\right) \text {, } \\ & f_{88}=-l_{2 t s} k_{2 t s} \delta_{a s} \exp \left(-\mathrm{i} k_{2 t s} n_{2 t s} H\right) \text {, } \\ & f_{89} \sim f_{8(12)}=-n_{2 r p n} k_{2 r p n} \delta_{a p n} \exp \left(\mathrm{i} k_{2 r p n} n_{2 r p n} H\right) \text {, } \\ & f_{8(13)}=-l_{2 r s} k_{2 r s} \delta_{a s} \exp \left(i k_{2 r s} n_{2 r s} H\right), \quad f_{91} \sim f_{93}=0 \text {, } \\ & f_{94} \sim f_{97}=\left(\bar{\lambda}+2 \mu n_{t p n}^2+D_1 \delta_{f p n}+D_2 \delta_{a p n}\right) k_{ {tpn }}^2-D_3 \delta_{T p n} \text {, } \\ & f_{98}=-2 \mu l_{2 t s} n_{2t s} k_{2t s}^2 \text {, } \\ & f_{99} \sim f_{9(12)}=\left(\bar{\lambda}+2 \mu n_{2 r p n}^2+D_1 \delta_{f p n}+D_2 \delta_{a p n}\right) k_{2 r p n}^2-D_3 \delta_{Tp n} \text {, } \\ & f_{9(13)}=2 \mu l_{2 r s} n_{2 r s} k_{2 r s}^2, \quad f_{10(1)} \sim f_{10(3)}=0 \text {, } \end{aligned}\\ &\begin{aligned} & f_{10(4)} \sim f_{10(7)}=2 \mu l_{t p n} n_{t p n} k_{t p n}^2,\\ & f_{10(8)}=\mu\left(1-2 l_{2 ts}^2\right) k_{2 ts}^2, \quad f_{10(9)} \sim f_{10(12)}=-2 \mu l_{2r p n} n_{2 r p n} k_{2 r p n}^2 \text {, } \\ & f_{10(13)}=\mu\left(1-2 l_{2 r s}^2\right) k_{2 r s}^2, \quad f_{11(1)}=f_{11(2)}=f_{11(3)}=0 \text {, } \\ & f_{11(4)} \sim f_{14(7)}=\left(B_1+B_2 \delta_{f p n}+B_3 \delta_{a p n}\right) k_{t p n}^2-B_4 \delta_{T_{p n}} \text {, } \\ & f_{11(8)}=0 \text {, } \\ & f_{11(9)} \sim f_{11(12)}=\left(B_1+B_2 \delta_{f p n}+B_3 \delta_{\text {apn }}\right) k_{2 r p n}^2-B_4 \delta_{Tp n} \text {, } \\ & f_{11(13)}=0, \quad f_{12(1)}=f_{12(2)}=f_{12(3)}=0 \text {, } \\ & f_{12(4)} \sim f_{12(7)}=\left(B_5+B_6 \delta_{f p n}+B_7 \delta_{a p n}\right) k_{t p n}^2-B_8 \delta_{Tp n} \text {, } \\ & f_{12(8)}=0 \text {, } \\ & f_{12(9)} \sim f_{12(12)}=\left(B_5+B_6 \delta_{f p n}+B_7 \delta_{a p n}\right) k_{2 r p n}^2-B_8 \delta_{Tp n} \text {, } \\ & f_{12(13)}=0, \quad f_{13(1)}=f_{13(2)}=f_{13(3)}=0 \text {, } \\ & f_{13(4)} \sim f_{13(7)}=n_{t p n}^2 k_{tpn}^2 \delta_{Tp n} \text {, } \\ & f_{13(8)}=0, \quad f_{13(9)} \sim f_{13(12)}=n_{2 r p n}^2 k_{2 r p n}^2 \delta_{Tp n}, \quad f_{13(13)}=0 \text {, } \\ & g_1=\left(\lambda^{\mathrm{e}}+2 \mu^{\mathrm{e}} n_{i p 1}^2\right) k_{i p 1}^2 \exp \left(-\mathrm{i} k_{i p 1} n_{i p 1} H\right)+3 K_b^{\mathrm{e}} \beta_T^{\mathrm{e}} \delta_{T p 1}^{\mathrm{e}} . \\ & \exp \left(-\mathrm{i} k_{i p 1} n_{i p 1} H\right), \quad g_2=2 \mu^{\mathrm{e}} l_{i p 1} n_{i p 1} k_{i p 1}^2 \exp \left(-\mathrm{i} k_{i p 1} n_{i p 1} H\right) \text {, } \\ & g_3=k_{i p 1} n_{i p 1} \exp \left(-\mathrm{i} k_{i p 1} n_{i p 1} H\right) \text {, } \\ & g_4=l_{i p 1} k_{i p 1} \exp \left(-\mathrm{i} k_{i p 1} n_{i p 1} H\right), \quad g_5=\delta_{T p 1}^{\mathrm{e}} \exp \left(-\mathrm{i} k_{i p 1} n_{i p1} H\right) \text {, } \\ & g_6=K^{\mathrm{e}} n_{i p 1} k_{i p 1} \delta_{T p 1}^{\mathrm{e}} \exp \left(-\mathrm{i} k_{i p 1} n_{i p 1} H\right), \quad g_7 \sim g_{13}=0 \text { 。 } \end{aligned}\\ &\text { 式中：} q=1,2 ; n=1,2,3,4 \text { 。 } \end{aligned}$

图 1 试验所用砂粒和粉粒的颗粒级配曲线

Figure 1. Grain-size distribution curves of sand and silt

下载: 全尺寸图片幻灯片

图 2 混合土细观特征观测流程

Figure 2. Observation process of micro-characteristics of mixed soil

下载: 全尺寸图片幻灯片

图 3 砂-黏混合土动孔压试验结果

Figure 3. Dynamic pore water pressures of sand-clay mixtures

下载: 全尺寸图片幻灯片

图 4 砂-黏混合土细观特征图像

Figure 4. Microscopic feature images of sand-clay mixtures

下载: 全尺寸图片幻灯片

图 5 砂-黏混合土颗粒接触状态示意图（e_s=0.9800）

Figure 5. Contact states of sand-clay mixtures(e_s=0.9800)

下载: 全尺寸图片幻灯片

图 6 砂-黏混合土颗粒接触状态示意图（e=0.8823）

Figure 6. Contact states of sand-clay mixtures(e=0.8823)

下载: 全尺寸图片幻灯片

图 7 砂-粉-黏混合土动孔压试验结果（e_s恒定）

Figure 7. Dynamic pore water pressures of sand-silt-clay mixtures (constant e_s)

下载: 全尺寸图片幻灯片

图 8 砂-粉-黏混合土动孔压试验结果（e恒定）

Figure 8. Dynamic pore water pressures of sand-silt-clay mixtures (constant e)

下载: 全尺寸图片幻灯片

图 9 砂-粉-黏混合土细观特征图像^[13]

Figure 9. Microscopic feature images of sand-silt-clay mixtures^[13]

下载: 全尺寸图片幻灯片

图 10 砂-粉-黏混合土颗粒接触状态示意图

Figure 10. Contact states of sand-silt-clay mixtures

下载: 全尺寸图片幻灯片

图 11 砂-粉-黏混合土颗粒接触状态示意图（e_s=0.9800）

Figure 11. Contact states of sand-silt-clay mixtures(e_s=0.9800)

下载: 全尺寸图片幻灯片

图 12 砂-粉-黏混合土颗粒接触状态示意图（e=0.8823）

Figure 12. Contact states of sand-silt-clay mixtures(e=0.8823)

下载: 全尺寸图片幻灯片

图 13 三模型对混合土动孔压的拟合结果

Figure 13. Fitting results of Booker, Baziar and Mod-Baziar models for dynamic pore water pressure of mixtures

下载: 全尺寸图片幻灯片

图 14 三模型对各个试样动孔压的拟合度

Figure 14. R-squared of Booker, Baziar and Mod-Baziar models for dynamic pore water pressure of mixtures

下载: 全尺寸图片幻灯片

图 15 饱和砂-黏、砂-粉-黏混合土动孔压发展

Figure 15. Development of dynamic pore water pressure for saturated sand-clay and sand-silt-clay mixtures

下载: 全尺寸图片幻灯片

图 16 Mod-Baziar模型中各个参数取值范围

Figure 16. Value ranges of each parameter in Mod-Baziar model

下载: 全尺寸图片幻灯片

表 1 试验设计方案

Table 1 Design schemes of dynamic tests

试验编号	细粒含量FC/%	细粒中黏、粉占比/%		黏粉比CS	固结压力σ₃/kPa	初始砂骨架孔隙比e_s	初始孔隙比e
试验编号	细粒含量FC/%	黏粒	粉粒	黏粉比CS	固结压力σ₃/kPa	初始砂骨架孔隙比e_s	初始孔隙比e
D1	0	－	－	－	200	0.9800	0.9800
D2	5	20	80	0.25	200	0.9800	0.8823
D3		50	50	1	200	0.9800	0.8823
D4		80	20	4	200	0.9800	0.8823
D5		100	0	－	200	0.9800	0.8823
D6	10	20	80	0.25	200	0.9800	0.7967
D7		50	50	1	200	0.9800	0.7967
D8		80	20	4	200	0.9800	0.7967
D9		100	0	－	200	0.9800	0.7967
D10	0	－	－	－	200	0.8823	0.8823
D11	10	20	80	0.25	200	1.0704	0.8823
D12		50	50	1	200	1.0704	0.8823
D13		80	20	4	200	1.0704	0.8823
D14		100	0	－	200	1.0704	0.8823

下载: 导出CSV

表 2 Mod-Baziar模型拟合参数取值范围

Table 2 Fitting parameter ranges of Mod-Baziar model

土体	模型参数
土体	θ	β	γ	ω
纯砂	1.5430~2.0362	-0.0057~-0.2760	0.1152 ~0.5637	0.4136~1.9222
砂-黏	1.3227~2.3576	-0.1089~-0.4709	0.0132~1.4852	-2.0889~2.7362
砂-粉-黏	1.2490~2.4263	-0.0510~-0.4252	-0.5369~1.8854	-2.8496~3.8353

下载: 导出CSV

参考文献(13)

[1]	SASAKI Y, TOWHATA I, MIYAMOTO K, et al. Reconnaissance report on damage in and around river levees caused by the 2011 off the Pacific coast of Tohoku earthquake[J]. Soils and Foundations, 2012, 52(5): 1016-1032. doi: 10.1016/j.sandf.2012.11.018
[2]	SEED H B, MARTIN P P, LYSMER J. Pore-water pressure changes during soil liquefaction[J]. Journal of the Geotechnical Engineering Division, 1976, 102(4): 323-346. doi: 10.1061/AJGEB6.0000258
[3]	BOOKER J R, RAHMAN M S, SEED H B. Gadflea: a computer program for the analysis of pore pressure generation and dissipation during cyclic or earthquake loading[R]. Berkeley, CA: Earthquake Engineering Research Center, University of California, 1976.
[4]	BAZIAR M, SHAHNAZARI H, SHARAFI H. A laboratory study on the pore pressure generation model for Firouzkooh silty sands using hollow torsional test[J]. International Journal of Civil Engineering, 2011, 9: 126-134.
[5]	许成顺, 王冰, 杜修力, 等. 循环加载频率对砂土液化模式的影响试验研究[J]. 土木工程学报, 2021, 54(11): 109-118. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC202111011.htm XU Chengshun, WANG Bing, DU Xiuli, et al. Experimental study on effect of cyclic loading frequency on liquefaction mode of sand[J]. China Civil Engineering Journal, 2021, 54(11): 109-118. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC202111011.htm
[6]	董青, 周正华, 苏杰, 等. 基于对数动骨架考虑可逆孔压的有效应力本构研究[J]. 岩土工程学报, 2020, 42(12): 2322-2329. doi: 10.11779/CJGE202012020 DONG Qing, ZHOU Zhenghua, SU Jie, et al. Constitutive model for effective stress based on logarithmic skeleton curve considering reversible pore pressure[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(12): 2322-2329. (in Chinese) doi: 10.11779/CJGE202012020
[7]	王桂萱, 桑野二郎, 竹村次朗. 循环荷载下砂质混合土孔隙水压力特性研究[J]. 岩土工程学报, 2004, 26(4): 541-545. doi: 10.3321/j.issn:1000-4548.2004.04.023 WANG Guixuan, JiRO Kuwano, JIRO Takemura. Study on excess pore water pressures of sands mixed with clays under cyclic loading[J]. Chinese Journal of Geotechnical Engineering, 2004, 26(4): 541-545. (in Chinese) doi: 10.3321/j.issn:1000-4548.2004.04.023
[8]	吴琪, 王路阳, 刘启菲, 等. 基于剪切应变特征的饱和珊瑚砂超静孔压发展模型试验研究[J]. 岩土工程学报, 2023, 45(10): 2091-2099. doi: 10.11779/CJGE20220956 WU Qi, WANG Luyang, LIU Qifei, et al. Experimental study on development model of excess pore pressure for saturated coral sand based on shear strain characteristics[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(10): 2091-2099. (in Chinese) doi: 10.11779/CJGE20220956
[9]	张伏光, 聂卓琛, 陈孟飞, 等. 不排水循环荷载条件下胶结砂土宏微观力学性质离散元模拟研究[J]. 岩土工程学报, 2021, 43(3): 456-464. doi: 10.11779/CJGE202103008 ZHANG Fuguang, NIE Zhuochen, CHEN Mengfei, et al. DEM analysis of macro- and micro-mechanical behaviors of cemented sand subjected to undrained cyclic loading[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(3): 456-464. (in Chinese) doi: 10.11779/CJGE202103008
[10]	王志华, 何健, 高洪梅, 等. 基于触变流体理论的可液化土体振动孔压模型[J]. 岩土工程学报, 2018, 40(12): 2332-2340. doi: 10.11779/CJGE201812023 WANG Zhihua, HE Jian, GAO Hongmei, et al. Dynamic pore water pressure model for liquefiable soils based on theory of thixotropic fluid[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(12): 2332-2340. (in Chinese) doi: 10.11779/CJGE201812023
[11]	土工试验方法标准: GB/T 50123—2019[S]. 北京: 中国计划出版, 2019. Standard for Soil Test Method: GB/T 50123—2019[S]. Beijing: China Planning Press, 2019. (in Chinese)
[12]	刘洋, 吴顺川, 周健. 循环荷载下砂土变形的细观数值模拟Ⅱ: 密砂试验结果[J]. 岩土工程学报, 2007, 29(11): 1676-1682. doi: 10.3321/j.issn:1000-4548.2007.11.014 LIU Yang, WU Shunchuan, ZHOU Jian. Micro-numerical simulation of cyclic biaxial test Ⅱ: results of dense sand[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(11): 1676-1682. (in Chinese) doi: 10.3321/j.issn:1000-4548.2007.11.014
[13]	LI T, TANG X W. Influences of low fines content and fines mixing ratio on the undrained static shear strength of sand-silt-clay mixtures[J]. European Journal of Environmental and Civil Engineering, 2022, 26(9): 3706-3728. doi: 10.1080/19648189.2020.1813206