搜索

x

留言板

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

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

球形复合柱表面波声子晶体的带隙特性仿真

谭自豪 孙小伟 宋婷 温晓东 刘禧萱 刘子江

引用本文:
Citation:

球形复合柱表面波声子晶体的带隙特性仿真

谭自豪, 孙小伟, 宋婷, 温晓东, 刘禧萱, 刘子江

Numerical simulation study on band gap characteristics of surface phononic crystal with spherical composite column

Tan Zi-Hao, Sun Xiao-Wei, Song Ting, Wen Xiao-Dong, Liu Xi-Xuan, Liu Zi-Jiang
PDF
HTML
导出引用
  • 设计了一种由镍球与环氧树脂垫层组成的复合柱沉积在铌酸锂基体上构成的表面波声子晶体结构, 采用有限元法计算了其能带结构和位移矢量场. 结果表明: 与具有相同晶格常数的倒圆锥形表面波声子晶体结构相比, 研究结构可以在更低的频率范围打开更宽的声表面波完全带隙, 且随着复合柱半径增大, 镍球体与压电基体的硬边界之间形成限制腔模, 相邻高阶带隙间存在能量的耦合以及振动模式的继承; 此外, 温度场的引入可以实现带隙的主动调控, 带隙频率范围随着温度升高向低频移动; 通过增加复合柱体的层数, 多振子结构与行波发生多极共振耦合, 可在高阶能带间打开完全带隙. 本文的研究结果为微米级表面波声子晶体结构在100 MHz以下频率范围的带隙特性优化提供了理论支持.
    In the study of acoustic characteristics of micro-scale surface phononic crystal, the band gap characteristics below 100 MHz need to be further optimized. In this work, a piezoelectric surface phononic crystal with a composite column composed of nickel balls and epoxy backing is proposed. The finite element method is used to calculate the band gap characteristics and displacement vector field of the model. The influence of column radius on the band structure is studied, and meanwhile, the effect of the multi-layer composite column structure on the band gap is discussed via increasing the number of elements in the composite column, while the reason for the opening of the high-order band gap is analyzed in detail by combining the vibration mode. Furthermore, the temperature adjustability of the band gap is further studied. The results show that the spherical composite column deposition structure can open a wider complete band gap of surface acoustic wave in a lower frequency range than the existing inverse conical surface phononic crystal structure with the same lattice constant (Hsu J C, Lin F S 2018 Jpn. J. Appl. Phys. 57 07LB01). The restricted cavity mode is easily formed between the hard boundaries with the increase of column radius, which provides a possible way for low-order vibration modes to open high-order band gaps. There exist mode inheritance and energy coupling between adjacent modes, which leads the band gap to flatten and anti-flatten. Moreover, the real-time adjustment of band gap frequency by external temperature field can be realized via introducing the temperature-sensitive material epoxy resin into the structure. The band gap frequency range can be effectively reduced by increasing the number of composite cylinder layers, while the multi-vibrator structure can generate multipole resonance coupling with traveling wave and finally open a complete band gap between high-order frequency bands. This work provides a theoretical reference for analyzing the low-frequency band gap mechanism of micron-scale surface phononic crystal.
      通信作者: 孙小伟, sunxw_lzjtu@yeah.net
    • 基金项目: 国家自然科学基金(批准号: 51562021)、甘肃省重点人才项目(批准号: 2020RCXM100)、甘肃省自然科学基金重点项目(批准号: 20JR5RA427, 20JR5RA211)、甘肃省高等学校创新基金项目(批准号: 2020A-039)和兰州市人才创新创业项目(批准号: 2020-RC-18)资助的课题
      Corresponding author: Sun Xiao-Wei, sunxw_lzjtu@yeah.net
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 51562021), the Key Talent Foundation of Gansu Province, China (Grant No. 2020RCXM100), the Key Program of the Natural Science Foundation of Gansu Province, China (Grant Nos. 20JR5RA427, 20JR5RA211), the Innovation Fund Project of Colleges and Universities in Gansu Province, China (Grant No. 2020A-039), and the Talent Innovation and Entrepreneurship Project of Lanzhou, China (Grant No. 2020-RC-18)
    [1]

    Yu N F, Genevet P, Kats M A, Aieta F, Tetienne J P, Capasso F, Gaburro Z 2011 Science 334 333Google Scholar

    [2]

    Ni X J, Emani N K, Kildishev A V, Boltasseva A, Shalaev V M 2012 Science 335 427Google Scholar

    [3]

    Grady N K, Heyes J E, Chowdhury D R, Zeng Y, Reiten M T, Azad A K, Taylor A J, Dalvit D A R, Chen H T 2013 Science 340 1304Google Scholar

    [4]

    Cummer S A, Christensen J, Alu A 2016 Nat. Rev. Mater. 1 16001Google Scholar

    [5]

    Xie Y B, Wang W Q, Chen H Y, Konneker A, Popa B-I, Cummer S A 2014 Nat. Commun. 5 5553Google Scholar

    [6]

    Xia J P, Zhang X T, Sun H X, Yuan S Q, Qian J, Ge Y 2018 Phys. Rev. Appl. 10 014016Google Scholar

    [7]

    丁昌林, 董仪宝, 赵晓鹏 2018 物理学报 67 194301Google Scholar

    Ding C L, Dong Y B, Zhao X P 2018 Acta Phys. Sin. 67 194301Google Scholar

    [8]

    Li Y, Assouar B M 2016 Appl. Phys. Lett. 108 063502Google Scholar

    [9]

    Chen C R, Du Z B, Hu G K, Yang J 2017 Appl. Phys. Lett. 110 221903Google Scholar

    [10]

    Wang X L, Luo X D, Zhao H, Huang Z Y 2018 Appl. Phys. Lett. 112 021901Google Scholar

    [11]

    Wu X X, Au-Yeung K Y, Li X, Roberts R C, Tian J X, Hu C D, Huang Y Z, Wang S X, Yang Z Y, Wen W J 2018 Appl. Phys. Lett. 112 103505Google Scholar

    [12]

    Lu Y J, Ge Y, Yuan S Q, Sun H X, Liu X J 2020 J. Phys. D: Appl. Phys. 53 015301Google Scholar

    [13]

    Pennec Y, Djafari-Rouhani B, Larabi H, Vasseur J O, Hladky-Hennion A C 2008 Phys. Rev. B 78 104105Google Scholar

    [14]

    Wu T T, Huang Z G, Tsai T C, Wu T C 2008 Appl. Phys. Lett. 93 111902Google Scholar

    [15]

    Wu T C, Wu T T, Hsu J C 2009 Phys. Rev. B 79 104306Google Scholar

    [16]

    Jin Y, Fernez N, Pennec Y, Bonello B, Moiseyenko R P, Hémon S, Pan Y D, Djafari-Rouhani B 2016 Phys. Rev. B 93 054109Google Scholar

    [17]

    Jin Y, EI Boudouti E H, Pennec Y, Bonello B, Moiseyenko R P, Hémon S, Pan Y D, Djafari-Rouhani B 2017 J. Phys. D 50 425304Google Scholar

    [18]

    Jin Y, Pennec Y, Pan Y D, Djafari-Rouhani B 2016 Crystals 6 64Google Scholar

    [19]

    曾伟, 王海涛, 田贵云, 胡国星, 汪文 2015 物理学报 64 134302Google Scholar

    Zeng W, Wang H T, Tian G Y, Hu G X, Wang W 2015 Acta Phys. Sin. 64 134302Google Scholar

    [20]

    Oudich M, Li Y 2017 J. Phys. D 50 315104Google Scholar

    [21]

    Oudich M, Djafari-Rouhani B, Bonello B, Pennec Y, Hemaidia S, Sarry F, Beyssen D 2018 Phys. Rev. Appl. 9 034013Google Scholar

    [22]

    周振凯, 韦利明, 丰杰 2013 物理学报 62 104601Google Scholar

    Zhou Z K, Wei L M, Feng J 2013 Acta Phys. Sin. 62 104601Google Scholar

    [23]

    钱莉荣, 杨保和 2013 物理学报 62 117701Google Scholar

    Qian L R, Yang B H 2013 Acta Phys. Sin. 62 117701Google Scholar

    [24]

    Benchabane S, Khelif A, Rauch J Y, Robert L, Laude V 2006 Phys. Rev. E 73 065601Google Scholar

    [25]

    Yudistira D, Boes A, Graczykowski B, Alzina F, Yeo L Y, Sotomayor Torres C M, Mitchell A 2016 Phys. Rev. B 94 094304Google Scholar

    [26]

    Ash B J, Worsfold S R, Vukusic P, Nash G R 2017 Nat. Commun. 8 174Google Scholar

    [27]

    Hsu J C, Lin F S 2018 Jpn. J. Appl. Phys. 57 07LB01Google Scholar

    [28]

    Coffy E, Euphrasie S, Addouche M, Vairac P, Khelif A 2017 Ultrasonics 78 51Google Scholar

    [29]

    Cheng Y, Liu X J, Wu D J 2011 J. Acoust. Soc. Am. 129 1157Google Scholar

    [30]

    Liu H, Huo S Y, Feng L Y, Huang H B, Chen J J 2019 Ultrasonics 94 227Google Scholar

    [31]

    Li Z, Li Y M, Kumar S, Lee H P 2019 J. Appl. Phys. 126 155102Google Scholar

    [32]

    郝娟, 周广刚, 马跃, 黄文奇, 张鹏, 卢贵武 2016 物理学报 65 113101Google Scholar

    Hao J, Zhou G G, Ma Y, Huang W Q, Zhang P, Lu G W 2016 Acta Phys. Sin. 65 113101Google Scholar

    [33]

    Yudistira D, Pennec Y, Rouhani B D, Dupont S, Laude V 2012 Appl. Phys. Lett. 100 061912Google Scholar

    [34]

    Zhang D B, Zhao J F, Bonello B, Li L B, Wei J X, Pan Y D, Zhong Z 2016 AIP Adv. 6 085021Google Scholar

    [35]

    Graczykowski B, Alzina F, Gomis-Bresco J, Sotomayor Torres C M 2016 J. Appl. Phys. 119 025308Google Scholar

    [36]

    Bian Z G, Zhang S, Zhou X L 2019 Mech. Adv. Mater. Struct. 28 1663321Google Scholar

  • 图 1  球形复合柱表面波声子晶体单胞结构及第一布里渊区示意图

    Fig. 1.  Schematic diagram of the unit cell and the first Brillouin zone of the surface phononic crystal with spherical composite column.

    图 2  用于计算所设计表面波声子晶体结构插入损失的半无限结构

    Fig. 2.  The semi-finite structure for calculating the insertion loss of designed surface phononic crystal.

    图 3  球形复合柱半径不同取值时所设计模型的能带结构以及对比模型的能带结构 (a)左图为柱体半径r = 0.3a时的能带结构图, 右图为对应的插入损失谱; (b)参考文献[26]中对比模型半径r = 0.32a时的能带结构图; (c)柱体半径r = 0.4a时的能带结构; (d)柱体半径r = 0.5a时的能带结构

    Fig. 3.  Band structures of the designed model with different radius of spherical composite column and the comparison model: (a) Designed model with cylinder radius r = 0.3a in the left and its insertion loss spectrum in the right; (b) comparison model with cylinder radius r = 0.32a from Ref. [26]; (c) designed model with cylinder radius r = 0.4a and (d) r = 0.5a, respectively.

    图 4  (a)—(j)分别为图3(a)AF点附近的振动模态. 红色代表振动部分, 振动位移的大小如图例所示

    Fig. 4.  (a)–(j) are the vibration modes at the marked points AF in Fig. 3, respectively. The vibration part corresponds to the red and the magnitude of vibration displacement is shown in the legend.

    图 5  r = 0.5a时, 所设计表面波声子晶体的第7—10条能带的结构图

    Fig. 5.  The band structure of the seventh to tenth bands of designed surface phononic crystal when r = 0.5a.

    图 6  (a)—(r)分别为图5中AR标记点处的振动模态.

    Fig. 6.  (a)–(r) are the vibration modes at the marked points AR in Fig. 5, respectively.

    图 7  (a)参考文献[36]中随温度的变化对环氧树脂弹性模量的影响; (b) r = 0.4a时, 前6条能带随温度的变化

    Fig. 7.  (a) Effect of temperature change on elastic modulus of epoxy resin in the Ref. [36]; (b) change of the first six bands with temperature when r = 0.4a.

    图 8  r = 0.4a, h = 0.2a时, 多层球形复合柱表面波声子晶体的模型结构示意图及其能带结构 (a)两层球形复合柱模型结构示意图; (b)两层球形复合柱模型的能带结构图; (c)球形复合柱层数分别为三层、四层和五层时的能带结构及其模型示意图

    Fig. 8.  Model structures and band structures of the surface phononic crystals with multi-layer spherical composite columns when r = 0.4a and h = 0.2a: (a) Schematic of the two-layer spherical composite column model structure; (b) band structures of the designed model with two-layer spherical composite column; (c) band structures of the designed model with three, four and five layers, respectively.

    图 9  两层球形复合柱模型的能带结构中k = X点处的振动模态

    Fig. 9.  Vibration modes at the points k = X of the band structure of the designed model with two-layer.

    表 1  所设计模型的弹性垫层材料参数

    Table 1.  Material parameters of elastic cushion of the designed model.

    材料密度
    ρ/(kg·m–3)
    杨氏模量
    E/(1010 Pa)
    剪切模量
    μ/(1010 Pa)
    环氧树脂11800.4350.159
    下载: 导出CSV

    表 2  所设计模型的基体与散射体材料参数

    Table 2.  Material parameters of matrix and scatterer of the designed model.

    材料密度
    ρ/(kg·m–3)
    弹性常数
    cij/(1010 N·m2)
    压电常数
    e/(C·m–2)
    介电常数
    ε/(10–11 F·m–1)
    c11c12c13c14c33c44e15e22e31e33ε11ε33
    铌酸锂470020.35.37.50.924.56.0 3.72.50.21.3 39.025.7
    890529.813.48.2
    下载: 导出CSV
  • [1]

    Yu N F, Genevet P, Kats M A, Aieta F, Tetienne J P, Capasso F, Gaburro Z 2011 Science 334 333Google Scholar

    [2]

    Ni X J, Emani N K, Kildishev A V, Boltasseva A, Shalaev V M 2012 Science 335 427Google Scholar

    [3]

    Grady N K, Heyes J E, Chowdhury D R, Zeng Y, Reiten M T, Azad A K, Taylor A J, Dalvit D A R, Chen H T 2013 Science 340 1304Google Scholar

    [4]

    Cummer S A, Christensen J, Alu A 2016 Nat. Rev. Mater. 1 16001Google Scholar

    [5]

    Xie Y B, Wang W Q, Chen H Y, Konneker A, Popa B-I, Cummer S A 2014 Nat. Commun. 5 5553Google Scholar

    [6]

    Xia J P, Zhang X T, Sun H X, Yuan S Q, Qian J, Ge Y 2018 Phys. Rev. Appl. 10 014016Google Scholar

    [7]

    丁昌林, 董仪宝, 赵晓鹏 2018 物理学报 67 194301Google Scholar

    Ding C L, Dong Y B, Zhao X P 2018 Acta Phys. Sin. 67 194301Google Scholar

    [8]

    Li Y, Assouar B M 2016 Appl. Phys. Lett. 108 063502Google Scholar

    [9]

    Chen C R, Du Z B, Hu G K, Yang J 2017 Appl. Phys. Lett. 110 221903Google Scholar

    [10]

    Wang X L, Luo X D, Zhao H, Huang Z Y 2018 Appl. Phys. Lett. 112 021901Google Scholar

    [11]

    Wu X X, Au-Yeung K Y, Li X, Roberts R C, Tian J X, Hu C D, Huang Y Z, Wang S X, Yang Z Y, Wen W J 2018 Appl. Phys. Lett. 112 103505Google Scholar

    [12]

    Lu Y J, Ge Y, Yuan S Q, Sun H X, Liu X J 2020 J. Phys. D: Appl. Phys. 53 015301Google Scholar

    [13]

    Pennec Y, Djafari-Rouhani B, Larabi H, Vasseur J O, Hladky-Hennion A C 2008 Phys. Rev. B 78 104105Google Scholar

    [14]

    Wu T T, Huang Z G, Tsai T C, Wu T C 2008 Appl. Phys. Lett. 93 111902Google Scholar

    [15]

    Wu T C, Wu T T, Hsu J C 2009 Phys. Rev. B 79 104306Google Scholar

    [16]

    Jin Y, Fernez N, Pennec Y, Bonello B, Moiseyenko R P, Hémon S, Pan Y D, Djafari-Rouhani B 2016 Phys. Rev. B 93 054109Google Scholar

    [17]

    Jin Y, EI Boudouti E H, Pennec Y, Bonello B, Moiseyenko R P, Hémon S, Pan Y D, Djafari-Rouhani B 2017 J. Phys. D 50 425304Google Scholar

    [18]

    Jin Y, Pennec Y, Pan Y D, Djafari-Rouhani B 2016 Crystals 6 64Google Scholar

    [19]

    曾伟, 王海涛, 田贵云, 胡国星, 汪文 2015 物理学报 64 134302Google Scholar

    Zeng W, Wang H T, Tian G Y, Hu G X, Wang W 2015 Acta Phys. Sin. 64 134302Google Scholar

    [20]

    Oudich M, Li Y 2017 J. Phys. D 50 315104Google Scholar

    [21]

    Oudich M, Djafari-Rouhani B, Bonello B, Pennec Y, Hemaidia S, Sarry F, Beyssen D 2018 Phys. Rev. Appl. 9 034013Google Scholar

    [22]

    周振凯, 韦利明, 丰杰 2013 物理学报 62 104601Google Scholar

    Zhou Z K, Wei L M, Feng J 2013 Acta Phys. Sin. 62 104601Google Scholar

    [23]

    钱莉荣, 杨保和 2013 物理学报 62 117701Google Scholar

    Qian L R, Yang B H 2013 Acta Phys. Sin. 62 117701Google Scholar

    [24]

    Benchabane S, Khelif A, Rauch J Y, Robert L, Laude V 2006 Phys. Rev. E 73 065601Google Scholar

    [25]

    Yudistira D, Boes A, Graczykowski B, Alzina F, Yeo L Y, Sotomayor Torres C M, Mitchell A 2016 Phys. Rev. B 94 094304Google Scholar

    [26]

    Ash B J, Worsfold S R, Vukusic P, Nash G R 2017 Nat. Commun. 8 174Google Scholar

    [27]

    Hsu J C, Lin F S 2018 Jpn. J. Appl. Phys. 57 07LB01Google Scholar

    [28]

    Coffy E, Euphrasie S, Addouche M, Vairac P, Khelif A 2017 Ultrasonics 78 51Google Scholar

    [29]

    Cheng Y, Liu X J, Wu D J 2011 J. Acoust. Soc. Am. 129 1157Google Scholar

    [30]

    Liu H, Huo S Y, Feng L Y, Huang H B, Chen J J 2019 Ultrasonics 94 227Google Scholar

    [31]

    Li Z, Li Y M, Kumar S, Lee H P 2019 J. Appl. Phys. 126 155102Google Scholar

    [32]

    郝娟, 周广刚, 马跃, 黄文奇, 张鹏, 卢贵武 2016 物理学报 65 113101Google Scholar

    Hao J, Zhou G G, Ma Y, Huang W Q, Zhang P, Lu G W 2016 Acta Phys. Sin. 65 113101Google Scholar

    [33]

    Yudistira D, Pennec Y, Rouhani B D, Dupont S, Laude V 2012 Appl. Phys. Lett. 100 061912Google Scholar

    [34]

    Zhang D B, Zhao J F, Bonello B, Li L B, Wei J X, Pan Y D, Zhong Z 2016 AIP Adv. 6 085021Google Scholar

    [35]

    Graczykowski B, Alzina F, Gomis-Bresco J, Sotomayor Torres C M 2016 J. Appl. Phys. 119 025308Google Scholar

    [36]

    Bian Z G, Zhang S, Zhou X L 2019 Mech. Adv. Mater. Struct. 28 1663321Google Scholar

  • [1] 黄铭贤, 胡文彬, 白飞明. 声表面波-自旋波耦合及磁声非互易性器件. 物理学报, 2024, 73(15): 158501. doi: 10.7498/aps.73.20240462
    [2] 韩东海, 张广军, 赵静波, 姚宏. 新型Helmholtz型声子晶体的低频带隙及隔声特性. 物理学报, 2022, 71(11): 114301. doi: 10.7498/aps.71.20211932
    [3] 郝娟, 周广刚, 马跃, 黄文奇, 张鹏, 卢贵武. 高温条件下Ga3PO7晶体热学及声表面波性质的理论研究. 物理学报, 2016, 65(11): 113101. doi: 10.7498/aps.65.113101
    [4] 曾伟, 王海涛, 田贵云, 胡国星, 汪文. 研究激光激发的声表面波与材料近表面缺陷的振荡效应. 物理学报, 2015, 64(13): 134302. doi: 10.7498/aps.64.134302
    [5] 陈阿丽, 梁同利, 汪越胜. 二维8重固-流型准周期声子晶体带隙特性研究. 物理学报, 2014, 63(3): 036101. doi: 10.7498/aps.63.036101
    [6] 张思文, 吴九汇. 局域共振复合单元声子晶体结构的低频带隙特性研究. 物理学报, 2013, 62(13): 134302. doi: 10.7498/aps.62.134302
    [7] 高东宝, 曾新吾, 周泽民, 田章福. 一维亥姆霍兹共振腔声子晶体中缺陷模式的实验研究. 物理学报, 2013, 62(9): 094304. doi: 10.7498/aps.62.094304
    [8] 周振凯, 韦利明, 丰杰. ZnO/Diamond/Si结构中声表面波传播特性分析. 物理学报, 2013, 62(10): 104601. doi: 10.7498/aps.62.104601
    [9] 丁红星, 沈中华, 李加, 祝雪丰, 倪晓武. 复合兰姆波声子晶体中超宽部分禁带. 物理学报, 2012, 61(19): 196301. doi: 10.7498/aps.61.196301
    [10] 胡家光, 徐文, 肖宜明, 张丫丫. 晶格中心插入体的对称性及取向对二维声子晶体带隙的影响. 物理学报, 2012, 61(23): 234302. doi: 10.7498/aps.61.234302
    [11] 高国钦, 马守林, 金峰, 金东范, 卢天健. 声波在二维固/流声子晶体中的禁带特性研究. 物理学报, 2010, 59(1): 393-400. doi: 10.7498/aps.59.393
    [12] 陈圣兵, 韩小云, 郁殿龙, 温激鸿. 不同压电分流电路对声子晶体梁带隙的影响. 物理学报, 2010, 59(1): 387-392. doi: 10.7498/aps.59.387
    [13] 肖 夏, 尤学一, 姚素英. 表征超大规模集成电路互连纳米薄膜硬度特性的声表面波的频散特性. 物理学报, 2007, 56(4): 2428-2433. doi: 10.7498/aps.56.2428
    [14] 钟会林, 吴福根, 姚立宁. 遗传算法在二维声子晶体带隙优化中的应用. 物理学报, 2006, 55(1): 275-280. doi: 10.7498/aps.55.275
    [15] 杨 光, P. V. Santos. 声表面波对GaAs(110)量子阱发光特性的调制. 物理学报, 2006, 55(8): 4327-4331. doi: 10.7498/aps.55.4327
    [16] 蔡 力, 韩小云. 二维声子晶体带结构的多散射分析及解耦模式. 物理学报, 2006, 55(11): 5866-5871. doi: 10.7498/aps.55.5866
    [17] 赵 芳, 苑立波. 二维复式格子声子晶体带隙结构特性. 物理学报, 2005, 54(10): 4511-4516. doi: 10.7498/aps.54.4511
    [18] 温激鸿, 王 刚, 刘耀宗, 郁殿龙. 基于集中质量法的一维声子晶体弹性波带隙计算. 物理学报, 2004, 53(10): 3384-3388. doi: 10.7498/aps.53.3384
    [19] 王 刚, 温激鸿, 韩小云, 赵宏刚. 二维声子晶体带隙计算中的时域有限差分方法. 物理学报, 2003, 52(8): 1943-1947. doi: 10.7498/aps.52.1943
    [20] 齐共金, 杨盛良, 白书欣, 赵 恂. 基于平面波算法的二维声子晶体带结构的研究. 物理学报, 2003, 52(3): 668-671. doi: 10.7498/aps.52.668
计量
  • 文章访问数:  5205
  • PDF下载量:  82
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-24
  • 修回日期:  2021-03-05
  • 上网日期:  2021-07-15
  • 刊出日期:  2021-07-20

/

返回文章
返回