高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

水辅助激光无重铸层钻孔Al2O3陶瓷实验研究

周翔 段军 陈航 张菲 白克强 韩小花

引用本文:
Citation:

水辅助激光无重铸层钻孔Al2O3陶瓷实验研究

    作者简介: 周翔(1992-), 男, 硕士研究生, 现主要从事硬脆性材料激光精密制造方面的研究.
    通讯作者: 段军, duans@mail.hust.edu.cn
  • 基金项目:

    国家自然科学基金资助项目 51475182

    国家自然科学基金资助项目 51675205

  • 中图分类号: TN249

Experimental study about water-assisted laser drill on Al2O3 ceramics without recast layer

    Corresponding author: DUAN Jun, duans@mail.hust.edu.cn ;
  • CLC number: TN249

  • 摘要: 为了解决传统加工过程中重铸层无法消除的问题,采用超快皮秒激光创新性地在水介质中对Al2O3陶瓷进行皮秒激光钻孔实验,并与空气中钻孔结果进行对比,研究了皮秒激光主要参量如单脉冲能量、扫描次数等对陶瓷微孔的孔径、锥度和重铸层厚度的影响规律,并分析讨论不同介质中皮秒激光与Al2O3陶瓷相互作用机理及材料去除机制。结果表明,在水介质中激光钻孔直径增加约35μm、微孔锥度降低至1.04°并可获得无重铸层钻孔效果;激光作用过程中水的存在会引起空泡作用、吸收作用和运输作用,有效防止了去蚀材料二次黏附,消除了重铸层和降低了微孔锥度,提升了微孔质量。该研究阐述了水辅助激光钻孔的具体影响状况并加深了对水辅助的影响机理理解。
  • Figure 1.  Schematic diagram of micro-machining system by picosecond laser

    Figure 2.  Test method for micro taper

    Figure 3.  Measurement method of recast-layer thickness

    Figure 4.  Relationship between absorption and wavelength in pure water

    Figure 5.  Effect of laser parameters on hole diameter and taper of alumina ceramics

    a—hole diameter in air b—hole diameter in water c—taper in air d—taper in water

    Figure 6.  Micrograph of hole

    a~d—drilling in air e~h—drilling in water

    Figure 7.  Scanning electron microscope picture of the drilled sample

    a, b—in air c, d—in water

    Table 1.  Physical parameters of soda-lime glass

    properties reference value
    melting point 2072℃
    boiling point 2977℃
    density 3.72g/cm3
    heat capacity 880J/(kg·K)
    elasticity modulus 74GPa
    shear modulus 300GPa
    coefficient of thermal expansion 8.2×10-6K-1
    下载: 导出CSV
  • [1]

    ZHU M X, XU X, LI Y M, et al. Research status on metal-toughed alumina ceramics[J]. Chinese Ceramics, 2012, 148(12):10-13(in Chinese). 
    [2]

    ZHANG F L, FENG D J, SHI J C, et al. Rotary ultrasonic machining technology of ceramic material[J]. Electromachining & Mould, 2001(1):1-5(in Chinese). 
    [3]

    QI L T. Ultra-short pulsed laser microfabrication technology[M]. Harbin:Harbin Engineering University Press, 2012:22-25(in Chin-ese).
    [4]

    SCHULZ W, EPPELT U, POPRAWE R. Review on laser drilling Ⅰ. Fundamentals, modeling, and simulation[J]. Journal of Laser Applications, 2013, 25(1):12006. doi: 10.2351/1.4773837
    [5]

    ZHANG F, DUAN J, ZENG X, et al. UV laser microprocessing and post chemical etching on ultrathin Al2O3 ceramic substrate[J]. Journal of the European Ceramic Society, 2011, 31(9):1631-1639. doi: 10.1016/j.jeurceramsoc.2011.03.034
    [6]

    ADELMANN B, HELLMANN R. Rapid micro hole laser drilling in ceramic substrates using single mode fiber laser[J]. Journal of Materials Processing Technology, 2015, 221(1):80-86. 
    [7]

    HSU J, LIN W, CHANG Y, et al. Continuous-wave laser drilling assisted by intermittent gas jets[J]. The International Journal of Advanced Manufacturing Technology, 2015, 79(1/4):449-459. 
    [8]

    ZHANG F. Study on UV laser microprocessing technology and mechanism for electronic materials[D]. Wuhan: Huazhong University of Science and Technology, 2012: 87-102(in Chinese).
    [9]

    ZHANG G, LIANG Y. Study on laser drilling alumina ceramic of recast layer research[J]. Journal of ShenYang Institute of Technology, 2001, 20(2):1-4(in Chinese).
    [10]

    LI C, LEE S, NIKUMB S. Femtosecond laser drilling of alumina wafers[J]. Journal of Electronic Materials, 2009, 38(9):20116-20121. 
    [11]

    WANG X C, ZHENG H Y, CHU P L, et al. Femtosecond laser drilling of alumina ceramic substrates[J]. Applied Physics, 2010, A101(2):271-278. 
    [12]

    KONG L R, ZHANG F, DUAN J, et al. Research of water-assisted laser etching of alumina ceramics[J]. Laser Technology, 2014, 38(3):330-334(in Chinese). 
    [13]

    TSAI C, LI C. Investigation of underwater laser drilling for brittle substrates[J]. Journal of Materials Processing Technology, 2009, 209(6):2838-2846. doi: 10.1016/j.jmatprotec.2008.06.057
    [14]

    KRUUSING A. Underwater and water-assisted laser processing:Part 1-general features, steam cleaning and shock processing[J]. Optics and Lasers in Engineering, 2004, 41(2):307-327. doi: 10.1016/S0143-8166(02)00142-2
    [15]

    LIU Y Zh, ZHU Y T, ZHU G X, et al. Study of underwater spot spread of collimated laser beam[J]. Chinese Journal of Lasers, 2001, 28(7):617-620(in Chinese). 
  • [1] 孔令瑞张菲段军罗瑞峰曾晓雁 . 水辅助激光刻蚀氧化铝陶瓷的研究. 激光技术, 2014, 38(3): 330-334. doi: 10.7510/jgjs.issn.1001-3806.2014.03.010
    [2] 王成曾晓雁 . Al2O3陶瓷的激光三维雕刻实验研究. 激光技术, 2007, 31(1): 18-21.
    [3] 周敏魏昕谢小柱胡伟 . 355nm紫外激光抛光Al2O3陶瓷工艺的研究. 激光技术, 2014, 38(4): 556-560. doi: 10.7510/jgjs.issn.1001-3806.2014.04.024
    [4] 郭商勇陈涛刘世炳 . 提高准分子激光打孔质量的方法研究. 激光技术, 2006, 30(6): 625-627.
    [5] 谢兵兵袁根福 . 水射流激光复合刻蚀陶瓷的研究. 激光技术, 2016, 40(5): 762-766. doi: 10.7510/jgjs.issn.1001-3806.2016.05.030
    [6] 龙芋宏熊良才史铁林 . 水辅助准分子激光微加工硅的实验研究. 激光技术, 2006, 30(6): 567-569.
    [7] 刘鑫龙芋宏鲍家定刘清原毛建冬 . 基于水辅助激光加工的水层流动特性的研究. 激光技术, 2017, 41(3): 442-446. doi: 10.7510/jgjs.issn.1001-3806.2017.03.027
    [8] 杨元政刘志国刘正义庄育智 . 等离子喷涂Al2O3陶瓷涂层激光熔化深度的研究. 激光技术, 2000, 24(5): 312-317.
    [9] 王昆孙桂芳张尉冯爱新袁国定 . 激光合金化Mn-Al2O3和Mn-Al2O3-NiWC涂层的磨蚀性能研究. 激光技术, 2015, 39(2): 274-279. doi: 10.7510/jgjs.issn.1001-3806.2015.02.027
    [10] 张光辉黄宇星黄平周辽焦辉龙芋宏 . 水导激光技术中水-光耦合传能规律研究. 激光技术, 2022, 46(6): 749-754. doi: 10.7510/jgjs.issn.1001-3806.2022.06.006
    [11] 庞继伟王超蔡玉奎 . 玻璃材料激光加工技术的研究进展. 激光技术, 2021, 45(4): 417-428. doi: 10.7510/jgjs.issn.1001-3806.2021.04.003
    [12] 李策赵培娥彭涛冯力天周杰罗雄周鼎富 . 3维激光测风雷达技术研究. 激光技术, 2017, 41(5): 703-707. doi: 10.7510/jgjs.issn.1001-3806.2017.05.017
    [13] 杨航刘晓东李君豪卢敏健 . 3维曲面振镜激光加工FPLSCM算法实践. 激光技术, 2019, 43(4): 482-487. doi: 10.7510/jgjs.issn.1001-3806.2019.04.009
    [14] 段军 . 激光微加工磁盘——激光毛化技术现状与发展. 激光技术, 2006, 30(5): 490-493.
    [15] 宋新华修腾飞金湘中袁江宋斌 . 激光辅助加热搅拌摩擦焊3维流场数值模拟. 激光技术, 2016, 40(3): 353-357. doi: 10.7510/jgjs.issn.1001-3806.2016.03.011
    [16] 刘泽宇魏昕谢小柱华显刚洪继伟 . 紫外激光加工陶瓷刀具表面微织构的实验研究. 激光技术, 2016, 40(4): 550-554. doi: 10.7510/jgjs.issn.1001-3806.2016.04.020
    [17] 张冲王冠刘赞丰张雅文 . 激光微加工对Ti6Al4V表面形貌及润湿性影响的研究. 激光技术, 2021, 45(1): 31-36. doi: 10.7510/jgjs.issn.1001-3806.2021.01.006
    [18] 李东华邓磊敏段军曾晓雁游欣易 . 激光大尺度3维动态聚焦扫描加工系统研究. 激光技术, 2016, 40(4): 466-471. doi: 10.7510/jgjs.issn.1001-3806.2016.04.002
    [19] 叶昌庚闫平欧攀巩马理 . 基于CO2激光的双包层光纤端帽熔接实验研究. 激光技术, 2007, 31(5): 456-458.
    [20] 吴超杰陈培锋王英夏兵兵 . 连续CO2激光在线膏药群孔加工系统的研究. 激光技术, 2016, 40(3): 422-425. doi: 10.7510/jgjs.issn.1001-3806.2016.03.026
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  6397
  • HTML全文浏览量:  3919
  • PDF下载量:  296
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-05-17
  • 录用日期:  2017-06-27
  • 刊出日期:  2018-03-25

水辅助激光无重铸层钻孔Al2O3陶瓷实验研究

    通讯作者: 段军, duans@mail.hust.edu.cn
    作者简介: 周翔(1992-), 男, 硕士研究生, 现主要从事硬脆性材料激光精密制造方面的研究
  • 华中科技大学 武汉光电国家实验室, 武汉 430074
基金项目:  国家自然科学基金资助项目 51475182国家自然科学基金资助项目 51675205

摘要: 为了解决传统加工过程中重铸层无法消除的问题,采用超快皮秒激光创新性地在水介质中对Al2O3陶瓷进行皮秒激光钻孔实验,并与空气中钻孔结果进行对比,研究了皮秒激光主要参量如单脉冲能量、扫描次数等对陶瓷微孔的孔径、锥度和重铸层厚度的影响规律,并分析讨论不同介质中皮秒激光与Al2O3陶瓷相互作用机理及材料去除机制。结果表明,在水介质中激光钻孔直径增加约35μm、微孔锥度降低至1.04°并可获得无重铸层钻孔效果;激光作用过程中水的存在会引起空泡作用、吸收作用和运输作用,有效防止了去蚀材料二次黏附,消除了重铸层和降低了微孔锥度,提升了微孔质量。该研究阐述了水辅助激光钻孔的具体影响状况并加深了对水辅助的影响机理理解。

English Abstract

    • Al2O3陶瓷是目前世界上应用最为广泛的陶瓷材料之一,在生物、电子等领域已成为不可或缺的材料[1]。传统的机械加工、电火花加工和超声波加工Al2O3陶瓷等方法,都存在无法钻微小孔、钻孔效率低和工件磨损严重等问题[2]。自从激光器诞生以来,激光的种类和利用激光进行材料加工已经获得蓬勃发展[3]。激光应用领域逐步扩大到各种金属和非金属板材的加工,解决了许多传统切削加工无法解决的问题[4]。微孔激光加工技术是激光加工技术的一个新兴和重要领域,随着科技和社会生产的迅速发展,一方面给激光钻孔提出了各种各样更高的要求;另一方面技术的进步和先进设备的出现也使得高效率、高质量钻孔成为可能[5]。目前微孔加工领域主要存在的问题是:在加工微孔过程中,不仅要使材料易于加工,还应努力提高微孔的钻孔质量,减少微孔锥度和消除重铸层、微裂纹等问题。

      目前,国内外学者相继采用不活泼气体辅助激光加工[6-7]、激光钻孔样品化学再腐蚀[8-9]和超短脉冲激光加工[10-11]等手段来去除重铸层,但效果均不完美,未能完全消除重铸层,主要原因在于上述手段不能完全避免热效应的存在,依然有熔融物重新附着于原材料表面。而水辅助激光钻孔由于激光诱导液体空化作用、水的冷却作用和水的流动作用能够在一定程度上降低光热效应和抑制去蚀材料的二次黏附。曾经有学者[12-13]在纳秒、毫秒等激光器下利用水辅助激光切割陶瓷,虽未能完全消除热效应和重铸层,但加工质量明显优于空气中激光加工。

      本文中在前人的基础上采用水辅助超短脉冲激光对Al2O3陶瓷进行钻微孔研究,对比水和空气介质中钻孔的直径、微孔锥度和重铸层厚度等质量效果,研究水辅助激光钻孔的影响机理。

    • 实验中所用Al2O3陶瓷3维尺寸为28mm×21mm×625μm,Al2O3的质量分数为0.96,能带宽度约为9.1eV,常用物理性能见表 1

      Table 1.  Physical parameters of soda-lime glass

      properties reference value
      melting point 2072℃
      boiling point 2977℃
      density 3.72g/cm3
      heat capacity 880J/(kg·K)
      elasticity modulus 74GPa
      shear modulus 300GPa
      coefficient of thermal expansion 8.2×10-6K-1
    • 实验中采用的激光加工设备为多功能皮秒激光微加工系统,主要由多波段全固态皮秒激光器、微加工光路系统、精密机械系统、定位系统和计算机控制系统构成。激光器是德国Edgewave公司生产的Nd:YVO4皮秒激光器,输出波长1064nm,光束质量因子M2≤1.3。激光输出平均功率为0W~80W,脉冲宽度10ps,脉冲重复频率为200kHz~20MHz,出口光斑直径3.0mm,经透镜聚焦后的焦点光斑直径约为20.0μm。使用振镜扫描配合加工平台X-Y-Z 3维运动方式,可实现对70mm×70mm范围内任意复杂图形快速加工,最高扫描速度为9600mm/s。加工系统示意图如图 1所示。

      Figure 1.  Schematic diagram of micro-machining system by picosecond laser

    • 实验研究激光主要参量(单脉冲能量和扫描次数)对微孔直径、锥度和重铸层厚度影响规律以及皮秒激光在不同介质中与Al2O3陶瓷的作用机理。钻微孔目标直径D=200μm,钻孔方式采取单层多次同心圆填充旋切钻孔法,填充间距固定为10μm。激光聚焦平面位于Al2O3陶瓷表面上,焦点处有效光斑中心位于同心圆上。在空气介质中加工只需将Al2O3陶瓷样品放在加工平台上直接按设定参量进行加工,而水辅助激光钻孔则需要支撑和固定装置将Al2O3陶瓷固定于水下。

      微孔锥度是陶瓷微孔在材料厚度方向上的倾角θ,如图 2所示。锥度T计算公式如下:

      $ T\left( \theta \right){\rm{ = }}\arctan \left( {\frac{{{D_{{\rm{ent}}}} - {D_{{\rm{ext}}}}}}{{2t}}} \right) $

      (1)

      Figure 2.  Test method for micro taper

      式中,θ为倾角,单位是°;DentDext分别为正、反面孔径,单位均是μm;t为Al2O3陶瓷厚度,单位是μm。

      重铸层主要由于材料熔化重新冷却固化在内壁上,如图 3所示。计算公式如下:

      $ {\zeta _{{\rm{HAZ}}}} = \frac{{{D_{{\rm{HAZ}}}} - {D_0}}}{2} $

      (2)

      Figure 3.  Measurement method of recast-layer thickness

      式中,ζHAZ是重铸层厚度,单位是μm;DHAZ是包括重铸层在内微孔直径,单位是μm;D0是材料完全去除时微孔直径,单位是μm。

    • 激光束在水中传输时,水对激光的吸收作用与散射作用同时存在。激光束能量的衰减主要是纯净水与水中杂质对激光的吸收引起,对激光束的扩散影响不大。激光在纯净水中的传输特性如图 4所示。在水辅助激光钻孔过程中,水能够吸收一部分激光能量并降低加工效率。水吸收的激光能量可以由Beer-Lambert定律[14]计算得到,如下式所示:

      $ {I_x}\left( \lambda \right) = {I_0}\left( \lambda \right)\exp \left( {\frac{{ - x}}{L}} \right) $

      (3)

      Figure 4.  Relationship between absorption and wavelength in pure water

      式中,I0(λ)是入射前激光初始辐照度;Ix(λ)是在液体中传输路程为x后的激光辐照度;L是激光在水中的吸收长度,即激光被完全吸收穿过的溶液长度。

      因此可以很容易看出, 波长λ=1064nm激光对应纯净水吸收长度L≈15mm。水深h=1mm时,吸收率δh/L=6.7%;水深h=2mm时, 吸收率δ≈13.3%。而散射作用主要使激光束扩散、光斑变大,对激光束能量衰减影响较弱。深度为x处的激光光斑大小S计算公式[15]如下:

      $ S\left( x \right) = {S_0} = \exp \left( {\rho \left\langle {{\sigma _{\rm{s}}}} \right\rangle x} \right) $

      (4)

      式中,ρσs〉是水质参量,S0是初始光斑直径,σs是平均散射截面。

    • 本文中对空气和水不同介质中激光加工参量对陶瓷微孔的锥度和密封性影响规律进行了研究,主要的研究参量是激光单脉冲能量和扫描次数。实验中设定的参量变化范围如下:激光单脉冲能量为75.0μJ~150.0μJ;空气介质中的扫描次数为50次~200次;水介质中的扫描次数为400次~1000次。脉冲重复频率恒定为200kHz;扫描速率恒定为200mm/s;水辅助时陶瓷表面离水面深度为2mm。微孔直径和锥度的数据曲线如图 5所示。

      Figure 5.  Effect of laser parameters on hole diameter and taper of alumina ceramics

      在上述钻孔激光参量条件下,微孔直径随着单脉冲能量和扫描次数的增加呈现增加趋势,锥度随着单脉冲能量和扫描次数的增加呈现下降趋势。水辅助激光加工孔径增大主要是水对激光散射作用致使激光光斑增加。已知空气中有效光斑直径d=20μm,中心位于最外圈同心圆上,设定钻孔直径Dent=200μm,空气介质中激光钻孔正面孔径在Dent,1≈220μm,水辅助激光钻孔正面孔径稳定在Dent,2≈250μm,则散射作用引起孔径增大,Δ2=Dent,2-(Dent+d)=30μm,从而计算出水辅助激光钻孔在最后阶段有效光斑直径d′=d+$\frac{{{\Delta _2}}}{2} $=35μm,根据(4)式推算出水介质中钻孔最后阶段对应的水质参量ρσs〉=0.28。空气介质中激光钻孔锥度分布范围为2.6°~3.6°,如图 5c所示,在单脉冲能量为150μJ和扫描次数为200次时能获得最低锥度(2.6°)。水介质中激光钻孔锥度分布范围为0.8°~2.0°,如图 5d所示,在单脉冲能量为150μJ和扫描次数为1000次时能获得最低锥度(0.8°)。

      图 6中是皮秒激光分别在空气和水介质中对陶瓷钻微孔样品的光学显微镜剖面微观形貌图。图 6a~图 6d是在空气介质中激光单脉冲能量范围在75μJ~150μJ,增量为25μJ,扫描次数50次所获得的微孔,在微孔顶部有明显的热影响区,存在较大的锥度。图 6e~图 6h是水介质中激光单脉冲能量范围在75μJ~150μJ,增量为25μJ,扫描次数600次所获得的微孔。随着单脉冲能量的增加,水辅助钻孔的孔型逐渐向圆柱状转变,锥度明显小于空气中激光钻孔,并且空气中激光钻孔在微孔顶部有明显的热影响区,水辅助激光钻孔则没出现热影响区,内壁颜色与陶瓷原表面相同。

      Figure 6.  Micrograph of hole

      空气介质中皮秒激光钻孔内壁重铸层厚度变化范围为3.5μm~8.0μm,重铸层厚度随着单脉冲能量和扫描次数增加而增加,在单脉冲能量为75μJ和扫描次数为50次时重铸层厚度为3.5μm。图 7是扫描电镜分别对空气和水介质中的激光钻孔样品进行微观形貌观察。所选用的激光参量如下:激光单脉冲能量为100μJ,空气介质中扫描次数为50次,水介质中扫描次数为600次。在空气介质中激光钻的微孔可以明显地发现, 孔内壁表面附着一层细小颗粒,即为在激光钻孔过程中熔融物重新凝结而形成的重铸层,颗粒细小,与陶瓷原表面相比重铸层有泛白倾向,但无微裂纹产生,如图 7a图 7b所示。

      Figure 7.  Scanning electron microscope picture of the drilled sample

      水介质中皮秒激光钻孔在低单脉冲能量和低扫描次数下,微孔表面无重铸层、微裂纹和发黑变质区,孔壁清洁度高、无激光烧蚀痕迹,如图 7c图 7d所示。在水介质中进行激光钻孔加工时,由于水的比热容远大于陶瓷材料,因此水的冷却作用使得陶瓷材料表面温度很难达到其熔点,陶瓷难以发生相变,因此无发黑变质和微裂纹现象。水辅助激光钻孔参量与重铸层厚度关系为:当单脉冲能量Q≤100μJ时,在扫描次数为400次~1000次范围内无重铸层产生;当单脉冲能量Q=125μJ时,在扫描次数为,在扫描次数为800次~1000次范围内出现重铸层,最大厚度为0.82μm;当单脉冲能量Q=150μJ,在扫描次数为600次~1000次范围内出现重铸层,最大厚度为1.15μm。

      水辅助激光钻孔过程中,主要存在水的冷却作用、空泡作用和运输作用。水的高比热容能够传递和吸收激光能量,能够初步降低热影响区和减少重铸层。空泡作用是皮秒激光加工时由于超高的峰值功率密度,超过物质击穿阈值产生液体空化,形成大量空泡,溃灭时产生的冲击波使熔融物能与基体脱离。同时水的存在会阻碍熔融物无法二次黏附在基体表面,同时微孔内外压强差,使得水不断流动将残渣送出陶瓷孔外,从而降低微孔锥度和去除重铸层。

    • 本文中采用水辅助脉宽10ps的激光加工系统和同心圆填充旋切钻孔技术,实现了Al2O3陶瓷无重铸层激光钻孔,并对微孔的孔径、锥度、重铸层厚度和微观形貌观察。

      (1) 水辅助激光钻孔时,孔径、锥度和重铸层厚度随激光参量的变化规律与空气中直接钻孔的变化趋势基本一致,但水辅助钻孔的孔径比空气中钻孔大,锥度和重铸层厚度比空气中钻孔小,在低单脉冲能量和低扫描速度下,能够实现完全无重铸层。

      (2) 水辅助激光钻孔时,孔径的扩大主要是由于水对激光束的散射作用引起的。本文中激光光斑直径因水中悬浮物扩大约为35μm,增大率约为75%,对应的水质参量约为0.28。

      (3) 水辅助激光钻孔在合适的参量能够实现无重铸层,小锥度高质量钻孔,具体参量是:脉冲重复频率200kHz,单脉冲能量125μJ和扫描次数600次,能够实现锥度为1.04°无重铸层激光钻孔。

      (4) 水辅助激光钻孔时,由于水的冷却作用使得陶瓷表面温度难以达到其熔点,表面基本看不到激光钻孔痕迹。皮秒激光钻孔过程中由于具有高功率密度激光,会引起液体空化现象,空泡溃灭对孔壁产生强冲击作用,使得重铸层脱离基体,并且通过水的流动作用将悬浮物运出微孔内部,从而彻底实现无重铸层钻孔。

参考文献 (15)

目录

    /

    返回文章
    返回