STELLA销售存储模拟模型的建立

1、STELLA销售存储模拟模型的建、/ :五、 STELLA 销售存储模拟模型的建立 下面 结合实际操作介绍如何利用 STELLA 建立并模拟 一个存储管理系统的模型，说明怎样使用 STELLA 软件建模，通过学习，了解 STELLA 的某些功 能。然而，要在一个模型中详细介绍 STELLA 的 全部功能是不可能的，如有需要请查阅参考文献 2， 3。研究存储模型，必须要有仓库，现在就建 立 一个。在 STELLA 环境下，选择累积量，将 其 命名为“ kucun 。本模型中为研究物资的流 动， 任一时刻的库存量是有进库和出库决定 的，进 库由订货决策决定， 出库由需求量决定。 目前模 型中的库存

2、还无法改变，因为它还没有 流进和流 出。库存的流进和流出为进货和售出量，选择 流 率项，使它们分别流进和流出库存，命名为 “jinhuo “ shouchuliang 。现在，我们把注 意力集中到影响进货和售出的逻辑关系上。进货 决策是怎样作出的呢？先要有一个每次进货量和平安库存量，然后比较库存量和平安库存量，女口 果库存量少于平安库存量，那么作出订货决策对销 售的影响因素是需求量。销售量是否一定 等于需 求量呢？实际情况是当产生一个需求量 时，首先 比较需求量和库存量，当需求量小于 库存量时， 销售量等于需求量；当需求量大于 库存量时，销 售量等于库存量。所以三者之间 也要建立联系选 择相应

3、元素作连线的以下图：用前一节中讲述的方法消灭问号。在实际工作中，货物延迟是普遍现象，在2、. "1 ArA上述模型中参加延迟，模型如以下图:延迟的添加，这里需要特殊说明一下：在STELLA环境下选择累积量， 建成一个累积变 量 并命名为“ del 表示延迟双击该变量，出现 如下的对话窗:单项选择按钮“ Reservoir 是累积变量，即我们 最常用的变量。 “Conveyor 为输送带，具 有延迟的行为，是本例中我们用到的一种， “ Queun “ Oven 分别是排队和处理箱，本例 中暂不用，有兴趣的话可以探索一下它的使 用方法 我们将单项选择按钮选中“ Conveyor 将随后出

4、现的 对话窗中的“ Transit Time ：栏中数值改为 “ 4，意为货物自订货日期始 4 天后到货， “Inflow Limit ：改为 60，“ Capacity ： 也 改为 60，意为货物“ ITIAL(del) 的初始值定 为“ 0，最终的对话窗为以下图：最后单击“ OK 退出，定义完毕。到货延迟期的大小对平均库存量和缺货损 失都有影响，对不同的到货延迟期，应制定不 同 的平安库存量和进货量。下面三个图为同一 库 存策略在不同延迟期下的库存状态，一个将支付过多的保管费，个将承受缺货损失。阵件量"< * * * «* 旨*if!几f1 1* - * ,*l

5、，*fi* 1I14-V-411 IiZ : 'i(, iILuI14lJ VII人A. J Ilip|-.|fhLin1(1IhlXI .1 11bi 'I% ?1iJWinlII11|nk j1 . 1 1 1 11II1F1 i 1 1j* -riA * 9- A 1 , a I *1 1X a *. .» L h J*11ItIIIfi1厂 HIi V1 11BI11 11 1 ;、 |1 1f11I11E耳r11tV h14tI1H|I1y: d ii!* ri <i n I m.i -ii Bidi I Ii d ia m | fMr1 n 仃 n

6、J 6 iI7 1 LI UU ir0.0025 DO 第二天刮赏50JC7i00；紀娥ICO DO库存筆阵存量1I41750 00I加如1L0 W25 OD第囚真刮赏75,00WO W模拟天数ID 00呻YDOO25.00第六天到橫75 OO10D0D模拟天数对销售量的随机性模拟是通过计算机产生 随机数来实现的，可调用相应的函数来完成。对 某些商品其日销售量的随机分布可近似看作 是 正态分布。对一些中小型农机具其日销售量 是 偏态的。而对一些大型农机具如联合收割机 其销售量的随机分布具有特殊的形式，需单独考虑。综合考虑有关因素，建立存储模拟模型如 图示：模拟过程可简述如下：模拟开始，先检查

7、这 一天是否为预定到货期，如果是那么原有库存 量增 加进货量，并把预定到货量清零。如果这 一天不 是预定到货期，那么库存量不变。然后检 查预期库 存量实际库存量加预订到货量是 否小于平安 库存量，如果是那么需要重新订货， 这时需要增加 一项订货费用，并按到货延迟期预订到货期。按给定的分布规律，由计算机产生随机数 模 拟随机需求量，并将需求量与库存量比较， 假设需求量小于库存量，那么库存量减去需求量,21: 13 2003Table-2-25并对剩余库存计算一次得保管费；假设需求 量等 于库存量，那么库存量为零；假设需求量大于库 存 量，那么库存量为零， 同时计算相应的缺货损 失。在一定模拟时间

8、内，日复一日进行模拟， 最 后将销售收入减去订货总费用，库存保管费 用以 及由于缺货造成的损失即得总收益。通过 比较不 同策略的模拟结果，以效益最好的一个 方案作为 决策方案。模拟结果可通过图表的形式显示如下：日期012345678910111213141516数台2164.00.3912.008213.00.757.00 f1342.00.205495.00.847.0012012.00.6311.003498.00.9612.007071.00.675.008635.00.148003.00.771.005870.00台40.0036.0024.0021.0014.0012.007.006

9、0.0048.0037.0029.0017.0016.0011.006.0063.0062.00.24820.55台台0.000.000.000.0060.000.000.000.000.000.000.000.000.000.000.000.000.000.000.000.0060.000.000.000.000.000.000.000.000.0060.000.000.000.000.0017 392.796.0062.000.000.0018593.958.0056.000.000.0019554.337.0048.000.000.0020643.329.0041.000.000.004

10、8608.458.0040.000.000.0049213.404.0032.0060.000.00Final 571.52 28.00a：D ODd9 DO<30 DO0 oo<i nr2475IS9 0035 DO n血 1207575 DOGraph 126 00m；恥21STELLA 方程:cost(t) = cost(t - dt) + (c) * dt INIT cost =0INFLOWS:c = stock*unit_costdel(t) = del(t - dt) + (order - r) * dtINIT del = 0TRANSIT TIME = 4INFL

11、OW LIMIT = 60CAPACITY = 60INFLOWS:order if(stock<safety_AND(t=1)then(order_volume)else(0)OUTFLOWS: r = CONVEYOR OUTFLOW expend(t) = expend(t - dt) + (e) * dt INIT expend = 0INFLOWS:if(order>0)then(1000)else(0)loss(t) = loss(t - dt) + (l) * dt INIT loss =0 INFLOWS: l if(stock<require)then(re

12、quire-stock)*4000)e lse(0)profit(t) = profit(t - dt) + (income - e- c) * dtINIT profit = 0INFLOWS:income = sell*priceOUTFLOWS:eif(order>0)then(1000)else(0)c = stock*unit_cost stock(t)= stock(t - dt) + (r - sell) * dt INIT stock =40INFLOWS:CONVEYORrOUTFLOWOUTFLOWS:sellIF(stock>=require)THEN(req

13、uire)else(stock)order_volume = 60 price = 6680 randam = RANDOM(0,1000) require=if(randam<=39)then(0)else( if(r andam<=59)then(1)else(if(ran dam<=129) then(2)else(if(randam<=209)then( 3)el se(if(randam<=229)then(4)els e(if(randam<=379)then(5)else(if(r andam<=469)then(6)else(if(ra

14、nda m<=559)then(7)else(if (randam<=619)then(8)else(if(ran dam<=679)then(9)else(if(randam<=759)then(10) else(if(randam<=819)then(11)else( if(randam< =869)then(12)else(if(r andam<=899)then(13)else(if(rand am<=919)then(14)else(if(randam <=929)then(15)else(if (randam<=949)then(16)e