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

相对密实度对沉贯中吸力基础桶壁–砂土界面力学特性的影响

李大勇, 侯新宇, 张雨坤, 高玉峰

李大勇, 侯新宇, 张雨坤, 高玉峰. 相对密实度对沉贯中吸力基础桶壁–砂土界面力学特性的影响[J]. 岩土工程学报, 2022, 44(9): 1598-1607. DOI: 10.11779/CJGE202209004
引用本文: 李大勇, 侯新宇, 张雨坤, 高玉峰. 相对密实度对沉贯中吸力基础桶壁–砂土界面力学特性的影响[J]. 岩土工程学报, 2022, 44(9): 1598-1607. DOI: 10.11779/CJGE202209004
LI Da-yong, HOU Xin-yu, ZHANG Yu-kun, GAO Yu-feng. Effects of relative densities on mechanical characteristics of interface between sand and suction caisson during penetration[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(9): 1598-1607. DOI: 10.11779/CJGE202209004
Citation: LI Da-yong, HOU Xin-yu, ZHANG Yu-kun, GAO Yu-feng. Effects of relative densities on mechanical characteristics of interface between sand and suction caisson during penetration[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(9): 1598-1607. DOI: 10.11779/CJGE202209004

相对密实度对沉贯中吸力基础桶壁–砂土界面力学特性的影响  English Version

基金项目: 

国家自然科学基金重点项目 51879044

国家自然科学基金重点项目 51808325

国家自然科学基金重点项目 51639002

山东科技大学科研创新团队项目 2015TDJH104

详细信息
    作者简介:

    李大勇(1971—),男,教授,从事海洋岩土工程和土力学理论等方面的研究。E-mail: ldy@upc.edu.cn

    通讯作者:

    张雨坤, E-mail: philc007@163.com

  • 中图分类号: TU470

Effects of relative densities on mechanical characteristics of interface between sand and suction caisson during penetration

  • 摘要: 吸力基础是一种系泊海洋浮动平台及固定海上风电塔架的基础形式。吸力基础在砂土中吸力贯入时,沉贯阻力主要由内外壁摩阻力组成。吸力基础沉贯过程中基础–土体界面力学特性决定侧壁摩阻力的大小。开展界面剪切试验,研究了相对密实度及沉贯深度对吸力基础–砂土界面沉贯过程中相互作用的影响。结果表明:吸力基础–砂土界面剪应力随剪切位移增加,呈现出软化特征;随砂土密实程度的增加,界面摩擦系数、界面峰值强度和相对应的剪切位移逐渐增大;但上述参数随沉贯深度的增加先减小后逐渐增大。研究工况条件下,界面摩擦系数量值处于0.2~0.34。根据试验结果,得到了不同沉贯深度处基础桶壁摩阻力,基础侧壁摩阻力–贯入深度之间关系可用指数函数表示;内外壁摩阻力随砂土密实程度增加逐渐增大,外壁摩阻力增长幅度较内壁摩阻力更加明显,外壁摩阻力最大值是相同沉贯深度处内壁摩阻力的2.1倍。界面剪切试验中,砂土发生明显挤出现象,砂土挤出主要发生在剪切初期,对应剪切位移较小,可用于揭示吸力基础沉贯过程中,由于桶壁贯入造成的砂土置换侧向挤出现象。界面剪切过程中砂土颗粒发生破碎,填充至砂土颗粒间孔隙中,导致土体试样高度下降,砂土呈现出“剪缩”现象。发现界面摩擦角与土体内摩擦角之比随土体密实程度的增加逐渐降低。研究内容有助于准确得到沉贯摩阻力和沉贯所需吸力,确保吸力基础沉贯至预定位置。
    Abstract: Suction caissons are popular for mooring oil and gas platforms and fixing wind turbines. The penetration resistance is mainly the friction between the inner and outer walls during suction penetration in sand. During installation of the suction caisson, the frictional forces between the caisson wall and the soil are dominated by the shear behaviors of the interface between the caisson wall and the soil. A series of interface shear tests with the GDS interface shearing instrument are conducted to investigate the influences of the initial relative densities and penetration depth on the interfacial shear characteristics of soil-suction caisson. The test results show that the shear stress–shear displacement relationship curve of interface under different relative densities exhibits strain softening, and interfacial shear strength and the corresponding shear displacement increase with the increasing relative densities. It is found that the friction coefficient of interface increases with the increase of the relative density, and firstly increases and then decreases with the increase of the penetration depth. Under the experimental programs of tests, the friction coefficient of interface ranges from 0.2 to 0.34, and the friction resistance at different depths can be calculated in terms of the test results, indicating the friction resistance between the inner and outer walls increases exponentially with the increase of the penetration depth. The friction resistance between the inner and outer walls increases with the increase of the relative density, and the increase rate of friction of the outer wall is greater than that of the inner wall. When the penetration completes, the friction of the outer wall is 2.12 times that of the inner wall. During shearing, sand particles are found to be squeezed out, and they are mainly squeezed out at the initial stage of the testing, i.e., the corresponding shear displacement is relatively small. This phenomenon indicates that during installation of the suction caisson, the sand particles are squeezed out because of the replacement between the sand volume and the volume of suction caisson wall. The sand particles broken in the process of testing causes the sand volume to decrease. Besides, it is found that the ratio of the interfacial friction angle to the internal friction angle decreases with the increase of the relative density. To obtain the accurate installation resistance and the required suction to install the suction caisson, it is important to investigate the shear behaviors of the interface between the suction wall and the sand.
  • 扩体预制桩由内芯预制桩和外围扩体材料固结体组成,一般可通过就地搅拌水泥土或机械扩孔并向孔内灌注扩体材料浆体,再植入内芯同轴预制桩的工艺成桩[1]。近年来,随着厂拌高质量扩体浆材的利用和成桩工艺的创新[2],扩体预制桩的桩身质量和承载性能大幅提高,开始被广泛应用于高层建筑、桩板路桥等承载力和沉降控制要求较高的建构筑物桩基础。

    与单一材料桩相比,扩体预制桩竖向承载时涉及到内芯预制桩、外围扩体和周围土3种不同材质的相互作用,荷载传递机制十分复杂[3-4]。然而,目前有关扩体预制桩承载机制的理论分析研究极为有限。考虑芯桩-扩体-周围土协同工作对承载性能的影响,Wang等[5]在假定芯桩与外围水泥土界面固定黏结的基础上,基于荷载传递法建立了劲性搅拌桩荷载传递简化计算模型;俞建霖等[6]针对刚性基础下砼芯水泥土桩复合地基,给出了砼芯、外围水泥土、周围土三者的受力平衡微分方程;结合引扩孔灌浆扩体预制桩工艺特点,张浩等[7]在考虑芯桩−扩体材料−周围土径向相互作用的基础上,提出了扩体预制桩整体复合模量表达式,并给出了全截面竖向受压模式下扩体桩的承载变形解析计算方法;考虑芯桩-水泥土-周围土界面的非线性力学特征,Zhu等[8]通过三折线模型建立了刚性基础下劲性搅拌桩的荷载传递理论分析模型,极大完善了扩体预制桩荷载传递机制的分析理论。然而,这些研究多集中于桩顶全截面受压工况,即桩顶处芯桩与扩体近似等沉变形受荷,较适用于建筑筏板、桥台等刚性基础下的桩基工况。

    近年来,随着扩体预制桩在桩板式无土路基结构、高桩承台结构中的广泛应用,芯桩单独承载的扩体预制桩荷载传递机制研究不容忽视。既有试验研究表明[9],芯桩单独承压下扩体预制桩中芯桩-扩体-土的荷载传递机制与全截面受压下显著不同。顾士坦等[10]基于Mindlin弹性解和剪滞原理,仅考虑芯桩-外围水泥土的相互作用,给出了芯桩单独受力时劲性搅拌桩芯桩荷载传递规律的理论分析方法;任连伟等[11]、Zhou等[12]基于荷载传递法,分别给出了高喷插芯组合桩和静钻根植桩芯桩受力时荷载传递简化计算方法,但芯桩-水泥土界面完全线弹性剪切模型的假定,无法全面反映此类扩体预制桩的工作状态。

    据此,本文针对芯桩单独承载工况下全长扩体预制桩中芯桩-扩体-周围土的荷载传递特点,在充分考虑芯桩-扩体、扩体-土界面非线性相互作用特性以及两界面剪切性能发挥交互影响的基础上,构建了全长扩体预制桩芯桩承载荷载传递计算分析模型,并参考工程实例,对全长扩体预制桩芯桩承载荷载传递机制及其影响因素进行了分析。

    图 1为典型全长扩体预制桩芯桩单独承载工况下桩土作用示意图。荷载P直接作用于内芯预制桩顶部,芯桩首先向下发生相对外围扩体材料的变形,并在芯桩-扩体界面产生剪切作用;外围扩体材料发生剪切变形,进一步将荷载传递至周围土。同时随着荷载向下传递,桩端处发生向下刺入变形趋势,并受到地基土的反力作用。因此芯桩承载工况下上部荷载是经由芯桩-扩体-周围土相互作用传递至周围土的。芯桩承载扩体预制桩的承载变形机制涉及到芯桩-扩体、扩体-周围土内外两个相互作用界面的协同工作。扩体预制桩整体受到周围土侧阻力τs、芯桩端阻力Rp和扩体端阻力Rc,内部则存在芯桩与扩体材料的剪切作用力τc

    图  1  承载变形机制示意图
    Figure  1.  Diagrammatic sketch of bearing mechanism

    图 1扩体预制桩的承载变形特征,可得芯桩竖向变形s(z)的表达式:

    s(z)=spc(z)+sc(z)+scs(z)+ss(z)=sp(z)+ss(z) (1)

    式中:spc(z),scs(z)为深度z处芯桩-扩体界面和扩体-周围土界面的相对位移;sc(z),ss(z)为z处扩体材料和周围土的剪切变形;sp(z)为芯桩与土的相对位移。

    既有研究表明[8, 13],理想弹塑性模型能较好的反映水泥土、水泥砂浆等水泥基胶凝材料与土之间的荷载传递关系,因此将扩体-周围土界面模型简化为理想弹塑性模型,如图 2(a)所示,即

    τs={ksscs (scs<sfs)τfs(scssfs) (2)
    图  2  界面模型示意图
    Figure  2.  Diagrammatic sketch of interface model

    式中:τsscs为扩体-土界面侧阻力和相对位移;ks为扩体-土剪切刚度;τfssfs界面极限侧阻力和极限相对位移。

    芯桩与水泥土等水泥基胶凝材料的剪切试验[14]表明,在较小相对位移下界面剪切应力即表现出较高水平,随着相对位移的增大,界面应力达到极限值后会出现迅速降低,界面呈脆性破坏模式,如图 2(b)所示。因此,芯桩-扩体界面不宜采用完全线弹性剪切模型,可简化为图 2(b)理想弹性-破坏模型,即

    τc={kcspc(spc<sfc)τfc (spcsfc) (3)

    式中:τcspc为芯桩-扩体界面侧阻力和相对位移;kc为芯桩-扩体界面剪切刚度,由试验得τc-spc曲线的直线段斜率确定;τfcsfc为界面极限侧阻力和相对位移;τfc为界面残余侧阻力,与扩体材料水泥掺量、龄期等有关,约为界面极限侧阻力的(0.057~0.153)倍[14]

    据此,分析时作如下基本假定:①内芯预制桩和外围扩体材料均为线弹性体;②针对目前以水泥土、水泥砂浆为扩体材料的扩体预制桩,扩体剪切刚度远大于周围土,且厚度相对有限[1-2],受荷时扩体剪切变形较小,忽略不计;③假定芯桩-扩体荷载传递为理想弹性-破坏模型,扩体-土荷载传递为理想弹塑性模型。

    图 3所示,取深度z处“扩体预制桩-土”单元进行受力分析,根据单元体静力平衡条件可得

    图  3  单元体受力变形示意图
    Figure  3.  Diagrammatic sketch of forces and deformations on elements
    dPp(z)=τc(z)up dz, dPc(z)=τc(z)up dzτs(z)uc dz} (4)

    式中:Pp(z),Pc(z)为z处芯桩和扩体的截面轴力;τc(z),τs(z)为z处芯桩-扩体、扩体-土界面的侧阻力;upuc为芯桩、扩体外侧周长。

    若受力变形后芯桩单元体和扩体单元体的压缩量分别为ΔpΔc,根据图 3(b)变形模式可得

    dspc(z)=ΔcΔpdscs(z)=Δc} (5)

    进而,基于胡克定律由式(4),(5)联合,可得芯桩-扩体、扩体-周围土的界面微分控制方程:

    d2spc(z)dz2=(upEcAc+upEpAp)τc(z)ucEcAcτs(z),d2scs(z)dz2=ucEcAcτs(z)upEcAcτc(z)} (6)

    式中:EpEc为芯桩、扩体的弹性模量;ApAc为芯桩、扩体的横截面积。

    芯桩承载扩体预制桩与土相互作用模式见图 4

    图  4  扩体预制桩-土相互作用模式
    Figure  4.  Reamed precast pile-soil interaction model

    随着桩顶荷载增加,扩体-土界面和芯桩-扩体界面将由完全弹性状态、从浅而深逐渐进入塑性或破坏状态。设桩长L,扩体-土界面塑性区深度L1,芯桩-扩体界面破坏区深度L2。一般情况下,由于芯桩-扩体界面剪切强度比扩体-土界面要大,工程实践中为充分发挥芯桩的承载性能,往往会通过调配扩体材料强度和厚度来保证芯桩-扩体共同作用,避免芯桩-扩体界面提前破坏[1],此时L1往往大于L2,如图 4(a)

    当外围扩体强度较小或厚度较薄等特殊情况下,芯桩单独承载桩顶附近芯桩和扩体应力差异显著,芯桩-扩体界面可能先出现破坏,而后扩体-土界面逐渐进入塑性状态,桩-土相互作用为图 4(b)所示模式。

    (1)L2L1工况

    根据芯桩-扩体、扩体-土两界面相互作用状态的不同,将扩体预制桩划分为3个区域,见图 4(a)

    a)当L1 < zL时,芯桩-扩体界面、扩体-土界面均处于弹性阶段,由式(2),(3),(6)可得

    d2spc(z)dz2=αspc(z)+βscs(z),d2scs(z)dz2=γspc(z)+δscs(z)} (7)

    式中:α=upkc/EcAc+upkc/EpApβ=ucks/EcAcγ=upkc/EcAcδ=ucks/EcAc

    求解方程,可得芯桩-扩体的相对位移spc1(z)、扩体-周围土的相对位移scs1(z)

    spcl(z)=1γ(C1a21ea1z+C2a22ea2z+C3a23ea3z+C4a24ea4z)δγ(C1ea1z+C2ea2z+C3ea3z+C4ea4z),scs1(z)=C1ea1z+C2ea2z+C3ea3z+C4ea4z} (8)

    式中:C1C2C3C4为方程待定系数;

    a1=[(δ+α)+(δ+α)24(δαγβ)]/2
    a2=[(δ+α)+(δ+α)24(δαγβ)]/2
    a3=[(δ+α)(δ+α)24(δαγβ)]/2
    a4=[(δ+α)(δ+α)24(δαγβ)]/2

    由式(5),(8)可得芯桩轴力和扩体轴力:

    Pp1(z)=EpAp(δγγ)(C1a1ea1z+C2a2ea2z+C3a3ea3z+C4a4ea4z)EpApγ(C1a13ea1z+C2a23ea2z+C3a33ea3z+C4a43ea4z)Pc1(z)=EcAc(C1a1ea1z+C2a2ea2z+C3a3ea3z+C4a4ea4z)} (9)

    此时,由式(1)可得芯桩竖向变形s1(z)

    s1(z)=sp1(z)+ss1(z) = spc1(z)+scs1(z)+ss1(z) (10)

    式中:ss1(z)=τs(z)rcln(rm/rc)/Gsrc为扩体外壁半径;Gs为桩周土剪切模量;rm为影响半径[15]

    b)当L2 < zL1时,扩体-周围土界面处于塑性状态,芯桩-扩体界面处于弹性阶段,

    d2spc(z)dz2=αspc(z)+κd2scs(z)dz2=γspc(z)κ} (11)

    式中:κ=ucτfs/EcAc

    求解方程,可得z处芯桩-扩体相对位移spc2(z)、扩体-周围土相对位移scs2(z)

    spc2(z)=D1eαz+D2eαzκαscs2(z)=γD1αeαz+γD2αeαzκ2(γα+1)z2+D3z+D4} (12)

    式中:D1D2D3D4为方程待定系数。

    进而,芯桩轴力Pp2(z)和扩体轴力Pc2(z)

    Pp2(z)=EpAp[γ+αα(D1eαzD2eαz)κ(γα+1)z+D3]Pc2(z)=EcAc[γD1αeαzγD2αeαzκ(γα+1)z+D3]} (13)

    此时,芯桩竖向变形s2(z)

    s2(z)=sp2(z)+ss2(z) = spc2(z)+scs2(z)+ss2(z) (14)

    式中:ss2(z)=τsfrcln(rmrc)/τsfrcln(rmrc)GsGs

    c)当0≤zL2时,扩体-周围土界面处于塑性状态,芯桩-扩体界面处于破坏阶段,

    d2spc(z)dz2=g, d2scs(z)dz2=f} (15)

    式中:g=(upEcAc+upEpAp)τfcucτfsEcAcf=ucτfsEcAcupτfcEcAc

    求解方程,可得z处芯桩-扩体相对位移spc3(z)、扩体-周围土相对位移scs3(z)

    spc3(z)=12gz2+F1z+F2scs3(z)=12fz2+F3z+F4} (16)

    式中:F1F2F3F4为方程待定系数。

    进而,芯桩轴力Pp3(z)和扩体轴力Pc3(z)

    Pp3(z)=EpAp(gz+fz+F1+F3)Pc3(z)=EcAc(fz+F3)} (17)

    此时,芯桩竖向变形s3(z)

    s3(z)=spc3(z)+scs3(z)+ss3(z) (18)

    式中:ss2(z)=τsfrcln(rm/rc)/Gs

    d)边界与连续条件

    根据扩体桩受力特点,桩顶边界条件为

    Pp(z)|z=0=PPc(z)|z=0=0} (19)

    桩端边界条件为

    Pp(z)|z=L=kLApsp(z)|z=LPc(z)|z=L=kLAcscs(z)|z=L} (20)

    式中:kL为桩端土刚度系数,基于Randolph等[16]刚性体压入弹性半空间解答,kL=4GL/4GL[πrc(1ν)][πrc(1ν)]GL为桩端土体剪切模量;ν为桩端土泊松比。

    桩身受力变形连续条件:

    spc2(z)|z=L1=spc1(z)|z=L1+scs2(z)|z=L1=scs1(z)|z=L1+Pp2(z)|z=L1=Pp1(z)|z=L1+Pc2(z)|z=L1=Pc1(z)|z=L1+scs2(z)|z=L1=γD1αeαL1+γD2αeαL1 κ2(γα+1)L12+D3L1+D4=sfs} (21)
    spc3(z)|z=L2=spc2(z)|z=L2+scs3(z)|z=L2=scs2(z)|z=L2+Pp3(z)|z=L2=Pp2(z)|z=L2+Pc3(z)|z=L2=Pc2(z)|z=L2+spc3(z)|z=L2=12gL22+F1L2+F2=sfc} (22)

    综上,针对L2L1工况,根据边界和连续条件(19)~(22),可得关于C1~C4D1~D4F1~F4L1、L2的14个方程,求解未知量并回代入式(8)~(10)、(12)~(14)、(16)~(18),即可得扩体桩桩身轴力(Pp(z),Pc(z))与位移(spc(z),scs(z),s(z))。

    考虑不同计算条件影响:①当桩顶荷载相对较小时,全桩长范围芯桩-扩体、扩体-土界面均处于弹性状态时(L1=0,L2=0),上述计算可退化为关于C1~C4共4个未知量的求解问题,仅需将式(8)、(9)代入边界条件(19)、(20)即可得出4个求解方程,进而可得该工况下桩身受力响应;②若随着荷载增加,扩体-土界面开始部分出现塑性区、但芯桩-扩体界面全桩长处于弹性阶段时(0 < L1 < LL2=0),上述计算退化为关于C1~C4D1~D4L1共9个未知量的求解问题,此时将式(8),(9),(12),(13)代入边界和连续条件(19)~(21)即可得出9个求解方程,进而可得该工况下桩身受力响应;③随着桩顶荷载的增加,芯桩-扩体界面也开始出现脆性破坏(0 < L1 < L,0 < L2 < L1),即进入式(7)~(18)计算工况;④当扩体-土界面全部进入塑性区而芯桩-扩体界面局部出现破坏时(L1=L,0 < L2 < L1),上述计算退化为关于D1~D4F1~F4L2共9个未知量的求解问题,此时将式(12),(13),(16),(17)代入边界和连续条件(19)~(21)即可得出9个求解方程;⑤当扩体-土界面全部进入塑性区且芯桩-扩体界面全部破坏时(L1= L2=L),上述计算退化为关于F1~F4共4个未知量的求解问题,将式(16),(17)代入边界条件(19),(20)即可得出4个求解方程;⑥若实际工况是非均质层状土,可通过桩长范围内各土层参数的加权平均进行简化计算;当然也可通过桩身离散的矩阵系数传递法,采用界面相对位移迭代计算的方法近似求解。

    (2)L2 > L1工况

    图 4(b)所示,当芯桩-扩体破坏区深度L2大于外部扩体-土塑性区深度L1时,沿桩长同样可以将扩体预制桩划分为3个区域:①L2 < zL时,芯桩-扩体界面、扩体-土界面均处于弹性阶段;②L1 < zL2时,芯桩-扩体界面处于破坏状态,而外部扩体-土界面仍处于弹性状态;③0 < zL1时,扩体-土界面处于塑性状态,芯桩-扩体界面处于破坏阶段。

    其中,在L2 < zL和0 < zL1区域,桩-土相互作用状态与前述2.3节中扩体桩上、下两段区域桩-土相互作用状态相同,芯桩-扩体、扩体-土两界面的微分控制方程与式(7),(15)相同,求解方法也相同,在此不再赘述。

    L1 < zL2时,由式(2),(3),(6)可得此时芯桩-扩体和扩体-土界面微分方程:

    d2spc(z)/dz2=ξd2scs(z)/dz2=δscs(z)ζ} (23)

    式中:ξ=upτfc/EpApζ=upτfc/EcAc

    求解方程,可得该区域z处芯桩位移spc2(z)、扩体-周围土相对位移scs2(z)

    spc2(z)=12ξz2+J1z+J2zscs2(z)=J3eζz+J4eζz+ζ/ξ} (24)

    式中:J1J2J3J4为方程待定系数。

    进而,芯桩轴力Pp2(z)、扩体轴力Pc2(z)和芯桩竖向变形s2(z)

    Pp2(z)=EpAp[ξz+J1]Pc2(z)=EcAc[J3eζz+J4eζz]} (25)
    s2(z)=spc2(z)+scs2(z)+ss2(z) (26)

    式中:s2(z)=τs(z)rcln(rm/rc)/Gs

    据此,可对L2 > L1工况下扩体预制桩荷载传递进行计算分析。其中,桩顶和桩端边界条件同式(19),(20);桩身受力变形连续条件为

    spc2(z)|z=L2=spc1(z)|z=L2+scs2(z)|z=L2=scs1(z)|z=L2+Pp2(z)|z=L2=Pp1(z)|z=L2+Pc2(z)|z=L2=Pc1(z)|z=L1+spc2(z)|z=L2=sfc} (27)
    spc3(z)|z=L1=spc2(z)|z=L1+scs3(z)|z=L1=scs2(z)|z=L1+Pp3(z)|z=L1=Pp2(z)|z=L1+Pc3(z)|z=L1=Pc3(z)|z=L1+scs3(z)|z=L1=sfs} (28)

    式中:spc1(z)scs2(z)Pp1Pc1(z)分别为L2 < zL芯桩和扩体的变形和轴力;spc2(z)scs2(z)Pp2(z)Pc2(z)分别为L1 < zL2芯桩和扩体的变形和轴力;spc3(z)scs3(z)Pp3(z)Pc3(z)分别为0 < zL1芯桩和扩体的变形和轴力。

    考虑不同计算条件,采用与2.3节相同的求解方法,可得L2 > L1工况下扩体桩桩身轴力与位移。

    水泥土扩体预制桩现场试验参见文献[17]。试桩桩长10 m,芯桩为PHC-400AB(95)管桩,外围为厚300 mm水泥土桩,水泥掺量18%,水泥土90 d抗压强度≥1.6 MPa。计算中,C80管桩弹性模量38 GPa,水泥土弹性模量800 MPa。针对界面荷载传递模型参数取值,根据文献[14]混凝土-水泥土界面极限侧阻力约抗压强度的0.064~0.259,临界相对位移约1.23~2.40 mm,界面残余侧阻力约为极限侧阻力的0.057~0.153,取芯桩-水泥土界面极限侧阻力为250 kPa,临界相对位移sfc=2 mm,残余侧阻力τfc=25 kPa。基于剪切位移法[15],以芯桩-水泥土极限剪应力估算外围水泥土剪切变形约为1.4×10-4 m。可见,虽然外围水泥土包裹厚度达300 mm,其剪切变形仍较小,可忽略不计。水泥土-周围土界面荷载传递参数取值,以文献[13]水泥土桩-土试验结果作为参考,取极限侧阻力τfs=65 kPa,ks=20 kPa/mm。依据桩底土体的性质,结合文献[17]现场载荷试验桩端阻力和桩端沉降量,反演计算取桩端土刚度系数kL为1000 kPa/mm。

    采用本文方法对现场试桩进行计算,如图 56。各级荷载作用下,桩顶沉降计算值大于实测值,芯桩轴力计算值整体上与实测值相比偏小,这是由于实际芯桩后植入过程中会对水泥土和周围土产生挤密作用,使界面黏结强度增大[2],桩体整体承载性能提升,但计算分析中界面参数未考虑该部分的影响。

    图  5  实测与计算沉降-荷载曲线对比
    Figure  5.  Comparison between measured and calculated curves of displacement-load
    图  6  芯桩轴力与桩身压缩量的实测与计算结果对比
    Figure  6.  Comparison of measured and calculated axial force of core pile and settlement compression of pile

    从芯桩轴力分布规律来看,芯桩轴力沿深度逐渐减小,且随着上部荷载增加,该变化趋势越来越明显。各级荷载下桩端轴力约占桩顶总荷载18.4%~21.3%,说明水泥土扩体预制桩整体上表现为端承摩擦桩的承载特性。进一步对比芯桩轴力计算结果与现场实测结果,两者在桩顶、桩端附近差异性较小,这主要因为理论计算中桩端刚度由现场实测反演所得,且桩顶荷载边界与现场试桩相同。在桩身中部范围芯桩轴力存在差异,且实测桩身轴力沿桩身分布有波动性,这是由于实际桩周水泥土就地搅拌时均匀性不能得到很好保证,从而影响内部芯桩受力。参考文献[17]对试桩沉降压缩量进行计算分析,如图 6(b)表 1所示,可以看出本文方法也能较好反映桩身压缩变形特征。

    表  1  试桩沉降压缩量对比分析
    Table  1.  Comparative analysis of settlement compression of pile
    计算分析 桩顶荷
    载/kN
    桩顶
    沉降/mm
    桩身压
    缩/mm
    桩端沉
    降/mm
    桩身压缩
    沉降比/%
    试桩[17] 2200 10.23 2.82 7.41 27.57
    本文方法 2200 9.83 2.73 7.10 27.76
    下载: 导出CSV 
    | 显示表格

    德上高速合肥-枞阳段地处安徽中南部,部分路段采用了“桩柱一体式”桩板路基结构。考虑沿线基岩较浅且上覆沉积物工程性能较好,传统打入式预制桩沉桩困难,因此采用了引扩孔灌浆植桩工艺。芯桩采用PHC-500AB(100)管桩,外围采用厂拌充盈水泥砂浆,28 d抗压强度 > 15 MPa,桩长10 m。为揭示水泥砂浆扩体预制桩的工程特性,开展扩孔直径分别为700,800,1000,1200 mm(对应扩体厚度100,150,250,350 mm)的静载试验,试验区土层性质见表 2

    表  2  土层主要物理力学指标
    Table  2.  Physico-mechanical indexes of soil layers
    土层 h/
    m
    γ/
    (kN·m-3)
    c/
    kPa
    ϕ/
    (°)
    Es1-2/
    MPa
    fsk/
    kPa
    ②粉土 6.1 19.3 20.6 10.2 7.50 140
    ②1粉质黏土 2.4 19.1 34.6 12.2 8.59 180
    ③粗砂 2.8 250
    ⑦1全风化角砾岩 51.9 12.5 350
    下载: 导出CSV 
    | 显示表格

    为了获取芯桩-砂浆、砂浆-土界面荷载传递参数,采用自研界面剪切装置(图 7)开展界面剪切试验,试验方法同直接剪切试验,每组试验3个试样,分别考虑100,200,300 kPa 3种不同的法向应力,参考《土工试验方法标准》(GB/T50123—2019),每个法向力采用一次轻轻施加。受制于剪切试样尺寸(ϕ61.8 mm×20 mm)较小,芯桩混凝土采用同强度C80灌浆料替代,所用土样和水泥砂浆皆取自工程现场。

    图  7  界面剪切试验示意图
    Figure  7.  Diagrammatic sketch of interface shear tests

    制作水泥砂浆-土试样时,为保证试样土重度与工程实际的一致性,首先参考地勘报告调整土样含水率,并以勘察报告土体重度为基准,根据下剪切盒体积,计算每个试样所需土的重量,将其填入下剪切盒至剪切界面高度,然后在其上剪切盒内浇入水泥砂浆。制作芯桩-水泥砂浆试样时,首先将替代灌浆料浇入模具中,待标准养护28 d后再在其上浇入水泥砂浆。制样完成后,所有试样标准养护14 d后开展界面剪切试验,部分成果如图 8所示。其中,散点曲线为试验测试结果,黑色实折线为试验数据拟合模型曲线。

    图  8  界面剪切试验结果
    Figure  8.  Results of interface shear tests

    依据界面剪切试验,取芯桩-砂浆界面极限侧阻力368 kPa,临界相对位移2.3 mm,残余侧阻力74 kPa;取砂浆-土界面极限摩阻力90 kPa,临界相对位移3 mm。计算中,内芯预制桩弹性模量38 GPa,扩体砂浆弹性模量7 GPa,参考《公路桥涵地基与基础设计规范》(JTG3363—2019)桩端岩石地基竖向抗力系数取值,桩端全风化角砾岩刚度系数kL取2500 kPa/mm。

    图 9为不同扩径工况扩体预制桩荷载-位移曲线。随着扩径直径增加,扩体预制桩的承载性能大体呈增大变化趋势,当扩径直径由700 mm增加到1000 mm时,同级荷载桩顶沉降降低约44.3%~59.5%;而扩径的进一步增加(1200 mm),桩顶沉降变化量降低。

    图  9  实测值与计算值对比
    Figure  9.  Comparison between measured and calculated values

    采用本文方法计算的荷载-位移曲线与实测结果具有较好的一致性,同级荷载作用下计算桩顶沉降相较实测值整体偏大。当扩径为700 mm时,各级荷载水平下桩顶沉降平均误差率6.7%;随着扩径直径增大,计算值误差率随之增大。当扩径直径由700 mm增加至800,1000,1200 mm时,桩顶沉降平均误差率分别增加至24.8%,34.6%,31.1%。分析其原因,上述计算分析未考虑芯桩后植入施工的扩挤影响,尤其是后植入芯桩对孔内砂浆挤密及砂浆渗滤增强桩-土界面剪切性能的作用效应,同时简化理论模型中忽略外围扩体剪切变形的基本假定也对计算结果有影响。考虑相关技术规程和工程实践中[1]为了充分发挥芯桩高承载力的特性和合理控制孔内灌浆后植入桩的挤土效应,推荐扩孔直径宜为芯桩直径的1.5倍,即扩体厚度多控制在100~200 mm,因此本文计算方法适用于常规扩体预制桩的工程计算分析。

    基于上述实例,进一步计算分析复合桩身荷载传递机制。以扩径直径700 mm(内芯预制桩500 mm)水泥砂浆扩体预制桩为例,图 10所示为计算的桩身轴力分布曲线。不同荷载水平下芯桩轴力沿桩身呈逐渐减小的变化趋势,随着桩顶荷载水平的增大,在上部0~3 m内芯桩轴力减小幅度相对较大,3~8 m内芯桩轴力减小幅度变小且相对稳定,而桩端附近(8~10 m)芯桩轴力减小幅度再次变大。

    图  10  扩体预制桩桩身轴力
    Figure  10.  Axial forces of reamed precast pile

    与之对应,外围水泥砂浆截面轴力沿桩身呈逐渐增大变化趋势,在0~3 m内外围砂浆轴力快速增加,说明芯桩单独承载工况下,桩顶附近芯桩和外围砂浆的相互作用较强;随着深度增加,外围砂浆轴力增加幅度降低且相对稳定,与芯桩趋于共同作用;桩端(8~10 m)附近,由于芯桩与外围砂浆刚度的显著差异,以及桩端角砾岩较高桩端反力作用,芯桩相对外围砂浆“下刺入”过程中,两者相互作用增强,外围砂浆轴力呈明显增大的变化趋势,如图 10(b)所示。

    图 11所示为水泥砂浆扩体预制桩内、外两界面阻力分布曲线。内侧芯桩-水泥砂浆界面阻力沿桩身呈两端大、中间小的分布规律。这主要是由于在桩顶、桩端附近芯桩和外围砂浆的刚度差异导致两者界面相对位移趋势明显,引起较强烈的界面相互作用。其中桩顶荷载4000 kN时,芯桩-水泥砂浆界面在桩顶0.2 m范围内出现剪切破坏,界面阻力降为残余侧阻力。

    图  11  扩体预制桩界面阻力
    Figure  11.  Interfacial resistances of reamed precast pile

    对于外侧水泥砂浆-土界面,界面阻力沿桩身整体呈逐渐减小分布。当桩顶荷载较大时,水泥砂浆-土界面由上而下逐渐进入塑性状态,界面阻力达到极限侧阻力。随着桩顶荷载增加,界面塑性区长度增大,桩顶荷载3000,3500,4000 kN时,对应水泥砂浆-土界面塑性区长度分别为3.8,6.0,7.5 m。

    进一步分析不同扩孔直径对扩体桩承载特性的影响,以桩顶荷载3000 kN工况为例,如图 12所示。随深度增加芯桩轴力逐渐减小,桩端处芯桩轴力占总截面轴力约56.1%,45.9%,33.3%,26.5%,扩体直径越大,桩端处芯桩轴力占比越小;砂浆截面轴力随深度增加而逐渐增大,占总截面轴力比例也逐渐增大,桩端处砂浆截面轴力占比约43.9%,54.1%,66.7%,73.5%。外围水泥砂浆对芯桩具有“扩径”作用。

    图  12  不同扩径情况下桩身轴力
    Figure  12.  Axial forces under different reaming diameters

    当扩体直径相对较大(1200 mm)时,在桩顶(1 m)附近芯桩-砂浆的荷载传递作用出现明显降低现象,这主要是因为,随着扩体截面增大,其单位面积应力水平降低,单独承载芯桩与外围砂浆的相对位移势能增强,两者共同作用逐渐减低。由此可见,虽然扩体直径增加可以提升复合桩身的整体承载性能,但较厚外围扩体也可能制约内部芯桩承载能力的发挥。工程实践中宜注意芯桩和外围扩体的工作协同性。针对水泥砂浆扩体预制桩,以工程现场试验和本计算算例分析来看,对于直径500 mm的芯桩,外围砂浆扩体直径不宜大于1000 mm。

    考虑外围扩体材料性能的影响,计算对比扩体弹性模量为0.8,4,7,10 GPa时扩体预制桩桩身内力,如图 13。随着扩体材料模量增大,芯桩轴力逐渐减小而外围扩体截面轴力逐渐增大。当扩体模量由较低的0.8 GPa(如水泥土),提升至7~10 GPa(如水泥砂浆)时,芯桩轴力由总截面轴力的93%,降低至77%~80%;而外围扩体截面轴力由总截面轴力的6%,提升至19%~23%,说明外围扩体材料刚度的改善,可提升内芯预制桩向周围传递荷载的效率。

    图  13  不同扩体材料情况下桩身轴力
    Figure  13.  Axial forces under different peripheral materials

    由扩体截面轴力分布曲线可以看出,当扩体采用刚度相对较低(0.8 GPa)的水泥土时,芯桩与扩体的显著相互作用仅发生在桩顶(1 m)和桩端(2 m)附近有限区域内;随着扩体刚度增加,芯桩与扩体的显著相互作用区域逐渐增大,两者共同作用效应增强。在工程实践中可通过改善扩体材料性能,增强内芯预制桩与外围扩体的共同作用。

    (1)根据芯桩承载全长扩体预制桩的受力变形特点,综合考虑芯桩-扩体、扩体-周围土内外两界面剪切特性的发挥与耦合作用,推导建立了芯桩承载扩体预制桩的荷载传递理论计算方法,工程案例分析验证了该方法的可靠性。

    (2)外围扩体可改善桩身整体承载性能,当扩体直径(芯桩直径500 mm)由700 mm增加至1000 mm时,桩顶沉降降低约44.3%~59.5%;外围扩体可将芯桩荷载有效传递至周围土体;竖向荷载作用下,随着深度增加,芯桩轴力逐渐减小,而外围扩体轴力逐渐增大,对芯桩具有“扩径”作用。

    (3)扩体刚度的增加可提升芯桩荷载传递效率,当扩体弹性模量由0.8 GPa增加至7~10 GPa时,芯桩轴力占总截面轴力比例由93%降低至77%~80%,而扩体截面轴力占比由6%增加至19%~23%,芯桩与扩体共同作用效应增强。

  • 图  1   GDS界面剪切仪

    Figure  1.   GDS interface shearing instrument

    图  2   轴向压杆

    Figure  2.   Axial pressure bar

    图  3   直剪仪

    Figure  3.   Direct shear apparatus

    图  4   颗粒级配曲线

    Figure  4.   Grain-size distribution curve of particles

    图  5   界面剪应力–剪切位移关系曲线

    Figure  5.   Shear stress–shear displacement relationships of interface under different relative densities

    图  6   砂土剪应力–剪切位移关系曲线

    Figure  6.   Shear stress–shear displacement curves of saturated sand under different relative densities

    图  7   不同相对密实度下界面竖向位移–剪切位移关系

    Figure  7.   Vertical displacement–shear displacement relationship curves of interface under various relative densities

    图  8   剪切试验前后颗粒级配曲线对比

    Figure  8.   Comparison between grain-size distribution curves of particles before and after shear test

    图  9   砂土挤出

    Figure  9.   Sand to be squeezed out

    图  10   不同初始相对密实度下砂土界面抗剪强度包络线

    Figure  10.   Shear strength envelopes of sand interface under different initial relative densities

    图  11   界面摩擦角–内摩擦角曲线

    Figure  11.   Relationship between interfacial friction angle and internal friction angle

    图  12   界面摩擦系数变化规律

    Figure  12.   Variation laws of friction coefficient of interface

    图  13   沉贯摩阻力变化规律

    Figure  13.   Variation of friction resistance with penetration depth in sand with different densities

    表  1   试验用砂相关物理参数

    Table  1   Relevant physical parameters of sand for tests

    试验用砂 Cu Cc D50/mm emax emin
    细砂 2.33 0.98 0.18 0.980 0.560
    下载: 导出CSV

    表  2   法向应力与沉贯深度对应关系

    Table  2   Normal stresses under various penetration depths

    土压力/kPa 10 20 30 50 100 150 200
    对应外壁沉贯深度/m 0.35 0.60 0.83 1.30 2.50 3.70 5.00
    对应内壁沉贯深度/m 0.55 1.06 1.56 2.58 5.00
    下载: 导出CSV

    表  3   试验工况

    Table  3   Experimental programs

    剪切类型 相对密实度 法向应力/kPa
    界面剪切 0.3, 0.5, 0.7, 0.9, 1.0 10, 20, 30, 50, 100, 150, 200
    直剪 0.3, 0.5, 0.7, 0.9, 1.0
    下载: 导出CSV

    表  4   抗剪强度–法向应力线性回归方程

    Table  4   Shear strength–linear regression equations for vertical pressure

    相对密实度 a b c R2
    1.0 0.00021 0.3050 -0.1190 0.997
    0.9 0.00015 0.2871 -0.1090 0.999
    0.7 0.00027 0.2885 -0.6047 0.999
    0.5 0.00060 0.2741 -0.8873 0.999
    0.3 0.00018 0.2206 -0.1780 0.997
    下载: 导出CSV

    表  5   不同相对密实度下砂土抗剪强度指标

    Table  5   Parameters of sand under different relative densities

    Dr 界面摩擦角δ/(°) 内摩擦角φ/(°) m
    1.0 18.12 54.07 0.335
    0.9 17.23 50.89 0.338
    0.7 16.13 47.20 0.342
    0.5 15.51 43.50 0.357
    0.3 14.05 39.15 0.359
    下载: 导出CSV
  • [1] 李大勇, 张雨坤, 高玉峰, 等. 中粗砂中吸力锚的负压沉贯模型试验研究[J]. 岩土工程学报, 2012, 34(12): 2277–2283. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201212022.htm

    LI Da-yong, ZHANG Yu-kun, GAO Yu-feng, et al. Model tests on penetration of suction anchors in medium-coarse sand[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(12): 2277–2283. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201212022.htm

    [2] 王胤, 朱兴运, 杨庆. 考虑砂土渗透性变化的吸力锚沉贯及土塞特性研究[J]. 岩土工程学报, 2019, 41(1): 184–190. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201901025.htm

    WANG Yin, ZHU Xing-yun, YANG Qing. Installation of suction caissons and formation of soil plug considering variation of permeability of sand[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(1): 184–190. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201901025.htm

    [3] 丁红岩, 张浦阳. 海上吸力锚负压下沉渗流场的特性分析[J]. 海洋技术, 2003, 22(4): 44–48. doi: 10.3969/j.issn.1003-2029.2003.04.011

    DING Hong-yan, ZHANG Pu-yang. Suction penetration seepage field's characteristics of suction anchor[J]. Ocean Technology, 2003, 22(4): 44–48. (in Chinese) doi: 10.3969/j.issn.1003-2029.2003.04.011

    [4]

    HOULSBY G T, BYRNE B W. Design procedures for installation of suction caissons in sand[J]. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering, 2005, 158(3): 135–144. doi: 10.1680/geng.2005.158.3.135

    [5]

    HOULSBY G T, BYRNE B W. Design procedures for installation of suction caissons in clay and other materials[J]. Proceedings of the Institution of Civil Engineers- Geotechnical Engineering, 2005, 158(2): 75–82. doi: 10.1680/geng.2005.158.2.75

    [6]

    MEHRAVAR M, HARIRECHE O, FARAMARZI A, et al. Modelling the variation of suction pressure during caisson installation in sand using FLAC3D[J]. Ships and Offshore Structures, 2017, 12(7): 893–899. doi: 10.1080/17445302.2015.1051311

    [7] 李大勇, 吴宇旗, 张雨坤, 等. 砂土中桶形基础吸力值的设定范围[J]. 岩土力学, 2017, 38(4): 985–992, 1002. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201704010.htm

    LI Da-yong, WU Yu-qi, ZHANG Yu-kun, et al. Determination of suction range for penetration of suction caissons in sand[J]. Rock and Soil Mechanics, 2017, 38(4): 985–992, 1002. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201704010.htm

    [8] 李大勇, 吴宇旗, 张雨坤. 砂土中吸力式桶形基础竖向承载性状研究[J]. 山东科技大学学报(自然科学版), 2016, 35(3): 33–39. doi: 10.3969/j.issn.1672-3767.2016.03.005

    LI Da-yong, WU Yu-qi, ZHANG Yu-kun. Vertical bearing capacity of suction caisson foundation in sand[J]. Journal of Shandong University of Science and Technology (Natural Science), 2016, 35(3): 33–39. (in Chinese) doi: 10.3969/j.issn.1672-3767.2016.03.005

    [9]

    CHEN F, LIAN J J, WANG H J, et al. Large-scale experimental investigation of the installation of suction caissons in silt sand[J]. Applied Ocean Research, 2016, 60: 109–120. doi: 10.1016/j.apor.2016.09.004

    [10]

    HARIRECHE O, MEHRAVAR M, ALANI A M. Suction caisson installation in sand with isotropic permeability varying with depth[J]. Applied Ocean Research, 2013, 43: 256–263. doi: 10.1016/j.apor.2013.10.008

    [11] 闫澍旺, 霍知亮, 楚剑, 等. 黏土中桶形基础负压下沉阻力及土塞发展试验[J]. 天津大学学报(自然科学与工程技术版), 2016, 49(10): 1027–1033. https://www.cnki.com.cn/Article/CJFDTOTAL-TJDX201610004.htm

    YAN Shu-wang, HUO Zhi-liang, CHU Jian, et al. Experiment on penetration resistance and soil plug development during suction caisson penetration in soft clay[J]. Journal of Tianjin University (Science and Technology), 2016, 49(10): 1027–1033. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TJDX201610004.htm

    [12]

    LIU R, MA W G, QI Y, et al. Experimental studies on the drag reduction effect of bucket foundation installation under suction pressure in sand[J]. Ships and Offshore Structures, 2019, 14(5): 421–431. doi: 10.1080/17445302.2018.1500188

    [13]

    NARDELLI A, CACCIARI P P, FUTAI M M. Sand-concrete interface response: the role of surface texture and confinement conditions[J]. Soils and Foundations, 2019, 59(6): 1675–1694. doi: 10.1016/j.sandf.2019.05.013

    [14]

    TIWARI B, AL-ADHADH A R. Influence of relative density on static soil–structure frictional resistance of dry and saturated sand[J]. Geotechnical and Geological Engineering, 2014, 32(2): 411–427. doi: 10.1007/s10706-013-9723-6

    [15]

    TEJCHMAN J, WU W. Experimental and numerical study of sand-steel interfaces[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1995, 19(8): 513–536. doi: 10.1002/nag.1610190803

    [16]

    JITSANGIAM P, PRA-AI S, BOULON M, et al. Characterization of a soil–rough structure interface using direct shear tests with varying cyclic amplitude and loading sequences under a large cyclic testing cycle condition[J]. Acta Geotechnica, 2021, 17(5): 1829–1845.

    [17] 郭聚坤, 寇海磊, 许泓霖, 等. 桩–海洋黏土界面剪切性状试验研究[J]. 长江科学院院报, 2019, 36(4): 104–108, 117. https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201904024.htm

    GUO Ju-kun, KOU Hai-lei, XU Hong-lin, et al. Experimental study on shear behaviors of interface between pile and marine clay[J]. Journal of Yangtze River Scientific Research Institute, 2019, 36(4): 104–108, 117. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201904024.htm

    [18] 祁越, 刘润, 练继建. 无黏性土中筒型基础负压下沉模型试验[J]. 岩土力学, 2018, 39(1): 139–150. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801018.htm

    QI Yue, LIU Run, LIAN Ji-jian. Model test of bucket foundation suction installation in cohesionless soil[J]. Rock and Soil Mechanics, 2018, 39(1): 139–150. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801018.htm

    [19] 袁中立, 秦延龙, 唐海燕. 浅海桶形基础平台[M]. 北京: 石油工业出版社, 2010.

    YUAN Zhong-li, QIN Yan-long, TANG Hai-yan. Bucket Foundation Platform in Shallow Sea[M]. Beijing: Petroleum Industry Press, 2010. (in Chinese)

  • 期刊类型引用(1)

    1. 朱晟,陆雪妮,汤洪洁,孟涛,路德任. 大粒径砂砾料的缩尺试验方法及应用. 水利规划与设计. 2022(12): 5-10 . 百度学术

    其他类型引用(1)

图(13)  /  表(5)
计量
  • 文章访问数:  204
  • HTML全文浏览量:  44
  • PDF下载量:  68
  • 被引次数: 2
出版历程
  • 收稿日期:  2021-04-13
  • 网络出版日期:  2022-09-22
  • 刊出日期:  2022-08-31

目录

/

返回文章
返回