CN108279647A - 一种两阶段装配流水车间调度方法 - Google Patents
一种两阶段装配流水车间调度方法 Download PDFInfo
- Publication number
- CN108279647A CN108279647A CN201711275686.5A CN201711275686A CN108279647A CN 108279647 A CN108279647 A CN 108279647A CN 201711275686 A CN201711275686 A CN 201711275686A CN 108279647 A CN108279647 A CN 108279647A
- Authority
- CN
- China
- Prior art keywords
- workpiece
- working position
- machine
- cell receptor
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012545 processing Methods 0.000 claims abstract description 55
- 238000000034 method Methods 0.000 claims abstract description 49
- 230000008569 process Effects 0.000 claims abstract description 35
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 27
- 210000000987 immune system Anatomy 0.000 claims abstract description 9
- 238000013459 approach Methods 0.000 claims abstract description 4
- 230000001360 synchronised effect Effects 0.000 claims description 88
- 108091008875 B cell receptors Proteins 0.000 claims description 75
- 108020003175 receptors Proteins 0.000 claims description 68
- 108091008874 T cell receptors Proteins 0.000 claims description 26
- 102000016266 T-Cell Antigen Receptors Human genes 0.000 claims description 26
- 210000003719 b-lymphocyte Anatomy 0.000 claims description 25
- 238000004519 manufacturing process Methods 0.000 claims description 19
- 210000004027 cell Anatomy 0.000 claims description 18
- 210000001744 T-lymphocyte Anatomy 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 17
- 230000003993 interaction Effects 0.000 claims description 14
- 239000000427 antigen Substances 0.000 claims description 8
- 102000036639 antigens Human genes 0.000 claims description 8
- 108091007433 antigens Proteins 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 8
- 230000035772 mutation Effects 0.000 claims description 8
- 108090000623 proteins and genes Proteins 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 6
- 230000006978 adaptation Effects 0.000 claims description 5
- 238000003780 insertion Methods 0.000 claims description 5
- 230000037431 insertion Effects 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 5
- 230000028327 secretion Effects 0.000 claims description 5
- 230000003466 anti-cipated effect Effects 0.000 claims description 3
- 230000008260 defense mechanism Effects 0.000 claims 1
- 230000008707 rearrangement Effects 0.000 claims 1
- 102000005962 receptors Human genes 0.000 description 55
- 238000011160 research Methods 0.000 description 6
- 241000728904 Iais Species 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 210000003563 lymphoid tissue Anatomy 0.000 description 3
- 238000003754 machining Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 102000006306 Antigen Receptors Human genes 0.000 description 1
- 108010083359 Antigen Receptors Proteins 0.000 description 1
- 108060003951 Immunoglobulin Proteins 0.000 description 1
- 230000033289 adaptive immune response Effects 0.000 description 1
- 230000004721 adaptive immunity Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005194 fractionation Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 102000018358 immunoglobulin Human genes 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 210000005210 lymphoid organ Anatomy 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 244000052769 pathogen Species 0.000 description 1
- 230000001717 pathogenic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000001541 thymus gland Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41865—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by job scheduling, process planning, material flow
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32252—Scheduling production, machining, job shop
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- General Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
Abstract
本发明属于车间调度技术领域，具体的说是一种两阶段装配流水车间调度方法。该方法包括以下步骤：步骤一、建立不同交货期窗口约束下的两阶段装配流水车间调度的目标函数；步骤二、确定调度过程的约束条件；步骤三、采用基于细胞受体的人工免疫系统算法求解所述目标函数，经过迭代运算，将满足停止条件的解解码后作为最终的调度方案。本发明是一种基于细胞受体的人工免疫系统算法获得工作加工最优调度的两阶段装配流水车间调度方法，解决了现有技术在求解不同交货期窗口约束下的两阶段装配流水车间调度问题时无法求解大规模问题、收敛精度低、收敛速度慢的问题。
Description
技术领域
本发明属于车间调度技术领域，具体的说是一种两阶段装配流水车间调度方法。
背景技术
两阶段装配流水车间调度问题(Two-stage Assembly Flowshop SchedulingProblem,TSAFSP)是一种车间调度问题，广泛存在于各种生产过程中，如计算机制造和消防车生产等。经典TSAFSP包括制造与装配两个阶段，其中制造阶段有m台并行同速机，装配阶段只有一台装配机。一个工件的加工工艺包含m+1项操作，其中m项操作在制造阶段进行，而后经过装配操作将零部件组装成最终工件。当制造阶段机器的数量为2时，该问题就被证实为NP-hard问题，所以TSAFSP是一个NP-hard问题。经典TSAFSP调度问题对实际问题做了较大的简化，然而实际生产环境要复杂的多，出于理论研究和工程实际的需要，研究者提出了各种TSAFSP调度问题模型进行研究。首先，制造阶段每项操作加工处理的零件可以是多个同类型的零件，如果将每项操作理解为一个批次的话，每项操作可以拆分为多个小的批次，并在不同的并行同速机上进行加工，这更贴近实际生产情况。其次，由于零件类型的不同及拆分的处理，准备时间应当从加工时间中分离出来，在每项操作前都会有相应的准备时间。另外，之前大多数的研究成果在考虑工件交货期约束时，将其视为一个时间点。然而，在制造业生产中，工作更倾向于被要求在一个时间窗口内完成。工作早于这个时间窗口的起始时间完工的话，会产生库存成本；而晚于这个时间窗口结束时间完工的话，会产生延迟费。如何综合考虑上述实际生产因素，科学地制定车间调度方案，对于降低生产管理成本、提高机器利用率、和满足客户个性化的收货时间窗需求起着至关重要的作用。
两阶段装配流水车间调度方法的研究最初主要应用动态规划、分支定界等精确算法，但这些方法通常计算时间较长，不适用于求解大规模问题。今年来，随着计算智能方法的发展，越来越多的智能优化算法，如变领域搜索算法(International Journal ofProduction Research,2014,52(19):5626-5639)、粒子群优化算法(Computers&Operations Research,2009,36(9):2682-2689)、模拟退火算法(Advances inEngineering Software,2010,41(10):1238-1243)被应用于求解两阶段装配流水车间调度问题。但这些算法由于收敛能力不足，对于上述更接近实际环境的两阶段装配流水车间问题，不能在规定的时间内求得满意解。
发明内容
本发明提供了一种基于细胞受体的人工免疫系统算法(RAIS)获得工作加工最优调度的两阶段装配流水车间调度方法，解决了现有技术在求解不同交货期窗口约束下的两阶段装配流水车间调度问题时无法求解大规模问题、收敛精度低、收敛速度慢的问题。
本发明技术方案结合附图说明如下：
一种两阶段装配流水车间调度方法，该方法包括以下步骤：
步骤一、建立不同交货期窗口约束下的两阶段装配流水车间调度的目标函数；
步骤二、确定调度过程的约束条件；
步骤三、采用基于细胞受体的人工免疫系统算法即RAIS算法求解所述目标函数，经过迭代运算，将满足停止条件的解解码后作为最终的调度方案。
步骤一中所述的目标函数如下：
i为工件号，n为工件个数；
αi为工件i的单位提前惩罚成本，βi为工件i的单位滞后惩罚成本；
Ei为工件i的提前时间，Ti为工件i的滞后时间；
Z为提前与滞后惩罚成本加权和。
步骤二中所述的约束条件包括：
零件在并行同速机上不同位置约束：
零件拆分处理约束：
M×xjtk≥qjtk j,t＝1,2,...,h,k＝1,2,...m (6)
qjtk-1≥M×(xjtk-1) j,t＝1,2,...,h,k＝1,2,...m (7)
同一并行同速机不同位置上的零件完工时间约束：
C0k＝0 k＝1,2,...,m (8)
工件在装配机不同位置约束：
装配机不同位置上的工件完工时间约束：
Cr≥Fi+Ai-M×(1-yir) i,r＝1,2,...,n (13)
工件提前成本与滞后成本约束：
Ei≥ei-Cr-M×(1-yir) i,r＝1,2,...,n (15)
Ti≥Cr-di-M×(1-yir) i,r＝1,2,...,n (16)
Ei≥0 i＝1,2,...,n (17)
Ti≥0 i＝1,2,...,n (18)
上述公式(2)～(18)中，存在以下决策标量：
r为装配机上工件加工的位置号；
m为并行同速机个数；
k为并行同速机号；
h为总的零件种类数；
j为零部件种类号；
t为并行同速机上零件加工的位置号；
Ji为工件i的零件种类集；
s1j为加工零件j前的准备时间；
q1j为零件j的总个数；
p1j为零件j的单个加工时间；
s2i为加工工件i前的准备时间；
Ai为工件i的装配时间；
ei为工件i交货期窗口的起始时间；
di为工件i交货期窗口的结束时间；
Ctk为并行同速机k位置t上的零件的完工时间；
C(t-1)k为并行同速机k位置t-1上的零件的完工时间；
C0k为并行同速机k位置0上的完工时间即机器未进行加工状态下的时间点；
Fi为工件i进行装配操作的最早开始时间；
Cr为装配机位置r上的工件的完工时间；
Cr-1为装配机位置r-1上的工件的完工时间；
qjtk为并行同速机k位置t上的零件j的个数；
M为一个足够大的正整数；
在上述公式(2)～(18)中，存在以下决策变量：
所述的步骤三的具体方法如下：
用RAIS求解不同交货期窗口约束下的TSAFSP包括编码、解码和适应度值计算的操作；
31)编码时，每个受体采用基于工件序列的编码，记作(π1,π2,...,πn)，即加工位置i处的工件为πi；
32)解码时，每个受体的解码步骤如下：
步骤3.2.1：依据公式(2)～(10)，采用拆分法确定每台并行同速机上零件种类、数量和完工时间；
a.设定i＝1,转b；
b.计算其中CT0k＝0，如果i≤n，则计算并对TPj按从大到小排序得转c；如果i＝n+1，则停止；
c.设定t＝1，MCTk＝0(k＝1,2,...,m)转d；
d.令w为min(MCT1,MCT2,...,MCTm)的机器号，如果q1t′＝0，则令t＝t+1转d；如果同时满足q1t′＞0，则q′tw＝q′1t，q1t′＝0，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw，令t＝t+1转d；如果同时满足 q1t′＞0，则q1t′＝q1t′-qtw′，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw，转d；如果同时满足q1t′＞0，则qtw′＝q1t′，q1t′＝0，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw令t＝t+1转d；如果且q1t′＝0，则i＝i+1转b；
步骤3.2.1中存在以下决策标量：
i为工件加工位置号，j为零件号，k为并行同速机器号；
为加工位置号为i的工件的零件在制造阶段的完工时间下界；
CTik为加工位置号为i的工件的零件在并行同速机k上的最大完工时间；
CTiw为加工位置号为i的工件的零件在并行同速机w上的最大完工时间；
CT(i-1)k为加工位置号为i-1的工件的零件在并行同速机k上的最大完工时间；
为装配工件πi所需的零件种类数；
为工件πi的零件j的不拆分总处理时间；
为装配工件πi所需的种零件中，第t大的不拆分总处理时间；
MCTk为并行同速机k的完工时间；
MCTw为并行同速机w的完工时间；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的总个数；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件在并行同速机w上的个数；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的准备时间；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的单个加工时间；
为对数a向上取整；
步骤3.2.2：依据公式(11)～(14)，采用空闲时间插入法确定装配机上工件的完工时间：
a.设定i＝1，ACT＝0，转b；
b.判断i≤n，如果不满足，设定i＝n，ξr＝0(r＝1,2,...,n)，l＝1，FJl＝LJl＝i转c；否则判断满足则ACT＝Ci，不满足则ACT＝Ci；令i＝i+1转b；
c.判断i≥1，如果不满足，则停止；否则判断i＝n，如果满足，则ξi＝1，否则判断如果满足，则ξi＝ξi+1，FJl＝i，否则ξi＝ξi+1+1，l＝l+1，FJl＝LJl＝i转d；
d.令σl＝{r|ξr＝ξi,r＝1,2,...,n}，若则若则若则ηr＝0；计算判断如果满足，则i＝i-1转c，否则转e；
e.若则若则Δ1r＝M；若则判断ξi＝1，如果满足，则Δ2l＝M，如果不满足，转f；
f.若则Cr＝Cr+Δ2l，ξr＝ξr-1，FJl-1＝FJl，l＝l-1转c；若则转c；
CTi1为加工位置号为i的工件的零件在并行同速机1上的最大完工时间；
CTi2为加工位置号为i的工件的零件在并行同速机2上的最大完工时间；
CTim为加工位置号为i的工件的零件在并行同速机m上的最大完工时间；
Ci为装配机位置i上的工件的完工时间；
Ci+1为装配机位置i+1上的工件的完工时间；
为加工工件πi前的准备时间；
为加工工件πi+1前的准备时间；
为加工工件前的准备时间；
为工件πi的装配时间；
为工件πi+1的装配时间；
为工件的装配时间；
ACT为装配机的完工时间；
ξi为装配机上加工位置号为i的工件的标记号；
ξi+1为装配机上加工位置号为i+1的工件的标记号；
ξr为装配机上加工位置号为r的工件的标记号；
l为集合σl的标记号；
FJl为标记号为l的工件集中最早进行装配的工件的加工位置号；
FJl-1为标记号为l-1的工件集中最早进行装配的工件的加工位置号；
LJl为标记号为l的工件集中最晚进行装配的工件的加工位置号；
为装配机位置LJl上的工件的完工时间；
为装配机位置LJl-1上的工件的完工时间；
为工件πr的单位提前惩罚成本；
为工件πr的单位滞后惩罚成本；
ηi为装配机上加工位置号为i的工件的完工时间增加一个单位造成的成本减少值；
ηr为装配机上加工位置号为r的工件的完工时间增加一个单位造成的成本减少值；
为工件πr交货期窗口的起始时间；
为工件πr交货期窗口的结束时间；
σl为标记号为l的工件集；
Δ1r为工件在不改变其完工状态(提前、滞后或如期完工)时，完工时间允许的最大增加值；
Δ2l为集合σl中的工件在不影响其他工件集中的工作时，完工时间允许的最大增加值。
步骤3.2.3：依据公式(1)和(15)～(18)计算目标函数；
33)由于优化目标为提前和滞后惩罚成本加权和最小，采用目标函数的倒数作为适应度函数对每个受体进行评价；
34)RAIS算法求解不同交货期窗口约束下的TSAFSP的步骤如下：
步骤3.4.1：随机初始化NB个B细胞受体和NT个T细胞受体种群，转步骤3.4.2；
步骤3.4.2：采用选择和插入的方式对B细胞受体和T细胞受体种群进行基因再排列：首先将当前适应度值最高的受体作为标准受体，然后从其他受体的工件序列中随机选择R个加工位置号，将这R个加工位置号上对应的工件分别插入到其在标准受体中的加工位置上，若这个新的工件序列解码的适应度值大于被选择的受体适应度值，则更新被选择的受体，否则被选择的受体保持不变，转步骤3.4.3；
步骤3.4.3：比对B细胞受体与当前适应度值最高的T细胞受体在相同加工位置处的工件，如果存在工件相同，则B细胞与T细胞发生同源相互作用，转步骤4.4.4，否则B细胞与T细胞不发生同源相互作用，转步骤3.4.6；
步骤3.4.4：采用部分随机排序的方式对B细胞受体进行细胞突变：对比B细胞受体与当前适应度值最高的T细胞受体在相同加工位置处的工件，工件相同情况下组成的加工位置号集合为工件不相同情况下组成的加工位置号集合为δ，集合中加工位置号上工件保持不变，对集合δ中加工位置号上的工作重新随机排序生成新的集合δ′，若解码的适应度值大于被选择的B细胞受体适应度值，则更新被选择的B细胞受体，否则被选择的B细胞受体保持不变，转步骤3.4.5；
步骤3.4.5：采用交换、插入、和颠倒的领域搜索机制对B细胞受体进行抗体转换：由于B细胞在T细胞受体的帮助下可以大量分化为分泌IgG、IgA、和IgE形态抗体的浆细胞进而消灭抗原，所以对被选择的B细胞受体，分别以1/3的概率进行交换、插入、或颠倒搜索进而达到局部最优，每个被选择的B细胞受体共进行I次抗体转换，若抗体转换后B细胞受体的适应度值提升，则更新被选择的B细胞受体，否则不变，转步骤3.4.7；
步骤3.4.6：采用插入和交换的领域搜索机制对B细胞受体进行IgM抗体转换：由于B细胞在不与T细胞发生同源相互作用时，只是表达IgM抗体，提供一个快速的早期防御机制，所以被选择的B细胞受体先进行插入，再做交换，若IgM抗体的表达后B细胞受体的适应度值提升，则更新被选择的B细胞受体，否则不变，转步骤3.4.7；
步骤3.4.7：当前适应度值最高的受体保持不变，其他受体重新随机生成，转3.4.8；
步骤3.4.8：判断是否满足预先设定的最长CPU运行时间Tmax，如果满足，将当前最好的解作为最优解，得到最终的调度方案并停止，如果不满足转步骤3.4.2。
本发明的有益效果为：
1.本发明采用拆分法对制造阶段零部件进行加工安排，可以提高解的质量；
2.本发明采用空闲时间插入法对装配阶段工件进行安排，可以提高解的质量；
3.本发明引入T细胞对B细胞的帮助过程，通过判断它们之间是否发生同源相互作用，将B细胞的增值分化分为两个方向：1.先进行细胞突变，再抗体转换；2.只表达IgM抗体，在保证解的收敛的速度的同时提高了收敛精度。
附图说明
图1为本发明中基于细胞受体的人工免疫系统算法的实现流程示意图；
图2为4工件基因再排列示意图；
图3为4工件细胞突变示意图；
图4为4工件9零件实例最优调度甘特图。
具体实施方式
参阅图1，一种两阶段装配流水车间调度方法，该方法包括以下步骤：
步骤一、建立不同交货期窗口约束下的两阶段装配流水车间调度的目标函数；
步骤一中所述的目标函数如下：
i为工件号，n为工件个数；
αi为工件i的单位提前惩罚成本，βi为工件i的单位滞后惩罚成本；
Ei为工件i的提前时间，Ti为工件i的滞后时间；
Z为提前与滞后惩罚成本加权和。
步骤二、确定调度过程的约束条件；
步骤二中所述的约束条件包括：
零件在并行同速机上不同位置约束：
零件拆分处理约束：
M×xjtk≥qjtk j,t＝1,2,...,h,k＝1,2,...m (6)
qjtk-1≥M×(xjtk-1) j,t＝1,2,...,h,k＝1,2,...m (7)
同一并行同速机不同位置上的零件完工时间约束：
C0k＝0 k＝1,2,...,m (8)
工件在装配机不同位置约束：
装配机不同位置上的工件完工时间约束：
Cr≥Fi+Ai-M×(1-yir) i,r＝1,2,...,n (13)
工件提前成本与滞后成本约束：
Ei≥ei-Cr-M×(1-yir) i,r＝1,2,...,n (15)
Ti≥Cr-di-M×(1-yir) i,r＝1,2,...,n (16)
Ei≥0 i＝1,2,...,n (17)
Ti≥0 i＝1,2,...,n (18)
上述公式(2)～(18)中，存在以下决策标量：
r为装配机上工件加工的位置号；
m为并行同速机个数；
k为并行同速机号；
h为总的零件种类数；
j为零部件种类号；
t为并行同速机上零件加工的位置号；
Ji为工件i的零件种类集；
s1j为加工零件j前的准备时间；
q1j为零件j的总个数；
p1j为零件j的单个加工时间；
s2i为加工工件i前的准备时间；
Ai为工件i的装配时间；
ei为工件i交货期窗口的起始时间；
di为工件i交货期窗口的结束时间；
Ctk为并行同速机k位置t上的零件的完工时间；
C(t-1)k为并行同速机k位置t-1上的零件的完工时间；
C0k为并行同速机k位置0上的完工时间即机器未进行加工状态下的时间点；
Fi为工件i进行装配操作的最早开始时间；
Cr为装配机位置r上的工件的完工时间；
Cr-1为装配机位置r-1上的工件的完工时间；
qjtk为并行同速机k位置t上的零件j的个数；
M为一个足够大的正整数；
在上述公式(2)～(18)中，存在以下决策变量：
步骤三、采用基于细胞受体的人工免疫系统算法即RAIS算法求解所述目标函数，经过迭代运算，将满足停止条件的解解码后作为最终的调度方案；
基于细胞受体的人工免疫系统算法的设计来源于自然免疫系统中的适应性免疫应答过程。在适应性免疫应答过程中，B细胞与T细胞共同发现并消灭侵入人体的病原菌，其中B细胞受体和T细胞受体起到关键性作用，它们的发展过程如下：
a：B细胞和T细胞在最初淋巴器官内经过基因再排列产生相应的抗原受体，而后B细胞会离开骨髓进入周边淋巴组织，T细胞会离开胸腺进入第二淋巴组织；
b：B细胞利用其表面受体蛋白质抓住胞外抗原后进入附近第二淋巴组织，在那里与T细胞受体接触并比对它们辨识的抗原构造；
c：如果B细胞受体与T细胞受体辨识的抗原相同，即它们发生同源相互作用，则B细胞会在T细胞的诱导下发生细胞突变，带有高亲和力受体的B细胞存活下来并分化为分泌抗体的浆细胞，B细胞受体从而以可溶型蛋白质-抗体的形式存在，包括IgG、IgA、和IgE共三种形态的抗体；如果B细胞受体与T细胞受体辨识的抗原不同，则B细胞只表达IgM抗体；
d：清除亲和力低的B细胞和T细胞。
所述的步骤三的具体方法如下：
用RAIS求解不同交货期窗口约束下的TSAFSP包括编码、解码和适应度值计算的操作；
31)编码时，每个受体采用基于工件序列的编码，记作(π1,π2,...,πn)，即加工位置i处的工件为πi；
32)解码时，每个受体的解码步骤如下：
步骤3.2.1：依据公式(2)～(10)，采用拆分法确定每台并行同速机上零件种类、数量和完工时间；
a.设定i＝1,转b；
b.计算其中CT0k＝0，如果i≤n，则计算并对TPj按从大到小排序得转c；如果i＝n+1，则停止；
c.设定t＝1，MCTk＝0(k＝1,2,...,m)转d；
d.令w为min(MCT1,MCT2,...,MCTm)的机器号，如果q1t′＝0，则令t＝t+1转d；如果同时满足q1t′＞0，则q′tw＝q′1t，q1t′＝0，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw，令t＝t+1转d；如果同时满足 q1t′＞0，则q1t′＝q1t′-qtw′，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw，转d；如果同时满足q1t′＞0，则qtw′＝q1t′，q1t′＝0，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw令t＝t+1转d；如果且q1t′＝0，则i＝i+1转b；
步骤3.2.1中存在以下决策标量：
i为工件加工位置号，j为零件号，k为并行同速机器号；
为加工位置号为i的工件的零件在制造阶段的完工时间下界；
CTik为加工位置号为i的工件的零件在并行同速机k上的最大完工时间；
CTiw为加工位置号为i的工件的零件在并行同速机w上的最大完工时间；
CT(i-1)k为加工位置号为i-1的工件的零件在并行同速机k上的最大完工时间；
为装配工件πi所需的零件种类数；
为工件πi的零件j的不拆分总处理时间；
为装配工件πi所需的种零件中，第t大的不拆分总处理时间；
MCTk为并行同速机k的完工时间；
MCTw为并行同速机w的完工时间；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的总个数；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件在并行同速机w上的个数；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的准备时间；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的单个加工时间；
为对数a向上取整；
步骤3.2.2：依据公式(11)～(14)，采用空闲时间插入法确定装配机上工件的完工时间：
a.设定i＝1，ACT＝0，转b；
b.判断i≤n，如果不满足，设定i＝n，ξr＝0(r＝1,2,...,n)，l＝1，FJl＝LJl＝i转c；否则判断满足则ACT＝Ci，不满足则ACT＝Ci；令i＝i+1转b；
c.判断i≥1，如果不满足，则停止；否则判断i＝n，如果满足，则ξi＝1，否则判断如果满足，则ξi＝ξi+1，FJl＝i，否则ξi＝ξi+1+1，l＝l+1，FJl＝LJl＝i转d；
d.令σl＝{r|ξr＝ξi,r＝1,2,...,n}，若则若则若则ηr＝0；计算判断如果满足，则i＝i-1转c，否则转e；
e.若则若则Δ1r＝M；若则判断ξi＝1，如果满足，则Δ2l＝M，如果不满足，转f；
f.若则Cr＝Cr+Δ2l，ξr＝ξr-1，FJl-1＝FJl，l＝l-1转c；若则转c；
步骤3.2.2中存在以下决策标量：
CTi1为加工位置号为i的工件的零件在并行同速机1上的最大完工时间；
CTi2为加工位置号为i的工件的零件在并行同速机2上的最大完工时间；
CTim为加工位置号为i的工件的零件在并行同速机m上的最大完工时间；
Ci为装配机位置i上的工件的完工时间；
Ci+1为装配机位置i+1上的工件的完工时间；
为加工工件πi前的准备时间；
为加工工件πi+1前的准备时间；
为加工工件前的准备时间；
为工件πi的装配时间；
为工件πi+1的装配时间；
为工件的装配时间；
ACT为装配机的完工时间；
ξi为装配机上加工位置号为i的工件的标记号；
ξi+1为装配机上加工位置号为i+1的工件的标记号；
ξr为装配机上加工位置号为r的工件的标记号；
l为集合σl的标记号；
FJl为标记号为l的工件集中最早进行装配的工件的加工位置号；
FJl-1为标记号为l-1的工件集中最早进行装配的工件的加工位置号；
LJl为标记号为l的工件集中最晚进行装配的工件的加工位置号；
为装配机位置LJl上的工件的完工时间；
为装配机位置LJl-1上的工件的完工时间；
为工件πr的单位提前惩罚成本；
为工件πr的单位滞后惩罚成本；
ηi为装配机上加工位置号为i的工件的完工时间增加一个单位造成的成本减少值；
ηr为装配机上加工位置号为r的工件的完工时间增加一个单位造成的成本减少值；
为工件πr交货期窗口的起始时间；
为工件πr交货期窗口的结束时间；
σl为标记号为l的工件集；
Δ1r为工件在不改变其完工状态(提前、滞后或如期完工)时，完工时间允许的最大增加值；
Δ2l为集合σl中的工件在不影响其他工件集中的工作时，完工时间允许的最大增加值。
步骤3.2.3：依据公式(1)和(15)～(18)计算目标函数；
33)由于优化目标为提前和滞后惩罚成本加权和最小，采用目标函数的倒数作为适应度函数对每个受体进行评价；
34)RAIS算法求解不同交货期窗口约束下的TSAFSP的步骤如下：
步骤3.4.1：随机初始化NB个B细胞受体和NT个T细胞受体种群，转步骤3.4.2；
步骤3.4.2：采用选择和插入的方式对B细胞受体和T细胞受体种群进行基因再排列：首先将当前适应度值最高的受体作为标准受体，然后从其他受体的工件序列中随机选择R个加工位置号，将这R个加工位置号上对应的工件分别插入到其在标准受体中的加工位置上，若这个新的工件序列解码的适应度值大于被选择的受体适应度值，则更新被选择的受体，否则被选择的受体保持不变，转步骤3.4.3；
步骤3.4.3：比对B细胞受体与当前适应度值最高的T细胞受体在相同加工位置处的工件，如果存在工件相同，则B细胞与T细胞发生同源相互作用，转步骤4.4.4，否则B细胞与T细胞不发生同源相互作用，转步骤3.4.6；
步骤3.4.4：采用部分随机排序的方式对B细胞受体进行细胞突变：对比B细胞受体与当前适应度值最高的T细胞受体在相同加工位置处的工件，工件相同情况下组成的加工位置号集合为工件不相同情况下组成的加工位置号集合为δ，集合中加工位置号上工件保持不变，对集合δ中加工位置号上的工作重新随机排序生成新的集合δ′，若解码的适应度值大于被选择的B细胞受体适应度值，则更新被选择的B细胞受体，否则被选择的B细胞受体保持不变，转步骤3.4.5；
步骤3.4.5：采用交换、插入、和颠倒的领域搜索机制对B细胞受体进行抗体转换：由于B细胞在T细胞受体的帮助下可以大量分化为分泌IgG、IgA、和IgE形态抗体的浆细胞进而消灭抗原，所以对被选择的B细胞受体，分别以1/3的概率进行交换、插入、或颠倒搜索进而达到局部最优，每个被选择的B细胞受体共进行I次抗体转换，若抗体转换后B细胞受体的适应度值提升，则更新被选择的B细胞受体，否则不变，转步骤3.4.7；
步骤3.4.6：采用插入和交换的领域搜索机制对B细胞受体进行IgM抗体转换：由于B细胞在不与T细胞发生同源相互作用时，只是表达IgM抗体，提供一个快速的早期防御机制，所以被选择的B细胞受体先进行插入，再做交换，若IgM抗体的表达后B细胞受体的适应度值提升，则更新被选择的B细胞受体，否则不变，转步骤3.4.7；
步骤3.4.7：当前适应度值最高的受体保持不变，其他受体重新随机生成，转3.4.8；
步骤3.4.8：判断是否满足预先设定的最长CPU运行时间Tmax，如果满足，将当前最好的解作为最优解，得到最终的调度方案并停止，如果不满足转步骤3.4.2。
实施例
本实施例是将RAIS算法应用到两阶段装配流水车间调度中，该车间生产4种工件，分别为工件A,B,C,D,其中装配工件A需要零件1,2，工件B需要零件3,4，工件C需要零件5,6,7，工件D需要零件8,9，工件的准备时间s2i、加工时间Ai、单位提前惩罚成本αi、单位滞后惩罚成本βi、交货期窗口的起始时间ei、和交货期窗口的结束时间di，零件的准备时间s1j、总个数q1j、和单个加工时间p1j等加工参数如表1所示，且该车间有3台并行同速机M1,M2,M3，一台装配机A。
表1工件与零件加工相关参数表
RAIS算法的解决不同交货期窗口约束下的两阶段装配流水车间调度方法，具体包括以下步骤：
(1)：对两阶段装配流水车间调度问题进行数学描述，并确立评价指标：提前与滞后惩罚成本加权和Z；
定义如下决策标量：
i为工件号；
r为装配机上工件加工的位置号；
n为工件个数；
j为零件种类号；
t为并行同速机上零件加工的位置号；
h为总的零件种类数；
Ji为工件i的零件种类集；
k为并行同速机号；
m为并行同速机个数；
s1j为零件j零件的准备时间；
q1j为零件j的总个数；
p1j为零件j的单个加工时间；
s2i为加工工件i的准备时间；
Ai为工件i的装配时间；
ei为工件i交货期窗口的起始时间；
di为工件i交货期窗口的结束时间；
αi为工件i的单位提前惩罚成本；
βi为工件i的单位滞后惩罚成本；
Ei为工件i的提前时间；
Ti为工件i的滞后时间；
Ctk为并行同速机k位置t上的零件的完工时间；
C(t-1)k为并行同速机k位置t-1上的零件的完工时间；
C0k为并行同速机k位置0上的完工时间(即机器未进行加工状态下的时间点)；
Fi为工件i进行装配操作的最早开始时间；
Cr为装配机位置r上的工件的完工时间；
Cr-1为装配机位置r-1上的工件的完工时间；
qjtk为并行同速机k位置t上的零件j的个数；
M为一个足够大的正整数；
定义如下决策变量：
(2)：建立目标函数；
(3)：建立约束条件；
M×xjtk≥qjtk j,t＝1,2,...,h,k＝1,2,...m (6)
qjtk-1≥M×(xjtk-1) j,t＝1,2,...,h,k＝1,2,...m (7)
C0k＝0 k＝1,2,...,m (8)
Cr≥Fi+Ai-M×(1-yir)i,r＝1,2,...,n (11)
Ei≥ei-Cr-M×(1-yir) i,r＝1,2,...,n (15)
Ti≥Cr-di-M×(1-yir) i,r＝1,2,...,n (16)
Ei≥0 i＝1,2,...,n (17)
Ti≥0 i＝1,2,...,n (18)
上述表达式中：式(2)表式一台并行同速机上的一个加工位置上的零件种类小于等于1；式(3)表示一种类型的零件在一台并行同速机上的加工位置个数小于等于1；式(4)表示如果同一台并行同速机某一位置上没有安排零件加工，其后置位上也不会有零件加工；式(5)表示一种类型零件的总个数等于所有并行同速机上该类型零件个数之和；式(6)和式(7)表示如果零件j在并行同速机k位置t上加工，则qjtk＞0，否则qjtk＝0；式(8)表示并行同速机位置0上的完工时间为0；式(9)表示同一台平行同速机上后一个加工位置上的零件的完工时间大于等于该零件的准备时间、该零件单个加工时间与该零件在该机器上加工个数的乘积、和其前一个加工位置上的零件的完工时间之和；式(10)和式(11)表示工件只有在其所需的零件都加工完后才能开始装配操作；式(12)表示装配机上同一加工位置上只能有一个工件；式(13)表示一个工件在装配机上只能对应一个加工位置；式(14)表示装配机上后一个加工位置上的工件的完工时间大于等于该工件的准备时间、装配时间、和其后一个加工位置上的工件的完工时间之和；式(15)～(18)表示每个工件的提前惩罚成本与滞后惩罚成本。
(4)：针对上述两阶段装配流水车间调度问题的目标函数和约束条件，设计RAIS算法得到调度方案；
(4.1)编码时，每个受体采用基于工件序列的编码，如(A,B,C,D)表示工件的加工次序是A,B,C,D，即工件A在1号位置，B在2号位置，C在3号位置，D在4号位置；
(4.2)解码时，每个受体的解码步骤如下：
步骤4.2.1：依据公式(2)～(10)，采用拆分法确定每台并行同速机上零件种类、数量和完工时间；
a.设定i＝1,转b；
b.计算其中CT0k＝0。如果i≤n，则计算并对TPj按从大到小排序得转c；如果i＝n+1，则停止；
c.设定t＝1，MCTk＝0(k＝1,2,...,m)转d；
d.令w为min(MCT1,MCT2,...,MCTm)的机器号。如果q1t′＝0，则令t＝t+1转d；如果同时满足q1t′＞0，则q′tw＝q′1t，q1t′＝0，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw，令t＝t+1转d；如果同时满足 q1t′＞0，则q1t′＝q1t′-qtw′，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw，转d；如果同时满足q1t′＞0，则qtw′＝q1t′，q1t′＝0，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw令t＝t+1转d；如果且q1t′＝0，则i＝i+1转b；
步骤4.2.1中存在以下决策标量：
i为工件加工位置号，j为零件号，k为并行同速机器号；
为加工位置号为i的工件的零件在制造阶段的完工时间下界；
CTik为加工位置号为i的工件的零件在并行同速机k上的最大完工时间；
CTiw为加工位置号为i的工件的零件在并行同速机w上的最大完工时间；
CT(i-1)k为加工位置号为i-1的工件的零件在并行同速机k上的最大完工时间；
为装配工件πi所需的零件种类数；
为工件πi的零件j的不拆分总处理时间；
为装配工件πi所需的种零件中，第t大的不拆分总处理时间；
MCTk为并行同速机k的完工时间；
MCTw为并行同速机w的完工时间；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的总个数；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件在并行同速机w上的个数；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的准备时间；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的单个加工时间；
为对数a向上取整。
步骤4.2.2：依据公式(11)～(14)，采用空闲时间插入法确定装配机上工件的完工时间：
a.设定i＝1，ACT＝0，转b；
b.判断i≤n，如果不满足，设定i＝n，ξr＝0(r＝1,2,...,n)，l＝1，FJl＝LJl＝i转c；否则判断满足则ACT＝Ci，不满足则ACT＝Ci；令i＝i+1转b；
c.判断i≥1，如果不满足，则停止；否则判断i＝n，如果满足，则ξi＝1，否则判断如果满足，则ξi＝ξi+1，FJl＝i，否则ξi＝ξi+1+1，l＝l+1，FJl＝LJl＝i转d；
d.令σl＝{r|ξr＝ξi,r＝1,2,...,n}，若则若则若则ηr＝0；计算判断如果满足，则i＝i-1转c，否则转e；
e.若则若则Δ1r＝M；若则判断ξi＝1，如果满足，则Δ2l＝M，如果不满足，转f；
f.若则Cr＝Cr+Δ2l，ξr＝ξr-1，FJl-1＝FJl，l＝l-1转c；若则转c；
步骤4.2.2中存在以下决策标量：
ACT为装配机的完工时间；
ξr为装配机上加工位置号为r的工件的标记号；
l为集合σl的标记号；
FJl为标记号为l的工件集中最早进行装配的工件的加工位置号；
LJl为标记号为l的工件集中最晚进行装配的工件的加工位置号；
ηi为装配机上加工位置号为i的工件的完工时间增加一个单位造成的成本减少值；
σl为标记号为l的工件集；
Δ1r为工件在不改变其完工状态(提前、滞后或如期完工)时，完工时间允许的最大增加值；
Δ2l为集合σl中的工件在不影响其他工件集中的工作时，完工时间允许的最大增加值。
(4.3)由于优化目标为提前和滞后惩罚成本加权和最小，采用目标函数的倒数作为适应度函数对每个受体进行评价；
(4.4)如图1，RAIS算法求解不同交货期窗口约束下的TSAFSP的步骤如下：
步骤4.4.1：设置参数NB,NT,R,I,随机初始化NB个B细胞受体和NT个T细胞受体种群，计算受体群的适应度值，转步骤4.4.2；
步骤4.4.2：采用选择和插入的方式对B细胞受体和T细胞受体种群进行基因再排列：首先将当前适应度值最高的受体作为标准受体，然后从其他受体的工件序列中随机选择R个加工位置号，将这R个加工位置号上对应的工件分别插入到其在标准受体中的加工位置上，若这个新的工件序列解码的适应度值大于被选择的受体适应度值，则更新被选择的受体，否则被选择的受体保持不变，转步骤4.4.3；
例如：参阅图2，R1表示当前适应度值最高的受体，R2为基因再排列的受体，R3为新的工件序列I。受体R1为(A,B,C,D),受体R2为(A,C,D,B),从受体R2中随机选取加工位置2，由于该位置处的工件C在受体R1中位于加工位置3，所以将受体R2中的工件C插入到位置3，得到新的工件序列R1为(A,D,C,B)。
步骤4.4.3：比对B细胞受体与当前适应度值最高的T细胞受体在相同加工位置处的工件，如果存在工件相同，则B细胞与T细胞发生同源相互作用，转步骤4.4.4，否则B细胞与T细胞不发生同源相互作用，转步骤4.4.6；
步骤4.4.4：采用部分随机排序的方式对B细胞受体进行细胞突变：对比B细胞受体与当前适应度值最高的T细胞受体在相同加工位置处的工件，工件相同情况下组成的加工位置号集合为工件不相同情况下组成的加工位置号集合为δ，集合中加工位置号上工件保持不变，对集合δ中加工位置号上的工作重新随机排序生成新的集合δ′，若解码的适应度值大于被选择的B细胞受体适应度值，则更新被选择的B细胞受体，否则被选择的B细胞受体保持不变，转步骤4.4.5；
例如：参阅图3，R4表示当前适应度值最高的T细胞受体，R5表示细胞突变的B细胞受体，R6表示新的工件序列II。受体R4为(A,B,C,D)，受体R5为(A,C,D,B)，与受体R4对比，受体R5在相同位置上相同的工件是A，相同位置上不同的工件是C,D,B，对工件C,D,B重新随机排序生成序列B,D,C，则新的工件序列R6为(A,B,D,C)。
步骤4.4.5：采用交换、插入、和颠倒的领域搜索机制对B细胞受体进行抗体转换：由于B细胞在T细胞受体的帮助下可以大量分化为分泌IgG、IgA、和IgE形态抗体的浆细胞进而消灭抗原，所以对被选择的B细胞受体，分别以1/3的概率进行交换、插入、或颠倒搜索进而达到局部最优，每个被选择的B细胞受体进行I次抗体转换，若抗体转换后B细胞受体的适应度值提升，则更新被选择的B细胞受体，否则不变，转步骤4.4.7；
步骤4.4.6：采用插入和交换的领域搜索机制对B细胞受体进行IgM抗体转换：由于B细胞在不与T细胞发生同源相互作用时，只是表达IgM抗体，提供一个快速的早期防御机制，所以被选择的B细胞受体先进行插入，再做交换，若IgM抗体的表达后B细胞受体的适应度值提升，则更新被选择的B细胞受体，否则不变，转步骤4.4.7；
步骤4.4.7：当前适应度值最高的受体保持不变，其他受体重新随机生成，转4.4.8；
步骤4.4.8：判断是否满足预先设定的最大CPU运行时间Tmax，如果满足，将当前最好的解作为最优解，得到最终的调度方案并停止，如果不满足转步骤4.4.2；
在当前的B细胞受体中选取目标函数最优的B细胞受体为(D,B,C,A)，其对应的目标函数值为19，具体的调度方案表如表2和表3所示：
表2并行同速机上各零件的种类、个数、开始时间、和完工时间
表3装配机上工件的开始时间和完工时间
最优调度甘特图如图4所示。
在本实施例中，将本方法与基于免疫球蛋白的人工免疫(IAIS)算法(AppliedSoft Computing,13(8),3729–3736)、变领域搜索(VNS)算法(European Journal ofOperational Research,130(3),449–467)进行性能比较，验证所提出的RAIS算法的有效性，实验环境采用C++语言编程，在CPU为Intel(R)Core(TM)i5-3210M 2.50GHz，内存为4.00GB的计算机上进行仿真试验。
加工参数设置：工件个数n∈{10,50,100}，并行同速机个数m∈{2,3,4}，零件的总个数q1j＝U[2,5]，零件的单个加工时间p1j＝U[1,6]，工件的装配时间Ai＝U[1,100]，准备时间与加工时间给出3个水平的比值：10％，20％，和30％，对应与这些比值，准备时间s2i＝U[1,10],U[1,20],或U[1,30]，装配工件所需的零件类型数hi＝U[2,5]。工件的单位提前惩罚成本αi＝U[1,10]，工件的单位滞后惩罚成本βi＝U[1,10]，交货期窗口的设置，先确定交货期窗口的中心DC＝max(0,U(PP(1-TF-RDD/2),PP(1-TF+RDD/2)))，其中滞后因素TF∈{0.2,0.6,1.0}，交货期窗口相对范围RDD∈{0.2,0.6,1.0}，PP为该车间调度问题最大完工时间的下界， 工件交货期窗口的起始时间ei＝max(0,DC-DC×H/100)；工件交货期窗口的结束时间di＝max(0,DC+DC×H/100)，其中交货期窗口宽度H＝U[1,10]。
RAIS算法参数设置如下：B细胞受体种群大小NB＝10，T细胞受体种群大小NT＝5，基因再排列选择工序号个数R＝n/5，抗体转换次数I＝5，最长CPU运行时间Tmax＝0.2×n秒。工件个数、并行同速机个数、准备时间与加工时间的比值、滞后因素、和交货期窗口相对范围的不同组合构成3×3×3×3×3＝243种不同的问题规模，为了消除随机因素对计算结果的影响，每种问题规模下产生5个计算算例，共产生1215个算例，每个实例运行10次，并报告平均计算结果。如表4所示，其中以偏差百分比PD来表示算法的性能。
其中，ZH表示由相关算法(VNS、IAIS、和RAIS)求得的目标函数值，Zbest表示三种算法求得的目标函数最小值。
表4针对1215个算例的运行结果比较
从表中可以看出，随着工件个数的增加，VNS、IAIS、和RAIS的偏差百分比呈增长趋势；VNS、IAIS、和RAIS的平均偏差百分比分别为188.36，6.86，和0.08，表明RAIS比VNS，IAIS有更好的求解性能。
Claims (4)
1.一种两阶段装配流水车间调度方法，其特征在于，该方法包括以下步骤：
步骤一、建立不同交货期窗口约束下的两阶段装配流水车间调度的目标函数；
步骤二、确定调度过程的约束条件；
步骤三、采用基于细胞受体的人工免疫系统算法即RAIS算法求解所述目标函数，经过迭代运算，将满足停止条件的解解码后作为最终的调度方案。
2.根据权利要求1所述的一种两阶段装配流水车间调度方法，其特征在于，步骤一中所述的目标函数如下：
i为工件号，n为工件个数；
αi为工件i的单位提前惩罚成本，βi为工件i的单位滞后惩罚成本；
Ei为工件i的提前时间，Ti为工件i的滞后时间；
Z为提前与滞后惩罚成本加权和。
3.根据权利要求1所述的一种两阶段装配流水车间调度方法，其特征在于，步骤二中所述的约束条件包括：
零件在并行同速机上不同位置约束：
零件拆分处理约束：
M×xjtk≥qjtk j,t＝1,2,...,h,k＝1,2,...m (6)
qjtk-1≥M×(xjtk-1) j,t＝1,2,...,h,k＝1,2,...m (7)
同一并行同速机不同位置上的零件完工时间约束：
C0k＝0 k＝1,2,...,m (8)
工件在装配机不同位置约束：
装配机不同位置上的工件完工时间约束：
Cr≥Fi+Ai-M×(1-yir) i,r＝1,2,...,n (13)
工件提前成本与滞后成本约束：
Ei≥ei-Cr-M×(1-yir) i,r＝1,2,...,n (15)
Ti≥Cr-di-M×(1-yir) i,r＝1,2,...,n (16)
Ei≥0 i＝1,2,...,n (17)
Ti≥0 i＝1,2,...,n (18)
上述公式(2)～(18)中，存在以下决策标量：
r为装配机上工件加工的位置号；
m为并行同速机个数；
k为并行同速机号；
h为总的零件种类数；
j为零部件种类号；
t为并行同速机上零件加工的位置号；
Ji为工件i的零件种类集；
s1j为加工零件j前的准备时间；
q1j为零件j的总个数；
p1j为零件j的单个加工时间；
s2i为加工工件i前的准备时间；
Ai为工件i的装配时间；
ei为工件i交货期窗口的起始时间；
di为工件i交货期窗口的结束时间；
Ctk为并行同速机k位置t上的零件的完工时间；
C(t-1)k为并行同速机k位置t-1上的零件的完工时间；
C0k为并行同速机k位置0上的完工时间即机器未进行加工状态下的时间点；
Fi为工件i进行装配操作的最早开始时间；
Cr为装配机位置r上的工件的完工时间；
Cr-1为装配机位置r-1上的工件的完工时间；
qjtk为并行同速机k位置t上的零件j的个数；
M为一个足够大的正整数；
在上述公式(2)～(18)中，存在以下决策变量：
。
4.根据权利要求1所述的一种两阶段装配流水车间调度方法，其特征在于，所述的步骤三的具体方法如下：
用RAIS求解不同交货期窗口约束下的TSAFSP包括编码、解码和适应度值计算的操作；
31)编码时，每个受体采用基于工件序列的编码，记作(π1,π2,...,πn)，即加工位置i处的工件为πi；
32)解码时，每个受体的解码步骤如下：
步骤3.2.1：依据公式(2)～(10)，采用拆分法确定每台并行同速机上零件种类、数量和完工时间；
a.设定i＝1,转b；
b.计算其中CT0k＝0，如果i≤n，则计算并对TPj按从大到小排序得转c；如果i＝n+1，则停止；
c.设定t＝1，MCTk＝0(k＝1,2,...,m)转d；
d.令w为min(MCT1,MCT2,...,MCTm)的机器号，如果q1t′＝0，则令t＝t+1转d；如果同时满足q1t′＞0，则qtw′＝q1′t，q1t′＝0，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw，令t＝t+1转d；如果同时满足 q1t′＞0，则q1t′＝q1t′-qtw′，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw，转d；如果同时满足q1t′＞0，则qtw′＝q1t′，q1t′＝0，CTiw＝MCTw+s1t′+p1t′×qtw′，MCTw＝CTiw令t＝t+1转d；如果且q1t′＝0，则i＝i+1转b；
步骤3.2.1中存在以下决策标量：
i为工件加工位置号，j为零件号，k为并行同速机器号；
为加工位置号为i的工件的零件在制造阶段的完工时间下界；
CTik为加工位置号为i的工件的零件在并行同速机k上的最大完工时间；
CTiw为加工位置号为i的工件的零件在并行同速机w上的最大完工时间；
CT(i-1)k为加工位置号为i-1的工件的零件在并行同速机k上的最大完工时间；
为装配工件πi所需的零件种类数；
为工件πi的零件j的不拆分总处理时间；
为装配工件πi所需的种零件中，第t大的不拆分总处理时间；
MCTk为并行同速机k的完工时间；
MCTw为并行同速机w的完工时间；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的总个数；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件在并行同速机w上的个数；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的准备时间；
为装配工件πi所需的种零件中，不拆分总处理时间第t大的零件的单个加工时间；
为对数a向上取整；
步骤3.2.2：依据公式(11)～(14)，采用空闲时间插入法确定装配机上工件的完工时间：
a.设定i＝1，ACT＝0，转b；
b.判断i≤n，如果不满足，设定i＝n，ξr＝0(r＝1,2,...,n)，l＝1，FJl＝LJl＝i转c；否则判断满足则ACT＝Ci，不满足则ACT＝Ci；令i＝i+1转b；
c.判断i≥1，如果不满足，则停止；否则判断i＝n，如果满足，则ξi＝1，否则判断如果满足，则ξi＝ξi+1，FJl＝i，否则ξi＝ξi+1+1，l＝l+1，FJl＝LJl＝i转d；
d.令σl＝{r|ξr＝ξi,r＝1,2,...,n}，若则若则若则ηr＝0；计算判断如果满足，则i＝i-1转c，否则转e；
e.若则若则Δ1r＝M；若则判断ξi＝1，如果满足，则Δ2l＝M，如果不满足，转f；
f.若则Cr＝Cr+Δ2l，ξr＝ξr-1，FJl-1＝FJl，l＝l-1转c；若则i＝i-1转c；
步骤3.2.2中存在以下决策标量：
CTi1为加工位置号为i的工件的零件在并行同速机1上的最大完工时间；
CTi2为加工位置号为i的工件的零件在并行同速机2上的最大完工时间；
CTim为加工位置号为i的工件的零件在并行同速机m上的最大完工时间；
Ci为装配机位置i上的工件的完工时间；
Ci+1为装配机位置i+1上的工件的完工时间；
为加工工件πi前的准备时间；
为加工工件πi+1前的准备时间；
为加工工件前的准备时间；
为工件πi的装配时间；
为工件πi+1的装配时间；
为工件的装配时间；
ACT为装配机的完工时间；
ξi为装配机上加工位置号为i的工件的标记号；
ξi+1为装配机上加工位置号为i+1的工件的标记号；
ξr为装配机上加工位置号为r的工件的标记号；
l为集合σl的标记号；
FJl为标记号为l的工件集中最早进行装配的工件的加工位置号；
FJl-1为标记号为l-1的工件集中最早进行装配的工件的加工位置号；
LJl为标记号为l的工件集中最晚进行装配的工件的加工位置号；
为装配机位置LJl上的工件的完工时间；
为装配机位置LJl-1上的工件的完工时间；
为工件πr的单位提前惩罚成本；
为工件πr的单位滞后惩罚成本；
ηi为装配机上加工位置号为i的工件的完工时间增加一个单位造成的成本减少值；
ηr为装配机上加工位置号为r的工件的完工时间增加一个单位造成的成本减少值；
为工件πr交货期窗口的起始时间；
为工件πr交货期窗口的结束时间；
σl为标记号为l的工件集；
Δ1r为工件在不改变其完工状态(提前、滞后或如期完工)时，完工时间允许的最大增加值；
Δ2l为集合σl中的工件在不影响其他工件集中的工作时，完工时间允许的最大增加值。
步骤3.2.3：依据公式(1)和(15)～(18)计算目标函数；
33)由于优化目标为提前和滞后惩罚成本加权和最小，采用目标函数的倒数作为适应度函数对每个受体进行评价；
34)RAIS算法求解不同交货期窗口约束下的TSAFSP的步骤如下：
步骤3.4.1：随机初始化NB个B细胞受体和NT个T细胞受体种群，转步骤3.4.2；
步骤3.4.2：采用选择和插入的方式对B细胞受体和T细胞受体种群进行基因再排列：首先将当前适应度值最高的受体作为标准受体，然后从其他受体的工件序列中随机选择R个加工位置号，将这R个加工位置号上对应的工件分别插入到其在标准受体中的加工位置上，若这个新的工件序列解码的适应度值大于被选择的受体适应度值，则更新被选择的受体，否则被选择的受体保持不变，转步骤3.4.3；
步骤3.4.3：比对B细胞受体与当前适应度值最高的T细胞受体在相同加工位置处的工件，如果存在工件相同，则B细胞与T细胞发生同源相互作用，转步骤4.4.4，否则B细胞与T细胞不发生同源相互作用，转步骤3.4.6；
步骤3.4.4：采用部分随机排序的方式对B细胞受体进行细胞突变：对比B细胞受体与当前适应度值最高的T细胞受体在相同加工位置处的工件，工件相同情况下组成的加工位置号集合为工件不相同情况下组成的加工位置号集合为δ，集合中加工位置号上工件保持不变，对集合δ中加工位置号上的工作重新随机排序生成新的集合δ′，若解码的适应度值大于被选择的B细胞受体适应度值，则更新被选择的B细胞受体，否则被选择的B细胞受体保持不变，转步骤3.4.5；
步骤3.4.5：采用交换、插入、和颠倒的领域搜索机制对B细胞受体进行抗体转换：由于B细胞在T细胞受体的帮助下可以大量分化为分泌IgG、IgA、和IgE形态抗体的浆细胞进而消灭抗原，所以对被选择的B细胞受体，分别以1/3的概率进行交换、插入、或颠倒搜索进而达到局部最优，每个被选择的B细胞受体共进行I次抗体转换，若抗体转换后B细胞受体的适应度值提升，则更新被选择的B细胞受体，否则不变，转步骤3.4.7；
步骤3.4.6：采用插入和交换的领域搜索机制对B细胞受体进行IgM抗体转换：由于B细胞在不与T细胞发生同源相互作用时，只是表达IgM抗体，提供一个快速的早期防御机制，所以被选择的B细胞受体先进行插入，再做交换，若IgM抗体的表达后B细胞受体的适应度值提升，则更新被选择的B细胞受体，否则不变，转步骤3.4.7；
步骤3.4.7：当前适应度值最高的受体保持不变，其他受体重新随机生成，转3.4.8；
步骤3.4.8：判断是否满足预先设定的最长CPU运行时间Tmax，如果满足，将当前最好的解作为最优解，得到最终的调度方案并停止，如果不满足转步骤3.4.2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711275686.5A CN108279647B (zh) | 2017-12-06 | 2017-12-06 | 一种两阶段装配流水车间调度方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711275686.5A CN108279647B (zh) | 2017-12-06 | 2017-12-06 | 一种两阶段装配流水车间调度方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108279647A true CN108279647A (zh) | 2018-07-13 |
| CN108279647B CN108279647B (zh) | 2019-04-23 |
Family
ID=62801324
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711275686.5A Expired - Fee Related CN108279647B (zh) | 2017-12-06 | 2017-12-06 | 一种两阶段装配流水车间调度方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108279647B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109325680A (zh) * | 2018-09-12 | 2019-02-12 | 山东大学 | 考虑学习效应的最小化延误总时间的单机批次排序方法 |
| CN109816262A (zh) * | 2019-01-31 | 2019-05-28 | 贵州大学 | 采用改良免疫遗传算法的柔性作业车间调度方法 |
| CN111144710A (zh) * | 2019-12-06 | 2020-05-12 | 重庆大学 | 一种可持续性混合流水车间的构建和动态调度方法 |
| CN112650187A (zh) * | 2021-01-22 | 2021-04-13 | 北京理工大学 | 一种车间调度方法、装置和系统 |
| CN114691327A (zh) * | 2022-03-23 | 2022-07-01 | 华南理工大学 | 面向两阶段任务调度的多目标群智能优化方法与系统 |
| CN117852840A (zh) * | 2024-03-07 | 2024-04-09 | 长春工业大学 | 一种基于多目标差分进化可变子批的柔性车间调度方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102929263A (zh) * | 2012-11-16 | 2013-02-13 | 北京理工大学 | 一种混合流水车间调度方法 |
| CN103116813A (zh) * | 2013-02-25 | 2013-05-22 | 西北工业大学 | 一种实时信息驱动的生产设备端制造活动智能导航方法 |
| CN104123587A (zh) * | 2014-04-04 | 2014-10-29 | 张涛 | Mto-mts管理模式下钢铁生产合同计划和多级库存匹配优化方法 |
| CN104392124A (zh) * | 2014-11-18 | 2015-03-04 | 广东工业大学 | 基于st的启发式算法的三阶段柔性流水车间调度方法 |
-
2017
- 2017-12-06 CN CN201711275686.5A patent/CN108279647B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102929263A (zh) * | 2012-11-16 | 2013-02-13 | 北京理工大学 | 一种混合流水车间调度方法 |
| CN103116813A (zh) * | 2013-02-25 | 2013-05-22 | 西北工业大学 | 一种实时信息驱动的生产设备端制造活动智能导航方法 |
| CN104123587A (zh) * | 2014-04-04 | 2014-10-29 | 张涛 | Mto-mts管理模式下钢铁生产合同计划和多级库存匹配优化方法 |
| CN104392124A (zh) * | 2014-11-18 | 2015-03-04 | 广东工业大学 | 基于st的启发式算法的三阶段柔性流水车间调度方法 |
Non-Patent Citations (1)
| Title |
|---|
| TSUI-PING CHUNG，HENG SUN，CHING-JONG LIAO: "Two new approaches for a two-stage hybrid flowshop problem with a single batch processing machine under waiting time constraint", 《COMPUTERS & INDUSTRIAL ENGINEERING》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109325680A (zh) * | 2018-09-12 | 2019-02-12 | 山东大学 | 考虑学习效应的最小化延误总时间的单机批次排序方法 |
| CN109816262A (zh) * | 2019-01-31 | 2019-05-28 | 贵州大学 | 采用改良免疫遗传算法的柔性作业车间调度方法 |
| CN109816262B (zh) * | 2019-01-31 | 2023-04-28 | 贵州大学 | 采用改良免疫遗传算法的柔性作业车间调度方法 |
| CN111144710A (zh) * | 2019-12-06 | 2020-05-12 | 重庆大学 | 一种可持续性混合流水车间的构建和动态调度方法 |
| CN111144710B (zh) * | 2019-12-06 | 2023-04-07 | 重庆大学 | 一种可持续性混合流水车间的构建和动态调度方法 |
| CN112650187A (zh) * | 2021-01-22 | 2021-04-13 | 北京理工大学 | 一种车间调度方法、装置和系统 |
| CN114691327A (zh) * | 2022-03-23 | 2022-07-01 | 华南理工大学 | 面向两阶段任务调度的多目标群智能优化方法与系统 |
| CN117852840A (zh) * | 2024-03-07 | 2024-04-09 | 长春工业大学 | 一种基于多目标差分进化可变子批的柔性车间调度方法 |
| CN117852840B (zh) * | 2024-03-07 | 2024-05-03 | 长春工业大学 | 一种基于多目标差分进化可变子批的柔性车间调度方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108279647B (zh) | 2019-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108279647B (zh) | 一种两阶段装配流水车间调度方法 | |
| Hyun et al. | A genetic algorithm for multiple objective sequencing problems in mixed model assembly lines | |
| Brintrup et al. | Ergonomic chair design by fusing qualitative and quantitative criteria using interactive genetic algorithms | |
| CN104751297A (zh) | 一种混流生产线产能分配方法 | |
| CN103272853A (zh) | 一种冷连轧各机架压下量与轧制速度的设定装置及方法 | |
| CN102411306A (zh) | 一种基于蜂社会自组织模型的混流装配生产调度控制方法 | |
| CN112862380B (zh) | 一种基于混合算法的项目型产品装配车间人员调度方法、设备及存储介质 | |
| CN107400935B (zh) | 基于改进elm的熔融纺丝工艺的调节方法 | |
| Badri et al. | A multi-criteria decision-making approach to solve the product mix problem with interval parameters based on the theory of constraints | |
| CN104463414B (zh) | 一种分时电价下热轧生产经济负荷调度方法 | |
| CN109934364A (zh) | 一种基于遗传模拟退火算法的刀具准时化配送方法 | |
| CN108681789A (zh) | 一种云制造服务优化方法 | |
| Lin et al. | A hybrid EA for reactive flexible job-shop scheduling | |
| CN104199820A (zh) | 云平台MapReduce工作流调度优化方法 | |
| CN107346469A (zh) | 云制造环境下多地运输多目标综合调度方法 | |
| CN110135752B (zh) | 一种具有切换时间的成套订单的调度方法 | |
| CN103761385A (zh) | 一种多热源环状管网的优化设计方法 | |
| CN110490446A (zh) | 一种基于改进遗传算法的模块化工艺重组方法 | |
| Löppenberg et al. | Self optimisation and automatic code generation by evolutionary algorithms in PLC based controlling processes | |
| Mu | Disruption Management of Flexible Job Shop Scheduling Considering Behavior Perception and Machine Fault Based on Improved NSGA-II Algorithm. | |
| CN110648037A (zh) | 一种整车生产评价方法及装置 | |
| LU502930B1 (en) | Flexible Job Shop Scheduling Method with Transportation Time and Adjustment Time | |
| CN116774657A (zh) | 一种基于鲁棒优化的再制造车间动态调度方法 | |
| Han et al. | Multiple rules decision-based DE solution for the earliness-tardiness case of hybrid flow-shop scheduling problem | |
| CN115169798A (zh) | 一种带准备时间分布式柔性作业车间调度方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190423 Termination date: 20191206 |