The objective function at Equation

The objective function at Equation 8.31 aims to minimize the total cost TC of maintenance for all periods in the planning horizon. Constraints at Equations 8.32 and 8.33 respectively balance the workforce size and the maintenance workload between adjacent periods. Constraints at Equations 8.34 and 8.35 respectively relate regular time and overtime work hours to the number of regular maintenance workers in each period. Constraints at Equation 8.36 ensure that the number of assigned subcontract work hours does not exceed the available limit in each period.

Example 8.10: Maintenance workload for the next five months is 2,500, 1,500, 1800, 2800 and 2200 man-hours. This workload can be met by employees on regular time at a cost of $10 per hour, employees on overtime at a cost of $15 per hour, or subcontractors at a cost of $18 per hour. The initial workforce size is 10 employees. Each employee works for 150 regular time hours and a maximum of 60 overtime work hours per month. Maximum capacity of subcontract workers is 200 h per month. Early maintenance costs $8 per hour per month, while late maintenance costs $14 per hour per month. For each employee, hiring cost is $800 and firing cost is $1000. Assuming zero starting and ending backlog, model and solve this capacity planning problem using mathematical programming.
The integer programming model is given by
T
min TC 10Rt 15Ot 18St 8At 14Bt 800H t 1000Ft
t 1
subject to
W1 = 10 + H1 – F1
Wt = Wt – 1 + Ht – Ft, t = 2, . . , 5
A1 – B1 = R1 + O1 + S1 – 2500
A2 – B2 = A1 – B1 + R2 + O2 + S2 – 1500 A3 – B3 = A2 – B2 + R3 + O3 + S3 – 1800 A4 – B4 = A3 – B3 + R4 + O4 + S4 – 2800 0 = A4 – B4 + R5 + O5 + S5 – 2200
Rt = 150Wt, t = 1, . . , 5
Ot 60Wt, t = 1, . . , 5
St 200, t = 1, . . , 5
The optimum solution of the above model was obtained by the optimization software package LINDO. The minimum total cost TC is $120,920. Decision variables with non-zero values are shown in Table 8.11.

Table 8.11. Integer programming optimal solution of Example 8.10

Month t 1 2 3 4 5
Workforce size Wt 11 11 12 15 15
Regular time hours Rt 1650 1650 1800 2250 2250
Overtime hours Ot 660 0 0 500 0
Subcontract hours St 40 0 0 0 0
Backlogged hours Bt 150 0 0 50 0
Hired employees Ht 1 0 1 3 0

182 H.K. Al-Fares and S.O. Duffuaa


8.9 Stochastic Techniques for Capacity Planning

Stochastic models for capacity planning consider various uncertainties ever present in real-life maintenance systems. Uncertainties in maintenance surround both maintenance workload or demand (i.e., timing and severity of equipment failure) and maintenance capacity (i.e., availability and effectiveness of maintenance resources). Usually, uncertainties are represented by probability distributions with specified values of the means and variances. Stochastic models for maintenance capacity planning include queuing models, simulation models, and stochastic programming. Stochastic programming models are mathematical programming models similar to the deterministic models discussed in the previous section, except that some of their elements are probabilistic. Although these models have been used for maintenance capacity planning (e.g., Duffuaa and Al-Sultan, 1999), they are beyond the scope of this chapter, and thus will not be discussed further. In the remainder of this section, queuing theory models and computer simulation models are presented.

8.9.1 Queuing Models

Queuing models deal with systems in which customers arrive at a service facility, join a queue, wait for service, get service, and finally depart from the facility. Queuing theory is used to determine performance measures of the given system, such as average queue length, average waiting time, and average facility utilization (Taha, 2003). In addition, queuing models can be used for cost optimization by minimizing the sum of the cost of customer waiting and the cost of providing service. In applying queuing theory to maintenance systems, the maintenance jobs or required maintenance tasks are considered as the customers, and maintenance resources such as manpower and equipment are considered as the servers.
Queuing systems differ from each other in terms of several important characteristics. To define clearly the characteristics of the given queuing situation, a standard notation (Taha, 2003) is used in the following format:
(a/b/c):(d/e/f)

where


a = customer inter-arrival time distribution
b = service time (or customer departure) distribution
c = number of parallel servers
d = queue discipline, i.e., order or priority of serving customers
e = maximum number of customers allowed in the system (queue plus service)
f = size of the total potential customer population


Standard symbols are used to represent individual elements of the above notation (symbols a and b). Arrival and service distributions (symbols a and b) are represented by the symbols M (Markovian or Poisson), D (deterministic or constant), E (Erlang or Gamma), and G (general). The queue discipline (symbol d) is represented by the symbols: FCFS (first come, first served), LCFS (last come, first served), SIRO (service in random order), and GD (general discipline). The

Maintenance Forecasting and Capacity Planning 183


symbol M corresponds to the exponential or Poisson distributions. If the inter- arrival time is exponential, then the number of arrivals during a specific period is Poisson. These complementary distributions have a significant role in queuing theory because they have the Markovian (or forgetfulness) property, which makes them completely random In order to introduce specific queuing models for maintenance capacity planning, the following notation is defined:
n = number of customers in the system (queue plus service)
n = customer arrival rate with n customers in the system
n = customer departure rate with n customers in the system
 = server utilization = n /n
pn = probability of n customers in the system
Ls = expected number of customers in the system Lq = expected number of customers in the queue Ws = expected waiting time in the system
Wq = expected waiting time in the queue

Waiting time and the number of customers are directly related by Little’s Law, one of the most fundamental formulas in queuing theory:
Ls = eff Ws, or Lq = eff Wq (8.38)

where


eff = effective customer arrival rate at the system

Most queuing models are applicable to maintenance capacity planning. Two of
these models are presented below, namely the (M/M/c):(GD//) system and the (M/M/R) (GD/k/k) system.

8.9.1.1 The (M/M/c):(GD//) System
This queuing system has Markovian inter-arrival and service times, c parallel servers (repairmen), and general service disciplines. Since there are no limits on the number of customers in the system, then = eff. Defining = /, the steady- state performance measures for this system are given by
c1

Lq  p0
(c 1)!(c )2

(8.39)


where

Ls = Lq +  (8.40)

⎪⎧c1  n
p0 ⎨ 

1
c ⎪⎫
⎬,

1


(8.41)

⎪⎩n0 n!

c!(1 / c) ⎪⎭ c

The expected number for waiting time in the queue Wq and expected total time in the system Ws are respectively obtained by dividing Lq and Ls by .

The above model can be used in maintenance capacity planning to determine the optimum number of servers c (maintenance workers). In this case, the objective would be to minimize the total cost TC of waiting (i.e., cost of equipment downtime) plus the cost of providing maintenance (i.e., cost of maintenance workers). For example, this objective can be expressed as follows:

184 H.K. Al-Fares and S.O. Duffuaa





where

min TC(c) = cM c + cW Ls(c) (8.42)

cM = cost of maintenance workers per employee
cW = cost of waiting time in the queue

It should be noted that Equation 8.42 is only a typical example of a relevant
objective in maintenance capacity planning. Several alternative objective functions are possible; for instance, c could be replaced by , while Ls could be replaced by Lq, Ws, or Wq.

Example 8.11: A maintenance department repairs a large number of identical machines. Average time between failures is 2 h and 40 min, and average repair time is 5 h; both are exponentially distributed. The hourly labor cost is $15 per maintenance employee, while the hourly cost of downtime is $40 per waiting machine. Use queuing theory to determine the optimum number of maintenance employees.

 = 1/2.6667 = 0.375
 = 1/5 = 0.2
 = 0.375 /0.2 = 1.875
Since /c = 1.875/c < 1, then c > 1.875, or c 2
For c = 2, the average number of waiting machines Ls(2) and associated total cost
TC(2) are calculated by Equations 8.39–8.42 as follows:

p0 (2) ⎨1.875 

1.8752

1
⎫
⎬  0.03226

⎩n0 n!

2!(1 1.875 / 2) ⎭ 31


Lq (2) 


1.875


21

( 1 ) 421.875 13.60887

(2 1)!(2 1.875) 2 31 31
Ls(2) = 13.60887 + 1.875 = 15.48387
TC(2) = 15(2) + 40(15.48387) = 649.35
For c 3, the average number of waiting machines Ls(c) and associated total cost TC(c) are similarly calculated by Equations 8.39–8.42. Because TC(c) is convex, we should start with c = 2 and increment c by one employee at a time until the total cost TC(c) begins to increase. The calculations are summarized in Table 8.12, showing that the optimum number of maintenance employees is equal to 4.

Table 8.12. Queuing model solution of Example 8.11

c p0(c) Ls(c) TC(c)
2 0.03226 15.48387 649.35
3 0.13223 2.52066 145.83
4 0.14924 2.00265 140.11
5 0.15255 1.90328 151.13

Maintenance Forecasting and Capacity Planning 185


8.9.1.2 The (M/M/R):(GD/K/K) System
This queuing system is called the machine repair or machine servicing model. It has Markovian inter-arrival and service times, R parallel servers (repairmen), and a general service discipline.

Kết quả (Việt) 1:

Hàm mục tiêu tại phương trình 8.31 nhằm mục đích giảm thiểu tổng chi phí TC bảo trì trong tất cả thời gian ở chân trời lập kế hoạch. Các hạn chế tại 8,32 cho phương trình và 8.33 tương ứng cân bằng lực lượng lao động kích thước và khối lượng công việc bảo trì giữa thời kỳ lân cận. Các hạn chế tại phương trình 8.34 và 8,35 tương ứng liên quan thường xuyên thời gian và thời gian bù giờ làm việc giờ để số lượng công nhân bảo trì thường xuyên trong từng thời kỳ. Các hạn chế tại 8.36 phương trình đảm bảo rằng số lượng subcontract được chỉ định giờ làm việc không vượt quá giới hạn có sẵn trong từng thời kỳ.Ví dụ 8,10: Khối lượng công việc bảo trì cho năm tháng tiếp theo là 2.500, 1.500, 1800, 2800 và 2200 giờ. Khối lượng công việc này có thể được đáp ứng bởi các nhân viên thường xuyên thời gian chi phí $10 cho mỗi giờ, nhân viên trên giờ làm thêm chi phí $15 cho mỗi giờ, hoặc nhà thầu phụ chi phí $ 18 cho giờ. Kích thước của lực lượng lao động ban đầu là 10 nhân viên. Mỗi nhân viên làm việc cho 150 thường xuyên thời gian giờ và tối đa là 60 giờ làm việc thêm giờ mỗi tháng. Sức chứa tối đa của người lao động subcontract là h 200 mỗi tháng. Bảo trì đầu chi phí $8 mỗi giờ mỗi tháng, trong khi chi phí bảo trì cuối $14 cho giờ mỗi tháng. Đối với mỗi nhân viên, thuê chi phí là $800 và bắn chi phí là $1000. Giả sử không bắt đầu và kết thúc backlog, mô hình và giải quyết vấn đề lập kế hoạch khả năng này bằng cách sử dụng toán học lập trình.Mô hình lập trình số nguyên được cho bởiTMin TC 10Rt 15Ot 18St 8At 14Bt 800H t 1000Ftt 1tùy thuộc vàoW1 = 10 + H1-F1Wt = Wt-1 + Ht-Ft, t = 2,. . , 5A1-B1 = R1 O1 ++ S1-2500A2-B2 = A1-B1 + R2 + O2 + S2-1500 A3-B3 = A2-B2 R3 + O3 ++ S3-1800 A4-B4 = A3-B3 R4 + O4 ++ S4-2800 0 = A4-B4 R5 + O5 + S5-2200RT = 150Wt, t = 1,. . , 5OT 60Wt, t = 1,. . , 5St 200, t = 1,. . , 5Các giải pháp tối ưu của các mô hình ở trên thu được bằng cách tối ưu hóa phần mềm gói LINDO. Tối thiểu tổng chi phí TC là $120,920. Các biến quyết định với giá trị không được hiển thị trong bảng 8.11.8.11 bảng. Số nguyên lập trình các giải pháp tối ưu của ví dụ 8,10Tháng t 1 2 3 4 5Lực lượng lao động kích thước Wt 11 11 12 15 15Thường xuyên thời gian giờ Rt 1650 1650 1800 2250 2250Thêm giờ giờ Ot 660 0 0 500 0Subcontract giờ St 40 0 0 0 0Trệ giờ Bt 150 0 0 50 0Thuê nhân viên Ht 1 0 1 3 0 182 H.K. Al-giá vé và so Duffuaa8.9 các kỹ thuật ngẫu nhiên lập kế hoạch nâng cao năngCác mô hình ngẫu nhiên năng lực lập kế hoạch cho xem xét sự không chắc chắn khác nhau đã từng có trong cuộc sống thực bảo trì hệ thống. Sự không chắc chắn trong duy trì bao quanh cả hai khối lượng công việc bảo trì hoặc nhu cầu (tức là, thời gian và mức độ nghiêm trọng của thiết bị thất bại) và khả năng bảo dưỡng (tức là, tính khả dụng và hiệu quả của tài nguyên bảo trì). Thông thường, sự không chắc chắn được đại diện bởi phân bố xác suất với các giá trị được chỉ định của các phương tiện và chênh lệch. Các mô hình ngẫu nhiên để bảo trì năng lực lập kế hoạch bao gồm xếp hàng các mô hình, mô hình mô phỏng, và lập trình ngẫu nhiên. Mô hình lập trình ngẫu nhiên là toán học lập trình mô hình tương tự như các mô hình xác định, thảo luận trong phần trước, ngoại trừ rằng một số yếu tố của họ là xác suất. Mặc dù những mô hình đã được sử dụng để bảo trì năng lực lập kế hoạch (ví dụ như, Duffuaa và Al-Sultan, 1999), họ vượt ra ngoài phạm vi của chương này, và do đó sẽ không được thảo luận thêm. Trong phần còn lại của phần này, xếp hàng các mô hình lý thuyết và mô hình mô phỏng máy tính cũng được trình bày trong phòng.8.9.1 mô hình xếp hàngMô hình xếp hàng đối phó với hệ thống trong đó khách hàng đến một cơ sở dịch vụ, tham gia một hàng đợi, chờ đợi cho dịch vụ, nhận được dịch vụ, và cuối cùng rời khỏi cơ sở. Xếp hàng lý thuyết được sử dụng để xác định các biện pháp hiệu suất của hệ thống nhất định, chẳng hạn như chiều dài hàng đợi trung bình, Trung bình là chờ đợi thời gian, và sử dụng trung bình cơ sở (Tran, 2003). Ngoài ra, xếp hàng các mô hình có thể được sử dụng để tối ưu hóa chi phí bằng cách giảm thiểu tổng chi phí của khách hàng chờ đợi và chi phí của việc cung cấp dịch vụ. Trong việc áp dụng lý thuyết xếp hàng để bảo trì hệ thống, các công việc bảo trì hoặc nhiệm vụ yêu cầu bảo trì được coi là các khách hàng, và duy trì nguồn tài nguyên như nguồn nhân lực và thiết bị được coi là các máy chủ.Hệ thống xếp hàng khác nhau trong điều khoản của một số đặc điểm quan trọng. Để xác định rõ ràng các đặc tính của tình huống nhất định xếp hàng, một ký hiệu chuẩn (Tran, 2003) được sử dụng trong các định dạng sau:(a/b/c):(d/e/f) nơi một = khách hàng thời gian liên đến phân phốib = thời gian phục vụ (hoặc khách hàng khởi hành) phân phốic = số của song song máy chủd = kỷ luật hàng đợi, tức là, đơn đặt hàng hoặc ưu tiên phục vụ khách hànge = số tối đa của khách hàng cho phép trong hệ thống (hàng đợi cộng với dịch vụ)f = kích thước tổng dân số khách hàng tiềm năng Biểu tượng tiêu chuẩn được sử dụng để đại diện cho các yếu tố cá nhân của các ký hiệu trên (biểu tượng một và b). Đến và dịch vụ phân phối (biểu tượng một và b) được đại diện bởi biểu tượng M (Markovian hoặc Poisson), D (xác định hoặc không đổi), E (Erlang hoặc Gamma) và G (tổng quát). Kỷ luật hàng đợi (biểu tượng d) được đại diện bởi các biểu tượng: FCFS (lần đầu tiên đến, lần đầu tiên phục vụ), LCFS (cuối đến, lần đầu tiên phục vụ), SIRO (dịch vụ theo thứ tự ngẫu nhiên), và GD (kỷ luật chung). Các Bảo trì dự báo và năng lực lập kế hoạch 183biểu tượng M tương ứng với các mũ hoặc phân phối Poisson. Nếu thời gian đến inter là mũ, sau đó số lượng khách đến trong một khoảng thời gian cụ thể là Poisson. Các bản phân phối bổ sung có một vai trò quan trọng trong lý thuyết xếp hàng vì họ có tài sản Markovian (hoặc forgetfulness), mà làm cho chúng hoàn toàn ngẫu nhiên để giới thiệu mô hình xếp hàng cụ thể cho khả năng bảo trì có kế hoạch, các ký hiệu sau được định nghĩa:n = số lượng khách hàng trong hệ thống (hàng đợi cộng với dịch vụ)n = tỷ lệ khách hàng đến với n khách hàng trong hệ thốngn = tỷ lệ khởi hành khách hàng với khách hàng n trong hệ thống = máy chủ sử dụng = n /nPN = xác suất của n khách trong hệ thốngLs = số dự kiến của khách hàng trong hệ thống Lq = số dự kiến của các khách hàng trong hàng đợi Ws = thời gian chờ đợi dự kiến trong hệ thốngWQ = dự kiến thời gian chờ đợi trong hàng đợiThời gian chờ đợi và số lượng khách hàng trực tiếp liên quan đến ít của pháp luật, một trong những công thức cơ bản nhất trong lý thuyết xếp hàng:Ls = eff Ws, hoặc Lq = eff Wq (8,38) nơi EFF = tỷ lệ đến khách hàng có hiệu quả tại hệ thống Đặt xếp hàng các mô hình được áp dụng để bảo trì năng lực lập kế hoạch. Hai trong sốCác mô hình này được trình bày dưới đây, cụ thể là các (hệ thống M/M/c):(GD//) và hệ thống (M/M/R) (GD/k/k).8.9.1.1 (M/M/c):(GD//) hệ thốngHệ thống xếp hàng này đã Markovian liên đến và dịch vụ thời gian, c máy chủ song song (repairmen), và các ngành dịch vụ tổng hợp. Kể từ khi không có không có giới hạn về số lượng khách trong hệ thống, sau đó  = eff. Xác định  = /, các biện pháp hiệu suất trạng thái ổn định cho hệ thống này được cho bởic1 LQ  p0(c 1)! (c ) 2 (8,39) nơi Ls = Lq +  (8,40) ⎪⎧c1  nP0 ⎨ đột 1c ⎪⎫⎬, 1 (8.41) ⎪⎩n0 n! c! (1  / c) ⎪⎭ c Số dự kiến cho thời gian chờ đợi trong hàng đợi Wq và dự kiến tổng số thời gian trong hệ thống Ws tương ứng được thu được bằng cách chia Lq và Ls bởi .Các mô hình ở trên có thể được sử dụng trong bảo trì năng lực lập kế hoạch để xác định số máy chủ c (bảo trì công nhân), tối ưu. Trong trường hợp này, mục tiêu sẽ là để giảm thiểu tổng chi phí TC chờ đợi (tức là, các chi phí của thiết bị downtime) cộng với chi phí của việc cung cấp bảo dưỡng (tức là, các chi phí của công nhân bảo trì). Ví dụ, mục tiêu này có thể được thể hiện như sau: 184 H.K. Al-giá vé và so Duffuaa nơi Min TC(c) = cM c + cW Ls(c) (8.42)cM = chi phí bảo trì công nhân cho mỗi nhân viêncW = chi phí thời gian chờ đợi trong hàng đợi Cần lưu ý rằng phương trình 8.42 là chỉ là một ví dụ điển hình của bản có liên quanmục tiêu trong bảo trì năng lực lập kế hoạch. Một số chức năng thay thế mục tiêu là có thể; Ví dụ, c có thể được thay thế bởi , trong khi Ls có thể được thay thế bởi Lq, Ws hoặc Wq.Ví dụ 8.11: Một bộ phận bảo trì sửa chữa một số lớn các máy giống hệt nhau. Thời gian trung bình giữa thất bại là 2 giờ 40 phút, và thời gian trung bình sửa chữa là 5 h; cả hai được phân phối theo cấp số nhân. Chi phí lao động theo giờ là $15 cho mỗi nhân viên bảo trì, trong khi chi phí theo giờ của thời gian chết là $40 cho mỗi chờ máy. Sử dụng lý thuyết xếp hàng để xác định số nhân viên bảo trì, tối ưu. = 1/2.6667 = 0.375 = 1/5 = 0,2 = 0.375 /0.2 = 1.875Kể từ /c = 1.875/c < 1, sau đó c > 1.875, hoặc c 2Cho c = 2, số trung bình của chờ đợi máy Ls(2) và liên quan đến tổng chi phíTC(2) được tính bằng phương trình 8,39-8.42 như sau: P0 (2) ⎨1.875 đột 1.8752 1⎫⎬  0.03226 ⎩n0 n! 2! (1 1.875 / 2) ⎭ 31 LQ (2)  1.875 21 (1) 421.875 13.60887 (2 1)! (2 1.875) 2 31 31LS(2) = 13.60887 + 1.875 = 15.48387TC(2) = 15(2) + 40(15.48387) = 649.35Cho c 3, số trung bình của chờ đợi máy Ls(c) và liên quan đến tổng chi phí TC(c) tương tự như vậy được tính bằng phương trình 8,39-8.42. Bởi vì TC(c) là lồi, chúng ta nên bắt đầu với c = 2 và tăng c bởi một nhân viên tại một thời gian cho đến khi tổng chi phí TC(c) bắt đầu tăng. Các tính toán được tóm tắt trong bảng 8,12, Hiển thị số nhân viên bảo trì, tối ưu là tương đương với 4.Bảng 8,12. Giải pháp mô hình xếp hàng của ví dụ 8.11c p0(c) Ls(c) TC(c)2 0.03226 15.48387 649.353 0.13223 2.52066 145.834 0.14924 2.00265 140.115 0.15255 1.90328 151.13 Bảo trì dự báo và năng lực lập kế hoạch 1858.9.1.2 (M/M/R):(GD/K/K) hệ thốngHệ thống xếp hàng này được gọi là sửa chữa máy hoặc máy phục vụ mô hình. Đô thị này có Markovian liên đến và thời gian dịch vụ, R song song máy chủ (repairmen), và một kỷ luật dịch vụ tổng hợp.

Kết quả (Việt) 2:

Hàm mục tiêu ở phương trình 8.31 nhằm mục đích giảm thiểu tổng chi phí TC bảo trì cho tất cả các giai đoạn trong chu kỳ kế hoạch. Ngại ở phương trình 8.32 và 8.33 tương ứng cân bằng kích thước lực lượng lao động và khối lượng công việc bảo trì giữa các kỳ kế cận. Ngại ở phương trình 8.34 và 8.35 tương ứng liên quan thời gian và công việc làm thêm giờ thường xuyên với số lượng công nhân bảo trì thường xuyên trong từng thời kỳ. Ngại ở phương trình 8.36 đảm bảo rằng số giờ làm việc được giao thầu lại không vượt quá giới hạn có sẵn trong từng thời kỳ. Ví dụ 8.10: Bảo trì khối lượng công việc cho những năm tháng tiếp theo là 2.500, 1.500, 1800, 2800 và 2200 người đàn ông giờ. Khối lượng công việc này có thể được đáp ứng bởi nhân viên trên thời gian thường xuyên với chi phí $ 10 mỗi giờ, nhân viên làm thêm giờ với chi phí $ 15 một giờ, hoặc nhà thầu phụ với giá $ 18 mỗi giờ. Kích thước ban đầu lực lượng lao động là 10 nhân viên. Mỗi nhân viên làm việc cho 150 giờ thời gian thường xuyên và tối đa là 60 giờ làm thêm mỗi tháng. Công suất tối đa của người lao động hợp đồng phụ là 200 h mỗi tháng. Bảo dưỡng đầu chi phí $ 8 mỗi giờ mỗi tháng, trong khi bảo trì muộn chi phí $ 14 mỗi giờ mỗi tháng. Đối với mỗi người lao động, tuyển dụng chi phí là $ 800 và chi phí sa thải là $ 1000. Giả sử không khởi đầu và kết thúc công việc tồn đọng, mô hình và giải quyết vấn đề quy hoạch năng lực này sử dụng lập trình toán học. Các mô hình lập trình số nguyên được đưa ra bởi T min TC 10Rt 15Ot 18St 8At 14Bt 800H t  1000Ft t 1 chịu W1 = 10 + H1 - F1 Wt = Wt - 1 + Ht - Ft, t = 2,. . , 5 A1 - B1 = R1 + O1 + S1 - 2500 A2 - B2 = A1 - B1 + R2 + O2 + S2 - 1500 A3 - B3 = A2 - B2 + R3 + O3 + S3 - 1800 A4 - B4 = A3 - B3 + R4 + O4 + S4 - 2800 0 = A4 - B4 + R5 + O5 + S5 - 2200 Rt = 150Wt, t = 1,. . , 5 Ot 60Wt, t = 1,. . , 5 St 200, t = 1,. . 5 Giải pháp tối ưu của các mô hình trên được thu thập bởi các LINDO gói phần mềm tối ưu hóa. Các TC tổng chi phí tối thiểu là $ 120,920. Quyết định biến với các giá trị khác không được thể hiện trong Bảng 8.11. Bảng 8.11. Lập trình Integer giải pháp tối ưu của Ví dụ 8.10 Tháng t 1 2 3 4 5 kích thước Workforce Wt 11 11 12 15 15 giờ thời gian thường xuyên Rt 1650 1650 1800 2250 2250 Giờ làm thêm Ot 660 0 0 500 0 giờ đồng phụ St 40 0 0 0 0 giờ Backlogged Bt 150 0 0 50 0 nhân viên thuê HT 1 0 1 3 0 182 HK Al-Giá vé và SO Duffuaa Kỹ thuật 8.9 Stochastic cho năng lực lập kế hoạch mô hình Stochastic cho kế hoạch năng lực xem xét những bất ổn khác nhau bao giờ có mặt trong các hệ thống bảo trì thực tế cuộc sống. Không chắc chắn trong việc bảo trì bao quanh cả khối lượng công việc bảo dưỡng hay nhu cầu (tức là, thời gian và mức độ nghiêm trọng của thiết bị thất bại) và khả năng bảo trì (tức là, sẵn sàng và hiệu quả của các nguồn lực bảo trì). Thông thường, sự không chắc chắn được đại diện bởi phân bố xác suất với giá trị quy định của các phương tiện và phương sai. Mô hình ngẫu nhiên cho kế hoạch năng lực bảo dưỡng bao gồm các mô hình xếp hàng, mô hình mô phỏng, và lập trình ngẫu nhiên. Mô hình lập trình Stochastic là mô hình lập trình toán học tương tự như mô hình xác định được thảo luận trong các phần trước, ngoại trừ một số yếu tố của họ là xác suất. Mặc dù các mô hình này đã được sử dụng để lập kế hoạch bảo trì công suất (ví dụ, Duffuaa và Al-Sultan, 1999), họ là vượt ra ngoài phạm vi của chương này, và do đó sẽ không được thảo luận thêm. Trong phần còn lại của mục này, xếp hàng mô hình lý thuyết và mô hình máy tính mô phỏng được trình bày. 8.9.1 Queuing Models Xếp hàng mô hình thỏa thuận với các hệ thống mà trong đó khách hàng đến một cơ sở dịch vụ, tham gia vào một hàng đợi, chờ đợi cho dịch vụ, nhận được dịch vụ, và cuối cùng khởi hành từ cơ sở. Lý thuyết hàng đợi được sử dụng để xác định các biện pháp thực hiện của hệ thống nhất định, chẳng hạn như độ dài trung bình hàng đợi, thời gian chờ đợi trung bình, và sử dụng cơ sở trung bình (Taha, 2003). Ngoài ra, các mô hình xếp hàng có thể được sử dụng để tối ưu hóa chi phí bằng cách giảm thiểu các khoản chi phí của khách hàng chờ đợi và chi phí cung cấp dịch vụ. Khi áp dụng lý thuyết xếp hàng để bảo trì hệ thống, các công việc bảo trì hoặc các công việc bảo dưỡng cần thiết được coi là khách hàng, và các nguồn lực bảo trì như nhân lực và trang thiết bị được coi là các máy chủ. Hệ thống xếp hàng khác nhau về một số đặc điểm quan trọng. Để xác định rõ đặc điểm tình hình xếp hàng đưa ra, một ký hiệu chuẩn (Taha, 2003) được sử dụng trong các định dạng sau: (a / b / c) :( d / e / f) nơi a = khách hàng liên đến phân phối thời gian b = thời gian dịch vụ (hoặc khởi hành của khách hàng) phân phối c = số lượng máy chủ song song d = xếp hàng kỷ luật, tức là, trật tự hoặc ưu tiên phục vụ khách hàng e = số lượng tối đa của khách hàng cho phép trong hệ thống (hàng đợi cộng với dịch vụ) f = kích thước của tổng số dân số khách hàng tiềm năng các biểu tượng tiêu chuẩn được sử dụng để đại diện cho yếu tố cá nhân của các ký hiệu trên (ký hiệu a và b). Đến và phân phối dịch vụ (ký hiệu a, b) được đại diện bởi các biểu tượng M (Markovian hoặc Poisson), D (xác định hoặc không đổi), E (Erlang hoặc Gamma) và G (nói chung). Kỷ luật hàng đợi (ký hiệu d) được đại diện bởi các biểu tượng: FCFS (đầu tiên đến trước được phục vụ), LCFS (cuối cùng đến, đầu tiên phục vụ), SIRO (dịch vụ trong thứ tự ngẫu nhiên), và GD (kỷ luật chung). Việc bảo trì dự báo và năng lực lập kế hoạch 183 biểu tượng M tương ứng với sự phân bố hàm mũ hay Poisson. Nếu thời gian đến tế là hàm mũ, sau đó số lượng khách đến trong một khoảng thời gian cụ thể là Poisson. Các bản phân phối bổ sung có vai trò quan trọng trong lý thuyết hàng đợi, vì họ có Markovian (hay quên) tài sản, mà làm cho họ hoàn toàn ngẫu nhiên Nhằm giới thiệu các mô hình xếp hàng cụ thể cho kế hoạch năng lực bảo trì, các ký hiệu sau đây được định nghĩa: n = số lượng khách hàng trong hệ thống (hàng đợi cộng với dịch vụ) n = khách hàng tốc độ đến với n khách hàng trong hệ thống n = tỷ lệ khách hàng rời với n khách hàng trong hệ thống  = sử dụng máy chủ = n / n pn = xác suất của n khách hàng trong hệ thống Ls = dự kiến số lượng khách hàng trong hệ thống Lq = số lượng dự kiến của khách hàng trong Ws queue = dự kiến thời gian chờ đợi trong hệ thống WQ = dự kiến thời gian chờ đợi trong hàng đợi thời gian và số lượng khách hàng có liên quan trực tiếp bởi Luật nhỏ, một trong những chờ các công thức cơ bản nhất trong lý thuyết hàng đợi: Ls = eff Ws, hoặc Lq = eff WQ (8.38), nơi eff = hiệu quả tốc độ đến khách hàng tại hệ thống các Hầu hết các mô hình xếp hàng được áp dụng để lập kế hoạch bảo trì công suất. Hai trong số các mô hình này được trình bày dưới đây, cụ thể là (M / M / c) :( GD / k / k hệ thống (GD /  / ) hệ thống và (M / M / R)). 8.9.1.1 Các (M / M / c) :( GD /  / ) Hệ thống xếp hàng này có liên đến và dịch vụ lần Markovian, c máy chủ song song (thợ sửa chữa), và các ngành dịch vụ nói chung. Vì không có giới hạn về số lượng khách hàng trong hệ thống, sau đó  = eff. Định nghĩa  =  / , các biện pháp thực hiện nhà nước steady- cho hệ thống này được đưa ra bởi c1 Lq  p0 (c 1)! (C ) 2 (8.39), nơi Ls = Lq +  (8.40) ⎪⎧c1  n p0  ⎨ 1 c ⎪⎫ ⎬, 1 (8,41) ⎪⎩n0 n! c! (1   / c) ⎪⎭ c Số lượng dự kiến cho thời gian chờ đợi trong hàng đợi WQ và dự kiến tổng thời gian trong Ws hệ thống tương ứng được thu được bằng cách chia Lq và Ls bởi . Các mô hình trên có thể được sử dụng trong kế hoạch năng lực bảo trì để xác định tối ưu số lượng máy chủ c (công nhân bảo trì). Trong trường hợp này, mục tiêu là để giảm thiểu tổng chi phí TC chờ đợi (tức là, chi phí thiết bị downtime) cộng với các chi phí cung cấp bảo trì (tức là, chi phí của công nhân bảo trì). Ví dụ, mục tiêu này có thể được thể hiện như sau: 184 HK Al-Giá vé và SO Duffuaa nơi min TC (c) = cM c + CW Ls (c) (8.42) cM = chi phí của công nhân bảo trì cho mỗi nhân viên Cw = chi phí chờ đợi thời gian trong hàng đợi Cần lưu ý rằng phương trình 8.42 chỉ là một ví dụ điển hình của một liên quan trong kế hoạch năng lực bảo trì. Một số chức năng mục tiêu thay thế có thể xảy ra; Ví dụ, c có thể được thay thế bằng , trong khi Ls có thể được thay thế bởi Lq, Ws, hoặc WQ. Ví dụ 8.11: Một bảo trì sửa chữa bộ phận một số lượng lớn các máy giống hệt nhau. Thời gian trung bình giữa các sự cố là 2 h và 40 phút, và thời gian sửa chữa trung bình là 5 h; cả hai đều được phân theo cấp số nhân. Các chi phí lao động theo giờ là $ 15 cho mỗi nhân viên bảo trì, trong khi chi phí theo giờ của thời gian chết là $ 40 cho mỗi máy chờ đợi. Sử dụng lý thuyết xếp hàng để xác định số lượng tối ưu của các nhân viên bảo trì.  = 1 / 2,6667 = 0,375  = 1/5 = 0,2  = 0,375 /0.2 = 1,875 Từ  / c = 1,875 / c <1, sau đó c> 1,875, hoặc c 2 Đối với c = 2, số lượng trung bình của chờ đợi máy Ls (2) và chi phí liên quan tổng TC (2) được tính bằng phương trình 8,39-8,42 như sau: p0 (2) ⎨1.875  1,8752 1 ⎫ ⎬  0.03226 ⎩n0 n! 2! (1 1.875 / 2) ⎭ 31 Lq (2)  1,875 21 (1) 421.875 13.60887 (2 1) ! (2 1.875) 2 31 31 Ls (2) = 13,60887 + 1,875 = 15,48387 TC (2) = 15 (2) + 40 (15,48387) = 649,35 Đối với c 3, số con trung bình của chờ đợi máy Ls ( c) và liên TC tổng chi phí (c) được tính tương tự của phương trình 8,39-8,42. Bởi vì TC (c) là lồi, chúng ta nên bắt đầu với c = 2 và tăng c bởi một nhân viên tại một thời gian cho đến khi tổng chi phí TC (c) bắt đầu tăng. Các tính toán được tóm tắt trong bảng 8.12, cho thấy rằng số lượng tối ưu của các nhân viên bảo trì bằng 4. Bảng 8.12. Xếp hàng giải pháp mô hình của Ví dụ 8.11 c p0 (c) Ls (c) TC (c) 2 0,03226 15,48387 649,35 3 0,13223 2,52066 145,83 4 0,14924 2,00265 140,11 5 0,15255 1,90328 151,13 bảo trì dự báo và năng lực lập kế hoạch 185 8.9.1.2 Các (M / M / R) :( GD / / K) Hệ thống K hệ thống sắp hàng này được gọi là sửa chữa máy hoặc máy mô hình phục vụ. Nó có Markovian liên đến và dịch vụ lần, R server song song (thợ sửa chữa), và một kỷ luật dịch vụ tổng hợp.