The Scientific World Journal

Volume 2014 (2014), Article ID 967140, 7 pages

http://dx.doi.org/10.1155/2014/967140

Research Article

The Study of the Optimal Parameter Settings in a Hospital Supply Chain System in Taiwan

Hung-Chang Liao,1,2 Meng-Hao Chen,1,2 and Ya-huei Wang3,4

¹Department of Health Services Administration, Chung Shan Medical University, No. 110, Section 1, Jian-Koa N. Road, Taichung 402, Taiwan

²Department of Medical Management, Chung Shan Medical University Hospital, No. 110, Section 1, Jian-Koa N. Road, Taichung 402, Taiwan

³Department of Applied Foreign Languages, Chung Shan Medical University, No. 110, Section 1, Jian-Koa N. Road, Taichung 402, Taiwan

⁴Department of Medical Education, Chung Shan Medical University Hospital, No. 110, Section 1, Jian-Koa N. Road, Taichung 402, Taiwan

Received 21 January 2014; Revised 29 July 2014; Accepted 4 August 2014; Published 27 August 2014

Academic Editor: Gonzalo Guillen-Gosalbez

Copyright © 2014 Hung-Chang Liao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This study proposed the optimal parameter settings for the hospital supply chain system (HSCS) when either the total system cost (TSC) or patient safety level (PSL) (or both simultaneously) was considered as the measure of the HSCS’s performance. Four parameters were considered in the HSCS: safety stock, maximum inventory level, transportation capacity, and the reliability of the HSCS. A full-factor experimental design was used to simulate an HSCS for the purpose of collecting data. The response surface method (RSM) was used to construct the regression model, and a genetic algorithm (GA) was applied to obtain the optimal parameter settings for the HSCS. The results show that the best method of obtaining the optimal parameter settings for the HSCS is the simultaneous consideration of both the TSC and the PSL to measure performance. Also, the results of sensitivity analysis based on the optimal parameter settings were used to derive adjustable strategies for the decision-makers.

1. Introduction

An optimal supply chain (SC) is a strategy plan that covers the network of suppliers, hospital’s manufacturing factors, and customers. An optimal SC will decrease the costs of material flow from suppliers to customers. Hence, most of the literature focuses on SC costs. García-Dastugue and Lambert [1] established an Internet-based system in order to integrate business processes using the flow of information to improve business practices in the SC. The big advantage of the Internet-based system was the opportunity it offered to purchase some items at a lower price. Wan et al. [2] used a simulation-based optimization framework for the manufacturing processes associated with SC. They used two small examples to simulate a three-stage SC for the function of total cost and to obtain the optimal inventory levels. Their results showed that the framework can accommodate chance constraints and scales up well. Yao et al. [3] explored the optimization of the ordering process in considering the minimization of a collaborative SC’s total costs by using a vendor-managed inventory. van der Vlist et al. [4] revised Yao et al.’s model [3] by considering redistributing risk, coordinating replenishment for multiple buyers, sharing downstream demand data, a choice of stock positioning, and giving suppliers the latitude to change delivery times and the authority to adjust delivery quantities in response to inventory developments. Thus, the use of the SC model could optimize total costs when coordinated inventory management is implemented. Frota Neto et al. [5] developed a sustainable logistics network and discussed how to evaluate the effects of the trade-off between the costs of the logistics network and its environmental impact. Ganga and Carpinetti [6] proposed a SC model and used fuzzy logic to forecast its performance as to cost, responsiveness, reliability, and so forth. They showed that fuzzy logic is a feasible technique that is helpful in managing SC performance.

Most studies of SCs are in the industrial field; only a few have focused on the hospital supply chain system (HSCS). Schut and van Bergeijk [7] indicated that hospitals in most developing countries have high medical costs and suggested that they should go through a centralized government agency to assist purchasing in order to decrease the general price level of pharmaceuticals. Shah [8] researched a strategy to optimize SC in the pharmaceutical industry and concluded that the pharmaceutical industry needs to balance future customer demand with the production planning necessitated by clinical and competitor activity. Hence, effective capacity utilization planning and robust infrastructure investment are important key issues. Lapierre and Ruiz [9] studied an innovative approach for helping hospitals improve their scheduling logistics by coordinating procurement and distribution operations while respecting inventory capacities. Ghandforoush and Sen [10] presented a decision support system for platelet production SC with flexible scheduling. They suggested that, to meet daily demand, the decision support system should implement a superior production and mobile assignment schedule.

On the other hand, in order to plan a SC with mathematical precision, many studies used an experimental design to construct the model. Holweg et al. [11] used the Taguchi method and a simulation of a multitier SC to investigate the impact of scheduling activities. They concluded with a set of recommendations on how to improve current vehicle supply systems by forecasting scheduling systems. Delavar et al. [12] proposed a SC to coordinate production scheduling and air transportation. A Taguchi experimental design was applied to obtain the parameter values that would best improve their performance. Tiwari et al. [13] addressed the problem of an integrated SC design, and in order to ensure high service levels, a novel algorithm combined the Taguchi technique with an Artificial Immune System that was used to solve the problem. Yang et al. [14] improved the robustness of SC information-sharing strategies using a hybrid Taguchi method. Their results showed that when customer demand is uncertain, e-shopping has the most robust performance.

Based on the above literature review and significance, the purpose of this study is to design and plan a robust HSCS. This study used a simulated HSCS to analyze the different factors involved and to decide which parameters would affect the performance of the HSCS. These parameters were analyzed and used to create a robust design for the HSCS.

The simulated HSCS used in this study was based on computation technology, which allowed for the construction of a dynamic system that included information acquisition, processing, and management. The computation technology included a full-factor experimental design and response surface method (), and a genetic algorithm () was applied here to obtain the optimal parameter settings for the HSCS. Three responses, including total system cost (), patient safety level (), and overall performance (), were considered in the simulated HSCS. This study also obtained the optimal parameter settings for the robust HSCS and explored the sensitivity analysis of the HSCS’s to derive the adjustable strategies. In addition, different weight combinations for exploring the trade-off effects of and were discussed.

The innovation of this paper is the use of the RSM and GA method to create the HSCS and to obtain the robust parameters’ setting. Based on the literatures, the authors cannot find the research using RSM to design the HSCS. So this is the knowledge gap this work aims to cover.

2. The Problem of HSCS Parameter Settings

In this paper, three regional hospitals formed an HSCS alliance that included suppliers and a centralized purchasing center (CPC). Centralized purchasing allows large orders to be placed, which may result in a discounted price for each hospital, thus giving them a competitive advantage by reducing each participating hospital’s total costs. However, the parameters used to coordinate the HSCS for centralized purchasing and distribution are very important. With robust parameter settings, the HSCS could attempt to minimize the and maximize the . In order to obtain robust HSCS parameter settings, the scenario for the HSCS is described below.

Three hospitals forecast the demand for clinical masks for each season based on past data. The demand random variable (uncertain demand in each period for hospital ) is the normal probability distribution , in which is the mean () and is the variance (). Here, the for Hospital is , for Hospital is , and for Hospital is . Because the for each hospital is independent of that of the others, the CPC collects the total demand from the three hospitals to calculate the random variable of the total ordering quantity and aggregates the three hospitals’ normal probability distribution . After aggregating the forecast demands, the CPC submits the total order quantity to a pharmaceutical company. The pharmaceutical company then delivers quantities according to the hospitals’ orders, meeting each hospital’s needs. The evaluation indexes of the HSCS’s performances are TSC and PSL. is procurement cost + inventory cost + transportation cost. The formulation is as follows:

is the total system cost of the HSCS. is Hospital , , or . is the period from the 1st to the 12th month. is the distributed quantity for period . is the procurement cost for each unit in period ; it is set to normal probability distribution . is the discount rate for the total ordering quantity; it is set to normal probability distribution . is the inventory level for period . is the rate of the holding cost per unit for period ; it is set to normal probability distribution . is the number of times orders are transported from the pharmaceutical company to the hospitals. is the pharmaceutical company’s shipping cost per time to the hospitals. When the transportation capacity is 500 units, it is set to normal probability distribution . When the transportation capacity is 230 units, it is set to normal probability distribution . When the transportation capacity is 100 units, it is set to normal probability distribution .

In addition, the HSCS’s is broken down as follows. Because there is a gap between forecasted demand and actual demand, the material flow of the HSCS [15–18] may result in a shortage and by extension affect patient safety. Hence, avoiding a shortage is another aspect of an HSCS’s performance. Service level can thus be a feasible index to measure quantity shortages. In this paper, the service level is the index of [19].

The is defined as “1 − the shortage of quantity/the actual demand,” formulated as follows:

where is the shortage quantity for period and is the actual demand for period .

The inventory and shortage level are defined in (4):

The following equation shows the quantity , which is subject to the safety stock () and the maximal level of inventory ():

From (1)~(4), we obtain the decision variable, .

The HSCS’s performances— and —will result in a trade-off effect because, by obtaining larger quantities, will improve because of a decrease in the probability of an inventory shortage; however, the inventory will increase .

The safety stock, maximum inventory level, transportation capacity, and HSCS reliability were selected as the parameters for the robust HSCS. Also, the four parameters affect the and . Generally, an increase/decrease of safety stock or maximum inventory level will cause an increase/decrease of and . An increase of transportation capacity will cause a decrease of transportation times; it will also affect the stock quantity and later affect the and the . The reliability of supply chains is defined as the ability to perform the promised service of customer expectations [20, 21]. Here, the performances include and , which are based on the distributed quantity. Hence, the reliability of HSCS is defined as the ratio of quantities supply in time for different suppliers. Table 1 shows the different levels of these parameters.

Table 1: The parameter levels in this study.

Parameters
Value

Level 1
Level 2
Level 3

A: safety stock
70
80
100

B: maximum inventory level
400
500
600

C: transportation capacity
100
230
500

D: the reliability of HSCS
99%
97%
95%

To solve this HSCS problem, this paper used a full-factor experimental design in which 81 () combinations of parameter levels were simulated. Each combination was simulated 1000 times in order to reduce the variation of the means of the and the ; the means of the and the were then calculated for each combination. In addition, because the and the are simultaneously considered here, to help the decision-makers to set the robust parameters’ level, it is necessary to transfer the two responses to one multiresponse objective. Hence, in order to integrate and into one , the ideal function was adopted here.

The normalized mean of for each combination () is defined as , where is the minimization of . Thus, the larger the the better, and the value is between 0 and 1.

The normalized mean of for each combination () is defined as , where is the maximization of . Thus, a larger is again better, and the value is between 0 and 1.

The is defined by the ideal function; hence, . Thus, the larger the the better, and the value is between 0 and 1.

The response surface method (RSM) [22] was applied here in order to obtain the regression model using the parameters and the responses—, , and . Because the RSM is a model building technique using statistical experimental design and least square error fitting, it can be used to approximate a response function in terms of predictor variables. The approximation of the response functions for the available data is set as follows. The linear first-order polynomial approximation of the response function is

The quadratic second-order polynomial approximation of the response function is

where , , , and are tuning parameters; is the error item; is the number of model parameters.

The Appendix shows the nonlinear regression formulations of for , for , and for . The nonlinear regression formulations are significant for (, ), (, ), and (, ). Also, the explained variances () are higher for (0.996), (0.999), and (0.993).

Furthermore, in order to obtain the optimal parameter settings, the GA was used to find the solution. A common method of terminating a GA is to test (after a specified number of generations) the quality and the convergence at the global optimum of the best members of the population against the tested problem definition. The roulette wheel approach was used as the selection procedure. The population size was set as 50, the mutation rate was set as 0.06, and the crossover rate was set as 0.5. The mutation rate and the crossover rate controlled the expected number of chromosomes to mate and the number of genes to mutate, respectively [23, 24]. The stopping condition of the GA procedure was set at 1,000 iterations or when the change in the previous 100 iterations was less than 1 percent. Parameters and were continuous variables, and Parameters and were discrete variables.

Table 2 shows the optimal parameter settings. For , the safety stock was set as 100 (level 3), the maximum inventory level was set as 400 (level 1), the transportation capacity was set as 100 (level 1), and the reliability of the HSCS was set as 95 percent (level 3); the was 1065479. For , the safety stock was set as 100 (level 3), the maximum inventory level was set as 550 (level 2.498), the transportation capacity was set as 500 (level 3), and the reliability of the HSCS was set as 99 percent (level 1); was 0.981. For the , the safety stock was set as 78 (level 1.810), the maximum inventory level was set as 400 (level 1), the transportation capacity was set as 100 (level 1), and the reliability of the HSCS was set as 99 percent (level 1); the was 0.962. To compare the results of decision behavior, when the decision-maker considers that the performance of the HSCS should depend on , the would be 1076531.743 and the would be 0.918. In this situation, the (1076531.743) is 1.037 percent higher than the derived (1065478.998) and the (0.918) is 3.146 percent higher than the derived (0.890) when the decision-maker considers that the HSCS performance should depend on . Additionally, when the decision-maker considers that the performance of the HSCS should depend on the , the would be 1065478.998 and the would be 0.890. In this situation, the (1065478.998) is 13.144 percent lower than the derived (1205528.075) and the (0.890) is 10.225 percent lower than the derived (0.981) when the decision-maker considers that the HSCS performance should depend on the . Finally, when the decision-maker considers that the HSCS performance should depend on the , the would be 1205528.075 and the would be 0.981. In this situation, the (1205528.075) is 11.983 percent higher than the derived (1076531.743) and the (0.981) is 6.863 percent higher than the derived (0.918) when the decision-maker considers that the performance of the HSCS should depend on . Hence, from the above comparison of decision behavior results, the HSCS’s for the parameter settings is the superior choice because the is 1076531.743, which is 1.037 percent higher than the derived when the decision-maker considers that the HSCS performance should depend on the . The is 0.918, which is 6.754 percent lower than the derived when the decision-maker considers that the HSCS performance should depend on the .

Table 2: The optimal parameter settings for  , , and .

Response

Total system cost ()
3
1
1
3
0.998

Patient safety level ()
3
2.498
3
1
1.000

Overperformance ()
1.810
1
1
1
0.962

3. Sensitivity Analysis

The sensitivity analysis was based on the most robust parameter settings when the performance of the HSCS depends on . The safety stock is from 70 (level 1) units to 100 (level 3) units, and the other parameter settings are fixed. The results (see Figure 1) show that the curve decreases from 0.993 to 0.970, which means that the increased. The curve increases from 0.929 to 0.945, which means that the increased. The curve slightly increases from 0.960 to 0.962 (the safety stock level is 1.810, i.e., 78 units) and then slightly decreases to 0.958. In addition, the maximum inventory level ranges from 400 units (level 1) to 600 units (level 3), and the other parameter settings are fixed. Figure 2 shows that the curve decreases from 0.988 to 0.905 (maximum inventory level is 2.4, i.e., 540 units), which means that the increased. The curve subsequently increases to 0.920, which means that the decreased. The curve increases from 0.936 to 0.994 (maximum inventory level is 2.4, i.e., 540 units), which means that the increased. It then decreases to 0.989, which means that the likewise decreased. The curve decreases slightly from 0.962 to 0.952.

Figure 1: The changed safety stock levels affecting , , and .

Figure 2: The changed maximum inventory levels affecting , , and .

When the transportation capacity increased from level 1 and level 2 to level 3, the results (see Table 3) show that decreases, which means that the increased. decreases slightly, which means that the decreased slightly. decreases. Table 4 shows that when the reliability of the HSCS increased from level 1 and level 2 to level 3, the in level 1 is 0.988, in level 2 is 0.984, and in level 3 is 0.987. decreases, which means that the similarly decreased. decreases slightly.

Table 3: The transportation capacity levels affecting , , and .

Transportation capacity

Level
Units

3
100
0.850
0.934
0.890

2
230
0.967
0.936
0.951

1
500
0.988
0.936
0.962

Table 4: The reliability of the HSCS levels affecting , , and .

The reliability of the HSCS

Level
Reliability

3
95%
0.987
0.921
0.953

2
97%
0.984
0.926
0.955

1
99%
0.988
0.936
0.962

In addition, to compare the trade-off effect, when different weight combinations are formed, where ; the results are shown in Figure 3 and as follows.

(1)When is from 0.1 to 0.4, increases from 0.902 to 0.912 and decreases from 1.007 to 1.002. The parameters’ setting for the safety stock is level 3, the maximum inventory level is from 2.453 to 2.084, the transportation capacity is level 1, and the reliability of HSCS is level 1.(2)When is from 0.4 to 0.5, immediately increases from 0.912 to 0.989 and immediately decreases from 1.00 to 0.934. The parameters’ setting for the safety stock level is 1.615, the maximum inventory level is level 1, the transportation capacity is level 1, and the reliability of HSCS is level 1.(3)When is from 0.5 to 0.8, increases from 0.989 to 0.994 and decreases from 0.934 to 0.929. The parameters’ setting for the safety stock level is from 1.615 to 1.00, the maximum inventory level is level 1, the transportation capacity is level 1, and the reliability of HSCS is level 1.(4)When is from 0.9 to 1.0, the is 0.998, and is 0.907. The parameters’ setting for the safety stock is level 3, the maximum inventory level is level 1, the transportation capacity is level 1, and the reliability of HSCS is level 1.(5)When is set to as the ideal function, , and the safety stock is level 1.810, the maximum inventory level is level 1, the transportation capacity is level 1, and the reliability of HSCS is level 1. This situation is near the situation when is set to ; in the (). is set between 0.4 to 0.5, and is between 0.5 to 0.6. It means that the ideal function weights and are set between 0.4 to 0.5 and between 0.5 to 0.6. Hence, when is less than the ideal function weight, the performance is better than . When is larger than the ideal function weight, the performance is better than .

Figure 3: The changed weights affecting , , and .

4. Conclusion

In this paper, the authors proposed a method of setting parameters to design a robust HSCS. The full-factor experimental design was used to obtain simulation data for the HSCS. The RSM was applied here to obtain a regression model that shows how the parameter settings affect the HSCS’s performance. Also, GA was used to search for the solution to the problem of identifying optimal parameter settings. The results show that the best choice for the HSCS parameter settings is when the decision-maker simultaneously takes the and the into consideration; in this scenario, the is only 1.037 percent higher than the derived when the decision-maker considers that the HSCS performance should depend on the . The is only 6.754 percent lower than the derived when the decision-maker considers that the HSCS performance should depend on the . Also, the adjustable strategies coming from the results of sensitivity analysis show that when the is considered and the HSCS parameters’ setting is set in the optimization (the safety stock: 78 units (level 1.810); the maximum inventory level: 400 units (level 1); the transportation capacity: 100 (level 1); the reliability of the HSCS: 99 percent (level 1)), if the safety stock is adjusted increasing/decreasing 1 unit, the will be changed (increasing/decreasing) 0.077%, and will be changed (increasing/decreasing) 0.053%. If the maximum inventory level is adjusted increasing from 400 units to 540 units, the will increase 0.593% per unit and will increase 0.414% per unit. In addition, if the maximum inventory level is adjusted increasing from 540 units to 600 units, the will decrease 0.250% per unit and will decrease 0.083% per unit. If the transportation capacity is adjusted to level 2 or level 3, the will increase 2.1% or 13.8%, and will decrease 0.000% or 0.002%. If the reliability of the HSCS is adjusted to level 2 or level 3, the will increase 0.04% or 0.01%, and will decrease 0.010% or 0.015%. Hence, the decision-makers can evaluate the changed and to adjust the parameters’ level as adjustable strategies when dealing with practical problems. The concept of adjustable strategies will play an important role in future research on HSCS.

In addition, to compare the trade-off effect between and , the results show that when the is less than the ideal function weight, the will be stronger to affect the than the . When the is larger than the ideal function weight, the will be stronger to affect than the . Future study may consider using any MOO (multiobjective optimization) method (epsilon constraint, normal boundary intersection) to obtain the optimal parameter settings for the HSCS and to evaluate the performance of the HSCS.

Appendix

Regression Formulations

Consider

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

S. García-Dastugue and D. M. Lambert, “Internet-enabled coordination in the supply chain,” Industrial Marketing Management, vol. 32, no. 3, pp. 251–263, 2003.

X. Wan, J. F. Pekny, and G. V. Reklaitis, “Simulation-based optimization with surrogate models: application to supply chain management,” Computers and Chemical Engineering, vol. 29, no. 6, pp. 1317–1328, 2005.

Y. Yao, P. T. Evers, and M. E. Dresner, “Supply chain integration in vendor-managed inventory,” Decision Support Systems, vol. 43, no. 2, pp. 663–674, 2007.

P. van der Vlist, R. Kuik, and B. Verheijen, “Note on supply chain integration in vendor-managed inventory,” Decision Support Systems, vol. 44, no. 1, pp. 360–365, 2007.

J. Q. Frota Neto, J. M. Bloemhof-Ruwaard, J. A. E. E. van Nunen, and E. van Heck, “Designing and evaluating sustainable logistics networks,” International Journal of Production Economics, vol. 111, no. 2, pp. 195–208, 2008.

G. M. D. Ganga and L. C. R. Carpinetti, “A fuzzy logic approach to supply chain performance management,” International Journal of Production Economics, vol. 134, no. 1, pp. 177–187, 2011.

F. T. Schut and P. A. G. van Bergeijk, “International price discrimination: the pharmaceutical industry,” World Development, vol. 14, no. 9, pp. 1141–1150, 1986.

N. Shah, “Pharmaceutical supply chains: key issues and strategies for optimisation,” Computers and Chemical Engineering, vol. 28, no. 6-7, pp. 929–941, 2004.

S. D. Lapierre and A. B. Ruiz, “Scheduling logistic activities to improve hospital supply systems,” Computers & Operations Research, vol. 34, no. 3, pp. 624–641, 2007.

P. Ghandforoush and T. K. Sen, “A DSS to manage platelet production supply chain for regional blood centers,” Decision Support Systems, vol. 50, no. 1, pp. 32–42, 2010.

M. Holweg, S. M. Disney, P. Hines, and M. M. Naim, “Towards responsive vehicle supply: a simulation-based investigation into automotive scheduling systems,” Journal of Operations Management, vol. 23, no. 5, pp. 507–530, 2005.

M. R. Delavar, M. Hajiaghaei-Keshteli, and S. Molla-Alizadeh-Zavardehi, “Genetic algorithms for coordinated scheduling of production and air transportation,” Expert Systems with Applications, vol. 37, no. 12, pp. 8255–8266, 2010.

M. K. Tiwari, N. Raghavendra, S. Agrawal, and S. K. Goyal, “A Hybrid Taguchi-Immune approach to optimize an integrated supply chain design problem with multiple shipping,” European Journal of Operational Research, vol. 203, no. 1, pp. 95–106, 2010.

T. Yang, Y. Wen, and F. Wang, “Evaluation of robustness of supply chain information-sharing strategies using a hybrid Taguchi and multiple criteria decision-making method,” International Journal of Production Economics, vol. 134, no. 2, pp. 458–466, 2011.

C. E. Schmidt and T. Bottoni, “Improving medication safety and patient care in the emergency department,” Journal of Emergency Nursing, vol. 29, no. 1, pp. 12–16, 2003.

J. B. Dimick, J. D. Birkmeyer, and G. R. Upchurch Jr., “Measuring surgical quality: what's the role of provider volume?” World Journal of Surgery, vol. 29, no. 10, pp. 1217–1221, 2005.

N. Menachemi and R. G. Brooks, “Reviewing the benefits and costs of electronic health records and associated patient safety technologies,” Journal of Medical Systems, vol. 30, no. 3, pp. 159–168, 2006.

A. Rogers, E. Jones, and D. Oleynikov, “Radio frequency identification (RFID) applied to surgical sponges,” Surgical Endoscopy, vol. 21, no. 7, pp. 1235–1237, 2007.

H. Liao, “Using GA to dispatch the obtaining quantity for minimizing the total cost based on consideration of patient safety in HSCM,” Expert Systems with Applications, vol. 36, no. 8, pp. 11358–11362, 2009.

H.-P. Wiendahl, G. von Cieminski, and C. Begemann, “A systematic approach for ensuring the logistic process reliability of supply chains,” CIRP Annals—Manufacturing Technology, vol. 52, no. 1, pp. 375–380, 2003.

I. Koçoǧlu, S. Z. Imamoǧlu, H. Ince, and H. Keskin, “The effect of supply chain integration on information sharing: enhancing the supply chain performance,” in Proceedings of the 7th International Strategic Management Conference, vol. 24, pp. 1630–1649, July 2011.

G. E. P. Box and K. B. Wilson, “On the experimental attainment of optimum conditions,” Journal of the Royal Statistical Society B, vol. 13, no. 1, pp. 1–45, 1951.

M. Gen and R. Cheng, Genetic Algorithms and Engineering Design, John Wiley & Sons, New York, NY, USA, 1997.

Z. Khan, B. Prasad, and T. Singh, “Machining condition optimization by genetic algorithms and simulated annealing,” Computers and Operations Research, vol. 24, no. 7, pp. 647–657, 1997.