1.一种柔性作业车间插单动态调度优化方法，其特征在于，建立柔性作业车间调度问题数学模型，采用实验仿真方式获得合理分批选择策略，同时根据对典型算例的仿真计算，给出分批数量的推荐值，其次基于工序、机器、分批后各子批的分单数量的三层基因染色体，采用粒子群算法与遗传算法结合的混合算法，形成改进的优化技术，进行问题优化求解；具体步骤如下：首先建立柔性作业车间多目标批量调度的数学模型，包括2个步骤：

第一、确定柔性作业车间调度问题描述和相关假设；

通过对柔性车间动态调度问题的研究，假设m台不同的机床集合{M1,M2,…,Mm}和n批待加工(含插单)工件集合{J1',J'2,…,J'n}，每批工件有不同的数量大小，并且包含多个工序，同时，每道工序都有可供选择的加工机床，且各工序加工时间随着机床不同而变化，为了减少拖延期，将某批次Je工件的工序分为z子批次，各个子批次在执行其每个工序时都视为不可分割的整体；根据现场实际情况，提出以下假设：

1)资源约束：每个资源同时只能服务一批工件；

2)顺序约束：每批工件只能由一台机器处理；

3)不同批次工件的加工顺序无限制；

第二、构建多目标批量调度的数学模型；

在以上假设条件下，并且根据生产要求，以最大完工时间最小和拖延期最短为优化目标，将子批订单和未分批的订单组合成新任务集合{J1,J2,…,Jn'}，在m台设备上进行合理安排；

其目标函数f为

minf＝xf1+(1-x)f2  (1)

f1＝max Ti  (2)

约束方程如公式(4)-(8)

当Xijk＝Xi(j-1)m＝1,k＝m  (4)

当Xijk＝Xi(j-1)m＝1,k≠m  (5)

当Xqgk＝Xijk＝1,Rijqgk＝1  (6)

当Xijk＝1  (7)

其中：：i的集合为i∈(1,n')，q的集合为q∈(1,n')，n'为工件分批后的总批次，且n'＝n+z-1，j∈(1,li)，g∈(1,lq)，li、lq分别为第i、q批工件的工序数；k,k'∈(1,m)，m为机床数量；Ji、Jq分别为第i、q批工件，n为工件分批前的总批次； 分别为第i批工件Ji、第q批工件Jq的第j、g道工序在设备k的上的单件加工时间； 分别为Ji、Jq工件第j、g道工序在设备k的上的批量开始时间；Qi、Qq分别为第i、q批工件的件数；Ti为工件Ji最后一道加工工序的完工时间；Ci为工件Ji加工工序的完工时间；Di为工件Ji最后一道加工工序的完工时间；Li为工件Ji的单位拖期；ωi为Li的惩罚系数或称惩罚因子；f1为所有批次的完工时间，f2为所有批次的总拖期时间；Sijk、Eijk分别表示工件Ji的工序在设备k上的开始时间和完成时间；Ei(j-1)m表示工件Ji的前一道工序的完成时间；Eqgk表示工件Jq的工序在设备k上的完成时间；x为优化目标函数权重；两个决策变量为：同理Xqgk为决策变量，即当Jq的工序g在设备k上加工时为1，否则为0；Xi(j-1)k'为决策变量，即当Ji的前一道工序在设备k'上加工时为1，否则为0。

式(3)、(4)表示任意加工批次必须在上一道工序完工后才能开始下一道工序；式(5)表示同一台设备不能同时加工两道工序；式(6)表示任意工序的完工时间与开始时间的差值不能小于其加工所需的时间；式(7)表示任意工序唯一由一台设备独立完成；

其次通过实验方式得到分批选择策略和数量，包括以下2个步骤：

第三、设计实验，获取分批选择策略

在三个合理的选择中选择最适合的选择策略，三个策略分别为：(a)选择单件加工时间最长的订单；(b)选择工件数量最多的订单；(c)选择加工时间最长的订单；通过实验仿真得出第三种方案最优；

第四、通过实验仿真方法即通过最优分单策略获得分批的数量

分别设置u+2、u+1、u三种不同的预设批数量，其中u为所有工序中能够同时选择机床的最大数，通过反复试验运算得到最优的分批数量应该小于等于u，且允许某一批数量为零；

最后，设置混合算法的初始化参数，产生离散组合问题的初始解种群，采用遗传算法进行优化求解，求解过程中，子批种群的更新方法采用粒子群算法更新原则，包括以下2个步骤：第五、设置以下初始化参数：种群数量、遗传代数、交叉率、变异率和学习因子，并以工序编码为基础，设计三层染色体，初始化种群；

初始参数按照经验选择，三层染色体则包含工序、机器、子批数量三层，同时保证子批数量总和为一定值，然后进行种群初始化，工序与机器基因通过随机方式产生，子批数量基因由于需要满足数值总和等于分批订单的工件数，则需要随机生成前几个批次的数量，最后一批数量采用总数减去前几批的和得到，同时若前几批批量和等于总数量，则剩下批次设置为零；

第六、利用遗传算法优化求解，并结合粒子群算法，改进子批数量每代的个体种群更新方式；

子批数量的更新，采用粒子群算法的方式，通过公式(9)、(10)更新分批数量粒子，使分批数量能够向着最优方向变异，并且通过控制速度矢量和为零以保证分批数量总和不变，同时加快了分批数量向最优解的收敛速度；

式中，d表示D维空间中的第d个变量，d＝1,2,…,D；vid和pid分别表示第i个粒子在第d维的速度和位置分量； 和 分别表示第i个粒子上一代第d维的速度和更新后的粒子速度，和 分别表示第i个粒子上一代第d维的位置和更新后的粒子位置；k和k+1为迭代次数；

Pid为第i个粒子迭代过程中最优适应度值对应的第d维位置分量；Pgd为全局最优适应度值对应的第d维位置分量；c1,c2为学习因子，为非负的加速度常量，其取值决定了粒子向Pid和Pgd的转移速度变化，通常都取2.0；r1为0-1的随机数，它与c1的乘积决定粒子向个体最优位置移动的概率；r2为0-1的随机数，它与c2的乘积决定粒子向全局最优位置移动的概率；w为非负的惯性因子；

其他基因采用传统的遗传选择、交叉、变异操作，循环更新个体，计算目标函数值，分配适应度函数值，直到满足迭代次数或是满足终止条件。

2.根据权利要求1所述的一种柔性作业车间插单动态调度优化方法，其特征在于：所述最优分单策略是通过程序化方式实现自动选择：首先选择批次工序加工时间最长的批次，如果存在相同的批次，则计算相同批次的所有工序加工时间，取加工时间长的批次，如果仍存在相同批次，则产生0-1的随机数，小于0.5选择前者，大于0.5选择后者。

3.根据权利要求1所述的一种柔性作业车间插单动态调度优化方法，其特征在于：所述三层染色体中的工序、机器、子批数量分别为3道工序3台机床和3批工件，其调度优化方法如下：首先假设要分单的工件是第3批工件，第3批工件的总数量为W,则分单的数量最大时等于同时选择的最大机床数量3，以保证分单后存在最佳情况，即各子批订单在不同机床上同时开工；三层基因染色体具体的对应关系如下，，工序1、2分别代表第一、第二批工件，

3.1、3.2、3.3分别代表第三批工件的第二、第三子批，加工机床为其在机器序列中对应位置的机床号，此外，工序序列中同一工件号出现的次数代表工件的加工工序，即第几次出现则为第几道工序；分单数量中A1A2A3A4分别代表不同的第三批的子批1,2,3,4的分单数量，同时各子批的数量和等于第三批工件总数，而且子批次中数量可以为零；

工序[2 1 3.1 1 3.2 3.3 1 2 3.1 3.2 3.2 3.3 3.1 3.3 2]机器[2 1 3 1 3 2 1 2 3 1 2 2 3 1 3]

子批数量[A1 A2 A3 A4]。

4.根据权利要求1所述的一种柔性作业车间插单动态调度优化方法，其特征在于：采用遗传算法随机产生初始种群，再与粒子群算法结合的方式，利用式(9)(10)更新子批数量的种群个体，以提高算法的收敛速度和稳定性。