1. Introduction
Fused-filament fabrication (FFF), also known as fused deposition modeling or material extrusion additive manufacturing, is one of the most widely used processes in additive manufacturing for fabricating three-dimensional (3D) plastic structures [
1]. In FFF, deposited strands are stacked together to form a plastic structure (
Figure 1). The bonding strength between the stacked strands is one of the important factors that determine the strength of the fabricated plastic product. In addition to the temperatures of the filaments [
2], the bonding width between the cross-sections of two neighboring strands also affects the bonding strength (
Figure 1). A large bonding width leads to a large bonding area, and thus a high bonding strength. As such, it is important to know the relation between the bonding width and processing parameters (such as extrusion and print speed).
In FFF, a solid filament is melted in an extruder, extruded from its nozzle, and then deposition is carried out on a build platform (
Figure 2a). Several analytical models have been developed to understand the melting processes and melt flows inside extruders. A case was considered in [
3], where the filament melted as soon as it entered the extruder. A limiting case was investigated in [
4], where a solid filament only began to melt when it reached near the end of the extruder [
4]. The effects of feed rate on the temperature profile and feed force were considered in [
5]. The upper bound of feed rates was determined in [
6] that did not cause jamming during the extrusion.
During the deposition, the cross-section of an extrudate is deformed when the extrudate contacts the build platform (
Figure 2a). This deformation decides the bonding width. Previously, extensive studies have been carried out on the deformation of an object in its contact region with another object, such as a solid sphere or cylinder pressed against a solid substrate [
7], water drop resting on a flat surface [
8], water drop rolling on an inclined surface [
8], drops of n-heptane [
9] and water [
10] impacting solid substrates, and molten particles of tin [
11], alumina [
12], and lead [
13] colliding with cold surfaces.
Recently, other researchers have done much work to specifically explore the deposition process in FFF, through numerical simulation [
14,
15,
16,
17] and theoretical analysis [
18]. They have explored the effects of the speed ratio [
14,
15], viscosity [
15], viscoelasticity [
16], deposition flows [
17], and interlayer contact pressures [
18] on the dimensions of the deposited strands for particular combinations of processing parameters. As indicated in [
19], a simple method is needed to calculate the geometries of the deposited filament for FFF. For this purpose, an empirical formula was proposed in [
19] to predict the widths of the deposited strands. Naturally, a rigorously derived formula, which involves less empirical hypotheses, should give a better prediction.
The liquids considered in [
8,
9,
10,
11,
12,
13] have dynamic viscosities smaller than 1 Pa s [
20]. Therefore, the viscous effect is not dominant over effects such as gravity [
8], inertia [
9], and surface tension [
10]. In contrast, the filament materials commonly used in FFF, such as acrylonitrile butadiene styrene (ABS), polylactic acid (PLA), and polylactic acid-polyhydroxybutyrate copolymer (PLA-PHB), have much higher viscosities at the temperatures that they are extruded [
14,
15,
21,
22]. The viscosities further increase with the decrease in temperature during the cooling of the extrudates. Consequently, the deformations of these extrudates should be more affected by viscous forces. The analysis of the corresponding deformation may be simplified by comparing the viscous effect with the other effects, such as those of surface tension, inertia, and gravity. This comparison may enable the identification of both significantly large and negligibly small factors on the deformation. Accordingly, it is used in the present study to help analyze the deformation of the deposited strands in FFF.
Based on this analysis, our aim is to derive simple theoretical relations, which can be used to quickly predict the geometries of the deposited filaments for different combinations of processing parameters, such as extrusion speed, print speed and nozzle height. Particularly, they can be applied to determine the combination of the processing parameters that may lead to a large bonding width.
2. Materials and Experimental Methods
Experiments are performed using ABS to validate the derived relations. Single strands are printed with white ABS filaments (MakerBot Company, Brooklyn, USA), using a Creality Ender Pro 3D printer (Creality 3D Technology Company, Shenzhen, China). The nozzle employed in the printer has an opening with a diameter of 0.4 mm. Temperature of printing platform and nozzle are maintained at 90 °C and 235 °C throughout the prints. The platform is kept at a temperature lower than glass-transition temperature for ABS as of 104.5 °C [
2], to avoid the effect of the bed temperature on the strands’ cross-sections.
A cuboidal wire model with a length of 120 mm is designed in software SOLIDWORKS and exported as STL file. Simplified 3D is used as the slicing software to import the STL file of the computer aid design model, and also to generate G-code embodying all printing parameters. The print speed (), which refers to the horizonal speed of the nozzle head, is set as 10 mm/s in the software. The extrusion speed (), which is the speed that the polymer melt exits the nozzle, is calculated according to the law of mass conservation: = Q/Ad, where Q denotes the amount of extruded material per unit time and Ad is the measured cross-sectional area of the deposited filament.
The deformations of extrudates are considered in four special cases, which will be discussed in
Section 4. In each case, the layer height input in the printer setting defines the nozzle height. As designed, a printed strand is 120 mm long. Its length is kept long enough to ensure that it has a uniform extrusion, primarily focusing on the mid-section. The strand is cut at the middle using a brand-new razor blade with one stroke. The cutting is square to the printing direction to expose its cross-section for analysis. No deformation is observed along the cutting edge that is big enough to alter the experimental result. Three repetitions are performed for each test and their results are used to calculate standard deviation. Average value of all measurements is taken as experimental value.
A Dino-Lite Pro Digital Microscope (AnMo Electronics Corporation, Taipei, Taiwan) is utilized to take a cross-sectional picture along the axial direction of the strand (
Figure 3d,
Figure 4d,
Figure 5e and
Figure 6c1–d4).
The contour of a strand’s cross-section is used to evaluate its area from the image captured via software AutoCAD. The height of each strand is measured manually with a digital caliper, which has a measurement error of 1 µm. This height is also used to calibrate the image. At the time of measuring a strand, it has already cooled down to a room temperature. As such, it is rigid, and has negligible elastic deformation during the measurement. Accordingly, the measurement error induced by the elastic deformation is negligible, in comparison with the cross-sectional dimensions of a strand, such as height and width, which are all over 100 µm.
A high-speed camera, Fastec TS3 (Fastec Imaging Company, San Diego, USA), is installed to record the deposition processes at 500 fps.
Figure 3b and
Figure 5b are extracted from the recorded videos. The corresponding images are calibrated according to the nozzle height.
3. Model
To release the residual stress induced during an extrusion, an extrudate may start to swell at the nozzle exit. This expansion is referred to as the die-swell effect. The swell ratio and terminal swell distance are parameters used to characterize this effect. The former parameter is the ratio between the diameters of the nozzle opening and swelled extrudate, while the latter is the distance where the maximum swell diameter is attained. These two parameters depend on the type of the used polymer, nozzle length/diameter ratio, and shear rate in the nozzle [
23,
24,
25,
26].
By default, the solid filaments currently used in FFF, such as those used in commercial 3D printers, have the same diameter of 1.75 mm. The nozzles normally have circular openings, with diameters in the range of 0.2 to 0.8 mm. Accordingly, the extrudate is considered to have a circular cross-section. The length/diameter ratio of the nozzle varies from one to five. The feed rate of a solid filament is in the order of 1 mm s
−1 [
1]. According to the law of conservation of mass, an extrudate near the nozzle exit has a speed in the order of 10 mm s
−1, as shown in the printing of ABS filaments [
2]. The corresponding shear rate in the nozzle is in the order of 10 s
−1.
Using a typical FFF extruder with this shear rate, we measured the swell ratios and terminal swell distances of extruded ABS filaments. The swell ratios varied between 1.05 and 1.10, while the terminal swell distances were smaller than 20% of the nozzle diameter. According to these data, even when the die-swell effect is considered, the extruded strand in FFF may have a diameter of the same order of magnitude as that of the nozzle opening before its deposition on the substrate. We denote this diameter as
D. Owing to the relaxation of the restraint of the nozzle and negligible air drag on its surface, the speed of the extrudate becomes uniform across its cross-section after it travels the terminal swell distance [
24,
27]. As
D is in the same order of magnitude as that of the nozzle opening, this speed is in the same order of magnitude as
.
Consider that the extrudate falls in a steady shape onto a stationary platform. We set up a rectangular
x − y coordinate system, as defined in
Figure 2a.
H denotes the distance between nozzle head and build platform,
s is the arc length measured from the exit of the nozzle,
ψ(
s) is the inclination of the centerline to the horizontal, and
is the speed of the extrudate (
Figure 2a). At the exit of the nozzle,
ψ (0)
and
=
.
The polymer melt is modelled as a power-law non-Newtonian fluid [
3,
21,
28,
29,
30,
31]. During the deposition process, an extrudate may be compressed. We assume that the elongation viscosity,
µ, has a power-law relation under the generalized Newtonian fluid (GNF) framework. Equation (1).
where
L and
m are the extensional consistency and thinning (or thickening) indices, respectively.
L changes with temperature. When
m = 1, Equation (1) is reduced to the constitutive relation for a Newtonian fluid. Let
K and
n denote the shear consistency and thinning (or thickening) indices, respectively. As in [
32], for simplicity, we also assume that
m = n. Define the Trouton ratio,
Tr =
L/K. Within the GNF framework,
Tr = 3, indicating that
L is in the same order as
K. For PLA,
n = 0.433, and
K = 3.54 × 10
4 Pa s
n at 170 °C (Table I of [
32]). In the case of ABS,
n = 0.32, and
K = 1.04 × 10
4 Pa s
n at 230 °C [
2]. 170 °C and 230 °C are extrusion temperatures that have been used to print PLA [
32] and ABS [
2], respectively.
According to the above discussion and material properties of the filaments used in FFF, Equation (2) [
1,
2,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33].
where
and
are surface tension and mass density of the polymer melt, respectively, and
L is estimated at the extrusion temperature.
The bonding is considered to occur mainly above the melting or glass-transition temperature. Within 2 s, the ABS melt was cooled from the extrusion temperature of 260 °C to its glass-transition temperature of 104.5 °C [
2,
34]. To ensure a larger bonding time, the travel time of the extrudate from the nozzle exit to the substrate has to be considerably smaller than 2 s. In addition, two results of [
2] are important. First, the cooling rate of the ABS extrudate is 10 °C per 0.1 s for the first second. Second, the viscosity of the ABS melt rapidly increases with the decrease in temperature. Therefore, to reduce the increase in the viscosity as well as to increase the bonding time, it is reasonable to limit the above travel time to values smaller than 0.1 s. Considering Equation (2), this time limit implies that the nozzle height should be smaller than 1 mm. In other words,
mm and
is less than 2.
We consider that, when a strand just leaves the nozzle, it has a circular cross-section. Let
and
, respectively, denote the compressed depth and flat width of the deposited strand’s bottom surface (
Figure 2b).
and
are the counterparts of
and
at the top surface of the deposited strand.
is the smaller value between
and
, and it is the bonding width between two stacked filaments.
An extruded strand may impact a solid surface in a manner similar to that a viscous thread [
27,
33,
34,
35,
36,
37,
38] or sheet [
38] falls from a nozzle onto a horizontal surface. The previous studies focused on the shapes and dynamics of the falling threads or sheets, but not on their deposition deformations. In addition, golden syrup was the liquid having the highest viscosity in the previous studies [
27,
33,
34,
35,
36,
37,
38]. Its viscosity is in the order of 10 Pa s [
27], which is at least three orders of magnitude lower than those of its counterparts in FFF. This difference also indicates that the viscous effect should be more dominant in FFF.
Considering these results and previous modelling studies on the viscous threads and sheets [
27,
33,
34,
35,
36,
37,
38], particularly that of [
27], in the present study we investigate the deformations of the extruded strands at their interfaces with the substrates. As will be detailed in the rest of this section, as well as in
Section 4, we begin this investigation with the theoretical model of [
27], which was originally developed for a Newtonian fluid. We extend it to consider a non-Newtonian fluid, and then simplify it by mainly considering the viscous effect. Finally, with the aid of the derived formulas and according to the balance of linear momentum and also the balance of energy, we explore the deposition deformation, obtaining theoretical relations for predicting the geometries of the deposited strands.
The following axial and transverse equations of motion were derived in [
27] for a viscous strand based on the extensional-flow approximation, Equations (3) and (4):
where
F is the axial force acting over the cross-section of the suspended strand. The viscous resistance to stretching or compression was considered in these two equations, while the viscous resistance to shear deformations was neglected. The shear deformations might be caused by the bending and air drag, which were neglected in [
27]. It was also indicated in [
27] that the neglect of the strand’s bending precluded the buckling mechanism. Since the nozzle height considered in the present study is not large, the buckling does not occur. As such, the bending deformation is not large. Meanwhile, the drag force of the surrounding air is negligible, due to the small density of air and relatively low extrusion speed. Therefore, Equations (3) and (4) still apply to the present study to consider the shape and dynamics of an extrudate before it falls on the build platform.
In addition, during the collision and molding processes that will be considered later, the strand mainly experiences compression. Accordingly, the viscous resistance to stretching or compression is considered, while the viscous resistance to shear deformations is still neglected.
The expression of
F was also derived in [
27] for a Newtonian fluid. Following its derivation in [
27] while considering that
µ has the expression of Equation (1), we obtain an expression of
F for a power-law non-Newtonian fluid, Equation (5):
According to the relations in Equation (2) and the value of
H, we obtain the following scaling relations, Equation (6):
All numbers in Equation (6) have the same unit of 1/sm. In the first relation, as m > 0, 20m > 1. The second, third, and fourth relations in Equation (6) represent the relative effects of the surface tension, inertia, and gravity with the viscous force, respectively. According to these relations, these three effects can be neglected in comparison with the viscous effect.
In addition, at a scaling level, we obtain, Equation (7):
Using Equations (5) and (7), both Equations (3) and (4) can be scaled to Equation (8):
Considering Equation (6), the second term in Equation (8) is considerably larger than the other three terms, which indicates that these three terms can be neglected in this equation. Consequently, Equations (3) and (4) are simplified to Equations (9) and (10)
The solution to the set of Equations (5), (9), and (10) can be obtained by, Equations (11)–(14):
If
in case
we obtain
where
c is a constant.
ABC and
JBCI are, respectively, the side and cross-sectional views of the compressed portion of the extrudate at its contact region with the substrate (
Figure 2c). Terms
Fx and
Fy denote the reaction forces applied by the substrate to the extrudate along the
x and
y directions, respectively,
is the distance between
A and
B, while
and
are the axial forces at the two ends of the extrudate. The two ends may have different temperatures and thus different consistency indices. Therefore, when Equation (5) is used to determine
and
, the corresponding values of
L may be different. As discussed in
Section 2, the time limit ensures that the temperature difference between the two ends of the extrudate is small, such that the consistency indices at these ends have the same order of magnitude. The corresponding consistency indices are not distinguished thereafter.
As what was done in [
39], considering the balance of linear momentum on the portion of the extrudate from the nozzle to the end of the contact region (
Figure 2c), Equations (15) and (16):
where
G denotes the weight of this portion of the extrudate.
Further, we analyze the balance of energy during the collision. In this process, the bottom of the extrudate’s cross-section is compressed. Therefore, the work is mainly carried out by
Fy along the vertical direction. It is, Equation (17):
The viscous dissipation is, Equation (18) [
9]
where
denotes the dissipation function,
is the characteristic volume of the extrudate’s compressed portion, and
is the characteristic time of formation of the compressed portion. At the contact region, the extrudate is impacted mainly by the stage along the
y direction. Therefore, we focus only on the viscous dissipation along this direction. Consequently, we obtain, Equations (19)–(21):
tc is estimated by the travel time of the extrudate through the contact region,
4. Results and Discussions
We consider that the deposition includes two processes: collision and molding. During the collision, only the bottom portion of the extrudate has contact with the substrate. It may be compressed. However, the top portion is not, since it does not have direct contact with the substrate during the collision. After the collision, the extrudate goes through the molding process. The substrate and nozzle serve as bottom and top parts of the mold, respectively. In the molding process, both top and bottom portions of the extrudate may be compressed.
We analyze two cases according to the relation between
and
(
Table 1). In Cases I and II,
and
, respectively. We investigate two subcases in either case. In the first subcase, after the collision, the thickness of the extrudate is less than
H. However, it is larger than
H in the second subcase. Therefore, in the first subcase, the extrudate only collides with the sustrate, while it is not molded. In the second subcase, it is further molded after the collision.
If the extrudate is not compressed by either substrate or nozzle during its deposition, its cross-sections are still circular in shape. According to the law of conservation of mass, the corresponding diameter,
, meets the following relation, Equation (22):
4.1. Case I.1
In this subcase, we assume that, Equation (23)
As indicated in [
27], before the extrudate collides with the substrate, it has a catenary shape (
Figure 3a). This is validated by our experimental result (
Figure 3b). Thus, Relation (11) can be used. Therefore, the extrudate has a constant speed of
before the collision, which implies that
. Using Equation (5), we obtain, Equations (24)–(26):
Subsequently, using Equations (6) and (24), according to Equation (16)
Furthermore, using Equations (6) and (25), according to Equation (17),
Along the
y direction, the maximum compressed size of the extrudate is
. The corresponding characteristic speed is
. Accordingly,
=
which, by Equation (21), yields, Equations (27)
Using Equation (27), according to Equation (19), Equation (28)
In addition, as what was done in [
7], according to the geometric relation in the region of
JBIC (
Figure 2c), Equation (29)
Using Equations (19)–(21) and (28), according to Equation (18), Equation (30)
Considering Equations (26) and (30), the balance of
W and
E leads to, Equation (31)
The combination of Equations (29) and (31) yields,
During the collision, as the top portion of the extrudate is not directly exposed to external forces, Equations (33) and (34):
According to Equations (31) and (33), Equation (34):
Equations (31)–(34) indicate that, during the collision, the extrudate’s cross-section experiences negligible deformation. As such, after the collision, the extrudate still has circular cross-sections, whose diameters are
(
Figure 3b). This explains a result published in [
15,
22]. For fast printing with a large nozzle height, it was found numerically in [
15] and experimentally in [
22] that the PLA extrudate’s cross-section was almost cylindrical. In their case,
and
. Accordingly, it belongs to Case I.1. Consequently, according to our results, the deposited strand should have almost circular cross-sections. The derived results for this subcase are also validated by experimental results. As shown in
Figure 3c1–c5, when
H was fixed to be 0.5 mm and
/
ranged from 0.6 to 1.0, the deposited extrudates have approximately circular cross-sections.
4.2. Case I.2
In this subcase, we assume that, Equation (36)
During the collision, although there is negligible compression, the speed at the centerline increases from
to
, and the diameter of the extrudate is changed from
D to
. If Equation (36) is satisfied, after the collision, the extrudate has a thickness larger than the nozzle height and it should go through a molding process (
Figure 4a,b). The extrudate is approximately axisymmetric before the molding. Owing to this symmetry, the top and bottom portions should be compressed by the same degree during the molding, which implies that. Equations (37) and (38):
As in [
15,
22], the compressed cross-section is assumed to have an oblong shape, which is a flat cuboid with rounded edges (
Figure 4b). Subsequently, by geometric analysis and with the assistance of Equation (22), Equations (39) and (40)
Consequently,
, and Equation (41)
The relations in Equations (39) and (40) are validated by experiments. The values of
that are predicted using Equation (39) have a good agreement with the corresponding experimental results (
Figure 4c). As further observed from
Figure 4d1–f4,
approximately equals
. The maximum difference is 16%, which appears in the sample shown in
Figure 4f4.
4.3. Case II.1
In this subcase, Relation (23) is also assumed to hold true. However, different from that in Case I.1, the extrudate drops along the vertical direction, and then has a sharp turn as it lands on the contact region (
Figure 5a). This has been validated by the experimental result (
Figure 5b).
Since Relation (23) is satisfied, the top portion of the extrudate is not directly exposed to external forces during the collision. As in Case I.1, the top portion should still have an approximately circular profile. However, the bottom portion may be compressed (
Figure 5c). Next, the focus is on finding its compressed width and depth.
In the vertical portion of the extrudate, Equation (13) is satisfied and
decreases from
to
. Consequently, according to Equation (5), Equation (42):
In the curved part of the extrudate, Equation (14) is satisfied. Therefore, in this part, the speed at the centerline remains
, which indicates that
Using Equation (5), we obtain, Equation (43):
Subsequently, using Equations (6) and (42), according to Equations (15) and (16), Equations (44) and (45):
Equations (42)–(45) indicate that
and
are the two main forces applied on the suspended portion of the extruded strand. To obtain the balance of moments on this portion, the lateral area of the collision and cross-sectional area of the extrudate, which are the areas where
and
apply, should be approximately equal
D. Therefore, we obtain Equation (46):
The combination of Equations (17) and (45) yields, Equation (47)
In addition, as in Case I.1, can also be expressed by Equation (27). According to Equations (19), (46), and (27), ø can be expressed by Equation (28), while E can be expressed by Equation (30).
Using Equations (29) and (46), according to Equations (47) and (30), the balance of
W and
E leads to Equation (48)
As shown in
Figure 5d,e1–e5, this equation matches well with experimental tests. In calculating the relation between
and
using Equation (48),
m was set to be 0.32 [
2],
was 0.4 mm, and
H was 0.5 mm.
According to Equation (29), Equation (49)
Since the extrusion speed is higher than the print speed, the substrate applies a force during the collision to slow down the extrudate, making its speed match with the print speed. This force is represented by . As observed from Equation (45), it is not negligible, and compresses the bottom of the extrudate.
As in Case I.1, since the top portion of the extrudate is not directly exposed to external forces during the collision, Equations (33)–(35) are also valid in Case II.1.
In addition, to validate Equation (45),
was also determined, using a model that was originally applied for bi-directional rectilinear compression molding of polymers between two parallel plates [
40]. In the compression model, the polymer melt was assumed to be a Newtonian fluid. According to this model, Equation (50)
where
is the volume being compressed,
is the clamping speed, and
h is the gap between the two plates. In our case,
is the volume of the compressed portion of the extrudate (
Figure 2c). Thus, we have, Equation (51)
Using Equations (51), (52), (21) and (49), according to Equation (50), Equation (53)
It is identical to Equation (45), when m is set to be 1 in Equation (45) for a Newtonian fluid. As such, the compression model also validates our derived result.
4.4. Case II.2
In this subcase (
Figure 6), we assume that Equation (54)
In Equation (54), denotes the depth that the extrudate’s bottom is compressed during the collision. It equals of Case II.1, and also has the expression of Equation (49).
If Equation (54) is valid, the collided extrudate still has a thickness larger than the nozzle height. Consequently, a molding process follows the collision. The bottom portion of the extrudate is compressed during the collision, while the top portion is not. Accordingly, when the impacted extrudate is subjected to molding, the top and bottom portions of the extrudate are not geometrically symmetric. Consequently, they may be compressed by different degrees, leading to different flat areas at the top and bottom of the extrudate. Given that they are compressed by the same degree during the molding, we may still have
. For example, if the collided extrudate has a thickness just slightly larger than the nozzle height, both of its top and bottom portions are not compressed much during the molding process. As such, after the molding,
and
. According to the above discussions, in general,
. Consequently,
, which is validated by our experimental results. As shown in
Figure 6c1–d4, the difference between
and
ranges from 20 to 40%, except for the sample shown in
Figure 6c3, which has a 4% difference between these two values.
A simple geometric analysis shows that the thickness of the molded extrudate should be equal to the nozzle height, which implies that, Equation (55)
Accordingly, either
or
is smaller than
Considering Equation (29),
or
is smaller than
. With the assistance of Equation (22), this result indicates that, Equation (56)