Processing math: 100%
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

对《“矢量和法”抗滑稳定计算模型的力学概念错误》一文的回应

郭明伟, 李春光, 刘艳章, 吴振君, 邓琴

郭明伟, 李春光, 刘艳章, 吴振君, 邓琴. 对《“矢量和法”抗滑稳定计算模型的力学概念错误》一文的回应[J]. 岩土工程学报, 2022, 44(9): 1760-1762. DOI: 10.11779/CJGE202209022
引用本文: 郭明伟, 李春光, 刘艳章, 吴振君, 邓琴. 对《“矢量和法”抗滑稳定计算模型的力学概念错误》一文的回应[J]. 岩土工程学报, 2022, 44(9): 1760-1762. DOI: 10.11779/CJGE202209022
GUO Ming-wei, LI Chun-guang, LIU Yan-zhang, WU Zhen-jun, DENG Qin. Reply to "Mechanical concept errors in anti-sliding stability computational model of vector sum method"[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(9): 1760-1762. DOI: 10.11779/CJGE202209022
Citation: GUO Ming-wei, LI Chun-guang, LIU Yan-zhang, WU Zhen-jun, DENG Qin. Reply to "Mechanical concept errors in anti-sliding stability computational model of vector sum method"[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(9): 1760-1762. DOI: 10.11779/CJGE202209022

对《“矢量和法”抗滑稳定计算模型的力学概念错误》一文的回应  English Version

详细信息
    作者简介:

    郭明伟(1981—),男,河南商丘人,副研究员,主要从事边坡工程、深基础工程及计算岩土力学方面的研究。E-mail:mwguo@whrsm.ac.cn

  • 中图分类号: TU431

Reply to "Mechanical concept errors in anti-sliding stability computational model of vector sum method"

  • 成都市主城八区均下伏中等风化泥岩,埋深一般为-0.5~-35m,层厚十几米到几十米不等[1-2],工程性质差,承载力估算困难。随着成渝地区双城经济圈建设的深度推进,成都地区已开建或筹建的高层、超高层建筑的基础持力层均涉及中等风化泥岩,研究该类岩体的承载能力具有一定的工程意义。

    江苏[3]、广东[4]、重庆[5]等地在实践中逐渐修正了规范[6]对软岩地基承载力的取值建议,提出了不同岩石单轴抗压强度分区所对应的折减系数及承载力;长沙、兰州等地以岩基载荷试验为主,辅以点荷载试验、旁压试验等,提出了泥岩承载力的界定标准[7-9];并通过不同类型试验结果的对比,探讨了影响地基承载力取值的因素[10-11],指出确定地基承载力时应考虑岩体的物理力学性质、成因、年代、地下水及其变化情况、建筑物对沉降及差异沉降的要求等[12]。成都地标DB51/T5026—2001[13]建议中等风化泥岩地基承载力特征值为500~1000 kPa,也可根据室内天然单轴抗压强度折减法确定地基承载力,工程实践表明规范条文存在两点不足。一方面是建议的承载力取值偏低,成都地区中风化泥岩极限承载力值试验值普遍分布在2000~6000 kPa,且主要集中在4000~6000 kPa[14],按照安全系数3考虑,其承载力特征值亦大于1000 kPa;另一方面,折减法确定承载力合宜性差,折减系数ψ的取值无既定条款可循,往往带有过多的人为因素,即使对同一类型、同一条件、同一环境的地基岩体,不同的设计人员也会取不同的ψ值,有时甚至出入较大,在没有原位试验的情况下,所提出的承载力缺乏依据,带有一定的随意性,对工程的安全性没有切实保障。在成都地区更为合理的取值研究至今鲜有突破,受到不少诟病。

    综上,本文通过对24个场地约100余处原位平板载荷试验与同条件岩石单轴抗压强度试验工程案例进行对比分析,讨论了成都地区中等风化泥岩计算公式中ψ值取值及其与承载力特征值的相关关系。进一步,结合成都高新区在建的某超高层建筑项目,采用原位平板载荷试验(9组)联合钻孔旁压试验法(17孔)针对性的探讨了该类岩体承载力特性,并分析了承压板尺寸等因素对承载能力发挥程度的影响。以上研究旨在探寻单轴抗压强度或旁压试验结果作为承载力设计依据的可能。

    成都地区红层属于盆西小区[15],是最具代表性的近水平红层分布区。区内中等风化泥岩以J3p侏罗系蓬莱镇组上统和K2g白垩系中统灌口组为主(见图1),其中J3p泥岩大范围的分布在成都天府新区高楼山—秦皇寺中央商务区—倒石桥一带,其他区域均为K2g泥岩,从城区东南方至西北方埋深为10~30 m,深度逐渐加深。

    图  1  成都地区中等风化泥岩分布略图
    Figure  1.  Distribution of moderately weathered mudstone in Chengdu

    区内中等风化泥岩呈紫红色,层理清晰,风化裂隙较发育,巨厚层构造,整体结构。锤击易碎、声哑,用镐难挖掘,岩芯钻方可钻进。岩芯采取率达90%以上(图2)。通过对岩芯的X衍射和矿物成份分析,可知区内泥岩及其风化物含有超50%的SiO2、Al2O3、Fe2O3,其不同的组合构成了以蒙脱石、伊利石为主的片状黏土矿物。表1为岩体物理力学参数统计表。

    图  2  成都地区典型中等风化泥岩岩芯
    Figure  2.  Typical moderately weathered mudstone core in Chengdu
    表  1  中等风化泥岩物理力学性质统计表
    Table  1.  Physical and mechanical properties of rock
    统计指标天然密度/(g·cm-3)单轴抗压强度/MPa黏聚力/MPa内摩擦角/(°)
    天然状态饱和状态天然状态饱和状态天然状态饱和状态
    最大值2.629.678.201.20.8540.138.6
    最小值2.302.711.560.40.2333.731.4
    平均值2.486.003.880.80.4937.935.5
    样本容量141513928172817
    下载: 导出CSV 
    | 显示表格

    目前成都地区确定岩石地基承载力特征值主要为岩石单轴抗压强度折减法和原位平板岩基载荷试验法:①第一种方法简便但计算结果与实际地基承载力有较大出入;②第二种计算方法得到的结果普遍被认为是真实的,但费时费力。

    笔者整理了24个工程实例中约100个中等风化泥岩的原位平板载荷试验点的相关测试数据,开展了2种方法确定地基承载力的对比分析,结果见表2,试验点覆盖了成都主城八区,具有较好的代表性。需要说明的是,文中所示意的原位平板载荷试验点绝大多数未达到明显的极限荷载,因此岩基的承载力特征值确定标准:①极限荷载fu/3与比例界限荷载的较小值;②未达到极限状态取比例界限荷载。

    表  2  岩基载荷试验与单轴抗压强度试验对比
    Table  2.  Comparison between load tests and uniaxial compressive strength tests
    工程名称工程地点frk/MPafa/kPaψ停载位移/mm点数
    隆鑫成华区4.0020440.515.4203
    润都科技成华区4.0019810.496.7603
    蓉华·上林成华区3.0518410.603.9908
    17730.589.730
    18610.618.680
    科创中心成华区5.0022630.4518.3708
    东方荟锦江区5.0021570.433.6404
    东大街项目锦江区3.0719000.623.5203
    仁恒置地广场锦江区3.4319900.582
    航天科技大厦锦江区6.0028800.483
    西部国金锦江区1.488430.5710.4103
    融御大楼锦江区7.3030100.410.5623
    海桐三期锦江区7.9029590.378.3004
    滨江30号地块高新区6.8027400.4010.7003
    半岛城邦高新区4.5020210.456.6009
    龙湖世纪城高新区2.7918000.6533.1908
    川大科技园双流区3.9218940.485.2903
    水电·云立方龙泉驿5.9824000.405.5003
    视高天府公园天府新区6.4824000.377.3203
    学府城天府新区4.9922900.461.4204
    22660.452.7202
    仁中风华一期天府新区5.6424000.435.9806
    驿都水岸丽园天府新区5.5825240.4517.5001
    5.0021620.4318.6302
    东方希望天祥武侯区4.3920000.501.1603
    海洋乐园武侯区6.5525000.384.2606
    中邮金融中心武侯区5.7424000.422.2003
    田园晶品新都区3.116970.553.7206
    注:frk为岩石天然单轴抗压强度标准值;fa为承载力特征值;ψ为折减系数。
    下载: 导出CSV 
    | 显示表格

    表2中,成都主城八区中等风化泥岩地基承载力特征值试验值70%以上集中在1500~2500 kPa,较规范DB51T5026建议最大值普遍大1.5倍以上,统计原位平板载荷试验与同条件岩石单轴抗压强度试验结果得出折减系数分布在0.37~0.70,总体值大于0.5(表3);综合可见,折减系数呈现出随单轴抗压强度增加而递减的趋势。

    表  3  折减系数和承载力特征值建议值
    Table  3.  Recommended values of reduction coefficient and characteristic value of bearing capacity
    frk/MPaψfa/kPa说明
    <20.7500~1400符号意义同表2,岩石单轴抗压强度为天然状态
    2~40.5~0.71400~2000
    4~60.4~0.52000~2400
    6~80.37~0.42400~3000
    >80.37>3000
    下载: 导出CSV 
    | 显示表格

    表3所列的统计结果对于岩石单轴抗压强度小于2 MPa和大于8 MPa的岩体,由于没有试验样本,可分别取统计数据中折减系数最高区段和最低区段的ψ值,即取0.7,0.37,承载力特征值最小值按现行规范[12]最小值500 kPa考虑。同时,上述统计结果多数是在未加载到岩基破坏的试验条件下得到的,实际的折减系数比该统计值还要大。但是,仍要注意的是,表3中岩石单轴抗压强度为天然状态,若当基岩浸水软化时,需要适当考虑折减系数减小。

    进一步依托成都高新区在建的某超高层建筑项目,探讨不同承压板直径对研究区中等风化泥岩地基承载力取值的影响。

    该超高层建筑位于成都市南部天府新区秦皇寺中央商务区。场地地貌单元属宽缓浅丘,为剥蚀型浅丘陵地貌,建筑基础底板持力层为侏罗系蓬莱镇组(J3p)中等风化泥岩。泥岩天然密度2.48 g/cm3、自然状态下单轴抗压强度平均值6.00 MPa、饱和状态单轴抗压强度平均值3.88 MPa。三轴天然抗剪强度参数平均值内摩擦角为37.9°,黏聚力为0.8 MPa。同时,根据钻孔波速测得岩体波速为2200~3200 m/s;室内岩块波速平均值为2370 m/s。

    承压板载荷试验在场地开挖的3口深井平硐中进行(编号分别为SJ01,SJ02,SJ03)。试验点点位见图3所示。

    图  3  试验点位示意图
    Figure  3.  Schematic diagram of test sites

    图3中,SJ01深30 m(高程459.06 m)、SJ02深36 m(高程451.15 m)、SJ03深39 m(高程447.57 m),井底平硐长×宽×高为2 m×2 m×8 m。每一平硐内分别开展承压板直径R为300,500,800 mm的载荷试验各1组,共9组。载荷试验操作过程严格参考规范[16],开挖至目标层后对试验点进行粗砂找平,放置承压板、千斤顶、百分表等试验设备。检查无误后采用慢速维持加载法加载,加载分8~12级,每级加载后连续三次读数之差均不大于0.01 mm可施加下一级荷载;当某级荷载下24 h变形速率不能达到相对稳定标准或本级荷载的变形量大于前级荷载变形量的5倍或承压板周边的岩土出现明显隆起或径向裂缝持续发展等现象终止加载。试验现场照片见图4

    图  4  岩基平板载荷试验过程照片
    Figure  4.  Photos of bearing plate tests

    试验数据绘制ps曲线,曲线起始直线段终点的前一级荷载为比例界限荷载;终载前一级荷载为极限荷载。将极限荷载除以3的安全系数,所得值与对应于比例界限的荷载相比较,取小值定为承载力特征值。

    典型ps曲线见图5所示。

    图  5  ps曲线(以02深井300 mm承压板为例)
    Figure  5.  ps curves (deep well No. 2, plate diameter 300 mm)

    (1)曲线整体上表现出变形随着荷载的增大而逐渐增大。

    (2)以SJ02中直径300 mm载荷试验结果为例,试验荷载分12级加载,试验结束荷载为7800 kPa。加载结束时,位移持续增加且不能稳定,累计变形量为4.182 mm。从曲线上均可直接读出比例界限荷载和极限荷载。曲线第一个拐点(2400 kPa)对应的前一级荷载(2100 kPa)为比例界限荷载;极限荷载为曲线第二个拐点(7800 kPa)对应的前一级荷载(6300 kPa)。极限荷载的1/3与比例界限荷载对比取小值为承载力特征值,即为2100 kPa。

    承载板直径与承载力特征值之间的关系见图6所示。其他试验点试验结果统计见表4所示。

    图  6  承载板直径与承载力特征值的关系
    Figure  6.  Plate diameter vs. bearing capacity
    表  4  研究区中等风化泥岩承载力特征值试验值统计表(承压板载荷试验)
    Table  4.  In-suit test values of bearing capacity of moderately weathered mudstone in study area
    试验点号
    竖井编号
    承压板直径D/mm原位试验极限值/3对应荷载fa/kPa岩体完整性指数Kv
    最大加载量/kPa比例界限对应荷载/kPa终止条件裂缝情况
    SJ01
    (459.06 m)
    30090003000趋近试验装置系统极限24000.46
    50096002400位移持续增大,不能稳定24000.46
    80063002700沉降增大,大于前一级沉降量的5倍23000.43
    SJ02
    (451.15 m)
    30078002400变形不能保持稳定21000.39
    500108002400变形不能保持稳定20000.39
    800117002700沉降增大,大于前一级沉降量的2倍21000.39
    SJ03
    (447.54 m)
    300118003600趋近试验装置系统极限28000.49
    50096003000趋近试验装置系统极限30000.50
    800108003600支硐上壁无法提供足够反力27000.43
    下载: 导出CSV 
    | 显示表格

    对于同一平硐,中等风化泥岩地基承载力特征值差异小于10%。具体来说,SJ01深井泥岩承载力特征值为2300~2400 kPa;SJ02深井泥岩承载力特征值为2000~2100 kPa;SJ03深井泥岩承载力特征值为2700~3000 kPa。说明,相似性状的泥岩,载荷板尺寸对于其地基承载力性能影响程度不大(见图7所示),可推测基础宽度对于同一性状的中等风化泥岩地基承载力的影响不明显。

    图  7  旁压试验过程照片
    Figure  7.  Photos of pressuremeter tests

    不同平硐内试验结果略有差异,与岩体个体差异尤其是岩体完整性和矿物组分不无关系,试验表明:三口平硐内泥岩岩体完整性指数分别为0.43~0.46,0.39,0.43~0.50;黏土矿物含量分别为29%~30%,36%~47%,17%~25%。上述差异导致不同平硐内中等风化泥岩地基承载性能出现差异。另外,岩石地基承载性能差异亦是试验环境条件、操作人为误差等因素综合影响所造成,需进一步详细甄别和探讨。

    上述超高层建筑项目中,3口深井平硐中进行承压板载荷试验同时亦在300 mm压板试验点位附近各取3个试样进行了室内单轴抗压强度试验(试验参照规范[16]),对修正折减系数法进行验证验证。根据(DB51/T5026—2001)[13]和(GB50007—2011)[6]中确定承载力特征值fa公式为

    fa=ψfrk, (1)

    式中,fa为岩石地基承载力特征值(kPa),frk为岩石饱和单轴抗压强度标准值(kPa),软岩取为天然单轴抗压强度标准值),ψ为折减系数,由当地经验确定。

    结合前述研究,对式(1)参数进行修正,得

    fa=ψfrk, (2)

    式中,fa为岩石地基承载力特征值(kPa),frk为岩石饱和单轴抗压强度标准值(kPa,软岩取为天然单轴抗压强度标准值),ψ为修正折减系数(按表3选取)。

    试验及验证结果见表5所示。

    表  5  修正折减系数法验证结果
    Table  5.  Test results of modified reduction coefficient method
    竖井编号编号fa/kPafrk/MPa平均frk/MPaψ计算值差异/%
    SJ01124004.95.960.4023840.6
    25.8
    37.19
    SJ02121003.4740.5020004.7
    23.78
    34.39
    SJ03128006.897.10.3927691.1
    28.31
    38.41
    下载: 导出CSV 
    | 显示表格

    表5中计算结果可见,本文所提出的折减系数修正法所估算的岩基承载力特征值与原位载荷试验结果相近,差异在5%以内,效果较好,可以作为承载力计算的依据。根据单轴抗压强度分区不同其修正系数可按表3选取,取值区间可采用线性内插来确定。

    旁压试验其实质是横向载荷试验。软岩中旁压临塑荷载既与岩石单轴抗压强度有关又与岩体的裂隙发育程度有关,利用旁压试验能更好地反映工程岩体的实际强度[10-11]。近年来,在长沙等地得到广泛应用,然而其在成都地区中等风化泥岩中的适宜性仍需进一步探讨。在上述项目场地中14口勘察钻孔+3口对比钻孔中旁压试验(见图5),测试深度15~30 m不等。其中,3口对比钻孔TL05,JK14,JK12紧邻SJ01,SJ02,SJ03使试验结果具有可对比性。

    采用PM-2B型预钻式旁压仪进行试验(见图7),主机尺寸23 cm×36 cm×90 cm,最大压力6 MPa,试验操作严格参考规范[17,18],根据试验结果绘制PV曲线,进一步确定初始压力P0和临塑压力Pf。两个压力差值为承载力特征值fak

    旁压试验结果见表6所示。

    表  6  旁压试验结果统计表
    Table  6.  Results of pressuremeter tests
    钻孔编号测试深度/mPf/kPaP0/kPafa/kPa钻孔波速/(m·s-1)岩体完整性指数Kv
    JK0428.53002.4625.22377.224650.380
    JK0321.01399.5350.01049.518730.231
    JK1020.01294.2310.5983.720230.256
    JK1213.01495.5311.21184.321600.292
    JK1220.02670.0698.51971.522040.304
    JK1428.03753.0850.62902.429120.507
    TL0126.53790.2849.62940.624450.414
    TL0318.54085.11120.32964.827540.530
    TL0528.02828.4475.62352.825420.447
    TL1120.53715.9850.62865.330130.469
    TL1316.53342.3942.32400.025550.452
    TL1323.02209.1487.21721.922410.348
    TL1417.52300.5422.11878.426130.473
    TL1614.52489.0475.62013.421780.296
    TL1620.52737.9475.62262.328240.498
    TL1823.02590.6432.22158.425230.441
    TL2223.02856.7511.22345.529300.415
    TL2524.52707.0530.12176.923470.381
    TL2526.04773.5681.24092.330940.663
    TL3523.03032.5625.22407.326590.490
    TL3718.02494.4550.01944.426900.501
    下载: 导出CSV 
    | 显示表格

    表6中可见,中等风化泥岩地基承载力特征值测试结果最小值为983.7 kPa,最大值为4092.3 kPa,范围较大,这与每一测试孔测试时周围环境、成孔质量、岩体质量、测试深度等有一定关系。因测试时环境和成孔质量属于随机无量化指标的影响因素。进一步对比表4,6中深井及其对比钻孔试验结果,以探究旁压试验适宜性。以深井平硐内的承压板试验结果为基准。对比结果见图8所示。

    图  8  试验结果对比图
    Figure  8.  Comparison between test results

    对比结果显示,各测点的泥岩均为中等风化泥岩,测试深度、岩体完整性均相似,两种方法测试结果极为相近。尤其是深井SJ01和钻孔TL05,测试深度分别为30 m和28 m,岩体完整性指数均为0.46左右,得到的地基承载力特征值2300 kPa范围,可以认为采用承压板载荷试验、钻孔旁压试验确定的泥岩地基承载力特征值fa可互为印证。

    本文通过对24个场地约100个原位平板载荷试验与同条件岩石单轴抗压强度试验工程案例进行对比分析,讨论了成都地区中等风化泥岩计算公式中折减系数取值。进一步,依托成都高新区在建的某超高层建筑项目,探讨了修正折减系数法和旁压试验结果作为设计依据的建议。

    (1)成都主城八区中等风化泥岩地基承载力特征值主要集中在1500~2500 kPa。

    (2)统计原位平板载荷试验与同条件岩石单轴抗压强度试验结果得出折减系数分布在0.37~0.70,总体值大于0.5。

    (3)本文所提出的折减系数修正法所估算的岩基承载力特征值在实际工程中实用性较强,其估算结果与原位载荷试验结果差异在5%以内,效果较好,可以作为承载力计算的依据。

    (4)对于相似性状的泥岩,泥岩地基承载力受承压板直径影响不显著;且承压板载荷试验、钻孔旁压试验确定的地基承载力特征值可互为印证;

    (5)上述成果,可用于该地区同类岩石地基参考,但是岩土体物理力学参数间的相关性取决于样本容量,对上述相关内容的细化有待更多的试验对比和检验及完善。

  • 图  1   矢量和法安全系数公式中定义的矢量

    Figure  1.   Sketch of the vector of the safety factor in the vector sum method

    图  2   该算例不同方向角

    Figure  2.   Different sliding angles of the slope example

    图  3   极限平衡M-P法计算模型

    Figure  3.   The calculating model of limit equilibrium M-P method

    图  4   两种方法安全系数对比图

    Figure  4.   Comparison of safety factors with both methods

    图  5   矢量和法与M-P法结果相对误差

    Figure  5.   Relative error of safety factors between the vector sum method and the limit equilibrium M-P method

    表  1   标准考题EX1C安全系数计算结果汇总表

    Table  1   The calculating results of the slope example EX1C

    计算方法 安全系数 下滑方向/(°)
    极限平衡法 推荐裁判答案Donald 1.3900
    SSA(Baker) 1.3900
    STAB(陈祖煜) 1.3850
    EMU 1.3900
    Fredlund 1.4060
    平均值 1.3922
    孙建生法 孙建生 1.4311 45.533
    本文方法 矢量和法 1.3939 23.1938
    下载: 导出CSV
  • [1] 孙建生. "矢量和法"抗滑稳定计算模型的力学概念错误[J]. 岩土工程学报, 2021, 43(5): 975–980. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202105028.htm

    SUN Jian-sheng. Mechanical concept errors in anti-sliding stability computational model of "vector sum method"[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(5): 975–980. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202105028.htm

    [2] 李广信, 张丙印, 于玉贞, 等. 土力学[M]. 2版. 北京: 清华大学出版社, 2018.

    LI Guang-xin, ZHANG Bing-yin, YU Yu-zhen, et al. Soil Mechanics[M]. 2nd ed. Beijing: Tsinghua university Press, 2018. (in Chinese)

    [3] 陈祖煜. 土质边坡稳定分析: 原理·方法·程序[M]. 北京: 中国水利水电出版社, 2003.

    CHEN Zu-yu. Soil Slope Stability Analysis[M]. Beijing: China Water Power Press, 2003. (in Chinese)

    [4]

    DUNCAN J M. State of the art: limit equilibrium and finite-element analysis of slopes[J]. Journal of Geotechnical and Geoenvironmental Engineering, 1996, 122(7): 577–596. doi: 10.1061/(ASCE)0733-9410(1996)122:7(577)

    [5]

    GUO M W, LI C G, WANG S L, et al. Vector-sum method for 2D slope stability analysis considering vector characteristics of force[J]. International Journal of Geomechanics, 2019, 19(6): 04019058. doi: 10.1061/(ASCE)GM.1943-5622.0001436

    [6] 孙建生. 基于有限元应力载荷宏观刚性滑裂面极限平衡抗滑稳定计算理论[J]. 岩石力学与工程学报, 2018, 37(4): 862–875. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201804009.htm

    SUN Jian-sheng. Theoretical model of stability calculation of macro-rigid sliding planes with FEM stress load based on limit equilibrium[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(4): 862–875. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201804009.htm

  • 期刊类型引用(0)

    其他类型引用(4)

图(5)  /  表(1)
计量
  • 文章访问数:  173
  • HTML全文浏览量:  30
  • PDF下载量:  138
  • 被引次数: 4
出版历程
  • 收稿日期:  2022-03-01
  • 网络出版日期:  2022-09-22
  • 刊出日期:  2022-08-31

目录

/

返回文章
返回