1. Introduction
Additive manufacturing (AM), also known as 3D printing (3DP), rapid prototyping (RP), or solid freedom (SF), is a process of creating 3D model components by gradually joining materials layer by layer. Although it was first introduced three decades ago, it is still gaining a lot of attention from researchers and investors worldwide. Generally, a design can be generated using computer-aided design (CAD) software with complex 3D shapes and exported as a surface tessellation language (STL) file to be compatible with 3D Printer software. Then the slicer software slices the 3D complex model as multiple 2D layers and the resulting toolpath file is sent to the machine for printing. The software also generates a simulation view of the printing process, allowing the designers to verify each printing layer and calculate the estimated time and extruded filaments. The 3DP of polymeric parts prefers thermoplastic polymer materials, such as acrylonitrile butadiene styrene (ABS), polylactic acid (PLA), and polycarbonate (PC), as well as thermosetting polymer materials, including epoxy resins [
1].
Since 3DP was first introduced, various 3DP techniques have been developed to accomplish different criteria. Notable current technologies are fused deposition modelling (FDM), powder bed and inkjet head 3DP, stereolithography (SLA), and 3D plotting/direct-write. FDM printing is amongst the most common printing techniques and thermoplastic polymers, such as PC, ABS, and PLA, are often selected as the preferred materials for printing due to their low melting temperatures. Before printing, the quality of the part is controlled in the slicer software by carefully selecting the best setting parameters, such as layer thickness, skirt, support, raster angle, fan speed, printing speed, etc. [
1]. In the FDM process, which is the focus of the current study, the filament is fed into a heated liquefier with the assistance of a gearbox. The heat block melts the filaments to a suitable temperature for different thermoplastic materials and feeds them into the nozzle. The nozzle head moves to controlled positions to place molten plastic filaments to form the first layer. After one layer is formed, the nozzle head moves upward and repeats the deposition process to form the next layer [
2,
3]. Some issues related to FDM printing are void formation during the manufacturing of the parts, layer delamination of the printed polymers, limited choice of thermoplastic polymer materials due to the nature of the method, and material anisotropy [
1,
4]. Amongst the common issues in FDM, anisotropy of material properties of 3D-printed parts, which can be divided into three categories (mechanical, electrical, and thermal anisotropies), was reported as the most significant issue [
5,
6,
7]. In the studies by Türk et al. and Mohan et al. [
8,
9], the low-layer adhesion and air voids produced by the printing method were confirmed to be the main reasons for the low mechanical, electrical, and thermal properties of 3D-printed parts printed in transverse build orientation.
According to the review conducted by Zohdi, et al. [
6], adding additives [
10,
11], build orientation [
12,
13], and bed and nozzle temperatures [
14,
15] are the top three control factors and parameters that can significantly contribute to mechanical anisotropy. Out of those, the effect of build orientation on the mechanical properties of the printed parts has gained a lot of attention from researchers [
6,
12,
13]. By evaluating the tensile properties of the printed parts, it has been found that the build orientation can affect the mechanical properties of the printed parts and can greatly influence the manufacturing time and the related costs. Layer delamination and air voids from the microstructural analysis were identified as the main reasons for the low strength results and, thus, can be identified as contributing factors in the mechanical anisotropy of printed parts via FDM. In one study conducted by Dul, et al., comparing 3D-printed parts with the compression moulded parts for ABS/Graphene affirmed the effect of build orientation on the mechanical properties of the fabricated samples [
16]. In addition to the build orientation influence, material type has an enormous influence on the mechanical properties of the printed parts. In our preliminary study by Zohdi et al., the mechanical anisotropy in parts 3D printed with high-impact polystyrene (HIPS) and ABS was investigated [
7]. This study revealed that the samples printed in two different build orientations of longitudinal and transverse with HIPS showed a negligible anisotropy degree compared to ABS samples. SEM images revealed a lower-layer adhesion for the parts made of ABS polymers compared to the HIPS polymers. Moreover, relative tensile strength values of 3D-printed HIPS samples compared to the values for the mould-injected samples prove that by choosing a proper polymer and printing at certain print parameters, almost the same mechanical properties as those of the mould-injected samples could be achieved.
Apart from the mechanical anisotropy, the effects of different printing parameters (infill density, temperature, and layer thickness) on printed parts have also been investigated in literature [
17,
18,
19]. For example, infill density was reported to influence the tensile strength and dynamic mechanical properties of conductive ABS/zinc oxide (CABS/ZnO) composites printed by FDM. Infill density profoundly improved the tensile strength and Young’s modulus properties and decreased the elongation at break. The noticeable increase in dynamic modulus happened when the infill density changed from 50% to 100% [
20]. A change in raster angle from 45°/−45° to 90°/0° showed an 8 per cent increase in mechanical strength for ABS [
19]. In another study by Rajpurohit, et al. [
21], the highest tensile strength was found to be at the 0° raster angle for PLA printed by FDM. The bonding strength was found to vary with the change in build orientation and level of layer thickness. With lower layer thickness, edgewise configuration indicated the highest bonding strength, while with higher layer thickness, flatwise configuration proved to have the highest bonding strength [
22]. As stated before, among different types of thermoplastic polymers for 3DP, this article focuses on the parts made of ABS polymer. ABS is a widely used polymer, mostly used in packaging, toys, bottles, housewares, electronic appliances, and light-duty industrial components because of its good rigidity and ease of colouring and processing. However, its applications are limited because of its relatively poor impact strength, heat deflection, and flame retardancy [
23].
To compare the experimental results with modelled values, different multiscale modelling methods on CAD programs were referenced by different researchers. For example, sequential multiscale modelling was employed by He, et al. [
24] to transfer the effective properties of 3D-braided composites from microscale to mesoscale and from mesoscale to macroscale. One typical method to execute multiscale modelling is constructing a 3D representative volume element (RVE). MSC Digimat-FE software was employed for this method to build 3D microstructure RVE models of aluminium nanocomposites and perform microstructural deformation analysis. The simulation results indicated great agreement between the tensile property of the created RVE models and the experimental values [
25]. In another study, MSC Digimat-FE software was also used to construct a 3D microstructural model of BaTiO
3 to compute the effective elastic constants and epoxy adhesive composites [
26].
For extracting the material properties of composites in modelling, the rule of mixtures (ROM), also known as the Voigt model, can be acquired to compute the overall properties based on the properties of each region and its volume fraction [
27]. With the analysis of scanning electron microscopy (SEM) images, the volume fraction can be reduced to the area fraction of the cross-section by assuming the constant thickness of testing samples [
28,
29]. In a study by Deng, et al., the successful application of the ROM to determine the elastic modulus and tensile strength of carbon-fibre-reinforced nylon (CFRN) fabricated by FDM was recorded [
28]. Several versions of ROM were employed to validate the experimental results of sisal fibre-reinforced polystyrene composites [
27]. For predicting longitudinal properties, Voigt and Halpin–Tsai models showed agreement in the equations. However, the Halpin–Tsai model suggested a correction factor for the transverse properties, while the original Voigt model was modified and replaced by the Reuss model to enhance the accuracy. Applying the developed models to derive the properties of the samples printed in longitudinal build orientation exhibited good agreement with experimental results, while higher discrepancies were observed for transverse properties, especially for lower fibre loadings [
27].
In this study, a systematic multiscale modelling and simulation process is developed to model the 3D-printed parts made of ABS polymer at the macro- and mesolength scales, respectively. Firstly, the macroscale analysis is devised using FE modelling as the reverse engineering method and it is employed to investigate the effects of infill density at different build orientations. Secondly, the mesoscale analysis is developed by devising two methods to extract the material properties of the ABS polymers. The first RVE-based numerical method investigates the effect of the process parameters, while the second Rule-of-Mixtures-based analytical method evaluates differences in material properties at two build orientations and aims to highlight the effect of mechanical anisotropy. Lastly, the numerical and analytical results were compared with the experimental results and a modified rule of mixtures was proposed. The whole paper is outlined in the following sections:
Section 2 focuses on the development of the experimental and numerical analysis procedure;
Section 3 presents the obtained results and conducts the discussion; and
Section 4 draws conclusions and outlooks for future work.
2. Experimental and Numerical Analysis
2.1. Sample Design and Experimental Procedure
Polymer pellets of ABS were used to extrude filaments with an outer diameter of 1.75 ± 0.07 mm. Pellets of ABS polymer grade PA747 (C
8H
8·C
4H
6·C
3H
3N)
n were purchased from the ChiMei Corporation, Taiwan. The dogbone samples were designed according to the ASTM D638 standard [
30]. The standard geometry and dimensions of selected Type V were employed to create the CAD model of the dogbone samples using SolidWorks.
Samples were fabricated in batches of five with a Prusa i3 MK3S+ printer and its embedded slicer (
Figure 1). Five replications for each parameter were prepared to minimise the error percentage. Four sets of samples were produced to investigate the effects of different printing parameters, including infill density, layer thickness, raster angle, and temperature. For each printing parameter, the samples were fabricated at two different build orientations: longitudinal and transverse. For the infill density set, other machine factors were set as default, while infill density was modified to 40%, 60%, 80%, and 100%, respectively. Next, the layer thickness set was created with the same method by changing the layer thickness to 0.15 mm, 0.25 mm, and 0.35 mm, respectively. Next, samples in the raster angle set were produced at ±45° to compare with the original 0°/90° raster angle. Finally, the printing temperature was varied from 220 °C to 280 °C with an increment of 20 °C to complete the temperature set.
An Instron 3365 machine with a 5-kN load cell was employed to extract the mechanical properties of produced samples. From the raw data of the five samples from each batch, the average ultimate tensile strength and Young’s modulus were calculated as a benchmark for the simulation process.
SEM images were collected using a Hitachi Flex SEM 1000 II machine. The fractured surfaces of the dogbone samples from mechanical testings were flattened and prepared with the assistance of a microtome instrument prior to imaging. Samples were first submerged into liquid nitrogen for at least 20 min and then immediately prepared using a microtome instrument. Images were collected at various magnification levels to fit the purpose of modelling.
2.2. Macroscale Analysis
The macroscale model and analysis were developed for the tensile specimens to perform static structural analysis using ANSYS Workbench and MSC Digimat Software. Its three-step workflow is displayed in
Figure 2.
The CAD model of the tensile specimens was created using SolidWorks and then imported into ANSYS Workbench to devise a standard static structural analysis with appropriate boundary conditions, according to the ASTM D638 standard. Meanwhile, the toolpath containing moving tools’ locations in 3D space was generated using the Prusa slicer and the material data were prepared with the mean-field theory using MSC Digimat-MF. Next, all those data were further used as inputs in the FE model developed using MSC Digimat-RP to simulate the 3DP process of the tensile specimens. ANSYS Solver was used to perform finite element analysis (FEA) via the interface between MSC Digimat and ANSYS. Finally, the results were viewed and displayed on the structural model in ANSYS Workbench for final data processing to extract the macroscopic failure indicator (MFI) values and then tensile strength.
Since the parts printed with ABS polymer had transverse isotropic behaviours, the material’s failure was modelled using the Tsai–Hill criteria available in MSC Digimat. The Tsai–Hill criteria require the input parameters as axial tensile strength (X), in-plane tensile strength (Y), and transverse shear strength (S) to assign an MFI value to each element. Those parameters can be extracted from the properties of samples printed at 0/0 raster angle–longitudinal build orientation and 90/90 raster angle–transverse build orientation. However, most currently available printers, including the employed printers for this project, could not accomplish the necessary settings for the above raster angles. Hence, the macroscale analysis was used to apply the reverse engineering method to extract the required parameters and match the tensile strength obtained for the sample manufactured at 100% infill density. The tensile strength was the maximum equivalent von-Mises stress when the MFI values in a location of the specimen exceeded 1.0. For calibration purposes that required multiple trials, this process was conducted at an element size of 1 mm.
In FEA, the result accuracy is strongly influenced by the element size. With a smaller element size, more elements are generated, thus, improving the accuracy of the results. Nonetheless, it also increases the computation time and might exceed the hardware limitation. Therefore, a mesh convergence test was conducted to derive the appropriate element size, guaranteeing accuracy and saving computation time for the macroscale analysis. For the convergence test, the maximum equivalent von-Mises stress was selected as the scrutinised criterion.
In 3DP, the printing time can be reduced at lower infill densities, but the 3D-printed materials experience more air voids. Air void is a crucial defect affecting the material’s mechanical properties and the software must recognise the infill density effect to successfully simulate the properties of 3D-printed parts [
6,
7,
31]. Therefore, the final aim of the macroscale analysis was to investigate the software’s capabilities in detecting the infill density effect. The toolpath files that contained the infill density information were replaced accordingly for each infill density ranging from 40% to 100% with an increment of 10%.
2.3. Mesoscale Analysis
The process of mesoscale modelling adopted the SEM images to create the representative volume elements (RVEs) and perform the FE analysis.
Figure 3 shows the SEM image of ABS-fabricated samples printed with 100%, 80%, 60%, and 40% infill density and longitudinal orientation, respectively. In these images, Regions 1 and 3, framed in blue, were identified as perimeter regions, while Region 2 represents the infill, which is framed in yellow. With the information from the SEM images, two methods were developed to study the material properties (tensile strength) of 3D-printed samples in MSC Digimat-FE. The results can be generalised to compare with experimental data by assuming the ultimate tensile strength occurs in the middle (gauge length) region of the dogbone specimens.
From the SEM image, one typical intrinsic element in each region was selected to measure the maximum width (W), maximum height (H), edge’s width (W
1), and edge’s height (H
1) to define the shape of the intrinsic element depicted in
Figure 4a. In the first method, the CAD model of the entire structure as shown in
Figure 4d, including perimeter and infill regions, was designed and analysed using MSC Digimat FE. The second method utilised the printing parameters, including nozzle diameter (0.4 mm) and layer thickness (0.2 mm), and manipulated the available models, the aligned and sparse models, to replicate the perimeter and infill regions as shown in
Figure 4b and
Figure 4c, respectively. After that, the tensile strength of the entire structure was computed using ROM, as shown in Equation (1). In this equation,
are the tensile strength of the perimeter region, infill region, and the resultant model combining all three regions, respectively. By selecting the constant thickness for all regions, the volume fraction
of the perimeter region to the resultant model was reduced to only define area fraction from the SEM images [
28].
With the same limitations described for the macroscale analysis in mind, the original material data for all methods were pre-calibrated to match the mechanical properties of the samples printed at 100% of infill density in the first method. Then the same material file was used to perform simulations for other sets to maintain consistency. Then, in each set (for instance: infill density, layer thickness, etc.), the properties of longitudinal and transverse orientations were extracted by modifying the loading directions onto the RVE models.
The first method was capable of generating the results for all four sets of parameters, including infill density, layer thickness, and temperature of the ABS polymer, when the selected build orientation was set as longitudinal. The FE model was automatically meshed with the non-conforming tetra method with local mesh refinements. The generated stress–strain curve at the macroscopic level was recorded to extract the mechanical properties. For the second method, a mesh convergence test was conducted, and an appropriate element size of 0.025 mm was applied to generate data for the infill density at both longitudinal and transverse orientations.
To the best of the author’s knowledge, no studies previously reported the use of ROM for the same purpose in the second method. Hence, the original ROM can be modified to enhance its accuracy. Based on the original ROM and the obtained data, an additional term was included to adjust the ROM suitably for each percentage of the infill density. Equations (2) and (3) were used previously in other studies, mostly focusing on composite and nanocomposite materials [
32]. However, so far, these equations were not used to describe the issues associated with the ROM for the polymer parts printed with AM. In these equations, additional functions of f(x) and g(x) were inserted into the original ROM to predict the material properties of the longitudinal and the transverse samples, respectively. In both additional functions, the variable x is the percentage of infill density. For instance, if the sample is fabricated at 80% infill density, the value of x is 0.8. For this study, the evaluated function was set in a quartic form, a polynomial of degree four with the same form as P(x) in Equation (4), and the coefficients of the functions
were derived in MATLAB using the collected data. The outcome of this task was expected to assist future modellings of the samples fabricated at different infill densities with the same pure polymer and enhance the accuracy of the original model.
4. Conclusions
In this study, the systematic FE-based multiscale numerical modelling and simulation processes, including the FE-based model at the macroscale and two methods at the mesoscale, were successfully developed for the ABS polymer samples printed using the FDM method.
For the proposed macroscale analysis, the following conclusions can be drawn:
With the 2.7% difference between the numerical and experimental tensile strength of transverse samples, the macroscale analysis revealed that the reverse engineering method can be used to determine material data, considering the effects of build orientation and infill density, which provides a close fit with those experimental data for the longitudinally built samples.
The longitudinal build orientation exhibits linear properties related to the change in the infill density, while the transverse build orientation does not show any specific trends due to its unpredicted behaviours.
At the mesoscale, the mesoscale analysis of the polymeric samples extracted the following research findings:
With the process of reproducing the entire RVE structure, including the perimeter and infill regions, the tensile strength results using the first method showed close relations with experimental data of longitudinally built samples.
Further, with less than 10% error percentages, the first method was proved to be the best to replicate those unexpected behaviours, considering printing parameter sets, e.g., layer thickness, raster angle, and temperature.
By employing ROM, the calculated tensile strengths using the second method have close relations with those experimental data at higher infill densities and significant deviations at lower infill densities.
For the transversely built samples, the modified model of ROM was proposed using parameters to improve the prediction accuracy.
Future work can be attempted to reproduce more accurate intrinsic models by obtaining all the SEM images of the cross-sections’ front, right, and top views and, thus, will wholly demonstrate the natural defections caused by the 3DP process. Furthermore, the devised macroscale and mesoscale modelling and analysis procedures can also be employed to model the nanocomposites’ material properties.