Zonal damage information and critical failure identification of CT images of rock under triaxial compression
-
摘要: 岩石三轴压缩试验过程中虽然损伤不断积累,但宏观破坏前其CT图像却没有明显的变化,因此在峰值前仅通过CT数难以准确分析其损伤破坏信息。为了识别岩石的损伤信息及临界破坏特征,将轴压过程中的砂岩CT图像进行分区,使用位置线上灰度的变化代表分区内灰度的变化。CT图像中某一位置线上灰度值的变化曲线是一种起伏、粗糙的曲线,利用描述粗糙曲线的分形Rd指标定量分析砂岩损伤规律,并基于CT图像上全部位置线的Rd值总和来定义损伤变量,用其表征砂岩的损伤程度。计算结果表明荷载水平在0.84左右时,Rd值波动大,损伤变量加速增大,表明损伤进入了加速发展阶段,岩石处于损坏破坏的临界状态。由此可以合理地预测岩石的损伤破坏。Abstract: The CT image and value do not change obviously before the macro-failure under triaxial compression experiments on rock during damage accumulation, so it is difficult to analyze accurately the damage information only through the change of CT number. In order to identify the damage information and critical failure characteristics of sandstone, the CT image is partitioned, and the change of grey level on the position line is used to represent the change of the gray level in the subarea. The grey value curve is a rough curve representing grey change of a certain location online in the CT image. It is used to describe quantitatively the rules of sandstone damage by fractal Rd index, the damage variable is defined as the sum of Rd value of all the position line in the CT image, and it reflects the damage degree of rock. The calculated results show that when the load is about 0.84, the value of Rd has obvious fluctuations and the damage variable increases rapidly, which denotes that the damage has entered a phase of accelerated development, and the rock is in a critical damage state. The proposed method can be used to predict reasonably the damage and failure of rock.
-
Keywords:
- sandstone /
- CT image /
- fractal index /
- critical damage information
-
0. 引言
中国是世界上膨胀土分布面积最广的国家之一,膨胀土分布面积超过10×104 km2,在膨胀土地区进行工程建设常常遇到“逢堑必滑,无堤不塌”的问题,对工程的危害是极大的,因此,进行膨胀土高边坡加固对确保工程安全稳定至关重要[1-2]。目前,国内外对膨胀土边坡的加固措施进行了深入的研究,常用的措施有边坡压脚、坡面防护、抗滑桩、挡墙、锚杆(索)、土钉等,针对不同的破坏特征和破坏模式,选用对应的加固措施是行之有效的,然而传统的加固方法也存在一定的不足,如抗滑桩加固方案工程造价高、圬工量大、施工周期长、施工场地大;锚(杆)索、预应力锚索等加固方案受注浆体或外部结构物凝期要求,施工工期长,质量不易保证;土钉(锚)等存在加固深度有限、锚固力不宜控制等问题。
针对膨胀土高边坡传统加固方法的不足,本文以南水北调中线工程淅川段膨胀土高渠坡抢险加固工程为背景,对滑坡失稳情况和原因进行分析,对变形体滑裂面位置和方向进行了推测;根据滑裂面的变形特征,提出采用新型伞型锚技术快速加固该滑坡的设计方案和施工技术要求,并通过伞型锚加固后的安全监测、加固工期等,验证了伞型锚技术在该滑坡除险加固工程中的技术优势。
1. 新型伞型锚技术
新型伞型锚主要由锚头、连接杆、锁定装置和承压板组成,见图 1所示,其中锚头是伞型锚的核心部件,主要组成部件包括锚板、击入头、滑移连杆、承压板等,见图 2所示。其加固原理[3-4]为钻机成孔后将伞型锚的锚头下至孔底,通过张拉设备施加张拉力进行张拉,通过系列装置使锚钉板呈伞形打开并嵌入两侧土体,达到张开角度要求后,通过限位装置使得锚头与两侧土体形成整体,以提供足够的锚固力,伞型锚锚力增长过程典型曲线见图 3所示。
由图 3可知:伞型锚是一种预应力锚,锚力增长经过了锚固力聚集、锚固力急剧积累和锚固力调整积累3个阶段。与传统的锚杆锚固技术相比,具有锚固力大、可即时提供锚力、不需注浆(永久加固时可通过注浆进行防腐处理)、施工简单快速、对环境影响小等优点。
2. 项目概况及滑坡原因分析
2.1 渠坡概况
该边坡为南水北调中线工程淅川段挖方膨胀土渠坡,最大坡高35.0 m,主要地层为第四系中更新统(al-plQ2)粉质黏土、黏土以及钙质结核粉质黏土组成。渠道原设计方案:第一级渠坡(过水断面)坡比1︰3.0,一级马道宽5 m为检修平台,一级马道以上每隔6 m设置一级马道,马道宽为2 m,二至四级渠坡坡比为1︰2.5,渠坡典型断面及地质分层情况见图 4所示,该段渠坡最上面两级边坡自2016年7月开始出现裂缝、变形。
2.2 滑坡变形情况分析
为确定边坡裂缝深度及边坡渗水情况,在渠坡坡顶、四级边坡1/3坡高及二、三级马道近坡脚位置进行了探槽开挖,探槽情况见表 1所示。
表 1 探坑情况表Table 1. Situations of exploring pits编号 位置 描述 WK1 距坡肩10 m 无裂缝,土体潮湿,无明水 WK2 距坡肩5 m 有裂缝,延伸至坑底23 cm,坑内干燥 WK3 距3级马道2.5 m 有裂缝,延伸至坑底10 cm, 坑底潮湿 WK4 4级边坡坡脚 无裂缝,无明显滑动迹象,有渗水 WK5 3级边坡坡脚 坑内无裂缝,坑底有少量渗水 现场查勘和探槽发现:渠坡二级马道土体下沉,压顶板与土体脱空,二级马道排水沟内侧沟壁与底板出现脱空现象,排水沟壁持续缩窄,缩窄速率为1 cm/月,目测有向渠道内位移迹象;坡顶距坡肩5 m位置和三级渠坡坡脚裂缝发育,且裂缝长度和宽度呈增大趋势,三级边坡坡脚有水渗出。
加固前的测斜发现:一级马道处渠坡深层水平位移几乎为零;二级马道渠坡深层水平位移值较小,加固前最大值为2.2 mm,拐点发生在坡面下1~2 m;三级马道渠坡深层水平位移值较大,加固前最大值达19.5 mm,且位移变化速率呈增大趋势,拐点发生在坡面下6~8 m处。
根据现场查勘、探槽、侧斜数据可判断,该渠坡三级边坡内部存在近似水平的滑裂面,且边坡变形在持续发生,无收敛趋势,因此,为确保渠道边坡安全稳定,必须进行堤坡的加固处理。
2.3 膨胀土边坡破坏模式及滑裂面推测
为了合理、准确确定膨胀土边坡加固措施和加固方案,确保边坡加固效果,在设计前应先进行边坡破坏模式分析,推测边坡潜在的滑裂面。
根据文献[5]膨胀土边坡的滑动模式主要有浅层滑动和深层滑动两种。程永辉等[6]通过离心模型试验揭示了浅层滑动的关键因素是膨胀变形,认为土体含水率变化会导致膨胀土体产生膨胀变形,膨胀变形引起边坡土体应力重分布,进而产生的顺坡向剪应力,产生浅层滑动;龚壁卫等[7]、胡波等[8]通过大量现场试验、大型静力模型试验、土工试验结合数值分析揭示了深层滑动的关键因素是裂隙面空间分布及其强度,认为膨胀土的裂隙面是具有一定厚度或充填物的软弱层面,其强度远远低于周围土体,成为膨胀土体中力学性能薄弱的一环,一旦地层走向合适,边坡开挖形成后,将产生塑性破坏,并随着变形发展,最终裂隙逐渐贯通,产生沿裂隙面的滑动。
根据二级马道和三级马道测斜实测数据,深层水平位移拐点均发生在二级马道以下1~2 m位置,根据勘察资料,所处地层为含铁锰石斑块的黏土层,偶见姜石,有存在裂隙面的可能,初步判断该渠坡滑动模式为裂隙控制的深层滑坡,结合渠坡坡面变形和裂缝分布情况,初步推测滑裂面位置如图 4所示。
3. 除险加固方案
3.1 加固方案
根据推测滑裂面,经方案比选,拟采用伞型锚加排水管的处理措施进行渠坡加固。根据类似工程经验,共布置伞型锚三排,单根伞型锚设计锚固力为120 kN,锚固方向与水平面角度为35°。第一排伞型锚沿坡面(下同)距二级马道1 m,锚固长度为15 m;第二排伞型锚距二级马道4 m,锚固长度为15 m;第三排伞型锚距二级马道7 m,锚固长度为20 m。三排伞型锚呈矩形布置,沿渠道水流方向间距为3 m。排水管布置为现有出水明显部位间距10 m,其它部位间距20 m。经稳定分析,加固后边坡的稳定系数能够满足规范要求。
3.2 施工技术要求
伞型锚现场施工仅用9 d时间(1月9日—17日),共施工伞型锚121根,施工速率为13~14根/天;渠坡加固长度120 m,加固速率为沿坡肩方向13.3 m/d,现场加固速度非常快,施工完成后伞型锚即可马上受力,及时阻止了渠坡变形进一步发生。
主要施工技术要求:
(1)钻孔首次采用履带潜孔钻车,该钻车在坡面上一定倾角钻孔非常方便、速度快,效果好。
(2)该段渠坡为永久加固工程,现场施工时,应对锚头和连接杆进行防腐处理,本工程采用灌注水泥浆对伞型锚及锚杆进行防腐处理。
(3)伞型锚张拉施工应按照相关规程规范[9-10]要求进行,张拉过程中记录拉拔力、锚杆上拔位移等,张拉至设计荷载且位移稳定后即可停止张拉。
(4)伞型锚锚固力锁定值应根据地层条件、变形要求确定[11]。
3.3 加固效果分析
深层水平位移采用测斜观测,布置测斜管5孔,孔深25 m,其中一级、二级马道各一孔,三级马道3孔,监测时间为2016年10月18日—2017年5月10日;锚杆拉力采用钢筋计监测,布置钢筋计9支,布置在3个断面,每个断面上、中、下三排锚各安装1支,监测时间为2017年1月17日—5月10日。测斜监测结果可知,截止2017年5月10日,一级马道和二级马道处测斜管监测坡体水平位移值很小,可认为渠道一级边坡和二级边坡坡体是稳定的。三级马道处测斜管监测坡体最大水平位移值26.5 mm,典型测斜曲线见图 5所示,测斜最大值变化曲线见图 6所示,测斜最大值变化速率图见7所示。
由图 6,7可知,加固完成时,已发生的累计位移值为19.5 mm;加固前深层水平位移最大值变形速率呈逐渐增大趋势,最大变形速率为0.65 mm/d,发生在1月8日;伞型锚加固完成后,变形速率逐渐趋缓,截止5月10日,变形速率为0.01 mm/d,可认为渠坡变形基本稳定。
锚力随时间变化曲线见图 8所示。由图 8可知,下排锚提供锚固力最大、中排锚次之、上排锚最小。下排锚单根锚固力43.8~93.4 kN,平均值69.6 kN,与初值相比锚固力增加98.9%;中排单根锚固力39.3~50.4 kN,平均值44.1 kN,与初值相比锚固力增加50.3%;上排锚提供锚固力最小,单根锚固力17.5~31.1 kN,平均值22.3 kN,与初值相比锚固力增加10.2%。由锚力变化曲线可知:下、中排锚锚固力先呈逐渐增大趋势,之后锚固力趋于定值,最大值发生在4月20日左右;上排伞型锚锚固力变化幅度不大,加固后基本趋于定值。结合测斜观测数据可知,锚力的逐渐增大阻止了渠坡变形的发生,当渠坡变形趋于稳定时,锚固力也逐渐趋于定值,说明伞型锚加固有效阻止了渠坡变形的进一步发生,加固效果良好。
4. 结论
(1)伞型锚技术是一种新型预应力锚固技术,具有锚固力大、施工便捷、施工速度快、可即时发挥锚固作用、对环境无污染等优点,可广泛应用在边坡加固、基坑支护、基坑抗浮、应急抢险等工程中,应用前景广阔。
(2)滑裂面的准确推测对加固方案的确定至关重要,本边坡综合考虑了坡面裂缝分布、现场探槽及坡体变形情况,结合地质资料,推测了滑裂面位置和滑动方向,加固后的边坡稳定分析和实测数据验证了推测滑裂面的合理性。
(3)伞型锚加固后,现场实测数据显示:边坡最大水平位移变形速率呈减缓趋势,变形逐渐趋于收敛;伞型锚的锚固力随水平位移增大逐渐增大,锚固力的增大又阻止了边坡变形的发生,当渠坡变形趋于稳定时,锚固力也逐渐趋于定值,监测结果表明采用伞型锚加固效果良好。
(4)伞型锚技术在南水北调中线工程淅川段膨胀土高渠坡抢险加固工程中的成功应用,充分验证了伞型锚在边坡加固工程中的技术可行性;快速的现场施工和施工完成即可发挥锚固力的特点,再次验证了伞型锚技术在膨胀土高边坡除险加固工程中的技术优势。
-
表 1 加载过程中的H值和SD值
Table 1 Values of H and SD during loading
扫描号 轴向应力/MPa 轴向应变/% 第20扫描层 第30扫描层 H SD H SD 1 0 0 2207.4 37.6 2192.0 40.9 2 5.82 0.25 2208.8 40.7 2191.7 39.9 3 12.71 0.43 2208.8 38.9 2193.9 40.4 4 24.53 0.76 2209.2 40.0 2192.4 41.6 5 32.14 1.01 2207.1 44.2 2190.5 45.7 6 38.24 1.26 2204.2 50.2 2182.7 48.3 7 25.61 1.67 2177.6 123.3 2168.4 70.7 8 24.36 1.83 2172.5 142.3 2165.9 79.8 表 2 第20扫描层各位置线上Rd值
Table 2 Values of Rd of different position lines
扫描次序 荷载水平/% 轴向应力/MPa Rd 15° 45° 75° 105° 135° 165° 总和 1 0 0 94.05 98.36 100.15 94.77 94.19 95.25 576.77 2 0.15 5.82 90.10 94.52 93.78 93.20 93.03 90.20 554.83 3 0.33 12.72 85.34 95.27 90.81 86.61 92.91 89.70 540.64 4 0.64 24.54 85.97 98.44 91.60 87.90 93.43 89.29 546.63 5 0.84 32.14 81.45 85.23 89.90 85.99 84.26 84.48 511.31 6 1.00 38.24 86.65 88.37 91.69 84.99 90.73 93.89 536.32 表 3 第30扫描层各位置线上Rd值
Table 3 Value of Rd of different position lines
扫描次序 荷载水平/% 轴向应力/MPa Rd 15° 45° 75° 105° 135° 165° 总和 1 0 0 98.26 105.05 99.05 100.75 99.77 99.39 602.27 2 0.15 5.82 93.52 98.68 94.29 95.20 91.28 98.03 571.00 3 0.33 12.72 95.07 85.71 90.34 94.70 88.61 96.91 551.34 4 0.64 24.54 95.27 96.50 90.87 94.29 91.21 94.02 562.16 5 0.84 32.14 85.23 81.14 86.35 88.48 89.80 85.26 516.26 6 1.00 38.24 88.18 82.67 73.86 78.20 90.21 87.12 500.24 表 4 第20,30扫描层的损伤变量值
Table 4 Values of damage variables for 20th and 30th scanning layers
扫描次序 荷载水平/% 损伤变量D dj 第20层 第30 第20层 第30层 第1次 0 0 0 576.77 602.27 第2次 0.15 0.23 0.27 554.83 571.00 第3次 0.33 0.38 0.49 540.64 551.34 第4次 0.64 0.40 0.52 546.63 562.16 第5次 0.84 0.84 0.86 511.31 516.26 第6次 1.00 1 1 536.32 500.24 -
[1] 李树刚, 陈高峰, 双海清, 等. 加载速率和初始损伤对砂岩能量演化影响的试验研究[J]. 采矿与安全工程学报, 2019, 36(2): 373-380. https://www.cnki.com.cn/Article/CJFDTOTAL-KSYL201902022.htm LI Shu-gang, CHEN Gao-feng, SHANG Hai-qing, et al. Experimental study on effect of loading rate and initial damage on energy evolution of sandstone[J]. Journal of Mining & Safety Engineering, 2019, 36(2): 373-380. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KSYL201902022.htm
[2] MEES F, SWENNEN R, VAN Geet M, et al. Applications of X-ray computed tomography in the geosciences[J]. Geological Society, London, Special Publications, 2003, 215(1): 1-6. doi: 10.1144/GSL.SP.2003.215.01.01
[3] MARTÍNEZ-MARTÍNEZ J, FUSI N, GALIANA-MERINO J J, et al. Ultrasonic and X-ray computed tomography characterization of progressive fracture damage in low-porous carbonate rocks[J]. Engineering Geology, 2016, 200: 47-57. doi: 10.1016/j.enggeo.2015.11.009
[4] 杨更社, 刘慧. 基于CT 图像处理技术的岩石损伤特性研究[J]. 煤炭学报, 2007, 32(5): 463-468. doi: 10.3321/j.issn:0253-9993.2007.05.004 YANG Geng-she, LIU Hui. Study on the rock damage characteristics based on the technique of CT image processing[J]. Journal of China Coal Society, 2007, 32(5): 463-468. (in Chinese) doi: 10.3321/j.issn:0253-9993.2007.05.004
[5] 杨更社, 孙钧, 谢定义. 岩石材料损伤变量与CT数间的关系分析[J]. 力学与实践, 1998, 20(4): 47-49. https://www.cnki.com.cn/Article/CJFDTOTAL-LXYS804.015.htm YANG Geng-she, SUN Jun, XIE Ding-yi. Analysis of the relation between damage variable and CT value of rock material[J]. Mechanics in Engineering, 1998, 20(4): 47-49. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-LXYS804.015.htm
[6] 张国凯, 李海波, 王明洋, 等. 岩石单轴压缩下损伤表征及演化规律对比研究[J]. 岩土工程学报, 2019, 41(6): 1074-1082. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201906013.htm ZHANG Guo-kai, LI Hai-bo, WANG Ming-yang, et al. Comparative study on damage characterization and damage evolution of rock under uniaxial compression[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(6): 1074-1082. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201906013.htm
[7] 任建喜, 惠兴田. 裂隙岩石单轴压缩损伤扩展细观机理CT分析初探[J]. 岩土力学, 2005, 26(增刊1): 18-23. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2005S1010.htm REN Jian-xi, HUI Xing-tian. Primary study on meso-damage propagation mechanism of cracked-sandstone using computerized tomography under uniaxial compression[J]. Rock and Soil Mechanics, 2005, 26(S1): 18-23. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2005S1010.htm
[8] REN J X, LUO Y, LIU W G. Application of computerized topography testing technology on studying rock failure mechanism under loading and unloading[J]. Journal of Glaciology and Geocryology, 2002, 24(5): 672-675. doi: 10.3969/j.issn.1000-0240.2002.05.036
[9] YU Q, YANG S, RANJITH P G, et al. Numerical modeling of jointed rock under compressive loading using X-ray computerized tomography[J]. Rock Mechanics and Rock Engineering, 2016, 49(3): 877-891. doi: 10.1007/s00603-015-0800-4
[10] 曹文贵, 林星涛, 张超, 等. 基于非线性动态强度准则的岩石动态变形过程统计损伤模拟方法[J]. 岩石力学与工程学报, 2017, 36(4): 794-802. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201704003.htm CAO Wen-gui, LIN Xing-tao, ZHANG Chao, et al. A statistical damage simulation method of dynamic deformation process for rocks based on nonlinear dynamic strength criterion[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(4): 794-802. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201704003.htm
[11] 丁卫华, 仵彦卿, 蒲毅彬. 受力岩石密度损伤增量及其数字图像[J]. 西安理工大学学报, 2000, 16(1): 61-64. https://www.cnki.com.cn/Article/CJFDTOTAL-XALD200001013.htm DING Wei-hua, WU Yan-qing, PU Yi-bin. The density damage increment and its digital image of rock in compression[J]. Journal of Xi'an University of Technology, 2000, 16(1): 61-64. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XALD200001013.htm
[12] 丁卫华, 仵彦卿, 蒲毅彬. 低应变率下岩石内部裂纹演化的X射线CT方法[J]. 岩石力学与工程学报, 2003, 22(11): 1793-1797. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200311007.htm DING Wei-hua, WU Yan-qing, PU Yi-bin. X-ray CT approach on rock-interior crack evolution under low strain rate[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(11): 1793-1797. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200311007.htm
[13] 党发宁, 方建银, 丁卫华. 基于CT的混凝土试样静动力单轴拉伸破坏裂纹分形特征比较研究[J]. 岩石力学与工程学报, 2015(增刊1): 2922-2928. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2015S1041.htm DANG Fa-ning, FANG Jian-yin, DING Wei-hua. Fractal comparison research of fracture of concrete samples under static and dynamic uniaxial tensile using CT[J]. Chinese Journal of Rock Mechanics and Engineering, 2015(S1): 2922-2928. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2015S1041.htm
[14] 易成, 张亮, 陈忠辉, 等. 一种新的描述粗糙表面形貌尺度分维参数的研究[J]. 中国矿业大学学报, 2007, 36(1): 75-80. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD200701015.htm YI Cheng, ZHANG Liang, CHEN Zhong-hui, et al. A novel description of roughness surface with a modified fractal index Rd[J]. Journal of China University of Mining & Technology, 2007, 36(1): 75-80. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD200701015.htm
[15] 李夕兵, 翁磊, 谢晓锋, 等. 动静载荷作用下含孔洞硬岩损伤演化的核磁共振特性试验研究[J]. 岩石力学与工程学报, 2015, 34(10): 1985-1993. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201510005.htm LI Xi-bing, WENG Lei, XIE Xiao-feng, et al. Study on the degradation of hard rock with a pre-existing opening under static-dynamic loadings using nuclear magnetic resonance technique[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(10): 1985-1993. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201510005.htm
[16] 曾鹏, 刘阳军, 纪洪广, 等. 单轴压缩下粗砂岩临界破坏的多频段声发射耦合判据和前兆识别特征[J]. 岩土工程学报, 2017, 39(3): 509-517. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201703021.htm ZENG Peng, LIU Yang-jun, JI Hong-guang, et al. Coupling criteria and precursor identification characteristics of multi-band acoustic emission of gritstone fracture under uniaxial compression[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(3): 509-517. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201703021.htm
[17] HOUNSFIELD G N. Computerized transverse axial scanning (tomography) part 1: description of system[J]. The British Journal of Radiology, 1973, 46(552): 1016-1022.
-
期刊类型引用(2)
1. 张远庆,陈勇,王世梅,王力. 岸坡渗流潜蚀模型试验系统变革研究. 三峡大学学报(自然科学版). 2025(02): 48-54 . 百度学术
2. 何健健,蒋希豪,汪玉冰. 离心模型试验中温度及孔隙率对砂土渗透系数的影响研究(英文). Journal of Zhejiang University-Science A(Applied Physics & Engineering). 2025(03): 177-194 . 百度学术
其他类型引用(1)