An in-line rheometer has been incorporated into a fused deposition modeling printer for the first time by designing a modified nozzle with a custom pressure transducer and a thermocouple for measuring the processed melt temperature. Additionally, volumetric flow rates and shear rates were monitored by counting the stepper motor pulses as well as the pulses from a custom filament encoder to account for filament slippage and skipped motor steps. The incorporation of the sensors and the design and development of the in-line rheometer are described; and pressures, temperatures, and viscosities within the 3D printing nozzle are presented. The in-line rheometer was validated against traditional, off-line rotational rheology and capillary rheology measurements by analyzing two polymeric materials: polycarbonate and high-impact polystyrene. A variety of rheological corrections were considered for the in-line rheometer, including entrance effects, non-Newtonian corrections, shear heating, pressure effects, and temperature fluctuations/inaccuracies. Excellent agreement was obtained between the in-line and off-line rheometers after applying the most critical corrections, which were found to be entrance effects, non-Newtonian corrections, and temperature inaccuracies. After applying the appropriate corrections, the in-line rheometer provides an accurate viscosity measurement that can be used for real-time monitoring and process control.

Fused deposition modeling (FDM®), also known as fused filament fabrication (FFF), is a 3D printing technique that deposits strands of a molten polymer that solidify into the final, desired shape as described in previous works [1–4]. FDM is an appealing process with low capital investment that allows users to quickly and easily create complex designs unattainable by traditional manufacturing processes. Despite the benefits and widespread use of FDM, there are many opportunities for improvement, including part strength [1–9], resolution [10–15], and material availability [16–18].1,2 Many of these problems could be improved with advances to FDM process control. Mani et al. from the National Institute of Standards and Technology (NIST) call out the need for more real-time measurements and in-line analysis of additive manufacturing to develop process understanding and predictive models that can be used for real-time quality control [19]. Turner and Gold propose that regulating flow rate through process monitoring and control algorithms is one of the greatest opportunities for improving FDM resolution [10].

A first step to implementing the control algorithms proposed within the NIST paper and article by Turner and Gold is the development of appropriate on-line sensors; Huang et al. similarly state that one of the largest gaps in the field of additive manufacturing is the lack of sensors capable of providing relevant, real-time data [20]. The only sensors and feedback controls in most FDM printers are proportional-integral-derivative temperature controllers for the thermistors of the nozzle and build platform. The motion, positioning, and extrusion are controlled with open-loop stepper motors in over 95% of entry-level FFF printers [21]. While open-loop stepper motors offer many benefits, they are prone to skipping steps where the lack of feedback control results in inaccurate positioning and speed [22]. Weiss et al., therefore, included an encoder on the X and Y gantries to monitor the trajectory and accuracy of the positional motors, they improved the accuracy by using a PI controller, and they recommend incorporating more closed-loop motion control into FDM printers [21].

To combat some of these control and resolution problems, high-end FDM companies experimentally optimize their print conditions for their filament materials by applying an empirical program during the acceleration phases of a print [23]. While effective for a limited set of materials, this experimental optimization makes it very difficult and time-consuming to introduce new materials. To combat resolution problems due to dimensional inaccuracies of the filament, the high-end printer companies have described many forms of area measurement sensors that could provide feedback to the drive system to correct for feed-rate errors [24]. Greeff and Schilling took flow rate control one step further by not only measuring filament diameter but also monitoring the extrusion speed and the motor gear speed through visualization and image analysis to detect motor errors and slippage in addition to filament inaccuracies; they incorporated a control algorithm to decrease the deposition speed when slippage occurred [15].

New in-line sensors, like the ones described by Greeff and Schilling, are critical for understanding the processing mechanisms of FDM and confirming theoretical, predictive models. Sbriglia et al. incorporated sensors into the FDM process by stopping the printer in the middle of a job to add accelerometers, which inform users about the frequency response of the system [25]. Bragg sensors and thermocouples were added using a similar approach to monitor strain and temperature during the FDM process [26], but such methods of adding sensors into a part cannot be used for real-time control. The most common type of sensor studied in the additive manufacturing literature has been the nondestructive infrared (IR) camera for monitoring material temperature evolution during processing. The heating and cooling of the powder bed has been studied for selective laser sintering [19,27], and the initial extrudate temperature and the subsequent cooling process have been studied for FDM [8,28–30]. IR measurements have become popular in research because the transient cooling process of the FDM extrudate has been shown to be a critical parameter for predicting the strength of FDM parts [1,5,8,9,31–33]. Similarly, Phan et al. used an IR camera to compare against predicted melt temperatures within the nozzle. By utilizing off-line rheological measurements as well as the power supplied to the stepper motor that drives the filament feeding gears, they predicted melt temperatures that showed similar trends to those measured by the IR camera [34].

The original research that is the focus of this journal article adds to the development of in-line sensors for FDM by incorporating a filament feed-rate encoder, a variety of thermocouples, and a custom pressure transducer to develop an in-line rheometer. Rheology is the most critical material property for most polymer melt and liquid processing methods, including extrusion, injection molding, coating, and FDM. Rheology influences the temperature of the melt through viscous dissipation [35–38], determines the material output and flow rate [35–37,39], determines the shear rate and velocity profiles along the material flow path [35,39,40], creates the pressure profiles that determine the force and torque requirements of the process [10,15,36,37,39,41], and can influence the final material shape and strength [9,32,42–44].

Due to the practical and theoretical importance of rheology, on-line and in-line rheometers have been developed and extensively studied for many polymer processing techniques, including extrusion [45–47], injection molding [48], and coating processing [49]. Dealy describes on-line rheometers as taking measurements on melt that has been redirected away from the processing melt stream while in-line rheometers, such as the in-line FDM rheometer proposed in this article, offer the benefit of measuring viscosity directly in the process line with minimal delay in response [50]. On-line (or in-line) rheometers are typically validated by applying the appropriate rheological corrections and comparing the on-line measurements against off-line rheological measurements (a similar approach is also used in this article to verify the FDM in-line rheometer). For example, Köpplmayr et al. have developed methods for simultaneous on-line shear viscosity and on-line extensional viscosity on an extruder. They compared the on-line measurements to rotational and capillary rheology after applying Rabinowitsch corrections to the on-line and capillary data, and they found the on-line extensional viscosity compared well with off-line Sentmanat extensional rheometer (SER) data [47]. Coates et al. incorporated in-line pressure transducers to obtain accurate viscosity measurements for injection molding, which they recommended using, to implement closed-loop controls [48].

Similar but more complex rheological corrections have been applied to better understand the flow through capillary rheometers. Moshe et al. analyzed the transient behavior of capillary rheometers by accounting for pressure measurement errors, end effects, viscous dissipation, pressure effects, and compressibility [51]. Ansari et al. account for pressure effects, viscous dissipation, compressibility, and wall slip to model the pressure drop through the coni-cylindrical section of a capillary rheometer for low density polyethylene [52]. Different rheological corrections are critical for obtaining accurate viscosity measurements, depending on the system of interest. Despite all the on-line rheological studies for other polymer processes, no on-line or in-line rheometer has previously been developed for FDM printing.

Huang et al. state the importance of understanding the relationship between viscosity, in-line measurements, and additive manufacturing processes [20]. Similarly, a variety of other sources stress the importance of rheology to the FDM process [5,6,10,40,41,53–55]; nonetheless, there has not been any systematic study comparing material rheology to FDM outputs, such as mechanical properties, accuracy, shrinkage, roughness, or resolution. A variety of works have calculated pressure predictions through the FDM nozzle for a given viscosity profile [40,41,56–58]. Gilmer et al. model the relationship between the viscosity profile and FDM pressure and relate these to the practical problems of filament buckling and melt backflow [56]. Ramanath et al. have performed simulations to predict temperature and pressure profiles during the FDM process [40], but their analysis only used a power-law model, neglected viscous dissipation, and was never validated against real, in-line measurements. Despite the criticality of rheology, in-line rheometers have yet to be incorporated or studied on FDM to confirm the theoretical calculations or to study the influence of rheology on final properties.

To the authors' knowledge, this is the first work that describes the development of an in-line rheometer for the FDM process. This article describes the design of the custom pressure transducer and custom nozzle required for the in-line rheometer. The performance of the rheometer is validated against off-line rheological measurements as well as with an in-line comparison to a capillary rheometer. Finally, a variety of rheological corrections are considered and discussed, including end effects, non-Newtonian flow, viscous dissipation, pressure effects, and temperature corrections, which are described in Sec. III. While viscous dissipation was analyzed, it was found to be negligible, so viscous dissipation corrections were not applied to the data and are only discussed in the supplementary material [74]. Although a similar work attempted to measure the melt pressure based on the power sent to the feeding motor of an FDM printer [34], this approach is limited due to the nonlinear relationships between power, torque, and speed [22]; the approach was unable to be utilized as an accurate rheometer, so the original work presented in this article has developed a direct measurement of melt pressure, which is able to provide very accurate viscosity measurements.

An FDM printer (Lulzbot Taz 6 v 2.1, Loveland, CO) with a single extruder is used in this research. The printer is provided with a 0.5 mm diameter, brass nozzle as shown in Fig. 1. This printer model can achieve a maximum nozzle temperature of 298 °C, a maximum bed temperature of 120 °C, a maximum print speed of 200 mm/s (12 000 mm/min), a minimum layer height of 0.05 mm, and a maximum layer height of 0.5 mm. The printer is designed to process 3 mm diameter filaments, and its maximum print area is 280 × 280 × 250 mm.

FIG. 1.

Dimensions (mm) of the original nozzle.

FIG. 1.

Dimensions (mm) of the original nozzle.

Close modal

A new nozzle system with multiple custom parts was produced to create the in-line rheometer on the printer hot end. First, a new nozzle was designed [Fig. 2(a)] with a pressure port hole for measuring the melt pressure. A load transfer column [Fig. 2(b)] was designed to be inserted into the pressure port within the nozzle [Fig. 2(c)] so that it is in direct contact with the melt. The load transfer column contains a thermocouple for measuring the melt temperature [Figs. 2(c) and 2(d)], and it transfers the load from the melt to a custom load transfer clamp, which is shown in Fig. 2(e). The clamp transfers the melt load to a 10-lb load cell (voltage from the load cell is then converted into melt pressure) as shown in the assembly of Fig. 2(f). The load cell is rigidly fixed to the nozzle system using a 1.5 in. long 6–32 screw as shown in Fig. 2(f).

FIG. 2.

In-line rheometer nozzle design: (a) custom nozzle, (b) load transfer column, [(c) and (d)] load transfer column plus thermocouple inserted into the nozzle pressure port, (e) custom clamps for the load cell, and (f) full assembly.

FIG. 2.

In-line rheometer nozzle design: (a) custom nozzle, (b) load transfer column, [(c) and (d)] load transfer column plus thermocouple inserted into the nozzle pressure port, (e) custom clamps for the load cell, and (f) full assembly.

Close modal

The nozzles and load cell clamps were all produced from brass castings of lost wax patterns produced by a 3D printing process (Shapeways, New York, NY). The nozzle was designed to match the original printer nozzle as closely as possible in terms of both dimensions (Fig. 1) and materials (brass). To achieve the necessary accuracy of the critical dimensions [shown in Figs. 2(a), (c), and (d)], the custom nozzles were finished by drilling, reaming, and tapping with a 3-axis mill (EZ Trak DX, Bridgeport, Elmira, NY). The 10-lb load cell was an LCL-010 from Omega (Norwalk, CT), and all thermocouples used in this work were 2 m long, type K, glass braid insulated, 30 gauge thermocouples from Omega (Norwalk, CT). The melt thermocouple was attached to the load transfer column using a silver solder (SSQ-6, Muggy Weld, Olympia, WA) as shown in Figs. 2(c) and 2(d). Additional information about the design and assembly of all the in-line rheometer parts is provided in the supplementary material [74].

The final dimensions of one of the custom nozzles are shown in Fig. 3. It is critical to accurately measure the nozzle dimensions highlighted in Fig. 3 because the length and diameter of the nozzle tip are used to calculate the shear rate and shear stress [Eqs. (1) and (2)], and the area of the pressure port hole is used to convert the load cell force to a pressure (described in Sec. II E). The diameters of the nozzle tip and pressure port were initially measured via optical microscopy while the nozzle length was obtained by measuring the length of calibration pins inserted into the nozzle. These three dimensions were verified using x-ray computed tomography (Phoenix v|tome|x s, General Electric, Boston, MA) as shown in Fig. 3, which confirmed the manual measurements as well as the diameter uniformity of the bores. X-ray scans were run at 210 kV and 180 mA, creating images with an 8 μm voxel size.

FIG. 3.

Nozzle dimensions measured by x-ray computed tomography.

FIG. 3.

Nozzle dimensions measured by x-ray computed tomography.

Close modal

Once designed and connected, the load cell voltage was collected in real time by transferring the voltage through a signal conditioner to a data acquisition device (DAQ), which transfers the data to a personal computer. Additionally, to measure the filament extrusion rate, the DAQ was counting the pulses from both the stepper motor that drives the filament feeding gears and an encoder in direct contact with the filament. Finally, four thermocouples (including the melt thermocouple shown in Fig. 2) were measuring temperature data, which were collected through a temperature input module. The system of controllers, sensors, conditioners, and data acquisition devices is sketched in Fig. 4. The DAQ is an NI USB-6210 from National Instruments (Austin, TX) with eight differential, 16-bit analog input channels, four digital input channels, four digital output channels, and two counters. The DAQ collected voltage data from the load cell at its fastest sampling rate of 250 000 Hz.

FIG. 4.

System of sensors, devices, and connections for in-line rheology and process monitoring.

FIG. 4.

System of sensors, devices, and connections for in-line rheology and process monitoring.

Close modal

The temperature input module is an NI USB-9162 from National Instruments and converts the thermocouple voltage into temperature. An intrusive thermocouple monitors the apparent polymer melt temperature (Fig. 2), which is of critical importance as later discussed. Another thermocouple is in contact with the load cell and secured between the load transfer clamp and the load cell clamp; it is used to apply an in-line temperature correction to the load cell data as described in Sec. II E. The other two thermocouples do not directly influence the in-line rheometer, but they are important for monitoring the 3D printing process. The third thermocouple is suspended in the environment, parallel to the nozzle but located about 2 in. away from the nozzle to monitor the environmental temperature during a print. The last thermocouple is in contact with the top of the build plate as shown in Fig. 4 to monitor and confirm the accuracy of the build plate temperature.

The printer's filament drive is a NEMA 17 stepper motor (type MS17HD6P4150-1, Moons' Industries, Shanghai, China); it steps in increments of 830 pulses/mm of filament. A pulse is sent to the DAQ counter from the FDM microcontroller every time the motor is commanded to step, providing a measure for the commanded extrusion rate, but not necessarily the actual extrusion rate since the filament may slip in the drive during high pressure applications. Therefore, a filament displacement encoder was mounted to the printer to verify the volumetric flow rate (and subsequently, the shear rate) and account for skipped motor steps or filament slippage. The filament displacement encoder was provided by Mosaic Manufacturing (Toronto, Canada). An encoder housing and mounting system was 3D printed and attached to the printer frame so that the encoder is in constant contact with the 3 mm filament as depicted in Fig. 4. The filament encoder provides 7.53 pulses/mm of filament, which is not high enough resolution for high-speed measurements but successfully provides confirmation of the flow rate and shear rate through a moving average filter. It was experimentally determined that the encoder measures shear rates within ±5% of those measured by the stepper pulses when neither the motor is skipping steps nor the filament is slipping. The accuracies of the volumetric flow rates calculated from both the stepper motor and the encoder were verified by measuring the mass of printed material as described in the supplementary material [74]. The testing showed that both stepper motor and encoder were accurate with roughly the same tolerances as the cross-sectional area of the acquired filaments (about 4%). The consistency of the volumetric flow rate including variations in filament diameter is one of the more significant sources of error in both viscosity characterization and FDM process control as described in the supplementary material [74].

Matlab R2015b (Mathworks, Cambridge, MA) receives all the data from both the DAQ and temperature input module, then processes and saves the data. The data are collected using the Data Acquisition Toolbox and are analyzed using the Signal Processing Toolbox as well as the Statistics and Machine Learning Toolboxes. Temperature values are collected for all four thermocouples from the temperature input module every 0.5 s. Data are collected from the DAQ buffer every 4 ms, during which the DAQ has collected 1000 signals from each channel. The 1000 voltage values from the load cell were averaged together to provide a single voltage measurement every 4 ms, which was found to be an effective filtering method known as oversampling [22]. A shear rate is calculated from both the stepper motor and encoder counters as described in Sec. III. A low pass Butterworth filter was applied to the stepper motor shear rate data to reduce high frequency noise. The filter was a first-order Butterworth filter with a normalized cutoff frequency (cutoff frequency divided by sampling rate) of 0.01, which was tuned to prevent overshoot while maintaining acceptable time delays. Additional filtering is applied to set the shear rate to zero when it falls below a critical value of 10 s−1 (the printer was observed to have difficulty in maintaining shear rates below 20 s−1).

The load cell outputs an analog voltage from 0 to 10 V, so a calibration curve was created to convert voltage into pressure. A calibration fixture was cast in brass; the design is the same as the load transfer clamp shown in Fig. 2(e), but the 0.5 mm recess is replaced with a 3.1 mm diameter column with a platform with two holes at the top. A mass is hung from the calibration fixture via a nylon string threaded through the two holes as shown in Fig. 5 such that the nylon does not contact either side of the load cell. The mass was converted to a pressure using the area of the pressure port hole (dimensions shown in Fig. 3). The calibration was carried out by hanging the following masses from the calibration fixture: 0, 0.541, 1.127, 1.668, and 2.287 kg. The procedure with all five masses was repeated four times at four temperatures ranging from 21 to 100 °C. Temperature was recorded with the thermocouple which was clamped against the load cell, and the 20 data points were used to create a least-squares fit for the calibration curve for pressure as a function of voltage and temperature. The coefficient of determination (R2) was greater than 0.99, and an example of a pressure versus voltage calibration curve is presented in Fig. 8.

FIG. 5.

Voltage-to-force calibration for the FDM® Rheometer strain gage.

FIG. 5.

Voltage-to-force calibration for the FDM® Rheometer strain gage.

Close modal

Both polycarbonate (PC) and high-impact polystyrene (HIPS) filaments were analyzed on the in-line rheometer. The filaments were PC Transparent from Ultimaker (Cambridge, MA) and HIPS Natural from eSun (Shenzhen, China), which were both chosen based on their reported diameter tolerances: 2.85 ± 0.05 mm for PC and 2.95 ± 0.05 mm for HIPS. Any diameter deviations will result in a shear rate error proportional to the squared diameter error as indicated by Eqs. (1) and (2), so these filaments were chosen because they had some of the best tolerances on the market. The glass transition temperatures of the PC and HIPS filaments are 112 and 98 °C, respectively, as measured by differential scanning calorimetry (Q2000, TA Instruments, New Castle, DE) with a 20 °C/min temperature ramp rate.

The validation of the in-line FDM rheometer (FDMRheo) was performed by printing into open space, meaning that the nozzle was stationary and raised about 30 cm above the build platform to avoid any exit pressure effects that could build up between the nozzle and the platform. Nine shear rates were measured at three temperatures for each material as shown in Table I to account for shear rate and temperature dependencies. Each shear rate sweep was repeated twice, and results from both tests are reported in this article. In order to achieve the melt temperatures listed in Table I, the nozzle temperatures were set to higher temperatures as shown in Tables SIV–SVII of the supplementary material [74] (for example, the nozzle temperature of the printer was set to 259 °C to achieve a polymer melt temperature of 250 °C within the nozzle); the nozzle settings were determined using the intrusive melt thermocouple. The measurements were repeated with two nozzle lengths (both with the same nozzle diameter and conical entrance) to account for end effects as described in Sec. III C. The corrections were all applied to the results from the 1.72 mm length nozzle because this nozzle more closely resembles the standard printer nozzle.

TABLE I.

Materials and conditions tested on the in-line rheometer.

ConditionMaterialNozzle length (mm)Melt temperature (°C)Shear rate (s−1)
PC 1.72 250 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 1.72 270 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 1.72 290 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 3.74 250 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 3.74 270 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 3.74 290 30, 46, 65, 100, 140, 210, 300, 460, 650 
HIPS 1.72 225 30, 46, 65, 100, 140, 210, 300, 460, 650 
HIPS 1.72 250 30, 46, 65, 100, 140, 210, 300, 460, 650 
HIPS 1.72 275 30, 46, 65, 100, 140, 210, 300, 460, 650 
10 HIPS 3.74 225 30, 46, 65, 100, 140, 210, 300, 460, 650 
11 HIPS 3.74 250 30, 46, 65, 100, 140, 210, 300, 460, 650 
12 HIPS 3.74 275 30, 46, 65, 100, 140, 210, 300, 460, 650 
ConditionMaterialNozzle length (mm)Melt temperature (°C)Shear rate (s−1)
PC 1.72 250 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 1.72 270 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 1.72 290 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 3.74 250 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 3.74 270 30, 46, 65, 100, 140, 210, 300, 460, 650 
PC 3.74 290 30, 46, 65, 100, 140, 210, 300, 460, 650 
HIPS 1.72 225 30, 46, 65, 100, 140, 210, 300, 460, 650 
HIPS 1.72 250 30, 46, 65, 100, 140, 210, 300, 460, 650 
HIPS 1.72 275 30, 46, 65, 100, 140, 210, 300, 460, 650 
10 HIPS 3.74 225 30, 46, 65, 100, 140, 210, 300, 460, 650 
11 HIPS 3.74 250 30, 46, 65, 100, 140, 210, 300, 460, 650 
12 HIPS 3.74 275 30, 46, 65, 100, 140, 210, 300, 460, 650 

To allow enough time for equilibration, each shear rate was applied for approximately 25 s before stepping up to the next shear rate. A 5 s delay was applied between each shear rate to observe the full pressure equilibration process as shown in the example in Fig. 6 for all shear rates from 30 to 650 s−1 for HIPS at 225 °C.

FIG. 6.

Pressure and shear rate curves for HIPS at 225 °C.

FIG. 6.

Pressure and shear rate curves for HIPS at 225 °C.

Close modal
FIG. 7.

CapRheo/FDMRheo joint setup for verifying FDMRheo pressure measurements against a CapRheo pressure transducer.

FIG. 7.

CapRheo/FDMRheo joint setup for verifying FDMRheo pressure measurements against a CapRheo pressure transducer.

Close modal

Both rotational rheology (RotRheo) and capillary rheology (CapRheo) were performed for both PC and HIPS at the same temperatures described in Table I to validate the performance of the FDMRheo. A DHR1 rotational rheometer (TA Instruments, New Castle, DE) and a Rosand RH10 capillary rheometer (Malvern Panalytical, Westborough, MA) were used. Frequency sweeps from 460 to 1 rad/s with three points per decade on a logarithmic scale were performed on the RotRheo. The Cox–Merz rule was assumed to apply [59]. A strain of 4% was used for PC, and a strain of 2% was used for HIPS to test within the linear viscoelastic region as determined by strain sweeps. Due to the high frequencies explored in this study, an inertia correction was applied to minimize the effects of inertia, since the raw phase angle was above 100° for most of the data above 100 rad/s. The accuracy of the correction was validated by testing a 200 Pa · s Newtonian polybutene viscosity standard (Cannon Instrument Company, State College, PA), which was confirmed to maintain its Newtonian viscosity up to 460 rad/s.

Shear rate sweeps were done with three points per decade on a logarithmic scale for both materials on the CapRheo. PC was characterized with a 20 mm length die with a 1 mm diameter from 100 to 1000 s−1 at all three PC temperatures listed in Table I as well as with a 20 mm length die with a 2 mm diameter from 10 to 460 s−1 at 270 °C to observe the Newtonian plateau. HIPS was characterized with a 20 mm length die with a 2 mm diameter from 6.5 to 1000 s−1 at all three HIPS temperatures listed in Table I. The Rabinowitsch and pressure effect corrections (see Secs. III D and III E) were applied to the CapRheo results, but the Bagley correction was found to be negligible for the CapRheo experiments. All measurements on both RotRheo and CapRheo were repeated twice, and both specimens are reported in this article.

The sensor data must be converted into an apparent shear rate and a shear stress at the wall to calculate the apparent viscosity. First, a volumetric flow rate, Q, is calculated based on the stepper motor pulses (due to the much higher resolution, the stepper counts were used by default rather than the filament encoder counts unless there was a shear rate discrepancy >5%, indicating skipped steps or filament slip) as described in the following equation:

Q=V˙×A=vζ×πRfil2,
(1)

where V˙ is the velocity of the filament, A is the area of the filament, v is the measure of pulses per second, ζ is the pulses per mm (830 pulses/mm for the stepper motor), and Rfil is the radius of the filament.

The uncorrected apparent shear rate, Γ, is then calculated by the following equation [35,39]:

Γ=4QπRnoz3,
(2)

where Rnoz is the radius of the nozzle tip. Note that Γ is corrected for non-Newtonian effects into the corrected shear rate, γ˙, as later described in Sec. III D. Shear stress at the wall, τw, is calculated based on the following equation [35]:

τw=ΔP×Rnoz2Lnoz,
(3)

where ΔP is the pressure drop measured by the load cell and Lnoz is the length of the nozzle tip. Apparent viscosity, ηapp, is then calculated using an average Γ and an average ΔP for each shear rate step as shown in Eq. (4). Since the pressure and shear rate take time to equilibrate, the data for each shear rate step were averaged after the first, transient 25% of each step, so that only the equilibrated data were considered.

ηapp=τwΓ.
(4)

A rheological model was fitted to all the rheology data collected by each rheometer for both PC and HIPS to interpolate and extrapolate the viscosity, η, to other temperatures and shear rates. Therefore, the Cross-Williams–Landel–Ferry (WLF) model was chosen, which combines the shear rate dependent Cross model as shown in Eq. (5) and the temperature dependent WLF model for predicting the zero-shear viscosity, η0, as shown in Eq. (6) [39,51].

η(γ˙,T)=η01+η0γ˙τ1n,
(5)
η0(T)=D1expA1(TT)A2+(TT),
(6)

where τ is the critical shear stress at which the viscosity profile moves from Newtonian to shear thinning, n is the power-law index (which is fit to experimental data), T* is a reference temperature (the glass transition temperature is used in this article), and A1, A2, and D1 are fit to experimental data. In addition to comparing the Cross-WLF curves produced by the various rheometers (Fig. 10), the Cross-WLF models are also used to correct for pressure effects and temperature inaccuracies as described in Sec. III F.

To account for end effects (end effects include entrance effects from the coni-cylindrical taper prior to the nozzle tip as well as exit effects), the Bagley correction [35,60] was applied by measuring the viscosity with the two nozzles listed in Table I: the standard nozzle has an L/R of 6.49 while the longer nozzle has an L/R of 14.11. The entrance pressures were calculated by extrapolating the pressure to a zero L/R. The entrance pressures were then subtracted from the raw pressure data of the standard nozzle to calculate the corrected pressure. Entrance pressure is a function of temperature and shear rate, so the Bagley correction was applied across all the measurements as shown in Table I. Entrance pressure was found to be very significant, ranging from 5% to 35% of the pressure drop for PC and 8–45% of the pressure for HIPS, making entrance pressure the most critical of all the rheological corrections explored in this work.

The shear thinning nature of the polymers causes the cross-sectional velocity profile to vary from the parabolic flow pattern of Newtonian liquids that is assumed in Eq. (2). To account for this, the Rabinowitsch correction was applied according to the following equation [58]:

γ˙=Γ34+14dlnΓdlnτw.
(7)

Equation (7) adjusts the cross-sectional flow field by considering a power-law model [8,53,58]. The viscosity is then corrected by using the following equation:

η=τwγ˙.
(8)

The high pressures seen in many plastics processing techniques cause a reduction in the free volume of the melt [61–63], causing an exponential increase in viscosity at elevated pressures [35,51,61,64]. Viscosity pressure dependency can be accounted for by applying a pressure shift factor, αp, using the Barus equation as shown in the following equation [52,65]:

αp=exp(β×p),
(9)

where β is the pressure coefficient and p is the pressure in the FDMRheo nozzle. The pressure shift factor can then be applied as shown in Eqs.(10) and (11) to obtain the zero-pressure viscosity, ηp0, and the corrected shear rate, γ˙p0, which provide the corrected values to compare off-line rheological results obtained at zero pressures.

ηp0=ηαp,
(10)
γ˙p0=αpγ˙.
(11)

The β values were taken to be typical values from the literature: β=0.03MPa for PC [64,65] and β=0.034MPa for acrylonitrile butadiene styrene (ABS) [65], which is assumed to be the same as HIPS. The calculated pressure corrections indicate that FDMRheo pressure effects are real but minimal.

Viscous dissipation was found to be negligible based on the simulations described in the supplementary material [74], but the melt temperature as monitored by the intrusive melt thermocouple was found to deviate from the desired temperature in certain scenarios. First of all, it was observed that the nozzle temperature setting on the printer was not an accurate representation of the polymeric melt as discovered from the intrusive melt thermocouple readings, so the nozzle settings were adjusted as shown in Tables SIV and SV of the supplementary material [74] to achieve the desired melt temperatures listed in Table I. Once the nozzle settings were tuned to achieve the desired melt temperatures, the real-time temperature in the FDMRheo was observed to deviate from the desired melt temperature by no more than 2 °C in most cases as shown in Tables SIV–SVII [74]; however, in one extreme case for PC at 290 °C, the melt temperature could only reach 281 °C due to the maximum temperature settings of the printer. The viscosity and shear rate values were therefore shifted to account for these temperature inaccuracies.

The temperature shift factor, α, is calculated based on the temperature dependent WLF model [Eq. (6)] by calculating the predicted viscosity at the measured melt temperature, η0(Tmeas), and the predicted viscosity at the desired melt temperature, η0(Tdes), as shown in the following equation:

α=η0(Tdes)η0(Tmeas).
(12)

The viscosity and shear rate are then corrected using the following equations:

ηCorrected=α×η,
(13)
γ˙Corrected=γ˙α.
(14)

As discussed later in the article, the temperature inaccuracies were relatively negligible in all cases other than the FDMRheo data for PC at 290 °C.

To validate the accuracy of the pressure measurements acquired by the FDMRheo, the nozzle and FDMRheo setup was removed from the FDM printer and attached to the CapRheo as shown in Fig. 7. The nozzle connected to the capillary rheometer had a nozzle length of 1.36 mm long (all other results in this report are based on a nozzle with a 1.72 mm length, and the 3.74 mm length was only used to apply the Bagley correction). Despite the slightly shorter length, the pressure is measured in the same exact way, making this a valid approach for verifying the FDMRheo pressure against the CapRheo pressure transducer.

PC was processed through the CapRheo/FDMRheo setup at 270 °C, and the test was repeated four times at shear rates ranging from 100 to 4000 s−1 depending on the test. The results of the CapRheo pressure versus FDMRheo voltage for all four tests are combined in Fig. 8(a) and compared against the FDMRheo calibration curve. The continuous results from one of the CapRheo/FDMRheo tests are shown in Fig. 8(b). The results from both Figs. 8(a) and 8(b) show very good agreement between the CapRheo and FDMRheo pressures, confirming the accuracy of the FDMRheo against the standard pressure transducer of the CapRheo. The slight differences between the CapRheo and FDMRheo data shown in Fig. 8 are all within the accuracy of the FDMRheo load cell (The multiple calibration curves obtained for the FDMRheo indicated a 95% confidence uncertainty window of ±0.12 MPa at 0 MPa, which increases to ±0.165 MPa at 5 MPa.). All the data from this experiment are tabulated in Table SII in the supplementary material [74].

All the viscosity versus shear rate data are compared between the three rheometers in Fig. 9, and representative viscosity values are compared between the three rheometers in Table II. All the corrections described in Sec. III have already been applied to the FDMRheo and CapRheo data shown in Fig. 9 and Table II (note that for Table II, the data are based on equations that have been fitted to the corrected data to compare viscosity values at the same shear rates, which provides a good comparison with coefficients of determination of >0.99 for all HIPS data and 0.87–0.98 for all PC data). The FDMRheo measurements show excellent agreement with both off-line rheometers. The FDMRheo data nearly identically match the CapRheo data for HIPS, and the FDMRheo data fall between the RotRheo and CapRheo data for PC. Additionally, the deviation between the FDMRheo and the off-line rheometers is on the same order and even smaller than the deviation between the rotational and capillary rheometers, further validating that the in-line rheometer is equally as effective as the off-line rheometers for providing representative viscosity curves.

FIG. 8.

(a) Pressures measured by CapRheo (data points) versus the FDMRheo pressure curve (continuous line) as a function of FDMRheo voltage; (b) pressure versus time at four shear rates as set by the CapRheo.

FIG. 8.

(a) Pressures measured by CapRheo (data points) versus the FDMRheo pressure curve (continuous line) as a function of FDMRheo voltage; (b) pressure versus time at four shear rates as set by the CapRheo.

Close modal
FIG. 9.

Viscosity versus shear rate at three temperatures for [(a)–(c)] HIPS and [(d)–(f)] PC.

FIG. 9.

Viscosity versus shear rate at three temperatures for [(a)–(c)] HIPS and [(d)–(f)] PC.

Close modal
FIG. 10.

Comparison of Cross-WLF models generated by the three rheometers for PC and HIPS analyzed at 260 °C.

FIG. 10.

Comparison of Cross-WLF models generated by the three rheometers for PC and HIPS analyzed at 260 °C.

Close modal
TABLE II.

Tabulated viscosity values including percent deviation between the various rheometers.

MaterialTemperature (°C)Shear rate (s−1)FDMRheo (Pa · s)RotRheo (Pa · s)CapRheo (Pa · s)FDMRheo to RotRheo (%)FDMRheo to CapRheo (%)CapRheo to RotRheo (%)
PC 250 46 459 613 328 −25 40 −46 
250 140 428 558 311 −23 38 −44 
250 460 337 407 259 −17 30 −36 
PC 270 46 282 288 249 −2 13 −14 
270 140 265 272 224 −3 18 −18 
270 460 216 223 157 −3 38 −30 
PC 290 46 182 150 134 21 36 −11 
290 140 154 145 118 31 −19 
290 460 129 129 103 26 −21 
HIPS 225 46 844 1100 828 −23 −25 
225 140 373 527 408 −29 −9 −23 
225 460 156 240 191 −35 −19 −20 
HIPS 250 46 493 730 538 −32 −8 −26 
250 140 242 376 271 −36 −11 −28 
250 460 113 185 130 −39 −13 −30 
HIPS 275 46 337 473 363 −29 −7 −23 
275 140 179 270 203 −34 −12 −25 
275 460 91 148 109 −39 −17 −26 
MaterialTemperature (°C)Shear rate (s−1)FDMRheo (Pa · s)RotRheo (Pa · s)CapRheo (Pa · s)FDMRheo to RotRheo (%)FDMRheo to CapRheo (%)CapRheo to RotRheo (%)
PC 250 46 459 613 328 −25 40 −46 
250 140 428 558 311 −23 38 −44 
250 460 337 407 259 −17 30 −36 
PC 270 46 282 288 249 −2 13 −14 
270 140 265 272 224 −3 18 −18 
270 460 216 223 157 −3 38 −30 
PC 290 46 182 150 134 21 36 −11 
290 140 154 145 118 31 −19 
290 460 129 129 103 26 −21 
HIPS 225 46 844 1100 828 −23 −25 
225 140 373 527 408 −29 −9 −23 
225 460 156 240 191 −35 −19 −20 
HIPS 250 46 493 730 538 −32 −8 −26 
250 140 242 376 271 −36 −11 −28 
250 460 113 185 130 −39 −13 −30 
HIPS 275 46 337 473 363 −29 −7 −23 
275 140 179 270 203 −34 −12 −25 
275 460 91 148 109 −39 −17 −26 

For a direct comparison between the final data between the three rheometers, the corrected viscosities at all corrected shear rates are shown in Table SIII in the supplementary material [74]. The raw data (temperatures, uncorrected pressures, uncorrected viscosities, etc.) for HIPS and PC at both nozzle lengths are tabulated in Tables SIV–SVII in the supplementary material [74]. Based on the high accuracy of the FDMRheo, it could potentially be combined with the interesting work from Phan et al. [34] to improve the accuracy of their pressure measurements, which could enable a more easily transferable in-line rheometer based on power measurements from the stepper motor that drives the filament gears.

A Cross-WLF model was fitted to the rheological data for the three rheometers using a least-squares fit. To compare the models, the viscosity was plotted over a very wide shear rate range at an interpolated temperature of 260 °C as shown in Fig. 10. The HIPS viscosity curves match well for all three rheometers; for PC, the FDMRheo and RotRheo data are in good agreement, but the CapRheo Cross-WLF model has an inaccurate profile as seen in Fig. 10 and indicated by the fitted Cross-WLF data in Table III. The major difficulty in fitting a Cross-WLF model to CapRheo experiments is that the measurements do not typically extend to low enough shear rates to observe the Newtonian plateau, which likely caused some of the error in this fit.

TABLE III.

Cross-WLF coefficients fitted to the PC and HIPS data for all three rheometers.

PC CapRheoPC RotRheoPC FDMRheoHIPS CapRheoHIPS RotRheoHIPS FDMRheo
D1 (Pa · s) 2.32 × 1010 5.71 × 1012 5.18 × 109 4.91 × 109 2.20 × 1012 1.18 × 1014 
A1 23.6 31.4 22.2 19.8 27.6 33.2 
τ* (Pa) 3.90 × 104 5.03 × 104 4.78 × 104 2.74 × 104 3.03 × 104 19.2 × 104 
n 0.594 0.149 0.161 0.305 0.293 0.277 
A2 (K) 51.6 51.6 51.6 51.6 51.6 51.6 
T* (K) 385 385 385 371 371 371 
PC CapRheoPC RotRheoPC FDMRheoHIPS CapRheoHIPS RotRheoHIPS FDMRheo
D1 (Pa · s) 2.32 × 1010 5.71 × 1012 5.18 × 109 4.91 × 109 2.20 × 1012 1.18 × 1014 
A1 23.6 31.4 22.2 19.8 27.6 33.2 
τ* (Pa) 3.90 × 104 5.03 × 104 4.78 × 104 2.74 × 104 3.03 × 104 19.2 × 104 
n 0.594 0.149 0.161 0.305 0.293 0.277 
A2 (K) 51.6 51.6 51.6 51.6 51.6 51.6 
T* (K) 385 385 385 371 371 371 

The Cross-WLF coefficients for both materials from the different rheometers are shown in Table III. A2 was kept constant at 51.6 K, and T* was fixed as the glass transition temperature, but D1, A1, τ*, and n were fitted to all the data from the corresponding rheometer and material. All the data, except for the PC CapRheo data, match relatively well. Note from Eq. (6) that D1 and A1 counter one another in terms of their influence on the zero-shear viscosity; for example, when observing the HIPS data, D1 is largest for the FDMRheo data, but A1 is also largest for the FDMRheo, which balances the large D1 value and provides similar viscosity results as the other rheometers as shown in Fig. 10. The fits for τ* and n match very closely for all the rheometers, except in the case of the PC CapRheo data, which is at least partially caused by the small shear rate range of the CapRheo, focusing on high shear rates. The CapRheo data itself were actually very accurate and comparable to the other results, as displayed in Fig. 9, indicating that the Cross-WLF fit would be acceptable if a second capillary die with a larger diameter was used to extend the measurements to low shear rates; however, the CapRheo data were collected at relatively high shear rates to match the most common conditions measured on a CapRheo. One benefit of the FDMRheo over CapRheo is that the FDMRheo enables much faster testing (particularly at low shear rates) because the melt volume is much smaller in the FDMRheo, enabling more rapid pressure equilibration at each shear rate, which can be observed by comparing the pressure curves in Fig. 6 (FDMRheo) and Fig. 8(b) (CapRheo). Additionally, the FDMRheo can collect more measurements in a single test because there is no limitation in material volume like there is in the barrel of the CapRheo. The CapRheo offers the benefit of easily testing raw material in the pellet form, but the FDMRheo is an appropriate rheometer for analyzing FDM filaments.

Applying all the corrections described in Sec. III enabled excellent agreement between the three rheometers as shown in Fig. 9, but it is important to understand which corrections are the most critical and which could potentially be neglected. The influence of each correction can be visualized in Fig. 11 where the raw data for PC were modified by applying the corrections in the following order: Bagley end effects → Rabinowitsch → Temperature inaccuracies → Pressure effects. PC at 270 °C [Fig. 11(a)] is taken as the representative case, and PC at 290 °C [Fig. 11(b)] is also examined because this sample had a much larger temperature inaccuracy relative to the other samples. The influence of the corrections on HIPS was similar to those of PC that are shown in Fig. 11(a).

FIG. 11.

Changes in apparent viscosity as the raw data are modified with Bagley, Rabinowitsch, temperature, and pressure corrections: PC at (a) 270 °C and (b) 290 °C.

FIG. 11.

Changes in apparent viscosity as the raw data are modified with Bagley, Rabinowitsch, temperature, and pressure corrections: PC at (a) 270 °C and (b) 290 °C.

Close modal

Accounting for the entrance effects with the Bagley correction is clearly the most critical correction that must always be applied to accurately measure viscosity. The Rabinowitsch correction consistently causes a slight shift in the data, and it is a very simple correction to apply if viscosity is measured at multiple shear rates. The temperature correction is typically negligible, but it becomes very significant in the case of PC at 290 °C, so it is necessary to monitor temperature accuracy with the intrusive melt thermocouple. The pressure correction is negligible except for at elevated shear rates (i.e., when the pressure is highest).

Wall slip, which typically occurs above a critical wall shear stress [52], was neglected due to the relatively low wall shear stresses observed in this work, ranging from 0.005 to 0.2 MPa. Moreover, the excellent agreement between the in-line and off-line rheometers indicates that slip is not a significant source of error for HIPS and PC at the shear rates studied. Other materials may begin to slip at elevated shear stresses, which can be monitored by looking for the onset of stick-slip, which will show up as fluctuations in the real-time pressure data; the stable pressure data (Fig. 6) further indicated that slip did not occur during any of the tests described in this article.

The following paragraphs are incorporated to illustrate how the FDMRheo addresses some of the shortcomings that are common among on-line and in-line rheometers. Dealy describes the most important requirements of process rheometers to be: negligible leakage from the rheometer, no process interference, fast sample renewal rate, and a robust design [50]. Dealy continues to mention other important considerations, including preventing friction interference and maintaining measurement accuracy by accounting for nonisothermal effects, temperature uniformity, shear rate dependency, and by avoiding material degradation and build-up [50]. The following paragraphs describe how well these requirements have been met with the proposed design.

The load transfer column was designed to fit into the pressure port with tight tolerances, but a slide fit was used to prevent any frictional forces from interfering with the measurement, which results in the possibility of melt leakage around the load transfer column. A series of 3–4 g parts were printed over a range of shear rates and temperatures to quantify the amount of leakage; during printing, negligible leakage was observed while printing most of the parts, but some leakage was observed during high pressure prints. When observed, the melt leakage was collected and weighed: the melt leakage was found to weigh between 0 and 14 mg, which corresponded to 0–0.4% of the weight of the actual part. These amounts of leakage have no influence on the process or part quality, but it means the nozzle needs to be cleaned after printing a handful of parts or after printing a large part. Leakage could be further minimized through higher precision machining or by making the load transfer column and nozzle out of the same materials for better thermal expansion matching. Another potential concern about melt leakage is the viscous resistance from the melt surrounding the load transfer column that could limit the velocity of the transfer column pressing into the load cell. Viscous calculations based on the Couette flow indicate that instantaneous pressure changes would be resisted, causing a time delay on the order of 1 s. Practically, this does not influence steady-state measurements and should not add additional time delays, considering the already existing delays that can be caused by natural melt compressibility and relaxation [51]. These response times of about 1 s are fast compared to commercial in-line and on-line rheometers, some of which can have response delays on the order of minutes [66].

The pressure port and load transfer column are observed to create negligible process interference. Although the pressure port creates a slight recess (the recess varies from about 0.2 to 1 mm, depending on the pressure in the melt stream) from the circular cross section of the melt stream [as shown in Fig. 2(c)], estimations using the Hagen–Poiseuille equation indicate that changes in the total pressure drop due to the pressure port are negligible. In contrast, Lodge has indicated the importance of flush mounted pressure transducers to avoid such errors [67]; however, Lodge's intended application was for pressure transducers mounted directly in the small diameter, high shear rate region of interest where any error in shear rate has a profound effect on the measurement. In this FDM application, the pressure port is mounted in the large diameter section above the critical 0.5 mm diameter tip as shown in Fig. 2(c); therefore, estimates indicate that changes in total pressure due to the recessed pressure port will be small, on the order of 0.15%. Similar reasoning can be applied to understand the influence of how the thermocouple may disrupt the flow through the process; the pressure drop caused by the thermocouple is estimated to be smaller than that caused by the recess in the pressure port. The negligible influence of the pressure port and thermocouple is further validated by the excellent agreement between the in-line and off-line rheological measurements after applying the Bagley correction.

Another potential concern caused by the recess from the pressure port is the ability for melt stagnation and degradation. Since collecting the data displayed in Fig. 9, 25 parts were printed with a total printing time of more than 21 h. The load transfer column was then removed, and the melt was purged: no discoloration of the melt was observed, indicating no visual degradation. Additionally, the change between polymers of different colors (white to black and back to white) was very responsive without a long tailing edge of residence time distribution and without any residual black or white material coming out during a later print. Next, the printer was run for an additional 10 h, and the shear rate sweeps were repeated for PC at all three temperatures both before and after the 10 h of printing. The results for PC at 250 °C are shown in Fig. 12. The results overlay nearly perfectly with no trend of either increasing or decreasing pressure, indicating that the measurement is very consistent without any build-up of the degraded polymer.

FIG. 12.

Comparison of melt pressures before and after printing parts for 10 h; the numbers indicate the average change in pressure between the “original measurement” and after printing for 10 h.

FIG. 12.

Comparison of melt pressures before and after printing parts for 10 h; the numbers indicate the average change in pressure between the “original measurement” and after printing for 10 h.

Close modal

Perhaps, the largest practical limitation of the design is its industrial robustness. The load cell protrudes away from the nozzle, occupying space around the nozzle and potentially conducting heat away from the nozzle. The load cell could be damaged if it hits against something or if excessive pressure is applied. While the design is excellent for research purposes and has provided consistent results, a more space efficient and durable design would be best for widespread use. These problems could be solved, and the design could become commercially viable by designing a smaller pressure transducer to be integrated directly into the nozzle.

Viscosity temperature dependence is very important for polymer melt in-line rheometers, and it was discovered that the in-line thermocouple is critical to the success of FDM in-line rheometers because the melt temperature setpoint is typically not the same as the actual melt temperature. As done in this paper, new materials should be analyzed at multiple temperatures to quantify the viscosity temperature dependence [for example, the WLF equation shown in Eq. (6) can be used] so that temperature errors can be corrected by shifting viscosity to a reference or desired temperature. Interestingly, it was observed that at high shear rates (∼500 s−1), melt temperatures begin to decrease with increasing shear rate as displayed in Tables SIV–SVII, indicating insufficient melting and nonuniform temperatures across the nozzle. The in-line rheometer is, therefore, an excellent tool for studying heat transfer and high-speed limitations during the FDM process. Although the high flow rate limitations were not significant for the shear rates explored in this paper, heat transfer limitations are of both practical and theoretical interest, so future work will use the in-line rheometer to further explore this concept while extruding at higher shear rates.

The other important temperature consideration is the fact that the nozzle is not directly heated; it is heated by conduction with the heating element in the melting section of the hot end. Therefore, the nozzle material is cooler at its exit tip than at its entrance (as validated through manual temperature measurements), which may explain why many printers use a short nozzle with a small length to radius (L/R) ratio in the nozzle tip. Short L/R ratios of about 7 and 15 were used in this paper to resemble a typical printer nozzle and to prevent the melt from cooling through the length of the nozzle. While shorter length nozzles are better for both temperature uniformity and lower pressure drops (lower pressure drops prevent inaccurate flow rates due to skipped steps from the stepper motor or filament slippage in the gears), they result in smaller lengths of fully developed flow and a greater influence of end effects as made apparent by the Bagley corrections in Fig. 11. Due to the excellent agreement between all the in-line and off-line data after the Bagley correction, the authors found that L/R ratios of 7 and 15 enable a successful in-line rheometer.

The Bagley correction must be obtained by using two nozzles of differing L/R ratios, which is practically time-consuming to do on a regular basis. This could be overcome by predicting the entrance pressure. Although entrance effects have been modeled via many approaches in the literature [47,57,68–73], the authors are not aware of any analytical expressions that can be applied in real time that are universally accepted for describing the entrance pressure for all polymer melts. Therefore, the Bagley procedure is the most universally accurate method for analyzing the entrance pressure. Real-time corrections can be applied by fitting a curve to the entrance pressures obtained from the Bagley correction. The Bagley procedure could therefore be performed once for each material of interest, and the correction could then be applied for all future testing. The supplementary material [74] describes how the Bagley entrance pressures can be fitted to an equation for real-time analysis and control.

By accounting for entrance pressure, the in-line rheometer showed excellent accuracy when printing into open space; however, exit pressures must be accounted for when printing actual parts. When printing with small layer heights, elevated pressures are observed due to the pressure drop experienced between the nozzle and the previous layer. Additional work is on-going to measure and predict this exit pressure. Once the exit pressures are characterized, the in-line rheometer can be used for accurate viscosity measurements throughout an entire print. The real-time viscosity analysis will enable a better understanding of the FDM processing window and the role of polymer rheology in process optimization. In-line rheology and pressure measurements can also be used for real-time process control and predictions of quality metrics, such as strength and part density.

The in-line FDM rheometer (FDMRheo) designed in this work has been found to provide very accurate viscosity measurements as confirmed against off-line rheometers. The FDMRheo can collect data across a wide range of temperatures and shear rates to generate a successful Cross-WLF model for analyzing continuous viscosity curves as a function of temperature, shear rate, and pressure. The sensor for measuring the filament feed rate as well as the thermocouple for measuring melt temperature were both critical for the deployment and accuracy of the in-line rheometer. Entrance effects were the most significant correction for obtaining an accurate viscosity, so the Bagley correction should be applied to allow the FDMRheo to be used for real-time process control of the FDM process. For example, a control scheme could be developed to optimize the printing speed while maintaining pressures and viscosities within the ideal processing window. The FDMRheo is suitable for analyzing the viscosity of new, 3D printable materials to more rapidly introduce new materials to the market; the vision is that the rheometer can enable automatic process optimization and quality assurance using physics-based models for weld fusion (i.e., interlayer strength), residual stress, print density, and shrinkage.

This work was funded jointly by Saint-Gobain Research North America and the University of Massachusetts, Lowell. The authors would also like to thank David Shackleford of Saint-Gobain for his insightful discussions on rheological flow, Steve Johnston of UMass Lowell for his recommendation of the CapRheo/FDMRheo experiments, and Brian Beauvais and Phil Desroches of Saint-Gobain for their support in performing the x-ray computed tomography experiments.

1

Kellner, T. (2014, June 23), “Fit to Print: New Plant Will Assemble World’s First Passenger Jet Engine With 3D Printed Fuel Nozzles, Next-Gen Materials—GE Reports,” http://www.gereports.com/post/80701924024/fit-to-print/ (accessed 19 March 2017).

2

Oak Ridge National Laboratory: Innovations in Manufacturing, “ORNL Revealed the 3D Printed Shelby Cobra—Next Generation of Manufacturing Takes on 50-Year-Old Icon,” http://web.ornl.gov/sci/manufacturing/media/news/detroit-show/ (accessed 19 March 2017).

1.
Coogan
,
T. J.
, and
D. O.
Kazmer
, “
Bond and part strength in fused deposition modeling
,”
Rapid Prototyping J.
23
,
414
422
(
2017
).
2.
Ahn
,
S.
,
M.
Montero
,
D.
Odell
,
S.
Roundy
, and
P. K.
Wright
, “
Anisotropic material properties of fused deposition modeling ABS
,”
Rapid Prototyping J.
8
,
248
257
(
2002
).
3.
Bellini
,
A.
, and
S.
Güçeri
, “
Mechanical characterization of parts fabricated using fused deposition modeling
,”
Rapid Prototyping J.
9
,
252
264
(
2003
).
4.
Bellehumeur
,
C.
,
L.
Li
,
Q.
Sun
, and
P.
Gu
, “
Modeling of bond formation between polymer filaments in the fused deposition modeling process
,”
J. Manuf. Process.
6
,
170
178
(
2004
).
5.
Sun
,
Q.
,
G. M.
Rizvi
,
C. T.
Bellehumeur
, and
P.
Gu
, “
Effect of processing conditions on the bonding quality of FDM polymer filaments
,”
Rapid Prototyping J.
14
,
72
80
(
2008
).
6.
Thomas
,
J.P.
, and
J.F.
Rodriguez
, “
Modeling the fracture strength between fused-deposition extruded rods
,” in
Solid Freeform Fabrication Symposium
(
University of Texas
,
Austin
,
TX
,
2000
), pp.
16
23
.
7.
Hill
,
N.
, and
M.
Haghi
, “
Deposition direction-dependent failure criteria for fused deposition modeling polycarbonate
,”
Rapid Prototyping J.
20
,
221
227
(
2014
).
8.
Seppala
,
J.
,
S. H.
Han
,
K. E.
Hillgartner
,
C. S.
Davis
, and
K. B.
Migler
, “
Weld formation during material extrusion additive manufacturing
,”
Soft Matter
13
,
6761
6769
(
2017
).
9.
Coogan
,
T. J.
, and
D. O.
Kazmer
, “
Healing simulation for bond strength prediction of FDM
,”
Rapid Prototyping J.
23
,
551
561
(
2017
).
10.
Turner
,
B. N.
, and
S. A.
Gold
, “
A review of melt extrusion additive manufacturing processes: II. Materials, dimensional accuracy, and surface roughness
,”
Rapid Prototyping J.
21
,
250
261
(
2015
).
11.
Ahn
,
D.
,
J.
Kweon
,
S.
Kwon
,
J.
Song
, and
S.
Lee
, “
Representation of surface roughness in fused deposition modeling
,”
J. Mater. Process. Technol.
209
,
5593
5600
(
2009
).
12.
Galantucci
,
L. M.
,
F.
Lavecchia
, and
G.
Percoco
, “
Experimental study aiming to enhance the surface finish of fused deposition modeled parts
,”
CIRP Ann. Manuf. Technol.
58
,
189
192
(
2009
).
13.
D’Amico
,
A. A.
,
A.
Debaie
, and
A. M.
Peterson
, “
Effect of layer thickness on irreversible thermal expansion and interlayer strength in fused deposition modeling
,”
Rapid Prototyping J.
23
,
943
953
(
2017
).
14.
Pennington
,
R. C.
,
N. L.
Hoekstra
, and
J. L.
Newcomer
, “
Significant factors in the dimensional accuracy of fused deposition modelling
,”
Proc. Inst. Mech. Eng. E J. Process Mech. Eng.
219
,
89
92
(
2005
).
15.
Greeff
,
G. P.
, and
M.
Schilling
, “
Closed loop control of slippage during filament transport in molten material extrusion
,”
Addit. Manuf.
14
,
31
38
(
2017
).
16.
Carneiro
,
O. S.
,
A. F.
Silva
, and
R.
Gomes
, “
Fused deposition modeling with polypropylene
,”
Mater. Des.
83
,
768
776
(
2015
).
17.
Chia
,
H. N.
, and
B. M.
Wu
, “
Recent advances in 3D printing of biomaterials
,”
J. Biol. Eng.
9
,
4
17
(
2015
).
18.
Serra
,
T.
,
J. A.
Planell
, and
M.
Navarro
, “
High-resolution PLA-based composite scaffolds via 3-D printing technology
,”
Acta Biomater.
9
,
5521
5530
(
2013
).
19.
Mani
,
M.
,
S.
Feng
,
B.
Lane
,
A.
Donmez
,
S.
Moylan
, and
R.
Fesperman
,
Measurement Science Needs for Real-Time Control of Additive Manufacturing Powder Bed Fusion Processes
(
U.S. Department of Commerce, National Institute of Standards and Technology
, Gaithersburg, MD,
2015
).
20.
Huang
,
Y.
,
M. C.
Leu
,
J.
Mazumder
, and
A.
Donmez
, “
Additive manufacturing: Current state, future potential, gaps and needs, and recommendations
,”
J. Manuf. Sci. Eng.
137
,
014001
(
2015
).
21.
Weiss
,
B.
,
D. W.
Storti
, and
M. A.
Ganter
, “
Low-cost closed-loop control of a 3D printer gantry
,”
Rapid Prototyping J.
21
,
482
490
(
2015
).
22.
Kazmer
,
D. O.
,
Plastics Manufacturing Systems Engineering
(
Hanser
,
Cincinnati
,
OH
,
2009
).
23.
Comb
,
J. W.
,
P. J.
Leavitt
, and
E.
Rapoport
, Velocity profiling in an extrusion apparatus, US patent No. 6,054,077 (25 April 2000).
24.
Zinniel
,
R. L.
, and
J. S.
Batchelder
, Volumetric feed control for flexible filament, US patent 6,085,957, (11 July 2000).
25.
Sbriglia
,
L. R.
,
A. M.
Baker
,
J. M.
Thompson
,
R. V.
Morgan
,
A. J.
Wachtor
, and
J. D.
Bernardin
,
Embedding sensors in FDM plastic parts during additive manufacturing
, in
Topics in Modal Analysis & Testing
(
Springer
,
Cham
,
2016
), pp.
205
214
.
26.
Kousiatza
,
C.
, and
D.
Karalekas
, “
In-situ monitoring of strain and temperature distributions during fused deposition modeling process
,”
Mater. Des.
97
,
400
406
(
2016
).
27.
Peyre
,
P.
,
Y.
Rouchausse
,
D.
Defauchy
, and
G.
Régnier
, “
Experimental and numerical analysis of the selective laser sintering (SLS) of PA12 and PEKK semi-crystalline polymers
,”
J. Mater. Proc. Technol.
225
,
326
336
(
2015
).
28.
Dinwiddie
,
R. B.
,
V.
Kunc
,
J. M.
Lindal
,
B.
Post
,
R. J.
Smith
,
L.
Love
, and
C. E.
Duty
,
Infrared imaging of the polymer 3D-printing process
, in
SPIE Sensing Technology+Applications
(
International Society for Optics and Photonics
,
San Diego
,
CA
,
2014
), p.
910502
.
29.
Seppala
,
J. E.
, and
K. D.
Migler
, “
Infrared thermography of welding zones produced by polymer extrusion additive manufacturing
,”
Addit. Manuf.
12
,
71
76
(
2016
).
30.
Wolszczak
,
P.
,
K.
Lygas
,
M.
Paszko
, and
R. A.
Wach
, “
Heat distribution in material during fused deposition modelling
,”
Rapid Prototyping J.
24
,
615
622
(
2018
).
31.
Davis
,
C. S.
,
K. E.
Hillgartner
,
S. H.
Han
, and
J. E.
Seppala
, “
Mechanical strength of welding zones produced by polymer extrusion additive manufacturing
,”
Addit. Manuf.
16
,
162
166
(
2017
).
32.
D’Amico
,
A.
, and
A. M.
Peterson
, “
An adaptable FEA simulation of material extrusion additive manufacturing heat transfer in 3D
,”
Addit. Manuf.
21
,
422
430
(
2018
).
33.
McIlroy
,
C.
, and
P.
Olmsted
, “
Disentanglement effects on the welding behaviour of polymer melts during the fused-filament-fabrication method for additive manufacturing
,”
Polymer
123
,
376
391
(
2017
).
34.
Phan
,
D. D.
,
Z. P.
Swain
, and
M. E.
Mackay
, “
Rheological and heat transfer effects in fused filament fabrication
,”
J. Rheol.
62
,
1097
1107
(
2018
).
35.
Tadmor
,
Z.
, and
C. G.
Gogos
,
Principles of Polymer Processing
(
John Wiley & Sons
,
Hoboken
,
NJ
,
2013
).
36.
Rao
,
N. S.
, and
N. R.
Schott
,
Understanding Plastics Engineering Calculations Hands-on Examples and Case Studies
(
Hanser
,
Cincinnati
,
OH
,
2012
).
37.
Rauwendaal
,
C.
,
Polymer Extrusion
(
Hanser
,
Cincinnati
,
OH
,
2014
).
38.
Fan
,
B.
,
D. O.
Kazmer
, and
R.
Nageri
, “
An analytical non-Newtonian and nonisothermal viscous flow simulation
,”
Polym. Plast. Technol. Eng.
45
,
429
438
(
2006
).
39.
Kazmer
,
D. O.
,
Injection Mold Design Engineering
(
Hanser
,
Cincinnati
,
OH
,
2016
).
40.
Ramanath
,
H. S.
,
C. K.
Chua
,
K. F.
Leong
, and
K. D.
Shah
, “
Melt flow behaviour of poly-ε-caprolactone in fused deposition modelling
,”
J. Mater. Sci. Mater. Med.
19
,
2541
2550
(
2008
).
41.
Bellini
,
A.
,
S.
Guceri
, and
M.
Bertoldi
, “
Liquefier dynamics in fused deposition
,”
J. Manuf. Sci. Eng.
126
,
237
246
(
2004
).
42.
Han
,
C. D.
,
M.
Charles
, and
W.
Philippoff
, “
Rheological implications of the exit pressure and die swell in steady capillary flow of polymer melts. I. The primary normal stress difference and the effect of L/D ratio on elastic properties
,”
Trans. Soc. Rheol.
14
,
393
408
(
1970
).
43.
McIlroy
,
C.
, and
P. D.
Olmsted
, “
Deformation of an amorphous polymer during the fused-filament-fabrication method for additive manufacturing
,”
J. Rheol.
61
,
379
397
(
2017
).
44.
Malguarnera
,
S. C.
, and
A.
Manisali
, “
The effects of processing parameters on the tensile properties of weld lines in injection molded thermoplastics
,”
Polym. Eng. Sci.
21
,
586
593
(
1981
).
45.
Covas
,
J. A.
,
J. M.
Nóbrega
, and
J. M.
Maia
, “
Rheological measurements along an extruder with an on-line capillary rheometer
,”
Polym. Test.
19
,
165
176
(
2000
).
46.
Raps
,
D.
,
T.
Köppl
,
L.
Heymann
, and
V.
Altstädt
, “
Rheological behaviour of a high-melt-strength polypropylene at elevated pressure and gas loading for foaming purposes
,”
Rheol. Acta
56
,
95
111
(
2017
).
47.
Köpplmayr
,
T.
,
H.
Luger
,
I.
Burzic
,
M. G.
Battisti
,
L.
Perko
,
W.
Friesenbichler
, and
J.
Miethlinger
, “
A novel online rheometer for elongational viscosity measurement of polymer melts
,”
Polym. Test.
50
,
208
215
(
2016
).
48.
Coates
,
P. D.
,
A. R.
Haynes
, and
R. G.
Speight
, “
In-line characterization of polymer deformation in melt and solid phase processing
,”
Polymer
35
,
3831
3843
(
1994
),
49.
Wiklund
,
J.
,
B.
Birkhofer
,
S. A.
Jeelani
,
M. T.
Stading
, and
E. J.
Windhab
, “
In-line rheometry of particulate suspensions by pulsed ultra-sound velocimetry combined with pressure difference method
,”
Appl. Rheol.
22
,
42232
(
2012
).
50.
Dealy
,
J. M.
, “
Challenges in process rheometry
,”
Rheol. Acta
29
,
519
522
(
1990
).
51.
Moshe
,
A.
,
D. O.
Kazmer
,
S. P.
Johnston
,
R. M.
Malloy
, and
S.
Kenig
, “
Analysis of variance in capillary rheometry
,”
Polym. Eng. Sci.
56
,
895
904
(
2016
).
52.
Ansari
,
M.
,
T.
Zisis
,
S. G.
Hatzikiriakos
, and
E.
Mitsoulis
, “
Capillary flow of low-density polyethylene
,”
Polym. Eng. Sci.
52
,
649
662
(
2012
).
53.
Nikzad
,
M.
,
S. H.
Masood
,
I.
Sbarski
, and
A. M.
Groth
, “
Rheological properties of a particulate-filled polymeric composite through fused deposition process
,”
Mater. Sci. Forum
654
,
2471
2474
(
2010
).
54.
Zhong
,
W.
,
F.
Li
,
Z.
Zhang
,
L.
Song
, and
Z.
Li
, “
Short fiber reinforced composites for fused deposition modeling
,”
Mater. Sci. Eng. A
301
,
125
130
(
2001
).
55.
Wang
,
X.
,
M.
Jiang
,
Z.
Zhou
,
J.
Gou
, and
D.
Hui
, “
3D printing of polymer matrix composites: A review and prospective
,”
Compos. B Eng.
110
,
442
458
(
2017
).
56.
Gilmer
,
E.L.
,
D.
Miller
,
C.A.
Chatham
,
C.
Zawaski
,
J. J.
Fallon
,
A.
Pekkanen
,
T. E.
Long
,
C. B.
Williams
, and
M. J.
Bortner
, “
Model analysis of feedstock behavior in fused filament fabrication: Enabling rapid materials screening
,”
Polymer
152
,
51
61
(
2018
).
57.
Mackay
,
M. E.
,
Z. R.
Swain
,
C. R.
Banbury
,
D. D.
Phan
, and
D. A.
Edwards
, “
The performance of the hot end in a plasticating 3D printer
,”
J. Rheol.
61
,
229
236
(
2017
).
58.
Turner
,
B. N.
,
R.
Strong
, and
S. A.
Gold
, “
A review of melt extrusion additive manufacturing processes: I. Process design and modeling
,”
Rapid Prototyping J.
20
,
192
204
(
2014
).
59.
Cox
,
W. P.
, and
E. H.
Merz
, “
Correlation of dynamic and steady flow viscosities
,”
J. Polym. Sci.
28
,
619
622
(
1958
).
60.
Bagley
,
E. B.
, “
End corrections in the capillary flow of polyethylene
,”
J. Appl. Phys.
28
,
624
627
(
1957
).
61.
Kim
,
Y. H.
, and
R. P.
Wool
, “
A theory of healing at a polymer-polymer interface
,”
Macromolecules
16
,
1115
1120
(
1983
).
62.
Wool
,
R. P.
,
B.-L.
Yuan
, and
O. J.
McGarel
, “
Welding of polymer interfaces
,”
Polym. Eng. Sci.
29
,
1340
1367
(
1989
).
63.
Dees
,
M.
,
M.
Mangnus
,
N.
Hermans
,
W.
Thaens
,
A.
Hanot
, and
P. V.
Puyvelde
, “
On the pressure correction of capillary melt rheology data
,”
Rheol. Acta
50
,
117
124
(
2011
).
64.
Sedlacek
,
T.
,
M.
Zatloukal
,
P.
Filip
,
A.
Boldizar
, and
P.
Saha
, “
On the effect of pressure on the shear and elongational viscosities of polymer melts
,”
Polym. Eng. Sci.
44
,
1328
1337
(
2004
).
65.
Aho
,
J.
, and
S.
Syrjälä
, “
Measurement of the pressure dependence of viscosity of polymer melts using a back pressure-regulated capillary rheometer
,”
J. Appl. Polym. Sci.
117
,
1076
1084
(
2010
).
66.
Dealy
,
J. M.
, “
Viscometers for online measurement and control
,”
Chem. Eng.
91
,
62
70
(
1984
).
67.
Lodge
,
A. S.
, Flush pressure transducers for measuring pressures in a flowing fluid, US patent 4141252A (2 February 1979).
68.
Kwon
,
T. H.
,
S. F.
Shen
, and
K. K.
Wang
, “
Pressure drop of polymeric melts in conical converging flow: Experiments and predictions
,”
Polym. Eng. Sci.
26
,
214
224
(
1986
).
69.
Cogswell
,
F. N.
, “
Converging flow of polymer melts in extrusion dies
,”
Polym. Eng. Sci.
12
,
64
73
(
1972
).
70.
Bersted
,
B. H.
, “
Refinement of the converging flow method of measuring extensional viscosity in polymers
,”
Polym. Eng. Sci.
33
,
1079
1083
(
1993
).
71.
Cogswell
,
F. N.
, “
Measuring the extensional rheology of polymer melts
,”
Trans. Soc. Rheol.
16
,
383
403
(
1972
).
72.
Padmanabhan
,
M.
, and
C. W.
Macosko
, “
Extensional viscosity from entrance pressure drop measurements
,”
Rheol. Acta
36
,
144
151
(
1997
).
73.
Binding
,
D. M.
, “
An approximate analysis for contraction and converging flows
,”
J. Non-Newtonian Fluid Mech.
27
,
173
189
(
1988
).
74.
See supplementary material at https://doi.org/10.1122/1.5054648 for additional details regarding the design of the in-line rheometer components, validation of the volumetric flow rates calculated by the pulses from the stepper motor and encoder, additional equations and analysis for calculating entrance pressure effects, and the shear heating analysis in the FDMRheo. The final, corrected viscosities are tabulated to compare the values from the three rheometers. The raw data of the equilibrated pressures, shear rates, viscosities, etc. are also tabulated in the supplementary material. Supplementary raw data, including the transient pressure, voltage, shear stress, shear rate, viscosity, and temperature, are provided in individual raw data files: HIPS 225C 1_72mmL, HIPS 250C 1_72mmL, HIPS 275C 1_72mmL, PC 250C 1_72mmL, PC 270C 1_72mmL, and PC 290C 1_72mmL. Note that these files are only for one of the two tests at each condition for the 1.72 mm length nozzle tip; the raw data provided contain data every 8 ms (data were collected every 4 ms) to reduce file sizes. Two sets of time data are included in each raw data file; “Time” refers to the time for the pressure, voltage, shear stress, shear rate, and viscosity data, while “Temp Time” refers to the time for the nozzle temperature data, which is collected at a slower rate.

Supplementary Material