使用Python的SimPy進(jìn)行制造仿真

pip install simpy

import simpywood_capacity = 1000initial_wood = 500dispatch_capacity = 500
class Guitar_Factory:    def __init__(self, env):        self.wood = simpy.Container(env, capacity = wood_capacity, init = initial_wood)        self.dispatch = simpy.Container(env ,capacity = dispatch_capacity, init = 0)

def body_maker(env, guitar_factory):    while True:        yield guitar_factory.wood.get(1)        body_time = 1        yield env.timeout(body_time)        yield guitar_factory.dispatch.put(1)
def neck_maker(env, guitar_factory):    while True:        yield guitar_factory.wood.get(1)        neck_time = 1        yield env.timeout(neck_time)        yield guitar_factory.dispatch.put(2)

hours = 8days = 5#定義仿真的時(shí)間長度total_time = hours * days
env = simpy.Environment()guitar_factory = Guitar_Factory(env)
body_maker_process = env.process(body_maker(env, guitar_factory))neck_maker_process = env.process(neck_maker(env, guitar_factory))print('仿真開始：')env.run(until = total_time)print(f'倉庫Dispatch中分別有%d 琴身和琴頸的庫存!' % guitar_factory.dispatch.level)print('仿真結(jié)束。')

electronic_capacity = 100initial_electronic = 100
pre_paint_capacity = 100post_paint_capacity = 200
class Guitar_Factory:    def __init__(self, env):        self.wood = simpy.Container(env, capacity = wood_capacity, init = initial_wood)        self.electronic = simpy.Container(env, capacity = electronic_capacity, init = initial_electronic)        self.pre_paint = simpy.Container(env, capacity = pre_paint_capacity, init = 0)        self.post_paint = simpy.Container(env, capacity = post_paint_capacity, init = 0)        self.dispatch = simpy.Container(env ,capacity = dispatch_capacity, init = 0)

def painter(env, guitar_factory):    while True:        yield guitar_factory.pre_paint.get(10)        paint_time = 4        yield env.timeout(paint_time)        yield guitar_factory.post_paint.put(10)
def assembler(env, guitar_factory):    while True:        yield guitar_factory.post_paint.get(2)        yield guitar_factory.electronic.get(1)        assembling_time = 1        yield env.timeout(assembling_time)        yield guitar_factory.dispatch.put(1)

env = simpy.Environment()guitar_factory = Guitar_Factory(env)
body_maker_process = env.process(body_maker(env, guitar_factory))neck_maker_process = env.process(neck_maker(env, guitar_factory))painter_process = env.process(painter(env, guitar_factory))assembler_process = env.process(assembler(env, guitar_factory))
print(f'仿真開始：')
env.run(until = total_time)
print(f'噴涂前有%d 琴身和琴頸準(zhǔn)備噴涂' % guitar_factory.pre_paint.level)print(f'噴涂后有 %d 琴身和琴頸準(zhǔn)備組裝' % guitar_factory.post_paint.level)print(f'有 %d 吉他成品!' % guitar_factory.dispatch.level)print(f'----------------------------------')print(f'仿真完成。')

import random
num_body = 2mean_body = 1std_body = 0.1
num_neck = 1mean_neck = 1std_neck = 0.2
num_paint = 1mean_paint = 4std_paint = 0.3
num_ensam = 4mean_ensam = 1std_ensam = 0.2

def body_maker(env, guitar_factory):    while True:        yield guitar_factory.wood.get(1)        body_time = random.gauss(mean_body, std_body)        yield env.timeout(body_time)        yield guitar_factory.pre_paint.put(1)
def neck_maker(env, guitar_factory):    while True:        yield guitar_factory.wood.get(1)        neck_time = random.gauss(mean_neck, std_neck)        yield env.timeout(neck_time)        yield guitar_factory.pre_paint.put(2)
def painter(env, guitar_factory):    while True:        yield guitar_factory.pre_paint.get(10)        paint_time = random.gauss(mean_paint, std_paint)        yield env.timeout(paint_time)        yield guitar_factory.post_paint.put(10)
def assembler(env, guitar_factory):    while True:        yield guitar_factory.post_paint.get(1)        yield guitar_factory.electronic.get(1)        assembling_time = max(random.gauss(mean_ensam, std_ensam), 1)        yield env.timeout(assembling_time)        yield guitar_factory.dispatch.put(1)之前的例子timeout的參數(shù)我們總是設(shè)定為一個(gè)固定值，現(xiàn)在我們將這幾個(gè)工序的生產(chǎn)時(shí)間定義為一個(gè)隨機(jī)數(shù)。
random.gauss(mean_body, std_body)

random.gauss(mean_body, std_body)

def body_maker_gen(env, guitar_factory):    for i in range(num_body):        env.process(body_maker(env, guitar_factory))        yield env.timeout(0)
body_gen = env.process(body_maker_gen(env, guitar_factory))

wood_critial_stock = (((8/mean_body) * num_body +(8/mean_neck) * num_neck) * 3) electronic_critical_stock = (8/mean_ensam) * num_ensam * 2

class Guitar_Factory:    def __init__(self, env):        self.wood = simpy.Container(env, capacity = wood_capacity, init = initial_wood)        self.wood_control = env.process(self.wood_stock_control(env))        self.electronic = simpy.Container(env, capacity = electronic_capacity, init = initial_electronic)        self.electronic_control = env.process(self.electronic_stock_control(env))        self.pre_paint = simpy.Container(env, capacity = pre_paint_capacity, init = 0)        self.post_paint = simpy.Container(env, capacity = post_paint_capacity, init = 0)        self.dispatch = simpy.Container(env ,capacity = dispatch_capacity, init = 0)

    def wood_stock_control(self, env):        yield env.timeout(0)        while True:            if self.wood.level <= wood_critial_stock:                print(f'在第{0}日 第{1}小時(shí)，木材庫存 ({2})低于預(yù)警庫存水平下 '.format(                    int(env.now/8), env.now % 8,self.wood.level))                print('聯(lián)系供應(yīng)商')                print('----------------------------------')                yield env.timeout(16)                print('在第{0}天 第{1}小時(shí)，木材送達(dá)'.format(int(env.now/8), env.now % 8))                yield self.wood.put(300)                print('當(dāng)前庫存是：{0}'.format(                    self.wood.level))                print('----------------------------------')                yield env.timeout(8)            else:                yield env.timeout(1)

guitars_made = 0
body_pre_paint_capacity = 60neck_pre_paint_capacity = 60body_post_paint_capacity = 120neck_post_paint_capacity = 120
class Guitar_Factory:    def __init__(self, env):        self.wood = simpy.Container(env, capacity = wood_capacity, init = initial_wood)        self.wood_control = env.process(self.wood_stock_control(env))        self.electronic = simpy.Container(env, capacity = electronic_capacity, init = initial_electronic)        self.electronic_control = env.process(self.electronic_stock_control(env))        self.body_pre_paint = simpy.Container(env, capacity = body_pre_paint_capacity, init = 0)        self.neck_pre_paint = simpy.Container(env, capacity = neck_pre_paint_capacity, init = 0)        self.body_post_paint = simpy.Container(env, capacity = body_post_paint_capacity, init = 0)        self.neck_post_paint = simpy.Container(env, capacity = neck_post_paint_capacity, init = 0)        self.dispatch = simpy.Container(env ,capacity = dispatch_capacity, init = 0)        self.dispatch_control = env.process(self.dispatch_guitars_control(env))

def body_maker(env, guitar_factory):    while True:        yield guitar_factory.wood.get(1)        body_time = random.gauss(mean_body, std_body)        yield env.timeout(body_time)        yield guitar_factory.body_pre_paint.put(1)

def painter(env, guitar_factory):    while True:        yield guitar_factory.body_pre_paint.get(5)        yield guitar_factory.neck_pre_paint.get(5)        paint_time = random.gauss(mean_paint, std_paint)        yield env.timeout(paint_time)        yield guitar_factory.body_post_paint.put(5)        yield guitar_factory.neck_post_paint.put(5)

def dispatch_guitars_control(self, env):        global guitars_made        yield env.timeout(0)        while True:            if self.dispatch.level >= 50:                print('成品庫存為：{0}, 在第{1}日 第{2}小時(shí) 聯(lián)系了商場(chǎng)取貨'.format(self.dispatch.level, int(env.now/8), env.now % 8))                print('----------------------------------')                yield env.timeout(4)                print('在第{0}日 第{1}小時(shí)，商場(chǎng)取走{2}吉他'.format(int(env.now/8), env.now % 8,self.dispatch.level))                guitars_made += self.dispatch.level                yield self.dispatch.get(self.dispatch.level)                print('----------------------------------')                yield env.timeout(8)            else:                yield env.timeout(1)

def assembler(env, guitar_factory):    while True:        yield guitar_factory.body_post_paint.get(1)        yield guitar_factory.neck_post_paint.get(1)        yield guitar_factory.electronic.get(1)        assembling_time = max(random.gauss(mean_ensam, std_ensam), 1)        yield env.timeout(assembling_time)        yield guitar_factory.dispatch.put(1)

print('當(dāng)前等待噴涂的琴身數(shù)量：{0} 和琴頸數(shù)量：{1}'.format(    guitar_factory.body_pre_paint.level, guitar_factory.neck_pre_paint.level))print('當(dāng)前等待組裝的琴身數(shù)量：{0} 和琴頸數(shù)量：{1}'.format(    guitar_factory.body_post_paint.level, guitar_factory.neck_post_paint.level))print(f'當(dāng)前成品庫存量：%d ' % guitar_factory.dispatch.level)print(f'----------------------------------')print('此周期的吉他總生產(chǎn)數(shù)量: {0}'.format(guitars_made + guitar_factory.dispatch.level))print(f'----------------------------------')print(f'仿真完成！')

class Guitar_Factory:    def __init__(self, env):        self.wood = simpy.Container(env, capacity = wood_capacity, init = initial_wood)        self.wood_control = env.process(self.wood_stock_control(env))        self.electronic = simpy.Container(env, capacity = electronic_capacity, init = initial_electronic)        self.electronic_control = env.process(self.electronic_stock_control(env))        self.body_pre_paint = simpy.Container(env, capacity = body_pre_paint_capacity, init = 0)        self.neck_pre_paint = simpy.Container(env, capacity = neck_pre_paint_capacity, init = 0)        self.body_post_paint = simpy.Container(env, capacity = body_post_paint_capacity, init = 0)        self.neck_post_paint = simpy.Container(env, capacity = neck_post_paint_capacity, init = 0)        self.dispatch = simpy.Container(env ,capacity = dispatch_capacity, init = 0)        self.dispatch_control = env.process(self.dispatch_guitars_control(env))        self.env_status_monitor = env.process(self.env_status(env))

def env_status(self, env):  global status  status = pd.DataFrame(columns = ["datetime", "dispatch_level",'wood','electronic','body_pre_paint','neck_pre_paint','body_post_paint','neck_post_paint'])  status[["datetime", "dispatch_level",'wood','electronic','body_pre_paint','neck_pre_paint','body_post_paint','neck_post_paint']] = status[["datetime", "dispatch_level",'wood','electronic','body_pre_paint','neck_pre_paint','body_post_paint','neck_post_paint']].astype(int)  while True:    im = plt.plot(status['datetime'], status['neck_pre_paint'],color='#4D9221')    ims.append(im)    plt.title('neck_pre_paint')    plt.pause(0.001)    print('{0}在第{1}日 第{2}小時(shí)，成品庫存量：{3}'.format(env.now,int(env.now/8), env.now % 8,self.dispatch.level))    tmp = {'datetime':env.now,         'dispatch_level':self.dispatch.level,         'wood':self.wood.level,         'electronic':self.electronic.level,         'body_pre_paint':self.body_pre_paint.level,         'neck_pre_paint':self.neck_pre_paint.level,         'body_post_paint':self.body_post_paint.level,         'neck_post_paint':self.neck_post_paint.level         }    status = status.append([tmp])
    yield env.timeout(1)