Abstract
A multitude of rapid prototyping (RP) systems and technologies have come up since the introduction of additive process. Owing to the enlarging number of these systems with distinctive efficacy, the problem of selecting an appropriate system for a particular requirement is a cumbersome task. Henceforth, this work comes up with a strategy based on multi-attribute decision making to select a most suitable RP system. The presence of subjectivity in decision making as well as the existence of imprecision from various sources emphasize the methods which must consider uncertainty and vagueness. A decision advisor based on uncertainty theories, including fuzzy analytical hierarchy process (FAHP) and grey relational analysis (GRA) has been introduced. It provides a comprehensive database comprising thirty nine commercially available RP systems. The evaluation attributes consisting of machine cost, accuracy, layer thickness, machine speed, material cost, net build size volume, machine weight, surface roughness, and material strength were utilized to characterize the different machines. The FAHP based on trapezoidal fuzzy number was implemented to determine the priority weights of various attributes, while the GRA was employed to realize the best RP system and technology. The authors believe that this system has the potential to transform into a fully developed RP selection system.
Keywords
Introduction
The rapid prototyping (RP) (also known as Additive Manufacturing, AM) due to its innumerable benefits has revolutionized the manufacturing sector. Since the inception of Stereolithography (SLA) the first commercial AM printer in 1987, different AM technologies have emerged and evolved over the years. The AM is no longer used only for visualization purpose or prototyping, but are also being used to fabricate fully functional products. Some of the primary benefits of AM include reduced lead time (most often by 75%), reduced cost (many times, the cost can be minimized by 50%), fabrication of highly complex shapes as well as high design flexibility [1]. The AM process exhibit applications in various manufacturing industries such as aerospace, automobiles, home appliances, electrical, defense and medical applications. According to the Wohler’s report of 2018, there were more than 135 AM companies in 2017 which produced and sold industrial AM system worldwide [2]. Moreover, the polymer based AM machines occupies the largest share among the global AM Industry [3]. In 2017, the growth of AM worldwide was $7.36 billion, that was 21% higher when compared to $6.063 billion in 2016 [4].
The AM machines depending on the type of material used, can be categorized into three systems, namely liquid, solid (wire, sheet) and powder. In addition, the AM systems can also be classified into different technologies (processes) such as SLA, Selective Laser Sintering (SLS), Fused Deposition Modeling (FDM), etc. Among them, the SLA belongs to liquid based system, SLS comes under powder and FDM belongs to solid based AM systems. Each of these processes (or technologies) have their own strength, utilities, applications and limitations. Each AM machine have several attributes (machine characteristics) such as machine speed, material cost, layer thickness, accuracy etc., which drive their selection. The speed, cost, weight and size are not the only factors that should be considered while choosing a particular AM technology. The other factors consisting of surface finish, material strength, accuracy, etc., also varies for each process. Certainly, the strength and weakness are inherent in any manufacturing process and AM is no exception. It is crucial to realize that not one AM machine can do everything or can be used in all applications.
Due to the rapid growth of AM technology, the selection of the most appropriate AM system or machine has become increasingly important. A large number of new technologies are being introduced each year. With such a wide range of AM systems, technologies and machine attributes, selecting an AM machine for a user or a manufacturer is a challenging task and it takes a considerable amount of time and require greater expertise. There is no formal evaluation procedure which is generally followed in the selection of AM technique [2]. Indeed, the selection of optimum or most suitable system for AM process depends on many factors such as accuracy, material strength, build speed, surface finish and other system related parameters. To select the most suitable AM machine for a specific application, one must understand the strength and limitations of each AM machine. A number of factors (selection criteria) such as accuracy, speed, layer thickness, surface finish etc., which can be quantitative or qualitative must be synthetically considered while evaluating and selecting an AM machine. With so many factors involved and different AM machines available, it is really difficult for the potential customer or user to differentiate among them. In addition, the potential customer also confronts further complexity in choosing AM machine, as a result of the lack of knowledge of new AM system and its application. Consequently, to administer the selection of AM system efficiently, which is a multi-attribute problem, the decision makers should employ Multi-Attribute Decision Making (MADM) approaches. The MADM method represents a means of formulating fundamental information and providing judgements in situations involving multifarious objectives [5]. In the course of time, numerous approaches including Elimination and Choice Expressing Reality (ELECTRE) [6], Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) [7], Preference Ranking Organization Method For Enrichment of Evaluations (PROMETHEE) [8], Data Envelopment Analysis (DEA) [9], VIKOR [10], Full Consistency Method (FUCOM) [11, 12], Multi-Attributive Border Approximation Area Comparison (MABAC) [13–15], Conjoint Analysis [16], Complex Proportional Assessment (COPRAS) [17, 18], etc., have been introduced. Henceforth, different MADM methods are available, which can be employed for various real decision problems. Certainly, no method is superior to the other because each possesses their own benefits and limitations in order to find an appropriate final solution [19]. It is generally difficult to differentiate between two MADM methods because they are based on different principles or they employ contrasting tactics, when dealing with diverse data sets [20, 21]. The selection of MADM is an extensive and a painstaking process which must take into regard the available information, decision maker’s role and all the decision process abilities, etc. [22]. For a given problem, the MADM approach can be selected either randomly or depending on the user’s knowledge and expertise with the particular technique or the ease of use of that technique [23, 24].
A number of studies were carried out based on MADM situations using different multicriteria decision support systems, benchmark studies and software tools for assisting AM users in selecting the most suitable AM machine. The first attempt in the selection of AM using a computer based system was done by Hornberger (1993) at Santa Clara University for education purpose [25]. In the other, relational databases were introduced by Campbell and Bernie (1996) [26]. Bauer et al introduced a software tool for AM selection that assisted the users in choosing the best combination of material and machines [27]. Philipson and Henderson (1997) from Arizona University introduced the AM advisor for selecting the best AM process based on relational database applying MS Access [28]. All these systems provided general and simple process selection information of AM process. Bibb et al. [29] and Masood and Soo [30] were the first to use the knowledge-based expert system for the selection of AM among availablecommercial AM machines, though this system lacked material, cost and build time information. Moreover, it was difficult to rank the most suitable AM among them using conditional statements. The web based AM selection techniques emerged in 2002. They consisted of a series of questions with a scaling factor for each question thus making it hard for the user in a single analysis [31]. Similarly, another web based AM selector was launched in 2005 with few questions and limitations in the ranking process [32]. Mahesh et al. (2005) [33] developed the first fuzzy decision making and benchmarking based on the five AM process by measuring several parts and building a database. The drawback of the system was that it did not adopt a generic approach and was applicable to specific household application. Byun and Lee [34] presented a modified TOPSIS method for AM selection process, but the fuzzy approach for assigning the machine attributes required relatively more computation. A total of six machine attributes comprising accuracy, surface roughness, part cost, strength, elongation and build time were used in their selection process. The Spanish Rapid Manufacturing Association (ASERM, 2006) launched an AM selector for the use in Additive prototype process and Rapid tooling methods. However the system was inadequate due to the absence of material information and costing [35]. In recent times, due to the prominence of grey theory in MADM, a number of approaches were introduced to assist in the RP selection process. For example, Mahapatra and Panda [36] proposed a benchmarking of AM system using Grey Relational Analysis (GRA) based on six AM machines and six machine characteristics such as dimensional accuracy, surface quality, part cost, material properties and build time. Similarly, Shende et al. [37] employed TOPSIS for ranking four AM processes: SLA, SLS, Three Dimensional Printing (3DP) and FDM based on seven attributes: accuracy, surface finish, tensile strength, elongation, part cost, process time, heat deflection temperature. In addition, Narayan panda et al. (2014) [38] developed an integrated Analytic Hierarchy Process (AHP) and fuzzy TOPSIS method for the selection of five AM processes: Laminated Object Manufacturing (LOM), SLS, 3DP, FDM and SLA based on five attributes: dimensional accuracy, surface quality, part cost, build time and material properties. Furthermore, Nitish et al. (2016) [39] compared two different AM system, namely FDM and an open source low cost AM machine for comparing dimensional performance using GRA method. Ivan peko et al. [40] also proposed an AHP, Fuzzy AHP (FAHP) and PROMETHEE method for AM selection based on four processes: 3DP, FDM, SLS, Photopolymer Jetting (PolyJet) and six attributes including dimensional accuracy, surface roughness, mechanical properties, process cost, process time and post-processing. Lately, different preference relations, including intuitionistic multiplicative preference relations [41], hesitant fuzzy preference relation [42, 43], etc., have also been utilized to derive priority weights in MADM. Although all the research works in literature utilized some logical process for AM selection. However, these works cannot be identified as the most comprehensive and best as far as the selection of AM process is concerned due to the following reasons. Most of the existing methodologies used in the selection of AM were incomplete owing to limited number of attributes and processes. All the methods were aimed at comparing the AM processes rather than selecting the best AM machine among the AM processes. Some AM selection process also used benchmarking approach which is a costly and lengthy procedure. Additionally, all the vendors do not agree or conform to benchmarking of AM process. Some authors also compared the polymer based system to metal based system which has different surface finish and accuracy. Most of the MADM employed in earlier works assumed accurate judgement and crisp evaluation.
Henceforth, there is a need of a comprehensive but simple and logical approach involving more processes and attributes for the selection of AM. This work is different from previous works due to its exhaustive and extensive nature. This study has incorporated three AM systems, six AM technologies and thirty nine different AM machines with nine machine attributes (criteria). It intends to provide a most appropriate AM machine among the three AM based systems, different technologies and all the available AM machines in the database. Verily, it is imperative that the selection of a particular machine from an inclusive database always involve significant imprecision. This imprecision may arise from various sources such as unquantifiable data, deficient information, personal biases, inaccessible knowledge, and limited expertise, especially on the part of the users [44, 45]. Generally, most of the stakeholders or the system users are likely to depend on their knowledge and personal skill during machine selection, which lead to highly unorganized and uncertain decisions [46]. As a result of diverse sources of uncertainty, combination of fuzzy and crisp data as well as abundance of attributes and machines, a decision support methodology based on uncertainty theories, FAHP and GRA have been introduced in this work for the selection of a particular RP system. The FAHP (extent analysis based on trapezoidal fuzzy numbers) has been utilized to estimate the weights of various attributes, while the GRA was executed to achieve the best RP system and technology. The FAHP and GRA have been chosen owing to their ease of use, wide applicability as well as their robust nature and well-established outcomes in a number of real life decision problems. The FAHP and GRA have leverage over the other methods because they are adequate and competent in quantifying the uncertainty of information as well as the vagueness of judgements, can handle innumerable criteria and alternatives. Fuzzy set theory and GRA are very powerful mathematical tools for modeling and control of uncertain systems, while AHP represents a robust and a flexible framework for solving complex decision problems. This system can be identified as a user-friendly decision advisor which permit partial information as well as qualitative and quantitative data. As a result of growing complexity in decision-making, most of previous MADM approaches have only considered limited criteria and alternatives. However, in this work, most of the known criteria associated with AM selection have been analyzed. Furthermore, multi-level classification of alternatives have been considered for determining the AM process, technology and machine according to the requirement of the users. This approach is adept to minimize subjective or cognitive errors by streamlining, segregating and correlating multitude of attributes. The application of FAHP would shape up the opinions of the experts as well as enhance the quality of gathered data, without restricting them to distinct terminology or parameter. The FAHP and GRA based integrated approach can result in flexible and sustainable solutions as well as will expand the prospects of their future applications.
Methodology
This section discusses the different steps of the adopted methodology in order to reach out the best RP solution.
Generation of database
The foremost step in this methodology was the classification of different RP machines depending on the principles and raw material used. This work has only considered polymer based AM technologies. The different polymer based AM machines can be categorized into three systems depending on the build material (or the raw material) used. The different systems were liquid based system, solid based system and the powder based system. Each of these AM systems consisted of different technologies or processes as shown in Table 1. The six AM technologies, two each for each AM systems, have been considered in this research study. The different technologies were SLA, PolyJet technology, FDM, Multi-Jet Printing (MJP) or Multi-Jet Machining (MJM), SLS and Inkjet Printing (INKJET), etc. These technologies were chosen as they are widely used in industries and manufactured by five major manufacturing firms: 3D systems, Stratasys, Z-Corporation, Object geometrics Ltd. and Electrical optical systems (EOS). The number of machines (or different models) chosen for each technology was based on their commercial availability in the market as well their popularity among the masses.
AM process classification
AM process classification
The selection of the most appropriate RP system becomes especially crucial because each machine or model possesses its benefits and limitations. For example, the dimensional accuracy and speed of SLA produced parts is good, but they require support structures and post processing operations [47]. Similarly, the PolyJet can print multiple materials simultaneously, however they also need support structures for overhanging parts and holes [48]. Moreover, the FDM produced parts are durable, functional and good in strength, while MJM offers the highest Z direction resolution and offers best combination of resolution and print speed [49, 50]. The build material used in SLS is in powder form and their fabricated parts are robust, durable and exhibits higher wear resistance as well as they do not require support material [51]. The INKJET originally developed by MIT in 1993 is capable of printing parts with multiple colors and offers faster speed when compared to other technologies of AM systems [52]. Therefore, it is very important to gather as much information as possible about the different machines and identify the best system which can satisfy most of the requirements of the users. The database was generated in four stages. The first stage involved the survey of different AM processes and technologies available in the market. This survey mainly focused on the high-end level, medium level and small machines of the major players in AM industry. The second stage consisted of the collection of data which took months as there were more than 2000 data (or specifications) files of AM machines on public domain. The third stage included the building of database using technical features (machine characteristics) of the AM machines available in the market. The database can be expanded when new AM technologies will be added with the passage of time. The database might not be exhaustive, but includes important players of AM Industry. The fourth stage comprised of the evaluation attributes needed in any MADM. The entire database can be summarized in a hierarchical structure with five levels as shown in Fig. 1. The level 1 represents the objective of selecting the best machine among thirty nine AM machines. The level 2 attributes consist of AM systems including liquid, solid and powder based material. The level 3 presents different AM technologies such as SLA, PolyJet, SLS, etc. The level 4 represents the machine characteristics for the thirty AM machines used in this study. The level 5 consist of thirty nine different AM machines assigned to these technologies.

Hierarchy of AM database.
The characteristic data (attributes) of thirty nine AM machines as shown in Table 2, was obtained from the Wohler’s reports, company catalogues that manufacture these AM machines, vendors and literature survey, thus making it highly reliable. The essential source of information was web-based or internet. A significant number of important attributes including the machine cost, error (or the accuracy), layer thickness, machine speed, material cost, net build size, machine weight, surface roughness and material strength were considered in this work.
Machine characteristics (attributes) data
A group of RP users or experts with different requirements were identified from academia and industries. Their requirements were gathered and a decision matrix was prepared using Saaty’s intensity scale as shown in Table 3. The different users or experts and their requirements can be described as follows.
Saaty’s intensity importance scale [53]
Saaty’s intensity importance scale [53]
User 1 – Industrial personnel who needed a solid based system, especially MJM. The crucial machine characteristics for this user were accuracy, machine speed, surface finish, material strength in the same order.
User 2 – Academician who required a solid based system, such as FDM for demonstration purposes. The machine cost, machine weight, surface finish and material strength were important attributes to accomplish his objectives.
User 3 – A user from a manufacturing company who was looking for a liquid based system, particularly SLA with higher surface finish, lower machine weight and enhanced material strength.
User 4 – A researcher who was keen on liquid based system (exclusively PolyJet) for experimentation. This user required a machine with lower material and machine cost, lesser machine weight and a higher net build size.
The three decision matrices for each level (machine characteristics, AM systems and AM technologies) were obtained as the outcome of this step. Three matrices included decisions regarding machine characteristics, AM technologies and AM systems. In order to model unbiased and logical behavior, the consistency analysis of the judgements was mandatory. Therefore, before processing further, the data obtained from different users was validated to conform their consistency [53]. To simplify the computation, the consistency test proposed by Saaty [54] which made use of crisp values was implemented. The Saaty’s consistency test was acceptable to ascertain that the judgements were rational and dependable to a reasonable extent, even when the fuzzy set theory was integrated with AHP. Indeed, the same procedure was also adopted by Amarasingha and Piantanakulchai [55], Radionovs and Užga-Rebrovs [56] as well as Meixner [57] in their analysis to verify the consistency of acquired judgements. A random consistency index (RI) as shown in Table 4 was employed to test the consistency of decision making. A consistency ratio, CR which compares between consistency index (CI) and RI has to be defined as follows.
Random consistency index
The λmax indicates the maximum eigenvalue of the decision matrix and n represents the number of attributes. If CR is less than or equal to 10%, the inconsistency is acceptable else the users have to revise their decision matrix.
In this research, the uncertainty and vagueness of expert’s opinions as well as of data acquisition have been managed through fuzzy sets theory proposed by Zadeh [58]. The FAHP based on the concepts of fuzzy set theory and hierarchical structure analysis can methodically deal with the RP machine selection problem [59, 60]. Indeed, the FAHP is an extension of a typical AHP technique into fuzzy environment through fuzzy numbers [61]. The original transformation of standard AHP into FAHP was proposed by Van Laarhoven and Pedrycz [62], Buckley [63] and Chang [64]. They represented the Saaty’s relative preferences by using fuzzy numbers with triangular membership functions. However, the crisp values in this work have been modeled using trapezoidal membership function owing to its superior performance as compared to triangular membership functions [65]. The different steps used in the estimation of fuzzy weights can be discussed as follows.
The decision matrices obtained for different users were transformed into trapezoidal fuzzy numbers in this step. The membership function of trapezoidal fuzzy number designated by F and defined as (l, m, n, u) can be described as follows [66].
In order to achieve this transformation, the crisp values based on Saaty scale were first transformed into triangular fuzzy numbers and then to the trapezoidal fuzzy numbers. The triangular fuzzy numbers were converted into trapezoid numbers by keeping the upper and lower bounds of the fuzzy numbers constant and extending the center to a certain range [67]. For example, the crisp value of x was converted to triangular fuzzy number (a, b, c) as a = x - 1; b = x and c = x + 1. Subsequently, the triangular fuzzy numbers were reshaped into trapezoidal fuzzy number (l, m, n, u) as l = a; u = c; m = l + 0.5 and n = u - 0.5.
In this study, priori weights of the performance measures were computed by using extent analysis of Chang [64]. The extent analysis for computation of priori weights can be described as follows [68, 69]. Let X = {x1, x2, x3, …… , x
n
} represent the set of objects and G = {g1, g2, g3, …… , g
n
} be an objective set. According to extent analysis approach, every object has to be subjected to extent analysis for each objective of the problem. As a result, m extent analysis values for each object can be obtained as follows.
Where,
The synthetic fuzzy values with respect to ith object can be computed by employing the Equation (1).
Since, FAHP utilizes four values (trapezoidal fuzzy number) for the evaluation of the particular attribute. Therefore, to compute
The value of
Suppose M1 = (l1, m1, n1, u1) and M2 = (l2, m2, n2, u2) represents two trapezoidal fuzzy numbers. The degree of possibility of M2 = (l2, m2, n2, u2) ⩾ M1 = (l1, m1, n1, u1) is dependent on conditions in Equation (5).
The degree of possibility of a fuzzy number greater than k fuzzy numbers can be defined as
Suppose that
The weight vector can be expressed as follows.
, where A i (i = 1, 2 3 4, …, n) are n elements. To obtain the weight vectors for individual elements, the normalization has to be carried out. The weight vectors are normalized to get the normalized vectors.
w = (d (A1) , d (A2) , d (A3) , d (A4) , …, d (A n )) T , where w is a vector of non-fuzzy numbers or crisp values.
In this work, three sets of weights (w1, w2 and w3) were computed using the above discussed FAHP as shown in Fig. 2. The first set of weights (w1) was computed using the decision matrix acquired for different attributes. These weights (represented by local weights 1) were used to determine the best machine within the respective technologies (SLA, PolyJet, FDM, MJM, INKJET and SLS). For example, the best SLA machine, the best MJM machine, etc. Consequently, the second set of weights (w2) were determined by using the decision matrix of different technologies. In order to estimate the most suitable AM system (liquid, solid or powder based system) depending on the requirements, local weights 2 (=local weights 1 * w2) were utilized as shown in Fig. 2. Finally, the decision matrix of AM systems was employed to compute the third set of weights (w3). Then the global weights (=local weights 2* w3) were determined and these were implemented to achieve the most appropriate RP machine.

Representation of weights at various levels of the hierarchy.
As the weights were computed for different users, therefore, the individual weights were combined using the geometric mean (GM) approach as follows in Equation (7).
Where, d is the number of users or experts.
The GRA has been employed in this work to rank the various alternatives in the RP selection problem. The notion of grey theory was first introduced by Prof. Deng from the grey set in combination with theory of space and control theory [70]. The objective behind utilizing grey theory is its ability to take into account the ambiguity and uncertainty associated with user decisions and data collection. The execution of GRA can be carried out in the following manner [71, 72].
Suppose
Using the comparability sequences, a reference sequence was determined. After data processing was performed, grey relational coefficient (GRC) were estimated using the preprocessed sequences. The GRC was computed using the following equation Equation (8).
Where Δ
oi
(k) represented the deviation sequence of the reference sequence
Where β
k
represented the weighting value of the kth performance characteristic, and
The ranking was carried out for the individual users as well as for an imagined user whose requirements was the combination of the individual four users.
The objective of sensitivity analysis was to ascertain the robustness and stability of ranking acquired through the proposed approach. It can be defined as a technique to understand the variation of input values on the final model results [73]. Certainly, it was critical to examine the deviations in the outcome parameters of given model due to deviations in input parameter values [74, 75]. It suggests that the model is reliable and trustworthy if the model solution is not very sensitive to changes in input. The most commonly used sensitive analysis is accomplished by varying the weights of the performance attributes and analyzing the change in results [76]. In earlier studies also, it has been shown that the ranks of the alternatives significantly depend on the values of the weight coefficients of the attributes [77]. Therefore, in this work, sensitivity analysis depending on changes in weight coefficients was employed to validate the model and ensure the reliability of results.
In sensitivity analysis, the first step was the arbitrary selection of a criterion. It was followed by a percentage change in the selected criterion (more or less). Consequently, the computation of weights for remaining criteria was carried out using Equation (10) [80].
Where,
w i : Original weights for attribute i
weight by a percentage change for attribute i
w n : Original weight for attribute n
In order to establish the robustness of the developed approach and to examine the similarities of ranks obtained by varying the weights, Kendall’s coefficient (Z) of concordance was utilized [81, 82]. The Kendall’s coefficient depicts the similarities in the ordering of ranked quantities and its values ranges from 0 to 1. For example, the value of 1 represents a perfect match between the various ranking orders. It suggests that closer the value of Z to 1, the better will be the robustness. The value of Z can be computed by using Equations (11–13).
Where, r i j= ranking scenario i provides to alternative j, m = number of scenarios and k = number of alternatives.
Altogether, four schemes depending on weight percentage (20%, 30%, 40% and 50%) were established for each user (decision maker). In each scheme, there were six scenarios based on variously selected attribute for percentage change. The different schemes and scenarios can be realized in Table 5. As an example, the scheme 2 and scenario 3 represents a 30% decrease in weight for surface finish (weights of remaining attributes were computed using Equation (10)), 30% increase in weight for SLS (weights of remaining technologies were calculated using Equation (10)) and 30% increase in weight for POWDER (weights of remaining systems were estimated using Equation (10)).
Scenarios used in sensitivity analysis
The database established for the RP selection problem consisted of the type of AM system (liquid, solid, powder), principle or the technology (SLA, PolyJet, FDM, MJM, 3DP, SLS), manufacturer, machine name and the machine code. As shown in Table 6, the database also included information about specifications of various machines. In the last column of Table 6, the value of three represented weak material strength, six expressed the strong strength and nine was used to show the very strong material strength.
Database developed for the selection of best RP machine
Database developed for the selection of best RP machine
In the consecutive step, the decision matrices were prepared for various users. The users were requested to measure their requirements on a scale of 1 to 9 (Saaty’s scale) needed in their particular applications. As result, three decision matrices were obtained for each user. The Table 7 represents the decision matrices for the user 1.
Decision matrices for the user 1: (a) RP system; (b) RP technology; (c) Machine attributes
The consistency was assessed for each user and the results are presented in Table 8. The outcome form consistency analysis implied that the value of the CR in each decision matrix did not exceed the maximum permissible value of 10%, therefore, these opinions from different users was acceptable in consequent computations.
Results obtained using consistency analysis
Subsequently, the crisp values gathered within each decision matrix for every user were transformed into trapezoidal fuzzy numbers. The trapezoidal fuzzy numbers for different intensities of importance are shown in Table 9.
Conversion of Saaty’s importance scale to trapezoidal fuzzy numbers
After transformation into trapezoidal fuzzy numbers, Chang’s extent analysis was employed to compute the weights for individual users. The different calculation steps adopted in Chang’s method are discussed as follows. The decision matrix for machine attributes by user 1 is utilized as an example in this calculation.
Fuzzy pairwise comparison matrix
S1= (5.218, 6.3773, 8.911, 10.45) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
S2= (38, 41.5, 48.5, 52) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
S3= (8.0933, 9.744, 13.257, 15.283) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
S4= (29.333, 32.9, 40.167, 44) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
S5= (3.329, 3.9953, 5.544, 6.593) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
S6= (11.9503, 14.09, 18.575, 21.083) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
S7= (2.44, 2.6064, 3.155, 3.704) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
S8= (22.583, 25.686, 32.067, 35.5) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
S9= (16.7833, 19.4082, 24.853, 27.833) ⊗ (1/216.446, 1/195.029, 1/156.3072, 1/137.7299)
V (S1 ⩾ S2) = 0.3804, V (S1 ⩾ S3) = 0.9378, V (S1 ⩾ S4) = 0.4717, V (S1 ⩾ S5) = 1, V (S1 ⩾ S6) = 0.8838, V (S1 ⩾ S c 7) = 1, V (S1 ⩾ S8) = 0.5855, V (S1 ⩾ S9) = 0.7229
V (S2 ⩾ S1) = 1, V (S2 ⩾ S3) = 1, V (S2 ⩾ S4) = 1, V (S2 ⩾ S5) = 1, V (S2 ⩾ S6) = 1, V (S1 ⩾ S7) = 1, V (S2 ⩾ S8) = 1, V (S2 ⩾ S9) = 1
V (S3 ⩾ S1) = 1, V (S3 ⩾ S2) = 0.5534, V (S3 ⩾ S4) = 0.6597, V (S3 ⩾ S5) = 1, V (S3 ⩾ S6) = 0.9244, V (S3 ⩾ S7) = 1, V (S3 ⩾ S8) = 0.7822, V (S3 ⩾ S9) = 0.922
V (S4 ⩾ S1) = 1, V (S4 ⩾ S2)=0.911, V (S4 ⩾ S3) = 1, V (S4 ⩾ S5) = 1, V (S4 ⩾ S6) = 1, V (S4 ⩾ S7) = 1, V (S4 ⩾ S8) = 1, V (S4 ⩾ S9) = 1
V (S5 ⩾ S1) = 1, V (S5 ⩾ S2) = 0.2064, V (S5 ⩾ S3) = 0.83, V (S5 ⩾ S4) = 0.2735, V (S5 ⩾ S6) = 0.6432, V (S5 ⩾ S7) = 1, V (S5 ⩾ S8) = 0.3677, V (S5 ⩾ S9) = 0.4893
V (S6 ⩾ S1) = 1, V (S6 ⩾ S2) = 0.7141, V (S6 ⩾ S3) = 1, V (S6 ⩾ S4) = 0.8272, V (S6 ⩾ S5) = 1, V (S6 ⩾ S7) = 1, V (S6 ⩾ S8) = 0.95, V (S6 ⩾ S9) = 0.9162
V (S7 ⩾ S1) = 0.755, V (S7 ⩾ S2) = 0.0487, V (S7geqslantS3) = 0.5368, V (S7 ⩾ S4) = 0.0856, V (S7 ⩾ S5) = 0.9929, V (S7 ⩾ S6) = 0.3659, V (S7 ⩾ S8) = 0.1502, V (S7 ⩾ S9) = 0.2404
V (S8 ⩾ S1) = 1, V (S8 ⩾ S2) = 0.9824, V (S8 ⩾ S3) = 1, V (S8 ⩾ S4) = 0.9037, V (S8 ⩾ S5) = 1, V (S8 ⩾ S6) = 1, V (S8 ⩾ S7) = 1, V (S8 ⩾ S9) = 1
V (S9 ⩾ S1) = 1, V (S9 ⩾ S2) = 0.8576, V (S9 ⩾ S3) = 1, V (S9 ⩾ S4) = 0.9713, V (S9 ⩾ S5) = 1, V (S9 ⩾ S6) = 1, V (S9 ⩾ S7) = 1, V (S9 ⩾ S8) = 0.9109
d(Machine Cost) = min (0.3804, 0.9378, 0.4717,1, 0.8838, 1, 0.5855, 0.7229) = 0.3804
d(Accuracy) = min (1, 1, 1, 1, 1, 1, 1, 1) = 1
d(Min Layer Thickness) = min (1, 0.5534, 0.6597, 1, 0.9244, 1, 0.7822, 0.922) = 0.5534
d(Machine Speed) = min (1, 0.911, 1, 1, 1, 1, 1, 1) = 0.911
d(Material Cost) = min (1, 0.2064, 0.83, 0.2735, 0.6432, 1, 0.3677, 0.4893) = 0.2064
d(Net Build Vol.) = min (1, 0.7141, 1, 0.8272, 1, 1, 0.95, 0.9162) = 0.7141
d(Machine Weight) = min (0.755, 0.0487, 0.5368, 0.0856, 0.9929, 0.3659, 0.1502, 0.2404) = 0.0487
d(Surface finish) = min (1, 0.9824, 1, 0.9037, 1, 1, 1, 1) = 0.9037
d(Material Strength) = min (1, 0.8576, 1, 0.9713, 1, 1, 1, 0.9109) = 0.8576
The weight vector can now be expressed as w = (0.3806, 1, 0.5534, 0.911, 0.2064, 0.7141, 0.0487, 0.9037, 0.8576)
After normalization, the weight vectors can be represented as w = (0.0682, 0.1794, 0.0993, 0.1634, 0.0370, 0.1281, 0.0087, 0.1621, 0.1538)
Similarly, the remaining weights for other decision matrices were computed. The individual weights were also combined and normalized as shown in Table 10 using GM method. The local weights 1, local weights 2 and global weights were also estimated as discussed in the previous section.
Once the weights were computed, the GRA was implemented to determine the best alternative for different users as well for the imagined user whose requirement was the combination of all the individual users. The implementation of GRA commenced with data normalization to obtain comparability sequences, subsequently deviation sequences and GRC. As an example, the calculation steps of GRC for Machine Cost (LSLA3DS1M1) are shown below. Since, the objective is the minimization of machine cost, therefore, the smaller the better condition will be utilized in normalization.
Max = 1200000; Min = 10900; Given value = 220700
Similarly, the remaining data were normalized depending on the conditions of maximization or minimization. Consequently, the deviation sequences were computed as follows.
Δ
oi
(k) = |1 - 0.8236|
Δ oi (k) =0.1764
The GRC was computed using Equation (8) and by incorporating the following values.
∈=0.5; Δmin = 0; Δmax = 1; Δ oi (k) =0.1764
GRC=(0 + 0.5*(1))/ (0.1764+(0.5*1)) GRC = 0.7392
For user 1, machine cost (machine attribute), SLA (technology) and Liquid (System), the following weights were obtained (please refer Fig. 2).
Local weight 1 = 0.0682
Local weight 2 = 0.0682*0.0768 = 0.005238
Local weight 3 = 0.005283*0.0838 = 0.0004389
Similarly, the GRC s for other combinations of machines, technologies and systems were carried out. The Tables 11–13 represent the comparability sequences, deviation sequences and GRC respectively.
Weights computed using FAHP
Comparability after data normalization using GRA
Deviation Sequences after data pre-processing
After estimating the GRC, the GRG and the respective ranking were obtained for different users and the combined user utilizing their respective weights. Using Equation (9), the values of GRG were computed to rank various machines, systems and technology. The GRG (k = 1) for machine cost, SLA and Liquid combination is as follows.
GRG1 (based on different technologies)
=0.0682*0.7392
=0.0504
GRG1 (based on build material)
=0.005238*0.7392
=0.003872
GRG1 (overall)
=0.0004389*0.7392
=0.0003245
In order to rank the various alternatives, the inclusive GRG was estimated by summation of GRGs for different attributes (k = 1, 2, ... , n) using Equation (9). For example,
GRG (within AM technology) = (0.0504 + 0.1748 + 0.0517 + 0.0557 + 0.0204 + 0.0432 + 0.0078 + 0.1430 + 0.0513) = 0.5983.
GRG (within AM system) = (0.0039 + 0.0134 + 0.0040 + 0.0043 + 0.0016 + 0.0033 + 0.0006 + 0.0110 + 0.0039) = 0.0460.
GRG (overall machine) = (0.000325 + 0.001126 + 0.000333 + 0.000358 + 0.000132 + 0.000278 + 0.000050 + 0.000921 + 0.000330) = 0.003853.
Similarly, the GRG for other machines, technologies and systems were achieved. The Tables 14–16 represents the GRG and the respective ranking of machines in a particular technology, AM systems and the most appropriate RP technology respectively for the user 1.
GRC computed based on the deviation sequences
GRG and ranking of the machines within various technologies
GRG and ranking of AM systems (build material classification)
Based on the analysis using FAHP and GRA, the machines or systems as shown in Table 17 can be suggested to various users and the fictitious combined user. It can be realized that the integrated approach based on FAHP and GRA is able to satisfy all the users. For instance, PROJECT 5000 from 3D Systems, which is a solid based system based on MJM technology has been recommended to the User 1 by the RP selection system. Moreover, for a fictitious user with combined requirements, a machine with model DIMENSION 768 SERIES (solid based system and FDM technology) can be suggested.
GRG and ranking to determine the most appropriate RP system
Similarly, the DIMENSION 768 SERIES from STRATASYS (solid based FDM system), I PRO 9000 XL from 3D systems (liquid based system, SLA) and 24 PERSONNEL 3D PRINTER from OBJECT GEOMETRICS (liquid based system, PolyJet) were recommended to the user 2, user 3 and user 4 respectively.
The outcomes from sensitivity analysis as shown in Table 18 implies that the ranking order established using FAHP and GRA, is robust and stable. It can be attributed to the fact that the value Z for each user was closer to 1. To compute the values of Z, the ranking order obtained for six scenarios in each scheme and primarily established ranking order (in this work) were analyzed for their similarirties. The higher values of Z suggest that the ranking order established in this work is reliable and is not very sensitive to changes in weight.
Appropriate outcomes of the RP selection system
Outcomes obtained in sensitivity analysis
A decision advisor based on FAHP and GRA has been utilized to methodically deal with an unstructured and a complex RP selection problem. This method has the competence to consider the unclearness of human thinking and expertly solve MCDM problems. This approach has the potential and the capability to overcome any inconsistency and uncertainty with increase in the number of attributes. Certainly, the decisions made by the developed RP selection advisor successfully convinced every user and adequately satisfied all their requirements. This research is distinctive owing to its comprehensive database, categorization of RP selection procedure into three main levels, application of trapezoidal fuzzy numbers in extent analysis technique (FAHP), conversion of crisp values into trapezoidal fuzzy numbers and the utilization of GRA to rank various alternatives. As a result of its unique aspects, this method is easy to use and can yield rational, competent and sustainable results. Since, the outcomes in this work have been computed keeping in mind the available uncertainty in decision making and ambiguity in data collection. Therefore, it can be asserted through sensitivity analysis that the weights and the ranking were robust and effective. Moreover, many attributes have been considered in this investigation. Indeed, the complicated problems can be represented conveniently with this decision advisor which is simple and flexible in nature as well as that can exploit both objective and subjective information. The adaptability and resilience of the final outcome can also be investigated through sensitivity analysis. The utilization of FAHP weights and GRA ranking made the selection procedure adaptable, reliable and more realistic. This model is especially useful in RP selection or other similar problems, where the extent of uncertainty are expected to vary due to the dynamic and competitive market in conjunction with fickle customer demands. However, this selection approach may become tedious, computationally expensive and time consuming with increase in the number of hierarchical levels and the pairwise comparisons. Moreover, since the effectiveness of this decision approach depends on the judgement of users, it may sometime tends to overestimate the ranking process. In addition, the consistency of fuzzy pairwise comparison matrix with trapezoidal fuzzy numbers should be studied comprehensively. As the problem of RP selection has not been dealt with efficiently and effectively due to its complex nature, this work intends to aid RP users in selecting the most appropriate system. This work would inevitably help the vendors and industrial firms as a guideline in selecting the most appropriate RP system. Apart of RP selection, the proposed method can also be employed to other managerial applications such as project selection, appropriate supplier identification, hospitality industry, infrastructure management, etc. On account of continuously evolving RP market and dynamic customer demands, the database has to be updated periodically with addition of new technologies and machines. Therefore, in the future work, more attributes and alternatives will be considered to develop a more effective RP selection system. The authors also plan to solve the presented RP selection problem using other decision models, such as the method based on unbalanced hesitant fuzzy linguistic term sets [78], computational model based on extended linguistic hierarchies [79], etc.
Footnotes
Acknowledgments
The authors extend their appreciation to the Deanship of Scientific Research at King Saud University for funding this work through research group number RG-1440-034.
