CN104573846B - 一种基于ca‑pso混合优化算法的多态性作业车间布局优化方法 - Google Patents
一种基于ca‑pso混合优化算法的多态性作业车间布局优化方法 Download PDFInfo
- Publication number
- CN104573846B CN104573846B CN201410741912.4A CN201410741912A CN104573846B CN 104573846 B CN104573846 B CN 104573846B CN 201410741912 A CN201410741912 A CN 201410741912A CN 104573846 B CN104573846 B CN 104573846B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- cellular
- spacing
- msubsup
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/043—Optimisation of two dimensional placement, e.g. cutting of clothes or wood
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
Abstract
一种基于CA‑PSO混合优化算法的多态性作业车间布局优化方法，包括以下步骤：1)建车间的机器元胞、辅设元胞、通道元胞和间距元胞，2)粒子基于2N维向量的编码，3)通过SLP方法产生初始粒子，4)计算粒子的约束惩罚函数值，5)计算粒子的适应度值，6)设定迭代次数t，加速系数c1、c2，惯性系数w，7)元胞通过跟踪适应度值来更新自己的状态，每个粒子按照位置变换式进行优化，最终得到最优粒子。通过本方法所得到的车间布局方案可有效降低车间物流搬运总费用。
Description
技术领域
本发明涉及一种基于元胞机(Cellular Automata，以下简称CA)和改进粒子群算法(Particle Swarm Optimization，以下简称PSO)的混合优化算法的多态性作业车间布局优化方法，属于企业管理工程与信息技术领域。
背景技术
随着市场竞争的日益激烈，一类复杂产品大规模定制企业和产业集群下大型代工企业，其生产过程在订单、设备、工艺、批量和瓶颈等层次上越来越凸显形态和状态的多样性，布局的柔性、鲁棒性和动态性需求逐步增强，此类企业车间布局新问题被称为“多态性作业车间布局”。
目前应用于车间布局研究的智能优化算法有遗传算法(GA)、模拟退火算法(SA)、禁忌搜索算法(TS)等，但是将这些算法运用于求解多态性作业车间布局问题的时候，存在着计算量大、收敛慢、效率低、寻优能力不足的问题。
发明内容
为克服现有技术在解决多态性作业车间布局存在的不足，本发明提出一种基于CA-PSO混合优化算法的多态性作业车间布局优化方法，使其在解决多态性作业车间布局时能够寻求最佳的设备布置，从而得到比现有布局方法更适应的优化方法。
本发明解决其技术问题所采用的技术方案是：
一种基于CA-PSO混合优化算法的多态性作业车间布局优化方法，所述优化方法包括以下步骤：
1.一种基于CA-PSO混合优化算法的多态性作业车间布局优化方法，其特征在于：所述布局优化方法包含以下步骤：
(1)建立多态性作业车间布局的元胞机模型。对车间(L×W，L为车间长度，W为车间宽度)当中的机器、辅助设备、通道、机器间的间距进行抽象化，建立机器元胞、辅设元胞、通道元胞和间距元胞如下：
机器元胞t时刻的状态属性表示为：
式中：
mp——机器元胞所在网格中的位置(i，j)；
mn——机器的名称；
mi——机器布置所需要的间距；
ml——机器长度；
mw——机器宽度；
ms——机器布置状态，ms∈{0,1}，0表示横向摆放，1表示纵向摆放。
辅设元胞t时刻的状态属性表示为：
式中：
ap——辅设元胞所在网格中的位置(i，j)；
an——辅助设备的名称；
ai——辅助设备布置所需要的间距；
al——辅助设备长度；
aw——辅助设备宽度；
as——辅助设备布置状态，as∈{0,1}，0表示横向摆放，1表示纵向摆放。通道元胞t时刻的状态属性表示为：
式中：
ps——通道元胞所在网格中的位置(i，j)；
pm——主通道大小；
ps——次通道大小；
pa——通道元胞属性，pa∈{0,1}，0表示主通道，1表示次通道。
间距元胞t时刻的状态属性表示为：
式中：
Ss——间距元胞所在网格中的位置(i，j)；
Sf——间距元胞在网格(i+1，j)方向的间距，Sf∈{al,aw,pm,ps}；
Sb——间距元胞在网格(i-1，j)方向的间距，Sf∈{al,aw,pm,ps}；
Sl——间距元胞在网格(i，j+1)方向的间距，Sf∈{al,aw,pm,ps}；
Sr——间距元胞在网格(i，j-1)方向的间距，Sf∈{al,aw,pm,ps}。
(2)编码。假设车间布局元胞机模型有N个元胞(包括机器元胞、辅设元胞、通道元胞和间距元胞)，把每个布局候选方案表示为一个粒子，则每个粒子的编码是一个2N维的向量：
P＝(x1,x2,…,xN，y1,y2,…,yN) (5)
式中，前N维表示每个元胞的x坐标，后N维表示每个元胞的y坐标。
每个粒子的飞行方向和速率也是一个2N维的向量：
V＝(vx1,vx2,…,vxN，yx1,yx2,…,yxN) (6)
式中，前N维表示每个元胞的x方向的移动速度，后N维表示每个元胞的y方向的移动速度。
(3)产生初始粒子。对车间待布置的元胞通过使用SLP法(System LayoutPlanning，系统布置设计)进行分析计算，得到初步的布局方案，作为初始粒子。
(4)按照式(7)计算粒子的约束惩罚值。
式中，δ1、δ2为惩罚系数；如果惩罚函数值为0时，则该粒子满足约束条件。
(5)按照式(8)计算粒子的适应度值。
式中，Q＝{1，2，…，q}表示车间单位生产周期内所加工的零件类型集合；表示零件q生产过程在元胞n和元胞m之间每单位距离的物料搬运费用；表示零件q生产过程在元胞n和元胞m之间物料搬运频率；Dnm表示元胞n和元胞m之间的距离。
(6)设定迭代次数t，加速系数c1、c2，惯性系数w。
迭代次数t作为该优化方法的终止条件；加速系数c1、c2是PSO算法中最重要的参数之一。它反映粒子在寻优过程中受个体和全局信息影响的程度，它对平衡飞行粒子的局部搜索能力和全局搜索能力起到至关重要的作用。加速系数的设定在不同情况下会产生不同结果，是针对某一特定问题而设定的；惯性系数w也是PSO算法中最重要的参数之一，它的大小决定了粒子飞行的新速度对原来速度继承保留的程度，本发明中采用线性递减的惯性系数，通过不断降低原有速度所包含的梯度信息，以减小局部极值在优化初期对算法的影响。
(7)优化。在优化过程中，元胞通过跟踪适应度值来更新自己的状态，每个粒子的位置按式(9)～式(12)进行变换，对适应度值进行比较，挑选出其中最优的粒子，即最优的车间布局方案。
Vxi(t+1)＝w·Vxi(t)+c1·rand()·[pxi(t)-xi(t)]+c2·rand()·[pxg(t)-xi(t)]
i＝1,2…,N； (9)
Vyi(t+1)＝w·Vyi(t)+c1·rand()·[pyi(t)-yi(t)]+c2·rand()·[pyg(t)-yi(t)]
i＝1,2,…,N； (10)
xi(t+1)＝xi(t)+Vxi(t+1) i＝1,2,…N； (11)
yi(t+1)＝yi(t)+Vyi(t+1) i＝1,2,…N； (12)
式(9)～式(12)中：
t——迭代次数；
c1、c2——加速系数；
rand()——均匀分布在(0，1)间的随机数；
w——惯性系数；
pxi——第i个元胞经历的x方向的最好位置；
pyi——第i个元胞经历的y方向的最好位置；
pxg——所有元胞经历过的x方向的最好位置；
pyg——所有元胞经历过的y方向的最好位置。
本发明的技术构思：由于多态性作业车间待布置的各种资源繁多，而且车间几何约束条件、物流关系和功能关系复杂，因此对车间布局的智能优化算法要求更高，要求能够解决大规模的布局问题，在保证不存在解早熟和陷入局部最优解的前提下，能够快速地寻找到全局最优解。
本发明正是利用CA-PSO算法强大的计算能力来解决多态性作业车间问题。元胞机是一个时间与空间都离散的动力系统，由有限个状态的元胞组成的离散空间，具有同步更新的特点；粒子群算法是一个无约束的优化算法，具有高效的计算的能力。CA-PSO算法思想：通过元胞机来抽象构建多态性作业车间的各资源元胞，再通过粒子群算法寻求全局最优解，在寻优过程中，各资源元胞状态同步更新。
本发明的实质效果体现在：1.CA-PSO优化算法具有收敛速度快，更容易趋向全局最优解的特点；2.基于CA-PSO优化算法所得到的多态性作业车间布局方案能够有效降低车间的物流费用，提高车间面积利用率。
附图说明
图1为多态性作业车间元胞机模型图
图2为本发明实施例车间布局图
具体实施方式
下面结合附图和实施例，具体说明本发明的内容。
某汽车尾管生产车间面积为35m×30m，待布置的资源有线切割机、冲床等见表1所示；各机器间的物流量见表2。机器的最小间距mi＝0.9m；辅设的最小间距ai＝0.8m；各设备之间物料搬运均采用手推车，搬运单价为1.5×10-4元/米×件。
表1 设备的尺寸和数量
表2 设备间物流量 (单位：件)
1)构建机器元胞45个、辅设元胞11个、通道元胞30个和间距元胞56个。
机器元胞t时刻的状态属性表示为：
式中：
mp——机器元胞所在网格中的位置(i，j)；
mn——机器的名称；
mi——机器布置所需要的间距；
ml——机器长度；
mw——机器宽度；
ms——机器布置状态，(ms∈{0,1})，0表示横向摆放，1表示纵向摆放。
辅设元胞t时刻的状态属性表示为：
式中：
ap——辅设元胞所在网格中的位置(i，j)；
an——辅助设备的名称；
ai——辅助设备布置所需要的间距；
al——辅助设备长度；
aw——辅助设备宽度；
as——辅助设备布置状态，(as∈{0,1})，0表示横向摆放，1表示纵向摆放。
通道元胞t时刻的状态属性表示为
式中：
ps——通道元胞所在网格中的位置(i，j)；
pm——主通道大小；
ps——次通道大小；
pa——通道元胞属性，(pa∈{0,1})，0表示主通道，1表示次通道。
间距元胞t时刻的状态属性表示为：
式中：
Ss——间距元胞所在网格中的位置(i，j)；
Sf——间距元胞在网格(i+1，j)方向的间距，Sf∈{al,aw,pm,ps}；
Sb——间距元胞在网格(i-1，j)方向的间距，Sf∈{al,aw,pm,ps}；
Sl——间距元胞在网格(i，j+1)方向的间距，Sf∈{al,aw,pm,ps}；
Sr——间距元胞在网格(i，j-1)方向的间距，Sf∈{al,aw,pm,ps}。
2)编码。车间布局元胞机模型有142个元胞，把每个布局候选方案表示为一个粒子，则每个粒子的编码是一个284维的向量：
P＝(x1,x2,…,x142，y1,y2,…,y142) (5)
式中，前142维表示每个元胞的x坐标，后142维表示每个元胞的y坐标。
每个粒子的飞行方向和速率也是一个284维的向量：
V＝(vx1,vx2,…,vx142，yx1,yx2,…,yx142) (6)
式中，前142维表示每个元胞的x方向的移动速度，后142维表示每个元胞的y方向的移动速度。
3)产生初始粒子。对车间待布置的元胞通过使用SLP法(System LayoutPlanning，系统布置设计)进行分析计算，得到初步的布局方案，作为初始粒子。
4)按照式(7)计算粒子的约束惩罚函数值。
取惩罚系数δ1、δ2为0.5。
5)按照式(8)计算粒子的适应度值。
式中，Q＝{1，2，…，q}表示车间单位生产周期内所加工的零件类型集合；表示零件q生产过程在元胞n和元胞m之间每单位距离的物料搬运费用；表示零件q生产过程在元胞n和元胞m之间物料搬运频率；Dnm表示元胞n和元胞m之间的距离。
6)设定迭代次数t＝100，加速系数c1＝1.2、c2＝1.2，惯性系数w＝0.5。
7)优化。在优化过程中，元胞通过跟踪适应度值来更新自己的状态，每个粒子的位置按式(9)～式(12)进行变换，对适应度值进行比较，挑选出其中最优的粒子。
Vxi(t+1)＝w·Vxi(t)+c1·rand()·[pxi(t)-xi(t)]+c2·rand()·[pxg(t)-xi(t)]
i＝1,2…,N； (9)
Vyi(t+1)＝w·Vyi(t)+c1·rand()·[pyi(t)-yi(t)]+c2·rand()·[pyg(t)-yi(t)]
i＝1,2,…,N； (10)
xi(t+1)＝xi(t)+Vxi(t+1) i＝1,2,…N； (11)
yi(t+1)＝yi(t)+Vyi(t+1) i＝1,2,…N； (12)
式(9)～式(12)中：
t——迭代次数；
c1、c2——加速系数；
rand()——均匀分布在(0，1)间的随机数；
w——惯性系数；
pxi——第i个元胞经历的x方向的最好位置；
pyi——第i个元胞经历的y方向的最好位置；
pxg——所有元胞经历过的x方向的最好位置；
pyg——所有元胞经历过的y方向的最好位置。
经上述具体实施步骤后，本实施例车间布局图如图2所示，计算得到的总物流费用为9860000元，而本车间未采用本布局优化方法前的总物流费用1250000元。对于待布置资源繁多的多态性作业车间，应用本发明方法所得到的布局方案可以明显降低物流总费用。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (1)
1.一种基于CA-PSO混合优化算法的多态性作业车间布局优化方法，其特征在于：所述布局优化方法包含以下步骤：
(1)建立多态性作业车间布局的元胞机模型；对车间当中的机器、辅助设备、通道、机器间的间距进行抽象化，车间尺寸是L×W，L为车间长度，W为车间宽度，建立机器元胞、辅设元胞、通道元胞和间距元胞如下：
机器元胞t时刻的状态属性表示为：
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>t</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mi>p</mi>
</msub>
<mo>,</mo>
<msub>
<mi>m</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>m</mi>
<mi>l</mi>
</msub>
<mo>,</mo>
<msub>
<mi>m</mi>
<mi>w</mi>
</msub>
<mo>,</mo>
<msub>
<mi>m</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：
mp——机器元胞所在网格中的位置(i，j)；
mn——机器的名称；
mi——机器布置所需要的间距；
ml——机器长度；
mw——机器宽度；
ms——机器布置状态，ms∈{0,1}，0表示横向摆放，1表示纵向摆放；
辅设元胞t时刻的状态属性表示为：
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>t</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>p</mi>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>l</mi>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>w</mi>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：
ap——辅设元胞所在网格中的位置(i，j)；
an——辅助设备的名称；
ai——辅助设备布置所需要的间距；
al——辅助设备长度；
aw——辅助设备宽度；
as——辅助设备布置状态，as∈{0,1}，0表示横向摆放，1表示纵向摆放；
通道元胞t时刻的状态属性表示为：
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>t</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：
ps——通道元胞所在网格中的位置(i，j)；
pm——主通道大小；
ps——次通道大小；
pa——通道元胞属性，pa∈{0,1}，0表示主通道，1表示次通道；
间距元胞t时刻的状态属性表示为：
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>t</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mi>S</mi>
<mi>f</mi>
</msub>
<mo>,</mo>
<msub>
<mi>S</mi>
<mi>b</mi>
</msub>
<mo>,</mo>
<msub>
<mi>S</mi>
<mi>l</mi>
</msub>
<mo>,</mo>
<msub>
<mi>S</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：
Ss——间距元胞所在网格中的位置(i，j)；
Sf——间距元胞在网格(i+1，j)方向的间距，Sf∈{al,aw,pm,ps}；
Sb——间距元胞在网格(i-1，j)方向的间距，Sf∈{al,aw,pm,ps}；
Sl——间距元胞在网格(i，j+1)方向的间距，Sf∈{al,aw,pm,ps}；
Sr——间距元胞在网格(i，j-1)方向的间距，Sf∈{al,aw,pm,ps}；
(2)编码；假设车间布局元胞机模型有N个元胞，包括机器元胞、辅设元胞、通道元胞和间距元胞，把每个布局候选方案表示为一个粒子，则每个粒子的编码是一个2N维的向量：
P＝(x1,x2,…,xN，y1,y2,…,yN) (5)
式中，前N维表示每个元胞的x坐标，后N维表示每个元胞的y坐标；
每个粒子的飞行方向和速率也是一个2N维的向量：
V＝(vx1,vx2,…,vxN，yx1,yx2,…,yxN) (6)
式中，前N维表示每个元胞的x方向的移动速度，后N维表示每个元胞的y方向的移动速度；
(3)产生初始粒子；对车间待布置的元胞通过使用SLP法进行分析计算，得到初步的布局方案，作为初始粒子；
(4)按照式(7)计算粒子的约束惩罚值；
<mrow>
<mi>P</mi>
<mo>=</mo>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mo>&times;</mo>
<mo>&lsqb;</mo>
<mi>L</mi>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>m</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mi>l</mi>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>a</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>l</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>&times;</mo>
<mo>&lsqb;</mo>
<mi>W</mi>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>m</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mi>w</mi>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>a</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>w</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，δ1、δ2为惩罚系数；如果惩罚函数值为0时，则该粒子满足约束条件；
(5)按照式(8)计算粒子的适应度值；
<mrow>
<mi>F</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</msubsup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
<mi>q</mi>
</msubsup>
<msubsup>
<mi>F</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
<mi>q</mi>
</msubsup>
<msub>
<mi>D</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，Q＝{1，2，…，q}表示车间单位生产周期内所加工的零件类型集合；表示零件q生产过程在元胞n和元胞m之间每单位距离的物料搬运费用；表示零件q生产过程在元胞n和元胞m之间物料搬运频率；Dnm表示元胞n和元胞m之间的距离；
(6)设定迭代次数t，加速系数c1、c2，惯性系数w；
(7)优化；在优化过程中，元胞通过跟踪适应度值来更新自己的状态，每个粒子的位置按式(9)～式(12)进行变换，对适应度值进行比较，挑选出其中最优的粒子，即最优的车间布局方案；
Vxi(t+1)＝w·Vxi(t)+c1·rand()·[pxi(t)-xi(t)]+c2·rand()·[pxg(t)-xi(t)]i＝1,2…,N； (9)
Vyi(t+1)＝w·Vyi(t)+c1·rand()·[pyi(t)-yi(t)]+c2·rand()·[pyg(t)-yi(t)]i＝1,2,…,N； (10)
xi(t+1)＝xi(t)+Vxi(t+1) i＝1,2,…N； (11)
yi(t+1)＝yi(t)+Vyi(t+1) i＝1,2,…N； (12)
式(9)～式(12)中：
t——迭代次数；
c1、c2——加速系数；
rand()——均匀分布在(0，1)间的随机数；
w——惯性系数；
pxi——第i个元胞经历的x方向的最好位置；
pyi——第i个元胞经历的y方向的最好位置；
pxg——所有元胞经历过的x方向的最好位置；
pyg——所有元胞经历过的y方向的最好位置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410741912.4A CN104573846B (zh) | 2014-12-08 | 2014-12-08 | 一种基于ca‑pso混合优化算法的多态性作业车间布局优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410741912.4A CN104573846B (zh) | 2014-12-08 | 2014-12-08 | 一种基于ca‑pso混合优化算法的多态性作业车间布局优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104573846A CN104573846A (zh) | 2015-04-29 |
| CN104573846B true CN104573846B (zh) | 2017-09-26 |
Family
ID=53089861
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410741912.4A Active CN104573846B (zh) | 2014-12-08 | 2014-12-08 | 一种基于ca‑pso混合优化算法的多态性作业车间布局优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104573846B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107256434B (zh) * | 2017-04-24 | 2021-03-02 | 深圳市蜗牛窝科技有限公司 | 家居自动布局的方法 |
| CN109543896B (zh) * | 2018-11-16 | 2021-07-27 | 闽江学院 | 一种基于绝对值和函数最小化与粒子群优化的车间设备双行布局方法 |
| TWI828512B (zh) * | 2023-01-10 | 2024-01-01 | 力晶積成電子製造股份有限公司 | 搬運排程方法及裝置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102608916A (zh) * | 2012-02-15 | 2012-07-25 | 浙江工业大学 | 一种基于元胞机的大型零件柔性作业车间的动态调度方法 |
-
2014
- 2014-12-08 CN CN201410741912.4A patent/CN104573846B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102608916A (zh) * | 2012-02-15 | 2012-07-25 | 浙江工业大学 | 一种基于元胞机的大型零件柔性作业车间的动态调度方法 |
Non-Patent Citations (2)
| Title |
|---|
| "基于差分元胞多目标遗传算法的车间布局优化";张屹 等;《计算机集成制造系统》;20130415;第19卷(第4期);第727-734页 * |
| "基于粒子群优化算法的车间布局设计";郭源源 等;《计算机集成制造系统》;20120301;第18卷(第11期);第2476-2484页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104573846A (zh) | 2015-04-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105488562B (zh) | 一种基于多因子粒子群算法的不规则件排样方法 | |
| CN105243430B (zh) | 列车节能运行的目标速度曲线的优化方法 | |
| CN103020358B (zh) | 针对机械装备的可适应性动态设计平台构建方法 | |
| CN104573846B (zh) | 一种基于ca‑pso混合优化算法的多态性作业车间布局优化方法 | |
| CN102622468B (zh) | 基于并行计算的大规模集成电路通道布线系统 | |
| CN102073311A (zh) | 基于离散量子微粒群算法的机器零件加工流水线调度方法 | |
| CN103488537A (zh) | 一种数据抽取、转换和加载etl的执行方法及装置 | |
| CN105858043B (zh) | 一种升降机与穿梭车结合的仓储系统优化调度方法 | |
| CN104063745A (zh) | 基于改进粒子群算法的多路径规划方法 | |
| CN112446642A (zh) | 一种多天车调度优化方法和系统 | |
| CN105373668B (zh) | 芯片版图设计方法 | |
| CN105447284B (zh) | 一种工艺余量的添加方法 | |
| Chen et al. | Particle swarm optimization based on genetic operators for sensor-weapon-target assignment | |
| CN101739713A (zh) | 一种基于剖分的并行碰撞检测方法 | |
| CN102222353A (zh) | 一种基于二次b样条迭代的曲线绘制方法 | |
| CN102692491A (zh) | 基于分阶段禁忌寻优算法的土壤水分特征参数计算方法 | |
| Zhang et al. | Digital twin-based production simulation of discrete manufacturing shop-floor for onsite performance analysis | |
| CN103984985A (zh) | 一种实现定制家具继列分批调度优化的方法 | |
| CN104346441B (zh) | 一种配电网信息数据动态集成交互方法 | |
| Huang et al. | Obstacle-avoiding octagonal Steiner tree construction based on particle swarm optimization | |
| Gong | Application of artificial intelligence computer intelligent heuristic search algorithm | |
| Zhang et al. | An improved genetic algorithm for multi-objective flexible job-shop scheduling problem | |
| CN107577894B (zh) | 一种四边形单元实时二分式加密-减密方法 | |
| CN105867131A (zh) | 电动汽车的控制效率分配方法及装置 | |
| CN102194033A (zh) | 一种变型设计方法 |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20191206 Address after: 313000 No.5-8, Changhe Road, Huaxi street, Changxing County, Huzhou City, Zhejiang Province Patentee after: Zhejiang Changxing green battery technology Co.,Ltd. Address before: 313000 Room 1403, 14th Floor, Building B, Freeport, Headquarters 1188 District Fulu Road, Wuxing District, Huzhou City, Zhejiang Province Patentee before: Zhejiang creation Intellectual Property Service Co.,Ltd. Effective date of registration: 20191206 Address after: 313000 Room 1403, 14th Floor, Building B, Freeport, Headquarters 1188 District Fulu Road, Wuxing District, Huzhou City, Zhejiang Province Patentee after: Zhejiang creation Intellectual Property Service Co.,Ltd. Address before: 310014 Hangzhou city in the lower reaches of the city of Zhejiang Wang Road, No. 18 Patentee before: Zhejiang University of Technology |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201104 Address after: 313000 no.55-9 Huyan South Road, Shuanglin Town, Nanxun District, Huzhou City, Zhejiang Province (self declaration) Patentee after: Huzhou Duxin Technology Co.,Ltd. Address before: 313000 No.5-8, Changhe Road, Huaxi street, Changxing County, Huzhou City, Zhejiang Province Patentee before: Zhejiang Changxing green battery technology Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220106 Address after: 310000 room 205, building 2, shijiajun block, Xiacheng District, Hangzhou City, Zhejiang Province Patentee after: Hangzhou Hongren Technology Co.,Ltd. Address before: 313000 no.55-9 Huyan South Road, Shuanglin Town, Nanxun District, Huzhou City, Zhejiang Province (self declaration) Patentee before: Huzhou Duxin Technology Co.,Ltd. |
|
| TR01 | Transfer of patent right |