CN113467405B - 工业4.0需求驱动的分布式动态边缘云智能制造方法及系统 - Google Patents

工业4.0需求驱动的分布式动态边缘云智能制造方法及系统 Download PDF

Info

Publication number
CN113467405B
CN113467405B CN202110879349.7A CN202110879349A CN113467405B CN 113467405 B CN113467405 B CN 113467405B CN 202110879349 A CN202110879349 A CN 202110879349A CN 113467405 B CN113467405 B CN 113467405B
Authority
CN
China
Prior art keywords
subtasks
subtask
distributed
factory
cost
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
CN202110879349.7A
Other languages
English (en)
Other versions
CN113467405A (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.)
Fuzhou University
Original Assignee
Fuzhou 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 Fuzhou University filed Critical Fuzhou University
Priority to CN202110879349.7A priority Critical patent/CN113467405B/zh
Publication of CN113467405A publication Critical patent/CN113467405A/zh
Application granted granted Critical
Publication of CN113467405B publication Critical patent/CN113467405B/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/41875Total 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 quality surveillance of production
    • 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/32368Quality control
    • 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)
  • General Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种工业4.0需求驱动的分布式动态边缘云智能制造方法及系统。包括:客户向云端发布个性化产品需求;设计师获取需求信息,与客户沟通,确定最终设计;云端根据最终设计,将个性化产品分解为几个子任务;在决策时刻,边缘端根据云端发布的待分配子任务信息,判断当前工厂状态是否能够执行任一待分配子任务,并将结果返回云端;建立动态分布式任务分配模型;采用基于l p ‑Box ADMM的方法完成任务分配;待个性化产品的所有子任务完成,将产品打包并运输给客户。本发明以客户需求为驱动,让客户参与产品生产全生命周期并实时提出修改意见,提高了产品个性化程度,同时,采用分布式边缘云架构,减轻云端负担,保障数据隐私,避免单点故障,提高鲁棒性。

Description

工业4.0需求驱动的分布式动态边缘云智能制造方法及系统
技术领域
本发明涉及工业4.0智能制造方法及系统,更具体地,涉及一种工业4.0需求驱动的分布式动态边缘云智能制造方法及系统。
背景技术
工业4.0,也称作第四次工业革命,是由德国在2013年提出的。除此之外,美国的“先进制造伙伴计划”(AMP)、英国的“未来制造”以及其他政府计划都致力于通过实际行动逐步实现工业4.0。简单来说,工业4.0可以概括为将所有物理设备集成到网络中,旨在以高度自动化和灵活性实现大规模个性化生产的工业革命。大规模个性化生产是为了解决客户日益增长的个性化需求与个性化产品短缺之间的矛盾。大规模个性化生产以客户的需求为驱动,期望客户能够参与产品生产的全生命周期,并实时提出修改意见,提高产品个性化程度以及客户产品体验满意度。
以客户需求为驱动的大规模个性化生产是工业4.0的目标。近几年来,全球各地也提出了许多方法来实现大规模个性化生产。Yao等在考虑总完成时间、成本和绿色水平的情况下,让制造商通过执行分解的生产集来协作和安排供应商以完成个性化产品生产。Ogunsakin等提出了一种自组织柔性制造系统,其移动加工站可以基于BEEPOST算法动态调整位置以形成最佳资源布局。Leng等针对大规模个性化生产,考虑了动态生产流程管理的复杂性和产品灵活性,提出了一种基于信息物理社会系统(CPSS)的情境感知的社会物联网(SIoT)体系结构。除此之外,还有社会制造以及云制造等先进制造模式致力于实现大规模个性化生产。
虽然还有一部分研究专注于分布式架构,但要实现工业4.0需求驱动的大规模个性化生产,仍然存在一些问题。Jian等在云制造的基础上添加了边缘层,先在云端进行一次任务分解,再在边缘端将子任务分解为不可分割的原任务,最终生成原任务调度方案。Ma等提出的架构分为边缘侧和云侧,云侧将订单分解成子任务,边缘侧进行制造任务和制造资源的匹配,最终将任务分配给工厂,并完成工厂中机器的调度方案。这两种方法由于其分布式架构,不依赖于中心节点,不存在单点故障等问题,但是所提出的都是调度解决方案,缺乏实时性,不能很好的考虑机器故障以及客户实时提出产品修改意见等动态情况。
现存的用于工业4.0需求驱动的大规模个性化生产的方法普遍存在三个问题。其一,客户参与度低,比如社会制造,客户只能参与产品设计过程以及制造商选择过程,其他的生产过程都不可参与。其二,即使客户可以参与产品生产的全生命周期,现存的方法也不能快速响应客户实时提出的产品修改意见,导致产品修改错过最佳时间,增加了个性化产品成本,也降低了产品个性化程度。其三,现存的实现大规模个性化生产的方法大多是集中式架构,依赖于中心节点,数据隐私难以保障,容易造成单点故障,且随着客户和工厂数量的增加,数据处理会更加缓慢甚至超出云端计算能力。
针对上述工业4.0需求驱动的大规模个性化生产存在的问题,亟需工业4.0需求驱动的分布式动态边缘云智能制造方法及系统。
发明内容
本发明的目的在于提供一种工业4.0需求驱动的分布式动态边缘云智能制造方法及系统,该方法及系统能够提高客户产品满意度。
为实现上述目的,本发明的技术方案是:一种工业4.0需求驱动的分布式动态边缘云智能制造方法,包括以下步骤:
步骤S1、客户向云端发布个性化产品需求;
步骤S2、设计师获取需求信息,确定最终设计;
步骤S3、云端根据最终设计,将个性化产品分解为几个子任务;
步骤S4、在决策时刻,边缘端根据云端发布的待分配子任务信息,判断当前工厂状态是否能够执行任一待分配子任务,并将结果返回云端;
步骤S5、建立动态分布式任务分配模型;
步骤S6、采用基于lp-BoxADMM的方法完成任务分配;
步骤S7、待个性化产品的所有子任务完成,将产品打包并运输给客户。
在本发明一实施例中,所述步骤S5具体包括以下步骤:
步骤A1、以最大化决策时刻t所有工厂总利润为目标,建立目标函数;
步骤A2、根据场景实际情况,添加约束条件;
步骤A3、将模型转化成标准的求解最小化问题的分布式优化数学模型。
在本发明一实施例中,所述步骤A1具体实现如下:
动态分布式任务分配模型的目标是最大化决策时刻t所有工厂总利润,目标函数如下:
Figure BDA0003190989300000021
其中,i,j分别表示工厂和子任务的下标,n,m分别表示可用工厂和待分配子任务的数量,
Figure BDA0003190989300000022
表示工厂
Figure BDA0003190989300000023
执行子任务
Figure BDA0003190989300000024
能够获得的效用,
Figure BDA0003190989300000025
表示执行代价,
Figure BDA0003190989300000026
表示运输代价;
Figure BDA0003190989300000027
表示工厂
Figure BDA0003190989300000028
在决策时刻t是否有能力执行子任务
Figure BDA0003190989300000029
如果是的话,值为1,否则,值为0;
Figure BDA00031909893000000210
为待求解的决策变量,表示在决策时刻t,是否将子任务
Figure BDA0003190989300000031
分配给工厂
Figure BDA0003190989300000032
执行,是的话,值为1,否则,值为0;每个工厂的利润都可以表示为效用减去执行代价减去运输代价;
效用
Figure BDA0003190989300000033
受工厂生产能力匹配度以及工厂执行子任务所花费的时间长短的影响,执行每个子任务
Figure BDA0003190989300000034
都需要原材料
Figure BDA0003190989300000035
以及生产能力
Figure BDA0003190989300000036
其中q为需要的生产能力种类总数,生产需要的每种原材料
Figure BDA0003190989300000037
的数量表示为
Figure BDA0003190989300000038
每个工厂
Figure BDA0003190989300000039
在决策时刻所拥有的生产能力表示为
Figure BDA00031909893000000310
因此,定义
Figure BDA00031909893000000311
表示工厂
Figure BDA00031909893000000312
的第h种生产能力是否与子任务
Figure BDA00031909893000000313
执行所需要的第k种能力相同,如果是的话,值为1,否则为0;效用的计算还需要获取先前完成子任务
Figure BDA00031909893000000314
最快用时TTj以及子任务
Figure BDA00031909893000000315
最高售价pricej,并且,预估工厂
Figure BDA00031909893000000316
执行子任务
Figure BDA00031909893000000317
用时为
Figure BDA00031909893000000318
由此可得,效用
Figure BDA00031909893000000319
计算如下:
Figure BDA00031909893000000320
执行代价
Figure BDA00031909893000000321
包括机器设置成本
Figure BDA00031909893000000322
生产成本
Figure BDA00031909893000000323
以及原材料购买成本,每种原材料
Figure BDA00031909893000000324
的成本的计算包括两种情况,第一种情况是原材料
Figure BDA00031909893000000325
需要的数量
Figure BDA00031909893000000326
大于工厂库存的数量
Figure BDA00031909893000000327
那么原材料成本就是购买成本
Figure BDA00031909893000000328
加上储存成本
Figure BDA00031909893000000329
如果原材料数量足够生产,那么只需要储存成本
Figure BDA00031909893000000330
基于此,引入
Figure BDA00031909893000000331
来判断需要的原材料数量是否多于原材料库存数量,如果是的话,值为1,否则,值为0;因此,执行代价计算公式如下:
Figure BDA00031909893000000332
运输代价
Figure BDA00031909893000000333
包括子任务
Figure BDA00031909893000000334
运输到工厂
Figure BDA00031909893000000335
的代价
Figure BDA00031909893000000336
以及各种需要购买的原材料的运输代价,每种需要购买的原材料
Figure BDA00031909893000000337
的运输代价为购买数量
Figure BDA00031909893000000338
乘以单位运输代价
Figure BDA00031909893000000339
因此,运输代价
Figure BDA00031909893000000340
计算公式表示如下:
Figure BDA0003190989300000041
在本发明一实施例中,所述步骤A2具体实现如下:
动态分布式任务分配模型的约束包括:
1)
Figure BDA0003190989300000042
以及
Figure BDA0003190989300000043
的取值为0或1,表示如下:
Figure BDA0003190989300000044
Figure BDA0003190989300000045
2)子任务运输代价
Figure BDA0003190989300000046
和原材料运输代价
Figure BDA0003190989300000047
的取值如下:
Figure BDA0003190989300000048
Figure BDA0003190989300000049
运输代价取值为-1,表示没有能力运输,即子任务
Figure BDA00031909893000000410
或者购买的原材料
Figure BDA00031909893000000411
无法运输到工厂
Figure BDA00031909893000000412
运输代价为0,对于子任务
Figure BDA00031909893000000413
来说,存在三种情况,当前子任务
Figure BDA00031909893000000414
为第一个子任务、当前子任务就在工厂
Figure BDA00031909893000000415
无需运输,或者因为其他情况导致子任务运输代价为0;对于原材料来说,就是指原材料
Figure BDA00031909893000000416
运输到工厂
Figure BDA00031909893000000417
代价为0,其他情况下,运输代价为大于0的实数;
3)假定待分配子任务在决策时刻至少有一个工厂可以执行,因此,关于
Figure BDA00031909893000000418
的约束:
Figure BDA00031909893000000419
4)
Figure BDA00031909893000000420
用来表示工厂
Figure BDA00031909893000000421
在决策时刻t是否有能力执行子任务
Figure BDA00031909893000000422
如果有的话,值为1,否则,值为0,有如下约束会导致
Figure BDA00031909893000000423
Figure BDA00031909893000000424
Figure BDA00031909893000000425
Figure BDA00031909893000000426
Figure BDA00031909893000000427
这四个约束从上到下依次表示:
1)如果子任务
Figure BDA0003190989300000051
无法运输到工厂
Figure BDA0003190989300000052
Figure BDA0003190989300000053
那么是工厂
Figure BDA0003190989300000054
无法执行子任务
Figure BDA0003190989300000055
2)如果工厂
Figure BDA0003190989300000056
执行子任务
Figure BDA0003190989300000057
必须购买原材料
Figure BDA0003190989300000058
Figure BDA0003190989300000059
但是购买原材料
Figure BDA00031909893000000510
无法运输到工厂
Figure BDA00031909893000000511
那么,子任务
Figure BDA00031909893000000512
无法在工厂
Figure BDA00031909893000000513
执行;
3)如果工厂
Figure BDA00031909893000000514
能够匹配给子任务
Figure BDA00031909893000000515
的能力小于总能力种类的一半,那么,也视为工厂
Figure BDA00031909893000000516
没有能力执行子任务
Figure BDA00031909893000000517
4)如果工厂
Figure BDA00031909893000000518
所获得的利润小于0的话,那么也视为工厂
Figure BDA00031909893000000519
无法执行子任务
Figure BDA00031909893000000520
5)关于决策变量
Figure BDA00031909893000000521
有如下约束:
Figure BDA00031909893000000522
Figure BDA00031909893000000523
前者表示任何一个子任务都最多被一个工厂执行,后者表示任何一个工厂都最多只能执行一个子任务;
6)关于
Figure BDA00031909893000000524
的约束:
Figure BDA00031909893000000525
即如果工厂
Figure BDA00031909893000000526
执行子任务
Figure BDA00031909893000000527
需要购买原材料
Figure BDA00031909893000000528
那么,
Figure BDA00031909893000000529
取值为1,否则为0。
在本发明一实施例中,所述步骤A3具体实现如下:
Figure BDA00031909893000000530
定义
Figure BDA00031909893000000531
效用
Figure BDA00031909893000000532
执行代价
Figure BDA00031909893000000533
运输代价
Figure BDA00031909893000000534
以及工厂执行能力变量
Figure BDA00031909893000000535
都可在决策时刻t根据计算实时获得,因此,求解时只需要考虑与
Figure BDA00031909893000000536
相关的约束,则数学模型可表示为:
Figure BDA00031909893000000537
Figure BDA00031909893000000538
Figure BDA0003190989300000061
Figure BDA0003190989300000062
由于每个管理工厂
Figure BDA0003190989300000063
的边缘端只知道自身的目标函数fi(xi),并且也只能决定与自身相关的决策变量
Figure BDA0003190989300000064
的取值情况,因此这是一个分布式优化问题,该问题假设每个边缘端之间能够正常通信,并且能够获得来自其他边缘端刚更新过的决策变量
Figure BDA0003190989300000065
在本发明一实施例中,所述步骤S6具体包括以下步骤:
步骤B1、添加虚拟工厂或虚拟子任务,将不等式约束转化为等式约束;
步骤B2、将决策变量二值约束转化成盒约束与lp球约束的交集;
步骤B3、构造增广拉格朗日函数;
步骤B4、将问题分解为每个工厂要解决的子问题;
步骤B5、根据更新公式以及收敛条件进行分布式求解。
在本发明一实施例中,所述步骤B1具体实现如下:
令u=max{n,m},对于n<m的情况,添加m-n个虚拟工厂,这些虚拟工厂执行所有子任务的代价
Figure BDA0003190989300000066
置为0,同样,对于n>m的情况,添加n-m个虚拟子任务,每个工厂执行虚拟子任务的代价
Figure BDA0003190989300000067
置为0,那么,问题就转化为将u个任务分配给u个工厂的分布式任务分配问题,
Figure BDA0003190989300000068
数学模型转化为:
Figure BDA0003190989300000069
Figure BDA00031909893000000610
Figure BDA00031909893000000611
Figure BDA00031909893000000612
所述步骤B2具体实现如下:
由于x∈{0,1}n等价于
Figure BDA00031909893000000613
Figure BDA00031909893000000614
替换为盒约束与lp球约束的交集,引入辅助变量yi以及zi,将问题模型转换如下:
Figure BDA0003190989300000071
Figure BDA0003190989300000072
Figure BDA0003190989300000073
Figure BDA0003190989300000074
Figure BDA0003190989300000075
yi∈Sb
zi∈Sp
其中,Sb={yi|||yi||≤1},
Figure BDA0003190989300000076
在本发明一实施例中,所述步骤B3具体实现如下:
引入
Figure BDA0003190989300000077
以及
Figure BDA0003190989300000078
考虑步骤B2中约束,构造增广拉格朗日函数如下:
Figure BDA0003190989300000079
所述步骤B4具体实现如下:
将问题根据工厂分解成u个子问题,定义
Figure BDA00031909893000000710
即边缘端通信并计算得到的其他边缘端对于子任务
Figure BDA00031909893000000711
的决策变量的和,每个边缘端需要最小化的拉格朗日函数表示如下:
Figure BDA00031909893000000712
Figure BDA0003190989300000081
在本发明一实施例中,所述步骤B5具体包括以下步骤:
步骤B51、每个边缘端初始化
Figure BDA0003190989300000082
yi,zi,λi,ω,θ,σ以及ρi
步骤B52、在未达到收敛条件的情况下,每个工厂根据更新公式更新yi,zi
Figure BDA0003190989300000083
ω,θ,σ,λi以及ρi
a)所述收敛条件为
Figure BDA0003190989300000084
Figure BDA0003190989300000085
其中,stop为大于零的收敛指标,可根据实际情况取值;
b)在每个边缘端,yi,zi
Figure BDA0003190989300000086
ω,θ,σ,λi以及ρi依次更新:
1)首先,对于每一个yij,其更新公式如下:
Figure BDA0003190989300000087
2)关于zi,更新公式如下:
Figure BDA0003190989300000088
其中
Figure BDA0003190989300000089
3)在工厂
Figure BDA00031909893000000810
与决策变量
Figure BDA00031909893000000811
相关的拉格朗日函数表示如下:
Figure BDA00031909893000000812
引入
Figure BDA00031909893000000813
即为除了
Figure BDA00031909893000000814
以外的其他u-1个决策变量的和,将问题拆分成u个子问题,每个子问题表示如下:
Figure BDA00031909893000000815
可以看出上述式子是关于
Figure BDA0003190989300000091
的二次函数,要求取最小值,应该让导数等于0,由此可得如下式子:
Figure BDA0003190989300000092
因此,
Figure BDA0003190989300000093
的更新公式为:
Figure BDA0003190989300000094
4)ωj的更新公式如下:
Figure BDA0003190989300000095
5)θij的更新公式如下:
Figure BDA0003190989300000096
6)σij的更新公式如下:
Figure BDA0003190989300000097
7)λi的更新公式如下:
Figure BDA0003190989300000098
8)ρi的更新公式如下:
Figure BDA0003190989300000099
μ会根据实际问题以及收敛速度取值,有时,为防止ρi增长过快,还会为其设置上限;
步骤B53收敛后,每个边缘端可根据自己的决策变量知道自己要执行的子任务,子任务也会立刻运输到边缘端所管理的工厂以投入生产。
本发明还提供了一种工业4.0需求驱动的分布式动态边缘云智能制造系统,采用如权利要求1-9任一所述的工业4.0需求驱动的分布式动态边缘云智能制造方法,包括需求获取与设计模块、子任务分配模块、打包运输模块;具体为:
A、需求发布与设计模块:用于客户发布个性化需求以及设计师设计产品,包括需求发布模块、设计模块以及产品分解模块;首先,客户向云端发布个性化产品需求,其次,设计师通过云端获取需求信息,发布初步设计,然后,设计师与客户沟通,确定最终设计,最后,云端根据最终设计,将个性化产品分解为几个子任务以供后续生产制造;
B、子任务分配模块:用于在决策时刻将来自多个个性化产品的待分配子任务通过分布式方法分配给可用工厂,从而最大化所有可用工厂总利润,包括待分配子任务信息发布模块,工厂能力判断模块以及基于lp-BoxADMM的分布式子任务分配模块;在决策时刻,有多个子任务待分配,边缘端根据云端发布的待分配子任务信息,判断当前工厂状态是否能够执行任一待分配子任务,并将结果返回云端,云端收集到可用工厂信息后,采用动态分布式任务分配模型以及基于lp-BoxADMM的方法通过协调可用工厂完成任务分配;
C、打包运输模块:用于打包和运输未完成/完成的产品以及原材料,包括未完成子任务运输模块、完成子任务运输模块以及原材料运输模块;首先,所述模块可用于工厂之间运输未完成的产品以将当前待分配子任务运输到对应的工厂,其次,所述模块可用于打包和运输完成的产品到达客户指定的地址,最后,所述模块还用于运输执行子任务所需要购买的原材料;该模块是否有能力运输未完成的产品以及工厂亟需的原材料,也决定了工厂在决策时刻是否有能力执行待分配子任务。
相较于现有技术,本发明具有以下有益效果:本发明以客户需求为驱动,让客户参与产品生产全生命周期并实时提出修改意见,提高了产品个性化程度,同时,采用分布式边缘云架构,减轻云端负担,保障数据隐私,避免单点故障,提高鲁棒性。
附图说明
图1是本发明提出的工业4.0需求驱动的分布式动态边缘云智能制造模式示意图。
图2是本发明基于lp-BoxADMM的方法完成任务分配的流程图。
图3是本发明提出的工业4.0需求驱动的分布式动态边缘云智能制造系统结构图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
本发明一种工业4.0需求驱动的分布式动态边缘云智能制造方法,如图1所示,包括以下步骤:
步骤S1、客户向云端发布个性化产品需求;
步骤S2、设计师获取需求信息,与客户沟通,确定最终设计;
步骤S3、云端根据最终设计,将个性化产品分解为几个子任务;
步骤S4、在决策时刻,边缘端根据云端发布的待分配子任务信息,判断当前工厂状态是否能够执行任一待分配子任务,并将结果返回云端;
步骤S5、建立动态分布式任务分配模型;
步骤S6、采用基于lp-BoxADMM的方法完成任务分配;
步骤S7、待个性化产品的所有子任务完成,将产品打包并运输给客户。
所述步骤S5包括以下步骤:
步骤A1、以最大化决策时刻t所有工厂总利润为目标,建立目标函数;
步骤A2、根据所述场景实际情况,添加必要约束;
步骤A3、将模型转化成标准的求解最小化问题的分布式优化数学模型。
所述步骤A1包括以下步骤:
动态分布式任务分配模型的目标是最大化决策时刻t所有工厂总利润,目标函数如下:
Figure BDA0003190989300000111
其中,i,j分别表示工厂和子任务的下标,n,m分别表示可用工厂和待分配子任务的数量,
Figure BDA0003190989300000112
表示工厂
Figure BDA0003190989300000113
执行子任务
Figure BDA0003190989300000114
能够获得的效用,
Figure BDA0003190989300000115
表示执行代价,
Figure BDA0003190989300000116
表示运输代价。
Figure BDA0003190989300000117
表示工厂
Figure BDA0003190989300000118
在决策时刻t是否有能力执行子任务
Figure BDA0003190989300000119
如果是的话,值为1,否则,值为0。
Figure BDA00031909893000001110
为待求解的决策变量,表示在决策时刻t,是否将子任务
Figure BDA00031909893000001111
分配给工厂
Figure BDA00031909893000001112
执行,是的话,值为1,否则,值为0。
每个工厂的利润都可以表示为效用减去执行代价减去运输代价。因此,效用
Figure BDA00031909893000001113
执行代价
Figure BDA00031909893000001114
以及运输代价
Figure BDA00031909893000001115
的获取至关重要。
效用
Figure BDA00031909893000001116
受工厂生产能力匹配度以及工厂执行子任务所花费的时间长短的影响。执行每个子任务
Figure BDA00031909893000001117
都需要原材料
Figure BDA00031909893000001118
以及一定的生产能力
Figure BDA00031909893000001119
其中q为需要的生产能力种类总数。生产需要的每种原材料
Figure BDA00031909893000001120
的数量表示为
Figure BDA00031909893000001121
Figure BDA00031909893000001122
每个工厂
Figure BDA00031909893000001123
在决策时刻所拥有的生产能力表示为
Figure BDA00031909893000001124
因此,定义
Figure BDA00031909893000001125
表示工厂
Figure BDA00031909893000001126
的第h种生产能力是否与子任务
Figure BDA00031909893000001127
执行所需要的第k种能力相同,如果是的话,值为1,否则为0。不仅如此,效用的计算还需要获取先前完成子任务
Figure BDA00031909893000001128
最快用时TTj以及子任务
Figure BDA00031909893000001129
最高售价pricej,并且,预估工厂
Figure BDA00031909893000001130
执行子任务
Figure BDA00031909893000001131
用时为
Figure BDA00031909893000001132
由此可得,效用
Figure BDA00031909893000001133
计算如下:
Figure BDA0003190989300000121
执行代价
Figure BDA0003190989300000122
包括机器设置成本
Figure BDA0003190989300000123
生产成本
Figure BDA0003190989300000124
(单位时间生产成本乘以生产时间)以及原材料购买成本。每种原材料
Figure BDA0003190989300000125
的成本的计算包括两种情况,第一种情况是原材料
Figure BDA0003190989300000126
需要的数量
Figure BDA0003190989300000127
大于工厂库存的数量
Figure BDA0003190989300000128
那么原材料成本就是购买成本
Figure BDA0003190989300000129
(需要购买的数量乘以单位购买成本)加上储存成本
Figure BDA00031909893000001210
(原材料储存数量乘以原材料单位储存成本),如果原材料数量足够生产,那么只需要储存成本
Figure BDA00031909893000001211
基于此,引入
Figure BDA00031909893000001212
来判断需要的原材料数量是否多于原材料库存数量,如果是的话,值为1,否则,值为0。因此,执行代价计算公式如下:
Figure BDA00031909893000001213
运输代价
Figure BDA00031909893000001214
包括子任务
Figure BDA00031909893000001215
运输到工厂
Figure BDA00031909893000001216
的代价
Figure BDA00031909893000001217
以及各种需要购买的原材料的运输代价。每种需要购买的原材料
Figure BDA00031909893000001218
的运输代价为购买数量
Figure BDA00031909893000001219
乘以单位运输代价
Figure BDA00031909893000001220
因此,运输代价
Figure BDA00031909893000001221
计算公式表示如下:
Figure BDA00031909893000001222
所述步骤A2包括以下步骤:
动态分布式任务分配模型的约束主要从以下几个方面考虑:
1)
Figure BDA00031909893000001223
以及
Figure BDA00031909893000001224
的取值为0或1,表示如下:
Figure BDA00031909893000001225
Figure BDA00031909893000001226
2)此外,子任务运输代价
Figure BDA00031909893000001227
和原材料运输代价
Figure BDA00031909893000001228
的取值如下:
Figure BDA00031909893000001229
Figure BDA0003190989300000131
运输代价取值为-1,表示没有能力运输,即子任务
Figure BDA0003190989300000132
或者购买的原材料
Figure BDA0003190989300000133
无法运输到工厂
Figure BDA0003190989300000134
运输代价为0,对于子任务
Figure BDA0003190989300000135
来说,可能存在三种情况,当前子任务
Figure BDA0003190989300000136
为第一个子任务、当前子任务就在工厂
Figure BDA0003190989300000137
无需运输,或者因为其他情况导致子任务运输代价为0;对于原材料来说,就是指原材料
Figure BDA0003190989300000138
运输到工厂
Figure BDA0003190989300000139
代价为0。其他情况下,运输代价为大于0的实数。
3)假定待分配子任务在决策时刻至少有一个工厂可以执行,因此,就有了如下关于
Figure BDA00031909893000001310
的约束:
Figure BDA00031909893000001311
4)
Figure BDA00031909893000001312
用来表示工厂
Figure BDA00031909893000001313
在决策时刻t是否有能力执行子任务
Figure BDA00031909893000001314
如果有的话,值为1,否则,值为0。有如下约束会导致
Figure BDA00031909893000001315
Figure BDA00031909893000001316
Figure BDA00031909893000001317
Figure BDA00031909893000001318
Figure BDA00031909893000001319
这四个约束从上到下依次表示:
·如果子任务
Figure BDA00031909893000001320
无法运输到工厂
Figure BDA00031909893000001321
(即
Figure BDA00031909893000001322
),那么是工厂
Figure BDA00031909893000001323
无法执行子任务
Figure BDA00031909893000001324
·如果工厂
Figure BDA00031909893000001325
执行子任务
Figure BDA00031909893000001326
必须购买原材料
Figure BDA00031909893000001327
(即
Figure BDA00031909893000001328
),但是购买原材料
Figure BDA00031909893000001329
无法运输到工厂
Figure BDA00031909893000001330
那么,子任务
Figure BDA00031909893000001331
无法在工厂
Figure BDA00031909893000001332
执行
·如果工厂
Figure BDA00031909893000001333
能够匹配给子任务
Figure BDA00031909893000001334
的能力小于总能力种类的一半,那么,也视为工厂
Figure BDA00031909893000001335
没有能力执行子任务
Figure BDA00031909893000001336
·如果工厂
Figure BDA00031909893000001337
所获得的利润小于0的话,那么也视为工厂
Figure BDA00031909893000001338
无法执行子任务
Figure BDA00031909893000001339
5)关于决策变量
Figure BDA00031909893000001340
有如下约束:
Figure BDA0003190989300000141
Figure BDA0003190989300000142
前者表示任何一个子任务都最多被一个工厂执行,后者表示任何一个工厂都最多只能执行一个子任务。
6)最后,是关于
Figure BDA0003190989300000143
的约束:
Figure BDA0003190989300000144
即如果工厂
Figure BDA0003190989300000145
执行子任务
Figure BDA0003190989300000146
需要购买原材料
Figure BDA0003190989300000147
那么,
Figure BDA0003190989300000148
取值为1,否则为0。
所述步骤A3包括以下步骤:
Figure BDA0003190989300000149
定义
Figure BDA00031909893000001410
效用
Figure BDA00031909893000001411
执行代价
Figure BDA00031909893000001412
运输代价
Figure BDA00031909893000001413
以及工厂执行能力变量
Figure BDA00031909893000001414
都可在决策时刻t根据计算实时获得,因此,求解时只需要考虑与
Figure BDA00031909893000001415
相关的约束,则,数学模型可表示为:
Figure BDA00031909893000001416
Figure BDA00031909893000001417
Figure BDA00031909893000001418
Figure BDA00031909893000001419
由于每个管理工厂
Figure BDA00031909893000001420
的边缘端只知道自身的目标函数fi(xi),并且也只能决定与自身相关的决策变量
Figure BDA00031909893000001421
的取值情况,因此这是一个分布式优化问题,该问题假设每个边缘端之间能够正常通信,并且能够获得来自其他边缘端刚更新过的决策变量
Figure BDA00031909893000001422
所述步骤S6如图2所示,主要包括以下步骤:
步骤B1添加虚拟工厂或虚拟子任务,将不等式约束转化为等式约束;
步骤B2将决策变量二值约束转化成盒约束与lp球约束的交集;
步骤B3构造增广拉格朗日函数;
步骤B4将问题分解为每个工厂要解决的子问题;
步骤B5根据更新公式以及收敛条件进行分布式求解。
所述步骤B1包括以下步骤:
为了便于求解,令u=max{n,m},对于n<m的情况,添加m-n个虚拟工厂,这些虚拟工厂执行所有子任务的代价
Figure BDA0003190989300000151
置为0,同样,对于n>m的情况,添加n-m个虚拟子任务,每个工厂执行虚拟子任务的代价
Figure BDA0003190989300000152
置为0,那么,问题就转化为将u个任务(包括虚拟任务)分配给u个工厂(包括虚拟工厂)的分布式任务分配问题,
Figure BDA0003190989300000153
数学模型转化为:
Figure BDA0003190989300000154
Figure BDA0003190989300000155
Figure BDA0003190989300000156
Figure BDA0003190989300000157
所述步骤B2包括以下步骤:
由于x∈{0,1}n等价于
Figure BDA0003190989300000158
Figure BDA0003190989300000159
替换为盒约束与lp球约束的交集,引入辅助变量yi以及zi,将问题模型转换如下:
Figure BDA00031909893000001510
Figure BDA00031909893000001511
Figure BDA00031909893000001512
Figure BDA00031909893000001513
Figure BDA00031909893000001514
yi∈Sb
zi∈Sp
其中,Sb={yi|||yi||≤1},
Figure BDA0003190989300000161
所述步骤B3包括以下步骤:
引入
Figure BDA0003190989300000162
以及
Figure BDA0003190989300000163
考虑步骤B2中所述约束,构造增广拉格朗日函数如下:
Figure BDA0003190989300000164
所述步骤B4包括以下步骤:
将问题根据工厂(包括虚拟工厂)分解成u个子问题,定义
Figure BDA0003190989300000165
Figure BDA0003190989300000166
即边缘端通信并计算得到的其他边缘端对于子任务
Figure BDA0003190989300000167
的决策变量的和,每个边缘端需要最小化的拉格朗日函数表示如下:
Figure BDA0003190989300000168
所述步骤B5包括以下步骤:
步骤B51每个边缘端初始化
Figure BDA0003190989300000169
yi,zi,λi,ω,θ,σ以及ρi
步骤B52在未达到收敛条件的情况下,每个工厂根据更新公式更新yi,zi
Figure BDA00031909893000001610
ω,θ,σ,λi以及ρi
a)所述收敛条件为
Figure BDA0003190989300000171
Figure BDA0003190989300000172
其中,stop为大于零的收敛指标,可根据实际情况取值。
b)在每个边缘端,yi,zi
Figure BDA0003190989300000173
ω,θ,σ,λi以及ρi依次更新:
1)首先,对于每一个yij,其更新公式如下:
Figure BDA0003190989300000174
2)关于zi,更新公式如下:
Figure BDA0003190989300000175
其中
Figure BDA0003190989300000176
3)在工厂
Figure BDA0003190989300000177
与决策变量
Figure BDA0003190989300000178
相关的拉格朗日函数表示如下:
Figure BDA0003190989300000179
引入
Figure BDA00031909893000001710
即为除了
Figure BDA00031909893000001711
以外的其他u-1个决策变量的和,将问题拆分成u个子问题,每个子问题表示如下:
Figure BDA00031909893000001712
可以看出上述式子是关于
Figure BDA00031909893000001713
的二次函数,要求取最小值,应该让导数等于0,由此可得如下式子:
Figure BDA00031909893000001714
因此,
Figure BDA00031909893000001715
的更新公式为:
Figure BDA0003190989300000181
4)ωj的更新公式如下:
Figure BDA0003190989300000182
5)θij的更新公式如下:
Figure BDA0003190989300000183
6)σij的更新公式如下:
Figure BDA0003190989300000184
7)λi的更新公式如下:
Figure BDA0003190989300000185
8)ρi的更新公式如下:
Figure BDA0003190989300000186
μ一般会根据实际问题以及收敛速度取值,有时,为了防止ρi增长过快,还会为其设置上限。
步骤B53收敛后,每个边缘端可根据自己的决策变量知道自己要执行的子任务,子任务也会立刻运输到边缘端所管理的工厂以投入生产。步骤B53收敛后,每个边缘端根据决策变量得到该工厂分配到的子任务。
本发明还提供了一种工业4.0需求驱动的分布式动态边缘云智能制造系统,如图3所示,采用工业4.0需求驱动的分布式动态边缘云智能制造方法,包括需求获取与设计模块、子任务分配模块、打包运输模块;具体为
A.需求发布与设计模块:用于客户发布个性化需求以及设计师设计产品,包括需求发布模块、设计模块以及产品分解模块;首先,客户向云端发布个性化产品需求,其次,设计师通过云端获取需求信息,发布初步设计,然后,设计师与客户沟通,确定最终设计,最后,云端根据最终设计,将个性化产品分解为几个子任务以供后续生产制造。
B.子任务分配模块:用于在决策时刻将来自多个个性化产品的待分配子任务通过分布式方法分配给可用工厂,从而最大化所有可用工厂总利润,包括待分配子任务信息发布模块,工厂能力判断模块以及基于lp-Box ADMM的分布式子任务分配模块。在决策时刻,有多个子任务待分配,边缘端根据云端发布的待分配子任务信息,判断当前工厂状态是否能够执行任一待分配子任务,并将结果返回云端,云端收集到可用工厂信息后,采用动态分布式任务分配模型以及基于lp-Box ADMM的方法通过协调可用工厂完成任务分配。
C.打包运输模块:用于打包和运输未完成/完成的产品以及原材料,包括未完成子任务运输模块、完成子任务运输模块以及原材料运输模块;首先,所述模块可用于工厂之间运输未完成的产品以将当前待分配子任务运输到对应的工厂,其次,所述模块可用于打包和运输完成的产品到达客户指定的地址,最后,所述模块还用于运输执行子任务所需要购买的原材料。该模块是否有能力运输未完成的产品以及工厂亟需的原材料,也决定了工厂在决策时刻是否有能力执行待分配子任务。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。

Claims (8)

1.一种工业4.0需求驱动的分布式动态边缘云智能制造方法,其特征在于,包括以下步骤:
步骤S1、客户向云端发布个性化产品需求;
步骤S2、设计师获取需求信息,确定最终设计;
步骤S3、云端根据最终设计,将个性化产品分解为几个子任务;
步骤S4、在决策时刻,边缘端根据云端发布的待分配子任务信息,判断当前工厂状态是否能够执行任一待分配子任务,并将结果返回云端;
步骤S5、建立动态分布式任务分配模型;
步骤S6、采用基于lp-BoxADMM的方法完成任务分配;
步骤S7、待个性化产品的所有子任务完成,将产品打包并运输给客户;
所述步骤S5具体包括以下步骤:
步骤A1、以最大化决策时刻t所有工厂总利润为目标,建立目标函数;
步骤A2、根据场景实际情况,添加约束条件;
步骤A3、将模型转化成标准的求解最小化问题的分布式优化数学模型;
所述步骤A1具体实现如下:
动态分布式任务分配模型的目标是最大化决策时刻t所有工厂总利润,目标函数如下:
Figure FDA0003729795900000011
其中,i,j分别表示工厂和子任务的下标,n,m分别表示可用工厂和待分配子任务的数量,
Figure FDA0003729795900000012
表示工厂
Figure FDA0003729795900000013
执行子任务
Figure FDA0003729795900000014
能够获得的效用,
Figure FDA0003729795900000015
表示执行代价,
Figure FDA0003729795900000016
表示运输代价;
Figure FDA0003729795900000017
表示工厂
Figure FDA0003729795900000018
在决策时刻t是否有能力执行子任务
Figure FDA0003729795900000019
如果是的话,值为1,否则,值为0;
Figure FDA00037297959000000110
为待求解的决策变量,表示在决策时刻t,是否将子任务
Figure FDA00037297959000000111
分配给工厂
Figure FDA00037297959000000112
执行,是的话,值为1,否则,值为0;每个工厂的利润都可以表示为效用减去执行代价减去运输代价;
效用
Figure FDA00037297959000000113
受工厂生产能力匹配度以及工厂执行子任务所花费的时间长短的影响,执行每个子任务
Figure FDA00037297959000000114
都需要原材料
Figure FDA00037297959000000115
以及生产能力
Figure FDA00037297959000000116
其中q为需要的生产能力种类总数,生产需要的每种原材料
Figure FDA00037297959000000117
的数量表示为
Figure FDA00037297959000000118
每个工厂
Figure FDA00037297959000000119
在决策时刻所拥有的生产能力表示为
Figure FDA00037297959000000120
因此,定义
Figure FDA00037297959000000121
表示工厂
Figure FDA00037297959000000122
的第h种生产能力是否与子任务
Figure FDA00037297959000000123
执行所需要的第k种能力相同,如果是的话,值为1,否则为0;效用的计算还需要获取先前完成子任务
Figure FDA0003729795900000021
最快用时TTj以及子任务
Figure FDA0003729795900000022
最高售价pricej,并且,预估工厂
Figure FDA0003729795900000023
执行子任务
Figure FDA0003729795900000024
用时为
Figure FDA0003729795900000025
由此可得,效用
Figure FDA0003729795900000026
计算如下:
Figure FDA0003729795900000027
执行代价
Figure FDA0003729795900000028
包括机器设置成本
Figure FDA0003729795900000029
生产成本
Figure FDA00037297959000000210
以及原材料购买成本,每种原材料
Figure FDA00037297959000000211
的成本的计算包括两种情况,第一种情况是原材料
Figure FDA00037297959000000212
需要的数量
Figure FDA00037297959000000213
大于工厂库存的数量
Figure FDA00037297959000000214
那么原材料成本就是购买成本
Figure FDA00037297959000000215
加上储存成本
Figure FDA00037297959000000216
如果原材料数量足够生产,那么只需要储存成本
Figure FDA00037297959000000217
基于此,引入
Figure FDA00037297959000000218
来判断需要的原材料数量是否多于原材料库存数量,如果是的话,值为1,否则,值为0;因此,执行代价计算公式如下:
Figure FDA00037297959000000219
运输代价
Figure FDA00037297959000000220
包括子任务
Figure FDA00037297959000000221
运输到工厂
Figure FDA00037297959000000222
的代价
Figure FDA00037297959000000223
以及各种需要购买的原材料的运输代价,每种需要购买的原材料
Figure FDA00037297959000000224
的运输代价为购买数量
Figure FDA00037297959000000225
乘以单位运输代价
Figure FDA00037297959000000226
因此,运输代价
Figure FDA00037297959000000227
计算公式表示如下:
Figure FDA00037297959000000228
2.根据权利要求1所述的工业4.0需求驱动的分布式动态边缘云智能制造方法,其特征在于,所述步骤A2具体实现如下:
动态分布式任务分配模型的约束包括:
1)
Figure FDA00037297959000000229
以及
Figure FDA00037297959000000230
的取值为0或1,表示如下:
Figure FDA00037297959000000231
Figure FDA00037297959000000232
2)子任务运输代价
Figure FDA00037297959000000233
和单位运输代价
Figure FDA00037297959000000234
的取值如下:
Figure FDA00037297959000000235
Figure FDA0003729795900000031
运输代价取值为-1,表示没有能力运输,即子任务
Figure FDA0003729795900000032
或者购买的原材料
Figure FDA0003729795900000033
无法运输到工厂
Figure FDA0003729795900000034
运输代价为0,对于子任务
Figure FDA0003729795900000035
来说,存在三种情况,当前子任务
Figure FDA0003729795900000036
为第一个子任务、当前子任务就在工厂
Figure FDA0003729795900000037
无需运输,或者因为其他情况导致子任务运输代价为0;对于原材料来说,就是指原材料
Figure FDA0003729795900000038
运输到工厂
Figure FDA0003729795900000039
代价为0,其他情况下,运输代价为大于0的实数;
3)假定待分配子任务在决策时刻至少有一个工厂可以执行,因此,关于
Figure FDA00037297959000000310
的约束:
Figure FDA00037297959000000311
4)
Figure FDA00037297959000000312
用来表示工厂
Figure FDA00037297959000000313
在决策时刻t是否有能力执行子任务
Figure FDA00037297959000000314
如果有的话,值为1,否则,值为0,有如下约束会导致
Figure FDA00037297959000000315
Figure FDA00037297959000000316
Figure FDA00037297959000000317
Figure FDA00037297959000000318
Figure FDA00037297959000000319
这四个约束从上到下依次表示:
1)如果子任务
Figure FDA00037297959000000320
无法运输到工厂
Figure FDA00037297959000000321
Figure FDA00037297959000000322
那么是工厂
Figure FDA00037297959000000323
无法执行子任务
Figure FDA00037297959000000324
2)如果工厂
Figure FDA00037297959000000325
执行子任务
Figure FDA00037297959000000326
必须购买原材料
Figure FDA00037297959000000327
Figure FDA00037297959000000328
但是购买原材料
Figure FDA00037297959000000329
无法运输到工厂
Figure FDA00037297959000000330
那么,子任务
Figure FDA00037297959000000331
无法在工厂
Figure FDA00037297959000000332
执行;
3)如果工厂
Figure FDA00037297959000000333
能够匹配给子任务
Figure FDA00037297959000000334
的能力小于总能力种类的一半,那么,也视为工厂
Figure FDA00037297959000000335
没有能力执行子任务
Figure FDA00037297959000000336
4)如果工厂
Figure FDA00037297959000000337
所获得的利润小于0的话,那么也视为工厂
Figure FDA00037297959000000338
无法执行子任务
Figure FDA00037297959000000339
5)关于决策变量
Figure FDA00037297959000000340
有如下约束:
Figure FDA0003729795900000041
Figure FDA0003729795900000042
前者表示任何一个子任务都最多被一个工厂执行,后者表示任何一个工厂都最多只能执行一个子任务;
6)关于
Figure FDA0003729795900000043
的约束:
Figure FDA0003729795900000044
即如果工厂
Figure FDA0003729795900000045
执行子任务
Figure FDA0003729795900000046
需要购买原材料
Figure FDA0003729795900000047
那么,
Figure FDA0003729795900000048
取值为1,否则为0。
3.根据权利要求2所述的工业4.0需求驱动的分布式动态边缘云智能制造方法,其特征在于,所述步骤A3具体实现如下:
Figure FDA0003729795900000049
定义
Figure FDA00037297959000000410
效用
Figure FDA00037297959000000411
执行代价
Figure FDA00037297959000000412
运输代价
Figure FDA00037297959000000413
以及工厂执行能力变量
Figure FDA00037297959000000414
都可在决策时刻t根据计算实时获得,因此,求解时只需要考虑与
Figure FDA00037297959000000415
相关的约束,则数学模型可表示为:
Figure FDA00037297959000000416
Figure FDA00037297959000000417
Figure FDA00037297959000000418
Figure FDA00037297959000000419
由于每个管理工厂
Figure FDA00037297959000000420
的边缘端只知道自身的目标函数fi(xi),并且也只能决定与自身相关的决策变量
Figure FDA00037297959000000421
的取值情况,因此这是一个分布式优化问题,该问题假设每个边缘端之间能够正常通信,并且能够获得来自其他边缘端刚更新过的决策变量
Figure FDA00037297959000000422
4.根据权利要求3所述的工业4.0需求驱动的分布式动态边缘云智能制造方法,其特征在于,所述步骤S6具体包括以下步骤:
步骤B1、添加虚拟工厂或虚拟子任务,将不等式约束转化为等式约束;
步骤B2、将决策变量二值约束转化成盒约束与lp球约束的交集;
步骤B3、构造增广拉格朗日函数;
步骤B4、将问题分解为每个工厂要解决的子问题;
步骤B5、根据更新公式以及收敛条件进行分布式求解。
5.根据权利要求4所述的工业4.0需求驱动的分布式动态边缘云智能制造方法,其特征在于,所述步骤B1具体实现如下:
令u=max{n,m},对于n<m的情况,添加m-n个虚拟工厂,这些虚拟工厂执行所有子任务的代价
Figure FDA0003729795900000051
置为0,同样,对于n>m的情况,添加n-m个虚拟子任务,每个工厂执行虚拟子任务的代价
Figure FDA0003729795900000052
置为0,那么,问题就转化为将u个任务分配给u个工厂的分布式任务分配问题,
Figure FDA0003729795900000053
数学模型转化为:
Figure FDA0003729795900000054
Figure FDA0003729795900000055
Figure FDA0003729795900000056
Figure FDA0003729795900000057
所述步骤B2具体实现如下:
由于x∈{0,1}n等价于
Figure FDA0003729795900000058
Figure FDA0003729795900000059
替换为盒约束与lp球约束的交集,引入辅助变量yi以及zi,将问题模型转换如下:
Figure FDA00037297959000000510
Figure FDA00037297959000000511
Figure FDA00037297959000000512
Figure FDA00037297959000000513
Figure FDA0003729795900000061
yi∈Sb
zi∈Sp
其中,Sb={yi|||yi||≤1},
Figure FDA0003729795900000062
6.根据权利要求5所述的工业4.0需求驱动的分布式动态边缘云智能制造方法,其特征在于,所述步骤B3具体实现如下:
引入
Figure FDA0003729795900000063
以及
Figure FDA0003729795900000064
考虑步骤B2中约束,构造增广拉格朗日函数如下:
Figure FDA0003729795900000065
所述步骤B4具体实现如下:
将问题根据工厂分解成u个子问题,定义
Figure FDA0003729795900000066
即边缘端通信并计算得到的其他边缘端对于子任务
Figure FDA0003729795900000067
的决策变量的和,每个边缘端需要最小化的拉格朗日函数表示如下:
Figure FDA0003729795900000068
7.根据权利要求6所述的工业4.0需求驱动的分布式动态边缘云智能制造方法,其特征在于,所述步骤B5具体包括以下步骤:
步骤B51、每个边缘端初始化
Figure FDA0003729795900000069
yi,zi,λi,ω,
Figure FDA00037297959000000610
σ以及ρi
步骤B52、在未达到收敛条件的情况下,每个工厂根据更新公式更新yi,zi
Figure FDA0003729795900000071
ω,
Figure FDA0003729795900000072
σ,λi以及ρi
a)所述收敛条件为
Figure FDA0003729795900000073
Figure FDA0003729795900000074
其中,stop为大于零的收敛指标,可根据实际情况取值;
b)在每个边缘端,yi,zi
Figure FDA0003729795900000075
ω,
Figure FDA0003729795900000076
σ,λi以及ρi依次更新:
1)首先,对于每一个yij,其更新公式如下:
Figure FDA0003729795900000077
2)关于zi,更新公式如下:
Figure FDA0003729795900000078
其中
Figure FDA0003729795900000079
3)在工厂
Figure FDA00037297959000000710
与决策变量
Figure FDA00037297959000000711
相关的拉格朗日函数表示如下:
Figure FDA00037297959000000712
引入
Figure FDA00037297959000000713
即为除了
Figure FDA00037297959000000714
以外的其他u-1个决策变量的和,将问题拆分成u个子问题,每个子问题表示如下:
Figure FDA00037297959000000715
可以看出上述式子是关于
Figure FDA00037297959000000716
的二次函数,要求取最小值,应该让导数等于0,由此可得如下式子:
Figure FDA00037297959000000717
因此,
Figure FDA00037297959000000718
的更新公式为:
Figure FDA0003729795900000081
4)ωj的更新公式如下:
Figure FDA0003729795900000082
5)
Figure FDA0003729795900000083
的更新公式如下:
Figure FDA0003729795900000084
6)σij的更新公式如下:
Figure FDA0003729795900000085
7)λi的更新公式如下:
Figure FDA0003729795900000086
8)ρi的更新公式如下:
Figure FDA0003729795900000087
μ会根据实际问题以及收敛速度取值,有时,为防止ρi增长过快,还会为其设置上限;
步骤B53收敛后,每个边缘端可根据自己的决策变量知道自己要执行的子任务,子任务也会立刻运输到边缘端所管理的工厂以投入生产。
8.一种工业4.0需求驱动的分布式动态边缘云智能制造系统,其特征在于,采用如权利要求1-7任一所述的工业4.0需求驱动的分布式动态边缘云智能制造方法,包括需求获取与设计模块、子任务分配模块、打包运输模块;具体为:
A、需求发布与设计模块:用于客户发布个性化需求以及设计师设计产品,包括需求发布模块、设计模块以及产品分解模块;首先,客户向云端发布个性化产品需求,其次,设计师通过云端获取需求信息,发布初步设计,然后,设计师与客户沟通,确定最终设计,最后,云端根据最终设计,将个性化产品分解为几个子任务以供后续生产制造;
B、子任务分配模块:用于在决策时刻将来自多个个性化产品的待分配子任务通过分布式方法分配给可用工厂,从而最大化所有可用工厂总利润,包括待分配子任务信息发布模块,工厂能力判断模块以及基于lp-BoxADMM的分布式子任务分配模块;在决策时刻,有多个子任务待分配,边缘端根据云端发布的待分配子任务信息,判断当前工厂状态是否能够执行任一待分配子任务,并将结果返回云端,云端收集到可用工厂信息后,采用动态分布式任务分配模型以及基于lp-BoxADMM的方法通过协调可用工厂完成任务分配;
C、打包运输模块:用于打包和运输未完成/完成的产品以及原材料,包括未完成子任务运输模块、完成子任务运输模块以及原材料运输模块;首先,所述模块可用于工厂之间运输未完成的产品以将当前待分配子任务运输到对应的工厂,其次,所述模块可用于打包和运输完成的产品到达客户指定的地址,最后,所述模块还用于运输执行子任务所需要购买的原材料;该模块是否有能力运输未完成的产品以及工厂亟需的原材料,也决定了工厂在决策时刻是否有能力执行待分配子任务。
CN202110879349.7A 2021-08-02 2021-08-02 工业4.0需求驱动的分布式动态边缘云智能制造方法及系统 Active CN113467405B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110879349.7A CN113467405B (zh) 2021-08-02 2021-08-02 工业4.0需求驱动的分布式动态边缘云智能制造方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110879349.7A CN113467405B (zh) 2021-08-02 2021-08-02 工业4.0需求驱动的分布式动态边缘云智能制造方法及系统

Publications (2)

Publication Number Publication Date
CN113467405A CN113467405A (zh) 2021-10-01
CN113467405B true CN113467405B (zh) 2022-09-02

Family

ID=77883738

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110879349.7A Active CN113467405B (zh) 2021-08-02 2021-08-02 工业4.0需求驱动的分布式动态边缘云智能制造方法及系统

Country Status (1)

Country Link
CN (1) CN113467405B (zh)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110429999A (zh) * 2019-04-24 2019-11-08 西安电子科技大学 基于lp-Box ADMM算法的大规模MIMO检测方法
CN112465394B (zh) * 2020-12-09 2022-06-14 福州大学 一种用于工业4.0大规模个性化生产的动态云制造方法
CN112884241B (zh) * 2021-03-12 2023-09-01 重庆大学 基于智能Agent的云边协作制造任务调度方法

Also Published As

Publication number Publication date
CN113467405A (zh) 2021-10-01

Similar Documents

Publication Publication Date Title
CN110046865B (zh) 分布式库存调度方法
Ding et al. Hidden Markov model-based autonomous manufacturing task orchestration in smart shop floors
Zhou et al. Reinforcement learning with composite rewards for production scheduling in a smart factory
CN110084405B (zh) 吞吐柔性智能装配物流路径规划方法
CN105974891B (zh) 一种基于动态看板的模具生产过程自适应控制方法
CN111882215A (zh) 一种含有agv的个性化定制柔性作业车间调度方法
CN112053082A (zh) 基于云计算的物料管理系统
CN111950939A (zh) 基于智能制造的物料管理方法
CN111798097B (zh) 一种基于市场机制的自主移动机器人任务分配处理方法
CN105005503B (zh) 基于元胞自动机的云计算负载均衡任务调度方法
CN114503038A (zh) 使用状态矩阵自学习柔性制造系统的制造调度的方法及设备
CN114936783B (zh) 一种基于mmddpg算法的rgv小车调度方法及系统
CN116706917A (zh) 基于快速交替方向乘子法的智慧园区协同调控方法及系统
Hong et al. A dynamic demand-driven smart manufacturing for mass individualization production
CN113467405B (zh) 工业4.0需求驱动的分布式动态边缘云智能制造方法及系统
CN117557016A (zh) 一种基于深度强化学习的整车制造冲压资源调度方法
Gu et al. Dynamic scheduling mechanism for intelligent workshop with deep reinforcement learning method based on multi-agent system architecture
CN112465391B (zh) 一种基于博弈论的分布式智能工厂供应任务分配方法
CN108090655A (zh) 用于服装生产的智能分配工序方法、模块及系统
CN112465394B (zh) 一种用于工业4.0大规模个性化生产的动态云制造方法
CN112508478B (zh) 一种基于自组织自动导引车的柔性物流配送任务分配方法
CN116993063A (zh) 一种任务求解方法及其装置
CN106408130B (zh) 一种用于电力系统经济调度的分布式拍卖方法
CN118261332B (zh) 基于数字化实现工业制造的场景建模方法及系统
Gao et al. An Integrated Two-Stage Optimization Method for Job-Shop Bottleneck Planning and Scheduling

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