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

doi:10.7498/aps.70.20210165

Abstract

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.

Keywords:

Authors and contacts

1.
School of Mathematics and Physics, Lanzhou Jiaotong University, Lanzhou 730070, China
2.
Department of Physics, Lanzhou City University, Lanzhou 730070, China

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)

FullText HTML : translate this paragraph

References

[1]	Yu N F, Genevet P, Kats M A, Aieta F, Tetienne J P, Capasso F, Gaburro Z 2011 Science 334 333 Google Scholar
[2]	Ni X J, Emani N K, Kildishev A V, Boltasseva A, Shalaev V M 2012 Science 335 427 Google 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 1304 Google Scholar
[4]	Cummer S A, Christensen J, Alu A 2016 Nat. Rev. Mater. 1 16001 Google Scholar
[5]	Xie Y B, Wang W Q, Chen H Y, Konneker A, Popa B-I, Cummer S A 2014 Nat. Commun. 5 5553 Google Scholar
[6]	Xia J P, Zhang X T, Sun H X, Yuan S Q, Qian J, Ge Y 2018 Phys. Rev. Appl. 10 014016 Google Scholar
[7]	丁昌林, 董仪宝, 赵晓鹏 2018 物理学报 67 194301 Google Scholar Ding C L, Dong Y B, Zhao X P 2018 Acta Phys. Sin. 67 194301 Google Scholar
[8]	Li Y, Assouar B M 2016 Appl. Phys. Lett. 108 063502 Google Scholar
[9]	Chen C R, Du Z B, Hu G K, Yang J 2017 Appl. Phys. Lett. 110 221903 Google Scholar
[10]	Wang X L, Luo X D, Zhao H, Huang Z Y 2018 Appl. Phys. Lett. 112 021901 Google 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 103505 Google Scholar
[12]	Lu Y J, Ge Y, Yuan S Q, Sun H X, Liu X J 2020 J. Phys. D: Appl. Phys. 53 015301 Google Scholar
[13]	Pennec Y, Djafari-Rouhani B, Larabi H, Vasseur J O, Hladky-Hennion A C 2008 Phys. Rev. B 78 104105 Google Scholar
[14]	Wu T T, Huang Z G, Tsai T C, Wu T C 2008 Appl. Phys. Lett. 93 111902 Google Scholar
[15]	Wu T C, Wu T T, Hsu J C 2009 Phys. Rev. B 79 104306 Google 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 054109 Google 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 425304 Google Scholar
[18]	Jin Y, Pennec Y, Pan Y D, Djafari-Rouhani B 2016 Crystals 6 64 Google Scholar
[19]	曾伟, 王海涛, 田贵云, 胡国星, 汪文 2015 物理学报 64 134302 Google Scholar Zeng W, Wang H T, Tian G Y, Hu G X, Wang W 2015 Acta Phys. Sin. 64 134302 Google Scholar
[20]	Oudich M, Li Y 2017 J. Phys. D 50 315104 Google Scholar
[21]	Oudich M, Djafari-Rouhani B, Bonello B, Pennec Y, Hemaidia S, Sarry F, Beyssen D 2018 Phys. Rev. Appl. 9 034013 Google Scholar
[22]	周振凯, 韦利明, 丰杰 2013 物理学报 62 104601 Google Scholar Zhou Z K, Wei L M, Feng J 2013 Acta Phys. Sin. 62 104601 Google Scholar
[23]	钱莉荣, 杨保和 2013 物理学报 62 117701 Google Scholar Qian L R, Yang B H 2013 Acta Phys. Sin. 62 117701 Google Scholar
[24]	Benchabane S, Khelif A, Rauch J Y, Robert L, Laude V 2006 Phys. Rev. E 73 065601 Google Scholar
[25]	Yudistira D, Boes A, Graczykowski B, Alzina F, Yeo L Y, Sotomayor Torres C M, Mitchell A 2016 Phys. Rev. B 94 094304 Google Scholar
[26]	Ash B J, Worsfold S R, Vukusic P, Nash G R 2017 Nat. Commun. 8 174 Google Scholar
[27]	Hsu J C, Lin F S 2018 Jpn. J. Appl. Phys. 57 07LB01 Google Scholar
[28]	Coffy E, Euphrasie S, Addouche M, Vairac P, Khelif A 2017 Ultrasonics 78 51 Google Scholar
[29]	Cheng Y, Liu X J, Wu D J 2011 J. Acoust. Soc. Am. 129 1157 Google Scholar
[30]	Liu H, Huo S Y, Feng L Y, Huang H B, Chen J J 2019 Ultrasonics 94 227 Google Scholar
[31]	Li Z, Li Y M, Kumar S, Lee H P 2019 J. Appl. Phys. 126 155102 Google Scholar
[32]	郝娟, 周广刚, 马跃, 黄文奇, 张鹏, 卢贵武 2016 物理学报 65 113101 Google Scholar Hao J, Zhou G G, Ma Y, Huang W Q, Zhang P, Lu G W 2016 Acta Phys. Sin. 65 113101 Google Scholar
[33]	Yudistira D, Pennec Y, Rouhani B D, Dupont S, Laude V 2012 Appl. Phys. Lett. 100 061912 Google 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 085021 Google Scholar
[35]	Graczykowski B, Alzina F, Gomis-Bresco J, Sotomayor Torres C M 2016 J. Appl. Phys. 119 025308 Google Scholar
[36]	Bian Z G, Zhang S, Zhou X L 2019 Mech. Adv. Mater. Struct. 28 1663321 Google Scholar

Cited By

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

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

DownLoad: Full-Size Img PowerPoint

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

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

DownLoad: Full-Size Img PowerPoint

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

Figure 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.

DownLoad: Full-Size Img PowerPoint

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

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

DownLoad: Full-Size Img PowerPoint

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

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

DownLoad: Full-Size Img PowerPoint

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

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

DownLoad: Full-Size Img PowerPoint

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

Figure 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.

DownLoad: Full-Size Img PowerPoint

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

Figure 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.

DownLoad: Full-Size Img PowerPoint

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

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

DownLoad: Full-Size Img PowerPoint

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

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

材料	密度 ρ/(kg·m^–3)	杨氏模量 E/(10¹⁰ Pa)	剪切模量 μ/(10¹⁰ Pa)
环氧树脂	1180	0.435	0.159

DownLoad: CSV

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

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

材料	密度 ρ/(kg·m^–3)	弹性常数 c_ij/(10¹⁰ N·m²)						压电常数 e/(C·m^–2)				介电常数 ε/(10^–11 F·m^–1)
材料	密度 ρ/(kg·m^–3)	c₁₁	c₁₂	c₁₃	c₁₄	c₃₃	c₄₄	e₁₅	e₂₂	e₃₁	e₃₃	ε₁₁	ε₃₃
铌酸锂	4700	20.3	5.3	7.5	0.9	24.5	6.0	3.7	2.5	0.2	1.3	39.0	25.7
镍	8905	29.8	13.4				8.2

DownLoad: CSV

[1]	Yu N F, Genevet P, Kats M A, Aieta F, Tetienne J P, Capasso F, Gaburro Z 2011 Science 334 333 Google Scholar
[2]	Ni X J, Emani N K, Kildishev A V, Boltasseva A, Shalaev V M 2012 Science 335 427 Google 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 1304 Google Scholar
[4]	Cummer S A, Christensen J, Alu A 2016 Nat. Rev. Mater. 1 16001 Google Scholar
[5]	Xie Y B, Wang W Q, Chen H Y, Konneker A, Popa B-I, Cummer S A 2014 Nat. Commun. 5 5553 Google Scholar
[6]	Xia J P, Zhang X T, Sun H X, Yuan S Q, Qian J, Ge Y 2018 Phys. Rev. Appl. 10 014016 Google Scholar
[7]	丁昌林, 董仪宝, 赵晓鹏 2018 物理学报 67 194301 Google Scholar Ding C L, Dong Y B, Zhao X P 2018 Acta Phys. Sin. 67 194301 Google Scholar
[8]	Li Y, Assouar B M 2016 Appl. Phys. Lett. 108 063502 Google Scholar
[9]	Chen C R, Du Z B, Hu G K, Yang J 2017 Appl. Phys. Lett. 110 221903 Google Scholar
[10]	Wang X L, Luo X D, Zhao H, Huang Z Y 2018 Appl. Phys. Lett. 112 021901 Google 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 103505 Google Scholar
[12]	Lu Y J, Ge Y, Yuan S Q, Sun H X, Liu X J 2020 J. Phys. D: Appl. Phys. 53 015301 Google Scholar
[13]	Pennec Y, Djafari-Rouhani B, Larabi H, Vasseur J O, Hladky-Hennion A C 2008 Phys. Rev. B 78 104105 Google Scholar
[14]	Wu T T, Huang Z G, Tsai T C, Wu T C 2008 Appl. Phys. Lett. 93 111902 Google Scholar
[15]	Wu T C, Wu T T, Hsu J C 2009 Phys. Rev. B 79 104306 Google 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 054109 Google 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 425304 Google Scholar
[18]	Jin Y, Pennec Y, Pan Y D, Djafari-Rouhani B 2016 Crystals 6 64 Google Scholar
[19]	曾伟, 王海涛, 田贵云, 胡国星, 汪文 2015 物理学报 64 134302 Google Scholar Zeng W, Wang H T, Tian G Y, Hu G X, Wang W 2015 Acta Phys. Sin. 64 134302 Google Scholar
[20]	Oudich M, Li Y 2017 J. Phys. D 50 315104 Google Scholar
[21]	Oudich M, Djafari-Rouhani B, Bonello B, Pennec Y, Hemaidia S, Sarry F, Beyssen D 2018 Phys. Rev. Appl. 9 034013 Google Scholar
[22]	周振凯, 韦利明, 丰杰 2013 物理学报 62 104601 Google Scholar Zhou Z K, Wei L M, Feng J 2013 Acta Phys. Sin. 62 104601 Google Scholar
[23]	钱莉荣, 杨保和 2013 物理学报 62 117701 Google Scholar Qian L R, Yang B H 2013 Acta Phys. Sin. 62 117701 Google Scholar
[24]	Benchabane S, Khelif A, Rauch J Y, Robert L, Laude V 2006 Phys. Rev. E 73 065601 Google Scholar
[25]	Yudistira D, Boes A, Graczykowski B, Alzina F, Yeo L Y, Sotomayor Torres C M, Mitchell A 2016 Phys. Rev. B 94 094304 Google Scholar
[26]	Ash B J, Worsfold S R, Vukusic P, Nash G R 2017 Nat. Commun. 8 174 Google Scholar
[27]	Hsu J C, Lin F S 2018 Jpn. J. Appl. Phys. 57 07LB01 Google Scholar
[28]	Coffy E, Euphrasie S, Addouche M, Vairac P, Khelif A 2017 Ultrasonics 78 51 Google Scholar
[29]	Cheng Y, Liu X J, Wu D J 2011 J. Acoust. Soc. Am. 129 1157 Google Scholar
[30]	Liu H, Huo S Y, Feng L Y, Huang H B, Chen J J 2019 Ultrasonics 94 227 Google Scholar
[31]	Li Z, Li Y M, Kumar S, Lee H P 2019 J. Appl. Phys. 126 155102 Google Scholar
[32]	郝娟, 周广刚, 马跃, 黄文奇, 张鹏, 卢贵武 2016 物理学报 65 113101 Google Scholar Hao J, Zhou G G, Ma Y, Huang W Q, Zhang P, Lu G W 2016 Acta Phys. Sin. 65 113101 Google Scholar
[33]	Yudistira D, Pennec Y, Rouhani B D, Dupont S, Laude V 2012 Appl. Phys. Lett. 100 061912 Google 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 085021 Google Scholar
[35]	Graczykowski B, Alzina F, Gomis-Bresco J, Sotomayor Torres C M 2016 J. Appl. Phys. 119 025308 Google Scholar
[36]	Bian Z G, Zhang S, Zhou X L 2019 Mech. Adv. Mater. Struct. 28 1663321 Google Scholar

[1]	Han Dong-Hai, Zhang Guang-Jun, Zhao Jing-Bo, Yao Hong. Low-frequency bandgaps and sound isolation characteristics of a novel Helmholtz-type phononic crystal. Acta Physica Sinica, 2022, 71(11): 114301. doi: 10.7498/aps.71.20211932
[2]	Hao Juan, Zhou Guang-Gang, Ma Yue, Huang Wen-Qi, Zhang Peng, Lu Gui-Wu. Theoretical study on thermal and acoustic surface wave properties of Ga3PO7 crystal at high temperature. Acta Physica Sinica, 2016, 65(11): 113101. doi: 10.7498/aps.65.113101
[3]	Zeng Wei, Wang Hai-Tao, Tian Gui-Yun, Hu Guo-Xing, Wang Wen. Research on the oscillation effect of near-surface metal defect based on laser-generated acoustic surface wave. Acta Physica Sinica, 2015, 64(13): 134302. doi: 10.7498/aps.64.134302
[4]	Chen A-Li, Liang Tong-Li, Wang Yue-Sheng. Study on band gap properties of two-dimensional 8-fold quasi-periodic phononic crystals. Acta Physica Sinica, 2014, 63(3): 036101. doi: 10.7498/aps.63.036101
[5]	Qian Li-Rong, Yang Bao-He. Calculation and analysis of surface acoustic wave properties of ZnO film on diamond under different excitation conditions. Acta Physica Sinica, 2013, 62(11): 117701. doi: 10.7498/aps.62.117701
[6]	Zhang Si-Wen, Wu Jiu-Hui. Low-frequency band gaps in phononic crystals with composite locally resonant structures. Acta Physica Sinica, 2013, 62(13): 134302. doi: 10.7498/aps.62.134302
[7]	Gao Dong-Bao, Zeng Xin-Wu, Zhou Ze-Min, Tian Zhang-Fu. Experiments on defect mode of one-dimensional phononic crystal containing Helmholtz resonators. Acta Physica Sinica, 2013, 62(9): 094304. doi: 10.7498/aps.62.094304
[8]	Zhou Zhen-Kai, Wei Li-Ming, Feng Jie. Simulation of characteristics of ZnO/diamond/Si structure surface acoustic wave. Acta Physica Sinica, 2013, 62(10): 104601. doi: 10.7498/aps.62.104601
[9]	Ding Hong-Xing, Shen Zhong-Hua, Li Jia, Zhu Xue-Feng, Ni Xiao-Wu. Large partial band-gaps for Lamb waves in multiple phononic crystals thin plates. Acta Physica Sinica, 2012, 61(19): 196301. doi: 10.7498/aps.61.196301
[10]	Hu Jia-Guang, Xu Wen, Xiao Yi-Ming, Zhang Ya-Ya. The two-dimensional phononic crystal band gaps tuned by the symmetry and orientation of the additional rods in the center of unit cell. Acta Physica Sinica, 2012, 61(23): 234302. doi: 10.7498/aps.61.234302
[11]	Gao Guo-Qin, Ma Shou-Lin, Jin Feng, Kim Tong-Beum, Lu Tian-Jian. Acoustic band gaps in finite-sized two-dimensional solid/fluid phononic crystals. Acta Physica Sinica, 2010, 59(1): 393-400. doi: 10.7498/aps.59.393
[12]	Chen Sheng-Bing, Han Xiao-Yun, Yu Dian-Long, Wen Ji-Hong. Influences of different types of piezoelectric shunting circuits on band gaps of phononic beam. Acta Physica Sinica, 2010, 59(1): 387-392. doi: 10.7498/aps.59.387
[13]	Xiao Xia, You Xue-Yi, Yao Su-Ying. Dispersion feature in arbitrary direction of surface acoustic wave applied to property characterization of ultra-large-scale integrated circuit interconnect films. Acta Physica Sinica, 2007, 56(4): 2428-2433. doi: 10.7498/aps.56.2428
[14]	Zhong Hui-Lin, Wu Fu-Gen, Yao Li-Ning. Application of genetic algorithm in optimization of band gap of two-dimensional phononic crystals. Acta Physica Sinica, 2006, 55(1): 275-280. doi: 10.7498/aps.55.275
[15]	Yang Guang, P. V. Santos. Photoluminescence of GaAs(110) quantum wells modulated by surface acoustic waves. Acta Physica Sinica, 2006, 55(8): 4327-4331. doi: 10.7498/aps.55.4327
[16]	Cai Li, Han Xiao-Yun. Study of the band-structure and the uncoupled modes in two-dimensional phononic crystals with the multiple-scattering theory. Acta Physica Sinica, 2006, 55(11): 5866-5871. doi: 10.7498/aps.55.5866
[17]	Zhao Fang, Yuan Li-Bo. Characteristics of the band structure in two-dimensional phononic crystals with complex lattices. Acta Physica Sinica, 2005, 54(10): 4511-4516. doi: 10.7498/aps.54.4511
[18]	Wen Ji-Hong, Wang Gang, Liu Yao-Zong, Yu Dian-Long. Lumped-mass method on calculation of elastic band gaps of one-dimensional phononic crystals. Acta Physica Sinica, 2004, 53(10): 3384-3388. doi: 10.7498/aps.53.3384
[19]	Wang Gang, Wen Ji-Hong, Han Xiao-Yun, Zhao Hong-Gang. Finite difference time domain method for the study of band gap in two-dimensiona l phononic crystals. Acta Physica Sinica, 2003, 52(8): 1943-1947. doi: 10.7498/aps.52.1943
[20]	Qi Gong-Jin, Yang Sheng-Liang, Bai Shu-Xin, Zhao Xun. A study of the band structure in two-dimensional phononic crystals based on plane-wave algorithm. Acta Physica Sinica, 2003, 52(3): 668-671. doi: 10.7498/aps.52.668

Catalog

Vol. 70, Issue 14 - 2021-07-20

Metrics

Abstract views: 3902
PDF Downloads: 66
Cited By: 0

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

Search

Article

留言板