1. Introduction
Some treatments of damaged tissue require a scaffold that can mimic aspects of the missing tissue in an in vivo environment to support healing. Materials taken from the human body have superior biocompatibility compared to synthetic material [
1,
2]. Furthermore, they integrate into the body—eliminating the need for a second removal surgery [
3]. Therefore, the donated tissues taken from a cadaver (known as allografts) have been a preferred material for regenerating injured tissues such as bone [
2]. However, allografts are available in limited geometries. A method of fabricating allograft structures of arbitrary geometry from donated tissues while retaining favorable features would allow the fabrication of customized geometry to meet patient requirements.
Additive manufacturing (AM) has the capability of fabricating arbitrary geometries, but typical AM materials have significant limitations as scaffolds [
4]. An AM method that utilized donated tissue, while maintaining key biological functionality, could provide significantly improved treatments, but the process may alter the biological features of the material during fabrication. In this study, we investigated the adaptation of laser sintering (LS) to fuse a demineralized bone matrix (DBM)/polymer composite. Laser sintering (LS) is an AM technique that uses laser to selectively sinter the particles based on the 2D slices of a CAD model [
5]. Additional 2D layers are processed sequentially to build up a 3D physical object.
Good mechanical properties, excellent accuracy, and high build rates have made LS valuable in the industry. Expanding the range of materials that can be sintered is key to further growth [
6]. In this regard, one of the most promising directions is the development of composite material systems [
7,
8,
9]. This is of particular value in biomaterials where mechanical properties, composition, and geometry have a strong impact on functionality [
10]. The pores should have the characteristics to facilitate cell ingrowth in the actual body microenvironment. Work in this area was reviewed by Loh et al. [
11].
Efforts to apply laser sintering for tissue engineering applications are actively being explored and diverse methods have been demonstrated to create porous scaffolds with an interconnected 3D flow channel network [
11,
12]. As an example, calcium phosphate was processed by an intermediate polymer which was thermally removed during post-processing [
13]. Tan et al. [
14] conducted a study on producing a biocomposite of polyetheretherketone (PEEK) and hydroxyapatite (HA). They found that increased bed temperature resolved a delamination problem and enhanced the layer integrity. The degree of interparticle boding could be controlled by adjusting LS process parameters.
LS fabrication with DBM particulates has not been reported previously. However, a composite of DBM and polycaprolactone (PCL) and the mixture has been fabricated by fused deposition modeling (FDM) [
15] and large area projection sintering (LAPS) [
16]. Hung et al. fabricated DBM/PCL scaffolds containing up to 70 wt % DBM using FDM [
17]. The composite samples had greater cell proliferation compared to control samples of pure PCL. However, this approach embeds much of the DBM in PCL, reducing its biological impact and potentially altering its geometry. Alternatively, a mixture of DBM and PCL powders could be selectively fused using a lower PCL content to increase porosity and exposure of DBM surfaces to improve biological performance if the mechanical properties are sufficient. Previous research showed promising mechanical properties by using a projected image to fuse the DBM/PCL blend. However, the geometric accuracy was insufficient due to the poor absorption of the projected wavelength [
1].
Although the advantages of using a CO
2 laser to process DBM particles were described earlier, the focused beam could pose the risk of overheating this valuable graft donation. This paper investigates the feasibility of using a CO
2 laser to process DBM/PCL composites to achieve basic mechanical functionality. The biological viability of the process is not assessed here. DBM particles impose a strict limitation on the process conditions due to their sensitivity to temperature. Processing conditions must be carefully controlled to melt the PCL while not overheating the DBM (T
max ~70 °C). There are over 50 variables that play a role in laser-based additive manufacturing that fall under four major groups: (1) laser parameters and scan strategies; (2) powder properties; (3) spread layer quality; and (4) processing environment conditions [
18,
19,
20]. This paper will focus on laser parameters and scanning strategies.
2. Experimental Methods and Materials
2.1. Materials and Layer Formation Study
Laser sintering requires a thermal bonding agent. Typically, all or part of the raw material melts and densifies under surface energy forces when heated by the laser. Since DBM does not melt on heating, a secondary thermal bonding agent must be added to the mixture. Potential biopolymers were screened and polycaprolactone (PCL) was selected as the matrix for the DBM powder due to its low melting point and biocompatibility. For this project, the PCL was purchased from Polysciences Inc. The as-received powder was sieved and the portion <300 μm was used in this research. The Discovery from TA instruments was utilized to take a differential scanning calorimeter (DSC) measurement of the PCL, which showed that this PCL has a melting point of Tm ~63 °C, below the processing limits of the DBM [
21]. The molecular weight of the PCL was reported 50,000 g/mol and the shape of the PCL was irregular. An initial test was performed on pure PCL to make sure that the laser effectively interacts with PCL. Where the CO
2 laser beam traces the surface, local bonds form between PCL particles, as seen in
Figure 1. The PCL is a feasible bonding agent for the DBM/PCL mixture.
Donated cadaveric bone was processed at LifeLink Tissue Bank, where a proprietary clinical cleaning process was used to remove the blood, lipids, viable cells, and microorganisms. The cortical bone was further processed by grinding to create a coarse powder. The resultant DBM powder was sieved and the granules smaller than 250 μm were utilized. The processing and mechanical properties are sensitive to the fraction of PCL and DBM used in the blend. Prior work showed that a 45 wt % DBM and 55 wt % PCL mixture results in good mechanical properties (comparable to pure DBM) while maximizing DBM content [
1]. The same mixture was used in this study. Scanning was done with a Rabbit HX-SE-1290 system with a CO
2 laser tube that operates at 80 watts nominal power. According to the laser apparatus specifications, the beam spot size in the focal distance of 51 mm is 0.2 mm. In order to minimize the local heating of the samples, relatively low laser power, and slow scan speeds were studied. Generally, at such low powers, laser power values may not be accurate, but the results of this study indicate key tradeoffs and reflect the potential for using a low-cost laser cutter for additive fabrication.
The interaction between the laser beam and the material system was studied by fabricating single layer primitives that were 1 cm × 1 cm with a range of speed, laser power, and line spacing. For consistency, throughout this research, all experiments were performed by a bidirectional raster pattern (see
Figure 2) at room temperature (no bed preheat). The laser parameter setting started with a low laser power at which no bonding initiated. For each set of scanning conditions, the laser power was gradually raised until interparticle bonding emerged, marking the lower bound of processing power. This incrementing continued by 0.008 or 0.08 watts until the color change of DBM material occurred, identifying an upper bound of the processing power range.
This comprehensive study identified the parameter window over which the materials could be processed without visible damage. Successful processing required that the powder mixture fused together without curling out of plane. The layer cohesion must be at least sufficient to be handled without damage. It was observed that the layers have insufficient strength to survive handling directly after fusing. To enhance the cohesiveness of the layers, two different methods were examined. In the first approach, the line spacing was reduced to stitch the neighboring traces to each other more effectively, and in the second, the distance between the bed surface and the laser exit aperture was increased to enlarge the laser spot size, delivering energy to a larger area. The laser spot size estimation for the new laser-to-powder distance was performed by knowing the fact that, based on laser specifications and at 51 mm focal distance, the spot size was 0.2 mm. This spot size was found to be the same as the individual laser trace line widths on plywood at a speed of 150 mm/s and p = 5.36 W. With the same parameters, the spot size for 57 mm distance was estimated as 800 μm. For each case, the weight of the primitive layers was considered as an indicator of thickness and integrity of the layer.
Multilayer parts were fabricated by adding a layer spreading apparatus into the laser cutting system. A layering set-up was designed and built considering the space limitations inside the work envelope of the laser cutter (
Figure 3). To form a layer, the powder supply piston raises the material; at the same time, the build piston goes down to accommodate one new layer. Then, the counter-rotating roller with a rotational speed of 120 rpm and 50 mm/s traverse speed recoats the material over the build platform, and retracts to its initial place. The laser scanning is manually triggered to fuse one 2D layer. By repeating this process and putting down successive layers, the 3D geometry was created.
2.2. Sample Preparation
The spreading test results performed in our previous study [
1] showed that the mixture of DBM <250 μm and PCL <300 μm powders can spread uniform layers of 500 μm thickness. Therefore, the layer thickness was set at 500 μm for fabricating test coupons. In this study, two types of samples were fabricated: tensile bars for dimensional change and tensile strength assessment, and 5 mm cubes for density and porosimetry analysis. All these samples, regardless of the geometry, were printed by laser power = 5.36 W, line spacing = 100 μm, and scan speed = 150 mm/s.
In total, 9 tensile bars composed of 45 wt % DBM-PCL with a 500 μm layer thickness were printed. The CAD dimensions of the gauge section were in accordance with ASTM D638-10 standard with a 13 mm width and 6 mm thickness. The overall length of the samples was 64 mm. The tensile bar samples were printed in the x–y plane with the long direction in the rolling direction. The as-printed tensile bars showed sufficient mechanical properties for handling but were weaker than pure DBM strips. Thus, isothermal heat treatments were carried out in a convection oven at 65 °C, which is slightly above the measured melting point of the PCL. Each set of three samples was considered as a group to undergo different heat treatment durations: 3 h, 5 h, and 7 h. Three DBM strips were machined to have a 13 mm gauge width to be consistent with the printed samples for strength comparison [
1].
Figure 4 shows an example of composite tensile bars before and after 5 h heat treatment and the reference sample which was machined from DBM. The cube samples which were for density and porosimetry studies were heat-treated for 5 h.
The 5 mm cubes were manufactured using two DBM weight fractions: 45 wt% DBM and 55 wt% DBM. This allows identifying the role of DBM content in density and porosity of the scaffold. The printed composite cubes heat-treated at 65 °C for 5 h in post-processing and to have a reference sample, the DBM strips were was cut into the same dimension for the density and pore characterization experiments. The skeletal density was measured by using two different techniques for comparison: the helium pycnometer (UltraPyc 1200e) and mercury intrusion porosimetry (MIP), model PoreMaster 60, both from Anton Paar. Besides density, the MIP technique is capable of analyzing porosity. The first set of measurements were taken by the helium pycnometer, then the same samples were utilized in the MIP instrument to minimize the error that would have been caused by using new samples. The mercury penetration porosimeter correlates the penetrated volume of mercury and the penetration pressure for each sample. The pressure increases from 3.5 Kpa to 4.1 × 105 Kpa and the porosity is determined herewith.
3. Results and Discussion
In this section, we describe how the key process parameters impact the layer quality and DBM discoloration. The results include the first series of experiments that were conducted on the layer primitives to find the best combination of process settings to use for multi-layer parts. Tensile bars were fabricated to investigate heat treatment impact on dimensional deviation and mechanical strength. In addition to the tensile bars, 5 mm cubes were manufactured to study skeletal density and void volume fraction in the final parts.
3.1. Process Parameters Effect on DBM/PCL Primitive Layers
The layers are composed of the adjacent lines cured by the laser path. Initially, the stepover was set at 100 μm for which the lines laterally met each other and formed a layer. However, it was expected that reducing stepover distance by half should stitch the lines to each other more effectively, but
Figure 5 shows that the variation of stepover changes the operation window significantly at the higher scan speeds. The processing window for the 100 μm line spacing expands for higher vector speeds. However, layers made with a line spacing of 50 μm show little change in the size of the processing window with scan speed in the ranges tested. Delivering excessive energy into the material may burn the DBM in the powder bed. This is more critical at the end of the laser vector where the beam turns back to trace a fresh line in which case the previously delivered energy has not sufficiently been dissipated and overheating is more likely to happen. Running the process with small line gaps would necessitate using a tightly controlled laser power with minimal fluctuation.
The second parameter examined was the beam spot size. The spot sizes were increased from 200 μm to 800 μm by changing the laser-to-powder bed distance from 51 mm to 57 mm. The experiments were conducted for a combination of power and speeds (see
Figure 6). The expanded beam spot delivers energy into a wider trace of material and keeps the PCL softened for a longer period of time yielding better cohesion in the layer. The attenuation of the energy intensity can be compensated by operating at a higher laser power. The feasible range of powers increased significantly, but the size of the print features is limited to beam size (800 μm).
3.2. Printing Defects
The low-cost laser system is subject to errors in power output and power fluctuations during operation [
22]. Wider processing windows are preferable for more stable and consistent outcomes. In this study, we observed that the power fluctuations near the edges of the power window could drastically degrade layer quality (see
Figure 7). When the parts are sintered near the lower bound of the processing window, laser power variation may fail to cure some sections. The resultant layer is formed incompletely, which is a critical problem for layer integrity. Similarly, when the power is close to the upper bound, power fluctuations may overheat the DBM particles. These defects emerged at slower scan speeds (50 mm/s and 100 mm/s) with the smallest processing window. To minimize the risk of bone discoloration, a higher speed of scan (here 150 mm/s or 200 mm/s) and operating power with some margins from the bonding and burning thresholds are preferred for LS of the DBM/PCL material.
The edges of the layers also curl upward under some process conditions. The raised regions can cause part shifting and other problems when spreading a fresh layer of material. As seen in
Figure 8, layers processed at higher power levels have a significantly larger mass. These parts have greater handling strength, but the parts with higher mass curl more. For a fixed laser power and line spacing, the scanning speed is also found to impact both layer mass and layer flatness (see
Figure 9). The curling of the edges occurs when the top surface of the powder bed experiences higher temperatures than the depth and temperature gradient curls the layer upward as the top surface shrinks faster than the bottom [
23]. This phenomenon can also be caused by premature crystallization of polymers. This limits the use of higher laser power to improve layer integrity [
24].
Looking at both figures reveals that by increasing the laser power, more material is cured in depth of powder bed. This increases the layer weight, but after a threshold causes warpage. So, the best process setting for manufacturing DBM/PCL scaffolds at ambient temperature cannot be determined solely based upon samples’ weight (thickness) because severe curling occurs for excessively cohesive layers. The process setting has to meet a combination of DBM discoloration requirements, acceptable cohesion, plus minimal warpage in the layers.
3.3. Optimized Processing Parameters
By incorporating the observations of the single layer primitive study, it is concluded that for processing at room temperature and with no preheated powder bed, the energy input must be maintained below the threshold that causes layer distortion. Among the processed layers, the following set of variables created the most favorable quality of the layer primitive: power = 5.36 W, line spacing = 100 μm, and scan speed = 150 mm/s. This set of variables falls in the safe operational window to process DBM particulates and was utilized to manufacture all tensile bars and cubes for further characterization. Unfortunately, these process parameters produce limited strength. Therefore, the initial fusion of the material to create the scaffold geometry was followed by an isothermal heat treatment to strengthen the parts as described below.
3.4. Dimensional Accuracy
The heat treatment of the samples in post-processing ensures that the PCL has melted throughout the part and improves the consolidation level in the samples. But, the heating shrinks the parts. To evaluate shrinkage, the sample height, gauge section, and the overall length of the tensile bars were measured before and after post-processing [
25]. The shrinkage of the nine tensile bars that were printed was assessed directly after laser fusion is summarized in
Table 1.
The length of the tensile bar was very close to the target of 64 mm. The scanning was performed lengthwise along the tensile bar for scanning layers. The lateral heat conduction causes some material to stick to the edges of the print and make planar green dimensions turn out slightly greater than the CAD model. It should be noted that since the beam width is constant during printing, the amount of the neighboring material that sticks to the layer width and length is expected to be approximately similar. This dimensional increment results in a greater percentage error in the narrower part (gauge) of the sample. However, scan paths could be offset to compensate for this deviation. In addition, the beam movement inaccuracy, optical aberration, and beam scattering in the powder bed could be other sources of error [
26]. The shrinkage in the thickness is highest as porosity between layers may be higher and gravity promotes the rearrangements in this direction [
27].
The dimensional change of the samples from the desired CAD values is presented in
Figure 10 as a function of the curing time. The samples undergo minimal shrinkage during the 3 h heat treatment. Samples cured for 5 and 7 h underwent comparable shrinkage though the variation between samples was reduced at the longer treatment. The longer treatment would be preferred because the initial model can be compensated to offset for the mean shrinkage; however, the variation becomes the error. The length has the smallest shrinkage. This is probably because the continuous scan of the laser along the sample length creates fully fused lines that resist shrinkage during post-processing. In these samples, every layer is scanned in the same direction so that this impact is clearly visible. If the scanning direction were varied by layer, the shrinkage would likely be reduced, but the part may distort during post-processing.
3.5. Tensile Strength Assessment
The heat treatment caused anisotropic shrinkage that varied with the heating time and is expected to influence the mechanical properties of the cured samples. The tensile strength of the samples was measured using an MTS 858 universal hydraulic testing system. In order to acquire accurate data, the grip section of the samples must not deform during testing. Thus, a low viscosity epoxy was infiltrated into the pores in the grip sections of the tensile bars and hardened overnight at room temperature. The induced tensile forces were measured with a 900 N load-cell at a grip displacement rate of 1.5 mm/min. An exemplary set of stress–strain curves for a 7 h heat-treated sample and a machined DBM strip is presented in
Figure 11.
Looking at the tensile test results (
Figure 12) reveals that maximum tensile strength is linked to the curing time duration. The samples cured for 3 h have the least dimensional change and retain their weak mechanical property. The 5 h heat treatment increases tensile strength significantly as the molten PCL has more time to migrate into the DBM particles and increase adhesion. Longer heat treatment (7 h) slightly improves the tensile strength, which suggests that maximum tensile strength reaches a plateau beyond which further heating does not increase shrinkage. The tensile strength of the cancellous bone falls between that of the DBM/PCL composite samples with heat treatments of 3 h and 5 h. This suggests that the post-processing or the PCL composition could be tailored to match the mechanical properties of the natural DBM if desired.
Table 2 compares Young’s modulus (E) for DBM strips and the LS-fabricated parts with different post-processing times. It can be seen that the demineralized bone strips have the least Young’s modulus. For the produced samples, stiffness is tied to heat treatment duration. Stiffness increases significantly with increased curing time from 0.22 MPa after 3 h to 4.333 MPa after 7 h heat treatment. This could be due to the fact that longer heat treatments allow for more viscous sintering evolution or because the polymer chains have more time for crystallization, leading to a reduction in the extensibility of the polymer [
28].
3.6. Density and Porosity Characterization
The density and pore channels in scaffolds are of paramount importance to regenerate the bone. Load distribution and proper biological functionality are closely linked to the void fraction and architecture of the pore network in the bone. These properties typically vary with location in the body [
29]. The pore volume of bone and the scaffold samples were measured using helium pycnometry [
30] and the skeletal density was calculated. Using these values, the density of PCL was estimated. Composite samples were fabricated at two DBM mass fractions (45%, 55%) to assess the composition sensitivity.
Figure 13 demonstrates that the DBM fraction is proportionally related to the skeletal density of the final product. This confirms that there is minimal closed porosity in the manufactured samples that would reduce the measured density.
The MIP measurements, shown in
Figure 14, indicate that, for the 45 wt % DBM and 55 wt % DBM blends, there is a similar trend as well as pore distribution through the scaffolds. It is also evident that the maximum pore size in both samples is ~250 μm.
The density measurements and porosity summarized in
Table 3 provide additional insight into the part structures. The skeletal density of the DBM/PCL composite cubes is comparable for pycnometry and MIP. The MIP additionally shows that there is significant open porosity (>40%) to support bone ingrowth. This porosity increases for higher DBM fractions. The DBM samples were not successfully measured using MIP due to the presence of very large pores and the possible collapse of the structure under pressure.
It is evident that skeletal densities from gas pycnometry for 45 wt % DBM and 55 wt % DBM samples are is quite similar to what measured by MIP. The slight difference could result from the limitation of mercury to protrude the holes smaller than 3 nm, which might exist within the DBM particulates. On the other hand, helium can intrude almost all the pores thanks to the small molecular size. Another noteworthy point is that the higher content of DBM yields more porosity in the construct which may be the result of different DBM and PCL arrangements in the print and/or fine pores of the DBM particles.
The results indicate that CO2 laser is able to create a porous construct with DBM/PCL material and that the final mechanical properties, as well as porosity, highly depend on DBM fraction and post-processing conditions. These results could be extended by adjusting the composition ratios, using modified material systems such as adding biocompatible porogens or using more controlled powder morphology to adjust the balance between porosity and mechanical properties. In vitro and in vivo testing of the resultant parts to evaluate the biological features are preserved is another open question that awaits further study.