CN111352397B - 基于数学规划算法的自动制造系统状态稳健性检测方法 - Google Patents
基于数学规划算法的自动制造系统状态稳健性检测方法 Download PDFInfo
- Publication number
- CN111352397B CN111352397B CN202010051204.3A CN202010051204A CN111352397B CN 111352397 B CN111352397 B CN 111352397B CN 202010051204 A CN202010051204 A CN 202010051204A CN 111352397 B CN111352397 B CN 111352397B
- Authority
- CN
- China
- Prior art keywords
- transition
- unreliable
- state
- resources
- sequence
- 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.)
- Active
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 47
- 238000001514 detection method Methods 0.000 title claims description 5
- 230000007704 transition Effects 0.000 claims abstract description 130
- 238000000034 method Methods 0.000 claims abstract description 91
- 230000008569 process Effects 0.000 claims abstract description 72
- 230000014509 gene expression Effects 0.000 claims abstract description 24
- 230000005540 biological transmission Effects 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims description 27
- 238000013508 migration Methods 0.000 claims description 6
- 238000003745 diagnosis Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000005012 migration Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 3
- 101000879840 Homo sapiens Serglycin Proteins 0.000 description 3
- 102100037344 Serglycin Human genes 0.000 description 3
- 101000999079 Homo sapiens Radiation-inducible immediate-early gene IEX-1 Proteins 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
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/41885—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 modeling, simulation of the manufacturing system
-
- 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/32339—Object oriented modeling, design, analysis, implementation, simulation language
-
- 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)
- Manufacturing & Machinery (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Multi Processors (AREA)
- Debugging And Monitoring (AREA)
Abstract
本发明属于自动制造系统技术领域,公开了一种基于数学规划算法的自动制造系统状态稳健性检测方法,对于自动制造系统的Petri网模型,将含有不可靠资源的进程的变迁划分为两类:一是不可靠资源的前序变迁,二是不可靠资源的后序变迁;对不可靠资源损坏进行文字叙述到数学表达式的转化;列出资源损坏下的数学规划算法(MP)表达式,找出使能变迁发射序列;通过验证上述变迁发射序列是否满足不使用不可靠资源的其他进程不被阻塞,得到当前状态的稳健性。本发明首先比没有考虑不可靠资源的研究,更加接近于实际自动制造系统的生产情况;其次将原先对于含有不可靠资源的自动制造系统停留在文字叙述提升到了数学表达式上,为相应算法的提出实现了可能。
Description
技术领域
本发明属于自动制造系统技术领域,尤其涉及一种基于数学规划算法的自动制造系统状态稳健性检测方法。
背景技术
目前,最接近的现有技术:自动制造系统(AMS)是指在较少人为因素的干扰下,加工装置自动地将原材料加工并组装成产品。自动制造系统一般是由众多小规模的,局部的,交互的,异步的以及事件驱动的并发子系统构成。考虑到资源成本、系统体积,资源的数量将会是有限的。因此资源共享在自动制造系统中屡见不鲜,正因为各个子系统之间存在着资源共享,那么各子系统之间竞争有限的资源将会随时发生,可能伴随着一个产品加工进程的所有加工阶段,资源竞争将会导致不期望的死锁以及阻塞等等;其次,在实际的自动制造系统中,有限的资源还会面临着资源故障这一大大的潜在危害,资源故障在实际的生产中也是很常见的现象,比如电器失灵,信号丢失,工件损坏等等。
一个典型的自动制造系统通常包含多个加工进程,每个加工进程分别负责加工生产不同的工件,当一个进程由于系统资源被损坏时,其加工进程将不得不被终止,但是人们总是希望其他不使用该损坏资源的加工进程不受影响,继续循环往复地加工工件。然而实际的生产过程中并非如此,一个加工进程因为损坏的资源被终止,往往会阻塞其他加工进程。那么在实际的自动制造系统中,一个加工进程受资源故障的影响,其他进程是否被阻塞,我们需要找出一种有效的方式进行诊断。在现有的技术中,根据资源是否发生故障,将所有资源分为可靠资源和不可靠资源,采用Petri网对自动制造系统可进行建模,一个含有不可靠资源的加工进程是否会阻塞其他不使用该不可靠资源的加工进程,被定义为该系统的稳健性。给定任意一个含有不可靠资源的自动制造系统当前所处状态,如何快速诊断出该状态的稳健性,现有的研究技术中并未给出有效的方法。
综上所述,现有技术存在的问题是:现有技术没有考虑实际的自动制造系统含有不可靠资源这一常见现象,而不考虑不可靠资源这一因素研究的结论都是处于理想状态,并不能真是反应自动制造系统的实际生产情况;并对含有不可靠资源的自动制造系统的稳健性给出一些诊断方法,在处理过程中,首先方法比较保守,其次诊断的状态稳健性也比较保守,进而不能满足自动制造系统的许可性更大。
解决上述技术问题的难度:首先,对于采用Petri网模型对含有不可靠资源的自动制造系统进行建模,不可靠资源被损坏,其定义仅仅停留在文字叙述,那么从文字叙述到数学表达式的转化这是其中最关键的一步;其次确定不可靠资源损坏导致的状态稳健是否存在类似于死锁问题中的坏状态这样的中间状态;最后快速诊断系统当前所处状态是否稳健及满足最大许可性。
解决上述技术问题的意义:实际的自动制造系统因为存在不可靠资源导致系统阻塞的问题,对于不可靠资源被损坏这一定义的数学化,为解决含有不可靠资源的自动制造系统的稳健性提供了新思路,以及新的解决方案。新的解决方案在解决系统最大许可性方面实现了突破。
发明内容
针对现有技术存在的问题,本发明提供了一种基于数学规划算法的自动制造系统状态稳健性检测方法。
本发明是这样实现的,一种基于数学规划算法的自动制造系统状态稳健性检测方法,所述基于数学规划算法的自动制造系统状态稳健性检测方法包括:
第一步,将含有不可靠资源的进程的变迁划分为两类:一是不可靠资源的前序变迁,二是不可靠资源的后序变迁;
第二步,对不可靠资源损坏进行文字叙述到数学表达式的转化;
第三步,列出资源损坏下的数学规划算法(MP)表达式,找出使能变迁发射序列;
第四步,通过验证上述变迁发射序列是否满足不使用不可靠资源的其他进程不被阻塞,得到当前状态的稳健性。
进一步,所述第一步还包括:划分变迁类,对于使用不可靠资源的加工进程,按照加工阶段的次序,把该进程中的所有变迁划分为不可靠资源的前序变迁和不可靠资源的后序变迁;
(1)输入Petri网中含有不可靠资源的进程的所有变迁、库所以及不可靠资源;
(3)按照该进程工件的加工顺序,即从当前进程的闲置库所P0开始,托肯的流动方向,以不可靠资源为分界线,从闲置库所的后置变迁{P0}·到不可靠资源的后置变迁{PRu}·,其中包含的变迁均为不可靠资源的前序变迁,依次加入前序变迁集合;从不可靠资源的前置变迁·{PRu}到闲置库所的前置变迁·{P0},其中所包含的变迁均为不可靠资源的后序变迁,依次加入后续变迁集合。
进一步,所述第二步还包括:从数学上定义不可靠资源被损坏,得到不可靠资源被损坏的数学意义,即数学表达式表述;
(1)不可靠资源被损坏的定义为:工作站中的处理器无法正常工作,而缓冲区仍将正常存放工件,且资源损坏不会影响正在被加工的工件;如果在缓冲区中存在没被使用的空闲位置,那么已经完成上一阶段加工任务的工件就进入到该缓冲区;另外对于已经完成当前操作的工件,则可以正常进入到下一个阶段的加工中,但对于没有加工完成或者正在等待资源的工件,将一直保持请求状态,直到该资源被修复;
(4)数学上定义不可靠资源损坏为:不可靠资源被损坏,则不可靠资源的后序变迁集合都将禁止发射,即|Tback|=[0],对于不可靠资源的前序变迁集合以及其他进程的变迁则可以任意发射。
进一步,所述第三步还包括:以不可靠资源被损坏的数学表达式作为约束条件,得出含有不可靠资源的数学规划算法,在当前状态下找出满足条件的变迁发射序列;
(1)输入任意一个状态M,步骤一所得后序变迁集合Tback,所有库所,变迁集合以及整个Petri网的关联矩阵[N];
(3)数学规划算法表达式:
式(1)在当前状态下,存在一个变迁序列,使得系统回到又当前状态;式(2)变迁序列满足不可靠资源的后序变迁均禁止发射;式(3)变迁序列至少存在一个变迁发射,即当前状态不能静止不动。
进一步,所述第四步还包括:在当前状态下,生成各个不使用不可靠资源的进程的可达图,即局部网的可达图,在各个局部网的可达图中寻找是否存在该变迁发射序列,使得从当前状态发射该变迁序列又回到当前状态,诊断出当前状态的稳健性;
(2)初始化当前状态Mcurrent=M'0,标记为flag=false,初始状态已访问,则标记为true;
(3)若系统中还存在没有被搜索到的状态,则继续执行下面的过程,否则终止;
(4)选择一个状态为“false”,即没有被访问过的标识M;
(5)对于M下所有满足使能条件的变迁t,执行下面的操作:
1)激发t得到新的标识M';
2)从M到M'添加弧t;
(6)将M的状态标记为“true”后回到(2);
(7)生成各个不使用不可靠资源的进程的局部网可达图;
进一步,所述选择一个状态为“false”,即没有被访问过的标识M包括:
1)若已搜索过M,则将其标记为“true”,并开始搜索其它“false”标识;
2)若在M下,不存在可以被使能的变迁,则将M标记为“deadlock”。
本发明的另一目的在于提供一种接收用户输入程序存储介质,所存储的计算机程序使电子设备执行包括下列步骤:
第一步,将含有不可靠资源的进程的变迁划分为两类:一是不可靠资源的前序变迁,二是不可靠资源的后序变迁;
第二步,对不可靠资源损坏进行文字叙述到数学表达式的转化;
第三步,列出资源损坏下的数学规划算法(MP)表达式,找出使能变迁发射序列;
第四步,通过验证上述变迁发射序列是否满足不使用不可靠资源的其他进程不被阻塞,得到当前状态的稳健性。
本发明的另一目的在于提供一种存储在计算机可读介质上的计算机程序产品,包括计算机可读程序,供于电子装置上执行时,提供用户输入接口以实施所述的基于数学规划算法的自动制造系统状态稳健性检测方法。
本发明的另一目的在于提供一种所述基于数学规划算法的自动制造系统状态稳健性检测方法在自动制造系统中的应用。
综上所述,本发明的优点及积极效果为:本发明首先比没有考虑不可靠资源的研究,更加接近于实际自动制造系统的生产情况;其次将原先对于含有不可靠资源的自动制造系统停留在文字叙述提升到了数学表达式上,为相应算法的提出实现了可能。自动制造系统采用的Petri建立模型,一个加工进程因不可靠资源损坏,需要检测某一个状态的稳健性,即不可靠资源损坏是否影响其他加工进程被阻塞。
本发明对不可靠资源的损坏实现了数学表达式的定义,使得不可靠资源损坏有了清晰的数学逻辑;本发明提出了一种基于数学规划的算法,不用生成整个Petri网的全局可达图,再基于全局可达图诊断状态的稳健性,而是非常巧妙地诊断给定状态的稳健性,避免了生成全局可达图产生的状态爆炸问题。
附图说明
图1是本发明实施例提供的基于数学规划算法的自动制造系统状态稳健性检测方法流程图。
图2是本发明实施例提供的含有不可靠资源的S3PR网模型示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
针对现有技术存在的问题,本发明提供了一种基于数学规划算法的自动制造系统状态稳健性检测方法,下面结合附图对本发明作详细的描述。
如图1所示,本发明实施例提供的基于数学规划算法的自动制造系统状态稳健性检测方法包括以下步骤:
S101:将含有不可靠资源的进程的变迁划分为两类:一是不可靠资源的前序变迁,二是不可靠资源的后序变迁;
S102:对不可靠资源损坏进行文字叙述到数学表达式的转化;
S103:列出资源损坏下的数学规划算法(MP)表达式,找出使能变迁发射序列;
S104:通过验证上述变迁发射序列是否满足不使用不可靠资源的其他进程不被阻塞,得到当前状态的稳健性。
下面结合具体实施例对本发明的技术方案作进一步的描述。
首先介绍一下关于Petri网的一些定义:
定义1:一个简单顺序过程(S2P)是Petri网N=(PA∪{p0},T,F),这里:
(3)N是强连通的状态机;
(4)N的每一条回路包含库所p0。
定义2:含有资源的简单顺序过程(S2PR)是Petri网N=(PA∪{p0}∪PR,T,F),使得:
(1)由生成的子网是S2P;
定义3:一个S3PR可递归定义如下:
(1)一个S2PR是一个S3PR;
(4)T=T1∪T2,F=F1∪F2。
定义4:在自动制造系统中,根据资源是否会发生故障,将所有资源分为可靠资源和不可靠资源;进一步,对于一个加工进程,根据是否使用不可靠资源,将所有加工进程区分为两类,一类是需要使用不可靠资源的进程,一类是仅使用可靠资源的进程。
本发明实施例提供的基于数学规划算法的自动制造系统状态稳健性检测方法具体包括以下步骤:
步骤一:划分变迁类,对于使用不可靠资源的加工进程,按照加工阶段的次序,把该进程中的所有变迁划分为不可靠资源的前序变迁和不可靠资源的后序变迁;具体实现过程如下:
(1)输入该S3PR网中含有不可靠资源的进程的所有变迁T={ti,i∈9,10,11,12,13,14}、库所P={pi,i=8,9,10,11,12,13}以及不可靠资源PRu={p17};
(3)按照该进程工件的加工顺序,即从当前进程的闲置库所p0=p8开始,托肯的流动方向,以不可靠资源为分界线,从闲置库所的后置变迁到不可靠资源的后置变迁/>其中包含的变迁均为不可靠资源的前序变迁,依次加入前序变迁集合,即Tfront={t9,t10};从不可靠资源的前置变迁·p17到闲置库所的前置变迁·p8,其中所包含的变迁均为不可靠资源的后序变迁,依次加入后续变迁集合,即Tback={t11,t12,t13,t14}。
步骤二:从数学上定义不可靠资源被损坏,得到不可靠资源被损坏的数学意义,即数学表达式表述;具体实现过程如下:
(1)不可靠资源被损坏的定义:工作站中的处理器无法正常工作,而缓冲区仍将正常存放工件,且资源损坏不会影响正在被加工的工件;
(2)数学上定义不可靠资源损坏为:不可靠资源被损坏,则不可靠资源的后序变迁集合都将禁止发射,对于该S3PR网,即|Tback={t11,t12,t13,t14}|=[0,0,0,0],对于不可靠资源的前序变迁集合以及其他进程的变迁则可以任意发射,即TnotBack≠[0]。
步骤三:以不可靠资源被损坏的数学表达式作为约束条件,得出含有不可靠资源的数学规划算法,在当前状态下找出满足条件的变迁发射序列;具体实现过程如下:
(1)输入任意一个状态M=[6,0,0,0,0,0,0,4,1,1,0,0,0,1,1,1,0,1,0]T,步骤一所得后序变迁集合Tback={t11,t12,t13,t14},所有库所,变迁集合以及整个Petri网的关联矩阵[N];
(3)计算满足该数学规划算法表达式:
步骤四:在当前状态Mcurrent=M=[6,0,0,0,0,0,0,4,1,1,0,0,0,1,1,1,0,1,0]T下,生成各个不使用不可靠资源的进程的可达图,即局部网的可达图,在各个局部网的可达图中寻找是否存在该变迁发射序列,使得从当前状态发射该变迁序列又回到当前状态,从而诊断出当前状态的稳健性。具体实现过程如下:
(1)输入步骤三所得变迁序列的子序列各个不使用不可靠资源的其他进程的局部网数据,对于该S3PR网,包含一个不含有不可靠资源的加工进程,该局部网初始状态M'0=[6,0,0,0,0,0,0,1,1,1,1,0]T,前置集,后置集;
(2)初始化当前状态Mcurrent=M'0=[6,0,0,0,0,0,0,1,1,1,1,0]T,标记为flag=false,初始状态已访问,则标记为true;
(3)若系统中还存在没有被搜索到的状态,则继续执行下面的过程,否则终止;
(4)选择一个状态为“false”,即没有被访问过的标识M;
1)若已搜索过M,则将其标记为“true”,并开始搜索其它“false”标识;
2)若在M下,不存在可以被使能的变迁,则将M标记为“deadlock”;
(5)对于M下所有满足使能条件的变迁t,执行下面的操作:
1)激发t得到新的标识M';
2)从M到M'添加弧t;
(6)将M的状态标记为“true”后回到(2);
(7)生成不使用不可靠资源的进程的局部网可达图PRG1;
本发明中符号说明:S3PR表示带有资源的简单顺序过程系统,M0表示初始状态,Mcurrent表示当前状态,p表示库所,t表示变迁,deadlock表示死锁状态,PRG表示局部网的可达图。
该S3PR网结构如下表所示:
其中p17为不可靠资源,则Tfront={t9,t10},Tback={t11,t12,t13,t14},禁止后续变迁集发射,即|Tback={t11,t12,t13,t14}|=[0,0,0,0]。假设需要诊断的状态为M1=[6,0,0,0,0,0,0,4,1,1,0,0,0,1,1,1,0,1,0]T,M2=[6,0,0,0,0,0,0,5,0,1,0,0,0,1,1,1,0,1,1]T,由数学规划算法表达式:
其中:
(1)局部网初始状态M'0=[6,0,0,0,0,0,0,1,1,1,1,0]T,生成局部网可达图如下所示:
Initial State[6 0 0 0 0 0 0 1 1 1 1 0]
State nr:1
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:6 0 0 0 0 0 0 1 1 1 1 0
==t1==>s2
State nr:2
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 1 0 0 0 0 0 1 1 1 0 0
==t3==>s3
==t2==>s4
State nr:3
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 0 0 1 0 0 0 0 1 1 1 0
==t4==>s5
==t1==>s6
State nr:4
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 01 0 0 0 0 1 0 1 1 0
==t5==>s5
==t1==>s7
State nr:5
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 0 0 0 1 0 0 1 1 1 0 0
==t6==>s8
State nr:6
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 1 0 1 0 0 0 0 1 1 0 0
==t2==>s9
State nr:7
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 1 1 0 0 0 0 1 0 1 0 0
==t3==>s9
State nr:8
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 0 0 0 0 1 0 1 1 0 1 0
==t1==>s10
State nr:9
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 1 1 0 0 0 0 0 1 1 0
==t5==>s11
==t4==>s12
==t1==>s13
State nr:10
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 1 0 0 0 1 0 1 1 0 0 0
==t3==>s14
==t2==>s15
State nr:11
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 1 1 0 0 0 1 1 0 0
==t6==>s14
State nr:12
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 1 0 1 0 0 1 0 1 0 0
==t6==>s15
State nr:13
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 1 1 0 0 0 0 0 1 0 0
Deadlock States
State nr:14
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 1 0 1 0 0 1 0 1 0
==t4==>s16
==t1==>s17
State nr:15
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 1 0 0 1 0 1 0 0 1 0
==t5==>s16
==t1==>s18
State nr:16
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 0 1 1 0 1 1 0 0 0
Deadlock States
State nr:17
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 0 1 0 1 0 0 1 0 0 0
==t2==>s19
State nr:18
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 1 0 0 1 0 1 0 0 0 0
==t3==>s19
State nr:19
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 1 1 0 1 0 0 0 0 1 0
==t5==>s20
==t4==>s21
==t1==>s22
State nr:20
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 0 1 1 1 0 0 1 0 0 0
Deadlock States
State nr:21
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 1 0 1 1 0 1 0 0 0 0
Deadlock States
State nr:22
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:2 1 1 1 0 1 0 0 0 0 0 0
Deadlock States
Total states count:22
由可达图可知,首先由于p12中没有托肯,到达状态s8无法发射t7,紧接着发射t8,从而不能回到初始状态;其次同样由于p12中没有托肯,到达状态s8无法发射t7,紧接着发射t8,从而不能回到初始状态。故序列/> 均不能满足,可以诊断M1=[6,0,0,0,0,0,0,4,1,1,0,0,0,1,1,1,0,1,0]T为不稳健状态。
(2)局部网初始状态M'0=[6,0,0,0,0,0,0,1,1,1,1,1]T,生成局部网可达图如下所示:
Initial State[6 0 0 0 0 0 0 1 1 1 1 1]
State nr:1
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:6 0 0 0 0 0 0 1 1 1 1 1
==t1==>s2
State nr:2
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 1 0 0 0 0 0 1 1 1 0 1
==t3==>s3
==t2==>s4
State nr:3
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 0 0 1 0 0 0 0 1 1 1 1
==t4==>s5
==t1==>s6
State nr:4
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 0 1 0 0 0 0 1 0 1 1 1
==t5==>s5
==t1==>s7
State nr:5
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 0 0 0 1 0 0 1 1 1 0 1
==t6==>s8
State nr:6
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 1 0 1 0 0 0 0 1 1 0 1
==t2==>s9
State nr:7
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 1 1 0 0 0 0 1 0 1 0 1
==t3==>s9
State nr:8
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 0 0 0 0 1 0 1 1 0 1 1
==t7==>s10
==t1==>s11
State nr:9
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 1 1 0 0 0 0 0 1 1 1
==t5==>s12
==t4==>s13
==t1==>s14
State nr:10
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:5 0 0 0 0 0 1 1 1 1 1 0
==t8==>s1
==t1==>s15
State nr:11
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 1 0 0 0 1 0 1 1 0 0 1
==t7==>s15
==t3==>s16
==t2==>s17
State nr:12
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 1 1 0 0 0 1 1 0 1
==t6==>s16
State nr:13
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 1 0 1 0 0 1 0 1 0 1
==t6==>s17
State nr:14
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 1 1 0 0 0 0 0 1 0 1
Deadlock States
State nr:15
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 1 0 0 0 0 1 1 1 1 0 0
==t8==>s2
==t3==>s18
==t2==>s19
State nr:16
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 1 0 1 0 0 1 0 1 1
==t7==>s18
==t4==>s20
==t1==>s21
State nr:17
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 1 0 0 1 0 1 0 0 1 1
==t7==>s19
==t5==>s20
==t1==>s22
State nr:18
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 1 0 0 1 0 1 1 1 0
==t8==>s3
==t4==>s23
==t1==>s24
State nr:19
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 1 0 0 0 1 1 0 1 1 0
==t8==>s4
==t5==>s23
==t1==>s25
State nr:20
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 0 1 1 0 1 1 0 0 1
==t7==>s23
State nr:21
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 0 1 0 1 0 0 1 0 0 1
==t7==>s24
==t2==>s26
State nr:22
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 1 0 0 1 0 1 0 0 0 1
==t7==>s25
==t3==>s26
State nr:23
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 0 1 0 1 1 1 1 0 0
==t8==>s5
==t6==>s27
State nr:24
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 0 1 0 0 1 0 1 1 0 0
==t8==>s6
==t2==>s28
State nr:25
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 1 0 0 0 1 1 0 1 0 0
==t8==>s7
==t3==>s28
State nr:26
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 1 1 0 1 0 0 0 0 1 1
==t7==>s28
==t5==>s29
==t4==>s30
==t1==>s31
State nr:27
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:4 0 0 0 0 1 1 1 1 0 1 0
==t8==>s8
==t1==>s32
State nr:28
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 1 1 0 0 1 0 0 1 1 0
==t8==>s9
==t5==>s33
==t4==>s34
==t1==>s35
State nr:29
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 0 1 1 1 0 0 1 0 0 1
==t7==>s33
State nr:30
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 1 0 1 1 0 1 0 0 0 1
==t7==>s34
State nr:31
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:2 1 1 1 0 1 0 0 0 0 0 1
==t7==>s35
State nr:32
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 1 0 0 0 1 1 1 1 0 0 0
==t8==>s11
==t3==>s36
==t2==>s37
State nr:33
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 0 1 1 0 1 0 1 1 0 0
==t8==>s12
==t6==>s36
State nr:34
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 1 0 1 0 1 1 0 1 0 0
==t8==>s13
==t6==>s37
State nr:35
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:2 1 1 1 0 0 1 0 0 1 0 0
==t8==>s14
State nr:36
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 0 1 0 1 1 0 1 0 1 0
==t8==>s16
==t4==>s38
==t1==>s39
State nr:37
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 1 0 0 1 1 1 0 0 1 0
==t8==>s17
==t5==>s38
==t1==>s40
State nr:38
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:3 0 0 0 1 1 1 1 1 0 0 0
==t8==>s20
State nr:39
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:2 1 0 1 0 1 1 0 1 0 0 0
==t8==>s21
==t2==>s41
State nr:40
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:2 1 1 0 0 1 1 1 0 0 0 0
==t8==>s22
==t3==>s41
State nr:41
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:2 0 1 1 0 1 1 0 0 0 1 0
==t8==>s26
==t5==>s42
==t4==>s43
==t1==>s44
State nr:42
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:2 0 0 1 1 1 1 0 1 0 00
==t8==>s29
State nr:43
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:2 0 1 0 1 1 1 1 0 0 0 0
==t8==>s30
State nr:44
p.nr:1 2 3 4 5 6 7 8 9 10 11 12
toks:1 1 1 1 0 1 1 0 0 0 0 0
==t8==>s31
Total states count:44
下面进一步给出具体实施例:
其中p14为不可靠资源,则Tfront={t1},Tback={t2,t3,t4,t5,t6},禁止后续变迁集发射,即|Tback={t2,t3,t4,t5,t6}|=[0,0,0,0,0]。任意给定一个状态M=[7,0,0,0,1,0,6,2,0,0,2,0,0,1,1]T,由数学规划算法,可得到所有发射变迁序列则变迁序列的子序列/>
不使用不可靠资源的进程(局部网)初始状态M'0=[6,2,0,0,2,0,0]T,即任意给定状态的局部网的托肯分布,生成局部网可达图如下所示:
Initial State[6 2 0 0 2 0 0]
State nr:1
p.nr:1 2 3 4 5 6 7
toks:6 2 0 0 2 0 0
Deadlock States
Total states count:1
从局部网可达图可知,其初始状态就为死锁状态,没有使能变迁可以发射,故不存在满足的发射变迁使得局部网,即不使用不可靠资源的进程,从初始状态经过变迁序列发射回到初始状态,故可以诊断任意给定的状态M=[7,0,0,0,1,0,6,2,0,0,2,0,0,1,1]T为不稳健状态。
应当注意,本发明的实施方式可以通过硬件、软件或者软件和硬件的结合来实现。硬件部分可以利用专用逻辑来实现;软件部分可以存储在存储器中,由适当的指令执行系统,例如微处理器或者专用设计硬件来执行。本领域的普通技术人员可以理解上述的设备和方法,可以使用计算机可执行指令或包含在处理器控制代码中来实现,例如在诸如磁盘、CD或DVD-ROM的载体介质、诸如只读存储器(固件)的可编程的存储器或者诸如光学或电子信号载体的数据载体上提供了这样的代码。本发明的设备及其模块可以由诸如超大规模集成电路或门阵列、诸如逻辑芯片、晶体管等的半导体、或者诸如现场可编程门阵列、可编程逻辑设备等的可编程硬件设备的硬件电路实现,也可以用由各种类型的处理器执行的软件实现,也可以由上述硬件电路和软件的结合例如固件来实现。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于数学规划算法的自动制造系统状态稳健性检测方法,其特征在于,所述基于数学规划算法的自动制造系统状态稳健性检测方法包括:
第一步,将含有不可靠资源的进程的变迁划分为两类:一是不可靠资源的前序变迁,二是不可靠资源的后序变迁;
第二步,对不可靠资源损坏进行文字叙述到数学表达式的转化;
第三步,列出资源损坏下的数学规划算法表达式,找出使能变迁发射序列;
第四步,通过验证上述变迁发射序列是否满足不使用不可靠资源的其他进程不被阻塞,得到当前状态的稳健性;
所述第一步还包括:划分变迁类,对于使用不可靠资源的加工进程,按照加工阶段的次序,把该进程中的所有变迁划分为不可靠资源的前序变迁和不可靠资源的后序变迁;
(1)输入Petri网中含有不可靠资源的进程的所有变迁、库所以及不可靠资源;
(3)按照该进程工件的加工顺序,即从当前进程的闲置库所P0开始,托肯的流动方向,以不可靠资源为分界线,从闲置库所的后置变迁{P0}·到不可靠资源的后置变迁其中包含的变迁均为不可靠资源的前序变迁,依次加入前序变迁集合;从不可靠资源的前置变迁到闲置库所的前置变迁·{P0},其中所包含的变迁均为不可靠资源的后序变迁,依次加入后序变迁集合;
所述第四步还包括:在当前状态下,生成各个不使用不可靠资源的进程的可达图,即局部网的可达图,在各个局部网的可达图中寻找是否存在发射该变迁发射序列,使得从当前状态发射该变迁序列又回到当前状态,诊断出当前状态的稳健性;
(2)初始化当前状态Mcurrent=M'0,标记为flag=false,初始状态已访问,则标记为true;
(3)若系统中还存在没有被搜索到的状态,则继续执行下面的过程,否则终止;
(4)选择一个状态为“false”,即没有被访问过的标识M;
(5)对于M下所有满足使能条件的变迁t,执行下面的操作:
1)激发t得到新的标识M';
2)从M到M'添加弧t;
(6)将M的状态标记为“true”后回到(2);
(7)生成各个不使用不可靠资源的进程的局部网可达图;
所述选择一个状态为“false”,即没有被访问过的标识M包括:
1)若已搜索过M,则将其标记为“true”,并开始搜索其它“false”标识;
2)若在M下,不存在可以被使能的变迁,则将M标记为“deadlock”。
2.如权利要求1所述的基于数学规划算法的自动制造系统状态稳健性检测方法,其特征在于,所述第二步还包括:从数学上定义不可靠资源被损坏,得到不可靠资源被损坏的数学意义,即数学表达式表述;
(1)不可靠资源被损坏的定义为:工作站中的处理器无法正常工作,而缓冲区仍将正常存放工件,且资源损坏不会影响正在被加工的工件;如果在缓冲区中存在没被使用的空闲位置,那么已经完成上一阶段加工任务的工件就进入到该缓冲区;另外对于已经完成当前操作的工件,则可以正常进入到下一个阶段的加工中,但对于没有加工完成或者正在等待资源的工件,将一直保持请求状态,直到该资源被修复;
(2)数学上定义不可靠资源损坏为:不可靠资源被损坏,则不可靠资源的后序变迁集合都将禁止发射,即|Tback|=[0],对于不可靠资源的前序变迁集合以及其他进程的变迁则可以任意发射。
3.如权利要求1所述的基于数学规划算法的自动制造系统状态稳健性检测方法,其特征在于,所述第三步还包括:以不可靠资源被损坏的数学表达式作为约束条件,得出含有不可靠资源的数学规划算法,在当前状态下找出满足条件的变迁发射序列;
(1)输入任意一个状态M,步骤一所得后序变迁集合Tback,所有库所,变迁集合以及整个Petri网的关联矩阵[N];
(3)数学规划算法表达式:
式(1)在当前状态下,存在一个变迁序列,使得系统回到又当前状态;式(2)变迁序列满足不可靠资源的后序变迁均禁止发射;式(3)变迁序列至少存在一个变迁发射,即当前状态不能静止不动。
4.一种接收用户输入程序存储介质,其特征在于,所存储的计算机程序使电子设备执行权利要求1~3任意一项所述方法的步骤。
5.一种如权利要求1~3任意一项所述基于数学规划算法的自动制造系统状态稳健性检测方法在自动制造系统中的应用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010051204.3A CN111352397B (zh) | 2020-01-17 | 2020-01-17 | 基于数学规划算法的自动制造系统状态稳健性检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010051204.3A CN111352397B (zh) | 2020-01-17 | 2020-01-17 | 基于数学规划算法的自动制造系统状态稳健性检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111352397A CN111352397A (zh) | 2020-06-30 |
CN111352397B true CN111352397B (zh) | 2023-03-31 |
Family
ID=71195799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010051204.3A Active CN111352397B (zh) | 2020-01-17 | 2020-01-17 | 基于数学规划算法的自动制造系统状态稳健性检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111352397B (zh) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7120699B2 (en) * | 2001-09-20 | 2006-10-10 | Ricoh Company, Ltd. | Document controlled workflow systems and methods |
US7627459B2 (en) * | 2004-01-08 | 2009-12-01 | Georgia Tech Research Corporation | Systems and methods for reliability and performability assessment |
CN102622510A (zh) * | 2012-01-31 | 2012-08-01 | 龚波 | 一种软件缺陷量化管理系统和方法 |
CN106200575B (zh) * | 2016-07-07 | 2018-07-20 | 西安电子科技大学 | 一种基于Petri网的自动制造系统的稳健性控制方法 |
CN108919645B (zh) * | 2018-07-09 | 2021-05-11 | 西安电子科技大学 | 一种存在不可控行为的自动制造系统的稳健性控制方法 |
CN108919644B (zh) * | 2018-07-09 | 2021-05-11 | 西安电子科技大学 | 存在不可观行为的自动制造系统的稳健性控制方法 |
CN110209118B (zh) * | 2019-05-22 | 2022-04-01 | 西安电子科技大学 | 具有使用多种类多数量资源和柔性路径特征的基于Petri网的稳健性控制方法 |
-
2020
- 2020-01-17 CN CN202010051204.3A patent/CN111352397B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111352397A (zh) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112740184B (zh) | 用于确定性地报告软件系统中的原因和影响的方法 | |
US9336123B1 (en) | Method and system for automatically establishing a component description format (CDF) debugging environment | |
US10241852B2 (en) | Automated qualification of a safety critical system | |
JP4404438B2 (ja) | 合成及びレイアウトのオーバーヘッドが少ないsatソルバー・アーキテクチャーを得るための方法及び装置 | |
CN110489812B (zh) | 多层层级网表处理方法、装置、计算机设备及存储介质 | |
Reveliotis | Accommodating FMS operational contingencies through routing flexibility | |
WO2014200551A1 (en) | Identifying the introduction of a software failure | |
Jannach et al. | Parallelized hitting set computation for model-based diagnosis | |
Hsiao et al. | Asyncclock: Scalable inference of asynchronous event causality | |
CN111352397B (zh) | 基于数学规划算法的自动制造系统状态稳健性检测方法 | |
US20240045787A1 (en) | Code inspection method under weak memory ordering architecture and corresponding device | |
US10007746B1 (en) | Method and system for generalized next-state-directed constrained random simulation | |
CN112631925A (zh) | 一种单变量原子违背缺陷的检测方法 | |
CN111124694A (zh) | 一种基于petri网的可达图的死锁检测和解决方法 | |
Bowie | Applications of graph theory in computer systems | |
Elsayed et al. | An optimal transition-based recovery policy for controlling deadlock within flexible manufacturing systems using graph technique | |
US7260803B2 (en) | Incremental dummy metal insertions | |
JP6037976B2 (ja) | ソースコード生成装置 | |
US8627301B2 (en) | Concurrent management of adaptive programs | |
Khalaf et al. | Revisiting the behavior of fault and compensation handlers in WS-BPEL | |
CN104145249A (zh) | 调试方法和计算机程序产品 | |
Yeung et al. | Modelling and model checking suspendible business processes via statechart diagrams and CSP | |
Kaiser et al. | BDD complexity reduction by component fault trees | |
CN117907812B (zh) | 电路检测方法及装置、电子设备、存储介质、程序产品 | |
US20170373925A1 (en) | Analyzing contention data and following resource blockers to find root causes of computer problems |
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 |