Concurrent optimization of process parameters and product design variables for near net shape manufacturing processes

This paper presents a new systematic approach to the optimization of both design and manufacturing variables across a multi-step production process. The approach assumes a generic manufacturing process in which an initial near net shape (NNS) process is followed by a limited number of finishing operations. In this context the optimisation problem becomes a multi-variable problem in which the aim is to optimize by minimizing cost (or time) and improving technological performances (e.g. turning force). To enable such computation a methodology, named conditional design optimization (CoDeO) is proposed which allows the modelling and simultaneous optimization of process parameters and product design (geometric variables), using single or multi-criteria optimization strategies. After investigation of CoDeO’s requirements, evolutionary algorithms, in particular Genetic Algorithms, are identified as the most suitable for overall NNS manufacturing chain optimization The CoDeO methodology is tested using an industrial case study that details a process chain composed of casting and machining processes. For the specific case study presented the optimized process resulted in cost savings of 22% (corresponding to equivalent machining time savings) and a 10% component weight reduction.


Introduction
The global desire to reduce energy, material consumption and emissions continuously push companies to adopt manufacturing technologies that reduce wastages and maintain productivity. In response to these drivers near net shape (NNS) manufacturing processes have emerged as a powerful tool for achieving such savings.
NNS manufacturing is a relative, rather than absolute, propriety of a process chain that minimizes raw material utilization and finishing machining operations (Altan and Allen 1990) in comparison with the alternatives (Marini et al. 2017a).
Optimization of such processes must take a broad view and control both process parameters and design variables in order to minimize the raw material utilization and the time of machining operations. The determination of optimal, or acceptable near-optimal, solutions by a suitable optimization technique has been recognized as critical but difficult task for researchers and practitioners (Mukherjee and Ray 2006).
However, in industrial enterprises and academic literature, product design and manufacturing process optimization are linked by predominately qualitative " Design-for-X" methodologies rather than the explicit setting of dimensions and process variables. In contrast, this paper presents a "offline" methodology for selecting the best configuration of design and process parameters in a NNS supply chain, which integrates both product and process design.

Define an optimization methodology 2. Quantify benefits by implement it on a case study
The proposed methodology, known as conditional design optimization (CoDeO), is limited to a NNS manufacturing chain: where, a NNS chain is composed of a primary shaping process and a limited number of machining steps.

Methodology
In other contexts the idea of concurrent optimization between process parameters, manufacturing systems and design have been explored using different methodologies.
In mechatronic, quantitative approaches have been adopted for developing "Design for Control" methodologies for concurrent optimization. Authors, such as Bi and Zhang (2001), Li et al. (2001), Rampersad (1995) and Ouyang et al. (2004), develop concurrent mathematical models able to control robot mechanics and concurrently optimize robot dynamics and product design. In comparison to these approaches, an "online" (i.e. real-time) process and product optimization of flexible systems is beyond the scope of this paper.
Elsewhere concurrent optimization methods have been extended to applications in operations management and management science. Qualitative and quantitative methodologies, such as Zhang et al. (2019) and Ha and Porteus (1995), have been developed for modelling the links between manufacturing and design operations, allowing their optimal scheduling and design. However, such a systemic and holistic managerial view of the enterprise systems (Zhang and Wang 2016;Giachetti 1998) is beyond of the scope of this work.
In the context of NNS manufacturing, the machining steps gives an opportunity to formulate a generic approach to optimize NNS processes in accordance with production requirements and processes dynamics. In the previous investigations, NNS processes have been investigated on case-bycase basis, for example as in Gupta (2019).
Product design variation is limited by design constraints and working conditions, while process parameters are limited by technological feasibility and consequently strictly process dependant. Understanding their correlations across the process chain is fundamental for a holistic optimization of both process and product.

Roadmap
The next section reviews the literature on both numerical and design optimization. After the literature survey, the problem formulation is defined (section "Concurrent optimization of near net shape processes") and the proposed methodology (illustrated by a validation case study) is presented. The case study (section "Case study") requires optimisation of a blank (produced by centrifugal casting) and finishing machining steps (done by turning). To optimize the turning process, the literature is also reviewed to identify viable optimization's algorithm settings. Section "Case study results" presents the case study results, which suggest an overall cost saving of 22%. Finally, section "Discussion" discusses the limitations of the CoDeO methodology before section "Conclusions" draws some conclusions.

Literature survey: process and product design optimization for manufacturing processes
Process optimization is a trade-off between maximizing the time performances and the process by considering the process dynamics. For this reason, the optimization model needs to take into consideration the process mechanics. In the literature, researchers introduce process models (e.g. in turning process: turning power, cutting force, obtainable roughness equations) into optimization algorithm in different ways: using these multi-variable equations (that often take the form of non-linear functions) as an objective function itself (i.e. having a multi-criteria optimization and therefore a Pareto solution), as non-linear constraints or in a combination of the previous.
The modelling of NNS processes depends on the optimization algorithm selected. In a multi-variable and multicriteria optimization, that NNS methods require, heuristics can be generally defined as search techniques that are able to provide, among several alternatives, the most effective solution for achieving a certain goal. The search based algorithm provides a rule, or set of rules, for identifying an acceptable solution, or solutions, at a certain (or acceptable) computational cost. Heuristic-based search techniques find viable application cases whenever conventional techniques fail to achieve results or their computational cost is too high (e.g. a high dimensional search space with several local optima). They have been widely used for combinatorial process optimization, providing near-optimal solution, or solutions. In these cases, researchers and practitioners prefer alternative cost effective near-optimal (or approximate) solution(s) rather than exact ones, as it may be extremely difficult to find the precise, optimal point in higher dimension, and multimodal search space (Mukherjee and Ray 2006).
Furthermore, optimization problem complexity has increased as the resolution of process measurement method have also improved. The availability of better measurement technology has resulted in discrete and continuous search spaces with multi-modal, differentiable as well as non-differentiable, objective functions or responses. The next section briefly reviews the main computational approaches researchers have adopted for multi-variable optimization in the NNS manufacturing and, generally, in primary shaping and machining process optimization. Generally, the different approaches to process optimization can be classified as conventional and non-conventional. While conventional techniques provide a local optimal solution, non-conventional techniques, which are based on extrinsic models or objective functions, are only approximations that aim to provide near-optimal condition. Conventional techniques can be classified in two categories: experimental and iterative mathematical search techniques. The first includes statistical design of experiment (e.g. Taguchi method) and response surface design methodology (RSM). On the other hand, iterative mathematical search techniques include linear programming (LP), non-linear programming (NLP), and dynamic programming (DP) algorithms. Non-conventional techniques include the meta-heuristic search based algorithm. Three main types of metaheuristic search based algorithm are applied in process parameters optimization: genetic algorithm (GA), simulated annealing (SA), and tabu search (TS).
Evolutionary algorithms have been found particularity effective in multi-objective optimization (Zitzler et al. 2000;Fonseca and Fleming 1995;Srinivas and Deb 1994;Bäck and Schwefel 1993) and, in particular, in manufacturing sector (Pierreval and Tautou 1997;Fonseca and Fleming 1998).
In the literature, evolutionary algorithms have been applied to the optimization of several different production systems and manufacturing processes. For example, evolutionary algorithms have been applied to optimize the performances of production (e.g. Aydoğan et al. 2019) and assembly lines (e.g. Minghai et al. 2019;Saif et al. 2019). Optimization techniques have also been widely applied to scheduling manufacturing problems, for example for taking into account random events (Fleury et al. 1999).
In particular, several authors apply and modify GAs for optimising conventional and non-conventional machining process. For example, general multi-criteria approaches are developed and used for optimizing turning, wire cutting and laser cutting by Rao et al. (2018). Dedicated evolutionary algorithms are developed and tested on specific processes such as water-jet (Rao et al. 2019) or grinding .
Evolutionary algorithms have been used for concurrently optimizing the product design and the correspondent manufacturing cost. In particular, authors ,such as Ramesh et al. (2009), Al-Ansary andDeiab (1997) and Huang et al. (2005), focus on the optimal allocation of dimensional tolerances for product cost minimization using GAs and other evolutionary methods.
Several authors use evolutionary techniques for "on-line" process optimization. For example, Del Prete et al. (2019) use genetic programming for maximizing turning productivity (maximising removal rate) and reducing defect rate.  compare Taguchi approach to non-linear methods in the context of on-line process control of injection moulding, assessing the superiority in performances of the latter. Similarly, Zainal et al. (2016) compare swarm optimization with Taguchi techniques for optimizing milling process.
Evolutionary optimization techniques have also been applied to design optimization in the early conceptual design phase, as in Burnell et al. (1991),  and Lee and Wang (1992).

Metaheuristic algorithms for manufacturing process and design optimization: application to near net shape processes
Regarding NNS process, few authors attempt to associate the NNS concept to process chain optimization using evolutionary algorithms. Generally, the different investigation approaches are related to the process to be investigated. For example, Caporalli et al. (1998) investigate hot forging process, using numerical and experimental methodologies. Investigating two intercepting manufacturing lines (one producing a forging tool, another producing a workpiece though machining, annealing and forging). Their optimization model minimized the flow stress of every hot forging operation, using FEM simulation. Denkena et al. (2011) develop analytical and derived experimental models, minimizing a global cost function. The developed expression of the product quality is dependent on the tolerances achievable in machining and forging process. This formula is used concurrently with the total cost and leading time formulas for a multi-criteria optimization using Genetic Algorithms. Davidson et al. (2008) optimize a flow forming process using experimental trials, maximizing process quality (achievable tolerances) and optimizing microstructure using a Taguchi array (Design of Experiments). The optimization target is related to the process dynamics in the experimental approach. Sardinas et al. (2006) use a micro GA for optimizing a turning process by two objective functions: turning time and tool life. The algorithm uses two different populations: one static and one dynamic. The first is processed with the classic GA operators, creating the second population.
In an effort to improve the resolution and accuracy of machining process researchers have reported the application of many different computational models and algorithms. The optimization of process parameters is the key for economically efficient machining operations (Khan et al. 1997) because metal cutting is still a significant step in most manufacturing processes.
However these established processes need to respond to the increasing demand for quality, and so consequently optimization methods in metal cutting are vital tool for continuous improvement of process output quality. But finding optimal machining conditions is a complex problem and consequently metaheuristic searches (particularly the evolutionary algorithms) have been developed as "problem specific" solver, guiding, or modifying heuristics to produce optimal solutions. For example, Mukherjee and Ray (2006) provide the general applications advantages and disadvantages of evolutionary algorithm usage. Khan et al. (1997)) benchmark evolutionary algorithms, in particular GA and SA and gradients method, using different test function. GA and SA show consistency in converging to the global minima and are not influenced by the dimensions of the input vectors (i.e. variables number), although they increase the convergence time in comparison with differential methods. Regarding the differences between the specific models, SA gives high precision and the code can be run longer to get higher precision. For GA, the precision is proportional and limited by the precision used for representing each variable (i.e. number of bits). GA and SA advantages over gradient methods include the higher total number of function and constrains evaluable. Although, the higher functions evaluations per run (i.e. giving a high convergence time) is a disadvantage that make them unsuitable for real-time optimization.
In conclusion, GAs are most frequently used for single and multi-objective optimization of large (i.e. high number of variables) and complex problems (Deb 2001), because they: -Do not require continuous or convex design space (as non-gradient based) -Are able to explore large search space and have a high chance of avoiding local optima (due to is probabilistic and not deterministic nature) -Provide multiple near-optimal solution.
-Are able to solve multiple objectives and non-linear response function problems, in both discrete and continuous cases.
Therefore, the GA is generally preferred when near-optimal condition(s) rather than exact optimal conditions are cost effective and acceptable for use by the manufacturers (Mukherjee and Ray 2006).

Concurrent optimization of near net shape processes
Problem formulation A generic NNS manufacturing chain can be modelled as a primary shaping process followed by a number of machining steps (Fig. 1). The primary shaping component can be either a forming, casting or additive layer manufacturing process. The industrial need is to have a simultaneous optimization of design and process parameters, achieving: -Product geometric requirements (e.g. assembly constraints) -Quantitative product requirements (e.g. surface roughness, component weight) -Process quality (e.g. defect rate, tool wear, forging force) -Production targets (e.g. product cost, production rate, lead-time) The process parameters to consider are: where N is the number of machining operations required, n is the number of near net shape process variables, k is the number of machining variables. Referring to Fig. 1, the design variables to consider are where V is the number of design variables.
Conditional design optimization (CoDeO) Figure 2 shows a schematic representation of the CoDeO methodology. The optimization method is associated with the primary shaping and machining process models. In general, for optimizing the product design and process parameters, it is necessary to have a cost model which is a function of both. The models needed for this optimization are: 1. Cost model of the primary shaping process, dependent on process parameters (e.g. forging force or intermediate geometry) and/or product design variables (e.g. diameter) 2. Cost model of the machining process dependent on process parameters and/or product design variables 3. Support models (i.e. linear or non-linear) for the primary shaping or machining processes (e.g. production rate target or turning force constraint) The presence of support models is dependent on the nature of the processes, production targets and product requirements. A model can investigate the product quality (i.e. if it can be quantified by a function) or represent the technological aspect (e.g. forging force) to optimize (or keep within acceptable limits). By adding these as objectives function, the optimization changes from a single to a multi-criteria process. If these criteria can be added into other forms (e.g. non-linear constraints), the optimization can retain its single scope and optimising the cost by giving boundaries to the variable optimization.
These models allow the exploration of iterative systematic product and process design variation that can be optimized concurrently by minimizing the process cost (as main Referring to the Fig. 1 schematic, a machining process (milling or tuning) parameters can generally be synthesized in Eq. (1) where, a is the cutting depth, v c is the cutting speed F is the feed rate, N is number of passes.
Using geometrical constrains and machining parameters some of the intermediate machining geometries can be written as a function of the final and primary shaped geometries, by the number of machining steps (Eq. (2)).
This allows the number of machining step N to be expressed as a function of the design variables.
Expressing the machining steps number as in Eq.
(2), the total machining time T M , and, consequently, the total machining cost C M , can be written as Eq. (3).
In case of single criteria optimization, the objective function OBJ i can be written as in Eq. (4) where, c M is the hourly machining cost and C NNS the cost model for the NNS process. In case of multi-criteria optimization, the objective functions can be formulated as in Eq. (5) The cost model for the NNS process C NNS needs to be built individually, depending on the process nature and data available. Generally the cost is function of input and output geometries and process parameters. Marini et al. (2017a) present a framework for technological and economic models building for near net shape process.
Non-linear constrains can be added to the objective functions, limiting the feasible solutions to those which match process limit capabilities (e.g. turning power constraint), quality targets (e.g. surface roughness constrain or product weight) and production targets (e.g. machining removal rate).
NNS process variables and models are dependent on the process nature, which could be established from literature or experimental investigation.
The design variables selection is similarly case dependant and can be influenced by several factors. The main constraints are usually geometrical or connected to key features for influencing component functionality or manufacturing costs however their definition is difficult to standardize.
The machining process variables are usually limited to the feed rate, cutting speed and cutting depth.

Case study
CoDeO has been applied to the optimization of the manufacturing of a control valve's cage, a complex component manufactured from 420 Stainless Steel. The schematic of a manufacturing chain for producing valve's component is showed in Fig. 3.
A centrifugal casting process and turning process are used for making the component. For brevity of presentation, the finishing turning and drilling processes are excluded from the optimization, although the finishing allowance is used for modeling the machining process (described in "Appendix B").
Variables definition: search space setting The optimization considers processes cost so consequently these need to be modelled into equations, optimizing both process parameters and design variables Fig. 4 displays the variables to be optimized and their collocation in the manufacturing chain and final/semi-finished products designs. Nine variables (defined in Table 1) have been considered for the optimization problem, as in Fig. 4: four process parameters and five design variables.
Using the centrifugal casting process parameters selection and modelling described in  and included in the cost equation (9), the diameter of the centrifugal castings internal cavity ( ID m ) is set in accordance to mould available and minimum required molten material for spinning. The final component dimensional variables ( OD, ID, OD b , l f , L f ) are informed by current assembly constraints and design targets. The cutting speed ( v c ), depth of cut (a) and feed rate(F) ranges are defined using (Kalpakjian and Schimd 2009) for medium and high C steels and ceramic coated carbide tool. This definitions include the values currently employed in the manufacturing process. The machine tool spindle speed (n) is also used as a variable: its range derived from the cutting speed ranges and the minimum and maximum diameters.
To determine appropriate GA functions to support the optimization the literature was reviewed.
Literature survey: GA functions and settings for turing process optimization Many different applications for turning process optimization through metaheuristic search based algorithms have been reported in literature, in particular for GA. Cus and Balic (2003) use a GA search base algorithm for optimizing the turning process of steel. Saravanan et al. (2003) use a lower number of variables (i.e. feed rate and cutting speed) for testing turning process optimization through GA and SA. Amiolemhen and Ibhadode (2004) differentiate between roughing and finishing operations, introducing different unit production cost fitness function and use different constrains for every analysed operation. Sardinas et al. (2006) use a micro GA for optimizing a turning process by two objective functions: turning time and tool life. The  first is processes with the classic GA operators, creating the second. This new individuals are added to the static ones in the next iteration, meanwhile the other population remain unchanged, even if they violate the constraints (measured by the violation grade by an unfeasibility index). Yildiz and Ozturk (2006) propose an integration of Taguchi robust design optimization (i.e. using S/N ratio and ANOVA) and GA. The authors test the computational power of single and multi-criteria optimization through this hybrid algorithm, comparing it with previous hybrid GA in literature. Yang and Natarajan (2010) compare a multi-objective differential evolution algorithm and non-dominated sorting genetic algorithm for the turning optimization of a particular tool (steel and tungsten carbide). The differential evolution (DE) algorithm acts on GA by narrowing down its search space during each step. D'Addona and Teti (2013) use the MAT-LAB GA optimization toolbox for optimizing the turning process of cast steel blank with an HSS tool, comparing two crossover operators (single point and two-points) and mutation coefficients. Tables 2 and 3 summarize the approaches proposed for optimization of single pass turning process. Most of the authors use only the main turning process parameters as variables (feed rate, cutting speed and depth of cut). Similarly, all the authors used the unit production cost as fitness function (or as a part of it). Similar settings are applied also in the crossover and mutation fractions. All the authors include linear and non-linear constraints for limiting the individuals' feasibility.
For multi-passes turning, the researchers tend to use iterative cycles, optimizing each pass individually, usually, having every pass interacting with the previous and next ones. So essentially, the authors are applying single objective optimization, using the unit process cost (Onwubolu and Kumalo 2001;Yildiz 2013;Alberti and Perrone 1999;Wang et al. 2002) or customized linear function combining tool life and cost (Schrader 2003).
Objective functions Having determined that the multicriteria optimization method will be a GA (section "Concurrent optimization of near net shape processes"), two objective functions must be defined. The first objective function (Eq. (7)) is the total cost (i.e. including centrifugal casting and machining cost). The second fitness function is the turning force (Eq. (8)), developed by Sardinas et al. (2006). The two functions are: 1. Total cost model (7) FF 1 = C CCB + T tu c m  Sardinas et al. (2006) Centrifugal casting model The first part of Eq. (7) refers to the centrifugal casting blank cost ( C CCB ), which is defined by the Eq. (9) With, OD m the mould outer diameter and L f mould length.
The selection of OD m and ID m from the target casting dimension, and the analytical relationship with ID m are detailed in .K c and K v constants value, depending on the centrifugal casting operation, considered 5.56 and 3.17 respectively for the selected material and mould external diameter .
Turning model The second part of Eq. (7) refers turning cost, which is proportional to the turning time ( t i T ) and the cost per hour ( c m ). The total turning cost is dependent on the machining strategy, showed in Fig. 5). The turning process is divided into three main operations, which correspond to three different machining time per pass ( T Tu ). The total turning time can be written as: With, t EX T , external turning time; t EX bT : external bottom turning time; t IN T , internal turning time (Fig. 5).
The turning process has been analysed as a helicoidal motion, resembling the cutting trajectory, instead of a linear one (i.e. classic turning formulation), resembling the tool motion. The aim of this new formulation is to gather the process parameters, diameters and number of passes in a (10) The derivation of the turning time formula is showed in "Appendix A".The generic helix turning time can be written as in Eq. 11.

T Tu = t IN T + t EX T + t EX bT
In this way, the model takes into consideration the turning process dynamics when optimizing the process parameters and the design variables. The solution is consequently a trade-off between process time minimization (thus turning cost) and the cutting force minimization. Non-linear constraints According to the literature, nonlinear constrains have been add to the GA, limiting the feasible individuals to matching the turning process limit capabilities (i.e. turning power), quality targets (i.e. surface roughness) and production targets (machining removal rate). Turning power (Eq. (12) (2010))are used as non-linear inequality constraints for the machining parameters. The machining removal rate constrain (Eq. (14)) ensure the required production volume by a machining parameters feasible setting (i.e. previous constraints and the "traded-off" by the turning force fitness function). The non-linear constraints are showed as follows.

Turning power ( W Tu )
With, maximum turning Power ( W Max ) of 82 kW, correspondent to the target machine capabilities 2. Surface roughness ( Ra Tu ) With, minimum machining removal rate higher ( MRR Min ) higher than the current one ( 1.39 10 5 mm 3 ∕min ). The GA settings are summarized in Table 4 Case study results The two objective functions (Equations (7) and (8)) are concurrent and the solution form a Pareto frontier. Accordingly to the machining parameters increase, the cutting force increase, remaining into feasible ranges. The surface roughness improvements cause the turning power to increase. The largest improvements are seen in the machining removal rate and, so, the machining time and cost decreases as a consequence.
Regarding the optimized variables, it is possible to notice a general reduction of the final dimensions, with the exception of the bottom section length, which increases. The machining parameters can be set higher than those currently applied ones (except from the feed rate that remaining unchanged). The following observations can be made 1. Design parameters vary in the first generations before converging over an optimal set of variables (constant for all the last generation individuals). 2. Casting process parameter behave similarly to the design parameters (connected by geometrical relationships). 3. Feed rate reach an optimal constant value. 4. Depth of cut and spindle speed vary in every solution. 5. Total cost varies and remain coherent with current total cost 6. Cutting force and turning power vary congruently with machining parameters 7. Surface roughness varies with the machining parameters, satisfying the target requirements. 8. Machining removal rate varies with the machining parameters, although it remains sensibly higher than the current rate (average of 5.28 10 5 mm 3 ∕min). 9. Similarly, machining time varies with machining parameters, remaining lower than current value (75 min on average). Tables 5 and 6 show the GA results, displaying optimized set of variables (fittest individuals) and resulting constraints, objectives functions and performances. Considering the average results of the optimal variable sets, the optimized geometry and process parameters adoption show an improvement in performances. Optimized machining parameters gives up to a 9.5 times higher removal rate. Using centrifugal cast blank, the machining time decreases between 8.3% and 21.6%. The latter process configuration has been applied by the manufacturers and resulted in a cost saving of 22% (The Weir Group PLC 2015). Considering the selected configuration (Table 7), the improvement of process performances (i.e. interaction between process parameters and design variables) combined with the adoption of optimal design variables resulted in a final weight reduction of 10%. This 19 kg material wastage reduction resulted in an estimated 2.8 ton scrap reduction per year (The Weir Group PLC 2015).

Discussion
Many reported NNS analytical approaches aim to optimize cost (Denkena et al. 2011;Sardinas et al. 2006), while the empirical experimental approaches tend to maximize process quality (Caporalli et al. 1998;Davidson et al. 2008). CoDeO optimize concurrently the manufacturing chain's cost and process feasibility, depending on analysed process characteristics and global cost model. For example, turning force has been used as trade-off (i.e. Pareto optimization) for the multi-criteria optimization of turning process using Genetic Algorithm. Similarly to Sardinas et al. (2006), The minimization of turning force makes the turning operation more feasible, as the minimization of costs tended to increase the machining parameters.
CoDeO is not able to use tool and equipment as design variable (unlike Caporalli et al. 1998's expert system that optimize the forming die shape and characteristics). On the other hand, CoDeO can optimize simultaneously process parameters and design variables.
In contrast with other NNS approaches (Caporalli et al. 1998;Denkena et al. 2011), CoDeO optimization can be adapted to different NNS processes.
In comparison with other concurrent optimization approaches, CoDeO achieves similar optimization targets to the methodologies used for robotic systems (Bi and Zhang 2001;Li et al. 2001;Rampersad 1995;Ouyang et al. 2004) in terms of process parameters and design (both subject to design constraints). CoDeO is not an "offline" methodology and its scope is limited to the concurrent optimization of the product design and process parameters of NNS manufacturing chains. In order words it is unable to influence the manufacturing control system ("on-line"). While the Design for Control methodologies are designed for the flexible manufacturing systems, CoDeO application is limited to NNS processes. However, CoDeO is able to consider a broader range of primary shaping processes (from sand casting to additive manufacturing) and different modelling methods (e.g. experimentally derived models). Although CoDeO has flexibility in the manufacturing processes and product designs considered. Its search space is limited in comparison to optimization methods that integrate both design and operations management Ha and Porteus 1995;Zhang and Wang 2016). Other researchers have tried to generalize the modelling of entire manufacturing ecosystem . Their frameworks take a global approach to the enterprise and supply chain, including the product design and operational phases.
While the resulting models are more flexible and holistic, their implementation is more complex as an additional phase of abstraction is required for utilizing their models in practice. In contrast, CoDeO limits its scope to a particular manufacturing sub-system of the whole production system (Giachetti 1998), although modelling it with higher accuracy.
CoDeO can be integrated with this general methodology in order to provide mathematical modelling and concurrent optimization of NNS manufacturing chains.

Conclusions
The application of the CoDeO methodology to the case study application resulted in significant saving in material and cost. The holistic approach over manufacturing optimization enable the creation of a reliable model for optimizing both economic and technological sides of a NNS manufacturing chain. CoDeO can be mathematically integrated with other models for generating a holistic production models and concurrent optimization.
Similarly to other general models, the methodology is not suitable for non-expert users, because of the choices that need to be made during the models development and algorithm setting phases. Variable and targets definition are critical for the optimization, however, in the context of a global NNS process chain, the general model developed defines the key outcomes (i.e. objective functions definition) and limit possible search space (i.e. variable selection).Further details of the methodology can be found in Marini (2020). Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.

Appendix A: Helix-turning formulation
This appendix details the representation of the helicoidal turning ar developed to support the optimization step. In Figs. 6 and 7, the turning helix has been displayed. In point O, the global frame of the component is composed by the axis x, y , z and the unit vectors ( i, j, k ). T is the point of contact of the tool. The local frame is composed of three unit vectors:c , tangent to the circumference, to which corresponds the cutting speed vector v c ; t , tangent to the helix trajectory, to which corresponds the helix speed vector v H ; a , parallel to the workpiece axis and coincident to the tool motion axis (liner trajectory), to which corresponds the feed speed vector v H . The helix equation parametric equation can be expressed as in (15) With, ∈ [0, k ] ; x, y ∈ [0, r] ; z ∈ [0, L].
The main features of the helix are: c = p 2 ; Helix pitch, p; Slope, c p ; Number of turns, T 2 ; Torsion, = c r 2 +c 2 ; Curvature, k = r r 2 +c 2 . The velocity of a point on the helix can be written as (16) Velocity vector can be written as in (17) If the motion is uniform, it is possible to write (18) Arc Length: the total helix arc length can be derived integrating the velocity (17) as in (19) and (20) Classic turning formulation The classic formulation use the following inputs: -F, feed rate (mm/lap) n, spindle speed (lap/min) a, depth of cut (mm) l, working engagement (mm) The feed velocity (21), cutting speed (22) and cutting time (23) classic formulations can be written as follows Helix equivalence For the helix the speed can be written as (24) Having the total number of turns T , the length (20) can be rewritten as (25) Given, z = c , with z = L , it is possible to write (26) The cutting time t c must be the same for both the formulations. Having L = l and t c = l Fn , the cutting time can be also written as (27) The hypothesis t c (helix) = t c (classic) can be dimensionally checked as in (28) Therefore, the helix pitch can be written as (29) Velocity hypothesis The helix speed and the cutting and feed rate speed must be complementary, so (30) can be formulated.

Appendix B: Turning model formulation
Referring to Fig. 8

Proof Machined Diameters (20 mm allowances from centrifugal casting blank)
For a single turning operation, the machining parameters can be considered as constant, so the following hypotheses result generally valid for the turning process parameters. In addition, in all three cases, the turning length can be considered constant ( L i = cost = L ) for a single turning operation. Therefore (47) can be written Therefore, using the turning helix formula, the total turning time can be written as (48) Distinction needs to be made for the external and internal turning.
In external turning case(D i = OD i ) Although with a = constant it becomes The last sum can be written as Therefore Therefore Similar conclusions can be drafted for the internal turning case ( D i = (ID i )) The general inner diameter after i-passes can be written as Therefore, dually to the outer diameter case Therefore the total external turning time formula become With N e = round OD 0 −OD f a Similarly the total internal turning time In the examined case, the total turning time can be written as (49) OD i = OD 0 − (a 1 + a 2 + ⋯ + a i ) (50) (57) 1 + 2 N I (ID + 10) + N I (NI+1)a 2 2 F 2 + 2 N I (ID + 10)