Project on Coordinating Demand Fulfillment With Supply Across A Dynamic Supply Chain

Description
A supply chain, as opposed to supply chain management, is a set of organizations directly linked by one or more upstream and downstream flows of products, services, finances, or information from a source to a customer.

ABSTRACT Title of Dissertation: COORDINATING DEMAND FULFILLMENT WITH SUPPLY ACROSS A DYNAMIC SUPPLY CHAIN Maomao Chen, Doctor of Philosophy, 2006 Dissertation directed by: Professor Michael Ball, Decision & Information Technologies Department

Today, technology enables companies to extend their reach in managing the supply chain and operating it in a coordinated fashion from raw materials to end consumers. Order promising and order fulfillment have become key supply chain capabilities which help companies win repeat business by promising orders competitively and reliably. In this dissertation, we study two issues related to moving a company from an Available to Promise (ATP) philosophy to a Profitable to Promise (PTP) philosophy: pseudo order promising and coordinating demand fulfillment with supply. To address the first issue, a single time period analytical ATP model for n confirmed customer orders and m pseudo orders is presented by considering both material constraints and production capacity constraints. At the outset, some analytical properties of the optimal policies are derived and then a particular customer promising scheme that depends on the ratio between customer service level and profit changes is presented. To tackle the second issue, we create a mathematical programming model and explore two cases: a deterministic demand curve or stochastic demand. A simple, yet generic optimal solution structure is derived and a

series of numerical studies and sensitivity analyses are carried out to investigate the impact of different factors on profit and fulfilled demand quantity. Further, the firm’s optimal response to a one-time-period discount offered by the supplier of a key component is studied. Unlike most models of this type in the literature, which define variables in terms of single arc flows, we employ path variables to directly identify and manipulate profitable and non-profitable products. Numerical experiments based on Toshiba’s global notebook supply chain are conducted. In addition, we present an analytical model to explore balanced supply. Implementation of these policies can reduce response time and improve demand fulfillment; further, the structure of the policies and our related analysis can give managers broad insight into this general decision-making environment.

COORDINATING DEMAND FULFILLMENT WITH SUPPLY ACROSS A GLOBAL SUPPLY CHAIN

by Maomao Chen

Dissertation submitted to the Faculty of the Graduate School of the University of Maryland, College Park, in partial fulfillment of the requirements for the degree of Doctor of Philosophy 2006

Advisory Committee: Professor Michael Ball, Chair/Advisor Professor Martin Dresner Professor Satyandra K. Gupta Professor Gilvan Souza Professor Paul Zantek

© Copyright by Maomao Chen 2006

Dedicated to Yanmei Zhao, my lovely wife Alina M. Chen, my daughter Bangying Chen and Yunju Yan, my parents

ii

ACKNOWLEDGEMENTS

I would like to express my gratitude to my advisor, Professor Michael Ball, for his support, patience, and encouragement throughout my graduated studies. It is not often that one finds an advisor that always finds the time for listening to the little problems and roadblocks that unavoidably crop up in the course of performing research. His technical and editorial advice was essential to the completion of this dissertation and has taught me innumerable lessons and insights on the workings of academic research in general.

I am also grateful to Professor Gilvan C. Souza and Dr. Zhenying Zhao for their invaluable advice, support and thought-provoking ideas.

I would like to express my deep appreciation to Dr. Dresner, Dr. Gupta and Dr. Zantek for being my committee members and providing insightful comments on the dissertation.

iii

TABLE OF CONTENTS List of Tables………………………………………………………………......viii List of Figures ………………………………………………………….…........ix Chapter 1: Introduction 1

1.1 Research Motivation …………………………………………………....2 1.1.1 1.1.2 Business Driving Forces ………………………………………. .....2 ATP Trend ………………………………………………………..4 ………………………………………………….....7

1.2 Research Questions

1.3 Contributions ……………………………………………………….……10 1.3.1 Research Contributions ………………………………………………10 1.3.2 Managerial Implications………………………………………………12 1.4 Overview of Chapters …………………………………………...……..….13

Chapter 2: Literature Review

15

2.1 Current ATP Research……………………………………………………...15 2.2 Stochastic ATP – Uncertainty……………………………………………...16 2.3 Coordinating Demand Fulfillment with Supply …………………………...18 2.4 Revenue Management ……………………………………………....……...21

Chapter 3: Pseudo Order Consideration in ATP

25

3.1 ATP Framework……………………………………………………..……..25
iv

3.2 Problem Formulation and Model …………………………………………27 3.2.1 Notation……………………………………………………………..28 3.2.2 Formulation……………………………………………………….....30 3.2.3 Model Analysis…………………………………………….………..32 3.3 Implementation Rule………………………………………………………38 3.4 Experimental Study and Results…………………………………………..39 3.5 Multiple Pseudo Orders in ATP …………………………………………..46 3.5.1 Notation ……………………………………………………………...46 3.5.2 Formulation………………………………………………………......48 3.5.3 Model Analysis…………………………………………….………...50 3.6 Remarks …………………………………………………………………...54

Chapter 4: Coordinating Demand Fulfillment with Supply

57

4.1 Analytical Model with Deterministic Demand Curve……………………...57 4.1.1 Model & Formulation………………………………………………...57 4.1.2 Model Analysis……………………………………………………….60 4.1.3 Numerical Experiments………………………………………………63

4.2 Analytical Model with Stochastic Demand ..………………………………..72 4.2.1 Model & Formulation………………………………………………....73 4.2.2 Model Analysis………………………………………………………..76

v

4.2.3 Algorithm ……………………………………………………………..77 4.2.3 Numerical Experiments ………………………………………………80 4.3 Analytical Model for Balanced Supply ..……………………………………85 4.3.1 Model & Formulation………………………………………………....87 4.3.2 Model Analysis ……………………………………………………….90 4.4 Available Resource Analysis…………………………………………….......91 4.4.1 Model Background……………………………………………….……92 4.4.2 Toshiba Global Supply Chain…………………………………………94 4.4.3 Model Formulation…………………………………………….………96 4.4.4 Result Analysis………………………………………………………..102 4.4.5 Remarks………………………………………………………...……...110

Chapter 5: Conclusions

114

Appendices Appendix 1…………………………………………………………………………117 Appendix 2……………………………………………………………………........118 Appendix 3……………………………………………………………………........120 Appendix 4…………………………………………………………….…………...122 Appendix 5……………………………………………….……………….…..........124 Appendix 6…............................................................................................................125 Appendix 7…...........................................................................................................126

vi

Appendix 8…...........................................................................................................128 Appendix 9…...........................................................................................................130 Appendix 10….........................................................................................................131 Appendix 11….........................................................................................................132 Appendix 12….........................................................................................................133 Appendix 13….........................................................................................................134 Appendix 14….........................................................................................................135 Appendix 15….........................................................................................................136 Appendix 16…....................................................................................................... .140 Appendix 17 ………………………………………………………………………142 Bibliography 143

vii

LIST OF TABLES 3.1 Effect of Pseudo Order’s Price Change…………………….. …………………41 3.2 Effect of Pseudo Order’s Variance……………………………………………..43 3.3 Effects of Service Level on Order Promising…………………………………..45 4.1 Experiment Coefficients ……………………………………………………….61 4.2 Q3 vs. Total Profit ……………………………………………………………...64 4.3 Q2 , Q3 vs. Total Profit ………………………………………….. …………....65 4.4 Q1 , Q3 vs. Total Profit ………………………………………………………....65 4.5 Q1 , Q2 , Q3 vs. Total Profit……………………………………………………..66 4.6 Effect of Transportation Cost …………………………………………………..67 4.7 Effect of Maximal Price …………………………………………………….......68 4.8 Effect of Price Elasticity ………………………………………………………..69 4.9 Experiment Coefficients ………………………………………………………..79 4.10 Effect of Demand’s Mean ……………………………………………………..80 4.11 Effect of Demand’s Stand Deviation …………………………………… ……81 4.12 Effect of Transportation Cost …………………………………………… …….83 4.13 Experimental Design for First Research Question …………………………....102 4.14 T–Statistic for Multiple Regressions Where Dependent Variables are Profit and Purchased Quantity and Independent Variables are Factors ………………………103

viii

LIST OF FIGURES 1.1 ATP Decision Space……………………………………………………………..6 1.2 Illustration of Path Variable……………………………………………………..12 3.1 Pseudo Order Decision Flow Chart …………………………………………….40 3 2. Committed Order as a function of pseudo order’s Price……………….…….....42 3.3. Total Benefit as a function of pseudo order’s Price…………………….…........42 3.4 Committed Order Quantity as a Function of Pseudo Order’s Variance……........43 3.5. Total Benefit as a Function of Pseudo Order’s Variance……………….……....44 3.6. Total Committed Quantity of Confirmed Order ………………………………..45 4.1 Coordination in Supply Chain (1) …….………………………..………………56 4.2 Coordination Algorithm (1)……………………………………………………..62 4.3 Coordination Algorithm (2)………………………………………………….….63 4.4 Q3 vs. Total Profit …………………………………………………….. ………64 4.5 Q2 , Q3 vs. Total Profit …………………………………………………….. ..….65 4.6 Q1 , Q3 vs. Total Profit …………………………………………………………..66 4.7 Q1 , Q2 , Q3 vs. Total Profit …………………………………………………..…66 4.8 Effect of Transportation Cost …………………………………………………...67 4.9 Effect of Maximal Price …………………………………………………………68 4.10 Effect of Price Elasticity ……………………………………………………….69 4.11 Supply Chain Path (2) …………………………………………………………72 4.12 Coordination Algorithm (3) ……………………………………………………77

ix

4.13 Coordination Algorithm (4) ……………………………………………………78 4.14 Effect of Demand’s Mean …………………………………………………….80 4.15 Effect of Demand’s Stand Deviation ………………………………………….81 4.16 Effect of Transportation Cost ………………………………………………….83 4.17 Supply Balance ………………………………………………………………...86 4.18 Toshiba Global Supply Chain …………………………………………………94 4.19 Effect of Holding Cost, Capacity, Fill Rate and Discount on Profit ………….105 4.20 Effect of Holding Cost, Capacity, Fill Rate, and Discount on the Purchased Quantity …………………………………………………………………………....106 4.21 Effect of Capacity Increase in the Production of Profitable Products ……..…106 4.22 Profit vs. Price Discount …………………………………………………...…109 4.23 Purchased Component Quantity vs. Price Discount ……………………….…109

x

Chapter 1 Introduction

Today’s technology enables companies to extend their reach in managing the supply chain and operating it in a coordinated fashion from purchasing raw material to fulfilling end consumers’ demands. Traditional cost and profit based supply chain strategies are no longer sufficient in the present competitive business environment. Leading companies are creating synchronized supply chains that are driven by market needs and, in essence, are moving the supply chain closer to the customer. As a result, demand fulfillment capabilities have become the key to the competitive strategies of many companies. Available to Promise (ATP) directly links customer orders with enterprise resources to achieve supply chain optimization. ATP had its origins in the late 1980’s with Manufacturing Resource Planning (MRP II). Traditionally, the ATP function provides a response to customer order requests based on resource availability by checking the uncommitted portion of a company’s inventory and planned production, maintained in the master schedule to support customer order promising (Ball, Chen, 2002). Supply Chain Management (SCM) introduces processes and systems to generate an ATP that is feasible and optimal with respect to resource constraints (Ervolina, 2001). Since this ATP strategy has the ability to optimize resource utilization through complicated material and process constraints, it is also referred as “advanced” ATP (Chen and Ball, 2000). Due to the complexity of ATP, only a very limited number of papers present quantitative models to support ATP

1

(Ball and Chen, 2001). One objective of this dissertation is to introduce the analytical model to deduce the generic rules that can provide managers with useful insights into the optimal policies for improving demand fulfillment.

1.1 Research Motivation
Our research is motivated both by business needs and gaps prior ATP research.

1.1.1 Business Driving Forces
Global competition and widely adopted e-commerce business models have imposed tremendous pressure on product and service providers to get closer to their customers. At the same time end consumers are increasingly knowledgeable and demanding. Supply chains are confronting the essential challenges in the current customer-centric business environment: real-time responsiveness, uncertain customer orders, globally dispersed locations and diminishing profit margin. As the front-end of a supply chain, order management must treat these challenges as the diving forces to gain the advantage. Detailed business transaction information has become accessible in real–time mode or near real-time modes throughout the supply chain, providing the possibility of real–time order management and optimization. Meanwhile, broad application of ecommence technology has challenged demand fulfillment manner and created the needs for new order promising styles. Customer order response time has become critical to customer satisfaction, especially when a real-time customer response is

2

required. If a company does not meet customer expectation in “real-time”, customers may look towards their competitors while waiting for their order promise. In addition, by responding in real-time to the customer, manufacturers and suppliers could also better collaborate to present jointly constructed campaigns to end-customers and therefore provide both the manufacturers and suppliers with unparalleled means of promising and earning new business (Zweben,1996). As the number of customer orders increase, batch ATP becomes inefficient; and rule-based decision mechanism becomes a requirement for achieving real time response. The solutions from the analytical model in this dissertation provide such a mechanism. Uncertainty is another challenge for ATP. Uncertainties across a supply chain generally come from inaccurate forecasting. Under severe competition, companies have to offer customers more flexibility canceling orders. It has become common for some customer orders to not show up or for customers to make changes that require "what-if" problem solving around cancellations, substitutions or reshuffling of orders. According to Fisher (1997), customer uncertainty is inherent in order promising and has a considerable impact on the supply chain structure. Similarly, uncertainty can be from supply side such as, e.g., delayed delivery of raw material or factory shutdowns. The company can hedge against uncertainty with excess inventory or excess capacity but this results in high inventory cost and capacity waste. The effective exploitation of uncertainty allows the supply chain to better calibrate service levels to meet the needs of various customer segments, as well as to reduce costs. Therefore: considering uncertainty in the ATP decision-making process is necessary to provide a greater degree of stability, continuity and predictability in the customer base in order

3

to earn the more business. We are going to study the uncertainty from the customers’ perspective by introducing the pseudo order. Today’s supply chains are continually increasing in complexity. Low level, even negative profit margins for a single product may make sense in the context of a large, complex supply chain. With improved communications and increased competition, consumers have been provided with more choice, while most competitors are similar in product performance, quality and price. As consumers expect new products, better quality, and shorter lead times at a reasonable price, strategic use of non-profitable products is not unusual. (Bhattacharjee, 2000). In facing this predicament, companies are showing enthusiasm in discovering how they can better use profitable products to serve their customers. From the supply chain’s perspective, this means products have to reach the customers from the right supplier. Supply chain diversity ranging from globally dispersed manufacturers, distribution centers and sales subsidiaries with different production cost, capacities, capabilities and lead-times for different products demands identifying the profitable path for effective order promising, capacity utilization and production smoothness. Overcoming the diminishing profit margin and achieving resource allocation efficiency stimulates us to perform path analyses for the global supply chain.

1.1.2 ATP Trend As the name suggests, the ATP function provides information regarding resource availability to promise delivery, in response to a customer order request. Conventional ATP quantity is a row under the Master Production Schedule (MPS),

4

and is responsible for keeping track of the uncommitted portion of current and future available finished goods. Unlike conventional ATP, which assigns existing inventory or pre–planned production capacity, advanced ATP refers to a systematic process for making best use of available resources including raw materials, work–in–process, and production and distribution capacity, in addition to finished goods. An increasingly dynamic and customer-centric environment is heightening the requirements in which companies perform order promising and fulfillment. However, even the most expensive and complex commercial ATP currently available typically promises orders on an incremental, first-come, first-served basis, and as such, have some obvious drawbacks: ? ? They do not consider the opportunity cost associated with committing supply to a particular order; for example, promising supply to a lower-margin order may preclude that supply from going to a higher-margin order that has yet to come in; ? They do not attempt to maximize the potential revenue of each order; ? They do not distinguish the products from different paths. In addition, the current advanced-ATP research generally focuses on two elements: quantity and due-date. Quantity quoting gives the customer flexibility often seen in the supply contract. Due-date quoting gives the customer a time buffer in which the order has to be delivered. We know that the primary purpose of the ATP function is to provide a response to customer orders. There are two levels of response: customer response space -- one involves the direct response to the customer and the most fundamental decision related to an order is whether or not to accept the order and if accepted, its committed delivery date and quantity; product processing space --

5

this involves the underlying activities required of the production and distribution systems to carry out the customer commitment. We can see that the decision based on quantity and delivery date elements is only in customer response space. Therefore, we present an additional element of ATP: product path (see the Figure 1.1). In addition, as we described before, uncertainty resulting from order cancellations is a critical factor in ATP decisions regardless of the decision space. However, only a few stochastic push-based ATP models have been built, and little work related to stochastic pull-based ATP has been completed. Here, push-based ATP models are designed to allocate available resources for promising future customer demands, and pull-based ATP performs dynamic resource allocation in direct response to actual customer orders (Chen, Zhao and Ball, 2001). According to this definition, our research falls within the pull-based ATP domain. Incorporating path analysis and uncertainty into our model allows us to set and manage customer expectations with accurate supply availability and build a responsive, agile and truly customer-centric supply chain.
Figure 1.1: ATP Decision Space
ATP Decision Space

Uncertainty

Quantity Due Date Customer Response

Product Path Supplier

Space

Product Processing

Space

6

1.2 Research Questions
Available-to-Promise (ATP) applications originated as a means for controlling the allocation of finished goods inventory and improving the quality of delivery promises to customers. It has since developed into a major operational tool that supports the management of customer demands, safety stocks, production efficiency and the available resource. ATP demonstrates the tremendous synergistic opportunities available within integrated manufacturing planning and control systems.

Unfortunately, though easily understood by most users and characterized by some researchers, many companies defer the development and/or implementation of ATP due to the shortage of the efficient ATP systems to clearly and visibly link the external commitments to the supporting manufacturing plan. These are the questions we would address at the strategic level: Question 1: Under what conditions do some of the commonly used ATP rules perform well? Are there any other appropriate rules? How can the model parameters be effectively set? We believe that our rule-based ATP solution is the answer to this question. Rulebased results can be obtained from analytical models and their solution can be easily implemented and deployed in Decision Support Systems (DSS).

Our previous discussion has clearly shown the value of explicitly including product path analyses in addressing ATP decision. It also shows the importance of

7

including model of uncertainty in such analysis. Accomplishing these objectives is not simple and is our major focus. Profitability and customer service are two fundamental drivers in determining a company’s performance. Of course, without profitable orders, a business cannot survive. By putting customers at the center of the supply chain, and using information about customer needs to drive it, companies can lower costs, boost revenues and greatly increase customer satisfaction. A clear understanding of customer profitability is critical, because it enables the organization to differentiate the various level of service it provides to various customer segments according to their needs and value to the company. It is important for different customers to get the service that is most appropriate for their needs and the company's profitability. A comprehensive view of customer profitability and customer service lets companies focus resources where they will do the most good in terms of strengthening key customer relationships and bolstering top-line growth. This leads us to the question below: Question 2: How can uncertainty be incorporated into ATP optimization models? We are interested in answering Question 2 by developing an analytical ATP model, which generates a set of business rules to guide ATP execution. We

incorporate both profit and customer service in the objective and include consideration of the pseudo orders. Here, we use a simple pseudo order to aggregate all potential future customer orders and the uncertainties surrounding them. The simple solution structure as well as the empirical result is provided. We devote Chapter 3 to this study.

8

It has been generally recognized that the coordination of demand fulfillment and purchasing is of critical importance to marketing and product managers, because this translates into increased customer satisfaction and cost. To counteract price erosion and the accompanying reduction in profit margin, manufacturers need to align production and logistics planning with end-sales to choose the right path. In addition, to avoid vulnerability in market competition and reduce the risk, the company also needs the right path. When we say “product path”, we refer to the path from the supplier where the raw materials are purchased to manufacturer where products are produced, through the distribution center, to the sales locations where the products are ultimately sold to the customers. Obviously, as the cost varies from every path, we need to ask: Question 3: What kind of strategies and models produce effective demand fulfillment through the right supplier/supply chain paths? We build constrained non-linear integer programming models in Chapter 4 to decide optimal supplying quantity of the products from each supplier to meet end customer’s demand so that the profit of the company is maximized.

We also investigate the dynamic management of the available resource in accordance with the order promising, for example, if the supplier offers price discount on the component /raw material, how should the company respond to such situation? We think the resource and order management (sales) interact with each other, and should be managed in such a way.

9

Question 4: How should companies coordinate the raw material and end product discounts? We use the Toshiba global supply chain as a case study, and introduce “path variables” to build a MIP (Mixed Integer Programming) model. Numerical results are presented in Chapter 4.2. We introduce “path variables” to directly determine the unit cost of all products, and identify, for example, the percentage of demand that is met by those cost effective paths. It is necessary to employ “path variables” to obtain such information and produce the appropriate decisions. This analysis leads to a question raised above: when can a discount on raw materials generate more profit, and when should a complementary discount on sales prices be offered to stimulate demand?

1.3 Contributions
1.3.1 Research Contributions

This dissertation provides several contributions to ATP research: Production pooling considerations in ATP models: Production cost is particularly significant in the ATP research as it is always a major factor in affecting resource allocation. The shapes of the production cost functions depend on many issues (Ghali, 2003), such as industry, the length of the time horizon, capacity, and even the product life cycle. It is recognized that production cost curves, with current technological change, become concave. We build ATP models that evaluate the impact of such changes.

10

Supply chain collaboration to achieve better demand fulfillment: We build our models to reflect global supply chain goals, while addressing demand fulfillment so as to enable sellers, distributors, manufacturers and suppliers to easily satisfy the endcustomer and to help them collaborate in sales, marketing and service initiatives. We should clarify, however, that our current research is related to classical revenue management but has certain differences as well. Revenue management encompasses all practices of discriminatory pricing used to enhance delivery reliability and maximize the profit generated from the resources. The key question facing us is how to allocate the resources shared by the various products, which some RM models address. In addition, demand fulfillment involves issues such as inventory and the resulting holding cost that are not covered in revenue management.

Employing path variables to directly identify and manipulate profitable and non-profitable products: Unlike most models in the global supply chain literature, which define variables in terms of flows along a single arc in the network and use flow balancing constraints at nodes, we employ path variables, which provide location–specific cost information and directly identify profitable and non–profitable products (see Figure 1.2). In Figure 1.2, xij and y jk are traditional flow variables for the product in the network. Demand is represented by d k at node k. Note that this choice of variables does not capture the unit costs for products sold at node k (for example, it is not possible to determine whether a particular product at node k comes from node i or i’.) By defining a path variable zijk , one can directly determine the unit

11

cost of all products, and identify, for example, the percentage of demand that is met by profitable products. In our two–year research project with Toshiba, we found this to be very important managerial information, as it serves as an aid for other strategic decisions such as product line offering at different locations. Although this modeling approach certainly increases the number of decision variables (for example, a typical model we study here has over a million variables and a similar number of constraints), we find that the resulting models, even with real–world data, are manageable with solution times being around 5 minutes using typical computing environments. Figure 1.2: Illustration of Path Variables k y jk i’

dk

j

xij
i

zijk

j’

k’

1.3.2 Managerial Implications We would like to show the practitioners that our research can also serve as a strategic weapon along several dimensions: ePromise Capability: Through our easy-to-implement analytical ATP solutions, a company can provide a so-called ePromise capability, which supports real-time resource allocation based on actual availability and dynamic order requests, in accordance with the company's business objectives. Such capabilities enable a

12

company to automate order-promising decisions based on electronic information, which is essential for an e-commerce environment. Enhanced Revenue and Service: Embedding the pseudo (future) order into confirmed orders to deal with uncertainty and finding the “right supplier” to serve the customer based on varying production status allow companies to identify customers’ real demand and provide a clear understanding of their internal capabilities, thus enabling managers to enhance service and profit simultaneously. Buy Smart – How to Benefit from Recession: Our model also identifies the best quantity of the raw materials/components to buy from suppliers when they offer price discounts. This helps understand how to coordinate raw material discounts with the end sales discounts to achieve the best resource utilization. Such a proactive approach to managing procurement can make a substantial difference during a recession, and can help managers capitalize on future opportunities. As a result, the proactive purchaser is able to take advantage of the recession and shape the supply chain to their long-term advantage.

1.4 Overview of Chapters
The remaining chapters in this dissertation are organized as follows. Chapter 2 gives the literature review. In Chapter 3, we introduce an analytical model that considers multiple confirmed orders, multiple pseudo orders and production pooling. This model can evaluate the usefulness of responding to customer enquiries in real time. In Chapter 4, we create mathematical programming models for two scenarios: a

13

company facing a certain demand curve or uncertain demand. These are strategic level models that provide insight into what quantity levels to purchase from multiple suppliers based on several cost factors. We also develop a mixed integer programming model that addresses how to coordinate raw material discounts offered by s supplier with end-consumers policies. In Chapter 5, we conclude the dissertation by summarizing the results and the contributions.

14

Chapter 2 Literature Review
2.1 Current ATP Research
Traditional ATP systems are based on the Master Production Schedule, which is derived from the aggregate production plan, detailed end item forecasts, and existing inventory and orders (Vollman, 1992). Thus, raw materials and production capacity constraints are taken into account in the MPS to the extent that they were previously considered in the firm’s aggregate production plan––an infeasible MPS is only detected after a more detailed resource planning, later in the planning process. In a differentiated product portfolio, detailed item forecasts can be highly inaccurate, unexpected demand events are more frequent, and developing a feasible MPS is more challenging thus compromising the availability of reliable ATP information. It is not surprising to find several papers (eB2x 2000 ) , (Fordyce and Sullivan, 1999), (Lee and Billingtion, 1995), (Robinson and Dilts, 1999), and (Zweben, 1996) that discuss the need for advanced ATP systems, which provide order promising capabilities based on current capacity and inventory conditions within the firm’s supply chain. Interestingly, there are relatively few papers that address quantitative models for order promising. Taylor (1999) introduces a heuristic that keeps track of traditional ATP quantities to generate feasible due dates for order promising. Kilger (2000) also proposes a search heuristic to promise orders, motivated by yield management algorithms used in the airline industry (see below). Ervolina (2001) presents models

15

developed at IBM for resource allocation in a CTO production environment. Moses (2002) investigates highly scalable methods for real–time ATP that are applicable to discrete BTO environments facing dynamic order arrivals, focusing on production scheduling. At an operational level, Yongjin (2002) discusses the relationship between the performance of dynamic vehicle routing algorithms and online ordering in conditions of demand saturation––where demand exceeds service capacity. Chen and Ball (2001) provide mixed–integer programming (MIP) formulations for order promising and due–date quoting, taking into account existing inventories for raw materials, components and finished goods, production capacities, and a flexible bill of materials (BoM) environment (where customers can select different suppliers for the same raw material). Their models address a static situation, computing ATP

quantities for orders in a batching interval––a batching interval is the time window over which customer orders are collected before the ATP function is executed to schedule their production ––which is an input to their models. These models

maximize profit for a batching interval only, without consideration of future demands. We propose to address the stochastic and dynamic nature of the problem. More importantly, we address the ATP decision not just from customer response space, but also from product processing space that includes the product path analysis and uncertainty. In the following sections, we are going to review these topics.

2.2 Stochastic ATP - Uncertainty
In our analytical stochastic model we analyze how to allocate a resource to optimally promise customer orders when the future orders are uncertain. The problem

16

of allocating scarce capacity to customer orders to promise for future deliveries can be viewed from another perspective when customers’ sensitivity to lead times (or their due dates) is significantly different. In this environment, the firm may design a menu of price and lead–time combinations to segment the market (F.Chen, 2001). Once the menu is designed, the firm may use operational policies based on resource rationing to maximize profit. But ATP models and systems are essentially different from inventory models and systems. Inventory focused systems with different demand classes, different margins, and with different stock–rationing policies, for example Kaplan (1969) and Topkis (1968), have dynamic replenishment. There are n demand classes, and the penalty for not satisfying demand depends on the demand class (for example, one demand class pays a higher price, or can be met at a lower cost, thus having a higher priority for the firm). A significant stream of research exists on inventory rationing, depending on how assumptions of review period, demand distribution, and unsatisfied demand are handled (Cattani 2002), (Deshpande and Donhoue, 2001), and (Ha,1997). These models, however, do not consider capacity limitations in the replenishment decision, and therefore are fundamentally different from our dynamic ATP models. A body of literature closely related to this topic is newsvendor-like problems, which build inventory model to trade off order placing cost and holding cost. The work of Agrawal and Seshadri (1999) considers a single-period inventory model in which a risk-averse retailer faces uncertain customer demand and makes a purchasing-order-quantity and a selling-price decision with the objective of maximizing expected utility. They analyze how price affects the distribution of the

17

demand and in turn how the quantity to be committed should be determined. The paper by Cattani (2000) examines a single-period stochastic inventory problem where

N distinct kinds of demand can be satisfied with a single kind of product. But it
assumes that priority is hierarchical -- all demand for a higher priority product is met before meeting demand for the next lower priority product. Also, we address the demand pooling effect from economies of the scale, but in Cattani’s paper (2000) the author studies pooling effect from negatively correlated demands.

2.3 Coordinating Demand Fulfillment with Supply
In order to coordinate demand fulfillment with supply, we build models that identify profitable paths from different suppliers and to different sales subsidiaries. We provide both an analytical model and a global supply chain MIP model to investigate that issue, since they involve interactions among several individual factors as well as thousands of products and product locations. Of the three levels of planning in a supply chain––strategic, tactical, and operational (see, e.g., Vidal and Goetschalckx, 1997) ––our models address primarily tactical decisions, that is, production and distribution decisions that span a maximum of four months, such as production and transportation choices that maximize profits, subject to capacity and other constraints. Thus, our models neither address strategic decisions such as facility location nor operational decisions such as daily production scheduling at plants. There are several streams of literature that are relevant to our research: global supply chains, supply chain coordination mechanisms, inventory models with pricing, and

18

product line offering; we review each stream separately and discuss our study accordingly. Past literature addresses issues such as supplier–buyer coordination, for example, an ordering policy that minimizes supply chain costs for both parties (for a review see Thomas and Griffin, 1996)). Most of these models, however, are stylized extensions of the basic EOQ model or the Wagner–Within algorithm, focusing on a single product under independent and deterministic demand. Closely related is the literature on supply chain contracts (e.g. Bassok and Anupindi 1997, Urban 2000, Chen and Krass 2001, Serel, Dada and Moskowitz 2001), where, under uncertain demand for a single product, the buyer commits to a minimum cumulative procurement quantity

over a long–term planning horizon in exchange for price discounts. Our model
differs from the past literature in various aspects: in the analytical model, we consider demand curves that vary by sales locations; in the resource analysis of Toshiba global supply chain, we consider a situation where a firm orders a component periodically (there are no fixed ordering costs), however, the supplier offers a one–time price discount. The literature on inventory and pricing is also relevant and extensive––where the firm decides, in addition to order quantity, on the price of the product (which influences demand). For a review of single–period models, readers can see Petruzzi and Dada (1999), and for multi–period models, see Federgruen and Heching (1999), Chen, Federgruen and Zheng (2001) and Zhao and Wang (2002). Unlike our work, again, this stream of research assumes a single product with unlimited capacity.

19

Cohen and Lee (1989) argue that differences between supply–chain planning within a single country and for a global network include the existence of duties, tariffs, tax rates across countries, currency exchange rates, multiple transportation modes, and local content rules, among others. There is a considerable body of research in strategic production–distribution models for global supply chains, and the reader could refer to Thomas and Griffin (1996), Vidal and Goetschalckx (1997), Goetschalckx, Vidal and Dogan (2002) for detailed literature reviews. In addition to production and distribution decisions, this body of literature addresses the more strategic problem of network design, which is usually formulated as a mixed–integer programming (MIP) or a non–linear programming (NLP) model. The global nature of the problem may require careful modeling of transfer pricing (e.g. Vidal and Goetschalckx 2001), and exchange rates (e.g. Huchzermeier and Cohen 1996). Our research is also related to the product line offering question––which products are profitable and should be offered at each location. The literature on the design of a product line that maximizes profitability focuses primarily on marketing issues, such as the interactions of a set of products, given their relative utilities and prices, in the market place (for a review see Yano and Dobson, 1998). A few papers consider the manufacturing and/or inventory implications of product line breadth (e.g., Van Ryzin and Mahajan 1999, Smith and Agrawal 2000, Morgan, Daniels and Kouvelis 2001), such as variable production cost, holding and setup costs, but these papers do not consider the complex interactions of sourcing, manufacturing, and distribution in global, capacitated, supply chains, where the profitability of a product can be different, depending on its supplier or path in the supply chain and the location where

20

it is sold. Summarily, our analysis bridges the gap between sales’ product profitability and supplier’s variety.

2.4 Revenue Management
Finally, we discuss the revenue management literature, which clearly has strong relevance to our work. Most revenue management models assume fixed (or “almost fixed”) resource availability, (e.g., airline seats) and balance resource allocation among multiple demand classes (e.g., fair segment of price-sensitive customers). A common way to model the airline booking process is to model it as a sequential decision problem over a fixed time period, in which one decides whether each request for a ticket should be accepted or rejected. The classical example is that of customers traveling for leisure and those traveling on business. The former group typically books in advance and is more price-sensitive, whereas the latter behaves in the opposite way. Airline companies attempt to sell as many seats as possible to high-fare paying customers and at the same time avoid the potential loss resulting from unsold seats. In most cases, rejecting an early (and lower-fare) request saves the seat for a later (and higher-fare) booking, but at the same time that creates the risk of flying with empty seats. On the other hand, accepting early requests raises the percentage of occupation but creates the risk of rejecting a future high-fare request because of the constraints on capacity. The airline booking problem was first addressed by Littlewood in 1972, when he proposed what is now known as the “Littlewood Rule”. Roughly speaking, the rule — proposed for a two class model — says that low-fare bookings should be accepted as long as their revenue value exceeds the expected
21

revenue of future full fare bookings. This basic idea was subsequently extended to multiple classes (Belaboba, 1990). Later, it was shown that, under certain conditions, it is optimal to accept a request only if its fare level is higher or equal to the difference between the expected total revenues from the current time to the end when respectively rejecting and accepting the request. This rule immediately leads to the question “How to evaluate or approximate the expected total revenue from the current time until the end of booking?” However, the drawback of solving it as a sequential decision problem is also clear in that the booking policy is only locally optimized and it cannot guarantee global optimality. Glover et al. were perhaps the first to describe a network revenue management problem in airlines. By assuming that passenger demands are deterministic, they focus on the network aspects of the model (e.g., using network flow theory) rather than on the stochastic aspect of customer arrivals. Dror, Trudeau, and Ladany propose a similar network model, again with deterministic demand. The proposed improvements allow for cancellations, which often happens in the real booking process. Booking methods based on linear programming were thoroughly investigated by Williamson (Williamson 1992). The basic models take stochastic demand into account only through expected values, thus yielding a deterministic program that can be easily solved. The major drawback of the approach above is that it ignores any distributional information about the demand. Later many other industries also applied these techniques to control their perishable or even non-perishable assets. Weatherford and Bodily (1992) not only

22

propose the new name, Perishable-Asset Revenue Management (PARM), but also provide a comprehensive taxonomy and research overview of the field. They identify fourteen important elements for defining revenue management problems. Although most of these elements are airline-orientated, many ATP problems share the similar characteristics: particularly the last three modeling-related elements: bumping procedure (for handling “overbooking”), asset control mechanism (for resource reservation), and decision rule (for resource allocation). More recently, McGill and Van Ryzin (1999) classify over 190 research papers into four groups: 1) forecasting, 2) overbooking research, 3) seat inventory control, and 4) pricing. The papers in the third group are more relevant to ATP models discussed here. For example, Weatherford and Bodily (1995) present a generic multiple-class PARM allocation problem. They first study a simplified two-class problem without diversion. The problem assumes that there are two demand classes, full-price and discount, share the fixed available capacity of q0 units and that no full-price customer would pay less than their willingness to pay (i.e., the full price). The purpose is to determine the number of units that should be allocated to discount-price customers before the number of full-price customer is realized. The authors further extended the problem to allow diversion in the multiple-class setting. Sen and Zhang (1999) worked on a similar but more complicated problem by treating the initial availability as the decision variable and model the problem as a newsboy problem with multiple demand classes. To some degree, our work can be seen as the extension of Littlewood twoclass model under the special business setting. Yet due to the characteristics of ATP, the holding cost and customer service level, which normally is not in the scope of
23

revenue management, are taken into consideration and plays a critical role in our model (see Chapter 3.2).

24

Chapter 3 Pseudo Order Consideration in Available to Promise (ATP)
3.1 ATP Framework
The fundamental decisions ATP models must address are: 1) which orders to accept 2) the committed quantity for accepted orders. A sophisticated approach to carry out ATP functionality, introduced by Chen et al. (Chen et. al., 2002), is to employ large-scale mixed–integer-programming (MIP) models. Other researchers have also developed ATP models like allocated ATP (a-ATP) and capacity ATP (cATP) to support ATP decision-making process. This model-based approach for ATP execution can make effective use of resource flexibility and generate reliable orderpromising results. It is indeed efficient in some complicated business environments to support resource allocation and rescheduling. However, it usually takes more execution time to solve these models compared to conventional simple finishedproduct level ATP search results. As described in Chapter 1, real-time response is becoming a requirement based on customer service. Moreover, large numbers of customer orders with both accurate information and inaccurate configurations may come simultaneously in the e-business environment and/or large number of customer service channels. With this in mind, the MIP-based ATP mechanism may not be suitable because of its heavy computation. In contrast, analytical ATP models, which are based on simple rules and principles, can provide effective mechanisms for order promising solutions by trading

25

off multiple business objectives under resource constraints. Another benefit of using analytical models to solve ATP problems is their capability to tackle uncertainty. Unconfirmed customer orders can capture real-life customer’s inquiry and order cancellation, which is difficult for MIP types of models to handle. Undoubtedly, a further advantage of analytical models is that they offer generic solutions that don’t require extensive experiments. In this chapter, an analytical ATP model will be introduced for orderpromising and fulfillment decisions based on consideration of both profit and customer service. Instead of putting the customer service level in the constraints, we include it in the objective function. This reflects the trend toward pushing the service levels as high as possible. Other feature of this analytical ATP model is: A pseudo order with stochastic characteristics is considered with other confirmed orders to represent uncertain customer inquires and order cancellations. Since one pseudo order may have a higher profit margin but also uncertainty, it will have an impact on the commitment of the confirmed orders. It’s worth mentioning here that the uncertainties in ATP are mainly caused by three factors: demand, lead-time and raw material purchasing price. According to Weber’s survey (Weber and Current 2000), the effect of lead time is only 10-20% of effect of demand on a company’s total profit. On the other hand, the uncertainty from purchase price can be compromised by supply contracts as analyzed in our model. Thus we believe that using a pseudo order to incorporate the current orders and future orders together not only reflects price uncertainty in rapidly changing environment, but also captures the Achilles’ heel.

26

We shall see that our objective function is neither convex nor concave, but a d.c.

function, i.e. a function that can be represented as the difference of two convex
(concave) functions (Horst, 1993). The problem of maximizing a d.c. function under linear constraints is a nonconvex global optimization problem, which may have multiple local minima with substantially different values. Such multiextremal problems cannot be solved by standard methods of nonlinear programming which can at best locate a local minimum. Outer approximation methods along with branch and bound methods for finding a global minimum have been suggested in (Tuy, 1987). However, most of these methods are able to solve limited size problem instances. This should not be surprising, since the problem is known to be NP-hard, see e.g. (Pardalos, 1984). Therefore, simultaneous consideration of the uncertain demand makes the problem more general, and also more difficult. Fortunately, we have derived some rule-based solutions which are presented later in this section.

3.2 Problem Formulation & Model
The problem under consideration is a single period, single product, multi-order ATP model. This model consists of N confirmed customer orders, which are assumed to be deterministic, and one pseudo order, which is stochastic. The fundamental decision in the model is to determine promised order quantities for the confirmed customer orders and a reserved quantity for the pseudo order by considering both production impact and material limitations. The objective of including the pseudo order is to anticipate near-term future customer orders based on customer inquiry

27

information. The resultant model is a constrained non-linear stochastic programming problem. In the model, we use “total benefit” as the objective function, which is defined as the sum of weighted expected profit and customer service level. This reflects common practice in most order fulfillment and optimization processes. The relative weights of expected profit and customer service level can be adjusted in the model to reflect their importance in specific business settings. Meanwhile, we assume the committed quantity of the confirmed order can never exceed the requested quantity. For the stochastic pseudo order, holding costs may be incurred if the committed quantity is greater than the specified quantity. Below is the notation that will be used in the formulation.

3.2.1 Notation and Remarks on Function Properties
1,2, L , N } be the index of a set of the confirmed customer orders. For all Let I = {
i ? I , let qi and ri be the requested quantity and sales price of the confirmed order i ,

respectively; ? i is a weighted constant of the customer service level for the confirmed order i . For the pseudo order, let u be the pseudo order quantity, which is a random variable with known probability density function (PDF) as f (?) , and p , h the unit sales price and unit holding cost, respectively. For the order promising decision, we consider two kinds of resources: material availability and production capacity. Here, we assume the production cost is a convex function of the quantity produced. A convex production cost exhibits non-increasing

28

returns to scale. Basically, two factors can lead to that. One is overload or overtime. As an example, consider a factory with a regular workforce. If demand is beyond the capacity of the regular workforce, management has to employ overtime at a higher cost, and, if needed, it may subcontract production at an even higher cost. Convex production costs therefore are incurred and reflect diseconomies of scale (Galeotti and Maccini 2004). The other factor that leads to convex production cost is quality. Consumers have heterogeneous willingness to pay for quality, and the unit cost as a function of product quality is technology specific. Based on a distinct engineering principle, for a given production technology, the unit production cost tends to rise more rapidly as quality increases, and an increasing, convex cost function effectively captures such decreasing returns (Rochet and Chone 1998). We let m(?) be the

convex production cost function. Another kind of resource for order promising is material availability. Let g (?) be the material cost function, which is a concave function of material quantity used to reflect the economies of scale. Since ATP is used to serve customers, both total profit and customer service level should be considered vital performance criteria. In this paper, we employ order fillrate to model customer service level. The order fill-rate is defined as the quantity committed over the quantity requested for a customer order. Thus, the “customer service benefit” is a function of the order fill rate. Let s (?) be the customer service benefit function. It should be a concave function since most companies only pursue a certain higher level of customer service and penalize seriously very lower fill rates.

29

To formulate the committed quantity of the customer orders, we define the following decision variables. For all i ? I , let: X i equals the committed quantity of the confirmed order i ;

Y equals the committed quantity of the pseudo order.
3.2.2 Model Formulation

Based on the previous notation, the total revenue, which comes from both confirmed customer orders and the pseudo order, can be formulated as

? (r X
i i?I

i

) + pMin(Y , u ) ,

(3.1)

Where the first term is the confirmed order’s revenue, and the second term is the pseudo order’s revenue. We assume that the production cost and material cost are both measured in terms of a unit of product. Then, the material used is proportioned to the number of products committed to customers. The material cost and production cost is given as equation (3.2) and (3.3), respectively. g (? X i + Y )
i?I

(3.2) (3.3)

m( ? X i + Y )
i?I

The holding cost incurred by the pseudo order will be: h × Max(Y ? u,0) . Hence, the total expected profit is: EP = (3.4)

?r X
i i?I i?I

i

+ E ( pMin(Y , u ) - g (? X i + Y )
i?I

- m(? X i + Y ) - h × E ( Max(Y ? u,0))

(3.5)

30

Substituting the expected revenue and holding cost of the pseudo order into the PDF function in the above equation yields: EP =

? ri X i + p(? uf (u )du + ? Y f (u )du ) ? g (? X i + Y )
i?I 0 Y i?I

Y

?

? m(? X i + Y ) ? h ? (Y ? u ) f (u )du
i?I 0

Y

(3.6)

We note that the unit profit margin of an individual customer order may be negative like most discount sales. The reason to promise such orders is due to customer service level considerations. For the ith customer order, the customer service level, specifically, order fill rate, is X i q i . Hence, the total customer service benefit from all customer orders can be represented as below:

? ? s( q
i i?I

Xi
i

)

(3.7)

in which ? i is the customer service benefit weight. The value of ? i indicates the importance of customer service in comparison to profit for the ith customer order. Therefore, the objective function, which is defined as the “benefit function” since it includes both profit and customer service level, can be written as: TB =

?r X
i i?I

i

+ p ( ? uf (u )du + ? Y f (u )du ) ? g (? X i + Y )
0

Y

?

Y

i?I

i ? m(? X i + Y ) ? h ? (Y ? u ) f (u )du + ? ? i s ( X qi )

Y

(3.8)

i?I

0

i?I

subject to constraints described as following: qi ? X i ? 0 (3.9)

31

Xi ? 0
Y ?0

(3.10) (3.11)

Constraints (3.9), (3.10) and (3.11) are order-promising limitations and nonnegativity. They are similar to the widely accepted concept of “booking limits” or “protection level” in the revenue management area. Using the Lagrange method, the first order condition of the problem (3.8)-(3.11) can be given as follows for all i ? I .

ri ? g ' (? X i + Y ) ? m' (? X i + Y ) +
i?I i?I

?i
qi

i s ?( X qi ) ? ? i + ? i = 0

(3.12) (3.13) (3.14)
(3.15) (3.16) (3.17) (3.18) (3.19)

p(1 ? F (Y )) ? g ' (? X i + Y ) ? m' (? X i + Y ) ? hF (Y ) + w = 0
i?I i?I

?i ( q i ? X i ) = 0
? × Xi = 0 ?i ? 0 ?i ? 0
w×Y = 0

w?0

Where ?i , ? i and w are Lagrange multipliers for constraints (3.9), (3.10) and (3.11), respectively.
3.2.3 Model Analysis

32

As stated earlier the objective function is neither convex nor concave, but a d.c.

function. We present the following Theorem to assure the existence of optimal
solution for the problem (3.8)-(3.11). The proof is shown in Appendix 1.

Theorem 3.1: The objective function defined by (3.8) will be strictly concave on the

feasible region defined by problem (3.9)-(3.11) if,
M> Max( LU ,?(h + p ) K ) n +1 (3.20)

Where U = Max s ' ' ( x) ; L = Min
0 ? x ?1
1? i ? n

?i
q
2 i

; K = Min f ( y ) with y 0 = F ?1 (
0? y ? y0

p ) ; and p+h

M = Min [ g ' ' ( z ) + m' ' ( z )] with z 0 = y 0 + ? q i .
0? z ? z0
i =1

n

Remark 3.1: There exists one and only one optimal solution for the nonlinear

stochastic problem (3.8)-(3.11) if g (?) is a concave function, s (?) is a strictly concave function, and
g" (?) + m' ' (?) > 0

(3.21)

Proof: see Appendix 2.

One can observe that condition (3.21) is the special case of condition (3.20). Theorem 3.1 gives conditions for the existence of an optimal solution for the problem (3.8)–(3.11). Condition (3.20) basically states that the total cost (production cost plus material cost) function should be “convex enough” to compensate the concaveness of projected customer service level functions.

33

Obviously, the problem (3.8)–(3.11) is a multiple constrained newsvendor problem, which has been shown by Lau (1995) to be very difficult to solve, and not much work has been done on capacitated systems (Tayur, 1998). From an order promising point of view, we are more interested in the solution structure, rather than the solution itself, of the problem since the solution structure can provide insights and guidance for optimal order promising. We present the following Theorem to illustrate the structure of the solutions for the problem (3.8)-(3.11).

Theorem 3.2: For ?i, j ? I , if the following condition holds

[ Z ij ? s' (0)][ Z ij ? s' (1)] ? 0 , where Z ij is defined by Equation (3.23) with ? ij =
? ? 1, ? Z ij = ? ri ? r j , ? ? ? ij if? ij = 0 Otherwise

(3.22)

?i
qi

?

?j
qj

(3.23)

then there exist two points k , k ? ? I (k ' ? k ) , the optimal solutions of problem (3.8)(3.11) have the following structure:
Xi = 0

for all i < k ' ,

X i = qi for all i > k , X i* for all k ? ? i ? k is solved from the following equations: s' (
n n X i* q ) = i [ g ' (? X l + Y ) + m' (? X l + Y ) ? ri ] ?i qk l =1 l =1

(3.24)

34

F (Y ) =

p ? g ' (? X i + Y ) ? m' (? X i + Y )
i =1 i =1

n

n

p+h

(3.25)

Proof: see Appendix 3.

Theorem 3.2 provides the structure of the optimal solution. Based on this Theorem we can see that some customer orders should be one hundred percent committed, and some others should not committed at all when the conditions above are hold. The Z ij variable can be interpreted as the ratio between profit and customer service level, and Theorem 3.2 just states how a particular customer order promising scheme should be adopted depending on that ratio.
Remark 3.2: If the following condition holds, ?i ? I ,

?i
qi

s ' (0) + ri <

? i +1
qi +1

s ' (1) + ri +1 ,

(3.26)

then the optimal solution for problem (3.8)–(3.11) is as given in Theorem 3.2 with
k = k?.

Proof: see Appendix 4.

From this remark, one can observe that the optimal order commitment policy is characterized by a single key order, which represents a threshold between orders whose quantities are either 0 or qi . That quantity can be found by simply searching the linear solution space I. This gives us a simplified policy to make the optimal order promising for customers.

35

In real life situations, one company may only concentrate on a single criteria, choosing between customer service level and profit in order promising practices. In such case, condition (3.22)-(3.23) in Theorem 3.2 can be further simplified as the following:
?

In the case when product sales price plays a dominant role in order promising practice, condition (3.22)-(3.23) becomes (3.27) in deriving optimal solution to problem (3.8)-(3.11) in Theorem 3.2.
c ri +1 ? ri > , l
i where c = s ' (0) Max( ? qi ) l ? I

(3.27)

This solution structure can be easily checked and allows sales personnel to determine how the orders should be committed. When l = 1 , then there exists one and only one solution.
?

In the case when customer service level plays a dominant role in order promising practice, condition (3.22)-(3.23) would be simplified to (3.28) in Theorem 3.2.

?i
qi

l

? <

? i +1
qi +1

where ? =

s' (0) s' (1)

(3.28)

This is symmetric to condition (3.27).

Remark: Extension of Littlewood Model

We can also explain the work above from the perspective of revenue management. If we aggregate all the confirmed orders into one class (since they are all deterministic and order promising can be easily carried out by ordering the price

36

from high to low), and the problem represented by (3.8-3.11) becomes tradeoff of pseudo order class and confirmed order class. This is a two-class revenue management problem except we take customer service level into consideration. Thus, it can be viewed as an extension of the Littlewood two-class model under a special business setting. For example, if we manipulate this equation and treat Y as the optimal booking of (3.25)
p ? g ' (? X i + Y ) ? m' (? X i + Y )
i =1 i =1 n n

F (Y ) =

p+h

Then we have
h + g ' (? X i + Y ) + m' (? X i + Y )
i =1 i =1 n n

F (Y ) = 1 ?

p+h p2 (the Littlewood two-class model ) p1
n n

F (Y ) = 1 ?

where p 2 = h + g ' (? X i + Y ) + m' (? X i + Y ) and p1 = p + h
i =1 i =1

This result is in the form of Littlewood model. Belobaba (1987a) heuristically extends Littlewood’s rule to multiple fare classes and introduces the term Extended Marginal Seat Revenue (EMSR) for the general approach. The EMSR method does not produce optimal booking limits except in the two-fare case, and Robinson (1995) shows that, for more general demand distributions, the EMSR method can produce arbitrarily poor results. There has been very little published research on joint capacity allocation/pricing decisions in the revenue management context. Weatherford (1994)

37

presents a formulation of the simultaneous pricing/allocation decision that assumes normally distributed demands, and models mean demand as a linear function of price. The corresponding expressions for total revenue as a function of both price and allocation are extremely complex, but no structural results are obtained. However, our study not only obtains the structural results but also shows the impact of holding costs on both pseudo order and confirmed order quantities. More importantly, our model takes service level into consideration and allows not only 100% or 0% commitment but also portion commitment for each single order, unlike most discussions in the field of Airline or Yield Management. Formulas (3.26) (3.27) and (3.28) describe the rule of the accepting or rejecting order, which is determined by different order parameters such as price and SL. We can see the order’s “marginal benefit”, characterized in Remark 3.2, includes both price and SL; the “marginal revenue” as defined in the Littlewood rule and extended EMSR (Expected Marginal Seat Revenue) model in other Revenue Management research, includes only one parameter (price). ). Thus, our work can be viewed as extending prior revenue management research.

3.3 Implementation Rule
The analysis above provides decision makers the effective rule-based mechanism to implement and deploy critical decisions in real-time ATP systems. We summarize and illustrate it in Figure 3.1.

38

3.4 Experimental Study and Results
Our goal in developing analytical ATP model is to develop real-time strategies to support order-promising process. In order to gain strategic insight from the model numerically, we design and carry out the following experiments.

Figure 3.1 Pseudo order decision flow chart

Material economies of Scale.

No

Yes

Theorem 3.1 or g’’(.)+m’’(.) >0

Company Level

Yes

Order Level

Condition (3.22)

Customized Goal

Yes
Optimal Solution Structure

39

The cost and service level functions are assumed to be as following: s ( X i qi ) = ( X i qi + ? )? m(? X i + Y ) = k1 (? X i + Y )?
i? I i? I

(0 < ? < 1)
(0 < ? < 1) (? > 1)

g ((? X i + Y )) = k2 (? X i + Y )?
i? I i? I

f (Y ) = N ( ? , ? )
where the coefficients are given as below:

? i = 800 h = 3
? = 0 .9 l = 1 ,
k1 = 43 ? = 0.8 k 2 = 0.002

?=2

Four confirmed orders and one pseudo order are considered in this experiment with specifications as below:
r1 = 5 r2 = 8 r3 = 12 r4 = 18

p = 20

q1 = q 2 = q3 = q 4 = 400 , ? = 500 , ? = 70 By using MathCAD software, we have the solution: ( X 1 , X 2 , X 3 , X 4 , Y ) = (0, 0, 383.5, 400, 619.2) We can observe that the solution has the exact same structure as Theorem 3.2 states. Now let’s see some interesting findings from the experiments. 1) The pseudo order’s price sensitivity analysis:

40

Table 3.1: Effect of Pseudo Order’s Price Change
Y's Price 18 19 20 21 22 X1 0 0 0 0 0 X2 0 0 0 0 0 X3 387.2 385.8 383.5 381 380 X4 400 400 400 400 400 Y 614 617 619 621 623 SUM(X) 787.2 785.8 783.5 781 780 TB 6391 6860 7330 7802 8276

Figure 3.2 Committed Order as a function of pseudo order’s Price
Y's Price vs. SUM(X),Y

800 750
SUM(X), Y

700 650 600 18 19 20 Y's Price SUM(X) 21
Y

22

41

Figure 3.3 Total Benefit as a Function of Pseudo Order’s Price
Price vs. Total Benefit

8500 8000 7500 7000 6500 6000 5500 5000 18 19 20
Price

TB

21

22
TB

Figure 3.2 and Figure 3.3 shows the quantity of committed orders and total benefit as a function of pseudo order’s price respectively. We observe that as the pseudo order’s price increases, both the commitment quantity of the pseudo order and total benefit increases. These results imply that even with embedded uncertainty, the pseudo order’s price is still a decisive factor in determining the commitment quantity and total benefit. The commitment quantity of the total confirmed orders decreases as the pseudo order’s price increases. This is because they are sharing the same materials whose supply is tight. As in our case, the material cost function m(?) is a increasing convex function. So as the pseudo order’s price increases, the confirmed orders become less economically attractive. 2) Pseudo order’s variance sensitivity analysis:

42

Table 3.2 Effect of Pseudo Order’s Variance
Variance of Y 60 70 80 90 100 X1 0 0 0 0 0 X2 0 0 0 0 0 X3 394.5 383.5 373 360.6 352.5 X4 400 400 400 400 400 Y 606.7 619.2 630 641.5 651.5 SUM(X) 794.5 783.5 773 760.6 752.5 TB 7639 7330 7031 6741 6458

Figure 3.4 Committed Order Quantity as a Function of Pseudo Order’s Variance

850 800
SUM(X), Y

750 700 650 600 550 60 70
Y's Variance

80

90
SUM(X)

100
Y

43

Figure 3.5 Total Benefit as a Function of Pseudo Order’s Variance
Variance vs. Total Benefit

8000 7500 7000 6500 6000 60 70 80
Variance TB

90

100
TB

Figure 3.4 shows the committed order quantity as a function of the pseudo order’s variance. As variance increases, the total commitment quantity of confirmed order decreases, and the commitment quantity of the pseudo order increases. Note that the reason that more pseudo orders are committed is to offset the larger uncertainty, not to generate more profit, as evidenced by Figure 3.5, the total benefit decreases when variance and commitment quantity of the pseudo order increase. This is an interesting finding and reveals that uncertainty of the forecasted order will not only lead to the bullwhip effect via information distortion, but ultimately will hurt the efficiency of the supply chain in the form of excess raw material inventory, misguided production schedules, missed target orders, and poor customer service levels. Predictably, more and more companies are adopting the various cutting-edge technologies like Electronic Data Interchange (EDI) to improve their forecasting ability and mitigate such uncertainty. 3) The service level sensitivity analysis:

44

Table 3.3
Beta 750 760 770 780 790 800

Effects of Service Level on Order Promising
X1 0 0 0 0 0 0 X2 0 0 0 0 0 0 X3 352.4 359 366 373 380 387 X4 400 400 400 400 400 400 Y 614.8 614.7 614.6 614.4 614.3 614.2 SUM(X) 752.4 759 766 773 780 787 TB 6108 6164 6221 6278 6334 6391

Figure 3.6 Total Committed Quantity of Confirmed Order as a Function of Service Level

beta vs. SUM(X)

790 780 770
SUM(X)

760 750 740 730 750 760 770
beta

780

790

800

SUM(X)

Figure 3.6 shows that the total commitment quantity of confirmed order increases with service level’s weight. This reflects the significance of service level, which, like price, can be a driving force of commitment quantity. This observation allows us to integrate the parameter beta, the service level’s weight, into a company’s decision support system, providing a way to differentiate customers, e.g. in customer

45

relationship management (CRM) system. This provides an approach to offer superb service levels to the most valuable customers, develop strategies for unprofitable customers, and differentiate proactively the handling of different needs-based customer categories. Simply put, look for individual solutions rather than mass solutions.

3.5 Multiple Pseudo Orders in ATP
We further consider the case that captures more facts in real-life: a company not only has multiple confirmed orders, but also has multiple pseudo orders to fulfill. The stochastic characteristics associated with pseudo orders represent uncertain customer inquires and order cancellations. Particularly, we are interested to see what commitment decisions should be made if those pseudo orders’ profit function is concave of the order quantity. This is very true in most of the mass-production environments such as lot-by-lot manufacturing process, or pallet-by-pallet transportation among semi-final and final assembly factories.

3.5.1 Notation

The problem under consideration is a single period, single product, multi-order ATP model. The model consists of N confirmed customer orders, which are assumed to be deterministic; and M pseudo orders, which are stochastic. The fundamental decision in the model is to determine promised order quantities for the confirmed customer orders and the reserved quantity for the pseudo orders by considering both

46

production impact and material limitations. The objective of including the pseudo order is to anticipate near-term future customer orders based on customer inquiry information. The resultant model is a constrained non-linear stochastic programming problem. In the model, we use “total benefit” as the objective function, which is defined as the sum of weighted expected profit and customer service level. This reflects common practice in most order fulfillment and optimization processes. The relative weights of expected profit and customer service level can be adjusted in the model to reflect their importance in specific business settings. Meanwhile, we assume the committed quantity of the confirmed order can never exceed the requested quantity. For the stochastic pseudo orders, holding costs may be incurred if the committed quantities are greater than the specified quantity.
1,2, L , N } be the index of a set of the confirmed customer orders. For all Let I = {
i ? I , let qi and ri be the requested quantity and sales price of the confirmed order i ,

respectively; ? i a weighted constant of the customer service level for the confirmed order i . Let J = { 1,2, L , M } be the index of a set of the pseudo orders, p j and h j the unit sales price and unit holding cost of the pseudo order j, respectively. u is the pseudo order quantity, which is a random variable with known probability density function (PDF) as f (?) . As described in section 3.2.1, let g (?) be the concave material cost function and m(?) be the convex production cost function. Also, let s(?) be the customer service benefit function. We model this function as a concave function, which would reflect a

47

natural business perspective of increasing the per unit penalty for customer service deviations as the level of deviation increases (lower fill-rate are associated with higher deviations). To formulate the committed quantity of the customer orders, we define the following decision variables. For all i ? I , let: X i equals the committed quantity of the confirmed order i ; Y j equals the committed quantity of the pseudo order j ;

3.5.2 Model Formulation

Based on the previous notation, the total revenue, which comes from both confirmed customer orders and the pseudo orders, can be formulated as

? (r X
i i?I

i

) + ? p j E min(Yj, u ) ,
j?J

(3.29)

where the first term is the confirmed order’s revenue and the second term is the pseudo order’s revenue. We assume that the production cost and material cost are both measured in terms of a unit of product. Then, the material used is proportional to the number of products committed to customers. The material cost and production cost are given by equations (3.30) and (3.31), respectively.
g (? X i + ? Y j )
i?I j?J

(3.30) (3.31)

m( ? X i + ? Y j )
i?I j?J

The holding cost incurred by the pseudo orders will be:

48

? h E max(Y
j j?J

j

? u,0) .

(3.32)

Hence, the total expected profit is: TP = ? ri X i + ? p j E min(Y j , u ) ? g (? X i + ? Y j )
i?I j?J i?I j?J

? m(? X i + ? Y j ) ? ? h j E max(Y j ? u ,0)
i?I j?J j?J

(3.33)

We note that the unit profit margin of an individual customer order could be negative like most discount sales. The reason to promise such orders is due to customer service level considerations. For the ith customer order, the customer service level, specifically, order fill rate, is X i q . Hence, the total customer service benefit from all customer orders can be represented as below:

? ? s( q
i i?I

Xi
i

)

(3.34)

in which ? i is the customer service benefit weight. The value of ? i indicates the importance of customer service in comparison to profit for the ith customer order. Therefore, the objective function, which is defined as the “reward function” since it includes both profit and customer service level, can be written as: TB =

?r X + ? p
i i i?I j?J
i?I

j

E min(Y j , u ) ? g (? X i + ? Y j )
i?I j?J

i ? m(? X i + ? Y j ) ? ? h j E max(Y j ? u ,0) + ? ? i s ( X qi )

(3.35)

j?J

j?J

i?I

subject to constraints described as following: qi ? X i ? 0 Xi ? 0 (3.36) (3.37)

49

Yj ? 0

(3.38)

Constraints (3.36), (3.37) and (3.38) are order-promising limitations and nonnegativity. Using the Lagrange method, the first order condition of the problem (3.35)-(3.38) can be given as follows for all i ? I .
ri ? g ' (? X i + ? Y j ) ? m' (? X i + ? Y j ) +
i?I j?J i?I j?J

?i
qi

i s ?( X qi ) ? ? i + ? i = 0

(3.39) (3.40) (3.41) (3.42) (3.43) (3.44) (3.45) (3.46)

p j (1 ? F (Y j )) ? g ' (? X i + ? Y j ) ? m' (? X i + ? Y j ) ? h j F (Y j ) + w j = 0
i?I j?J i?I j?J

?i ( q i ? X i ) = 0

? × Xi = 0 ?i ? 0 ?i ? 0
wj ×Yj = 0 wj ? 0

Where ?i , ? i and w are Lagrange multipliers for constraints (3.36), (3.37) and (3.38), respectively.
3.5.3 Model Analysis Lemma 3.5.1 : The profit function of pseudo orders is concave in terms of order

quantity. Proof: See Appendix 5.

50

As stated earlier the objective function (3.35) is neither convex nor concave, but a d.c. function. We present the following Theorem to assure the existence of optimal solution for the problem (3.35)-(3.38). The proof is shown in Appendix 6.

Theorem 3.3: The objective function defined by (3.35) will be strictly concave on the

feasible region defined by problem (3.36)-(3.38) if, M > Max( LU ,? K ) n+m (3.47)

Where U = Max s ' ' ( x) ; L = Min
0 ? x ?1 1? i ? n

?i
q i2

; K = Min ( p j + h j ) K j ; Min K j = Min f ( y )
1? j ? m 0? y ? y j
n

with y j = F ?1 (

pj p j + hj

) ; and M = Min [ g ' ' ( z ) + m' ' ( z )] with z 0 = y 0 + ? q i .
0? z ? z0
i =1

Corollary 3.5.2: There exists one and only one optimal solution for the nonlinear

stochastic problem (3.35)-(3.38) if g (?) is a concave function, s (?) is a strictly concave function, and

g" (?) + m' ' (?) > 0
Proof: Obvious from Theorem 3.3.

(3.48)

One can observe that condition (3.48) is the special case of condition (3.47). Theorem 3.3 gives conditions for the existence of an optimal solution. Condition (3.48) basically states that the total cost (production cost plus material cost) function should be convex enough to compensate for the concavity of projected customer service level functions. This is very much true when material supply is tight.

51

From an order promising point of view, we are more interested in the solution structure rather than the solution itself, since the solution structure can provide insights and guidance for optimal order promising. We present the following Theorem to illustrate the structure of the solutions for the problem (3.35)-(3.38).

Theorem 3.4:

For confirmed orders, we have: Let Z ij = ri ? r j and ? ij =

?i
qi

?

?j
qj

. If we assume, ?i, j ? I ,

[ Z ij + ? ij s ' (0)][ Z ij + ? ij s ' (1)] ? 0 . Then there exist two points k , k ? ? I (k ' ? k ) such that the optimal solution
?0; ? X i = ? X i* ; ?q ; ? i ? i < k' ?k ' ? i ? k , Where X i* solves equation ?i > k

s' (

n m n m X i* q ) = i [ g ' (? X l + ? Y j ) + m' (? X l + ? Y j ) ? ri ] qk ?i l =1 j =1 l =1 j =1

(3.49)

For pseudo orders: we have: If we assume p1 ? p 2 ? L ? p m , then there exists k ?? ? I such that

? ?0; Yj = ? * ? ?Y j ;

? j < k ?? ? j ? k ??

, where Y j* solves equation

52

F (Y j* ) =

p j ? g ' (? X i + ? Yl ) ? m' (? X i + ? Yl )
i =1 l =1 i =1 l =1

n

m

n

m

p j + hj

(3.50)

Proof: see Appendix 7.

Theorem 3.4 provides a simple structure for the optimal solution for both confirmed orders and pseudo orders: some customer orders should be one hundred percent committed, and some others should not be committed at all. More specifically, confirmed orders’ commitment quantities are determined by K and K’ and pseudo orders’ commitment quantities are determined by K. Therefore we only need to compute K, K’ for confirmed orders and K” for pseudo orders. The interpretation for Z ij and ? i are the ratio between profit and customer service level, and Theorem 3.4 just states how a particular customer order promising scheme should be adopted depending on that ratio.

Corollary 3.5.3: For ?i ? I , if the following condition holds,

?i
qi

s ' (0) + ri <

? i +1
qi +1

s ' (1) + ri +1 ,

(3.51)

then the optimal solution for problem (3.35)–(3.38) would be as given in Theorem 3.4 with k = k ? . Moreover, we have k ?? ? arg min{ j | p j > g ?(q * ) + m?(q * ), j ? J } , where q* =

i = k +1

?q

n

i

.

53

Proof: see Appendix 8.

From this Corollary, one can observe that the optimal order commitment policy is decided by only one order commitment, and the others becomes either 0 or qi if the condition (3.51) is satisfied. That quantity can be found by simply searching the linear solution space I+J. It gives us a simplified policy to make the optimal order promising for customers.

3.6 Remarks
It is clear to see that demand uncertainty will not only lead to the bullwhip effect via information distortion, but ultimately will hurt the efficiency of a supply chain in the form of excess raw material inventory, misguided production schedules, missed target orders, and poor customer service level. Our analytical ATP model, as the core of customer-driven order management, indeed provides the manager with useful insights to conquer the resulting problems. Our models can provide guidance in both customer-service dominant and profit dominant business environments. There are two other types of approaches to pseudo orders promising: a traditional simulation-based approach and the more recent “scenario generation” approach. The simulation-based approach for promising the delivery of future orders is based on dynamic buffer adjustment coupled with forecasting the amount of buffer required. The primary objective of that is to frame the problem and suggest methods of analysis. Since this approach is simulation-based, optimal policies or solutions are

54

not obtained; simulation run times can also grow very large, e.g. from tens of minutes to several days to run one data set. Grant and Moses (2002) applied this approach but have to deal with an issue that inherently exists in this methodology: multiple yet conflicting performance measures, such as average lateness to promise, variance of the lateness to promise, and constraining the probability of a failed promise. We can see, all of these ask for more rather than less judgment calls from decision-makers. We should also note that systems based on simulation are not only
extremely sensitive to the distribution of the process times but also the congestion (arrivals of orders) of the system. Meixell and Chen develop a “scenario generation”

approach for efficiently incorporating uncertainty about future demand into an ATP system, so that demand scenarios with their associated probability distributions may be initially established and then continually revised. They first define demand as a random variable and then estimate the parameters of its distribution to define the uncertainty associated with demand for future order. From this distribution they derive scenario probabilities based on expected values and variances associated with the demand model for each order type, and finally to assign values for each of the branches of the scenario tree. They use “scenario generation” methods to convert stochastic pseudo orders into a pre-defined order types so that a dynamic linear model, instead of stochastic programming model, can be built. Further, in this approach, only a procedure for generating demand scenarios that is compatible with a stochastic mixed-integer-programming (MIP) model is developed, neither optimal solution structures nor optimal order promising policies are presented.

55

In contrast to these two approaches, we derive a new stochastic programmingbased customer order promising scheme along with efficient algorithms to implement it. These rule-based algorithms can be solved in real-time, and provide decision makers a clear answer: accept or reject certain orders based on a pre-set ratio of profit and customer service level without having to provide estimates of order arrivals.

56

Chapter 4 Coordinating Demand Fulfillment with Supply
In the current supply chains, suppliers, manufacturers and retailers are globally dispersed. Material cost differs among suppliers; production cost varies across manufacturers, and transportation cost depends on locations. On the other hand, as customers’ demand for the end products varies across the sales locations and is highly correlated to the sales price, the revenues are also different among paths. Consequently, the same product might be profitable along one path but unprofitable along another path. The choice of supplier or manufacturer can provide location– specific cost information and directly identify profitable and non–profitable products. In this section, we will answer two related questions: 1) How to coordinate demand fulfillment with supplies through the right path? 2) How to manage the available resource through the coordination of the raw material discount and end sales based on path analysis? First, we build and solve the analytical models, then provide the results of numerical experiments.

4.1 Analytical Model with Deterministic Demand Curve
4.1.1 Model & Formulation

The problem under consideration below is a single time period, single product, multiple-tier supply chain model. To simplify the deriving process, we assume the model consists of one manufacturer and n sales locations. The manufacturer buys raw

57

materials from the supplier, and after assembling, sells finished products to endcustomers across all sales locations. In most situations, the supplier will offer price discounts to buyers if their purchase quantity increases and therefore total purchase cost is concave. For every sales location we assume the demand curve is different. A supply chain path is from supplier to sales location. (See Figure 4.1). The fundamental decision in the model is to derive optimal purchase quantity of the raw materials and sales quantity of the end products in every location to maximize the profit of the company. This is a constrained non-linear programming problem.
Figure 4.1 Coordination in Supply Chain (1)
Supplier

f (Q)

Li

pi=ai-biQi

Sales Location

Below is the notation which will be used in the model formulation.
Notation:

Let I = { 1,2, L , N } be the index set of the sales locations. For all i ? I , let the path unit cost to be li and the demand curve be: pi = ai ? bi Qi (We use this format just to simplify the following deductive process. Of course, it is equivalent to the

58

widely used format: Qi =

ai 1 ? pi ). pi is the sales price. Let the purchase cost bi bi

function be f (?). It is a concave function of the purchased quantity. So the decision variables are:

Q : the raw material purchase quantity.
Qi : the quantity of finished product sold at location i.

We assume one raw material unit per finished product unit, so that:

?Q
i =1

n

i

=Q

Model Formulation:

Based on the previous notation, the total revenue can be formulated as

? p Q =? (a
i i i?I i?I

i

? bi Qi )Qi

(4.1)

The raw material purchase cost and the path cost is given by expressions (4.2) and (4.3), respectively.
f (? Qi )
i?I

(4.2) (4.3)

?l Q
i i?I

i

Hence, the total profit is:

TP =

? (a
i?I

i

? bi Qi )Qi - f (? Qi ) i?I

?l Q
i i?I

i

(4.4)

We reorganize the formula and have the objective function:

59

TP =

? (a
i?I

i

? li ? bi Qi )Qi - f (? Qi )
i?I

(4.5)

Subject to constraints described as following:
Qi ? 0

(4.6)

From the Lagrange method, the first order condition of problem (4.5)-(4.6) can be given as following: for all i ? I .

?i = ai ? li ? 2bi Qi ? f ' (? Qi ) ? 0
i?I

?i Qi = 0
Where ?i is Lagrange multiplier for constraint (4.6).

(4.7)

4.1.2 Model Analysis Theorem 4.1: Let (ai ? l i ) be in increasing order of i. There exists a point k ? I , the

optimal solutions of problem (4.5)–(4.6) have the following structure:
Qi = 0
a i ? l i ? f ' (? Qi )
i?I

for all i ? k ,

Qi =

2bi

for all i > k .

Proof: see the Appendix 9. We have observed that f ' (? Qi ) is the marginal purchase cost, and we can
i?I

clearly see how the material price discounts would have an impact on the end sales Qi based on the numerical experiments we conducted in Chapter 4.4. However, Theorem 4.1 reveals in close-form how the discount offered by a supplier can have an impact

60

on profit and purchased quantity. It will not only change the Qi ( i > k ), but it will also change the Qi ( i ? k ). It gives us a formula, i.e. ai ? li ? f ' (? Qi ) , to identify
i?I

"profitable" and "non-profitable" demand, and only "profitable" demand will be fulfilled and "non-profitable" demand will not be fulfilled -- a conclusion which is different from the one shown in Chapter 4.4 because we don't have a "minimum fill rate" constraint in this case. Furthermore, since material cost function f (?) is a

concave function, i. e. f ' (?) < 0 , the following applies: when the supplier offers discount, f ' (?) decreases and Qi increases, which is in line with the one in Chapter 4.4. The other mechanism being investigated in Chapter 4.4 is offering of discount on profitable products, or the equivalent of the following: If we reduce bi , then from Theorem 4.1 Qi will increase, which is also in line with conclusion of Chapter 4.4: Offering discount for profitable products will increase total purchased component quantity.

From this Theorem we can have the following steps to derive the optimal Qi . Algorithm: Let ?1 =
n ai ? li 1 , ? = . ? ? 2 i = k +1 2bi i = k +1 2bi n

Beginning Step 1: Let k = 1 Step 2: Derive the MC by solving the equation: f ' ( ?1 ? ? 2 ? MC ) = MC

61

Step 3: Check if the MC satisfy the assumption:
a1 ? l1 ? ... ? a k ? l k ? MC < a k +1 ? l k +1 ? ... ? a n ? l n

(4.8)

If Yes, go to step 4; otherwise go to step 5. Step 4: Check if the Hessian matrix at this point is negative definitive. If so, store this solution, go to step 5; otherwise, go to step 5 without storing solution. Step 5: If k < n , let k = k + 1 go to step 2. If k = n , go to step 6; Step 6: Choose the solution with maximal objective value from the stored, this is the optimal solution. End. (See Figure 4.2) Since f (?) is a concave function, we have f ' ' (?) ? 0 . Further, if the f (?) satisfy:
? 2 1 ? i?I bi < f ' ' (?) < 0 , we can have the following conclusion:

Remark 4.1: Let the (ai ? l i ) in increasing order of i. If the f (?) satisfy:
? 2 1 ? i?I bi < f ' ' (?) < 0 , there exists one and only one point k ? I , the optimal solutions

of problem (5)–(6) have the following structure:
Qi = 0
a i ? l i ? f ' (? Qi )
i?I

for all i ? k ,

Qi =

2bi

for all i > k .

Proof: see the Appendix 10. With Remark 4.1, we can simplify the solution process as following:

62

Beginning Step 1: Let k = 1 Step 2: Derive the MC by solving the equation: f ' ( ?1 ? ? 2 ? MC ) = MC , Step 3: Check if the MC satisfies the assumption:
a1 ? l1 ? ... ? a k ? l k ? MC < a k +1 ? l k +1 ? ... ? a n ? l n

If it does, stop, this is the optimal solution. Otherwise, let k = k + 1 go to step 2. End. (See Figure 4.3)
4.1.3 Numerical Experiments

Based on the analysis above, the most important factors in our model that determine the structure of optimal solution to fulfill demand would be: maximal sales price ai and unit transportation cost li . The price elasticity bi and purchase cost function f (?) would only influence the optimal value of demand fulfillment. In the following experiments we show how the structure of optimal solution is in line with Theorem 4.1; then we run multiple scenarios for some sensitivity analysis. Assuming we have 3 sales locations, i.e. i =3. With demand curve pi = ai ? bi Qi , the coefficients are shown in Table 4.1 below:
Table 4.1: Experiment Coefficients Parameters/ Locations Maximal Price: ( ai ) Price Elasticity ( bi ) Transportation Cost ( li ) 1 4 1.8 0.12 2 6 2 0.11 3 10 1.9 0.15

63

Thus (ai ? l i ) is in increasing order of

i.

Further, let the purchase cost

function f (Q) = 6 Q , which is an increasing concave function.

Figure 4.2 Coordination Algorithm (1)
Let k = 1

f ' (?1 ? ?2 ? MC) = MC

k = k +1
(4.8) Hold

No

Yes
Hessian Matrix is

No

negative definitive
Yes

Store the solution

No

k=n
Yes

OUTPUT Optimal Solution Figure 4.2 Coordination Algorithm (1) (the one with the maximal objective value)

64

Figure 4.3 Coordination Algorithm (2)
Let k = 1

f ' (?1 ? ?2 ? MC) = MC

k = k +1
(4.8) Hold
Yes

No No

k=n

Yes

OUTPUT Optimal Solution

First, we solve the problem analytically (4.5 - 4.6) by Theorem 4.1 and algorithm shown in Figure 4.2, the optimal solution is: (Q1 , Q2, Q3 ) = (0, 0, 2.05), the objective function value (total profit) is 3.62. Now, let us examine the numerical experiments to determine where our optimal solution stands in terms of the possible different solution structures: (0, 0, Q3 ), (0, Q2 , Q3 ), ( Q1 ,0, Q3 ), ( Q1 , Q2 , Q3 ), which are shown in Table 4.2 - 4.5 and the corresponding Figure 4.4 – 4.7 , respectively. Figure 4.4 clearly shows that as Q3 increases, total profit increases, but it begins to decline beyond a certain value: 2.05. This occurs because, in order to have more demand, the charged price has to be at a lower level. Note that after sales quantity reaches 2.05, the additional profit can not offset the lost margin and increased transportation cost anymore. Figure 4.5 – 4.7 show that if we increase Q1 or Q2 , or both from 0 to some value without following Theorem 4.1, the total profit would indeed decrease. The
65

reason that the additional fulfilled demand quantity can not overcome the loss in revenues is because (ai ? li ) is too low to be profitable. In fact, under such situation the more Q1 or Q2 sold, the less profit generated. Overall, these experiments numerically demonstrate the implication of the solution structure in the Theorem and Remark 4.1. Price elasticity and purchase cost function only affect the optimal solution’s value. However, more importantly, the factors that determine the fundamental structure of the optimal solution is the maximum price ai and transportation unit cost li . Now we examine some other interesting findings from further experiments. First, we
Table 4.2: Q3 vs. Total Profit

Q1
0 0 0 0 0 0

Q2
0 0 0 0 0 0 0

Q3
1.9 1.95 2 2.05 2.1 2.15 2.2

SUM (Q)
1.9 1.95 2 2.05 2.1 2.15 2.2

TP
3.59 3.6 3.61 3.62 3.61 3.6 3.57

Figure 4.4: Q3 vs. Total Profit
3.8 3.7 3.6 3.5 1.9 1.95 2 2.05
Q3 TP

Q3 vs. TP

2.1

2.15

2.2

66

Table 4.3: Q2 , Q3 vs. Total Profit:
Q1
0 0 0 0 0 0 0

Q2
0 1.05 1.1 1.15 1.2 1.25 1.3

Q3
2.05 2.05 2.05 2.05 2.05 2.05 2.05

SUM (Q)
2.05 3.1 3.15 3.2 3.25 3.3 3.35

TP
3.62 3.52 3.42 3.3 3.18 1.58 1.45

Figure 4.5: Q2 , Q3 vs. Total Profit
(Q2,Q3) vs. TP
4 3 2 1 0 2.05 3.1 3.15 3.2 Q2+Q3 TP 3.25 3.3 3.35

Table 4.4: Q1 , Q3 vs. Total Profit:

Q1
0 1.05 1.1 1.15 1.2 1.25 1.3

Q2
0 0 0 0 0 0 0

Q3
2.05 2.05 2.05 2.05 2.05 2.05 2.05

SUM (Q)
2.05 3.1 3.15 3.2 3.25 3.3 3.35

TP
3.62 1.63 1.45 1.26 1.06 0.85 0.63

67

Figure 4.6: Q1 , Q3 vs. Total Profit
4 3 2 1 0 2.05 3.1 3.15 3.2
Q1+Q3 TP

(Q1,Q3) vs. TP

3.25

3.3

3.35

Table 4.5: Q1 , Q2 , Q3 vs. Total Profit:

Q1
0 1.05 1.1 1.15 1.2 1.25 1.3

Q2
0 1 1 1 1 1 1

Q3
2.05 2.05 2.05 2.05 2.05 2.05 2.05

SUM (Q)
2.05 4.1 4.15 4.2 4.25 4.3 4.35

TP
3.62 3.61 3.61 3.6 3.59 3.56 3.52

Figure 4.7: Q1 , Q2 , Q3 vs. Total Profit
(Q1,Q2,Q3) vs. TP
3.7 3.6 3.5 3.4 2.05 4.1 4.15 4.2 Q1+Q2+Q3 TP 4.25 4.3 4.35

examine the effect of transportation cost l3 . Not surprisingly, with transportation cost increasing, both fulfilled demand quantity and total profit decrease, and its effect on

68

profit is more significant since it results in both a decrease in fulfilled demand and an increase in cost, as shown in Figure 4.8. Note that as long as the (ai ? li ) is increasing in the order of i, and the condition in the Remark 4.1 holds, the solution still has Q1 and Q2 at 0.
Table 4.6: Effect of Transportation Cost l3 ( l1 = 0.12, l 2 = 0.11)

l3
0.06 0.08 0.1 0.12 0.15 0.18 0.2 0.22

Q1
0 0 0 0 0 0 0 0

Q2
0 0 0 0 0 0 0 0

Q3
2.068 2.062 2.058 2.053 2.05 2.03 2.028 2.025

TP
3.8 3.76 3.72 3.69 3.62 3.57 3.51 3.44

Figure 4.8: Effect of Transportation Cost l3
4 3.5 3 2.5 2 1.5 0.06 0.08 0.1 0.12 L3
TP Q

0.15

0.18

0.2

0.22

Next, we look at the impact of the maximal unit price that a firm can charge on demand quantity and profit by varying a3 from 8 to 14. These are shown in Table 4.7 and Figure 4.9. They display both Q and total profit increase as maximal price a3 increases. This is not surprising; as matter of the fact, maximal price reflects the

69

customers’ desire and affordability of the product. Total profit increases more significantly than Q due to the dual effects of raising both Q and price. Again, we should notice that (a3 ? l 3 ) in this scenario is always larger than

(a1 ? l1 ) and (a 2 ? l 2 ) , and in optimal solutions Q1 and Q2 remain at 0, which is exactly

as Theorem 4.1 points out.
Table 4.7: Effect of Maximal Price a3 ( a1 = 4, a 2 = 6 )

a3
8 9 10 11 12 13 14

Q1
0 0 0 0 0 0 0

Q2
0 0 0 0 0 0 0

Q3
1.4 1.8 2.05 2.35 2.65 2.96 3.2

TP
0.17 1.72 3.62 5.81 8.29 11.07 14.13

Figure 4.9: Effect of Maximal Price a3
15 10 5 0 8 9 10 11 A3 Q TP 12 13 14

Table 4.8 and Figure 4.10 show the optimal solution of fulfilled demand quantity and total profit as a function of price elasticity. For various levels as price elasticity,
b3 , increases, both Q and total profit decrease, and profit decreases more steeply due

to the same dual effects as mentioned in maximal price discussion. Further, if we

70

compare Figure 4.8, 4.9 and 4.10, we can see that price-related coefficients, i. e., maximal price a3 and price elasticity b3 , have a more significant effect on profit than transportation cost, does as evidenced by the “flat” curve in Figure 4.8 vs. “steep” curve in Figure 4.9 and 4.10. This occurs because when end customers are buying something, they are always more sensitive to visible price than transparent transportation cost. In addition, in most cases, price accounts for a larger portion of the cost than the transportation cost itself. Thus the firm should carefully determine the quantity of demand to fulfill as demonstrated in Theorem 4.1 and shown in the results above. Fulfilling more demand doesn’t necessarily mean more profit is secured. Sometimes it may even hurt the business.
Table 4.8: Effect of Price Elasticity’s b3 ( b1 =1.8, b2 = 2 )

b3
1.6 1.7 1.8 1.9 2 2.1 2.2

Q1
0 0 0 0 0 0 0

Q2
0 0 0 0 0 0 0

Q3
2.49 2.32 2.18 2.05 1.92 1.82 1.72

TP
5.14 4.56 4.06 3.62 3.22 2.88 2.56

Figure 4.10: Effect of Price Elasticity’s b3
5.5 5 4.5 4 3.5 3 2.5 2 1.5 1.6 1.7 1.8 1.9 B3 Q TP 2 2.1 2.2

71

4.2 Analytical Model with Stochastic Demand
In Section 4.1, we assumed that the demand is price-dependent with different demand curves at each sales location. Therefore no inventory cost is considered. This scenario is true as stated in the theory of monopolistic competition (Chamberlain and Robinson 1954), where the demand curve of a monopoly firm is downward sloping and the inventory will not be an issue when price can be used as a decision tool. However, if one company did not posses the monopolistic power, e.g., under fierce global competition, it has to face the uncertainty of customer demand. In such cases, inventory holding to support sales becomes significant. Effectively managing and minimizing the inventory holding cost while satisfying customer demands could provide a competitive advantage to any firm in the marketplace. Besides the inventory cost, we also take into account the lost sales (supply shortage) penalty that would be incurred when demand can not be met. We further assume that the uncertain demand on each path can have different distributions. The fundamental decision in this constrained stochastic model is still to derive an optimal purchase quantity for the raw materials and a sales quantity for the end products in every location to maximize the profit of the company. A body of literature closely related to this topic is newsvendor-like problems, which are to find the order quantity which maximizes the expected profit in a singleperiod probabilistic demand framework. By assuming that the retailers’ demand obeys a normal or lognormal distribution and that the retailers place orders according to the Newsvendor Rule, researchers derive the necessary and sufficient conditions

72

for the optimal solution of production size (Khouja, 2000; Z.Weng 2003). Our model differs from theirs on: 1) It is more generalized and all demands are equally treated; (Unlike Cattni (2000)’s paper, we don’t need hierarchical priority in the model.) 2) it can be applied to any unknown-but-bounded disturbances; 3) with both manufacturers and retailers in the picture, it is a constrained stochastic problem. We also would like to point out that this model is different from a multi-inventory system. Most papers that have discussed multi-inventory systems have assumed that the retailer’s ordering policy is only related to the demand and supplier but unrelated to the other retailers. Some supply-chain and inventory models use the following twoechelon symmetric-information and deterministic gaming structure: a

‘‘manufacturer’’ wholesales a product to a ‘‘retailer,’’ who in turn retails it to the consumer. The retail market demand varies either with the retail price according to a deterministic ‘‘demand function’’ that is known to both the manufacturer and the retailer or with the known distribution. The manufacturer is a Stackelberg leader and the retailer is a Stackelberg follower (from Franco Blanchini etc. 2004). Thus they completely separate the “market channel” from “manufacturing channel” to make decisions. Under the approach described, in the next section, we derive an optimal promising policy from a centralized supply chain’s point, in which the company considers the “market channel” and “manufacturing channel” together as a path.

4.2.1 Model & Formulation

The problem under consideration is a single time period, single product, twoechelon supply chain system. As before, we assume there is one supplier and n sales

73

locations. The manufacturer buys raw materials from the supplier, and after assembling, sells finished products to end-customers via sales locations. In practical situations, the supplier usually offers price discounts to manufactures if their purchase quantity increases and therefore total purchase cost is concave (the cost function is denoted by g (?) ). The supply chain path is the path from one manufacturer to sales
1,2, L , N } be the index set of the sales locations. For location (Figure 4.11). Let I = {

all i ? I , let the path unit cost be li and the demand of each market as u i , which is a random variable with known probability density function (PDF) denoted by f (?) and
c. d. f. denoted by F (?) . pi is the sales price in each market. We assume one raw

material unit per finished product unit. Our decision variables are:

Q : the raw material purchase quantity.
S i : the quantity of finished product sold at location i.
Figure 4.11 Supply Chain Path (2)
Supplier g (Q)

li

Price: pi Demand: ui Holding : hi

Sales Location

74

Formulation:

Revenue:

? pi min(S i , ui ) =
i

? pi ( ? ui f (ui )dui + ? S i f (ui )dui )
i

Si

?

0

Si

Purchase Cost: g (Q) , Transportation Cost:

g (?) is a concave function

?l S
i i

i

Holding Costs on sales location:

? hi × max(S i ? ui ,0) =
i

? [hi ? (S i ? ui ) f (u i )dui ]
i

Si

0

Shortage Cost:

? d i × max(ui ? S i ,0) = ?[di ? (ui ? Si ) f (ui )dui ]
i i Si

?i

We would like to: Maximize: Profit = Revenue - Purchase Cost - Transportation Cost - Holding Costs – Shortage cost So the problem now is: Max:

? p (? u
i i

Si

i

0

f (u i )du i + ? S i f (u i )du i ) - g (Q) Si

?

? l S - ? [h ? ( S
i i i i i

Si

i

? u i ) f (u i )du i ] -

0

? [d i ? (ui ? S i ) f (ui )dui ]
i Si

?

s. t.

Q ? ? Si ? 0
i

Si ? 0 , Q?0

It is equivalent to:

75

Min: s. t.

? ( pi + hi + d i ) ? (S i ? u i ) f (ui )du i + g (? S i ) ? ? ( pi + d i ? li )S i
i

Si

(4.9) (4.10)

0

i

i

Si ? 0

4.2.2 Model Analysis

From the Lagrange method, the first order condition of this problem can be given as the following: Kuhn-Tucker condition for all i ? I .

?i = ( pi + d i + hi )Wi ' ( S i ) + g ' (? S i ) ? ( pi + d i ? li ) ? 0
i

(4.11)

?i S i = 0 ?i ? 0
Wi ( S i ) is defined by Ei [ S i ? u i ] = ? ( S i ? u i ) f (u i )du i . Assume Wi (?) is everywhere
+
0
Si

differentiable, we have: Wi ' ( S i ) = F ( S i ) . It’s easy to show: Wi (0) = Wi (0) = 0 , and
Wi (?) is convex (see Appendix 11).

'

Theorem 4.2: Let the ( pi + d i ? li ) in increasing order of i. There exists a point
k ? I such that the optimal solutions of problem (4.9)–(4.10) has the following

structure:
Si = 0
( p i + d i ? l i ) ? g ' (? S i )
i?I

for all i ? k ,

(4.12)

F (S i ) =

p i + d i + hi

for all i > k .

(4.13)

Proof: see Appendix 12.
76

We can further show that an upper bound exists for Q (see Appendix 13), which reveals that, even with a discount on raw material, it would not be foolish to limit the total quantity purchased. Unsurprisingly, the transportation cost li , shortage penalty
d i and sales price pi have an impact on the decision. Also, we notice that g ' (? S i )
i?I

can be interpreted as the marginal purchase cost. From Theorem 4.2 we can easily see it is a decisive factor of the end sales S i . Another interesting finding is that inventory cost hi has no effect at the first stage where we determine whether the finished product should be sold at one location or not. It only plays a role at the second stage where we determine the sold quantity S i (for i > k ).
4.2.3 Algorithm:

BEGIN: Step 1: Let k = 1 Step 2: Derive the S i by solving the equations (4.12)-(4.13). Step 3: Check if the g ' (? S i ) satisfy the assumption:
i

p1 + d1 ? l1 ? ... ? p k + d k ? l k ? g ' (? S i ) < p k +1 + d k =1 ? l k +1 ? ... ? p n + d n ? l n
i

(4.14) If so, go to step 4; otherwise go to step 5. Step 4: Check if the Hessian matrix at this point is negative definitive. If so, store this solution, go to step 5; otherwise, go to step 5 without storing. Step 5: If k < n , let k = k + 1 go to step 2. If k = n , go to step 6;

77

Step 6: Choose the solution with maximal objective value from the stored, this is the optimal solution. END. (see the Figure 4.12)
Remark 4.2: Let the ( pi + d i ? li ) be increasing i. If 1 > UF × UG ? i ( p i + hi + d i )

where UF and UG is the upper bound of f (?) and ? g ' ' (?) respectively, then there exists a unique k ? I for which the optimal solutions of problem (4.9) – (4.11) have the following structure:
Si = 0
( p i + d i ? l i ) ? g ' (? S i )
i?I

for all i ? k ,

F (S i ) =

p i + d i + hi

for all i > k .

Proof: see the Appendix 14. Using Remark 4.2, we can state the following simplified algorithm:
Algorithm 2:

BEGIN: Step 1: Let k = 1 Step 2: Derive the S i by solving the equations (4.12)-(4.13). Step 3: Check if the g ' (? S i ) satisfy the assumption:
i

p1 + d1 ? l1 ? ... ? pk + d k ? lk ? g ' (? Si ) < pk +1 + d k =1 ? lk +1 ? ... ? pn + d n ? ln
i

78

If so, Stop; Otherwise, let k = k + 1 go to step 2. END. (see the Figure 4.13)

Let k = 1

Solve Equations (4.12-4.13)

k = k +1

pk + d k ? l k ? g ( ? S i )
No Yes

Hessian

Matrix

is

No

negative definitive

Yes

Store the solution

No

k=n
Yes

OUTPUT Optimal Solution (the one with the maximal objective value)
function

Figure 4.12. Coordination Algorithm (3)

79

Let k = 1

Solve Equations (4.12-4.13)

k = k +1

(4.14) Hold
No No Yes

Yes

k=n

OUTPUT Optimal Solution

Figure 4.13. Coordination Algorithm (4)

4.2.4 Numerical Experiments

Based on the analysis above, the most important factors in our model that determine the structure of optimal solution to fulfill demand are: sales price pi , unit transportation cost li , and unit shortage cost d i . Holding cost hi and purchase cost function f (?) would only influence the optimal value but not the structure. In the following experiments we show how the structure of optimal solution is in line with Theorem 4.2; then we run multiple scenarios for some sensitivity analysis. Assuming we have 3 sales locations, i.e. i = 3. The demand of each market u i is a random variable with known probability density function (PDF), f (?) , and c. d. f. , F (?) . Assume f (?) is the widely-used lognormal distribution with E( u ) = 300, and std( u ) =

80

60. In addition, let the purchase cost function f (Q) = 6 Q , which is an increasing concave function. The other coefficients are shown in Table 4.9 below:
Table 4.9: Experiment Coefficients Parameters/ Locations Sales Price: ( pi ) Transportation Cost ( li ) Holding Cost ( hi ) Shortage Cost ( d i )

1 4 0.7 0.2 0.5

2 8 0.8 0.25 0.8

3 10 0.9 0.4 1

First, we solve the problem analytically (4.9 - 4.11) by Theorem 4.2 and using the algorithm shown in Figure 4.12 yields the optimal solution: ( S1 , S 2, S 3 ) = (0, 0, 198), with objective function value (total profit) 762. Like the experiments shown in Table 4.2 - 4.5, the solution structure has exactly the same form as Theorem 4.2 and Remark 4.2, since in our case ( pi + d i ? l i ) is in increasing order of i. In addition, it is confirmed numerically that the purchase cost function and holding cost only affect the optimal solution’s value. Most importantly, the sales price pi , transportation cost li and shortage cost d i determine the fundamental structure of the optimal solution. Now let us examine some interesting findings from the following experiments. First, we investigate the effect of varying the demand mean. Unsurprisingly, as the mean of demand increases, both fulfilled demand quantity and total profit increases, its effect on profit is more significant since it results in both sales increase and purchase cost decrease due to economies of scale, as shown in Figure 4.14. Note that as long as ( pi + d i ? l i ) is increasing in the order of i, and the condition of Remark 4.2 holds. S1 and S 2 will remain at 0.
81

Table 4.10: Effect of Demand’s Mean
AVG
297 300 317 328 380

(STD = 60)

S1
0 0 0 0 0

S2
0 0 0 0 0

S3
170 198 205 224 296

Q = sum( S i )
170 198 205 224 296

TP
684 762 886 932 1013

Figure 4.14: Effect of Demand’s Mean
1300 1100 900 700 500 300 100 297 300 317
TP

328
S3, Q

380

Mean of Demand

Next, we look at the impact of demand’s standard deviation on fulfilled demand quantity and profit by varying std from 30 to 150. These are shown in Table 4.11 and Figure 4.15. They show that the fulfilled demand quantity increases, while the total profit decreases as the standard deviation of demand increases. This can be explained as follows: large variance in demand patterns create havoc for a company, which in turn must sell more quantity to offset it. Even so, the generated revenues are still not enough to overcome the loss due to overstocking (holding cost) and lost sales (shortage cost). Thus, accurate demand forecasts are a business imperative. After all, the better you can match supply with actual demand, the more streamlined your business operations will be. Thus, a company needs the capability to track forecast accuracy in real-time, allowing management to take immediate action to eliminate forecast errors and improve the overall demand forecasting process.

82

Again, we should notice that ( p3 + d 3 ? l 3 ) in this scenario is always larger than
( p1 + d 1 ? l1 ) and ( p 2 + d 2 ? l 2 ) , and in optimal solutions S1 and S 2 stay at 0, which is

exactly what Theorem 4.1 points out.
Table 4.11: Effect of Demand’s Stand Deviation (AVG = 300)
STD
30 60 90 120 150

S1
0 0 0 0 0

S2
0 0 0 0 0

S3
191 198 216 237 263

Q = sum( S i )
191 198 216 237 263

TP
778 762 746 722 688

Figure 4.15: Effect of Demand’s Stand Deviation
900 700 500 300 100 30 60 90 STD of Demand TP S3,Q 120 150

For various levels of transportation cost l3 , as l3 increases (Table 4.12), both Q and total profit monotonically decrease. This result is obvious as the company is increasingly more costly now, profits decrease more steeply (Figure 4.16) due to the dual effects from increased cost and reduced fulfilled quantity. We also notice that
S 2 changes from 0 to 142 when l3 increases to 3.2. This should not surprise us as,

according to Theorem 4.2, the break point for ( p3 + d 3 ? l3 ) in this scenario is l3 = 3.0. In addition, even the optimal value of S 2 changes, the solution structure still observes the formula (4.12) and (4.13).
83

Table 4.12: Effect of Transportation Cost ( l1 = 0.7, l 2 = 0.8)

l3
0.5 0.9 1.5 2.5 3.2

S1
0 0 0 0 0

S2
0 0 0 0 142

S3
207 198 184 169 0

Q = sum( S i )
207 198 184 169 142

TP
778 762 746 722 688

Figure 4.16: Effect of Transportation Cost
800 600 400 200 0 0.5 0.9 1.5 2.5 3.2 L3 (transportation cost) TP S3 S2 Q

It’s also interesting to study the effect of price ( pi ) on sales price and total sales quantity, and we have:
Remark 4.3: If there exists optimal solution S i* for problem (4.9)-(4.10), then Q *
(Q =
*

?S
i

* i

) increases in pi .

Proof: See appendix 15. This result can be intuitively explained as follows: sales locations are only tied together by a concave production function so if something induces more sales at one location then this will drive down unit costs and cause an increase in profitability at other locations.

84

4.3 Analytical Model for Balanced Supply
Competition becomes stiffer and margins get smaller. Globalization, more demanding customers and shortened product lifecycles are challenges in today’s competition. Companies are continuously forced to improve their performance in order to create value-added (VA) to customers from day to day, if they want to remain profitable. As companies continuously seek to provide their products and services to customers faster, cheaper, and better than their competitors, they realized that they cannot do it without considering the “product path”. From supply chain’s perspective, this means that products have to reach the customers from the right supplier. Supply chain diversity ranges from globally dispersed manufacturers, distribution centers and sales subsidiaries with different production cost, capacities, capabilities and lead-times for different products. Therefore identifying the right supplier/path is crucial for effective order promising, capacity utilization and production smoothness. Overcoming the diminishing profit margin and achieving the resource allocation efficiency stimulates us to do path analysis for the global supply chain. Again, when we say product path, we refer to the path from the supplier where the raw materials are purchased to manufacturer where products are produced, through the distribution center (D.C.), and to the sales locations where the products are finally sold to the customers.

85

In previous section we stated that, to counteract price erosion and the resulting reduction in profit margin, manufacturers need to align the production and logistics planning with the end-sales to choose the right path. Especially, in such a risky and uncertainty environment, from the purchasing of raw material to manufacturing, to the D.C, and finally to the customer, companies have to use different suppliers in order to avoid vulnerability in market competition and reduce risk (L.Hunter 2004). This brings us our major concern: the “balanced supply”. To meet customer demand, decision-makers on the supply side include suppliers, manufacturers and distributors. They are concerned not only with profit maximization but also with risk minimization (A. Nagruney 2004). At the same time, developing more responsive strategies requires multiple suppliers to respond to different supply chain drivers. One of the time-based competition strategies proposed by D. Kritchanchaj (1999). focuses on how to improve flexibility and responsiveness of business processes to meet customer requirements. Organizations need to ensure that they continually monitor the changing demands of customers and then attempt to meet their customer's expectations in order to defend their market position against competitors. The capability of responding quickly to customers’ demand is a key business process. But the bottom line is the responsiveness of different suppliers varies. Hence, meeting customer demands with different suppliers instead of adopting only one supplier, thus has become a prerequisite for business survival in the face of market globalization where rapidly changing business environments and seemingly insatiable customer expectations have become the norm. This is especially true for companies that incur the burden of high logistics costs. For example, Amazon.com has at least 2 fulfillment

86

centers to cover any customer zone in order to satisfy its trademark “shipping guarantee policy”. Achieving this feat during Christmas season is paramount for any retailer. All these point us in one direction: reducing the risk through balanced supply.

4.3.1 Model & Formulation

The problem under consideration below is a single time period, single product, multiple-tier supply chain model. In particular, we consider one focal company with multiple suppliers, which provide the material to manufacturing sites / distribution sites. The manufacturing sites are involved in the production of a homogeneous product which is then shipped to distributors, who, in turn, ship the product to end customers. To simplify the deriving process, the supply chain path is illustrated to be from one supplier to one sales location. (See Figure 4.17). The fundamental decision in the model is to decide the optimal supply of products from each supplier to meet the end customer’s demand so that the profit of the company is maximized, while
supply balance is also achieved. This problem is modeled as a constrained non-linear

integer programming problem.
4.3.1.1 Inputs
• • •

Set F of Suppliers, I = n. In general elements of I will be identified by i. Customer demand Q : a given quantity. Unit cost Ci : This represents linear cost such as materials needed

87

Figure 4.17 Supply Balance

Unit cost: Ci

Supplier

Path i

Customer Demand : Q



Cost function f:

This cost function represents nonlinear costs such as

production costs and transportation costs, which include shipping and handling fees. It is concave function to reflect economies of the scale.


Supply balance reward Bi : With more suppliers to commit the customer demand, the company has more flexibility to quickly respond to changes and more alternatives to overcome the uncertainties, and thus reduce the risk. We use Bi to represent such benefits in the model.



Supply Lower bound li : the minimum quantity that each supplier should provide in one time period, so that path shutdown can be avoided and meanwhile, production smoothness, TL transportation and the economies of the scale can be achieved.

4.3.1.2 Decision Variables

xi : the quantity committed by path i. yj = 1 if supplier i is chosen to fulfill demand , and yj = 0 otherwise.

88

In other words: yi = 1 if xi >0,
yi = 0 if xi = 0.
4.3.1.3 Model

Minimize Subject to:

? ( f (x ) + C x ) ? ? B y
i i i i i i

i

(4.15)

?x
i

i

=Q

(4.16)

yi = 1 if xi >0, yi = 0 if xi = 0. xi ? l i if xi >0

(4.17) (4.18)

We can see, (4.17) and (4.18) are equivalent to: if y i = 1, then xi ? li ; if y i = 0, then xi = 0 , which can be expressed as:
xi ? M ? y i xi ? l i ? y i xi ? 0 y i = 0,1

and the model can be simplified as: Min subject to:

? ( f (x ) + C x ) ? ? B y
i i i i i i

i

(4.19)

?x
i

i

=Q

(4.20)

89

xi ? M ? y i xi ? l i ? y i xi ? 0 y i = 0,1

(4.21) (4.22) (4.23) (4.24)

We will always assume that Q > ? l i , all parameters are positive, and f (?) is
i

concave.

4.3.2 Model Analysis

Without loss of generality, we assume: C1 < C 2 < ..... < C n , the following results are based on assumption that: C 2 ? C1 > f ' (0) ? f ' (Q) .
Theorem 4.3: For any optimal solution X*, x1 ? l1 .

Proof: see Appendix 16.

This Theorem tells us 1) that the supplier with the lowest unit cost will definitely be allocated at least the lower-bound quantity 2) it only depends on the unit cost, the nonlinear shipping and handling cost doesn’t play any role here if

C 2 ? C1 > f ' (0) ? f ' (Q) .

From Theorem 4.3, we can further have:
Theorem 4.4: For any optimal solution X*, we have: xi ? li . (Since li is lower

bound, this is equivalent to say: for all i ? 2 : either xi = 0 or xi = li )

90

Proof: see Appendix 17.

These Theorems give us a clear and simple solution. It shows us that the first supplier is used at l1 and for all the other suppliers, the solution space is not infinite but a very limited number of possibilities.

From Theorem 4.3 and Theorem 4.4, we can see that, without any complicated assumptions and constraints involved, this model provides a simple structure relative to demand fulfillment and balanced utilization of the supply chain. The model’s conclusion can also be broadly applied to other similar situations.

4.4 Available Resource Analysis
Identifying profitable paths is also a powerful weapon coordinating an available resource such as component inventory or capacity, with order management (sales subsidiaries) policies because it enables a company to know which path is profitable and therefore can improve the customer service at least cost. The path in the real-life global supply chain is difficult to handle analytically since they involve interactions among several individual factors as well as thousands of products and product locations. To reduce complexities and improve insight, we create an aggregate level model, which means sales subsidiaries, rather than the individual orders, are the unit of attention.

91

In the research below we study the firm’s optimal response to a price discount on a component (raw material) under two different strategies. In the first strategy, the firm does not pass along the discount to its customers (sales subsidiaries); the firm simply coordinates sales among the different products and subsidiaries to minimize the financial impact of non–profitable products. In the second strategy, the firm offers price discounts in different sales subsidiaries to increase the demand for profitable products. We carried out experiments for the two strategies based on a mathematical programming model, built around Toshiba’s global notebook supply chain. Model constraints include, among others, material constraints, bill–of–

materials, capacity and transportation constraints, and a constraint on minimum fill rate (service level constraint). Unlike most models of this type in the literature, which define variables in terms of single arc flows, we employ path variables, which allow for direct direction identification and manipulation of profitable and non–profitable products.
4.4.1 Model background

Facing fierce competition, many companies have to cope with two conflicting goals: one is to maximize total profit; the other one is to maintain a high level of service, usually measured by fill–rate. For strategic reasons, a firm may set a lower bound on its fill rate, be it global (across products and locations), or local for a given location. This constraint usually forces a company to sell non–profitable products. Meantime, due to the different path involved, similar products might be profitable in one path (location) but non-profitable in the other. In this study, we consider n

92

products, which are variations of a generic product, and can be thought of as configurations for a notebook PC. The products are manufactured in different

factories at different costs, and sold in different sales subsidiaries, also called sales locations, at different prices. Temporary price discounts by suppliers are a relatively common phenomenon–– these occur as a result of supply–demand imbalances (e.g., production overruns, poor forecasting, etc.), competition among suppliers, retooling, etc. (Tersine and Barman 1995). In this research we consider the issue of coordinating available resource with sales at different locations in a global, capacitated, sourcing–manufacturing– distribution supply chain under such temporary price discount. This research is based on a two–year project conducted with Toshiba Corporation. All products have a common and critical component, for which the supplier offers a price discount in one time period. Customers are price–sensitive and demand for a product is assumed to decrease linearly with its price. We define a sales subsidiary’s fill rate as total committed quantity divided by its total demand. There are different minimum fill rate

requirements across different sales subsidiaries. This minimum fill rate requirement may induce sales products that are non–profitable in that location. We propose two coordinating mechanisms to improve the resource utilization and reduce the financial impact of non–profitable products. The first mechanism simply buys a large amount of components at the discounted price to use them in non–profitable products for as long as possible––the length of time is related to inventory holding costs, and capacity considerations along the chain (transportation, production, etc.). An optimal

93

procurement quantity trades off gains from the discounted price against additional holding costs, taking into account the various capacity constraints. The second mechanism reduces the quantity of non–profitable products needed to meet the minimum service level constraint by increasing sales of profitable products via a price discount. Note that offering price discounts for profitable products has the potential of increasing profits due to two effects: increased demand, and decreased sales of non–profitable products, however, too large a price discount can ultimately hurt profits by decreasing total revenue. We should point out here, without the product path to identify profitable product, such coordinating mechanism is impossible to be carried out.
4.4.2 Toshiba Supply Chain

In this section we model Toshiba’s notebook (PC) global supply chain (Figure 4.18), which comprises four final assembly and testing (FAT) factories and six sales subsidiaries, including locations in Japan, The Philippines, United States and Germany. Toshiba buys PC components, including CPU, hard disk drive, keyboard, LCD, CD–ROM, and DVD–ROM, directly from suppliers, and transports the motherboard subassembly from subassembly factories; these are shipped to the FAT factories for assembly, and finished PCs are finally shipped to the different sales subsidiaries. The motherboard components are bought from suppliers, and all subassembly factories can produce all types of motherboards. There are over 3,500 different product models offered across the different sales subsidiaries. Although our

94

model is based on Toshiba Corporation, it is clear that the structure of its supply chain is similar to many other firms operating in similar industries.

Suppliers

FAT

Sales

Sales Australia Japan FAT Supplier Singapore Philippine FAT Supplier Japan USA FAT Supplier German Sales Japan Sales USA Subsidiaries In Europe Sales Singapore Sale Canada

Customer Customer

Customer

Customer Customer Customer

Germany FAT

Figure 4.18 Toshiba Global Supply Chain

Our sourcing–production–distribution planning model spans a time horizon of 13 weeks. Unmet demand for a sales subsidiary at any period cannot be backordered; consequently committed quantity is equal or lower than demand. The supplier

network is not considered in this model and thus we ignore holding costs at the suppliers. The integrated supply chain model is a multi–period, multi–echelon and multi–product MIP model (integer variables are necessary because of minimum lot size requirements), which is presented in §4.4.3 From a modeling perspective, our global supply chain model is rather straightforward, except for a choice in the decision variables, which allows us to answer the research questions posed in §4.4.1. As explained in the Chapter 1, unlike

95

most models in the global supply chain literature, which define variables in terms of flows along a single arc in the network and use flow balancing constraints at nodes, we employ path variables, which provide location–specific cost information and directly identify profitable and non–profitable products In our two–year research project with Toshiba, we found this to be very important managerial information, as it serves as an aid for other strategic decisions such as product line offering at different locations. Although this modeling approach certainly increases the number of

decision variables (for example, a typical model we study here has over a million variables and a similar number of constraints), we find that the resulting models, even with real–world data, are manageable with solution times being around 5 minutes using typical computing environments.

4.4.3 Model Formulation Notation

Index Use (Indices and Index Sets)
s ? S Sales subsidiaries

f ? F Final assembly and testing (FAT) factories
l ? L Subassembly factories. i?I

Product (notebook PC models)

j ? J Components sourced directly from suppliers j ? J Motherboard (sourced from subassembly factories)
k ? K Motherboard components.

96

Data (Lower Case Letters)

?

Tangible profit weight (objective function tradeoff parameter)

d is (t ) Demand for product i ? I at sales subsidiary s ? S at time t = 1,2,L , T pis (t ) Sales price for product i ? I at sales subsidiary s ? S at time t = 1,2,L , T
a jf (t ) Quantity of component j ? J available at factory f ? F at time t = 1,2,L, T

bij

Quantity of component j ? J needed to produce a unit of product i ? I ; similarly for bi j and b jk

?s

Minimum fill rate for sales subsidiary s ? S (total quantity committed divided by total demand)

q fs (t ) Transportation capacity between factory f ? F and sales subsidiary s ? S at

time t = 1,2,L , T ; similarly for q lf (t )
c f (t ) Maximum production rate at factory f ? F , notebook PC units, at time t = 1,2,L , T ; c l (t ) is similarly defined for motherboard units zif

Minimum lot–size for product i ? I at factory f ? F ; z lj is similarly defined Transportation lead–time between factory f ? F and sales subsidiary s ? S ;

? fs

? lf is similarly defined ?if
Production lead–time for product i ? I at factory f ? F ; ? lj is similarly defined

97

rj f

Unit purchasing cost of component
t = 1,2,L , T ; rkl is similarly defined

j?J

at factory

f ?F

at time

his

Unit holding cost for product i ? I at sales subsidiary s ? S ; similarly defined

hkl , h jf are

vif

Production cost per unit for product i ? I at factory f ? F ; v lj is similarly defined

w fs

Transportation cost per unit between FAT factory f ? F and sales subsidiary
s ? S , including taxes or duties; wlf is similarly defined.

Path Decision Variables (capital letters) Dilfs (t ) Quantity of product i ? I produced at factory f ? F by using motherboard

from sub–assembly factory l ? L and shipped to sales subsidiary s ? S at time t = 1,2,L , T

Auxiliary Variables (Capital Letters) M is (t ) Demand commitment for product i ? I at sales subsidiary s ? S at time

t = 1,2,L, T
Qi f (t ) Quantity of product i ? I produced at factory f ? F at time t = 1,2,L, T ;
Q lj (t ) is similarly defined

98

Z i f (t ) Binary variable (= 1, if product i ? I is produced at factory f ? F at time t = 1,2,L , T ; 0 otherwise); Z lj (t ) is similarly defined H is (t ) Inventory on hand for product i ? I at sales subsidiary s ? S at time t = 1,2,L , T ( H is (0) is known and given);
H jf (t ) and H kl (t ) are similarly

defined
C is (t ) Total cost for product i ? I sold at sales subsidiary s ? S at time t = 1,2,L, T Pi s (t ) Total profit for product i ? I sold at sales subsidiary s ? S at time

t = 1,2,L, T
% (t ) C
Total PC component and motherboard component holding cost at time
t = 1,2,L , T

The objective function maximizes profit––revenue from promised orders minus transportation, production, duty, component, and inventory costs. Constraints

include: demand commitment and fill rate constraints, inventory balance constraints at factories and sales subsidiaries, minimum lot size, and production and transportation capacity constraints. Objective function: maximize profit max Subject to: Demand commitment and minimum fill rate constraints
M is (t ) ? dis (t ) i ? I , 1 ? t ? T ,
s ? ? Pi (t ) T s

t =1 i

(4.25)

(4.26)

99

i?I

s s s ? M i (t ) ? ? ? di (t ) s ? S , 1 ? t ? T , i?I

(4.27)

Inventory balance constraints at sales subsidiaries
H is (t ) = H is (t ? 1) + ? Dilfs (t ? ? fs ) ? M is (t )
l?L , f ?F

i ? I , s ? S ,1 ? t ? T

(4.28)

Flow conservation constraints at factories
Qi f (t ? ? if ) = ? Dilfs (t ) i ? I , f ? F ,1 ? t ? T
l?L , s?S

(4.29)

Inventory balance constraints for components at factories
H jf (t ) = H jf (t ? 1) + a jf (t ) ? ? Qi f (t ) ? bij j ? J , f ? F ,1 ? t ? T
i?I f ? Qi (t ) ? bi j =

(4.30) (4.31)

i?I

i?I , s?S ,l?L

?

Dilfs (t ? ? lf ) ? bi j

j ? J , f ? F ,1 ? t ? T

Flow conservation constraints at subassembly factories
Q lj (t ? ? lj ) =
i?I , s?S , f ?F

?

Dilfs (t ) ? bi j

j ? J , l ? L,1 ? t ? T

(4.32)

Inventory balance constraints for raw material at sub–assembly factories
l H kl (t ) = H kl (t ? 1) + ak (t ) ? ? Q lj (t ) ? b jk k ? K , l ? L,1 ? t ? T j ?J

(4.33)

Minimum lot size constraints at factories
Qi f (t ) ? Z i f (t ) ? zif i ? I, f ? F,

(4.34) (4.35)

Qi f (t ) ? Z i f (t ) ? N i ? I , f ? F ,

Minimum lot size constraints at subassembly factories
Q lj (t ) ? Z lj (t ) ? z lj j ? J , l ? L, Q lj (t ) ? Z lj (t ) ? N j ? J , l ? L,

(4.36) (4.37)

Production and transportation capacity constraints

100

i?I

f f ? Qi (t ) ? c (t ) f ? F ,1 ? t ? T

(4.38) (4.39) (4.40) (4.41)

j ?J

l l ? Q j (t ) ? c (t ) l ? L,1 ? t ? T

i?I ,l?L

lfs fs ? Di (t ) ? q (t )

f ? F , s ? S ,1 ? t ? T l ? L, f ? F ,1 ? t ? T

i?I , s?L , j ?J

?

Dilfs (t ) ? bij ? q lf (t )

Cost and profit calculation constraints
?vif + w fs + ? bij ? rjf + ? j?J ? ? s s s lfs Ci (t ) = hi ? H i (t ) + ? Di (t ) l l? ? f ?F ,l?L ? bi j ? v j + ? bi j ? b jk ? rk ? ? j ?J , k?K ? j ?J ? i ? I , s ? S ,1 ? t ? T ,

(4.42)

% (t ) = C

j?J , f ?F

f f l l ? h j ? H j (t ) + ? hk (t ) ? H k (t ) k?K ,l?L

1? t ? T,

(4.43) (4.44)

% (t ) / d s (t ) i ? I , s ? S ,1 ? t ? T , Pi s (t ) = pis (t ) ? M is (t ) ? Cis (t ) ? C i

Integrality and Non–negativity
Z i f (t ) ? {0,1}, Z lj (t ) ? {0,1} ,

(4.45)

Cis (t ) ? 0, Qi f (t ) ? 0, Q lj (t ) ? 0, Dilfs (t ) ? 0, H i f (t ) ? 0, H jf (t ), H kl (t ) ? 0, Pi s (t ) ? 0

(4.46)

We now elaborate on how (purchasing) component costs are handled in the model. Materials availability, given by a jf (t ) ’s, are parameters (inputs) to the model. Consequently, one would be inclined to treat purchasing costs as “sunk” costs; in this manner they would not influence the optimal solution––given by the path decision variables Dilfs (t ) . In the computation of profits (4.42)– (4.44), however, we explicitly incorporate these purchasing costs as part of the unit cost for the path decision
101

variables, and thus they affect the firm’s optimal solution. We justify this approach by noting that the problem has a long (13–period) planning horizon, and therefore the firm may hold components in inventory for future use in any given period.
4.4.4 Result Analysis

In the fist coordinating mechanism, we don’t differentiate the profitable or nonprofitable product and we shall see how the available resource would affect the profit and interact each other under such scenario. Then in the second strategy, we will treat the profitable product differently to get optimal response. The model is solved with the MPL 4.11, with a Cplex 7.0 solver (Cplex 1998) operating on an IBM server with a Pentium III processor, under Windows NT (4.0). The computer RAM is 1,047,960. A typical model has 1.16 million variables and 810,000 constraints; a typical running time is 7 – 10 minutes.
4.4.4.1

To address the first research question, we design a numerical study based on data collected at Toshiba (Figure 4.18), however, we disguise the actual magnitude of profits and demands for obvious reasons. Based on the literature review in Chapter 2, the most important factors in our model that influence the firm’s optimal response to a one–time discount by a critical supplier would be: the magnitude of the discount rate, the minimum fill rate (service level), production and transportation capacity, and the holding cost. Accordingly, we design a full–factorial experimental design with three levels for each of these factors, where one of the levels represent the current value used at Toshiba; the other two levels represent alternative scenarios that are

102

considered plausible by the company. We thus have a 34 factorial design, which is shown in Table 4.13 and explained below.
Table 4.13: Experimental Design for First Research Question
Parameter Discount Fill Rate Capacity Holding Cost Levels Low (LO) 0% 0.65 0.65 10% Medium (MID) 30% 0.75 0.85 30% High (HI) 50% 0.85 1.00 60%

Discount is defined as a % reduction in the nominal price of a critical component from a major supplier (in this case, CPU), from no reduction to a 50% reduction. The firm has a contractual obligation to purchase “normal” amounts of components from its suppliers throughout the planning horizon; these “normal” amounts are the parameters a jf (t ) . The discount mentioned above is given only to the marginal amount above the “normal” amount. We set three values for a common minimum fill rate ? s across all sales subsidiaries s: 0.65, 0.75 and 0.85. We consider three

scenarios for capacity as follows: we take existing values for the capacity parameters
c f (t ) and q fs (t ) at Toshiba, and consider three multipliers; the relative proportions

among these three multipliers are the capacity ratios shown in Table 4.13. (we disguise which one is the current capacity at Toshiba). Finally, we consider three values for annual holding cost, expressed as a percentage of the product’s price at the sales subsidiary ( his ) or as a percentage of the component’s price from the supplier ( h jf and hkl ): 10%, 30% and 60%; these simulate various levels of obsolescence in this high–velocity environment.

103

We report key summary measures of the model’s optimal solution for this research question: total profit, and total quantity of components purchased (the component for which the supplier offers a price discount). Given the full factorial study design, an effective way to measure the sensitivity of profit and purchased quantity (dependent variables) to the factors (Table 4.13, independent variables) is to compute t–statistics for the respective coefficients in multiple regressions for each dependent variable as a function of all independent variables (Wagner 1995). These statistics are reported in Table 4.14.
Table 4.14: T–Statistic for Multiple Regressions Where Dependent Variables are Profit and Purchased Quantity and Independent Variables are Factors (n = 81)
Factor Discount Fill Rate Capacity Holding Cost Profit (R2 = 0.78) 0.4 –11.4 12.3 1.1 Purchased Quantity (R2 = 0.95) 2.6 39.1 10.8 –4.2

Not surprisingly, both fill rate and capacity have a strong effect on both profit and purchased quantity, as evidenced by the high values of the t–statistic. The effect of the discount rate on profit is non–significant (t–value = 0.4) since purchasing cost for the component represents only a small fraction of total profit; the same is true for holding cost. Note that the discount is offered only in one period, and total profit is defined as the sum over the entire planning horizon of 13 periods. The effect of the discount rate on purchased quantity is relatively mild compared to the effects of the minimum fill rate and capacity. There are several reasons for this surprising mild effect. There are production and transportation capacities in the supply chain, as well as constraints for other components. In addition, the discount is offered only for units

104

above the minimum purchased quantity specified in the contract with the supplier. Thus, the discounted components, without a change in the final product’s price, does not result in a demand increase for the final product, and is therefore mostly used to increase the fraction of demand that is met with profitable products. The effect of unit holding cost on profit is also non–significant (Table 4.14); this occurs because holding costs comprise only a fraction of total profit over the 13– period planning horizon. The effect of holding cost on the purchased quantity is mild compared to the effect of the minimum fill rate and the capacity, which is a surprising result, considering the lot–sizing literature. For example, if the unit holding cost increases six–fold, the EOQ formula posits that the optimal lot size would decrease by 60% (= 1 – 1 6 ); our model’s optimal solution recommends instead that the total quantity of purchased components decreases by an average of only 7.3% (Figure 4.19). Overall, these are interesting results because they point out weaknesses of simple models in the literature––where assumptions of a single product, uncapacitated supply chains, and no service level constraints are frequently used––in dealing with the complexities of global supply chains, where these assumptions do not hold, and which significantly influence a firm’s optimal sourcing–production– distribution decision. Next, we look at the impact of each factor on profit and purchased quantity, by taking an average across experiments, for each factor level. These are shown in Figure 4.19 and Figure 4.20 for profit and purchased quantity, respectively. The relationship between (minimum) fill rate and profit (Figure 4.19) displays the usual

105

shape of profit–service trade–off curves—that is, increasing minimum fill rates decreases profits at an increasing rate. Also, increasing capacity beyond a certain level (here, mid level) does not result in significant profit increases (Figure 4.19) or in additional purchased components (Figure 4.20) and this can be explained as follows. Additional capacity beyond the low level is used to produce profitable products, however, additional capacity beyond the middle level would be used to produce non– profitable products beyond the required service level, which is not economically attractive and thus does not happen (Figure 4.21); the purchased quantity for non– profitable products remains unchanged with the capacity level at 267). As discussed before, holding cost and discount do not significantly affect profits, as evidenced by the respective “flat” curves.
Figure 4.19: Effect of Holding Cost, Capacity, Fill Rate and Discount on Profit
Profit (averaged across experiments)
2.40E+08 2.35E+08 2.30E+08 2.25E+08 2.20E+08 2.15E+08 2.10E+08 2.05E+08 2.00E+08 LO MID HI

Level of Input

Holding Cost

Capacity

Fill Rate

Discount

Figure 4.20: Effect of Holding Cost, Capacity, Fill Rate, and Discount on the Purchased Quantity

106

Purchased quantity (average across experiments)
90,000 80,000 70,000 60,000 50,000 40,000 30,000 LO MID HI

Level of Input

Holding Cost

Capacity

Fill Rate

Discount

Figure 4.21: Effect of Capacity Increase in the Production of Profitable Products (Average Across Experiments Where Minimum Fill Rate = 0.85)
Purchased Component Quantity for Profitable Products
1700 1650 1600 1550 1500 1450 LO MID HI

Capacity

4.4.4.2

We now investigate the possibility of passing on part of the component savings to the customers by offering price discounts for profitable products at sales subsidiaries, which increases their demand. As demand for profitable products increases, the

107

amount of non–profitable products needed to meet the minimum fill rate requirement decreases (we note that the total demand at a sales subsidiary s is composed of several individual demand quantities, d is (t ) .). For ease of presentation, we assume that all
profitable products at all locations are offered the same price discount ?, and for

simplicity we assume the same (linear) demand curve across products and locations:

% s (t ) = d s (t ) ? (1 + ? ? ? ) , where ? is a price elasticity parameter, and we use the “~” to d i i
differentiate between the “discounted” demand and the base demand data dis (t ) . Offering discounts for profitable products potentially increases profit due to two effects: increased overall demand, and decreased sales of non–profitable products, however, too large a price discount may hurt profits by decreasing total revenue. We use the same base data from §4.4.4.1; in addition, we fix the capacity ratio at its high level 1.0, the holding cost at its low level 10%, and the minimum fill rate at a high level of 0.9; we vary the price discount ? from 5% to 20%. The price discount is initially offered in period 2 (of our 13–period model), and we study four scenarios, where we discount the product for 1, 2, 3 and 4 consecutive periods. We consider three levels for the demand elasticity parameter ?: 0.5, 1 and 10 (these values indicate that a 10% price discount results in a demand increase of 5%, 10% and 100%, respectively.) Across all numerical examples, the CPU supplier offers a discount of 25% on the price of the component, where, as described in §4.4.4.1, the discount is offered only to the marginal amount above the “normal” amount, or contractual obligation, a jf (t ) . In contrast, the price discount ? offered to the final product

108

customers is applied to all products sold in that period and location (that is, there is no price discrimination across customers in a location.) The curves showing profit as a function of price discount can be seen in Figure 4.22, for various levels of demand elasticity ?, when the price discount is offered 4 periods. Note that for ? = 0.5––an inelastic demand curve––profit decreases

monotonically with the discount rate, and thus the firm should not offer a price discount. Note that profits decrease steeply for discount rates above 10% when ? = 0.5; this is because the additional demand generated is not large enough to overcome the loss in revenues, and several products may become non–profitable. The optimal discount rate is around 10% for ? = 1 and 10. To meet the additional demand generated by the price discount, the firm increases the quantity of components purchased at the discounted rate, as shown in Figure 4.23 (note that Figure 4.23 shows the total quantity of the component purchased at the discounted price not the total quantity of the component used, which is higher).

109

Figure 4.22: Profit vs. Price Discount
Profit (Price Discount Offered in 4 Periods)
2.80E+08 2.60E+08 2.40E+08 2.20E+08 2.00E+08 1.80E+08 1.60E+08 1.40E+08 1.20E+08 1.00E+08 0% 5% 10% 15% 20%

Discount Rate ? ? = 10 ?=1 ? = 0.5

Figure 4.23: Purchased Component Quantity vs. Price Discount
Purchased Quantity (Price Discount Offered in 4 Periods)
800,000 700,000 600,000 500,000 400,000 300,000 200,000 100,000 0 0% 5% 10% 15% 20%

Discount Rate ? ? = 10 ?=1 ? = 0.5

4.4.5 Remarks

We have studied the problem of responding to a supplier’s one–time price discount for a critical component that is common to several products on a global

110

supply chain, using real data from Toshiba Corporation.

We use an MIP–type

formulation for a 13–week planning horizon, but for implementation purposes we simplify it to a linear program. Our model is a tactical planning model; we neither consider strategic issues such as facility location nor tactical issues such as production scheduling. Our model’s primary decision variables are flows of materials in the supply chain along paths (as opposed to flows along single arcs), which isolates critical information regarding profitable and non–profitable products across locations—products that are profitable in one location may not be profitable in another location. This information serves as an input to the more strategic decision of product line offering. Our model maximizes profit subject to a critical service level constraint—a minimum fill rate, which may be different across locations. Because of the minimum fill rate, the firm may need to sell non–profitable products. Thus, our research studies one of many weapons a firm has to mitigate the effect of non– profitable products on the firm’s profitability. Temporary price discounts by suppliers are a relatively common phenomenon. We consider the issue of coordinating a temporary price discount for a critical component with pricing and sales at different locations in a global, capacitated, sourcing–manufacturing–distribution supply chain. Based on data from Toshiba, we perform a numerical study where we vary critical parameters on the supply chain— capacity, discount rate, holding cost, number of discounted periods, and minimum fill rate—and analyze their impact on the firm’s profitability. We find that capacity critically impact profits, however, adding capacity beyond a certain level, without changes in the firm’s pricing strategy, only impacts the production of non–profitable

111

products, which does not result in profit improvement if the firm is above the minimum fill rate threshold. The relationship between minimum fill rate and profit has the usual shape of a profit–service trade–off curve. We also point out the

weaknesses of analytical models in dealing with the complexities of decision–making in global supply chains. For example, we show how a six–fold increase in holding cost for a component only results in a 7% increase in the optimal lot size, whereas the EOQ model predicts a 60% increase. This occurs because the assumptions of

uncapacitated supply chains—including production and transportation—as well as a single product type, typical of analytical models, do not hold in complex supply chains. We find that passing on temporary component price discounts to final product customers is attractive if the demand curve is “elastic” enough such that a price discount significantly stimulates demand. Also, a price discount needs to be offered for several periods (weeks in this case) for it to be effective and profitable. For example, it is not optimal to offer price discount on the final product during only one period under any scenario analyzed in this paper. We find that the relationship between price discount and the number of periods when the discount is offered is non monotonic, requiring careful optimization. Although our study was conducted with real data from a global company, we caution on some of its limitations. The primary limitation regards modeling the demand curve. Essentially, we have assumed demands to be independent across products—discounts offered in one product would not impact demand for other products. Also, our analysis considers competitive issues only on a limited basis,

112

through the estimated demand curve––for example, we assume that if the firm offers no discounts on its final product, even after receiving a discount on a critical component by a key supplier, then its demand does not change. In reality, it is likely that the supplier offers discounts to several manufacturers; thus the demand curve should include price dynamics by all other competitive products in the marketplace; this, however, significantly complicates modeling and parameter estimation for the demand curves.

113

Chapter 5 Conclusions
Global competition and e-commerce have imposed tremendous pressure on product and service providers to get closer to the customers. At the same time end consumers are becoming increasingly knowledgeable and demanding. Consequently, in the current customer-centric business environment, supply chains must be designed to accommodate for real-time responsiveness, uncertain customer orders, globally dispersed locations and diminishing profit margin. Therefore, managing supply chains in today's distributed manufacturing environment has become more complex. To remain competitive in today's global marketplace, organizations must streamline their supply chains. The practice of coordinating design, procurement, flow of goods, services, information and finances, beginning from raw material flows, parts supplier, manufacturer, distributor, retailer, and finally to consumer requires synchronized planning and execution. It is of critical importance to understand how an efficient and effective supply chain is affected by order promising and order fulfillment strategies. In this dissertation, we study two issues related to moving a company from an Available to Promise (ATP) philosophy to a Profitable to Promise (PTP) philosophy: pseudo order promising and coordinating demand fulfillment with supply. To address the first issue, a single time period analytical ATP model for n confirmed customer orders and m pseudo orders is presented by considering both

114

material constraints and production capacity constraints. At the outset, some analytical properties of the optimal policies are derived and then a particular customer order promising scheme is derived based on the ratio between customer service level and profit. Algorithms presented to solve this problem provide decision makers the effective rule-based mechanism to implement and deploy critical decisions in realtime ATP systems. To tackle the second issue, we explore two cases: a deterministic demand curve or stochastic demand. In the first case, a constrained non-linear programming model is developed. For the second case a stochastic constrained model is formulated to determine both the quantity of raw materials to purchase from suppliers and demand fulfillment levels for each end product in every location. A simple, yet generic optimal solution structure is derived and a series of numerical studies and sensitivity analysis are carried out to investigate the impact of different factors on profit and fulfilled demand quantity. The objective of this analysis is to understand the implication of factors like price elasticity, path cost, shortage cost and holding cost, along with their interaction effects within the firm's supply chain. Further, we present an analytical model to explore balanced supply. Implementation of the resulting generic rules reduces response time and provides managers with insight into the optimal policies that improve demand fulfillment and simultaneously reduce risk through balanced supply. Further, the firm’s optimal response to a one-time-period discount offered by the supplier of a key component is studied under two different strategies: a) Not passing along the discount to its customers b) Offering price discounts to increase the demand
115

for profitable products. Unlike most models of this type in the literature, which define variables in terms of single arc flow, this model employs path variables to directly identify and manipulate profitable and non-profitable products. Numerical experiments based on Toshiba’s global notebook supply chain are conducted. Based on the results, an interesting relationship among capacity, fill rate and profit is observed. Furthermore, it is found that passing on temporary component price discounts to end product customers is attractive if the demand curve is “elastic” enough so that a price discount significantly stimulates demand.

116

Appendices
Appendix 1:

Proof of Theorem 3.1: It suffices to show that the Hessian matrix of the objective function, Hn+1×n+1, is negative definite. Let En+1×n+1 be the matrix with each entry of 1. Then we have H=diag(
?1 ? X ? X X s' ' ( 21 ) , 2 s' ' ( 22 ) ,…, n s' ' ( 2n ) , -(p+h)f(Y))-( g' ' (z) + m' ' (z) )E 2 2 2 q1 q1 q 2 q2 qn qn
n

where z= ? X i + Y .
i =1

Let ?=(?1,…, ?n+1)T,
? T H? =?
i =1 n n +1 X ?i 2 2 s' ' ( i )? i ? (p + h )f (Y)? n +1 ? [g' ' (z) + m' ' (z)](? ? i ) 2 qi qi i =1 n +1 i =1 n 2 2

? ?M (? ? i ) 2 + LU? ? i ? (p + h )K? n +1
i =1

If M > 0, RHS ? 0 and RHS = 0 only if ? =0. Hence H will be negative definite. If M ? 0, we have,
? T H? ? ?M (n + 1)? ? i + LU? ? i ? (p + h )K? n +1
2 2 i =1 i =1 n 2 n +1 n 2

= [LU ? M (n + 1)]? ? i ? [(p + h )K + M (n + 1)]? n +1
i =1

2

then, by the assumption of the theorem, H is still negative definitive. Q.E.D.

117

Appendix 2:

Proof of Remark 3.1: The second derivative of TB at X i and Y :
n n ? X ? 2 TB = ? ' ' ( + ) ? ' ' ( g X Y m X i + Y ) + i2 s ' ' ( i ) ? ? i 2 qi ? Xi qi i =1 i =1

n n ? 2TB = ? pf ( Y ) ? g ' ' ( X + Y ) ? m ' ' ( X i + Y ) ? hf (Y ) ? ? i ? 2Y i =1 i =1

n n ? 2TB = ? g ' ' ( ? X i + Y ) ? m' ' ( ? X i + Y ) ?X i ?X j i =1 i =1

n n ? 2TB = ? g ' ' ( ? X i + Y ) ? m' ' ( ? X i + Y ) ?X i ?Y i =1 i =1

n n X ? ? 2TB = ? ' ' ( + ) ? ' ' ( g X Y m X i + Y ) + i2 s ' ' ( i ) ? 0 ? ? i 2 qi ? Xi qi i =1 i =1 n n ? 2TB = ? pf ( Y ) ? g ' ' ( X + Y ) ? m ' ' ( X i + Y ) ? hf (Y ) ? 0 ? ? i ? 2Y i =1 i =1

Now we need to prove the Hessian matrix to be negative definitive. Hessian Matrix:

118

? ? 2TB ? 2 ? ? X1 ? ? 2TB ? ? ?X 2 ?X 1 . ? ? ? 2TB ? ? ?X 2n ?X 1 ? ? TB ? ?Y?X 1 ?

? 2TB ?X 1?X 2 ? 2TB ?2 X 2 . ? 2TB ?X n ?X 2 ? 2TB ?Y?X 2

.. .. . . .
n

? 2TB ?X 1?X n ? 2TB ?X 2 ?X n . ? 2TB ?X n ?X n ? 2TB ?Y?X n

? 2TB ?X 1?Y ? 2TB ?X 2 ?Y . ? 2TB ?X n ?Y ? 2TB ? 2Y

? ? ? ? ? ? ( n +1)×( n +1) ? = HM ? ? ? ? ? ?
n

For simplicity, let g (?) = g (? X i + Y ) , and m(?) = m(? X i + Y )
i =1 i =1

We can decompose the HM ( n +1)×( n +1) to be: HM ( n +1)×( n +1) = H 1 + H 2 + H 3 ? ?1 ? 2 ? q1 ? ? where : H 1 = ? ? ? ? ? ? s' ' ( 0 . 0 0 X1 ) q12 0 0 ? ? ? ? . 0 ? ? ,which is a . . ? Xn ?n ? . s' ' ( 2 ) 2 ? qn qn . ? ( p + h) f (Y ) ? ? . 0

?2
q
2 2

s' ' ( 0 0

X2 ) 0 2 q2 . . . 0

semi-negative definitive.

Q g ' ' () + m' ' () > 0
H 2 = ? g ' ' (?)(I )

H 3 = ? m' ' (?)(I ),

? H 2 + H 3 = ?[ g ' ' (?) + m' ' (?)]( I ) , which is a negative definitive.
HM ( n +1)×( n +1) = H 1 + H 2 + H 3 ,

? HM ( n +1)×( n +1) is negative definitive. Q.E.D.

119

Appendix 3:

Proof of Theorem 3.2: i) Let the sequences { ri +

?i
qi

s ' (1) } and { ri +

?i
qi

s ' (0) } be non-decreasingly ordered.

By [ Z ij ? s ' (0)][ Z ij ? s ' (1)] ? 0 , we can prove for any i, j ri +

?i
qi

s ' (1) < r j +

?j
qj

s ' (1)

ri +

?j ?i s' (0) ? r j + s' (0) and qi qj

ri +

?j ?i s' (1) > r j + s' (1) qi qj

ri +

?j ?i s' (0) ? r j + s' (0) qi qj

Then, if necessary, we can adjust the order of elements in the above two sequences to make sure we have simultaneously for any i,
ri +
?i ? ? ? s' (1) ? ri +1 + i +1 s' (1) and ri + i s' (0) ? ri +1 + i +1 s' (0) qi q i +1 qi q i +1
n * n *

ii) For simplicity, let f * = g' (? X i + Y * ) + m' (? X i + Y * ) .
i =1 i =1

Then there exist unique k and k’ such that
rk +
?k ? ? ? s' (1) ? f * < rk +1 + k +1 s' (1) and rk ' + k ' s' (0) < f * ? rk '+1 + k '+1 s' (0) qk q k +1 q k' q k '+1

By s’(0)>s’(1)>0, it is trivial to prove k ? k’ iii) For all i ? k + 1 : ri ? f (?) +

?
qi

s' (

Xi ) = ?i ? ? i > 0 , qi

i.e.: ?i > 0 , from ?i (qi ? X i ) = 0 , we have: X i = qi , ( i ? k + 1 )

120

For all i ? k ' : ri ? f (?) +

?i
qi

s' (

Xi ? ) < ri ? f (?) + i s ' (0) < 0 qi qi

?

?i ? ? i = ri ? f (?) +

?
qi

s' (

Xi ) < 0, qi

i.e. ? i > 0 , from: ? × X i = 0 , we have: X i = 0, ( i ? k ' ). Q.E.D.

121

Appendix 4:

Proof of Remark 3.2. Let {Xi*} and Y* be an arbitrary solution to (3.24)-(3.25). From
?i ? s' (0) + ri < i +1 s' (1) + ri +1 , we can have: qi q i +1
s' (1) < r1 +

r1 +

?1
q1

?1
q1

s' (0) < r2 +
n

?2
q2
*

s' (1) < r2 +
n

?2
q2

s ' (0) < .... < rn +
*

?n
qn

s' (1) < rn +

?n
qn

s ' ( 0)

For simplicity, let f (?) = g ' (? X i + Y ) + m' (? X i + Y )
i =1 i =1

Without the loss of the generality, there are two cases: (A): if rk +

?k
qk

s ' (0) ? f (?) ? rk +1 +

? k +1
q k +1

s ' (1) Xi ) = ?i ? ? i > 0 , qi

then: 1)for i ? k + 1 : ri ? f (?) +

?
qi

s' (

i.e.: ?i > 0 , from ?i (qi ? X i ) = 0 , we have: X i = qi , ( i ? k + 1 ) 2) For all i ? k :
ri ? f (?) +

?i
qi

s' (

Xi ? ) < ri ? f (?) + i s ' (0) < 0 qi qi

?

?i ? ? i = ri ? f (?) +

?
qi

s' (

Xi ) < 0, qi

i.e. ? i > 0 , from: ? × X i = 0 , we have: X i = 0, ( i < k ) (B) if: rk +
?k ? s' (1) < g ' (?) + m' (?) < rk + k s' (0) qk qk

then: 1) for i ? k + 1 : ri ? f (?) +

?
qi

s' (

Xi ) = ?i ? ? i > 0 , qi

122

i.e.: ?i > 0 , from ?i (qi ? X i ) = 0 , we have: X i = qi , ( i ? k + 1 ) 2) For all i ? k ? 1 :
ri ? f (?) +
n n Xi ? ) < ri ? g ' (? X i + Y ) ? m' (? X i + Y ) + k s ' (1) < 0 ? qi qi qk i =1 i =1 X ? ?i ? ? i = ri ? f (?) + i s ' ( i ) < 0 , qi qi

?i

s' (

i.e. ? i > 0 , from: ? × X i = 0 , we have: X i = 0, ( i < k ) for i =k, need to be determined by (3.24) and (3.25). Combining both (A) and (B), the conclusion is proven. Q.E.D.

123

Appendix 5:

Proof of Corollary 3.5.1. After separating the confirmed orders and pseudo orders, the profit function of pseudo orders is: PP = ? p j E min(Y j , u ) ? ? h j E max(Y j ? u,0)
j?J j?J

Its first order derivative would be: PP’ = ? ( p j (1 ? F (Y j ) ? h j F (Y j ))
j?J

PP’ = ? ( p j ? ( p j + h j ) F (Y j )) ? 0
j?J

And its second order derivative is: PP’’ = ? (? p j f (Y j ) ? h j f (Y j )) , which is:
j?J

PP’’ = ? ? ( p j + h j ) f (Y j )) ? 0
j?J

So we can see, the first order of that function is greater than or equal to 0; and its second order is smaller than or equal to 0; thus the Profit function of Pseudo orders is concave in terms of order quantity. Q. E. D.

124

Appendix 6:

Proof of Theorem 3.3: It suffices to show that the Hessian matrix of the objective function, Hn+m×n+m, is negative definite. Let En+m×n+m be the matrix with each entry of 1. Then we have H=diag(
?1 ? X ? X X s' ' ( 21 ) , 2 s' ' ( 22 ) ,…, n s' ' ( 2n ) , -(pj+hj)f(Yj))-( g' ' (z) + m' ' (z) )E 2 2 2 q1 q1 q 2 q2 qn qn
n

where z= ? X i + ? Y j .
i =1 j?J

Let ?=(?1,…, ?n+m)T,

? T H? n ?i
=?
i =1

qi

s' ' (

n+m Xi 2 2 )? i ? ? ( p j + h j ) f (Y j )? n + j ? [ g ' ' ( z ) + m' ' ( z )]( ?? i ) 2 qi j?J i =1 n 2 2

? ? M ( ?? i ) 2 + LU ?? i ? ? ( p j + h j )K j? n + j
i =1 i =1 j?J

n+ m

If M > 0, RHS ? 0 and RHS = 0 only if ? =0. Hence H will be negative definite. If M ? 0, we have,

? T H?
? ? M (n + m) ?? i + LU ?? i ? ? ( p j + h j )K j? n + j
2 2 i =1 i =1 j?J n+m n 2

? [ LU ? M (n + m)]?? i ? ? [K + M (n + m)]? n + j
2 i =1 j?J

n

2

then, by the assumption of the Theorem, H is still negative definitive. Q.E.D.

125

Appendix 7:

Proof of Theorem 3.4: We first prove part i) Let the sequences { ri +

?i
qi

s ' (1) } and { ri +

?i
qi

s ' (0) } be non-decreasingly ordered.

By [ Z ij + ? ij s ' (0)][Z ij + ? ij s ' (1)] ? 0 , we can prove for any i, j
ri +

?i
qi

s ' (1) < r j +

?j
qj

s ' (1)

ri +

?j ?i s' (0) ? r j + s' (0) and qi qj

ri +

?j ?i s' (1) > r j + s' (1) qi qj

ri +

?j ?i s' (0) ? r j + s' (0) qi qj

Then, if necessary, we can adjust the order of elements in the above two sequences to make sure we have simultaneously for any i,
ri + ?i ? ? ? s' (1) ? ri +1 + i +1 s' (1) and ri + i s' (0) ? ri +1 + i +1 s' (0) qi q i +1 qi q i +1
n * m * n * m *

For simplicity, let f * = g ' (? X i + ? Y j ) + m' (? X i + ? Y j ) .
i =1 j =1 i =1 j =1

Then there exist unique k and k’ such that
rk +

?k
qk

s ' (1) ? f * < rk +1 +

? k +1
q k +1

s ' (1) and rk ' +

?k' ? s' (0) < f * ? rk ' +1 + k '+1 s' (0) q k' q k ' +1

By s’(0)>s’(1)>0, it is trivial to prove k ? k’ For all i ? k + 1 : ri ? f * +

?
qi

s' (

Xi ) = ?i ? ? i > 0 , qi

i.e.: ?i > 0 , from ?i (qi ? X i ) = 0 , we have: X i = qi , ( i ? k + 1 ) For all i ? k ' : ri ? f * +

?i
qi

s' (

? Xi ) < ri ? f * + i s ' (0) < 0 qi qi

126

?

?i ? ? i = ri ? f * +

?
qi

s' (

Xi ) < 0, qi

i.e. ? i > 0 , from: ? × X i = 0 , we have: X i = 0, ( i ? k ' ).

Now we will show ii) From (3-47), we have F (Y j ) =
p j + wj ? f * p j + hj

, where wj>=0. Obviously when

p j < f * , wj has to be positive which then implies Yj=0 by (3-47). On the other hand,

if p j < f * , we have Yj>0 and solves equation (3-50). Apply similar arguments to the case where p j = f * , we still get Yj=0. We complete part (ii) by rephrasing our arguments above. Q.E.D.

127

Appendix 8:

Proof of Corollary 3.5.3. Let {Xi*} and {Yj}* be an arbitrary solution to (3-49)-(3-50). From
?i ? s' (0) + ri < i +1 s' (1) + ri +1 , we can have: qi q i +1
s' (1) < r1 +

r1 +

?1
q1

?1
q1

s' (0) < r2 +
n

?2
q2
*

s' (1) < r2 +

?2
q2

s ' (0) < .... < rn +
n *

?n
qn

s' (1) < rn +

?n
qn

s ' ( 0)

For simplicity, let f (?) = g ' (? X i + ? Y j ) + m' (? X i + ? Y j )
i =1 j? i =1 j?

Without the loss of the generality, there are two cases: (A): if rk +

?k
qk

s ' (0) ? f (?) ? rk +1 +

? k +1
q k +1

s ' (1) Xi ) = ?i ? ? i > 0 , qi

then: 1)for i ? k + 1 : ri ? f (?) +

?
qi

s' (

i.e.: ?i > 0 , from ?i (qi ? X i ) = 0 , we have: X i = qi , ( i ? k + 1 ) 2) For all i ? k :
ri ? f (?) +

?i
qi

s' (

? Xi ) < ri ? f (?) + i s ' (0) < 0 qi qi ?
qi s' ( Xi ) < 0, qi

?

?i ? ? i = ri ? f (?) +

i.e. ? i > 0 , from: ? × X i = 0 , we have: X i = 0, ( i < k ) (B) if: rk +
?k ? s' (1) < g ' (?) + m' (?) < rk + k s' (0) qk qk

then: 1) for i ? k + 1 : ri ? f (?) +

?
qi

s' (

Xi ) = ?i ? ? i > 0 , qi

128

i.e.: ?i > 0 , from ?i (qi ? X i ) = 0 , we have: X i = qi , ( i ? k + 1 ) 2) For all i ? k ? 1 :
ri ? f (?) +
n n ? Xi * * ) < ri ? g ' (? X i + ? Y j ) ? m' (? X i + ? Y j ) + k s ' (1) < 0 ? qi qk qi i =1 j? i =1 j? ? X ?i ? ? i = ri ? f (?) + i s ' ( i ) < 0 , qi qi

?i

s' (

i.e. ? i > 0 , from: ? × X i = 0 , we have: X i = 0, ( i < k ) for i =k, need to be determined by (3-49) and (3-50). Combining both (A) and (B), the first part is proven. Now we know that

?X
i?I

i

? q * From part ii) of Theorem 3.4, Yj>0 only if

p j > g ' ( ? X i + ? Y j ) + m' (? X i + ? Y j )
i?I j?J i?I j?J

Since

g ' ' ( ? X i + ? Y j ) + m' ' ( ? X i + ? Y j ) > 0 ,
i?I j?J i?I j?J

we

know

g ' (? X i + ? Y j ) + m' (? X i + ? Y j ) is increasing.
i?I j?J i?I j?J

Thus we have
p j > g ' (? X i + ? Y j ) + m' (? X i + ? Y j ) > g ' (q*) + m' (q*)
i?I j?J i?I j?J

when Yj > 0. This is equivalent to the second part of the corollary. Q.E.D.

129

Appendix 9:

Proof of Theorem 4.1. Let the sequences (ai ? li ) non-decreasingly ordered. Without loss of generality, assume we have:
a1 ? l1 ? ... ? a k ? l k ? f ' (? Qi ) < a k +1 ? l k +1 ? ... ? a n ? l n
i?I

then for all i ? k , the formula (4.7) is strict negative, so (4.8) must hold, and we have
Qi = 0 ,

for all i ? k ; for all i > k , formula (4.7) = 0, so we have:

Qi =

a i ? l i ? f ' (? Qi )
i?I

2bi

for all i > k .

Q.E.D.

130

Appendix 10:

Proof of Remark 4.1. It suffices to show that the Hessian matrix of the objective function (4.5) is negative definitive.
?b1 ?0 ? H = ?2 × ? 0 ? ?0 ? ?0 0 b2 0 0 0 0 0 0 0 0? 0? ? 0 ? ? f ' ' (? Qi )I × I T ? i?I 0? bn ? ? < f ' ' (?) < 0 , then: X T HX < 0 .

. 0 0 bn ?1 0 0 2

We need to show: if ?

?b
i?I

1
i

X T HX = ?2? bi xi ? f ' ' (? Qi )(? xi )
2 2

i?I

i?I

i?I

According to Cauchy formula, we have: (? xi ) = ( ?
2

1
bi

i?I

i?I

? bi ? xi ) 2 ? (?
i?I
2

1
bi

)(? bi xi ) , then we easily have:
2

i?I

? 2? bi xi ? f ' ' (? Qi )(? xi ) ? 0
2

i?I

i?I

i?I

Q.E.D.

131

Appendix 11:

Proof of Remark 4.1
Wi (0) = ? ( S i ? u i ) f (u i )du i = 0
0 0

Wi ' (0) = F (0) = 0

and: Q Wi ' ( S i ) = F ( S i ) ? Wi " ( S i ) = f ( S i ) > 0 So the Wi (?) is a convex function. Q.E.D.

132

Appendix 12 :

Proof of Theorem 4.2. Let the sequences ( pi + d i ? l i ) non-decreasingly ordered. Without loss of generality, assume we have:
p1 + d1 ? l1 ? ... ? p k + d k ? l k ? g ' (? S i ) < p k +1 + d k +1 ? l k +1 ? ... ? p n + d n ? l n
i?I

then for all i ? k , the formula (4.12) is strict positive; to let (4.13) hold, we have
S i = 0 for all i ? k ;

On the other hand, for all i > k , when g ' (? Si ) < pk +1 + d k +1 ? lk +1 ? ... ? pn + d n ? ln ,
i?I

we have: Wi ( S i ) > 0 . We also know: Wi (?) is convex and Wi (0) = 0 , ? Si > 0 .
Q S i > 0 and ?i S i = 0

'

'

? ?i = ( pi + d i + hi )Wi ( S i ) + g ' (? S i ) ? ( pi + d i ? li ) = 0
'

i

so we have: Wi ' ( S i ) = F ( S i ) = Q.E.D.

( p i + d i ? l i ) ? g ' (? S i )
i?I

pi + d i + hi

for all i > k .

133

Appendix 13.

Q Wi (?) is convex function, ? Wi ' (?) is a increasing function which goes from 0

when S i =0 to 1 when S i = infinity
Q Wi ' (0) = 0 , and
Wi ' ( S i ) < ( pi + d i ? li ) <1 pi + d i + hi

? S i has upper bound. ? Q = ? S i has upper bound.
i

Q.E.D.

134

Appendix 14:

Proof of Remark 4.2. It suffices to show that the Hessian matrix of the objective function

?(p
i

i

+ hi + d i ) ? ( S i ? u i ) f (u i )du i + g (? S i ) ? ? ( pi + d i ? l i ) S i
0

Si

i

i

is positive definitive. We have: ?( p1 + h1 + d1 ) f ( S1 ) ? 0 ( p 2 + h2 ? H =? 0 ? 0 ? ? 0 ?
i?I

0 + d 2 ) f (S 2 ) 0 0 0

0 0 0 0 . 0 0 ( p n ?1 + hn ?1 + d n ?1 ) f ( S n ?1 ) 0 0 ( p n + hn

+ g ' ' (? S i )I × I T

0 ? ? 0 ? ? 0 ? 0 ? + d n ) f ( S n )? ?

Now we need to show: if

? (p
i

i

1 > UF × UG , then: X T HX > 0 ; where UF + hi + d i )

and UG is the upper bound of f (?) and ? g ' ' (?) respectively.
X T HX = ? ( pi + hi + d i ) f ( S i )xi + g ' ' (? S i )(? xi )
2 2

i?I

i?I

i?I

According to Cauchy formula, we have:
(? xi ) = (?
2 i?I i?I

1 ( pi + hi + d i ) f ( S i ) 1
i?I

? ( p i + hi + d i ) f ( S i ) ? xi ) 2

? (?
i?I

( pi + hi + d i ) f ( S i )

)(? ( pi + hi + d i ) f ( S i ) xi )
2

,

then we easily have:

?( p
i?I i?I

i

+ hi + d i ) f ( S i )xi > ? g ' ' (? S i )(? xi )
2 2 i?I i?I 2 2

? ? ( pi + hi + d i ) f ( S i )xi + g ' ' (? S i )(? xi ) > 0
i?I i?I

Q.E.D.

135

Appendix 15:

Proof of Remark 4.3. Proof: If there exists optimal solution for problem (4.9)-(4.10), then Hessian Matrix of the objective function

?( p
i

i

+ hi + d i ) ? ( S i ? u ) f (u )du + g (? S i ) ? ? ( pi + d i ? li ) S i
0 i i

Si

is positive definite, and its Hessian matrix can be written as:

? ( p1 + h1 + d1 ) f ( S1 ) ? 0 ( p 2 + h2 ? ? H= . ? 0 ? ? 0 ?

0 0 . + d 2 ) f (S 2 ) 0 . . 0 0

? ? ? ?+ . . . ? . . . ? ? 0 . ( p n + hn + d n ) f ( S n ) ?

0 0

?1 ? ?1 g ' ' (? S i ) ? . ? i ?1 ? ?1

1 1 . 1 1

1 1 . . 1

. . . . .

1? ? 1? .? ? 1. ? ? 1?

H being positive definite means, for any X = ( X 1 , X 2 ,... X n ) ? 0 , XHX T > 0 .

We have:
XHX T = ? ( pi + hi + d i ) f ( S i ) X i2 + g ' ' (? S i )(? X i ) 2
i i i

Using Cauchy inequality, we have:
136

(?
i

1 )(? ( pi + hi + d i ) f ( S i ) X i2 ) ? (? X i ) 2 ( pi + hi + d i ) f ( S i ) i i

So:
XHX T ? ( 1 1 ? i ( p i + hi + d i ) f ( S i ) + g ' ' (? S i ))(? X i ) 2
i i

And Positive Definite implies:

1 1 ? i ( p i + hi + d i ) f ( S i )

+ g ' ' (? S i ) > 0
i

In Theorem 4.2, we have optimal solution:
Si = 0

for all i ? k ,

(4.12)
( p i + d i ? l i ) ? g ' (? S i )
i?I

F (S i ) =

p i + d i + hi

for all i > k .

(4.13) So from (4.13), we have:

S i = F (1 ?

?1

hi + li + g ' (? S i )
i?I

pi + d i + hi

)

137

Q = ? Si =
i

i?K *

?F

?1

(1 ?

hi + li + g ' (Q) ) pi + d i + hi

K * is the set of i such that (4.13) holds.

Without loss of generality, we assume 1 ? K * , Now we look at

?Q : ?p1

Define function R (Q) as: R(Q) =

i?K *

? F ?1 (1 ?

hi + li + g ' (Q) )?Q pi + d i + hi

?R ?p ?Q ?R ?R = ? 1 , and let’s look as signs of and : ?R ?p1 ?p1 ?Q ?Q
h1 + l1 + g ' (Q) ( p1 + d1 + h1 ) 2 ?R = ?0 hi + li + g ' (Q) ?p1 ?1 ) f ( F (1 ? p i + d i + hi g ' ' (Q) pi + d i + hi ?1 hi + li + g ' (Q) ?1 ) f ( F (1 ? pi + d i + hi ?

?R =? ?Q i?K *

= ??

i?K *

g ' ' (Q) ? 1, f ( S i )( pi + d i + hi )

? ??

g ' ' (Q) ?1 1? i ? n f ( S i )( p i + d i + hi )

<0

138

?R ?p ?Q Therefore, = ? 1 ? 0. ?R ?p1 ?Q so Q = ? S i is increasing of pi .
i

Q.E.D.

139

Appendix 16.

Proof of Theorem 4.3. Proof: suppose the optimal solution is ( x1 , x 2 ,....x n ) , where x1 = 0 . We will construct a new solution ( x1 ' , x 2 ' ,....x n ' ) . If we can show this new solution is better than the old one, we complete our proof by contraposition. The construction of new solution is as follows: For 2 ? i ? n , we have: xi ' = 0 if xi = 0 ; xi ' = li if xi ? li , For i = 1 , we have: x1 ' = Q ?

1<i ? n

? x '.
i

Since we assume Q > ? li , it is trivial to show that x1 ' ? l1 .
i

Now we compare the objective value of the two solutions. Let z and z ' as the objective values associated with old solution and new solution respectively. We only need to show that: z ' < z . We have: z ? z ' = ? ( f (xi ) ? f ( xi ' ) + C i ( xi ? xi ' )) + B1 , (note the reason we have
i

time B1 is because x1 = 0 and x1 ' ? l1 ) We denote ? i = xi ? xi ' , and we have:

??
i

i

= 0,

?1 < 0 , )

? i ? 0 for 1 ? i ? n

140

Applying Mean-Value Theorem to the term f ( xi ) ? f ( xi ' ) yields:
z ? z ' = ? (? i ? i + C i ? i ) + B1
i

where ? i = Plugging
z ? z' =

f ( xi ) ? f ( xi ' )

?i
i

and 0 < f ' (Q) ? ? i ? f ' (0) .

??
i i

= 0 into the above formula, we have: ? C1 + ? i ? ?1 )? i + B1 ?

1<i ? n

? (C

1<i ? n

? (C

2

? C1 + f ' (Q) ? f ' (0))? i + B1 ? B1 > 0

Q.E.D.

141

Appendix 17

Proof of Theorem 4.4. suppose some xi , say x 2 > l 2 , we will construct a new solution which turns out to be better than the old one. Now we know x1 ? l1 from Theorem 1. To construct a new solution, we keep all xi except x1 and x 2 unchanged in the new solution. However, x 2 ' = l 2 , and x1 ' = x1 + x 2 ? l 2 . It is trivial to show the new solution is feasible. So we have:
z ? z ' = f ( x1 ) + f ( x 2 ) ? f ( x1 ' ) ? f ( x 2 ' ) + C1? 1 + C 2? 2

where 0 < ? 2 = ?? 1 , still applying Mean-Value Theorem, we have:
z ? z ' = (C 2 ? C1 + ? 2 ? ?1 )? 2 > 0

Q.E.D.

142

BIBLIOGRAPHY

[1] Agrawal, V. and Seshadri, S, 1999, “Impact of Uncertainty and Risk Aversion on Price and Order Quantity in the Newsvendor Problem”, Manufacturing &
Service Operations Management, 2, 410-423.

[2] Ardalan A., 1995, "A Comparative Analysis of Approaches for Determining Optimal Price and Order Quantity When a Sale Increases Demand", European
Journal of Operational Research, 84, 416–430.

[3] Balakrishnan, A. and J. Geunes, 2000, "Requirements Planning with Substitutions: Exploring Bill–of–Materials Flexibility in Production Planning",
Manufacturing & Service Operations Management, 2, 166–185.

[4] Ball, M. O., S. Boyson, C.–Y. Chen, E. Rabinovich, L. Raschid, V. Sambamurthy and V. Zadorozhny, 2001, "Scalable Supply Chain Infrastructures and Their Impact on Organizations", Proceedings of the 2001 NSF Design and
Manufacturing Research, Tampa.

[5] Ball, M. O., S. Boyson, L. Raschid and V. Sambamurthy, 2000, "Scalable Supply Chain Infrastructures: Models and Analysis", Proceedings of the 2000 NSF Design
& Manufacturing Research Conference, Vancouver.

[6] Ball, M.O., Chen, C-Y. and Zhao, Z., 2003, “Material Compatibility Constraints for Make-To-Order Production Planning”, Operations Research Letters , 31 , 420428.

143

[7] Bassok, Y. and R. Anupindi, 1997, "Analysis of Supply Contracts With Total Minimum Commitment", IIE Transactions, 29, 373–381. [8] Baumol, W.J. and Wolf, P, 1958, “A Warehouse Location Problem”, Operations
Research, 6, 252-8.

[9] Belobaba, Peter, 1987, "Airline Yield Management: An Overview of Seat Inventory Control", Transportation Science, 21, 63-73. [10] Belobaba, Peter,1989, ``Application of a Probabilistic Decision Model to Airline Seat Inventory Control''. Operations Research, 37, 183-197. [11] Bhattacharjee, S. and R. Ramesh, 2000, "A Multi–Period Profit Maximizing Model for Retail Supply Chain Management: An Integration of Demand and Supply–Side Mechanisms", European Journal of Operational Research, 122, 584–601. [12] Bitran, G. and S. Mondschein, 1995, "An Application of Yield Management to the Hotel Industry Considering Multiple Day Stays", Operations Research, 43, 427- 448.
[13] Bitran, G., R. Caldentey, S. Mondschein, 1998, "Coordinating Clearance Markdown Sales of Seasonal Products in Retail Chains", Operations

Research,
46, 609--624.

[14] Blanchini, Franco, Stefano Miani and Franca Rinaldi, 2004, “Guaranteed Cost Control for Multi-inventory Systems with Uncertain Demand”, Automatica,

144

40, 213-223 [15] Cattani, K, 1998, “Demand Pooling Effects of a Universal Product when Demand is from Distinct Markets”, European Journal of Operations
Research, 135, 582-601.

[16] Cattani, K. and G. C. Souza (2002), "Inventory Rationing and Shipment Flexibility Alternatives for Direct Market Firms," Production and Operations
Management, 11, 4.

[17] Chen, C-Y. , Zhao, Z. and Ball, M.O., 2002, “A Model for Batch Available-toPromise”, Production and Operations Management, 11, 424 - 440 [18] Chen, C-Y. , Zhao, Z. and Ball, M.O, 2001, “Quantity and Due Date Promising Available to Promise”, Information Systems Frontiers, 3 , 477-488. [19] Chen, F., A. Federgruen and Y.–S. Zheng, 2001, "Near–Optimal Pricing and Replenishment Strategies for a Retail/Distribution System," Operations
Research, 49, 839–853.

[20] Chen, F. Y. and D. Krass, 2001, "Analysis of Supply Contracts with Minimum Total Order Quantity Commitments and Non–Stationary Demands," European
Journal of Operational Research, 131, 302–233.

[21] Chen, F., 2001, "Market Segmentation, Advanced Demand Information, and Supply Chain Performance," Manufacturing & Service Operations Management, 3, 53–67. [22] Cheng, F. and M.Ettl, 2002, “Inventory-service Optimization in Configure-ToOrder Systems”, Manufacturing & Service Operations Management, 4, 114126

145

[23] Cohen, M. and H. L. Lee, 1989, "Resource Deployment Analysis of Global Manufacturing and Distribution Networks," Journal of Manufacturing and
Operations Management, 2, 81–104.

[24] Cplex, 1998, "Using The Cplex Callable Library", Version 6.0. Incline Village, NV, ILOG, Inc. [25] Das, S. and L.Abdel-Malek, 2005, “Modeling the flexibility of order quantities and lead-times in supply chains”, International Journal of Production
Economics, 85, 171-181

[26] Deshpande, V., M. Cohen and K. Donohue, 2003, "A Threshold Inventory Rationing Policy for Service Differentiated Demand Classes," Management
Science, 49, 683-703

[27] eB2x, I., 2000, "Why Promising Orders? How Strategic Order Promising Helps Supply Chains Compete in the New Economy," Online, http://www.eb2x.com. [28] Federgruen, A. and A. Heching, 1999, "Combined Pricing and Inventory Control under Uncertainty," Operations Research, 47, 454–475. [29] Feldman, E., Lehrer, F.A. and Ray, T.L., 1966, “Warehouse Location Under Continuous Economies of Scale”, Management Science, 12, 670-84.
[30] Feng, Y., G. Gallego., 2000, "Perishable Asset Revenue Management with Markovian Time Dependent Demand Intensities", Management

Science, 46,
941--956.

146

[31] Fordyce, K. and G. Sullivan, 1999, "Supply Chain Management, Decision Technology, and E–Business Information Technology at IBM Microelectronics,"
MicroNews (http://chips.ibm.com/micronews), 5, 4.

[32] Galeotti, Marzio and Louis J. Maccini, Fabio Schiantarelli, 2005, “Inventories, Employment and Hours”, Journal of Monetary Economics,138, 72-85
[33] Gallego, G., G. van Ryzin, 1997, "A Multi-product, Multi-resource Pricing Problem and Its Applications to Network Yield Management",

Operations Research, 45, 24--41.
[34] Ghali, M., 2003, "Production-planning Horizon, Production Smoothing, and Convexity of the Cost Functions", International Journal of Production Economics, 81, 67 –74. [35] Godfrey, G. and Powell, W., 2002, "An Adaptive, Approximate Method for Stochastic Dynamic Programs with Applications to Inventory and Distribution Problems", Transportation Science 36, 21-39. [36] Goetschalckx, M., C. J. Vidal and K. Dogan, 2002, "Modeling and Design of Global Logistics Systems: A Review of Integrated Strategic and Tactical Model and Design Algorithms," European Journal of Operational Research, 143, 1–18. [37] Goh, M. and M. Sharafali, 2002, "Price–Dependent Inventory Models with Discount Offers Random Times", Production and Operations Management, 11, 139–156.

147

[38] Grant, Hank and Scott Moses, Dave Goldsman, 2002, "Using Simulation to
Adjust Buffer In Order Promising’’, Proceedings of the 2002 Winter Simulation Conference

[39] Ha, A., 1997, "Inventory Rationing in a Make–to–Stock Production System with Several Demand Classes and Lost Sales", Management Science, 43, 1093–1103. [40] Hall, George, 2002, "Non-Convex Costs and Capital Utilization: A Study of Production Scheduling at Automobile Assembly Plants", Journal of Monetary
Economics, 47, 113-128.

[41] Horst, R., Tuy, H., 1993, Global Optimization, second edition. Springer, Berlin. [42] Huchzermeier, A. and M. Cohen, 1996, "Valuing Operational Flexibility under Exchange Rate Risk", Operations Research, 44, 100–113. [43] Hunter, Lisa M., 2004, “A Classification of Business-to-Business Buying Decisions: Risk Importance and Probability as a Framework for e-business Benefits”, Industrial Marketing Management, 33,145-154. [44] Johnson, L. and D. Montgomery, 1974, Operations Research in Production
Planning, Scheduling and Control, New York, NY, John Wiley & Sons, Inc.

[45] Kaplan, A., 1969, "Stock Rationing", Management Science, 5, 260–267. [46] Kelle, P. and E. Silver, 1990, "Safety Stock Reduction by Order Splitting",
Naval Research Logistics, 37, 725–743.

[47] Khouja, Moutaz J., 2000, “Optimal Ordering, Discounting, and Pricing in the Single-period Problem”, International Journal of Production Economics, 65, 201-216

148

[48] Kilger, C. and L. Schneeweiss, 2000, "Demand Fulfillment and ATP. Supply Chain Management and Advanced Planning: Concepts, Models, Software and Case Studies", H. Stadtler and C. Kilger. Berlin, Germany, Springer, 135–148. [49] Kritchanchaj, D. and B.L. MacCarthy, 1999, “Responsiveness of the Order Fulfillment Process", International Journal Of Operations & Production
Management, 19, 145-154

[50] Lau, H and A.Lau, 1995, “Multiple-Product, Multiple-Constraint Newsboy Problem: Formulation, Application, and Solution”, Journal of Operations
Management 13, 153-162.

[51] Lau, H. and L. Zhao, 1993, "Optimal Ordering Policies with Two Suppliers When Lead Times and Demands Are All Stochastic", European Journal of
Operational Research, 68, 120–133.

[52] Lau, H and Hon-Shiang Lau, 2003, “Some Two-echelon Supply-chain Games: Improving from Deterministic-symmetric-information to Stochastic-asymmetricInformation Models”, European Journal of Operational Research, 147, 530-548. [53] Lee, H. L. and C. Billington, 1995, "The Evolution of Supply Chain Management Models", Interfaces, 25, 43–63. [54] Lin, Chinho and Yihsu Lin, 2005, “The Production Size and Inventory Policy for A Manufacturer in a Two-echelon Inventory Model”, Computers & Operations
Research, 32, 1181-1196

[55] Littlewood, K., 1972, ``Forecasting and Control of Passenger Bookings'',
AGIFORS Symposium Proceedings, American Airlines, New York, 95--117.

[56] Lovell, C. Michael, 2003, “Optimal Lot Size, Inventories, Prices and JIT under

149

Monopolistic Competition”, International Journal of Production Economics, 81, 59-66. [57] Meixell, M and Chien-Yu Chen, 2004, “A Scenario-Based Bayesian Forecasting Engine for Stochastic Available-to-Promise”, 15th Annual Production and
Operations Management Conference

[58] Moinzadeh, K.,1997, "Replenishment and Stock Policies for Inventory Systems With Random Deal Offerings," Management Science, 43, 334–342. [59] Morgan, L. O., R. L. Daniels and P. Kouvelis, 2001, "Marketing/Manufacturing Trade–Offs in Product Line Management", IIE Transactions, 33, 949–962. [60] Moses, S., H. Grant, L. Gruenwald and S. Pulat, 2002, "Real–Time Due Date Promising by Build–to–Order Environments", in Proceedings of the 2002 NSF
Design and Manufacturing Research, San Juan, Puerto Rico.

[61] Nagruney, Annua and Jose Cruz, 2005, “Supply Chain Networks, Electronic Commerce, and Supply Side and Demand Side Risk”, European Journal Of
Operations Research, 26, 120-142

[62] Pardalos, P.M., Rosen, J.B., 1984, "Constrained Global Optimization: Algorithms and Applications", Computer Science, 11, 268, [63] Petruzzi, N. C. and M. Dada, 1999, "Pricing and the Newsvendor Problem: a Review with Extensions," Operations Research, 47, 183–194. [64] Ramasesh, R. V., J. Ord and J. Hayya, 1993, "Note: Dual Sourcing with Nonidentical Suppliers," Naval Research Logistics, 40, 279–288. [65] Ramasesh, R. V., J. Ord, J. Hayya and A. Pan, 1991, "Sole Versus Dual Sourcing In Stochastic Lead–Time (S,Q) Inventory Models," Management Science, 37,

150

428–443. [66] Robinson, A. and D. Dilts, 1999, "OR & ERP: A Match for the New Millennium?," OR/MS Today, 26, 30–35. [67] Ryan, S., 2002, "Planning for and Responding to Changing Demand", in
Proceedings of the 2002 NSF Design and Manufacturing Research Conference,

San Juan, Puerto Rico.
[68] Serel, D. A., M. Dada and H. Moskowitz, 2001, "Sourcing Decisions with

Capacity Reservation Contracts," European Journal of Operational Research, 131, 635–648. [69] Smith, S. and N. Agrawal, 2000, "Management of Multi–Item Retail Inventory Systems with Demand Substitution," Operations Research, 48, 1, 50–64. [70] Souza, G., Z. Zhao, M. Chen, and M. Ball, 2004, “Coordinating Sales and Raw Material Discounts in a Global Supply Chain,” Production and Operations
Management, 13 , 34-45.

[71] Subramanian, J., S. Stidham and C. Lautenbacher, 1999, "Airline Yield Management with Overbooking, Cancellations, and No–Shows," Transportation
Science, 33, 147–167.

[72] Taylor, S. and G. Plenert, 1999, "Finite Capacity Promising," Production &
Inventory Management, 40, 3, 50–56.

[73] Tayur, S. and Ganeshan R., 1998, “Quantityative Models for Supply Chain
Management”, Kluwer Academic Publishers.

[74] Tersine, R. J. and S. Barman, 1995, "Economic Purchasing Strategies for Temporary Price Discounts," European Journal of Operational Research, 80,
151

328–343 [75] Thomas, D. J. and P. M. Griffin, 1996, "Coordinated Supply Chain Management," European Journal of Operational Research, 94, 1–15. [76] Topkis, D., 1968, "Optimal Ordering and Rationing Policies in a Nonstationary Dynamic Inventory Model with N Classes", Management Science, 15, 160– 176. [77] Urban, T., 2000, "Supply Contracts with Periodic, Stationary Commitment,"
Production and Operations Management, 9, 400–413.

[78] Van Ryzin, G. and S. Mahajan, 1999, "On the Relationship Between Inventory Costs and Variety Benefits in Retail Assortments," Management Science, 45, 1496–1509. [79] Vidal, C. J. and M. Goetschalckx, 1997, "Strategic Production–Distribution Models: A Critical Review with Emphasis on Global Supply Chain Models,"
European Journal of Operational Research, 98, 1–18.

[80] Vidal, C. J. and M. Goetschalckx, 2001, "A Global Supply Chain Model with Transfer Pricing and Transportation Cost Allocation," European Journal of
Operational Research, 129, 134–158.

[81] Vollman, T. E., W. L. Berry and D. C. Whybark, 1992, Manufacturing Planning
and Control Systems, Homewood, IL, Business One Irwin.

[82] Wagner, H.,1995, "Global Sensitivity Analysis," Operations Research, 43, 948– 969.
[83] L.R. Weatherford, S.E. Bodily and P.E. Pfeifer., 1993, "Modeling the Customer

152

Arrival Process and Comparing Decision Rules in Perishable Asset Revenue Management Situations", Transportation Science 27, 239-251.

[84] Weber, Charles A.; Current, John; Desai, Anand, 2000, “An Optimization Approach to Determining the Number of Vendors to Employ”, Supply Chain
Management: An International Journal, 5, 90-98.

[85] Weng, Z. Kevin, 1995, “Coordinating Order Quantities Between the Manufacturer and the Buyer: A Generalized Newsvendor Model”, European
Journal of Operational Research, 41, 1509-1522.
[86] Williamson, E. L. 1992, "Airline Network Seat Control", Ph.D. thesis, MIT, Cambridge, Massachusetts.

[87] Yano, C. and G. Dobson, 1998, "Profit–Optimizing Product Line Design, Selection and Pricing with Manufacturing Cost Considerations," Product
Variety Management: Research Advances, T. H. Ho and C. Tang. Norwell,

MA, Kluwer, 145–176. [88] Yongjin, K., M. Hani and J. Patrick, 2002, "Dynamic Truckload Truck Routing and Scheduling in over–Saturated Demand Situations," Proceedings of the 2002
NSF Design and Manufacturing Research Conference, San Juan, Puerto Rico.

[89] Youssef, M.,1996, " An Iterative Uncapacitated Production-Distribution Problem under Concave Cost Function", International Journal of Operations &
Production Management, 16, 18-27.

153

[90] Zhao, L. and H. Lau, 1992, "Reducing Inventory Costs and Choosing Suppliers with Order Splitting," Journal of the Operational Research Society, 43, 10031017. [91] Zhao.W and Y. Wang, 2002, "Coordination of Joint Pricing–Production Decisions in a Supply Chain," IIE Transactions, 34, 701–715 [92] Zhao, Z., M. Ball and M. Kotake, 2005, “Optimization-Based Available-toPromise and Multi-Stage Resource Availability”, Annals of Operations
Research, 135, 65 – 85.

[93] Zweben, M., 1996, "Delivering on Every Promise: Real–Time Sales Order Promising of Supply Chain Resources," APICS: The Performance Advantage.

154



doc_621432060.pdf
 

Attachments

Back
Top