CN111352397B - 基于数学规划算法的自动制造系统状态稳健性检测方法 - Google Patents

基于数学规划算法的自动制造系统状态稳健性检测方法 Download PDF

Info

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
Application number
CN202010051204.3A
Other languages
English (en)
Other versions
CN111352397A (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN202010051204.3A priority Critical patent/CN111352397B/zh
Publication of CN111352397A publication Critical patent/CN111352397A/zh
Application granted granted Critical
Publication of CN111352397B publication Critical patent/CN111352397B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/418Total 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/41885Total 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
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/32Operator till task planning
    • G05B2219/32339Object oriented modeling, design, analysis, implementation, simulation language
    • 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/02Total 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网中含有不可靠资源的进程的所有变迁、库所以及不可靠资源;
(2)初始化数据,前序变迁集合
Figure GDA0004027178420000031
后序变迁集合/>
Figure GDA0004027178420000032
(3)按照该进程工件的加工顺序,即从当前进程的闲置库所P0开始,托肯的流动方向,以不可靠资源为分界线,从闲置库所的后置变迁{P0}·到不可靠资源的后置变迁{PRu}·,其中包含的变迁均为不可靠资源的前序变迁,依次加入前序变迁集合;从不可靠资源的前置变迁·{PRu}到闲置库所的前置变迁·{P0},其中所包含的变迁均为不可靠资源的后序变迁,依次加入后续变迁集合。
进一步,所述第二步还包括:从数学上定义不可靠资源被损坏,得到不可靠资源被损坏的数学意义,即数学表达式表述;
(1)不可靠资源被损坏的定义为:工作站中的处理器无法正常工作,而缓冲区仍将正常存放工件,且资源损坏不会影响正在被加工的工件;如果在缓冲区中存在没被使用的空闲位置,那么已经完成上一阶段加工任务的工件就进入到该缓冲区;另外对于已经完成当前操作的工件,则可以正常进入到下一个阶段的加工中,但对于没有加工完成或者正在等待资源的工件,将一直保持请求状态,直到该资源被修复;
(4)数学上定义不可靠资源损坏为:不可靠资源被损坏,则不可靠资源的后序变迁集合都将禁止发射,即|Tback|=[0],对于不可靠资源的前序变迁集合以及其他进程的变迁则可以任意发射。
进一步,所述第三步还包括:以不可靠资源被损坏的数学表达式作为约束条件,得出含有不可靠资源的数学规划算法,在当前状态下找出满足条件的变迁发射序列;
(1)输入任意一个状态M,步骤一所得后序变迁集合Tback,所有库所,变迁集合以及整个Petri网的关联矩阵[N];
(2)初始化当前状态Mcurrent=M,发射变迁序列
Figure GDA0004027178420000041
/>
(3)数学规划算法表达式:
Figure GDA0004027178420000042
Figure GDA0004027178420000043
Figure GDA0004027178420000044
式(1)在当前状态下,存在一个变迁序列,使得系统回到又当前状态;式(2)变迁序列满足不可靠资源的后序变迁均禁止发射;式(3)变迁序列至少存在一个变迁发射,即当前状态不能静止不动。
进一步,所述第四步还包括:在当前状态下,生成各个不使用不可靠资源的进程的可达图,即局部网的可达图,在各个局部网的可达图中寻找是否存在该变迁发射序列,使得从当前状态发射该变迁序列又回到当前状态,诊断出当前状态的稳健性;
(1)输入第三步所得变迁序列的子序列
Figure GDA0004027178420000045
各个不使用不可靠资源的其他进程的局部网数据,包含局部网初始状态M'0,即需要判断的状态下各进程中托肯的分布,前置集,后置集;
(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)生成各个不使用不可靠资源的进程的局部网可达图;
(8)根据第三步所得变迁序列的子序列
Figure GDA0004027178420000051
找出各可达图中均是否存在该发射变迁序列,使得回到M'0,若均不存在,则需要诊断的状态为不稳健的。
进一步,所述选择一个状态为“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),这里:
(1)
Figure GDA00040271784200000712
称为工序库所的集合;
(2)
Figure GDA00040271784200000713
称为闲置进程库所或简称闲置库所;
(3)N是强连通的状态机;
(4)N的每一条回路包含库所p0
定义2:含有资源的简单顺序过程(S2PR)是Petri网N=(PA∪{p0}∪PR,T,F),使得:
(1)由生成的子网是S2P;
Figure GDA0004027178420000071
Figure GDA0004027178420000072
Figure GDA0004027178420000073
Figure GDA0004027178420000074
定义3:一个S3PR可递归定义如下:
(1)一个S2PR是一个S3PR;
(2)令
Figure GDA0004027178420000075
是两个S3PR,使得
Figure GDA0004027178420000076
那么,由N1和N2通过PC复合而成的Petri网/>
Figure GDA0004027178420000077
仍然是S3PR,定义为:
Figure GDA0004027178420000078
Figure GDA0004027178420000079
Figure GDA00040271784200000710
(4)T=T1∪T2,F=F1∪F2
由n个S2PR网N1,N2,…,Nn复合而成的S3PR网N记为
Figure GDA00040271784200000711
定义4:在自动制造系统中,根据资源是否会发生故障,将所有资源分为可靠资源和不可靠资源;进一步,对于一个加工进程,根据是否使用不可靠资源,将所有加工进程区分为两类,一类是需要使用不可靠资源的进程,一类是仅使用可靠资源的进程。
本发明实施例提供的基于数学规划算法的自动制造系统状态稳健性检测方法具体包括以下步骤:
步骤一:划分变迁类,对于使用不可靠资源的加工进程,按照加工阶段的次序,把该进程中的所有变迁划分为不可靠资源的前序变迁和不可靠资源的后序变迁;具体实现过程如下:
(1)输入该S3PR网中含有不可靠资源的进程的所有变迁T={ti,i∈9,10,11,12,13,14}、库所P={pi,i=8,9,10,11,12,13}以及不可靠资源PRu={p17};
(2)初始化数据,前序变迁集合
Figure GDA0004027178420000081
后序变迁集合/>
Figure GDA0004027178420000082
(3)按照该进程工件的加工顺序,即从当前进程的闲置库所p0=p8开始,托肯的流动方向,以不可靠资源为分界线,从闲置库所的后置变迁
Figure GDA0004027178420000083
到不可靠资源的后置变迁/>
Figure GDA0004027178420000084
其中包含的变迁均为不可靠资源的前序变迁,依次加入前序变迁集合,即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];
(2)初始化当前状态Mcurrent=M=[6,0,0,0,0,0,0,4,1,1,0,0,0,1,1,1,0,1,0]T,发射变迁序列
Figure GDA0004027178420000091
(3)计算满足该数学规划算法表达式:
Figure GDA0004027178420000092
Figure GDA0004027178420000093
Figure GDA0004027178420000094
得到所有发射变迁序列
Figure GDA0004027178420000095
步骤四:在当前状态Mcurrent=M=[6,0,0,0,0,0,0,4,1,1,0,0,0,1,1,1,0,1,0]T下,生成各个不使用不可靠资源的进程的可达图,即局部网的可达图,在各个局部网的可达图中寻找是否存在该变迁发射序列,使得从当前状态发射该变迁序列又回到当前状态,从而诊断出当前状态的稳健性。具体实现过程如下:
(1)输入步骤三所得变迁序列的子序列
Figure GDA0004027178420000096
各个不使用不可靠资源的其他进程的局部网数据,对于该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;
(8)根据步骤三所得变迁序列的子序列
Figure GDA0004027178420000102
找出可达图PRG1中是否存在该发射变迁序列,使得回到M'0=[6,0,0,0,0,0,0,1,1,1,1,0]T,若不存在,则需要诊断的状态为不稳健的。
本发明中符号说明:S3PR表示带有资源的简单顺序过程系统,M0表示初始状态,Mcurrent表示当前状态,p表示库所,t表示变迁,deadlock表示死锁状态,PRG表示局部网的可达图。
该S3PR网结构如下表所示:
Figure GDA0004027178420000101
/>
Figure GDA0004027178420000111
其中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,由数学规划算法表达式:
Figure GDA0004027178420000112
Figure GDA0004027178420000113
Figure GDA0004027178420000114
其中:
Figure GDA0004027178420000121
均可得到所有发射变迁序列
Figure GDA0004027178420000122
Figure GDA0004027178420000123
实质上/>
Figure GDA0004027178420000124
和/>
Figure GDA0004027178420000125
为该Petri网的T不变式。
则变迁序列的子序列
Figure GDA0004027178420000126
(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
由可达图可知,首先
Figure GDA0004027178420000161
由于p12中没有托肯,到达状态s8无法发射t7,紧接着发射t8,从而不能回到初始状态;其次
Figure GDA0004027178420000162
同样由于p12中没有托肯,到达状态s8无法发射t7,紧接着发射t8,从而不能回到初始状态。故序列/>
Figure GDA0004027178420000163
Figure GDA0004027178420000164
均不能满足,可以诊断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
由可达图可知,首先
Figure GDA0004027178420000241
从而能够找到一条回路回到初始状态;其次
Figure GDA0004027178420000242
同样能回到初始状态。故序列
Figure GDA0004027178420000243
均能满足,可以诊断M1=[6,0,0,0,0,0,0,5,0,1,0,0,0,1,1,1,0,1,1]T为稳健状态。
下面进一步给出具体实施例:
Figure GDA0004027178420000244
/>
Figure GDA0004027178420000251
其中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,由数学规划算法,可得到所有发射变迁序列
Figure GDA0004027178420000252
则变迁序列的子序列/>
Figure GDA0004027178420000253
不使用不可靠资源的进程(局部网)初始状态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
从局部网可达图可知,其初始状态就为死锁状态,没有使能变迁可以发射,故不存在满足
Figure GDA0004027178420000254
的发射变迁使得局部网,即不使用不可靠资源的进程,从初始状态经过变迁序列发射回到初始状态,故可以诊断任意给定的状态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网中含有不可靠资源的进程的所有变迁、库所以及不可靠资源;
(2)初始化数据,前序变迁集合
Figure FDA0004083607790000011
后序变迁集合/>
Figure FDA0004083607790000012
(3)按照该进程工件的加工顺序,即从当前进程的闲置库所P0开始,托肯的流动方向,以不可靠资源为分界线,从闲置库所的后置变迁{P0}·到不可靠资源的后置变迁
Figure FDA0004083607790000013
其中包含的变迁均为不可靠资源的前序变迁,依次加入前序变迁集合;从不可靠资源的前置变迁
Figure FDA0004083607790000014
到闲置库所的前置变迁·{P0},其中所包含的变迁均为不可靠资源的后序变迁,依次加入后序变迁集合;
所述第四步还包括:在当前状态下,生成各个不使用不可靠资源的进程的可达图,即局部网的可达图,在各个局部网的可达图中寻找是否存在发射该变迁发射序列,使得从当前状态发射该变迁序列又回到当前状态,诊断出当前状态的稳健性;
(1)输入第三步所得变迁序列的子序列
Figure FDA0004083607790000015
各个不使用不可靠资源的其他进程的局部网数据,包含局部网初始状态M'0,即需要判断的状态下各进程中托肯的分布,前置集,后置集;
(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)生成各个不使用不可靠资源的进程的局部网可达图;
(8)根据第三步所得变迁序列的子序列
Figure FDA0004083607790000021
找出各可达图中均是否存在该变迁发射序列,使得回到M'0,若均不存在,则需要诊断的状态为不稳健的;
所述选择一个状态为“false”,即没有被访问过的标识M包括:
1)若已搜索过M,则将其标记为“true”,并开始搜索其它“false”标识;
2)若在M下,不存在可以被使能的变迁,则将M标记为“deadlock”。
2.如权利要求1所述的基于数学规划算法的自动制造系统状态稳健性检测方法,其特征在于,所述第二步还包括:从数学上定义不可靠资源被损坏,得到不可靠资源被损坏的数学意义,即数学表达式表述;
(1)不可靠资源被损坏的定义为:工作站中的处理器无法正常工作,而缓冲区仍将正常存放工件,且资源损坏不会影响正在被加工的工件;如果在缓冲区中存在没被使用的空闲位置,那么已经完成上一阶段加工任务的工件就进入到该缓冲区;另外对于已经完成当前操作的工件,则可以正常进入到下一个阶段的加工中,但对于没有加工完成或者正在等待资源的工件,将一直保持请求状态,直到该资源被修复;
(2)数学上定义不可靠资源损坏为:不可靠资源被损坏,则不可靠资源的后序变迁集合都将禁止发射,即|Tback|=[0],对于不可靠资源的前序变迁集合以及其他进程的变迁则可以任意发射。
3.如权利要求1所述的基于数学规划算法的自动制造系统状态稳健性检测方法,其特征在于,所述第三步还包括:以不可靠资源被损坏的数学表达式作为约束条件,得出含有不可靠资源的数学规划算法,在当前状态下找出满足条件的变迁发射序列;
(1)输入任意一个状态M,步骤一所得后序变迁集合Tback,所有库所,变迁集合以及整个Petri网的关联矩阵[N];
(2)初始化当前状态Mcurrent=M,发射变迁序列
Figure FDA0004083607790000031
(3)数学规划算法表达式:
Figure FDA0004083607790000032
Figure FDA0004083607790000033
Figure FDA0004083607790000034
式(1)在当前状态下,存在一个变迁序列,使得系统回到又当前状态;式(2)变迁序列满足不可靠资源的后序变迁均禁止发射;式(3)变迁序列至少存在一个变迁发射,即当前状态不能静止不动。
4.一种接收用户输入程序存储介质,其特征在于,所存储的计算机程序使电子设备执行权利要求1~3任意一项所述方法的步骤。
5.一种如权利要求1~3任意一项所述基于数学规划算法的自动制造系统状态稳健性检测方法在自动制造系统中的应用。
CN202010051204.3A 2020-01-17 2020-01-17 基于数学规划算法的自动制造系统状态稳健性检测方法 Active CN111352397B (zh)

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)

* Cited by examiner, † Cited by third party
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网的稳健性控制方法

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