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

面膜堆石坝不同品种土工膜力学特性

张宪雷, 马仲阳, 吴云云

张宪雷, 马仲阳, 吴云云. 面膜堆石坝不同品种土工膜力学特性[J]. 岩土工程学报, 2023, 45(5): 940-952. DOI: 10.11779/CJGE20220151
引用本文: 张宪雷, 马仲阳, 吴云云. 面膜堆石坝不同品种土工膜力学特性[J]. 岩土工程学报, 2023, 45(5): 940-952. DOI: 10.11779/CJGE20220151
ZHANG Xianlei, MA Zhongyang, WU Yunyun. Mechanical properties of different geomembranes in membrane-faced rockfill dam[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(5): 940-952. DOI: 10.11779/CJGE20220151
Citation: ZHANG Xianlei, MA Zhongyang, WU Yunyun. Mechanical properties of different geomembranes in membrane-faced rockfill dam[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(5): 940-952. DOI: 10.11779/CJGE20220151

面膜堆石坝不同品种土工膜力学特性  English Version

基金项目: 

国家自然科学基金项目 51609087

国家自然科学基金项目 51709114

国家重点研发计划基金项目 2018YFC0406803

详细信息
    作者简介:

    张宪雷(1984—),男,博士,副教授,主要从事膜防渗结构方面研究。E-mail: zhangxianlei@ncwu.edu.cn

  • 中图分类号: TV641.4

Mechanical properties of different geomembranes in membrane-faced rockfill dam

  • 摘要: 针对面膜堆石坝主体防渗材料常用的几个品种土工膜力学性能,展开了土工膜单向拉伸、液胀、耐局部划痕损伤和适应颗粒垫层变形能力的试验性研究。试验分析成果表明:1.0 mm厚的膜材中PVC-HX膜的杨氏模量、应力应变线弹性区间和适应变形能力方面具有一定优越性;划痕长度、方向和深度占膜厚比值相同的条件下,TPO膜耐局部划痕损伤能力较强,PVC-HX膜仍展现良好的耐损伤极限强度和极限延伸率,能够满足膜防渗结构力学性能要求;厚度大于0.8 mm的HDPE/PE膜刚度偏大,适应变形和耐局部划痕损伤能力偏低,不适用于高水头压力下多孔隙介质垫层的膜防渗结构。试验成果分析表明膜厚大于1.0 mm的PVC-HX膜和TPO膜在高面膜堆石坝膜防渗结构中应用具有潜在优势,工程设计厚度不大于0.8 mm的土工膜可优先选择HDPE/PE膜。
    Abstract: Aiming at the mechanical properties of several geomembranes (GMs) commonly used as the main impervious materials of membrane-faced rockfill dam (MFRD), the uniaxial tensile tests, liquid expansion tests, local scratch resistance tests and capability to accommodate deformation tests are carried out. The test results show that compared to other GMs, the PVC-HX with the thickness of 1.0 mm is superior in the Young's modulus, linear elastic range of stress-strain and capability to accommodate deformation. Under the same scratch length, direction and ratio of depth to thickness, the TPO has strong resistance to the local scratch damage. However, the PVC-HX behaves satisfactorily in the ultimate strength against damage and the ultimate elongation, which can meet the mechanical performance requirements of impervious structures of the GMs. The HDPE/PE with thickness greater than 0.8 mm has high stiffness, low capability to accommodate deformation and local scratch damage, therefore, it is not suitable for the impervious structures of the GMs under high pressure and porous medium cushion. The comprehensive results indicate that the PVC-HX and TPO with thickness greater than 1.0 mm have potential advantages in the impervious structures of high MFRDs, and the HDPE/PE is recommended when the thickness is less than 0.8 mm required in engineering design.
  • 边坡、隧道、路基和地下工程的失稳垮塌机制复杂,岩土体是这些典型岩土工程灾变的主要赋存介质,深入掌握土的本构特性对工程安全分析十分必要。Rendulic[1]和Henkel[2]针对饱和正常固结黏土和具有相同先期固结压力的超固结土的三轴试验研究发现,不管是固结排水试验还是固结不排水试验,破坏前土体有效应力和孔隙比间存在唯一性关系,土破坏时孔隙比和强度间存在唯一性关系,且结果与试验类型及应力路径等无关。这就意味着土的抗剪强度与有效围压和土的孔隙比(反映了土的密实程度)之间存在着唯一性关系。

    在上述研究基础上,Schofield等[3]、Muir Wood等[4]提出了土的剑桥模型(Cambridge Model Cam- clay),形成了经典的临界状态土力学理论,解释了黏土的临界土力学行为。经典剑桥模型是基于黏土试验数据总结提炼的,难以描述砂土在临界状态下的应力软化和剪胀行为。砂土的经典剪胀模型可反映土密实状态和围压等物理量对无黏性土摩擦强度的影响[5]。为完善临界土力学理论描述无黏性土,Jefferies[6]引入了状态量ψψ=eecse为土的孔隙比,ecs为土临界状态下的孔隙比),假定应力增量与应变增量方向相同、应变增量方向与屈服面垂直,提出了一种砂土的修正剑桥模型,即Nor-Sand模型。

    不管是黏土还是砂土,根据临界状态理论,土体能够在临界状态下发生无限变形,应力和孔隙比恒定,在常规三轴应力条件下,其数学表达式可表示为

    pεq=0 qεq=0 εvεq=0 (1)

    式中:p=(σ1+2σ3/3q=σ1σ3εq为剪应变,可表示为εq=2(ε1ε3)/3。且当ψ0时,εq。在p-q-v三维空间坐标系下(其中v=1+e),临界土力学理论定义了物态边界面(罗斯柯面,State Boundary Surface),即以NCL(normal consolidation line,正常固结线)和CSL(critical state line,临界状态线)为边界的空间曲面,可定量化描述临界状态下土平均应力p、抗剪强度q与孔隙比e之间的关系。除了土临界状态理论所描述的土的抗剪强度与有效围压、密实度间的相关性,大量土力学试验发现,加载速率与土的不排水强度存在正相关关系(如图1[7],一般的三轴试验会持续数十分钟或数小时,试验发现,低加载速率下测得的土体强度比高加载速率下小得多,最大可相差50%。尽管这种加载效应得到了广泛认识,其物理机理尚未明确。

    图  1  应变率对不排水强度的影响[7],其中su为三轴压缩试验的不排水剪切强度,以剪切速率1%/h的强度为参考值,˙ϵ为剪切速率
    Figure  1.  Influences of strain rate on undrained strength [7] (su represents the undrained shear strength from triaxial compression tests, normalized by that at a shear rate of 1%/hour, and ˙ϵ is shear rate)

    采用堆石料和无黏性土在不同围压下开展三轴试验,试样开始发生应变硬化达到峰值强度,低围压下随后转为应变软化发展到临界状态[8-9]。总结试验数据发现[7],峰值摩擦角和临界状态摩擦角随着围压增大而减小,临界状态的摩擦角符合下列公式:

    φcr=φ0Δφlg(σ3pa) (2)

    式中:pa为与围压σ3量纲相同的大气压力;φ0Δφ为试验参数。

    近百年土力学领域大量科研工作,揭示了土的内摩擦角(摩擦特性)与围压(外荷载大小)、孔隙比(体变特征)、加载速率(变形速率)等因素存在相关性,基于试验结果建立了各自变量相关性的经验公式。这种相关性需满足普适性的物理定律,即应该存在深层次的颗粒动力统一机制。考虑到颗粒物质是自然土体的实际载体,本研究从颗粒物理的视角出发,基于颗粒力学理论提出准静态颗粒惯性数的概念,基于大量三轴试验数据,建立临界状态下砂土摩擦特性、体变特征、剪切速率间的标度律,发展土力学的临界状态理论。从另一角度看,以大量既有三轴试验数据为支撑建立准静态颗粒的标度律,可为经典流态μ(I)颗粒流变准则的准静态修正以及颗粒材料统一准静态-流态本构法则的构建提供可行方案,也能够为深入剖析岩土灾变机理创造条件。

    颗粒体的快速平板剪切试验是定量研究颗粒体在快速剪切状态下应力-应变行为的有效手段,Baran等[10]开展了三维颗粒体剪切试验,试验中对顶板施加恒定的压力P,在两板之间的颗粒(颗粒密度为ρs、粒径为d)以固定剪切速率˙γ受剪,上板可以垂向自由移动,试验中测出极限状态时恒定的体积分数ϕ和剪切应力τ。进而可以计算广义摩擦系数(friction coefficient)μ=τ/P。试验发现在稳态变形状态(steady-state,对应土的临界状态),可通过量纲分析确定颗粒的本构关系主要特征[11],该体系的性质受颗粒重分布微观时间尺度tmicro与剪切速率定义的宏观时间尺度的比值控制。颗粒体系从一个稳定的定向排列变形到下一个稳态(如图2),过程满足牛顿第二定律,即

    md2zdt2=Fz (3)
    图  2  颗粒受剪变形分析概化模型,惯性数为微观颗粒重分布时间尺度(左)和宏观变形时间尺度(右)比值
    Figure  2.  Schematic of granular system under shear (inertia number is defined as time scale ratio between microscopic rearrangement (left) and macroscopic deformation rate (right))

    考虑到mρsd3d2z/dt2d/t2以及FzPd2,代入式(3)化简后,可得围压约束条件下颗粒重新排列时间t=tmicro,即

    tmicro = d/P/ρs (4)

    又由相邻错动的两层颗粒相对变形速度Δu=˙γd,相邻颗粒间隙的距离为单倍粒径d(如图2),可得到特定颗粒从某粒间隙移动到相邻间隙的宏观时间为

    tmacro=dΔu=d˙γd=˙γ1 (5)

    式中:˙γ为剪应变速率。将微观时间尺度与宏观变形时间尺度作比,可综合表征颗粒变形的剧烈程度,可定义无量纲数I,又称为惯性数[12-13]

    I=tmicrotmacro=|˙γ|dP/ρs (6)

    根据以上量纲分析的结果,并结合斜面上颗粒流试验[14]、三维颗粒体剪切数值试验[15]、二维颗粒体数值试验[16]、环剪数值试验[17],GDR[18]、Forterre等[19]定义了流态颗粒的广义摩擦系数μ、颗粒体积分数ϕ与惯性数I之间的相关性,拟合这些数据点可得到被广泛认可的μ(I)ϕ(I)颗粒流变准则:

    μ(I)=μs+μ2μsI0/I0II+1 (7)
    ϕ=ϕmax+(ϕminϕmax)I (8)

    式中:μsμ2分别为广义摩擦系数的临界值(critical value)和极限值(limiting value),I0为试验参数。μsμ2I0与自然土颗粒的尺寸参数(长度、宽度、面积、周长)和形貌特征(长宽比、磨圆度、粗糙度)相关。试验参数ϕmaxϕmin分别为颗粒体积分数的最大和最小值,对应颗粒体系最密实和最疏松两种状态。

    目前关于μϕI间相关性在数值计算、室内试验等方面的探究多是针对颗粒流态(I>10-3)状态下开展的。μ(I)和ϕ(I)颗粒流变准则揭示了颗粒摩擦强度μ,颗粒体积分数ϕ与外力作用P、剪应变速率˙γ、颗粒材料特性(dρs)的相关性。虽然μ(I)ϕ(I)流变准则对流态(I>10-3)颗粒动力行为的刻画已得到广泛认可,但在分析准静态变形颗粒(I<10-3)时失效。为描述解析颗粒物质的准静态力学行为,Kamrin等[20]、Henann等[21]在塑性力学框架下,引入颗粒关联长度的概念修正塑性势函数,提出非局部流变准则;Bouzid等[22-23]通过I的空间梯度项表示其耗散作用,直接对摩擦系数进行修正,即非局部流变准则。这些修正方法较为复杂,与唯象的土力学理论存在应用鸿沟。

    考虑到经典土力学三轴压缩试验中土样变形多处于准静态(I<10-3),近百年土力学研究积累了大量试验的三轴试验数据,本研究基于这些试验数据发展了颗粒惯性数概念,分析初始孔隙比、剪切速率、围压和粒径、中主应力等因素对颗粒摩擦特性和体变规律的影响,探讨颗粒准静态变形动力特征的物理标度律。

    土的三轴试验是土力学中最常用的经典试验方法,可完整地反映试样受力变形直到破坏的全过程,用作土体强度、应力应变关系分析。用三轴试验开展土的受力分析,与经典流态颗粒μ(I),ϕ(I)流变准则受力变形状态有所差别,惯性数有所不同,本节着重分析三轴试验中的颗粒惯性数。

    经典μ(I),ϕ(I)颗粒流变准则是基于单剪物理和数值试验方法取得的,试样处于二维纯剪状态,即颗粒以均一的剪切速率发生剪切变形。为进一步描述复杂三维颗粒动力特征,假设颗粒在运动过程中保持体积不变(即不可压缩),基于偏应力张量与应变率张量共轴性假定(即线性相关),Jop等[24]将标量形式的μ(I)局部流变准则扩展成张量形式,推广了该准则在连续介质模型中的应用:

    σ=Pδ+τ, (9)
    τ=η(|˙γ|,P)˙γ (10)
    η(|˙γ|,P)=μ(I)P|˙γ| (11)

    式中:στ˙γ分别为应力张量、偏应力张量和应变率张量;δ为单位张量;|˙γ|=(12˙γij˙γji)12为应变率第二不变量。

    在常规三轴试验过程中,施加围压后通过施加轴向力产生偏应力使得土样发生变形破坏。相对密实的土样会出现明显的剪切带,此时变形主要集中在剪切带区域,即土样发生局部剪切破坏。相对疏松土样的变形则比较均匀,无剪切带出现。需要说明的是,本文所取三轴试验数据均来自于局部剪切型破坏。对于三轴试验中的三维受剪土样建立物理模型,考虑到试样横向变形相对轴向较小,式(6)中定义应变率第二不变量主要受轴向变形速率控制,即为试样的剪切速率v与圆柱试样轴向尺寸H之比(即˙ε=v/H)。此外,取试验围压定义P(即P=σ3),天然土样的颗粒粒径d取平均粒径d50。根据式(6),可以得到三轴试验中砂土的惯性数为

    I=|˙γ|dP/ρs = |˙ε|d50σ3/ρs (12)

    土力学和颗粒物理领域常用内摩擦角或摩擦系数来描述土和颗粒的内摩擦特性。而通过试验得到的内摩擦角与颗粒材料摩擦系数间的关系与试验方法息息相关。例如,在重力驱动颗粒流动的斜槽试验中[14],摩擦系数可定义为临界滑动状态下平面倾斜角度的正切值,这个临界角度类似于准静态砂土休止角,可以用于评价砂土的摩擦属性。

    临界状态理论[3]定义了在大变形条件下,材料达到一种临界状态(如图3所示),其特征为体积分数和剪应力(对应摩擦特性)均保持不变。由于试验数据均来源于三轴试验,因此本文中的摩擦系数均采用临界状态内摩擦角φcs进行计算,即μ=tanφcs。此外,为与颗粒物理描述表达保持一致,本研究将土力学里面的孔隙比e通过下式换算为颗粒体积分数ϕ

    ϕ=VsVs+Vv=11+e (13)
    图  3  密实和疏松土受剪变形的临界状态
    Figure  3.  The critical state of dense and loose soil under shear

    式中:Vs为土体中土颗粒的体积;Vv为土体中孔隙的体积。

    本研究选取代表性无黏性颗粒材料的三轴试验数据进行分析。Alshibli等[25]采用长宽比、磨圆度、粗糙度等形貌特征差别较大的3种砂在4种围压下进行了三轴排水压缩试验,部分试验物理指标如表1所示。每种试验材料设置了3种不同的密实度(松散L、中密M、密实D),中密和密实试样利用Al-Shibli等[26]的空气淋滤装置制备,松散的样品使用ASTM D4254(ASTM 2016b)等程序制备(手拿一个带有13 mm喷嘴的漏斗,从25 mm的落点高度沉积砂子)。圆柱形试验试样直径为71.1 mm,高度为142.2 mm,偏移应力以0.5 mm/min的恒定位移速率施加。虽然3种砂的Gs=2.65和d50=0.36 mm参数一致,但是电镜图像发现3种材料的长宽比、磨圆度、粗糙度等差别很大。Ottawa砂、玻璃砂和Columbia grout 砂的μ-ln(I)关系如图4所示。

    表  1  试验材料的物理指标
    Table  1.  Physical properties of test materials
    材料 Gs d50/mm 密实状态 e
    F-35 Ottawa 砂 2.65 0.36 疏松 0.733
    中等 0.614
    致密 0.494
    #1玻璃砂 2.65 0.36 疏松 0.928
    中等 0.778
    致密 0.638
    GS#40 Columbia grout砂 2.65 0.36 疏松 0.921
    中等 0.762
    致密 0.639
    下载: 导出CSV 
    | 显示表格

    图 4可以看出,3种砂的临界状态摩擦系数μ均与惯性数和颗粒的密实度相关。对于不同类型的土,土越密实摩擦系数越大。如果将采用颗粒物理中初始颗粒体积分数ϕ0来描述颗粒的密实程度,可以得到颗粒准静态变形摩擦特性μ、初始密实状态ϕ0和颗粒惯性数I间的标度律,符合线性关系,如图5示。

    图  4  颗粒松密程度对摩擦系数μ与ln(I)关系的影响
    Figure  4.  Definition of internal friction angle at critical state
    图  5  不同密实状态土摩擦系数μ与准静态惯性数Q标度律
    Figure  5.  Scaling laws between friction coefficient μ and Q for soils at different density states

    不同于流态颗粒运动状态下摩擦系数μ与惯性数I之间的一一对应关系,颗粒在高围压、准静态下,摩擦系数μϕ0[ln(I)+α]之间存在线性一一对应关系。其中,对于3种试验材料,α的数值均为α=30,这可能是因为3种材料的比重和中值粒径等参数相近。ϕ0这一变量的相关性表明颗粒体系的摩擦特性受其密实状态影响,土颗粒越密实内摩擦特征越显著。另一变量颗粒惯性数I表征了颗粒形变剧烈程度客观规律。将发现的规律看成是考虑密实程度的准静态颗粒惯性数,定义它为Q=ϕ0[ln(I)+α]。

    临界土力学理论认为,在剪切的最终阶段,只有剪应变还持续增加,而土颗粒所受的有效应力和体积趋于不变。流态颗粒μ(I)和ϕ(I)流变准则认为摩擦系数μ分别与惯性数Iϕ一一对应,对应的是达到临界状态后受力的变形特征。经典的颗粒μ(I),ϕ(I)流变准则在描述I>10-3的密实均匀颗粒流时效果显著,但对颗粒准静态I<10-3动力特性描述时失效。发现的新标度律,揭示了颗粒准静态受力变形的本构法则。颗粒达到临界状态后产生流态化变形,ϕ(I)流变准则表明在固定的I条件下颗粒体积分数ϕ是定值,故颗粒的摩擦特性μ只反映出与颗粒惯性数I的对应关系。而在达到流态前的准静态变形状态,颗粒的摩擦特性还受体积分数ϕ0的影响:ϕ0越大,准静态惯性数Q=ϕ0 [ln(I)+α]越大,如图6示。也就意味着,颗粒在达到临界状态前发生准静态变形,对于同一种土,在同样的受力变形条件下(颗粒惯性数I相同,即外部围压(正应力)P和变形速率|˙γ|相同),更密实的土颗粒(更大的颗粒体积分数ϕ0)表现出更强的摩擦强度特性,这与常理认知一致。

    对线性关系进行标定,可得下列关系表达:

    μ=ξQ (14)

    其中,拟合参数见表2所示。可以发现,ξ近似等于0.11。

    图  6  更密实的土颗粒表现出更强的摩擦强度示意Q1<Q2
    Figure  6.  Denser soils exhibiting greater frictional strength, Q1<Q2

    从另一个角度来分析上述标度率所表征的体积分数对摩擦特性的影响,临界土力学定义了临界状态比Mc

    q=Mcp (15)

    这个固定的临界状态比与摩擦系数之间的关系为

    Mc=6μ/(3μ) (16)

    在固定的临界状态比(摩擦系数)条件下,正应力ps(=P)与颗粒体积分数的关系为[27-28]

    ps=aϕϕrlpϕrcpϕ (17)

    式中:ϕrlp为疏松颗粒临界状态下的体积分数;ϕrcp为密实颗粒临界状态下的体积分数;a为参数。

    根据莫尔库仑准则,无黏性颗粒强度τ可以表示为

    τ=μps (18)

    以莫尔库仑准则为纽带,提炼标度律所反映的体积分数与摩擦系数μ之间正相关的关系(P不变),与临界土力学给出的体积分数与正应力P间的关系是一致的,都反应了体积分数对颗粒强度正向影响。总的来讲,文章所分析的试验规律和新揭示的物理规律,把颗粒不同方面的性质,包括颗粒的强度μ、粒径d、颗粒密度ρs、密实状态ϕ0,正应力P和对应的变形速率|˙γ|有机的结合在一起,反应了颗粒重要的受力变形物理规律。

    表  2  3种材料三轴试验数据获取标度律的拟合参数
    Table  2.  Fitting parameters of scaling law for obtaining triaxial test data for three types of materials
    材料 ξ
    F-35 Ottawa砂 0.1021
    #1 玻璃砂 0.1173
    GS#40 Columbia grout砂 0.1179
    下载: 导出CSV 
    | 显示表格

    自然状态下的土体多处于三向应力状态,由于上述分析中p=(σ1+2σ3)/3只考虑了σ1σ3 的影响,需进一步考虑主应力σ2的影响。Xiao等[29]采用120 mm×60 mm×120 mm试样开展真三轴固结排水试验,研究了在较低应力水平下中主应力对堆石颗粒材料强度影响。引入毕肖普常数b描述中主应力大小:

    b=σ2σ3σ1σ3 (19)

    选取不同的毕肖普常数0,0.25,0.5,0.75,1开展独立试验。每次试验中,剪切过程中保持小主应力σ3b值不变,加载速率为恒定的1 mm/min。数据处理结果如图7所示。可以发现,在相同的试验条件下,摩擦系数μ随着b值的增大而减小。

    图  7  不同毕肖普常数bμ与ln(I)相关性影响
    Figure  7.  Influences of intermediate principal stress on relationship between μ and ln(I)

    土的三维莫尔库仑强度准则的表达形式为

    psinφ13q[13sinθsinφ+cosθ]+ccosφ=0 (20)

    其中应力罗德角θb之间符合下列关系

    tanθ=2b13 (21)

    认为土达到临界状态与弹塑性材料达到屈服准则进入塑性区的力学特征一致,也就是认为临界内摩擦角满足屈服准则。对于砂土c=0,将μ=tanφcs代入上式,可得

    pμq(13sinθμ+cosθ3)=0 (22)

    土达到临界状态时,根据剑桥模型得

    Mc=μ13sinθμ+cosθ3 = 1sinθ3+cosθ3μ (23)

    这里可以定义一个考虑应力罗德角(即中主应力影响)的摩擦系数μβ。考虑到Mc=6μβ/(3μβ),其中μβ=tanϕβ(cs)代入公式可得

    μβ=12(sinθ3+cosθ3μ + 16) (24)

    该式子定义了应力罗德角对临界状态内摩擦强度的影响。

    为定量描述中主应力(应力罗德角θ或毕肖普常数b)的影响,定义一个新的无量纲数“中主应力数”β=μ/μβ

    β=2μ(sinθ3+cosθ3μ + 16) (25)

    引入新的无量纲数“中主应力数”用于试验数据分析,发现μ+η(β-1)与准静态惯性数Q=ϕ0[ln(I)+α]之间呈线性相关关系,如图8所示。常规三轴试验中b=0,则中主应力数β=1,此时标度律回归到不考虑中主应力的标度律(即式(14))。

    图  8  通过中主应力数β反映中主应力对标度律的影响
    Figure  8.  Influences of intermediate principal stress on scaling law through intermediate principal stress number β

    上节通过建立准静态颗粒惯性数Q,建立标度率描述临界状态砂土的摩擦特性,本节探讨临界状态体积分数与土摩擦系数μ、颗粒惯性数I、以及初始体积分数ϕ0间的临界状态标度律。选取姜景山等[30]、Xiao等[31-32]的固结排水常规三轴压缩试验(试验中使用土样的物理参数如表3所示),试验发现土达到临界状态后,初始孔隙比、剪切速率、颗粒粒径对临界状态体积分数有直接的影响。导入试验数据,可以初步得到ϕ-ln(I)相关关系,如图9所示。可以看到,初始体积分数越小(初始孔隙比越大),临界状态的体积分数也越小。

    表  3  主要试验参数
    Table  3.  Primary experimental parameters
    初始孔隙比 Gs d50/mm 剪切速率
    0.360 2.72 14.39 0.6 mm/min 姜景山 等[30]
    0.333 2.72 14.39 0.6 mm/min
    0.308 2.72 14.39 0.6 mm/min
    0.283 2.72 14.39 0.6 mm/min
    0.189 2.69 26.48 1 mm/min Xiao 等[31]
    0.224 2.69 26.48 1 mm/min
    0.285 2.69 26.48 1 mm/min
    0.317 2.69 26.48 1 mm/min
    0.189 2.69 23.48 1 mm/min Xiao 等[32]
    0.224 2.69 23.48 1 mm/min
    0.285 2.69 23.48 1 mm/min
    0.317 2.69 23.48 1 mm/min
    下载: 导出CSV 
    | 显示表格
    图  9  初始孔隙比对ϕ与ln(I)关系的影响
    Figure  9.  Effect of initial void ratio on the relationship between ϕ and .ln(I)

    引入本文建立的准静态颗粒惯性数Q,可以得到准静态变形颗粒在临界状态下,颗粒体积分数ϕ与初始ϕ0和颗粒惯性数I间的标度律,如图10示。发现临界状态的颗粒体积分数与准静态惯性数Q= ϕ0[ln(I)+α]呈线性关系,可知准静态颗粒惯性数Q对揭示临界状态砂土的体积特征同样关键。

    图  10  不同初始孔隙比条件下ϕQ间的标度律
    Figure  10.  Scaling between ϕ and Q for sands with different initial void ratios

    本研究从土的颗粒物理本源出发,采用3种无黏性土的试验数据获取标度律的数理表达,发现临界状态摩擦系数μ值与惯性数I、颗粒初始体积分数ϕ0相关:土体越密实,临界状态摩擦系数越大。为定量描述土力学试验所揭示的土体强度与加载速率的相关性,引入了初始颗粒体积分数ϕ0来描述颗粒的初始密实程度,提出了考虑初始密实程度的准静态颗粒惯性数Q=ϕ0[ln(I)+α],这样通过ϕ0可以反映密实状态对颗粒体系的摩擦特性的影响,同时发现摩擦系数μ与准静态颗粒惯性数Q之间满足μ=ξQ的简单线性关系。

    经典的三轴试验还发现,体积分数越小(初始孔隙比越大),临界状态的颗粒体积分数也越小。为揭示临界状态土体的体变规律,采用本研究提出的准静态颗粒惯性数Q,得到了准静态变形的临界状态砂土,颗粒体积分数ϕ与初始体积分数ϕ0和颗粒惯性数I间的相关性。

    需要说明的是,自然土形貌特征等性质复杂,本文在建模分析先是将土颗粒抽象概化为圆形颗粒,开展无量纲数的数理推导,在此建模分析中粒径、压力、密度被认为是最主要的控制因素。土颗粒的重排列时间实际还与颗粒的其他尺寸参数、形貌特征等因素相关。

    此外,本研究基于多组试验数据定义的准静态颗粒惯性数Q=ϕ0[ln(I)+α]中,拟合参数值为α=30,过程中没有特意调整拟合参数,参数α的归一可能源自3种试验颗粒材料的基本性质相似。参数α对土类型和性质的依赖性,仍待进一步补充研究。考虑到研究所采用的三轴试验经典颗粒惯性数I>e -30,因此准静态颗粒惯性数Q恒大于零,符合定义的物理意义。对于I<e -30的情况,需要采用对应条件的试验数据拟合得到新的拟合参数α值。

  • 图  1   单向拉伸试验示意图

    Figure  1.   Diagram of uniaxial tensile tests

    图  2   应力-应变关系曲线

    Figure  2.   Relationship between stress and strain

    图  3   试验用土工膜应力应变关系曲线

    Figure  3.   Stress-strain curves of GMs used in tests

    图  4   杨氏模量/屈服应变与膜厚的关系

    Figure  4.   Relationship between Young's modulus/yield strain and thickness

    图  5   液胀试验破坏试样

    Figure  5.   Failure specimens in liquid expansion tests

    图  6   压力强度-液胀体积关系曲线

    Figure  6.   Relationship between pressure intensity and liquid expansion volume

    图  7   液胀模型拟合分析

    Figure  7.   Fitting analysis of liquid expansion model

    图  8   数值分析模型

    Figure  8.   Numerical analysis model

    图  9   拉伸强度-延伸率关系曲线

    Figure  9.   Relationship between tensile strength and elongation

    图  10   损伤试样拉伸强度-延伸率关系曲线

    Figure  10.   Relationship between tensile strength and elongation of damaged specimens

    图  11   泄压后土工膜变形图

    Figure  11.   Deformations of GMs after pressure relief

    表  1   土工膜主要参数表

    Table  1   Typical parameters of GMs

    技术指标 土工膜名称
    HDPE PE TPO PVC-HX PVC-CAPI
    平均厚度/ mm 0.5 1.0 1.2 1.5 2.0 0.5 0.8 1.0 3.65 1.0 2.0 2.85 3.65
    单位面积质量/
    (g·cm-2)
    0.082 0.096 0.114 0.135 0.194 0.074 0.066 0.1 0.608 0.14 0.278 0.354 0.42
    断裂强度/MPa 4.31 14.97 13.28 13.79 7.74
    屈服强度/MPa 1.94 1.75 1.64 1.62 2.8 1.77 1.27 1.0 1.32 4.09 4.51 3.01 1.75
    断裂延伸率/
    %
    144 166 185 256 257
    屈服延伸率/
    %
    0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 2.0 33.0 34.0 35.0 38.0
    注:技术指标依据SL235—2012《土工合成材料测试规程》测定;“—”表示未能在拉伸过程中发生损坏。
    下载: 导出CSV

    表  2   土工膜泊松比

    Table  2   Poisson's ratios of GMs

    名称 HDPE PE PVC-HX PVC-CAPI TPO
    厚度h/mm 0.50 1.00 1.20 1.50 2.00 0.50 0.80 1.00 1.00 2.00 2.85 3.65 3.65
    μy 0.53 0.42 0.46 0.60 0.44 0.56 0.50 0.43 0.44 0.42 0.51 0.54 0.51
    μz 0.47 0.58 0.54 0.40 0.56 0.44 0.50 0.57 0.56 0.58 0.49 0.46 0.49
    下载: 导出CSV

    表  3   应力-应变特征点参数

    Table  3   Parameters of characteristic points in stress-strain curves

    材料 厚度/mm εA/% σA/MPa εB/% σB/MPa
    HDPE 0.5 0.1 0.83 110 39.39
    1.0 0.1 0.85 108 40.37
    1.2 0.1 0.96 112 55.65
    1.5 0.1 0.93 109 58.43
    2.0 0.1 1.19 111 60.15
    PE 0.5 0.1 1.78 95 26.99
    0.8 0.1 0.69 90 26.10
    1.0 0.1 0.59 87 27.04
    PVC-HX 1.0 26.0 5.53 82 36.20
    2.0 28.0 5.35 85 32.38
    2.85 33.0 5.25 105 52.15
    PVC-CAPI 3.65 23.0 1.99 112 36.06
    TPO 3.65 2.0 1.35 124 28.84
    下载: 导出CSV

    表  4   杨氏模量成果表

    Table  4   Values of Young's modulus

    材料 厚度h/mm 线性区间/% 杨氏模量E/MPa R2
    HDPE 0.50 [0, 0.10] 837.11 0.9998
    1.00 [0, 0.10] 853.57 0.9978
    1.20 [0, 0.10] 874.82 0.9990
    1.50 [0, 0.10] 883.60 0.9999
    2.00 [0, 0.10] 1509.41 0.9996
    PE 0.50 [0, 0.10] 1228.00 0.9970
    0.80 [0, 0.10] 689.07 0.9998
    1.00 [0, 0.10] 545.22 0.9996
    PVC-HX 1.00 [0, 26.0] 21.27 0.9960
    2.00 [0, 28.0] 21.11 0.9960
    2.85 [0, 33.0] 19.05 0.9950
    PVC-CAPI 3.65 [0, 23.0] 15.70 0.9950
    TPO 3.65 [0, 2.0] 49.57 0.9960
    下载: 导出CSV

    表  5   液胀强度特征值

    Table  5   Characteristic values of liquid expansion strength

    材料
    名称
    厚度
    h/ mm
    峰值强度
    Pm/MPa
    峰值强度对应体积
    VP/(10-5 m³)
    最大液胀体积
    Vm/(10-5 m³)
    平台期范围
    V/(10-5 m³)
    HDPE 0.50 1.100 2.55 2.98 2.60~2.83
    1.00 2.274 3.68 3.68 0.00
    1.20 2.482 3.75 3.75 0.00
    1.50 2.789 3.75 3.75 0.00
    2.00 5.432 5.45 5.45 0.00
    PE 0.50 0.829 2.27 3.12 2.55~2.97
    0.80 1.476 2.20 3.12 2.76~3.12
    1.00 1.596 2.98 2.98 0.00
    PVC-HX 1.00 1.236 4.68 4.68 0.00
    2.00 2.095 4.75 4.75 0.00
    2.85 2.862 4.25 4.25 0.00
    PVC-CAPI 3.65 2.837 4.25 4.25 0.00
    TPO 3.65 1.659 3.26 5.53 3.75~5.38
    下载: 导出CSV

    表  6   液胀模型拟合参数

    Table  6   Fitting parameters of liquid expansion model

    名称 h/
    mm
    Pm/
    MPa
    XC/
    (10-5 m³)
    W/
    (10-5 m³)
    平均误差/百分比/%
    HDPE 0.50 1.100 2.55 0.63 0.0813/7.39
    1.00 2.274 3.68 0.82 0.0044/0.19
    1.20 2.482 3.75 0.94 0.0000/0.00
    1.50 2.789 3.75 1.08 0.0868/3.11
    2.00 5.432 5.45 1.75 0.1406/2.59
    PE 0.50 0.829 2.27 0.54 0.0429/5.18
    0.80 1.476 2.48 0.61 0.0025/0.17
    1.00 1.896 2.98 0.74 0.0015/0.08
    PVC-HX 1.00 1.236 4.68 1.22 0.0031/0.25
    2.00 2.095 4.75 1.25 0.0015/0.07
    2.85 2.862 4.25 1.40 0.0197/0.69
    PVC-CAPI 3.65 2.837 4.25 1.24 0.0596/2.10
    TPO 3.65 1.659 3.26 0.76 0.0069/0.42
    下载: 导出CSV

    表  7   液胀强度

    Table  7   Liquid expansion strengths

    材料
    名称
    厚度
    h/ mm
    计算液胀
    强度P/MPa
    极限液胀
    强度Pm/MPa
    允许液胀
    强度Pa/MPa
    HDPE 0.50 0.42 1.100 0.44
    1.00 0.15 2.274 0.91
    1.20 0.20 2.482 0.99
    1.50 0.28 2.789 1.12
    2.00 0.39 5.432 2.17
    PE 0.50 0.36 0.829 0.33
    0.80 0.16 1.476 0.59
    1.00 0.30 1.596 0.64
    PVC-HX 1.00 0.22 1.236 0.49
    2.00 0.28 2.095 0.84
    2.85 0.45 2.862 1.14
    PVC-CAPI 3.65 0.32 2.837 1.13
    TPO 3.65 0.29 1.659 0.66
    下载: 导出CSV

    表  8   损伤试样拉伸强度/延伸率表

    Table  8   Tensile strengths/elongations of damaged specimens

    材料 厚度/mm 横向 斜向
    极限延伸率/% 极限拉
    伸强度/(N·mm-1)
    极限
    延伸率/%
    极限拉
    伸强度/(N·mm-1)
    HDPE 0.50 8.80 5.13 25.20 5.20
    1.00 77.27 9.55 104.53 10.53
    1.20 78.20 10.05 110.81 10.58
    PE 0.50 10.80 3.68 14.85 3.95
    0.80 83.94 6.50 203.88 6.89
    1.00 96.73 8.62 254.67 22.62
    PVC-HX 1.00 156.23 10.41 225.14 13.82
    2.00 189.30 16.66 269.16 19.47
    2.85 102.80 12.80 209.40 19.90
    PVC-CAPI 3.65 110.20 9.40 249.30 21.00
    下载: 导出CSV

    表  9   土工膜变形量统计表

    Table  9   Deformations of GMs

    名称 厚度/ mm 初始变形量/mm 12 h后变形量/ mm 回复量/ mm 回复率/%
    HDPE 0.50 8.20 5.60 2.60 31.71
    1.00 7.40 3.90 3.50 47.30
    1.20 7.10 2.70 4.40 61.97
    1.50 6.80 2.30 4.50 66.18
    2.00 6.40 1.70 4.70 73.44
    PE 0.50 8.80 7.70 1.10 12.50
    0.80 7.80 5.90 1.90 24.36
    1.00 7.70 4.70 3.00 38.96
    PVC-HX 1.00 7.30 1.90 5.40 73.97
    2.00 5.10 0.20 4.90 96.08
    2.85 4.35 0.15 4.20 96.55
    PVC-CAPI 3.65 5.65 0.05 5.60 99.12
    TPO 3.65 3.45 0.65 2.80 81.16
    下载: 导出CSV
  • [1] 束一鸣. 我国水库大坝土工膜防渗工程进展[J]. 水利水电科技进展, 2015, 35(5): 20-26. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201505005.htm

    SHU Yiming. Progress in geomembrane barriers for seepage prevention in reservoirs and dams in China[J]. Advances in Science and Technology of Water Resources, 2015, 35(5): 20-26. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201505005.htm

    [2] 顾淦臣. 土工薄膜在坝工建设中的应用[J]. 水力发电, 1985, 11(10): 43-50. https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD198510012.htm

    GU Ganchen. Application of geomembrane in dam construction[J]. Water Power, 1985, 11(10): 43-50. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD198510012.htm

    [3] 束一鸣, 吴海民, 姜晓桢. 高面膜堆石坝发展的需求与关键技术: 高面膜堆石坝关键技术(一)[J]. 水利水电科技进展, 2015, 35(1): 1-9.

    SHU Yiming, WU Haimin, JIANG Xiaozhen. Need and key technology on development of high membrane faced rockfill dam: key technology of high membrane faced rockfill dam(Ⅰ)[J]. Advances in Science and Technology of Water Resources, 2015, 35(1): 1-9. (in Chinese)

    [4] 束一鸣, 吴海民, 姜晓桢, 等. 高面膜堆石坝周边的夹具效应机制与消除设计方法: 高面膜堆石坝关键技术(二)[J]. 水利水电科技进展, 2015, 35(1): 10-15. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201501003.htm

    SHU Yiming, WU Haimin, JIANG Xiaozhen, et al. Mechanism of anchoring influence at perimeter of high membrane faced rockfill dam and eliminating approach: key technology of high membrane faced rockfill dam(Ⅱ)[J]. Advances in Science and Technology of Water Resources, 2015, 35(1): 10-15. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201501003.htm

    [5] 吴海民, 束一鸣, 姜晓桢, 等. 高面膜堆石坝运行状态下土工膜双向拉伸力学特性: 高面膜堆石坝关键技术(三)[J]. 水利水电科技进展, 2015, 35(1): 16-22.

    WU Haimin, SHU Yiming, JIANG Xiaozhen, et al. Biaxial tensile mechanical property of geomembrane used as high membrane faced rockfill dam: key technology of high membrane faced rockfill dam(Ⅲ)[J]. Advances in Science and Technology of Water Resources, 2015, 35(1): 16-22. (in Chinese)

    [6] 吴海民, 束一鸣, 滕兆明, 等. 高堆石坝面防渗土工膜锚固区夹具效应破坏模型试验[J]. 岩土工程学报, 2016, 38(S1): 30-36. doi: 10.11779/CJGE2016S1006

    WU Haimin, SHU Yiming, TENG Zhaoming, et al. Model tests on failure properties of geomembrane anchorage due to clamping effect in surface barrier of high rock-fill dam[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(S1): 30-36. (in Chinese) doi: 10.11779/CJGE2016S1006

    [7] 张宪雷, 刘云锋, 顾克, 等. 高面膜土石坝防渗结构中土工膜弯折(褶皱)试验研究[J]. 岩土工程学报, 2019, 41(8): 1555-1561. doi: 10.11779/CJGE201908021

    ZHANG Xianlei, LIU Yunfeng, GU Ke, et al. Experimental study on geomembrane bending (folding) in anti-seepage structure of membrane-faced rockfill dam[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(8): 1555-1561. (in Chinese) doi: 10.11779/CJGE201908021

    [8] 宁宇, 喻建清, 崔留杰. 软岩堆石高坝土工膜防渗技术[J]. 水力发电, 2016, 42(5): 62-67, 105. https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD201605020.htm

    NING Yu, YU Jianqing, CUI Liujie. Anti-seepage of geomembrane for high soft rock filling dam[J]. Water Power, 2016, 42(5): 62-67, 105. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD201605020.htm

    [9]

    ICOLD. Geomembrane Sealing Systems for Dams: Design Principles and Return of Experience (Bulletin 135)[R]. Paris: International Commission on Large Dams, 2010.

    [10]

    SCUERO A, VASCHETTI G. Geomembrane sealing systems for dams: ICOLD bulletin 135[J]. Innovative Infrastructure Solutions, 2017, 2(1): 29.

    [11]

    CAZZUFFI D, MORACI N, CALVARANO L S, et al. The influence of vertical effective stress and geogrid length on interface behavior under pullout conditions[J]. Geosynthetics, 2014, 32(2): 40-50.

    [12] 束一鸣, 吴海民, 姜晓桢. 中国水库大坝土工膜防渗技术进展[J]. 岩土工程学报, 2016, 38(增刊1): 1-9.

    SHU Yiming, WU Haimin, JIANG Xiaozhen. The development of anti-seepage technology with geomembrane on reservoirs and dams in China[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(S1): 1-9. (in Chinese)

    [13] 中华人民共和国水利部. 碾压式土石坝设计规范: SL 274—2020[S]. 北京: 中国水利水电出版社, 2020.

    Ministry of Water Resources of the People's Republic of China. Design code for rolled earth-rock fill dams: SL 274—2020[S]. Beijing: China Water & Power Press, 2020. (in Chinese)

    [14] 国家能源局. 碾压式土石坝施工组织设计规范: NB/T 35062—2015[S]. 北京: 中国水利水电出版社, 2016.

    National Energy Bureau of the People's Republic of China. Design Code for Construction Planning of Rolled Earth-Rock Dam: NB/T 35062—2015[S]. Beijing: China Water & Power Press, 2016. (in Chinese)

    [15]

    CHEVALIER L, CALLOCH S, HILD F. Digital image correlation used to analyze the multiaxial behavior of rubber-like materials[J]. European Journal of Mechanics, 2001, 20(2): 169-187.

    [16]

    AYDILEK A H, GULER M, EDIL T B. Use of image analysis in determination of strain distribution during geosynthetic tensile testing[J]. Journal of Computing in Civil Engineering, 2004, 18(1): 65-74.

    [17]

    GÓRSZYK J, MALICKI K. Digital image correlation method in monitoring deformation during geogrid testing[J]. Fibres and Textiles in Eastern Europe, 2019, 27(2(134)): 84-90.

    [18]

    MERRY S M, BRAY J D. Geomembrane response in the wide strip tension test[J]. Geosynthetics International, 1996, 3(4): 517-536.

    [19]

    Standard Practice for Determination of 2% Secant Modulus for Polyethylene Geomembranes: ASTM D5323-92[S]. 2018.

    [20]

    ZHANG X L, MA Z Y, WU Y Y, et al. Response of mechanical properties of polyvinyl chloride geomembrane to ambient temperature in axial tension[J]. Applied Sciences, 2021, 11(22): 10864.

    [21] 束一鸣, 叶乃虎. LDPE土工膜液胀极限荷载的工程仿真实验[J]. 水利水电科技进展, 2003, 23(5): 1-3, 69. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD200305001.htm

    SHU Yiming, YE Naihu. Simulation test for Mullen burst limit load of LDPE geomembrane[J]. Advances in Science and Technology of Water Resources, 2003, 23(5): 1-3, 69. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD200305001.htm

    [22] 李旺林, 刘占磊, 孟祥涛, 等. 土工膜环向约束气胀变形试验研究[J]. 岩土工程学报, 2016, 38(6): 1147-1151. doi: 10.11779/CJGE201606023

    LI Wanglin, LIU Zhanlei, MENG Xiangtao, et al. Experimental study on air expansion deformation of geomembrane under ring-restrained conditions[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(6): 1147-1151. (in Chinese) doi: 10.11779/CJGE201606023

    [23]

    XU F, LI W L, LIU Z L, et al. Study of factors that influence geomembrane air expansion deformation under ring-restrained conditions[J]. Geotextiles and Geomembranes, 2017, 45(3): 178-183.

    [24] 束一鸣, 陆忠民, 侯晋芳, 等. 土工合成材料防渗排水防护设计施工指南[M]. 北京: 中国水利水电出版社, 2020.

    SHU Yiming, LU Zhongmin, HOU Jinfang. Guide for Design and Construction of Seepage Control and Drainage Protection of Geosynthetics[M]. Beijing: China Water & Power Press, 2020. (in Chinese)

    [25] 徐光明, 章为民, 彭功勋. HDPE膜的力学特性受损伤影响初步研究[J]. 河海大学学报(自然科学版), 2004, 32(1): 76-80.

    XU Guangming, ZHANG Weimin, PENG Gongxun. Effect of damage on mechanical behavior of HDPE geomembrane[J]. Journal of Hehai University (Natural Sciences), 2004, 32(1): 76-80. (in Chinese)

    [26] 岑威钧, 波兰汗·开肯, 刘凤茹. 缺陷和损伤土工膜的拉伸特性试验[J]. 水利水电科技进展, 2021, 41(1): 69-73. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD202101012.htm

    CEN Weijun, BOLANHAN Kaiken, LIU Fengru. Experimental study on tensile properties of defective and damaged geomembranes[J]. Advances in Science and Technology of Water Resources, 2021, 41(1): 69-73. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD202101012.htm

    [27] 吴海民, 束一鸣, 滕兆明, 等. 高堆石坝面膜防渗体非散粒体垫层工程特性试验: 高面膜堆石坝关键技术(五)[J]. 水利水电科技进展, 2015, 35(1): 29-36, 66. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201501007.htm

    WU Haimin, SHU Yiming, TENG Zhaoming, et al. Experimental research on engineering property of non-granular cushion for high membrane faced rockfill dam: key technology of high membrane faced rockfill dam(Ⅴ)[J]. Advances in Science and Technology of Water Resources, 2015, 35(1): 29-36, 66. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201501007.htm

    [28] 姜晓桢, 束一鸣. 基于最大熵原理的土工膜与垫层颗粒间接触力随机状态分析[J]. 岩土工程学报, 2016, 38(增刊1): 49-55. doi: 10.11779/CJGE2016S1009

    JIANG Xiaozhen, SHU Yiming. Stochastic analysis of contact force between geomembrane and particle underlayer based on principle of maximum entropy[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(S1): 49-55. (in Chinese) doi: 10.11779/CJGE2016S1009

图(11)  /  表(9)
计量
  • 文章访问数:  273
  • HTML全文浏览量:  66
  • PDF下载量:  60
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-02-10
  • 网络出版日期:  2023-05-18
  • 刊出日期:  2023-04-30

目录

/

返回文章
返回