IE 503 - Operations Analysis: Narayan - Rangaraj@iitb - Ac.in
IE 503 - Operations Analysis: Narayan - Rangaraj@iitb - Ac.in
IE 503 - Operations Analysis: Narayan - Rangaraj@iitb - Ac.in
Any cheating will result in grade penalty and reporting to the Dean AP’s office
No tutorials. Assigned reading, some assignments. TAs and I will be available at some
times for discussion on some topics, to be announced on Friday
TAs will do some grading of questions, but will not answer cribs about evaluation –
that is done only by me
Background material and linkages with
other courses (IEOR students)
Modeling exercises in lab courses (IE 507)
– Especially formulation of decision models
Dynamic Programming
– Shortest Path Models
B.Tech/BS/B.Des - DD students
You may have done works visits, internships,
projects etc. in some industry
No pre-requisites
1) In steady state, how many hostel rooms does the institute require,
assume all durations are deterministic?
2) Suppose c) students take 5-7 years and d) students take 4-6 years, how
will you plan capacity? You will have to assume something about the
uncertain duration.
3) For a)-d), past data is known regarding proportion of acceptance and
also male female split (which is different year to year). How do you plan
hostel accommodation for males and females separately, if desired?
Assignment 0 – problem 2
On a (linear) stretch of beach, a vendor intends to set up a stall to sell food. Customers
are uniformly spread over the entire length L of the beach. Propose different criterion
on the basis of which the vendor could decide the location of the stall.
a) If the vendor wants to minimize total or average distance travelled by customers
where would she locate?
b) Customers accesses stall with probability p inversely proportional to distance to
stall. If we want to set up stall to maximize the expected number of customers that we
will get, where should we put it up?
c) If we can set up two stalls, where should we put them up?
d) If the two vendors are competing, how would they position themselves? Argue that
if one vendor is at location a in the interval [0,L], what would be the best location for
the other vendor, and thereby find the best location chosen by both vendors.
e) If the two vendors co-operate, then decide where they should set up.
X = (x, y)
• d(X, Pi) = (x - a i ) 2 + (y - bi )2
X = (x, y)
• Squared Euclidean
• Chebyshev distance (L∞ norm) Pi = (ai, bi)
• Minkowski distance
'
• d(X, Pi)= ∑ 𝑥 − 𝑎 % & (
Distance measures - continued
• Rectilinear/Manhattan/Mannheim/Chandigarh
• Facility allocation in Grid like cities and warehouses (Look at maps of a few cities)
• Euclidian distances (as the crow flies)
• Squared euclidean also used
• Building of cell-phone towers
• Chebyshev distance
• Used to calculate the time taken for movement of overhead crane (crane can move on
the x and y axes at the same time but at the same speed along each axis)
• Network distance
Mannheim, in
Baden-Württemberg
Germany
From
Openstreetmap.org
Chandigarh
From
Openstreetmap.org
Example
• 5 demand locations on the plane, e.g. machines on a shop-floor
• P1 = (1,1), P2 = (6,2), P3 = (2,8), P4 = (3,6) and P5 = (8,4)
• New general purpose machine to be located on the shop-floor
• Number of demands (trips per day) from new machine to existing
machines
• 10, 20, 25, 20 and 25
• Cost of movement proportional to distance
• Where should new machine be located?
• In general, given locations (𝑎𝑖 , 𝑏𝑖 ) and weights 𝑤𝑖, where should we
locate the facility to minimize sum of weighted distance to all 𝑖’s?
Solution concept for different distance measures
If distances are rectilinear distances, then solution by linear programming (LP)
or other algos
Basic idea: Min |𝑥| is equivalent to Min 𝑧, s.t. 𝑧 ≥ 𝑥, 𝑧 ≥ −𝑥
We can carefully replace all absolute value terms with such a construction
Simplified algorithms available which solve this LP
BU SIN ESS
From dw.com Why are German chemical plants located near big cities?
(Deutsche Welle) à A fatal explosion at a chemical complex near Cologne has raised concerns about the park's proximity to one of Germany's
biggest cities. Chemical plants located near major cities are a common sight in the country.
The chemical industry in the Rhineland — a region of the German state of North Rhine-Westphalia — which includes the chemical park in
Leverkusen, is one of the largest hubs for the industry in Europe. The region has over 260 chemical companies employing more than 70,000
people. The firms are often clustered together in so-called chemical parks.
https://www.dw.com/en/why-are-german-chemical-plants-located-near-big-cities/a-58685101
From Wesseling in the south to Dormagen in the north: Global companies such as Bayer, ExxonMobil Chemical, Ineos, Covestro,
Lyondellbasell and Lanxess are located within a few kilometers of Cologne, a city with a population of over 1 million people.
To understand how a major city and a chemical park can exist in such close proximity, one has to look back at the long history of the Rhineland
as an industrial powerhouse.
3
(4) Objective, to minimize sum of
weighted distance to all nodes
A 5
B
(5) (2) 4
Example: If on arc A-B, x units from A,
6 Weighted distance =
D (5)*(x) + (2)*(5-x) +
(3) (4)*(8-x) + (3)*(11-x)
Thegooglemaplink(h
ttps://goo.gl/maps/UZ6cSeb1F6mpnMoe8)showsthemapofthe
centreofMannheimtowninGermany.Ithasagrid-likestructure(thisareaiscalledthe
Quadrate).
Therearemanyrestaurantsofvariouscuisinesinthetownandoffersomeparticularkindsof
service(dine-in,takeawayandhomedelivery)tothecustomers.
Afewoftheserestaurantsareconsideredforfurtheranalysisandthetablebelowdescribes
theirnameandexpectedno.ofweeklyordersbelongingtothecategoryofhomedeliveryin
thisarea.
Restaurantname Weeklyordersatlunchtime
Lenok’sBurger 22
RoyalDonuts 13
LittleIndia 9
Subway 25
SurinKuche 7
PontePizzeria 18
Supposethatthecustomersareuniformlydistributedovertimeandasingledelivery
operatorinthetownproposestoofferacommonservicetofulfilthedemand.Constructan
appropriategridstructureandlocatetheaboverestaurantstocarryoutthefirstcutanalysis.
Usingthedataandinformationgivenabove,answerthefollowing:
a)Iftheoperatorwantstominimizetotaloraveragedistancetravelled,wherewouldhe
locatehimself?
Trythegriddistancealongtheroadsandtheresultingsolutionandalsothecentreofgravity
solutionandcheckhowdifferenttheyare.
b)Iftherearetwooperatorslocatedsuchthatthedemandofarestaurantisfulfilledbythe
operatornearesttoit,thencommentwheretheyshouldlocatethemselves.Notethatweare
nottalkingaboutdynamicallyallocatingorderstotheoperators,butwedonotdecideapriori
thatsoandsorestaurantwillbeservedbysoandsooperator.Aftertheone-time
calculation,theallocationisfixed.
More on location and
location-allocation problems
Location on networks
• p-median problem, to locate p facilities at nodes in a network, with the
objective of minimizing weighted distance to n given demand or service
locations (note that demand is at discrete locations on a network)
• For a large network can be tough, nCp possible combinations of locations
• Each combination will result in an allocation of demand nodes and then the
objective needs to be computed
• Generally speaking, no efficient way available to do a much simpler
computation (NP hard)
• Many good heuristics available
Location along arcs
• What if we can locate facilities along arcs (if we have the choice)?
• Check for simple networks (one arc, then two, then a triangle) that if
a location inside an arc is optimal, then one of the (neighbouring?)
nodes is also optimal
• Check on all (weighted) demands met through one end point of arc
and through the other end point
• Summary: It is enough to look at solutions at nodes
• Solutions on arcs may still be optimal
• For practical reasons, location at some nodes may not be possible
Try the following
For this network, also try a 1 median solution which is as far away as possible from any of the nodes, i.e. if x is the
chosen location, then min(d(x,A), d(x,B), d(x,C)) is as large as possible. This is the obnoxious facility location problem. Is
this also guaranteed to have a solution at a node? What about min(max d(x,A), d(x,B), d(x,C)) - the center problem?
Median problems on networks – Hakimi’s
theorem
• A network is given and there are demands at nodes to be satisfied
• Demand at node j is hj
• If facilities are set up at points Xk = x1, x2, …, xk anywhere on the network, then every
demand point j will access the facility that is closest to it, call this distance as d(j,Xk),
i.e. distance from j to set Xk
• Distance to a point on an arc is measured through the end points of the arc (i.e. all movements
are on the network of arcs)
• The weighted accessibility function is then ∑ ℎ# 𝑑(𝑗, 𝑋) )
• The k median problem on the network is to find the set Xk which minimizes the
weighted accessibility function
• Hakimi’s theorem says that a solution to the k median problem on a network can be
found among the nodes of the network
• Verify from earlier problem that the 1-median problem always has a solution at a node
• Our beach problem had continuous spatial demand (there was no concept of a node)
1 and 2 median problems on networks -
example
https://web.mit.edu/urban_or_book/www/book/chapter6/6.5.2.html
e.g.
In the 1-median case, cost if facility is located
at D is 43
Packaging
Milling
Assembly
Drilling
Packaging
Milling
Assembly
Drilling
Example 1 Example 2
T D T *
T P P M A P
P T A * D M
A M One possible
Available Another possible
D A area allocation allocation
M D For each allocation of departments, the associated
item/people cost can be computed using the movement
Available One possible Another possible information
area allocation allocation
Appropriate distance/cost measure needed
For specific example, pl see note Could be multiple arrangements where the same cost is
accompanying lecture slides incurred, based on relative positioning of departments
Linear and U-shaped layouts
• If we have n facilities to locate in a product oriented layout, the
relative location is decided by the product flow (i.e. the sequence of
operations will determine the relative location)
• So the default layout is linear, in order of the operation sequence
• If the available area has a particular shape, can lay out the facilities in
a U shaped or other convenient shape, respecting the available area
• U shaped cells are preferred in manual systems because there is
better interaction between workforce members and movement
between stations (in case of re-work and any trouble-shooting) is
more convenient
The following simple layout problem has to be modelled.
A four-way intersection has 2 bus routes, number 61 and number 62, which cross each other at
the intersection (61 goes east <—> west and 62 goes north <—> south at the intersection).
Assume that buses drive on the left of the road (as in India). There is space to create bus stops
either before or after the junction, so the 61 E-W bus can stop at location J before the junction
or M after the junction. We would like to locate the bus stops so as to maximise the benefit to
passengers who wish to change buses (for example, as few people as possible should cross the
road for bus change).
North
G H
F I
West East
M J
L K
\ South
The buses are every ten minutes and staggered so that transfers both ways can take place.
Assume we have data about how many passengers generally change from 61 E-W to 62 N-S and
all other combinations. Here is the estimate of number of people per bus who wish to change
(the traffic is symmetric, when viewed for the entire day, ie. N-S people who wish to change to
E-W will be the same as W-E changing to S-N, etc.):
N-S S-N E-W W-E
N-S -- -- 8 12
S-N -- -- 5 10
E-W 10 12 -- --
W-E 5 8 -- --
For example, 5 people per S-N bus wish to transfer to an E-W bus on average.
Propose a model to answer this question in general.
Layout problem : Linear assignment and Quadratic assignment
Data: The owner of a job shop has to locate four machines: (1) a punch press, (2) a
grinder, (3) a lathe, and (4) and a welding machine. Four locations are available for
this purpose - labeled (A), (B), (C) and (D).
The foreman has estimated the costs of material movement to and from the new
machines from the existing facilities in Rs per hour as follows
Questions:
b. Write out a mathematical expression for the optimization problem that achieves
this.
See if you can solve this. How will the formulation change if you have one
more location E available as a possibility?
c. From an enumeration point of view, how many possible ways are there to assign
facilities to locations in this sort of problem?
d. How do you think the costs of locating a facility at a particular site (A, B, etc.
are computed)? Here the foreman has given you this data, but you should know
how to get this.
More difficult problem:
A variant is the following. In addition to the movement between (1), (2), (3) and
(4) and existing facilities, there is movement between (1), … (4) among the new
machines also (e.g. between (1) and (3) and other pairs). This also has to be
included in the decision.
Now suppose costs of item movement between the new facilities (1), (2), (3) and
(4) are as follows (in Rs per unit distance per hour – this factors in the
requirements per hour of jobs moving to and fro these facilities):
Given this and the earlier data, formulate an appropriate cost function and
optimization problem and then determine the optimum location of the new
facilities (1), (2), (3) and (4).
Sample problems for IE 503 : set 1
Answer: Median location : co-ordinates (10,6) - keep adjusting one of the co-ordinates at a time
till half the weight is exceeded for the first time.
Notes: Check that the problem is separately solvable in x and y co-ordinates. Check for
alternate solutions. Check that for some range of values of the extreme co-ordinates of the
existing departments, the median solution would remain the same.
Note: The constructive approach to the 1-median problem works, but is not really extendable to
p medians.
Answer: Node B.
The obnoxious facility location problem is NOT guaranteed to have a solution at nodes.
Typically such problems are not solved as weighted objectives (all demands are equally
important, but weighted versions are also possible). Single arc network with d = 2 and equal
weights on demands will give the unique solution on the mid-point of the arc.
6. Modeling for location on networks: Lecture 4, slide 6 (these models are not unique
and could vary depend on the notation used). Please read if interested.
Use notation y_i = 1 if facility at i and y_i =0 if not, and x_ij = 1 if demand at j is (fully) supplied
from i and x_ij = 0 if not.
Can also try with 0 <= x_ij <= 1 representing the fraction of j’s demand met from i.
Answer :
∑ 𝑦𝑖 ≤ 𝑝
𝑖
● Fixed cost of facilities that are selected,
∑ 𝑥𝑖𝑗 ≥ 1 ∀ j
𝑖
∑ 𝑑𝑗𝑥𝑖𝑗 ≤ 𝐾𝑖 ∀ i
𝑗
∑ 𝑥𝑖𝑗 ≤ 𝑘 ∀ i
𝑗
Layout Decisions as
Location/Assignment Problems
Example
Owner of a shop has to locate 4 machines:
• (1) A Punch Machine
• (2) A Grinder
• (3) A Lathe
• (4) A welding Machine
• 4 locations are available for this purpose labelled as : A, B, C, D
• Locate machines so that the total cost of material transfer to and from
the new machines is as small as possible
The Cost matrix for movement is
Machine | Location A B C D
• Let us try assigning low cost assignments first: gives (4) at A, then
(1) at B, then (3) and C and (2) at D
𝑠𝑢𝑏𝑗𝑒𝑐𝑡 𝑡𝑜 𝑐𝑜𝑛𝑠𝑡𝑟𝑎𝑖𝑛𝑡𝑠
∑C∈N 𝑥CD = 1 ∀ j ∈ L
∑D∈M 𝑥CD = 1 ∀ i ∈ 𝑀
𝑥CD ∈ {0,1}
Relaxing the constraints 𝑥𝑖𝑗 ∈ {0,1} to 0 ≤ 𝑥𝑖𝑗 ≤ 1 and solving
as a Linear Programme (which can be done efficiently) works
in this case
Here the foreman has given you this data, but you should
know how to get this
More general layout problems
• Suppose the distance between two locations is given
𝑠𝑢𝑏𝑗𝑒𝑐𝑡 𝑡𝑜 𝑐𝑜𝑛𝑠𝑡𝑟𝑎𝑖𝑛𝑡𝑠
∑C∈M 𝑥CY = 1 ∀ 𝑘 ∈ 𝑀
∑Y∈N 𝑥CY = 1 ∀ i ∈ 𝐿
𝑥CY ∈ 0,1
Product Variety Large variety High variety Low variety Very low or no variety
Flow of Product Unique for every Unique for every batch Sequenced Sequenced
product
Type of Customer One at a time Few individual customers Mass market Mass market
Worker skills required Experts, Craftsmen Wide range of skills Limited range of skills Equipment Monitors
Production System Long term project Intermittent job shops Flow lines, Assembly Process Industry
lines
Advantages Custom work, Latest Flexibility, Quality Efficiency, Speed, Low cost Highly efficient, Large
Technology Capacity, Ease of Control
Examples Construction, ship Machine Shops, Print Automobiles, Paint, Chemicals, Food
building, Aircraft Shops, Bakeries, Televisions, Computers, stuff
Manufacturing Education Fast food
Main issue in product oriented layout
● Relative locations of facilities already determined (exact orientation can be
decided – often U shaped)
● Most often used for high volume production of standardized items
● Main issue is to balance the throughput of each stage
● Idea of bottleneck resource
● Smooth material flow through every stage, no accumulation between stages
● With multiple operations at each stage/machine/work station, for example,
assembly lines, then we need to do line balancing
Main issue in process layout
● Many different technologies/departments/facilities
● Relative locations to be determined
● Multiple products each with a routing through the factory
● Depending on volumes of each product, we may decide on relative locations
to minimize total cost of item movement
● With several possible locations, we may think of it as an Assignment problem,
i.e. the decision is “what goes where”
○ If relevant costs can be computed once we decide on the location of one facility, then it is a
linear assignment problem – easy to solve
○ If we need to compute costs only after deciding on pairs of locations, it is a quadratic
assignment problem – much more difficult to solve, for large cases
Automobile - robotic assembly line
Source: indiamart.com
Automobile - manual assembly line
Source: toyotavn.com
A scooter production line
Boeing assembly line
Boeing introduced a faster layer of automation into its 737 factory, adding a new
robotic system to drill holes in the main beams inside each wing
Source: assemblymag.com
Planning of Assembly Lines
Instance of Product-
oriented Layouts
Flows and bottlenecks
Machine A Machine B Machine C Machine D
Bottleneck
facility
Flow
Shifting bottlenecks
• Bottleneck is a process that limits (for now) the output of a system
• Could be machines, departments or other facilities
• Because facility addition is sometimes chunky and because everyone likes to have more
capacity to work with
Shifting bottlenecks
24 items an 40 items an 12 items an 30 items an
hour hour hour hour
Bottle neck
Flow
12 items an
24 items an 40 items an hour 30 items an
hour hour 20 items an hour
hour
Bottle neck
Flow
Assembly lines – an important example of
product oriented layouts
• Many operations in assembly
• Subassemblies
• Components
• Manual/semi-automated/automated/robotic
• Sequence of operations is important
• Example 1: Some fastening of internal components has to be done before
external cover attached
• Example 2: All operations on one side to be done before turning around
• Want to have unidirectional flow of the assembled product
• Save time
• Avoid confusion (what operations are done and what not done)
Work station
• A physical location where one or more operations are performed
• Fixtures have to be available to do all these operations
• All components/sub-assemblies have to be stored and made available
• If manual, then workers are trained and able to perform all the
operations
• Even within a workstation, there could be a sequence of operations
• All predecessor operations should have been done in prior
workstations – so that there is no back-flow of items
• Total work content of a work station is the sum of all operation times
Planning of assembly lines
• Operations have:
• Processing time
• Predecessor activities which need to be completed before
that operation
• Exercise: The same procedure can be used to find the smallest cycle time for a
given number of workstations, by successively trying different cycle times.
Example from Nahmias – 12 operations, total work content = 70 mins
Operation Immediate Processing time Ranked positional
predecessor weight
1 - 12 70
2 1 6 58
3 2 6 31
4 2 2 27
5 2 2 20
6 2 12 29
7 3,4 7 25
8 7 5 18
9 5 1 18
10 9,6 4 17
11 8,10 6 13
12 11 7 7
In the example, suppose the desired number of
products/day manufactured is 32
• The quality of the solution is measured by the total idle time for
all the workstations in every cycle.
• For the given operation times, the smallest theoretical number
of workstations may or may not be achievable
Actual operation
• The smoothest operation of such an assembly line is a
transfer line, with products moving from one station
to the other every cycle time minutes.
• In steady state, there is no pile up of inventory
between workstations and material is ideally pulled
into the system at the first station when needed
• Within a workstation, operations can be done by the
operator in any sequence that respects
the precedence constraints
Example of Assembly Line
Balancing (from Nahmias)
Operation list for (old) computer assembly
1. Drill holes in metal casing and mount brackets to 7. Mount hard disk drive. Attach hard disk controller
hold disk drives (12) and power supply to drives (7)
2. Attach motherboard to casing (6) 8. Set switch settings on the motherboard for specific
configuration of system (5)
3. Mount power supply and attach to motherboard (6) 9. Attach monitor to graphics board prior to running
system diagnostics (1)
4. Place main processor and memory chips on 10. Run system diagnostics (4)
motherboard (2)
5. Plug in graphics card (2) 11. Seal the casing (6)
6. Mount disk drives and attach drive controller and 12. Attach company logo and pack system for shipping
power supply to drives (12) to customer (7)
For example, when operations 5,6 and 9 (in some order) are being carried
out on assembly C in the third workstation, operations 2,3 and 4 (in some
order) will be carried out on assembly B in the second workstation
Location Problem
D Subway 25 (4,1)
1a) Objective is to minimize the total distance (rectilinear distance in this case) that needs to
be travelled to fulfil all the demand. Therefore,
i) Rank the coordinates of supply locations in ascending order and accumulate the weights.
ii) The optimal value of the coordinate is found by dividing the total cumulative weight by 2
and identifying the first location at which the cumulative weight exceeds this value.
∑𝑤𝑖*𝑥𝑖 ∑𝑤𝑖*𝑦𝑖
ℎ = 𝑖
= 4. 94 , 𝑘 = 𝑖
= 3. 19
∑𝑤𝑖 ∑𝑤𝑖
𝑖 𝑖
i) Divide the demand location into clusters (equal to no. of operators, 2 in this case), where
each cluster is served by one of the operators.
ii) The best possible location for both the operators can be calculated separately using the
same approach used in part (a).
iii) Check the feasibility i.e., whether the restaurant is being served by the nearest operator
or not. Calculate the total distance.
iv) The process is repeated by making changes in the cluster formation until we find the
optimal solution i.e., the solution with least distance.
Optimal solution:
We aim to minimize the number of road crossings by the people who wish to change the bus
Objective Function Value : Total Number of Road crossings
So, no improvement in objective value in this step so our final Locations for bus
stops are at G,F,K,J with total number of road crossings = 68
Subjected to constraints:
d) Machines are involved in machining of various products. To calculate the cost for
installation of a particular machine at a location, we first see for how many products
this machine is used, then we see the if the machine is at a particular location what is
its distance from all other machines, then total cost of machine for one product type
is calculated by multiplying quantity with the sum of distance from which product is
coming and to which it is going after getting machined from the particular machine.
2nd part: where cost of material movement from one machine to another machine is also
considered.
- Mathematical model for this: Similar to above one, but some additional cost that to be
considered is as follows:
As the objective function has two variables in the same expression, therefore it becomes a
quadratic mathematical model which is hard to solve.
Approach to solve: heuristic based solution can be used to solve this problem.
Improvement based heuristics :
➔ Start with some layout (assignment of facilities)
➔ Check on neighbouring solutions (e.g. interchanging the location of two
facilities)
➔ Accept improvements and continue till no more improvement possible
Quantity Length
2 2
4 3
3 4
1 6
2 7
First Fit Decreasing Method : Fit the item of large size first
Bin 1 7,4
Bin 2 7, 4
Bin 3 6, 4
Bin 4 3, 3, 3, 3
Bin 5 2, 2
First Fit Increasing Method: Fit the item of small size first
Bin 1 2, 2, 3, 3
Bin 2 3, 3, 4
Bin 3 4, 4
Bin 4 6
Bin 5 7
Bin 6 7
Optimal Solution
Bin 1 7 , 3, 2
Bin 2 7, 3, 2
Bin 3 4, 4, 4
Bin 4 6, 3, 3