This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Research ArticleA Seed-Based Plant Propagation AlgorithmThe Feeding Station Model
Muhammad Sulaiman12 and Abdellah Salhi1
1Department of Mathematical Sciences University of Essex Colchester CO4 3SQ UK2Department of Mathematics Abdul Wali Khan University Mardan Khyber Pakhtunkhwa Pakistan
Correspondence should be addressed to Muhammad Sulaiman sulaiman513yahoocouk
Received 25 December 2014 Revised 10 February 2015 Accepted 10 February 2015
Academic Editor Xinyu Li
Copyright copy 2015 M Sulaiman and A Salhi This is an open access article distributed under the Creative Commons AttributionLicense which permits unrestricted use distribution and reproduction in any medium provided the original work is properlycited
The seasonal production of fruit and seeds is akin to opening a feeding station such as a restaurant Agents coming to feed onthe fruit are like customers attending the restaurant they arrive at a certain rate and get served at a certain rate following someappropriate processes The same applies to birds and animals visiting and feeding on ripe fruit produced by plants such as thestrawberry plantThis phenomenon underpins the seed dispersion of the plants Modelling it as a queuing process results in a seed-based searchoptimisation algorithm This variant of the Plant Propagation Algorithm is described analysed tested on nontrivialproblems and compared with well established algorithms The results are included
1 Introduction
Plants have evolved a variety of ways to propagate Propaga-tion with seeds is perhaps the most common of them all andone which takes advantage of all sorts of agents ranging fromwind to water birds and animals In [1] a Plant PropagationAlgorithm based on the way the strawberry plant propagatesusing runners has been introduced Here we consider thecase where the strawberry plant uses seeds to propagate
Plants rely heavily on the dispersion of their seeds tocolonise new territories and to improve their survival [2 3]There are a lot of studies and models of seed dispersionparticularly for trees [2ndash6] Dispersion by wind and ballisticmeans is probably the most studied of all approaches [7ndash9] However in the case of the strawberry plant given theway the seeds stick to the surface of the fruit (Figure 1(a))[10] dispersion by wind or mechanical means is very limitedAnimals however and birds in particular are the ideal agentsfor dispersion [2 3 11 12] in this case
There are many biologically inspired optimization algo-rithms in the literature [13 14] The Flower PollinationAlgorithm (FPA) is inspired by the pollination of flowersthrough different agents [8] the swarm data clustering algo-rithm is inspired by pollination by bees [15] Particle SwarmOptimization (PSO) is inspired by the foraging behavior of
groups of animals and insects [16 17] the Artificial BeeColony (ABC) simulates the foraging behavior of honey bees[18 19] the Firefly algorithm is inspired by the flashingfireflies when trying to attract a mate [20 21] the SocialSpider Optimization (SSO) algorithm is inspired by thecooperative behavior of social spiders [22] The list couldeasily be extended
The Plant Propagation Algorithm (PPA) also known asthe strawberry algorithm was inspired by the way plants andspecifically the strawberry plants propagate using runners[1 23] The attraction of PPA is that it can be implementedeasily for all sorts of optimization problems Moreover it hasfew algorithm specific arbitrary parameters It follows theprinciple that plants in good spots with plenty of nutrientswill send many short runners They send few long runnerswhen in nutrient poor spots With long runners PPA triesto explore the search space while short runners enable it toexploit the solution space well In this paper we investigatean alternative PPAwhich is entirely based on the propagationby seeds of the strawberry plant Because of the periodicnature of fruit and seed production it amounts to setting up afeeding station for the attention of potential seed-dispersingagents [24] Hence the feeding station model used here andthe resulting Seed-Based Plant Propagation Algorithm orSbPPA
Hindawi Publishing Corporatione Scientific World JournalVolume 2015 Article ID 904364 16 pageshttpdxdoiorg1011552015904364
2 The Scientific World Journal
(a) Strawberry fruit with seeds (b) Strawberry flower (c) A strawberry eaten by bird(s)
(d) A bird eating strawberries (e) Strawberry plants showing runners and seededfruit
Figure 1 Strawberry plant propagation through seed dispersion [25ndash28]
SbPPA is tested on both unconstrained and constrainedbenchmark problems also used in [22 29 30] Experimentalresults are presented in Tables 4ndash7 in terms of best meanworst and standard deviation for all algorithms The paperis organised as follows In Section 2 we briefly introducethe feeding station model representing strawberry plants infruit and the main characteristics of the paths followed bydifferent agents that disperse the seeds Section 3 presentsthe SbPPA in pseudocode form The experimental settingsresults and convergence graphs for different problems aregiven in Section 4
2 Aspects of the Feeding Station Model
Some animals and plants depend on each other to conservetheir species [31]Thusmany plants require for effective seed
dispersal the visits of frugivorous birds or animals accordingto a certain distribution [2 3 32 33]
Seed dispersal by different agents is also called ldquoseedshadowrdquo [32] this shows the abundance of seeds spreadglobally or locally around parent plants Here a queuingmodel is used which in the context of a strawberry feedingstation model involves two parts
(1) the quantity of fruit or seeds available to agents whichimplies the rate at which the agents will visit theplants
(2) a probability density function that tells us about theservice rate with which the agents are served by theplants
The model estimates the quantity of seeds that is spreadlocally compared to that dispersed globally [34ndash38] Thereare two aspects that need to be balanced exploitation which
The Scientific World Journal 3
is represented by the dispersal of seeds around the plantsand exploration which ensures that the search space is wellcovered
Agents arrive at plants in a random process Assumethat at most one agent arrives to the plants in any unit oftime (orderliness condition) It is further supposed that theprobability of arrivals of agents to the plants remains thesame for a particular period of timeThis period correspondsto when the plants are in fruit and during which timethe number of visitors is stable (stationarity condition)Furthermore it is assumed that the arrival of one agent doesnot affect the rest of arrivals (independence)
With these assumptions in mind the arrival of agentsto plants follows a Poisson process [39 40] which can beformally described as follows Let 1198831015840 be the random variablerepresenting the number of arrivals per unit of time 119905 Thenthe probability of 119896 arrivals over 119905 is
119875 (1198831015840= 119896) =
(120582119905)119896119890minus120582119905
119896 (1)
where 120582 denotes the mean arrival rate of agents per time unit119905 On the other hand the time taken by agents in successfullyeating fruit and leaving to disperse its seeds in other wordsthe service time for agents is expressed by a random variablewhich follows the exponential probability distribution [41]This can be expressed as follows
where 120583 is the average number of agents that can feed at time119905 Let us assume that the arrival rate of agents is less than thefruits available on all plants per unit of time therefore 120582 lt 120583
We assume that the system is in steady state Let119860 denotethe average number of agents in the strawberry field (somealready eating and the rest waiting to feed) and119860
119902the average
number of agents waiting to get the chance to feed If wedenote the average number of agents eating fruits by120582120583 thenby Littlersquos formula [42] we have
119860 = 119860119902+120582
120583 (3)
Since the plant needs to maximise dispersion this isequivalent to having a large 119860
119902in (3) Therefore from this
equation we need to solve the following problem
Maximize 119860119902= 119860 minus
120582
120583
subject to 1198921(120582 120583) = 120582 lt 120583 + 1
120582 gt 0 120583 gt 0
(4)
where 119860 = 10 which represents the population size in theimplementationThe simple limits on the variables are 0 lt 120582120583 le 100 The optimum solution to this particular problem is120582 = 11 120583 = 01 and 119860
119902= 1
Frugivores may travel far away from the plants and hencewill disperse the seeds far and wide This feeding behaviourtypically follows a Levy distribution [43ndash45] In the followingwe present some basic facts about it
21 Levy Distribution The Levy distribution is a probabilitydensity distribution for random variables Here the randomvariables represent the directions of flights of arbitrary birdsThis function ranges over real numbers in the domainrepresented by the problem search space
The flight lengths of the agents served by the plants followa heavy tailed power law distribution [14] represented by
119871 (119904) sim |119904|minus1minus120573
(5)
where 119871(119904) denotes the Levy distribution with index 120573 isin
(0 2) Levy flights are unique arbitrary excursions whose steplengths are drawn from (5) An alternative form of Levydistribution is [14]
In terms of the Fourier transform [14] the limiting value of119871(119904) can be written as
lim119904rarrinfin
119871 (119904) =120572120573Γ (120573) sin (1205871205732)
120587 |119904|1+120573
(8)
where Γ(120573) is the Gamma function [46] defined by
Γ (120573) = int
infin
0
119909120573minus1
119890minus119909119889119909 (9)
The steps 119871(119904) are generated by Mantegnarsquos algorithm [14]This algorithm ensures that the behaviour of Levy flights issymmetric and stable as shown in Figure 3(b)
3 Strawberry Plant Propagation AlgorithmThe Feeding Station Model
We assume that the arrival of different agents (birds andanimals) to the plants to feed is according to the Poissondistribution [40] As per the solution of problem (4) themean arrival rate is 120582 = 11 and NP = 10 is the size ofthe agents population Let 119896 = 1 2 119860 be the possiblenumbers of agents visiting the plants per unit time Withthese assumptions the graphic representation of (1) results inFigure 2
As already stated it is essential in this algorithm tobalance exploration and exploitation To this end we choosea threshold value of the Poisson probability that dictateshow much exploration and exploitation are done duringthe search The probability Poiss(120582) lt 005 means that
4 The Scientific World JournalPr
obab
ilitie
s Po
iss(120582)
Arrival rate of agents per unit time
025
02
015
01
005
00 2 4 6 8 10
Figure 2 Distribution of agents arriving at strawberry plants to eatfruit and disperse seeds
exploitation is covered In this case (10) below is used whichhelps the algorithm to search locally
119909lowast
119894119895=
119909119894119895+ 120585119895(119909119894119895minus 119909119897119895) if PR le 08 119895 = 1 2 119899
119894 119897 = 1 2 NP 119894 = 119897
119909119894119895 Otherwise
(10)
where PR denotes the rate of dispersion of the seeds locallyaround SP 119909lowast
119894119895and 119909
119894119895isin [119886119895119887119895] are the 119895th coordinates of
the seeds 119883lowast119894and119883
119894 respectively 119886
119895and 119887119895are the 119895th lower
and upper bounds defining the search space of the problemand 120585119895isin [minus1 1] The indices 119897 and 119894 are mutually exclusive
On the other hand if Poiss(120582) ge 005 then global disper-sion of seeds becomes more prominent This is implementedby using the following equation
119909lowast
119894119895=
119909119894119895+ 119871119894(119909119894119895minus 120579119895) if PR le 08 120579
119895isin [119886119895119887119895]
119894 = 1 2 NP
119895 = 1 2 119899
119909119894119895 Otherwise
(11)
where 119871119894is a step drawn from the Levy distribution [14] and
120579119895is a random coordinate within the search space Equations
(10) and (11) perturb the current solution the results of whichcan be seen in Figures 3(a) and 3(b) respectively
As mentioned in Algorithm 1 we first collect the bestsolutions from the first NP trial runs to form a populationof potentially good solutions denoted by popbest The conver-gence rate of SbPPA is shown in Figures 4 and 5 for differenttest problems used in our experiments (see Appendices) Thestatistics values best worst mean and standard deviation arecalculated based on popbest
The seed-based propagation process of SP can be repre-sented in the following steps
Sear
ch sp
ace
Spring Design Problem
times105
15
10
5
0
minus5
minus10
minus150 2 4 6
(a) Perturbations by (10)Se
arch
spac
eSpring Design Problem
times105
40
20
0
minus20
minus40
minus60
minus800 5 10 15
(b) Perturbations by (11)
Figure 3Overall performance of SbPPAonSpringDesignProblem
(1) The dispersal of seeds in the neighbourhood of theSP as shown in Figure 1(e) is carried out either byfruits fallen from strawberry plants after they becomeripe or by agents The step lengths for this phase arecalculated using (10)
(2) Seeds are spread globally through agents as shownin Figures 1(c) and 1(d) The step lengths for thesetravelling agents are drawn from the Levy distribution[14]
(3) The probabilities Poiss(120582) that a certain number 119896 ofagents will arrive to SP to eat fruits and disperse it
The Scientific World Journal 5
(1) NPlarr Population size 119903 larr Counter of trial runs MaxExp larrMaximum experiments(2) for 119903 = 1 MaxExp do(3) if 119903 le NP then(4) Create a random population of seeds pop = 119883
119894| 119894 = 1 2 NP
using (12) and collect the best solutions from each trial run in popbest(5) Evaluate the population pop(6) end if(7) while 119903 gt NP do(8) Use updated population popbest(9) end while(10) while (the stopping criteria is not satisfied) do(11) for 119894 = 1 to NP do(12) if Poiss(120582)
119894ge 005 then ⊳ (Global or local seed dispersion)
(13) for 119895 = 1 to 119899 do ⊳ (119899 is number of dimensions)(14) if rand le PR then ⊳ (PR = Perturbation Rate)(15) Update the current entry according to (11)(16) end if(17) end for(18) else(19) for 119895 = 1 to 119899 do(20) if rand le PR then(21) Update the current entry according to (10)(22) end if(23) end for(24) end if(25) end for(26) Update current best(27) end while(28) Return Updated population and global best solution(29) end for
isin [119886119895119887119895] is the 119895th coordinate of solution 119883
119894 119886119895
and 119887119895are the 119895th coordinates of the bounds describing the
search space of the problem and 120578119895isin (0 1) This means that
119883119894= [119909119894119895] for 119895 = 1 119899 represents the position of the 119895th
seed in population pop
4 Experimental Settings and Discussion
In our experiments we tested SbPPA against some recentlydeveloped algorithms and somewell established and standardones Our set of test problems includes benchmark con-strained and unconstrained optimization problems [22 3048 49] The results are compared in terms of statistics (bestworst mean and standard deviation) for solutions obtained
by SbPPA ABC [18 50] PSO [51] FF [21] HPA [29] SSO-C [22] Classical Evolutionary Programming (CEP) [30] andFast Evolutionary Programming (FEP) [30] The detaileddescriptions of these problems are given in Appendices
In Tables 4 and 7 the significance of results is shown interms of wintieloss (see Table 2 in [52]) according to thefollowing notations
(i) (+) when SbPPA is better(ii) (asymp) when the results are approximately the same as
those obtained with SbPPA(iii) (minus) when SbPPA is worse
Moreover in Tables 5 and 6 the significance of resultsobtained with SbPPA is highlighted
41 Parameter Settings The parameter settings are given inTables 1ndash3
5 Conclusion
In this paper a new metaheuristic referred to as the Seed-Based Plant Propagation Algorithm (SbPPA) [47] has beenproposed Plants have evolved a variety of ways to propagatePropagation through seeds is perhaps the most commonof them all and one which takes advantage of all sorts of
agents ranging from wind to water birds and animals Thestrawberry plant uses both runners and seeds to propagateHere we consider the propagation through seeds that thestrawberry plant has evolved to design an efficient optimiza-tion algorithm
To capture the dispersal process we adopt a queuingapproach which given the extent of fruit produced indicatesthe extent of seeds dispersed and hence the effectivenessof the searchoptimization algorithm based on this processLooking at the random process of agents using the plants(feeding station) it is reasonable to assume that it is of thePoisson type On the other hand the time taken by agentsin successfully eating fruit and leaving to disperse its seedsin other words the service time for agents is expressed by arandom variable which follows the exponential probabilitydistribution To this end we choose a threshold value ofthe Poisson probability that dictates how much explorationand exploitation are done during the search An alternativestrategy has been adopted here This strategy consists inmaking sure that the initial population is as good as theuser can afford it to be by using best solutions found sofar The effects of this strategy on convergence are shownthrough convergence plots of Figures 4 and 5 for some of thesolved problems SbPPA is easy to implement as it requiresless arbitrary parameter settings than other algorithms Thesuccess rate of SbPPA increases as it gets its population of bestsolutions It has been implemented for both unconstrainedand constrained optimization problems Its performancecompared to that of other algorithms points to SbPPA asbeing superior
Appendices
A Unconstrained GlobalOptimization Problems
See Tables 8 and 9
B Set of Constrained Global OptimizationProblems Used in Our Experiments
B1 CP1 Consider the following
Min 119891 (119909) = 5
4
sum
119889=1
119909119889minus 5
4
sum
119889=1
1199092
119889minus
13
sum
119889=5
119909119889
subject to 1198921(119909) = 2119909
1+ 21199092+ 11990910+ 11990911minus 10 le 0
1198922(119909) = 2119909
1+ 21199093+ 11990910+ 11990912minus 10 le 0
1198923(119909) = 2119909
2+ 21199093+ 11990911+ 11990912minus 10 le 0
1198924(119909) = minus8119909
1+ 11990910le 0
1198925(119909) = minus8119909
2+ 11990911le 0
1198926(119909) = minus8119909
3+ 11990912le 0
1198927(119909) = minus2119909
4minus 1199095+ 11990910le 0
1198928(119909) = minus2119909
6minus 1199097+ 11990911le 0
1198929(119909) = minus2119909
8minus 1199099+ 11990912le 0
(B1)
where bounds are 0 le 119909119894le 1 (119894 = 1 9 13) 0 le
119909119894le 100 (119894 = 10 11 12) The global optimum is at 119909lowast =
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
(a) Strawberry fruit with seeds (b) Strawberry flower (c) A strawberry eaten by bird(s)
(d) A bird eating strawberries (e) Strawberry plants showing runners and seededfruit
Figure 1 Strawberry plant propagation through seed dispersion [25ndash28]
SbPPA is tested on both unconstrained and constrainedbenchmark problems also used in [22 29 30] Experimentalresults are presented in Tables 4ndash7 in terms of best meanworst and standard deviation for all algorithms The paperis organised as follows In Section 2 we briefly introducethe feeding station model representing strawberry plants infruit and the main characteristics of the paths followed bydifferent agents that disperse the seeds Section 3 presentsthe SbPPA in pseudocode form The experimental settingsresults and convergence graphs for different problems aregiven in Section 4
2 Aspects of the Feeding Station Model
Some animals and plants depend on each other to conservetheir species [31]Thusmany plants require for effective seed
dispersal the visits of frugivorous birds or animals accordingto a certain distribution [2 3 32 33]
Seed dispersal by different agents is also called ldquoseedshadowrdquo [32] this shows the abundance of seeds spreadglobally or locally around parent plants Here a queuingmodel is used which in the context of a strawberry feedingstation model involves two parts
(1) the quantity of fruit or seeds available to agents whichimplies the rate at which the agents will visit theplants
(2) a probability density function that tells us about theservice rate with which the agents are served by theplants
The model estimates the quantity of seeds that is spreadlocally compared to that dispersed globally [34ndash38] Thereare two aspects that need to be balanced exploitation which
The Scientific World Journal 3
is represented by the dispersal of seeds around the plantsand exploration which ensures that the search space is wellcovered
Agents arrive at plants in a random process Assumethat at most one agent arrives to the plants in any unit oftime (orderliness condition) It is further supposed that theprobability of arrivals of agents to the plants remains thesame for a particular period of timeThis period correspondsto when the plants are in fruit and during which timethe number of visitors is stable (stationarity condition)Furthermore it is assumed that the arrival of one agent doesnot affect the rest of arrivals (independence)
With these assumptions in mind the arrival of agentsto plants follows a Poisson process [39 40] which can beformally described as follows Let 1198831015840 be the random variablerepresenting the number of arrivals per unit of time 119905 Thenthe probability of 119896 arrivals over 119905 is
119875 (1198831015840= 119896) =
(120582119905)119896119890minus120582119905
119896 (1)
where 120582 denotes the mean arrival rate of agents per time unit119905 On the other hand the time taken by agents in successfullyeating fruit and leaving to disperse its seeds in other wordsthe service time for agents is expressed by a random variablewhich follows the exponential probability distribution [41]This can be expressed as follows
where 120583 is the average number of agents that can feed at time119905 Let us assume that the arrival rate of agents is less than thefruits available on all plants per unit of time therefore 120582 lt 120583
We assume that the system is in steady state Let119860 denotethe average number of agents in the strawberry field (somealready eating and the rest waiting to feed) and119860
119902the average
number of agents waiting to get the chance to feed If wedenote the average number of agents eating fruits by120582120583 thenby Littlersquos formula [42] we have
119860 = 119860119902+120582
120583 (3)
Since the plant needs to maximise dispersion this isequivalent to having a large 119860
119902in (3) Therefore from this
equation we need to solve the following problem
Maximize 119860119902= 119860 minus
120582
120583
subject to 1198921(120582 120583) = 120582 lt 120583 + 1
120582 gt 0 120583 gt 0
(4)
where 119860 = 10 which represents the population size in theimplementationThe simple limits on the variables are 0 lt 120582120583 le 100 The optimum solution to this particular problem is120582 = 11 120583 = 01 and 119860
119902= 1
Frugivores may travel far away from the plants and hencewill disperse the seeds far and wide This feeding behaviourtypically follows a Levy distribution [43ndash45] In the followingwe present some basic facts about it
21 Levy Distribution The Levy distribution is a probabilitydensity distribution for random variables Here the randomvariables represent the directions of flights of arbitrary birdsThis function ranges over real numbers in the domainrepresented by the problem search space
The flight lengths of the agents served by the plants followa heavy tailed power law distribution [14] represented by
119871 (119904) sim |119904|minus1minus120573
(5)
where 119871(119904) denotes the Levy distribution with index 120573 isin
(0 2) Levy flights are unique arbitrary excursions whose steplengths are drawn from (5) An alternative form of Levydistribution is [14]
In terms of the Fourier transform [14] the limiting value of119871(119904) can be written as
lim119904rarrinfin
119871 (119904) =120572120573Γ (120573) sin (1205871205732)
120587 |119904|1+120573
(8)
where Γ(120573) is the Gamma function [46] defined by
Γ (120573) = int
infin
0
119909120573minus1
119890minus119909119889119909 (9)
The steps 119871(119904) are generated by Mantegnarsquos algorithm [14]This algorithm ensures that the behaviour of Levy flights issymmetric and stable as shown in Figure 3(b)
3 Strawberry Plant Propagation AlgorithmThe Feeding Station Model
We assume that the arrival of different agents (birds andanimals) to the plants to feed is according to the Poissondistribution [40] As per the solution of problem (4) themean arrival rate is 120582 = 11 and NP = 10 is the size ofthe agents population Let 119896 = 1 2 119860 be the possiblenumbers of agents visiting the plants per unit time Withthese assumptions the graphic representation of (1) results inFigure 2
As already stated it is essential in this algorithm tobalance exploration and exploitation To this end we choosea threshold value of the Poisson probability that dictateshow much exploration and exploitation are done duringthe search The probability Poiss(120582) lt 005 means that
4 The Scientific World JournalPr
obab
ilitie
s Po
iss(120582)
Arrival rate of agents per unit time
025
02
015
01
005
00 2 4 6 8 10
Figure 2 Distribution of agents arriving at strawberry plants to eatfruit and disperse seeds
exploitation is covered In this case (10) below is used whichhelps the algorithm to search locally
119909lowast
119894119895=
119909119894119895+ 120585119895(119909119894119895minus 119909119897119895) if PR le 08 119895 = 1 2 119899
119894 119897 = 1 2 NP 119894 = 119897
119909119894119895 Otherwise
(10)
where PR denotes the rate of dispersion of the seeds locallyaround SP 119909lowast
119894119895and 119909
119894119895isin [119886119895119887119895] are the 119895th coordinates of
the seeds 119883lowast119894and119883
119894 respectively 119886
119895and 119887119895are the 119895th lower
and upper bounds defining the search space of the problemand 120585119895isin [minus1 1] The indices 119897 and 119894 are mutually exclusive
On the other hand if Poiss(120582) ge 005 then global disper-sion of seeds becomes more prominent This is implementedby using the following equation
119909lowast
119894119895=
119909119894119895+ 119871119894(119909119894119895minus 120579119895) if PR le 08 120579
119895isin [119886119895119887119895]
119894 = 1 2 NP
119895 = 1 2 119899
119909119894119895 Otherwise
(11)
where 119871119894is a step drawn from the Levy distribution [14] and
120579119895is a random coordinate within the search space Equations
(10) and (11) perturb the current solution the results of whichcan be seen in Figures 3(a) and 3(b) respectively
As mentioned in Algorithm 1 we first collect the bestsolutions from the first NP trial runs to form a populationof potentially good solutions denoted by popbest The conver-gence rate of SbPPA is shown in Figures 4 and 5 for differenttest problems used in our experiments (see Appendices) Thestatistics values best worst mean and standard deviation arecalculated based on popbest
The seed-based propagation process of SP can be repre-sented in the following steps
Sear
ch sp
ace
Spring Design Problem
times105
15
10
5
0
minus5
minus10
minus150 2 4 6
(a) Perturbations by (10)Se
arch
spac
eSpring Design Problem
times105
40
20
0
minus20
minus40
minus60
minus800 5 10 15
(b) Perturbations by (11)
Figure 3Overall performance of SbPPAonSpringDesignProblem
(1) The dispersal of seeds in the neighbourhood of theSP as shown in Figure 1(e) is carried out either byfruits fallen from strawberry plants after they becomeripe or by agents The step lengths for this phase arecalculated using (10)
(2) Seeds are spread globally through agents as shownin Figures 1(c) and 1(d) The step lengths for thesetravelling agents are drawn from the Levy distribution[14]
(3) The probabilities Poiss(120582) that a certain number 119896 ofagents will arrive to SP to eat fruits and disperse it
The Scientific World Journal 5
(1) NPlarr Population size 119903 larr Counter of trial runs MaxExp larrMaximum experiments(2) for 119903 = 1 MaxExp do(3) if 119903 le NP then(4) Create a random population of seeds pop = 119883
119894| 119894 = 1 2 NP
using (12) and collect the best solutions from each trial run in popbest(5) Evaluate the population pop(6) end if(7) while 119903 gt NP do(8) Use updated population popbest(9) end while(10) while (the stopping criteria is not satisfied) do(11) for 119894 = 1 to NP do(12) if Poiss(120582)
119894ge 005 then ⊳ (Global or local seed dispersion)
(13) for 119895 = 1 to 119899 do ⊳ (119899 is number of dimensions)(14) if rand le PR then ⊳ (PR = Perturbation Rate)(15) Update the current entry according to (11)(16) end if(17) end for(18) else(19) for 119895 = 1 to 119899 do(20) if rand le PR then(21) Update the current entry according to (10)(22) end if(23) end for(24) end if(25) end for(26) Update current best(27) end while(28) Return Updated population and global best solution(29) end for
isin [119886119895119887119895] is the 119895th coordinate of solution 119883
119894 119886119895
and 119887119895are the 119895th coordinates of the bounds describing the
search space of the problem and 120578119895isin (0 1) This means that
119883119894= [119909119894119895] for 119895 = 1 119899 represents the position of the 119895th
seed in population pop
4 Experimental Settings and Discussion
In our experiments we tested SbPPA against some recentlydeveloped algorithms and somewell established and standardones Our set of test problems includes benchmark con-strained and unconstrained optimization problems [22 3048 49] The results are compared in terms of statistics (bestworst mean and standard deviation) for solutions obtained
by SbPPA ABC [18 50] PSO [51] FF [21] HPA [29] SSO-C [22] Classical Evolutionary Programming (CEP) [30] andFast Evolutionary Programming (FEP) [30] The detaileddescriptions of these problems are given in Appendices
In Tables 4 and 7 the significance of results is shown interms of wintieloss (see Table 2 in [52]) according to thefollowing notations
(i) (+) when SbPPA is better(ii) (asymp) when the results are approximately the same as
those obtained with SbPPA(iii) (minus) when SbPPA is worse
Moreover in Tables 5 and 6 the significance of resultsobtained with SbPPA is highlighted
41 Parameter Settings The parameter settings are given inTables 1ndash3
5 Conclusion
In this paper a new metaheuristic referred to as the Seed-Based Plant Propagation Algorithm (SbPPA) [47] has beenproposed Plants have evolved a variety of ways to propagatePropagation through seeds is perhaps the most commonof them all and one which takes advantage of all sorts of
agents ranging from wind to water birds and animals Thestrawberry plant uses both runners and seeds to propagateHere we consider the propagation through seeds that thestrawberry plant has evolved to design an efficient optimiza-tion algorithm
To capture the dispersal process we adopt a queuingapproach which given the extent of fruit produced indicatesthe extent of seeds dispersed and hence the effectivenessof the searchoptimization algorithm based on this processLooking at the random process of agents using the plants(feeding station) it is reasonable to assume that it is of thePoisson type On the other hand the time taken by agentsin successfully eating fruit and leaving to disperse its seedsin other words the service time for agents is expressed by arandom variable which follows the exponential probabilitydistribution To this end we choose a threshold value ofthe Poisson probability that dictates how much explorationand exploitation are done during the search An alternativestrategy has been adopted here This strategy consists inmaking sure that the initial population is as good as theuser can afford it to be by using best solutions found sofar The effects of this strategy on convergence are shownthrough convergence plots of Figures 4 and 5 for some of thesolved problems SbPPA is easy to implement as it requiresless arbitrary parameter settings than other algorithms Thesuccess rate of SbPPA increases as it gets its population of bestsolutions It has been implemented for both unconstrainedand constrained optimization problems Its performancecompared to that of other algorithms points to SbPPA asbeing superior
Appendices
A Unconstrained GlobalOptimization Problems
See Tables 8 and 9
B Set of Constrained Global OptimizationProblems Used in Our Experiments
B1 CP1 Consider the following
Min 119891 (119909) = 5
4
sum
119889=1
119909119889minus 5
4
sum
119889=1
1199092
119889minus
13
sum
119889=5
119909119889
subject to 1198921(119909) = 2119909
1+ 21199092+ 11990910+ 11990911minus 10 le 0
1198922(119909) = 2119909
1+ 21199093+ 11990910+ 11990912minus 10 le 0
1198923(119909) = 2119909
2+ 21199093+ 11990911+ 11990912minus 10 le 0
1198924(119909) = minus8119909
1+ 11990910le 0
1198925(119909) = minus8119909
2+ 11990911le 0
1198926(119909) = minus8119909
3+ 11990912le 0
1198927(119909) = minus2119909
4minus 1199095+ 11990910le 0
1198928(119909) = minus2119909
6minus 1199097+ 11990911le 0
1198929(119909) = minus2119909
8minus 1199099+ 11990912le 0
(B1)
where bounds are 0 le 119909119894le 1 (119894 = 1 9 13) 0 le
119909119894le 100 (119894 = 10 11 12) The global optimum is at 119909lowast =
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
is represented by the dispersal of seeds around the plantsand exploration which ensures that the search space is wellcovered
Agents arrive at plants in a random process Assumethat at most one agent arrives to the plants in any unit oftime (orderliness condition) It is further supposed that theprobability of arrivals of agents to the plants remains thesame for a particular period of timeThis period correspondsto when the plants are in fruit and during which timethe number of visitors is stable (stationarity condition)Furthermore it is assumed that the arrival of one agent doesnot affect the rest of arrivals (independence)
With these assumptions in mind the arrival of agentsto plants follows a Poisson process [39 40] which can beformally described as follows Let 1198831015840 be the random variablerepresenting the number of arrivals per unit of time 119905 Thenthe probability of 119896 arrivals over 119905 is
119875 (1198831015840= 119896) =
(120582119905)119896119890minus120582119905
119896 (1)
where 120582 denotes the mean arrival rate of agents per time unit119905 On the other hand the time taken by agents in successfullyeating fruit and leaving to disperse its seeds in other wordsthe service time for agents is expressed by a random variablewhich follows the exponential probability distribution [41]This can be expressed as follows
where 120583 is the average number of agents that can feed at time119905 Let us assume that the arrival rate of agents is less than thefruits available on all plants per unit of time therefore 120582 lt 120583
We assume that the system is in steady state Let119860 denotethe average number of agents in the strawberry field (somealready eating and the rest waiting to feed) and119860
119902the average
number of agents waiting to get the chance to feed If wedenote the average number of agents eating fruits by120582120583 thenby Littlersquos formula [42] we have
119860 = 119860119902+120582
120583 (3)
Since the plant needs to maximise dispersion this isequivalent to having a large 119860
119902in (3) Therefore from this
equation we need to solve the following problem
Maximize 119860119902= 119860 minus
120582
120583
subject to 1198921(120582 120583) = 120582 lt 120583 + 1
120582 gt 0 120583 gt 0
(4)
where 119860 = 10 which represents the population size in theimplementationThe simple limits on the variables are 0 lt 120582120583 le 100 The optimum solution to this particular problem is120582 = 11 120583 = 01 and 119860
119902= 1
Frugivores may travel far away from the plants and hencewill disperse the seeds far and wide This feeding behaviourtypically follows a Levy distribution [43ndash45] In the followingwe present some basic facts about it
21 Levy Distribution The Levy distribution is a probabilitydensity distribution for random variables Here the randomvariables represent the directions of flights of arbitrary birdsThis function ranges over real numbers in the domainrepresented by the problem search space
The flight lengths of the agents served by the plants followa heavy tailed power law distribution [14] represented by
119871 (119904) sim |119904|minus1minus120573
(5)
where 119871(119904) denotes the Levy distribution with index 120573 isin
(0 2) Levy flights are unique arbitrary excursions whose steplengths are drawn from (5) An alternative form of Levydistribution is [14]
In terms of the Fourier transform [14] the limiting value of119871(119904) can be written as
lim119904rarrinfin
119871 (119904) =120572120573Γ (120573) sin (1205871205732)
120587 |119904|1+120573
(8)
where Γ(120573) is the Gamma function [46] defined by
Γ (120573) = int
infin
0
119909120573minus1
119890minus119909119889119909 (9)
The steps 119871(119904) are generated by Mantegnarsquos algorithm [14]This algorithm ensures that the behaviour of Levy flights issymmetric and stable as shown in Figure 3(b)
3 Strawberry Plant Propagation AlgorithmThe Feeding Station Model
We assume that the arrival of different agents (birds andanimals) to the plants to feed is according to the Poissondistribution [40] As per the solution of problem (4) themean arrival rate is 120582 = 11 and NP = 10 is the size ofthe agents population Let 119896 = 1 2 119860 be the possiblenumbers of agents visiting the plants per unit time Withthese assumptions the graphic representation of (1) results inFigure 2
As already stated it is essential in this algorithm tobalance exploration and exploitation To this end we choosea threshold value of the Poisson probability that dictateshow much exploration and exploitation are done duringthe search The probability Poiss(120582) lt 005 means that
4 The Scientific World JournalPr
obab
ilitie
s Po
iss(120582)
Arrival rate of agents per unit time
025
02
015
01
005
00 2 4 6 8 10
Figure 2 Distribution of agents arriving at strawberry plants to eatfruit and disperse seeds
exploitation is covered In this case (10) below is used whichhelps the algorithm to search locally
119909lowast
119894119895=
119909119894119895+ 120585119895(119909119894119895minus 119909119897119895) if PR le 08 119895 = 1 2 119899
119894 119897 = 1 2 NP 119894 = 119897
119909119894119895 Otherwise
(10)
where PR denotes the rate of dispersion of the seeds locallyaround SP 119909lowast
119894119895and 119909
119894119895isin [119886119895119887119895] are the 119895th coordinates of
the seeds 119883lowast119894and119883
119894 respectively 119886
119895and 119887119895are the 119895th lower
and upper bounds defining the search space of the problemand 120585119895isin [minus1 1] The indices 119897 and 119894 are mutually exclusive
On the other hand if Poiss(120582) ge 005 then global disper-sion of seeds becomes more prominent This is implementedby using the following equation
119909lowast
119894119895=
119909119894119895+ 119871119894(119909119894119895minus 120579119895) if PR le 08 120579
119895isin [119886119895119887119895]
119894 = 1 2 NP
119895 = 1 2 119899
119909119894119895 Otherwise
(11)
where 119871119894is a step drawn from the Levy distribution [14] and
120579119895is a random coordinate within the search space Equations
(10) and (11) perturb the current solution the results of whichcan be seen in Figures 3(a) and 3(b) respectively
As mentioned in Algorithm 1 we first collect the bestsolutions from the first NP trial runs to form a populationof potentially good solutions denoted by popbest The conver-gence rate of SbPPA is shown in Figures 4 and 5 for differenttest problems used in our experiments (see Appendices) Thestatistics values best worst mean and standard deviation arecalculated based on popbest
The seed-based propagation process of SP can be repre-sented in the following steps
Sear
ch sp
ace
Spring Design Problem
times105
15
10
5
0
minus5
minus10
minus150 2 4 6
(a) Perturbations by (10)Se
arch
spac
eSpring Design Problem
times105
40
20
0
minus20
minus40
minus60
minus800 5 10 15
(b) Perturbations by (11)
Figure 3Overall performance of SbPPAonSpringDesignProblem
(1) The dispersal of seeds in the neighbourhood of theSP as shown in Figure 1(e) is carried out either byfruits fallen from strawberry plants after they becomeripe or by agents The step lengths for this phase arecalculated using (10)
(2) Seeds are spread globally through agents as shownin Figures 1(c) and 1(d) The step lengths for thesetravelling agents are drawn from the Levy distribution[14]
(3) The probabilities Poiss(120582) that a certain number 119896 ofagents will arrive to SP to eat fruits and disperse it
The Scientific World Journal 5
(1) NPlarr Population size 119903 larr Counter of trial runs MaxExp larrMaximum experiments(2) for 119903 = 1 MaxExp do(3) if 119903 le NP then(4) Create a random population of seeds pop = 119883
119894| 119894 = 1 2 NP
using (12) and collect the best solutions from each trial run in popbest(5) Evaluate the population pop(6) end if(7) while 119903 gt NP do(8) Use updated population popbest(9) end while(10) while (the stopping criteria is not satisfied) do(11) for 119894 = 1 to NP do(12) if Poiss(120582)
119894ge 005 then ⊳ (Global or local seed dispersion)
(13) for 119895 = 1 to 119899 do ⊳ (119899 is number of dimensions)(14) if rand le PR then ⊳ (PR = Perturbation Rate)(15) Update the current entry according to (11)(16) end if(17) end for(18) else(19) for 119895 = 1 to 119899 do(20) if rand le PR then(21) Update the current entry according to (10)(22) end if(23) end for(24) end if(25) end for(26) Update current best(27) end while(28) Return Updated population and global best solution(29) end for
isin [119886119895119887119895] is the 119895th coordinate of solution 119883
119894 119886119895
and 119887119895are the 119895th coordinates of the bounds describing the
search space of the problem and 120578119895isin (0 1) This means that
119883119894= [119909119894119895] for 119895 = 1 119899 represents the position of the 119895th
seed in population pop
4 Experimental Settings and Discussion
In our experiments we tested SbPPA against some recentlydeveloped algorithms and somewell established and standardones Our set of test problems includes benchmark con-strained and unconstrained optimization problems [22 3048 49] The results are compared in terms of statistics (bestworst mean and standard deviation) for solutions obtained
by SbPPA ABC [18 50] PSO [51] FF [21] HPA [29] SSO-C [22] Classical Evolutionary Programming (CEP) [30] andFast Evolutionary Programming (FEP) [30] The detaileddescriptions of these problems are given in Appendices
In Tables 4 and 7 the significance of results is shown interms of wintieloss (see Table 2 in [52]) according to thefollowing notations
(i) (+) when SbPPA is better(ii) (asymp) when the results are approximately the same as
those obtained with SbPPA(iii) (minus) when SbPPA is worse
Moreover in Tables 5 and 6 the significance of resultsobtained with SbPPA is highlighted
41 Parameter Settings The parameter settings are given inTables 1ndash3
5 Conclusion
In this paper a new metaheuristic referred to as the Seed-Based Plant Propagation Algorithm (SbPPA) [47] has beenproposed Plants have evolved a variety of ways to propagatePropagation through seeds is perhaps the most commonof them all and one which takes advantage of all sorts of
agents ranging from wind to water birds and animals Thestrawberry plant uses both runners and seeds to propagateHere we consider the propagation through seeds that thestrawberry plant has evolved to design an efficient optimiza-tion algorithm
To capture the dispersal process we adopt a queuingapproach which given the extent of fruit produced indicatesthe extent of seeds dispersed and hence the effectivenessof the searchoptimization algorithm based on this processLooking at the random process of agents using the plants(feeding station) it is reasonable to assume that it is of thePoisson type On the other hand the time taken by agentsin successfully eating fruit and leaving to disperse its seedsin other words the service time for agents is expressed by arandom variable which follows the exponential probabilitydistribution To this end we choose a threshold value ofthe Poisson probability that dictates how much explorationand exploitation are done during the search An alternativestrategy has been adopted here This strategy consists inmaking sure that the initial population is as good as theuser can afford it to be by using best solutions found sofar The effects of this strategy on convergence are shownthrough convergence plots of Figures 4 and 5 for some of thesolved problems SbPPA is easy to implement as it requiresless arbitrary parameter settings than other algorithms Thesuccess rate of SbPPA increases as it gets its population of bestsolutions It has been implemented for both unconstrainedand constrained optimization problems Its performancecompared to that of other algorithms points to SbPPA asbeing superior
Appendices
A Unconstrained GlobalOptimization Problems
See Tables 8 and 9
B Set of Constrained Global OptimizationProblems Used in Our Experiments
B1 CP1 Consider the following
Min 119891 (119909) = 5
4
sum
119889=1
119909119889minus 5
4
sum
119889=1
1199092
119889minus
13
sum
119889=5
119909119889
subject to 1198921(119909) = 2119909
1+ 21199092+ 11990910+ 11990911minus 10 le 0
1198922(119909) = 2119909
1+ 21199093+ 11990910+ 11990912minus 10 le 0
1198923(119909) = 2119909
2+ 21199093+ 11990911+ 11990912minus 10 le 0
1198924(119909) = minus8119909
1+ 11990910le 0
1198925(119909) = minus8119909
2+ 11990911le 0
1198926(119909) = minus8119909
3+ 11990912le 0
1198927(119909) = minus2119909
4minus 1199095+ 11990910le 0
1198928(119909) = minus2119909
6minus 1199097+ 11990911le 0
1198929(119909) = minus2119909
8minus 1199099+ 11990912le 0
(B1)
where bounds are 0 le 119909119894le 1 (119894 = 1 9 13) 0 le
119909119894le 100 (119894 = 10 11 12) The global optimum is at 119909lowast =
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
Figure 2 Distribution of agents arriving at strawberry plants to eatfruit and disperse seeds
exploitation is covered In this case (10) below is used whichhelps the algorithm to search locally
119909lowast
119894119895=
119909119894119895+ 120585119895(119909119894119895minus 119909119897119895) if PR le 08 119895 = 1 2 119899
119894 119897 = 1 2 NP 119894 = 119897
119909119894119895 Otherwise
(10)
where PR denotes the rate of dispersion of the seeds locallyaround SP 119909lowast
119894119895and 119909
119894119895isin [119886119895119887119895] are the 119895th coordinates of
the seeds 119883lowast119894and119883
119894 respectively 119886
119895and 119887119895are the 119895th lower
and upper bounds defining the search space of the problemand 120585119895isin [minus1 1] The indices 119897 and 119894 are mutually exclusive
On the other hand if Poiss(120582) ge 005 then global disper-sion of seeds becomes more prominent This is implementedby using the following equation
119909lowast
119894119895=
119909119894119895+ 119871119894(119909119894119895minus 120579119895) if PR le 08 120579
119895isin [119886119895119887119895]
119894 = 1 2 NP
119895 = 1 2 119899
119909119894119895 Otherwise
(11)
where 119871119894is a step drawn from the Levy distribution [14] and
120579119895is a random coordinate within the search space Equations
(10) and (11) perturb the current solution the results of whichcan be seen in Figures 3(a) and 3(b) respectively
As mentioned in Algorithm 1 we first collect the bestsolutions from the first NP trial runs to form a populationof potentially good solutions denoted by popbest The conver-gence rate of SbPPA is shown in Figures 4 and 5 for differenttest problems used in our experiments (see Appendices) Thestatistics values best worst mean and standard deviation arecalculated based on popbest
The seed-based propagation process of SP can be repre-sented in the following steps
Sear
ch sp
ace
Spring Design Problem
times105
15
10
5
0
minus5
minus10
minus150 2 4 6
(a) Perturbations by (10)Se
arch
spac
eSpring Design Problem
times105
40
20
0
minus20
minus40
minus60
minus800 5 10 15
(b) Perturbations by (11)
Figure 3Overall performance of SbPPAonSpringDesignProblem
(1) The dispersal of seeds in the neighbourhood of theSP as shown in Figure 1(e) is carried out either byfruits fallen from strawberry plants after they becomeripe or by agents The step lengths for this phase arecalculated using (10)
(2) Seeds are spread globally through agents as shownin Figures 1(c) and 1(d) The step lengths for thesetravelling agents are drawn from the Levy distribution[14]
(3) The probabilities Poiss(120582) that a certain number 119896 ofagents will arrive to SP to eat fruits and disperse it
The Scientific World Journal 5
(1) NPlarr Population size 119903 larr Counter of trial runs MaxExp larrMaximum experiments(2) for 119903 = 1 MaxExp do(3) if 119903 le NP then(4) Create a random population of seeds pop = 119883
119894| 119894 = 1 2 NP
using (12) and collect the best solutions from each trial run in popbest(5) Evaluate the population pop(6) end if(7) while 119903 gt NP do(8) Use updated population popbest(9) end while(10) while (the stopping criteria is not satisfied) do(11) for 119894 = 1 to NP do(12) if Poiss(120582)
119894ge 005 then ⊳ (Global or local seed dispersion)
(13) for 119895 = 1 to 119899 do ⊳ (119899 is number of dimensions)(14) if rand le PR then ⊳ (PR = Perturbation Rate)(15) Update the current entry according to (11)(16) end if(17) end for(18) else(19) for 119895 = 1 to 119899 do(20) if rand le PR then(21) Update the current entry according to (10)(22) end if(23) end for(24) end if(25) end for(26) Update current best(27) end while(28) Return Updated population and global best solution(29) end for
isin [119886119895119887119895] is the 119895th coordinate of solution 119883
119894 119886119895
and 119887119895are the 119895th coordinates of the bounds describing the
search space of the problem and 120578119895isin (0 1) This means that
119883119894= [119909119894119895] for 119895 = 1 119899 represents the position of the 119895th
seed in population pop
4 Experimental Settings and Discussion
In our experiments we tested SbPPA against some recentlydeveloped algorithms and somewell established and standardones Our set of test problems includes benchmark con-strained and unconstrained optimization problems [22 3048 49] The results are compared in terms of statistics (bestworst mean and standard deviation) for solutions obtained
by SbPPA ABC [18 50] PSO [51] FF [21] HPA [29] SSO-C [22] Classical Evolutionary Programming (CEP) [30] andFast Evolutionary Programming (FEP) [30] The detaileddescriptions of these problems are given in Appendices
In Tables 4 and 7 the significance of results is shown interms of wintieloss (see Table 2 in [52]) according to thefollowing notations
(i) (+) when SbPPA is better(ii) (asymp) when the results are approximately the same as
those obtained with SbPPA(iii) (minus) when SbPPA is worse
Moreover in Tables 5 and 6 the significance of resultsobtained with SbPPA is highlighted
41 Parameter Settings The parameter settings are given inTables 1ndash3
5 Conclusion
In this paper a new metaheuristic referred to as the Seed-Based Plant Propagation Algorithm (SbPPA) [47] has beenproposed Plants have evolved a variety of ways to propagatePropagation through seeds is perhaps the most commonof them all and one which takes advantage of all sorts of
agents ranging from wind to water birds and animals Thestrawberry plant uses both runners and seeds to propagateHere we consider the propagation through seeds that thestrawberry plant has evolved to design an efficient optimiza-tion algorithm
To capture the dispersal process we adopt a queuingapproach which given the extent of fruit produced indicatesthe extent of seeds dispersed and hence the effectivenessof the searchoptimization algorithm based on this processLooking at the random process of agents using the plants(feeding station) it is reasonable to assume that it is of thePoisson type On the other hand the time taken by agentsin successfully eating fruit and leaving to disperse its seedsin other words the service time for agents is expressed by arandom variable which follows the exponential probabilitydistribution To this end we choose a threshold value ofthe Poisson probability that dictates how much explorationand exploitation are done during the search An alternativestrategy has been adopted here This strategy consists inmaking sure that the initial population is as good as theuser can afford it to be by using best solutions found sofar The effects of this strategy on convergence are shownthrough convergence plots of Figures 4 and 5 for some of thesolved problems SbPPA is easy to implement as it requiresless arbitrary parameter settings than other algorithms Thesuccess rate of SbPPA increases as it gets its population of bestsolutions It has been implemented for both unconstrainedand constrained optimization problems Its performancecompared to that of other algorithms points to SbPPA asbeing superior
Appendices
A Unconstrained GlobalOptimization Problems
See Tables 8 and 9
B Set of Constrained Global OptimizationProblems Used in Our Experiments
B1 CP1 Consider the following
Min 119891 (119909) = 5
4
sum
119889=1
119909119889minus 5
4
sum
119889=1
1199092
119889minus
13
sum
119889=5
119909119889
subject to 1198921(119909) = 2119909
1+ 21199092+ 11990910+ 11990911minus 10 le 0
1198922(119909) = 2119909
1+ 21199093+ 11990910+ 11990912minus 10 le 0
1198923(119909) = 2119909
2+ 21199093+ 11990911+ 11990912minus 10 le 0
1198924(119909) = minus8119909
1+ 11990910le 0
1198925(119909) = minus8119909
2+ 11990911le 0
1198926(119909) = minus8119909
3+ 11990912le 0
1198927(119909) = minus2119909
4minus 1199095+ 11990910le 0
1198928(119909) = minus2119909
6minus 1199097+ 11990911le 0
1198929(119909) = minus2119909
8minus 1199099+ 11990912le 0
(B1)
where bounds are 0 le 119909119894le 1 (119894 = 1 9 13) 0 le
119909119894le 100 (119894 = 10 11 12) The global optimum is at 119909lowast =
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
(1) NPlarr Population size 119903 larr Counter of trial runs MaxExp larrMaximum experiments(2) for 119903 = 1 MaxExp do(3) if 119903 le NP then(4) Create a random population of seeds pop = 119883
119894| 119894 = 1 2 NP
using (12) and collect the best solutions from each trial run in popbest(5) Evaluate the population pop(6) end if(7) while 119903 gt NP do(8) Use updated population popbest(9) end while(10) while (the stopping criteria is not satisfied) do(11) for 119894 = 1 to NP do(12) if Poiss(120582)
119894ge 005 then ⊳ (Global or local seed dispersion)
(13) for 119895 = 1 to 119899 do ⊳ (119899 is number of dimensions)(14) if rand le PR then ⊳ (PR = Perturbation Rate)(15) Update the current entry according to (11)(16) end if(17) end for(18) else(19) for 119895 = 1 to 119899 do(20) if rand le PR then(21) Update the current entry according to (10)(22) end if(23) end for(24) end if(25) end for(26) Update current best(27) end while(28) Return Updated population and global best solution(29) end for
isin [119886119895119887119895] is the 119895th coordinate of solution 119883
119894 119886119895
and 119887119895are the 119895th coordinates of the bounds describing the
search space of the problem and 120578119895isin (0 1) This means that
119883119894= [119909119894119895] for 119895 = 1 119899 represents the position of the 119895th
seed in population pop
4 Experimental Settings and Discussion
In our experiments we tested SbPPA against some recentlydeveloped algorithms and somewell established and standardones Our set of test problems includes benchmark con-strained and unconstrained optimization problems [22 3048 49] The results are compared in terms of statistics (bestworst mean and standard deviation) for solutions obtained
by SbPPA ABC [18 50] PSO [51] FF [21] HPA [29] SSO-C [22] Classical Evolutionary Programming (CEP) [30] andFast Evolutionary Programming (FEP) [30] The detaileddescriptions of these problems are given in Appendices
In Tables 4 and 7 the significance of results is shown interms of wintieloss (see Table 2 in [52]) according to thefollowing notations
(i) (+) when SbPPA is better(ii) (asymp) when the results are approximately the same as
those obtained with SbPPA(iii) (minus) when SbPPA is worse
Moreover in Tables 5 and 6 the significance of resultsobtained with SbPPA is highlighted
41 Parameter Settings The parameter settings are given inTables 1ndash3
5 Conclusion
In this paper a new metaheuristic referred to as the Seed-Based Plant Propagation Algorithm (SbPPA) [47] has beenproposed Plants have evolved a variety of ways to propagatePropagation through seeds is perhaps the most commonof them all and one which takes advantage of all sorts of
agents ranging from wind to water birds and animals Thestrawberry plant uses both runners and seeds to propagateHere we consider the propagation through seeds that thestrawberry plant has evolved to design an efficient optimiza-tion algorithm
To capture the dispersal process we adopt a queuingapproach which given the extent of fruit produced indicatesthe extent of seeds dispersed and hence the effectivenessof the searchoptimization algorithm based on this processLooking at the random process of agents using the plants(feeding station) it is reasonable to assume that it is of thePoisson type On the other hand the time taken by agentsin successfully eating fruit and leaving to disperse its seedsin other words the service time for agents is expressed by arandom variable which follows the exponential probabilitydistribution To this end we choose a threshold value ofthe Poisson probability that dictates how much explorationand exploitation are done during the search An alternativestrategy has been adopted here This strategy consists inmaking sure that the initial population is as good as theuser can afford it to be by using best solutions found sofar The effects of this strategy on convergence are shownthrough convergence plots of Figures 4 and 5 for some of thesolved problems SbPPA is easy to implement as it requiresless arbitrary parameter settings than other algorithms Thesuccess rate of SbPPA increases as it gets its population of bestsolutions It has been implemented for both unconstrainedand constrained optimization problems Its performancecompared to that of other algorithms points to SbPPA asbeing superior
Appendices
A Unconstrained GlobalOptimization Problems
See Tables 8 and 9
B Set of Constrained Global OptimizationProblems Used in Our Experiments
B1 CP1 Consider the following
Min 119891 (119909) = 5
4
sum
119889=1
119909119889minus 5
4
sum
119889=1
1199092
119889minus
13
sum
119889=5
119909119889
subject to 1198921(119909) = 2119909
1+ 21199092+ 11990910+ 11990911minus 10 le 0
1198922(119909) = 2119909
1+ 21199093+ 11990910+ 11990912minus 10 le 0
1198923(119909) = 2119909
2+ 21199093+ 11990911+ 11990912minus 10 le 0
1198924(119909) = minus8119909
1+ 11990910le 0
1198925(119909) = minus8119909
2+ 11990911le 0
1198926(119909) = minus8119909
3+ 11990912le 0
1198927(119909) = minus2119909
4minus 1199095+ 11990910le 0
1198928(119909) = minus2119909
6minus 1199097+ 11990911le 0
1198929(119909) = minus2119909
8minus 1199099+ 11990912le 0
(B1)
where bounds are 0 le 119909119894le 1 (119894 = 1 9 13) 0 le
119909119894le 100 (119894 = 10 11 12) The global optimum is at 119909lowast =
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
agents ranging from wind to water birds and animals Thestrawberry plant uses both runners and seeds to propagateHere we consider the propagation through seeds that thestrawberry plant has evolved to design an efficient optimiza-tion algorithm
To capture the dispersal process we adopt a queuingapproach which given the extent of fruit produced indicatesthe extent of seeds dispersed and hence the effectivenessof the searchoptimization algorithm based on this processLooking at the random process of agents using the plants(feeding station) it is reasonable to assume that it is of thePoisson type On the other hand the time taken by agentsin successfully eating fruit and leaving to disperse its seedsin other words the service time for agents is expressed by arandom variable which follows the exponential probabilitydistribution To this end we choose a threshold value ofthe Poisson probability that dictates how much explorationand exploitation are done during the search An alternativestrategy has been adopted here This strategy consists inmaking sure that the initial population is as good as theuser can afford it to be by using best solutions found sofar The effects of this strategy on convergence are shownthrough convergence plots of Figures 4 and 5 for some of thesolved problems SbPPA is easy to implement as it requiresless arbitrary parameter settings than other algorithms Thesuccess rate of SbPPA increases as it gets its population of bestsolutions It has been implemented for both unconstrainedand constrained optimization problems Its performancecompared to that of other algorithms points to SbPPA asbeing superior
Appendices
A Unconstrained GlobalOptimization Problems
See Tables 8 and 9
B Set of Constrained Global OptimizationProblems Used in Our Experiments
B1 CP1 Consider the following
Min 119891 (119909) = 5
4
sum
119889=1
119909119889minus 5
4
sum
119889=1
1199092
119889minus
13
sum
119889=5
119909119889
subject to 1198921(119909) = 2119909
1+ 21199092+ 11990910+ 11990911minus 10 le 0
1198922(119909) = 2119909
1+ 21199093+ 11990910+ 11990912minus 10 le 0
1198923(119909) = 2119909
2+ 21199093+ 11990911+ 11990912minus 10 le 0
1198924(119909) = minus8119909
1+ 11990910le 0
1198925(119909) = minus8119909
2+ 11990911le 0
1198926(119909) = minus8119909
3+ 11990912le 0
1198927(119909) = minus2119909
4minus 1199095+ 11990910le 0
1198928(119909) = minus2119909
6minus 1199097+ 11990911le 0
1198929(119909) = minus2119909
8minus 1199099+ 11990912le 0
(B1)
where bounds are 0 le 119909119894le 1 (119894 = 1 9 13) 0 le
119909119894le 100 (119894 = 10 11 12) The global optimum is at 119909lowast =
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
agents ranging from wind to water birds and animals Thestrawberry plant uses both runners and seeds to propagateHere we consider the propagation through seeds that thestrawberry plant has evolved to design an efficient optimiza-tion algorithm
To capture the dispersal process we adopt a queuingapproach which given the extent of fruit produced indicatesthe extent of seeds dispersed and hence the effectivenessof the searchoptimization algorithm based on this processLooking at the random process of agents using the plants(feeding station) it is reasonable to assume that it is of thePoisson type On the other hand the time taken by agentsin successfully eating fruit and leaving to disperse its seedsin other words the service time for agents is expressed by arandom variable which follows the exponential probabilitydistribution To this end we choose a threshold value ofthe Poisson probability that dictates how much explorationand exploitation are done during the search An alternativestrategy has been adopted here This strategy consists inmaking sure that the initial population is as good as theuser can afford it to be by using best solutions found sofar The effects of this strategy on convergence are shownthrough convergence plots of Figures 4 and 5 for some of thesolved problems SbPPA is easy to implement as it requiresless arbitrary parameter settings than other algorithms Thesuccess rate of SbPPA increases as it gets its population of bestsolutions It has been implemented for both unconstrainedand constrained optimization problems Its performancecompared to that of other algorithms points to SbPPA asbeing superior
Appendices
A Unconstrained GlobalOptimization Problems
See Tables 8 and 9
B Set of Constrained Global OptimizationProblems Used in Our Experiments
B1 CP1 Consider the following
Min 119891 (119909) = 5
4
sum
119889=1
119909119889minus 5
4
sum
119889=1
1199092
119889minus
13
sum
119889=5
119909119889
subject to 1198921(119909) = 2119909
1+ 21199092+ 11990910+ 11990911minus 10 le 0
1198922(119909) = 2119909
1+ 21199093+ 11990910+ 11990912minus 10 le 0
1198923(119909) = 2119909
2+ 21199093+ 11990911+ 11990912minus 10 le 0
1198924(119909) = minus8119909
1+ 11990910le 0
1198925(119909) = minus8119909
2+ 11990911le 0
1198926(119909) = minus8119909
3+ 11990912le 0
1198927(119909) = minus2119909
4minus 1199095+ 11990910le 0
1198928(119909) = minus2119909
6minus 1199097+ 11990911le 0
1198929(119909) = minus2119909
8minus 1199099+ 11990912le 0
(B1)
where bounds are 0 le 119909119894le 1 (119894 = 1 9 13) 0 le
119909119894le 100 (119894 = 10 11 12) The global optimum is at 119909lowast =
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
agents ranging from wind to water birds and animals Thestrawberry plant uses both runners and seeds to propagateHere we consider the propagation through seeds that thestrawberry plant has evolved to design an efficient optimiza-tion algorithm
To capture the dispersal process we adopt a queuingapproach which given the extent of fruit produced indicatesthe extent of seeds dispersed and hence the effectivenessof the searchoptimization algorithm based on this processLooking at the random process of agents using the plants(feeding station) it is reasonable to assume that it is of thePoisson type On the other hand the time taken by agentsin successfully eating fruit and leaving to disperse its seedsin other words the service time for agents is expressed by arandom variable which follows the exponential probabilitydistribution To this end we choose a threshold value ofthe Poisson probability that dictates how much explorationand exploitation are done during the search An alternativestrategy has been adopted here This strategy consists inmaking sure that the initial population is as good as theuser can afford it to be by using best solutions found sofar The effects of this strategy on convergence are shownthrough convergence plots of Figures 4 and 5 for some of thesolved problems SbPPA is easy to implement as it requiresless arbitrary parameter settings than other algorithms Thesuccess rate of SbPPA increases as it gets its population of bestsolutions It has been implemented for both unconstrainedand constrained optimization problems Its performancecompared to that of other algorithms points to SbPPA asbeing superior
Appendices
A Unconstrained GlobalOptimization Problems
See Tables 8 and 9
B Set of Constrained Global OptimizationProblems Used in Our Experiments
B1 CP1 Consider the following
Min 119891 (119909) = 5
4
sum
119889=1
119909119889minus 5
4
sum
119889=1
1199092
119889minus
13
sum
119889=5
119909119889
subject to 1198921(119909) = 2119909
1+ 21199092+ 11990910+ 11990911minus 10 le 0
1198922(119909) = 2119909
1+ 21199093+ 11990910+ 11990912minus 10 le 0
1198923(119909) = 2119909
2+ 21199093+ 11990911+ 11990912minus 10 le 0
1198924(119909) = minus8119909
1+ 11990910le 0
1198925(119909) = minus8119909
2+ 11990911le 0
1198926(119909) = minus8119909
3+ 11990912le 0
1198927(119909) = minus2119909
4minus 1199095+ 11990910le 0
1198928(119909) = minus2119909
6minus 1199097+ 11990911le 0
1198929(119909) = minus2119909
8minus 1199099+ 11990912le 0
(B1)
where bounds are 0 le 119909119894le 1 (119894 = 1 9 13) 0 le
119909119894le 100 (119894 = 10 11 12) The global optimum is at 119909lowast =
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
B6 Welded Beam Design Optimisation The welded beamdesign is a standard test problem for constrained designoptimisation [55 56] There are four design variables thewidth 119908 and length 119871 of the welded area and the depth119889 and thickness ℎ of the main beam The objective is tominimise the overall fabrication cost under the appropriateconstraints of shear stress 120591 bending stress 120590 buckling load119875 and maximum end deflection 120575 The optimization modelis summarized as follows where 119909119879 = (119908 119871 119889 ℎ)
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
B7 Speed Reducer Design Optimization The problem ofdesigning a speed reducer [57] is a standard test problemIt consists of the design variables as face width 119909
1 module
of teeth 1199092 number of teeth on pinion 119909
3 length of the first
shaft between bearings 1199094 length of the second shaft between
bearings 1199095 diameter of the first shaft 119909
6 and diameter
of the first shaft 1199097(all variables are continuous except 119909
3
that is integer) The weight of the speed reducer is to beminimized subject to constraints on bending stress of the gearteeth surface stress transverse deflections of the shafts andstresses in the shaft [55]Themathematical formulation of theproblem where 119909119879 = (119909
1 1199092 1199093 1199094 1199095 1199096 1199097) is as follows
Minimise 119891 (119909) = 0785411990911199092
2
sdot (333331199092
3+ 149334119909
3430934)
minus 15081199091(1199092
6+ 1199093
7)
+ 74777 (1199093
6+ 1199093
7)
+ 07854 (11990941199092
6+ 11990951199092
7)
subject to 1198921(119909) =
27
11990911199092
21199093
minus 1 le 0
1198922(119909) =
3975
11990911199092
21199092
3
minus 1 le 0
1198923(119909) =
1931199093
4
119909211990931199094
6
minus 1 le 0
1198924(119909) =
1931199093
5
119909211990931199094
7
minus 1 le 0
1198925(119909) =
10
1101199093
6
radic(7450119909
4
11990921199093
)
2
+ 169 times 106
minus 1 le 0
1198926(119909) =
10
851199093
7
radic(7450119909
5
11990921199093
)
2
+ 1575 times 106
minus 1 le 0
1198927(119909) =
11990921199093
40minus 1 le 0
1198928(119909) =
51199092
1199091
minus 1 le 0
1198929(119909) =
1199091
121199092
minus 1 le 0
11989210(119909) =
151199096+ 19
1199094
minus 1 le 0
11989211(119909) =
111199097+ 19
1199095
minus 1 le 0
(B8)
The simple limits on the design variables are 26 le 1199091le 36
07 le 1199092le 08 17 le 119909
3le 28 73 le 119909
4le 83 78 le 119909
5le 83
29 le 1199096le 39 and 50 le 119909
7le 55
B8 Spring Design Optimisation The main objective of thisproblem [58 59] is to minimize the weight of a ten-sioncompression string subject to constraints of minimumdeflection shear stress surge frequency and limits on outsidediameter and on design variables There are three designvariables the wire diameter 119909
1 the mean coil diameter 119909
2
and the number of active coils 1199093[55] The mathematical
formulation of this problem where 119909119879 = (1199091 1199092 1199093) is as
follows
Minimize 119891 (119909) = (1199093+ 2) 119909
21199092
1
subject to 1198921(119909) = 1 minus
1199093
21199093
7 1781199094
1
le 0
1198922(119909) =
41199092
2minus 11990911199092
12 566 (11990921199093
1) minus 1199094
1
+1
5 1081199092
1
minus 1 le 0
1198923(119909) = 1 minus
140451199091
1199092
21199093
le 0
1198924(119909) =
1199092+ 1199091
15minus 1 le 0
(B9)
The simple limits on the design variables are 005 le 1199091le 20
025 le 1199092le 13 and 20 le 119909
3le 150
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper
The Scientific World Journal 15
Acknowledgments
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
The authors are grateful to anonymous reviewers for theirvaluable reviews and constructive criticism on earlier ver-sion of this paper This work is supported by Abdul WaliKhan University Mardan Pakistan Grant no F16-5PampDAWKUM238
References
[1] A Salhi and E S Fraga ldquoNature-inspired optimisationapproaches and the new plant propagation algorithmrdquo in Pro-ceedings of the International Conference on Numerical Analysisand Optimization (ICeMATH rsquo11) pp K2-1ndashK2-8 YogyakartaIndonesia 2011
[2] C M Herrera and O Pellmyr Plant Animal Interactions AnEvolutionary Approach John Wiley amp Sons 2009
[3] C M Herrera ldquoSeed dispersal by vertebratesrdquo in Plant-AnimalInteractions An Evolutionary Approach pp 185ndash208 2002
[4] W G Abrahamson and T N Taylor Plant-Animal InteractionsMcGraw Hill 1989
[5] A N Andersen and R W Braithwaite ldquoPlant-animal inter-actionsrdquo in Landscape and Vegetation Ecology of the KakaduRegion Northern Australia pp 137ndash154 Springer 1996
[6] J P Bryant ldquoPlant-animal interactionsrdquo Environmental Ento-mology vol 19 no 4 pp 1169ndash1170 1990
[7] B J GloverUnderstanding Flowers and Flowering an IntegratedApproach Oxford University Press Oxford UK 2007
[8] X-S Yang ldquoFlower pollination algorithm for global optimiza-tionrdquo in Unconventional Computation and Natural Computa-tion pp 240ndash249 Springer 2012
[9] X-S YangMKaramanoglu andXHe ldquoMulti-objective floweralgorithm for optimizationrdquo in Proceedings of the 13th AnnualInternational Conference on Computational Science (ICCS rsquo13)vol 18 pp 861ndash868 June 2013
[10] H J du Plessis R J Brand C Glyn-Woods and M AGoedhart ldquoEfficient genetic transformation of strawberry (Fra-garia x Ananassa Duch) cultivar selektardquo in III InternationalSymposium on In Vitro Culture andHorticultural Breeding ISHSActa Horticulturae 447 pp 289ndash294 ISHS 1996
[11] L W Krefting and E I Roe ldquoThe role of some birds andmammals in seed germinationrdquo Ecological Monographs vol 19no 3 pp 269ndash286 1949
[12] D G Wenny and D J Levey ldquoDirected seed dispersal bybellbirds in a tropical cloud forestrdquo Proceedings of the NationalAcademy of Sciences of the United States of America vol 95 no11 pp 6204ndash6207 1998
[14] X-S Yang Nature-Inspired Metaheuristic Algorithms LuniverPress Beckington UK 2011
[15] M Kazemian Y Ramezani C Lucas and B Moshiri ldquoSwarmclustering based onflowers pollination by artificial beesrdquo Studiesin Computational Intelligence vol 34 pp 191ndash202 2006
[16] R Eberhart and J Kennedy ldquoA new optimizer using particleswarm theoryrdquo in Proceedings of the 6th International Sympo-sium on Micro Machine and Human Science (MHS rsquo95) pp 39ndash43 IEEE Nagoya Japan October 1995
[17] M Clerc Particle Swarm Optimization vol 93 John Wiley ampSons 2010
[18] D Karaboga ldquoAn idea based on honey bee swarm for numericaloptimizationrdquo Tech Rep TR06 Erciyes University Press Kay-seri Turkey 2005
[19] D Karaboga and B Basturk ldquoOn the performance of artificialbee colony (ABC) algorithmrdquo Applied Soft Computing Journalvol 8 no 1 pp 687ndash697 2008
[20] X-S Yang ldquoFirefly algorithm stochastic test functions anddesign optimisationrdquo International Journal of Bio-Inspired Com-putation vol 2 no 2 pp 78ndash84 2010
[21] A H Gandomi X-S Yang and A H Alavi ldquoMixed variablestructural optimization using Firefly Algorithmrdquo Computersand Structures vol 89 no 23-24 pp 2325ndash2336 2011
[22] E Cuevas and M Cienfuegos ldquoA new algorithm inspired inthe behavior of the social-spider for constrained optimizationrdquoExpert Systems with Applications vol 41 no 2 pp 412ndash4252014
[23] M Sulaiman A Salhi B I Selamoglu and O B KirikchildquoA plant propagation algorithm for constrained engineeringoptimisation problemsrdquoMathematical Problems in Engineeringvol 2014 Article ID 627416 10 pages 2014
[24] J L Tellerıa A Ramırez and J Perez-Tris ldquoConservationof seed-dispersing migrant birds in Mediterranean habitatsshedding light on patterns to preserve processesrdquo BiologicalConservation vol 124 no 4 pp 493ndash502 2005
[25] Wikipedia Contributors ldquoStrawberryrdquo 2015 httpbitly17REoNP
[26] Ruth Strawberries eaten by animals 2012 httpbitly1zZg0jS[27] Anisognathus ldquoBlue Winged Mountain-Tanager eating tree
strawberries at San isidro Lodgerdquo 2012 httpbitly1we3thV[28] lifeisfullineedexchangecom ldquoLooking for nice strawberries
to pick among the tiny ones rotten ones and thistlesrdquo 2014httpslifeisfullwordpresscompage5
[29] M S Kıran and M Gunduz ldquoA recombination-basedhybridization of particle swarm optimization and artificialbee colony algorithm for continuous optimization problemsrdquoApplied Soft Computing Journal vol 13 no 4 pp 2188ndash22032013
[30] X Yao Y Liu and G Lin ldquoEvolutionary programming madefasterrdquo IEEE Transactions on Evolutionary Computation vol 3no 2 pp 82ndash102 1999
[31] N E Stork and C H C Lyal ldquoExtinction or lsquoco-extinctionrsquoratesrdquo Nature vol 366 no 6453 p 307 1993
[32] P Jordano ldquoFruits and frugivoryrdquo in Seeds The Ecology ofRegeneration in Plant Communities vol 2 pp 125ndash166 CABIWallingford UK 2000
[33] M Debussche and P Isenmann ldquoBird-dispersed seed rain andseedling establishment in patchy Mediterranean vegetationrdquoOikos vol 69 no 3 pp 414ndash426 1994
[34] D H Janzen ldquoHerbivores and the number of tree species intropical forestsrdquoThe American Naturalist vol 104 no 940 pp501ndash528 1970
[35] S A Levin ldquoPopulation dynamic models in heterogeneousenvironmentsrdquo Annual Review of Ecology and Systematics vol7 no 1 pp 287ndash310 1976
[36] S A H Geritz T J de Jong and P G L Klinkhamer ldquoTheefficacy of dispersal in relation to safe site area and seedproductionrdquo Oecologia vol 62 no 2 pp 219ndash221 1984
[37] S A Levin D Cohen and A Hastings ldquoDispersal strategiesin patchy environmentsrdquoTheoretical Population Biology vol 26no 2 pp 165ndash191 1984
16 The Scientific World Journal
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985
[38] C K Augspurger and S E Franson ldquoWind dispersal of artificialfruits varying in mass area and morphologyrdquo Ecology vol 68no 1 pp 27ndash42 1987
[39] R B Cooper Introduction to Queueing Theory 1972[40] J A Lawrence and B A Pasternack Applied Management
Science Wiley New York NY USA 2002[41] A H S Ang andWH Tang ldquoProbability concepts in engineer-
ingrdquo Planning vol 1 no 4 pp 1ndash3 2004[42] J D C Little ldquoA proof for the queuing formula 119871 = 120582119882rdquo
Operations Research vol 9 no 3 pp 383ndash387 1961[43] D WThompson On Growth and Form Courier 1942[44] K S van Houtan S L Pimm J M Halley R O Bierregaard Jr
and T E Lovejoy ldquoDispersal of Amazonian birds in continuousand fragmented forestrdquo Ecology Letters vol 10 no 3 pp 219ndash229 2007
[45] A M Reynolds and M A Frye ldquoFree-flight odor tracking inDrosophila is consistent with an optimal intermittent scale-freesearchrdquo PLoS ONE vol 2 no 4 article e354 2007
[46] Wikipedia Contributors Gamma function 2015 httpbitly1w6scza
[47] M Sulaiman and A Salhi ldquoA seed-based plant propagationalgorithm the feeding station modelrdquo in Proceedings of the 5thInternational Conference on Metaheuristics and Nature InspiredComputing Marrakesh Morocco October 2014 httpmeta-2014sciencesconforg40158
[48] P N Suganthan N Hansen J J Liang et al ldquoProblemdefinitions and evaluation criteria for the cec 2005 specialsession on real-parameter optimizationrdquo Tech Rep NanyangTechnological University Singapore 2005
[49] J J Liang T Runarsson E Mezura-Montes et al ldquoProblemdefinitions and evaluation criteria for the CEC 2006 specialsession on constrained real-parameter optimizationrdquo Journal ofApplied Mechanics vol 41 p 8 2006
[50] D Karaboga and B Akay ldquoAmodified artificial bee colony (abc)algorithm for constrained optimization problemsrdquo Applied SoftComputing Journal vol 11 no 3 pp 3021ndash3031 2011
[51] Q He and LWang ldquoA hybrid particle swarm optimization witha feasibility-based rule for constrained optimizationrdquo AppliedMathematics and Computation vol 186 no 2 pp 1407ndash14222007
[52] H-Y Wang X-J Ding Q-C Cheng and F-H Chen ldquoAnimproved isomap for visualization and classification of multiplemanifoldsrdquo inNeural Information Processing M Lee A HiroseZ-G Hou and R M Kil Eds vol 8227 of Lecture Notes inComputer Science pp 1ndash12 Springer Berlin Germany 2013
[53] D B Fogel System Identification Through Simulated EvolutionAMachine Learning Approach to Modeling Ginn Press 1991
[54] T Back and H-P Schwefel ldquoAn overview of evolutionary algo-rithms for parameter optimizationrdquo Evolutionary Computationvol 1 no 1 pp 1ndash23 1993
[55] L C Cagnina S C Esquivel and C A C Coello ldquoSolvingengineering optimization problemswith the simple constrainedparticle swarm optimizerrdquo Informatica vol 32 no 3 pp 319ndash326 2008
[56] X-S Yang and S Deb ldquoEngineering optimisation by cuckoosearchrdquo International Journal of Mathematical Modelling andNumerical Optimisation vol 1 no 4 pp 330ndash343 2010
[57] J Golinski ldquoAn adaptive optimization system applied tomachine synthesisrdquoMechanism and MachineTheory vol 8 no4 pp 419ndash436 1973
[58] J S Arora Introduction to Optimum Design Academic Press2004
[59] A D Belegundu and J S Arora ldquoA study of mathematicalprogramming methods for structural optimization I TheoryrdquoInternational Journal for Numerical Methods in Engineering vol21 no 9 pp 1583ndash1599 1985