Next Article in Journal
Investigation of Welding Parameters of Dissimilar Weld of SS316 and ASTM A36 Joint Using a Grey-Based Taguchi Optimization Approach
Previous Article in Journal
Load Introduction Specimen Design for the Mechanical Characterisation of Lattice Structures under Tensile Loading
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanical Performance over Energy Expenditure in MEX 3D Printing of Polycarbonate: A Multiparametric Optimization with the Aid of Robust Experimental Design

by
Nectarios Vidakis
1,*,
Markos Petousis
1,
Constantine N. David
2,
Dimitrios Sagris
2,
Nikolaos Mountakis
1 and
Emmanuel Karapidakis
3
1
Department of Mechanical Engineering, Hellenic Mediterranean University, 71410 Heraklion, Greece
2
Department of Mechanical Engineering, International Hellenic University, 62124 Serres, Greece
3
Electrical and Computer Engineering Department, Hellenic Mediterranean University, 71410 Heraklion, Greece
*
Author to whom correspondence should be addressed.
J. Manuf. Mater. Process. 2023, 7(1), 38; https://doi.org/10.3390/jmmp7010038
Submission received: 26 December 2022 / Revised: 23 January 2023 / Accepted: 28 January 2023 / Published: 1 February 2023

Abstract

:
Sustainability and energy efficiency of additive manufacturing (AM) is an up-to-date industrial request. Likewise, the claim for 3D-printed parts with capable mechanical strength remains robust, especially for polymers that are considered high-performance ones, such as polycarbonates in material extrusion (MEX). This paper explains the impact of seven generic control parameters (raster deposition angle; orientation angle; layer thickness; infill density; nozzle temperature; bed temperature; and printing speed) on the energy consumption and compressive performance of PC in MEX AM. To meet this goal, a three-level L27 Taguchi experimental design was exploited. Each experimental run included five replicas (compressive specimens after the ASTM D695-02a standard), summating 135 experiments. The printing time and the power consumption were stopwatch-derived, whereas the compressive metrics were obtained by compressive tests. Layer thickness and infill density were ranked the first and second most significant factors in energy consumption. Additionally, the infill density and the orientation angle were proved as the most influential factors on the compressive strength. Lastly, quadratic regression model (QRM) equations for each response metric versus the seven control parameters were determined and evaluated. Hereby, the optimum compromise between energy efficiency and compressive strength is attainable, a tool holding excessive scientific and engineering worth.

1. Introduction

Sustainability is lately a key issue for all aspects of society (government, enterprises, and public) for both products and processes [1]. AM is considered a process with better sustainability characteristics compared to traditional manufacturing processes due to its energy consumption and reduced material waste [1,2]. Additionally, the sustainability aspects of AM are gaining attention in the literature since it has become one of the important tools facilitating the manufacturing of products [3,4]. The effect of environmental factors [5,6,7], the recycling of the materials [8,9], and the process parameters [10] on the sustainability of the AM process have been investigated. One of the factors affecting the sustainability of a process is energy consumption [11]. The energy efficiency of MEX polymer processing has been investigated in the literature [12]; however, research is not yet extensive. The energy demands in MEX 3D printing of high-performance polymers, such as PEEK, have been investigated with a full factorial design, considering two 3D printing parameters, i.e., layer thickness and printing speed [13]. Interpretive structural modeling (ISM) has also been employed for the analysis of parameters affecting energy consumption in MEX 3D printing [11]. The effect of parameters, such as the infill pattern, on energy consumption, has also been investigated [14] for materials such as the ABS polymer [15]. Energy consumption has also been studied using machine learning techniques [16]. The effect of six 3D printing parameters on the energy consumption of parts built by ABS polymer with the MEX 3D printing process has been analyzed and optimized with statistical modeling tools [17]. Energy modeling techniques have also been applied [18]. The effect of energy consumption on parameters such as the required 3D printing time [19] and geometric accuracy [20] has been studied.
PC is a high-performance polymer [21], with tensile strength higher than 62 MPa, even when recycled, and a glass transition of about 121 °C [22], which makes it suitable for applications requiring high strength and high thermal loadings [23]. Therefore, it is used in safety parts, such as bulletproof glass in automobiles, helmets, and bumpers, but also electronics in capacitors [21,24,25]. In composite form, it has been used in electromagnetic shielding [26,27] and in construction [28,29] to improve the resistance to fire of materials [30]. PC blends are popular in biomedical devices and aerospace parts [21,29,31] and exhibit high energy-absorbing properties [32]. In 3D printing, it has been applied in aeronautical [33] and biomedical applications [34,35,36]. Still, the use and the research on the PC polymer in AM, overall, are not so extensive. Its mechanical properties under tensile [37] and shear forces have been investigated [38,39]. The effect of 3D printing parameters on the mechanical response of parts made with the PC polymer under tensile loading has been reported. Statistical tools have been employed for the analysis and optimization of the results [40]. Such an approach, by employing modeling tools, statistical and regression analysis, or even neural networks for the analysis and optimization of the performance of the 3D-printed parts, is often used in the literature [41,42,43,44].
The effect of 3D printing parameters on the response of PC 3D-printed parts under different strain rates in tensile tests has also been reported [45], and the impact strength of such parts as well [46]. In MEX 3D printing, the mechanical performance of recycled PC polymer has been reported, promoting a sustainable character to the polymer [23]. Quality characteristics related to the dimensional accuracy of 3D-printed parts built with the PC polymer have been studied [47]. The mechanical properties of PC composites in 3D printing have been investigated for carbon-based composites [48,49], and nanocomposites with titanium carbide [50], aluminum nitride [51], titanium nitride [52], silicon carbide [53], silica [54], and cellulose nanofibers [55]. PC/ABS blends in 3D printing have also been presented aiming to enhance the mechanical properties and the printability of the materials [56,57,58,59].
For the analysis of the performance of the PC polymer in MEX 3D printing, statistical tools, such as the Taguchi experiment design, have been employed [40,47,60,61,62], along with simulation tools [63]. The PC polymer has also been applied in hybrid additive manufacturing processes to further expand its applicability. MEX 3D printing has been combined with friction stir welding for the production of large-size parts [64] and with CO2 laser cutting processes for the improvement in surface quality characteristics [65]. The compressive strength of the PC polymer has been studied in the pure form [66], recycled form [23] in 3D printing [67,68], and specifically for aeronautical [33] and medical [34] applications. Additionally, the effect of strain rate on the compressive strength of 3D-printed parts has been reported [69]. Still, there is a gap in the literature in the thorough analysis of the effects of the MEX 3D printing parameters on the compression response of the PC polymer. Additionally, the energy demands when 3D printing PC parts with the MEX process have not yet been addressed in the literature, and especially a study is missing on how the energy requirements are affected by the 3D printing parameters.
In the current study, the effect of seven 3D printing parameters, i.e., raster deposition angle, orientation angle, layer thickness, infill density, nozzle temperature, bed temperature, and printing speed, on the energy consumption and the response under compressive loading of parts made of PC with MEX 3D printing is reported. The 3D printing settings were selected to be generic and machine independent. Ten metrics related to energy consumption, the compression properties, and the quality of the 3D printing parts are evaluated with statistical and regression modeling tools.
Despite the compressive specimens’ increased printing time, volume, and weight (when compared to tensile specimens), this work focused on this mechanical property for reasons explained further in the text. Compression test data in multi-parametric experiments are limited in the literature, probably due to the increased effort required for the preparation of the samples with the 3D printing process. Considering that the compressive loads are the most frequent in the functional 3D-printed parts life cycle [70], this shortfall of data is remarkable. The increased mass and volume of the compression test specimens facilitate the monitoring and documentation of parameters in this work, such as the printing time, weight, and energy consumption, leading to more accurate results. Additionally, data for the mechanical performance of the PC polymer in MEX 3D printing are overall limited in the literature.
Equations as functions of the 3D printing settings (control parameters) were compiled, and their reliability was verified with confirmation experiments. The energy-related metrics were monitored during the MEX 3D printing process. The PC polymer was extruded from raw material to filament form, and specimens were 3D printed with this filament to be tested for their performance under compressive load. For the preparation of the specimens and performing the compression tests, the ASTM D695-02a international standard was followed. We also examined 3D-printed specimens for their morphological characteristics before and after the compression experiments. To the authors’ best knowledge, no research so far simultaneously studies seven 3D printing parameters for their effect on energy consumption and the thorough analysis of the compressive behavior of PC MEX 3D-printed parts. The prediction models provide valuable information about the expected performance of the parts, indicating direct industrial merit. The statistical analysis showed that each parameter affects the metrics studied in this work differently, making the modeling procedure a requisite. Infill density was the dominant parameter regarding compression strength, while layer thickness was the dominant parameter regarding energy consumption.

2. Materials and Methods

2.1. Methodology for Sample Preparation and Testing

Figure 1 presents snapshots of the steps followed in this work for the preparation of the samples and the experiments conducted. The PC polymer was sourced from Styron Europe GmbH (Trinseo Europe GmbH, Horgen, Switzerland) in pellet form (grade EMERGE 8430–15, a tensile strength of 70.0 MPa, and a density of 1.20 g/cm3). Pellets were initially dried (Figure 1a) for 14 h at 60 °C. Then a 3devo Precision (3devo, Utrecht, The Netherlands) MEX filament extruder was used for the production of 1.75 mm filament compatible with the MEX 3D printing process (Figure 1b, the first three heating zones operated at 200 °C, while the fourth heating zone operated at 240 °C, and screw speed was set at 4.8 rpm). The thermal properties of the specific PC grade were evaluated with TGA (Figure 2a) and DSC (Figure 2b) to ensure that the extrusion temperatures in this work do not cause any degradation in the PC polymer.
The produced filament was dried (Figure 1c, 60 °C for 4 h). All the filaments were produced, and then the 3D printing process followed, which was not conducted on the same day. After the production of the filament, it was not stored in a vacuum case or a laboratory oven, so it might have had moisture, as polymeric materials attract moisture easily. This moisture significantly affects the 3D printing process. So, in order to ensure that there was no humidity trapped in the filament before it was used in the 3D printing process, the filament was dried initially before being used in MEX 3D printing process.
With this filament, specimens were MEX 3D printed (Figure 1d) on an Intamsys, Funmat HT (Intamsys, Shanghai, China) MEX 3D printer. Five specimens were fabricated for each different set of 3D printing parameters. The 3D printing settings are shown in Figure 3, along with the specimens’ dimensions, in accordance with the ASTM D695-02a standard. During the 3D printing process, the printing time and the consumed energy were monitored with the Rigol DM3058E device (RIGOL Technologies, Shanghai, China) (Figure 1e). The produced specimens were inspected with an optical microscope Kern OKO 1 with a 5MP ODC 832 camera (KERN, Balingen, Germany) (Figure 1f). It should be noted that in the 3D printer used (Intamsys, Funmat HT), parts are built inside a chamber. In this chamber, the filament is dried before its use. The table in which the part is built is heated (80–120 °C in the current work), and the nozzle in the 3D printing head also has a high temperature (260–300 °C in the current work). So, any humidity present in the 3D printing chamber is expected to be negligible during the 3D printing process. It should be noted also that no delamination was observed during the 3D printing process, attributed to the robustness and the geometry of the compression test specimen, which does not have thin walls or low thickness.
Finally, the specimens were experimentally tested in compression tests to determine the corresponding mechanical properties (Figure 1g). The compression experiments were conducted on a testing machine type Instron KN1200 (Instron Corp., Norwood, MA, USA). The load cell is build-in into the machine, and it has a capacity of 1200 N and an accuracy of ±24 N. Following the compression test standard, the elongation speed in the tests was set at 1.3 mm/min. It should be mentioned that the experiments were conducted according to the ASTM D695-02a international standard for testing polymeric materials under compression loading. No specific load was applied to the specimens. The load was increased from zero (0) up to the value at which each specimen failed. When the specimen failed, the experiment for the specific specimen was terminated. So, the maximum load for each specimen differs.

2.2. Energy Indicators

Regarding energy consumption, it can be distinguished into three main phases: i machine startup, ii 3D printing process, and iii machine shutdown and can be calculated utilizing the equations below [17]:
E total = E thermal + E motion + E auxiliary
where,
E thermal = E heating + E cooling
E motion is the consumed energy by the motors of the 3D printer and
E auxiliary = E startup + E steadystate + E shutdown
of the electronics and the remaining parts of the 3D printer.
The specific printing energy index is calculated by the following equation:
SPE = EPC w MJ / g
The specific printing power index is calculated by the following equation:
SPP = EPC PT w 10 3 kW / g
where energy printing consumption (EPC) represents the energy used by the 3D printer ( E total ), w the actual weight of each specimen, and PT the actual printing time for each experimental run.

2.3. Design of Experiment (DOE), Regression Analysis (ANOVA)

In this work, seven generic, continuous (not categorial) machine-independent 3D printing settings, i.e., raster deposition angle (RDA, deg), orientation angle (ORA, deg), layer thickness (LT, mm), infill density (ID, %), nozzle temperature (NT, °C), bed temperature (BT, °C) and printing speed (PS, mm/min), were the control parameters of the robust design approach followed [71]. Each control parameter had three levels, and five replicas were run for each different case studied in this work. The control parameters and their levels were selected following the literature review for the PC polymer in MEX 3D printing.
An L27 array was formed, which means 135 experiments were carried out for the modeling process. The Taguchi design of experiments was followed, as the corresponding full factorial design would require 5 × 37 experiments to be carried out. The effect of these control parameters on ten different metrics (response parameters) related to the 3D-printed parts’ energy consumption during the MEX 3D printing process (printing time—s the sample weight—g; EPC—MJ; SPE—MJ/g; and SPP—KW/g), the mechanical properties under compressive loading (compressive strength—MPa; compression modulus of elasticity—MPa; and compression toughness—MJ/m3), and the quality metrics of the parts (deviation of the samples’ cross-sectional area with the nominal cross-sectional area—Area2Nom; and volume deviation compared with the nominal volume of the specimens—Volume2Nom) was evaluated. Regression analysis provided prediction models as functions of the control parameters for each one of the response parameters. The prediction models were verified with confirmation runs, also carried out with five replicas each.

3. Results and Discussion

3.1. Examination of the Morphological Characteristics of the Samples and Their Behavior during the Compression Test

Figure 4 depicts top surface images acquired with optical microscopy from the samples. Images were taken on samples built with various control parameter values. Due to the various 3D printing parameters utilized in each sample, the changes in the 3D printing structure can be observed. An error-free 3D printing structure is shown in every image, demonstrating that the control parameter levels were suitable for producing the PC samples using MEX 3D printing. Please refer to this work’s Supplementary Material, where a microscope image from each of this work’s runs is displayed.
Compression tests are performed to ascertain the compressive strength of the 3D-printed specimens fabricated under various printing parameters, as well as to study and analyze the way they are prone to fail. Figure 5 illustrates the specimens subjected to a compression test according to the ASTM D695 standard. Based on the conducted compression tests, a classification of the resulting compressive failure modes can be established. The failure modes are grouped into three main failure types and can be further distinguished into subcategories. The first one pertains to kinking failure, either single or multiple (runs 1, 2, 3, 4, 5, 6). The kinking plane is formed at an angle of 45° to the specimen’s longitudinal axis, where the maximum shear stress occurs. The second failure type, characterized as buckling failure, consists of three discrete failure modes, i.e., lateral buckling (runs 7, 8, 9), buckling with interlaminar delamination (runs 19, 20, 21), and buckling accompanied by interfacial debonding (runs 25, 26, 27). Finally, the third failure type observed concerns shear failure comprising two distinct secondary modes, i.e., sliding shear (runs 10, 11, 12, 14) with a brittle fracture or successive degradation (runs 13, 15), and flexural shear failure (runs 16, 17, 18).
To explain the failure mechanism of the specimens, each type of failure mode is related to the 3D printing settings. It was noticed that the failure mode of the specimens is primarily dependent on the printing orientation angle (ORA). For a printing angle of 0o, the specimens exhibit flexural failure either of kinking type or lateral buckling type, demonstrating maximum compressive strength in comparison to other printing angle values, which is even enhanced proportionally by the infill density (ID). Regarding a printing orientation angle of 45°, it is obvious that it leads to shear failure since sliding at the shear plane of 45° angle is triggered. Shearing is happening in the direction perpendicular to the filament strands deposited on the shear plane.
A printing angle of 90° leads to flexural failure either with interlaminar delamination in the case of the raster deposition angle (RDA) being equal to 0° or with interfacial debonding between adjacent filaments for RDA equal to 90°. The raster deposition angle (RDA) plays a supplemental role by modifying the principal failure mechanism. For example, the RDA angle of 45° generally alters the specimen’s failure, either by the shearing deformation with respect to the longitudinal axis or by preventing the separation of layers due to high interlaminar stresses acting on the interface between two adjacent layers [72,73,74,75,76]. Finally, the test results showed that the ultimate compressive resistance of the 3D-printed specimens is consistently related to the infill density (ID). The ultimate compressive resistance is defined as the maximum compressive load the specimen is capable of withstanding during the compression test. That means that by the ultimate compression load, a rupture occurs, indicating the initialization of the material failure. Regardless of the type of failure mode taking place, the infill density (ID) of the specimen enhances the compressive strength. With respect to the contribution of the other 3D printing settings, i.e., PS, BT, NT, on the failure mode of the specimens, the results prove that these parameters do not notably affect the compressive failure mode.
Figure 6 visualizes the methodology of investigating the 3D-printed PLA specimens’ failure caused by applying compressive load. The particular case presented here concerns the shear failure mode. The test specimen is placed between the plates of the compression tool, taking care to align its longitudinal axis with the center line of the plunger and to ensure that the ends of the specimen are parallel to the plates of the compression tool. Subsequently, the compression load is applied by axial movement of the plunger at the standard testing speed of 1.3 mm/min. While the specimen is exposed to a compressive load and gradually undergoes deformation, shear failure is initiated depending on the printed structure of the specimen. The compressive load is applied by the Instron KN1200 testing machine, which is suitable for these types of tests and was used in this work to conduct the compressive tests until the specimen completely fails. To explore the failure mechanism, the fractured surfaces are examined by microscopy. The micrograph images of two representative specimen failures (runs 10 and 22) reveal the sliding marks clearly formed on the PC strands being stressed under shear loading. The fracture plane where the shear failure occurs has an inclination of 45° to the loading axis due to the fact that on this plane, the maximum shear stresses are induced. In both cases presented here, mobilization of sliding fracture is facilitated from the specimens’ structure due to the intrinsic shear plane built in.

3.2. Design of Experiment, Experimental Results, and Statistical Analysis

The control parameter levels are shown in Table 1, whereas Table 2 and Table 3 present for each run carried out in this study, the average response parameters values and their deviation, as they were determined in the corresponding experiments by the authors of this work within the context of this specific study. More specifically, Table 2 depicts the weight (g), printing time (s), compressive strength (MPa), compressive modulus of elasticity (MPa), and compressive toughness (MJ/m3) response parameter results, and Table 3 depicts the results for EPC (MJ), SPE (MJ/g), SPP (kW/g), area to nominal (%), and volume to nominal (%). Please refer to this work’s Supplementary Material, where an optical microscope image from each of this work’s runs is displayed.
In Figure 7, by analyzing the experimental results of this work, considering the printing time (s), the weight (g), the compressive strength (MPa), and the EPC (MJ), box plots were produced:
  • For the printing time (s), only an LT of 0.3 mm shows a compact response. All the other parameters and levels show a scatter response, indicating a strong influence on the printing time (s) response parameter.
  • For the part weight (g), an ID of 60% shows a compact response. All the other parameters and levels show a scatter response, indicating a strong influence on the part weight (g) response parameter.
  • For the compressive strength (MPa), an ID of 60% and 100%, ORA 0 deg, and RDA 0 deg show a compact response. All the other parameters and levels show a scatter response, indicating a strong influence on the compressive strength (MPa) response parameter.
  • For the EPC (MJ), LT of 0.2 mm and 0.3 mm, ID of 60%, and ORA 0 deg show a compact response. All the other parameters and levels show a scatter response, indicating a strong influence on the EPC (MJ) response parameter.
MEP plots for the response are presented in Figure 8 and Figure 9, and the Supplementary Material of this work:
  • For the printing time (s), LT (mm) is the most critical parameter (rank No. 1), and then PS (mm/s). The increase in both leads to a decrease in the printing time (s). The increase in ID increases the printing time (s). The median value of 45 deg for the ORA increases the printing time (s), while low and high values lead to reduced printing time (s) values. The remaining control parameters (BT, NT, and RDA) do not significantly affect the printing time (s) response parameter, with RDA being the least important control parameter.
  • For the part weight (g), ID is the most important control parameter (rank No. 1). The increase in ID increases the part weight (g). The rank No. 2 control parameter is ORA, with the increase in the control parameter decreasing the part weight (g). The remaining control parameters (PS, LT, RDA, NT, and BT) do not significantly affect the part weight (g) response parameter, with BT being the least important control parameter.
  • For the compressive strength (MPa), ID (%) is the most critical parameter (rank No. 1), and then ORA (deg). The increase in ID leads to an increase in compressive strength (MPa). Higher compressive strength (MPa) strength values are achieved with low ORA and RDA values. The increase in LT (mm) decreases compressive strength (MPa). The remaining control parameters (PS, NT, and BT) do not significantly affect the compressive strength (MPa) response parameter, with PS being the least important control parameter.
  • For the EPC (MJ), LT (mm) is the most critical parameter (rank No. 1), and then ID (%). Higher LT (mm) values decrease the EPC (MJ) values, while low ID (%) values achieve the same effect. The increase in ORA (deg) increases the EPC (MJ). For the RDA (deg) control parameter, the median values resulted in higher EPC (MJ) values. Higher PS (mm/s) values decrease the EPC (MJ). Lower BT (°C) values also decrease the EPC (MJ). Only the NT (°C) control parameter had no significant effect on this metric, and it was at the same time the least important control parameter for the EPC (MJ) response parameter.
MEPs are not sufficient to depict the interaction between the control parameters, and thus, interaction plots were formed (Figure 10, additional interaction plots are available in the Supplementary Material of this work) to depict such connections. For the compressive strength (MPa) and the EPC (MJ), the interaction plot findings are depicted in Table 4.

3.3. Regression Analysis

The reduced quadratic regression model (RQRM) for each response is calculated:
Y k = a k + i = 1 n b i , k x i + i = 1 n c i , k x i 2 + e k
where k represents the quality output (e.g., weight, printing time, compressive strength, compressive modulus of elasticity, compressive toughness, EPC, SPE), a is a constant value, b are the coefficients of the linear terms, c are the coefficients of the quadratic terms, e is the error, and xi refers to the seven (n = 7) control parameters, i.e., the orientation angle, raster deposition angle, layer thickness, infill density, printing speed, nozzle temperature, and bed temperature.
The quadratic regression model (QRM) for each response is calculated:
Y k = a i , k + i = 1 n b i , k x i + i = 1 n c i , k x i 2 + i j d ij , k x i x j + e k
where k represents the response output (e.g., SPP, area to nominal, volume to nominal), a is a constant value, b refers to the coefficients of the linear terms, c refers to the coefficients of the quadratic terms, d refers to the coefficients of the two-way interaction terms, e is the error, and xi refers to the seven (n = 7) control parameters, i.e., the orientation angle, raster deposition angle, layer thickness, infill density, printing speed, nozzle temperature, and bed temperature.
The regression analysis for the metrics evaluated in this work is presented in corresponding tables provided in the Supplementary Material of this work (one for each metric). The corresponding predicting models formed for each metric are shown below (Equations (8)–(17)):
  • Weight (g): the F-value is 60.09 (>4), and the P-value is almost zero. The regression values are higher than 84.20%, indicating that model (8) is sufficient for the prediction of this specific metric.
  • Printing time (s): the F-value is 47.63 (>4), and the P-value is almost zero. The regression values are higher than 80.70%, indicating that the model (9) is sufficient for the prediction of this specific metric.
  • Compression strength (MPa): the F-value is 110.32 (>4), and the P-value is almost zero. The regression values are higher than 90.88%, indicating that the model (10) is sufficient for the prediction of this specific metric.
  • Compression modulus of elasticity (MPa): the F-value is 111.95 (>4), and the P-value is almost zero. The regression values are higher than 91.00%, indicating that the model (11) is sufficient for the prediction of this specific metric.
  • Compression toughness (MJ/m3): the F-value is 103.93 (>4), and the P-value is almost zero. The regression values are higher than 90.36%, indicating that the model (12) is sufficient for the prediction of this specific metric.
  • EPC (MJ): the F-value is 61.28 (>4), and the P-value is almost zero. The regression values are higher than 84.47%, indicating that the model (13) is sufficient for the prediction of this specific metric.
  • SPE (MJ/g): the F-value is 48.28 (>4), and the P-value is almost zero. The regression values are higher than 80.92%, indicating that the model (14) is sufficient for the prediction of this specific metric.
  • SPP (kW/g): the F-value is 5.03 (>4), and the P-value is almost zero. The regression values are higher than 24.56% (15). These results are marginal, and the prediction model accuracy is expected to be low for this specific metric.
  • Area2Nom (%): the F-value is 20.67 (>4), and the P-value is almost zero. The regression values are higher than 69.76% (16). These results are lower compared to the other metrics (except SPP). Although they are highly acceptable, these results indicate lower reliability in the prediction of the specific metric.
  • Vol2Nom (%): the F-value is 18.00 (>4), and the P-value is almost zero. The regression values are higher than 66.52% (17). These results are lower compared to the other metrics (except SPP). Although they are highly acceptable, these results indicate lower reliability in the prediction of the specific metric.
Weight = 23.2 0.00605 × ORA 0.00951 × RDA + 2.15 × LT + 0.0312 × ID + 0.0894 × PS + 0.179 × NT + 0.0428 × BT 0.000046 × ORA 2 + 0.000124 × RDA 2 10.60 × LT 2 + 0.000188 × ID 2 0.001133 × PS 2 0.000327 × NT 2 0.000223 × BT 2
PrintingTime = 5241 + 77.82 × ORA + 9.30 × RDA 71541 × LT + 72.0 × ID 123.4 × PS + 191 × NT 195.1 × BT 0.8477 × ORA 2 0.1187 × RDA 2 + 129,172 × LT 2 0.233 × ID 2 + 0.739 × PS 2 0.360 × NT 2 + 1.032 × BT 2
sB = 362 0.7102 × ORA + 0.1001 × RDA 30.4 × LT + 0.525 × ID 0.063 × PS 2.94 × NT + 1.246 × BT + 0.005843 × ORA 2 0.002439 × RDA 2 + 17.0 × LT 2 + 0.00118 × ID 2 + 0.00010 × PS 2 + 0.00531 × NT 2 0.00587 × BT 2
E = 7277 19.037 × ORA 1.651 × RDA 3219 × LT + 45.90 × ID 0.08 × PS 69.3 × NT + 36.56 × BT + 0.16808 × ORA 2 0.01567 × RDA 2 + 7537 × LT 2 0.2284 × ID 2 0.0058 × PS 2 + 0.1237 × NT 2 0.1777 × BT 2
Toughness = 50.5 0.04132 × ORA + 0.05584 × RDA + 24.67 × LT 0.1919 × ID 0.0065 × PS 0.349 × NT + 0.1129 × BT + 0.000175 × ORA 2 0.000633 × RDA 2 70.2 × LT 2 + 0.001955 × ID 2 + 0.000047 × PS 2 + 0.000622 × NT 2 0.000519 × BT 2
EPC = 18.3 + 0.01296 × ORA + 0.03477 × RDA 44.70 × LT 0.0316 × ID 0.0203 × PS + 0.156 × NT + 0.0564 × BT 0.000037 × ORA 2 0.000390 × RDA 2 + 85.72 × LT 2 + 0.000400 × ID 2 + 0.000048 × PS 2 0.000286 × NT 2 0.000218 × BT 2
SPE = 1.45 + 0.001552 × ORA + 0.004148 × RDA 5.659 × LT 0.00392 × ID 0.00621 × PS + 0.0130 × NT + 0.00990 × BT 0.000004 × ORA 2 0.000047 × RDA 2 + 11.09 × LT 2 + 0.000043 × ID 2 + 0.000049 × PS 2 0.000023 × NT 2 0.000042 × BT 2
SPP = 0.564 0.000419 × ORA 0.00063 × RDA + 1.342 × LT + 0.00022 × ID + 0.000721 × PS + 0.00284 × NT + 0.00131 × BT + 0.000007 × ORA 2 0.000006 × RDA 2 + 0.835 × LT 2 0.000000 × ID 2 0.000002 × PS 2 0.000003 × NT 2 0.000007 × BT 2 + 0.000007 × RDA × PS 0.000000 × RDA × NT + 0.000010 × RDA × BT 0.00267 × LT × PS 0.00578 × LT × NT
Area 2 Nom = 0.26 + 0.000576 × ORA + 0.01594 × RDA 0.616 × LT + 0.00512 × ID + 0.00429 × PS 0.00432 × NT + 0.01101 × BT 0.000005 × ORA 2 0.000010 × RDA 2 1.531 × LT 2 0.000025 × ID 2 0.000037 × PS 2 + 0.000010 × NT 2 0.000040 × BT 2 0.000020 × RDA × PS 0.000034 × RDA × NT 0.000042 × RDA × BT + 0.00226 × LT × PS + 0.00452 × LT × NT
Volume 2 Nom = 0.25 + 0.000988 × ORA + 0.01664 × RDA 0.90 × LT + 0.00527 × ID + 0.00470 × PS 0.00452 × NT + 0.01131 × BT 0.000009 × ORA 2 0.000011 × RDA 2 1.578 × LT 2 0.000026 × ID 2 0.000041 × PS 2 + 0.000010 × NT 2 0.000041 × BT 2 0.000021 × RDA × PS 0.000035 × RDA × NT 0.000045 × RDA × BT + 0.00231 × LT × PS + 0.00559 × LT × NT
To identify the statistically important parameters for the metrics of this work, Pareto charts were formed, and they are presented in Figure 11 and Figure 12 below. Additional Pareto charts are available in the Supplementary Material of this work. In the Pareto charts, the parameters that cross the 1.98 margin are the statistically important parameters for each metric. On the right side of each Pareto chart, a graph shows a comparison between the experimental and the predicted values for each specific metric. In each such graph, the mean absolute percentage error (MAPE) [77] and the Durbin–Watson factor, which is a measurement indicator of the autocorrelation in the residuals, are presented. In Figure 11 and Figure 12 below, the following can be observed:
  • Figure 11a (printing time—s): the statistically important parameters are ORA, ORA2, LT, LT2, PS, PS2, NT, and NT2. The MAPE is 22.99%, which is an acceptable result. The Durbin–Watson factor is 1.53, showing a positive autocorrelation of the prediction residuals.
  • Figure 11b (part weight—g): the statistically important parameters are RDA, RDA2, PS, and PS2. The MAPE is 4.47%, which is a very acceptable result, verifying the reliability of the model. The Durbin–Watson factor is 0.96, showing a positive autocorrelation of the prediction residuals.
  • Figure 12a (compressive strength—MPa): the statistically important parameters are ORA, ORA2, RDA2, NT, NT2, BT, and BT2. The MAPE is 8.65%, which is a very acceptable result, verifying the reliability of the model. The Durbin–Watson factor is 1.18, showing a positive autocorrelation of the prediction residuals.
  • Figure 12b (EPC-MJ): statistically important parameters are ORA, RDA, RDA2, LT, and LT2. The MAPE is 35.86%, which is an acceptable result. The Durbin–Watson factor is 1.42, showing a positive autocorrelation of the prediction residuals.
Figure 13 shows, for some critical metrics of this work, three-dimensional surface graphs of the control parameters (rank No. 1 and 2) vs. the corresponding metric.

3.4. Confirmation Experiments

To verify the regression analysis and the prediction models forms, two confirmation runs were conducted, with five replicas each (runs 28 and 29). Table 5 presents the levels for the control parameters in each run. The experimental results for the metrics of this work in each run (average value and deviation) are presented in Table 6 and Table 7 (the confirmation runs results are analytically accessible in the Supplementary Material of this work). The accuracy of the prediction models in the conformation runs is depicted in Table 8, in which the actual and the predicted values for the metrics of this work are presented along with the deviation (error) between them. An eminent accuracy of the prediction models was found for run 28, with 5.62% and 8.52% deviation between the actual and the calculated compression strength (MPa) and EPC (MJ) values, respectively. For run 29, the accuracy of compression strength (MPa) was not as high as in run 28 (18.03%), while for the EPC (MJ) metric, the model failed to predict the expected value. This indicates that the control parameter levels were outside the limits for the specific metric. The acceptable limits of the prediction models were not within the scope of the current work. It should be noted that the compression strength and the EPC are the two main metrics related to the aims of this work. Therefore, they were the metrics that were more critical to be evaluated in the confirmation tests. Confirmation tests were conducted with the more suitable control parameter levels for each metric. Since the optimum parameter levels for each metric studied in this work differ, the extent of this work would be significantly increased if more metrics were considered in the confirmation tests.

4. Conclusions

In this study, the energy consumption, compressive strength, and two quality metrics of parts made with the high-performance PC polymer with MEX 3D printing were investigated compared to seven continuous, machine-independent, generic 3D printing settings. The study of energy consumption is directly related to the sustainability of the MEX 3D printing process, which is already a critical aspect of AM. Additionally, a thorough investigation was conducted on the performance of MEX 3D-printed parts when subjected to compressive load. The two quality parameters studied are related to the dimensional accuracy of the MEX 3D-printed parts, and it was shown how it is affected by seven different 3D printing settings when varying their values within limits determined by the literature review conducted in this work.
DOE and regression analysis ranked the importance of each 3D printing setting on the energy, mechanical properties under compression loading, and dimensional accuracy metrics studied in this work. Additionally, prediction models were compiled for each one of the ten metrics considered in this work. The accuracy of the prediction models was sufficient for direct industrial use in all metrics except the SPP, which showed reduced accuracy. A reduced accuracy of about 60% was also found for the two quality metrics in this work, i.e., actual to nominal for the samples’ cross-section area and volume. The validity of these prediction models was verified with two confirmation runs. In the confirmation runs, the predicted values for the compression strength were very accurate in both runs. In the second confirmation run, the model failed to predict the EPC value, showing that the prediction models function properly within limits for the control parameter levels.
No set of 3D printing settings optimizes both the energy and the mechanical properties metrics, indicating that priorities should be made for the optimization of the desirable metrics, with the models that can predict the expected values for the remaining metrics. Infill density was the most important parameter affecting the compression strength of the PC MEX 3D-printed parts. It was ranked No. 1 in importance, and high ID values (i.e., 100%) are recommended to achieve increased compression strength values. Low ORA values (i.e., 0 deg) also contribute to this result. Printing speed was the least important parameter for the compression strength metric. The layer thickness was the most important parameter (ranked as No. 1) affecting the energy consumption of the PC MEX 3D-printed parts. High values of LT (i.e., 0.3 mm) significantly reduce the energy consumption when 3D printing PC parts with the MEX process. Infill density also significantly affects energy consumption, and low values (i.e., 60%) are recommended when the energy requirements need to be minimum. Nozzle temperature was the least important parameter for this metric. To broaden the range of applicability for the current study’s findings, future work may expand on it by examining control parameter values in a new value range.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jmmp7010038/s1.

Author Contributions

Conceptualization, N.V. and M.P.; methodology, N.V. and M.P.; software, N.M., C.N.D. and D.S.; validation, N.V., D.S. and E.K.; formal analysis, N.M., C.N.D. and E.K.; investigation, D.S.; resources, N.V. and C.N.D.; data curation, N.M., C.N.D. and E.K.; writing—original draft preparation, M.P.; writing—review and editing, M.P.; visualization, N.V.; supervision, N.V.; project administration, M.P.; funding acquisition, N.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

3DP3D Printing
ABSAcrylonitrile Butadiene Styrene
AMAdditive Manufacturing
ANOVAAnalysis of Variances
BTBed Temperature
DFDegrees of Freedom
DOEDesign of Experiment
DSCDifferential Scanning Calorimetry
ETensile Modulus of Elasticity
EPCEnergy Printing Consumption
FFFFused Filament Fabrication
IDInfill Density
LTLayer Thickness
MEPMain Effect Plot
MEXMaterial Extrusion
NTNozzle Temperature
ORAOrientation Angle
PAPolyamide
PCPolycarbonate
PEEKPolyether-ether-ketone
PLAPolylactic Acid
PTPrinting Time
PSPrinting Speed
RDARaster Deposition Angle
QRMQuadratic Regression Model
RQRMReduced Quadratic Regression Model
sBCompression strength
SEMScanning Electron Microscopy
SPESpecific Printing Energy
SPPSpecific Printing Power
TgGlass Transition Temperature
TGAThermogravimetric Analysis

References

  1. Freitas, D.; Almeida, H.A.; Bártolo, H.; Bártolo, P.J. Sustainability in extrusion-based additive manufacturing technologies. Prog. Addit. Manuf. 2016, 1, 65–78. [Google Scholar] [CrossRef]
  2. Liu, Z.; Jiang, Q.; Zhang, Y.; Li, T.; Zhang, H.C. Sustainability of 3D printing: A critical review and recommendations. In Proceedings of the ASME 2016 11th International Manufacturing Science and Engineering Conference, Blacksburg, VI, USA, 27 June–1 July 2016; Volume 2, pp. 1–8. [Google Scholar] [CrossRef]
  3. Malik, A.; Ul Haq, M.I.; Raina, A.; Gupta, K. 3D printing towards implementing Industry 4.0: Sustainability aspects, barriers and challenges. Ind. Rob. 2022, 49, 491–511. [Google Scholar] [CrossRef]
  4. Jiang, J.; Fu, Y.F. A short survey of sustainable material extrusion additive manufacturing. Aust. J. Mech. Eng. 2020, 1–10. [Google Scholar] [CrossRef]
  5. Dey, A.; Eagle, I.N.R.; Yodo, N. A review on filament materials for fused filament fabrication. J. Manuf. Mater. Process. 2021, 5, 69. [Google Scholar] [CrossRef]
  6. Suárez, L.; Domínguez, M. Sustainability and environmental impact of fused deposition modelling (FDM) technologies. Int. J. Adv. Manuf. Technol. 2020, 106, 1267–1279. [Google Scholar] [CrossRef]
  7. Khosravani, M.R.; Reinicke, T. On the environmental impacts of 3D printing technology. Appl. Mater. Today 2020, 20, 100689. [Google Scholar] [CrossRef]
  8. Stoof, D.; Pickering, K. Sustainable composite fused deposition modelling filament using recycled pre-consumer polypropylene. Compos. Part B Eng. 2018, 135, 110–118. [Google Scholar] [CrossRef]
  9. Fico, D.; Rizzo, D.; Casciaro, R.; Corcione, C.E. A Review of Polymer-Based Materials for Fused Filament Recycled Materials. Polymers 2022, 14, 465. [Google Scholar] [CrossRef]
  10. Kumar, S.; Kazancoglu, Y.; Deniz, M.; Top, N.; Sahin, I. Optimizing fused deposition modelling parameters based on the design for additive manufacturing to enhance product sustainability. Comput. Ind. 2023, 145, 103833. [Google Scholar] [CrossRef]
  11. Ramesh, P.; Vinodh, S. Analysis of factors influencing energy consumption of material extrusion-based additive manufacturing using interpretive structural modelling. Rapid Prototyp. J. 2021, 27, 1363–1377. [Google Scholar] [CrossRef]
  12. Abeykoon, C.; McMillan, A.; Nguyen, B.K. Energy efficiency in extrusion-related polymer processing: A review of state of the art and potential efficiency improvements. Renew. Sustain. Energy Rev. 2021, 147, 111219. [Google Scholar] [CrossRef]
  13. Kim, K.; Noh, H.; Park, K.; Jeon, H.W.; Lim, S. Characterization of power demand and energy consumption for fused filament fabrication using CFR-PEEK. Rapid Prototyp. J. 2022, 28, 1394–1406. [Google Scholar] [CrossRef]
  14. Hassan, M.R.; Jeon, H.W.; Kim, G.; Park, K. The effects of infill patterns and infill percentages on energy consumption in fused filament fabrication using CFR-PEEK. Rapid Prototyp. J. 2021, 27, 1886–1899. [Google Scholar] [CrossRef]
  15. De Bernardez, L.; Campana, G.; Mele, M.; Sanguineti, J.; Sandre, C.; Mur, S.M. Effects of infill patterns on part performances and energy consumption in acrylonitrile butadiene styrene fused filament fabrication via industrial-grade machine. Prog. Addit. Manuf. 2022, 1–13. [Google Scholar] [CrossRef]
  16. El youbi El idrissi, M.A.; Laaouina, L.; Jeghal, A.; Tairi, H.; Zaki, M. Energy Consumption Prediction for Fused Deposition Modelling 3D Printing Using Machine Learning. Appl. Syst. Innov. 2022, 5, 86. [Google Scholar] [CrossRef]
  17. Vidakis, N.; Kechagias, J.D.; Petousis, M.; Vakouftsi, F.; Mountakis, N. The effects of FFF 3D printing parameters on energy consumption. Mater. Manuf. Process. 2022, 1–18. [Google Scholar] [CrossRef]
  18. Peng, T. Energy Modelling for FDM 3D Printing from a Life Cycle Perspective. Int. J. Manuf. Res. 2017, 11, 1. [Google Scholar] [CrossRef]
  19. Hassan, M.R.; Noh, H.; Park, K.; Jeon, H.W. Simulating energy consumption based on material addition rates for material extrusion of CFR-PEEK: A trade-off between energy costs and cycle time. Int. J. Adv. Manuf. Technol. 2022, 120, 4597–4616. [Google Scholar] [CrossRef]
  20. Tian, W.; Ma, J.; Alizadeh, M. Energy consumption optimization with geometric accuracy consideration for fused filament fabrication processes. Int. J. Adv. Manuf. Technol. 2019, 103, 3223–3233. [Google Scholar] [CrossRef]
  21. Kausar, A. A review of filled and pristine polycarbonate blends and their applications. J. Plast. Film Sheeting 2018, 34, 60–97. [Google Scholar] [CrossRef]
  22. Park, S.J.; Lee, J.E.; Lee, H.B.; Park, J.; Lee, N.K.; Son, Y.; Park, S.H. 3D printing of bio-based polycarbonate and its potential applications in ecofriendly indoor manufacturing. Addit. Manuf. 2020, 31, 100974. [Google Scholar] [CrossRef]
  23. Reich, M.J.; Woern, A.L.; Tanikella, N.G.; Pearce, J.M. Mechanical properties and applications of recycled polycarbonate particle material extrusion-based additive manufacturing. Materials 2019, 12, 1642. [Google Scholar] [CrossRef] [PubMed]
  24. Mallakpour, S.; Hussain, C.M. Handbook of Consumer Nanoproducts. Handb. Consum. Nanoproducts 2022, 257–281. [Google Scholar] [CrossRef]
  25. Pedrosa, P.; Alves, E.; Barradas, N.P.; Fiedler, P.; Haueisen, J.; Vaz, F.; Fonseca, C. TiNx coated polycarbonate for bio-electrode applications. Corros. Sci. 2012, 56, 49–57. [Google Scholar] [CrossRef]
  26. Pande, S.; Chaudhary, A.; Patel, D.; Singh, B.P.; Mathur, R.B. Mechanical and electrical properties of multiwall carbon nanotube/polycarbonate composites for electrostatic discharge and electromagnetic interference shielding applications. RSC Adv. 2014, 4, 13839–13849. [Google Scholar] [CrossRef]
  27. dos Anjos, E.G.R.; Moura, N.K.; Antonelli, E.; Baldan, M.R.; Gomes, N.A.S.; Braga, N.F.; Santos, A.P.; Rezende, M.C.; Pessan, L.A.; Passador, F.R. Role of adding carbon nanotubes in the electric and electromagnetic shielding behaviors of three different types of graphene in hybrid nanocomposites. J. Thermoplast. Compos. Mater. 2022, 1–27. [Google Scholar] [CrossRef]
  28. Chen, F.; Bao, Y.; Zhang, J.; Yang, M.; Ruan, M.; Feng, W.; Jiang, Y.; Li, M.; Chen, Y. Comparative study on the mechanical and thermal properties of polycarbonate composites reinforced by KH570/SA/SDBS modified wollastonite fibers. Polym. Compos. 2022, 43, 8125–8135. [Google Scholar] [CrossRef]
  29. Wijerathne, D.; Gong, Y.; Afroj, S.; Karim, N.; Abeykoon, C. Mechanical and thermal properties of graphene nanoplatelets-reinforced recycled polycarbonate composites. Int. J. Light. Mater. Manuf. 2023, 6, 117–128. [Google Scholar] [CrossRef]
  30. Zhang, J.; Koubaa, A.; Xing, D.; Wang, H.; Tao, Y.; Wang, X.M.; Li, P. Fire Behavior and Failure Model of Multilayered Wood Flour/HDPE/Polycarbonate Composites with a Sandwich Structure. Polymers 2022, 14, 2833. [Google Scholar] [CrossRef]
  31. Blanco, I.; Cicala, G.; Ognibene, G.; Rapisarda, M.; Recca, A. Thermal properties of polyetherimide/polycarbonate blends for advanced applications. Polym. Degrad. Stab. 2018, 154, 234–238. [Google Scholar] [CrossRef]
  32. Kim, M.K.; Kim, H.I.; Nam, J.D.; Suhr, J. Polyamide-nylon 6 particulate polycarbonate composites with outstanding energy-absorbing properties. Polymer 2022, 254, 125082. [Google Scholar] [CrossRef]
  33. Costa, P.; Dios, J.R.; Cardoso, J.; Campo, J.J.; Tubio, C.R.; Gon alves, B.F.; Castro, N.; Lanceros-Méndez, S. Polycarbonate based multifunctional self-sensing 2D and 3D printed structures for aeronautic applications. Smart Mater. Struct. 2021, 30, 085032. [Google Scholar] [CrossRef]
  34. Alaboodi, A.S.; Sivasankaran, S. Experimental design and investigation on the mechanical behavior of novel 3D printed biocompatibility polycarbonate scaffolds for medical applications. J. Manuf. Process. 2018, 35, 479–491. [Google Scholar] [CrossRef]
  35. Gómez-Gras, G.; Abad, M.D.; Pérez, M.A. Mechanical performance of 3D-printed biocompatible polycarbonate for biomechanical applications. Polymers 2021, 13, 3669. [Google Scholar] [CrossRef] [PubMed]
  36. Tichy, A.; Simkova, M.; Schweiger, J.; Bradna, P.; Güth, J.-F. Release of Bisphenol A from Milled and 3D-Printed Dental Polycarbonate Materials. Materials 2021, 14, 5868. [Google Scholar] [CrossRef] [PubMed]
  37. Bahar, A.; Belhabib, S.; Guessasma, S.; Benmahiddine, F.; Hamami, A.E.; Belarbi, R. Mechanical and Thermal Properties of 3D Printed Polycarbonate. Energies 2022, 15, 3686. [Google Scholar] [CrossRef]
  38. Cantrell, J.T.; Rohde, S.; Damiani, D.; Gurnani, R.; DiSandro, L.; Anton, J.; Young, A.; Jerez, A.; Steinbach, D.; Kroese, C.; et al. Experimental characterization of the mechanical properties of 3D-printed ABS and polycarbonate parts. Rapid Prototyp. J. 2017, 23, 811–824. [Google Scholar] [CrossRef]
  39. Rohde, S.; Cantrell, J.; Jerez, A.; Kroese, C.; Damiani, D.; Gurnani, R.; DiSandro, L.; Anton, J.; Young, A.; Steinbach, D.; et al. Experimental Characterization of the Shear Properties of 3D–Printed ABS and Polycarbonate Parts. Exp. Mech. 2018, 58, 871–884. [Google Scholar] [CrossRef]
  40. Vidakis, N.; Petousis, M.; Kechagias, J.D. A comprehensive investigation of the 3D printing parameters’ effects on the mechanical response of polycarbonate in fused filament fabrication. Prog. Addit. Manuf. 2022, 7, 713–722. [Google Scholar] [CrossRef]
  41. Vidakis, N.; Petousis, M.; Kechagias, J.D. Parameter effects and process modelling of Polyamide 12 3D-printed parts strength and toughness. Mater. Manuf. Process. 2022, 37, 1358–1369. [Google Scholar] [CrossRef]
  42. Kechagias, J.D.; Vidakis, N.; Petousis, M. Parameter effects and process modeling of FFF-TPU mechanical response. Mater. Manuf. Process. 2021, 1–11. [Google Scholar] [CrossRef]
  43. Kechagias, J.D.; Vidakis, N.; Petousis, M.; Mountakis, N. A multi-parametric process evaluation of the mechanical response of PLA in FFF 3D printing. Mater. Manuf. Process. 2022, 1–13. [Google Scholar] [CrossRef]
  44. Vidakis, N.; Petousis, M.; Mountakis, N.; Maravelakis, E.; Zaoutsos, S.; Kechagias, J.D. Mechanical response assessment of antibacterial PA12/TiO2 3D printed parts: Parameters optimization through artificial neural networks modeling. Int. J. Adv. Manuf. Technol. 2022, 121, 785–803. [Google Scholar] [CrossRef] [PubMed]
  45. Vidakis, N.; Petousis, M.; Korlos, A.; Velidakis, E.; Mountakis, N.; Charou, C.; Myftari, A. Strain Rate Sensitivity of Polycarbonate and Thermoplastic Polyurethane for Various 3D Printing Temperatures and Layer Heights. Polymers 2021, 13, 2752. [Google Scholar] [CrossRef]
  46. Vries, H.D.; Engelen, R.; Janssen, E. Impact Strength of 3d-Printed Polycarbonate. Facta Univ. Ser. Electron. Energetics 2020, 33, 105–117. [Google Scholar] [CrossRef]
  47. Maurya, N.K.; Rastogi, V.; Singh, P. Investigation of dimensional accuracy and international tolerance grades of 3D printed polycarbonate parts. Mater. Today Proc. 2019, 25, 537–543. [Google Scholar] [CrossRef]
  48. Gupta, A.; Fidan, I.; Hasanov, S.; Nasirov, A. Processing, mechanical characterization, and micrography of 3D-printed short carbon fiber reinforced polycarbonate polymer matrix composite material. Int. J. Adv. Manuf. Technol. 2020, 107, 3185–3205. [Google Scholar] [CrossRef]
  49. Jahangir, M.N.; Billah, K.M.M.; Lin, Y.; Roberson, D.A.; Wicker, R.B.; Espalin, D. Reinforcement of material extrusion 3D printed polycarbonate using continuous carbon fiber. Addit. Manuf. 2019, 28, 354–364. [Google Scholar] [CrossRef]
  50. Vidakis, N.; Petousis, M.; Grammatikos, S.; Papadakis, V.; Korlos, A.; Mountakis, N. High Performance Polycarbonate Nanocomposites Mechanically Boosted with Titanium Carbide in Material Extrusion Additive Manufacturing. Nanomaterials 2022, 12, 1068. [Google Scholar] [CrossRef]
  51. Vidakis, N.; Petousis, M.; Mangelis, P.; Maravelakis, E.; Mountakis, N.; Papadakis, V.; Neonaki, M.; Thomadaki, G. Thermomechanical Response of Polycarbonate/Aluminum Nitride Nanocomposites in Material Extrusion Additive Manufacturing. Materials 2022, 15, 8806. [Google Scholar] [CrossRef]
  52. Vidakis, N.; Petousis, M.; Mountakis, N.; Grammatikos, S.; Papadakis, V.; Kechagias, J.D.; Das, S.C. On the thermal and mechanical performance of Polycarbonate/Titanium Nitride nanocomposites in Material Extrusion Additive Manufacturing. Compos. Part C Open Access 2022, 8, 100291. [Google Scholar] [CrossRef]
  53. Petousis, M.; Vidakis, N.; Mountakis, N.; Grammatikos, S.; Papadakis, V.; David, C.N.; Moutsopoulou, A.; Das, S.C. Silicon Carbide Nanoparticles as a Mechanical Boosting Agent in Material Extrusion 3D-Printed Polycarbonate. Polymers 2022, 14, 3492. [Google Scholar] [CrossRef]
  54. Kodali, D.; Umerah, C.O.; Idrees, M.O.; Jeelani, S.; Rangari, V.K. Fabrication and characterization of polycarbonate-silica filaments for 3D printing applications. J. Compos. Mater. 2021, 55, 4575–4584. [Google Scholar] [CrossRef]
  55. Vidakis, N.; Petousis, M.; Velidakis, E.; Spiridaki, M.; Kechagias, J. Mechanical Performance of Fused Filament Fabricated and 3D-Printed Polycarbonate Polymer and Polycarbonate/Cellulose Nanofiber Nanocomposites. Fibers 2021, 9, 74. [Google Scholar] [CrossRef]
  56. Verma, N.; Aiswarya, S.; Banerjee, S.S. Development of material extrusion 3D printable ABS/PC polymer blends: Influence of styrene–isoprene–styrene copolymer on printability and mechanical properties. Polym. Technol. Mater. 2022, 1–14. [Google Scholar] [CrossRef]
  57. Kumar, P.; Gupta, P.; Singh, I. Performance Analysis of Acrylonitrile–Butadiene–Styrene–Polycarbonate Polymer Blend Filament for Fused Deposition Modeling Printing Using Hybrid Artificial Intelligence Algorithms. J. Mater. Eng. Perform. 2022, 1–14. [Google Scholar] [CrossRef]
  58. Yap, Y.L.; Toh, W.; Koneru, R.; Lin, K.; Yeoh, K.M.; Lim, C.M.; Lee, J.S.; Plemping, N.A.; Lin, R.; Ng, T.Y.; et al. A non-destructive experimental-cum-numerical methodology for the characterization of 3D-printed materials—Polycarbonate-acrylonitrile butadiene styrene (PC-ABS). Mech. Mater. 2019, 132, 121–133. [Google Scholar] [CrossRef]
  59. Kumar, M.; Ramakrishnan, R.; Omarbekova, A. 3D printed polycarbonate reinforced acrylonitrile–butadiene–styrene composites: Composition effects on mechanical properties, micro-structure and void formation study. J. Mech. Sci. Technol. 2019, 33, 5219–5226. [Google Scholar] [CrossRef]
  60. Sumalatha, M.; Rao, J.N.M.; Reddy, B.S. Optimization Of Process Parameters In 3d Printing-Fused Deposition Modeling Using Taguchi Method. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1112, 012009. [Google Scholar] [CrossRef]
  61. Salokhe, O.A.; Shaikh, A.M. Study of Fused Deposition Modeling Process Parameters for Polycarbonate/Acrylonitrile Butadiene Styrene Blend Material using Taguchi Method. Int. Res. J. Eng. Technol. 2019, 4244–4249. [Google Scholar]
  62. Domingo-Espin, M.; Borros, S.; Agullo, N.; Garcia-Granada, A.A.; Reyes, G. Influence of Building Parameters on the Dynamic Mechanical Properties of Polycarbonate Fused Deposition Modeling Parts. 3D Print. Addit. Manuf. 2014, 1, 70–77. [Google Scholar] [CrossRef]
  63. Domingo-Espin, M.; Puigoriol-Forcada, J.M.; Garcia-Granada, A.-A.; Llumà, J.; Borros, S.; Reyes, G. Mechanical property characterization and simulation of fused deposition modeling Polycarbonate parts. Mater. Des. 2015, 83, 670–677. [Google Scholar] [CrossRef]
  64. Vidakis, N.; Petousis, M.; Mountakis, N.; Kechagias, J.D. Optimization of Friction Stir Welding Parameters in Hybrid Additive Manufacturing: Weldability of 3D-Printed Poly (methyl methacrylate) Plates. J. Manuf. Mater. Process. 2022, 6, 77. [Google Scholar] [CrossRef]
  65. Kechagias, J.D.; Vidakis, N.; Ninikas, K.; Petousis, M.; Vaxevanidis, N.M. Hybrid 3D printing of multifunctional polylactic acid/carbon black nanocomposites made with material extrusion and post-processed with CO2 laser cutting. Int. J. Adv. Manuf. Technol. 2022, 124, 1843–1861. [Google Scholar] [CrossRef]
  66. Blumenthal, W.R. Influence of Temperature and Strain Rate on the Compressive Behavior of PMMA and Polycarbonate Polymers. Int. J. Solids Struct. 2003, 665, 665–668. [Google Scholar] [CrossRef]
  67. Li, S.; Liu, Z.; Shim, V.P.W.; Guo, Y.; Sun, Z.; Li, X.; Wang, Z. In-plane compression of 3D-printed self-similar hierarchical honeycombs–Static and dynamic analysis. Thin-Walled Struct. 2020, 157, 106990. [Google Scholar] [CrossRef]
  68. Ramírez-Revilla, S.; Camacho-Valencia, D.; Gonzales-Condori, E.G.; Márquez, G. Evaluation and comparison of the degradability and compressive and tensile properties of 3D printing polymeric materials: PLA, PETG, PC, and ASA. MRS Commun. 2022, 1–8. [Google Scholar] [CrossRef]
  69. Siviour, C.R.; Walley, S.M.; Proud, W.G.; Field, J.E. The high strain rate compressive behaviour of polycarbonate and polyvinylidene difluoride. Polymer 2005, 46, 12546–12555. [Google Scholar] [CrossRef]
  70. Vidakis, N.; Petousis, M.; Vairis, A.; Savvakis, K.; Maniadi, A. On the compressive behavior of an FDM Steward Platform part. J. Comput. Des. Eng. 2017, 4, 339–346. [Google Scholar] [CrossRef]
  71. Phadke, M.S. Quality Engineering Using Robust Design, 1st ed.; Prentice Hall PTR: Hillsdale, NJ, USA, 1995; ISBN 0137451679. [Google Scholar]
  72. Capote, G.A.M.; Rudolph, N.M.; Osswald, P.V.; Osswald, T.A. Failure surface development for ABS fused filament fabrication parts. Addit. Manuf. 2019, 28, 169–175. [Google Scholar] [CrossRef]
  73. Mustafa, A.; Aloyaydi, B.; Subbarayan, S.; Al-Mufadi, F.A. Mechanical properties enhancement in composite material structures of poly-lactic acid/epoxy/milled glass fibers prepared by fused filament fabrication and solution casting. Polym. Compos. 2021, 42, 6847–6866. [Google Scholar] [CrossRef]
  74. Daniel, I.M. Failure mechanisms in thick composites under compressive loading. Compos. Part B Eng. 1996, 27, 543–552. [Google Scholar] [CrossRef]
  75. Sun, Q.; Guo, H.; Zhou, G.; Meng, Z.; Chen, Z. Experimental and computational analysis of failure mechanisms in unidirectional carbon fi ber reinforced polymer laminates under longitudinal compression loading. Compos. Struct. 2018, 203, 335–348. [Google Scholar] [CrossRef]
  76. Lemanski, S.L.; Wang, J.; Sutcliffe, M.P.F.; Potter, K.D.; Wisnom, M.R. Composites: Part A Modelling failure of composite specimens with defects under compression loading. Compos. Part A 2013, 48, 26–36. [Google Scholar] [CrossRef] [Green Version]
  77. Swamidass, P.M. (Ed.) MAPE (Mean Absolute Percentage Error)Mean Absolute Percentage Error (Mape) BT-Encyclopedia of Production and Manufacturing Management; Springer: Boston, MA, USA, 2000; p. 462. ISBN 978-1-4020-0612-8. [Google Scholar]
Figure 1. The experimental procedure followed in this work: (a) raw materials drying, (b) filament extrusion, (c) filament drying, (d) 3D printing of the samples, (e) energy monitoring during the 3D printing process, (f) sample inspection in the microscope, (g) compression test of the 3D printing samples.
Figure 1. The experimental procedure followed in this work: (a) raw materials drying, (b) filament extrusion, (c) filament drying, (d) 3D printing of the samples, (e) energy monitoring during the 3D printing process, (f) sample inspection in the microscope, (g) compression test of the 3D printing samples.
Jmmp 07 00038 g001
Figure 2. For the PC grade studied here, (a) weight loss vs. temperature graph (°C) (TGA), (b) endo graph (DSC).
Figure 2. For the PC grade studied here, (a) weight loss vs. temperature graph (°C) (TGA), (b) endo graph (DSC).
Jmmp 07 00038 g002
Figure 3. Graphic representation of the 3D printing process, defining the 3D printing parameters studied in this work. On the right side of the figure, the 3D printing parameter values and the geometry of the compression test sample fabricated following the ASTM D695 standard are presented.
Figure 3. Graphic representation of the 3D printing process, defining the 3D printing parameters studied in this work. On the right side of the figure, the 3D printing parameter values and the geometry of the compression test sample fabricated following the ASTM D695 standard are presented.
Jmmp 07 00038 g003
Figure 4. Images from the optical microscope of 3D-printed specimens, built with the various settings studied in this work. The 3D printing orientation angle that each sample was built on is presented on the left side in each case.
Figure 4. Images from the optical microscope of 3D-printed specimens, built with the various settings studied in this work. The 3D printing orientation angle that each sample was built on is presented on the left side in each case.
Jmmp 07 00038 g004
Figure 5. Specimens in the compression test at the failure stage, depicting their different behavior for the various 3D printing sets of settings.
Figure 5. Specimens in the compression test at the failure stage, depicting their different behavior for the various 3D printing sets of settings.
Jmmp 07 00038 g005
Figure 6. (a) Compression test of a 3D-printed specimen until failure (Run 24), (b) graphical representation of the compression test and microscopic inspection of the fracture surface, (c,d) micrograph images of the fracture surface of the upper and the lower part of two representative specimens (Run 10, 22) subjected to shear failure.
Figure 6. (a) Compression test of a 3D-printed specimen until failure (Run 24), (b) graphical representation of the compression test and microscopic inspection of the fracture surface, (c,d) micrograph images of the fracture surface of the upper and the lower part of two representative specimens (Run 10, 22) subjected to shear failure.
Jmmp 07 00038 g006
Figure 7. Box plots for (a) printing time (s) vs. LT (mm), PS (mm/s), and ORA (deg); (b) part weight (g) vs. ID (%), ORA (deg), and PS (mm/min); (c) compressive strength (MPa) vs. ID (%), ORA (deg), and RDA (deg); (d) EPC (MJ) vs. LT (mm), ID (%), and ORA (deg).
Figure 7. Box plots for (a) printing time (s) vs. LT (mm), PS (mm/s), and ORA (deg); (b) part weight (g) vs. ID (%), ORA (deg), and PS (mm/min); (c) compressive strength (MPa) vs. ID (%), ORA (deg), and RDA (deg); (d) EPC (MJ) vs. LT (mm), ID (%), and ORA (deg).
Jmmp 07 00038 g007
Figure 8. MEP for the printing time (s) and the part weight (g) vs. the control parameters of this work.
Figure 8. MEP for the printing time (s) and the part weight (g) vs. the control parameters of this work.
Jmmp 07 00038 g008
Figure 9. MEP for the compressive strength (MPa) and the energy (MJ) vs. the control parameters of this work.
Figure 9. MEP for the compressive strength (MPa) and the energy (MJ) vs. the control parameters of this work.
Jmmp 07 00038 g009
Figure 10. Interaction plots for the energy (MJ) and the compressive strength (MPa) vs. the control parameters of this work.
Figure 10. Interaction plots for the energy (MJ) and the compressive strength (MPa) vs. the control parameters of this work.
Jmmp 07 00038 g010
Figure 11. Pareto chart for (a) the printing time (s) and (b) the part weight (g). The statistically significant parameters are those crossing the threshold value of 1.98. On the right side, the experimental vs. the predicted values are presented.
Figure 11. Pareto chart for (a) the printing time (s) and (b) the part weight (g). The statistically significant parameters are those crossing the threshold value of 1.98. On the right side, the experimental vs. the predicted values are presented.
Jmmp 07 00038 g011
Figure 12. Pareto chart for (a) compressive strength (MPa) and the (b) energy (MJ). The statistically significant parameters are those crossing the threshold value of 1.98. On the right side, the experimental vs. the predicted values are presented.
Figure 12. Pareto chart for (a) compressive strength (MPa) and the (b) energy (MJ). The statistically significant parameters are those crossing the threshold value of 1.98. On the right side, the experimental vs. the predicted values are presented.
Jmmp 07 00038 g012
Figure 13. Three-dimensional surface graphs: (a) part weight (g) vs. ID (%) and ORA (deg), (b) compression strength (MPa) vs. ID (%) and ORA (deg), (c) energy (MJ) vs. LT (mm) and ID (%), (d) printing time (s) vs. LT (mm) and PS (mm/s), (e) compression strength (MPa) vs. RDA (deg) and LT (mm), (f) energy (MJ) vs. ORA (deg) and RDA (deg).
Figure 13. Three-dimensional surface graphs: (a) part weight (g) vs. ID (%) and ORA (deg), (b) compression strength (MPa) vs. ID (%) and ORA (deg), (c) energy (MJ) vs. LT (mm) and ID (%), (d) printing time (s) vs. LT (mm) and PS (mm/s), (e) compression strength (MPa) vs. RDA (deg) and LT (mm), (f) energy (MJ) vs. ORA (deg) and RDA (deg).
Jmmp 07 00038 g013
Table 1. Taguchi L27 design: control parameters and levels.
Table 1. Taguchi L27 design: control parameters and levels.
RunORARDALTIDPSNTBT
1000.1602026080
2000.16040280100
3000.16060300120
40450.2802026080
50450.28040280100
60450.28060300120
70900.31002026080
80900.310040280100
90900.310060300120
104500.210020280120
114500.21004030080
124500.210060260100
1345450.36020280120
1445450.3604030080
1545450.36060260100
1645900.18020280120
1745900.1804030080
1845900.18060260100
199000.38020300100
209000.38040260120
219000.3806028080
2290450.110020300100
2390450.110040260120
2490450.11006028080
2590900.26020300100
2690900.26040260120
2790900.2606028080
Table 2. Mean average values and standard deviations of measured responses for weight, printing time, compressive strength (sB), compressive modulus of elasticity (E), and compressive toughness.
Table 2. Mean average values and standard deviations of measured responses for weight, printing time, compressive strength (sB), compressive modulus of elasticity (E), and compressive toughness.
RunWeight (g)Printing Time (s)sB [MPa]E [MPa]Toughness [MJ/m3]
16.83 ± 0.217318.80 ± 1656.7547.65 ± 3.961062.93 ± 71.464.27 ± 0.14
28.33 ± 0.214326.00 ± 973.5650.72 ± 0.441132.37 ± 13.404.66 ± 0.13
36.70 ± 0.063340.20 ± 721.7451.63 ± 1.001118.10 ± 26.885.22 ± 0.49
48.33 ± 0.183621.80 ± 722.5558.88 ± 0.381191.34 ± 12.067.75 ± 0.16
58.46 ± 0.412557.00 ± 474.8660.77 ± 0.611196.78 ± 17.217.93 ± 0.13
67.68 ± 0.361868.20 ± 367.6762.84 ± 1.261155.62 ± 45.748.15 ± 0.37
79.52 ± 0.103431.20 ± 705.0964.28 ± 2.571143.33 ± 38.629.15 ± 0.39
89.47 ± 0.041953.00 ± 369.3962.77 ± 0.701166.96 ± 54.748.72 ± 0.13
99.42 ± 0.181475.00 ± 309.7362.14 ± 0.971177.49 ± 13.828.45 ± 0.27
109.55 ± 0.236906.20 ± 1201.2355.84 ± 2.71886.82 ± 95.628.34 ± 0.43
119.45 ± 0.054348.00 ± 815.6755.10 ± 1.66825.61 ± 70.707.94 ± 1.63
128.81 ± 0.073422.00 ± 691.6556.68 ± 0.98887.99 ± 36.128.95 ± 0.14
135.80 ± 0.044458.00 ± 595.6424.75 ± 0.59410.29 ± 30.553.76 ± 0.09
146.61 ± 0.062399.00 ± 429.4526.26 ± 0.44437.88 ± 39.414.16 ± 0.10
156.57 ± 0.092095.20 ± 466.2223.14 ± 1.89475.80 ± 34.233.37 ± 0.62
168.15 ± 0.449484.00 ± 170.2134.03 ± 4.25523.85 ± 91.064.85 ± 0.89
178.30 ± 0.086921.00 ± 1519.4829.35 ± 1.67556.02 ± 81.624.22 ± 0.05
188.21 ± 0.175927.80 ± 1173.1435.40 ± 0.34693.81 ± 38.125.17 ± 0.23
197.41 ± 0.042745.60 ± 559.1654.33 ± 4.881142.79 ± 23.903.65 ± 0.63
207.20 ± 0.131202.00 ± 238.0538.68 ± 9.24954.44 ± 122.893.54 ± 1.45
216.68 ± 0.101893.00 ± 331.6434.44 ± 6.39873.71 ± 111.672.92 ± 0.89
227.54 ± 0.307859.00 ± 1649.0265.23 ± 3.561138.30 ± 99.489.38 ± 0.25
239.63 ± 0.438291.80 ± 1830.2067.08 ± 3.001107.30 ± 53.799.06 ± 0.69
248.75 ± 0.403954.00 ± 802.0552.23 ± 6.07817.24 ± 157.847.00 ± 0.82
256.62 ± 0.162300.80 ± 432.6216.59 ± 1.55331.94 ± 20.632.29 ± 0.15
266.15 ± 0.091891.80 ± 387.2320.57 ± 2.06376.71 ± 15.252.53 ± 0.14
276.47 ± 0.10996.00 ± 193.8123.25 ± 1.34435.20 ± 127.423.20 ± 0.22
Table 3. Mean average values and standard deviations of measured responses for EPC, SPE, SPP, area to nominal, and volume to nominal.
Table 3. Mean average values and standard deviations of measured responses for EPC, SPE, SPP, area to nominal, and volume to nominal.
RunEPC (MJ)SPE (MJ/g)SPP (kW/g)Area 2 Nom [%]Volume 2 Nom [%]
11.627 ± 0.3090.238 ± 0.0450.035 ± 0.01271.32 ± 0.7468.18 ± 0.68
21.822 ± 0.3780.220 ± 0.0490.052 ± 0.01296.98 ± 0.4796.45 ± 0.55
31.217 ± 0.3370.182 ± 0.0510.056 ± 0.019101.09 ± 0.18100.35 ± 0.29
41.008 ± 0.2650.121 ± 0.0330.035 ± 0.01399.57 ± 0.2399.48 ± 0.39
50.929 ± 0.2460.110 ± 0.0270.044 ± 0.013102.18 ± 0.49100.43 ± 0.64
60.828 ± 0.2440.108 ± 0.0330.062 ± 0.03199.45 ± 0.2997.80 ± 0.25
70.857 ± 0.2470.090 ± 0.0260.027 ± 0.010103.24 ± 0.46101.01 ± 0.64
80.835 ± 0.2410.088 ± 0.0260.047 ± 0.01798.37 ± 0.7596.71 ± 0.77
90.540 ± 0.1670.057 ± 0.0180.040 ± 0.013104.14 ± 0.67103.21 ± 0.62
101.692 ± 0.4600.177 ± 0.0470.027 ± 0.01197.65 ± 0.4497.34 ± 0.55
111.116 ± 0.2230.118 ± 0.0230.027 ± 0.00599.64 ± 0.4998.99 ± 0.53
121.584 ± 0.3450.180 ± 0.0390.053 ± 0.01298.33 ± 0.3698.20 ± 0.42
130.965 ± 0.1790.166 ± 0.0310.037 ± 0.00296.78 ± 0.5596.67 ± 0.58
140.547 ± 0.1580.083 ± 0.0240.036 ± 0.01297.90 ± 1.2298.21 ± 1.27
150.756 ± 0.1820.115 ± 0.0270.057 ± 0.019100.54 ± 1.27100.36 ± 1.27
164.003 ± 0.3660.494 ± 0.0700.052 ± 0.00799.19 ± 1.2398.59 ± 1.21
171.842 ± 0.5800.222 ± 0.0720.034 ± 0.01498.21 ± 0.6798.31 ± 0.68
181.728 ± 0.3980.210 ± 0.0480.038 ± 0.015101.52 ± 0.90101.16 ± 0.82
191.332 ± 0.3950.180 ± 0.0530.069 ± 0.02999.42 ± 0.2399.24 ± 0.32
200.612 ± 0.1820.085 ± 0.0250.074 ± 0.02898.52 ± 0.6297.94 ± 0.53
210.540 ± 0.1630.081 ± 0.0260.044 ± 0.01891.13 ± 0.3490.63 ± 0.35
224.997 ± 0.8440.663 ± 0.1120.086 ± 0.01498.01 ± 0.5197.49 ± 0.49
235.136 ± 0.3120.535 ± 0.0560.068 ± 0.019103.09 ± 0.61103.16 ± 0.59
243.353 ± 0.8740.383 ± 0.0980.101 ± 0.03297.73 ± 0.6297.48 ± 0.71
250.324 ± 0.0990.049 ± 0.0160.022 ± 0.00997.54 ± 0.2097.10 ± 0.18
260.828 ± 0.2230.134 ± 0.0350.075 ± 0.03397.22 ± 0.7296.84 ± 0.79
270.331 ± 0.0640.051 ± 0.0100.054 ± 0.018100.79 ± 0.49100.58 ± 0.48
Table 4. Interaction plot findings.
Table 4. Interaction plot findings.
MetricsCompressive Strength (MPa)EPC (MJ)
Control ParametersSynergisticAntagonisticSynergisticAntagonistic
ORAPS, NT, and BTRDA, LT, and IDPS, NT, and BTRDA, LT, and ID
RDAPS, NT, and BTLT, ID, and ORAPS, NT, and BTLT, ID, and ORA
LTPS, NT, and BTID, RDA, and ORAPS, NT, and BTID, RDA, and ORA
IDPS, NT, and BTLT, RDA, and ORAPS, NT, LT, and BTRDA and ORA
PSID, LT, RDA, and ORANT and BTID, LT, RDA, and ORANT and BT
NTID, LT, RDA, and ORAPS and BTID, LT, RDA, and ORAPS and BT
Table 5. Control parameters for the confirmation runs.
Table 5. Control parameters for the confirmation runs.
RunORARDALTIDPSNTBT
28020.90.110020300106.3
290900.26606030080
Table 6. Mean average values and standard deviations of measured responses for weight, printing time, compressive strength, compressive modulus of elasticity, and compressive toughness for the confirmation runs.
Table 6. Mean average values and standard deviations of measured responses for weight, printing time, compressive strength, compressive modulus of elasticity, and compressive toughness for the confirmation runs.
RunWeight (g)Printing Time (s)sB [MPa]E [MPa]Toughness [MJ/m3]
289.52 ± 0.1110,854.60 ± 497.5080.72 ± 2.041369.03 ± 46.0710.23 ± 0.60
297.27 ± 0.241158.00 ± 144.1942.46 ± 1.061068.52 ± 41.694.73 ± 0.10
Table 7. Mean average values and standard deviations of measured responses for EPC, SPE, SPP, area to nominal, and volume to nominal for the confirmation runs.
Table 7. Mean average values and standard deviations of measured responses for EPC, SPE, SPP, area to nominal, and volume to nominal for the confirmation runs.
RunEPC (MJ)SPE (MJ/g)SPP (kW/g)Area 2 Nom [%]Volume 2 Nom [%]
284.032 ± 0.4050.423 ± 0.0380.039 ± 0.00291.55 ± 0.5691.09 ± 0.70
290.518 ± 0.0750.071 ± 0.0100.062 ± 0.010106.99 ± 0.5187.90 ± 0.49
Table 8. Validation table.
Table 8. Validation table.
Run 2829
ActualsB (MPa)80.7242.46
EPC (MJ)4.030.52
PredictedsB (MPa)85.2634.80
EPC (MJ)3.69Vague
Absolute ErrorsB (%)5.6218.03
EPC (%)8.52Vague
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vidakis, N.; Petousis, M.; David, C.N.; Sagris, D.; Mountakis, N.; Karapidakis, E. Mechanical Performance over Energy Expenditure in MEX 3D Printing of Polycarbonate: A Multiparametric Optimization with the Aid of Robust Experimental Design. J. Manuf. Mater. Process. 2023, 7, 38. https://doi.org/10.3390/jmmp7010038

AMA Style

Vidakis N, Petousis M, David CN, Sagris D, Mountakis N, Karapidakis E. Mechanical Performance over Energy Expenditure in MEX 3D Printing of Polycarbonate: A Multiparametric Optimization with the Aid of Robust Experimental Design. Journal of Manufacturing and Materials Processing. 2023; 7(1):38. https://doi.org/10.3390/jmmp7010038

Chicago/Turabian Style

Vidakis, Nectarios, Markos Petousis, Constantine N. David, Dimitrios Sagris, Nikolaos Mountakis, and Emmanuel Karapidakis. 2023. "Mechanical Performance over Energy Expenditure in MEX 3D Printing of Polycarbonate: A Multiparametric Optimization with the Aid of Robust Experimental Design" Journal of Manufacturing and Materials Processing 7, no. 1: 38. https://doi.org/10.3390/jmmp7010038

APA Style

Vidakis, N., Petousis, M., David, C. N., Sagris, D., Mountakis, N., & Karapidakis, E. (2023). Mechanical Performance over Energy Expenditure in MEX 3D Printing of Polycarbonate: A Multiparametric Optimization with the Aid of Robust Experimental Design. Journal of Manufacturing and Materials Processing, 7(1), 38. https://doi.org/10.3390/jmmp7010038

Article Metrics

Back to TopTop