Classification of metal PBF-LB parts manufactured with different process parameters using resonant ultrasound spectroscopy

To face the challenges raised by the qualification of metallic additively manufactured (AM) complex shaped and rough finish parts, non-destructive testing (NDT) volumetric methods are required. X-ray computed tomography (XCT) is presently the favored technique; however, alternative methods are needed to overcome the requirement of technical skills and the high cost of the technique. XCT also has limitations regarding the size and density of parts. Here, we propose an easy to use, fast, and efficient global NDT volumetric method based on resonant ultrasound spectroscopy (RUS) which basic principle relies on the comparative analysis of natural resonant frequency spectra of similar parts from the same family, both of which vibrating as free as possible. The methods have already proven to have the ability to sort parts with defects from flawless parts. In the present study, we demonstrate that RUS can also segregate metallic parts manufactured with different AM system process parameters. Eleven sets of three parts were manufactured, using a metal laser-powder bed fusion process, with different wall thicknesses, laser powers, scanning speeds, and scanning strategies. These parts were tested by RUS and then analyzed using the Z-score statistical method. The AM process parameter changes clearly influenced the resonance responses of the parts, and thus, the method is able to classify the different groups of parts according to their process parameters. Hence, the RUS methods can provide industries convenient tools to not only identify defective parts but to also configure AM machine parameters according to the expected and desired material properties.


Introduction
The layer wise process to build parts in additive manufacturing (AM), powder bed fusion (PBF), material jetting (MJT), vat photopolymerisation (VPP), and binder jetting (BJT) categories of processes provides several advantages over conventional manufacturing technologies, which includes the possibility to fabricate more complex geometries, for instance, lattice or internal structures. Nevertheless, in metal AM, the as-built AM surface finish can be very rough, and like all other manufacturing methodologies, it cannot systematically produce parts that are imperfection free and compliant with required specifications. To check the integrity of these parts, before use in their intended applications, some inspection methodologies, preferably non-destructive, are required. However, the inspection of AM parts raises numerous challenges for non-destructive testing (NDT) which is used commonly for quality assurance of parts manufactured conventionally. Indeed, some of the imperfections are specific to AM processes [1], to which should be added to the complex geometry and the roughness of the parts. The requirements of the inspection methods are that they need to be non-destructive, but also volumetric to inspect internal structures, not limited by surface roughness and/or complex shapes, and they also need to be fast and inexpensive if possible. Among the existing volumetric NDT methods, the most powerful and accepted in terms of inspection capability of complex shapes and rough surface finish part is X-ray computed tomography (XCT) [2][3][4][5]. Generally, it has enough spatial resolution to locate and even evaluate the dimensions of most internal imperfections in AM parts [6]. However, size and density of parts can be a major obstacle to the inspection. The combination of a part's density and thickness can strain the practical limitations of the XCT source power required for X-rays to fully penetrate the part. Additionally, the XCT chamber needs to be large enough to contain the entire part. XCT is also a relatively high-cost, time-consuming inspection technique that requires skilled and certified technicians to perform the testing. In addition, the resulting scan files are very large, and thus, they are not easily transferred, stored, or viewable. These disadvantages do not allow all parts to be scanned systematically. Hence, alternative non-destructive volumetric methods are required. Among all of the available non-destructive technologies, resonant ultrasound spectroscopy (RUS) [7][8][9][10] gives the most valuable results-after XCT-to control metallic parts [8,11,12]. RUS inspections are global or whole-body methods able to sort parts based on the analysis of their natural resonant frequencies by comparison to the natural resonant frequencies of parts from the same family or to the modelled natural resonant frequencies. In addition to the sorting, the measurement of material elasticity constants (Young's modulus E, shear modulus G, Poisson's ratio ν, and anisotropy characteristics) of parts whose geometry is easily measurable such as bars, cylinders, tubes, and discs is possible. RUS is directly tied to the structural integrity of the parts, easy to use, fast, and exponentially affordable compared to XCT. They also can inspect any part size and shape and are insensitive to the inherent surface roughness, and they can detect material property variations that XCT cannot.
RUS is commonly used on simple geometrical specimens for defect detection or to measure the material elasticity constants [13][14][15][16][17][18][19]. A.S. WU et al. [20] performed a RUS parametric study into the effects of laser scanning pattern, power, speed, and build direction in powder bed fusion (PBF) AM on residual stress. Ibrahim et al. [12] used RUS to determine material properties and quality of lattice structures through the measurements of the modulus of elasticity. In our previous studies, we have evaluated RUS methods on complex AM parts containing different numbers of geometries added in their numerical design (i.e., in computed-aided design (CAD)) to simulate flaws, designated by ISO-TC261/ASTM-F42 as CAD seeded flaws, and manufactured with laserpowder bed fusion (PBF-LB) process [8,9]. We have also investigated RUS methods on PBF-LB lattice structures with different numbers of missing struts [10,11]. In both studies, the methods have given valuable results. The methods have been proven to be successful in the sorting of complex parts with CAD seeded flaws and/or missing struts from identical reference parts, supposedly flawless, but does the RUS method have the capability to classify parts according to  High laser power for contour 0.5 2 53 900 2 53 900 SCh2 SCh3 S1C1 One contour instead of two 0.5 2 53 900 1 80 900 S1C2 S1C3 manipulation of AM process parameters, representing as closely as possible, the potential naturally occurring flaws? In the previous study [8][9][10][11], the tested parts with CAD seeded flaws or missing struts, manufactured using PBF-LB, were compared to the reference parts, whereas in this study, we propose to classify different PBF-LB parts manufactured with different process parameters.
The paper is structured as follows: the first section describes the principle of RUS; the second section is dedicated to the presentation of the investigated metal PBF-LB parts manufactured with different process parameters; the third section reports on the experimental RUS tests and data statistical analysis of these parts. Finally, the last section presents the validation of the RUS results by comparison with experimental RUS tests with another RUS system and by XCT inspections. The RUS methods are whole-body comparative inspection approaches. Their basic principle relies on the excitation of the part under test to make it vibrate as free as possible to record the natural resonant frequency spectrum of its vibrational modes and then to compare this spectrum with those of parts from the same family [21] or to modelled natural resonant frequencies. According to ASTM E2001-18 [22], there are two types of RUS methods: "swept sine method" and "impulse excitation method (IEM)." This distinction is based on the way the physical excitation of the parts under test is generated. In both cases, it is a mechanical impulse of the parts, but, in the first case, with contact piezoelectric transducers, whereas in the second case, the excitation is generated with a mechanical impulse. Then, the frequency spectrum is monitored with contact piezoelectric transducers in the case of a "swept sine  SCh1 SCh2 SCh3 SH1 SH2 SH3 S1C1 S1C2 S1C3 All Parts σ(Hz) 13 18 11  19  16  16  21  19 24 σ mean (Hz) 17 method," while it is recorded with either a microphone or a piezoelectric transducer in the case of IEM methods. Two similar objects will have comparable resonant frequency spectra. A shift of a frequency or multiple frequencies of a tested part, compared to other parts from the same family, will be the indicator of a structural difference in the tested part: changes in its geometry, density, elasticity, or external/internal structural integrity (e.g., material alterations or imperfections and AM process deviations).

Additively manufactured (AM) investigated parts with different process parameters
To evaluate the influence of AM machine process parameters on the properties of additively manufactured parts, cylindershaped parts (Fig. 1), with different thicknesses, were fabricated implementing various laser powers, scanning speeds, and scanning strategies. These parts were then measured with a coordinate measuring machine (CMM) to assess geometrical  Fig. 1a.
Several cylinders were then manufactured (Fig. 1b), in Inconel 625, with a 3DSystems ProX DMP320 PBF-LB machine. The scanning strategy, for all cylinders, consisted of back and forth trajectories with a rotation of 65° between each layer and a contour step occurring after the filling step. The hatch space was set by default to be 100 µm and the layer thickness to be 60 µm. Supports were added below the bases to remove the parts easily from the build platform. The difference between the cylinders, in wall thicknesses, laser powers, scanning speeds, and scanning strategies, are synthesized in Table 1. Eleven sets of process parameters or wall thickness combinations were chosen, and three parts per combination were manufactured to check the AM process repeatability.
A set of parts with default machine parameters recommended by the PBF-LB machine manufacturer (Std) was produced. Then, the first combination of three sets was manufactured with various laser powers (P in Watt) and scanning speeds (v in mm/s) (Pb, Ph, and Pm), defined by their volumetric energy density (VED in J/mm 3 ): where S represents the surface of the laser beam which was measured with a beam profiling camera (in mm 2 ).
A second combination of three sets with three different thicknesses (Db, Dm, and Dh) was manufactured, and finally, a third combination of four sets with various scanning strategies (SH, SCb, SCh, and S1C) as specified in Table 1. A stripe scanning strategy and two contour steps were chosen for all cylinders except for the set referenced SH, for which a hatching scanning strategy was chosen, and for the set referenced S1C, there is only one contour step. The difference between hatching and stripe scanning strategies is illustrated in Fig. 2.

SH3
On the edge SH3 and SH1 or SH2 S1C3 On the edge S1C3 and S1C1 or S1C2 Fig. 9 RUS system from GrindoSonic used to test the investigated parts The parts were manufactured within the same batch with other parts. They were localized as illustrated in Fig. 3.

RUS tests
In this study, LNE tested the parts with its IEM system (RAM-UP11HK) developed by The Modal Shop (TMS) (Fig. 4a). The system is equipped with a connected modal hammer, enabling the impact force to be measured, and a microphone. The output of the microphone is converted into a frequency spectrum in the range from 100 Hz to 94 kHz using a high-speed analog to digital converter (24 bits) performing a fast Fourier transform (FFT). For these tests, to excite the cylinders, the parts were presented for impulse from a slide to the modal hammer such as shown in Fig. 4b. Each cylinder was tested and positioned similarly on the slide (reference at the bottom, cylinder perpendicular to the slide, launched from about halfway up the slide), in the frequency range of 100 Hz to 94 kHz, with the highest resolution of 1.465 Hz which gives an acquisition time of 0.683 s.
Each cylinder was tested ten separate times to allow for verification of repeatability. The overall averaged spectra for the parts with default machine parameters recommended by the manufacturer (Std1-3) is presented on Fig. 5. Table 2 gives the square root of the variances' average (σ) for the frequencies corresponding to the maximum amplitudes of the seven first peaks, for each part, as well as the square root of the average of these σ (σ mean ) to quantify the repeatability for separate tests and all tests.

RUS data statistical analysis
The data statistical analysis was performed by Weaver NDT using Z-score method to classify the different groups of parts within the overall parts by comparison to the set of parts with default process parameters recommended by the manufacturer (Std1-3). For each part, ten RUS tests were performed. The analysis implies the averages (avg) of the ten frequencies of the maxima for the seven first peaks in the whole frequency spectra. This selection is because lower frequency data are usually key to material property differences. This choice turned out to be relevant to classify all the cylinders as this is presented below.
The Z-score [23] was, thus, calculated from the below relation: The Z-score deviation was then plotted versus the Z-score average for all parts as presented in Figs. 6 and 7.
These plots are particularly evocative. The AM process parameter changes clearly affect the frequency position of the resonance peaks, and one can classify the different set of three parts according to their process parameters. As shown in Fig. 6, the key parameter is wall thickness, and then, as shown in Fig. 7, high laser power variation for the contour is very important, more important than the one contour instead of two. A low VED has more impact than a high VED, and if the VED is the same as the default parameters, but with higher laser power and increased scanning speed, the impact is quite important. The type of scanning strategies (hatch or stripe) has a lower impact on the integrity of the parts. Such plots clearly show how the AM process parameter changes influence the resonance responses of the parts.
Among the Pm, SH, and S1C series, dispersion within each set is observed as highlighted in Fig. 8. Pm1, SH3, and S1C3 are outliers (Table 3). This may be because these parts are all located on the edge of the AM platform (see Fig. 3 parts framed by a red square). To identify the origins of the differences in these responses, they were

Experimental RUS tests with another RUS system
To demonstrate that the resonant frequency spectra of the tested parts are not system or operator dependent, the parts were also tested with an IEM RUS system developed and commercialized by GrindoSonic, namely, the MK7 (Fig. 9). Both methods, from TMS and Grindo-Sonic, include a modal hammer and a microphone. However, the excitation of the part was different in the two cases. LNE presented the parts to the TMS system from a slide to the hammer, whereas GrindoSonic hit the parts, in the middle of the cylinder height, laying on a foam, with the hammer such as show in Fig. 9. Two tests were performed per part as the standard deviation over ten separate tests on a part was evaluated to be 0.093 Hz on the 8.5 kHz peak. Figure 10 shows the comparison of a RUS spectrum acquired by GrindoSonic and a RUS spectrum acquired by LNE on a part with default machine parameters recommended by the manufacturer (Std1) all over the frequency range. Most peaks below 50 kHz are located at the same frequency. However, LNE did not observe any peak around 25.4 kHz whereas GrindoSonic does. We attributed this dissimilarity to the difference in excitation between LNE and GrindoSonic, but this is still invalidated. The mode is possibly not activated by the free-falling process used by LNE. Above 50 kHz, the comparison was not possible as the range of sensitivity of the systems are different. Figure 11 presents a comparison, performed by GrindoSonic, of the frequency positions of the two first vibrational modes for all parts. This figure confirms the statistical analysis performed by Weaver NDT on the LNE data. The highest shifts from the Std set are observed for the Dm, Db, and Dh sets, then for Pb, SCh, and SCb. A difference is also observed between Pm1, SH3, and S1C3 and the rest of their respective set.

X-ray computed tomography (XCT) inspections
In order to find an explanation of the dispersion of RUS results within a set of similar parts with identical process parameters, the outlier as well as another part from the same set was scanned by Baker Hughes using their Waygate Technologies v|tome|x m 300/180 XCT system (Fig. 12). This 300 kV system is equipped with a micro focus source and a dyn|41-100 detector with a 100 µm pitch distance. The distance between the source and the detector (FDD, focal detector distance) is around 800 mm. The XCT parameters used to scan the parts of interest are given in Table 4, and Fig. 13 presents an XCT image of the whole cylinder.
The XCT unrolled images for Pm1/Pm2, SH3/SH1, and S1C3/S1C2 for three different layers on the cylinder element, in the center (defined as position 0 mm), and then on both borders (defined as position -0.15 and + 0.15 mm from the center) are presented in Table 5. This table also gives the porosity rate (in %) for each layer.
It can be observed, either on the images or on the porosity average and/or the standard deviation, that the porosity rate or the distribution is significantly different for all outliers compared to the part from the same family. A higher standard deviation indicates that the porosity is not evenly distributed among the parts, while a higher average indicates that the porosity rate is higher. Thus, the differences between the parts, produced with the same parameters, can be assigned to a difference in porosity levels and/or the uneven distribution of the porosity and therefore a difference in overall material modulus.

Conclusions
We have conducted an experimental and statistical study to investigate the feasibility of resonance ultrasound spectroscopy (RUS) to classify complex shaped and rough finished additively manufactured (AM) parts, according to their machine process parameters.
In the first section, the principle of RUS has been reported. In the second section, the AM parts have been detailed. Eleven sets of three similar cylinders with a base have been manufactured, in Inconel 625, using a laser-powder bed fusion process. Either the wall thicknesses, the laser powers, the scanning speeds, or the scanning strategies were modified in between these eleven sets. In the third section, the RUS experimental conditions have been described, the experimental results on the cylinders are provided, and their statistical analysis has been presented. The results are particularly conclusive. The method could successfully classify the different parts according to their AM process parameters. In the fourth section, the RUS results are validated using another system and operator, and the XCT inspection is documented. The other RUS system gave similar results. One single response difference was identified, and this was determined to be dependent on the method by which the part is excited. The XCT scans highlighted that the outlier within a same parameter series is attributed to location on the edge of the AM build platform and porosity distribution.
This study demonstrated that the IEM RUS methods provide qualitative results in identifying defective parts, but they also are successful in classifying parts according to their AM process parameters. The methods provide valuable benefits. They are non-destructive, volumetric, easy to use, very quick, and simple to operate, and no part preparation or fixturing is required. In addition, there is no restriction in size and shape for the test parts, and surface roughness does not prevent accurate and repeatable testing. This study has shown that RUS methods are very suitable tools for the AM industry, as an alternative to XCT, for routine quality inspections for medium to high volumes.