arxiv: v1 [cond-mat.soft] 27 Mar 2017

Transcription

1 Disentanglement Effects on Weling Behaviour of Polymer Melts uring the Fuse-Filament-Fabrication Metho for Aitive Manufacturing C. McIlroy, P.D. Olmste Department of Physics, an Institute for Soft Matter Synthesis an Metrology, Georgetown University, Washington DC, 57, USA arxiv:73.995v [con-mat.soft] 7 Mar 7 Abstract Although 3D printing has the potential to transform manufacturing processes, the strength of printe parts often oes not rival that of traitionally-manufacture parts. The fuse-filament fabrication metho involves melting a thermoplastic, followe by layer-by-layer extrusion of the molten viscoelastic material to fabricate a three-imensional object. The strength of the wels between layers is controlle by interiffusion an entanglement of the melt across the interface. However, iffusion slows own as the printe layer cools towars the glass transition temperature. Diffusion is also affecte by high shear rates in the nozzle, which significantly eform an isentangle the polymer microstructure prior to weling. In this paper, we moel non-isothermal polymer relaxation, entanglement recovery, an iffusion processes that occur post-extrusion to investigate the effects that typical printing conitions an amorphous non-crystalline) polymer rheology have on the ultimate wel structure. Although we fin the wel thickness to be of the orer of the polymer size, the structure of the wel is anisotropic an relatively isentangle; reuce mechanical strength at the wel is attribute to this lower egree of entanglement. Keywors: Fuse filament fabrication, polymer melt, weling, isentanglement, non-isothermal. Introuction Fuse filament fabrication FFF) [] has become an essential tool for the rapi fabrication of custom parts via aitive manufacturing. Although there are numerous avantages to this technique [], incluing ease of use, cost an flexibility, improving the strength of printe parts to rival that of traitionally-manufacture parts remains an unerlying issue. The most common printing materials are amorphous polymer melts such as linear polycarbonate PC) [3] an acrylonitrile butaiene styrene ABS) [], a melt containing rubber nano-particles that provie toughness even at low temperatures. FFF printers can also hanle semicrystalline polymers such as poly-lactic aci PLA) [4]. The printing process involves melting a soli filament of the printing material an extruing it through a nozzle. To fabricate a three-imensional object the melt is eposite layer-by-layer, as illustrate in Fig.. The key to ensuring the strength of the final printe part is successful interiffusion an re-entanglement of the polymer melt across the layer-layer interfaces. In general, the wel strength of a polymer-polymer interface grows sub-iffusively with weling time as t /4 until the bulk strength plateau is reache [5, 6, 7]. Several molecular mechanisms are propose to explain this scaling. Since the wel thickness arising from interpenetration epth also scales as t /4 until the raius of gyration R g is reache ue to polymer reptation, one suggeste mechanism is that this interpenetration epth etermines the wel strength [8, 7, 9]. Others suggest that the formation of briges between the surfaces is the key strengthening mechanism [, ]. Both approaches are motivate by the iea of entanglement formation across the interface []. Whilst some stuies assume a simple proportionality between the interpenetration istance an entanglement formation [3, 8,, 4], others assume a minimum interpenetration istance for an entanglement to form [9, 5]. Few experiments fin that iffusion by the raius of gyration R g is require to achieve bulk strength [6], but in many cases wels reach bulk strength in much shorter times [3]. In FFF the weling behaviour is essentially a thermallyriven iffusive process [7], an interiffusion is limite as the melt rapily cools towars its glass transition temperature [8]. In aition, large shear rates in the nozzle eform an align the polymer microstructure prior to weling; it is suggeste that this alignment can affect the iffusive behaviour at the wel line causing e-boning [9]. The eformation inuce by the FFF extrusion an eposition process, which involves a 9 o turn, has recently been investigate using a molecularly-aware moel for a non-crystalline polymer melt []. In that paper both stretch an orientation of the polymer are incorporate using the Rolie-Poly constitutive equation [] an the entanglement ensity is allowe to vary with the flow []. Preprint submitte to Polymer March 9, 7

2 a) soli filament R pinch rollers b) printe wall t P L P m P heate section T N = 5 o C U L z b P t P- R U N L wel region z L P- m P- b P- t P- L P H m P U L L P- m P- L P- m P- z L P- m P- buil plate z motion of buil plate b P- z Figure : Simple schematic of typical FFF process, as escribe in text. a) In the frame of the moving nozzle, the melt exits the nozzle at spee U N an the buil plate moves at spee U L in the ŷ-irection. The current printe layer is enote L p an the mile of the layer is enote m p. b) Printe wall geometry with sites at the top t, bottom b an mile m of each layer labelle. Wels occur at the interface between layers. See Appenix, Tables A., A. for typical moel parameters. Flow through the nozzle followe by eposition into an elliptically-shape layer inuces complex, non-axisymmetric polymer configurations, with the polymer microstructure varying ramatically from the top to the bottom of the printe layer. This eformation significantly isentangles the polymer melt via convective constraint release [3] in an inhomogeneous way. Due to this eformation impose by the FFF extrusion process, interiffusion oes not necessarily occur from an equilibrium state. Non-equilibrium molecular ynamics NEMD) calculations of a iffusion tensor for relatively short polymer melts uner both planar Couette an planar elongational flow show a significant enhancement of iffusion parallel to the flow irection [4]. Recently, moels that incorporate an anisotropic shear-rate-epenent friction coefficient, or mobility tensor, have been propose to successfully reprouce the ynamics of polymers uner shear [5, 6]. Suitably accounting for flow-inuce friction-reuction effects is also require to quantitatively moel uniaxial extensional ata [7]. The simple ynamical moel of polymers uner shear with an anisotropic mobility tenor of Uneyama et al. [5] is consistent with NEMD simulations an experimental ata; polymer segment alignment is suggeste as the main cause for the anisotropy of the iffusion tensor. However, this moel is not expecte to apply o flows other than planar shear, an oes not capture the anisotropic relaxation ynamics of aligne polymers after flow cessation. Furthermore, ue to the nature of the eposition flow [] polymers on either sie of the interface resie in ifferent eformation environments. Thus, a mutual iffusion mechanism shoul be consiere, in which the iffusion coefficient epens on the local composition of chain mobility [8, 9]. In particular, the theory of Kramer et al. [9], suggesting that mutual iffusion is controlle by the mobility of faster moving chains, can successfully escribe experimentally measure iffusion coefficients [3, 3, 3]. In this paper, we investigate the post-extrusion iffusive behaviour at the wel between two printe layers of an amorphous non-crystalline) polymer. We use the proceure evelope by McIlroy & Olmste [] to calculate the polymer conformation tensor an corresponing isentanglement that is inuce by the extrusion process. We then introuce a spatio-temporal temperature profile to examine how this eformation relaxes in the wel region between two layers whilst cooling. In particular, we stuy how the structure of the wel evolves, how entanglements recover an calculate an interpenetration istance that incorporates both anisotropic an mutual iffusion effects. Finally, we aress how these wel properties are affecte by molecular weight, nozzle shear rate an print temperature.. FFF Moel.. Ieal Extrusion Process The soli filament feestock is melte within the nozzle. Recently Mackay et al. evelope a moel that solves an approximate energy balance to correlate the maximum fee velocity with the print temperature T N [33]. In a frame moving with the nozzle, the melt exits the nozzle in steay state at mass-average spee U N. It is then eposite onto the buil surface, where the material must spee up an eform to make a 9 o turn. The buil

3 a) Stage : Printing Lp- b) Stage : Printing Lp z z -H tp provies wel region for next layer mp tp- mp- H bp tp- wel mp- bp- H -R R bp- x -R R c) Spatio-Temporal Temperature Profile at x= 5 TN H tp mp -z Tg 5 Ta bp tp- mp- -H bp t-twp-) s) -4 ) Temperature Evolution at Sites 35 t b temperature oc) x p- mp- t-twp) s) mp bp p- - tp- mp- 5 TN 5 Tg 5 Stage -4 Ta Stage - t-twp-) s) -4 - t-twp) s) Figure : Schematic illustrating a) first printing stage where Lp is eposite at time twp ) an b) secon printing stage where Lp is eposite on top of Lp- at time twp). Wel sites on either sie of the layer-layer interface are enote tp an bp, an layer mipoints are enote mp an mp. c) Spatio-temporal temperature profile of two printing stages an ) temperature evolution at sites inicate in a) an b). Print temperature is TN = 5o C an ambient temperature is Ta = 95o C. Weling between layers begins at time twp) an the layer-layer interface z = approaches the glass transition temperature Tg = 4o C at time tw g = s. For comparison with Ref. [8], we consier an extrusion process that eposits a single filament or layer) in the xy-plane. Subsequent filaments are eposite on top of the previously-eposite filament to create a vertical printe wall in the z-irection Fig. ). Due to this geometry we use the term layer to refer to a single eposite filament; note that some authors use layer to refer to the planar geometry of the buil. Figs. a an b illustrate these two key stages of FFF printing; layer Lp is eposite at time twp ), followe by layer Lp on top of Lp at time twp). During the secon stage the previously-printe layer heats up an a wel forms between the ajacent layers. The temperature profile T t, x) of the two layers rives plate moves horizontally in relation to the nozzle at massaverage spee UL. Assuming mass conservation an eformation into a cylinrical eposite filament, the two spees are relate by πr UN = πrh UL, ) where R is the nozzle raius an H is the thickness of the eposition. Typically H < R, so that the eposite filament is elliptically shape. The eposite filament then cools from the print temperature TN towars the glass transition temperature Tg. Typical moel parameters are etaile in the Appenix Table A.). 3

4 the weling across the layer-layer interface z = in the region of the polymer size ±R g ); we enote wel sites either sie of the interface to be t p at the top of L p ) an b p at the bottom of L p ), respectively. Inter-iffusion between the layers occurs until the wel temperature reaches T g at time t W g... Temperature Profile In typical FFF processes, the nozzle is fixe at temperature T N, the buil plate is hel at the ambient temperature T a usually just below T g ), an printing occurs within an oven. The temperature profile of the oven is inhomogeneous ue to the moving nozzle, which is continually accelerating an ecelerating accoring to the print geometry an generating complex air-flow patterns. Thus, heat flow through the layers an exchange with the air is a very complicate problem. Typically, the bounary conitions at the layer-layer an layer-air interfaces are etermine by a combination of convection, conuction an raiation. However, the heat transfer coefficient to escribe this cooling process is not well unerstoo. We neglect temperature variations in the x-irection an moel the one-imensional temperature profile T t, z) across two layers z [ H, H], where z = is the layerlayer interface an z = H is the layer-air interface) via the one-imensional heat equation T t = αt T ) z. ) This gives the temperature evolution through the centre x = ) of the two layers Fig. b). For polycarbonate, which is a typical printing material, the thermal iffusivity α has a linear temperature epenence [34] of the form αt ) = α BT ), 3) for constants α, B an T g = 4 o C. Thermal iffusivity changes by 3% between the melt at T N an the soli at T a see Appenix, Table A.). We solve Eqs. -3 via an explicit finite-ifferencing scheme with the following bounary conitions. We assume Dirichlet bounary conitions at the interfaces such that T t, z = ±H) = T a, 4) an the initial temperature profile is a step function, { T N, z >, T t =, z) = T a, z <. Although Eqs. -5 give a crue approximation for cooling, this thermal protocol mimics experimental infrare measurements of the surface temperature of a printe ABS wall [8], which has similar thermal iffusivity to polycarbonate. In particular, Fig. 3 shows that the oneimensional moel well escribes the infrare-measure 5) 4 temperature o C) 4 6 wel line z= exp m p exp m p- sim m p sim m p- T g time s) Figure 3: Temporal temperature evolution at layer mipoints m p an m p uring the secon stage of printing compare to experimental infrare measurements of the surface temperature of a printe ABS wall [8]; ABS has similar thermal iffusivity to polycarbonate, a smaller T g = o C, a smaller nozzle temperature T N = o C an equivalent ambient temperature T a = 95 o C. The temperature evolution of the wel line is ifficult to measure using infrare so only the preiction is plotte. temperature evolution at layer mipoints m p an m p uring the secon stage of printing. The temperature at the wel line is har to measure using infrare imaging ue to the curvature of the surface; in Ref. [8] it is etermine by taking an average of the temperature evolution at m p an m p. Arguably, the one-imensional calculation shown in Fig 3 gives a more accurate escription of the wel line temperature evolution, although complex cooling ynamics uring eposition are neglecte. The iealise wall geometry use here also oes not capture the effects of ajacent layers in the xy-plane on the temperature profile. Figs. c shows the spatio-temporal temperature profile T t, z) preicte by Eqs. -5 for polycarbonate properties see Table A.) - the material we will consier henceforthly as in Ref. []. Fig. shows the temporal temperature evolution of sites at the top t, mile m an bottom b of each layer. We use a log-linear scale for comparison with later results. We choose T N = 5 o C an T a = 95 o C, as in typical FFF systems. Weling at the interface z = is governe by the temperature evolution an material ynamics on either sie of the wel line at wel sites t p an b p. These sites are chosen such that t p b p = R g, 6) where R g nm is the raius of gyration of an iniviual polymer. In particular, Fig shows how layer-air interface at t p rapily cools below T g uring stage. At t wp) the next layer is eposite, creating a layer-layer interface, an instantly heats up wel site t p ; the secon wel site b p instantly cools by the same egree. During stage, the layer-layer interface cools much slower than the layer-air interface reaching T g in approximately t W g = s. We will stuy how this wel line temperature evolution affects the relaxation of the printing material at wel sites t p an

5 b p, an consequently the characteristics of the wel..3. Polymer Dynamics We escribe the polymer microstructure using a moifie version Rolie-Poly moel [] that inclues flow-inuce changes in the entanglement fraction, as in Ref. []. Essentially the Rolie-Poly moel is a variation of the stanar Doi-Ewars tube moel for linear entangle polymer networks, which approximates the more powerful but unwiely microscopic GLaMM moel [35] to provie a simple one-moe constitutive equation for the stress tensor see Appenix Appenix A for more etails). Since surrouning chains restrict transverse motion in a melt, a polymer chain is restricte to a tube-like region. This tube represents topological constraints ue to entanglements [36]. At equilibrium, the entanglement number of a melt is relate to the molecular weight M w via Z eq = M w M e, 7) where M e is the molecular weight between entanglements see Table A.). Motion of a chain along the contour of the tube is unhinere by topological constraints an is known as reptation. The polymer microstructure can be parametrise by a conformation tensor A = RR 3Rg, 8) for en-to-en vector R an raius of gyration R g, which satisfies the Rolie-Poly equation [] DA Dt = K A + A KT A I) τ T, γ) ) ) 3 tra A + β A I), τ R T ) tra 3 where D Dt = t + u ) is the material erivative for flui velocity u, K = αβ u α is the velocity graient tensor an tra enotes the trace of tensor A. The first term in Eq. 9 escribe how chains become stretche an oriente in the flow fiel, whereas the last two terms efine two relaxation mechanisms: reptation along the tube an Rouse relaxation of the tube stretch, respectively. The convective constraint release CCR) mechanism, where the motion of neighbouring tubes can release a topological constraint, is controlle by the parameter β. In fast flow conitions, CCR together with alignment can lea to a flow-inuce ecrease in the entanglement fraction ν = Z/Z eq. Thus, we incorporate flow-inuce isentanglement via the recent kinetic equation of Ianniruberto []: ν t = β K : A ) tra t tra) ν + ν ) T ). 9) 5 Entanglement loss can be moifie by changing β see Ref. []) an entanglements are regaine via curvilinear iffusion along the tube i.e. reptation. At equilibrium, the relaxation times epen on temperature via the typical Williams-Lanel-Ferry WLF) equation [37]: the reptation time governs the orientation of the tube: T ) = τ exp an the Rouse time R stretch: R T ) = τ R exp C T T ) T + C T ), ) governs the relaxation of the tube C T T ) T + C T ). ) Here C an C are the WLF constants, T is the reference temperature an τ an τ R are the reptation an Rouse time at T given by [38] τr = τe Zeq, τ = 3τ e Z 3 eq 3.38 Zeq Z eq.55 Zeq 3 ) 3a), 3b) respectively, where τe is the Rouse time of one entanglement segment at T see Table A. ). To incorporate the anisotropic nature of polymers in flow, the reptation time is moifie accoring to [] τ T, γ) = T ) + β K : A ) tra. 4) tra t where the temperature-epenence of the equilibrium reptation time is given by Eq. ). Thus, polymers that are more aligne an therefore partially isentangle) can relax faster at a rate proportional to CCR. The Rouse time in Eq. 9 oes not epen on the flow. Often the entanglement number is written as Z eq eq 3τ. 5) R However, this relation is strictly only true for Z eq >. Since printing materials typically have fewer entanglements, it is important to acknowlege that from Eq. 3b the reptation time scales as τ = 3 τ e Z 7/ eq, 6) in the range 6 < Z eq < 5, rather than Z 3 eq. This is consistent with experiments that observe a 3.4 scaling [36]. Flow through the nozzle is characterise by the equilibrium mass-average nozzle Weissenberg number by W i = U N R T N). 7) For W i >, the flow rate excees the characteristic relaxation time an can thus inuce significant eparture

6 from the equilibrium polymer configuration. Similarly the equivalent average Rouse Weissenberg number is given by an signifies stretching of the tube. 3. Moelling the Wel Region 3.. Initial Conition W i R = U N R R T N ), 8) Successful weling epens on the how the polymer melt interiffuses an re-entangles across the layer-layer interface. To unerstan the interiffusion process, we first must quantify the microstructure inuce by the extrusion an eposition process. This configuration will provie an initial conition to calculate the evolution of the polymer structure at the wel. During eposition the polymer melt must eform to make the 9 o turn an transform into the elliptical geometry. Ref. [] calculates the eformation A impose by the extrusion an eposition process using the Rolie-Poly moel Sec..3). The moel parameters are T N = 5 o C, W i =, Z eq = 37 an β =.3. In particular, Fig. 4a shows the stretch across the elliptical cross section of the eposit. There is a istinct graient in the stretch from the top to the bottom of the layer, with the bottom half becoming much more stretche ue to the stretching of the free surface uring the 9 o turn. Due to polymer realignment uring the eposition flow, the melt becomes significantly isentangle across the layer Fig. 4b). 3.. Dynamics at the Wel To etermine how the polymer microstructure evolves at the wel after eposition, we solve the moifie Rolie- Poly Eqs. 9 an uner zero flow conitions i.e. K = ). The initial conition is calculate by the proceure in Ref. [] e.g. Fig. 4). We then calculate the relaxation process at the wel sites t p an b p on either sie of the wel line uring the two stages of printing. During stage, t p is at a layer-air interface, whereas uring stage, t p an b p form a layer-layer interface Fig. 5). The temperature epenence of the reptation an Rouse times Eqs. an ) is etermine by the temperature profile at the wel calculate in Sec... The moel parameters are set to T N = 5 o C, Z eq = 37 an W i =, which are typical for polycarbonate printing material. The CCR parameter is set to β =.3 as in Ref. []. 4. Evolution of Disentangle Wel Structure Fig. 5 shows an elliptical representation of tensor A at wel sites t p an b p ; a sphere represents an uneforme polymer at equilibrium, whereas an ellipse represents a stretche an oriente polymer. The grey circles correspon to the equilibrium shape. We iscuss the effects of changing the print spee an the CCR parameter β on the relaxation ynamics in the Appenices Appenix B, Appenix C. Stage Fig. 5a): Layer L p exits the nozzle at temperature T N an is eposite at time t wp ). The leftmost ellipse represents the initial polymer configuration at wel site t p inuce by the eposition process. Since site t p is expose to the air it rapily cools to the ambient temperature T a Fig. ); the temperature of the free surface rops below T g in t F g =.5µs. Thus, the extrusioninuce eformation an corresponing isentanglement fraction o not have time fully relax an a non-equilibrium polymer configuration is locke into the wel site prior to the creation of the layer-layer interface. Stage Fig. 5b): Layer L p exits the nozzle at temperature T N an is eposite on top of the cool layer L p at time t wp). This creates a layer-layer interface between wel sites t p an b p. Again the left-most ellipses represent the initial polymer configuration at time t wp). Each site has a ifferent initial microstructure ue to ifferent egrees of eformation at the top an bottom of the layer uring eposition, as well as the thermal history. In particular, the initial ellipse at t p represents the eformation inuce at the top of a layer that has been frozen in by the cooling of the layer-air interface uring stage. On the other han, the initial ellipse at b p represents the eformation inuce at the bottom of a layer at print temperature T N, which is much larger ue to the greater stretch aroun the outer corner. Since the layer-layer interface cools much slower than the layer-air interface Fig ), the polymer has much longer to relax before T g is reache t F g t W g ). Thus we see relaxation of the polymer at t p. Similarly, the larger eformation at b p also relaxes with a similar temperature evolution, also reaching T g at time t W g. There is still insufficient time for the polymer to fully relax before the onset of the glass transition, so the wel region remain slightly anisotropic at t W g. Fig. 6 shows how the reptation time τ, the tube stretch tra 3, the principle shear component A yz an the entanglement fraction ν at the wel sites t p an b p evolve uring the two stages of printing; note that b p only exists once L p has been eposite uring the secon stage of printing. We iscuss these features in turn. Reptation Time Fig. 6a,b): After eposition of L p, the reptation time τ rapily iverges Fig 6a), since the temperature at the layer-air interface t p rops below T g in less than µs. When the layer-layer interface is forme at time t wp), heat transfer between the two layers causes the polymer at t p to instantaneously become more mobile an the reptation time becomes finite Fig 6b). Despite having a similar temperature evolution, the two wel sites t p an b p have ifferent reptation times ue to the ifferent egree of stretch inuce by the extrusion process Eq. 4). Deformation Relaxation Fig. 6c,,e,f): During stage, Fig. 6c,e shows that the initial tube stretch tra an shear eformation A yz at the free surface t p o not 6

7 8 H H 4 -z.5 -z x R c) tra-3 ν 4.5 -R tra)-3 b) Entanglements, ν ν 8 a) Stretch, tra)-3 -R x R.5.5 z/h.75 Figure 4: a) Initial stretch profile tra an b) initial isentanglement fraction ν inuce by the extrusion process across a printe layer. b) Quantitative plot of the stretch an ν along the z-axis x = ). The moel parameters are TN = 5o C, W i =, Zeq = 37 an β =.3. The polymer structure is highly stretche an isentangle at the bottom compare to the top of the layer. C). relax ue to the iverging reptation time an persist to the secon printing stage. Once Lp is eposite the increase temperature of the wel site allows relaxation of tra an Ayz Fig. 6,f). Since W ir >, linear relaxation of the eformation at both tp an bp oes not apply. Instea, the relaxation process is two-stage; the first relaxation moe is Rouselike an governe by τr. Once the tube length returns to the equilibrium value tra = 3, the usual reptation behaviour prevails an the reptation time τ becomes equivalent for both wel sites Fig. 6b), now epening only on the temperature evolution. Convective constraint release parametrise by β) only contributes to relaxation whilst tra > 3 an therefore only has a small effect on the relaxation ynamics uring the first Rouse relaxation moe see Appenix Appenix C). Fig 6 shows that the stretch of the tube at both tp an bp has sufficient time to relax prior to the onset of the glass transition at time tw g. In contrast, the relaxation of principle shear component Ayz Fig 6f) is arreste by the glass transition. This anisotropy is particularly prominent at site bp. Thus, a non-equilibrium polymer orientation becomes locke into the wel region at tw g so that the structure of the wel is slightly anisotropic. Recovery of Entanglements Fig. 6g,h): Finally, Fig. 6h shows how the entanglement fraction ν at both sites recovers towars unity uring stage accoring to Eq. ). Initially, the evolution of entanglements is governe by the tube stretch; whilst tra > 3, ν remains constant since entanglements cannot be gaine when the tube is retracting. Once the stretch has returne to equilibrium, entanglements recover at a rate etermine by the reptation time τeq, which epens only on the temperature evolution. Due to the arresting glass transition, entanglements are not able to fully recover an ν.5 at tw g for both sites tp an bp. Thus, the wel region is approximately 5% less entangle than the equilibrium material, presumably yieling a lower mechanical strength compare to the equilibrium material. Since β only affects entanglement recovery uring the Rouse relaxation moe, we see see similar ν at tw g for all β see Appenix Appenix 5. Wel Inter-Penetration Thickness We now consier iffusion ynamics at the two wel sites tp an bp locate at istances ±Rg from the wel line Eq. 6). First we consier isotropic iffusion that is spe up by the relaxation of the tube stretch accoring to Eq. 4. We then inclue the effects that an anisotropic environment has on the iffusion irection. Finally, we incorporate inhomogeneous iffusivity across the chain ue to the polymer iffusing into a ifferent eformation environment as it crosses the wel line. 5.. Isotropic Weling Approximation The curvilinear iffusion coefficient along the contour length of the tube is given by [36] Dc = kb T Nζ 9) where N is the number of Kuhn steps along the path, ζ is the friction coefficient an kb is the Boltzman constant. The Kuhn length b is the statistical length of a polymer segment an often represents chain stiffness [36]. The curvilinear istance ` travelle by a polymer chain in time t is then given by the Einstein relation for a oneimensional ranom walk [8] h` i = Dc t = Lc t. τ ) The tube contour length is efine by [36] Lc = N b, at for Kuhn length b an tube iameter p at = b Ne, ) ) where Ne is the number of Kuhn steps between entanglements; at equilibrium, Ne = N/Zeq. 7

8 a) Stage : Printing Lp- b) Stage : Printing Lp nozzle air z y nozzle Lp- twp-) air printe layer time tp- Lp Lp- twp) printe layers W tg time bp tp- Figure 5: Elliptical representation of polymer conformation tensor A in the yz-plane uring two stages of printing. a) Stage : Layer L p is eposite at time t wp ) an weling site t p at the layer-air interface rapily cools below T g an so that the eformation inuce by extrusion only relaxes slightly. b) Stage : Layer L p is eposite at time t wp) creating layer-layer interface between wel sites t p an b p an weling begins. Polymer relaxes until temperature of the wel reaches the glass transition at time t W g. Moel parameters are T N = 5 o C, W i =, Z eq = 37 an β =.3. The structure of the wel is anisotropic at t W g. τ tra)-3 A yz t p- a) b p tg F Stage : relaxation at layer-air interface c) e) Stage : relaxation at layer-layer interface W t g b) ) f) - g) h) ν t-t wp-) s) -4 - t-t wp) s) Figure 6: Time evolution of a,b) the reptation time τ given by Eq. 4) c,) the tube stretch tra 3, e,f) the principle shear component A yz, an g,h) the entanglement fraction ν at wel sites t p an b p uring two stages of printing. Dashe lines inicate the time at which the glass transition occurs at the free surface t F g.5µs an at the wel t W g s. Moel parameters are T N = 5 o C, W i =, Z eq = 37 an β =.3. The structure of the wel is anisotropic an isentangle at t W g ν.5). 8

9 a) Wel Site b p 3 χ χ zz t W g b) Wel Site t p- χ χ zz t W g R E χ/r g R g -4 - t-t wp) s) -4 - t-t wp) s) Figure 7: Inter-penetration istance calculate at wel sites a) b p an b) t p. Isotropic weling given by χ Eq. 4) an anisotropic weling given by χ zz Eq. 6). Dashe lines inicate R g an R E an t W g s. Moel parameters are T N = 5 o C, W i =, Z eq = 37 an β =.3. Inter-iffusion across the wel in z-irection) is asymmetric an slowe by the anisotropic environment inuce by extrusion. The D tube contour executes a ranom walk in 3Dspace, leaing to an interpenetration istance χ given by [8, ] χ = l / a T, ) / t 3) = Nb τ Thus, a polymer chain iffuses via a ouble ranom walk process. In the FFF problem the ecreasing temperature progressively slows the motion. Thus, the interpenetration istance is calculate via the integral ) χ t W /4 g = 36 R g t w τ T t), γt)) t, 4) where R g = Nb / 6 an the reptation time epens on both temperature an shear rate Eq. 4). For polycarbonate printe at T N = 5 o C, Fig. 7 shows the interpenetration epth for polymers locate at the wel: polymers locate at b p travel slightly further than polymers locate at t p ue to a smaller reptation time, which is a result of the increase tube stretch at this wel site. In particular, for our moel parameters χ R g at b p before the glass transition arrests iffusion at t W g. Usually, experiments fin that bulk strength is achieve once iffusion of the orer of R g has occurre [3] an molecular simulations suggest that only a few entanglement lengths are require []. Despite a iffusion istance greater than R g, we fin that orientations o not fully relax uring this time, thus creating an anisotropic wel structure Fig. 5). The polymer must iffuse its en-to-en istance R E = 6R g to fully escape its tube an relax to equilibrium Anisotropic Weling Approximation Eq. 4 oes not account for a preferre iffusion irection ue to the local anisotropic structure of the polymer melt. Similar to the work of Ilg & Kroger [6], we propose a time-epenent anisotropic iffusion tensor of the form D = D I + ηa I)), 5) where D = Nb /τ t). The anisotropy parameter was foun to be η /3 by molecular simulations of long chains [6]. The interpenetration istance across the interface is then given by χ zz R g = 36 ) t W /4 g t w τ t) + ηa zzt) )) t. 6) For polycarbonate printe at T N = 5 o C, Fig. 7 shows how the interpenetration epth across the wel line is reuce for anisotropic iffusion of polymers locate at wel sites t p an b p. Yet χ zz is also larger than R g, which accoring to molecular ynamics simulations [] suggests that bulk strength shoul be achieve in the wel region Mutual Diffusion Approximation When a polymer molecule iffuses across the interface, it is not only affecte by the anisotropy of its own structure, but also the anisotropy of the environment into which it iffuses. We consier polymers locate either sie of the interface of type t an b to signify the eformation at wel sites t p an b p, respectively. Similar to the work of Kramer et al. [9], we propose a mutual iffusion coefficient of the form D M t, z) = φt, z))d b t) + φt, z)d t t), 7)

10 a) φ φ φ t p- t p- t= -5 b p t= -3 b p t= - t p- b p - z/r g t p- t p- t= -4 b p t= - b p t p- b p - z/r g with/r g b) 6 4 t= with left with right with t W g -4 - t-t wp) s) Figure 8: a) Evolution of the volume fraction φ profile Eq. 9) in the region of the layer-layer interface z = ) an wel sites t p an b p. Arrows inicate the interfacial with on the left re) an on the right green) of the interface. b) Evolution of the interfacial with, an withs calculate to the left an right of the interface. Dashe line inicates R g an t W g s. Moel parameters are T N = 5 o C, W i =, Z eq = 37 an β =.3. The wel thickness is asymmetric. where D b t) = Nb τ b t) I + ηa b t) I) ), D t t) = Nb τ t t) I + ηat t) I)). 8a) 8b) The volume fraction occupie by type t polymers is enote by φ. Initially φ = at wel site t p an φ = at wel site b p. The volume fraction evolves accoring to φ t = z D M zzt, z) φ z ), 9) where Dzz M is the zz-component of the mutual iffusion tensor governe by Eq. 7. In this way a iffusing chain carries its mobility across the wel line an the iffusion coefficient epens on the local composition of mobility. For polycarbonate printe at T N = 5 o C, Fig. 8a shows the evolution of φ uring the relaxation process. The wel is forme ue to polymers iffusing across the interface. Assuming that a wel is forme in the region % < φ < 98%, Fig. 8b. shows the evolution of the interfacial with of the wel region. We highlight the asymmetric nature of iffusion across the interface by plotting the with calculate to the left an the right of the interface. Asymmetry arises ue to the ifferent egree of eformation at wel sites t p an b p. The Kramer moel [9] escribes the mutual iffusion of two ifferent molecular weight polymers [3]. However, since the FFF problem iscusse here involves a single molecular weight inter-iffusing between ifferent-mobility environments, a polymer chain will inherit the relaxation characteristics of the environment into which it is iffusing. Thus, a similar approach can be taken to calculate the mutual interpenetration istance of a molecule locate at wel site t p, namely χ M. That is χ M R g = 36 t W g t w DM zz t)t ) /4, 3) where the iffusion coefficient of this molecule is given by χ ) M D t + χ M D b, χ M < R g, D M t) = R g R g D b, χ M R g, 3) an χ M parametrises how far the polymer has penetrate. D t an D b are given by Eq. 8. A similar relation is use for molecules locate at wel site b p. In this way the mutual iffusion coefficient of a chain is given by a weighte average of the two environments on either sie of the interface that epens on the polymer s location. Compare to χ zz Eq. 6), Fig. 9a shows how iffusion of molecules initially locate at wel site b p is slowe own by iffusing into a slower moving environment for polycarbonate printe at T N = 5 o C. In contrast, the iffusion of molecules initially locate at t p is increase by iffusing into a faster-moving environment Fig. 9b). For this mutual iffusion moel, type t polymers ultimately iffuse further across the interface ue to the increase mobility of iffusing into a faster environment. Thus, espite type b polymers initially being more mobile ue to the egree of eformation, the type t polymers ultimately create a thicker interfacial with. This approach gives a final asymmetric interfacial with similar to that seen in

11 3 a) Wel Site b p t W g b) Wel Site t p- t W g R E χ/r g χ zz χ M with χ zz χ M with R g -4 - t-t wp) s) -4 - t-t wp) s) Figure 9: Time evolution of the interpenetration istance χ zz incorporating anisotropic iffusion Eq. 6), an the interpenetration istance χ M given by the mutual iffusion theory Eq. 3) for polymers initially locate at wel sites a) b p an b) t p. Dashe lines inicate R g an R E an t W g s. Distances are compare to the the interfacial with etermine by φ shown in Fig. 8b. Moel parameters aret N = 5 o C, W i =, Z eq = 37 an β =.3. The wel thickness is asymmetric, with material iffusing further into more mobile environment i.e. t p towars b p). Fig. 8b for the Kramer moel base on local composition Eq. 9). Incluing the effects of mutual iffusion also yiels χ M > R g at this print temperature. Thus, for this case the mechanical strength of the wel region is limite the molecular structure an entanglement fraction) of the wel region itself Fig 6), rather than the interpenetration epth of the polymer. 6. Controlling Wel Strength 6.. Final Wel Properties The mechanical properties of a material can be characterise by the plateau moulus G e [36] an the fracture toughness G c [39, 4]. Both properties are proportional to the molecular weight between entanglements, M e : G e M e, 3a) G c M ) e, 3b) qm w for q.6 [39]. Smaller M e results in a greater entanglement ensity an therefore increases both the bulk moulus an fracture toughness of a material. For polycarbonate, the entanglement molecular weight in etermine by M e =.56a T g/mol, 33) where a T = 37.9 Å is the tube iameter Eq. ) an the pre-factor accounts for bon angle, characteristic ratio an monomer weight [4]. At a wele interface the mechanical strength is etermine by how many segments of length a T cross the interface [4]. Thus, the mechanical strength of a wel is attribute to the wel thickness χ W, as well as the integrity of the entanglement network at the wel ν W. Bulk strength is expecte for ν W = an χ W /R g >. For certain printing conitions we have seen that, although the wel thickness excees R g, the ultimate structure of the wel is anisotropic with a weaker entanglement structure. Thus, reliable preiction of the wel properties from material properties an printing parameters is key to ensuring wel strength an avancing FFF technology. The FFF extrusion moel in [] assumes an iealise eposition process uring which there is no cooling or relaxation of the melt. In aition, the assume temperature profile in the nozzle neglects inhomogeneities ue to thermal iffusion an shear heating effects. Both of these factors may significantly affect the temperature an polymer structure at the wel site, an consequently affect the wel properties. Thus, it is important to valiate χ W against experimentally measure wel thicknesses, which can be challenging to measure reproucibly an reliably from FFF-printe parts. In this section we raw qualitative conclusions of the wel properties base on our moel. First we calculate the final wel entanglement number ν W an the interpenetration epth χ W at time t W g an examine how they vary with the equilibrium mass-average nozzle Weissenberg number Eq. 7), print temperature T N an molecular weight Z eq. Finally, we suggest how these wel properties can increase the entanglement molecular weight in the wel region an therefore ecrease the mechanical strength of a printe part.

12 . fully entangle a) 6 b) ν W.8.4 Z eq =37 Z eq =8 Z eq = -4 - Wi χ W /R g 4 anisotropic isotropic R E R g -4 - Wi. fully entangle c) 6 ) ν W T N χ W /R g 4 anisotropic isotropic R E R g T N Figure : Wel characteristics a,c) final wel entanglement fraction ν W an b,) final wel thickness χ W measure at t W g for a,b) a range of Weissenberg numbers W i at T N = 5 o C an c,) a range of print temperatures T N at the slow print spee U N = mm/s. For χ W < R E = 6R g the wel structure will be anisotropic. Equilibrium entanglement numbers are Z eq = 37, 8 an, similar to typical printing materials. 6.. Effect of Print Spee For a fixe print temperature T N an a range of typical entanglement numbers Z eq, Figs. a an b show wel properties ν W an χ W for increasing W i. We observe a slight ecrease in ν W with W i for the largest molecular weight Z eq = 37, whereas χ W is inepenent of the nozzle shear rate for all three Z eq shown Fig. b). Thus, the increase initial stretch impose by larger shear rates in the nozzle is not large enough to influence the wel thickness by reucing the reptation time. Melts with fewer entanglements are more mobile an are therefore able to iffuse further uring the weling process, creating a thicker, more entangle wel structure. For χ W < R E the wel structure will be anisotropic at the glass transition. Once the interpenetration epth surpasses R E we expect an isotropic wel structure Effect of Print Temperature In contrast to Weissenberg number, we fin that print temperature significantly affects weling behaviour Figs. c,). Both ν W an χ W increase as a function of T N since higher temperatures significantly spee up iffusion. Notably, even at low print temperature T N = o C the wel thickness continues to surpass R g, suggesting healing of the interface. However, the wel site remains almost fully isentangle in this case Fig. c), so that bulk strength is unattainable. This is a consequence of the elay in entanglement recovery ue to the relaxation of the tube stretch Fig. 6h). Moreover, since smaller molecular weights have shorter Rouse times, entanglement recovery begins earlier, which allows the formation of a more entangle wel. Only for T N 3 o C oes the wel becomes fully entangle for all Z eq Effect of Molecular Weight For fixe print spee U N = mm/s, Figs. a,b show how ν W an χ W vary with a broaer range of molecular weights. Since entanglements only recover once the stretch has relaxe, the final wel entanglement can be preicte from Eq. as follows: νw ν t = ν ep ν ν = ν T t)), 34a) t W g t w T t))t, 34b)

13 a) 6 b) T N =3 T N =5 T N = Eq.38 Eq.38 Eq.38 ν W. χ W /R g 8 4 T N =3 T N =5 T N = Eq.35 Eq.35 Eq.35 Eq Z eq Z eq Figure : Log-log plot of wel characteristics at t W g for a broa range of entanglement numbers Z eq, fixe print spee U N = mm/s an three print temperatures T N =, 5 an 3 o C: a) wel entanglement ν W with lines given Eq. 35. For T N = 5 o C we compare Eq. 35 to Eq. 37. b) The final wel thickness χ W with lines given by Eq. 38. The entanglement number for typical printing materials is in the range Z eq 4. where ν ep is the entanglement fraction at the wel site after eposition. Eq. 34 yiels ) ν W = ν ep Z eq )) exp t W g t w T t))t 35) weight at which iffusion is arreste uring relaxation of the tube stretch. In this regime the tube ynamics reuce to the reptation time via Eq. 4 an spee up iffusion. This effect results in a larger χ W than preicte by Eq. 38. an is plotte in Fig. a. The integral in Eq. 35 can be approximate as t W g t w T t))t C t W g T t, 36) N) t w where the constant C accounts for being approximately two orers of magnitue larger when entanglement recovery begins ue to cooling Fig. 6b). Thus, the final wel entanglement can be written as ) ν W ν ep Z eq )) exp C tw g T N), 37) an is also plotte in Fig. a for T N = 5 o C an C =.6. Z eq has two effects on entanglement recovery at the wel. First, larger Z eq increases the time scale so that larger molecular weights iffuse slower an therefore recover fewer entanglements. Secon, at a fixe print spee increasing Z eq yiels a larger Weissenberg number, which leas to greater isentanglement uring eposition i.e. smaller ν ep []. The wel interpenetration epth scales as χ W R g ) /4 T t)), 38) Zeq 7/8 since scales as Zeq 7/ see Eq. 6). The preicte χ W eviates from Eq. 38 at some Z eq epening on print temperature Fig. b). This critical Z eq efines the molecular Wel Fracture Toughness Mechanical strength at the wel is etermine in part by the molecular weight between entanglements Eq. 3). In the wel region we have M W e = M e ν W. 39) Thus, from Eq. 3b the fracture toughness at the wel is given by G W c M e W ) qm w qν W Z eq ). 4) Hence ecreasing ν W increases Me W an reuces the toughness of the wel. The equilibrium bulk strength G c is only achieve for ν W =. For polycarbonate, Fig. shows that G W c can ecrease by 5% epening on the entanglement number Z eq an the print temperature T N. In particular, for a prescribe print temperature, there exists a maximum molecular weight that can be printe whilst maintaining G W c = G c. For example, Z eq 4 can be printe at T N = 3 o C before strength is affecte by isentanglement, whereas only Z eq maintains bulk strength at lower temperature T N = o C. 7. Conclusion We have evelope a moel for the non-isothermal FFF weling process to test the effect of changing print spee,

14 G c wel /Gc bulk strength T N =3 T N =5 T N = Z eq Eq.4 Eq.4 Eq.4 theory against molecular ynamics simulations of isentangle melts, in particular to verify that ν < for χ > R E is a result of the initial elay in entanglement recovery ue to tube stretch. Practically, thicker an more entangle wels can be forme by increasing the print temperature or using a less entangle printing material, since both parameters significantly reuce the reptation time. We fin that the wel thickness is inepenent of the eformation inuce by ifferent nozzle shear rates, as is the final entanglement network for Z eq =. Consequently equivalent mechanical integrity is expecte across all print spees. Thus, the maximum print spee available can be exploite to increase prouctivity. Figure : Log-log plot of the fracture toughness at the wel G W c normalise by the equilibrium fracture toughness G c for fixe print spee U N = mm/s an three print temperatures T N =, 5 an 3 o C. The entanglement number for typical printing materials is in the range Z eq 4. Lines are given by Eq. 4. temperature an entanglement number on the ultimate weling characteristics of an non-crystalline polymer melt. It has previously been shown that the extrusion process can significantly eform an isentangle the polymer microstructure prior to weling []. After eposition the printe layer cools an this eformation relaxes via reptation, whilst inter-iffusing with the previously-printe layer. The temperature profile at the wel between the two layers is calculate by solving the heat equation in one imension. The cooling rate inhibits the total relaxation of the eformation inuce by printing, so that the ultimate structure of the wel is anisotropic an less entangle than the equilibrium material for typical printing conitions. Solving a iffusion equation that incorporates anisotropic an mutual iffusion yiels the thickness of the wel forme between the layers. The moel preicts that the wel thickness typically surpasses R g, but not quite enough to fully relax. Thus, mechanical strength shoul not be limite by interpenetration epth. However, espite sufficient wel thickness for bulk strength at the interface, entanglements o not have sufficient time to recover uring cooling; ν W is as low as 5% for typical printing conitions. Since a isentangle wel structure can significantly increase the entanglement molecular weight in the wel region, the mechanical properties at the wel may be significantly reuce. These finings suggest that isentanglement in the nozzle combine with the elay in entanglement recovery ue to the relaxation of the tube stretch is the key mechanism responsible for a reuce mechanical strength in the wel region. Although the theory of flow-inuce isentanglement has been compare to molecular ynamics simulations of planar shear flow [43], the way in which polymers recover entanglements in the event of flow cessation is yet to be aresse. It is crucial to benchmark this re-entanglement 4 Acknowlegements We thank Jonathan Seppala an Kalman Migler for avice an an enjoyable collaboration, as well as the National Institute for Stanars an Technology NIST), Georgetown University, an the Ives Founation for funing. Appenix A. Extrusion an Deposition Moel Here we summarise the FFF extrusion an eposition moel etaile in Ref. []. The printing material is heate to temperature T N an extrue through a nozzle of raius R at mass-average spee U N. Assuming a steay state the momentum balance is given by σ =, A.) for stress tensor σ. The total stress in the polymer melt comprises solvent an polymer contributions σ = pi + G e A I) + µ s K + K T ), A.) where p is the isotropic pressure an G e is the plateau moulus. For times shorter than τ e, Rouse moes corresponing to lengths shorter than M e contribute to a backgroun viscosity efine as [44] µ s = π G e R Z. eq A.3) The temperature profile is assume to be uniform across the nozzle an Eq. A. is solve alongsie the Rolie-Poly Eqs. 9) an ) to calculate the plug-like velocity profile, an the polymer eformation an isentanglement across the nozzle. Nozzle Weissenberg numbers see Eqs. 7) an 8)) for typical fast an slow print spees are quote in Table A.; the polymer is foun to stretch an orient in the nozzle epening on print spee. The material is then eposite into a layer of thickness H, which travels horizontally at spee U L, in a frame moving with the nozzle. During this eposition the material must spee up an eform to make a 9 o turn an

15 Table A.: Material parameters for a typical amorphous printing material, polycarbonate. Polycarbonate Properties Notation Value Units Glass-Transition Temperature T g 4 o C Thermal Diffusivity [34] at 5 o C) α.4 mm /s Eq. 3 α.6 mm /s Eq. 3 B.3 3 o C Molecular Weight M w 6 kda Entanglement Molecular Weight [4] M e.6 kda Plateau Moulus [4] G e.6 6 Pa Entanglement Time [36] at T ) τe s Reference Temperature T 6 o C WLF parameter C 3 - WLF parameter C 6 - Equilibrium Entanglement Number Z eq 37 - CCR parameter β.3 - Equilibrium Reptation Time at T N ).33 s Equilibrium Rouse Time at T N ) R.57 s Table A.: Moel printing parameters for two typical print spees; a fast an a slow case. Printing Parameters Notation Fast Case Slow Case Units Extrusion Temperature T N 5 5 o C Ambient Temperature T a o C Mean Extrusion Spee U N.75. m/s Mean Print Spee U L..3 m/s Nozzle Raius R.. mm Layer Thickness H.3.3 mm Reptation Weissenberg Number average) W i 3 - Rouse Weissenberg Number average) W i R transform from circular to elliptical geometry since typical H < R). In orer to conserve mass, U L > U N Eq. ). Rather than calculating the full flui mechanics, the shape of the eposition is prescribe ignoring ie swell) an the flui is avecte using local flux conservation. We assume a uniform temperature profile through out the eposition an neglect polymer relaxation uring this stage. Due to the assumption of zero polymer relaxation, the moel breaks own for slower print spees such that W i < H/R []. In this way the velocity profile u is prescribe only by the geometry of the eposition an the Rolie-Poly Eq. 9 reuces to Appenix u )A = K A + A + K T. A.4) B. Wel Structure an Print Spee Fig. B. shows the relaxation ynamics at the wel site b p for the two typical print spees given in Table A., corresponing to W i = an 3. Parameters are T N = 5 o, Z eq = 37 an β =.3. Time is scale by T N) in orer to highlight the effect of the cooling temperature profile an the case τ = with an initial polymer configuration impose by the fast print spee) is plotte to high- 5 light the effect that polymer stretch has on the relaxation process. The reptation time τ T, γ) an the Rouse time τ R T ) are plotte in Figs. B.a,b. Although the stretch inuce by printing significantly reuces the reptation time compare to Eq. 4), the ifference in the reptation for the two typical print spees is small. The Rouse time is inepenent of print spee Eq. ). Figs. B.c, shows the relaxation of the principle shear component A yz an the the tube stretch tra 3. The principle shear component A yz relaxes on a time scale t/ T N ) >, emonstrating how the cooling temperature profile inhibits the relaxation process. We see a similar two-moe relaxation Rouse followe by reptation) for both print spees. The stretch has sufficient time to relax prior to the glass transition but the polymer orientation remains out of equilibrium for both print spees. The structure is slightly closer to equilibrium for the fast-printing case ue to a slightly smaller reptation time at t wp). Fig. B.e shows how the entanglement fraction evolves at the wel for the two print cases. The wel region is approximately 5% less entangle than the bulk material an the final ν is similar for both printing spees. Finally, Fig. B.f shows the isotropic interpenetration istance χ calculate by Eq. 4). Since the reptation time is