Abstract
In this research, a novel fusion model which incorporates a type-1 fuzzy logic system (T1FLS) and the Dempster-Shafer theory (DST) is proposed to model the 3D fused deposition process. First, the experimental data, which are constructed based on the Taguchi L18 array for two filaments (i.e., Polyether-ether-ketone and Polyether-ketone-ketone), are used to develop T1FLSs that have different structures. Second, a fusion algorithm integrating fuzzy logic and the DST is presented to integrate the predicted values of the T1FLSs after analysing the behaviours of such systems in the space examined. The proposed model can integrate such systems in a way that can resolve possible conflicts among the models developed in the first stage and, as a result, it can improve the predictive performance. Validated on a set of experimental data, the proposed fusion model has improved the predictive performance with an average improvement of 25.6%. The fusion model developed is employed to anticipate the mechanical characteristics of a dental part produced by fused deposition.
Introduction
Additive manufacturing, or the so-called 3D printing and rapid manufacturing, is well established as an innovative manufacturing approach for complex manufacturing industries (Conner et al., 2014; Huang et al., 2013; Levy et al., 2003). In brief, such an approach involves the generation of complex 3D parts from pre-designed 3D structures by adopting material additive processes including, for instance, layer-by-layer or surface-by-surface ones (Patil et al., 2021). In addition to its considerable effect on innovation and the manufacturing and its related industries, additive manufacturing is considered to be a vital constituent of the Industry 4.0 revolution because of its (i) characteristics that facilitate low-volume but cost-effective manufacturing; (ii) ability to generate highly complex monolithic structures; and (iii) capability of handling changes and, thus, allowing flexible production of customized 3D printed parts using different colours and materials at small as well as large scales (Cano-Vicent et al., 2021; Doshi et al., 2021; Ivanova et al., 2014; Murr, 2015). Consequently, many research studies have been directed to investigate the various 3D printing techniques across a vast myriad of industries including, but not limited to, aerospace, pharmaceutics and tissue engineering (Alizadeh-Osgouei et al., 2021; Boesch et al., 2019; Parulski et al., 2021; Short, 2015). Various state-of-the-art techniques have been, in general, covered under the umbrella of additive manufacturing. Such techniques include fused deposition modelling (FDM), resin printing (i.e., stereolithography (SLA)) and selective laser sintering (SLS) (Crump, 1992; Deckard, 1989; Hull, 1986; Sachs et al., 1993). FDM, whose essence is the production of 3D printed parts by controllably depositing liquefied thermoplastic polymer layers, is considered to be one of the most common techniques in various applications (e.g., biomedical science) (Singh, Prakash et al., 2019). This is because of the ability of the FDM process to deal with different biomedical polymers (e.g., Polyether-ether-ketone (PEEK) and Polyether-ketone-ketone (PEKK)) (Singh, Babbar et al., 2019). Thus, many studies have concentrated on the FDM process and its various applications.
The various parts and polymers 3D printed using FDM for various applications, in particular, biomedical ones were examined in order to replace the aggravated autograft and allograft ones (Singh, Prakash et al., 2019). Polylactic acid (PLA) scaffolds, for example, were 3D printed by employing the FDM technique under different sets of FDM parameters, in order to examine their properties (Singh, Babbar et al., 2019). Furthermore and to enhance its quality attributes and imitate the human bone structure, PLA was amalgamated with polyvinyl alcohol, thermal-stimulus-based hydroxyapatite and nanohydroxyapatite and, then, used to produce scaffolds. The properties of such scaffolds (e.g., mechanical and rheological properties) were analysed and then compared to the PLA ones (Alizadeh-Osgouei et al., 2021; Esposito Corcione et al., 2017; Prakash et al., 2021; Singh et al., 2020; Song et al., 2018). In addition, various orthopaedic and dental 3D printed parts were produced using the FDM process by PEEK and carbon-fibre-reinforced PEEK, and their mechanical and biocompatibility characteristics were assessed (Dawood et al., 2015). Since various parameters can play a significant role in defining the fate of the final 3D printed parts, several research papers investigated these parameters. These parameters can be, in general, categorised into: (i) material-related parameters such as the mechanical and rheological characteristics of the polymers; (ii) process-related parameters such as humidity; and (iii) machine-related parameters such as print speed (Nasereddin et al., 2018; Parulski et al., 2021). To design and print 3D printed parts having the required quality attributes, such parameters need to be controlled and optimized. Therefore, a spectrum of experimental techniques has been proposed to control and optimize these parameters (Sheoran & Kumar, 2020). An analytical paradigm was, for instance, developed to assess the combinatory influences of the material- and machine-related parameters on the attributes of 3D printed specimens (Khan et al., 2021). Likewise, the effects of various parameters such as orientation, print speed and infill type on the attributes of 3D printed PLA parts were statistically analysed by determining the correlation coefficient and the P-value using the analysis of variance test (Ansari & Kamil, 2021; Carlier et al., 2019; Chacón et al., 2017; Ouhsti et al., 2018; Pandzic et al., 2019). Moreover, the impact of copper electroless plating on 3D printed specimens produced using the acrylonitrile butadiene styrene (ABS) polymer was also studied (Karthikeyan et al., 2025). Furthermore, digital imaging was utilized to show how the mechanical attributes of the 3D printed parts were affected by different thickness values (Xu et al., 2021). Likewise, the influences of four FDM parameters, namely, layer thickness, raster angle, infill density, and nozzle temperature on the mechanical properties were examined for PLA, ABS and Chlorinated Polyethylene (CPE) polymers (Alsa’di et al., 2024). In addition, the influence of the infill pattern on fatigue life of 3D printed parts produced using ABS was studied by examining various infill pattern geometries (Dudescu et al., 2023). It was proved by comparative evaluation of the determined fatigue life that infill pattern geometries had considerable effects on fatigue life (Dudescu et al., 2023).
Nowadays, there is a strong need to anticipate the quality characteristics of 3D printed specimens accurately, in particular, in biomedical science, pharmaceutics and tissue engineering. Therefore, several studies have been devoted to implementing and developing systems-engineering paradigms that can simulate the FDM process. For instance, various multi-criteria decision-making paradigms (e.g., VIKOR) were employed to determine the best combination of the FDM parameters (Deomore & Raykar, 2021; Raykar & D'Addona, 2020). In addition, the genetic algorithm was embedded in the structure of the artificial neural network to anticipate the strength of 3D printed specimens prepared using various polymers (Giri et al., 2021; Pazhamannil et al., 2021; Teharia et al., 2021; Yadav et al., 2020). Furthermore, a modelling structure based on conceptual modelling and neural networks was proposed to better understand and, consequently, model the FDM process. For instance, an artificial neural network was designed to represent the mechanical attributes of three polymers (i.e., PLA, ABS and CPE) as a function of four FDM parameters (Alsa’di et al., 2024). In addition, an interval type-2 fuzzy logic system (IT2FLS) was designed to mimic the FDM process and to handle the uncertainties in the measurements of the quality attributes of the 3D printed parts (AlAlaween, Alalawin et al., 2022). Such a system had the ability to predict the quality attributes of the 3D printed parts and provide a linguistic understanding of the relationships between the FDM parameters and the quality attributes. However, the performance of such a system was represented by the average predicted performance for all the areas in the space examined. In other words, the system behaviours at different areas in the spaces examined cannot be shown. In addition, a modelling paradigm based on basis functions was designed to mimic the FDM process and predict the attributes of the parts produced (AlAlaween, Abueed et al., 2023). Such a paradigm consisted of a number of radial basis functions that play an integral role in extracting the relationships between the FDM parameters and the attributes examined. However, such a paradigm cannot handle uncertainties and cannot provide users with linguistic representation of the relationships extracted. Furthermore, some research papers have focused on optimizing the FDM process and its different parameters (AlAlaween, AlAlawin et al., 2023; Lalegani Dezaki et al., 2021). For instance, the process was optimized by defining the different challenges faced while printing polymers, composites, geopolymers and novel materials, where machine parameters were investigated to optimize printing novel polymers for various applications (Lalegani Dezaki et al., 2021). Likewise, a right-first-time concept-based paradigm that integrated fuzzy logic and multi-objective swarm optimization was proposed to define the best FDM parameters that need to be employed to 3D print parts with predefined mechanical attributes with minimum waste and recycling ratios (AlAlaween, AlAlawin et al., 2023).
Although a huge body of research has concentrated on understanding and modelling the FDM process, and mapping the quality attributes of the 3D printed parts to the FDM parameters, there is an urge to develop a paradigm that can (i) accurately predict the quality attributes for a highly dimensional space resulted from the considerable number of the FDM parameters that affect these attributes; (ii) deal with a limited number of sparse data points because of the cost and time required to produce more 3D printed test specimens; and (iii) handle measurement uncertainties. Ascertaining all these features may not be possible with a single paradigm. Therefore, in this paper, a novel fusion model that incorporates the Dempster-Shafer theory (DST) and a type-1 fuzzy logic system (T1FLS) is presented to simulate the FDM process and study the effects of the process parameters on the mechanical characteristics of the parts produced using PEEK and PEKK. The fusion model is developed to (i) integrate various T1FLSs that have different structures developed to represent the FDM process, such systems can play an integral role in extracting all possible relationships between the process parameters and the mechanical characteristics; (ii) examine the behaviours of these systems in the space investigated instead of representing the predictive performance by its average values in the different areas of the space; (iii) resolve possible conflicts among the various T1FLSs; (iv) improve the modelling predictive performance in the areas with unacceptable performance by considering the systems that provide the best predictive performance; and (v) tackle measurement uncertainties by employing T1FLS and the DST. The paper is structured as follows. Section 2 briefly summarizes the equipment used and the experimental work. The development of the proposed fusion model and the mathematics behind the DST and the T1FLS are discussed in Section 3, whereas its results are summarized in Section 4. Conclusions are, then, listed in Section 5.
Experimental Work
Two types of polymers (i.e., engineering-grade PEEK and PEKK) filaments were studied in this paper. PEEK and PEKK are both semicrystalline thermoplastic high-temperature polymers that have relatively excellent mechanical as well as chemical resistance characteristics. Such polymers were provided by 3DXTECH (Michigan, USA). ASTM-D638 parts were produced using the FUNMAT HT 3D printer, as a functional-materials FDM printer, as shown in Figure 1, (INTAMSYS Technology Inc., Minneapolis, USA). Such a printer is supported by InstamSuite 3.6.2 program that was employed to convert the 3D designed structure of the ASTM-D638 parts to a GCODE format. The ASTM-D638 parts were printed using various sets of operating conditions. In addition to the three mechanical characteristics studied in this paper, the eight FDM parameters with the identified levels, which were studied for PEEK and PEKK, are listed in Table 1.

The FUNMAT HT 3D Printer.
The Fused Deposition Modelling Parameters and the Mechanical Characteristics (AlAlaween et al., 2022).
The Taguchi L18 array was utilized to understand the relationships and to model the FDM process. Therefore, 18 experiments were carried-out for each polymer. To ensure good repeatability, each experiment was repeated three times. Once they were printed, the specimens were removed from the printer glass plate. According to manufacturer recommendations, all specimens produced were annealed. The mechanical attributes were experimentally estimated. The micro-hardness was measured using the micro Vickers hardness (HTMV 2000M, echo LAB, Italy) with load and time values of 0.98 N and 15 s, respectively. Whereas the ultimate tensile strength and elongation were determined using Instron (SHIMADZU, USA) with load and speed values of 50 kN and 1 mm/min, respectively. The parameters of these apparatuses are summarized in Table 1. For each experiment, the average of the three repetitions was then estimated.
For the various sets of the FDM parameters, various fracture patterns occurred at different locations. It was apparent that the investigated parameters have various influences on the mechanical attributes examined in terms of the nature of the relationships and their strength, which were proved by calculating the correlation values that are listed in Table 2. Such coefficients are reasonable. In addition, some of the FDM parameters have different correlation values for PEEK and PEKK. To elucidate, the layer thickness has a negligible impact on the strength for PEEK, however it has a significant one on that for PEKK. It also has a considerable effect on the elongation for PEKK but a negligible one for PEEK. It is also worth noting that the relationships between the layer thickness and the strength for PEEK and PEKK are direct and indirect, respectively. Likewise, the relationships between the layer thickness and the elongation for PEEK and PEKK are direct and indirect, respectively. Moreover, different natures of the relationships can be seen between the FDM parameters and the mechanical attributes examined. For instance, the relationships between the print speed and specimens’ strength for PEEK and PEKK are indirect and direct, respectively. In addition, the layer thickness had direct and indirect effects on micro-hardness for PEEK and PEKK, respectively. It is also apparent that the layer orientation can considerably affect the strength and elongation. However, it has an inconsiderable impact on micro-hardness for both materials. For the infill pattern, the analysis of variance was used to study its influences. It was found that it significantly affected the mechanical attributes examined with P-values of less than 0.05. It was found that the P-values for the ultimate tensile strength, elongation and microhardness are 0.032, 0.030 and 0.027.
The Correlation Coefficients.
The advances in the computing power that have been witnessed recently make data-driven models an ideal way to represent and model complex processes, in particular, when physical models are difficult to derive and/or implement (AlAlaween, AlAlawin, Al-Durgham et al., 2021). This has been the main reason behind the extensive use of various data-driven models such as linear regression, artificial neural networks and fuzzy logic in various areas including, for example, pharmaceutics, medicine, automobile, supply chain and manufacturing (Alalaween, Alalawin, Mahfouf et al., 2021; AlAlaween, Faouri et al., 2022; AlAlawin et al., 2022). The core of these models, as the name indicates, rests with data in terms of the amount of data available as well as the distribution in the space investigated (AlAlaween et al., 2020; Alshafiee et al., 2019). To elucidate, having a good number of data points is of paramount importance to extract and simulate the inputs/outputs relationships of a process. In addition, a good distribution of the data in the space examined can prevent developing a biased paradigm that can perform satisfactorily in those areas where the amount of data is enough and unsatisfactorily elsewhere (AlAlaween, Mahfouf et al., 2021). Therefore, integrating several models that may have different structures can circumvent these challenges and help in extracting relationships that may not be extractable using a single model. However, integrating different models is not as simple as it may sound, as one needs to consider the behaviours of the models in the space examined and, thus, integrate them in a way that can improve the predictive performance (AlAlaween et al., 2017). Various paradigms can be employed to integrate different models. Information fusion, as a concept that simulates the human cognitive process, is considered to be one of the effective ways of integrating information from different sources in a way that improves the reliability of such information and, thus, allows decision makers to take optimal decisions (Frikha & Moalla, 2015).
Among the various algorithms presented in the related literature, the DST is considered to be an efficient one, this being due to its ability to consider imprecision and possible conflicts among various sources of information (AlAlaween et al., 2017). However, the DST deals only with uncertainties due to probabilities and due to lack of specifications. Therefore, and in order to tackle the uncertainties due to fuzziness, an information fusion algorithm based on both fuzzy logic and the DST is presented in this research, and incorporated with the T1FLS to (i) integrate various T1FLSs that have different structures; (ii) examine the different behaviours of the T1FLSs developed to simulate the FDM process in the space investigated; (iii) resolve possible conflicts among the various T1FLSs; and (iv) improve the modelling predictive performance, in particular when the data available are limited and/or sparse.
Figure 2 depicts the schematic structure of the proposed fuzzy fusion model that incorporates the DST and the T1FLS. The proposed structure consists of several stages. First, the experimental data are used to develop M T1FLSs that include different architectures (i.e., the number of fuzzy sets and fuzzy parameters). Such T1FLSs can capture the process behaviours represented by the data and act synergetically when simulating and extracting the possible behaviours of the process in the space examined. Second, based on the predicted outputs obtained from the T1FLSs developed, the fusion algorithm is implemented to integrate the predicted outputs of the T1FLSs. The behaviours of the T1FLSs that have different structures are evaluated in the space investigated by clustering the performance of these models into several clusters described linguistically (e.g., Unsatisfactory, Satisfactory and Excellent). For each predicted output, the membership degrees are then estimated with respect to these clusters. Then, the mass functions are assigned. Such a step is followed by combining the hypothesis which are the clusters indicating the performance of the T1FLSs. In this structure, a hypothesis of one T1FLS should be combined with another hypothesis that has an equal or better predictive performance value. In other words, a high level of conflict needs to be assumed between a hypothesis and a one with low predictive performance. This can improve the predictive performance. The membership degrees with respect to the combined hypotheses are then calculated and defuzzified to determine the final predicted outputs.

The Schematic Structure of the Fuzzy Fusion Model.
The mathematics behind the proposed fuzzy fusion model can be presented by including the mathematics behind the T1FLS and the DST based on fuzzy logic. Although the mathematics behind them have been publicised (AlAlaween et al., 2017; Bishop & Nasrabadi, 2006; Boudraa et al., 2004; Frikha & Moalla, 2015), in this research, the key developments are briefly summarized.
Nowadays, fuzzy systems have been used to simulate and represent complicated processes using data provided and/or expert knowledge, this being due to their ability to (i) represent the relationships among inputs and outputs; (ii) handle process uncertainties by using fuzzy sets; and (iii) expand the linguistic representation of a process by extracting fuzzy If/Then rules that can be employed to control such a process. Therefore, they have been utilized in many applications such as supply chain and pharmaceutics (Alalaween, Alalawin, Mahfouf et al., 2021). Among the various fuzzy systems in the related literature, the T1FLS is utilized in this research paper because of its simplicity.
Figure 3 depicts the representation of the T1FLS, which contains four steps: fuzzification, inference, extracting rules and defuzzification. First, the process inputs, which are commonly in the crisp/singleton form (x1, x2 … xn), are fuzzified to identify the fuzzy inputs (i.e.,

The T1FLS Structure.
Application-wise, a single output is needed, therefore, the output fuzzy set is defuzzified via the many techniques already available. In this research article, the centroid defuzzifier was utilized.
Information fusion is considered to be one of the main cognitive processes employed by human to integrate information from different sources. Such a process aims to realise efficient inferences that lead to making the best decision (Boudraa et al., 2004). The need for such a process stems from the fact that (i) information obtained from a single source may be limited; and (ii) information accuracy may not be as required (Maseleno et al., 2015). It has been utilized in different areas such as medical, marine and manufacturing to enhance the reliability and accuracy of information provided (Maseleno et al., 2015).
Various information fusion algorithms have been proposed in the literature. Such algorithms include Bayesian, fuzzy logic and the DST (Maseleno et al., 2015). Among the various algorithms presented, the DST is considered to be an efficient one, this being due to its ability to take into consideration imprecision and resolve possible conflicts among various sources of information (AlAlaween et al., 2017). However, the DST can only consider uncertainties due to probabilities and lack of specification (AlAlaween et al., 2017). However, in order to develop a reliable and effective fusion paradigm, three types of uncertainties due to probabilities, lack of specification and fuzziness need to be carefully considered (AlAlaween et al., 2017). However, the DST can only deal with the first two types of uncertainties. Therefore, in this research article, a paradigm that integrates the DST and fuzzy logic is presented, where the third type of uncertainties can be dealt with using fuzzy logic.
Dempster-Shafer Theory Based on Fuzzy Logic
In general, the DST, or the so-called the evidence theory, is a commonly used framework in the cognitive process for reasoning with uncertainties. Such a theory can only handle uncertainties due to probabilities and lack of specifications. In addition, determining the mass function of the hypothesis to be combined is considered to be a challenge in performing the DST (AlAlaween et al., 2017; Boudraa et al., 2004). It can be estimated using various approaches including distance or probabilities. In order to circumvent such a challenge and consider uncertainties due to fuzziness, the DST is integrated with fuzzy logic in this research to successfully identify the mass function.
The hypotheses are identified based on the performance of the predictive models in the space investigated. The performance of the developed models is measured using the error residuals. In this research paper, unsupervised clustering is utilized to identify such hypotheses by classifying the predictive performance values into three classes: Unsatisfactory, Satisfactory and Excellent. The input parameters and the error values for each data point from all the models developed are then classified. The membership degree for each point is calculated with respect to the three clusters as follows:
Fusion Model: Implementation and Results
In this research work and as depicted in Figure 2, 10 T1FLSs that have different structures, as previously stated, were developed for each output (i.e., a mechanical characteristic). Each T1FLS was established by partitioning the data into training (28) and testing (8) sets. The training and testing sets are generally utilized to extract the relationships and evaluate the T1FLS's generalization capabilities, respectively (Alalaween, Alalawin, Mahfouf et al., 2021). The training and testing sets of the 10 T1FLSs have the same numbers of data points. However, such points have different distributions in the space examined. Therefore, the various T1FLSs developed can simulate the input/output relationships. To develop a representative system for the FDM process, the FDM parameters needs to be carefully considered. For instance, material type and infill pattern were considered as singleton variables (i.e., crisp or discrete), whereas the remaining FDM parameters were considered as continuous variables. Various numbers of rules were examined for each T1FLS. The optimal one that was chosen was the one that provided the optimal predictive performance determined by the root mean square error (RMSE) (i.e., the minimum RMSE value). For a specific number of rules, the system's parameters were initially assigned via the hierarchical approach presented in (Zhang & Mahfouf, 2011). These parameters were then optimized during the training process by utilizing the steepest descent procedure that was embedded in the back-propagation paradigm (Mendel, 2001).

One of the 11 Rules used for the Ultimate Tensile Strength Rule (Materials (1) PEEK and (2) PEKK; Infill Pattern (1) Cubic and (2) Grid).

The Performance of a T1FLS for the Ultimate Tensile Strength with a 90% Confidence Level: (a) Training and (b) Testing Sets.
For the ultimate tensile strength, 11 rules, as the best number of rules, were used to develop the T1FLS. Figure 4 shows an example of one of these rules. The linguistic representation of such a rule is presented in Figure 4.

The Three Clusters of the Predictive Performance of the T1FLSs.

The Performance of a T1FLS in the Space Area of the Ultimate tensile Strength.

The Predictive Performance of the Proposed Fused Model for the Strength with a 90% Confidence band: (a) Training and (b) Testing Sets.
Figure 5 shows the predictive performance for one of the T1FLSs developed using 11 rules, as the best number of rules for this system. The RMSE values for training and testing sets are 9.65 MPa and 8.97 MPa, respectively, whereas the R2 values for the corresponding sets are 0.81 and 0.78, respectively. Such performance measures represent the average values of all points, in other words, they do not show how the model behaves in the space examined. To elucidate, it is apparent from Figure 5 that 14 data points in the training set out of 28 and 4 data points in the testing set out of 8 lay outside the 90% confidence interval. Furthermore, the data points in both sets are scattered around the best-fit line. Such behaviours were not detected by looking at both the RMSE and R2 values. Such behaviours can be attributed to (i) the highly nonlinear relationships between the FDM parameters and the ultimate tensile strength; (ii) the limited number of data points that are used to develop a T1FLS, as a data-driven model; and (iii) the number of input parameters (i.e., nine) examined which can result in the phenomenon of the “curse of dimensionality” when the data points are limited. Therefore, and based on the results of the 10 T1FLSs developed that have similar behaviours to the system presented in Figure 5 but different in the space examined, it was obvious that a single T1FLS cannot be used to simulate the FDM process and anticipate the mechanical characteristics using such limited data points. Therefore, there is a need for the new fusion model presented in Section 3.
Once the 10 T1FLSs were developed, their performance behaviours in the space were classified into the three abovementioned classes using the k-means clustering algorithm. Therefore, the FDM parameters and the residuals of these 10 T1FLSs were employed to identify the behaviours of the models (i.e., areas where the performance of a model is Unsatisfactory, Satisfactory and Excellent). The parameters of the identified clusters were estimated. Figure 6 depicts the membership functions for these three clusters with their parameters for the 10 T1FLSs developed.
Figure 7 shows how one of the developed T1FLSs behaved in an area for the strength. It indicates that the performance of such a T1FLS is Excellent when the material is PEEK, infill pattern is grid, layer thickness is high, infill density is high, print speed is high, number of shells is medium, cooling rate is small, layer orientation is 45o and raster width is medium.
Based on the identified clusters for the 10 T1FLSs, the defined membership was utilized to determine the mass functions, as described in Section 3. Such mass functions were then integrated as mathematically shown in Equation (4) in order to find the ones for the fusion model that integrated the 10 T1FLSs. Then, the membership function values of such a model were estimated by numerically solving Equation (3). Finally, the outputs were then determined by the height defuzzification. Figure 8 depicts the performance of the proposed fusion model for the strength, where the RMSE (training, testing) and R2 (training, testing) are [3.57, 6.99] MPa and [0.96, 0.96], respectively. The testing RMSE value is twice the training one. However, such a high value is due to that 6 data points that have values of 60 MPa or greater, whereas the majority of the data points in the training set have values that are less than 60 MPa. Therefore, such a high value is not due to overfitting. This can be proven by looking at the R2 values. It is apparent that the performance of such a newly model outperformed the single system with an average overall improvement of 20% in R2 for the ultimate tensile strength. It is also noticeable that all the points are within the confidence band. Furthermore, they are scattered neatly around the best-fit line.
In a similar manner, a T1FLS and the fusion model based on 10 T1FLSs were utilized to predict the elongation values and micro-hardness. The predictive performance for these models is listed in Table 3. Such a table shows that the predictive performance values of T1FLSs developed for both elongation and micro-hardness were not acceptable. Therefore, the fusion model was used to improve the predictive performance. The performance values of the fusion models developed for the two mechanical characteristics are better than those of the T1FLSs with average overall improvements of 27% and 30% in R2, respectively.
The Performance Values.
For comparison purposes, IT2FLS and the radial based integrated network presented in (AlAlaween, Alalawin et al., 2022) and (AlAlaween, Abueed et al., 2023), respectively, were employed to predict the three mechanical characteristics. The results obtained are also summarized in Table 3. It is noticeable that the fusion model outperformed these models with a considerable improvement. This can be attributed to (i) its structure that consists of several T1FLSs that can play an integral role in extracting and representing the relationships between the FDM parameters and the mechanical characteristics; (ii) its ability to analyse the behaviours of these systems in the different areas of the space examined; and (iii) its ability to integrate the predicted outputs in these areas in a way that can improve the predictive performance. In addition, the fusion model can provide a linguistic understanding of the relationships between the FDM parameters and the mechanical characteristics that can be utilized to control and understand the process.
The capabilities of the fusion paradigm in anticipating the three mechanical characteristics of the 3D-printed parts were demonstrated by predicting these attributes for a dental implant (i.e., MidFace Rim) presented in Figure 9. Such a part was produced by the same printer mentioned in the second section using PEEK. The printing parameters in terms of infill pattern and its density, thickness, printing speed, number of shells, cooling rate, printing orientation and raster width were cubic, 100%, 0.12 mm, 20 mm/s, 3, 0%, 0° and 0.4 mm, respectively. It is worth noting that measuring the strength and elongation for such a part, which is not printed according to the specifications of the standard test parts, is not as simple as it may seem and may not provide accurate measurements. Therefore, the new fusion model can be employed in this case and perhaps similar ones to anticipate these attributes. Therefore, the mechanical characteristics were predicted using the fusion models developed in Section 4.1. The anticipated values of the MidFace rim were 102.4 MPa, 19.5 and 29.3 for the strength, elongation and micro-hardness respectively. The average of measured values of the micro-hardness for the MidFace Rim is 30.5 which is close to the anticipated one (i.e., 29.3).

Midface Rim.
A novel fusion model based on incorporating a type-1 fuzzy logic system (T1FLS) and the Dempster-Shafer theory (DST) was presented in this research article to simulate the 3D fused deposition process. Such a model was developed in two stages. First, several T1FLSs having different structures were developed. Second, the fuzzy based DST was employed to integrate the predicted outputs of the T1FLSs after analysing their behaviours in the space examined. In addition to integrating and analysing the behaviours of the T1FLSs, such a model can resolve possible conflicts among these systems and, as a result, can lead to a better predictive performance. In summary, the fusion model presents a promising advancement not only in 3D printing and its techniques but also in other equally complex processes with highly dimensional spaces and a limited number of sparsely distributed data points. In the future, the fusion model can be incorporated with multi-objective optimization paradigms in order to identify the best 3D printing parameters that need to be employed to produce 3D printed specimens with predefined attributes. This incorporation can lead to minimizing waste and recycling ratio and reducing time-to-market.
Footnotes
Acknowledgment
The authors would like to thank The Royal Academy of Engineering (UK) and Industrial Scientific Research and Development Fund-The Higher Council for Science and Technology (Jordan) for the financial support. The authors would like to thank Abdul Hameed Shoman Foundation for the financial support.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by The Royal Academy of Engineering (UK) and Industrial Scientific Research and Development Fund-The Higher Council for Science and Technology (Jordan) for the financial support (Grant number TSP1134). This work was also supported by Abdul Hameed Shoman Foundation (Grant number 230800045).
Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
