With a non-invasive and real-time technique, nanoscopic polymer motion inside a polymer filament is imaged during 3D printing. Fine-tuning this motion is crucial for producing constructs with optimal performance and appearance. This method reaches the core of plastic layer fusion, thus offering insights into optimal printing conditions and material design criteria.
In recent times, 3D printing technology has revolutionized our ability to design and produce products, but optimizing the print quality can be challenging. The process of extrusion 3D printing involves pressuring molten material through a thin nozzle and depositing it onto previously extruded material. This method relies on bonding between the consecutive layers to create a strong and visually appealing final product. This is no easy task, as many parameters, such as the nozzle temperature, layer thickness, and printing speed, must be fine-tuned to achieve optimal results. In this study, a method for visualizing the polymer dynamics during extrusion is presented, giving insight into the layer bonding process. Using laser speckle imaging, the plastic flow and fusion can be resolved non-invasively, internally, and with high spatiotemporal resolution. This measurement, which is easy to perform, provides an in-depth understanding of the underlying mechanics influencing the final print quality. This methodology was tested with a range of cooling fan speeds, and the results showed increased polymer motion with lower fan speeds and, thus, explained the poor printing quality when the cooling fan was turned off. These findings show that this methodology allows for optimizing the printing settings and understanding the material behavior. This information can be used for the development and testing of novel printing materials or advanced slicing procedures. With this approach, a deeper understanding of extrusion can be built to take 3D printing to the next level.
The method of 3D printing is an additive manufacturing technique in which an object is manufactured layer by layer to form the desired shape. This method has a large and diverse user-base thanks to its versatility, affordability, and ease of use. Fused deposition modeling features a moving extruder (with a diameter of hundreds of microns to a couple of millimeters) to deposit molten plastic into the desired shape1. The extruded plastic should behave in a liquid-like manner for a certain duration to achieve good fusion with the previously printed plastic and form a strongly cohesive material. However, the plastic should cool down and solidify quickly after printing to prevent the plastic from flowing away from the print location and reducing the print quality. This delicate interplay between heating and cooling has been shown to directly underpin the balance between the mechanical strength and the geometric accuracy of the final 3D-printed object2. To achieve the optimal heating-cooling balance, the plastic is extruded at a temperature just above its melting temperature, and a fan head, attached to the printer, is used to cool the plastic down quickly. An in-depth understanding of the effects of printing temperatures and cooling speeds could provide the insights required for developing advanced slicing and printing protocols that maximize the mechanical or geometric results in the areas where they are most important. Efforts to gain more insight into these processes often rely on infrared (IR) imaging, which only visualizes the surface temperature3,4,5 and does not indicate the inner temperature of the plastic. Local heating beyond the melting transition drastically increases the polymer mobility and, thus, allows polymer entanglement between the old and new material. This temporally enhanced polymer motion is a requirement for the formation of the final cohesive material6,7, but IR imaging can only measure polymer motion indirectly through the surface temperature8,9. Translating the surface temperature to layer bonding thus necessitates precise knowledge of the core-surface temperature gradient and the associated complex polymer dynamics over a range of time and length scales. A direct measurement of the layer bonding (i.e., the polymer entanglement process) would allow the visualization of the mechanism underlying bulk material cohesion without a priori information or assumptions.
To gain an understanding of the spatial and temporal distribution of layer bonding, an imaging technique that directly quantifies the dynamics of the polymers that make up the plastic filament is employed in this work. This technique, laser speckle imaging (LSI), relies on interferometric light scattering to visualize nanoscopic motions, independent of the chemical composition. Depending on the optical properties of the sample, it can accurately measure several millimeters to centimeters into non-transparent materials10,11,12, unlike IR imaging, which reports only surface temperatures8,9. These attributes have recently made speckle-based methods popular in understanding dynamic processes in a plethora of materials, although they were originally developed for medical applications10,11,12. Recently, LSI has been used to gain insight into the behavior of advanced polymeric materials such as self-cleaning liquid crystal polymer networks13,14, as well as for predicting fracture in rubber15 and for studying self-healing materials16.
The feasibility of applying LSI to 3D printing was showcased in a previous article17, where a portable LSI setup with real-time analysis capabilities was presented, and it was shown that the deposition of molten plastic results in increased polymer motion multiple layers below the current layer. In the paper presented here, systematic research into the effects of the cooling fan speed on the degree of multi-layer bonding is performed. An improved plug-and-play version of the portable instrument is used that can be operated by users without optics or programming expertise. The speckle images are analyzed in real-time using Fourier transforms17, which visualize the amplitude of the speckle intensity fluctuations. This instrument has an additional brightfield camera that is aligned with the speckle camera so that the LSI motion maps can be overlayed with the brightfield images for easier interpretation without the brightfield light affecting the motion maps. The experimental approach presented in this article can be used to gain more insight into the melting, layer bonding, and solidification of extruded plastic during the 3D printing of challenging geometries and materials.
1. Setup and alignment of the LSI instrument with the 3D printer
2. Preparation of the 3D-printing design and G-code
Figure 1: Object design. A 3D view (left) and 2D view (right) from the side, front, and top of the object design. The grid represents 1.0 mm x 1.0 mm, with 1.0 cm x 1.0 cm in bold. The wall is 25 mm x 12 mm x 1.2 mm (width x height x depth), and the ridges have a width of 1.0 mm, have a depth of 0.4 mm, and are separated by 1.0 mm. The windows have a width of 1.0 mm and a height of 2.0 mm. The 3D design can be found in Supplementary Coding File 1. Please click here to view a larger version of this figure.
Property/setting | Value |
Filament | Polylactic acid (PLA), white |
Nozzle diameter | 0.4 mm |
Layer thickness | 0.2 mm |
Nozzle temperature | 210 °C |
Cooling fan speed | 100% |
Printing speed | 10 mm/s |
Travel speed | 10 mm/s |
Bed Temperature | 60 °C |
Table 1: The 3D print settings. The settings and printer properties used for slicing the object design. For the second experiment, the fan speed was manually changed to 0%.
3. Performing the experiment
A simple object was drawn as a test target for the experiments: a wall with ridges on the back, two windows, and a large hole (Figure 1). The object was sliced with the printer settings and properties listed in Table 1.
The LSI instrument was aligned with the 3D printer, and the experiment was performed. The user-friendly setup features an additional brightfield camera, which helps during alignment and allows for an easy comparison between the plastic extrusion and the measured polymer motion. The speckle and the brightfield cameras are both equipped with optical filters that prevent interference from the other channel. More technical details about the setup can be found in Supplementary File 1, and an explanation of the analysis routine is presented in Supplementary File 2. Highlights of the results of this experiment are shown in Figure 2, and the full movie can be found in Supplementary Movie 1. As shown before, the experiment can be performed just as well with a home-built instrument17.
Figure 2: Time-lapse of printing with a 100% cooling fan speed. Left: Brightfield, front-view image of the object when the printer is almost finished. The quality of the print looks good upon inspection; although the surface shows the layer lines, the overall designed geometry has been produced. Right: Four LSI snapshots from the white delineated region during the printing process; the blue arrows indicate the print head position at the time of the snapshot, as the LSI images do not correspond in time with the brightfield image. The lighter colors in each snapshot indicate increased polymer motion, which is observed in the most recently printed layers. Note that the region with enhanced motion (the welding zone) is multiple layers thick. The full detailed movie of the experiment is available in Supplementary Movie 1. Please click here to view a larger version of this figure.
Complementary to these results, the print was visually inspected; as expected for these commonly used polymer filaments and print settings, the quality was good. The designed geometry was, indeed, reproduced, and the surface was even, with a small line visible on each layer. With the LSI data, it was possible to gain in-depth insight into the printing process. The freshly extruded plastic was visible as highly mobile, and the mobility decreased gradually as it cooled down. The height of the area with high mobility (i.e., the welding zone) was four to five layers thick throughout the printing procedure, indicating a well-defined duration of layer fusion.
The experiment was repeated with the cooling fan speed manually adjusted to 0%. With this setting, the plastic did not cool down fast enough, which affected the print quality. Highlights of the results are shown in Figure 3, and the full detailed movie can be found in Supplementary Movie 2.
Figure 3: Time-lapse of printing with a 0% cooling fan speed. Left: Brightfield, front-view image of the object when the printer is almost finished. The visual quality of the print looks poor; the surface shows irregular layer lines and big blobs. Additionally, the overall designed geometry has been imperfectly reproduced; notably, the windows and holes are deformed. Right: Four LSI snapshots from the white delineated region during the printing process; the blue arrows indicate the print head position at the time of the snapshot, as the LSI images do not correspond in time with the brightfield image. The lighter colors in each snapshot indicate increased polymer motion, which can be observed throughout the whole object. The full detailed movie of the experiment is available in Supplementary Movie 2. Please click here to view a larger version of this figure.
In line with expectations, visual inspection of the 3D-printed construct indeed showed poor print quality. The layers were unevenly distributed, and the designed geometry was reproduced with deformations. A comparison of the brightfield images in Figure 2 and Figure 3 shows the major effect of the cooling fan on the surface quality and shape of the print result. The origin of this effect was determined by comparing the LSI results from Figure 2 and Figure 3. With a 100% cooling fan speed, enhanced polymer motion was observed in a region only a few layers below the extruded plastic. Therefore, each layer was liquified moderately a few times to achieve layer bonding without plastic flow. With the 0% cooling fan speed, enhanced polymer motion was observed through the whole object. Thus, each layer was liquified many times and extremely close to the freshly extruded plastic, resulting in a loss of geometric accuracy through plastic flow.
To gain a more quantitative view of the cooling fan effect in more moderate situations, the cooling fan speed was systematically varied. The object design was simplified to a wall of 25 mm x 12 mm x 0.8 mm (width x height x depth) without holes or ridges. The same print settings as in Table 1 were used. The experiment was performed 12 times, with cooling fan speeds of 0%, 20%, 40%, 60%, 80%, and 100%, each in duplicate. The resulting movies can be found in Supplementary Movies 3, Supplementary Movie 4, Supplementary Movie 5, Supplementary Movie 6, Supplementary Movie 7, and Supplementary Movie 8, as well as Supplementary Coding File 6, Supplementary Coding File 7, Supplementary Coding File 8, Supplementary Coding File 9, Supplementary Coding File 10, and Supplementary Coding File 11.
To quantitatively compare the welding zones for different fan speeds, advanced data analysis was performed on the LSI results. The goal of this data analysis was to obtain a height profile of the extent of the polymer motion in the welding zone. The associated fully commented MATLAB script can be found in Supplementary Coding File 4 and is described briefly. For every LSI image in the movie, a height profile is computed by taking the mean along the horizontal direction. The profiles of the images where the print head is in the ROI show a distinct peak around the welding zone. To exclusively select those profiles, only profiles with a peak above 8 dB are taken into account. Profiles in which this peak is too close to the edge of the ROI are also discarded. The peak positions of all profiles are subsequently aligned to give an average profile relative to the height at which the polymers are most mobile. The resulting profiles for the six different cooling fan speeds are plotted in Figure 4.
Figure 4: Height profiles for the systematic variation of the cooling fan speed. Left: The welding zone profiles for cooling fan speeds of 100% (black), 80% (blue), 60% (purple), 40% (red), 20% (orange), and 0% (yellow), obtained from the advanced data analysis script in Supplementary Coding File 4. The shaded region is the standard deviation between duplicate experiments. The right schematic explains the averaging procedure to obtain the profile of a typical LSI image. By aligning the maximums of the peaks of all the obtained profiles, the welding zone is obtained. The maximum of the welding zone (relative height = 0) is the height at which the polymers are most mobile. Full detailed LSI and brightfield movies of each experiment are available in Supplementary Movie 3, Supplementary Movie 4, Supplementary Movie 5, Supplementary Movie 6, Supplementary Movie 7,and Supplementary Movie 8. The objectprinted for this figure can be found in Supplementary Coding File 5, with corresponding G-code files in Supplementary Coding File 6. Please click here to view a larger version of this figure.
The welding zone profiles for 40%-100% cooling were nearly identical. The welding zone for 20% cooling had a shoulder reaching into several deeper layers. The welding zone for 0% cooling extended over the entire measured area. The height at which the polymers were most mobile lay in or slightly below the most recently printed layer. This phenomenon explains the presence of an LSI signal at positive relative heights, as there is printed material above the mobility peak. In all cases, the welding zone reached much deeper than the 0.2 mm layer thickness.
Supplementary File 1: LSI setup.xls. Hardware parameters of the LSI instrument used here. Please click here to download this File.
Supplementary File 2: LSI analysis.docx. Explanation of the conversion of the raw speckle images into LSI images. Please click here to download this File.
Supplementary Movie 1: LSI and brightfield movie of the experiment described in Figure 2. The movie is played at 12.5x real-time speed. The upper part is the LSI result, and the lower part is the synchronized brightfield view with the LSI ROI indicated. Please click here to download this Movie.
Supplementary Movie 2: LSI and brightfield movie of the experiment described in Figure 3. The movie is played at 12.5x real-time speed. The upper part is the LSI result, and the lower part is the synchronized brightfield view with the LSI ROI indicated. Please click here to download this Movie.
Supplementary Movie 3: The 100% cooling fan speed experiment described in Figure 4. The movie is played at 12.5x real-time speed. The upper part is the LSI result, and the lower part is the synchronized brightfield view with the LSI ROI indicated. Please click here to download this Movie.
Supplementary Movie 4: The 80% cooling fan speed experiment described in Figure 4. The movie is played at 12.5x real-time speed. The upper part is the LSI result, and the lower part is the synchronized brightfield view with the LSI ROI indicated. Please click here to download this Movie.
Supplementary Movie 5: The 60% cooling fan speed experiment described in Figure 4. The movie is played at 12.5x real-time speed. The upper part is the LSI result, and the lower part is the synchronized brightfield view with the LSI ROI indicated. Please click here to download this Movie.
Supplementary Movie 6: The 40% cooling fan speed experiment described in Figure 4. The movie is played at 12.5x real-time speed. The upper part is the LSI result, and the lower part is the synchronized brightfield view with the LSI ROI indicated. Please click here to download this Movie.
Supplementary Movie 7: The 20% cooling fan speed experiment described in Figure 4. The movie is played at 12.5x real-time speed. The upper part is the LSI result, and the lower part is the synchronized brightfield view with the LSI ROI indicated. Please click here to download this Movie.
Supplementary Movie 8: The 0% cooling fan speed experiment described in Figure 4. The movie is played at 12.5x real-time speed. The upper part is the LSI result, and the lower part is the synchronized brightfield view with the LSI ROI indicated. Please click here to download this Movie.
Supplementary Coding File 1: wall_with_holes.stl. The 3D design for the object described in Figure 1. Please click here to download this File.
Supplementary Coding File 2: wall_with_holes.gcode. The sliced object wall_with_holes.stl with the settings from Table 1. Please click here to download this File.
Supplementary Coding File 3: config.ini. The configuration file for the slicing software. Please click here to download this File.
Supplementary Coding File 4: AdvancedDataAnalysis_FanSpeed.m. The script to perform the advanced data analysis on the cooling fan sweep data and plot Figure 4. The script is fully commented. Please click here to download this File.
Supplementary Coding File 5: wall.stl. The 3D design of the object used to collect the data in Figure 4. Please click here to download this File.
Supplementary Coding File 6: wall_100%fan.gcode. The sliced object wall.stl with a 100% cooling fan speed. Please click here to download this File.
Supplementary Coding File 7: wall_80%fan.gcode. The sliced object wall.stl with an 80% cooling fan speed. Please click here to download this File.
Supplementary Coding File 8: wall_60%fan.gcode. The sliced object wall.stl with a 60% cooling fan speed. Please click here to download this File.
Supplementary Coding File 9: wall_40%fan.gcode. The sliced object wall.stl with a 40% cooling fan speed. Please click here to download this File.
Supplementary Coding File 10: wall_20%fan.gcode. The sliced object wall.stl with a 20% cooling fan speed. Please click here to download this File.
Supplementary Coding File 11: wall_0%fan.gcode. The sliced object wall.stl with a 0% cooling fan speed. Please click here to download this File.
The experiments and results described in this research show that LSI is an easily applied tool that allows for a deeper understanding of layer bonding during additive manufacturing. LSI allows for the direct measurement of the polymer motion, which must be finely tuned to form a coherent material by the interpenetration and subsequent entangling of the polymer chains. The most common alternative for measuring layer bonding in situ is infrared imaging3,4,5. This well-established method images the local surface temperature of the plastic8,9, which is an indirect measurement of the polymer motion inside the material. With hotter plastic, the motion is faster, and the bonding becomes stronger. However, the relation between temperature and motion is not linear, as the printing temperatures cross the melting and glass transition temperatures6,7. This non-trivial relationship can be directly observed in the LSI images; specifically, there is a sharp transition between the liquid-like top and the solid-like bottom regions, while the temperature gradient is expected to be far more gradual. Another drawback of IR imaging is that it measures the surface temperature only, while LSI measures polymer motion typically several millimeters deep inside the material.
Just like with IR imaging, this implementation of LSI is essentially a point-and-shoot method; it can be used in situ if the camera can be pointed at the region of interest. The versatile tripod and long working distance of 0.7 m give the freedom to use any available 3D printer. Crucially, LSI is sensitive to nanoscopic motions, and, thus, vibrations from the surroundings and the printing process itself have to be minimized17. For example, performing another task on the same table or slamming a door will cause interference. Therefore, one should carefully walk around the setup; however, room lights or airflow do not generally interfere with the process.
LSI gives detailed insight into the layer bonding process and can be applied as easily as IR imaging. We envision that LSI has great potential in aiding the development and understanding of advanced 3D printing methods. The cooling fan speed sweep shows a glimpse of what is possible by combining LSI with 3D printing. As discussed in the introduction, the optimal cooling speed is a balance between keeping the plastic molten for long enough to improve layer bonding but cooling it down fast enough to prevent flow. The 40%-100% cooling fan speed results were very similar; indeed, these fan speeds did not show any flow and produced a good surface quality. With the 0% cooling fan speed, the material started to flow away from the printed location, but ample layer bonding was observed in the LSI measurement. Based on our results, the 20% cooling fan speed could be optimal for achieving slightly improved layer bonding without compromising the surface quality. However, to draw conclusions that can be applied in practice, more cooling fan speeds between 0% and 40% must be assessed. It is also desirable to establish quantitative measures for the surface quality and material strength to gain an objective and complete view of the effects of polymer motion on the desired properties. With this addition, the approach could be made more powerful for evaluating creative 3D printing advances.
The exact settings that are chosen for the LSI analysis are not prone to critical errors as long as liquid-like plastic and solid-like plastic phases can be distinguished clearly. The polymer motion changes dramatically when crossing the melting and glass transition temperatures, so a broad range of LSI settings capture the contrast well. This can be easily tested with a test print of a straightforward object (e.g., a straight wall) with the 3D printer settings that are recommended by the material supplier. For more advanced LSI users, delving deeper into the frequency range can yield extra information, as different types of polymer motion can be quantitatively distinguished. For example, high-frequency polymer motion is associated with the highest temperatures, which are only present close to the printer head. Lower-frequency polymer motion is associated with moderate temperatures, which are present in a much larger area around the printer head and also for a much longer time17. Whether the degree of bonding for cumulative low-frequency polymer motion could be equal to that with short, high-frequency motion (e.g., with dynamic mechanical analysis) must be examined. Most other settings, like the colormap scaling, ROI, saving interval, and experiment length, are solely chosen to give a visually clear and appealing result. Regarding the 3D printing settings, there is also much freedom, as LSI allows the user to objectively assess the results of changing any of the settings. Notably, drastically changing the printing speed changes the interpretation of the LSI data. In this work, a slow printing and travel speed of 10 mm/s were used in order to capture multiple LSI images during one pass of the printer head. If a more common print speed of 60 mm/s for PLA were used, roughly one full layer would be printed per LSI image, and, thus, averaging within one layer would occur. If experimenting with high-end speeds such as 300 mm/s and faster, averaging over multiple layers would occur. Nevertheless, this is completely dependent on the exact print geometry and LSI settings and could be easily mitigated by an experienced LSI user through advanced machine design, adjusting the size of the field of view, or using a faster camera. Both approaches necessitate a more powerful laser, which, in combination with the reflective printer head, requires additional laser safety precautions. The relatively slow printing speed also has a positive influence on the layer bonding, as it was previously proven that heat transfer to the plastic increases with slower printing speeds5.
One possible new direction for this approach is the testing of novel materials; for example, LSI could be used to visualize the relevant transitions and objectively quantify the recommended printer settings that give a five-layer welding zone upon the application of the top layer. Another application could be to study the welding zone in specific situations where the print quality is not reliably good, such as for bridges, overhangs, or sharp corners. If the welding zone in difficult situations can be better understood, it should be possible to compensate in the G-code. It is already common practice to print the first layer hotter and slower than the rest of the layers to achieve good adhesion to the build plate18. We envision the use of similar dynamic G-code slicing where, for example, the fan cooling could be adjusted to produce corners or bridges. It should also be possible to print the outer wall material with a smoother finish and the rest of the material and infill rougher but stronger to maximize both the material strength and visual appearance.
This article has discussed the application of LSI to study the layer bonding process after plastic extrusion. The technique is excellent for this task, as it can visualize the underlying polymer motion without a priori assumptions in real-time during 3D printing. However, it does not give any information on the material cohesion, so additional testing will be required. The other drawbacks discussed are situational; the limited imaging speed of four LSI images per second can be increased with a bigger laser and additional laser safety measures, and the vibration sensitivity requires precautions or vibration reduction hardware. LSI can be performed with cheap and small digital cameras and lasers19,20, which allows for integration into virtually every 3D printer for live quality control and dynamic tuning of the printing parameters. However, it makes more sense to employ LSI to develop thorough knowledge of layer bonding during 3D printing. If this understanding is used to develop more advanced slicing software, every consumer 3D printer could benefit from the knowledge gained.
The authors have nothing to disclose.
The authors received no external funding.
3D-drawing sofware | Autodesk | TinkerCad | tinkercad.com |
3D-Printer | Prusa3D | Original Prusa i3 MK3S | |
Advanced data analysis software | MathWorks | MATLAB R2018b | |
Image viewing sofware | National Institutes of Health | ImageJ 1.47v | |
LSI instrument | NanoMoI | NanoMoi allround | company to be founded 2023 |
Polylactic acid (PLA) filament | REAL | filament white 1,75 mm PLA 1 kg | |
Slicing software | Prusa3D | PrusaSlicer-2.5.0 |