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

考虑超载和纵坡效应的隧道开挖面破坏分析

杨峰, 高益康, 秦傲韩, 宋智慧, 赵炼恒

杨峰, 高益康, 秦傲韩, 宋智慧, 赵炼恒. 考虑超载和纵坡效应的隧道开挖面破坏分析[J]. 岩土工程学报, 2024, 46(11): 2323-2332. DOI: 10.11779/CJGE20230733
引用本文: 杨峰, 高益康, 秦傲韩, 宋智慧, 赵炼恒. 考虑超载和纵坡效应的隧道开挖面破坏分析[J]. 岩土工程学报, 2024, 46(11): 2323-2332. DOI: 10.11779/CJGE20230733
YANG Feng, GAO Yikang, QIN Aohan, SONG Zhihui, ZHAO Lianheng. Failure analysis of tunnel face influenced by overload and longitudinal slope[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(11): 2323-2332. DOI: 10.11779/CJGE20230733
Citation: YANG Feng, GAO Yikang, QIN Aohan, SONG Zhihui, ZHAO Lianheng. Failure analysis of tunnel face influenced by overload and longitudinal slope[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(11): 2323-2332. DOI: 10.11779/CJGE20230733

考虑超载和纵坡效应的隧道开挖面破坏分析  English Version

基金项目: 

湖南省杰出青年科学基金项目 2021JJ10063

中国铁路广州局集团有限公司科技研究开发计划项目 2021K094-Z

详细信息
    作者简介:

    杨峰(1981—),男,博士,副教授,主要从事隧道与地下工程方面的教学和科研工作。E-mail: yf5754@csu.edu.cn

    通讯作者:

    赵炼恒, E-mail: zlh8076@163.com

  • 中图分类号: TU432

Failure analysis of tunnel face influenced by overload and longitudinal slope

  • 摘要: 针对地表超载和纵向坡度影响下的隧道开挖面稳定性问题,提出源于构造破坏机制的初始三角形网格更新策略,应用刚体平动运动单元上限分析(UB-RTME)方法开展非线性规划运算,获得地表临界超载系数σs/c图表与滑移线网破坏模式。结果表明:①以网格更新得到的UB-RTME滑移线网与Optum G2自适应加密网格,均可反映隧道开挖面破坏影响范围及承载关键区域,且UB-RTME以有效间断线及单元速度确定破坏范围,可揭示相对破坏范围S/D2定量数据及滑移线网破坏形态与影响因素的关联规律。②临界超载系数σs/c随纵向坡度角θ增大而显著降低,随埋深比C/D增大而增大,随重度系数γD/c增大而不断减小,随内摩擦角ϕ增大而非线性增大。③对于C/Dϕ数值较大的情况,网格密度将显著影响UB-RTME运算得到的σs/c数值解精度;网格密度增加后,σs/c计算结果精度显著提高,且体现破坏特征的网状滑移线变得更加光滑。
    Abstract: Aiming at the stability of excavation face of tunnels under the influence of surface overload and longitudinal slope, an initial triangular mesh updating strategy derived from the structural failure mechanism is proposed. The computation is conducted by using the upper bound method with the rigid translatory moving elements (UB-RTME) with nonlinear programming to generate stability charts for the load factor σs/c and obtain the slip line failure modes. The results show that the range affected by the failure modes and the major failure domain are given by both the UB-RTME slip line mesh and the Optum G2 adaptive encryption mesh. The range of failure area is determined by the effective discontinuities and element velocity characterized by the UB-RTME. The quantitative data of the failure area S/D2 and the relationship between the slip line failure modes and the influencing factors are provided. The load factor σs/c decreases significantly with the increase of the longitudinal slope angle θ, increases with the increasing buried depth ratio C/D. It decreases with the increase of the gravity coefficient γD/c while increases nonlinearly with the increase in the internal friction angle ϕ. For the case with large values of C/D and ϕ, the mesh density will significantly affect the accuracy of the σs/c numerical solution obtained by the UB-RTME operation. With the increase of the mesh density, the accuracy of the calculated results of σs/c is significantly improved, and the slip line reflecting the failure characteristics becomes smoother.
  • 随着环保理念的发展,众多学者尝试使用新型生物技术代替传统无机材料如硅酸盐水泥等进行土体加固[1-3]。生物聚合物土体加固(biopolymer reinforced soil,简称BPRS)是一种绿色环保的新型土体固化技术,近年来得到了国内外众多学者的广泛关注[4-6]。BPRS主要是利用微生物或藻类植物产生的胞外多糖与水和土体混合,饱和状态下以胶体形式填充于土体内部,干燥状态下其以纤维状或网状结构填充于土体孔隙中,粘结相邻土颗粒,以提高土体的物理力学性质[7]。Cabalar等[8]通过无侧限压缩试验研究了低掺量黄原胶对砂土的强度的影响。结果表明,1.5%含量黄原胶固化土的抗压强度提升了两倍以上。生物聚合物具有一定的水溶性,经加固处理的土壤遇水稳定性表现不佳。刘瑾等[9]通过CT扫描技术和力学测试研究了干湿循环条件下黄原胶加固土的劣化机制。结果表明:黄原胶能够提高土体的抗压强度和耐干湿循环效果,但当1.5%含量黄原胶固化黏土干湿循环4次时,抗压强度仍具有14.04%的损失。王天亮等[10]通过室内力学试验研究了黄原胶改良膨润土的力学性质。结果表明:经历3次干湿循环后,黄原胶改良土的无侧限抗压强度降低为未经干湿循环黄原胶改良土的0.3倍。Muguda等[11]通过浸渍试验和耐久性试验表明了黄原胶在短期浸水后表现出较上等耐久性与吸湿性。宋泽卓等[12]对干湿循环条件下黄原胶改良砂土的强度劣化特性进行了研究,但并未进一步提出改善水耐久性方案。

    对于海洋基础工程,海岸处的地基等建筑物在海水涨落、波浪侵蚀情况下会表现出较为显著的干湿循环特征,对钙质粉土地基等建筑物的稳定性产生影响。有学者在利用生物技术和材料加固海洋土方面开展了相关研究。Shi等[13]通过无侧限抗压试验和剪切波速测量试验研究了不同生物聚合物固化钙质粉土的加固效果。结果表明,经过生物聚合物加固的钙质粉土强度得到较大提升,UCS强度从未加固土的0.254 MPa提升至10.099 MPa,接近40倍。但经一次浸水试验后,除卡拉胶外的BPRS强度下降明显。改良BPRS的水稳性具有重要意义。

    植物脲酶诱导碳酸钙沉淀技术(EICP)利用大豆离心液中的脲酶诱导Ca2+沉淀,是一种环境友好型土体加固技术[14]。目前有学者将EICP与BPRS结合,进行土体加固。Arab等[15]通过无侧限抗压试验和3点弯曲试验研究了EICP和海藻酸钠联合胶结土的力学性能。结果表明,EICP和海藻酸钠联合固化生物砖的抗弯强度从未加固的0.76 MPa提升至2.09 MPa。Chiao- Lo等[16]的进一步研究表明共聚物结合EICP可显著地改善固化土的遇水强度,且当Ca2+的存在时,共聚物会发生非共价交联,EICP沉淀会均匀嵌入聚合物基体中形成双网络形态。Hamdan等[17]通过烧杯反应和水蒸汽试验研究了黄原胶氢化凝胶保水性对EICP沉淀的影响。其试验结果表明,黄原胶凝胶可以辅助EICP溶液保留反应基质,提高沉淀效率,减少氨气的产生。

    鉴于EICP和生物聚合物各自的特点,以及在土体加固方面的作用方式,本文利用EICP联合黄原胶(XG)对南海钙质粉土进行加固,通过无侧限压缩试验和剪切波速测试试验对不同黄原胶含量、不同EICP浓度和不同干湿循环次数下XG固化土和XG-EICP固化土的力学特性进行研究,并利用扫描电子显微镜对固化土的微观孔隙结构进行了分析,意在为海洋地基加固工程提供新的思路与方法。

    本试验中的钙质粉土取自中国南海地区。将所取土样干燥并利用筛机筛分去除大于1 mm的砂粒。通过标准击实试验得到的土壤最佳含水率为13.1%,最大干密度为1.94 g/cm3。土的塑限为15%,液限为27%[13]。对土样的物理性质进行测定,其性质如表 1所示。试验材料的XRD图谱如图 1(a)所示,经X射线衍射分析发现黄原胶样品中以碳氢氧化合物组成的高聚物为主,钙质粉土中以碳酸钙为主,存在少量有机物。细粒和粗粒的扫描电镜图像分别如图 1(b)(c)所示。

    表  1  土样物理性质[13]
    Table  1.  Physical properties of soil
    物理性质 数值
    D50 0.095 mm
    Cu 16.8
    Cc 0.05
    相对质量密度 2.77
    液限 27%
    塑限 15%
    最大干密度 1.94 g/cm3
    下载: 导出CSV 
    | 显示表格
    图  1  XRD图谱,钙质粉土形貌和钙质砂形貌
    Figure  1.  (a) Spectra of XRD; (b) Morphology of calcareous silt; (c) Morphology of calcareous sand

    本试验所用的生物聚合物为黄原胶。黄原胶是野油菜黄单胞菌发酵产生的胞外多糖,现被广泛应用于食品和化妆品行业。Freitas等[18]验证了黄原胶的安全性能。本试验的黄原胶购自河南春祥化工产品有限公司,白色-浅黄色粉末,吸水,在冷水中不易溶解,在温水中充分搅拌可形成透明黏稠状溶液。

    为了研究在不同干湿循环次数下XG固化土和XG-EICP固化土的物理力学劣化特征,对不同黄原胶含量和EICP浓度的固化土进行了质量损失及尺寸变化测量试验,剪切波速测试试验和无侧限压缩试验。试验前,烘干钙质粉土和黄原胶粉末,根据土的干密度和最优含水率进行设计,初始含水率为14%。黄原胶的含量分别为1%,3%,5%。

    EICP溶液通过黄豆粉末、尿素和无水氯化钙进行制备。称取黄豆粉末,按质量比1︰20加入去离子水,充分搅拌均匀。在离心速度和离心时间分别为2000转/min和10 min的条件下进行离心,取离心上清液,调节上清液脲酶电导率至100 μs/cm/min。尿素和无水氯化钙溶液体积比为1︰1,按1,1.5,2 mol/L配置尿素-氯化钙溶液。最后按体积比1︰1将离心上清液和尿素-氯化钙溶液进行充分混合,制得EICP溶液。本文中只对1%和5% XG含量试样联合EICP加固。

    试样采用模具直径为39.1 mm。称取180 g钙质粉土,并混合设计比例的黄原胶粉末,加入14%质量的水或EICP溶液拌匀后压样,制样干密度为1.843 g/cm3,试样高度为81 mm。为确保结果的可靠性,每组工况制备3个平行试样。将制备好的XG固化试样放入60℃的烘箱中干燥18 h,XG-EICP固化试样置于室温下反应12 h后进行相同的烘干操作。测量烘干后的试样的重量,尺寸和剪切波速。将试样置于静水中浸泡6 h,再放入60℃的烘箱中干燥18 h为1次干湿循环。本试验设置的干湿循环次数分别为0,4,8,12次。

    Shi等[19]的研究表明,弯曲元试验中采用Tip-tip的方式选取剪切波的传播时间具有较好的客观性。本试验利用一组弯曲元传感器,在激发幅值和频率分别为20 V和15 kHz的条件下进行剪切波速测量[20]。为了保证弯曲元件与试样的紧密接触,在仪器顶部放置一个500 g砝码,剪切波速计算原理如下:

    Vs=L(T2T1Δt)
    (1)

    式中:Vs为剪切波速;L为试样高度;T2为接收端波峰对应时间值;T1为发射端峰值对应时间;Δt为仪器系统误差所延迟时间。

    根据《土工试验方法标准》[21]进行无侧限压缩试验,轴向位移加载速率为2 mm/min。循环结束后,通过测量质量损失,尺寸变化、剪切波速和无侧限抗压强度变化来量化物理与力学劣化特征。试验流程如图 2所示。扫描电镜试验使用的是赛默飞世尔公司生产的Quattro S型环境电子显微镜。

    图  2  试样制备及试验流程
    Figure  2.  Sample preparation and test flow

    为进一步观察XG-EICP沉淀并验证XG与EICP混合反应机理,利用载玻片烧杯静置法收集沉淀进行SEM观察与EDS图谱分析。将100 mL烧杯、干净载玻片和透气封口膜高温灭菌,固定载玻片与烧杯中央,倒入50 mL EICP溶液和XG-EICP溶液于室温下反应12 h,如图 3所示。利用去离子水清洗载玻片,去除溶于水的有机物和Ca2+并在60℃的烘箱中烘干24 h后进行SEM测试和EDS图谱测试。

    图  3  溶液试验
    Figure  3.  Solution tests

    图 4(a)~(c)为不同干湿循环条件下试样受干湿循环侵蚀导致的质量损失情况。损失的质量与原试样质量之比为质量损失百分率。通过图 4(a)可以看出,随着干湿循环次数的增加,XG固化土的质量损失率逐渐增大。由于循环浸泡导致黄原胶体积膨胀析出和浸水溶解量增加,5%含量的XG固化土质量损失率大于低含量XG固化土。图 4(b)(c)表明,相同黄原胶浓度和干湿循环次数条件下,XG-EICP固化土的质量率损失低于XG固化土。对于1%含量的XG,随着EICP溶液浓度的提高,质量损失率不断降低。以上试验结果表明,在EICP作用下,黄原胶与Ca2+发生交联反应[22],溶解度变小,稳定性增强。同时,EICP生成的碳酸钙沉淀可能会胶结土体颗粒并填充孔隙,增强试样整体性,减少干湿循环对固化土的侵蚀。

    图  4  不同干湿循环次数下质量、尺寸与密度变化
    Figure  4.  Variation of mass, size and density under different wetting-drying cycles

    Cabalar等[23]通过膨胀和收缩试验表明,黄原胶含量增加会导致黏土出现膨胀现象。干湿循环作用下,黄原胶会随着水分浸润和干燥而不断地膨胀和收缩,因此研究不同循环次数下固化土的尺寸变化显得十分必要。图 4(d)为不同黄原胶浓度的XG固化土高度变化情况。由于XG含量的增大会导致整体弹性增加,拆模后试样会出现回弹,因此对于不同XG含量的试样初始高度存在一定差别。5%XG试样初始高度为83 mm而1%和3%XG试样的初始高度为81 mm左右。

    对于高度随干湿循环次数变化情况,1%和3% XG含量的试样高度变化幅度小于XG 5%的试样。当ECIP和XG共同加固试样时,如图 4(e)(f)所示,XG-EICP加固试样高度的降幅小于纯XG固化土试样,但EICP浓度的增加不会改变试样高度变化幅度和趋势,表明低浓度的EICP即可对XG加固试样起到稳定尺寸变化的作用。图 4(g)~(i)列出了不同EICP-XG固化土试样直径随干湿循环的变化情况。由图 4(g)可知,随着干湿循环次数的增加,XG加固试样的直径随干湿循环次数的增加而增加,其中,5%XG含量的试样增幅最为显著。这是由于,随着XG含量的增大,试样塑性增大,且初始试样较高,试样脱模后,XG固化土发生轴向收缩和横向膨胀的程度较大,导致直径增加较为显著。图 4(h)(i)对比了相同XG含量下不同EICP反应液浓度对径向变化的影响,可以看出,EICP对试样径向尺寸随干湿循环变化的影响较小。较高XG含量下,试样的径向膨胀性在干湿循环的作用下依旧显著。

    图 4(j)~(l)列出了不同工况下的试验密度随干湿循环变化情况。从图 4(j)可以得出,初始密度随着XG含量的增加而减小,黄原胶含量越高,固化土密度随干湿循环下降幅度越明显。这是由于黄原胶粉末密度小于钙质粉土的密度,随着XG含量和干湿循环次数的增加,固化试样体积膨胀愈显著,但试样质量在干湿侵蚀作用下不断降低,进而导致密度减小。图 4(k)4(l)表明,EICP反应会增大XG固化土的密度。EICP溶液会产生碳酸钙沉淀,2.0 mol/L浓度下最为显著。同时EICP会降低质量损失率,相同XG含量和干湿循环次数下,较高EICP浓度会得到较大试样质量,但对于5% XG试样,未经历干湿循环时2.0 mol/L试样的密度依旧最高,但随着干湿循环次数的增加,密度与其他加固试样趋于一致,表明较高XG浓度下,EICP对土体的作用发生了变化。

    图 5(a)为不同XG含量条件下,XG固化土在未经历干湿循环情况下的应力应变曲线。黄原胶含量会显著影响XG固化土的应力应变曲线。黄原胶含量越高,无侧限抗压强度越大,且达到峰值应力对应的应变值越大,表明试样延展性更好。5%XG试样的最大无侧限抗压强度为3.425 MPa,相较于未加固土的UCS强度,即0.208 MPa,增大约16倍。图 5(b)~(d)给出了不同XG含量试样在经历干湿循环作用下的强度变化情况,可以看到,随着干湿循环次数的增加,XG固化土都呈现不同程度的强度劣化现象。第4次干湿循环后,1%,3%,5% XG试样的最大UCS强度分别降低了0.342,0.79,0.47 MPa。强度降低幅度在4循环时不受XG含量的影响。同时,可以观察到,第四次循环时试样的峰值强度应变增大,表明试样延性增加。第8次循环后,1%,3%,5% XG含量的试样UCS峰值强度分别为1.13,1.95,2.01 MPa,相较未经干湿循环处理的试样,XG固化土强度降幅显著但仍高于未加固试样。峰值强度应变在1% XG含量下继续增大,但3%,5%情况下呈现减小的趋势。12次干湿循环后,试样强度低于初始值一半,峰值强度应变减小,试样脆性增加。

    图  5  不同干湿循环次数下固化土应力应变曲线
    Figure  5.  Stress-strain curves of treated soil under different wetting-drying cycles

    图 5(e)~(f)表示的是1%,5% XG与不同浓度EICP溶液联合加固试样未经历干湿循环的UCS强度情况。可以看出,联合EICP加固后,试样的强度得到明显提升。对于1% XG的试样,试样无侧限峰值强度从无EICP条件下的1.87 MPa提高至2.0 mol/L EICP条件下的3.77 MPa,强度提升一倍。同时,峰值强度应变向后移动,表明由于碳酸钙的生成,导致试样脆性增加。5% XG含量下EICP作用更加显著,2.0 mol/L对应的试样其峰值强度达到7.21 MPa,同样强度提升是未掺入EICP强度3.41 MPa接近一倍。以上可以看出,在未经历干湿循环时,EICP联合XG加固钙质粉土加固效果显著。图 5(g)~(i)为1% XG含量下不同EICP溶液浓度加固试验在干湿循环作用下的强度表现。对于1 mol/L EICP溶液,试样在经历干湿循环后强度出现劣化,但劣化的幅度较小,12次循环后UCS峰值强度为1.98 MPa,而未经干湿循环峰值强度为2.63 MPa。当EICP溶液浓度为1.5,2.0 mol/L时,试样强度随干湿循环劣化的程度更加显著。这是由于在Ca2+充足的条件下,XG会与Ca2+产生无机交联反应,生成强度较高且不溶于水的复合生物高分子聚合物沉淀[24-25]。其中,2.0 mol/L下试样12次干湿循环后强度降低近50%,但仍高于图 5(b)中未联合EICP加固的试样强度。此外,峰值强度应变随干湿循环次数增加不断增大,试样延性增加。图 5(j)~(l)为5% XG含量下的试样强度变化情况。可以看出,强度劣化随干湿循环次数的增加依旧显著,对比相同干湿循环次数但不同EICP浓度试样的强度发现,较高EICP浓度对应着较高的强度,表明在干湿循环作用下,EICP仍能提供试样的加固作用。

    为了更直观的展示XG和EICP的加固效果,图 6(a)~(c)绘制了加固试样的UCS峰值强度随干湿循环次数的变化情况。由于纯XG加固试样在经历12次干湿循环后的峰值强度较为接近,而初始强度差别较大,因此,高XG含量的试样峰值强度降低较快。对于1% XG含量的加固试样,峰值强度随干湿循环的的变化规律与EICP溶液浓度无明显的相关性,下降规律较为一致。然而对于5%含量的试样,4~8次干湿循环阶段纯在一个强度下降平缓阶段,在此阶段内,峰值强度变化较小,当干湿循环次数大于8时,强度降低幅度再次加大。

    图  6  不同干湿循环次数下峰值强度、割线模量、剪切波速与初始剪切模量变化
    Figure  6.  Variation of peak intensity, secant modulus, shear wave velocity and initial shear modulus under different wetting-drying cycles

    图 6(d)~(f)为试样的割线模量E50随干湿循环次数的变化情况。从图中可以看出,随着干湿循环次数的增加,1%含量的XG固化试样的平均刚度大于5%含量的XG固化试样,且随着EICP溶液浓度的提高,XG-EICP固化试样平均刚度增加。经历12次干湿循环后,XG固化试样和XG-EICP固化试样的平均刚度降幅较大,但相较于未处理的钙质粉土试样的刚度(图 7(d)仍具有较高的平均刚度。此外,对于XG-EICP联合加固,1.5,2.0 mol/L试样的割线模量大小相近,明显高于1.0 mol/L和无EICP参与的试样,且随干湿循环的变化规律也相近。

    图  7  含量固化土浸水侵蚀情况: 1%黄原胶
    Figure  7.  Immersion erosion of soil treated with xanthan gum of 1%

    图 6(g)~(i)为不同循环次数下,XG固化土和XG-EICP固化土的剪切波速Vs随着干湿循环次数的增加,试样的剪切波速逐渐减小。通过对比可以看出,低浓度的XG固化土具有更高的剪切波速,这可能是由于高浓度的黄原胶过量填充于钙质粉土的空隙中,导致钙质粉土颗粒间隔增大,黄原胶胶体阻隔剪切波传输,使得5%含量的XG固化土剪切波速减小。这与Shi等[13]和Ta等[26]的结论规律性一致。

    1%含量的XG固化土在经历12次干湿循环后的剪切波速减少率为7.6%,大于XG-EICP固化土。这是由于随着干湿循环次数的增加,固化土内部孔隙开始增大,钙质粉土颗粒间距增大,粒间连接减少。这与宋泽卓等[12]的结论一致。但是掺入EICP溶液后,EICP生成的碳酸钙沉淀使得XG-EICP固化土内部孔隙填充较紧密,浸水侵蚀得以缓解。由于EICP产生的碳酸钙沉淀降低了由黄原胶填充导致的剪切波传播阻隔,使得XG-EICP固化土在多次浸水循环条件下的剪切波速较纯XG固化土明显增大。同样,在5% XG含量,相同干湿循环次数下,XG-EICP固化土的剪切波速大多高于XG固化土,4次干湿循环次数后的波速情况认为是试验误差导致的。

    图 6(j)~(l)为不同循环次数下,XG固化土和XG-EICP固化土的初始剪切模量Gmax变化情况,Gmax可由剪切波速推到而出,具体公式可见文献[27]。随着干湿循环次数的增加,试样初始剪切模量Gmax逐渐减小,1%和5%含量的XG-EICP固化土变化规律与剪切波速类似。但是由于5%含量的XG试样由于体积膨胀和开裂导致密度减小,5%含量XG试样的初始剪切模量较其他试样更低。

    干湿循环侵蚀物理与力学测试结果表明,EICP的掺入为XG固化土提供了抗水侵蚀保护。但是Shi等[13]发现由于黄原胶具有水溶性,水分浸泡会导致高浓度的黄原胶膨胀析出或遇水溶解。在本试验中,通过循环浸水-烘干试验发现XG固化土会出现显著的裂纹。图 7中列出了1%含量纯XG固化土不同干湿循环次数下,试样裂缝开展情况。可以看出,对于纯XG加固试样,第8次循环时,试样上端部出现裂缝;第12次循环时,上端部裂缝加深,下端部出现裂纹。图 8中5%含量纯XG固化土浸水侵蚀情况表明,4次干湿循环时,试样两端部出现明显的膨胀,第8次干湿循环上下端部开始出现明显裂缝,并随干湿循环次数进一步加深。图 7中同样列出了1%黄原胶含量,不同浓度XG-EICP固化土随干湿循环次数增加,裂缝开展情况。可以看出,12次干湿循环几乎不会导致1%的XG-EICP试样产生裂缝。图 8中为5%黄原胶含量,不同浓度XG-EICP固化土的侵蚀情况。在第8次干湿循环时,1 mol/L的XG-EICP固化土端部出现鼓胀现象,同时伴随有裂缝开展,但是其深度和广度均小于5%含量XG固化土。

    图  8  含量固化土浸水侵蚀情况: 5%黄原胶
    Figure  8.  Immersion erosion of soil treated with xanthan gum of 5%

    图 9为5%黄原胶含量,不同浓度XG-EICP固化土下端部侵蚀情况图。从图 9可以看出,XG固化土的下端部随干湿循环次数增加,出现明显的坑洞与裂缝。第八次干湿循环时,2 mol/L固化土的端部裂纹明显少于其他试样,且随着EICP浓度的降低,裂纹愈发显著,表明EICP能够缓解试样开裂情况,但当干湿循环达到12次时,虽然较高EICP浓度下裂纹较少,但仍出现了显著的开裂情况。

    图  9  5%黄原胶含量固化土下端部浸水侵蚀情况
    Figure  9.  Bottom erosion of samples treated with xanthan gum of 5%

    图 10(a)为黄原胶粉末的电镜图,干燥状态下黄原胶粉末显现出单颗粒状形态。图 10(b)为无侧限压缩试验后裂缝位置处XG固化土的局部电镜图。钙质粉土颗粒间的黄原胶呈现纤维状结构,充斥于土壤孔隙间,连接不同土颗粒。部分颗粒之间由带状的胶体连接,起到加固土体的作用。图 10(c)为无侧限压缩试验后裂缝位置处XG-EICP固化土的局部电镜图,可以看出经过EICP溶液中Ca2+的无机交联改性,钙质粉土颗粒间的纤维状黄原胶减少,大部分黄原胶以凝胶状沉淀存在于土体孔隙中。黄原胶胶状沉淀中存在EICP产生的细颗粒碳酸钙沉淀,两者相互结合,这使得XG-EICP固化土强度和整体性更高。图 10(d)为连接体的局部放大图,可以看出凝胶的表面和内部存在大量EICP产生的细颗粒碳酸钙沉淀。

    图  10  XG和XG-EICP固化土扫描电镜图
    Figure  10.  SEM images of treated soil

    Bergmann等[24]将黄原胶胶体与钙离子进行室温下进行混合交联反应,生成了离子交换的黄原胶联合胶体。他们发现联合胶体的力学特性要显著优于纯黄原胶胶体。本文中,为了探究EICP和黄原胶之间的交联反应,对两者溶液混合后得沉淀物进行了微观分析。图 11(a)为EICP溶液沉淀的SEM图,EICP生成的碳酸钙晶体大小大致介于1~3 μm,这与Ahenkorah等[28]的结论一致。图 11(b)为EICP溶液沉淀的EDS图谱,试验结果表明:在喷金时间一致情况下,EICP溶液沉淀中Ca含量占全部元素的65.28%,C元素占全部元素的9.5%。这表明经过去离子水清洗过的EICP沉淀中Ca元素的含量远高于其他元素,与试验预期相符。图 11(c)为XG-EICP溶液沉淀的SEM图,可以看出掺入XG的EICP溶液中出现较大颗粒的XG-EICP沉淀。对其中一个颗粒进行EDS分析,如图 11(d)所示,可以看出其组成成分中C元素的含量明显高于纯EICP沉淀中C元素含量,且与图 2(a)中的XG的XRD图谱较为相近,由于拍摄电镜之前已将试样用去离子水反复冲洗,且在电镜图中未观测到XG胶状体,可以认为,所观测到的固体颗粒为XG-EICP混合溶液中XG与Ca2+反应生成的联合产物与EICP生成的碳酸钙混合体。根据图 3中的直观观测,这种混合体较纯XG胶体不易溶于水,对改善加固土体在干湿循环作用下的劣化具有正面作用。

    图  11  溶液沉淀扫描电镜图与EDS图谱
    Figure  11.  SEM images of solution precipitation and EDS spectrum

    通过无侧限压缩试验,剪切波速测试试验和干湿循环试验对不同黄原胶含量、不同EICP溶液浓度和不同干湿循环次数条件下的XG-EICP固化土的劣化特性进行研究,并利用扫描电子显微镜对无侧限压缩试验后裂缝位置处的XG-EICP固化土的微观孔隙结构进行了分析,由试验结果可以得到4点结论。

    (1)黄原胶与EICP技术联合应用对钙质粉土的耐水性具有显著改善。在相同次数的干湿循环条件下,XG-EICP固化土受浸泡侵蚀造成的质量损失明显小于XG固化土。胶体遇水膨胀造成的质量损失率随着黄原胶浓度的增加而增大。随着干湿循环次数的增加,XG-EICP固化土高度降低率,直径增长率和密度降低率均小于XG固化土。

    (2)黄原胶和EICP对钙质粉土的无侧限抗压强度有较大的影响。相较于未处理的钙质粉土,XG固化土强度有明显提升。掺入2.0 mol/L的EICP溶液时,最大抗压强度为7.21 MPa。

    (3)由于高浓度的黄原胶过量填充于钙质粉土的空隙,剪切波速随着XG固化土中黄原胶含量的增加而减小。随着干湿循环次数的增加,固化土内部孔隙增大,粒间连接减少,剪切波速降低。掺入EICP溶液后,碳酸钙沉淀使得固化土内部孔隙填充较紧密,剪切波速有所提升。经历12次干湿循环时,XG固化土和XG-EICP固化土的剪切波速降幅显著,高浓度XG固化土剪切波速和剪切模量的下降幅度大于低浓度XG固化土。

    (4)通过SEM和EDS分析可以得出,纯XG固化土经历干湿循环和烘干后其内部会出现针状黄原胶和带状黄原胶两种土颗粒连接形式,而XG-EICP固化土中针状黄原胶数量显著下降,XG和钙离子联合反应会生不易溶于水的沉淀产物,碳酸钙颗粒会嵌入进联合沉淀产物中,对改善加固土体在干湿循环作用下的劣化性质具有正面作用。

  • 图  1   考虑超载和纵坡效应的隧道开挖面稳定性分析模型

    Figure  1.   Stability analysis of model for tunnel face influenced by overload and longitudinal slope

    图  2   运动单元上限分析网格更新过程样例图(θ =0°, C/D =3, γD/c =0, ϕ =15°)

    Figure  2.   Sample charts of mesh updates with UB-RTME(θ =0°, C/D =3, γD/c =0, ϕ =15°)

    图  3   非结构化初始网格更新样例图(θ =0°, C/D =3, γD/c =0, ϕ =15°)

    Figure  3.   Sample charts of unstructured initial mesh update(θ =0°, C/D =3, γD/c =0, ϕ =15°)

    图  4   超载和纵坡影响下的隧道开挖面临界超载系数σs/c上限解曲线图

    Figure  4.   Curves of variation of upper bound solution of σs/c for tunnel face influenced by overload and longitudinal slope

    图  5   隧道开挖面滑移线网破坏模式及网格变形特征图(ϕ=15°, γD/c=2, C/D=3, θ=0°)

    Figure  5.   Slip line failure mechanisms and mesh deforming characteristics for tunnel face (ϕ=15°, γD/c=2, C/D=3, θ=0°)

    图  6   埋深比C/D影响下的隧道开挖面滑移线网破坏模式

    Figure  6.   Slip line failure mechanisms for tunnel face influenced by C/D

    图  7   重度系数γD/c和内摩擦角ϕ影响下的隧道开挖面滑移线网破坏模式

    Figure  7.   Slip line failure mechanisms of tunnel face influenced by γD/c and ϕ

    图  8   网格疏密影响下的隧道开挖面滑移线网破坏模式(θ=0°, ϕ=25°, C/D=2, γD/c=2)

    Figure  8.   Slip line failure mechanisms for tunnel face influenced by mesh density (θ=0°, ϕ=25°, C/D=2, γD/c=2)

    表  1   隧道开挖面稳定性分析计算参数

    Table  1   Dimensionless parameters for stability of tunnel face

    计算参数 参数取值
    C/D 1, 2, 3, 4
    γD/c 0, 1, 2, 3
    ϕ/(°) 10, 15, 20, 22.5, 25
    θ/(°) -10, 0, 10
    下载: 导出CSV

    表  2   运动单元上限分析网格更新操作

    Table  2   Mesh update operation of UB-RTME method

    序号 名称 网格原状 网格更新 操作说明
    1 移动节点 修改坐标调整节点位置
    2 新增节点 在公共边按长度比例新增节点,添加间断线
    3 单元缩减 合并单元缩减单元
    4 对角互换 相邻单元公共边互换
    5 单元加密 单元一分为四加密
    6 成层加密 沿指定方向成层加密
    下载: 导出CSV

    表  3   隧道开挖面破坏临界超载系数σs/c上限解

    Table  3   Upper bound solutions of critical load factor σs/c for failure of tunnel face

    C/D ϕ/(°) γD/c (θ=10°) γD/c (θ=0°) γD/c (θ=-10°)
    0 1 2 3 0 1 2 3 0 1 2 3
    10.0 7.47 5.50 3.50 1.47 8.17 6.25 4.29 2.29 9.01 7.16 5.25 3.30
    15.0 11.33 8.90 6.38 3.79 12.62 10.29 7.87 5.35 14.24 12.07 9.78 7.37
    1 20.0 19.22 15.93 12.49 8.87 21.93 18.89 15.65 12.20 25.60 22.93 20.01 16.84
    22.5 26.40 22.40 18.18 13.70 30.67 27.04 23.13 18.92 36.36 33.33 29.94 26.20
    25.0 38.06 32.99 27.58 21.80 44.90 40.42 35.52 30.19 54.60 51.11 47.08 42.54
    10.0 11.34 7.87 4.30 0.63 12.11 8.75 5.27 1.67 13.01 9.78 6.42 2.93
    15.0 19.45 14.92 10.16 5.13 21.15 16.89 12.35 7.53 23.16 19.25 15.03 10.50
    2 20.0 39.13 32.45 25.26 17.49 43.61 37.63 31.05 23.84 49.07 44.06 38.35 31.97
    22.5 60.05 51.35 41.89 31.56 67.97 60.44 52.02 42.67 77.78 71.94 65.06 57.17
    25.0 98.40 86.30 73.02 58.41 113.40 103.36 91.95 79.09 132.34 125.40 116.84 106.70
    10.0 14.43 9.44 4.27 -1.14 15.24 10.41 5.37 0.10 16.19 11.57 6.69 1.58
    15.0 26.89 20.20 13.06 5.36 28.91 22.66 15.81 8.56 31.22 25.53 19.27 12.44
    3 20.0 60.84 50.41 38.98 26.40 66.93 57.75 47.45 35.97 74.10 66.61 57.85 47.87
    22.5 100.59 86.36 70.61 53.12 112.27 100.30 86.61 71.11 126.21 117.33 106.52 93.82
    25.0 179.87 158.75 135.21 108.88 204.11 187.32 167.74 145.25 233.42 222.73 208.87 191.93
    10.0 17.06 10.55 3.72 -3.50 17.88 11.60 4.96 -2.06 18.90 12.86 6.43 -0.37
    15.0 33.88 24.97 15.36 4.84 36.15 27.85 18.79 8.84 38.77 31.24 22.87 13.62
    4 20.0 83.79 69.34 53.33 35.48 91.38 78.83 64.54 48.40 100.30 90.21 78.21 64.33
    22.5 146.57 126.17 103.33 77.66 161.81 144.98 125.43 103.01 180.38 168.30 153.21 135.18
    25.0 279.81 248.03 212.19 171.67 313.31 288.78 259.67 228.22 355.36 340.68 320.89 296.10
    下载: 导出CSV

    表  4   隧道开挖面破坏临界重度系数γcrD/c对比数据

    Table  4   Comparison of critical gravity load factor γcrD/c for tunnel face

    C/D ϕ/(°) 文献[12] 变形单元 文献[13] 运动单元 Optum G2结果 本文解
    γcrD/c σs/c
    1 10 4.00 4.04 3.97 4.00 0.255
    1 15 4.89 4.88 4.77 4.90 0.264
    1 20 5.90 5.89 5.73 5.99 0.247
    2 10 3.48 3.41 3.37 3.40 0.199
    2 15 4.39 4.36 4.28 4.40 0.180
    2 20 5.62 5.55 5.43 5.67 0.137
    3 10 3.05 2.99 2.96 3.00 0.102
    3 15 4.08 4.01 3.95 4.05 0.112
    3 20 5.37 5.34 5.24 5.45 0.193
    下载: 导出CSV

    表  5   隧道开挖面相对破坏范围S/D2计算结果

    Table  5   Computated results for failure area S/D2 of tunnel face

    C/D ϕ/(°) γD/c (θ=10°) γD/c (θ=0°) γD/c (θ=-10°)
    0 1 2 3 0 1 2 3 0 1 2 3
    1 10 2.17 2.05 1.93 1.81 2.36 2.21 2.06 1.93 2.54 2.35 2.19 2.04
    15 2.16 2.03 1.90 1.76 2.35 2.19 2.04 1.89 2.59 2.39 2.21 2.05
    20 2.15 2.04 1.92 1.79 2.35 2.21 2.08 1.94 2.57 2.41 2.25 2.11
    25 2.17 2.08 1.99 1.88 2.37 2.26 2.16 2.05 2.60 2.47 2.36 2.24
    2 10 5.70 5.23 4.78 4.29 6.17 5.60 5.08 4.56 6.73 6.05 5.47 4.90
    15 5.74 5.30 4.85 4.37 6.21 5.69 5.20 4.67 6.76 6.15 5.62 5.08
    20 5.78 5.44 5.08 4.69 6.24 5.84 5.45 5.05 6.81 6.35 5.92 5.52
    25 5.86 5.63 5.38 5.11 6.33 6.06 5.79 5.52 6.89 6.59 6.30 6.03
    3 10 10.97 9.85 8.75 7.47 11.77 10.46 9.26 7.97 12.68 11.15 9.87 8.57
    15 11.08 10.11 9.12 7.97 11.85 10.74 9.70 8.60 12.83 11.57 10.41 9.32
    20 11.18 10.50 9.79 8.99 11.96 11.20 10.44 9.66 12.94 12.05 11.24 10.46
    25 11.37 10.94 10.50 10.01 12.15 11.69 11.23 10.75 13.15 12.61 12.11 11.64
    4 10 17.99 15.89 13.79 10.86 19.22 16.78 14.51 11.77 20.53 17.75 15.40 12.81
    15 18.20 16.48 14.71 12.46 19.29 17.41 15.55 13.52 20.83 18.64 16.64 14.71
    20 18.40 17.27 16.08 14.71 19.62 18.35 17.11 15.81 21.01 19.59 18.27 17.02
    25 18.72 18.07 17.40 16.63 19.93 19.19 18.48 18.48 21.38 20.56 19.81 19.09
    下载: 导出CSV
  • [1] 严松宏, 李国良. 隧道与地下工程荷载计算的研究现状[J]. 兰州交通大学学报, 2020, 39(2): 1-7. doi: 10.3969/j.issn.1001-4373.2020.02.001

    YAN Songhong, LI Guoliang. Research status of the load calculation of tunnel and underground engineering[J]. Journal of Lanzhou Jiaotong University, 2020, 39(2): 1-7. (in Chinese) doi: 10.3969/j.issn.1001-4373.2020.02.001

    [2] 李修磊, 李金凤, 杨超. 黏土地层浅埋盾构隧道开挖面三维稳定性上限分析[J]. 铁道学报, 2021, 43(4): 166-174. doi: 10.3969/j.issn.1001-8360.2021.04.021

    LI Xiulei, LI Jinfeng, YANG Chao. Upper-bound solutions for three-dimensional excavation face stability of shallow shield tunnel in cohesive-frictional soils[J]. Journal of the China Railway Society, 2021, 43(4): 166-174. (in Chinese) doi: 10.3969/j.issn.1001-8360.2021.04.021

    [3] 张明告, 周顺华, 黄大维, 等. 地表超载对地铁盾构隧道的影响分析[J]. 岩土力学, 2016, 37(8): 2271-2278.

    ZHANG Minggao, ZHOU Shunhua, HUANG Dawei, et al. Analysis of influence of surface surcharge on subway shield tunnel under[J]. Rock and Soil Mechanics, 2016, 37(8): 2271-2278. (in Chinese)

    [4] 黄大维, 周顺华, 冯青松, 等. 地表均布超载作用下盾构隧道上覆土层竖向土压力转移分析[J]. 岩土力学, 2019, 40(6): 2213-2220.

    HUANG Dawei, ZHOU Shunhua, FENG Qingsong, et al. Analysis for vertical earth pressure transference on overlaying soils of shield tunnel under uniform surface surcharge[J]. Rock and Soil Mechanics, 2019, 40(6): 2213-2220. (in Chinese)

    [5] 曾明生, 王文法, 康海波, 等. 基于不同破坏模式的盾构隧道掌子面支护力上限分析[J]. 铁道科学与工程学报, 2022, 19(8): 2360-2368.

    ZENG Mingsheng, WANG Wenfa, KANG Haibo, et al. Upper limit analysis of support force on shield tunnel face based on different failure modes[J]. Journal of Railway Science and Engineering, 2022, 19(8): 2360-2368. (in Chinese)

    [6]

    HUANG M S, SONG C X. Upper-bound stability analysis of a plane strain heading in non-homogeneous clay[J]. Tunnelling and Underground Space Technology, 2013, 38: 213-223. doi: 10.1016/j.tust.2013.07.012

    [7] 黄茂松, 宋春霞, 王浩然. 基于上限定理的软土隧道开挖面稳定性分析[J]. 防灾减灾工程学报, 2014, 34(3): 330-335.

    HUANG Maosong, SONG Chunxia, WANG Haoran. Upper bound limit analysis for face stability of shield tunnel[J]. Journal of Disaster Prevention and Mitigation Engineering, 2014, 34(3): 330-335. (in Chinese)

    [8] 黄茂松, 宋春霞, 吕玺琳. 非均质黏土地基隧道环向开挖面稳定上限分析[J]. 岩土工程学报, 2013, 35(8): 1504-1512. http://cge.nhri.cn/article/id/15259

    HUANG Maosong, SONG Chunxia, LÜ Xilin. Upper bound analysis for stability of a circular tunnel in heterogeneous clay[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(8): 1504-1512. (in Chinese) http://cge.nhri.cn/article/id/15259

    [9]

    SURAPARB K, BOONCHAI U. Design equation for stability of a circular tunnel in anisotropic and heterogeneous clay[J]. Underground Space, 2022, 7(1): 76-93. doi: 10.1016/j.undsp.2021.05.003

    [10] 阳军生, 张箭, 杨峰. 浅埋隧道掌子面稳定性二维自适应上限有限元分析[J]. 岩土力学, 2015, 36(1): 257-264.

    YANG Junsheng, ZHANG Jian, YANG Feng. Stability analysis of shallow tunnel face using two-dimensional finite element upper bound solution with mesh adaptation[J]. Rock and Soil Mechanics, 2015, 36(1): 257-264. (in Chinese)

    [11]

    YANG F, ZHANG J, YANG J S, et al. Stability analysis of unlined elliptical tunnel using finite element upper-bound method with rigid translatory moving elements[J]. Tunnelling and Underground Space Technology Incorporating Trenchless Technology Research, 2015, 50: 13-22.

    [12]

    YANG F, ZHANG J, ZHAO L H, et al. Upper-bound finite element analysis of stability of tunnel face subjected to surcharge loading in cohesive-frictional soil[J]. KSCE Journal of Civil Engineering, 2016, 20(6): 2270-2279. doi: 10.1007/s12205-015-0067-z

    [13] 孙雁军, 阳军生, 罗静静, 等. 隧道工作面稳定性与滑移线网破坏模式研究[J]. 岩土工程学报, 2019, 41(7): 1374-1380. doi: 10.11779/CJGE201907024

    SUN Yanjun, YANG Junsheng, LUO Jingjing, et al. Stability and mesh-like collapse mechanism of tunnel face[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(7): 1374-1380. (in Chinese) doi: 10.11779/CJGE201907024

    [14] 杨峰, 何诗华, 吴遥杰, 等. 非均质黏土地层隧道开挖面稳定运动单元上限有限元分析[J]. 岩土力学, 2020, 41(4): 1412-1419, 1436.

    YANG Feng, HE Shihua, WU Yaojie, et al. Tunnel face stability analysis by the upper-bound finite element method with rigid translatory moving element in heterogeneous clay[J]. Rock and Soil Mechanics, 2020, 41(4): 1412-1419, 1436. (in Chinese)

    [15]

    WENG X L, SUN Y F, YAN B H, et al. Centrifuge testing and numerical modeling of tunnel face stability considering longitudinal slope angle and steady state seepage in soft clay[J]. Tunnelling and Underground Space Technology, 2020, 101: 103406. doi: 10.1016/j.tust.2020.103406

    [16]

    CHENG C, JIA P J, ZHAO W, et al. Experimental and analytical study of shield tunnel face in dense sand strata considering different longitudinal inclination[J]. Tunnelling and Underground Space Technology, 2021, 113: 103950. doi: 10.1016/j.tust.2021.103950

    [17] 马辉, 王飞, 刘泽挂, 等. 深埋特长大坡度斜井变形规律与影响因素分析[J]. 地下空间与工程学报, 2020, 16(增刊2): 950-956.

    MA Hui, WANG Fei, LIU Zegua, et al. Analysis of deformation law and influencing factors of deep-buried extra-long inclined shaft with large slope[J]. Chinese Journal of Underground Space and Engineering, 2020, 16(S2): 950-956. (in Chinese)

    [18]

    SHI X M, LIU B G, TANNANT D, et al. Influence of consolidation settlement on the stability of inclined TBM tunnels in a coal mine[J]. Tunnelling and Underground Space Technology, 2017, 69: 64-71. doi: 10.1016/j.tust.2017.06.013

    [19] 雷明锋, 彭立敏, 施成华, 等. 迎坡条件下盾构隧道开挖面极限支护力计算与分析[J]. 岩土工程学报, 2010, 32(3): 488-492. http://cge.nhri.cn/article/id/12413

    LEI Mingfeng, PENG Limin, SHI Chenghua, et al. Calculation and analysis of limit support force of shield tunnel excavation face under facing-slope conditions[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(3): 488-492. (in Chinese) http://cge.nhri.cn/article/id/12413

    [20] 周峻, 杨子松, 彭芳乐. 上坡条件下盾构开挖面极限支护压力研究[J]. 地下空间与工程学报, 2011, 7(5): 914-918. doi: 10.3969/j.issn.1673-0836.2011.05.016

    ZHOU Jun, YANG Zisong, PENG Fangle. A study on the support pressure limit of the excavation face of shield tunnel under upslope conditions[J]. Chinese Journal of Underground Space and Engineering, 2011, 7(5): 914-918. (in Chinese) doi: 10.3969/j.issn.1673-0836.2011.05.016

    [21] 赵智涛, 曹伍富, 王霆. 隧道爬坡开挖与水平开挖掌子面稳定性对比分析[J]. 现代隧道技术, 2018, 55(6): 94-100.

    ZHAO Zhitao, CAO Wufu, WANG Ting. Contrastive analysis of working face stability in cases of climbing excavation and horizontal excavation[J]. Modern Tunnelling Technology, 2018, 55(6): 94-100. (in Chinese)

    [22]

    ZHAO L H, LI D J, LI L, et al. Three-dimensional stability analysis of a longitudinally inclined shallow tunnel face[J]. Computers and Geotechnics, 2017, 87: 32-48. doi: 10.1016/j.compgeo.2017.01.015

    [23] 程诚, 赵文, 王迎超, 等. 密砂地层盾构隧道纵坡开挖面稳定性理论分析[J]. 中国公路学报, 2023, 36(4): 157-168. doi: 10.3969/j.issn.1001-7372.2023.04.014

    CHENG Cheng, ZHAO Wen, WANG Yingchao, et al. Theoretical analysis of longitudinally-inclined shield tunnel face stability in dense sand stratum[J]. China Journal of Highway and Transport, 2023, 36(4): 157-168. (in Chinese) doi: 10.3969/j.issn.1001-7372.2023.04.014

    [24] 张箭, 戚瑞宇, 宗晶瑶, 等. 盾构隧道环向开挖面破坏机制及剪胀效应研究[J]. 岩土力学, 2022, 43(7): 1833-1844.

    ZHANG Jian, QI Ruiyu, ZONG Jingyao, et al. Failure mechanism of shield tunnel circumferential excavation face and the influence of the dilatancy effect on the tunnel stability[J]. Rock and Soil Mechanics, 2022, 43(7): 1833-1844. (in Chinese)

图(8)  /  表(5)
计量
  • 文章访问数:  317
  • HTML全文浏览量:  64
  • PDF下载量:  62
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-07-31
  • 网络出版日期:  2024-04-23
  • 刊出日期:  2024-10-31

目录

/

返回文章
返回