连云港建设工程安全网站,网页设计师相关职业前景,大学生ppt免费模板网站,个人网页设计作品代码本系列文章旨在介绍 SimPy 在工业仿真中的应用。在物流行业/工厂制造业/餐饮服务业存在大量急需优化的场景#xff0c; 例如#xff1a;如何最优化快递分拣人员的排班表以满足双十一突发的快递件量如何估算餐厅在用餐高峰的排队时长估算特定工序下#xff0c;工厂生产所需要…本系列文章旨在介绍 SimPy 在工业仿真中的应用。在物流行业/工厂制造业/餐饮服务业存在大量急需优化的场景 例如如何最优化快递分拣人员的排班表以满足双十一突发的快递件量如何估算餐厅在用餐高峰的排队时长估算特定工序下工厂生产所需要的物料成本/人力成本/时间成本这类场景无法通过常规算法求出最优解, 但是我们可以通过大量业务实践中总结出一些接近的次优解。实际生产中随时调整厂房的生产线来试验最优解是非常昂贵的。引进仿真技术可以给业务研究员无限的自由度去调整验证不同的优化方案。仿真的成本无非是计算机的算力以及程序员编写业务逻辑的时间。行业上其实已经存在一些工业仿真软件。但这类仿真软件往往针对某些特定场景高度定制化数据埋点往往不全缺乏通用的数据库接口难以结合真实业务产生的数据进行仿真这样就失去与真实业务进行比较的可能。利用 SimPy 我们可以构建一套完全开源的仿真方案可以完全私有定制业务场景。利用 Python 强大的生态仿真数据从来源到输出分析可以衔接所有开源流行的数据分析框架。目前我们已经利用 SimPy 仿真模拟物流核心分拣业务结合 MySQLTableauPandasSpark 构建一整套完整的报表可视化分析体系已经能够成功应用于现代物流中为分拣业务提供持续优化改良方案。我们创造性地解决了一些原有软件仿真中欠缺的环节这些内容将会在接下来的文章中分享。作为本系列文章的开篇我们将简要地介绍 SimPy 框架的基本理念。官方资料SimPy 是一个基于标准 Python 以进程为基础的离散事件仿真框架。SimPy 中的进程是由 Python 生成器构成生成器的特性可以模拟具有主动性的物件比如客户、汽车、或者中介等等。SimPy也提供多种类的共享资源(shared resource)来描述拥挤点(比如服务器、收银台和隧道)。仿真运行速度非常快仿真中的模拟时间长短不影响仿真运行效率仿真中的模拟时间单位可以任意指定一秒、一年、一小时都是允许的。SimPy 安装SimPy 可以同时在 Python 2 (2.7)以及 Python 3(3.2)上运行。只要有pip轻松安装。“$ pip install simpy”手动安装 SimPy 也非常方便。提取存档打开存放 SimPy 的 terminal 窗口然后输入“$ python setup.py install”你可以选择性地运行 SimPy 测试文件以了解软件是否可行。前提是要安装 pytest 包。并在 SimPy 的安装路径下运行下列命令行“$ py.test --pyargs simpy”SimPy 核心概念SimPy 是离散事件驱动的仿真库。所有活动部件例如车辆、顾客,、即便是信息都可以用 process (进程) 来模拟。这些 process 存放在 environment (环境) 。所有 process 之间以及与environment 之间的互动通过 event (事件) 来进行.process 表达为 generators (生成器) 构建event(事件)并通过 yield 语句抛出事件。当一个进程抛出事件进程会被暂停直到事件被激活(triggered)。多个进程可以等待同一个事件。 SimPy 会按照这些进程抛出的事件激活的先后 来恢复进程。其实中最重要的一类事件是 Timeout 这类事件允许一段时间后再被激活 用来表达一个进程休眠或者保持当前的状态持续指定的一段时间。这类事件通过 Environment.timeout来调用。EnvironmentEnvironment 决定仿真的起点/终点 管理仿真元素之间的关联, 主要 API 有simpy.Environment.process - 添加仿真进程simpy.Environment.event - 创建事件simpy.Environment.timeout - 提供延时(timeout)事件simpy.Environment.until - 仿真结束的条件(时间或事件)simpy.Environment.run - 仿真启动样例代码说明 API:下面是来自官方文档的两个例子:第一个例子 描述如何定义一个进程 并添加到 env 内 简单展示启动仿真的代码结构第二个例子 描述一个汽车驾驶一段时间后停车充电 汽车驾驶进程和电池充电进程通过事件的激活来相互影响Example 1import simpy# 定义一个汽车进程def car(env):while True:print(Start parking at %d % env.now)parking_duration 5yield env.timeout(parking_duration) # 进程延时 5sprint(Start driving at %d % env.now)trip_duration 2yield env.timeout(trip_duration) # 延时 2s# 仿真启动env simpy.Environment() # 实例化环境env.process(car(env)) # 添加汽车进程env.run(until15) # 设定仿真结束条件, 这里是 15s 后停止Example 2from random import seed, randintseed(23)import simpyclass EV:def __init__(self, env):self.env envself.drive_proc env.process(self.drive(env))self.bat_ctrl_proc env.process(self.bat_ctrl(env))self.bat_ctrl_reactivate env.event()self.bat_ctrl_sleep env.event()def drive(self, env):驾驶进程while True:# 驾驶 20-40 分钟print(开始驾驶 时间: , env.now)yield env.timeout(randint(20, 40))print(停止驾驶 时间: , env.now)# 停车 1-6 小时print(开始停车 时间: , env.now)self.bat_ctrl_reactivate.succeed() # 激活充电事件self.bat_ctrl_reactivate env.event()yield env.timeout(randint(60, 360)) self.bat_ctrl_sleep # 停车时间和充电程序同时都满足print(结束停车 时间:, env.now)def bat_ctrl(self, env):电池充电进程while True:print(充电程序休眠 时间:, env.now)yield self.bat_ctrl_reactivate # 休眠直到充电事件被激活print(充电程序激活 时间:, env.now)yield env.timeout(randint(30, 90))print(充电程序结束 时间:, env.now)self.bat_ctrl_sleep.succeed()self.bat_ctrl_sleep env.event()def main():env simpy.Environment()ev EV(env)env.run(until300)if __name__ __main__:main()Resource 和 StoreResource/Store 也是另外一类重要的核心概念, 但凡仿真中涉及的人力资源以及工艺上的物料消耗都会抽象用 Resource 来表达, 主要的 method 是 request. Store 处理各种优先级的队列问题, 表现跟 queue 一致, 通过 method get / put 存放 itemStore - 抽象队列simpy.Store - 存取 item 遵循仿真时间上的先到后到simpy.PriorityStore - 存取 item 遵循仿真时间上的先到后到同时考虑人为添加的优先级simpy.FilterStore - 存取 item 遵循仿真时间上的先到后到, 同时队列中存在分类, 按照不同类别进行存取simpy.Container - 表达连续/不可分的物质, 包括液体/气体的存放, 存取的是一个 float 数值Resource - 抽象资源simpy.Resource - 表达人力资源或某种限制条件, 例如某个工序可调用的工人数, 可以调用的机器数simpy.PriorityResource - 兼容Resource的功能, 添加可以插队的功能, 高优先级的进程可以优先调用资源, 但只能是在前一个被服务的进程结束以后进行插队simpy.PreemptiveResource - 兼容Resource的功能, 添加可以插队的功能, 高优先级的进程可以打断正在被服务的进程进行插队样例代码说明 API:Example 3银行排队服务例子情景:一个柜台对客户进行服务, 服务耗时, 客户等候过长会离开柜台import randomimport simpyRANDOM_SEED 42NEW_CUSTOMERS 5 # 客户数INTERVAL_CUSTOMERS 10.0 # 客户到达的间距时间MIN_PATIENCE 1 # 客户等待时间, 最小MAX_PATIENCE 3 # 客户等待时间, 最大def source(env, number, interval, counter):进程用于生成客户for i in range(number):c customer(env, Customer%02d % i, counter, time_in_bank12.0)env.process(c)t random.expovariate(1.0 / interval)yield env.timeout(t)def customer(env, name, counter, time_in_bank):一个客户表达为一个协程, 客户到达, 被服务, 然后离开arrive env.nowprint(%7.4f %s: Here I am % (arrive, name))with counter.request() as req:patience random.uniform(MIN_PATIENCE, MAX_PATIENCE)# 等待柜员服务或者超出忍耐时间离开队伍results yield req | env.timeout(patience)wait env.now - arriveif req in results:# 到达柜台print(%7.4f %s: Waited %6.3f % (env.now, name, wait))tib random.expovariate(1.0 / time_in_bank)yield env.timeout(tib)print(%7.4f %s: Finished % (env.now, name))else:# 没有服务到位print(%7.4f %s: RENEGED after %6.3f % (env.now, name, wait))# Setup and start the simulationprint(Bank renege)random.seed(RANDOM_SEED)env simpy.Environment()# Start processes and runcounter simpy.Resource(env, capacity1)env.process(source(env, NEW_CUSTOMERS, INTERVAL_CUSTOMERS, counter))env.run()Example 4# python 3.6 with SimPy工厂工序和传送带情景:一个机器处理物件, 处理完毕后放上传送带, 传送带传送一段时间后到达下个一个机器设备.[last_q][machine1] ----[con_belt]---- [next_q][machine2]import simpyimport randomPROCESS_TIME 0.5 # 处理时间CON_BELT_TIME 3 # 传送带时间WORKER_NUM 2 # 每个机器的工人数/资源数MACHINE_NUM 2 # 机器数MEAN_TIME 0.2 # 平均每个物件的到达时间间距def con_belt_process(env,con_belt_time,package,next_q):模拟传送带的行为while True:print(f{round(env.now, 2)} - item: {package} - start moving )yield env.timeout(con_belt_time) # 传送带传送时间next_q.put(package)print(f{round(env.now, 2)} - item: {package} - end moving)env.exit()def machine(env: simpy.Environment,last_q: simpy.Store,next_q: simpy.Store,machine_id: str):模拟一个机器, 一个机器就可以同时处理多少物件 取决资源数(工人数)workers simpy.Resource(env, capacityWORKER_NUM)def process(item):模拟一个工人的工作进程with workers.request() as req:yield reqyield env.timeout(PROCESS_TIME)env.process(con_belt_process(env, CON_BELT_TIME, item, next_q))print(f{round(env.now, 2)} - item: {item} - machine: {machine_id} - processed)while True:item yield last_q.get()env.process(process(item))def generate_item(env,last_q: simpy.Store,item_num: int100):模拟物件的到达for i in range(item_num):print(f{round(env.now, 2)} - item: item_{i} - created)last_q.put(fitem_{i})t random.expovariate(1 / MEAN_TIME)yield env.timeout(round(t, 1))if __name__ __main__:# 实例环境env simpy.Environment()# 设备前的物件队列last_q simpy.Store(env)next_q simpy.Store(env)env.process(generate_item(env, last_q))for i in range(MACHINE_NUM):env.process(machine(env, last_q, next_q, machine_idfm_{i}))env.run()结语文章暂时结束文章主要通过代码来展示 SimPy 的仿真能力。接下来的文章计划是介绍如何构建一个基于时间动态的人力资源排班表来模拟不同时段人力的分布比如流水线上不同岗位在不同的时间段需求的人力资源是不均等我们可以通过调岗的形式来达到人力资源调优让人力资源在时间分布上最优。介绍如何一个队列同时实现时间先后优先和优先级上优先来模拟比如银行柜台客户排队 vip 进行插队的情景启动一个翻译 SimPy 官方文档的众包计划希望在国内降低大众学习 SimPy 的语言成本本文作者
