CN103336855B - A kind of two-dimentional irregular nesting method based on many subgroups particle cluster algorithm - Google Patents
A kind of two-dimentional irregular nesting method based on many subgroups particle cluster algorithm Download PDFInfo
- Publication number
- CN103336855B CN103336855B CN201310199780.2A CN201310199780A CN103336855B CN 103336855 B CN103336855 B CN 103336855B CN 201310199780 A CN201310199780 A CN 201310199780A CN 103336855 B CN103336855 B CN 103336855B
- Authority
- CN
- China
- Prior art keywords
- particle
- subgroup
- rand
- solution
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000001788 irregular Effects 0.000 title claims abstract description 15
- 239000000463 material Substances 0.000 claims abstract description 32
- 230000000694 effects Effects 0.000 claims description 4
- 239000000470 constituent Substances 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000002922 simulated annealing Methods 0.000 description 2
- 238000000137 annealing Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Abstract
A kind of two-dimentional irregular nesting method based on improved Particle Swarm Algorithm, comprise the following steps: the first step, the geometric figure of print and material is converted to a series of two-dimensional coordinate interval, then use the left searching algorithm in the heuristic end to judge the two dimension interval of print and material whether the overlapping print that moves is relative to the position in material; Second step, modified PSO search procedure: cross the method dividing multiple subgroup, when not changing the parameter tending to current optimum solution, add the impact of subgroup optimum solution on particle in subgroup, during the maximum iteration time set when the iterations of population reaches initial, particle stops iteration, obtains current globally optimal solution as final layout project.The invention provides a kind of there is good search capability while, search speed is fast, last solution is better, stock layout the is respond well two-dimentional irregular nesting method based on improved Particle Swarm Algorithm.
Description
Technical field
The present invention relates to Computer Aided Nesting technology, especially a kind of Nesting.
Background technology
At present, both at home and abroad the Nesting of irregular part is concentrated in didactic modern optimization method substantially, mainly contain simulated annealing, genetic algorithm, particle cluster algorithm (PSO) etc.In practice, because simulated annealing is by the impact of annealing speed, speed is fast, is easily absorbed in local extremum, and speed is then difficult to the needs that can not meet people slowly; Genetic algorithm is then often than being easier to " precocity ", and local search ability is poor; PSO algorithm has stronger local search ability, but is also easily absorbed in extremal region and is difficult to jump out.
The range of application of Cutting Stock Problem widely, as panel beating stock layout, glass tailor, cloth-cutting etc., because above-mentioned industry turnout is large, particularly some material is costly, and therefore, the raising of Nesting Algorithms efficiency can produce very large social benefit.
Summary of the invention
In order to overcome the deficiency being easily absorbed in extremal region, have impact on stock layout effect of existing Nesting, the invention provides a kind of there is good search capability while, search speed is fast, last solution is better, stock layout the is respond well two-dimentional irregular nesting method based on many subgroups particle cluster algorithm.
The technical solution adopted for the present invention to solve the technical problems is:
Based on a two-dimentional irregular nesting method for many subgroups particle cluster algorithm, comprise the following steps:
The first step, is converted to a series of two-dimensional coordinate by the geometric figure of print and material interval, then use the left searching algorithm in the heuristic end to judge the two dimension interval of print and material whether the overlapping print that moves is relative to the position in material;
Second step, modified PSO search procedure is as follows:
1) inverse of the height after entering using print is as fitness value, and fitness value 1/H is larger, enters effect better;
2) state that enters of each print mainly contains three kinds: enter order, the anglec of rotation and mirror image, described in enter the span 1 ~ n of order order, n is print sum; The span of anglec of rotation angle 0 ° ~ 360 °, whether symmetrical about y-axis mirror image mirror represents, 0 ~ 1,0 expression is not in relation to y-axis symmetry, and 1 represents symmetrical about y-axis;
3) using 2) in three parameters proposing as three elements of the elementary particle in constituent particle group, elementary particle described in random initializtion;
4) calculate the European geometric position of each particle, according to the order from small to large from initial point distance, all particles are divided into M subgroup, M < n;
5) each optimum configurations is as follows:
x
ij=〈order
ij,angle
ij,mirror
ij〉
The position vector of a jth particle in i-th subgroup;
v
ij=〈v_ord
ij,v_ang
ij,v_mir
ij〉
The velocity vector of a jth particle in i-th subgroup;
p
ij=〈p_ord
ij,p_ang
ij,p_mir
ij〉
The history optimum position vector of a jth particle in i-th subgroup;
psg
i=〈psg_ord
i,psg_ang
i,psg_mir
i〉
History optimum position, i-th subgroup vector;
p
g=〈pg
ord,pg
ang,pg
miri〉
Global history optimum position vector;
6) speed of particle and location updating formula as follows:
v
ij(d+1)=w×v
ij(d)+c
1×rand
1ij×[p
ij(d)-x
ij(d)]
+c
2×rand
2ij×[psg
i(d)-x
ij(d)]
+c
3×rand
3ij×[p
g(d)-x
ij(d)]
x
ij(d+1)=x
ij(d)+v
ij(d+1)
Wherein d is iterations, c
l, c
2, c
3represent the history optimum solution of trend particle own, the optimum solution of subgroup, the speeds control factor of globally optimal solution respectively, wherein c
3> c
2> c
l> 0, rand
1ij, rand
2ij, rand
3ijbe the random factor between 0 ~ 1, w is the value w(d of inertial factor, w) linear decrease with the increase of iterations:
w(d)=u-v×d/D
D is maximum iteration time, and the value of u, v meets;
7) more new historical optimum solution
After each particle carries out speed renewal, enter in material by the left searching algorithm in the heuristic end, calculate the fitness value F of particle, thus the history optimal location of more new particle itself, the optimal location of subgroup and the overall situation;
8), during the maximum iteration time set when the iterations of population reaches initial, particle stops iteration, obtains current globally optimal solution as final layout project; If do not reach maximum iteration time, continue to perform from step 6).
Further, in described step 7), when more excellent solution is not all found in certain subgroup in continuous n iterative process, judge that this subgroup has been absorbed in locally optimal solution, to the random replacement of described subgroup.
Further, in the described first step, it is as follows that the left searching algorithm in the heuristic end obtains process: print first enters from the lower left corner of material, print is below turned right and is entered successively, when print exceeds the right side of material, print is moved relative on material, and again search for the remaining space of material from the left side, so move in circles, until stop when print all enters or print overflows material top.
Technical conceive of the present invention is: by dividing the method for multiple subgroup, when not changing the parameter tending to current optimum solution, add the impact of subgroup optimum solution on particle in subgroup, local search ability is strengthened with this, and when subgroup is absorbed in local extremum, carry out eliminating to subgroup and reset, make colony have an opportunity to jump out local extremum.
Beneficial effect of the present invention is mainly manifested in: while having good search capability, search speed is fast, last solution is better, stock layout is respond well.
Accompanying drawing explanation
Fig. 1 is the program flow diagram of two-dimentional irregular nesting underlying algorithm.
Fig. 2 is the program flow diagram of modified PSO algorithm.
Fig. 3 is for certain cover fashion plate is based on stock layout design sketch of the present invention.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described.
With reference to Fig. 1 ~ Fig. 3, a kind of two-dimentional irregular nesting method based on many subgroups particle cluster algorithm, comprises the following steps:
The first step, the selection of two-dimentional irregular nesting underlying algorithm
For the typesetting of the two-dimentional irregular print of process, first the profile determining print is needed, to enter in process at print and to judge whether each print can be overlapping, the present invention selects the geometric figure of print and material to be converted to a series of two-dimensional coordinate interval, thus the complicacy departing from irregular geometric figures carries out stock layout, then use the heuristic end left searching algorithm (HBLS) judge the two dimension interval of print and material whether the overlapping print that moves is relative to the position in material, complete the underlying algorithm of two-dimentional irregular nesting with this, its flow process can refer to Fig. 1.
The basic ideas of HBLS algorithm are: print first enters from the lower left corner of material, print is below turned right and is entered successively, when print exceeds the right side of material, print is moved relative on material, and from the left side, again search for the remaining space of material, so move in circles, until stop when print all enters or print overflows material top.
Second step, modified PSO searching algorithm
1) print of layer with the end enters algorithm, the state parameter just entered by print by searching algorithm passes to underlying algorithm and processes, thus the fitness value F(the present invention calculating each individuality enter using print after the inverse (1/H) of height as fitness value, fitness value is larger, enters effect better).
2) state that enters of each print mainly contains three kinds: enter order (order, 1 ~ n, n are print sum), the anglec of rotation (angle, 0 ° ~ 360 °), mirror image (mirror, symmetrical about y-axis, 0 ~ 1)
3) using 2) in three parameters proposing as three elements of the elementary particle in constituent particle group, the sum of particle is set as 30 here, and these 30 particles of random initializtion.
4) calculate the European geometric position of each particle, according to the order from small to large from initial point distance, 30 particles are divided into 5 parts, every part of 6 particles, namely divide into 5 subgroups.Be to allow in original state each subgroup be distributed in zones of different as far as possible according to the object that European geometric distance divides, strengthen the ability of searching optimum of algorithm.
5) each optimum configurations is as follows:
x
ij=(order
ij,angle
ij,mirror
ij〉
The position vector of a jth particle in i-th subgroup
v
ij=〈v_ord
ij,v_ang
ij,v_mir
ij〉
The velocity vector of a jth particle in i-th subgroup
p
ij=〈p_ord
ij,p_ang
ij,p_mir
ij〉
The history optimum position vector of a jth particle in i-th subgroup
psg
i=〈psg_ord
i,psg_ang
i,psg_mir
i〉
History optimum position, i-th subgroup vector
p
g=<pg
ord,pg
ang,pg
miri>
Global history optimum position vector
6) speed of particle and location updating formula as follows:
v
ij(d+1)=w×v
ij(d)+c
1×rand
1ij,×[p
ij(d)-x
ij,(d)]
+c
2×rand
2ij×[psg
i(d)-x
ij(d)]
+c
3×rand
3ij×[p
g(d)-x
ij(d)]
x
ij(d+1)=x
ij(d)+v
ij(d+1)
Wherein d is iterations, it can thus be appreciated that the renewal of the every generation of particle is determined by the data of previous generation, c
l, c
2, c
3represent the history optimum solution of trend particle own, the optimum solution of subgroup, the speeds control factor of globally optimal solution respectively, wherein c
3> c
2> c
l> 0, rand
1ij, rand
2ij, rand
3ijbe the random factor between 0 ~ 1, w is inertial factor, and when w is larger, algorithm is to solution space extensive search; Otherwise algorithm is searched among a small circle, the value of w is linear decrease with the increase of iterations:
w(d)=u-v×d/D
D is maximum iteration time, and the value of u, v meets O < v < u < 1.
In classical PSO algorithm, the speed of particle more new formula is:
v
ij(d+1)=w×v
ij(d)+c′
1×rand
1ij×[p
ij(d)-x
ij(d)]
+c′
3×rand
3ij×[p
g(d)-x
ij(d)]
Contrast with the PSO algorithm parameter improved, have following relation:
c′
l=c
l+c
2
c′
3=c
3
Modified PSO algorithm adds the concept of subgroup, and in speed more new formula, add the impact of subgroup optimum solution, strengthens Local Search, but c
3do not change, thus ensure that algorithm convergence is not weakened in the trend of more excellent solution.
7) more new historical optimum solution
After each particle carries out speed renewal, enter in material by bottom Nesting Algorithms, calculate the fitness value F of particle, thus the history optimal location of more new particle itself, the optimal location of subgroup and the overall situation.
8) the superseded and reset mechanism of subgroup
The problem of locally optimal solution is easily absorbed in for PSO algorithm, eliminating and reset mechanism of the subgroup of subgroup is added in the present invention, when more excellent solution is not all found in certain subgroup in continuous n iterative process, so, think that this subgroup has been absorbed in locally optimal solution, now by the random replacement to subgroup, this subgroup can be made to have an opportunity to jump out locally optimal solution, thus enhance the ability of searching optimum of algorithm.
9), during the maximum iteration time set when the iterations of population reaches initial, particle stops iteration, obtains current globally optimal solution as final layout project.If do not reach maximum iteration time, continue to perform from step 6).
With reference to Fig. 3, stock layout design sketch of the present invention, (material width: 1200mmm, print number: 18, iterations: 600, stock layout height H: 901mm, stock utilization: 78.1%).
Claims (2)
1., based on a two-dimentional irregular nesting method for many subgroups particle cluster algorithm, it is characterized in that: comprise the following steps:
The first step, is converted to a series of two-dimensional coordinate by the geometric figure of print and material interval, and whether the overlapping print that moves is relative to the position in material for the two-dimensional coordinate interval of print and material then to use the left searching algorithm in the heuristic end to judge;
Second step, modified PSO search procedure is as follows:
1) inverse of the height after entering using print is as fitness value, and fitness value 1/H is larger, enters effect better;
2) state that enters of each print mainly contains three kinds: enter order, the anglec of rotation and mirror image, described in enter the span 1 ~ n of order order, n is print sum; The span of anglec of rotation angle 0 ° ~ 360 °, whether symmetrical about y-axis mirror image mirror represents, 0 ~ 1,0 expression is not in relation to y-axis symmetry, and 1 represents symmetrical about y-axis;
3) using 2) in three parameters proposing as three elements of the elementary particle in constituent particle group, elementary particle described in random initializtion;
4) calculate the European geometric position of each particle, according to the order from small to large from initial point distance, all particles are divided into M subgroup, M < n;
5) each optimum configurations is as follows:
x
ij=<order
ij,angle
ij,mirror
ij>
The position vector of a jth particle in i-th subgroup;
v
ij=<v_ord
ij,v_ang
ij,v_mir
ij>
The velocity vector of a jth particle in i-th subgroup;
p
ij=<p_ord
ij,p_ang
ij,p_mir
ij>
The history optimum position vector of a jth particle in i-th subgroup;
psg
i=<psg_ord
i,psg_ang
i,psg_mir
i>
History optimum position, i-th subgroup vector;
p
g=<pg
ord,pg
ang,pg
miri>
Global history optimum position vector;
6) speed of particle and location updating formula as follows:
v
ij(d+1)=w×v
ij(d)+c
1×rand
1ij×[p
ij(d)-x
ij(d)]
+c
2×rand
2ij×[psg
i(d)-x
ij(d)]
+c
3×rand
3ij×[p
g(d)-x
ij(d)]
x
ij(d+1)=x
ij(d)+v
ij(d+1)
Wherein d is iterations, c
1, c
2, c
3represent the history optimum solution of trend particle own, the optimum solution of subgroup, the speeds control factor of globally optimal solution respectively, wherein c
3> c
2> c
1> 0, rand
1ij, rand
2ij, rand
3ijbe the random factor between 0 ~ 1, w is inertial factor, and the value w (d) of w is linear decrease with the increase of iterations:
w(d)=u-v×d/D
D is maximum iteration time, and the value of u, v meets 0 < v < u < 1;
7) more new historical optimum solution
After each particle carries out speed renewal, enter in material by the left searching algorithm in the heuristic end, calculate the fitness value F of particle, thus the history optimal location of more new particle itself, the optimal location of subgroup and the overall situation;
When more excellent solution is not all found in certain subgroup in continuous n iterative process, judge that this subgroup has been absorbed in locally optimal solution, described subgroup is reset at random;
8), during the maximum iteration time set when the iterations of population reaches initial, particle stops iteration, obtains current globally optimal solution as final layout project; If do not reach maximum iteration time, continue from step 6) perform.
2. a kind of two-dimentional irregular nesting method based on many subgroups particle cluster algorithm as claimed in claim 1, it is characterized in that: in the described first step, it is as follows that the left searching algorithm in the heuristic end obtains process: print first enters from the lower left corner of material, print is below turned right and is entered successively, when print exceeds the right side of material, print is moved relative on material, and from the left side, again search for the remaining space of material, so move in circles, until stop when print all enters or print overflows material top.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310199780.2A CN103336855B (en) | 2013-05-24 | 2013-05-24 | A kind of two-dimentional irregular nesting method based on many subgroups particle cluster algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310199780.2A CN103336855B (en) | 2013-05-24 | 2013-05-24 | A kind of two-dimentional irregular nesting method based on many subgroups particle cluster algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103336855A CN103336855A (en) | 2013-10-02 |
CN103336855B true CN103336855B (en) | 2016-04-13 |
Family
ID=49245019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310199780.2A Active CN103336855B (en) | 2013-05-24 | 2013-05-24 | A kind of two-dimentional irregular nesting method based on many subgroups particle cluster algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103336855B (en) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678800A (en) * | 2013-12-09 | 2014-03-26 | 华南理工大学 | Method for stock layout of irregular parts on rectangular board |
CN105488562B (en) * | 2015-11-27 | 2017-11-10 | 浙江工业大学义乌科学技术研究院有限公司 | A kind of layout of polygon method based on multiple-factor particle cluster algorithm |
CN105740953B (en) * | 2016-01-29 | 2018-01-12 | 温州大学 | A kind of irregular nesting method based on Real-coded quantum evolutionary algorithm |
CN105956663A (en) * | 2016-03-04 | 2016-09-21 | 安徽工程大学 | Parallel stock layout optimization method of special shaped part blanking |
CN105989605A (en) * | 2016-03-04 | 2016-10-05 | 安徽工程大学 | Irregular part blanking layout positioning method |
CN105845427B (en) * | 2016-04-08 | 2017-08-29 | 辽宁石油化工大学 | Cross section of transformer core design method based on particle group optimizing |
CN106156534B (en) * | 2016-09-10 | 2019-01-29 | 天津大学 | A kind of metal plate batch Nesting of the special container based on Creo |
CN106571610B (en) * | 2016-10-29 | 2017-12-12 | 南京理工大学 | A kind of transmission line of electricity mechanical deicing order optimization method based on particle cluster algorithm |
CN106959608B (en) * | 2017-02-27 | 2019-10-01 | 同济大学 | A kind of water supply network leakage optimal control method based on cluster particle swarm algorithm |
CN109214557A (en) * | 2018-08-15 | 2019-01-15 | 浙江优尼家装饰科技有限公司 | A kind of irregular plate Nesting of two dimension |
CN111275243B (en) * | 2020-01-08 | 2020-09-29 | 广东工业大学 | Two-dimensional irregular layout method for single-specification plate |
CN111488686B (en) * | 2020-04-14 | 2023-07-11 | 常州纳捷机电科技有限公司 | Cross stripe cloth position matching layout method based on genetic algorithm |
CN113189553A (en) * | 2021-01-28 | 2021-07-30 | 北京理工大学 | Multi-modal optimization-based micro-maneuvering multi-target rapid detection method |
CN113032921B (en) * | 2021-03-16 | 2022-12-13 | 山东科技大学 | Layout algorithm based on parallel adaptive parameter cuckoo search and lowest horizontal line |
CN113762607B (en) * | 2021-08-26 | 2023-05-30 | 甘肃同兴智能科技发展有限责任公司 | Prediction method for carbon emission of power grid enterprise |
CN115629537A (en) * | 2022-12-22 | 2023-01-20 | 南京净环热冶金工程有限公司 | Heating furnace combustion control method and system based on subgroup improved particle swarm optimization PID |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102313989A (en) * | 2011-09-13 | 2012-01-11 | 山东理工大学 | A kind of method of single non-spherical lens aberration correction |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7636700B2 (en) * | 2004-02-03 | 2009-12-22 | Hrl Laboratories, Llc | Object recognition system incorporating swarming domain classifiers |
-
2013
- 2013-05-24 CN CN201310199780.2A patent/CN103336855B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102313989A (en) * | 2011-09-13 | 2012-01-11 | 山东理工大学 | A kind of method of single non-spherical lens aberration correction |
Non-Patent Citations (2)
Title |
---|
智能数控裁床的研究与开发;黄建江;《中国优秀硕士学位论文全文数据库》;20090315(第3期);论文第7页第2.1、2.2节,第11页第2.3节-第14页第2.3.3节,第18页第3.1节-第31页第4.2.3节 * |
纵向参数多子群粒子群算法的研究与应用;常彦伟;《中国博士学位论文全文数据库》;20110515(第5期);论文第24页第2.10.1节-第28页第2.10.6节 * |
Also Published As
Publication number | Publication date |
---|---|
CN103336855A (en) | 2013-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103336855B (en) | A kind of two-dimentional irregular nesting method based on many subgroups particle cluster algorithm | |
CN102520718B (en) | Physical modeling-based robot obstacle avoidance path planning method | |
CN106985395B (en) | Increasing material manufacturing method and device based on feature | |
CN103611646B (en) | Spray robot space path planing method | |
CN101769754A (en) | Quasi three-dimensional map-based mobile robot global path planning method | |
CN106054882A (en) | Robot obstacle avoidance method | |
JP2017502430A (en) | Topology optimization method of structure preservation of lattice structure for additive manufacturing | |
JP2008186440A (en) | Topology evolution optimization computing method for structural design | |
CN110398967A (en) | A kind of multirobot collaboration trace information processing method using discretization method | |
Janchiv et al. | Complete coverage path planning for multi-robots based on | |
CN106931970A (en) | Robot security's contexture by self air navigation aid in a kind of dynamic environment | |
CN102831648A (en) | Fixed node discharge curve approximation and curved surface mesh generation optimizing technology | |
CN104548597A (en) | Navigation grid automatic-generation method and navigation grid automatic-generation device | |
CN102129517A (en) | High-precision two-phase fluid interface capturing method | |
CN112621754A (en) | Design method for multi-robot-cooperated assembly line safety layout | |
CN104407613B (en) | Obstacle avoidance path smooth optimization method | |
CN115454062A (en) | Robot dynamic path planning method and system based on Betz curve | |
CN108391275A (en) | Sensor network cover control method based on three-dimensional directional sensing model and system | |
Zhang et al. | A coalition-structure’s generation method for solving cooperative computing problems in edge computing environments | |
CN105760549B (en) | Nearest Neighbor based on attribute graph model | |
CN101572014A (en) | Graphic fusion based method for automatically generating air traffic control airway route RVSM region | |
CN116009552A (en) | Path planning method, device, equipment and storage medium | |
CN101794534B (en) | Automatic typesetting and positioning method of statistical graphs in zonal statistical map | |
CN114995384A (en) | Improved HHO algorithm AGV path planning fusing neural network | |
CN110399679B (en) | Foundation position optimization design method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |