# Evolution of Computer-Aided Process Planning for Hybrid Additive/Subtractive Process

Authors: Osama Abdulhameed; Abdulrahman Al-Ahmari; Syed Hammad Mian; Abdulmajeed Dabwan; Hisham Alkhalefah
Journal: Advances in Materials Science and Engineering
 (2020)

---

## Abstract

The hybrid process, which integrates two or more different processes such as additive manufacturing and subtractive manufacturing, has gained appreciable considerations in recent years. This process exploits the benefits of individual processes while overcoming their limitations. Lately, the combination of additive, subtractive, and inspection methods is a valuable conglomeration, considering its potential to produce complicated components precisely. Certainly, computer-aided process plan (CAPP) provides a crucial link among different processes and is essential to avail the benefits of hybridization. However, a valuable process plan can only be achieved through the optimization of its different elements. Therefore, the objective of this work is the accomplishment of an optimized CAPP to fabricate parts in the shortest time employing the hybrid additive, subtractive, and inspection processes. In this work, mathematical models have been developed to optimize part orientation as well as minimize additive and subtractive times. Additionally, the genetic algorithm has been employed to obtain the best path with minimum inspection time. The feasibility and capability of the proposed approach as well as the optimized CAPP for the hybrid process have been demonstrated through a case study.


---

## Body

## 1. Introduction

The growth and development of the manufacturing business depend on the ingenious studies related to the optimized utilization of manufacturing processes, materials, and product design. The highly competitive market, as well as the requirements of complex and accurate parts at a reasonable price, has given rise to the concept of hybridization. The hybrid manufacturing (HM) can be described as the fabrication approach, wherein the combination of manufacturing processes is integrated. HM helps to enhance the performance of individual processes while providing a compounded advantage from these processes. It is crucial to figure out the principles, advantages, and shortcomings of discrete processes for realizing HM. In recent times, the hybridization of additive and subtractive manufacturing is the most demanding due to its numerous benefits.A subtractive manufacturing (SM) process such as computer numerical controlled (CNC) material removal machine can be integrated with the additive manufacturing (AM) process to overcome the constraints such as tool inaccessibility, higher error, degraded surface conditions, and lower building rate [1]. The combinations of AM and SM processes are increasingly being used to enhance manufacturing performance in terms of higher tool life, lower building time, effective machining with close tolerances, and enhanced surface finish. Although HM exhibits many benefits, its applications are limited as a result of the scarcity of computer-aided process plan (CAPP) approaches. Moreover, the few CAPPs which are existing have not been optimized to achieve effective and efficient HM. Indeed, the research works have so far focused on improving the individual processes instead of applying a comprehensive and universal perspective to optimize HM [2, 3]. Consequently, an effective and efficient CAPP is needed to be developed which can effectively shorten the fabrication time. As reported by many researchers, there is a significant lack of CAPP techniques, which can adequately utilize manufacturing resources to accomplish hybridization.This work develops a hybrid process integrating AM, SM, and the measurement process. A new process-planning system has been proposed for the production of prismatic parts in the shortest manufacturing time. The developed CAPP would provide enough flexibility to adapt to any modifications in production assets, machining settings, etc. on the working area as well as allow minimization of part production time. In this work, the CAPP has been optimized to improve the hybridization of AM, SM, and the inspection processes. The mathematical models have been developed to obtain the best part orientation as well as minimize the additive and subtractive times. Additionally, the genetic algorithm (GA) has been implemented to achieve the best path with minimum inspection time. A case study has also been discussed to depict the performance of the optimized CAPP in the hybrid process.


## 2. Literature Review

AM parts usually possess a poor surface quality; therefore, they are machined to produce the smooth surface finish. AM process can be integrated with the SM process to produce the finely finished parts. Largely, CNC milling is the frequently used procedure to process the near-net shapes obtained from AM. This is because the flat surface obtained using an uncomplicated rotary-tool polishing or grinding does not meet the norms needed for precise and accurate parts [4]. Therefore, advanced postprocessing using CNC polishing and milling, ultrasonic machining, glass blasting, etc. are required to produce finely finished parts [5]. There have been numerous efforts, where AM and SM have been combined to overcome the limitations of each other. Karunakaran et al. [6] explored the benefits of hybridization between CNC and AM. They combined a CNC machine and arc weld deposition to enhance the quality of the product and reduce production span and overall manufacturing cost. In the proposed setup, the layers were accumulated using metal inert-active gas welding and three-axis milling was employed to obtain the finely finished surfaces. LUMEX represented a single setup, which combined a fiber laser for metal sintering and a machining center for high accuracy and high-speed milling [7]. This system could fabricate complicated mold dies in lesser time and reduced production cost. Similarly, Kulkarni and Dutta [8] overcame the limitations of AM and SM processes through an integrated production setup. A hybrid process based on the integration of selective laser melting and CNC milling was employed by Kashouty et al. [9] to achieve desired accuracy and precision in the manufacture of an automobile tool insert. To comprehend the benefits of HM, Ren et al. [10] introduced a unified process-planning system considering a single setup combining AM and SM. It consisted of the extraction of information from the computer-aided design (CAD) model, optimization of the tool path, and collision detection algorithms. Liou et al. [11] and Ruan et al. [12] also developed a process-planning approach for laser sintering and CNC machining. Their approach was based on dividing the parts in nonuniform layer thickness and accomplishing collision-free tool paths.Few works that focused on optimizing the process plans can also be pointed out in the literature. Although, there has been a significant scarcity of optimized process plans for hybrid processes. Most of the efforts directed to the optimization of process plans were dedicated to individual processes or manufacturing systems. For instance, Shabaka and ElMaraghy [13] developed a GA-based approach to optimize the manufacturing cost of the process plans in reconfigurable manufacturing systems (RMSs). The researchers identified various parameters such as machine appointment, machine structure, operation order, operation cluster sequencing, selected tools, and tool approach directions (TADs) of the operations, etc., to obtain the mathematical models. Likewise, Jia et al. [14] utilized the object-oriented representation to accomplish manufacturing information, process decision insight, and decision procedure control knowledge. Some works that have concentrated on the optimization of process plans for the hybrid process can also be listed out. Ahsan et al. [15] developed an innovative two-step optimization methodology to figure out the appropriate build direction for the object and material accumulation direction for layers. They considered the lowest contour plurality, surface characteristics, build height, fabrication parameter, and layer contour concavity to overcome the limitations associated with fabrication and resource availability. The developed methodology was carried out on a case study, and the outcomes indicated a significantly reduced total build time in contrast to the total build time related to subjective build and deposition directions. Similarly, Um et al. [16] proposed a description method of HM and developed process plans at both macro- and microlevels. The stereolithography (STL) data, which are utilized to fabricate the part by AM, comprise a large number of errors, owing to several geometric processes. Therefore, the representation of solid models utilizing the Standard for the Exchange of Product model data (STEP) maintains the error as low as possible during geometric reasoning. The realization of geometric reasoning makes the process of planning automated with high accuracy. Moreover, Le et al. [17] demonstrated the utility of aggregating two technologies. The proposed approach aimed at reusing existent parts directly to produce final parts, evading the material recycling stage. This paper particularly aimed at designing the process plan for a composite AM and SM processes. Aboutaleb et al. [18] proposed an integrated framework to systematically accomplish suitable fused filament fabrication (FFF) process parameters, which resulted in optimal geometry integrity. The developed method was verified through a real-world case study. The outcome confirmed that optimal properties were attained more efficiently as compared to other existing methods. Reiher and Koch [19] considered part orientation as well as quality management guidelines to accomplish a suitable production plan. They showed numerous methods and measures to attain product optimization, information acquisition, decision-making, and shape identification for part optimization in AM. Udroiu and Nedelcu [20] proposed an optimization approach which aimed at three AM applications. The first application was the part orientation considering the minimum build time, minimum support structure, and surface accuracy. The second application was targeted at fitted validation of parts obtained using three-dimensional (3D) printing. The study by Byun and Lee [21] aimed at the recognition of optimal buildup direction for different rapid prototyping (RP) systems. The computation of surface quality considered the layering effect, building time and part cost, labor cost rate, material cost, etc. The developed algorithm was utilized to obtain the most appropriate buildup direction and generate an optimal process planning.Certainly, the research studies presented above have contributed immensely to the literature on hybrid techniques and their process planning. However, as shown in Table1, these works were adequate and acceptable, but they were not all-inclusive as well as limited in their utilization. Moreover, there has been a dearth of research works that have focused on the development of optimization approaches for process planning of hybrid additive and subtractive manufacturing methods. There has always been a requirement to develop an optimized CAPP system for the hybrid process that can adjust to volatile machining conditions. Therefore, in this work, an optimized CAPP has been proposed to fabricate parts in the shortest time utilizing the hybrid additive, subtractive, and inspection processes. The mathematical models have been employed to optimize part orientation as well as minimize additive and subtractive times. Besides, the GA has been utilized to obtain the best path with minimum inspection time.Table 1
Distinction of present contribution in comparison to recent related works.
Recent related research workCurrent workReferenceContribution and limitationsKamrani et al. [22]A framework was developed to generate an instruction file for guiding the CMM measurement process. This work primarily emphasized on developing a computer-aided process plan for the inspection process. The feature extraction, feature recognition, and probe approach directions were utilized for ease of access, setup planning, and inspection point’s distribution on the measurand. The techniques such as the clustering algorithm, graphical method, and ANNs were used to carry out the competent inspection. The proposed approach was effective and competent, but it did not deal with the manufacturing aspects. The stated process plan did not consolidate the manufacturing processes, such as additive and subtractive processes.This work is distinctive because it integrates additive, subtractive, and inspection processes. The proposed CAPP is holistic and comprehensive, which incorporates manufacturing as well as inspection information in the detailed process plan. The plan is flexible to adjust to any changes in production resources, machining, and inspection settings, etc.Homar and Pušavec [23]A processing time- and cost-oriented algorithm was reported to minimize the fraction of the component that has to be fabricated additively. A novel algorithm based on an optimized manufacturing operation sequence was utilized for feature recognition in hybrid processing (subtractive and additive). The authors presented an ingenious concept but did not recognize inspection, which has now become an integral link in the production process. Furthermore, this approach can be inefficient if the part is very complex and the product requires an additive process for the majority of its features.In this work, the inspection activity has been an indispensable element of the entire methodology. The additive, subtractive, and inspection processes have been synchronized to attain an efficient process plan. The established approach is formative and can be utilized for prismatic parts with any complexity.Al-Ahmari et al. [24]The relationship between the attribute and build orientation was studied in this work. An optimal part orientation system was established to improve part quality/accuracy in additive manufacturing. In this approach, the weighted sum function was employed to compute the part volume and time required for printing. This work was enlightening, but it was limited only to additive manufacturing and didn’t describe the proposed algorithm for its application in hybrid additive-subtractive-inspection process.The existent research incorporated and unified the developed mathematical models to enhance the effectiveness of the extensive process plan for the hybrid additive-subtractive-inspection process.Abdulhameed et al. [25]The concept of the multiagent system was proposed to develop a computer-aided process plan for the hybrid additive-subtractive-inspection process. It was based on a three-step methodology, namely, a static process plan, selection of optimum scenario, and dynamic process plan. The crux of this work was the application of agent theory in acquiring a proficient computer-aided process plan for the hybrid process. However, this work did not describe the optimization of individual phases. Although the inspection was an integral part of their methodology, the minimization of inspection time was not considered. Moreover, this work primarily focused on the software front and did not present the actual setup for hybrid machining.The present work has also introduced a computer-aided process plan for the hybrid additive-subtractive-inspection process, but it does not consider the agent theory. Instead, the idea of peaks is implemented to slice the entire part to subpart. Moreover, the optimization of each phase in the process planning of the hybrid process is the essence of this work. The genetic algorithm has been implemented in this study to minimize the inspection time, besides other optimization models. The setup designed for hybrid additive-subtractive-inspection process has also been referred to in this work.


## 3. Methodology

The methodology proposed in this work would enhance the product quality as well as reduce the manufacturing costs. The scheme explained in Figure1 represents specifically a generic overview of the hybrid additive, subtractive, and inspection system.Figure 1
Methodology of the hybrid additive, subtractive, and inspection processes.


### 3.1. Feature Extraction and Recognition

A rational consistent feature recognition approach has been employed to achieve feature recognition. This system possesses the competence to communicate with hybrid additive, subtractive, and inspection systems. The proposed methodology has been introduced for 3D prismatic parts designed by a solid modeling package (CATIA V5). The system considers a neutral file in STEP format as input and then interprets the information in the file to the manufacturing information. The ISO 10303 STEP Standards AP 203 E2 file has been utilized. The steps of the feature extraction and recognition methodology are as follows:Step 1. Obtain the geometry and topology items for the designed model using STEP file format.Step 2. Obtain topology items in each basic surface and determine its type.Step 3. Verify the feature’s presence in the basic surface depending on loops.Step 4. Determine feature type.Step 5. Determine the detailed features and obtain the related feature geometry parameters.Step 6. Determine the comprehensive machining data for each feature and the designed part.The sample output of a feature extraction file that contains a feature ID, feature name, feature face IDs, and feature dimensions is shown in Table2.Table 2
Feature recognition output [24].
Feature IDFeature nameFeature face IDsDimensions1Pocket throughF1, F2, F3, F4L =, W =, H =2HoleF5, F6H =, R =The rule-based approach as reported in [26] was adopted in this work to recognize the SM features. The following example represents the rule for recognizing the pocket through feature:(1)
For every 4 faces of the surface type plane (named as face 1, face 2, face 3, and face 4) as shown in Figure2(a).(2)
If concave edge count of the outer loop equals 2 for each face.2.1. If each consecutive pair of faces (F1F2, F2F3, F3F4, and F4F1) has one common concave edge of type Line and the pair of faces is perpendicular to each other.2.1.1. Pocket through found.2.1.2. Create a new pocket through object and add to the feature list.(3)
End ForFigure 2
Different faces in features: (a) pocket through; (b) slot blind feature.
(a)(b)For inspection, the data such as faces, edges, and vertices in the extraction file were classified and listed on the basis of their type and angles. Dimensional tolerances (DTs) were calculated from the direct limits of the extracted feature’s faces plus/minus tolerance values. For example, the feature length was calculated as the difference between F3 and F5, etc., as shown in Table3 and Figure 2(b).Table 3
Inspection data extracted for the slot blind feature.
Feature IDFeature nameFeature face IDsBlock face IDsDimensions1Slot blindF1, F2, F3, F4F5, F6Length = F3–F5, width = F1–F2, height = F4–F6


### 3.2. Part Orientation

The orientation rules for the improved surface accuracy depend on the relative orientation of various features/entities (in the extraction file) concerning the build direction. The extracted face orientation can be parallel, perpendicular, or at some angle. To satisfy better Geometric Dimensioning and Tolerancing (GD&T) [24], the objective functions of the optimal orientation of the part were modeled. The objective functions were the maximization of GD&T values and the minimization of production time. The automatic setup generation was achieved by gathering information from the feature extraction and recognition file as shown in Figure 3.Figure 3
Part orientation.In this diagram (Figure3), the setup class is consisted of get these methods, get all feature’s faces, get all part faces, remove interaction faces, generate frequency interaction edges, generate probe feature accessibility, and make feature clustering. The face class has consisted of get face ID, get face name, get face type, get face loop, get face edges, get face vertices, and get face direction.


### 3.3. Automatic Slicing Part to Subpart

The automated hybrid additive, subtractive, and inspection cannot be achieved without the automatic slicing of a part to a subpart. In this study, automatic slicing of a part to a subpart was achieved by dividing the CAD model into different CAD files using the CATVBA editor within CATIA V5. After slicing a part to a subpart, output data for every subpart was obtained in the STEP file format. In this automatic slicing, all data were divided automatically in an empty target file by running the code. The part peaks were the inputs to the slicing code. After that, the subpart area and volume parameters were calculated automatically and data were exported to an Excel datasheet. The exported data were used to calculate the total additive and subtractive area and volume of the part as shown in Figure4.Figure 4
Automatic slicing of a part to a subpart [24].


### 3.4. Simulation

In this work, a finite element analysis (FEA) framework (see Figure5) has also been developed for simulating the additive process. The objective of this analysis was to arrive at the constraints for process planning. The FEA approach in this work was established similar to the scheme and mathematical models reported in [27]. This analysis was subdivided across two phases. The thermal analysis, which solves the heat equations, was carried out to assess the temperature field. The resulting temperature field was then used in the mechanical analysis to determine von Mises stress and component deformation.Figure 5
Schematic diagram of FEA.The additive system employed in the current investigation was equipped with two extruders, which could be heated up to 600°C, and a nozzle of 0.2 mm of diameter was used. The bed temperature could reach up to 100°C. The different nozzle and bed temperatures, part orientations, extruder velocities, and layer thicknesses were applied to study the stress and deformations. The conditions were then identified to produce the most accurate part using additive manufacturing. This analysis also provided the regions (such as areas with high distortion) which could be prioritized for machining using the subtractive process. The FEA results were also utilized to define the sampling strategy required for inspection. For example, the highly deformed regions needed a higher number of points as compared to areas with low distortion.


### 3.5. Process Plan of Additive Operations

AM process-planning steps as shown in Figure6 can be interpreted as the layer-based manufacturing approach, which comprises of the resource requirements including time and travel distance. These aspects comprise the number of layers, shape and size, single or multiple contours in a slice, support material, functionality, accuracy, and surface characteristics [21]. AM process steps are evenly critical and possess a significant effect on the attributes of the produced part.Figure 6
Additive process plan steps.


### 3.6. Process Plan of Subtractive Operations

The subtractive process plan consists of the machining data of the designed part such as sequence of machining, type of operation, machine type, cutting tool, tool approach, and the volume of material removed for each extracted feature as shown in Table4.Table 4
Subtractive process plan.
Feature IDFeature nameFace IDsFeature dimensionFeature volumeFeature locationMachine operationCutting toolMachine toolSetups


### 3.7. Process Plan of Inspection Operations

The inspection planning can be summarized in Table5. The different attributes of the process plan table can be described as follows:(1)
Face identification (ID) number, it determines ID of the face in the feature extraction and recognition output.(2)
ID of the inspection operation, it depicts the kind of GD&T such as flatness, parallelism, and perpendicularity.(3)
Tolerance value, it permits designers to establish tolerance limits for the different part characteristics after identifying its function and correlation to mating parts.(4)
Tool employed, it identifies the tools (its type, geometry, etc.), such as contact systems, noncontact techniques, horizontal probe, vertical probe, and star probe, employed for inspection. The choice of the inspection sensor for data acquisition is critical for the entire hybrid process, both in terms of measurement accuracy and time. Certainly, the contact systems, such as the touch-trigger probe is needed when inspection accuracy is favored and time is not a constraint [28]. The scanning touch probe also constitutes a contact inspection system, but it is quite fast and comparatively less accurate than a touch-trigger probe. Similarly, the noncontact systems, namely, laser scanners are preferred when accuracy can be compromised for lesser measurement time. Therefore, this step selects the most appropriate tool to prevent any damage to the part or the probe as well as satisfying the inspection requirements. Moreover, this selection is based on the location, depth, dimensions, measurement cost, orientation, and accessibility of the features.(5)
Datum (reference) faces, they represent the reference points, lines, planes, and axis. The coordinate system and the part origin is generated using these reference features.(6)
Part orientation, it represents the appropriate orientation (which covers a maximum number of features) of the prismatic part, and the setup number is identified in the table.(7)
Number of touch point, it identifies the number of points required to measure the given feature. In this work, stratified and adaptive sampling techniques have been utilized. The stratified sampling is a technique in which the complete set of points is split into subsets (also called strata) [29, 30]. The adaptive sampling demonstrates that regions with higher deformations will receive more points and regions with lower deformation or distortion should receive fewer points [31]. From the FEA, it was found that there were more distortion in the regions near the edges, and they acquired more points than in the middle regions.(8)
Touch point, it provides the coordinate of the probing point. Each boundary of the rectangular face was defined by its boundary corner and using coordinates (x, y, and z) in the extraction module. These coordinates were utilized as input for the algorithm [22].(9)
Geometric inspection boundary–It shows the boundary of the testing face.Table 5
Inspection plan table.
Face IDID ofinspectionoperationTolerancevalueToolusedDatumface IDOrientation of partNo. of touch pointsCoordinates of touch pointsGeometric inspection boundary1GD&TTolerance valueTouch-trigger probeFace IDSetup no.n(Ui, Vi, Wi)—(X1, Y1, Z1)-(Ui+1,Vi+1,Wi+1)(X2, Y2, Z2)-(Ui+2,Vi+2,Vi+2)(X3, Y3, Z3)-.(X4, Y4, Z4)-.(Un,Vn,Vn)


## 3.1. Feature Extraction and Recognition

A rational consistent feature recognition approach has been employed to achieve feature recognition. This system possesses the competence to communicate with hybrid additive, subtractive, and inspection systems. The proposed methodology has been introduced for 3D prismatic parts designed by a solid modeling package (CATIA V5). The system considers a neutral file in STEP format as input and then interprets the information in the file to the manufacturing information. The ISO 10303 STEP Standards AP 203 E2 file has been utilized. The steps of the feature extraction and recognition methodology are as follows:Step 1. Obtain the geometry and topology items for the designed model using STEP file format.Step 2. Obtain topology items in each basic surface and determine its type.Step 3. Verify the feature’s presence in the basic surface depending on loops.Step 4. Determine feature type.Step 5. Determine the detailed features and obtain the related feature geometry parameters.Step 6. Determine the comprehensive machining data for each feature and the designed part.The sample output of a feature extraction file that contains a feature ID, feature name, feature face IDs, and feature dimensions is shown in Table2.Table 2
Feature recognition output [24].
Feature IDFeature nameFeature face IDsDimensions1Pocket throughF1, F2, F3, F4L =, W =, H =2HoleF5, F6H =, R =The rule-based approach as reported in [26] was adopted in this work to recognize the SM features. The following example represents the rule for recognizing the pocket through feature:(1)
For every 4 faces of the surface type plane (named as face 1, face 2, face 3, and face 4) as shown in Figure2(a).(2)
If concave edge count of the outer loop equals 2 for each face.2.1. If each consecutive pair of faces (F1F2, F2F3, F3F4, and F4F1) has one common concave edge of type Line and the pair of faces is perpendicular to each other.2.1.1. Pocket through found.2.1.2. Create a new pocket through object and add to the feature list.(3)
End ForFigure 2
Different faces in features: (a) pocket through; (b) slot blind feature.
(a)(b)For inspection, the data such as faces, edges, and vertices in the extraction file were classified and listed on the basis of their type and angles. Dimensional tolerances (DTs) were calculated from the direct limits of the extracted feature’s faces plus/minus tolerance values. For example, the feature length was calculated as the difference between F3 and F5, etc., as shown in Table3 and Figure 2(b).Table 3
Inspection data extracted for the slot blind feature.
Feature IDFeature nameFeature face IDsBlock face IDsDimensions1Slot blindF1, F2, F3, F4F5, F6Length = F3–F5, width = F1–F2, height = F4–F6


## 3.2. Part Orientation

The orientation rules for the improved surface accuracy depend on the relative orientation of various features/entities (in the extraction file) concerning the build direction. The extracted face orientation can be parallel, perpendicular, or at some angle. To satisfy better Geometric Dimensioning and Tolerancing (GD&T) [24], the objective functions of the optimal orientation of the part were modeled. The objective functions were the maximization of GD&T values and the minimization of production time. The automatic setup generation was achieved by gathering information from the feature extraction and recognition file as shown in Figure 3.Figure 3
Part orientation.In this diagram (Figure3), the setup class is consisted of get these methods, get all feature’s faces, get all part faces, remove interaction faces, generate frequency interaction edges, generate probe feature accessibility, and make feature clustering. The face class has consisted of get face ID, get face name, get face type, get face loop, get face edges, get face vertices, and get face direction.


## 3.3. Automatic Slicing Part to Subpart

The automated hybrid additive, subtractive, and inspection cannot be achieved without the automatic slicing of a part to a subpart. In this study, automatic slicing of a part to a subpart was achieved by dividing the CAD model into different CAD files using the CATVBA editor within CATIA V5. After slicing a part to a subpart, output data for every subpart was obtained in the STEP file format. In this automatic slicing, all data were divided automatically in an empty target file by running the code. The part peaks were the inputs to the slicing code. After that, the subpart area and volume parameters were calculated automatically and data were exported to an Excel datasheet. The exported data were used to calculate the total additive and subtractive area and volume of the part as shown in Figure4.Figure 4
Automatic slicing of a part to a subpart [24].


## 3.4. Simulation

In this work, a finite element analysis (FEA) framework (see Figure5) has also been developed for simulating the additive process. The objective of this analysis was to arrive at the constraints for process planning. The FEA approach in this work was established similar to the scheme and mathematical models reported in [27]. This analysis was subdivided across two phases. The thermal analysis, which solves the heat equations, was carried out to assess the temperature field. The resulting temperature field was then used in the mechanical analysis to determine von Mises stress and component deformation.Figure 5
Schematic diagram of FEA.The additive system employed in the current investigation was equipped with two extruders, which could be heated up to 600°C, and a nozzle of 0.2 mm of diameter was used. The bed temperature could reach up to 100°C. The different nozzle and bed temperatures, part orientations, extruder velocities, and layer thicknesses were applied to study the stress and deformations. The conditions were then identified to produce the most accurate part using additive manufacturing. This analysis also provided the regions (such as areas with high distortion) which could be prioritized for machining using the subtractive process. The FEA results were also utilized to define the sampling strategy required for inspection. For example, the highly deformed regions needed a higher number of points as compared to areas with low distortion.


## 3.5. Process Plan of Additive Operations

AM process-planning steps as shown in Figure6 can be interpreted as the layer-based manufacturing approach, which comprises of the resource requirements including time and travel distance. These aspects comprise the number of layers, shape and size, single or multiple contours in a slice, support material, functionality, accuracy, and surface characteristics [21]. AM process steps are evenly critical and possess a significant effect on the attributes of the produced part.Figure 6
Additive process plan steps.


## 3.6. Process Plan of Subtractive Operations

The subtractive process plan consists of the machining data of the designed part such as sequence of machining, type of operation, machine type, cutting tool, tool approach, and the volume of material removed for each extracted feature as shown in Table4.Table 4
Subtractive process plan.
Feature IDFeature nameFace IDsFeature dimensionFeature volumeFeature locationMachine operationCutting toolMachine toolSetups


## 3.7. Process Plan of Inspection Operations

The inspection planning can be summarized in Table5. The different attributes of the process plan table can be described as follows:(1)
Face identification (ID) number, it determines ID of the face in the feature extraction and recognition output.(2)
ID of the inspection operation, it depicts the kind of GD&T such as flatness, parallelism, and perpendicularity.(3)
Tolerance value, it permits designers to establish tolerance limits for the different part characteristics after identifying its function and correlation to mating parts.(4)
Tool employed, it identifies the tools (its type, geometry, etc.), such as contact systems, noncontact techniques, horizontal probe, vertical probe, and star probe, employed for inspection. The choice of the inspection sensor for data acquisition is critical for the entire hybrid process, both in terms of measurement accuracy and time. Certainly, the contact systems, such as the touch-trigger probe is needed when inspection accuracy is favored and time is not a constraint [28]. The scanning touch probe also constitutes a contact inspection system, but it is quite fast and comparatively less accurate than a touch-trigger probe. Similarly, the noncontact systems, namely, laser scanners are preferred when accuracy can be compromised for lesser measurement time. Therefore, this step selects the most appropriate tool to prevent any damage to the part or the probe as well as satisfying the inspection requirements. Moreover, this selection is based on the location, depth, dimensions, measurement cost, orientation, and accessibility of the features.(5)
Datum (reference) faces, they represent the reference points, lines, planes, and axis. The coordinate system and the part origin is generated using these reference features.(6)
Part orientation, it represents the appropriate orientation (which covers a maximum number of features) of the prismatic part, and the setup number is identified in the table.(7)
Number of touch point, it identifies the number of points required to measure the given feature. In this work, stratified and adaptive sampling techniques have been utilized. The stratified sampling is a technique in which the complete set of points is split into subsets (also called strata) [29, 30]. The adaptive sampling demonstrates that regions with higher deformations will receive more points and regions with lower deformation or distortion should receive fewer points [31]. From the FEA, it was found that there were more distortion in the regions near the edges, and they acquired more points than in the middle regions.(8)
Touch point, it provides the coordinate of the probing point. Each boundary of the rectangular face was defined by its boundary corner and using coordinates (x, y, and z) in the extraction module. These coordinates were utilized as input for the algorithm [22].(9)
Geometric inspection boundary–It shows the boundary of the testing face.Table 5
Inspection plan table.
Face IDID ofinspectionoperationTolerancevalueToolusedDatumface IDOrientation of partNo. of touch pointsCoordinates of touch pointsGeometric inspection boundary1GD&TTolerance valueTouch-trigger probeFace IDSetup no.n(Ui, Vi, Wi)—(X1, Y1, Z1)-(Ui+1,Vi+1,Wi+1)(X2, Y2, Z2)-(Ui+2,Vi+2,Vi+2)(X3, Y3, Z3)-.(X4, Y4, Z4)-.(Un,Vn,Vn)


## 4. Optimization of CAPP

This step aimed to reduce the manufacturing time and hence the fabrication cost of the part produced using HM. To optimize the HM, the mathematical models were employed to obtain the most suitable part orientation as well as minimize the additive and subtractive times. The inspection time was optimized through the minimization of the probe measurement path using GA.


### 4.1. Part Orientation

The optimal orientation can enhance part quality, minimize production time, and lessen support requirements for building the part. An appropriate part orientation for lesser support volume increases the building height and accordingly decreases the fabrication time. It suggests part orientation has to be optimized to minimize support volume as well as manufacturing time. Therefore, the optimal part setup is very critical to an efficient and effective CAPP for HM. There are several parameters such as machine parameters, part geometry, part height, layer thickness, and support volume, which control the part orientation. The primary operations to achieve the optimal orientation are described as follows:Operation 1: obtain the geometry and topology items for the designed object model from STEP file format and extract topology entities in each basic face and determine its type.Operation 2: calculate the GD&T (the first objective value) using the following model:(1)w1j=0θj=0°,90°,180°,−cosθ2otherwise,(2)w2j=0no support,−0.3need  support,(3)w2j=0no  overcure,−0.2need overcure,(4)Wj=1+w1j+w2j+w3jθj=0°,90°,180°,w1j+w2j+w3jθj=0°,90°U90°,180°,−1+w1jw1j≠w2j≠w3j≠0,(5)Qi=∑j=1nNijWj,(6)Nj=nFace∗AjATotal.Operation 3: ifMaxQi=Q¯,Q¯ is a set of equal maximum Qi, the following should be done:(a)
Calculate building time and support volume(b)
Otherwise, go to Step 4Operation 4: selectMaxQiOperation 5: stopwhere Aj represents the total face area for different faces in part orientation i, Nj represents the relative number, which can be calculated by dividing the area of the type (j) divided by the summation of all types, and Wj represents the weight assigned to each type depending on its angle.The goal of this part orientation model was the maximization of horizontal and perpendicular faces in thez-axis (build direction), maximization of base surface area as well as minimization of angular/inclined surfaces, overhanging area, and trapped volume. There can also be parts, where the same value of max Q could be obtained for more than one orientation. Therefore, to select the best orientation, additive time was also be computed for these orientations as shown in Section 4.1.1.


#### 4.1.1. Additive Time

The time of AM is an important factor, especially when different part orientations are available. The build time was calculated for different orientations to select the optimal one. The minimum build time indicates a minimum number of slices and a minimum volumetric error. The volumetric error is the volumetric difference between the CAD model and the AM part. Magic 18 software was used to generate the support structure of the part, and that was exported as an STL file [21]. Subsequently, the area and volume were computed.The build time depends on machine parameters, part geometry, part height, layer thickness, and the amount of support volume needed for building the part. Therefore, the build time was opted as an output variable in each orientation of the part. The total build time was calculated by taking the summation of the time needed to cover the part geometry (Tgeometry), the time needed to build support volume (Tsupport), and the machine time (Tm/c):(7)Total building timeTbuilding=Tgeometry+Tsupport+Tm/c.Tgeometry was calculated by dividing the part volume with machine build rate and was expressed as follows:(8)Tgeometry=Part volumeBuilding  rate of machine.Tsupport was calculated by dividing the amount of support volume (which was provided by the Magic 18 software as shown in Figure 7) with the build rate of the machine:(9)Tsupport=Support volumeBuilding rate  of machine.Figure 7
Support volume generated by Magic 18 software [24].Tm/c was computed by measuring the machine setup time and postprocessing time.


### 4.2. Inspection Time

The reduction of inspection time can be achieved by minimizing the cycle time which primarily depends on the inspection probe path. A GA-based approach was developed to attain the best measurement path for the given part. The GA was employed because it has shown in the past its relevance and benefits in addressing path-planning problems. There are many reasons for its superior performance in comparison with other metaheuristics. The GA is swift, easy to execute, and cost-effective in terms of computational resources. It yields the most stable results and can be used to solve challenging problems in a limited timeframe. It is far better than any other method considering the size of the search area. GA generally performs well because of its parallel search features, where several solutions are merged dynamically until an effective solution is obtained. Since the GA deploys the random-based search algorithms to analyze most of the state space. Therefore, the quality of the acquired solution gets enhanced. Besides, GA’s fewer parameter interactions with just three parameters (population size, crossover probability, and mutation probability) make it simpler, easier to implement, and avoid numerous local optimal solutions in the parameter domain to provide quality solutions. The subsequent section explains the implementation of the GA in this work. The touch points (or measurement points) were specified based on the objective function to ensure the shortest path between two measuring features, as shown in Figure8.Figure 8
Probe path planning based on an objective function.Step 1: random solutions (chromosomes)Initially, the measuring features were grouped based on a generated random solution (chromosome). Subsequently, the touch-point coordinates (input) of the sequenced features were rearranged according to the objective function for obtaining the shortest path between two sequenced features as shown in Figure9.Figure 9
(a) Measured part; (b) measurement path before GA; (c) measurement path after GA.
(a)(b)(c)Step 2: computation of the fitness function for each chromosomeThe fitness function can be defined as the total shortest path between two sequenced features. The objective function can be described as(10)∑k=1PMini=RRMinj=1m=1xR−xj2+yR−yj2+zR−zj2,where i is the first touch point of the first sequenced feature, j is the first touch point of the next sequenced feature, n is the number of touch point coordinates of the last sequenced feature, m is the number of touch point coordinates of the next sequenced feature, and P is the number of sequenced features.Step 3: crossover and mutationSingle-point and two-point crossover operators as shown in Figure10 were used. For the mutation operator, the random swap mutation operator was utilized.Figure 10
Crossover and mutation operators: (a) single-point crossover. (b) Two-point crossover. (c) Swap mutation.
(a)(b)(c)Step 4: stopping criterionThe stopping criterion employed in this study was the maximum number of iterations. This stopping criterion was adopted because it helped to control the processing time as well as eliminate any bias.


## 4.1. Part Orientation

The optimal orientation can enhance part quality, minimize production time, and lessen support requirements for building the part. An appropriate part orientation for lesser support volume increases the building height and accordingly decreases the fabrication time. It suggests part orientation has to be optimized to minimize support volume as well as manufacturing time. Therefore, the optimal part setup is very critical to an efficient and effective CAPP for HM. There are several parameters such as machine parameters, part geometry, part height, layer thickness, and support volume, which control the part orientation. The primary operations to achieve the optimal orientation are described as follows:Operation 1: obtain the geometry and topology items for the designed object model from STEP file format and extract topology entities in each basic face and determine its type.Operation 2: calculate the GD&T (the first objective value) using the following model:(1)w1j=0θj=0°,90°,180°,−cosθ2otherwise,(2)w2j=0no support,−0.3need  support,(3)w2j=0no  overcure,−0.2need overcure,(4)Wj=1+w1j+w2j+w3jθj=0°,90°,180°,w1j+w2j+w3jθj=0°,90°U90°,180°,−1+w1jw1j≠w2j≠w3j≠0,(5)Qi=∑j=1nNijWj,(6)Nj=nFace∗AjATotal.Operation 3: ifMaxQi=Q¯,Q¯ is a set of equal maximum Qi, the following should be done:(a)
Calculate building time and support volume(b)
Otherwise, go to Step 4Operation 4: selectMaxQiOperation 5: stopwhere Aj represents the total face area for different faces in part orientation i, Nj represents the relative number, which can be calculated by dividing the area of the type (j) divided by the summation of all types, and Wj represents the weight assigned to each type depending on its angle.The goal of this part orientation model was the maximization of horizontal and perpendicular faces in thez-axis (build direction), maximization of base surface area as well as minimization of angular/inclined surfaces, overhanging area, and trapped volume. There can also be parts, where the same value of max Q could be obtained for more than one orientation. Therefore, to select the best orientation, additive time was also be computed for these orientations as shown in Section 4.1.1.


### 4.1.1. Additive Time

The time of AM is an important factor, especially when different part orientations are available. The build time was calculated for different orientations to select the optimal one. The minimum build time indicates a minimum number of slices and a minimum volumetric error. The volumetric error is the volumetric difference between the CAD model and the AM part. Magic 18 software was used to generate the support structure of the part, and that was exported as an STL file [21]. Subsequently, the area and volume were computed.The build time depends on machine parameters, part geometry, part height, layer thickness, and the amount of support volume needed for building the part. Therefore, the build time was opted as an output variable in each orientation of the part. The total build time was calculated by taking the summation of the time needed to cover the part geometry (Tgeometry), the time needed to build support volume (Tsupport), and the machine time (Tm/c):(7)Total building timeTbuilding=Tgeometry+Tsupport+Tm/c.Tgeometry was calculated by dividing the part volume with machine build rate and was expressed as follows:(8)Tgeometry=Part volumeBuilding  rate of machine.Tsupport was calculated by dividing the amount of support volume (which was provided by the Magic 18 software as shown in Figure 7) with the build rate of the machine:(9)Tsupport=Support volumeBuilding rate  of machine.Figure 7
Support volume generated by Magic 18 software [24].Tm/c was computed by measuring the machine setup time and postprocessing time.


## 4.1.1. Additive Time

The time of AM is an important factor, especially when different part orientations are available. The build time was calculated for different orientations to select the optimal one. The minimum build time indicates a minimum number of slices and a minimum volumetric error. The volumetric error is the volumetric difference between the CAD model and the AM part. Magic 18 software was used to generate the support structure of the part, and that was exported as an STL file [21]. Subsequently, the area and volume were computed.The build time depends on machine parameters, part geometry, part height, layer thickness, and the amount of support volume needed for building the part. Therefore, the build time was opted as an output variable in each orientation of the part. The total build time was calculated by taking the summation of the time needed to cover the part geometry (Tgeometry), the time needed to build support volume (Tsupport), and the machine time (Tm/c):(7)Total building timeTbuilding=Tgeometry+Tsupport+Tm/c.Tgeometry was calculated by dividing the part volume with machine build rate and was expressed as follows:(8)Tgeometry=Part volumeBuilding  rate of machine.Tsupport was calculated by dividing the amount of support volume (which was provided by the Magic 18 software as shown in Figure 7) with the build rate of the machine:(9)Tsupport=Support volumeBuilding rate  of machine.Figure 7
Support volume generated by Magic 18 software [24].Tm/c was computed by measuring the machine setup time and postprocessing time.


## 4.2. Inspection Time

The reduction of inspection time can be achieved by minimizing the cycle time which primarily depends on the inspection probe path. A GA-based approach was developed to attain the best measurement path for the given part. The GA was employed because it has shown in the past its relevance and benefits in addressing path-planning problems. There are many reasons for its superior performance in comparison with other metaheuristics. The GA is swift, easy to execute, and cost-effective in terms of computational resources. It yields the most stable results and can be used to solve challenging problems in a limited timeframe. It is far better than any other method considering the size of the search area. GA generally performs well because of its parallel search features, where several solutions are merged dynamically until an effective solution is obtained. Since the GA deploys the random-based search algorithms to analyze most of the state space. Therefore, the quality of the acquired solution gets enhanced. Besides, GA’s fewer parameter interactions with just three parameters (population size, crossover probability, and mutation probability) make it simpler, easier to implement, and avoid numerous local optimal solutions in the parameter domain to provide quality solutions. The subsequent section explains the implementation of the GA in this work. The touch points (or measurement points) were specified based on the objective function to ensure the shortest path between two measuring features, as shown in Figure8.Figure 8
Probe path planning based on an objective function.Step 1: random solutions (chromosomes)Initially, the measuring features were grouped based on a generated random solution (chromosome). Subsequently, the touch-point coordinates (input) of the sequenced features were rearranged according to the objective function for obtaining the shortest path between two sequenced features as shown in Figure9.Figure 9
(a) Measured part; (b) measurement path before GA; (c) measurement path after GA.
(a)(b)(c)Step 2: computation of the fitness function for each chromosomeThe fitness function can be defined as the total shortest path between two sequenced features. The objective function can be described as(10)∑k=1PMini=RRMinj=1m=1xR−xj2+yR−yj2+zR−zj2,where i is the first touch point of the first sequenced feature, j is the first touch point of the next sequenced feature, n is the number of touch point coordinates of the last sequenced feature, m is the number of touch point coordinates of the next sequenced feature, and P is the number of sequenced features.Step 3: crossover and mutationSingle-point and two-point crossover operators as shown in Figure10 were used. For the mutation operator, the random swap mutation operator was utilized.Figure 10
Crossover and mutation operators: (a) single-point crossover. (b) Two-point crossover. (c) Swap mutation.
(a)(b)(c)Step 4: stopping criterionThe stopping criterion employed in this study was the maximum number of iterations. This stopping criterion was adopted because it helped to control the processing time as well as eliminate any bias.


## 5. Case Study

The proposed method can successfully be applied to any prismatic part. The established approach had in fact been implemented for many parts with prismatic features to validate its performance. However, to limit the length of the manuscript, the paper contains only one case study as shown in Figure11 to depict the working and the findings of the proposed method. The 3D solid model was designed in CATIA V5, and the GD&T data were obtained from the STEP AP 203 file utilizing an object-oriented technique. Seventy-two faces were obtained from the geometric information file, including their loops, edges, vertices, face direction, face type, and face area. Twelve manufacturing features were obtained and recognized, inclusive of their dimensions. The feature ID, feature name, feature face IDs, and feature dimensions are indexed in the output extraction table (Table 6). Table 7 represents the subtractive process plan obtained after feature extraction and recognition.Figure 11
Test part (solid body) with feature IDs.Table 6
Feature extraction and recognition.
Feature IDFeature nameFace IDsFeature dimensionFeature peak1Slot through2, 1, 3Height: 40, width: 30, length: 300–402Slot through5, 6, 4Height: 40, width: 30, length: 300–403STEP blind round corner9, 7, 10, 7Length: 30, width: 25, height: 10, radius: 530–404STEP blind round corner13, 12, 14, 12Length: 30, width: 25, height: 10, radius: 530–405STEP blind round corner16, 17, 18, 17Length: 30, width: 25, height: 10, radius: 530–406STEP blind round corner20, 19, 22, 19Length: 30, width: 25, height: 10, radius: 530–407Slot blind round corner27, 23, 24, 26, 24, 25Length: 35, width: 30, height: 15, radius: 525–408Pocket blind round corner35, 37, 34, 33, 30, 31, 32, 29, 36Length: 40, width: 30, height: 10, radius: 530–409Slot blind round corner42, 40, 43, 41, 43, 38Length: 35, width: 30, height: 15, radius: 525–4010Slot blind round corner48, 44, 45, 47, 45, 46Length: 35, width: 80, height: 12, radius: 50–1211Slot blind round corner54, 52, 51, 53, 51, 50Length: 35, width: 80, height: 12, radius: 50–1212Pocket blind round corner62, 64, 61, 60, 57, 58, 59, 56, 63Length: 40, width: 30, height: 10, radius: 50–1013Slot through67, 68, 66Length: 110 width: 10 height: 634–4014Slot through70, 69, 71Length: 110, width: 10, height: 634–4015Slot through76, 77, 75Length: 40, width: 5, height: 100–4016Slot through79, 80, 78Length: 40 width: 5 height: 100–4017Slot through86, 85, 87Length: 40, width: 5, height: 100–4018Slot through89, 88, 90Length: 40, width: 5, height: 100–4019Slot blind94, 95, 93, 96Length: 20, width: 10, height: 1510–2020Slot blind98, 99, 97, 100Length: 15, width: 10, height: 2010–2021Slot blind104, 105, 103, 106Length: 20, width: 10, height: 1510–2022Slot blind108, 109, 107, 110Length: 20, width: 10, height: 1510–2023Hole blind112, 113, 114Radius: 5, height: 1515–2524Hole blind116, 115, 117Radius: 5, height: 1515–2525Hole blind120, 119, 121Radius: 5, height: 1515–2526Hole blind122, 123, 124Radius: 5, height: 1515–25Table 7
Subtractive process plan.
Feature IDFeature nameFeature volumeMachine operationCutting toolMachine toolWork piece materialBrinell hardnessCutting tool materialMachine parameterFeature locationSetups7Slot blind round corner15709Slot millingCorner rounding milling cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, 1, 0), location ZX_FRONT; TAD (0, 0, −1), location XY_TOP1 (XY_TOP)9Slot blind round corner15709Slot millingCorner rounding milling cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, −1, 0), location ZX_BACK; TAD (0, 0, −1), location XY_TOP1 (XY_TOP)8Pocket blind round corner11946Pocket millingCorner rounding milling cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, 0, −1), location XY_TOP1 (XY_TOP)3STEP blind round corner7486Shoulder millingCorner rounding milling cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (1, 0, 0), location YZ_LEFT; TAD (0, 1, 0), location ZX_FRONT; TAD (0, 0, −1), location XY_TOP1 (XY_TOP)4STEP blind round corner7486Shoulder millingCorner rounding milling cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (−1, 0, 0), location YZ_RIGHT; TAD (0, 1, 0), location ZX_FRONT TAD (0, 0, −1), location XY_TOP1 (XY_TOP)5STEP blind round corner7486Shoulder millingCorner rounding milling cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (−1, 0, 0), location YZ_RIGHT; TAD (0, −1, 0), location ZX_BACK; TAD (0, 0, −1), location XY_TOP1 (XY_TOP)6STEP blind round corner7486Shoulder millingCorner rounding milling cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (1, 0, 0), location YZ_LEFT; TAD (0, −1, 0), location ZX_BACK; TAD (0, 0, −1), location XY_TOP1 (XY_TOP)13Slot through6600Slot millingEnd mill cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, 0, −1), location XY_TOP1 (XY_TOP)14Slot through6600Slot millingEnd mill cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, 0, −1), location XY_TOP1 (XY_TOP)23Hole blind1177.5DrillingTwist drillMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, −1, 0), location ZX_BACK2 (ZX_BACK)24Hole blind1177.5DrillingTwist drillMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, −1, 0), location ZX_BACK2 (ZX_BACK)25Hole blind1177.5DrillingTwist drillMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, 1, 0), location ZX_FRONT3 (ZX_FRONT)26Hole blind1177.5DrillingTwist drillMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, 1, 0), location ZX_FRONT3 (ZX_FRONT)2Slot through36000Slot millingEnd mill cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (−1, 0, 0), location YZ_RIGHT4 (YZ_RIGHT)15Slot through2000Slot millingEnd mill cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (−1, 0, 0), location YZ_RIGHT4 (YZ_RIGHT)16Slot through2000Slot millingEnd mill cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (−1, 0, 0), location YZ_RIGHT4 (YZ_RIGHT)1Slot through36000Slot millingEnd mill cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (1, 0, 0), location YZ_LEFT5 (YZ_LEFT)17Slot through2000Slot millingEnd mill cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (1, 0, 0), location YZ_LEFT5 (YZ_LEFT)18Slot through2000Slot millingEnd mill cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (1, 0, 0), location YZ_LEFT5 (YZ_LEFT)12Pocket blind round corner11946Pocket millingCorner rounding milling cutterMilling machineLow-carbon steel150–200HSS21–30 (m/min), 0.3 (mm/rev)TAD (0, 0, 1), location XY_BOTTOM6 (XY_BOTTOM)


### 5.1. Part Orientation

A comparison between the different orientations was examined depending on the GD&T and building time needed for the entire part [24]. Figure 12 represents the hybrid setup consisting of additive, subtractive, and inspection centers.Figure 12
Hybrid additive, subtractive, and inspection machine setup.There were four possible orientations of the part. These were (+x, −x) orientations and (+y, −y) orientations as shown in Figure 13.Figure 13
(a) Orientation no. 1; (b) orientation no. 2; (c) orientation no. 3; (d) orientation no. 4.
(a)(b)(c)(d)An overhanging surface is a type of prismatic geometric feature that warps easily if the inclination angle is smaller. It results in a staircase volume, which causes a dimensional error and poor surface roughness. The resulting data as shown in Table8 represent the total face area (Aj) for different faces in part orientation (i). To determine the best orientation, the relative number (Nj), which is the area of type (j) divided by the summation of all types, was calculated. Next, the weight (Wj) was assigned to each type depending on its angle. Next, Qi was calculated for each type by multiplying the relative number (Nj) of the face (j) with its weight. Finally, the Qi for each type was added up to get the final rating for the given orientation. This final rating for each orientation was then compared with the ratings obtained for other orientation to select the best orientation of the part.Table 8
GD&T objective function for orientations 1, 2, 3, and 4.
Type[1][2][3][4][5][6][7][8][9][10][11][12][13]Orient. 1Qi14.1249.7018.870−24.7800000−0.77−1.480Sum (Qi)55.655Orient. 2Qi40.7245.966.590−9.3700002.070.5900Sum (Qi)86.56Orient. 3Qi20.7238.806.5920−22.2400000.471.0300Sum (Qi)45.38Orient. 4Qi23.1344.243.760−22.9100000−0.207−1.110Sum (Qi)46.915It is obvious from the above table that orientation 2 had the maximumQi value, based on the output of the GD&T objective function.


### 5.2. Additive Time

The constant factors used in the model are shown in Table9.Table 9
The constant factors.
FactorValueUnitsBuilding rate30cm3/hBatch size2PartsScan speed25m/sDensity0.9g/cm3Material price60$/kgThe additive time for different orientation was also computed to further examine the best orientation. After computing the additive time, orientation 2 was found to have the least building time as shown in Table10.Table 10
Additive process plan.
Orient. 1Orient. 2Orient. 3Orient. 4Part building time (h)20 hr 45 min15 hr 25 min15 hr 59 min18 hr 35 minPart volume (mm3)281160281160281160281160Support volume (mm3)44020302603389033940Preparation time (h)0.50.50.50.5Postprocessing time (0.10 h/kg)41 min26 min28 min37 minTotal building time (h)22 hr 6 min16 hr 21 min16 hr 57 min19 hr 42 min


### 5.3. Requirements of Process Planning

The mandatory requirements of the process planning were obtained through FEA discussed in Section3.5. The FEA was performed using the commercial software ABAQUS (Abaqus v 6.14, Dassault Systems). The Visual Studio (v 2012) and compiler (Intel Parallel Studio XE 2013) were employed to connect the subroutine with the Abaqus. The material properties were adapted to conform to a commercial acrylonitrile butadiene styrene (ABS) filament as shown in Table 11.Table 11
Mechanical and thermal properties of the ABS at the ambient temperature.
PropertyValueDensity1040 kg/m3Young’s modulus1826 MPaPoisson’s ratio0.32Yield stress25 MPaThermal expansion coefficient9 × 10−5°C−1Specific heat1290 J/(°C. kg)Thermal conductivity0.15 W/(°C. m)Glass transition temperature95°CFigures14(a) and 14(b) demonstrate the orientation of the part with regard to the building platform and the layers of the sliced model, respectively. The meshing strategy consisted of Hex-dominated elements with a minimum size of 0.2 mm. This simulation involved approximately 92000 elements in total. The simulation was performed for five layers due to limited computational resources.Figure 14
(a) Layer orientation; (b) meshing strategy.
(a)(b)The simulation was performed in several conditions, involving different part orientations as well as temperatures. The lowest von Mises stress and deformation as shown in Figures15(a) and 15(b), respectively, were obtained under following conditions: part setup, orientation 2; layer thickness, 0.2 mm; extruder temperature, 240°C; bed temperature, 80°C; ambient temperature, 25°C; extruder velocity, 20 mm/s; and nozzle diameter, 0.2 mm. The maximum value of von Mises stress was 0.02875 MPa, and maximum deformation was 1.6 × 10−5 mm. These conditions also defined the constraints for the additive process. The additive process had to be executed in these conditions to get the best results and minimize the usage of subtractive machining. It was also found that the deformation was maximum at the edges; that is, it increased from the center of the part to its boundary. As a consequence, the sample size was optimized, where more points were assigned in regions near the edges.Figure 15
(a) Stress distribution; (b) resultant total deformation.
(a)(b)


### 5.4. Inspection Time

For the present case study, the maximum number of iterations of 2000 yielded the best result. With this maximum number of iterations, the length of the obtained probe path was 617.23 mm. As shown in Table12, the probe’s travel distance between two features decreased with an increasing number of iterations. However, it stabilized after 2000 generations, and no additional improvement in probe distance was observed with an increase in the number of iterations. Therefore, a maximum number of iterations of 2000 was selected as the stopping criterion of the proposed GA.Table 12
Application of the GA to determine the inspection path for a subpart.
No. of generationProbe traveling path (mm)Genetic algorithm sequencesCMM running time (minutes)CPU time (second)Fitness graphBefore GAAfter GABefore GAAfter GA501007.3745.425-6-3-2-8-9-1-10-7-46 min 32 s4 min 19 s131.13100674.747-4-5-6-3-2-8-9-1-10262.09200617.235-6-3-7-10-1-9-8-2-4521.88500617.234-2-8-9-1-10-7-3-6-5863.141000617.235-6-3-7-10-1-9-8-2-41706.72000617.235-6-3-7-10-1-9-8-2-43641.7The best sequence obtained from the GA was 5-6-3-7-10-1-9-8-2-4, which corresponded to the minimum probe travel distance of 617.23 mm among all features. Initially, when the measurement path was selected randomly, the traveling distance was 1007.3 mm and the coordinate measuring machine (CMM) running time was 6 min 32 s. However, after the application of GA, the travel path of the probe decreased by 38.72%, from 1007.3 to 617.23 mm, while the CMM operation time was significantly reduced from 6 min 32 s to 4 min 19 s. A time reduction of almost 33.93% was observed during the actual measurement of the part. Although the CPU time also increased with an increasing number of iterations, this increase was negligible as compared with the reduction in the CMM operation time.From the previous data, the best setup was the orientation no. 2, and depending on the peaks of all features, the whole part was divided to seven subparts within values as shown in Table13.Table 13
All possible peaks and subparts.
Sequence operationSubpartSubpart peak (or subpart peak height)Additive processesSubtractive processesMachine typeSubpart volume (ABS filament volume) (mm3)Building time (hr)Subpart and subtractive feature IDsSubtractive operationSubtractive area (mm2)Machine type10–10FDM machine66429.203.52ID 1, ID 2, ID 10, ID 11, ID 12, ID 15, ID 16, ID 17, ID 18Surface finish9000Milling machine210–12FDM machine13242.921.02ID 1, ID 2, ID 10, ID 11, ID 15, ID 16, ID 17, ID 18, ID 19, ID 20, ID 21, ID 22Surface finish1000Milling machine312–20FDM machine95243.814.16ID 1, ID 2, ID 15, ID 16, ID 17, ID 18, ID 19, ID 20, ID 21, ID 22, ID 23, ID 24, ID 25, ID 26Surface finish6000Milling machine420–25FDM machine64643.813.15ID 1, ID 2, ID 15, ID 16, ID 17, ID 18, ID 23, ID 24, ID 25, ID 26Surface finish5000Milling machine525–30FDM machine56607.302.45ID 1, ID 2, ID 7, ID 9, ID 15, ID 16, ID 17, ID 18Surface finish5000Milling machine630–34FDM machine28657.521.51ID 1, ID 2, ID 3, ID 4, ID 5, ID 6, ID 7, ID 8, ID 9, ID 15, ID 16, ID 17, ID 18Surface finish4000Milling machine734–40FDM machine29786.282.37ID 1, ID 2, ID 3, ID 4, ID 5, ID 6, ID 7, ID 8, ID 9, ID 13, ID 14, ID 15, ID 16, ID 17, ID 18Surface finish8000Milling machineThe dynamic process plan sequence acquired for the test part can be seen in Figure16. As depicted in Figure 16, the entire part was first divided into subparts based on different peaks. Each subpart was fabricated in an additive fashion. Then, the subtractive operation was employed to obtain accurate features depending on the quality requirement. The inspection was also carried out simultaneously to decide if the subpart needed a subtractive operation or not.Figure 16
Steps of manufacturing part.


## 5.1. Part Orientation

A comparison between the different orientations was examined depending on the GD&T and building time needed for the entire part [24]. Figure 12 represents the hybrid setup consisting of additive, subtractive, and inspection centers.Figure 12
Hybrid additive, subtractive, and inspection machine setup.There were four possible orientations of the part. These were (+x, −x) orientations and (+y, −y) orientations as shown in Figure 13.Figure 13
(a) Orientation no. 1; (b) orientation no. 2; (c) orientation no. 3; (d) orientation no. 4.
(a)(b)(c)(d)An overhanging surface is a type of prismatic geometric feature that warps easily if the inclination angle is smaller. It results in a staircase volume, which causes a dimensional error and poor surface roughness. The resulting data as shown in Table8 represent the total face area (Aj) for different faces in part orientation (i). To determine the best orientation, the relative number (Nj), which is the area of type (j) divided by the summation of all types, was calculated. Next, the weight (Wj) was assigned to each type depending on its angle. Next, Qi was calculated for each type by multiplying the relative number (Nj) of the face (j) with its weight. Finally, the Qi for each type was added up to get the final rating for the given orientation. This final rating for each orientation was then compared with the ratings obtained for other orientation to select the best orientation of the part.Table 8
GD&T objective function for orientations 1, 2, 3, and 4.
Type[1][2][3][4][5][6][7][8][9][10][11][12][13]Orient. 1Qi14.1249.7018.870−24.7800000−0.77−1.480Sum (Qi)55.655Orient. 2Qi40.7245.966.590−9.3700002.070.5900Sum (Qi)86.56Orient. 3Qi20.7238.806.5920−22.2400000.471.0300Sum (Qi)45.38Orient. 4Qi23.1344.243.760−22.9100000−0.207−1.110Sum (Qi)46.915It is obvious from the above table that orientation 2 had the maximumQi value, based on the output of the GD&T objective function.


## 5.2. Additive Time

The constant factors used in the model are shown in Table9.Table 9
The constant factors.
FactorValueUnitsBuilding rate30cm3/hBatch size2PartsScan speed25m/sDensity0.9g/cm3Material price60$/kgThe additive time for different orientation was also computed to further examine the best orientation. After computing the additive time, orientation 2 was found to have the least building time as shown in Table10.Table 10
Additive process plan.
Orient. 1Orient. 2Orient. 3Orient. 4Part building time (h)20 hr 45 min15 hr 25 min15 hr 59 min18 hr 35 minPart volume (mm3)281160281160281160281160Support volume (mm3)44020302603389033940Preparation time (h)0.50.50.50.5Postprocessing time (0.10 h/kg)41 min26 min28 min37 minTotal building time (h)22 hr 6 min16 hr 21 min16 hr 57 min19 hr 42 min


## 5.3. Requirements of Process Planning

The mandatory requirements of the process planning were obtained through FEA discussed in Section3.5. The FEA was performed using the commercial software ABAQUS (Abaqus v 6.14, Dassault Systems). The Visual Studio (v 2012) and compiler (Intel Parallel Studio XE 2013) were employed to connect the subroutine with the Abaqus. The material properties were adapted to conform to a commercial acrylonitrile butadiene styrene (ABS) filament as shown in Table 11.Table 11
Mechanical and thermal properties of the ABS at the ambient temperature.
PropertyValueDensity1040 kg/m3Young’s modulus1826 MPaPoisson’s ratio0.32Yield stress25 MPaThermal expansion coefficient9 × 10−5°C−1Specific heat1290 J/(°C. kg)Thermal conductivity0.15 W/(°C. m)Glass transition temperature95°CFigures14(a) and 14(b) demonstrate the orientation of the part with regard to the building platform and the layers of the sliced model, respectively. The meshing strategy consisted of Hex-dominated elements with a minimum size of 0.2 mm. This simulation involved approximately 92000 elements in total. The simulation was performed for five layers due to limited computational resources.Figure 14
(a) Layer orientation; (b) meshing strategy.
(a)(b)The simulation was performed in several conditions, involving different part orientations as well as temperatures. The lowest von Mises stress and deformation as shown in Figures15(a) and 15(b), respectively, were obtained under following conditions: part setup, orientation 2; layer thickness, 0.2 mm; extruder temperature, 240°C; bed temperature, 80°C; ambient temperature, 25°C; extruder velocity, 20 mm/s; and nozzle diameter, 0.2 mm. The maximum value of von Mises stress was 0.02875 MPa, and maximum deformation was 1.6 × 10−5 mm. These conditions also defined the constraints for the additive process. The additive process had to be executed in these conditions to get the best results and minimize the usage of subtractive machining. It was also found that the deformation was maximum at the edges; that is, it increased from the center of the part to its boundary. As a consequence, the sample size was optimized, where more points were assigned in regions near the edges.Figure 15
(a) Stress distribution; (b) resultant total deformation.
(a)(b)


## 5.4. Inspection Time

For the present case study, the maximum number of iterations of 2000 yielded the best result. With this maximum number of iterations, the length of the obtained probe path was 617.23 mm. As shown in Table12, the probe’s travel distance between two features decreased with an increasing number of iterations. However, it stabilized after 2000 generations, and no additional improvement in probe distance was observed with an increase in the number of iterations. Therefore, a maximum number of iterations of 2000 was selected as the stopping criterion of the proposed GA.Table 12
Application of the GA to determine the inspection path for a subpart.
No. of generationProbe traveling path (mm)Genetic algorithm sequencesCMM running time (minutes)CPU time (second)Fitness graphBefore GAAfter GABefore GAAfter GA501007.3745.425-6-3-2-8-9-1-10-7-46 min 32 s4 min 19 s131.13100674.747-4-5-6-3-2-8-9-1-10262.09200617.235-6-3-7-10-1-9-8-2-4521.88500617.234-2-8-9-1-10-7-3-6-5863.141000617.235-6-3-7-10-1-9-8-2-41706.72000617.235-6-3-7-10-1-9-8-2-43641.7The best sequence obtained from the GA was 5-6-3-7-10-1-9-8-2-4, which corresponded to the minimum probe travel distance of 617.23 mm among all features. Initially, when the measurement path was selected randomly, the traveling distance was 1007.3 mm and the coordinate measuring machine (CMM) running time was 6 min 32 s. However, after the application of GA, the travel path of the probe decreased by 38.72%, from 1007.3 to 617.23 mm, while the CMM operation time was significantly reduced from 6 min 32 s to 4 min 19 s. A time reduction of almost 33.93% was observed during the actual measurement of the part. Although the CPU time also increased with an increasing number of iterations, this increase was negligible as compared with the reduction in the CMM operation time.From the previous data, the best setup was the orientation no. 2, and depending on the peaks of all features, the whole part was divided to seven subparts within values as shown in Table13.Table 13
All possible peaks and subparts.
Sequence operationSubpartSubpart peak (or subpart peak height)Additive processesSubtractive processesMachine typeSubpart volume (ABS filament volume) (mm3)Building time (hr)Subpart and subtractive feature IDsSubtractive operationSubtractive area (mm2)Machine type10–10FDM machine66429.203.52ID 1, ID 2, ID 10, ID 11, ID 12, ID 15, ID 16, ID 17, ID 18Surface finish9000Milling machine210–12FDM machine13242.921.02ID 1, ID 2, ID 10, ID 11, ID 15, ID 16, ID 17, ID 18, ID 19, ID 20, ID 21, ID 22Surface finish1000Milling machine312–20FDM machine95243.814.16ID 1, ID 2, ID 15, ID 16, ID 17, ID 18, ID 19, ID 20, ID 21, ID 22, ID 23, ID 24, ID 25, ID 26Surface finish6000Milling machine420–25FDM machine64643.813.15ID 1, ID 2, ID 15, ID 16, ID 17, ID 18, ID 23, ID 24, ID 25, ID 26Surface finish5000Milling machine525–30FDM machine56607.302.45ID 1, ID 2, ID 7, ID 9, ID 15, ID 16, ID 17, ID 18Surface finish5000Milling machine630–34FDM machine28657.521.51ID 1, ID 2, ID 3, ID 4, ID 5, ID 6, ID 7, ID 8, ID 9, ID 15, ID 16, ID 17, ID 18Surface finish4000Milling machine734–40FDM machine29786.282.37ID 1, ID 2, ID 3, ID 4, ID 5, ID 6, ID 7, ID 8, ID 9, ID 13, ID 14, ID 15, ID 16, ID 17, ID 18Surface finish8000Milling machineThe dynamic process plan sequence acquired for the test part can be seen in Figure16. As depicted in Figure 16, the entire part was first divided into subparts based on different peaks. Each subpart was fabricated in an additive fashion. Then, the subtractive operation was employed to obtain accurate features depending on the quality requirement. The inspection was also carried out simultaneously to decide if the subpart needed a subtractive operation or not.Figure 16
Steps of manufacturing part.


## 6. Conclusion

An efficacious hybridization of additive, subtractive, and inspection processes can be accomplished through a vital link known as CAPP. Indeed, the most relevant and invaluable CAPP can be achieved through the optimization of its various components. Henceforth, in this work, different approaches were employed to obtain the most appropriate CAPP for the hybrid process, which could manufacture parts efficiently and effectively. Consequently, mathematical models were utilized to optimize part orientation and minimize additive and subtractive times. The GA was also used to achieve the shortest measurement time. An implementation of the introduced methodology has also been presented by using a test part case study for demonstrating its performance. Certainly, the optimized CAPP resulted in an adequate and the most competent hybrid process.This paper has indeed developed a significant understanding and knowledge regarding the CAPP of the hybrid process. It would facilitate the users to acquire a most convenient strategy, to enable accurate and efficient hybrid fabrication of a given component. The primary contribution of this work pertains to the development of an optimized CAPP for a combined additive, subtractive, and inspection processes, using mathematical models, rules, and metaheuristics. The generation of CAPP was initiated depending on the information of additive and subtractive processes, available resources, technological constraints, etc. The presented approach exhibited several benefits, including reduced fabrication time, minimal resource utilization, and lesser human interference. This approach also had a downside in addition to several benefits that will be explored in future research. This limitation was mainly correlated with the effective specification of the dimensional tolerances. This could be because the STEP file used to extract and recognize features did not provide adequate information about the number of points, their distribution, the evaluation algorithm, etc. Certainly, this problem becomes more complex with the presence of complicated features, such as cylinders and curved surfaces. The specification of dimensional tolerances is indeed a key obstacle in the development of an inspection process plan as part of the hybrid process. Furthermore, as this approach can only be used for prismatic parts, therefore the future work also aims to upgrade its application for cylindrical, round, and free-form components. A comprehensive FEA scheme involving both additive as well as subtractive manufacturing processes will also be established to further enhance the quality of process planning.



---

*Source: 7458239-2020-07-18.xml*