The Solution of Homework #2(Isye 6202)

The Solution of Homework #2(Isye 6202)

The solution of homework #2(ISyE 6202)

The Solution of Homework #2

#12.22

We have 6 product families and 40 storage bays that can be used for storing product families. For further reference, we number each location as follows:

1 / 2 / 3 / 4 / 5 / 6 / 7 / 8
9 / 10 / 11 / 12 / 13 / 14 / 15 / 16
17 / 18 / 19 / 20 / 21 / 22 / 23 / 24
25 / 26 / 27 / 28 / 29 / 30 / 31 / 32
33 / 34 / 35 / 36 / 37 / 38 / 39 / 40

Then the problem can be solved as follows:

Assumptions:

  • We consider single-command type of operation.
  • Distance is computed to/from the center of each location.
  • The activity associated with each location is distributed evenly between receiving and shipping (otherwise, we would have accumulation of material in the warehouse or “systematic” shortages).
  • Distance is measured according to the rectilinear metric.

Parameters:

  • I = the index set of sku’s: I = {1,2, …, 6}
  • J = the index set of locations: J = {1,2,…, 40}
  • THi = number of units of sku i handled per unit of time
  • Ni = number of storage locations allocated to sku i
  • dj = expected travel distance per unit load for location j

Decision variables:

  • Xij = 1 if location j is allocated to sku i; 0 otherwise for all iI and jJ.

From the problem statement, we have the following data with respect to each product family:

Product Family / Area / Load Rate(THi) / Ni / THi/ Ni / Rank of THi/ Ni
1 / 2,400 / 500 / 6 / 83.33 / 3
2 / 3,200 / 250 / 8 / 31.25 / 6
3 / 2,000 / 650 / 5 / 130.00 / 2
4 / 2,800 / 450 / 7 / 64.29 / 4
5 / 4,000 / 375 / 10 / 37.50 / 5
6 / 1,600 / 750 / 4 / 187.50 / 1

where, Ni =storage locations required by SKU I = Area to be occupied by SKU I / (20 ft * 20ft) and THi = throughput rate of SKU I.

Since , this problem corresponds to a balanced transportation problem and its formulation is:

Since the cost cij associated with variable xij has the product form of a term associated with SKU I (THi / Ni)and a term associated with location j (dj), We can get the optimal solution of the above problem using the following matching procedure:

  • Rank all the available storage locations in increasing distance of dj
  • Rank all sku’s in decreasing turns of THi/ Ni
  • Move down the two lists, assigning the next most highly ranked sku i to the next Ni locations

Now we can compute the expected distance for each location as follows:

j / 1 / 2 / 3 / 4 / 5 / 6 / 7 / 8 / 9 / 10 / 11 / 12 / 13 / 14 / 15 / 16 / 17 / 18 / 19 / 20
Lin,j / 40 / 60 / 80 / 100 / 120 / 140 / 160 / 180 / 20 / 40 / 60 / 80 / 100 / 120 / 140 / 160 / 20 / 40 / 60 / 80
Lj,out / 180 / 160 / 140 / 120 / 100 / 100 / 120 / 140 / 160 / 140 / 120 / 100 / 80 / 80 / 100 / 120 / 140 / 120 / 100 / 80
dj / 110 / 110 / 110 / 110 / 110 / 120 / 140 / 160 / 90 / 90 / 90 / 90 / 90 / 100 / 120 / 140 / 80 / 80 / 80 / 80
Rank / 4 / 4 / 4 / 4 / 4 / 5 / 7 / 8 / 2 / 2 / 2 / 2 / 2 / 3 / 5 / 7 / 1 / 1 / 1 / 1
j / 21 / 22 / 23 / 24 / 25 / 26 / 27 / 28 / 29 / 30 / 31 / 32 / 33 / 34 / 35 / 36 / 37 / 38 / 39 / 40
Lin,j / 100 / 120 / 140 / 160 / 40 / 60 / 80 / 100 / 120 / 140 / 160 / 180 / 60 / 80 / 100 / 120 / 140 / 160 / 180 / 200
Lj,out / 60 / 60 / 80 / 100 / 120 / 100 / 80 / 60 / 40 / 40 / 60 / 80 / 100 / 80 / 60 / 40 / 20 / 20 / 40 / 60
dj / 80 / 90 / 110 / 130 / 80 / 80 / 80 / 80 / 80 / 90 / 110 / 130 / 80 / 80 / 80 / 80 / 80 / 90 / 110 / 130
Rank / 1 / 2 / 4 / 6 / 1 / 1 / 1 / 1 / 1 / 2 / 4 / 6 / 1 / 1 / 1 / 1 / 1 / 2 / 4 / 6

where, Lin,j and Lj,out represent the distance from the location j to the receiving and shipping dock respectively and dj=0.5*(Lin,j + Lj,out). Using the actual layout of the 40 locations, we get the following ranking profile for them:

4 / 4 / 4 / 4 / 4 / 5 / 7 / 8
2 / 2 / 2 / 2 / 2 / 3 / 5 / 7
1 / 1 / 1 / 1 / 1 / 2 / 4 / 6
1 / 1 / 1 / 1 / 1 / 2 / 4 / 6
1 / 1 / 1 / 1 / 1 / 2 / 4 / 6

The rank of each sku has been already computed in the first table above. So, following the matching procedure outlined above, an optimal allocation is as follows:

5 / 5 / 5 / 5 / 5 / 2 / 2 / 2
4 / 4 / 4 / 4 / 4 / 5 / 2 / 2
6 / 6 / 3 / 1 / 1 / 4 / 5 / 2
6 / 3 / 3 / 1 / 1 / 4 / 5 / 2
6 / 3 / 3 / 1 / 1 / 5 / 5 / 2

Notice that there are several alternative optimal solutions. The one provided above seeks also to cluster all the locations assigned to each SKU in the same area of the warehouse, since this would facilitate the better monitoring and management of the warehouse activity.

# 12.24

We have 3 products (X, Y, and Z) and 44 storage bays that can be used for storing products. Again, we number locations as follows:

Assumptions:

  • We only consider single command case.
  • Distance is computed to/from the center of each location considering round trip.
  • In an aisle, a worker/vehicle moves along the center of the aisle.
  • Travel distances are computed as shortest paths through the depicted warehouse aisle structure.

Since , this problem corresponds to a balanced transportation problem and we can use the same definitions for parameters and variables and the same formulation as in the problem #12.22. We can compute the following:

Product / Area / Load Rate(THi) / Ni / THi/ Ni / Rank of THi/ Ni
X / 6,400 / 500 / 16 / 31.25 / 2
Y / 8,800 / 400 / 22 / 18.18 / 3
Z / 2,400 / 500 / 6 / 83.33 / 1

Now we can apply the same optimal solution procedure as in #12.22 and the computation of the expected distance for each location is as follows (Notice that in computing these distances, we have taken advantage of the symmetry of the warehouse layout with respect to the central horizontal axis).

j / j / Distance for stocking / Distance for shipping / Expected distance / Rank of dj
1 / 34 / (38+18)+(18+78+56+96)=304 / (248+38+14)*2=600 / 452 / 6
2 / 35 / (56+38+18)+(18+58+56+96)=340 / (156+56+58+18)*2=576 / 458 / 7
3 / 36 / (56+58+18)+(18+38+56+96)=340 / (156+56+38+18)*2=536 / 438 / 5
4 / 37 / (56+78+18)+(18+18+56+96)=340 / (156+56+18+18)*2=496 / 418 / 4
5 / 38 / (96+38+18)+(18+38+96)=304 / (156+38+18)*2=424 / 364 / 3
6 / 39 / (96+56+38+18)+(18+38+56+96)=416 / (156+56+38+18)+(18+114+56+4)=460 / 438 / 5
7 / 40 / (96+56+58+18)+(18+58+56+96)=456 / (156+56+58+18)+(18+94+56+4)=460 / 458 / 7
8 / 41 / (96+56+78+18)+(18+78+56+96)=496 / (156+56+78+18)+(18+74+56+4)=460 / 478 / 8
9 / 42 / (96+56+98+18)+(18+98+56+96)=536 / (156+56+98+18)+(18+54+56+4)=460 / 498 / 9
10 / 43 / (96+56+118+18)*2=576 / (156+56+118+18)+(18+34+56+4)=460 / 518 / 11
11 / 44 / (96+56+138+18)+(14+38+152+96)=608 / (156+56+138+18)+(14+38+4)=424 / 516 / 10
12 / 23 / (18+18)*2=72 / (234+18)*2=504 / 288 / 1
13 / 24 / (38+18)*2=112 / (214+18)*2=464 / 288 / 1
14 / 25 / (58+18)*2=152 / (194+18)*2=424 / 288 / 1
15 / 26 / (78+18)*2=192 / (174+18)*2=384 / 288 / 1
16 / 27 / (122+18)*2=280 / (138+18)*2=312 / 296 / 2
17 / 28 / (142+18)*2=320 / (118+18)*2=272 / 296 / 2
18 / 29 / (162+18)*2=360 / (98+18)*2=232 / 296 / 2
19 / 30 / (182+18)*2=400 / (78+18)*2=192 / 296 / 2
20 / 31 / (202+18)*2=440 / (58+18)*2=152 / 296 / 2
21 / 32 / (222+18)*2=480 / (38+18)*2=112 / 296 / 2
22 / 33 / (242+18)*2=520 / (18+18)*2=72 / 296 / 2

The location ranking based on the expected travel distances computed above is as follows:

So an optimal solution is:

The solution of homework #2(ISyE 6202)

# 5.3

We are considering dedicated storage policy for three products (1, 2, and 3). According to the problem, the (daily) storage requirement for each product is a discrete random variable distributed according to a uniform distribution, defined by the following probability mass functions:

Let Qi denote the number of storage locations allocated to product i (i=1,2,3). Then, we can formulate the problem as follows:

Since , we can compute the feasible range of each Qi:

Then we can solve the problem using Dynamic Programming, as we did in class. The optimal solution is Q1*=100, Q2*=150, and Q3*=200. (The detailed computation can be found in the attached tables.)

1)Q1=88 2) Q1=89 3)Q1=90 4) Q1=91

5)Q1=92 6) Q1=93 7)Q1=94 8) Q1=95

9)Q1=96 10) Q1=97 11)Q1=98 12) Q1=99

13) Q1=100

# 5.10

For this problem, we have the following conditions.

Assumptions:

  • No salvage value at the end of 500 periods
  • We do not consider operating cost of storage after construction (C1 = $ 0)

Parameters:

  • T = length of planning horizon (T=500)
  • C0 = construction cost per storage position with C0 = $63
  • C2 = leasing cost per storage position with C2 = $0.3
  • dt = demand of period t

Variables:

  • N = storage space to be constructed in order to minimize total costs over the 500 periods

Then we can formulate the problem as follows:

As it was shown in class, this is equivalent to the following optimization problem:

where

  • F = the number of periods with storage demand dt > N

This last problem indicates the following solution procedure:

Demand / 4700 / 4500 / 4200 / 4000 / 3600 / 3500 / 3400 / 3200 / 3000 / 2800 / 2750 / 2500
# of periods / 38 / 41 / 56 / 74 / 34 / 88 / 12 / 39 / 81 / 22 / 10 / 5
Cumulative # of periods / 38 / 79 / 135 / 209 / 243 / 331 / 343 / 382 / 463 / 485 / 495 / 500

The optimal number of storage locations to construct for this item is equal to the first demand value for which the corresponding partial sum exceeds 210; therefore, N* = 3600.

# 2.

Let yi (i=1,2,3) denote the number of new storage locations to be allocated to each SKU i.

Obviously, 0yi5 with y1+y2+y3=5 for all i=1,2,3. Also, let Si = Ni + yi, where Ni is the current number of locations allocated to SKU i. Then, our problem is formulated as follows:

Now we can apply Dynamic Programming for solving this problem. We can get Fi(Si) from the table of Poisson distribution (c.f., pg. 264 in Francis et. al.)

1) S3=24 2) S3=25 3) S3=26

4) S3=27 5) S3=28 6) S3=29

From the above calculation we can see that the optimal solution is S1*=14, S2*=20, S3*=25 and the maximal service level is 0.746. For comparison, the original service level, with N1=12, N2=18, N3=24, was F1(12)F2(18)F3(24) = 0.7916*0.8195*0.8432 = 0.547.

APPENDIX : Solution of # 5.3 considering the storage requirement distributions as continuous random variables.

Remember that the considered optimization problem is formulated as follows:

where, Fi is the cumulative distribution function of Xi.

is equivalent to . And since , we can compute that .------(1)

But since Fi(Qi) is a monotone increasing function with respect to Qi, and the problem is the minimization problem of the sum of continuous random variables(Qi’s), the optimal solution exists on the boundary points of the feasible region. Hence, the optimal solution satisfies the constraint as equality, and the above problem is equivalent to the following:

We can prove that optimality is attained on the boundary of the feasible region as follows: Suppose that Qi*(i=1,2,3) is an optimal solution that is not on the boundary; i.e.,. Then there exists a Qi*, let’s say Q1*, such that . But this contradicts the optimality of Qi*(i=1,2,3), i.e., the sum of Qi*’s is not minimal among the triplets of Qi’s satisfying the constraint.

Now since

we can obtain an equivalent two-variable minimization problem by solving explicitly constraint (2) for one of the problem variables, e.g., Q3. Then,

(3)

and substituting this expression in the objective we have:

s.t

50  Q1  100(4)

0  Q2  200(5)

100 200(6)

Next notice that

So z is a monotone decreasing function with respect to Q1 and for any value of Q2, it is minimized at Q1=100.

When setting Q1=100, we have

So z(100,Q2) is a convex function with respect to Q2, and it is minimized at Q2=122.47, since . But from condition (6), Q1=100 => Q2150, and since z is a monotone increasing function with respect to Q2 in [122.47, 200], it follows that: Q2=150. If we substitute Q1=100 and Q2=150 to equation (3), we obtain Q3=200. So the optimal solution is Q1*=100, Q2*=150 and Q3*=200.