Experiments and Finite Element (FE) calculations were performed to study the raster angle dependent fracture behaviour of acrylonitrile butadiene styrene (ABS) t hermoplastic processed via fused filament fabrication (FFF) additive manufacturing (AM). The fracture properties of 3D printed ABS were characterized based on the concept of essential work of fracture (EWF), utilizing double edge notched tension (DENT) spe cimens considering rectilinear infill patterns with different raster angles (0°, 90° and +45/ -45°). The measurements show ed that the resistance to fracture initiation of 3D printed ABS specimens is substantially higher for the printing direction perpendic ular to the crack plane (0° raster angle) as compared to that of the samples wherein the printing direction is parallel to the crack (90° raster angle), reporting EWF values of 7.24 kJ m-2 and 3.61 kJ m-2, respectively. A relatively high EWF value was also reported for the specimens with +45/ -45° raster angle (7.40 kJ m-2). Strain field analysis performed via digital image * Corresponding author : s.kumar@eng.oxon.org 2 correlation showed that connected plastic zones existed in the ligaments of the DENT specimens prior to the onset of fracture, and this was corroborated by SEM fractography which showed that fracture proceeded by a ductile mechanism involving void growth and coalescence followed by drawing and ductile tearing of fibrils. It was further shown that the raster angle dependent strength and fra cture properties of 3D printed ABS can be predicted with an acceptable accuracy by a relatively simple FE model considering the anisotropic elasticity and failure properties of FFF specimens. The findings of this study offer guidelines for fracture -resistant design of AM -enabled thermoplastics. Keywords: Additive manufacturing; 3D printing; fused filament fabrication; fracture toughness; thermoplastics 1. Introduction Additive manufacturing (AM), commonly known as 3D printing, has attracted the interest of the academic and industrial research community because of its ability to fabricate objects with complex geometries at relatively low cost and with high flexibility [1-9]. It is an emerging technology where structures are printed layer by layer with the help of computer aided design (CAD) models. Popular AM techniques include Fused Filament Fab rication (FFF) [4, 10 -12][13], Selective Laser Sintering (SLS) [14, 15] , Digital Light Processing (DLP) [16, 17] [18], and Stereolithography (SLA) [19-21][22-24]. FFF has several advantages including simplicity of operation, lower cost of fabrication compared to conventional methods, for e. g. injection moulding, fast processing and the ability to process a versatile range of thermoplastics with complex structures [25-30]. The feedstock used in FFF is usually a thermoplastic polymer -based filament of a specific diameter . As the printing occurs in a layer by layer fashion by local welding of discrete extruded layers , the printing direction (aka raster angle) and the choice of the infill pattern play an important role in the mechanical response of 3D printed thermoplastic structures processed via FFF . The layered structure of the FFF print forms discrete bonds between adjacent beads and typically includes a significant amount of interlayer porosity, which affects the mechanical properties of the printed parts [31], and results in build orientation and raster angle dependent fracture properties. 3 Several phenomenological studies have been reported in the literature to foster understanding of the process -property relationships associated with 3D printed thermoplastics [32-34]. Aliheidari et al. [35] studied the effect of process parameters (melt temperature, bed temperature, layer height and layer width) on the fracture resistance of acrylonitrile butadiene styrene (ABS) and found that the critical J-integral is strongly correlated with the process parameters through both the interlayer adhesion and the mesostructure. Sood et al. [36] used a surface optimization method and found that the tensile, flexural and impact strength of 3D printed ABS samples can be enhanced by altering the printing parameters such as layer height, width, air gap and orientation. Tymark et al. [37] optimized the layer height and printing orientation for enhanced tensile strength, strain and modulus of polylactic acid (PLA). Durgun et al. [38] reported that the out -of-plane orientation has a greater influence on the mechanical properties of 3D printed parts than in -plane orientation of layers. Gurrla et al. [39] suggested a simplified mathematical model to understand the neck growth between cylindrical filament layers. By using ultimate tensile strength values obtained for two different print directions, they established the simulation output with experimental data and concluded that neck growth between layers has a dominant effect on the strength of printed parts. Yang et al. [40] studied the effect of thermal processing conditions in FFF o n crystallinity and mechanical properties of PEEK samples in order to develop the relation among them. Arif et a l. [41, 42] evaluated the influence of the print -orientation and raster angle on the tensile, flexural , thermal and fracture behaviour of FFF PEEK . McLouth et al. [43] studied the effect of print orientation and raster pattern on fracture toughness of ABS matrix and analysed variations of the fracture toughness with mesostructure. They concluded that, as alignment of printing filament changed from 0 to 90° (parallel to perpendicular) with respect to crack plane, the fracture toughness increased by 54%. It was also observed that the printing orientation significantly affect the frac ture toughness. A decrease of 11% was observed in fracture toughness for printing orientation of 0/90° compared to +45/ -45°. Recently, Cuesta et al. [44] used t he essential work of fracture (EWF) approach to examine the mechanical behaviour and fracture properties of 3D printed polymeric materials, including fiber -reinforced plastics, and proposed a new miniature test specimen for measuring the fracture propertie s of polymers processed via AM. Moreover , Lorenzo -Bañuelos et al. [45] examined the effect of raster orientation on the fracture p roperties of thin polypropylene (PP) components processed via material extrusion AM, and found that the 4 printing orientation can have a significant effect on the essential work of fracture of the 3D printed PP. As the applications of FFF processes are exp anding, good knowledge and understanding of the relationship between the printing orientation and the mechanical performance is becoming more important from the process and product design point of view to achieve high quality prints. Although previous stud ies have provided useful guidelines on the process -property relationship of FFF enabled parts, there is a need to combine experiments and modelling to understand how the printing direction and build orientation affect the crack initiation resistance of 3D printed thermoplastics that undergo significant inelastic deformation before the onset of fracture. Here in, we study the fracture properties of ABS thermoplastics processed via FFF using the concept of essential work of fracture (EWF) [46-48]. The effect of the raster angle on the EWF and non - essential work of fracture (non -EWF) of 3D printed ABS are assessed via both experimental and numerical approaches. Tensile tests are performed on double edge notched tension (DENT) specimens printed with different raster angles (0°, 90° and +45/ -45°), and the essential and non - essential work of fracture is determined from the obtained measurements. Moreover, the evolution of the strain fields in the ligaments of the specimens are exam ined via digital image correlation (DIC), and scanning electron microscopy (SEM) images are obtained from the fracture surfaces to identify dominant failure mechanisms and correlate them with the results extracted from the macroscopic fracture tests. The e xperimental results are further corroborated with the predictions from a 2D plane -stress FE model, considering the anisotropic elastic response and failur e properties associated with the FFF print. 2. Experimental Methods 2.1 Sample fabrication via Fused Filament Fabrication 3D printing Acrylonitrile butadiene styrene (ABS) is a rubber toughened, amorphous polymer with exceptional printability. ABS filament feedstocks (diameter 1.7 5 mm) supplied by leapfrog MAXX PRO, Netherlands with a batch number A -22-003 was used for the preparation of double edge notch tension (DENT) specimens via FFF 3D printing [49] As shown in Fig . 1, CAD models of double DENT specimens (length h = 90 mm, width W = 30 mm, thickness t = ~1 mm) were prepared in 3D printable STL format and a slicing tool (Simplify3D) was used to convert them in to a G -code 5 which the 3D printer (Creator Pro Flash Forge ) can process and execute prints. Four different ligament lengths ( L = 4, 6, 8 or 12 mm) were considered, and for each choice of L, three different specimens were printed with raster angles of 0°, 90° and +45/ -45° (see Fig. 1), respectively, giving a total of 12 different specimens. The process parameters such as nozzle tip temperature (230 °C), bed temperature (105 °C), layer height (0.17 mm), infill density (100%), extrusion width (0.4 mm) and number of layers (6) were kept constant for all the samples [50] . To examine possible changes in the crystallinity of the polymer as a result of the 3D printing process, wide angle X -ray diffraction (WAXD) analysis was performed on ABS samples before and after 3D printing (see Section S1, Supplementary Inform ation). 2.2 Measurement of the elastic constants of 3D printed ABS It is known that the mechanical properties of components processed via FFF are anisotropic, and are dependent on the printing direction. Since all infill patterns in this work are rectilinear (see Fig. 1), we consider our 3D printed ABS as a transversely isotropic solid composed of vertically stacked layers of fused beads. To determine their elastic properties, we performed quasi -static tensile tests on dogbone specimens printed wit h the same raster angles as the DENT specimens (0o, 900 and +45/-45o), as described in Section S2 (Supplementary Information). The engineering stress vs. strain responses for specimens with the aforementioned raster angles are shown in Fig. S2 (Supplementary Information). The in -plane Young’s moduli with respect to the pri nting direction, E11, and the transverse direction, E22, were found by measuring the initial slopes of the stress -strain responses associated with the 0° and 90° samples, respectively, while the in -plane Poisson’s ratio, 12 22 11 /  =− , was mea sured from the slope of the transverse strain, 22 , vs. longitudinal strain, 11 , plot which was constructed using the DIC results for the 0° sample. The in -plane shear modulus, 12 12 12 / G= , was found by measuring the slope of the shear stress vs. shear strain responses deduced from the tensile test results for the +45/ -45o sample via 122P A= and 12 yy xx  =− (1) where, P is the applied load, A is the cross sectional area, and yy and xx represent, respectively, the strains in longitudinal and transverse direction of the specimen. Since the elastic properties in 6 the 2 -3 plane are considered isotropic, we may write 22 33EE= , 12 13GG= and 12 13= , while the Poison’s ratio 23 and shear modulus 23G are given by 21 23 32 12 121 1−==− and ()22 23 23 21EG=+ (2) where 21 12 22 11 /EE= . The elastic constants for the transversely isotropic 3D printed ABS are listed in Table 1. Table 1 Elastic constants for the transversely isotropic 3D -printed ABS 𝑬𝟏𝟏 (MPa) 𝑬𝟐𝟐 (MPa) 𝑬𝟑𝟑 (MPa) 𝑮𝟏𝟐 (MPa) 𝑮𝟏𝟑 (MPa) 𝑮𝟐𝟑 (MPa) 𝝂𝟏𝟐 𝝂𝟏𝟑 𝝂𝟐𝟑 2038 1787 1787 634 634 641 0.366 0.366 0.392 2.3 Essential Work of Fracture (EWF) assessment The EWF concept has been widely used to describe the fracture behaviour of ductile materials [51][48, 52, 53] and is employed here to quantify directional changes in the fracture properties of the 3D printed ABS. The EWF concept is based on the idea that the fracture process zone (FPZ) within the ligament of a DENT specimen can be divided into (i) the inner fracture process zone (IFPZ), in which energy is absorbed through the separation of atomic bonds and the creation of new surfaces, and (ii) the outer plastic deformation zone (OPDZ), where energy is absorbed primarily through plastic deformation processes. Thus, the total work of fracture, Wf, can be additively decomposed as f e pW W W=+ (3) where We is the essential work of fracture, representing the work needed to advance the crack through the material in the IFPZ, while Wp denotes the non -essential work of fracture associated with plastic deformation within the OPDZ. Assuming that We is proportional to the cross -sectional area of the ligament, AL = Lt, and that Wp scales with the volume of the material undergoing plastic 7 deformation, the above equation can be re -written in terms of the fracture energy or specific work of fracture (ene rgy per unit area), wf , as follows f e pw w w L=+ (4) where we is the fracture initiation energy, i.e. essential work of fracture (EWF), wp is the specific plastic work (energy per unit volume) required to yield the material around the ligament, i.e. non - essential work of fracture (non -EWF), L is the ligament length a nd β is a plastic zone shape factor. To determine wf and βwp experimentally, multiple fracture tests were performed on DENT specimens with different ligament lengths ( L = 4, 6, 8, and 12 mm) and raster angles (0°, 90° and +45/-45°) using a Zwick -Roell uni versal testing machine (2.5 kN load cell) with a constant crosshead speed of 1 mm/min. Five repeat tests for each sample were conducted to generate statistically consistent results. To enable strain field measurements via DIC, all DENT specimens were spray -coated with solvent -dispersed black dye prior to testing. Images of the specimen’s surface were recorded during the test using a CCD camera of 5.0 MP at 2 Hz. The obtained images were analysed using VIC -2D software and the nominal average strains in the loading direction, yy , and transverse direction, xx , were evaluated. Note that the average strains were measured in an area within the gauge length of the sample over a width slightly less than the specimen width to reduce the error in strain calculation at the boundaries. The total work of fracture, Wf , was evaluated for each specimen by calculating the area under the measured load -displacement curve. Then, the specific work of fracture, wf , was calculated by dividing Wf by the respective ligament area AL, and the obtained data was used to construct a wf vs. L plot. Since the EWF theory predicts a linear relationship between wf and L (see eq. (4)), the EWF and non -EWF can be determined by fitting a regression line to the wf vs. L plot, where the intercept and slope represent the we and βwp, respectively. Furthe rmore, t he fracture surface morphologies of the DENT specimens were characterized through scanning electron microscopy (SEM) using a Nova Nano SEM (50 series) operated at 10 kV. Each sample was gold sputter coated prior to the SEM examination. 3. Finite E lement Modelling To better understand the experimental results and to develop a simple tool for strength prediction of notched samples, finite element (FE) models were developed in ABAQUS 6.14. The 3D printed ABS is treated as a transversely isotropic solid with elastic pr operties chosen according to the 8 measurements given in Table 1. The criteria for damage initiation in the 3D printed ABS are based on Hashin’s theory [54], which was originally developed for fiber reinforced c omposite materials but can also be used for analysis of orthotropic materials [55] by using the definitions of the directions as described in equations (5 -8). In this model, the print direction is treated analogous to the fibre direction and the transverse to the print direction is treated as the transverse directi on. The Hashin damage theory is available in ABAQUS, and introduces four damage initiation criteria associated with different damage mechanisms [54]: Axial direction tension: 22 11 12 1T 12fTF   =+        for 110 (5) Axial direction compression: 2 11 1CfCF =  for 110 (6) Transverse direction tension: 2 11 1CmTF =  for 220 (7) Transverse direction compression: 22 2 2 22 22 12 13 13 2 12122C mC CF                  = + − +              for 220 (8) Here, 11 , 22 , 12 and 13 are the stress components of the effective stress tensor with respect to the material coordinate system (1,2,3) where 1 denotes the printing (or longitudinal) direction, 2 the transverse direction and 3 the through -thickness direction. The material streng th of the transversely isotropic material is represented by 1T , 1C , 2T , 2C , 12 and 13 where the subscripts ‘T’ a nd ‘C’ denote tension and compression, respectively, and α in Eq. (5) is a parameter that quantifies the effect of the in -plane shear stress on the damage initiation criterion for fiber tension. The values for the tensile strength in the longitudinal, and transverse direction were taken from the stress -strain curves (see Fig. S2, Supplementary Information) obtained for dogbone specimens with 0° and 90° raster angles, 1T39.6MPa= and 2T25.6MPa= , respectively, while the longitudinal shear strength was taken equal to half of the tensile strength measured for the specimen with +45/ -45° raster angle, 1231.2 / 2 15.6 MPa== . The compressive strength is assumed to be the same as the tensile streng th in both the longitudinal and transverse direction, 9 1T 1 C= and 2T 2 C= . The transverse shear strength 13 is assumed to be the same as in -plane shear strength, and α was set to 0. The effectiv e stress tensor is defined as () () ()1 1 0 0 0 1 1 0 0 0 1 1f m sd d d−  =− −σσ (9) where σ is the Cauchy stress tensor [34]. The parameters df, dm and ds are damage variables associated with axial, transverse and shear damage, respectively, and are used to compute the damaged material response as follows: ()()() ()() () ()1 21 1 12 2 21 1 1 0 11 1 1 0 0 0 1f f m f m m sd E d d E d d E d ED d GD − − −  = − − − − σ . (10) In the latter equation, G is the in -plane shear modulus and 12 21 1 (1 )(1 )fm D d d  = − − − . The fiber, matrix and shear damage variables are defined as 11 11for 0 for 0fT f fCd dd  = (11) 22 22for 0 for 0mT m mCddd  =  (12) ()()()() 1 1 1 1 1s fT fC mT mCd d d d d= − − − − − (13) where fTd , fCd , mTd and mCd are damage variables associated with the four modes of failure (see Eqs. (5) -(8)). Note that the latter damage variables take values between 0 (undamaged state) and 1 (fully damaged state) and are in itially set to 0 to obtain a linear elastic response prior to the initiation of damage, as seen from Eq. (10). Once a specific damage criterion (Eqs. (5) -(8)) is met, the corresponding damage variable is activated and evolves as a function of the deformati on state [56]. Here, the evolution equations for the damage variables are chosen such to obtain a linear stress softening response for each failure mode, ensuring that a certain amount of fracture energy is dissipated when the fully damaged state is reached (see [56] for furthe r information). The 10 fracture energies associated with longitudinal tension, * 1TG , and transverse tension, * 2TG , were taken equal to 90% of the total area under the stress -strain curve measured in uniaxial tension for the specimens with 0° and 90° raster angles, respectively, and the latter values were assumed to be identical for both tension and compression; hence ** 11TCGG= and ** 22TCGG= . The parameters chosen for the Hashin failure criteria and the progressive damage model are summarized in Table 2. The DENT specimen is modelled as 2D plane -stress continuum . The thickness of the samples is small (only ~1mm) and thus plane -stress assumptions were made. Plane -strain assumption is typically made if out -of-plane dimension is significantly larger than the in -plane dimensions and hence not valid in this case. Moreov er, the Hashin Failure criterion implementation in ABAQUS is limited to plane -stress elements. This is likely due to the fact that the implementation of the Hashin Failure criterion implementation in ABAQUS was primarily aimed at composite laminates (which like our experimental samples are quite thin). To mimic the loading conditions applied in the tests, Dirichlet boundary conditions were applied on the DENT specimen by constraining all DOFs of the nodes at one end, while imposing a linear ramp on the axia l displacements on all nodes at the opposite end. The characteristic element size in the finite element mesh was chosen as 0.25 mm throughout the specimen (see Fig. S3, Supplementary Information). The latter choice was ascertained by performing a sensitivi ty study on the L = 4 mm DENT specimen with 0° raster angle, which showed that further decreases in the mesh size led to less than 3 -5% change in the calculated strength. The element type chosen is CPS4R, which is a 4 -node bilinear plane stress quadrilater al element with reduced integration and hourglass control, and the total number of elements is about 42,000. Table 2 Parameters chosen for the progressive damage model. Description Symbol Value Unit Longitudinal tensile strength 1T 39.6 MPa Longitudinal compressive strength 1C 39.6 MPa Transverse tensile strength 2T 25.6 MPa Transverse compressive strength 2C 25.6 MPa In-plane shear strength 12 15.6 MPa Transverse shear strength 13 15.6 MPa 11 Shear stress factor Α 0 - Fracture energy for longitudinal tension * 1TG 0.917 J mm-2 Fracture energy for longitudinal compression * 1CG 0.917 J mm-2 Fracture energy for transverse tension * 2TG 0.509 J mm-2 Fracture energy for transverse compression * 2CG 0.509 J mm-2 4. Results and discussion 4.1 Load -displacement responses In Figs. 2a -c we plot the predicted (dashed lines) and measured (solid lines) load –displacement curves of DENT specimens with raster angles of 0°, 90° and +45/ -45°; the contours included in each fi gure relate to specimens with different ligament length L = 4, 6, 8 or 12 mm. The curves in Figs. 2a -c show an initial linear elastic response followed by a blunt peak and a steep drop in load, coinciding with the propagation of a crack through the IFPZ th at had formed in the ligament of the DENT specimen. For the 90° infill pattern, the ultimate loads are generally lower, and the peaks appear to be sharper than those associated with the 0° and +45/ -45° raster angles, suggesting that the 90° specimens fract ure in a more brittle manner. We also observe that the measured load - displacement curves generally follow a pattern of self -similarity with the total area under the curve, Wf, the peak load, Pmax, and the extension at break, Δ b, increasing with increasing ligament length, L. Note that the latter self -similarity of the load -displacement curves is an important criterion for the validity of the EWF method and ensures that the cracks propagate under similar conditions [57-59]. Similar trends are predicted numerically by the FE model which is capable of replicating some of the measured load -displacement curves with good accuracy, particularly those shown in Fig. 2a for the 0° raster angle and for specimens with lower ligament length. However, the model consistently unde r predicts the measured peak loads for the specimens corresponding to ligament length, L =12 mm. Further details on the predicted failure modes are provided in Fig. S3 (Supplementary Information), showing contour plots of the damage variables df, dm and ds at the onset of fast fracture. The predicted and measured failure (or peak) loads, Pmax, are plotted as function of the ligament length L in Figs. 2d -f for 0°, 90° and +45/ -45° raster angles, respectively. The FE model predicts a nearly linear relation between Pmax and L, due to the observed self - similarity in the load -displacement curves (see Figs. 2a -c), and similar trends are observed fo r the measurements. 12 4.2 Visualization of ligament strains via DIC Fig. 3a shows the measured load -displacement curve for a DENT specimen with L = 12 mm and 0° raster angle along with contour plots of the longitudinal strain yy in the vicinity of the ligament, as measured in situ using DIC. During the initial phase of the response, the ligament strains resemble the Mode -I strain field of a linear -elastic isotropic solid [60], as seen from frames I and II. A clo se examination of the strain field reveals that a small plastic zone exists at the tip of the notch where the strains exceed the material’s yield strain ( 1.9%Y= for the 0° raster angle). This is observed even at low values of displaceme nt (i.e. frames I and II), due to the stress concentration effect produced by the notch. As the applied load approaches the peak load, the intensity of the strain field in the ligament increases significantly and the strain starts to localize in the form o f two strip -shaped zones within the ligament (see frame IV), which are the precursors to the propagation of a straight crack through the ligament, as shown in Fig. S4 (Supplementary Information). Similar information is shown in Figs. 3b -c for the case of 9 0° and +45/ -45° raster angles (both with L = 12 mm). At small displacements, the DENT specimen with +45/ -45° raster angle (frames I and II in Fig. 3c) shows similar strain contours as those associated with the response of the 0° specimen (see Fig. 3a). How ever, the zone of localized strains (i.e. the red zone in the contour plot) starts to expand further outside the ligament in the +45/ -45° specimen as the peak load is approached (see frame IV in Fig. 3c); this is not observed for the 0° specimen. The latte r observation can be explained by the fact that the +45/ -45° specimen fractures in a zigzag pattern along the relatively weak interfaces between adjacent beads which are subjected to a combination of tensile and shear loads, as shown in Fig. S4 (Supplement ary Information). Due to the brittle behaviour of the DENT specimen with 90° raster angle, the strain field could not be clearly resolved within the strain ranges shown in Fig. 3b but will be further discussed below. In Fig. 4, we compare the strain fields of the three different raster angles for DENT specimens with L = 4 mm (shortest ligament) and L = 12 mm (longest ligament), respectively. Note that the strain fields in Fig. 4 were evaluated via DIC at the point of u ltimate failure. Here, an effort was made to illustrate changes in the shape and size of the plastic zones by selecting a strain range limited by the yield strain measured in uniaxial tension for each raster angle. Note that the yield strain is defined her e as the total strain where the plastic strain would reach 0.2% upon unloading and was determined from the stress -strain curves in Fig. S2 (Supplementary information) as 1.9%, 13 1.6% and 2.0% for specimens with raster angles of 0°, 90° and +45/ -45°, respecti vely. Hence, the red zones in the contour plots shown in Fig. 4 can be considered equal to the size of the plastic zone in the ligament. For DENT specimens with L = 4 mm, a strip -shaped plastic zone exists for all raster angles considered here, while for L = 12 mm, the plastic zones occupy significantly more material volume, particularly at the centre of the ligaments, due to the fact that the plastic work scales (theoretically) with the square of the ligament length, 2 pWL : (see discuss ion in section 2.3). For all choices of L and raster angles shown in Fig. 4, the material in the ligament has fully yielded before undergoing fracture, giving us confidence that the EWF concept is applicable for the DENT samples analysed in this study. How ever, for the 90° specimen with the largest ligament size L = 12 mm, we notice an elastic core zone at the ligament centre. While the existence of the an elastic core might alter the fracture behaviour of the specimen, the 2 pWL : scalin g, as required by the EWF concept, remains unaffected, since this core region is surrounded by plastically deformed material. 4.3 Essential work of fracture assessment Having analyzed the load -displacement curves and the strain fields induced in the DENT specimens, we proceed to evaluate the essential and non -essential work of fracture associated with each raster angle. The total work of fracture, Wf, is evaluated by cal culating the area under the load- displacement curves for each DENT specimen (see Fig. 2) while the specific work of fracture wf is obtained by dividing Wf by the ligament area: 0f fW Pd = and f fWwLt= . (14) Here, Δ is the displacement applied on the specimen boundary, P is the induced load and Δ f is the displacement at ultimate failure. In Figs. 5a -c, we plot the predicted and measured total work of fracture, Wf, as a function of the ligament length, L, for 0°, 90° and +45/ -45° raster angles, respectively, while the corresponding wf vs. L trends are presented in Figs. 5d -f. The trend lines in Figs. 5d -f represent eq. (4) and were obtained through least square fits to the predictions and measurements, respectively. It is found that that the linear fitting lines indeed provide a good descrip tion of the measured and predicted wf vs. L trends, in line with the EWF theory [61]. For the 0° and 90° raster angles, very good 14 agreement is reported between the predictions and measurements; however, larger discrepancies are observed for the specimens with +45/ -45° raster angle. The intercepts and slopes of the fitting lines in Figs. 5d -f were used to determine the essential work of fracture (EWF) and non -EWF, and are denoted as we and βwp, respectively, in accordance with Eq. (4). The predicted and measured EWF a nd non -EWF of specimens with 0°, 90° and +45/ -45° raster angle are listed in Table 3 and plotted in Fig. 6. The experimental data show that the 0° and +45/-45° specimens possess a comparable EWF of 7.24 and 7.40 kJ m-2, respectively, while the ones with 90 ° raster angle possess a significantly lower EWF of only 3.61 kJ m-2 (see Table 3 or Fig. 6a). Similar trends are observed experimentally for the non -EWF, reporting 0.71, 0.46 and 0.80 kJ m-3 for 0°, 90° and +45/ -45° specimens, respectively (see Table 3 or Fig. 6b). The FE model predicts a higher EWF for the 0° and 90° specimens, but under -predicts the EWF of the +45/-45° specimen (see Fig. 6a), in line with what observed from Figs. 5d -f. Furthermore, the FE predictions of the non -EWF lie below the measured data by a considerable margin for all choices of raster angles, as shown in Fig. 6b. Although the FE model is capable of predicting the total specific work of fracture, wf, of the 0° and 90° specimens with good accuracy (see Figs. 2d -f), it is unable to c orrectly decompose wf into essential ( we) and non -essential ( βwp) components (see Figs. 6a -b), reporting an over -prediction in EWF and an under -prediction in non -EWF. The latter observation can be explained by the fact that the formulation of the model does not permit a clear distinction between the plastic dissipation in the OFPZ and the energy required to propagate the crack through the IFPZ, which would require a different constitutive description. For the specimens with +45/ -45° raster angle, the FE model under -predicts the total specific work of fracture, wf (see F ig. 5f), which, in turn, causes the predicted values of EWF and non -EWF to lie below the corresponding measurements (see Fig. 6 ). In the load -displacement curves the slope of the predicted curve matches well with that of experimental curves for low displa cements (See Figure 2a –c), suggesting that the choice of orthotropic material behaviour is reasonable. The significant mismatch of predicted fracture toughness values with those of experiments is likely due to one of the following two reasons: (i) the dam age modelling capabilities in ABAQUS FEA, namely the Hashin Failure criterion is based on simplistic assumptions and (ii) the simplifications and assumptions involved in the choice of parameters used for the analysis (See Table 2, for example, the value of the compressive strength is assumed to be the same as tensile strength). Parametric study with different values of α and exploring more advanced failure theories like phase 15 field criteria is left to subsequent work. The deviations are primarily due to no n-repeatable microstructure (rather, defect -structure) of the specimens stemming from the non -uniform quality of the filament feedstocks, under the same processing conditions. We note that the EWF values listed in Table 3 are comparable to those reported b y Luna et al. [62] for compression molded ABS polymer which were found to range between 3.4 -4.7 kJ m -2 at room temperature. Table 3 Essential and non -essential work of fracture evaluated for the three different raster angles considered in this study. Raster angle EWF, ew (kJ m-2) Non-EWF, pw (kJ m-3) Experiment FEA Experiment FEA 0° 6.9 ± 1.38 9.58 0.74 ± 0.18 0.44 90° 3.46 ± 0.91 5.04 0.46 ± 0.11 0.35 +45/-45° 7.01 ± 3.12 5.38 0.89 ± 0.42 0.51 Fig. 7 presents SEM scans of the fracture surfaces of DENT specimens with 0°, +45/ -45° and 90° raster angles. All fracture surfaces are microscopically rough and show a combination of dimples and fibrillated structures. It is likely that the formation of t he dimples resulted from extensive cavitation inside the polymer particles which then coalesced to form larger cavities, as those visible in the micrographs, while the fibrillated structures appear as a result of continued stretching and tearing of the ligaments between the cavities [63]. The fibrillated structures are most pronounced in the specimen with 0° raster angle (see Fig. 7a) where individual fibrils reach a length of >10 μm. For a raster angle of 90°, however, the dimples on the fracture surface appear to be smaller and the formation of fibrils cannot be clearly observed which suggests that plastic flow was highly restricted and led to fractur e of the ligaments s oon after cavitation commenced, in line with the results obtained from our macroscopic mechanical tests (see Fig. S2, Supplementary Information). The information presented in this section allows concluding that the choice of raster angl e does not only affect the macroscopic fracture properties of the 3D printed ABS but can also alter the microscopic deformation mechanisms operative during the fracture process. 16 5. Conclusions In this study, we examined experimentally and numerically the e ffect of the raster angle on the essential and non -essential work of fracture of acrylonitrile butadiene styrene (ABS) thermoplastics processed via FFF AM. Tensile tests were conducted on DENT specimens with three different raster angles (0°, +45/ -45° or 9 0°) and their essential and non -essential work of fracture (EWF and non -EWF) were determined via least square fits to the experimental data. To examine the effect of the raster angle on the details of the ligament strain field, digital image correlation (D IC) analysis was performed in -situ, and contour plots of the longitudinal strains in the vicinity of the ligament were constructed. In addition, we developed a 2D FE model to predict the response of the 3D printed DENT specimens, taking into account the an isotropic elastic response and fracture properties associated with the FFF print. The measurements show that the resistance to fracture initiation of 3D printed DENT specimens is substantially higher when the printing direction is perpendicular to the crac k plane (0° raster angle) as compared to the samples where the printing direction is parallel to the crack (90° raster angle), reporting EWF values of 7.24 kJ m-2 and 3.61 kJ m-2 for the 0° and 90° infill patterns, respectively. A high EWF value was also r eported for the +45/ -45° raster angle (7.40 kJ m-2) where the beads of the FFF print are subjected to pronounced shear loading. The DIC analysis showed that connected plastic zones existed in the ligament of the DENT specimens prior to the onset of fractur e, and this was corroborated by SEM fractography which showed that fracture proceeded by a ductile mechanisms involving void growth and coalescence followed by drawing and ductile tearing of fibrils. The ultimate strength and fracture properties predicted by the FE model were found in good agreement with the measurements obtained for the different DENT specimens, while larger discrepancies were reported when the ligament length was L =12 mm. Future work should involve in -depth study of the effects of the different parameters of the Hashin criteria as well as to explore the efficacy of different failure criteria like the Tsai Wu or Puck criteria among others. The results presented in this study offer physical insight into the fracture mechanisms of 3D print ed ABS and provide useful guidelines for the design of fracture resistant thermoplastics processed via Fused Filament Fabrication. Acknowledgements SK would like to thank the University of Glasgow for the start -up grant [award no: 144690 -01]. AJ would like to thank Indian Institute of Technology Kharagpur for the ISIRD grant. 17 Authors’ contributions Experimentation: Pawan Verma; numerical modelling: Ja bir Ubaid, Atul Jain; writing (original draft preparation): Pawan Verma, Andreas Schiffer, S Kumar, Jabir Ubaid, Atul Jain, Emilio Martínez -Pañed a; writing (review and editing): Pawan Verma, Andreas Schiffer, S Kumar, Jabir Ubaid, Atul Jain, Emilio Martínez -Pañed a. Funding S Kumar would like to thank the University of Glasgow for the start -up grant [award no: 144690 -01] Availability of data and mat erial The authors confirm that material supporting the findings of this work is available within the main text and supplementary Information. Compliance with ethical standards Competing interests The authors declare that they have no competing interests. Ethical approval The article follows the guidelines of the Committee on Publication Ethics (COPE) and involves no studies on human or animal subjects. Consent to participate Not applicable. The a rticle involves no studies on humans. Consent to publish Not applicable. The article involves no studies on humans.