Loading [MathJax]/jax/output/SVG/jax.js
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

热效应作用下P波入射非饱和土自由场地地震响应研究

马强, 杨奕琪, 周凤玺, 邵生俊

马强, 杨奕琪, 周凤玺, 邵生俊. 热效应作用下P波入射非饱和土自由场地地震响应研究[J]. 岩土工程学报, 2025, 47(3): 569-579. DOI: 10.11779/CJGE20230841
引用本文: 马强, 杨奕琪, 周凤玺, 邵生俊. 热效应作用下P波入射非饱和土自由场地地震响应研究[J]. 岩土工程学报, 2025, 47(3): 569-579. DOI: 10.11779/CJGE20230841
MA Qiang, YANG Yiqi, ZHOU Fengxi, SHAO Shengjun. Seismic response of free-field earthquakes in unsaturated soils with P-wave incidence under thermal effects[J]. Chinese Journal of Geotechnical Engineering, 2025, 47(3): 569-579. DOI: 10.11779/CJGE20230841
Citation: MA Qiang, YANG Yiqi, ZHOU Fengxi, SHAO Shengjun. Seismic response of free-field earthquakes in unsaturated soils with P-wave incidence under thermal effects[J]. Chinese Journal of Geotechnical Engineering, 2025, 47(3): 569-579. DOI: 10.11779/CJGE20230841

热效应作用下P波入射非饱和土自由场地地震响应研究  English Version

基金项目: 

国家自然科学基金项目 52168053

青海省自然科学基金项目 2024-ZJ-922

详细信息
    作者简介:

    马强(1990—),男,博士,副教授,主要从事土动力学方面的研究工作。E-mail: maqiang0104@163.com

    通讯作者:

    周凤玺, E-mail: geolut@163.com

  • 中图分类号: TU43

Seismic response of free-field earthquakes in unsaturated soils with P-wave incidence under thermal effects

  • 摘要: 基于单相热弹性介质和非饱和多孔热弹性介质中波的传播理论,建立了平面P波入射下非饱和土自由场地模型,采用亥姆霍兹矢量分解原理,对热效应作用下非饱和土自由场地中的波场进行分析,获得了热效应作用下平面P波入射非饱和土自由场地地震响应的解析解答。通过数值计算,分析了热传导系数、介质温度、热膨胀系数等热物性参数对非饱和土自由场地地震响应所产生的影响规律。研究结果表明:考虑热效应和不考虑热效应两种理论模型下所得到的地表位移放大系数有着明显差异;水平位移放大系数随着热传导系数的增大而逐渐增大,但竖向位移放大系数几乎不发生变化;随着热膨胀系数的增大,水平位移放大系数逐渐增大,竖向位移放大系数逐渐减小;随着介质温度的增大,水平位移放大系数逐渐减小,竖向位移放大系数逐渐增大;此外,随着饱和度的增大,水平位移放大系数逐渐减小,竖向位移放大系数逐渐增大,说明土体中气相的存在对地表位移放大系数有着显著的影响。
    Abstract: Based on the theories of wave propagation in single-phase and unsaturated porous thermoelastic media, a free-field model for unsaturated soils under plane P-wave incidence is established, and the Helmholtz vector decomposition principle is employed to analyze the wave fields in the free field of unsaturated soils under the action of thermal effects, and analytical solutions for the seismic response of the free field of unsaturated soils under the action of thermal effects with plane P-wave incidence are obtained. The influences of thermal physical parameters such as thermal conductivity, medium temperature, and thermal expansion coefficient on the seismic response of a free field of unsaturated soils is analyzed through numerical calculations. The results show that the displacement magnification factors of the ground obtained under the two theoretical models considering and not considering thermal effects have apparent differences. The horizontal displacement magnification factor gradually increases with the increase of the thermal conductivity, but the vertical one is almost unchanged. The horizontal displacement magnification factor gradually increases with the increase of the coefficient of thermal expansion, while the vertical one decreases gradually. With the increase of the medium temperature, the horizontal displacement magnification factor decreases, and the vertical one gradually increases. In addition, with the increase of the saturation degree, the horizontal displacement magnification factor gradually decreases, and the vertical one gradually increases, which indicates that the existence of gas phase in the soils has a significant influence on the surface displacement magnification factor.
  • 场地土性质的改变会引起地震波的放大或者衰减,因此在进行地震响应的研究时要选择符合实际的场地模型。然而,以往对场地地震响应问题的研究,是将场地土简化为单相弹性地基或饱和土地基,且没有考虑热效应的影响。但是,在自然界中,非饱和土是土体更为常见的一种状态,且随着温室效应的加重,温度的变化越来越显著。此外,随着各国工业化和现代化进程的加快,以及热力学工程,化工工程和核废料管理等学科领域的大力发展,温度变化等热效应对多孔介质中波的传播特性的影响越来越受到众多学者的关注。由于热和力之间的相互耦合效应影响着波的形成和传播,并且多相介质之间的温度和应力-应变的耦合比经典情况下更复杂,这些因素叠加起来使得波的传播也更为复杂。因此,考虑热效应影响,建立一个能够反映实际场地土特性的非饱和土场地模型,并研究其地震地面运动具有重要意义。

    自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]采用了考虑孔隙流体与孔隙材料固相相互作用的非等温动力学模型,利用混合物理论分析了均匀平面P1波在非饱和多孔热弹性介质边界处的反射特性。Wei等[19]研究了各向同性、均匀的热弹性介质与多孔热弹性介质在平面界面上,由斜入射纵波所引起的反射和透射现象。Hou等[20]研究了在任意入射角下,两种不相似的饱和多孔热弹性介质界面上,非均匀平面波的反射和透射特性。Liu等[21]在热流-固体耦合的非饱和土波动控制方程的基础上,进一步研究了热弹性固体介质与非饱和多孔热弹性介质之间的平面界面上,由均匀平面P波入射所产生的反射和透射现象。

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

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

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

    μe2us+(λe+μe)(us)3KbβsT=ρe¨us, (1a)
    3KbβsT0(˙us+τqe¨us)+ρecse(˙T+τqe¨T)=Ke2(T+τθe˙T) (1b)

    式中:Kb=λe+2μe/3为体积模量;λeμe为单相介质的拉梅常量;βs为热膨胀系数;ρe为密度;cse为固相比热;τqeτθe分别为热通量相位延迟时间和温度梯度相位延迟时间;TT0分别为开尔文介质温度和初始温度;Ke为单相介质的热传导系数;us为固相位移矢量。

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

    us=ψs+×Hs (2)

    式中:ΨsHs分别为固相的标量和矢量势函数。

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

    ψs=Asexp[i(kpxωt)] (3a)
    Hs=Bsexp[i(ksxωt)] (3b)
    T=ATexp[i(kpxωt)] (3c)

    式中:AB分别为P波(包含T波)和S波的波幅系数;kpks分别为P波(包括T波)和S波的波数;ω表示频率。

    将式(2)代入式(1a),(1b),并结合式(3a)~(3c),可得热弹性介质中体波的特征方程:

    |k11k12k21k22|=0 (4a)
    |ρeω2μek2s|=0 (4b)

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

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

    μ2us+(ˉλ+μ)(us)+D1(uf)+D2(ua)+D3T=ρ¨us+ρ1¨uf+ρg¨ua, (5a)
    B1(us)+B2(uf)+B3(ua)+B4T=ρ1¨us+ϑi¨uf+vl˙uf, (5b)
    B5(us)+B6(uf)+B7(ua)+B4T=ρg¨us+ϑg¨ua+vg˙ua (5c)
    C1(˙us+τq¨us)+C2(˙uf+τq¨uf)+C3(˙ua+τq¨ua)+C4(˙T+τq¨T)=K2(T+τθ˙T) (5d)

    式中:λ=λ+D0λμ表示Lame常数;uf=nSl(ul-us)为土体孔隙中液相相对于固体骨架的相对位移;ua=nSg(ug-us)为孔隙中气相相对于固体骨架的相对位移,其中usulug为固相、液相和气相的绝对位移向量,n为土体的孔隙率;Sl为孔隙内液态水的饱和度、Sg为孔隙内气体的饱和度;ρ为非饱和多孔介质的总密度,并且满足ρ=ˉρs+ˉρl+ˉρg,其中ˉρsˉρlˉρg表示固-液-气三相的表观密度,且3个表观密度可以分别表示为ˉρs=(1-n)ρsˉρl=nSlρlˉρg= nSgρg,其中ρsρlρg分别为固相、液相和气相的密度;τq为热通量的相位延迟时间;K为热传导系数;τθ为温度梯度的相位延迟时间。此外,其他符号的表达式见文献[14]。

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

    uα=ψα+×Hα (6)

    式中:ψαHα为各相介质的势函数(α=s,f,a),s表示固相、f表示液相、a表示气相。

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

    ψα=Aαexp[i(kpxωt)] (7a)
    Hα=Bαexp[i(ksxωt)] (7b)
    T=ATexp[i(kpxωt)] (7c)

    式中:AαBαα=s,f,a)分别为各相中P波(包括T波)和S波对应的势函数幅值;kpks分别为压缩波(包括热波)和剪切波的复波数;ω为角频率。

    将式(6)代入式(5a)~(5d),并结合式(7a),(7b),(7c),可得到热弹性体波的特征方程如下:

    |b11b12b13b14b21b22b23b24b31b32b33b34b41b42b43b44|=0 (8a)
    |c11c12c13c21c22c23c31c32c33|=0 (8b)

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

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

    图  1  非饱和土自由场地简化模型
    Figure  1.  Simplified model for free site of unsaturated soils

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

    (1)单相热弹性介质中体波的势函数

    P 波: ψe=Asip1exp[ikip1(lip1xnip1zcip1t)]+2q=1As1rpqqexp[ik1rpq(l1rpqx+n1rpqqzc1rpqt)], (9a)
    TTe=Asip1δeTp1exp[ikip1 (lip1xniplzcip1t)]+2q=1As1rpqδeTpqexp[ik1rpq(l1rpqx+n1rpqzc1rpqt)], (9b)
    SHe=B1rsexp[ik1rs(l1rsx+n1rszc1rst)] (9c)

    (2)非饱和土层中体波的势函数

    P 波: ψαu=4n=1Aαtpnexp[iktpn(ltpnxntpnzctpnt)]+4n=1Aα2rpnexp[ik2rpn(l2rpnx+n2rpnzc2rpnt)], (10a)
    T 波: Tu=4n=1AstpnδTpnexp[iktpn(ltpnxntpnzctpnt)]+4n=1As2rpnδTpnexp[ik2rpn(l2rpnx+n2rpnzc2rpnt)], (10b)
    S 波: Hαu=Bα2tsexp[ik2ts(l2tsxn2tszc2tst)]+Bα2rsexp[ik2rs(l2rsx+n2rszc2rst)] (10c)

    式中:α取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,分别代表反射(透射)P1波、反射(透射)P2波、反射(透射)P3波和反射(透射)T波。kip1cip1分别表示入射P波的波数和波速;k1rpqk1rsc1rpβc1rs分别表示半无限基岩层中2种反射压缩波(P波、T波)和1种反射剪切波(S波)的波数以及对应的波速;ktpnk2ts分别表示非饱和土层中4种透射波(包括P波和T波)和1种透射S波的波数;ctpnc2ts分别表示非饱和土层中4种透射纵波(包括热波)和1种透射S波的波速;k2rpnk2rs分别表示非饱和土层中4种反射波(包括P波和T波)和1种反射S波的波数;c2rpnc2rs分别表示非饱和土层中4种反射波(包括P波和T波)和1种反射S波的波速;ln分别表示xz方向的方向矢量;i=1

    由Snell定理可得

    lip1kip1=l1rp1k1p1=l1rp2k1rp2=l1rsk1rs=ltp1ktp1=ltp2ktp2=ltp3ktp3=ltp4ktp4=l2tsk2ts (11)

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

    δeTpn=ATAs=k11k12(n=1,2 (12)

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

    δfpn=AfAs=d11d15d12d16d13d16d14d15 (n=1,2,3,4)  (13a)
    δapn=AaAs=d21d25d22d26d23d26d24d25 (n=1,2,3,4)  (13b)
    δTpn=ATAs=d31d35d32d36d33d36d34d35 (n=1,2,3,4)  (13c)
    δfs=BfBs=c12c22,δas=BaBs=c21c12c11c22c13c22  (13d)

    式中:元素d11~d36的表达式详见文献[21]。

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

    交界面上应力连续:

    σzz=σezz,σxz=σexz  (14a)

    交界面上位移连续:

    usz=uez,usx=uex  (14b)

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

    T=Te,KTz=KeTez  (14c)

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

    ufz=uaz=0  (14d)

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

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

    σzz=σxz=pw=pa=Tz=0 (15)

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

    uex=ψexHez (16a)
    uez=ψez+Hex (16b)
    σexz=μe(22ψexz+2Hex22Hez2) (16c)
    σezz=λe(2ψex2+2ψez2)+2μe(2ψez2+2Hexz)3KbβTT  (16d)

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

    uαx=ψαuxHαuz (17a)
    uαz=ψauz+Hαux (17b)
    σzz=ˉλ(2ψsux2+2ψsuz2)+2μ(2ψsuz2+2Hsuxz)+D1(2ψfux2+2ψfuz2)+D2(2ψaux2+2ψauz2)+D3Tu, (17c)
    σxz=μ(22ψsuxz+2Hsux22Hsuz2) (17d)
    pg=B5(2ψsux2+2ψsuz2)B6(2ψfux2+2ψfuz2)B7(2ψaux2+2ψauz2)B8Tu (17e)
    pw=B1(2ψsux2+2ψsuz2)B2(2ψfux2+2ψfuz2)B3(2ψaux2+2ψauz2)B4Tu (17f)

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

    FN=AipG (18)

    式中:

    N=[As1rp1,As1rp2,Bs1rs,Astp1,Astp2,Astp3,Astp4,Bs2ts,As2rp1,As2rp2,As2rp3,As2rp4,Bs2rs]T矩阵F和矩阵G中的元素在附录A中给出。

    波幅系数确定后,场地各点的应力和位移均可确定。将式(10a)~(10c)代入到水平方向位移ux和竖直方向位移uz的表达式中可得

    ux=Astp1iktp1ltp1+Astp2iktp2ltp2+Astp3ktp3ltp3+Astp4iktp4ltp4+As2rp1ik2rp1l2rp1+As2rp2ik2rp2l2rp2+As2r3ik2rp3l2rp3+As2rp4ik2rp4l2rp4+Bs2tsik2tsn2tsBs2rsik2rsn2rs, (19a)
    uz=Astp1iktp1ntp1Astp2iktp2ntp2Astp3iktp3ntp3Astp4ktp4ntp4+As2rp1ik2rp1n2rp1+As2rp2ik2rp2n2rp2+As2rp3k2rp3n2rp3+As2rp4ik2rp4n2rp4+Bs2tsik2tsl2ts+Bs2rsik2rsl2rs (19b)

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

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

    图  2  本文数值解与文献解的比较
    Figure  2.  Comparison between numerical solutions in this study and literature solutions

    随后通过数值计算,研究热效应作用下P波入射非饱和土自由场地的地震响应问题,具体分析了非饱和土层中的介质温度等热物性参数对场地地震响应所产生的影响,在数值计算中所取的热物性参数[21]表 12所示。本文分析中采用无量纲频率ω0,其中ω0=ω/ω1ω0=1,ω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 Ke/(J·s-1·m-1·K-1)
    量值 12.0×109 8.0×109 2700 4.0×10-4 2.5
    参数 Te/K τqe/s τθe/s T0/K cse/(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×109 2.8×109 2650 1.0×10-4 0.6
    参数 T/K τq/s τθ/s T/K cs/(J·kg-1·K-1)
    量值 293.2 2.0×10-7 1.5×10-7 300 1000
    参数 cl/(J·kg-1·K-1) cg/(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
    参数 Sl Ssat Sres n k/m2
    量值 0.6 1.0 0.05 0.4 1.0×10-10
    参数 ρl/(kg·m-3) Ks/kPa Pg/kPa m d
    量值 1000 3.5×107 101.3 0.5 2
    下载: 导出CSV 
    | 显示表格

    热传导系数表征了岩石的导热性,不仅是研究地壳热状态的重要参数,也被广泛应用于地震勘探和采矿等领域。为了研究不同热传导系数下和等温模型下地表位移放大系数随入射角的变化曲线,忽略热效应的影响,令波动控制方程中温度的变化量以及温度随时间的导数为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给出了不同饱和度下地表位移放大系数随入射角的变化曲线。首先从图 4可知,随着入射角的增大,水平位移放大系数和竖向位移放大系数均呈现先增加后减小的趋势,分别在θ=45°和θ=5°时达到峰值,另外当入射角θ=0°时,地表水平位移放大系数为0,而当入射角θ=90°时,地表水平位移和竖向位移放大系数均为0。其原因是:平面P波的能量主要沿着传播的方向传递,当平面P波沿着基岩和非饱和土层的分界面水平入射时(θ=90°),由于没有垂直于分界面的分量,非饱和土层中的反射和透射波消失[21],故非饱和土表面不会产生位移。其次,从图 4可以发现随着饱和度的增大,地表水平位移放大系数逐渐减小,竖向位移放大系数逐渐增大。其原因是:P2波是由液相的作用产生的,P3波是由孔隙液体和气体的压力差产生的,当土体接近完全饱和时,土体孔隙中液相的作用力达到最大,P2波的波速达到最大值,而随着气体的不断消散,P3波的波速减小到最小值。此外,入射P波的能量主要分配给P1波和S波,而多孔介质孔隙中存在的流体相会阻碍P1波和S波的传播,液相的阻碍能力强于气相[24]。因此这些因素的共同作用使得不同饱和度下的土体,由于气体含量的不同对场地地震地面运动的影响也有所不同。此外,本文所得的这一变化趋势与文献[25]所得到的等温条件下饱和度变化对地表位移放大系数的影响规律一致。由此可见,图 4所表现出的地表位移放大系数随饱和度变化的规律,再一次验证了本文计算结果的有效性,又说明非饱和土中气相的存在对地表位移放大系数的影响不可忽视。把场地土视为饱和土地基从而没有考虑土体饱和度对场地地震地面运动的影响,这与实际现象是有出入的。

    图  4  不同饱和度下位移放大系数随入射角度变化的曲线
    Figure  4.  Variation curves of displacement amplification coefficient with incident angle under different saturation

    热膨胀系数主要涉及描述固相应力应变关系的物理方程和描述热传导的热传导方程,为了分析热膨胀系数变化对地表位移放大系数的影响规律,图 5给出了不同热膨胀系数下位移放大系数随入射角的变化曲线。首先从图 5可知,地表水平位移放大系数和竖向位移放大系数随着入射角度的增大,均呈现先增加后减小的趋势。其次,从图 5可以看出,随着热膨胀系数的增大,水平位移放大系数逐渐增大,竖向位移放大系数也逐渐减小,且变化的幅值较大。Liu等[21]研究发现反射P波、反射S波以及透射P1波和透射S波所占据的能量较多,且当热膨胀系数βT增加1×10-4 K−1时,透射P1波的振幅比和能量比分别增加约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

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

    图  6  不同温度下位移放大系数随入射角度变化的曲线
    Figure  6.  Variation curve of displacement amplification coefficient with incident angle at different temperatures

    相对于等温条件下的多孔介质理论,在多孔热弹性介质理论下将会产生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波的小好几个数量级。因此,热通量相位延迟时间的变化对位移放大系数的影响较小。

    为了分析土层深度变化对土层应力的影响,图 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

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

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

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

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

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

    f11f12=[(λe+2μen21rpn)k21rpn3kebβeTδeTpq)]exp(ik1rpqn1rpqH),f13=2μel1rsn1rsk21rsexp(ik1rsn1rsH),
    f14f17=[(ˉλ+2μn2tpn+D1δfpn+D2δapn)k2tpnD3δTpn].exp(iktpnntpnH),f18=2μl2tsn2tsk22tsexp(ik2tsn2tsH),
    f19f1(12)=[(ˉλ+2μn22rpn+D1δ2rpn+D2δ2rpn)k22rpnD3δTpn].exp(ik2rpnn2rpnH),f1(13)=2μl2rsn2rsk22rsexp(ik2rsn2rsH)f21f22=2μel1rpqn1rpqk21rpqexp(ik1rpqn1rpqH),f23=μe(l21rsn21rs)k21rsexp(i1rsn1rsH),f24f27=2μltpnntpnk2tpnexp(iktpnntpnH),f28=μ(n22tsl22ts)k22tsexp(ik2tsn2tsH)f29f2(12)=2μl2rpnn2rpnk22rpnexp(ik2rpnn2rpnH)f2(13)=μ(n22rsl22rs)k22rsexp(ik2rsn2rsH)f31f32=n1rpqk1rpnexp(ik1rpqn1rpqH)f33=l1rsk1rsexp(ik1rsn1rsH)f34f37=ntpnktpnexp(iktpnntpnH),f43=n1rsk1rsexp(ik1rsn1rsH)f38=l2tsk2tsexp(ik2tsn2tsH)f39f3(12)=n2rpnk2rpnexp(i2rpnn2rpnH),f3(13)=l2rsk2rsexp(i2rsn2rsH)f41f42=l1rpqk1rpqexp(ik1rpqn1rpqH),f44f47=ltpnktpnexp(iktpnntpnH),f48=n2tsk2tsexp(ik2tsn2tsH)f49f4(12)=l2rpnk2rpnexp(ik2rpnn2rpnH),f4(13)=n2rsk2rsexp(ik2rsn2rsH)f51f52=δeTpqexp(ik1rpqn1rpqH),f54f57=δTpnexp(iktpnntpnH),f59f5(12)=δTpnexp(ik2rpnn2rpnH),f53=0,f58=0,f5(13)=0f61f62=Ken1rpqk1rpqδeTpqexp(ik1rpqn1rpqH),f63=0,f64f67=KntpnktpnδTpnexp(iktpnntpnH),f68=0,f6(13)=0f69f6(12)=Kn2rpnk2rpnδTpnexp(ik2rpnn2rpnH),f71=f72=f73=0f74f77=ntpnktpnδfpnexp(iktpnntpnH),
    f78=l2tsk2tsδfsexp(ik2tsn2tsH),f79f7(12)=n2rpnk2rpnδfpnexp(ik2rpnn2rpnH)f7(13)=l2rsk2rsδfsexp(ik2rsn2rsH),f81f83=0f84f87=ntpnktpnδapnexp(iktpnntpnH)f88=l2tsk2tsδasexp(ik2tsn2tsH)f89f8(12)=n2rpnk2rpnδapnexp(ik2rpnn2rpnH)f8(13)=l2rsk2rsδasexp(ik2rsn2rsH),f91f93=0f94f97=(ˉλ+2μn2tpn+D1δfpn+D2δapn)k2tpnD3δTpnf98=2μl2tsn2tsk22tsf99f9(12)=(ˉλ+2μn22rpn+D1δfpn+D2δapn)k22rpnD3δTpnf9(13)=2μl2rsn2rsk22rs,f10(1)f10(3)=0f10(4)f10(7)=2μltpnntpnk2tpn,f10(8)=μ(12l22ts)k22ts,f10(9)f10(12)=2μl2rpnn2rpnk22rpnf10(13)=μ(12l22rs)k22rs,f11(1)=f11(2)=f11(3)=0f11(4)f14(7)=(B1+B2δfpn+B3δapn)k2tpnB4δTpnf11(8)=0f11(9)f11(12)=(B1+B2δfpn+B3δapn )k22rpnB4δTpnf11(13)=0,f12(1)=f12(2)=f12(3)=0f12(4)f12(7)=(B5+B6δfpn+B7δapn)k2tpnB8δTpnf12(8)=0f12(9)f12(12)=(B5+B6δfpn+B7δapn)k22rpnB8δTpnf12(13)=0,f13(1)=f13(2)=f13(3)=0f13(4)f13(7)=n2tpnk2tpnδTpnf13(8)=0,f13(9)f13(12)=n22rpnk22rpnδTpn,f13(13)=0g1=(λe+2μen2ip1)k2ip1exp(ikip1nip1H)+3KebβeTδeTp1.exp(ikip1nip1H),g2=2μelip1nip1k2ip1exp(ikip1nip1H)g3=kip1nip1exp(ikip1nip1H)g4=lip1kip1exp(ikip1nip1H),g5=δeTp1exp(ikip1nip1H)g6=Kenip1kip1δeTp1exp(ikip1nip1H),g7g13=0 。  式中:q=1,2;n=1,2,3,4 。 
  • 图  1   非饱和土自由场地简化模型

    Figure  1.   Simplified model for free site of unsaturated soils

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    参数 λe/Pa μe/Pa ρe/(kg·m-3) βs/K-1 Ke/(J·s-1·m-1·K-1)
    量值 12.0×109 8.0×109 2700 4.0×10-4 2.5
    参数 Te/K τqe/s τθe/s T0/K cse/(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×109 2.8×109 2650 1.0×10-4 0.6
    参数 T/K τq/s τθ/s T/K cs/(J·kg-1·K-1)
    量值 293.2 2.0×10-7 1.5×10-7 300 1000
    参数 cl/(J·kg-1·K-1) cg/(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
    参数 Sl Ssat Sres n k/m2
    量值 0.6 1.0 0.05 0.4 1.0×10-10
    参数 ρl/(kg·m-3) Ks/kPa Pg/kPa m d
    量值 1000 3.5×107 101.3 0.5 2
    下载: 导出CSV
  • [1]

    BIOT M A. General solutions of the equations of elasticity and consolidation for a porous material[J]. Journal of Applied Mechanics, 1956, 23(1): 91-96. doi: 10.1115/1.4011213

    [2]

    LIU Z X, LIU J Q, MENG S B, et al. Diffraction of elastic waves by a fluid-filled crack in a fluid-saturated poroelastic half-space[J]. Geophysical Journal International, 2021, 225: 1530-1553. doi: 10.1093/gji/ggab043

    [3] 范凯祥, 申玉生, 闻毓民, 等. 平面Rayleigh波入射下饱和土中浅埋隧道复合式衬砌的动力响应[J]. 岩土工程学报, 2022, 44(3): 444-455. doi: 10.11779/CJGE202203006

    FAN Kaixiang, SHEN Yusheng, WEN Yumin, et al. Dynamic response of composite linings of shallowly buried tunnels in saturated soils subjected to incidence of plane Rayleigh waves[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(3): 444-455. (in Chinese) doi: 10.11779/CJGE202203006

    [4] 禹海涛, 王治坤, 刘中宪. SV波入射下均匀饱和地层渗透系数对深埋隧道的影响机制[J]. 岩土工程学报, 2022, 44(2): 201-211. doi: 10.11779/CJGE202202001

    YU Haitao, WANG Zhikun, LIU Zhongxian. Influence mechanism of permeability coefficient in homogeneously saturated strata on responses of deep tunnels under incidence of SV waves[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(2): 201-211. (in Chinese) doi: 10.11779/CJGE202202001

    [5] 巴振宁, 梁建文, 梅雄一. 斜入射平面SH波在层状饱和半空间中沉积谷地周围的三维散射[J]. 岩土工程学报, 2013, 35(3): 476-486. http://cge.nhri.cn/article/id/15001

    BA Zhenning, LIANG Jianwen, MEI Xiongyi. Three-dimensional scattering of obliquely incident plane SH waves in alluvial valley embedded in layered saturated half-space[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(3): 476-486. (in Chinese) http://cge.nhri.cn/article/id/15001

    [6] 邹新军, 杨紫健, 吴文兵. 非饱和土地基中端承桩对SH波的水平地震响应[J]. 岩土工程学报, 2024, 46(1): 72-80. doi: 10.11779/CJGE20221125

    ZOU Xinjun, YANG Zijian, WU Wenbing. Horizontal seismic response of end-bearing piles in unsaturated soil foundation under SH waves[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(1): 72-80. (in Chinese) doi: 10.11779/CJGE20221125

    [7]

    WEI C F, MURALEETHARAN K K. A continuum theory of porous media saturated by multiple immiscible fluids: Ⅰ linear poroelasticity[J]. International Journal of Engineering Science, 2002, 40(16): 1807-1833. doi: 10.1016/S0020-7225(02)00068-X

    [8]

    CHEN W Y, XIA T D, HU W T. A mixture theory analysis for the surface-wave propagation in an unsaturated porous medium[J]. International Journal of Solids and Structures, 2011, 48(16/17): 2402-2412.

    [9]

    CHEN W Y, XIA T D, CHEN W, et al. Propagation of plane P-waves at interface between elastic solid and unsaturated poroelastic medium[J]. Applied Mathematics and Mechanics, 2012, 33(7): 829-844.

    [10]

    TOMAR S K, ARORA A. Reflection and transmission of elastic waves at an elastic/porous solid saturated by two immiscible fluids[J]. International Journal of Solids and Structures, 2006, 43(7/8): 1991-2013.

    [11]

    LORD H W, SHULMAN Y. A generalized dynamical theory of thermoelasticity[J]. Journal of the Mechanics and Physics of Solids, 1967, 15(5): 299-309.

    [12]

    PECKER C, DERESIEWICZ H. Thermal effects on wave propagation in liquid-filled porous media[J]. Acta Mechanica, 1973, 16(1/2).

    [13]

    YOUSSEF H M. Theory of generalized porothermoelasticity[J]. International Journal of Rock Mechanics and Mining Sciences, 2007, 44(2): 222-227.

    [14] 柳鸿博, 周凤玺, 岳国栋, 等. 非饱和土中热弹性波的传播特性分析[J]. 岩土力学, 2020, 41(5): 1613-1624.

    LIU Hongbo, ZHOU Fengxi, YUE Guodong, et al. Analysis of propagation characteristics of thermoelastic waves in unsaturated soils[J]. Rock and Soil Mechanics, 2020, 41(5): 1613-1624. (in Chinese)

    [15] 郑荣跃, 刘干斌, 邓岳保, 等. SV波在饱和多孔热弹性介质平面界面上的反射[J]. 岩土工程学报, 2013, 35(增刊2): 839-843. http://cge.nhri.cn/article/id/15502

    ZHENG Rongyue, LIU Qianbin, DENG Yuebao, et al. Reflection of SV waves on plane interface of saturated porous thermoelastic media[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(S2): 839-843. (in Chinese) http://cge.nhri.cn/article/id/15502

    [16] 柳鸿博, 周凤玺, 郝磊超. 平面S波在饱和多孔热弹性介质边界的反射问题研究[J]. 地震工程学报, 2021, 43(1): 105-112.

    LIU Hongbo, ZHOU Fengxi, HAO Leichao. Reflection of plane S waves at the boundary of saturated porous thermo-elastic media[J]. China Earthquake Engineering Journal, 2021, 43(1): 105-112. (in Chinese)

    [17]

    ZHOU F X, ZHANG R L, LIU H B, et al. Reflection characteristics of plane-S-wave at the free boundary of unsaturated porothermoelastic media[J]. Journal of Thermal Stresses, 2020, 43(5): 579-593.

    [18]

    LIU H B, DAI G L, ZHOU F X, et al. A mixture theory analysis for reflection phenomenon of homogeneous plane-P1-wave at the boundary of unsaturated porothermoelastic media[J]. Geophysical Journal International, 2021, 228(2): 1237-1259.

    [19]

    WEI W, ZHENG R Y, LIU G B, et al. Reflection and refraction of P wave at the interface between thermoelastic and porous thermoelastic medium[J]. Transport in Porous Media, 2016, 113(1): 1-27.

    [20]

    HOU W T, FU L Y, CARCIONE J M, et al. Reflection and transmission of inhomogeneous plane waves in thermoelastic media[J]. Frontiers in Earth Science, 2022, 10: 850331.

    [21]

    LIU H B, DAI G L, ZHOU F X, et al. Propagation behavior of homogeneous plane-P1-wave at the interface between a thermoelastic solid medium and an unsaturated porothermoelastic medium[J]. European Physical Journal Plus, 2021, 136(11): 1163.

    [22] 李伟华, 郑洁. 饱和度对平面P波入射下自由场地地震反应的影响分析[J]. 岩土工程学报, 2017, 39(3): 427-435. doi: 10.11779/CJGE201703005

    LI Weihua, ZHENG Jie. Effects of saturation on free-field responses of site due to plane P-wave incidence[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(3): 427-435. (in Chinese) doi: 10.11779/CJGE201703005

    [23]

    LI W H, ZHENG J, TRIFUNAC M D. Saturation effects on ground motion of unsaturated soil layer-bedrock system excited by plane P and SV waves[J]. Soil Dynamics and Earthquake Engineering, 2018, 110: 159-172.

    [24]

    LIU H B, JIANG M, ZHOU F, et al. Attenuation characteristics of thermoelastic waves in unsaturated soil[J]. Arabian Journal of Geosciences, 2021, 14(18): 1-11.

    [25] 郑洁. 非饱和土自由场地地震地面运动的解析分析[D]. 北京: 北京交通大学, 2015.

    ZHENG Jie. Analytical Analysis of Earthquake Ground Motion in Unsaturated Soil Free Field[D]. Beijing: Beijing Jiaotong University, 2015. (in Chinese)

图(8)  /  表(2)
计量
  • 文章访问数:  311
  • HTML全文浏览量:  22
  • PDF下载量:  36
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-09-03
  • 网络出版日期:  2024-08-20
  • 刊出日期:  2025-02-28

目录

/

返回文章
返回