ISIJ International
Online ISSN : 1347-5460
Print ISSN : 0915-1559
ISSN-L : 0915-1559
Special Issue on "Frontier in Characterization of Materials and Processes for Steel Manufacturing"
Rapid Analysis of Steel Powder for 3D Printing Using Laser-Induced Breakdown Spectroscopy
Haorong GuoZhongqi FengMinchao Cui Yoshihiro DeguchiLiang TanDacheng ZhangChangfeng YaoDinghua Zhang
Author information
JOURNAL OPEN ACCESS FULL-TEXT HTML

2022 Volume 62 Issue 5 Pages 883-890

Details
Abstract

In recent years, metal additive manufacturing technology (also known as 3D printing) has achieved dramatic development and has gradually occupied an indispensable position in modern manufacturing. However, the purity and uniformity of the metal powder, which is the raw material for additive manufacturing, generate significant impacts on the performance of the product. Therefore, a rapid analysis method of the steel powders based on laser induced breakdown spectroscopy (LIBS) is investigated for the potential applications in 3D printing process. In the current experiment, two pretreatment methods of the novel adsorbing process and the traditional tableting process are employed for performing comparative study. Through the adsorbing process with conductive tape, the signal intensity of LIBS spectrum has been significantly enhanced. The crater topography demonstrates that the adsorption band samples can effectively overcome the shock wave effect of the laser induced plasma. The calibration curve of LIBS measurement indicates that the pellet samples can be easily affected by the plasma shock wave. Differently, the adsorption band samples are slightly influenced by the plasma shock wave. The determination coefficient and data stability of the calibration curve are better in the condition of adsorption band samples. Through the simulated quantitative analysis with totally 11 artificial steel powder samples, the analytical performance of the two sample pretreatment processes is compared. Manganese (Mn) is selected as the target element with the matrix element of iron (Fe). When compared with the pellet samples, the average value of Relative Error of Prediction (REP) is reduced from 7.73% to 1.36% and the average value of Relative Standard Deviation (RSD) is reduced from 20.34% to 6.22% for the adsorption band samples. Based on the experimental results in the current work, it can be verified that the adsorbing process is a promising sample pretreatment method for steel powders. Moreover, this provides an effective sample pretreatment method for LIBS analysis of steel powder for 3D printing.

1. Introduction

Recent theoretical developments have revealed that metal additive manufacturing technology is a component manufacturing process, which is also called 3D printing. The current technology has been influential in the field due to its simple manufacturing process, low production cost, as well as short research and development cycle. Additionally, it has the advantages of free design, mass customization, waste minimization, and rapid prototyping.1,2,3) Meanwhile, it has been extensively adopted in the field of steel industry,4) aerospace,5) shipbuilding,6) and automotive applications.7) At present, with the development of the additive manufacturing industry, new materials and devices continue to appear. In the future, it will move toward precision improvement, and cost reduction, etc.8) However, a challenging problem that arises in this domain refers to the quality stability of parts. Therefore, it is of importance to explore the composition uniformity of raw materials for the performance of 3D printing parts.9,10)

Although a whole range of different approaches to the problem are available, the uniformity of the composition of the metal powder, the contamination of impurities, the morphological characteristics of the powder, etc. will influence the stability of the 3D metal printing parts. Moreover, it is prone to fluctuations in powder composition during the milling process. On this basis, these constraints make the problem difficult to solve some problems such as element enrichment, segregation, and oxide inclusions in local areas. Thus, a possible solution to the problem at hand is the analysis and control of the raw material powder composition in 3D printing process, which is a critical point to improve the quality and stability of the parts. As we all know, there are several common kinds of powder analysis methods, such as X-ray Fluorescence Spectroscopy (XRF),11) X-ray Tomography (XRT),12) Inductively Coupled Plasmay (ICP),13) Scanning Electron Microscope (SEM),14) Energy Chromatography Element Analysis (EDS)15) and X-ray Fluorescence Spectroscopy (XFS).16)

Nevertheless, the application of laser-induced breakdown spectroscopy to analyze the chemical elements of materials is widely regarded as an efficient method.17) It is characterized by not being restricted by the sample state, multi-element simultaneous analysis, fast analysis speed, and little sample damage.18) Various papers have been published to report the application of rapid analysis of powder samples, such as coal powder,19,20) soil,21,22) plants,23) food,24) cosmetic powder,25) geology26) and medicine.27) Especially in metal powder analysis of LIBS, K. Novotný et al.28) have made use of embedding hard tungsten carbide into the resin, and pressing the powder material with a hydraulic press to prepare particles, which can quantitatively determine the element of niobium, titanium, tantalum, and cobalt. Marco Aurelio Sperança et al.29) have proposed a new method for quantifying Al, Cr, Fe, Mg, Mn, and Ni in nickel-containing minerals through LIBS. Firstly, solid mineral samples are mixed, which are ground grinding into a fine powder. Water in different proportions is used to make a slurry. Then, polyvinyl alcohol solution is added. Secondly, after being placed in an oven at 50°C for 2 hours, the sample fixed in the polymer film is subjected to be adopted for exploring chemical component, aiming to reduce signal fluctuations and improve the uniformity of the sample. V. N. Lednev et al.30) have used double-sided tape to attach loose metal powder. As a result, LIBS with a fast-real-time on-site analysis function can better match the automatic, fast, and direct-production characteristics of 3D printing.31)

The reason is that most of the preparation methods for powder samples is difficult to achieve rapid online measurement. Therefore, we present a new method for sample preparation to enhance the sensitivity and stability of signal detection and perform well in preventing powder from being washed away. In this study, it is shown that the diagram of the sample pretreatment processes and pictures of real products are given, respectively, the pellet sample and the adsorption band sample. The enhancement of LIBS spectrum signal is discussed by improving sample preparation method. The different crater topography formed by the tableting process and adsorbing process of the original powder after being excited by the laser can be observed. The trend of calibration curves was analyzed. A total of nine samples from S1 to S9 are cross-validated. The value of REP and RSD that affects precision and uncertainty measurement can be obtained. Finally, we find that all aspects of the adsorption band samples generate the better results. More importantly, the adsorption band is easier to realize continuous and rapid sample preparation process, and it is more convenient to apply to industrial sites. The results indicate that the new preparement method exerts a critical role in the rapid analysis of steel powder for 3D printing based on laser-induced breakdown spectroscopy.

2. Experimental Setup

2.1. Device Configuration

Figure 1 shows the schematic diagram of LIBS system which is employed to measure the steel powder samples in the current work. The pulsed laser beam is generated by a Nd: YAG nanosecond laser (LS-2134UTF, LOTIS TII Co. Ltd., Republic of Belarus). At Q-switched mode of the laser, the diameter of laser beam is 6 mm and the pulse width is 5.4 ns. In this experiment, the laser is operated at 1064 nm with the pulse energy of 20.5 mJ and the frequency of 10 Hz. Besides, the laser beam is reflected by a mirror (reflectivity >99% at 1064 nm), and is then passed through a long pass filter (>450 nm), and is finally concentrated on the sample by a quartz lens (focal length 200 mm). The plasma light is collected in the reverse direction. It is reflected by the long pass filter and focused on the end of an optical fiber by a quartz lens (focal length 120 mm). After transmission through the optical fiber, the plasma light is dispersed by a grating spectrometer (NP-250-2M, SOL Instrument Co. Ltd., Republic of Belarus), which is detected by an ICCD (iStar DH334T-18U-03, Andor Technology Co. Ltd., England). A grating with 3600 lines per millimeter was installed in the spectrometer. Thus, the detection system has achieved a spectral resolution of 0.0097 nm in the wavelength range of 400–410 nm. In addition, the trigger timing between the devices is controlled by a delay generator (Model DG645, Stanford Research System Co. Ltd., USA). The detection data of ICCD is transferred to a computer for storage and analysis. In this work, the parameters of LIBS acquisition are set as: delay time 3000 ns, gate width 1000 ns, 50 times accumulation and 5 repeated experiments under each condition.

Fig. 1.

Experimental setup. (Online version in color.)

2.2. Sample Pretreatment

In this work, the steel powder samples were made of pure iron (Fe) powder and pure manganese (Mn) powder. The Fe powder and Mn powder were provided by Beijing Research Institute for Nonferrous Metals, China. The certified purity of Fe powder reaches 99.995% and the certified purity of Mn powder is 99.99%. With the size between 0.5 μm and 10 μm, both the powders are spherical particles. The artificial steel powders were mixed by weighting a certain mass of pure Fe powder and pure Mn powder respectively. In the meanwhile, the powders were mixed in a ball mill for 2 hours to ensure uniform mixing of the two powders. Table 1 shows the number of the artificial steel powders and the corresponding Mn and Fe content. For example, the S6 sample was made by mixing 3 g of pure Mn powder and 97 g of pure Fe powder in the ball mill for 2 hours. Therefore, the theoretical value of Mn content in the S6 sample is 3%.

Table 1. Elemental composition of the artificial steel powders.
Sample No.Content (%)
MnFe
S00.0100.0
S10.599.5
S21.099.0
S31.598.5
S42.098.0
S52.597.5
S63.097.0
S73.596.5
S84.096.0
S94.595.5
S105.095.0

In this study, these artificial steel powders were further pretreated into different states for LIBS measurement. Figure 2 shows the diagram of two pretreatment processes employed in the current work. Specifically, the first method is called tableting process. The original steel powder sample is put into the mold and then pressed into the pellet sample through 2k N pressure. This method has already been extensively applied in the LIBS measurements of powder samples, such as soils32) and anthracites.33,34) The other method is named as adsorbing process. The conductive tape is employed in the proposed method. Conductive tape is usually used in the pretreatment of electron microscope samples. It is made by uniformly coating highly conductive adhesive on both faces of copper band (mostly the thickness of the adhesive is a few microns). The conductive tape can bind the powders to form a dense sample surface. Additionally, the conductive tape can conduct electric charge and avoid the charging and discharging effects of thin samples. In this work, the steel powders were sprayed on the one face of the conductive tape and the surface with steel powders was compacted by a cylindrical tool. Subsequently, the other face of the conductive tape was pasted on an aluminum band with the thickness of 2 mm. Finally, the steel powders pretreated by the adsorbing process are called adsorption band samples, as displayed in Fig. 2. To the best of our knowledge, the adsorption band method is the first time to be investigated for the LIBS measurement of steel powder. From the perspective of the LIBS application in 3D printing filed, the adsorption band method shows great potential in industrial application. In comparison with the tableting process, the adsorbing process uses simpler tools, which is conducive to the high efficiency and automation of sample pretreatment. It is significant to explore the effect of adsorbing process on the results of LIBS measurement of steel powders.

Fig. 2.

Diagram of the sample pretreatment processes. (Online version in color.)

3. Results and Discussion

3.1. Signal Enhancement

In this study, the steel powder samples were measured in the three states. The first state is directly spread on the sample stage as the original powder state, the second state is the pellet sample by the tableting process, and the third state is the adsorption band sample by the adsorbing process. For the sake of clarity, the photos of the three states are shown in Fig. 3. To be specific, Fig. 3(a) shows the original powder sample without any pretreatment. Figure 3(b) presents the original steel powder pressed into a tablet form a mold. Figure 3(c) exhibits the powder adsorbed by conductive tape and turns into the adsorption sample. Besides, the dense top layer can be made by the above-mentioned two sample preparation methods. Figure 4 shows the recorded spectra from different samples. Several sharp spectral lines can be distinguished from these spectra. At the same time, the related parameters of the spectral lines are listed in Table 2 according to the NIST database. It can be observed from Fig. 4 that the signals of the three samples are significantly different. In the three states of the sample, the iron spectrum at different wavelengths and the manganese spectrum all show different signal intensities. For example, LIBS can not obtain effective spectral signal in the powder state. However, useful signal lines can be obtained when the original powder sample goes through the tableting process or adsorbing process. In terms of the original powder, the overall signal is relatively flat and there is no obvious peak attribute to the noise with continuous background, which is impossible to be explored. In the state of the pellet sample, the isolated peak is distinct, and the spectral signal is enhanced compared to the original powder. Regarding the adsorption band, the spectral intensity has been remarkably improved compared with other sample preparation methods. Obviously, the peak of the spectral line is higher. The individual atomic spectrum lines can also be clearly distinguished.

Fig. 3.

Pictures of the samples prepared by different pretreatment processes. (a) As powder, (b) pellet and (c) adsorption band. (Online version in color.)

Fig. 4.

Measured spectra of different sample pretreatment processes. (Online version in color.)

Table 2. Parameters of the atomic emission lines and signal enhancement.
Wavelength (nm)Ei (cm−1)Ek (cm−1)Aki (107 s−1)IT/IP
Fe I 400.52412560.9337521.162.043.1
Fe I 404.58111976.2436686.188.623.1
Fe I 406.35912560.9337162.756.653.4
Fe I 407.17412968.5537521.167.643.6
Mn I 403.076024802.251.73.4
Mn I 403.307024788.051.653.4
Mn I 403.449024779.321.583.5

Symbol annotations: Ei: Energy of lower level, Ek: Energy of upper level, Aki: Transition probability of spectral line, IT/IP: Signal enhancement ratio between adsorption band samples and pellet samples.

In addition, the signal enhancement ratio between adsorption band samples and pellet samples (IT/IP) is also listed in Table 2. According to Fig. 4 and Table 2, the intensities of the spectral lines are enhanced 3–4 folds by the sample pretreatment with adsorbing process. It is significant to discuss the reason for the above experimental results. The reason for the significant difference in the spectral signal is considered to be the shock wave effect of plasma during LIBS excitation. The shock wave of plasma will disperse the powder from the focal area when the powder is in a loose state. As a result, in the original powder state, the laser beam cannot continue to interact with the sample materials, leading to a weaker spectral signal during LIBS measurement. To improve the signal, it is vital to make the powder sample possess a dense surface. According to the experimental results presented in Fig. 4, both the pellet method and adsorption band method can significantly improve the signal quality. However, the adsorption band method exhibits a stronger signal. This result indicates that the adsorption band sample prepared by the adsorbing process is more suitable for LIBS excitation. Through the adsorbing process, the powder particles are fixed on the conductive tape. Due to The metal powder is bound to the surface of the aluminum plate substrate through a layer of binder. The interaction process between the laser beam and the surface of the adsorption band sample is closer to the interaction process between the laser beam and the surface of the solid sample. The adsorption band method yields a dense surface for LIBS measurement which significantly enhances the signal intensity. In consequence, a dense surface generates a significant effect on LIBS measurement of powder sample, which is the reason accounting for the difference in signals obtained by the three sample preparation methods.

3.2. Crater Topography

To figure out the surface differences between the pellet samples and the adsorption band samples, a scanning electron microscope (SEM, model: TM4000) is employed to investigate the crater topography of LIBS. Figure 5 shows the SEM images of pellet sample and adsorption band sample. A definite difference in crater topography was observed from the samples surface of different pretreatment processes. SEM images of the pellet sample are illustrated in Figs. 5(a) and 5(b). SEM images of the adsorption band sample are shown in Figs. 5(c) and 5(d). In addition, the upper right corner of each figure is a photo which was taken by an optical microscope. According to Figs. 5(a) and 5(c), the surface of the pellet sample is loose and non-uniform, while the surface of the adsorption band sample is dense and uniform. At the same time, the photos of the optical microscope show that the surface of the pellet sample contains more holes when compared with the surface of the adsorption band sample. Based on parts b and d of Fig. 5, there are the laser craters of pellet sample and adsorption band sample. According to Fig. 5(b), the laser crater of the pellet sample is irregular in shape, and the edge of the laser crater is characterized by an obvious laser shock. By contrast, the laser crater of the adsorption band sample presents a relatively regular spot shape without obvious laser shock characteristic. The same characteristic of laser craters can also be observed in the optical micrographs, as shown in the right corner of Figs. 5(b) and 5(d). The reason for the difference in crater topography is considered to be the different surface states of the samples.

Fig. 5.

Images of pellet sample and adsorption band sample. (a) Fresh surface and (b) laser crater on pellt sample, (c) fresh surface and (d) laser crater on adsorption band sample. (Online version in color.)

Because continuous contact based on the electrostatic adsorption force between the powder and the substrate during the laser irradiation, the steel powders form a dense surface on the adsorption band sample. During the plasma excitation and evolution, the surface of the sample is less influenced by the shock wave. And the interaction process between laser beam and sample material is comparatively sufficient and stable. Therefore, the laser crater of the adsorption band sample is close to that of the solid sample. However, in the pellet sample, the steel powders are combined relying on the residual force of the tableting process. This combined state is unstable, which is not difficult to be destroyed by the plasma shock wave. As a result, the interaction process between the laser beam and the sample material tends to be affected by the shock wave. The laser crater shows irregular shape and is accompanied by material splash marks in the edge area. The obtained result demonstrates that the surface of the adsorption band sample is more suitable for LIBS measurement.

3.3. Calibration Curves

The calibration curves for the tableting and adsorbing sample pretreatment processes are compared in the current work. Figure 6 shows the calibration curves of LIBS measurement for steel powder in different states. The measurement for each sample was repeated 5 times. The data points are the mean values and the error bars are the standard deviations of 5 measurements. All of the data obtained from 10 samples were adopted for establishing the above calibration curves. Furthermore, Mn I 403.307 nm and Fe I 404.581 nm were selected to calculate the intensity ratio of LIBS measurement. As illustrated in Fig. 6, the determination coefficient (R2) of the calibration curve is 0.910 for the pellet samples, whereas it is 0.998 for the adsorption band samples. In this perspective, the trend of the data points is clearer in the condition of adsorption band samples. Based on Fig. 6(a), several data points deviate from the fitted calibration curve, which leads to a lower R2 value. It is considered that the LIBS measurement of pellet samples is affected by the plasma shock wave. Because the laser ablation process is accompanied by the splashing of steel powder, the plasma excitation and evolution process can fluctuate in the state of the pellet sample. Hence, the intensity ratio of LIBS measurement is also affected by the plasma shock wave. Figure 6(b) shows the calibration curve of adsorption band samples. Obviously, all the data points are near the fitted calibration curve. This result indicates that the LIBS quantitative analysis of the adsorption band samples is less affected by the plasma shock wave.

Fig. 6.

Calibration curves of LIBS measurement for steel powder. (a) Measurement with pellet samples and (b) measurement with adsorption band samples. Intensity ration selects the spectral line pair of Mn I 403.307/Fe I 404.58. (Online version in color.)

In addition, as for the error bar, the results of pellet samples are more discrete than adsorption samples. This result indicates that large irregular craters are formed on the surface of the pellet sample. There are a lot of sample splash particles around the crater, seriously affecting the detection results of LIBS spectral signals. Briefly, this phenomenon can also be well explained according to the crater topography, which shows that the degree of signal fluctuation is closely associated with laser ablation.

By comparing the slope of the calibration curve, it can be discovered that the slope of the pellet sample is greater than that of the adsorption band sample. As the content of manganese continues to increase, the trend of the slope gradually becomes saturated. Apparently, it is reported that LIBS has caused the effect of self-absorption, contributing to a decrease in the slope of the calibration curve.35,36,37) These effects include not only the resonance spectrum of manganese selected in this paper which is easily affected by the self-absorption effect, but also the matrix effect is closely related to the sample preparation process for powder samples. In consequence, the intensity and stability of the spectral signal obtained by the proposed method are weaker than others. Therefore, it is impossible to obtain a satisfactory calibration curve of the pellet sample. Meanwhile, it is not suitable for performing LIBS real-time online quantitative analysis.

3.4. Measurement Uncertainty and Precision

Using the cross-validation method to evaluate the prediction model of Mn element concentration,38) one of the sample(from S1 to S9) is taken as an unknown sample and the remaining samples are taken as calibrated samples to perform cross validation. The predictive test with S6 is taken as the tested sample and the other samples are taken as the calibrated samples, as shown in Fig. 7. Moreover, the predictive test for pellet samples is shown in Fig. 7(a) and the predictive test for adsorption band samples is shown in Fig. 7(b). R2 of the pellet samples model is 0.932, whereas R2 of the adsorption band samples model is 0.998. Thus, the data points have a better fit of adsorption band samples. To explain the observed activity, we might consider that the analysis capability of artificial steel samples is improved by employing the sample preparation method.

Fig. 7.

Predictive test with S6 as the tested sample and the other samples as the calibrated samples. (a) Measurement with pellet samples and (b) measurement with adsorption band samples. Intensity ration selects the spectral line pair of Mn I 403.307/Fe I 404.581. (Online version in color.)

To quantitatively compare measurement uncertainty and precision between tableting process and adsorbing process sample pretreatment method, the statistic factors like relative error of prediction (REP) and the relative standard deviation (RSD) were calculated.39) REP is used to evaluate the precision of measurement and RSD is applied to evaluate the uncertainty of measurement. Uncertainty is defined as the degree of agreement between the measured experimental value and the true value. Therefore, the REP was computed according to Eq. (1).40)   

%REP= 100% Y m 1 N i=0 N ( Y s - Y p ) 2 (1)
where Ym is the average reference concentration of the predicted sample, YP denotes the prediction concentration predicted using the established calibration model of the ith sample, Ys is the reference concentration of the ith sample and N represents the total number of samples. For a fine calibration model, REP should be as small as possible.

Precision describes the degree of coordination between single test data when repeated measurements. Thus, RSD was calculated using the following Equation:   

s= [ ( x i -M ) 2 n-1 ] 2 %RSD=100%×[ s M ] (2)
where xi is the data measured by LIBS, M is the average of this set of data, n is obtained by measuring parallel data and s is the standard deviation of the evaluation of the measurement precision.

Tables 3 and 4 show a comparison of RSD and REP for the pellet samples and the adsorption band samples. The REP values and RSD values of the adsorption band samples which are applied to represent uncertainty and precision are lower than the pellet samples. The smallest value of REP for the pellet sample is 0.1%, which is 0.3% for the adsorption band sample. On the contrary, the largest value of REP for them was 14.4% and 5.5% respectively. For the same reason, the smallest value of RSD for the pellet sample is 17.9%, which is 4.8% for the adsorption band sample. Meanwhile, the largest value of RSD for them was 24.8% and 8.7% respectively. Moreover, the average value of REP is reduced from 7.73% to 1.36% and the average value of RSD is decreased from 20.34% to 6.22% through adsorbing process. These results indicate that the improvement of measurement uncertainty and precision is attributed to the preparation method of the adsorption band process. Therefore, the result demonstrates that the use of adsorbing process is beneficial for better analytical results of manganese in steel powder.

Table 3. Measurement uncertainty and precision of pellet samples.
Tested sampleMeasured value (%)REP (%)RSD (%)
S10.54649.317.9
S20.90599.416.6
S31.716714.424.8
S42.14357.219.2
S52.50250.122.1
S62.79286.920.0
S73.82389.320.3
S83.84893.819.2
S94.91499.220.7

Table 4. Measurement uncertainty and precision of adsorption band samples.
Tested sampleMeasured value (%)REP (%)RSD (%)
S10.49850.36.3
S20.9960.44.8
S31.50650.45.0
S41.97091.58.1
S52.63765.56.0
S63.0080.38.7
S73.34840.45.5
S83.82880.45.1
S94.63463.06.5

4. Conclusion

To conclude, the performance of conventional and novel sample pretreatment methods of the analysis of manganese in steel powders was in this study. According to the experimental results, the following conclusions can be drawn in the following aspects.

(1) The LIBS spectra of pellet samples and adsorption band samples were compared. The intensities of Mn and Fe spectral lines are significantly enhanced by the sample pretreatment with the adsorbing process. 3–4 folds signal enhancement has been obtained by adopting the novel sample pretreatment method. Through the investigation into the topography of laser craters with SEM, the dense surface formed by the adsorbing process is contributed to the signal enhancement of LIBS measurement.

(2) The calibration curves of pellet samples and adsorption band samples, which were prepared for 11 artificial steel powders, were established based on the LIBS measurement results. The data points of the pellet samples slightly deviate from the calibration curve while the data points of adsorption band samples are located very close to the calibration curve. Although the determination coefficient (R2) for adsorption band samples is 0.998, it is only 0.910 for pellet samples. This result indicates that the LIBS measurement data about adsorption band samples remain more stable.

(3) The simulated quantitative analysis of steel powders demonstrates that the analytical performance of LIBS was improved by the adsorbing process. The mean REP value was reduced from 7.73% to 1.36%, and the mean RSD value was decreased from 20.34% to 6.22%. The pretreatment of steel powders with the adsorbing process can significantly enhance the spectral signal and facilitate the analysis of chemical elements. The mechanism of the improvement refers to that the dense surface of the adsorption band weakens the plasma shockwave effect during LIBS measurement.

Acknowledgement

The current work was supported by National Natural Science Foundation of China (No. 62005218), National Natural Science Foundation of China for key program (No. 91860206) and National Science and Technology Major Project (No. 2017-VII-0001-0094).

References
 
© 2022 The Iron and Steel Institute of Japan.

This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs license.
https://creativecommons.org/licenses/by-nc-nd/4.0/
feedback
Top