基于仿真优化的 MRO调度问题研究
□ 丁金想，栾世超，石晓磊，连福诗
(中国航空综合技术研究所，北京 100028)
【摘 要】 与传统的制造业相比，MRO(Maintenance,Repair and Overhaul，即航空维修，包括飞机维护、修理和翻修)系统面临更多的不确定性和复杂性，如拆分-修理-组装的三级结构、待维修件质量水平的不确定、物料匹配需求以及不确定工艺路线和工时，使得MRO系统的计划和调度问题更加难以处理。在对MRO调度问题详细分析的基础上，文中建立以最小化期望权重延误时间为目标的仿真模型，并应用调度规则、NEH和GA三种启发式算法对仿真模型进行求解，最后，通过基于实际背景的算例验证了模型的可行性并对比分析了多种仿真优化算法的特点，文中的结论也从不同角度给实际MRO企业调度管理者提供决策参考。
【关键词】 航空维修；调度；不确定工艺路线；不确定工时；仿真优化
1 引言
在航空产业链中，全球民用航空运输市场的发展带动了民用航空维修业(MRO)和民用航空制造业的发展，民用航空维修业的经济价值越来越受到关注和青睐。从行业类型来看，MRO企业属于服务行业，民航公司是其最大客户群，与民航运输业呈相互依存关系[1]。
对于MRO企业来说，除定期维修外，更多的维修工作存在不确定性，例如，飞机故障发生时间、原因、程度不可预测，导致维修的工作内容、时间和周期等无法按例行工作安排，因而MRO企业面临诸多不确定因素，调查显示，47%的MRO公司维修过程中面临着40%-59%不确定工作内容[2]。虽然近年来国内外航空维修市场发展迅速，但是生产管理却面临着很多难以有效解决的问题，如何有效地解决维修生产调度中的问题对于MRO企业提高生产效率、提升生产能力、优化库存管理以及提高客户满意度等方面具有重要意义。
2 国内外研究现状
MRO企业是再制造工业的重要组成部分，一些复杂的特点使得MRO计划与调度问题与传统制造业有很大不同，例如：如拆分-修理-组装的三级结构、待维修件质量水平的不确定、物料匹配需求以及不确定工艺路线和工时。MRO是生产计划与调度的一种新的研究领域，它的复杂性和不确定性引起学术界的广泛关注。
仿真优化是一种适合解决含有不确定性调度问题的有效方法，Guide[3]等人研究了再制造的调度策略，分解机制和优先调度规则被建立在仿真模型中，通过大量的对比实验证明了最早交期(EDD)的调度规则往往能够得到很不错的结果，而分解机制的不同对调度效果的影响比较小。Guide等人[4]以最短交货期为目标，研究了调度规则对再制造系统的订单执行策略的影响，结果表明，通过采用简单的调度规则能显著提高再制造系统的性能指标。Christoph[5]等人将研究目标定位于寻找适合修理车间job-shop运作的调度规则，首先假设使用调度规则可以提高调度结果，考虑到MRO系统的复杂特点，文章通过建立仿真模型，将准时交货率、在制品率、利用率以及生产周期作为衡量标准，结果表明，调度规则对MRO的调度效果的确有一定的影响，并且发现FIFO/Slack和ESD/Slack的规则组合会取得更好的效果。
Kim[6]等人首先分析了再制造系统的复杂特点，考虑到分解-修理-组装的三级结构，并将决策变量定义为待维修件在分解车间中被分解的顺序、子部件在修理车间中每个工作站前的修理顺序以及待维修件在组装车间被组装的顺序，文章建立基于仿真的整数规划模型，通过分析，认为该模型是NP-hard问题，文章提出基于调度规则的启发式算法、基于NEH的启发式算法以及IG(Iterated greedy)算法进行求解，并通过大量算例对比分析了不同启发式算法的优缺点。然而，模型中针对修理车间的建模是确定的，他们认为子部件在修理车间有固定的维修工艺路线，此外也没有考虑到工时的不确定性。
郝[1]等人在对航空维修企业调度问题国内外研究现状详细分析的基础上，建立以最小化期望权重延误时间为目标的混合整数线性规划模型，因为该问题是NP-hard问题，通过使用传统的优化方法很难求得最优解，他们提出基于多精度仿真模型的MO2TOS算法对该问题求解，并通过基于实际背景的算例验证了模型的可行性。然而，当问题规模增大时，由于解空间巨大，MO2TOS算法将不再适合求解该问题。本文的研究背景是在该研究基础上开展的，通过基于仿真优化的方法解决大规模问题情况下的MRO调度问题。
3 MRO调度问题建模
3.1 MRO问题描述
一般的MRO系统基本框架包含三个车间，为了简便，即使现实中每个车间可能包含多个，但我们考虑每种车间只含有一个，三个车间各自主要的工序内容如下：
根据水资源承载力多目标情景分析法的边界条件和基本原理，山塘水资源承载力的情景方案组合有108种，本文选取未来较可能出现的3种情景组合进行计算，变量如表2所示。
“嚼舌根”虽是一种劣习，但它却反映了当今社会老人受到忽视产生的畸形心理。其实，老人更需要存在感。子女应时时处处善待老人，家中大小事情请他们一起商量、做决定，经常带他们四处走走、给他们买礼物……让老人常“露面”、常“做主”，老人便会有“存在感”。 而老年人在享受生活时，也需要找寻一些具有正能量的兴趣爱好，多参加集体的公益活动，树立正确的人生观和价值观。
其中OP1和OP8分别表示分解车间和组装车间，剩余工序表示修理车间。假设调度开始时刻，n个待维修件已经可以被安排调度，所有待维修件在每个机器上的加工顺序是相同的，这个顺序也是我们要决策的，其中OP4和OP5各自拥有两个不同的机器处理不同的工作内容，并且各自拥有一个概率值，这体现出MRO系统在修理车间的工艺路线不确定性。基于实际的维修企业信息，并不是所有的工序的加工时间都是不确定的，例如，清洗工序。因此，我们只假设OP6和OP7有不确定的加工时间，其他工序加工时间是确定的。
模糊PID控制器采用单片机编程设计，由于MSP430单片机内部没有专用的浮点数处理器，因此在数据的处理过程中，浮点数的计算是通过特定的算法程序来实现，如果采用浮点数计算来进行数据处理，将消耗大量的CPU资源，同时数据的处理周期较长，影响其单片机的实时控制，因此在数据处理时应尽量少用实型数据计算处理。在实际设计中将浮点数的小数部分放大，在满足精度要求的基础上，尽可能采用整形数据来处理数据计算，也可以采用长整形来实现数据的处理（见图4）。
步骤7：把当前序列Seq作为NEH算法得到的最优调度顺序，并输出该顺序，算法结束。
待维修件首先在分解车间根据BOM结构分解形成若干子部件，每个子部件的检测后状态可能是修理、完好、报废三种状态中的一种，每种维修状态表示子部件在修理车间不同的维修工艺路线，因此，每个子部件的实际工艺路线是不知道的，直到待维修件完成分解和检测工序；此外，因为子部件损坏程度不同，即使对于两个相同的子部件，其在同一个工序上的加工时间也存在不确定性。
可见，整个维修过程受到很多约束条件，在我们的问题中，对于N个给定的待维修件，决策变量是待维修件在分解车间分解的先后顺序、拆分出的子部件在修理车间每个工序上的修理顺序、以及子部件在组装车间的组装顺序。优化目标是期望权重延误时间之和最小。
1.2.2 研究组 标本采集和测定方法同对照组。在新生儿疾病筛查管理中，应用网络技术，构建新生儿疾病筛查信息化管理平台[4-7]，对新生儿进行管理，具体方法如下：①管理平台的建立。利用互联网技术和数据库技术，助产单位使用条形码技术将新生儿基础信息、标本采集资料输入系统；②建立新生儿疾病筛查三级网络：由助产机构(采血单位)、市县新生儿疾病筛查分中心(辖区妇幼保健机构)、省新生儿疾病筛查中心组成，完善新生儿遗传代谢病筛查、确诊、治疗等流程，进一步落实新生儿疾病筛查管理工作；③筛查结果查询。通过系统记录，可以了解每一个新生儿疾病筛查进展工作，包括筛查结果、确诊情况和治疗情况等。
3.2 MRO调度问题的仿真模型
如上文所述，MRO问题不确定性因素多(不确定达到时间，不确定加工时间，不确定工艺路线等)，问题结构复杂，为了实现对该问题的求解，结合作者在一家航空维修企业的项目经历，为不失一般性，本文对所研究的这类MRO调度问题做出以下假设：
信息网络渠道是互动行为产生的重要动因。新媒体环境下，沟通交流离不开信息网络[9]。围绕公众号推文及线上直播教学模块，公众号粉丝群体建立和维护了互动交流渠道，通过这个渠道增加公众号热度，从而进一步提升平台的活跃度和影响关注度。运营方可以在今后的推文以及线上直播教学模块中，鼓励用户留言、点赞、转发、弹幕等行为，并通过抽奖或者评选最佳留言的方式，来回馈活跃用户；同时，通过开通微社区功能，为学者之间或名师与学者之间提供一个相互分享学习及比赛等经验的渠道；此外还可以增加跟踪报道国内大型体育舞蹈赛事的栏目，并在其中设置粉丝投票环节，从而充分提高用户互动性。
整个MRO系统只由一个分解车间、一个修理车间和一个组装车间构成；
分解车间和组装车间分别由一个工序组成，每种工序只包含一个机器；
修理车间有多种工序类型，每个子部件在修理车间的工艺路线是不确定的；
修理车间每个工序只包含一个机器；
每个拆分出的子部件都只能组装到原主部件上；
每个机器同一时间只能加工一个部件。
细纸筒一个，粗纸筒一个，纸圈两个，黑色套圈两个，黑色圆孔观察盖一个，彩色珠子若干，长方形镜片三个（也可以用能够反光的镜面纸，一般的快餐盒或礼品盒内都可以找到），彩色塑料片一片，透明塑料片一片，胶水、胶带、画笔，画纸若干。
首先我们基于C++建立以最小化期望权重延误时间为目标的仿真模型，示意图如图1所示。
图1 MRO系统结构图
修理车间：包含很多不同的工序和机器，一般认为是job-shop类型，子部件在该车间的工艺流程不确定；
4 基于仿真优化的 MRO调度问题算法研究
4.1 调度规则
调度规则算法是一种简单、易操作的启发式算法，因为其具有容易理解和方便实现的特点，在实际中有着广泛的应用。在我们问题中，所有待维修件或子部件在分解车间、修理车间和组装车间内所有工序上的加工顺序需要决策。每个工序前都可以通过指定一种规则来实现调度，根据现有文献中对调度规则的研究[6]，我们对常见并且有效的调度规则分为两类，如表1所示。其中，第一类调度规则可以考虑到待维修件在其维修工艺路线上所有工序的加工特点，而第二类调度规则主要针对某个特定工序。
表1 两类调度规则
对于所有待维修件在分解车间中分解工序前的加工顺序，我们采用第一类调度规则；对于子部件在修理车间以及组装车间中工序前的加工顺序，我们采用第二类调度规则。显然为了实现MRO企业维修现场的生产调度，选择何种调度规则组合是我们要研究的问题。
4.2 NEH算法
在构造型启发式算法中，Nawaz-Enscore-Ham(NEH)[7]被认为是计算以make-span为目标的排列flow-shop问题的一种非常高效的算法[6][8]。基于NEH算法的基本思路，Fernandez-Viagas[9]等人提出一种适用于以延误时间为目标的改进NEH算法。这种新的NEH算法与传统NEH算法有两方面不同：第一，以EDD规则得到的解作为NEH的初始解；第二，通过计算总延误时间来评价一个部分顺序。
加强对公众的宣传，改变传统的“租不如买”的旧思想，树立“轻占有、重使用”的新观念。由于目前融资租赁业务涉及的主要是大型固定资产，很多行业没有机会参与到租赁业务中，对融资租赁认识不足。首先，企业或政府提供平台，向东部地区或西方国家学习其先进的租赁经验。其次，融资租赁企业定期开展公司展览活动，向外部人员介绍融资租赁业的基本常识和发展趋势。另外，还可以政企合作，成立专门的咨询机构，对企业、政府进行融资租赁只是培训，加强政府对租赁行业的重视程度，拓宽企业参与融资租赁业务渠道。
对于一个给定的顺序，在NEH算法计算期望权重延误时间之和时，我们首先根据该加工顺序对应的待维修件的工时，计算得到每个待维修件在分解车间的完工时间，也即可以进入修理车间加工的时间。对于修理车间的任何一个工序，子部件的加工顺序通过调度规则给出，根据每个子部件的维修工艺路径和调度规则，可以计算得到每个子部件离开修理车间的时间；对于组装车间内工序前的加工顺序，我们根据子部件是否全部到达组装车间以及全部到达的时间来决定组装顺序。最后，我们通过计算
给出该顺序的期望权重延误时间之和。该算法的具体步骤如下：
因此，如何激发中药学专业本科生对于中医药文化的兴趣，提高其中医药文化底蕴、培养其中医药思维方式，是我们中医药教育工作者需要思考、研究和探索的重要课题。笔者所在中医药学院对此进行了深入探索与思考，近年来开展了针对新生的中医药文化游学活动，现介绍如下，以期与相关教育工作者探讨。
本文运用的治疗方法是国际公认的液氮冷冻治疗。具体的治疗步骤是结合患者的临床情况，有针对性的使用棉签或冷冻枪治疗，当疣体四周出现2 mm左右的白晕后停止冷冻。一次治疗给患者进行2～4个冻融，每20～30天之内只能对患者进行一次治疗，对每个患者最多进行四次治疗。
步骤1：通过上节调度规则中给出的某种第一类调度规则得到一个初始顺序，记为Seq0=(J1,J2,…,Jn)；
步骤2：设置k=2，从Seq0中选择前两个待维修件J1和J2，分别形成序列(J1,J2)和(J2,J1)，计算各自的tardiness，以较小tardiness的序列作为当前序列Seq，假设Seq=J1,J2,…,Jn；
其中，γ是负常数，std(Y )表示Y的标准差，εLMD为局部均值差异，LSD1和LSD2为局部方差差异。当曲线内部局部灰度方差大于外部时，即std(L M1(X ) )＞std(L M2(X ) )，可以认为曲线内部匀质性低于外部匀质性，此时，应增加曲线内部能量所占权重，来降低其对曲线演化的影响；同理std(L M1(X ) )＜std(L M2(X ) )时，增加曲线外部能量所占比重，当两者相差较小时，能量函数达到最优。综上所述，本文模型的总能量泛函被定义为：
步骤3：设置k=3，如果k≤N，转向步骤4；否则转向步骤7；
步骤4：从Seq0中选择第k个待维修件；
步骤5：将待维修件Jk分别插入到序列Seq中的k个可能位置，得到k个新的序列，如：{Jk,J1,J2,…,Jk-1}，{J1,Jk,J2,…,Jk-1},…,{J1,J2,…,Jk-1),Jk}。计算各自的tardiness，把最小tardiness的序列作为当前序列Seq；如果有多个相同的最小tardiness，则通过最小总空闲时间的选择机制进行选择。
步骤6：设置k=k+1，转向步骤3。
组装车间：组装和最后的检验。
分解车间：分解维修件、清洗以及检测；
4.3 遗传算法(GA)
NEH算法虽然效率高，但只是一个局部搜索的构造型启发式算法，为了获得更好的解，我们需要具有全局搜索特性的改进型启发式算法。在生产计划与调度领域，遗传算法(GA)已经被广泛应用于求解NP-hard组合优化问题[10][11][12]。基于生物进化的“优胜劣汰”原理，GA通过选择、交叉和变异操作来实现群体并行优化。作为一种通用的全局优化算法，GA近年来在各领域得到广泛应用，但大量研究表明GA存在易早熟、算法参数敏感等缺点，取得良好的性能需要依赖较大的种群并对算法参数作精心设计[13]。
在遗传算法中，我们直接通过顺序编号作为染色体编码。针对本文提出的问题，与NEH算法一样，通过GA求解待维修件在分解车间的加工顺序，对于修理车间和组装车间内工序前的加工顺序，我们通过调度规则得到加工顺序。
染色体的表示：因为本文模型是求解分解车间待维修件的加工顺序，因此可以直接用自然数顺序编码。假设当前有8个待维修件，编号分别为1,2,3,4,5,6,7,8。则某一排列(染色体)可为：86 2 3 7 5 1 4。
针对MRO仿真模型，我们通过NEH求解待维修件在分解车间的加工顺序，对于修理车间和组装车间内工序前的加工顺序，我们通过调度规则得到加工顺序。与一般NEH算法相比，本文提出的NEH算法会考虑不同调度规则得到的解作为NEH初始解，与修理车间和组装车间的调度规则的最佳组合。此外，在算法执行过程中，对于具有相同目标函数的两个顺序，我们通过使用最小总空闲时间[10]的选择机制进行选择，而非随机选择。
交叉操作方法：本文采用典型的PMX方法，假如有两个个体A和B，A：8 62 3 7 5 1 4，B：1 5 4 2 8 6 3 7。PMX 操作方法的基本过程如下：首先在A、B 染色体上随机选取一段，例如2,3,7 和4,2,8，用该段内的基因对应关系来决定一系列的交换，杂交变换后的新个体A和B分别为：7 64 2 8 5 1 3，1 5 2 3 7 6 4 8。
该算法的具体步骤如下：
步骤5：如果t>3*n，转向步骤11；否则，转向步骤6；
步骤2：设置t=0；通过随机产生的方式初始化种群S(0)；
步骤3：通过计算适应度函数来评价种群S(0)中所有个体；
∀i∈si
其中，
为顺序
的目标函数值。
超声波测距模块用来测量模块距离地面的距离d;MPU-9250模块用来测量拐杖运动的角速度w。通过大量的实验来模拟老人摔倒时的状况,发现当d>240 cm、w>5 rad/s且最终测量角度大于80°时,有99.2%的情况老人处于跌倒状态,将此作为判定老人跌倒的标志[9],同时控制GSM模块将报警信息发送至远程的手机监测软件,其软件流程如10所示。
在辽宁，监管部门要求银行将“公平信贷”原则贯穿信贷管理全程，从业务受理、客户评级、放款约束、贷后管理等方面，检查是否存在可能影响公平信贷的政策，取消不合理的条款、要求和做法。在重庆，监管部门划分机构、片区，组织机构负责人和一线客户经理，对全市6000余家民营工业企业进行全覆盖式走访、回访，对民企广泛关注的困难和问题拿出切实举措，目前已现场为民企解决各类融资问题30个，落实资金24.6亿元。在浙江，监管部门组织银行对全省授信10亿元以上企业进行筛选，探索建立优质企业“白名单”，对暂时遇到困难但仍有发展前景的企业，主动开展帮扶，助力企业提前做好债券兑付安排，推动市场化法治化债转股落地。
步骤4：计算选择概率；
步骤1：设置参数，根据文献对GA参数的研究，我们设置种群大小为2.5*n，迭代次数为3*n，交叉和变异概率分别为0.75和0.25。
步骤6：设置t=t+1；
步骤7：从种群S(t-1)选择新的种群S(t)；
步骤8：交叉操作，交叉操作使用PMX方法[10][14]；
步骤9：变异；
步骤10：评价新的种群S(t)；
步骤11：把当前种群中适应度值最高的顺序作为GA算法得到的调度顺序，输出该顺序，算法结束。
上面给出的GA算法的初始种群完全是随机产生的，本文对这种传统的GA算法进行改进，改进的GA算法的初始种群中有一个个体是NEH算法的解，并且有0.25*n个个体是以NEH的解为基础随机变异产生的，在GA迭代过程中，始终保证每代种群中最优秀的个体被保留下来。
4.4 算例设计与结果分析
为了研究不同算法在处理MRO调度问题上的效果，本节我们设计一系列算例进行验证，并对计算结果进行对比分析。这设计算例时，我们需要考虑以下四个变量：待维修件数量、每个待维修件可以拆分的子部件数量、每个子部件可能的工艺路线个数以及每个子部件所有可能工艺路线上的最大工序个数。
本文给出的算例中，待维修件数量可以为(6,10,15,30,50)，每个待维修件可以拆分出的子部件个数为(2,3)，对于每一个子部件，其可能的工艺路线为(2,3)，而对于每个子部件，其所有工艺路线中最大工序个数为(3,5,7)。根据这些参数，我们可以随机产生60个算例，为了表述方便，对其进行编号如表2所示：
表2 60个算例编号
在随机产生算例时，考虑到模型中工时和工艺路线是不确定的，所以，我们将每个工序的工时设定为一个三点分布，例(2,4,5)，而这三个值相应的概率从U(0,1)中随机产生，例(0.1,0.6,0.3)。其中U(a,b)表示一个在区间[a,b]内的均匀分布。具体来讲，基于一定的实际数据，每个待维修件在修理车间任一工序的工时从U(2,6)中产生，而在分解车间和组装车间任一工序的工时从U(1,3)中产生。
乌龟、兔子的同时比赛借助多线程变为了现实。实践生活所出现的多线程状况也比较常见，包括：很多飞机或火车售票口在相同时间卖票，各窗口表示独立线程。对于此，存在不同线程对相同资源进行竞争的问题，代表不同的窗口存在对相同票进行售卖的可能性。所以在处理过程中，就需要运用多线程同步的技术。
每个待维修件的交期(due date)从U(P(1-T-R/2),P(1-T+R/2))中随机产生[15]，其中T和R分别表示延迟参数和交期范围参数，并分别设置为0.37和0.26，P是make-span的下界，每个待维修件的延迟权重(λ值)从U(1,2)中随机产生。
本文所有算法程序都使用C++编码，并在英特尔酷睿i3，主频3.40GHz的CPU和4.00GB内存的64位win7电脑上进行运算。对于每个解，我们设置仿真次数为1000以消除噪音获得稳定目标函数值。
由于篇幅有限，不失一般性，我们只展示编号为(6,12,18,24,30,36,42,48,54,60)的算例的计算结果。为了比较不同算法的表现，两种测量标准在本文被采用，分别是平均相对性能比率(ARPR)和平均提高百分比(API)。
其中，Heuristici表示对于第i个算例，该启发式算法计算得到的目标函数值，Besti表示对于第i个算例，所有启发式算法得到的最好值。为了比较不同算法的计算效果，我们计算每个目标函数值的均值和方差，最优性比较是基于假设检验，下文所有结果的置信度都是95%。
调度规则：表1给出了调度规则算法的表现，因为第一类调度规则有9种，第二类调度规则有6中，所以一共具有54种组合。因篇幅限制，本文根据平均ARPR值的大小只给出最好3种和最差3种的组合。其中EDD-FIFO表示分解车间的调度规则是EDD，而修理和组装车间的调度规则是FIFO。
从表3可以看出，EDD-EDD、EDD-FIFO和STPT-FIFO三种调度规则普遍来说表现比较好，但是并不能说其中哪一种调度规则绝对比其他更好。随着问题规模增大，显然EDD-EDD调度规则逐渐更具有优势。此外，调度规则算法计算速度很快，即使对于包含50个待维修件的算例-60，其计算时间也在1秒之内，这也是调度规则能够在实际中广泛使用的原因之一。
表3 ARPR结果对于调度规则算法
NEH算法：NEH算法的计算时间要比调度规则大，对于算例60，NEH算法的计算时间约为8分钟。虽然NEH算法的计算比较耗时，但是其解的质量相对于调度规则也有一定的提高，表4给出NEH算法与同种调度规则相比解的API值。从表4可以看出，当把调度规则作为NEH的初始解时，NEH算法获得的解的质量有很大的提高，最大的提高值可以达到23.51%(如NEH-LTRT-FIFO)。显然，从解的质量方面来看，即使调度规则支持快速决策，但是解的质量仍有很大的提升空间。
表4 NEH算法与调度规则相比的 API值
GA算法：表5给出了GA算法、NEH算法的比较评价，表中的值都是与同种组合的调度规则相比得到的API。根据表3、4的结果，我们可以看出不管是NEH算法还是调度规则，EDD-FIFO和STPT-FIFO两种组合表现的都比较好。所以在验证GA算法表现时，我们只展示这两种组合分别作为GA的初始解时GA解的情况。在表5中，符号GA-EDD-FIFO中，EDD-FIFO表示将NEH-EDD-FIFO算法得到的解作为GA算法初始种群中的一个，与其它算法一样，FIFO表示修理车间和组装车间的调度规则是FIFO。
从表5中可以看出，对于给出的算例，GA-EDD-FIFO和GA-STPT-FIFO两种组合得到的解与相同组合的调度规则得到的解相比，其平均API值分别为21.16%和21.18%，而对于NEH-EDD-FIFO和NEH-STPT-FIFO两种组合，其平均API值分别为18.07%和16.57%。所以，与NEH算法相比，有初始解的GA算法可以得到质量更高的解。相比有初始解的GA算法，纯GA算法的平均API仍然大于NEH算法，但是随着待维修件数量增大，纯GA算法的API值远小于有初始解的GA算法，可见，纯GA算法此时效率较低，改进的GA算法表现优于纯GA算法。
表5 GA、 NEH算法分别与调度规则相比的 API值
同样我们可以注意到，对于同种调度规则组合，NEH算法对解质量的提高百分比大于GA算法对NEH算法解质量的提高百分比，这可能是因为GA算法得到的解已经接近最优解或者是GA算法并没有比NEH算法好太多。然而从算法计算时间方面来看，GA算法是非常耗时的，表5中GA算法最大计算时间设为3600秒。因此，实际应用中，当对解的质量要求很高时，GA算法才有应用价值。
5 结论
本文针对MRO调度问题的复杂性和不确定性构建仿真模型，并提出三种启发式算法，分别是调度规则、NEH算法和GA算法，其中NEH算法属于构造型启发式算法，GA属于改进型启发式算法。然后对三种启发式算法的计算流程进行详细介绍，本文最后给出若干计算算例，分别使用三种启发式算法进行计算，并将计算结果进行展示和对比分析。
通过大量算例可以得出结论：调度规则算法计算时间快，但计算结果较差，GA算法的计算结果最好，但是非常耗时，适用于对解的质量要求很高的情况下，而NEH算法在计算时间和求解质量方面居于调度规则和GA之间。
[参考文献]
[1]郝春锋，丁金想，栾世超.基于MO2TOS算法的航空维修业调度问题研究[J].工业工程,2017,6(20):56-64.
[2]Reményi C,Staudacher S.Systematic simulation based approach for the identification and implementation of a sche-duling rule in the aircraft engine maintenance[J].International Journal of Production Economics,2014,147:94-107.
[3]Guide V D R,Kraus M E,Srivastava R.Scheduling policies for remanufacturing[J].International Journal of Production Economics,1997,48(2):187-204.
[4]Guide V D R,Srivastava R.An evaluation of order release strategies in a remanufacturing environment[J].Computers & operations research,1997,24(1):37-47.
[5]Guide V D R,Souza G C,van der Laan E.Performance of static priority rules for shared facilities in a remanufacturing shop with disassembly and reassembly[G].European Journal of Operational Research,2005,164(2):341-353.
[6]Kim M G,Yu J M,Lee D H.Scheduling algorithms for remanufacturing systems with parallel flow-shop-type reprocessing lines[J].International Journal of Production Research,2015,53(6):1819-1831.
[7]Nawaz M,Enscore E E,Ham I.A heuristic algorithm for the m-machine,n-job flow-shop sequencing problem[J].Omega,1983,11(1):91-95.
[8]Dong X,Huang H,Chen P.An improved NEH-based heuristic for the permutation flowshop problem[J].Compu-ters & Operations Research,2008,35(12):3962-3968.
[9]Fernandez-Viagas V,Framinan J M.NEH-based heuristics for the permutation flowshop scheduling problem to minimise total tardiness[J].Computers & Operations Research,2015,60:27-36.
[10]Allahverdi A,Aydilek H.The two stage assembly flowshop scheduling problem to minimize total tardiness[G].Journal of Intelligent Manufacturing,2014,26(2):225-237.
[11]Bertel S,Billaut J C.A genetic algorithm for an industrial multiprocessor flow shop scheduling problem with recirculation[G].European Journal of Operational Research,2004,159(3):651-662.
[12]Huang J,Süer G A,Urs S B.Genetic algorithm for rotary machine scheduling with dependent processing times[J].Journal of Intelligent Manufacturing,2012,23(5):1931-1948.
[13]王凌,郑大钟.求解同顺序加工调度问题的一种改进遗传算法[J].系统工程理论与实践,2002,22(6):74-79.
[14]Murata T,Ishibuchi H.Positive and negative combination effects of crossover and mutation operators in sequencing problems[C].Evolutionary Computation,1996.,Proceedings of IEEE International Conference on.IEEE,170-175.
[15]Potts C N,Van Wassenhove L.A decomposition algorithm for the single machine total tardiness problem[J].Opera-tions Research Letters,1982,1(5):177-181.
Research on MRO Scheduling Problem Based on Simulation Optimization
□DING Jin -xiang ，LUAN Shi -chao ，SHI Xiao -lei ，LIAN Fu -shi
(AVIC China Aero -Polytechnology Establishment ,Beijing 100028,China )
【Abstract 】MRO is the abbreviation of aviation maintenance,namely,maintenance,repair and overhaul.Compared to the traditional manufacturing system,MRO enterprises face more uncertainty and complexity,such as,Disassembly-Repair-Assembly three-level structure、Material recovery uncertainty、Material matching requirements、Stochastic routings and Variable processing times,these characteristics make the planning and scheduling of MRO system more difficult to deal with.On the basis of a detailed analysis of the MRO scheduling problem,a simulation model to minimize the expected total weighted tardiness is presented and three typical simulation optimization methods are applied to solve this problem,namely,dispatching rule,NEH-based heuristics and genetic algorithm.Finally,the characteristics of different algorithms are analyzed by the numerical examples based on the practical background.Some decision-making suggestions for the scheduling manager of MRO enterprise are proposed.
【Key words 】aviation maintenance;scheduling;stochastic routings;variable processing times;simulation optimization
【中图分类号】 F251;TN207
【文献标识码】 B
【文章编号】 1674-4993(2019)05-0137-05
doi: 10.3969/j.issn.1674-4993.2019.05.049
【收稿日期】 2019-03-21
【作者简介】 丁金想(1992— ),男,安徽阜阳人,管理学硕士,研究方向:精益生产、运营管理、系统优化。栾世超(1982— ),男,山东高密人,管理学博士,研究方向:工业工程、精益生产、系统优化。