+Advanced Search
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

Evaluation of Comprehensive Performance of High Damping and High Strength As-Cast SiCp/Mg94Zn5Y1 Composites  PDF

  • Wan Diqing
  • Dong Shaoyun
  • Wang Houbin
  • Hu Jiajun
  • Xue Yandan
  • Han Guoliang
  • Kang Jie
  • Zeng Guanmei
  • Wang Yu
  • Tang Hao
  • Yang Fan
School of Materials Science and Engineering, East China Jiaotong University, Nanchang 330013, China

Updated:2022-12-07

  • Full Text
  • Figs & Tabs
  • References
  • Authors
  • About
CN CITE
OUTLINE

Abstract

SiCP/Mg94Zn5Y1 composites with 0.5wt%~2.0wt% SiCp were prepared by the casting method. The mechanical properties and damping capacities of the composites were investigated. The microstructure and phase components of the composites were analyzed via the scanning electron microscopy and X-ray diffraction. Results show that after the addition of SiCp into the matrix, the SiCp is evenly distributed in the matrix, which refines the microstructure of the composite. The SiCp/Mg94Zn5Y1 composites contain the α-Mg, I-phase (quasicrystal phase), and SiCp phase. The damping capacities and mechanical properties of the SiCp/Mg94Zn5Y1 composites were evaluated by the dynamic mechanical analyzer and an AG-X testing machine, respectively. The mechanical properties of the composites are better than those of the original Mg94Zn5Y1 alloy. The 1.0wt% SiCp/Mg94Zn5Y1 composite exhibits the compressive strength of 350 MPa. The damping properties of all the composites are much higher than those of the parent alloy. The optimal damping capacity is achieved when the composite contains 0.5wt% SiCp. Moreover, according to the efficiency coefficient method, the 1.0wt% SiCp/Mg94Zn5Y1 composite has the optimal comprehensive properties.

Science Press

Magnesium and its alloys have received considerable attention as the promising novel engineering materials due to their high specific strength, low density, and high damping capacities[

1-4]. However, the engineering applications of high damping magnesium alloys are restricted due to the discord between their damping and mechanical properties[5]. Adding suitable reinforcement phases into the Mg-based materials can improve the strength and maintain the high damping capacity. Because Mg-based composites also have the advantage of low density, the shorcoming of low stiffness can be overcome. Thus, the Mg-based composites can be used as lightweight nonferrous materials[6-8]. Zhang et al[9] prepared MgO/ZK60 nanocomposites by traditional powder metallurgy, and the composite strength was improved. Yu et al[10] prepared the AZ91D composites reinforced by Ti2AlC phase with good damping performance. Normally, the mechanical or damping properties of Mg matrix composites can be solely improved. Therefore, the novel Mg matrix composites with both good mechanical and damping properties are required.

The performance of Mg-based composites mainly depends on the matrix, reinforcement, and the bonding strength between the reinforcement and the matrix[

11]. The reinforcement phase should have solid load-bearing capacity, good wettability, thermodynamic stability, high hardness, and high elastic modulus[12]. SiC particles (SiCp) have the advantages of high strength, large modulus, good thermal conductivity, and good wettability with Mg alloy melt[13]. Prasad et al[14] reported that SiCp can refine the grains and provide good damping properties for the composites. According to Ref.[15], the damping capacity of SiCp/AZ61 composites is better than that of AZ61 alloy due to the SiCp addition.

The Mg-Zn-Y alloys with quasi-crystalline particles are potential high-strength magnesium alloys[

16-20], because their quasicrystal phase can improve the mechanical properties of the alloys. However, the damping capacity of the alloy is correspondingly reduced due to the formation of quasicrystal[21]. The synergistic effect of SiCp and quasicrystal phase (I-phase) in the composites can achieve the favorable coordination between the mechanical and damping properties. However, the properties of coexisting SiCp and quasicrystal phases in the as-cast Mg-based composites are rarely reported. In this research, the SiCp ceramic reinforcement phase was added into the Mg-Zn-Y matrix. The mechanical and damping properties of Mg-Zn-Y composites with SiCp and quasicrystal reinforcement phases were studied, and the strengthening and damping mechanisms were also investigated. Additionally, the mathematical models were used to evaluate the comprehensive properties of the materials.

1 Experiment

The raw Mg ingot (99.9% purity), Zn ingot (99.9% purity), and Mg-Y master alloy (25wt% Y) were used to prepare the Mg94Zn5Y1 matrix. The SiCp was selected as the reinforcement phase. Firstly, the pure Mg was placed into a crucible preheated to 300 °C and smelted in a resistance furnace at 720 °C. Once the Mg ingot was completely molten, the Mg-Y master alloy and pure Zn ingot were sequentially added to obtain Mg94Zn5Y1 alloy. Finally, the 0wt%, 0.5wt%, 1.0wt%, and 2.0wt% SiCp was separately added into the mixed melt. After stirring at 700 r/min for 90 s in Ar atmosphere, the melt was poured into a steel mold. The composites of different SiCp contents were produced.

The specimen microstructure was observed by the scanning electron microscope (SEM, Hitachi SU-8010) equipped with energy dispersive spectroscope (EDS). The phase constituent analysis was performed via X-ray diffraction (XRD, Shimadzu XRD-6100) with monochromatic Cu Kα radiation. In addition, the specimens were compressed at 0.5 mm·min-1 by the AG-X testing machine (Shimadzu corporation). The damping capacities of specimens were measured by the dynamic mechanical analyzer (DMA Q800, TA Instruments). The strain ranged from 1×10-4 to 4×10-3 under the conditions of room temperature and frequency of 1 Hz.

2 Results and Discussion

2.1 Microstructure

The microstructures of the as-cast Mg94Zn5Y1 alloy and different SiCp/Mg94Zn5Y1 composites are shown in Fig.1. It can be seen that the Mg94Zn5Y1 alloy is composed of the matrix (α-Mg) and a lamellar eutectic phase (α-Mg+I-phase). During the solidification of Mg94Zn5Y1 alloy, the primary α-Mg phase is firstly crystallized, and the eutectic structure (α-Mg+I-phase) is formed via the pseudo-eutectic reaction[

22,23]. Yi et al[23] indicated that the two-phase (α-Mg+I-phase) field exists in the Mg-Zn-Y ternary alloy system when the content ratio of Zn to Y is close to 5 or 6. Kim et al[24] revealed that the α-Mg is firstly crystallized from the melt to form a dendritic morphology, and the remaining liquid in the interdendritic region is solidified into the eutectic structure (α-Mg+I-phase) when the composition changes. Fig.1b shows the magnificent microstructure of I-phase, namely quasicrystal phase, in the Mg94Zn5Y1 composite. The eutectic phase can hinder the diffusion of Mg, Zn, and Y elements more effectively than the grainy I-phase does. Therefore, the eutectic phase can refine the α-Mg dendrites more effectively[25]. The SiCp addition affects the temperature distribution and solute distribution of the solid-liquid, thereby influencing the morphology and growth kinetics of the lamellar eutectic phase. After the SiCp is added into the matrix (Fig.1c~1e), the composites are composed of α-Mg dendrites, eutectic phase, and scattered particles. Fig.2b shows EDS Si element distribution in the 0.5wt%SiCp/Mg94Zn5Y1 composite. It can be seen that the Si element is relatively evenly distributed in the matrix. Fig.3 shows XRD patterns of different SiCp/Mg94Zn5Y1 composites. The diffraction peaks are mainly related to the α-Mg, I-phase, and SiCp phase. According to Fig.2 and Fig.3, the SiCp stably exists in the matrix and has a phase strengthening effect.

Fig.1  Microstructures of SiCp/Mg94Zn5Y1 composites with different SiCp contents: (a, b) 0wt%, (c) 0.5wt%, (d) 1.0wt%, and (e) 2.0wt%

Fig.2  SEM microstructure (a) and corresponding EDS Si element distribution (b) of 0.5wt% SiCp/Mg94Zn5Y1 composite

Fig.3  XRD patterns of different SiCp/Mg94Zn5Y1 composites

2.2 Compressive properties

The elastic modulus of SiCp and Mg-based matrix is different. During the compressive deformation process, the SiCp and Mg-based matrix are coordinatively deformed, i.e., the load is transferred from the soft Mg-based matrix to the hard SiCp reinforcement phase, which improves the mechanical properties of the composites. Particularly, the I-phase improves the mechanical properties of the composites. Fig.4 shows the compressive stress-strain curves of different SiCp/Mg94Zn5Y1 composites. After the SiCp is added into the matrix, the compressive properties of the SiCp/Mg94Zn5Y1 composites are significantly better than those of the Mg94Zn5Y1 alloy. According to Table 1, the optimal mechanical properties (compressive strength of 350 MPa) are achieved when the SiCp content is 1.0wt%, which is 12.5% higher than the compressive strength of the Mg94Zn5Y1 base alloy. With increasing the SiCp content, the compressive properties of composites are increased firstly and then decreased. When the SiCp content is 0.5wt% and 1.0wt%, the SiCp is evenly distributed in the matrix without apparent agglomeration. With further increasing the SiCp content, the secondary phase strengthening, dislocation strengthening, and other strengthening mechanisms are more obvious. However, when 2wt% SiCp is added into the matrix, a small amount of SiCp is agglomerated in the composite and the fine grain strengthening effect is weakened. Thus, the mechanical properties of the material are reduced.

Fig.4  Compressive stress-strain curves of different SiCp/Mg94Zn5Y1 composites

Table 1  Mechanical properties of different SiCp/Mg94Zn5Y1 composites
SiCp content/wt%00.51.02.0
Compressive strength/MPa 311 344 350 342
Compression ratio/% 0.110 0.155 0.166 0.152

The strengthening mechanism of SiCp/Mg94Zn5Y1 compo-sites is primarily related to the secondary phase strengthening, fine grain strengthening, and dislocation strengthening mechanisms. In the SiCp-reinforced Mg-Zn-Y composites, the thermal expansion coefficient of SiCp is much smaller than that of the Mg-Zn-Y alloy[

26,27], which causes the residual stress in the material. As a result, the plastic deformation occurs in the Mg-based composites, forming the high-density dislocations, and thereby causing the new dislocation strengthening. Furthermore, the α-Mg interdendritic micro-structure in the lamellar eutectic phases (α-Mg+I-phase) causes the eutectic phase strengthening effect[28]. The addition of SiCp has a significant effect on the grain refinement of the matrix alloy. Because the grains are refined, the applied force is increased to activate the sources of dislocations in adjacent grains. A greater external force is required when the fine grains produce the plastic deformation[29].

2.3 Damping capacities

The damping capacities of the SiCp/Mg94Zn5Y1 composites are shown in Fig.5. It can be seen that all the damping properties of the composites with SiCp addition are higher than those without the SiCp addition. The curves can be divided into the strain-amplitude-independent part and the strain-amplitude-dependent part. In the strain-amplitude-independent part, the damping value of the composite is close to that of the matrix alloy. At high strain amplitudes, the damping performance is firstly increased and then decreased with increasing the SiCp content. Furthermore, after SiCp is added to the matrix, the critical strain amplitude of all the composites materials is higher than that of the matrix alloy, which indicates that the SiCp addition improves the damping performance of the matrix alloy. In the strain-amplitude-dependent part, the optimal damping capacity is achieved when the SiCp content is 0.5wt%. When the SiCp content is 1.0wt%, different damping mechanisms are functioning well, and the damping value of the 1.0wt% SiCp/Mg94Zn5Y1 composite is close to that of the 0.5wt% SiCp/Mg94Zn5Y1 composite. The damping value Qd-1 of the SiCp/Mg94Zn5Y1 composites at the strain of 4×10-4 is shown in Fig.5b. It is known that when the damping capacity of the material is over 0.01, the material is a high damping material[

21]. It can be seen that the damping value of all the composites with SiCp addition is greater than 0.01, indicating that the SiCp/Mg94Zn5Y1 composites are the high-damping materials. In addition, the 0.5wt% SiCp/Mg94Zn5Y1 composite shows the highest damping value of 0.011 80.

Fig.5  Damping value Q-1-strain amplitude curves (a) and damping value Q-1 at strain of 4×10-4 (b) of different SiCp/Mg94Zn5Y1 composites

Generally, the damping mechanism of Mg-based alloys is related to the dislocations and dominated by the Granato-Lücke (G-L) theory[

30,31], which proves that the dislocations are pinned by strong pinners (dislocation nodes, secondary phases, and grain boundaries) and weak pinners (impurity atoms, vacancies, and disorders). At low strain amplitudes, the dislocations are reciprocating between the weak pinners, and the damping is independent of or weakly dependent on the strain amplitude. With increasing the strain amplitude, the dislocations break away from the weak pinners and slip between the strong pinners. The damping capacity of the composites is increased with increasing the strain amplitude. The dislocation damping (Qd-1), strain-independent partial damping capacity (Q0-1), and strain-related partial damping capacity (Qh-1) can be expressed by Eq.(1~3), respectively, as follows:

Qd-1=Q0-1+Qh-1 (1)
Q0-1=ρBLC4ω36Gb2 (2)
Qh-1=C1εexp(-C2ε) (3)
C1=ρFBLN36bELC2 C2=FBbELC (4)

where ρ is the density of a removable dislocation; B is the damping constant; LN and LC are the average distances between successive strong pinners and successive weak pinners, respectively; ω is the angular frequency; FB is binding force between the weak pinners and the dislocations; b is the Burgers vector; G is the shear modulus; E is the elastic modulus; ε is the strain.

The G-L model considers that the dislocation in the crystal is strongly pinned by the immobile defects, such as grain boundaries, precipitation phases, or dislocation network nodes, and the intermediate portion is pinned by impurity atoms. The dislocation segment length LN and the average length between the weak pinners LC are governed by the defects and impuri-ties. The schematic diagram of the changes in the dislocation length with increasing the external stress[

30,31] is shown in Fig.6. When the applied stress is zero, the line (LN) is pinned by the strong pinners (line A in Fig.6). With increasing the stress, the loops (LC) gradually leave the weak pinners and continue to increase until the breakaway stress is reached (evolution from line B to line D in Fig.6). With further increa-sing the stress, the dislocation is also increased (evolution from line E to line F in Fig.6). Then, the new closed disloca-tion loops form (line G in Fig.6). When the stress is removed, the line (LN) returns to its original position and is repined.

Fig.6  Schematic diagram of G-L dislocation pinning model

The Mg94Zn5Y1 composite has high residual stress around the SiCp during solidification, which results in a high density of dislocations. It is known from the G-L theory that Q0-1 is proportional to ρ at low strain amplitudes. Hence, in the low-strain-amplitude stage, the damping capabilities of the SiCp/Mg94Zn5Y1 composites are better than those of the Mg94Zn5Y1 alloy.

According to the G-L theory, the relationship of Qh-1 with C1 and C2 can be expressed by Eq.(5), as follows:

ln(Qh-1ε)=lnC1-C2/ε (5)

Fig.7 shows the G-L plots of the composites in the strain amplitude correlation stage, and Table 2 shows the values of C1 and C2 of the composites. It can be seen that the C1 value of the composites with SiCp addition is greater than that without SiCp addition, but the C2 value of the composites with SiCp addition is smaller than that without the SiCp addition. According to Eq.(3), the Qh-1 of the composites with SiCp addition is better than that without the SiCp addition. According to Eq.(4), the LC value of the composites with the SiCp addition is smaller than that without the SiCp addition. This is mainly because the added SiCp becomes the heterogeneous nucleation core to promote the formation of quasicrystal phases, which results in a decrease in the solute atoms in the matrix, i.e., the number of weak pinners is decreased. With increasing the SiCp content, the C1 value is gradually decreased. It can be seen from Eq.(4) that C1 is proportional to ρLN3. With increasing the content of the reinforcement phase and the secondary phase, the pinning effect of the dislocation becomes more obvious and the dislocation mobility is reduced.

Fig.7  G-L plots of the different SiCp/Mg94Zn5Y1 composites

Table 2  Values of C1 and C2 of different SiCp/Mg94Zn5Y1 composites
SiCp content/wt%C1C2/×102
0 0.0295 0.9125
0.5 0.0457 0.7854
1.0 0.0433 0.8047
2.0 0.0396 0.7377

Although various damping mechanisms are functioning in the SiCp/Mg94Zn5Y1 composite, the dislocation damping mechanism is in the dominant state. The intrinsic damping mechanism and interface damping mechanism are functioning coordinately. The intrinsic damping mechanism of SiCp is different from that of the matrix; the intrinsic damping mechanism of the lamellar eutectic phase is also different from that of the primary phase[

32]. The structure and properties of different phases in the SiCp/Mg94Zn5Y1 composites are different. Therefore, the phase interface and the interface between SiCp and the matrix both slip under the application of stress and dissipate the energy. The interface between the α-Mg phase and the lamellar eutectic phase also dissipates the energy under the cyclic loading.

2.4 Comprehensive performance

In order to determine the contribution of each parameter to the coordination of the mechanical and damping properties of the composites, the efficiency coefficient method is used to evaluate the comprehensive performance of composite materials. This method is based on the principle of multi-objective planning, which determines the satisfactory value (xih) and unallowable value (xis) for each evaluation index. The satisfactory value is the highest value of the indicator, and the unallowable value is the smallest value of the indicator. The degree of each index (xi) was calculated to obtain the satisfactory value according to the equation, and the efficiency coefficient of each index can also be obtained. Then the total efficiency coefficient can be calculated according to Eq.(6) and Eq.(7), as follows:

fi=xi-xisxih-xis (6)
F=i=1nfin (7)

where the subscript i indicates various indicators; xi is an index of i indicator; n is the total number of indexes; fi is the efficiency coefficient of each index; F is the total efficiency coefficient. In the determination of the importance of each evaluation index, the total efficiency coefficient is calculated according to Eq.(7).

In this research, the SiC content in composites has an important influence on the mechanical properties and damping properties of the composites. The comprehensive properties of the SiCp/Mg94Zn5Y1 composites are evaluated by the efficiency coefficient method, and the main evaluation indexes are the compressive strength and damping capacities (Q0-1 and Qh-1), which are equally important.

The efficiency coefficients of the SiCp/Mg94Zn5Y1 composites with different SiCp contents can be calculated according to Eq.(6) and Eq.(7), and the calculation results are shown in Table 3. It can be seen that the efficiency coefficient of the composites with SiCp addition is much higher than that without SiCp addition. The comprehensive properties of the

Table 3  Efficiency coefficient F of different SiCp/Mg94Zn5Y1 composites
SiCp content/wt%00.51.02.0
Efficiency coefficient, F 0.048 0.581 0.668 0.523

SiCp/Mg94Zn5Y1 composites are firstly increased and then decreased with increasing the SiCp content. The maximum efficiency coefficient is obtained with 1.0wt% SiCp content, i.e., the optimal synergy between the quasicrystal phase and the SiCp phases is achieved with 1.0wt% SiCp content. Hence, the optimal comprehensive properties of the SiCp/Mg94Zn5Y1 composites are achieved with 1.0wt% SiCp content.

3 Conclusions

1) The SiCp added into the Mg94Zn5Y1 matrix is distributed evenly in the matrix and refines the composite grains. The SiCp/Mg94Zn5Y1 composite is composed of the α-Mg phase, quasicrystal phase, and SiCp phase.

2) The mechanical properties of the SiCp/Mg94Zn5Y1 composites are firstly increased and then decreased with increasing the SiCp content. The mechanical properties of the composites with SiCp addition are better than those of the Mg94Zn5Y1 base alloy. The 1.0wt% SiCp/Mg94Zn5Y1 composite has the compressive strength of 350 MPa, which is 12.5% higher than that of the Mg94Zn5Y1 base alloy.

3) The damping properties of the SiCp/Mg94Zn5Y1 composites are slightly improved at the low strain amplitudes. At the strain of 4×10-4, the 0.5wt% SiCp/Mg94Zn5Y1 composite shows the highest damping value of 0.011 80, inferring the high-damping material. At high strain, with increasing the SiCp content, the damping performance is firstly increased and then decreased. The damping properties of all SiCp/Mg94Zn5Y1 composites are better than those of the Mg94Zn5Y1 base alloy, and the best damping capacity is achieved with 0.5wt% SiCp addition.

4) According to the efficiency coefficient method, the optimal comprehensive properties of the SiCp/Mg94Zn5Y1 composites are achieved with 1.0wt% SiCp addition.

References

1

Liu Enyang, Yu Sirong, Zhao Yan et al. Rare Metal Materials and Engineering[J], 2017, 46(11): 3298 (in Chinese) [Baidu Scholar] 

2

Zhai J Y, Song X Y, Xu A Y et al. Metals and Materials International[J], 2021, 27(6): 1458 [Baidu Scholar] 

3

Ramalingam V V, Ramasamy P, Kovukkal M D et al. Metals and Materials International[J], 2020, 26(4): 409 [Baidu Scholar] 

4

Agnew S R, Nie J F. Scripta Materialia[J], 2010, 63(7): 671 [Baidu Scholar] 

5

Anbuchezhiyan G, Mohan B, Sathyanarayanan D et al. Journal of Alloys and Compounds[J], 2017, 719: 125 [Baidu Scholar] 

6

Wang X J, Xu D K, Wu R Z et al. Journal of Materials Science & Technology[J], 2017, 34(2): 245 [Baidu Scholar] 

7

Nie K B, Guo Y C, Deng K K et al. Journal of Alloys and Compounds[J], 2019, 792: 267 [Baidu Scholar] 

8

Wang Xiaojun, Xiang Yeyang, Hu Xiaoshi et al. Acta Metallurgica Sinica[J], 2019, 55(1): 73 (in Chinese) [Baidu Scholar] 

9

Zhang Z Y, Guo Y H, Zhao Y T et al. Materials Characterization[J], 2019, 150: 229 [Baidu Scholar] 

10

Yu W B, Li X B, Vallet M et al. Mechanics of Materials[J], 2019, 129: 246 [Baidu Scholar] 

11

Essa Y E S, Fernadez-Saez J, Perez-Castellanos J L. Journal of Testing and Evaluation[J], 2003, 31(6): 449 [Baidu Scholar] 

12

Dieringa H. Fundamentals of Magnesium Alloy Metallurgy[M]. Cambridge: Woodhead Publishing, 2013: 317 [Baidu Scholar] 

13

Zhen Xu. Thesis for Master[D]. Xi'an: Xi'an University of Science and Technology, 2016 (in Chinese) [Baidu Scholar] 

14

Prasad D S, Shoba C. Journal of Materials Research and Technology[J], 2016, 5(2): 123 [Baidu Scholar] 

15

Yan H, Xin H, Hu Q. Advanced Materials Research[J], 2010, 123-125: 35 [Baidu Scholar] 

16

Bae D H, Lee M H, Kim K T et al. Journal of Alloys and Compounds[J], 2002, 342(1-2): 445 [Baidu Scholar] 

17

Yuan G Y, Amiya K, Kato H et al. Journal of Materials Research[J], 2004, 19(5): 1531 [Baidu Scholar] 

18

Singh A, Tsai A P, Nakamura M et al. Philosophical Magazine Letters[J], 2003, 83(9): 543 [Baidu Scholar] 

19

Lü Binjiang, Peng Jian, Han Wei et al. Rare Metal Materials and Engineering[J], 2014, 43(7): 1643 (in Chinese) [Baidu Scholar] 

20

Singh A, Watanabe M, Kato A et al. Materials Science and Engineering A[J], 2004, 385(1-2): 382 [Baidu Scholar] 

21

Ma Rong. Thesis for Master[D]. Wuhan: Huazhong University of Science and Technology, 2011 (in Chinese) [Baidu Scholar] 

22

Wan Diqing, Hu Yinglin, Li Zhumin et al. Rare Metal Materials and Engineering[J], 2019, 48(1): 71 [Baidu Scholar] 

23

Yi S, Park E S, Ok J B et al. Materials Science and Engineering A[J], 2001, 300(1-2): 312 [Baidu Scholar] 

24

Kim Y K, Kim W T, Kim D H. Science and Technology of Advanced Materials[J], 2014, 15(2): 24 801 [Baidu Scholar] 

25

Zhang Y B, Yu S R, Song Y L et al. Journal of Alloys and Compounds[J], 2008, 464(1-2): 575 [Baidu Scholar] 

26

Long Qiansheng, Wang Weiwei, Ren Guangxiao et al. Foundry Technology[J], 2016, 37(5): 848 (in Chinese) [Baidu Scholar] 

27

Huang S J, Ali A N. Journal of Materials Processing Technology[J], 2019, 272: 28 [Baidu Scholar] 

28

Ren L B, Zhou M Y, Lu T H et al. Materials Science and Engineering A[J], 2020, 770: 138 548 [Baidu Scholar] 

29

Wan Diqing, Hu Yinglin, Ye Shuting et al. Journal of Functional Materials[J], 2018, 49(5): 5035 (in Chinese) [Baidu Scholar] 

30

Granato A V, Lücke K. Journal of Applied Physics[J], 1956, [Baidu Scholar] 

27(7): 583 [Baidu Scholar] 

31

Granato A V, Lücke K. Journal of Applied Physics[J], 1956, [Baidu Scholar] 

27(7): 789 [Baidu Scholar] 

32

Wan D Q, Wang J C, Wang G F et al. Transactions of Nonferrous Metals Society of China[J], 2009, 19(1): 45 [Baidu Scholar]