OPRE 6366. SCM: 3. Aggregate Planning
OPRE 6366. SCM: 3. Aggregate Planning
Aggregate Planning
1. MPS
MPS stands for Master Production Schedule. It is a company’s plan of how many and when items will
be delivered to the customer. MPS includes deliveries that are fixed by the customer and promised by
the company. It can include deliveries that are not yet promised but that are still negotiated. It may also
include a forecast of customer demand. MPS covers a planning horizon of several months into the future.
Initially a bigger portion of MPS is constituted by the firm customer orders and promised deliveries. But
towards the end of the planning horizon, a larger portion of MPS is driven by demand forecasts.
Since firm orders cannot be altered by the customer and the company tends to hold delivery promises,
initial months of MPS is fairly stable. Because of this, the portion of MPS associated with those months is
said to be frozen. Outside the frozen horizon, orders can be modified relatively easily so MPS is flexible; see
Figure 1.
Volume
Forecasts
Firm
Orders
Firm Orders
Having fixed the demand during the frozen horizon, companies use linear and integer programming
formulations to make production/distribution plans. On the other hand, demands are random during
the flexible zone. There are ways of incorporating randomness into linear and integer programs or ways
to create models from scratch. What is important is to recognize the random demand during the flexible
horizon. Failing to do so, i.e., assuming pseudocertainty, will result in plans that are too adamant to adjust
for various demand scenarios that can materialize.
Demand forecasting is naturally an important ingredient in MPS construction, especially during the
flexible horizon. We will not discuss forecasting methods. We want to emphasize that there are non-
1
traditional forecasting models (Guerrero and Elizondo [3], and Bodily and Freeland [2]) where demand is
assumed to reveal itself in steps over time. Partial and earlier observations of demand is used to forecast
the future demand later on. For example, in one of the models it is assumed that the ratio of orders received
up to a certain time to the whole demand is approximately constant. After the proportionality constant is
estimated, forecasts are readily generated from the partially observed demands.
We also note that forecasts for immediate future are more accurate (in terms of less randomness or
lower variance) than those for the far future. Thus, planners should handle forecasts of future demands
with some suspicion. If possible, it is a smart strategy to wait for forecasts to become more accurate before
committing to meet demands. An example of this strategy is the postponement of product differentiation.
Demand forecasts are the drivers of the SC operations, once they are ready we can plan the operations as
discussed in the next section.
These strategies are extreme although we often resort to compromises among the extremes. To achieve
such compromises, we need to detail the strategies further, which can be done through a more quantitative
approach.
• The total demand for the product in month j is d j , for j = 1, 2, . . . , 12. These could either be targeted
values or be based on forecasts.
• The cost of producing each unit of the product in month j is c j (dollars), for j = 1, 2, . . . , 12. There is
no setup/fixed cost for production.
• The inventory holding cost per unit for month j is h j (dollars), for j = 1, 2, . . . , 12. These are incurred
at the end of each month.
• The production capacity for month j is m j , for j = 1, 2, . . . , 12.
The manager’s task is to generate a production schedule that minimizes the total production and inventory-
holding costs over this twelve-month planning horizon.
To facilitate the formulation of a linear program, the manager decides to make the following simplifying
assumptions for now:
2
1. There is no initial inventory at the beginning of the first month.
2. Units scheduled for production in month j are immediately available for delivery at the beginning of
that month. This means in effect that the production rate is infinite.
3. Shortage of the product is not allowed at the end of any month.
To understand things better, let us consider the first month. Suppose, for that month, the planned pro-
duction level equals 100 units and the demand, d1 , equals 60 units. Then, since the initial inventory is 0
(Assumption 1), the ending inventory level for the first month would be 0+100-60=40 units. Note that all
100 units are immediately available for delivery (Assumption 2); and that given d1 = 60, one must produce
no less than 60 units in the first month, to avoid shortage (Assumption 3). Suppose further that c1 = 15
and h1 = 3. Then, the total cost for the first month can be computed as: 15 · 100 + 3 · 40 = 1380 dollars.
At the start of the second month, there would be 40 units of the product in inventory, and the cor-
responding ending inventory can be computed similarly, based on the initial inventory, the scheduled
production level, and the total demand for that month. The same scheme is then repeated until the end of
the entire planning horizon.
• The total production cost for the entire planning horizon equals
12
∑ c j x j ≡ c1 x1 + c2 x2 + · · · + c12 x12 ,
j =1
where we have introduced the standard summation notation (“≡” means by definition).
3
• The total inventory-holding cost for the entire planning horizon equals
[ ] [ ] [ ]
12 j 1 2
∑ h j ∑ ( xk − dk ) ≡ h1 ∑ ( xk − dk ) + h2 ∑ ( xk − dk ) +...
j =1 k =1 k =1 k =1
[ ]
12
+h12 ∑ ( xk − dk )
k =1
= h1 [ x1 − d1 ] + h2 [( x1 − d1 ) + ( x2 − d2 )] + . . .
+h12 [( x1 − d1 ) + ( x2 − d2 ) + · · · + ( x12 − d12 )] .
Since our goal is to minimize the total production and inventory-holding costs, the objective function can
now be stated as [ ]
12 12 j
Min ∑ c j x j + ∑ h j ∑ ( xk − dk ) .
j =1 j =1 k =1
xj ≤ mj
for j = 1, 2, . . . , 12; and since shortage is not allowed (Assumption 3), we require
j
∑ ( xk − dk ) ≥ 0
k =1
for j = 1, 2, . . . , 12. This results in a set of 24 functional constraints. Of course, being production levels, the
x j ’s should be nonnegative.
3.4 LP Formulation
In summary, we have arrived at the following formulation:
[ ]
12 12 j
Min ∑ c j x j + ∑ h j ∑ ( xk − dk )
j =1 j =1 k =1
Subject to :
xj ≤ mj f or j = 1, 2, . . . , 12
j
∑ ( xk − dk ) ≥ 0 f or j = 1, 2, . . . , 12
k =1
xj ≥ 0 f or j = 1, 2, . . . , 12 .
This is a linear program with 12 decision variables, 24 functional constraints, and 12 nonnegativity con-
straints. In an actual implementation, we need to replace the c j ’s, the h j ’s, the d j ’s, and the m j ’s with
explicit numerical values.
4
3.5 An Alternative Formulation for Production Planning
In the above formulation, the expression for the total inventory-holding cost in the objective function in-
volves a nested sum, which is rather complicated. Notice that for any given j, the inner sum in that expres-
j
sion, ∑k=1 ( xk − dk ), is simply the ending inventory level for month j. This motivates the introduction of an
additional set of decision variables to represent the ending inventory levels. Specifically, let
then, the objective function can be rewritten in the following simpler-looking form:
12 12
Min ∑ cj xj + ∑ hj yj .
j =1 j =1
With these new variables, the no-shortage constraints also simplify to y j ≥ 0 for j = 1, 2, . . . , 12. However,
we now need to introduce a new set of constraints to “link” the x j ’s and the y j ’s together.
Consider the first month again. Denote the initial inventory level as y0 ; then, by assumption, we have
y0 = 0. Since the production level is x1 and the demand is d1 for this month, we have y1 = y0 + x1 − d1 .
Continuation of this argument shows that for j = 1 , 2, . . . , 12,
y j = y j −1 + x j − d j ;
and these relations should appear as constraints to ensure that the y j ’s indeed represent ending inventory
levels. We have, therefore, arrived at the following new formulation:
12 12
Min ∑ cj xj + ∑ hj yj
j =1 j =1
Subject to :
xj ≤ mj f or j = 1, 2, . . . , 12
yj = y j −1 + x j − d j f or j = 1, 2, . . . , 12
xj ≥ 0 f or j = 1, 2, . . . , 12
yj ≥ 0 f or j = 1, 2, . . . , 12.
which is a linear program with 24 decision variables, 24 functional constraints, and 24 nonnegative vari-
ables.
Although there are twice as many decision variables in the new formulation, both formulations have
the same number of functional constraints. We will show in a later section that the total amount of ef-
fort necessary to arrive at an optimal solution to a linear program depends primarily on the number of
functional constraints. In general, it is not uncommon to have several equivalent formulations of the same
problem.
3.6 Remarks
• If Assumption 1 is relaxed, so that the initial inventory level is not necessarily zero, we can simply
set y0 to whatever given value.
• In our formulation, we assumed that there is no production delay (Assumption 2). This assumption
can be easily relaxed. Suppose instead there is a production delay of one month; that is, the scheduled
production for month j, x j , is available only after a delay of one month, i.e., in month j + 1. Then,
5
in the alternative formulation, we can simply replace the constraint y j = y j−1 + x j − d j by y j =
y j−1 + x j−1 − d j (with x0 ≡ 0), for j = 1, 2, . . . , 12. Of course, for the first month, the given value of y0
must be no less than d1 ; otherwise, the resulting LP will not have any solution.
• Assumption 3 can also be relaxed. If shortages are allowed, we can let yt be negative as well. Then
we define positive part of y j as on-hand inventory Ij and negative part as backorder S j :
Ij = max{0, y j }, S j = max{0, −y j }
so that
y j = Ij − S j
while Ij , S j ≥ 0. We also need to introduce a backorder penalty cost of, say, p j per unit of shortage at
the end of month j.
Solution: Let y BP = 1 if Betsy is produced at Plano, 0 otherwise. Define binary variables y BR , yVP , yVR ,
yWP , yWR similarly. Let x BP be the number Betsy produced at Plano. Define variables x BR , xVP , xVR , xWP , xWR
similarly.
6
yVP = 1 ⇒ yVR = 0, y BR = 0, yWR = 0; yWP = 1 ⇒ yWR = 0, y BR = 0, yVR = 0.
x BP /50 + xVP /80 + xWP /60 ≤ 300; x BR /40 + xVR /60 + xWR /70 ≤ 200.
These constraints say that no more than available hours used for production.
x BP ≤ 300(50)y BP ; xVP ≤ 300(80)yVP ; xWP ≤ 300(60)yWP ;
x BR ≤ 200(40)y BR ; xVR ≤ 200(60)yVR ; xWR ≤ 200(70)yWR ;
These constraints say that production can take place if the factory is set up.
Capacity planning models attempt to answer this question by telling the size of capacity increments and
the time between the increments.
It is hard to forecast demand with long planning horizons. Remember the forecasting discussion and
the fact that forecast become more variable as they are made for further periods in to the future. Tradition-
ally this variability is studied by treating demands as random variables. Then the stochastic optimization
techniques become the venue of choice. Because of the complexity of these techniques, we will restrict our
scope to deterministic capacity planning. Even within the deterministic capacity planning stream, there are
many different approaches. We next look at what criteria could be useful to classify deterministic capacity
planning approaches:
1. Single vs. Multiple facilities. When the capacity of a facility is planned independent of other existing
facilities such that facilities cannot transfer products among themselves, the problem is said to be
single facility type.
2. Single vs. Multiple resource capacities. Almost all operations require multiple resources (various
machines, workforce, etc.), but one can adapt an OPT-like approach and focus only on the single
resource. Then the capacity of the bottleneck resource is planned and the remaining resources are
adjusted according to the bottleneck resource capacity. This is typically how single resource capacity
models are used. In the multiple resource case, each resource’s capacity is studied explicitly and
optimized jointly. In this case, the sequence of capacity increments (increment first resource A, second
resource B, third resource A, fourth resource C, . . . ) needs to be decided on as well as the size of
increments and the time between increments.
3. Single vs. Multiple demands. In aggregate planning, demands of several products are lumped to-
gether to obtain a single demand stream. This is of course an approximation. The exact way of
handling multiple demands is representing them as a multidimensional vector. Then the capacity is
planned to match demand in every dimension. Clearly handling multiple demand streams explicitly
is a big challenge and there are very few models in this category.
7
4. Expansion vs. Contraction. Many models consider only capacity expansions. Capacity contractions
are not financially important, especially when machine capacities are concerned. This is because, ma-
chines tend to fill their economic life (amortized value of 0 dollars) before contraction, so companies
do not care to deal with the issue of getting rid of these no-dollar-value machines. The situation is
different for labor, because as long as the labor is employed there is an associated cost. Therefore,
contraction needs to be studied more carefully when studying labor capacity.
5. Discrete vs. Continuous expansion times. The choice here is made for modelling convenience. De-
pending on the structure, sometimes discrete models are easier to handle, and other times continuous
ones. The quality of solution whether the expansion times are discrete or continuous does not change
much. Sometimes vendors accept/deliver orders on certain times, say every monday, such cases give
rise to discrete times. But these cases can be studied almost equivalently well with continuous times.
6. Discrete vs. Continuous capacity increments. Almost all capacities can be incremented only in dis-
crete sizes. We can not buy half a machine or tenth of a truck. Thus, it will not be wrong to say that
capacity comes in quantum. However, sometimes the scale of increment is so large that capacity can
be treated as continuous. This is generally the case in workforce planning. The difference between
say 15 workers and 15.5 is not much. The justification of treating capacity as a continuous variable is
the same as the justification of using continuous variables for intrinsically-discrete activities in linear
programming.
7. Capacity costs. In most of the models, resources are priced according to their type and the increment
size. Naturally 1 lathe costs different than 1 press. The more interesting concept is the economies of
scale: 2 lathes cost less than twice the price of 1 lathe. When there is economies of scale, costs are
represented as increasing and concave costs. There can also be fixed costs associated with capacity
expansions. For example, when a building is expanded a big portion of the cost is independent of the
size of the expansion and it is fixed.
8. Penalty for demand-capacity mismatch. Many models plans capacity such that demand is always
met. However sometimes it is more profitable to delay an expansion (to save from the opportunity
cost of the expansion) while falling short of the demand. In order to evaluate the cost of this option
properly, we must charge ourselves when shortages happen. Shortage costs are generally propor-
tional to the magnitude of shortages. Penalty costs for not meeting demands are similar to those of
inventory.
9. Single vs. Multiple decision makers. In competitive markets, companies pay attention to each other’s
capacity expansions to avoid too much slack capacity. Then, expansion decisions become interdepen-
dent across companies. There are models that use game-theoretic techniques to study multi-player
markets with or without information symmetries.
A Model of Constant Capacity Increments: We now study a simple capacity expansion model where
capacity increments are all constant and equal to x. We keep the deterministic demand assumption and let
demand be D (t) = µ + δt. Here µ is the demand at the present time t = 0 and δ is the rate of demand
increase. Note that such an increasing trend in the demand cannot happen forever but for long enough
periods in specific industries. Examples can include global PC demand in the 1990s, global cell phone
demand in the 2000s, electricity demand in developing countries in the 2010s.
We assume that capacity is determined in such a way that the demand is always met. This assump-
tion applies to industries where storage is not possible/economical, such as power (electricity generation)
industry. Since capacity increments of size x are depleted in x/δ time units, there will be a new capacity
increment every x/δ time units. Thus, x/δ is the time between successive expansions. See Figure 2.
8
Units
Capacity
x
Demand
D(t)= µ+δ t
x
x/ δ Time (t)
Let f ( x ) be the cost of an expansion of size x. In general, expansion costs have economies of scale and
f ( x ) is an increasing but a concave curve. Let r be the interest rate (say, about 5%), i.e., discount rate of
money. Then the present value of the cost of an expansion of size x made t time units later is exp(−rt) f ( x ).
There is a brief discussion of continuous discounting and ”exp” in the Appendix.
Let C ( x ) be the discounted cost of expansions of size x over an infinite horizon:
where term exp(−rkx/δ) f ( x ) is the present value of the k th expansion. We can write the present value of
the total cost as
∞
1
C ( x ) = f ( x ) ∑ (exp(−rx/δ))k = f ( x ) ,
k =0
1 − exp(−rx/δ)
where we have used the identity
∞
1
∑ ak = 1 − a for a < 1.
k =0
For illustration we choose f ( x ) = x0.5 dollars, r=5%, δ=1 unit per week, and draw C ( x ) in Figure 3.
The minimum cost is achieved around the optimal expansion size x ∗ = 30 units. That is every expansion
must be of the size equivalent to 30 weeks demand. The total discounted cost of implementing x ∗ is about
7 dollars. In general, we suggest that C ( x ) graphed and x ∗ is found from the graph. This will be a safe way
to find x ∗ especially when C ( x ) is not convex. The details of this simple capacity expansion model can be
found in [5].
An important assumption in the previous model was 100% demand satisfaction. Capacity is often ex-
pensive and companies may choose to pay for shortage costs rather than install capacity. Once shortage
9
Figure 3: Discounted expansion cost C ( x ) as x varies.
becomes substantial, capacity is installed. Right before capacity installation, companies may use subcon-
tracting to boost up their capacity. Inventories can be used to regulate capacity and demand mismatch.
Inventories can be built up when there is surplus capacity. These inventoried are later depleted when
demand grows above capacity. We summarize these observations in Figure 4.
Units Capacity
Demand
10
5. Stochastic Capacity Planning: The case of flexible capacity
Many manufacturing companies own several plants where they produce various products. But not all
products can be produced at all plants. In order to produce a product at a plant, special tooling may need
to be installed at the plant. A plant that can produce many products is called flexible. To be accurate, it
is said to have product variety flexibility (compare this concept to product mix and volume flexibility).
Flexible plants provide a richer set of options when allocating products to plants.
First consider a small example with 3 plants and 2 products A and B. Suppose that Plant 1 and 2 pro-
duces only A, Plant 3 produces only B. Neither of these plants are flexible. Suppose further that we retool
Plant 2 so that it can produce B as well. With this added flexibility, we may produce all A or all B’s at
Plant 2. However, most often the production mixture is such that some As are produced along side with
some Bs at Plant 2. The more interesting issue is that the actual mixture of production at Plant 2 can be
determined after observing the actual demand. For example when demand for A is high, Plant 2 produces
more As and fewer Bs, and vice versa. Note that in both cases, the total production at Plant 2 is stable
although individual product volumes vary. When individual products are substitutes for each other (i.e.
they have negatively correlated demand), variations in A and B demand will cancel each other, keeping
the total demand constant. Therefore, the resulting system is more manageable.
In order to appreciate the flexibility of Plant 2, we must be explicitly model the correlation between
product A and B demand. Note that correlations are ignored when decisions are based on expected rev-
enues / profits. As a result, flexibility will have limited value if only discounted expected cash flows are
studied after ignoring demand correlations. The flexibility will be most valuable when the correlation is
strong but negative. By now, it must be clear that demands must be treated as random variables.
As for the decision variables, we have the plant to product assignment binary variable yij and production
levels xij . Then the objective function immediately follows
11
Suppose that we know the demands in advance as D j = d j (Capital letters used for random variables,
small letters for observations from associated random variables). Then we can add the demand constraint
as
∑ xij ≤ d j for each product j
i
However demands are not known in advance in general. Product to plant assignment decisions yij are
strategical so they are made against uncertain demand. On the contrary, production levels xij are deter-
mined more frequently and we can safely assume that they are determined after observing the demand.
In summary, the chronology of our decision making and demand observation is as: Decide on yij , Observe
D j = d j and Decide on xij .
Due to the specific sequence of decisions and demand observation we adopted, yij must be indepen-
dent of demand while xij can depend on the demand. To make ideas concrete, suppose that demands are
likely to be either d1j or d2j . More simplistically, we may say that the observed demand depends on a coin
toss after choosing yij . If head comes up, demand will be d1j , otherwise d2j . However, whether head or
tail comes up, yij will be the same: yij cannot anticipate the outcome of the coin toss. This limitation is
very common in decision making under uncertainty and leads to Nonanticipatory decision variables, e.g.
yij . On the other hand, xij is chosen after the demand so it is not need to be nonanticapatory: For differ-
ent demand observations, different xij can be used. For example if head comes up, we use xij1 , otherwise xij2 .
Each outcome of the demand is named as a scenario. Suppose that scenarios are indexed with k and
use pk to denote the probability of scenario k. Let us see now how the formulation changes
xijk ≤ ri yij for each plant and product i, j and for each scenario k
∑ xijk ≤ dkj for each product j and for each scenario k
i
By introducing probabilities pk into the objective, objective function becomes the expected total contribu-
tion to margin minus the plant configuration cost. The meaning of constraints do not change much but they
must be enforced for each scenario. The linear program given above is an example of a Stochastic Linear
Program. Note that going from a deterministic setting to a stochastic one both the number of constraints
and variables have increased significantly. The biggest issue in Stochastic Linear Programming is to solve
large Linear Programs using decomposition techniques ([1]). The plant flexibility problem presented here
is a simplification of the model given by [4].
12
demands. Thus, configuration decisions are nonanticipatory with respect to demand. Here we have se-
quence of events: Configure (nonanticipatory decision) - Find out demand (uncertainty resolution) - Ship
(anticipatory decision).
When Toyota decides on a recall, the recall decision is nonanticipatory with respect to the root cause
of the problem. The root cause can be deformation of pedal springs as opposed to software problems or
humidity build up that increases friction in the gas pedal. Once the uncertainty around the root cause is
resolved, Toyota has to decide on how to fix the root cause. The decision regarding the way the root cause
is fixed is anticipatory with respect to the root cause. Here the sequence of events: Recall (nonanticipatory
decision) - Find root cause (uncertainty resolution) - Fix root cause (anticipatory decision).
This distinction between anticipatory/nonanticipatory decisions is very important. In management, it
is natural that only a few decisions must be made simultaneously. Many decisions are staged over time.
The decisions that are made earlier (before uncertainty resolution) are nonanticipatory.
A word of caution on etymology: If you look at the dictionary, you will see that “anticipate” comes from
anti and capare. Anti means before and capare means to take in Latin. To anticipate is to take something
“before” it happens. Both “predict” and “forecast” have the same concept. Predict comes from prae and
dicere, so it is about saying something “before” it happens. Forecast comes from fore and cast, so it is about
puting forth something “before” it happens. All these three verbs, anticipate, forecast and predict, have
the same construction and they all include the word “before”. However, anticipatory variables relate to
decision that happen “after” uncertainty resolution. This terminology can be confusing if you take the
word “before” in anticipate too literally. Perhaps it is best to ignore the word “before” in anticipate.
6. Exercises
1. Aggregated computing capacity: The SOM currently has about 100 instructors, each of whom is
given a laptop computer to use in the classroom. These laptops have course slides and associated
software (access, excel, etc.) used during a lecture. A recent proposal is advocating for hosting these
slides and software on a central computer. According to this proposal, the instructors will use a key-
board and a computer screen in the classroom to connect to the central computer to project their slides
and to use the software. Each of the laptops currently has 1.8 GHz processing capacity and 75 GB
storage capacity. Suppose that if 100 laptops are connected appropriately (with an efficient parallel
processing architecture) the total processing capacity will become 180 GHz while the total storage
capacity becomes 7500 GB. First note that the capacity demanded by each instructor is random. Then
provide a verbal argument for why the central computer can serve the instructors as well as 100 lap-
tops even it has less capacity than 180 GHz and 7500 GB. Use the concept that extreme values cancel
each other in aggregation.
2. Deterministic capacity expansion with infinite horizon: Replicate Figure 3 by changing the follow-
ing parameters (a and b are independent):
a) f ( x ) = x0.2 .
b) r=20%.
Discuss how x ∗ changes in each case. Intuitively explain if you expected these changes, why?
√
3. Capacity expansion with finite horizon: While drawing Figure 3, we had f ( x ) = x, r = 0.05, and
δ = 1. We keep the same parameters but would like to compute the discounted expansion cost only
over T = 100 time periods. Note that we have considered an infinite horizon in Figure 3, but now we
are switching to consider a finite horizon of T = 100.
a) Let N ( x ) be the number of expansions made up to and including time T. For example, if x = 25,
we would make N ( x = 25) = 5 expansions over the horizon. These expansions are made at times 0,
13
25, 50, 75 and 100. If x = 30, we make N ( x = 30) = 4 expansions over the horizon, at times 0, 30, 60
and 90. In general,
⌊ ⌋
T
N (x) = + 1,
x/δ
where notation ⌊z⌋ indicates rounding down the number z. For example, ⌊2.2⌋ = ⌊2.7⌋ = 2. Explain
the formula for N ( x ) in English. Specifically, explain why we round down Tδ/x and why we add 1.
b) The cost C T ( x ) over the finite horizon of T becomes
N ( x )−1
CT (x) = f (x) ∑ exp(−rkx/δ).
k =0
4. Location with exchange rate uncertainty: Consider the facility location formulation with fixed in-
frastructure cost f i for plant i.
Minimize ∑14 14 18
i =1 f i ai + ∑i =1 ∑ j=1 cij xij
s.t.
∑14
i =1 xij = d j
j=1 xij ≤ Ci ai
∑18
This formulation does not have the variable and parameter names that you are used to. You can
guess that ai denotes the indicator variable for opening up plant i.
a) What does Ci denote, how many plant locations and markets are there in this formulation?
b) Suppose that all of our plants and markets are in the same country except for plant 1. Since Plant
1 is in a foreign country, the production cost of producing there is random in terms of the currency of
our home country. This production cost randomness makes the production plus transportation costs
of the form c1j random. Despite this randomness, we need to decide on plants to operate now, next
we will see the exchange rates and then we will determine the production/transportation quantity
xij . Is the production/transportation decision anticipatory with respect to exchange rate random-
ness? Is plant location decision anticipatory with respect to exchange rate randomness?
c) Modify the formulation above to incorporate two exchange rate scenarios: low cost c1j l and high
h that happen with equal probabilities. Alter anticipatory and nonanticipatory variables appro-
cost c1j
priately.
14
5. In many industries, such as Semiconductor and Automobile, subcontracting is becoming a popular
way of capacity management. A strong motivation for subcontracting is usually the cost advantage.
a) Explain, how can subcontractors provide products at a lower cost than in-house manufacturing
costs?
b) Given this cost advantage, what might be the disadvantages of subcontracting?
c) Comment on whether the disadvantages in (b) are more pronounced in Semiconductor or Auto-
mobile industry?
z∗ ( advance) ≥ z∗ .
Hint: Start with the formulations, one for each scenario, eventually yielding the average objective
value z∗ ( advance), combine them into a single formulation. Note that the only difference between
the combined formulation with objective value z∗ ( advance) and the original Stochastic Linear Program-
ming formulation with objective value z∗ is a constraint (which?). Finally, use the fact that adding a
constraint to a maximization Linear Program cannot increase its optimal objective value.
8. Read the Stochastic Capacity Planning section. Consider two demand scenarios d1j and d2j given for
all products j such that d1j > d2j for all products. Suppose that we set D j = d1j and D j = d2j , and solve
the initial LP to obtain objective values z1∗ and z2∗ . How does z1∗ compare against z2∗ , why?
15
9. Recall the 2 product and 3 plant example in the Stochastic Capacity Planning section. We want to
measure the net benefit of reconfiguring plant 2 to produce B. Suggest a step by step framework
using the Stochastic Linear Program given in the section to measure the net benefit.
10. Read the Stochastic Capacity Planning section. Suppose that there is a current plant to product as-
signment given by y0ij , i.e. Plant i has the tooling to produce Product j. However, we do not know if
additional Plant configurations can increase the profit.
a) Let yij (y0ij ) be the assignment after additional configurations (if any) made starting from y0ij , modify
the Stochastic Linear Program to find optimal yij∗ (y0ij ).
b) Let z∗ (y0ij ) be the optimal solution if the initial assignment is y0ij . Mathematically argue that
z∗ (y0ij ) ≥ z∗ (ȳ0ij ) if y0ij and ȳ0ij are two initial configurations such that y0ij ≥ ȳ0ij for all plants i and
products j. Express this result in English as well.
11. Refer to the solved exercise dealing with PlaToy. Now suppose that different dolls can be produced at
different factories but each doll is produced at only one factory. For example, Betsy and Vicky can be
produced at Plano while Wendy is produced at Richardson. But it is not possible to produce any one
of the dolls at two locations. For example, Betsy cannot be produced both at Plano and at Richardson.
In addition, there is an additional $10000 administrative cost incurred at a factory if any one of the
new dolls are produced at the factory: If 2, 1, 0 factories are used the cost is $20,000, $10,000, $0.
Formulate an MILP.
12. When the number of students suddenly increase at the School of Management in the Spring term,
what instruction capacity management strategy can be used? Chase strategy or level strategy, or
should the instruction capacity be maintained significantly over the regular teaching demand levels?
Answer this question twice, once from the perspective of the dean of the school, and once from the
perspective of UT system administrators.
13. Download flexman question.xls, which includes the data pertaining to question 8.4 on p.236 of the
textbook. Read the question statement.
a) The question does not mention any raw material costs for the routers or switches. What assump-
tions do you need to make to conclude that these costs are sunk (irrelevant): i) No backlog, ii) No
subcontracting, iii) No layoffs, iv) No new hires, v) Constant raw material costs over the year.
b) Compute the monthly regular time and overtime capacity in the spreadsheet.
c) By respecting the capacities computed in b), find a production plan that minimizes the total rel-
evant cost while ending december with exactly the same level of inventories at the beginning of
January. Report the production quantities, provide a printout of your final spreadsheet with the
objectives, decision variables and constraints clearly written.
14. Consider the product-to-plant assignment of major automobile manufacturers in North America in
scaggregate.ppt.
a) Which company owns the closest car manufacturing plant to UTD and what does that plant pro-
duce?
b) On the US map, consider a rectangle whose corners are the states of North Dakota, Oklahoma,
Arizona, Montana (NOAM rectangle). How many car manufacturing plants are there in the NOAM
rectangle? (Hint: What can NOAM stand for other than the first letters of the states?) In 1-2 sentences,
explain what could be the underlying reason(s) behind the automobile manufacturing activity in the
NOAM rectangle?
c) Consider the Toyota Plants in the US from west to east: 1: Nummi, 2: Long Beach, 3: San Antonio,
4: Blue Springs, 5: Princeton, 6: LaFayette, 7: Georgetown. Also consider the products of these plants:
16
A: Corolla, B: Tacoma, C: Hino, D: Tundra, E: Highlander, F: Avalon, G: Camry, H: Solara, I: Sequoia,
J: Sienna. List the product to plant assignment variables yi,j for i ∈ {1, 2, . . . , 7} and j ∈ { A, B, . . . , J }
that currently have value of 1.
d) Count the number of brands/types of cars/trucks GM is producing in the US and compare that
number against 10 for Toyota. Although Toyota has larger sales, it seems to be producing fewer num-
ber of types of cars/trucks than GM does. Would you expect a decrease in the number products in
GM’s product portfolio in an economic downturn, explain?
8. Appendix: MRP/OPT/JIT
8.1 MRP
MRP stands for material requirements planning. The name is reminiscent of the limited capabilities old
MRPs which had little functionality beyond figuring out a schedule of raw material requirements. How-
ever, modern day MRP systems come along with various modules associated with many business functions
from scheduling to accounting and to inventory control.
17
The first step of MRP is to look at the MPS and decide which, how many and when components are
needed to meet MPS. This process is called MRP explosion. During the explosion MRP uses a network
representation of product assembly called Bill of Materials (BOM). Suppose TexBag is company in manu-
facturing bags, then its BOM may look like Figure 5. According to the BOM a bag is made of two main
parts: a body and a strap. A body is manufactured by sewing 12 zips on a leather body. A strap is made by
attaching two hooks to the ends of a leather strap. Next to each operation, we also have its manufacturing
lead time. Lead times should not be confused with actual processing times. Lead times include transporta-
tion and wait-in-queue times as well as processing times. For example, assembly of a strap to a body can
be done in 5-6 minutes but the lead time can be a day mainly because of waiting. In summary, BOM is a
diagram for representing parts used in assembly and manufacturing operations, and lead times of these
operations.
Bag, 1 day
Leather Leather
body (1) Zips (12) strap (1) Hooks (2)
Now suppose that we need to deliver 50 bags on friday morning. Then we must have 50 bodies and
50 straps ready by thursday morning, we obtain thursday morning by backing up by 1 day assembly lead
time. To have 50 bodies on thursday morning, we must have 50 leather bodies and 600 zips on tuesday
morning. To have 50 straps on thursday morning, we must have 50 straps and 100 hooks on wednesday
morning. One can continue by backing up in time by operation lead times while computing how many
parts are needed for a certain number of assemblies. This process is called MRP explosion.
Once MRP explosion reaches the raw material level (leaves of the BOM tree), a schedule of how many
parts are needed and at what dates becomes available. This schedule of raw material requirements is
passed to Procurement department which orders raw materials form upstream companies in the supply
chain. Raw material acquisition schedule becomes an MPS for those upstream companies. MRP explosion
information helps to plan operations so that correct amount and type of inputs are ready for operations
and enough time is allocated to operations. For discrete product manufacturers whose BOM has multiple
levels, MRP explosion must be computerized. Computerization not only expedites the process but also
disciplines it and avoids manual errors. Enterprise resource planning (ERP) applies MRP logic to SCs, over
various production facilities of a firm or several firms.
MRP is the first conceptual representation that realized the underlying network structure in complex
production systems. It provides a systematic and coordinated way to schedule a large number of items over
this network. MRP is basically a specialized database system, so all advantages of database systems apply
to MRP as well. MRP streamlines information; it quickly makes correct information available. Streamlining
can foster standardization and further product development. For example, BOM is relevant for product
developers because it usually is a starting point for design improvements. Many design improvements
18
come as elimination or standardization of various parts. Suppose that TexBag is also manufacturing suit-
cases and uses zips in the suitcases. If the design team can standardize zips (say their length) then the same
zip can be used in bags and suitcases. Such a standardization decreases the complexity of processes; one
fewer product to name, buy, store, transport.
A major drawback of MRP is its failure to incorporate capacity restrictions appropriately during the
planning process. MRP first makes the plan with explosion and then checks for if there is enough produc-
tion capacity. If the capacity is enough, the plan is accepted. Otherwise, it has to be tweaked. MRP does
not have the capability to allocate the scarce capacity to products.
Another drawback is self fulfilling prophecy of lead times in MRP. Determining assembly and produc-
tion lead times is a tricky issue. MRP primarily aims to keep the production going on while efficiency is a
taken as a secondary objective. Thus, it puts in ample slack time into lead times even for short operations.
Later on, these long lead times become standard and short operations really take longer than they should.
This can be explained with a industrial psychological point of view: Those expectations that are set low
actually decreases the performance.
As a summary, we conclude that MRP is not the ultimate solution for production planning. It must
be supported with some decision making mechanisms. Also MRP data and performance levels should be
monitored closely to avoid self fulfilling prophecies.
8.2 OPT
Optimized Production Technology (OPT) attempts to overcome MRP’s insensitivity to capacity. OPT
divides resources into two categories: bottleneck resources and nonbottleneck resources. Bottleneck re-
sources are those resources that limit the production volume severely. For example in a serial production
system, the machine with the smallest capacity is the bottleneck resource. Once bottleneck resources are
identified, aggregate planning is geared entirely for the bottlenecks. In other words, it is assumed there
is enough capacity at nonbottleneck resources and they are assigned infinite capacity. Naturally, this ap-
proach simplifies the planning process by focusing the attention on a few of the resources. The exact steps
of how one obtains an aggregate plan with OPT is propriety information. However, it is general under-
standing that OPT software uses heavily tested network based heuristics for planning.
The main disadvantage of OPT is the concept of shifting bottlenecks. When the production volume and
the mix of products are known, we can find the bottlenecks in a system. However, the aggregate planning
exercise is done at least over several months and the volume or the mix may change from one week to
another. Different volumes or mixes can lead to different bottlenecks. When that happens, we say that the
bottleneck is shifting, i.e., it differs from one week to another. It is not clear how OPT handles this situation
because it relies on a clearly identified stationary bottleneck.
OPT focuses on bottleneck machines and ignores nonbotllenecks during planning. Thus, OPT provides
a plan for a production system that approximates the actual production system. In order for an OPT plan to
work, it is necessary to have plenty of nonbottleneck resources. When the cost of nonbottleneck resources
are not small, the OPT plan may have a high cost. This restricts the usefulness of OPT for cases when
nonbottleneck resources are expensive.
Both MRP and OPT are centralized decision making tools but in SCs many parties take part. The
rigid logic of MRP and OPT cannot accommodate multiple decision making parties, each serving its own
objective. Even if all players agree on a single objective, centralized decision making leads to overwhelming
bureacracy: inputs and data must be transmitted up to a single decision maker and each decision must be
transmitted down to interested parties. These transmissions increase the likelihood of error occurrences.
Moreover, since the parties and the decision maker are substantially apart, the decision maker cannot detect
errors in the data. Similarly, the parties cannot detect the flaws in the decisions. Transmissions also delay
the decision making process, that is decision can be made much faster locally.
19
8.3 JIT
JIT (just-in-time) differs from MRP and OPT in its philosophy of decentralization. Decisions and improve-
ments are made locally. For example, JIT proponents argue that workers know enough about their job to
make correct decisions about it. Industrial psychologist also support JIT by saying that worker empower-
ment increases the worker performance. Naturally, passing decisions to workers, decreases the coordina-
tion among different group of workers.
The coordination among different units is achieved with Kanban system. Kanban is basically a work
order that generally travels upstream in the SC. It indicates what is to be done (produced) in what amount
(number of units). The Kanban’s are driven from MPS. The important distinction between JIT and MRP
or OPT is that no production takes place at a workcenter before a Kanban card reaches that work center.
Therefore, JIT works in the pull mode whereas MRP and OPT are more of push types.
Some of the JIT aspects are given below:
• Set up time/cost reduction by taking set up time out of the production cycle, modifying jigs and
fixtures, etc.
• Quality improvement. Concepts such as total quality management has come from these efforts.
• Improving supplier relations by reducing the number of suppliers and strengthening the communi-
cation and coordination with the remaining suppliers.
• Process improvement is encouraged. Especially factory workers are encouraged to report their sug-
gestions and prize money is given to good suggestions.
JIT philosophy is quite different from production planning and inventory control in that it attempts to
remove problems rather than planning around them. Quality specialists term this as a root cause approach.
They would like to find the real problem that causes the symptoms. For example, unreliable suppliers can
be a result of wrong choice of suppliers or the lack of coordination with the suppliers. In this case, JIT will
attempt to increase supplier reliability whereas an inventory control based approach would study the best
plan for the given levels of unreliability. Perhaps the success of JIT is due to its nature of modifying the
rules of the game rather than playing with these rules.
Discussion Questions:
1. Put the following industries in an order of the highest benefactors of MRP application to the lowest:
Pharmaceutical, Catering, Automobile, Semiconductor, Health. Briefly explain how you came up
with the order.
2. Remember the last time you have been to a Taco Bell restaurant. Do you think Taco Bell uses a
modular BOM, explain why by giving examples.
3. Suppose that we are TexBag’s supplier of Leather Straps and Hooks (see Figure 5). Further suppose
that the leather strap and hooks are our only products and we need 4 labor hours and 1 labor hour to
manufacture 1 leather strap and 1 hook respectively.
a) Suppose that TexBag is our sole customer, aggregate our two products into one generic product
so that we can handle labor force planning problem as a single product problem. How many labor
hours does the generic product require?
b) Suppose that we have another customer than TexBag called ArkBag which buys only leather straps.
Can we still stick with our answer to (a)? If not, do we have enough information to modify our answer
to (a)? If not, what additional information is needed?
20
References
[1] J.R. Birge and F. Louveaux (1997). Introduction to Stochastic Programming. Springer-Verlag, New York.
[2] S.E. Bodily and J.R. Freeland (1988). A simulation of techniques for forecasting shipments using firm
orders-to-date. Journal of Operational Research Society Vol.39 No.9: 833-846.
[3] V.M. Guerrero and J.A. Elizondo (1997). Forecasting a cumulative variable using its partially accumu-
lated data. Management Science Vol43. No.6: 879-889.
[4] W.C. Jordan and S.C. Graves (1991). Principles on the benefits of manufacturing process flexibility.
Technical report GMR-7310, General Motors Research Laboratories, Warren, MI.
[5] A.S. Manne (1961). Capacity expansion and probabilistic growth. Econometrica 29: 632-649.
21