例4.20轴承修理策略模拟模型之二
概率分布Page 1轴承寿命(小时)分布 延误时间(分钟)分布累计概率 轴承寿命 概率 累计概率 延误时间 概率0 1000 0.05 0 5 0.60.05 1100 0.12 0.6 10 0.30.17 1200 0.161 0.9 15 0.10.331 1300 0.1930.524 1400 0.1830.707 1500 0.1510.858 1600 0.0860.944 1700 0.0280.972 1800 0.0160.988 1900 0.0121AB C D E FGH I J K LM N1234567891011121314概率分布Page 2轴承寿命(小时) 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900概率 0.1 0.13 0.25 0.13 0.09 0.12 0.02 0.06 0.05 0.05修理工延误时间(分钟) 5 10 15概率 0.6 0.3 0.1O P Q R S T U V W X Y1234567891011121314当前轴承替换策略模拟Page 3每个轴承成本 100 1号轴承 2号轴承 3号轴承每分钟停机成本 5 寿命 累计寿命 延误 寿命 累计寿命 延误 寿命 累计寿命 延误每小时修理工工时成本 12 1 1500 1500 5 1900 1900 10 1400 1400 10换1个轴承的时间 (分) 20 2 1400 2900 5 1400 3300 10 1000 2400 53 1400 4300 5 1300 4600 5 1400 3800 10轴承更换个数 45 4 1100 5400 5 1300 5900 10 1400 5200 5延误时间(分) 350 5 1500 6900 10 1500 7400 10 1600 6800 15更换轴承停机时间 (分) 900 6 1300 8200 5 1700 9100 10 1100 7900 5轴承成本 4500 7 1400 9600 10 1500 10600 5 1700 9600 10延误时间成本 1750 8 1000 10600 10 1500 12100 5 1100 10700 5更换轴承停机成本 4500 9 1200 11800 5 1500 13600 15 1500 12200 5修理工人工成本 180 10 1300 13100 5 1300 14900 10 1200 13400 10总成本 10930 11 1000 14100 5 1800 16700 5 1500 14900 1512 1200 15300 10 1500 18200 10 1600 16500 513 1100 16400 5 1200 19400 5 1200 17700 1014 1700 18100 5 1400 20800 5 1100 18800 1515 1200 19300 10 1400 20200 516 1500 20800 51718192016 105 14 115 15 130更换次数 总延误时间当 前 轴 承 替 换 策 略 模 拟AB C D EFG H I J K L M N O P Q1234567891011121314151617181920212223242526当前轴承替换策略模拟Page 4in I5: =VLOOKUP(RAND(),LifeProbTable,2)in J5: =I5in K5: =VLOOKUP(RAND(),DelayProbTable,2)in I6: =IF(OR(J5=20000,I5=0),0,VLOOKUP(RAND(),LifeProbTable,2)in J6: =IF(I6=0,0,I6+J5)in K6: =IF(I6=0,0,VLOOKUP(RAND(),DelayProbTable,2)in I25: =COUNTIF(I5:I24,0) in K25: =SUM(K5:K24)in D8: =SUM(I25,L25,O25) in D10: =D8*D6 in D13:in D9: =SUM(K25,N25,Q25) in D11: =D8*D3 in D14:in D12: =D9*D4 in D15:RS T U V W X Y Z1234567891011121314151617181920212223242526当前轴承替换策略模拟Page 5=D10*D4=D10/60*D5=SUM(D11:D14)AA1234567891011121314151617181920212223242526新轴承替换策略模拟Page 6每个轴承成本 100 序 1号轴承 2号轴承 3号轴承 实际 累计 延误每分钟停机成本 5 号 可能寿命可能寿命可能寿命 寿命 寿命 分钟每小时修理工工时成本 12 1 1600 1400 1600 1400 1400 10同时换3个轴承的时间 (分) 40 2 1200 1500 1500 1200 2600 103 1200 1600 1000 1000 3600 10轴承更换个数 54 4 1200 1100 1500 1100 4700 10延误时间(分) 125 5 1000 1000 1100 1000 5700 5更换轴承停机时间 (分) 720 6 1500 1300 1100 1100 6800 5轴承成本 5400 7 1500 1800 1500 1500 8300 10延误时间成本 625 8 1300 1600 1000 1000 9300 5更换轴承停机成本 3600 9 1200 1000 1500 1000 10300 5修理工人工成本 144 10 1100 1100 1600 1100 11400 10总成本 9769 11 1800 1500 1500 1500 12900 512 1100 1500 1600 1100 14000 513 1100 1200 1600 1100 15100 514 1000 1300 1200 1000 16100 515 1700 1500 1700 1500 17600 516 1100 1200 1100 1100 18700 517 1000 1200 1300 1000 19700 1018 1400 1500 1100 1100 20800 51920更换次数 18总延误 125新 轴 承 替 换 策 略 模 拟AB C D E FG H I J K L M N OP1234567891011121314151617181920212223242526新轴承替换策略模拟Page 7in I5: =VLOOKUP(RAND(),LifeProbTable,2)in L5: =MIN(I5:K5)in M5: =L5in N5: =VLOOKUP(RAND(),DelayProbTable,2)in I6: =IF(OR($M5=20000,I5=0),0,VLOOKUP(RAND(),LifeProbTable,2)in M6: =IF(L6=0,0,L6+M5)in N6: =IF(L6=0,0,VLOOKUP(RAND(),DelayProbTable,2)in L25: =COUNTIF(L5:L24,0)in N25: =SUM(N5:N24)in D8: =L25*3 in D10: =L25*D6 in D13: =D10*D4in D9: =N25 in D11: =D8*D3 in D14: =D10/60*D5in D12: =D9*D4 in D15: =SUM(D11:D14)Q R S T U V W X Y1234567891011121314151617181920212223242526模拟记录Page 8当前策略 新策略当前策略总成本 10930 9,455 8,099 均值 9,719.6 7,901.8 7000新策略总成本 9769 9,804 7,731 标准差 268.9 248.1 71009,579 7,706 极小值 9131 7188 72009,829 8,199 极大值 10451 8449 7300已模拟次数 100 9,554 8,224 74009,604 7,831 外矩形 75009,455 7,781 9,131 7,188 760010,327 8,149 9,131 8,449 77009,729 7,731 10,451 8,449 78009,704 8,224 10,451 7,188 79009,903 7,906 9,131 7,188 80009,729 8,124 内矩形 81009,703 7,731 9,451 7,654 82009,654 7,781 9,451 8,150 83009,680 8,324 9,988 8,150 84009,306 7,731 9,988 7,654 85009,654 7,806 9,451 7,654 860010,277 7,188 中心折线 87009,131 7,731 9,720 5,000 88009,156 7,906 9,720 7,902 89009,679 7,756 7,000 7,902 90009,604 7,706 91009,231 8,199 92009,904 7,706 93009,555 7,731 940010,102 7,656 已模拟次数=100 95009,704 7,681 960010,053 7,731 97009,853 7,881 98009,455 7,706 99009,928 8,199 100009,853 8,249 101009,779 7,681 当前策略 新策略 102009,754 7,681 均值 9719.55 7901.80 103009,604 7,706 标准差 268.93 248.08 104009,853 7,706 105009,829 7,706 106009,903 7,781 107009,779 7,7069,654 7,7569,430 8,0999,779 8,14910,003 8,19910,451 7,7319,455 8,1999,480 7,88110,178 7,8319,331 7,7319,853 8,2241 0 0模 拟AB C D EFG H I JKL M N O PQR S123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051模拟记录Page 99,306 8,4499,430 8,1999,605 8,2499,829 7,6819,356 7,7819,605 7,7819,605 7,7069,580 7,7069,754 8,1749,531 7,8319,729 7,85610,054 8,2249,978 7,88110,128 8,1249,804 7,7319,629 7,8069,381 7,7069,954 8,1749,928 7,7569,256 8,2999,928 8,2499,255 7,8069,554 7,58110,053 8,1249,679 8,14910,178 7,8069,754 7,65610,028 8,1249,306 8,1999,580 8,0999,804 8,1999,579 7,7319,829 8,3249,878 8,1499,654 7,36310,127 7,6819,654 7,65610,327 7,6319,879 8,1749,679 7,7069,629 7,6319,555 7,6319,255 8,1499,579 8,1749,854 7,8319,728 7,7069,829 8,2249,704 7,7569,754 7,6819,480 8,24910,152 7,756AB C D EFG H I JKL M N O PQR S525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102模拟记录Page 10AB C D EFG H I JKL M N O PQR S103模拟记录Page 11新政策当前政策新政策当前政策 in N3: =IF($D$161,AVERAGE(H3:H102)in N5:0 0 0 0 in N4: =IF($D$16=1,STDEV(H3:H102)in N6:0 0 0 00 1 0 0.01 in N9: =COUNT($H:$H)0 0 0 00 1 0 0.01 in N9,N10,N13: =MIN($H$3:$H$102) in N15,N16,N19:0 0 0 0 in N11,N12: =MAX($H$3:$H$102) in N17,N18:0 1 0 0.01 in O9,O12,O13:=MIN($I$3:$I$102) in O15,O18,O19:0 12 0 0.12 in O10,O11: =MAX($I$3:$I$102) in O16,O17:0 34 0 0.340 12 0 0.12 in N23:N24: =$N$3 in O24:O25:0 2 0 0.020 3 0 0.030 21 0 0.210 10 0 0.10 2 0 0.020 1 0 0.010 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 02 0 0.02 04 0 0.04 06 0 0.06 08 0 0.08 010 0 0.1 017 0 0.17 015 0 0.15 015 0 0.15 08 0 0.08 05 0 0.05 06 0 0.06 01 0 0.01 02 0 0.02 01 0 0.01 00 0 0 00 0 0 01 1T U V W XY Z AA AB AC123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051模拟记录Page 12T U V W XY Z AA AB AC525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102模拟记录Page 13T U V W XY Z AA AB AC103模拟记录Page 14=MIN(H3:H202)=MAX(H3:H202)=IF($N$4=,$N$3-$N$4)=IF($N$4=,$N$3+$N$4)=IF($O$4=,$O$3-$O$4)=IF($O$4=,$O$3+$O$4)=$O$3AD123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051新老策略成本散布图Page 15均 值 加 减 标 准 差极 限 矩 形9 ,7 2 07 ,9 0 2新 老 策 略 成 本 散 布 图1 2 0 0 011 0 0 01 0 0 0 09 0 0 08 0 0 07 0 0 05 0 0 06 0 0 07 0 0 08 0 0 09 0 0 01 0 0 0 0当 前 策 略 成 本新策略成本已 模 拟 次 数 = 1 0 0当 前 策 略 新 策 略均 值9 7 1 9 .5 5 7 9 0 1 .8 0标 准 差2 6 8 .9 3 2 4 8 .0 81 0 0 次模 拟轴 承 成 本 = 1 0 0新老策略成本散布图Page 16当 前 策 略 新 策 略均 值9 7 1 9 .5 5 7 9 0 1 .8 0标 准 差2 6 8 .9 3 2 4 8 .0 8两 种 轴 承 替 换 策 略 下 总 成 本 的 样 本 分 布 直 方 图7 0 0 0 7 5 0 0 8 0 0 0 8 5 0 0 9 0 0 0 9 5 0 0 1 0 0 0 0 1 0 5 0 00 %1 0 %2 0 %3 0 %4 0 %总成 本当 前 政 策 新 政 策已 模 拟 次 数 = 1 0 0当 前 策 略 新 策 略均 值 9 7 1 9 .5 5 7 9 0 1 .8 0标 准 差 2 6 8 .9 3 2 4 8 .0 81 0 0 次模 拟