模拟退火算法在动态设施布置中的应用研究

【本文由PB创新网为您整理】摘 要 制造业必须以运营效率和快速的反应适应产品的复杂化和需求。探讨了对生产部门间制造设施的布置和重组，以使物料搬运和重组成本最小化，进而提出动态设施布置问题及资源的有效组合和配置，确保设施间的物流通畅和提高企业生产效率。
　　关键词 模拟退火 启发式 动态设施布置
　　静态设施布置问题（static facility layout problem，SFLP）被认为是解决设施布置的有效方法。资源（如机器、部门或者劳动力）的有效组合和配置，可以确保设施间的物流通畅和提高企业生产效率。
当设施间的物流量在布置范围内变化时，SFLP就成了动态。这就是由Rosenblatt首次提出的著名的动态设施布置问题（dynamic facility layout problem，DFLP）。
　　启发式算法成功发展之前，曾经用禁忌搜索技术、遗传算法等工具来求解大型组合优化问题。换言之，就是利用最速下降成对交换启发式（steepest-descent pairwise exchange heuristic）从初始解开始产生邻域解。通常这类启发式的时间效率不高，并且只收敛于局部最优。为了克服这些缺陷，本文提出用模拟退火启发式解决DFLP最优问题，用固体退火的思想来接受邻域解，以免陷入局部最优。
1 DFLP的基本思想
　　企业随着市场的变化而调整其设施布置，本文称其为柔性布置。DFLP就是以将来的可预测变化为基础的。预期的未来可以划分为很多区段，这些区段可以定义为周、月甚至是年。研究动态设施布置问题时，设每一区段的流量数据是可预测和连续的，则设施布置问题中的每个区段，可以用SFLP（QAP）进行解决。
　　DFLP的布置规划是以可预测未来为基础的一系列布置，每个布置规划跟每个区段有关。在布置过程中，在原有基础上对设施的移动而产生的成本称为再布置成本，设施的再布置可能导致产品的损失，还可能需要专业人员和专门的设备。因此，再布置成本由劳动力成本、设备成本和产品损耗成本组成；另外，对制造设施来说，在设施间还要对物料进行搬运以满足设施加工的需要，搬运物料所投入的成本，称为物料搬运成本。成本的大小由设施间物料的流量以及设施间的距离决定，它是决定布置是否合理的最重要的衡量标准，一般占总运作成本的20%～50%，占产品制造成本的15%～70%，这成为设施布置中需要考虑的重要指标。因此，布置规划的总成本由所有区段的物料搬运成本和与再布置成本之和组成。
　　布置规划的目的是使搬运成本与再布置成本之和达到最小。在这过程中，须重复交换部门间的位置以满足上述要求，当物料搬运成本大于再布置成本时，可以把DFLP看作一系列的SFLPs（QAPs）来求解。
　　附图所示是具有6个设施在3个区段的DFLP事例。在第一区段，设施1、2、3、4、5和6分别被安置在位置3、4、1、5、2和6，由于设施3和5在第2区段被分配到不同的位置（即分别在位置2和1），则其再布置成本就是把设施3移动到位置2所产生的成本与把设施5移动到位置1所产生的成本之和。另外，在布置中，由于在阶段2和在阶段3的布置是一样的，故在阶段3中没有再布置成本。
为了操作方便，可以对DFLP进行如下的假设：设施间的流量是动态而确定的（deterministic），设施的面积和位置的大小一致，布置类型为已知的（即如附图为2×3的布置）；部门之间的距离确定为一个单位。

　　关于DFLP问题的求解，本文在Urban提出的最速下降成对交换启发式解法之上，提出用一种通用的模拟退火算法和最速下降成对交换相结合的启发式算法来求解DFLP的最优化问题。
2 模拟退火算法
　　模拟退火（Simulated Annealing，SA）算法的思想最早是由N.Metropolis等人在1953年借鉴统计力学中物质退火方法而提出的。其思想观念来自固体的退火过程，加热固体至最高温使之溶化，冷却时，液体中原子的热运动渐渐减弱，随着温度的徐徐降低，原子运动渐趋有序，达到固体的最低能量状态或者基态。根据Metropolis准则，粒子在温度T时趋于平衡的概率为e■，其中E为温度T时的内能，ΔE为其改变量，k为Boltzmann常数。
　　1982年，Kinkpatrick等人首次用模拟退火算法解决组合优化问题，将内能E模拟为目标函数值f，温度T演化成控制参数t，即得到解组合优化问题的模拟退火算法。由初始解i和控制参数初值t开始，对当前解重复“产生新解-计算目标函数差-判断是否接受-接受或舍弃”的迭代，并逐步衰减t值，算法终止时的当前解即为所得近似最优解。
　　下面给出模拟退火算法的基本步骤：
　　（1）给定模型每一个参数变化范围，在这个范围内随机选择一个初始模型m0，并计算相应的目标函数值E(m0)。
　　（2）对当前模型进行扰动产生一个新模型m，计算相应的目标函数值E(m)，得到△E=E（m）-E（m0）。
　　（3）若ΔE<0，则新模型m被接受；若ΔE>0，则新模型按概率P=exp（-△E/T）进行接受，T为温度。当模型被接受时，置m0=m，E（m0）=E（m）。
　　（4）在温度T下，重复一定次数的扰动和接受过程，即重复步骤（2）、（3）。
　　（5）缓慢降低温度T。
　　（6）重复步骤（2）、（5），直至收敛条件满足为止。
3 DFLP中的模拟退火启发式解法


3.1 参数设置
　　（1）

[1] [2]  下一页