CN105279605B - 一种动态变化权重下的批加工机调度方法 - Google Patents

一种动态变化权重下的批加工机调度方法 Download PDF

Info

Publication number
CN105279605B
CN105279605B CN201510624138.3A CN201510624138A CN105279605B CN 105279605 B CN105279605 B CN 105279605B CN 201510624138 A CN201510624138 A CN 201510624138A CN 105279605 B CN105279605 B CN 105279605B
Authority
CN
China
Prior art keywords
indicate
buffer
product
processing machine
product family
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.)
Expired - Fee Related
Application number
CN201510624138.3A
Other languages
English (en)
Other versions
CN105279605A (zh
Inventor
贾文友
刘莉
贾昊瑞
刘涛
丁洋洋
陈开开
黄西西
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Anhui Polytechnic University
Original Assignee
Anhui Polytechnic University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Anhui Polytechnic University filed Critical Anhui Polytechnic University
Priority to CN201510624138.3A priority Critical patent/CN105279605B/zh
Publication of CN105279605A publication Critical patent/CN105279605A/zh
Application granted granted Critical
Publication of CN105279605B publication Critical patent/CN105279605B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • General Factory Administration (AREA)

Abstract

一种动态变化权重下的批加工机调度方法,包括产生触发事件,确定调度对象,计算动态变化权重,获取优先权最高的产品族,对空闲可用的批加工机派工和装载,调度终止条件判断;在极小化加权总完工时间的调度目标中,空闲可用批加工机缓冲器里待加工产品为研究对象;根据不同产品族的待加工产品数量的动态变化情况等计算动态变化权重;根据动态权重确定极小化加权总完工时间的目标下优先权最高的产品族;按“满批原则”组批;“先到先服务规则”排序的第一组批被派工和装载;通过滚动循环执行,直至满足调度终止条件。它实现方便,调度精度高,利于半导体等有批加工机的行业推广应用。

Description

一种动态变化权重下的批加工机调度方法
所属技术领域
本发明涉及生产线的批加工机调度,尤其涉及一种动态变化权重下的批加工机调度方法。
背景技术
批加工是指在不超过加工设备的最大加工能力时,一次性可以加工多个待加工产品的加工操作。此加工设备称作批加工机;每次实际加工的多个未待加工产品称为一个批。如,在半导体芯片生产线前段的炉管区,该区的加工设备属于的批加工机,每次实际的加工的多个产品称作一个批。批加工机具有价格昂贵,加工时间长,产品具有不同到达时间,需要不断更新组批、排序,以及不同产品族不能混合组批等特性,是生产系统中的瓶颈机台。故批加工机调度制约着生产系统的整体绩效。批加工机的优化调度控制是一个非多项式的难问题,一直是学术界与应用界的研究热点。
有许多调度目标可以衡量批加工机的调度效率,其中众多的调度目标的衡量式常有权重的参数,如在极小化加权总完工时间的调度目标中,就包含权重。
目前,对于权重的参数的常规赋值方法是根据生产经验对不同产品族的待加工产品赋予不同定值。但是,每当批加工机空闲可用时,因不同产品族的待加工产品数量变化、相同产品族内的待加工产品到达时间不同等动态变化情况,以及具有等待时间限制等约束要求,根据生产经验给权重所赋予的定值不能反映待加工产品的这些动态变化情况和约束要求,致使批加工机优化调度过程偏离实际情况,使得优化调度解的精度降低。
发明内容
为了克服现有因不同产品族的待加工产品数量变化、相同产品族内的待加工产品到达时间不同等动态变化情况,以及具有等待时间限制等约束要求,根据生产经验给权重所赋予的定值不能反映待加工产品的这些动态变化情况和约束要求,致使批加工机优化调度过程偏离实际情况,使得优化调度解的精度降低等缺陷,本发明提供一种动态变化权重下的批加工机调度方法。
本发明具体解决其技术问题所采用的技术解决方案是:
一种动态变化权重下的批加工机调度方法,以当批加工机空闲可用时其缓冲器里不同产品族的待加工产品为研究对象。
动态变化权重包括不同产品族的待加工产品数量的动态变化情况,相同产品族内的待加工产品到达时间不同的动态变化情况,以及具有等待时间限制的约束要求,动态变化权重计算模型如下
其中
qf≥(tf)max,当qf=(tf)max时,
nf≤n,当nf=n时,
这里F表示缓冲器里产品族的种类总数量;μf表示缓冲器里产品族f的动态变化权重,其中f表示对应着F中第f个产品族;wf表示缓冲器里产品族f的固有权重;qf表示缓冲器里产品族f在缓冲器里的等待约束限制时间;tf表示缓冲器里产品族f在缓冲器里的已等待时间;ti表示缓冲器里产品族f中待加工工件i在缓冲器里的已等待时间,其中i表示对应着产品族f中第i个待加工工件;n表示缓冲器里所有产品族的待加工工件总数量;nf表示缓冲器里产品族f中包含待加工工件数量;(tf)max表示缓冲器里产品族f中所有待加工工件在缓冲器里的最大已等待时间;表示对ti从i=1到i=nf进行求最大值运算;“”表示μf计算值与qf,(tf)max,wf,n和nf具体数值相关;“qf≥(tf)max,当qf=(tf)max时,”表示必须满足qf≥(tf)max,且规定当qf=(tf)max时,取值无穷大,即qf=(tf)max所对应的产品族f具有无穷大动态变化权重值;“nf≤n,当nf=n时,”表示必须满足nf≤n,且规定当nf=n时,取值无穷大,即nf=n所对应的产品族f具有无穷大动态变化权重值;表示任取第f个产品族都属于缓冲器里产品族的种类总数量F中。
“满批原则”组批定义是按除最后一个批中的待加工工件数量大于或等于批加工机的最小加工限制对应待加工工件数量外,其余各批中待加工工件数量必须等于批加工机的最大加工能力对应待加工工件数量,“满批原则”组批模型如下
这里hf表示缓冲器里产品族f可组批的数量;Bmax表示批加工机的最大加工能力对应待加工工件数量;Bmin表示批加工机的最小加工限制对应待加工工件数量;表示返回不大于nf除以Bmax的商的最大整数;xf表示0-1变量;MOD(nf,Bmax)表示nf除以Bmax的余数;“”表示hf计算值与和xf具体数值相关;“”用于确定xf取值,当满足MOD(nf,Bmax)≥Bmin时,xf取值1,否则xf取值0;“”表示任取第f个产品族都属于缓冲器里产品族的种类总数量F中。
动态变化权重下的批加工机调度方法,批加工机调度步骤如下:
步骤1,产生触发事件:一台批加工机处于空闲可用状态;
步骤2,确定调度对象:步骤1中的空闲可用的批加工机的缓冲器里不同产品族的待加工产品为调度对象;
步骤3,计算动态变化权重:根据动态变化权重计算模型进行计算;
步骤4,获取优先权最高的产品族:根据不同产品族的动态变化权重确定在确定的调度目标下,建立和运行数学模型;
步骤5,对空闲可用的批加工机派工和装载:对所获取的优先权最高的产品族中的不同待加工工件根据“满批规则”组批,在确定的调度目标下进行批排序,优先级最高批中的待加工工件被派工和装载到步骤1的空闲可用的批加工机;
步骤6,调度终止条件判断:如果不符合所规定的整个实时调度任务的调度终止条件,进行滚动执行,跳转到步骤1;否则动态变化权重下的批加工机调度方法终止。
调度终止条件是所有的待调度工件已被装载到批处理机上或制造系统的实时调度总的运行时间域结束。
本发明的有益效果是,采用动态变化权重下的批加工机调度方法,通过计算和使用动态改变权重,解决了目前因不同产品族的待加工产品数量变化、相同产品族内的待加工产品到达时间不同等动态变化情况,以及具有等待时间限制等约束要求,根据生产经验给权重所赋予的定值不能反映待加工产品的这些动态变化情况和约束要求,致使批加工机优化调度过程偏离实际情况,使得优化调度解的精度降低等问题。它概念简单,实现方便,调度精度高,适应半导体等具有批加工机的行业实际优化调度需求。
附图说明
图1为本发明动态变化权重下的批加工机调度方法的算法流程图。
具体实施方式
下面结合附图和实施例对本发明专利作进一步地说明:
本发明所提供的一种动态变化权重下的批加工机调度方法,包括产生触发事件,确定调度对象,计算动态变化权重,获取优先权最高的产品族,对空闲可用的批加工机派工和装载,调度终止条件判断。
动态变化权重包括不同产品族的待加工产品数量的动态变化情况,相同产品族内的待加工产品到达时间不同的动态变化情况,以及具有等待时间限制的约束要求,动态变化权重计算模型如下
其中
qf≥(tf)max,当qf=(tf)max时,
nf≤n,当nf=n时,
这里F表示缓冲器里产品族的种类总数量;μf表示缓冲器里产品族f的动态变化权重,其中f表示对应着F中第f个产品族;wf表示缓冲器里产品族f的固有权重;qf表示缓冲器里产品族f在缓冲器里的等待约束限制时间;tf表示缓冲器里产品族f在缓冲器里的已等待时间;ti表示缓冲器里产品族f中待加工工件i在缓冲器里的已等待时间,其中i表示对应着产品族f中第i个待加工工件;n表示缓冲器里所有产品族的待加工工件总数量;nf表示缓冲器里产品族f中包含待加工工件数量;(tf)max表示缓冲器里产品族f中所有待加工工件在缓冲器里的最大已等待时间;表示对ti从i=1到i=nf进行求最大值运算;“”表示μf计算值与qf,(tf)max,wf,n和nf具体数值相关;“qf≥(tf)max,当qf=(tf)max时,”表示必须满足qf≥(tf)max,且规定当qf=(tf)max时,取值无穷大,即qf=(tf)max所对应的产品族f具有无穷大动态变化权重值;“nf≤n,当nf=n时,”表示必须满足nf≤n,且规定当nf=n时,取值无穷大,即nf=n所对应的产品族f具有无穷大动态变化权重值;“”表示任取第f个产品族都属于缓冲器里产品族的种类总数量F中。
“满批原则”组批定义是按除最后一个批中的待加工工件数量大于或等于批加工机的最小加工限制对应待加工工件数量外,其余各批中待加工工件数量必须等于批加工机的最大加工能力对应待加工工件数量,“满批原则”组批模型如下
这里hf表示缓冲器里产品族f可组批的数量;Bmax表示批加工机的最大加工能力对应待加工工件数量;Bmin表示批加工机的最小加工限制对应待加工工件数量;表示返回不大于nf除以Bmax的商的最大整数;xf表示0-1变量;MOD(nf,Bmax)表示nf除以Bmax的余数;“”表示hf计算值与和xf具体数值相关;“”用于确定xf取值,当满足MOD(nf,Bmax)≥Bmin时,xf取值1,否则xf取值0;“”表示任取第f个产品族都属于缓冲器里产品族的种类总数量F中。
一种动态变化权重下的批加工机调度方法,在极小化加权总完工时间的调度目标中,以当批加工机空闲可用时其缓冲器里不同产品族的待加工产品为研究对象,批加工机调度步骤如下:
步骤1,产生触发事件:一台批加工机处于空闲可用状态;
步骤2,确定调度对象:步骤1中的空闲可用的批加工机的缓冲器里不同产品族的待加工产品为调度对象;
步骤3,计算动态变化权重:根据动态变化权重计算模型进行计算;
步骤4,获取优先权最高的产品族:根据不同产品族的动态变化权重确定在极小化加权总完工时间的调度目标下,建立和运行下述数学模型
极小化目标
约束条件
Cf=s(f)×Pf×hf
s(f)∈{1,2,…,F}
如果f≠l;那么s(f)≠s(l)
这里Cf表示缓冲器里产品族f的完工时间;“”表示缓冲器里所有产品族的加权总完工时间;优化序列位置为一维数组,表示为:s(1),s(2),…,s(f),…,s(F);Pf表示产品族f在批加工机上的加工时间;l表示对应着F中第l个产品族;“Cf=s(f)×Pf×hf”表示hf计算值与s(f),Pf和hf具体数值乘积相关;“s(f)∈{1,2,…,F}”表示优化序列位置一维数组中优化序列位置总数量等于缓冲器里产品族的种类总数量F;“如果f≠l;那么s(f)≠s(l)”表示不同的第l个产品族和第f个产品族所对应优化序列位置一维数组中优化序列位置s(l)和优化序列位置s(f)是不同的;“”表示任取第f个产品族都属于缓冲器里产品族的种类总数量F中;“”表示任取第l个产品族都属于缓冲器里产品族的种类总数量F中;
步骤5,对空闲可用的批加工机派工和装载:对所获取的优先权最高的产品族中的不同待加工工件的按其到达缓冲器的时间不增的“先到先服务规则”排序;根据“满批规则”组建该排序工件的第一个批;所组建的第一个批的待加工工件被派工和装载到步骤1的空闲可用的批加工机;
步骤6,调度终止条件判断:如果不符合所规定的整个实时调度任务的调度终止条件,进行滚动执行,跳转到步骤1;否则极小化加权总完工时间的调度目标下动态变化权重下的批加工机调度方法终止。
在极小化加权总完工时间的调度目标下,运行动态变化权重下的批加工机调度方法的输出是极小化加权总完工时间的目标函数值和优化序列位置一维数组值。
调度终止条件是所有的待调度工件已被装载到批处理机上或制造系统的实时调度总的运行时间域结束。

Claims (4)

1.一种动态变化权重下的批加工机调度方法,包括批加工机调度步骤:步骤1,产生触发事件:一台批加工机处于空闲可用状态;步骤2,确定调度对象:步骤1中的空闲可用的批加工机的缓冲器里不同产品族的待加工产品为调度对象;步骤3,计算动态变化权重:根据动态变化权重计算模型进行计算;步骤4,获取优先权最高的产品族:根据不同产品族的动态变化权重确定在极小化加权总完工时间的调度目标下,建立和运行数学模型;步骤5,对空闲可用的批加工机派工和装载:对所获取的优先权最高的产品族中的不同待加工工件,按其到达缓冲器的时间不增的“先到先服务规则”排序;根据“满批规则”组建该排序工件的第一个批;所组建的第一个批的待加工工件被派工和装载到步骤1的空闲可用的批加工机;步骤6,调度终止条件判断:如果不符合所规定的整个实时调度任务的调度终止条件,进行滚动执行,跳转到步骤1;否则极小化加权总完工时间的调度目标下动态变化权重下的批加工机调度方法终止,其特征在于,所述动态变化权重包括不同产品族的待加工产品数量的动态变化情况,相同产品族内的待加工产品到达时间不同的动态变化情况,以及具有等待时间限制的约束要求,动态变化权重计算模型如下
其中
qf≥(tf)max,当qf=(tf)max时,
nf≤n,当nf=n时,
这里F表示缓冲器里产品族的种类总数量;μf表示缓冲器里产品族f的动态变化权重,其中f表示对应着F中第f个产品族;wf表示缓冲器里产品族f的固有权重;qf表示缓冲器里产品族f在缓冲器里的等待约束限制时间;tf表示缓冲器里产品族f在缓冲器里的已等待时间;ti表示缓冲器里产品族f中待加工工件i在缓冲器里的已等待时间,其中i表示对应着产品族f中第i个待加工工件;n表示缓冲器里所有产品族的待加工工件总数量;nf表示缓冲器里产品族f中包含待加工工件数量;(tf)max表示缓冲器里产品族f中所有待加工工件在缓冲器里的最大已等待时间;表示对ti从i=1到i=nf进行求最大值运算;“”表示μf计算值与qf,(tf)max,wf,n和nf具体数值相关;“qf≥(tf)max,当qf=(tf)max时,”表示必须满足qf≥(tf)max,且规定当qf=(tf)max时,取值无穷大;“nf≤n,当nf=n时,”表示必须满足nf≤n,且规定当nf=n时,取值无穷大;“”表示任取第f个产品族都属于缓冲器里产品族的种类总数量F中。
2.根据权利要求1所述的一种动态变化权重下的批加工机调度方法,其特征在于,所述“满批原则”组批是按除最后一个批中的待加工工件数量大于或等于批加工机的最小加工限制对应待加工工件数量外,其余各批中待加工工件数量必须等于批加工机的最大加工能力对应待加工工件数量,“满批原则”组批模型如下
这里hf表示缓冲器里产品族f可组批的数量;Bmax表示批加工机的最大加工能力对应待加工工件数量;Bmin表示批加工机的最小加工限制对应待加工工件数量;表示返回不大于nf除以Bmax的商的最大整数;xf表示0-1变量;MOD(nf,Bmax)表示nf除以Bmax的余数;“”表示hf计算值与和xf具体数值相关;“”确定xf取值,当满足MOD(nf,Bmax)≥Bmin时,xf取值1,否则xf取值0;“”表示任取第f个产品族都属于缓冲器里产品族的种类总数量F中。
3.根据权利要求1所述的一种动态变化权重下的批加工机调度方法,其特征在于,所述不同产品族的动态变化权重确定在极小化加权总完工时间的调度目标下,建立和运行数学模型如下
极小化目标
约束条件
Cf=s(f)×Pf×hf
s(f)∈{1,2,…,F}
如果f≠l;那么s(f)≠s(l)
这里hf表示缓冲器里产品族f可组批的数量;Cf表示缓冲器里产品族f的完工时间;“”表示缓冲器里所有产品族的加权总完工时间;优化序列位置为一维数组,表示为:s(1),s(2),…,s(f),…,s(F);Pf表示产品族f在批加工机上的加工时间;l表示对应着F中第l个产品族;“Cf=s(f)×Pf×hf”表示hf计算值与s(f),Pf和hf具体数值乘积相关;“s(f)∈{1,2,…,F}”表示优化序列位置一维数组中优化序列位置总数量等于缓冲器里产品族的种类总数量F;“如果f≠l;那么s(f)≠s(l)”表示不同的第l个产品族和第f个产品族所对应优化序列位置一维数组中优化序列位置s(l)和优化序列位置s(f)是不同的;“”表示任取第f个产品族都属于缓冲器里产品族的种类总数量F中;“”表示任取第l个产品族都属于缓冲器里产品族的种类总数量F中。
4.根据权利要求1所述的一种动态变化权重下的批加工机调度方法,其特征在于,所述调度终止条件是所有的待调度工件已被装载到批处理机上或制造系统的实时调度总的运行时间域结束。
CN201510624138.3A 2015-09-25 2015-09-25 一种动态变化权重下的批加工机调度方法 Expired - Fee Related CN105279605B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510624138.3A CN105279605B (zh) 2015-09-25 2015-09-25 一种动态变化权重下的批加工机调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510624138.3A CN105279605B (zh) 2015-09-25 2015-09-25 一种动态变化权重下的批加工机调度方法

Publications (2)

Publication Number Publication Date
CN105279605A CN105279605A (zh) 2016-01-27
CN105279605B true CN105279605B (zh) 2018-10-12

Family

ID=55148577

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510624138.3A Expired - Fee Related CN105279605B (zh) 2015-09-25 2015-09-25 一种动态变化权重下的批加工机调度方法

Country Status (1)

Country Link
CN (1) CN105279605B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113267609B (zh) * 2021-05-27 2022-09-13 重庆钢铁股份有限公司 高炉喷吹用煤质量评价方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102129242A (zh) * 2011-04-12 2011-07-20 上海大学 基于两层混合智能优化的批处理生产过程产品质量控制方法
CN103136631A (zh) * 2013-01-22 2013-06-05 安徽工程大学 基于产品族排序实现最小总完工时间的批调度方法
CN104460590A (zh) * 2014-10-23 2015-03-25 北京化工大学 一种半导体生产线多产品工件合并方法
JP5982187B2 (ja) * 2012-06-13 2016-08-31 株式会社エヌ・ティ・ティ・データ 作業計画管理装置、作業計画管理方法、および作業計画管理プログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102129242A (zh) * 2011-04-12 2011-07-20 上海大学 基于两层混合智能优化的批处理生产过程产品质量控制方法
JP5982187B2 (ja) * 2012-06-13 2016-08-31 株式会社エヌ・ティ・ティ・データ 作業計画管理装置、作業計画管理方法、および作業計画管理プログラム
CN103136631A (zh) * 2013-01-22 2013-06-05 安徽工程大学 基于产品族排序实现最小总完工时间的批调度方法
CN104460590A (zh) * 2014-10-23 2015-03-25 北京化工大学 一种半导体生产线多产品工件合并方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
批处理机生产与成批配送的协调调度研究;宫华等;《沈阳理工大学学报》;20130630;第32卷(第3期);全文 *
面向产品族优化时间窗下可重入批处理机调度;贾文友等;《机械工程学报》;20150630;全文 *

Also Published As

Publication number Publication date
CN105279605A (zh) 2016-01-27

Similar Documents

Publication Publication Date Title
CN102110015B (zh) 一种任务进度和任务处理时间的确定方法及装置
CN104504540A (zh) 基于多阶段智能优化算法的动态柔性作业车间调度控制方法
CN108171372B (zh) 一种存在准备时间的多品种生产分批调度方法
CN103136631B (zh) 基于产品族排序实现最小总完工时间的批调度方法
CN111356988B (zh) 用于调度半导体加工的方法和系统
CN101604409B (zh) 应用于多重入复杂制造系统的调度方法
CN104156266B (zh) 确定实时任务或事件可调度性测试最小区间的方法
Singholi et al. Impact of manufacturing flexibility on FMS performance: a simulation study
CN105279605B (zh) 一种动态变化权重下的批加工机调度方法
CN111857081A (zh) 基于Q-learning强化学习的芯片封装测试生产线性能控制方法
CN105204433B (zh) 一种随机加工间隔内机床状态切换的机床节能运行方法
CN110806737B (zh) 一种基于最小能耗及最小时间的生产线设备数量优化方法
CN103199040A (zh) 一种暖机控制方法
CN106327053B (zh) 一种基于多模式集合的纺织工艺推荐模型的构建方法
CN108829036A (zh) 一种金属铸坯切削成形加工过程的优化调度方法
CN103116809B (zh) 面向产品族排序的批处理机的调度装置及方法
CN113848904A (zh) 一种基于惩罚能耗优化多移动机器人任务分配的方法
CN111290360B (zh) 一种铸造生产线多目标优化方法
CN108873850A (zh) 一种机械制造自动化生产过程的优化调度方法
TW202004642A (zh) 用於改善生產週期的方法
CN115081763A (zh) 一种数控车床车削过程能耗预测系统
CN110750074B (zh) 基于物理仿真的数控加工优化系统及方法
CN106033217B (zh) 派工系统及派工方法
CN113641101A (zh) 一种基于数值仿真的多渠池控制参数寻优算法
Monostori et al. Multipurpose modelling and optimisation of production processes and process chains by combining machine learning and search techniques

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20181012

Termination date: 20190925