+高级检索
网刊加载中。。。

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

确定继续浏览么?

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

In-situ Spectral Monitoring of Microstructure During Laser Melting Deposition: A Case Study of CoCrMoW Alloy  PDF

  • Xia Yelin 1
  • Huang Zhaozhen 1
  • Chen Hanning 1,2
  • Liang Xiaodan 2
  • Shi Chuan 3
1. School of Mechanical Engineering, Tiangong University, Tianjin 300387, China; 2. School of Computer Science and Technology, Tiangong University, Tianjin 300387, China; 3. Institute for Applied Materials, Karlsruhe Institute of Technology, Karlsruhe 12, 76131, Germany

Updated:2021-11-25

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

Abstract

The relationship between the optical emission spectroscopy and the microstructure of CoCrMoW coatings by laser melting deposition was investigated. The relationship between the primary dendrite spacing and microhardness of CoCrMoW coatings at different laser powers of spectral signal was studied. A new spectral signal index, the integral area, was proposed, and the electron temperature was calculated from four discrete Cr I spectral lines. The results show that with increasing the laser power from 400 W to 1000 W, the average primary dendrite spacing is increased from 3.426 μm to 7.420 μm, and the microhardness HV0.2 is reduced from 3461 MPa to 3095 MPa. The integral area and electron temperature of the spectral is also increased with increasing the laser power. The primary dendrite spacing of the coating is positively linear-related with the spectral signal, while the microhardness of the coatings is negatively linear-related with the spectral signal. In this research, compared with the electron temperature, the integral area shows a better potential for the prediction of the primary dendrite spacing and microhardness of coatings.

Science Press

CoCr alloy has been used in aero-engine, turbine blade[

1,2], and bioimplants, such as denture and artificial joints[3,4], due to its high-temperature strength, high wear resistance, excellent corrosion resistance, and biocompatibility. Traditional manufacturing methods used in metal implants, including casting, cutting, and hot pressing, can lead to high defect, low material utilization, a long production cycle, and high cost[5,6]. In recent years, the laser-based additive manufacturing (AM) methods, such as selective laser melting (SLM) and laser melting deposition (LMD), have attracted wide attention, because they can overcome many shortcomings of traditional manufacturing methods and offer a way for feasible manufacturing of parts with complex structure. Thus, the laser-based AM has gradually become one of the most important methods for repairment and fabrication of objects with 3D structure[7].

It is known that the mechanical properties of CoCr alloys, such as tensile properties, wear resistance, and microhardness, are determined by their microstructure. For example, Qian et al[

8] reported that the Vickers hardness and tensile strength of CoCrMo alloy fabricated by SLM are better than those made by casting. The improvement is caused by the heterogeneous multiscale (from macroscale to nanoscale) textures formed by SLM. Lee et al[9] found that the microstructures (crystal orientation) determine the anisotropic mechanical properties of SLMed CoCrW alloys. Hedberg et al[10] found that the SLMed CoCrMo specimens show better corrosion resistance than the cast ones do due to the rapid cooling and fine cellular structure. Ram et al[11] found that the wear resistance of laser cladded CoCrMo alloy is affected by the number and morphology of cemented carbide particles in the microstructure. In general, almost all these metallographic specimens need to be carefully prepared (machined and etched) to characterize the microstructure. Ideally, the manufactured specimen should be non-destructive and prepared without complex preparation procedures.

Recently, in order to improve the quality and reliability of parts manufactured by laser-based AM, much attention is paid on real-time monitoring techniques[

12,13]. For example, the optical emission spectroscopy is applied to real-time monitoring using the plasma produced by the interaction between the laser and metal powder. Song[14] and Shin[15] et al established the relationship between plasma parameters and elemental composition by analyzing element spectral lines, which achieved the real-time detection of elements the LMD process. Lednev et al[16,17] showed the in-situ quantitative analysis of key elements by laser-induced breakdown spectroscopy in LMD. Ya et al[18] studied the spectral signal of the onset and extent of metallic bonding in LMD, and suggested that the electron temperature shows a better correlation with the dilution. Lough et al[19] developed a set of equipment for in-situ spectral detection during SLM and found that the temperature of the vapor plume is correlated with the melting pool size. Chen et al[20,21] established the LMD spectral diagnosis system to investigate the relationship between processing parameters and plasma radiation intensity and the related influence on forming quality. In brief, the spectral detection is mainly used for determination of element composition and surface quality analysis in AM. However, it is more important to establish a relationship among processing, structure, and mechanical properties.

In this research, an in-situ spectral monitoring system was established for manufacturing CoCrMoW alloy using LMD. The spectra of the melting pool at different laser powers were analyzed. The relationship among the spectral signal of the melting pool, primary dendrite spacing, and microhardness at different processing parameters using the integral area of spectral curve and electron temperature was established.

1 Experiment

In this study, CoCrMoW powders with size ranging from 45 μm to 106 μm were used. Fig.1 shows the morphology of CoCrMoW alloy powers observed by scanning electron microscope (SEM). The elements of CoCrMoW alloy included Co, Cr, Mo, W, and Si, as shown in Table 1. CoCrMoW alloy powders were deposited on the Q235 steel substrate with a dimension of 210 mm×60 mm×10 mm.

Fig.1 SEM morphology of CoCrMoW alloy powder

Table 1 Chemical composition of CoCrMoW alloy power (wt%)
CrMoSiWCo
24.52 4.82 0.61 5.12 Bal.

In this study, a beam of continuous wave (CW) fiber laser (wavelength of 1064 nm) was used for LMD. The powder was fed into the melting pool through a coaxial nozzle using a RC-PGF-D powder feeder (Raycham, China). In the multi-track cladding experiment, the laser power was increased from 400 W to 1000 W with a step size of 200 W at a scanning speed of 8 mm/s, and the spot diameter was 2 mm. An overlap ratio of 50% and a powder feeding rate of 25 g/min were selected.

The spectra produced during laser cladding were collected by a flame miniature spectrometer (Ocean Optics, Flame-S) with an entrance slit of 25 μm. The spectrometer used a linear silicon charge transfer device array (Sony ILX 511) with 2048 pixels and an integration time of 10 ms. The detector was fixed on the laser head with an angle of about 80°, i.e., about 10° from the substrate, as indicated in Fig.2. The detector was focused on the melting pool. The spectrum was converged by the focusing lens onto the optical fiber with a diameter of 600 μm, then guided into the spectrometer, and finally converted into the required digital signal. OceanView software was used to control the parameters of spectrum acquisition and subsequent data processing.

Fig.2 Schematic diagram of experiment device

Metallographic specimens were prepared to study the relationship between the spectral signal and the microstructure of the coating. SEM (ZEISS Sigma 300) and X-ray diffraction (XRD, D/MAX-2500, Cu target, 40 kV, 140 mA) were used to analyze the microstructure and phases of the coating. Microhardness of the coating was measured by HV-1000 Vickers hardness tester at load of 0.2 kg for loading time of 10 s. The collected spectra were analyzed offline to identify the discrete spectra, and then the integral area and electron temperature signals were used for further study.

2 Results and Discussion

2.1 Microstructure of multi-track coating

G=2KT-T02εP (1)
R=Vscosθ (2)

where K is the thermal conductivity of alloy, T is liquid temperature of alloy, T0 is initial temperature of substrate, ε is absorption efficiency of laser, P is laser power, Vs is laser scanning speed, θ is the included angle between Vs and R. During LMD, the cooling rate GR determines the grain size of solidified structure with an inversely proportional relationship. It is also known that the morphology of the solidified structure is affected by G/R. With the decrease of G/R, the degree of constitutional supercooling is increased, and the solidification mode changes from planar to cellular and columnar crystal, and finally to equiaxed crystal[

22]. The primary dendrite spacing δ reflects the extent of refinement of the microstructure, which is related to G and R. It can be expressed by Eq.(3), as follows[23]:

δ=aGR (3)

where a is a coefficient.

It can also be seen from Fig.3a1~3a4 that the cracks appear at 400 and 600 W, but disappear at 800 and 1000 W. This may be due to the high-temperature gradient G; when the laser power is low, it results in the high residual stress of the coating and then the crack[

24].

Fig.3 shows SEM images of CoCrMoW coatings prepared at different laser powers with a constant scanning speed of 8 mm/s. It is obvious that the average coating thickness is increased from 0.37 mm to 1.23 mm with increasing the laser power from 400 W to 1000 W, as indicated in Fig.3a1~3a4. Due to the rapid melting and solidification of powders in LMD, the microstructure can be determined by constitutional supercooling theory[

22], indicating that the liquid temperature is lower than the liquidus temperature, which is related to temperature gradient G and solidification rate R. It can be expressed by Eq.(1) and Eq.(2), as follows[23]:

At the bottom part of coating, as shown in Fig.3b1~3b4, the heat is mainly dissipated through the substrate. Columnar dendrite crystals grow from the bottom of the melting pool because of the large temperature gradient G and a small extent of constitutional supercooling in this area. The growth direction of columnar dendrite crystals is mostly along the temperature gradient direction which is perpendicular to the boundary of the melting pool[

25].

The average primary dendrite spacing is determined by the number of grains passing through a straight line of a certain length, as shown in Fig.4. When the laser power is 400 W, the average primary dendrite spacing of columnar dendrite crystals in bottom part of coating is 4.538 μm, as listed in Table 2. With increasing the laser power to 600, 800, and 1000 W, the average primary dendrite spacing in bottom part of coating is 4.969, 5.856, and 7.020 μm, respectively. This change can be explained by the fact that the cooling rate GR is decreased with increasing the laser power, leading to a coarser microstructure. Meanwhile, the primary dendrite spacing is increased with decreasing the GR, according to Eq.(3).

Table 2 Primary dendrite spacing of CoCrMoW coatings at different laser powers (μm)
RegionLaser power/W
4006008001000
Bottom 4.538 4.969 5.856 7.020
Middle 2.939 5.176 5.355 7.014
Top 2.801 4.544 5.176 8.225
Average value 3.426 4.896 5.462 7.420

Fig.3c1~3c4 show the microstructure of middle part of coatings. When the laser power is 400 and 600 W, the microstructures are mainly equiaxed crystals with a small number of columnar dendrite crystals. This is because the G/R value decreases rapidly from the bottom part to the middle part of the melting pool. Correspondingly, the columnar dendrite crystals are transformed into equiaxed crystals[

26]. When the power is 800 and 1000 W, uniform equiaxed crystals are formed in the middle of the coating, and the grain size of equiaxed crystal continues to grow. For equiaxed crystals, the primary dendrite spacing can be regarded as their grain size, so it is observed that the average primary dendrite spacing is increased from 2.939 μm to 7.014 μm with increasing the laser power from 400 W to 1000 W, as shown in Table 2.

Fig.3d1~3d4 show the microstructures of the top part of coatings. It can be seen that almost all the microstructures in this region are equiaxed crystals at different laser powers. At the top part of coatings, the main heat dissipation mode is heat transfer inside the coating and heat convection on the surface. Therefore, the cooling rate and solidification rate in this region is fast and fine equiaxed grains are obtained. The average primary dendrite spacing of equiaxed crystals is increased significantly from 2.801 μm to 8.225 μm with increasing the laser power from 400 W to 1000 W, as shown in Table 2. In general, the average primary dendrite spacing is increased from 3.426 μm to 7.420 μm with increasing the laser power from 400 W to 1000 W.

2.2 Phase analysis

The phase component of CoCrMoW coatings at different powers is analyzed by XRD. The influence of laser power on phase component is presented in Fig.5a. The XRD patterns show the presence of γ and ε phases with different grain orientations[

27-30], and the change of laser power has little effect on the types of the main phases. Fig.5b shows that the peak position of the phase shifts with the increase of laser power. When the power is 400, 600, 800, and 1000 W, the peak position is 43.62°, 43.84°, 43.56°, and 43.78°, respectively. This may be due to the shrinkage of lattice anisotropy and lattice distortion caused by residual stress, which makes the diffraction peak shift to a small angle[31]. However, when the residual stress is released, the diffraction peak shifts to a large angle at laser power of 400 and 600 W. The diffraction angle is related to the inter-plane distance, according to Bragg's Law[32], as expressed by Eq.(4):

2dsinθ=

(n=1,2,3…) (4)

where d is the inter-plane distance, θ is half scattering angle, and λ is the X-ray wavelength. The inter-plane distance shows negative correlation with the diffraction angle. Consequently, the inter-plane distance of CoCrMoW coatings decreases at first, then increases, and decreases finally.

Meanwhile, the full-width half-maximum (FWHM) of the diffraction peaks is decreased with increasing the laser power from 400 W to 1000W, as shown in Fig.5b. According to the Scherrer formula[

33], as expressed by Eq.(5):

L=Kλβcosθ (5)

where L is the average grain size, K is a constant related to grain shape, and β is FWHM. Therefore, FWHM is inversely proportional to the average grain size, and the average grain size is decreased with increasing the laser power.

2.3 Microhardness

Fig.6a shows the microhardness of the coating at different positions, and a clearly increasing microhardness can be observed in the coating. The average microhardness HV0.2 is reduced from 3461 MPa to 3095 MPa with increasing the laser power from 400 W to 1000 W. When the laser power increases from 400 W to 1000 W, the grain size increases as the solidification rate decreases, as shown in Fig.3. According to the Hall-Petch relationship[

34], it is expected that the microhardness gradually decreases, as shown in Fig.6b. Interestingly, the microhardness decreases almost linearly with a constant decrement of about 4% after increasing the laser power by 200 W each time. The result of a linear relationship can be a good sign for accurate prediction.

2.4 Spectra of melting pool

2.4.1 Integral area of spectral curve

The original spectra of CoCrMoW alloys obtained during LMD at different laser powers are shown in Fig.7a. Under the laser cladding condition, the major component of the optical emission is Planck's radiation of the melt pool surface, and thus a possible reason for the continuous background is the thermal radiation from the surface. The spectral lines are produced by the transition of electrons from high energy level m to low energy level n[

35]. The spectral intensity Imn can be described as follows[18]:

Imn=NmAmnhυmn (6)

where Nm is the population of the upper state m (cm-3), Amn is the transition probability, h is Planck's constant, and υmn is the frequency of a photon. Assuming that the plasma is at local thermal equilibrium, then the population of the upper state can be expressed by Boltzmann law[

18], as expressed by Eq.(7):

NmN=gmZTeexp-EmkTe (7)

where N is the population of ground state particles (cm-3), Te is the electron temperature of the plasma (K), Z(Te) is the partition function at Te, gm is the degeneracy of upper state m, Em is the energy of upper state m (cm-1), and k is the Boltzmann constant.

It is known that in thermal radiation field, Planck's formula describes the relationship between the emissivity and frequency of electromagnetic radiation emitted from a blackbody. Stefan Boltzmann's law states the total energy radiated by a blackbody surface per unit area across all wavelengths per unit time, which can be obtained by integrating Planck's formula over the entire spectrum[

36].

The spectral intensity Imn in the whole monitored wavelength range is integrated without background subtraction, while the continuous background radiation is treated as noise which needs to be deducted in normal spectral analysis. The total integral area IA of the spectra can be obtained, as shown in Fig.7b, which is regarded as the total radiation intensity in this range. With increasing the laser power from 400 W to 1000 W, IA is increased almost linearly from 8.392×105 to 14.879×105. It can be seen that the spectral continuum is the main factor for the increase of the integral area, which is closely related to the temperature of plasma.

2.4.2 Electron temperature

Electron temperature is an important parameter to characterize plasma produced by the interaction between laser and metal powder, which reflects the thermal input of laser. Electron temperature can be calculated by the relative intensities of discrete spectral lines by the Boltzmann plot method[

18]. In the calculation, the continuous spectrum should be eliminated to reduce the error. Thus, the iterative stripping method[37] was used for background treatment.

By comparing the intensity Ii of a wavelength i and the average value of two adjacent sampling points in the spectrum, the rapid changing information in the spectrum can be stripped[

37]. The flow chart of the background subtraction algorithm using iterative stripping is shown in Fig.8. Then, by substituting Eq.(7) into Eq.(6), after some mathematical treatments, Eq.(6) can be expressed as follows:

lnImnλmngmAmnlnNhcZTe-1kTeEm (8)

where λmn is the wavelength of the spectral line (nm), c is the light speed (m/s). Eq.(8) shows a linear relationship between ln(Imnλmn/gmAmn) and Em, where ln(Nhc/Z(Te)) is a constant.

Fig.8 Algorithm flow chart of iterative stripping

Thus, the electron temperature Te can be estimated from the inverse of the slope -1/(kTe) fitted from the calculated data of multiple spectral lines, as shown in Fig.9a. The parameters used for calculation in Eq.(8) can be found in National Institute of Standards and Technology (NIST) database. To ensure the calculation accuracy, the energy level difference between the two spectral lines used for calculation should be as large as possible. Table 3 shows the spectral data of four selected Cr I lines at 400 W because some parameters of other lines required in Eq.(8) are not specified in NIST database. The calculated electron temperatures in Fig.9b are increased almost linearly from 2286.206 K to 2463.742 K with increasing the laser power from 400 W to 1000 W.

Table 3 NIST data of four selected Cr I lines at 400 W
λmn/nmAmn/s-1gmEm/cm-1
534.044 69 1.45×107 3 46 448.629
544.240 69 3.90×106 5 45 966.360
569.474 04 1.40×107 7 48 661.508
571.981 57 4.96×105 9 41 782.131

2.4.3 Spectral intensity

Fig.10 shows the intensity of spectral lines at different laser powers. Three spectral lines with the wavelength of 520.017, 580.022, and 605.513 nm were selected because of their appropriate intensities, FWHM, as well as less interference. It can be clearly seen that the intensities of three chosen spectral lines are increased with increasing the laser power. The reason for this phenomenon can be explained by Eq.(6) and Eq.(7). Once the spectral line is determined, its intensity is mainly related to the population of ground state particles N and the electron temperature Te. When the laser power increases and other processing parameters remain unchanged, there is no clear change of the population of ground state particles N. Thus, the increasing electron temperature, as shown in Fig.9b, is the main factor in increasing the spectral line intensity[

20].

2.5 Relationship between spectrum signal and microstruc-ture

The relationship between the spectral related information (integral area IA and average electron temperature Te) and microstructures (average primary dendrite spacing δ and average microhardness) of the CoCrMoW coatings are shown in Fig.11. It can be clearly seen that all the four fitting curves have linear relationships. Moreover, the integral area and electron temperature are positively correlated with the primary dendrite spacing δ of CoCrMoW coatings, as shown in Fig.11a and 11c; while they are negatively correlated with the microhardness of CoCrMoW coatings, as shown in Fig.11b and 11d. This is because the microstructure and integral area are closely related to the melting pool. The microstructure is formed after the melting pool solidifies, and the spectral signal radiated from the melting pool is determined by the state of the melting pool which is related to the heat input of laser. The increase of heat input leads to the increase of electron temperature, which promotes the collisional excitation of particles, increases the number of excited state particles, radiates more photons through energy level transition, and finally leads to the increase of the spectral intensity radiated from the melting pool. Thus, the larger integral area IA and higher electron temperature Te can directly reflect the increasing heat input, which leads to the grain coarsening and reduces the microhardness of the coatings. The results show that the linearity of relationship corresponding to the integral area is more obvious than that corresponding to electronic temperature due to the larger R-Square, namely higher accuracy.

Therefore, the spectral signals can be used to predict the average primary dendrite spacing δ and average micro-hardness of CoCrMoW coatings processed by LMD. In addition, the integral area can directly reflect the heat input determined by the laser power without complex calculations. Thus, it can also be used to monitor LMD after calibration, and to avoid the inaccurate output power caused by the damage of optical lens or laser attenuation.

3 Conclusions

1) The integral area of spectral curve proposed in this study and the electron temperature are increased almost linearly with increasing the laser power. However, the microhardness of CoCrMoW coating is decreased with increasing the laser power due to the coarsening of microstructure.

2) The spectral analyses indicate that the integral area show good capability for predicting the average primary dendrite spacing and the average microhardness of CoCrMoW coating. The results suggest that the spectral signal can be used to predict the microstructure (primary dendrite spacing and microhardness) of the coating.

References

1

Sun S H, Koizumi Y, Kurosu S et al. Acta Mater[J], 2014, 64: 154 [Baidu Scholar

2

Barekat M, Shoja Razavi R, Ghasemi A. J Laser Appl[J], 2016, 28(4): 42 005 [Baidu Scholar

3

Santecchia E, Gatto A, Bassoli E et al. J Alloy Compd[J], 2019, 797: 652 [Baidu Scholar

4

Lu Y J, Guo S, Yang Y et al. J Alloy Compd[J], 2018, 730: 552 [Baidu Scholar

5

Zhang M K, Yang Y Q, Song C H et al. J Alloy Compd[J], 2018, 750: 878 [Baidu Scholar

6

Liverani E, Fortunato A, Leardini A et al. Mater Des[J], 2016, 106: 60 [Baidu Scholar

7

Ahmadi S M, Hedayati R, Li Y et al. Acta Biomater[J], 2018, 65: 292 [Baidu Scholar

8

Qian B, Saeidi K, Kvetková L et al. Dent Mater[J], 2015, [Baidu Scholar

31(12): 1435 [Baidu Scholar

9

Lee H W, Jung K H, Hwang S K et al. Mater Sci Eng A[J], 2019, 749: 65 [Baidu Scholar

10

Hedberg Y S, Qian B, Shen Z J et al. Dent Mater[J], 2014, [Baidu Scholar

30(5): 525 [Baidu Scholar

11

Ram G D J, Esplin C K, Stucker B E. J Mater Sci Mater Med[J], 2008, 19(5): 2105 [Baidu Scholar

12

Song L J, Huang W K, Han X et al. IEEE Trans Ind Electron[J], 2016, 64(1): 633 [Baidu Scholar

13

Stutzman C B, Nassar A R, Reutzel E W. Addit Manuf[J], 2018, 21: 333 [Baidu Scholar

14

Song L, Mazumder J. IEEE Sens J[J], 2011, 12(5): 958 [Baidu Scholar

15

Shin J, Mazumder J. Opt Laser Technol[J], 2018, 106: 40 [Baidu Scholar

16

Lednev V N, Sdvizhenskii P A, Asyutin R D et al. Opt Express[J], 2019, 27(4): 4612 [Baidu Scholar

17

Lednev V N, Sdvizhenskii P A, Asyutin R D et al. Addit Manuf[J], 2019, 25: 64 [Baidu Scholar

18

Ya W, Konuk A R, Aarts R et al. J Mater Process Technol[J], 2015, 220: 276 [Baidu Scholar

19

Lough C S, Escano L I, Qu M et al. Solid Freeform Fabrication[C]. Austin: University of Texas, 2018: 2192 [Baidu Scholar

20

Chen B, Yao Y Z, Tan C W et al. Int J Adv Manuf Technol[J], 2018, 96(9-12): 4231 [Baidu Scholar

21

Chen B, Yao Y Z, Tan C W et al. Transactions on Intelligent Welding Manufacturing[M]. Singapore: Springer Singapore, 2018: 60 [Baidu Scholar

22

Kou S. Welding Metallurgy[M]. New York: Wiley, 2002 [Baidu Scholar

23

Huang Y J, Zeng X Y, Hu Q W et al. Appl Surf Sci[J], 2009, [Baidu Scholar

255(7): 3940 [Baidu Scholar

24

Stutzman C B, Nassar A R, Reutzel E W. Addit Manuf[J], 2018, 21: 333 [Baidu Scholar

25

Wei H L, Elmer J W, DebRoy T. Acta Mater[J], 2017, 133: 10 [Baidu Scholar

26

Hemmati I, Ocelík V, De Hosson J T M. J Mater Sci[J], 2011, [Baidu Scholar

46(10): 3405 [Baidu Scholar

27

Limmahakhun S, Oloyede A, Sitthiseripratip K et al. Mater Des[J], 2017, 114: 633 [Baidu Scholar

28

Barekat M, Razavi R S, Ghasemi A. J Mater Eng Perform[J], 2017, 26(7): 3226 [Baidu Scholar

29

Takaichi A, Nakamoto T, Joko N et al. J Mech Behav Biomed Mater[J], 2013, 21: 67 [Baidu Scholar

30

Béreš M, Silva C C, Sarvezuk P W C et al. Mater Sci Eng A[J], 2018, 714: 36 [Baidu Scholar

31

Ferreira N G, Abramof E, Corat E J et al. Carbon[J], 2003, [Baidu Scholar

41(6): 1301 [Baidu Scholar

32

Wang Z Q, Denlinger E, Michaleris P et al. Mater Des[J], 2017, 113: 169 [Baidu Scholar

33

Vladescu A, Dinu M, Braic M et al. Ceram Int[J], 2015, 41(6): 8051 [Baidu Scholar

34

Wei X F, Li W J, Liang B J et al. Tribol Int[J], 2016, 97: 212 [Baidu Scholar

35

Ancona A, Spagnolo V, Lugarà P M et al. Appl Opt[J], 2001, [Baidu Scholar

40(33): 6019 [Baidu Scholar

36

Kawai J, Ishii H. Spectrochim Acta Part B: At Spectrosc[J], 2005, 60(12): 1586 [Baidu Scholar

37

Wang Jingge, Li Xinzhong, Li Hehe et al. Spectrosc Spectr Anal[J], 2018, 38(1): 276 (in Chinese) [Baidu Scholar