Open Access
Issue
Int. J. Simul. Multisci. Des. Optim.
Volume 5, 2014
Article Number A11
Number of page(s) 10
DOI https://doi.org/10.1051/smdo/2013008
Published online 10 February 2014

© V. Savsani, Published by EDP Sciences, 2014

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Stiffeners are required to improve the load carrying capacity of structures. Stiffeners help to achieve lightweight and robust design of the structure. For this purpose they have wide use in structural engineering domain. Specially, stiffened plates are used in critical and sensitive structures such as in aircrafts, ship hulls and box girders in which safety and a perfect design is crucial. Buckling is the one of the most complex phenomenon that is inevitable for heavily axially loaded stiffened plate structures. The wide area of application of such structures is an aircraft industry where light, stiff and highly resistant structures are required. It is possible to manufacture aircraft elements as thin panels reinforced by stiffeners. The proper stiffeners arrangement in a domain of the structure decides about the effectiveness of the construction or about the effectiveness of reinforcement. Optimal properties of structures can be searched using optimization tools which are discussed in the later part of this paper.

For steel structures a general total cost function KT can be defined in the following form [19] K T = K m + K f + K t + K e , $$ {K}_{\mathrm{T}}={K}_{\mathrm{m}}+{K}_{\mathrm{f}}+{K}_{\mathrm{t}}+{K}_{\mathrm{e}}, $$(1)where Km is the material cost of structural members, Kf is the fabrication cost, Kt is the cost of transporting the fabricated pieces to the construction field, and Ke is the erection cost. Many studies are reported in literature for the cost optimization of different structures.

Anderson and Chong [4] presented the minimum cost design of homogeneous and hybrid stiffened steel plate girders according to the AISC code [2]. Douty [13] described the minimum cost design of three different types of bolted and welded connections used in steel frames based on the AISC specifications [1] The cost function is presented in terms of the connection variables, such as the diameter of the bolts, flange plate width and thickness, shear plate length and thickness, and the leg size of the fillet weld for connecting the shear plate to the column flange and for the flange plate moment connection. Lorenz [25] discussed the minimum cost design of composite beams based on the AISC load and resistance factor design (LRFD) code [3]. Farkas [15] presents closed-form solutions for optimum cost values of the cross-sectional variables for simply supported welded box girders subjected to a uniformly distributed load and constraints on bending stress, local flange and web buckling, shear fatigue for longitudinal fillet welds, and deflection. Farkas and Jarmai [16] presented the minimum cost design of laterally loaded welded cellular steel plates using three different approaches: the backtracking method, the hill-climbing method, and feasible sequential quadratic programming. Bhatti [7] presented the minimum cost design of simply supported composite I-shaped steel beams with concrete slabs subjected to a uniformly distributed load, and strength, deflection, and vibration constraints using the Lagrange multiplier approach. Farkas [18] derived a differential equation for the radial shell deformation using the bending theory of cylindrical shells. The solution of this equation gives an approximate formula for the maximum radial deformation. Comparing this formula with the limiting deformation given by ECCS allows calculating the required shell thickness. Another thickness can be obtained using ECCS buckling formulae. Using an interpolation, a thickness can be calculated, which fulfills the constraints. Jarmai et al. [21] demonstrated the optimal design of a cylindrical orthogonally stiffened shell member of an offshore fixed platform truss, loaded by axial compression and external pressure. Ring stiffeners of welded box section and stringers of halved rolled I-section were used. The design variables considered in the optimization were the shell thickness as well as the dimensions and numbers of stiffeners. The design constraints relate to the shell, panel ring and panel stringer buckling, as well as manufacturing limitations. The cost function includes the cost of material, forming of plate elements into cylindrical shape, welding and painting. In the optimization a number of relatively new mathematical optimization methods such as leap-frog – LFOPC, Dynamic-Q, ETOPC, and Particle swarm optimization – PSO were used. Jarmai [20] presented the single and multiobjective optimization of a welded stringer-stiffened cylindrical steel shell. A column fixed at the bottom and free on the top was constructed of stringer-stiffened cylindrical shell and loaded by axial compression as well as by horizontal force acting on the top. Halved rolled I-section stringers were welded outside of the shell by longitudinal fillet welds. Design Variables taken were the shell thickness as well as dimension and number of stringers. Three different objective functions were considered viz. material cost, forming cost to form the shell elements into the cylindrical shape, assembly and welding, and painting cost. The original PSO algorithm was modified to handle multiobjective optimization techniques and to find discrete values of design variables. Luis et al. [26] presented the optimization of a uniaxially compressed stiffened plate subjected to static and fatigue loading. The design variables were the thickness of the base plate, the number and dimension of stiffeners of the orthogonally stiffened plate. Cost optimization was carried out considering material, assembly, welding and painting costs. Farkas and Jarmai [14] presented the cost minimization of a cantilever stub column of square box cross section with orthogonally stiffened side plates. Design variables include the thickness and width of the side plates as well as the dimensions and numbers of stiffeners. PSO was implemented for the cost optimization. Marcin et al. [27] presented the optimal configuration of stiffening ribs in the steel welded I-beams and columns. The configuration with diagonal ribs is proved better than traditional ones. Farkas et al. [17] used refined version of the original Snyman-Fatti (SF) global continuous optimization algorithm to the optimal design of welded square stiffened plates. Square plates of square symmetry subjected to uniformly distributed normal static loads, supported at four corners, and stiffened by a square symmetrical orthogonal grid of ribs was investigated. Halved rolled I-section stiffeners are used welded to the base plate by double fillet welds. The cost function includes material, welding as well as painting costs, and is formulated according to the fabrication sequence. Design variables include base plate thickness as well as the dimensions of the edge and internal stiffeners. Constraints on stress in the base plate and in stiffeners, as well as on deflection of edge stiffeners and of internal stiffeners were considered. PSO algorithm was applied to confirm the results given by the SF algorithm.

It has been observed from the review of research work on design optimization of welded structures that the recently used advanced optimization technique is PSO. PSO requires many controlling parameters such as inertia weight, maximum velocity, and learning factors which influence the performance of the algorithm. Therefore the efforts are continuing to use more different optimization algorithms, which are more powerful, robust and able to provide accurate solution. This paper is intended to apply one such optimization algorithms known as artificial bee colony (ABC) [22] optimization technique. In the present work, an effort is made to verify if any improvement in the solution is possible by employing ABC optimization technique to two different cost optimization problems for welded structures reported in the literature. Moreover ABC optimization technique is the most recently developed technique. It is not so far tried for the complex optimization problems for structural design. The next section gives a brief description of ABC optimization techniques.

2 Artificial bee colony (ABC) techniques

Artificial bee colony (ABC) Algorithm is an optimization algorithm based on the intelligent foraging behaviour of honey bee swarm. The colony of artificial bees consists of three groups of bees: employed bees, onlookers and scouts [5, 22].

The colony of the artificial bees is divided into two groups, first half of the colony consists of the employed artificial bees and the second half includes the onlooker bees. Scout bees are the employed bee whose food source has been abandoned. In ABC algorithm, the position of a food source represents a possible solution to the optimization problem (value of design variables) and the nectar amount of a food source corresponds to the quality of the associated solution (value of objective function). At the first step, the ABC generates a randomly distributed initial population Pinitial of N solutions, where N denotes the size of population. Each solution xi is an S-dimensional vector where S is the number of optimization parameters (design variables). After initialization, the population of the solutions is subjected to repeated cycles, C = 1, 2, …, G, of the search processes of the employed bees, the onlooker bees and scout bees. An employed bee which is equal to N/2 produces a modification on the solution in her memory depending on the local information. If the objective function value (fitness) of the new solution is higher than that of the previous one, the bee memorizes the new position and forgets the old one. Otherwise she keeps the position of the previous one in her memory. After all employed bees complete the search process; they share the nectar information of the food sources and their position information with the onlooker bees on the dance area. An onlooker bee evaluates the fitness information taken from all employed bees and chooses a food source with a probability related to its fitness value. An onlooker bee also produces a new solution and it memorizes the new position if its fitness value is better than the previous position. An artificial onlooker bee chooses a food source depending on the probability value associated with that food source, pi, calculated by the following expression p i = F i n = 1 N b F n , $$ {p}_i=\frac{{F}_i}{\sum_{n=1}^{{N}_b}{F}_n}, $$(2)where Fi is the fitness value of the solution i which is proportional to the nectar amount of the food source in the position i and Nb is the number of food sources which is equal to the number of employed bees. In order to produce a candidate food position from the old one in memory, the ABC uses the following expression (3): v ij = x ij + R ij ( x ij - x kj ) , $$ {v}_{{ij}}={x}_{{ij}}+{R}_{{ij}}\left({x}_{{ij}}-{x}_{{kj}}\right), $$(3)where k ∈ {1, 2, …, N} and j ∈ {1, 2, …, D} are randomly chosen indexes. Although k is determined randomly, it has to be different from i. Rij is a random number between (−1, 1). It controls the production of neighbour food sources around xij and represents the comparison of two solutions. As can be seen from (3), as the difference between the parameters of the xij and xkj decreases, the perturbation on the position xij gets decrease, too. Thus, as the search approaches to the optimum solution in the search space, the step length is adaptively reduced.

If the position of the food source cannot be improved for some predetermined number of cycles than that food source is abandoned. The abandoned food source is replaced with a new food source by the scouts. In ABC, this is simulated by producing a position randomly and replacing it with the abandoned one. The value of predetermined number of cycles is an important control parameter of the ABC algorithm, which is called “limit” for abandonment. The value of limit is generally taken as Number of employed bees*S [23, 24].

Assume that the abandoned source is xi and j ∈ {1, 2, …, D}, then the scout discovers a new food source to be replaced with xi. This operation can be defined as in (4) x i j = x min j + rand ( 0,1 ) ( x max j - x min j ) . $$ {x}_i^j={x}_{\mathrm{min}}^j+\mathrm{rand}\left(\mathrm{0,1}\right)\left({x}_{\mathrm{max}}^j-{x}_{\mathrm{min}}^j\right). $$(4)

It is clear from the above explanation that there are three control parameters used in the ABC: The number of food sources which is equal to the number of employed or onlooker bees (N), the value of limit, the maximum cycle number (G).

As suggested by Karaboga [22] Rij is uniformly distributed random number between –1 and 1. In this paper ABC is modified by using normally distributed random number with mean 0 and variance 1 instead of uniformly distributed random number between –1 and 1. The distributions of both the random numbers generated using MATLAB with 1e7 random points are shown in Figures 1 and 2 respectively.

thumbnail Figure 1.

Normal distribution of 1e7 random numbers with mean 0 and variance 1.

thumbnail Figure 2.

Uniformly distributed random number between –1 and 1 for 1e7 random points.

Another modification is made in equation (3) by replacing (xij − xik) by (xbestj − xij). Here xbestj is the best solution from the population. The reason is that in the earlier formula the solution tries to move towards any random solution (xik). There is no guarantee for the xik to be better than xij. So solution can move towards worst solution also, which may require more computational time to reach the optimum solution. By replacing xik with xbestj, the solution will always move towards the best solution in each iteration, which will lead to global optimum solution in less computational effort. Moreover effort is made to increase the convergence rate of the original ABC algorithm. Many experiments were conducted to change equation (3) so that algorithm converges to optimum solution in less number of generations. In this paper original ABC will be referred as ABC_1 and modified ABC as ABC_2. Experiments shows that by using ABC_2 solutions improves by approximately 87% and by using ABC_1 solutions improves by 39% in successive first five generations which indicates the fast convergence rate of modified ABC. Convergence rate for the benchmark Sphere function ( f ( X ) = i = 1 n x i 2 , x [ - 100   100 ] ,   i = 20 ) $ f(X)=\sum_{i=1}^n{x}_i^2,x\in \left[-100\enspace 100\right],\enspace i=20)$ for five different runs for both the variants of ABC is shown in Figures 3 and 4. From Figures 3 and 4 also it can be noted that modified ABC have better convergence than classical ABC.

thumbnail Figure 3.

Variation of solution for the sphere function for classical ABC for five different runs.

thumbnail Figure 4.

Variation of solution for the sphere function for modified ABC for five different runs.

The algorithm discussed above takes real number in the aspect of representation solution and updates themselves with certain rules based on the behavior of algorithm. In this algorithm, process looks only for the best solution and tends to converge to the best solution. In the implementation process, solutions are randomly generated at the beginning and usually violate the system constraints resulting in infeasible solutions. Therefore, the handling of system constraints, particularly nonlinear equation constraints, and the measurement and evaluation of infeasible solutions is very important. To cope with constrained problems with evolutionary computation, various approaches such as rejection of infeasible individuals, repair of infeasible individuals, replacement of individuals by their repaired versions, and penalty function methods can be adopted. Among them, the penalty function method is particularly promising as evidenced by recent developments [12] and the same is adapted in the present work. Penalty function changes the constrained optimization problem into unconstrained optimization problem. For example consider an optimization problem as Minimize :   f ( X ) ,   Subjected   to :   g i ( X ) 0 . $$ \mathrm{Minimize}:\enspace f(X),\enspace \mathrm{Subjected}\enspace \mathrm{to}:\enspace {g}_i(X)\ge 0. $$(5)

This problem is converted in unconstrained form as Minimize :   f ( X ) ,   if   g i ( X ) 0 ,   else , $$ \mathrm{Minimize}:\enspace f(X),\enspace \mathrm{if}\enspace {g}_i(X)\ge 0,\enspace \mathrm{else}, $$(6) Minimize :   f ( X ) + R i = 1 n g i ( X ) , $$ \mathrm{Minimize}:\enspace f(X)+R\sum_{i=1}^n{g}_i(X), $$(7)where R is very large number and n is number of constraints.

The next section presents examples to demonstrate and validate the considered optimization techniques for the cost optimization of welded structures.

3. Examples

To check the performance of ABC some classical benchmark problems from Karaboga and Basturk [24] were experimented. In experiments, Schaffer function has two parameters (−100 ≤ X ≤ 100), while Sphere (−100 ≤ X ≤ 100), Griewank (−600 ≤ X ≤ 600), Rastrigin (−5.12 ≤ X ≤ 5.12) and Rosenbrock (−50 ≤ X ≤ 50) functions have 20 parameters each. Moreover some engineering benchmark problems were also taken for the investigation. They include pressure vessel problem [32] consisting of four design variable and four constraints, Spring weight minimization problem [6]) consisting of three design variable and four constraints, and speed reducer problem [28] consisting of 7 design variables and 11 constraints. As this paper focuses on the optimization of welded structures two special problems of welded structure optimization is taken from the literature and is given in detail as follows. More detailed information can be obtained from the mentioned reference. Experiments were conducted for different population size and number of generations. In this paper N50G200 indicates population size of 50 with number of generation of 200. Results were obtained for 50 different runs.

3.1 Example – 1 welded beam

The objective is to design a welded beam for minimum cost. There are four design variables as shown in Figure 5: h (x1), l (x2), t (x3) and b (x4). Design is subjected to the constraints on shear stress (τ), bending stress in the beam (σ), buckling load on the bar (Pc), end deflection of the beam (δ), and side constraints.

thumbnail Figure 5.

Welded beam design (Example 1).

The problem can be stated as follows: Minimize : f ( x ) = 1.10471 x 1 2 x 2 + 0.04811 x 3 x 4 ( 14.0 + x 2 ) . $$ \mathrm{Minimize}:f(x)=1.10471{x}_1^2{x}_2+0.04811{x}_3{x}_4\left(14.0+{x}_2\right). $$(8)

Subject to: g 1 ( x ) = τ ( x ) - τ max 0 , $$ {g}_1(x)=\tau (x)-{\tau }_{\mathrm{max}}\le 0, $$(9) g 2 ( x ) = σ ( x ) - σ max 0 , $$ {g}_2(x)=\sigma (x)-{\sigma }_{\mathrm{max}}\le 0, $$(10) g 3 ( x ) = x 1 - x 4 0 , $$ {g}_3(x)={x}_1-{x}_4\le 0, $$(11) g 4 ( x ) = 0.10471 x 1 2 + 0.04811 x 3 x 4 ( 14.0 + x 2 ) - 5.0 0 , $$ {g}_4(x)=0.10471{x}_1^2+0.04811{x}_3{x}_4(14.0+{x}_2)-5.0\le 0, $$(12) g 5 ( x ) = 0.125 - x 1 0 , $$ {g}_5(x)=0.125-{x}_1\le 0, $$(13) g 6 ( x ) = δ ( x ) - δ max 0 , $$ {g}_6(x)=\mathrm{\delta }(x)-{\mathrm{\delta }}_{\mathrm{max}}\le 0, $$(14) g 7 ( x ) = P - P c ( x ) 0 , $$ {g}_7(x)=P-{P}_c(x)\le 0, $$(15)

where τ ( x ) = ( τ ) 2 + 2 τ τ '' x 2 2 R + ( τ '' ) 2 , $$ \tau (x)=\sqrt{({\tau }^\mathrm{\prime}{)}^2+2{\tau }^\mathrm{\prime}{\tau }^\mathrm{\prime\prime}\frac{{x}_2}{2R}+({\tau }^\mathrm{\prime\prime}{)}^2}, $$(16) τ = P 2 x 1 x 2 , $$ {\tau }^\mathrm{\prime}=\frac{P}{\sqrt{2}{x}_1{x}_2}, $$(17) τ '' = MR J , $$ {\tau }^\mathrm{\prime\prime}=\frac{{MR}}{J}, $$(18) M = P ( L + x 2 2 ) , $$ M=P\left(L+\frac{{x}_2}{2}\right), $$(19) R = x 2 2 4 + ( x 1 + x 3 2 ) 2 , $$ R=\sqrt{\frac{{x}_2^2}{4}+{\left(\frac{{x}_1+{x}_3}{2}\right)}^2}, $$(20) J = 2 [ 2 x 1 x 2 { x 2 2 12 + ( x 1 + x 3 2 ) 2 } ] , $$ J=2\left[\sqrt{2}{x}_1{x}_2\left\{\frac{{x}_2^2}{12}+{\left(\frac{{x}_1+{x}_3}{2}\right)}^2\right\}\right], $$(21) σ ( x ) = 6 PL x 4 x 3 2 , $$ \sigma (x)=\frac{6{PL}}{{x}_4{x}_3^2}, $$(22) δ ( x ) = 4 P L 3 E x 3 3 x 4 , $$ \mathrm{\delta }(x)=\frac{4P{L}^3}{E{x}_3^3{x}_4}, $$(23) P c ( x ) = 4.013 E x 3 2 x 4 6 36 L 2 ( 1 - x 3 2 L E 4 G ) P = 6000 lb , L = 14 in , E = 30 e 6 psi , G = 12 e 6 psi , τ max = 13600 psi , σ max = 30000 psi , δ max = 0.25 in 0.1 x 1 2.0,0.1 x 2 10.0,0.1 x 3 10.0,0.1 x 4 2.0 $$ \begin{array}{l}{P}_c(x)=\frac{4.013E\sqrt{\frac{{x}_3^2{x}_4^6}{36}}}{{L}^2}\left(1-\frac{{x}_3}{2L}\sqrt{\frac{E}{4G}}\right)\\ P=6000{lb},L=14{in},E=30e6{psi},G=12e6{psi},\\ {\tau }_{\mathrm{max}}=13600{psi},{\sigma }_{\mathrm{max}}=30000{psi},{\delta }_{\mathrm{max}}=0.25{in}\\ 0.1\le {x}_1\le \mathrm{2.0,0.1}\le {x}_2\le \mathrm{10.0,0.1}\le {x}_3\le \mathrm{10.0,0.1}\le {x}_4\le 2.0\end{array} $$(24)

3.2 Example-2 orthogonally stiffened cylindrical shell

Objective is to optimize cost of a cylindrical shell member that is orthogonally stiffened by using ring stiffeners of box cross section and stringers of halved I-section Figures 6 and 7.

thumbnail Figure 6.

Orthogonally stiffened cylindrical shell with stringer and ring stiffener acted by compression and external pressure [21].

thumbnail Figure 7.

Halved rolled I – section as a stringer [21].

The cost function includes material cost, manufacturing cost and painting cost.

Cost function is given by [21] K = K M + i K Fi + K P , $$ K={K}_M+\sum_i{K}_{{Fi}}+{K}_P, $$(25)where KM is the material cost, KF is the fabrication cost and KP is the painting cost.

The manufacturing cost is derived from the manufacturing sequence which is explained as follows:

(KF0) Cylindrical shapes are formed from plate elements of 3 m length (Ls). K F 0 = 5 k F θ e μ . $$ {K}_{F0}=5{k}_F\theta {e}^{\mu }. $$(26)

(KF1) 3 m shell segments are welded from two curved plate elements with two butt welds using gas metal arc welding. K F 1 = 5 k F ( θ 2 ρ V 1 + ( 1.3 ) ( 0.1520 X 10 - 3 t 1.9358 ) ( 2 L s ) . $$ {K}_{F1}=5{k}_F(\theta \sqrt{2\rho {V}_1}+\left(1.3\right)\left(0.1520\mathrm{X}{10}^{-3}{t}^{1.9358}\right)\left(2{L}_s\right). $$(27)

(KF2) 15 m un-stiffened shell is welded from five shell segments with four circumferential butt welds using gas metal arc welding. K F 2 = k F ( θ 25 ρ V 1 + ( 1.3 ) ( 0.1520 X 10 - 3 t 1.9358 ) ( 8 R π ) . $$ {K}_{F2}={k}_F(\theta \sqrt{25\rho {V}_1}+\left(1.3\right)\left(0.1520\mathrm{X}{10}^{-3}{t}^{1.9358}\right)\left(8R\pi \right). $$(28)

(KF3) nr ring stiffeners are welded from three plate elements with two fillet welds each using shielded metal arc welding. K F 3 = n r k F [ 3 3 ρ V R + 1.3 X 0.3394 X 10 - 3 a wr 2 4 π ( R - h r ) ] . $$ {K}_{F3}={n}_r{k}_F\left[3\sqrt{3\rho {V}_R}+1.3X0.3394X{10}^{-3}{a}_{{wr}}^24\pi \left(R-{h}_r\right)\right]. $$(29)

(KF4) nr ring stiffeners are welded into the whole shell with 2 nr circulferential fillet welds using shielded metal arc welding. K F 4 = k F [ 3 ( n r + 1 ) ρ ( 5 V 1 + n r V R ) + 1.3 X 0.3394 X 10 - 3 a wr 2 n r 4 π R ] . $$ {K}_{F4}={k}_F\left[3\sqrt{\left({n}_r+1\right)\rho \left(5{V}_1+{n}_r{V}_R\right)}+1.3X0.3394X{10}^{-3}{a}_{{wr}}^2{n}_r4{\pi R}\right]. $$(30)

(KF5) ns stringers are welded into the shell with 2 ns fillet welds using shielded metal arc welding. K F 5 = k F [ 3 ( n r + n s + 1 ) ρ ( 5 V 1 + n r V R + n s A s L ) + 1.3 X 0.3394 X 10 - 3 a ws 2 n s 2 L ] . $$ {K}_{F5}={k}_F\left[3\sqrt{\left({n}_r+{n}_s+1\right)\rho \left(5{V}_1+{n}_r{V}_R+{n}_s{A}_sL\right)}+1.3X0.3394X{10}^{-3}{a}_{{ws}}^2{n}_s2L\right]. $$(31)

The material cost is given by K M = k M 1 5 ρ V 1 + k M 1 ρ n r V R + K M 2 ρ n s A s L . $$ {K}_M={k}_{M1}5\rho {V}_1+{k}_{M1}\rho {n}_r{V}_R+{K}_{M2}\rho {n}_s{A}_s{L}. $$(32)where V1 is the volume of the shell segment and VR is the volume of ring stiffener.

The painting cost is given by K P = k P [ 2 L + 2 ( L - n r h r ) + 2 n r π h r ( R - h r ) + 4 π R h r ( R - h r 2 ) + n s L ( h 1 + 2 b ) ] . $$ {K}_P={k}_P\left[2{R\pi L}+2{R\pi }\left(L-{n}_r{h}_r\right)+2{n}_r\pi {h}_r\left(R-{h}_r\right)+4{\pi R}{h}_r\left(R-\frac{{h}_r}{2}\right)+{n}_sL\left({h}_1+2b\right)\right]. $$(33)

Shell buckling constraint g1(X) associated with the design is given by [11] σ e = σ a 2   - σ a σ p + σ p 2   f y 1 / 1 + λ s 4 , $$ {\sigma }_e=\sqrt{{\sigma }_a^{2\enspace }-{\sigma }_a{\sigma }_p+{\sigma }_p^{2\enspace }}\le {f}_{y1}/\sqrt{1+{\lambda }_s^4}, $$(34)where σe is the equivalent stress according to Det Norske Veritas rules. And σ a = N F / 2 t e $ {\sigma }_a={N}_F/2{R\pi }{t}_e$; t e = t + A ss $ {t}_e=t+{A}_{{ss}}$; s = 2 / n s $ s=2{R\pi }/{n}_s$; σ p = p F R / t ( 1 + α ) $ {\sigma }_p={p}_FR/t(1+\alpha )$; α = A R / L e 0 t $ \alpha ={A}_R/{L}_{e0}t$; L e 0 = min ( L r , L er ) = 1.56 Rt $ {L}_{e0}=\mathrm{min}({L}_r,{L}_{{er}})=1.56\sqrt{{Rt}}$; L r = L n r - 1 $ {L}_r=\frac{L}{{n}_r}-1$; λ s 2 = f y 1 σ e ( σ a σ Eas + σ p σ Eps ) $ {\lambda }_s^2=\frac{{f}_{y1}}{{\sigma }_e}\left(\frac{{\sigma }_a}{{\sigma }_{{Eas}}}+\frac{{\sigma }_p}{{\sigma }_{{Eps}}}\right)$; σ Eas = C as π 2 E 12 ( 1 - γ 2 ) ( t s ) 2 $ {\sigma }_{{Eas}}={C}_{{as}}\frac{{\pi }^2E}{12\left(1-{\gamma }^2\right)}{\left(\frac{t}{s}\right)}^2$; C as = Ψ as 1 + ( ρ as ξ as Ψ as ) 2 $ {C}_{\mathrm{as}}={\mathrm{\Psi }}_{{as}}\sqrt{1+{\left(\frac{{\rho }_{{as}}{\xi }_{{as}}}{{\mathrm{\Psi }}_{{as}}}\right)}^2}$; Ψ as = 4 $ {\mathrm{\Psi }}_{{as}}=4$; Z as = s 2 Rt 1 - v 2 $ {Z}_{{as}}=\frac{{s}^2}{{Rt}}\sqrt{1-{v}^2}$; ξ as = 0.702 Z as $ {\xi }_{{as}}=0.702{Z}_{{as}}$; ρ as = 0.5 ( 1 + R 150 t ) - 0.5 $ {\rho }_{{as}}=0.5{\left(1+\frac{R}{150t}\right)}^{-0.5}$; σ Eps = C ps π 2 E 10.92 ( t s ) 2 $ {\sigma }_{{Eps}}={C}_{{ps}}\frac{{\pi }^2E}{10.92}{\left(\frac{t}{s}\right)}^2$; C ps = Ψ ps 1 + ( ρ ps ξ ps Ψ ps ) 2 $ {C}_{{ps}}={\mathrm{\Psi }}_{{ps}}\sqrt{1+{\left(\frac{{\rho }_{{ps}}{\xi }_{{ps}}}{{\mathrm{\Psi }}_{{ps}}}\right)}^2}$; ρ ps = 0 . $ {\rho }_{{ps}}=0.$; ξ ps = 1.04 s L r Z ps $ {\xi }_{{ps}}=1.04\frac{s}{{L}_r}\sqrt{{Z}_{{ps}}}$; Z ps = Z as $ {Z}_{{ps}}={Z}_{{as}}$; Ψ ps = [ 1 + ( s L r ) 2 ] 2 $ {\mathrm{\Psi }}_{{ps}}={\left[1+{\left(\frac{s}{{L}_r}\right)}^2\right]}^2$.

Stringer buckling constraint g2(X) is given by σ e f y 1 1 + λ p 4 , $$ {\sigma }_e\le \frac{{f}_{y1}}{\sqrt{1+{\mathrm{\lambda }}_p^4}}, $$(35)where, λ p 2 = f y 1 σ e ( σ a σ Eap + σ p σ Epp ) $ {\lambda }_p^2=\frac{{f}_{y1}}{{\sigma }_e}\left(\frac{{\sigma }_a}{{\sigma }_{{Eap}}}+\frac{{\sigma }_p}{{\sigma }_{{Epp}}}\right)$; σ Eap = C ap π 2 E 10.92 ( t L r ) 2 $ {\sigma }_{{Eap}}={C}_{{ap}}\frac{{\pi }^2E}{10.92}{\left(\frac{t}{{L}_r}\right)}^2$; C ap = Ψ ap 1 + ( ρ ap ξ ap Ψ ap ) 2 $ {C}_{{ap}}={\mathrm{\Psi }}_{{ap}}\sqrt{1+{\left(\frac{{\rho }_{{ap}}{\xi }_{{ap}}}{{\mathrm{\Psi }}_{{ap}}}\right)}^2}$; ρ ap = 0.5 $ {\rho }_{{ap}}=0.5$; ξ ap = 0.702 Z ap $ {\xi }_{{ap}}=0.702{Z}_{{ap}}$; Z ap = L r 2 Rt 0.9539 $ {Z}_{{ap}}=\frac{{L}_r^2}{{Rt}}0.9539$; Ψ ap = 1 + γ s 1 + A s s e t $ {\mathrm{\Psi }}_{{ap}}=\frac{1+{\gamma }_s}{1+\frac{{A}_s}{{s}_et}}$; γ s = 10.92 I sef s t 3 $ {\gamma }_s=10.92\frac{{I}_{{sef}}}{s{t}^3}$; s E = 1.9 t E f y $ {s}_E=1.9t\sqrt{\frac{E}{{f}_y}}$; Z G = h 1 2 t w ( h 1 4 + t 2 ) + b t f ( h 1 + t + t f 2 ) s e t + b t f + h 1 t w 2 $ {Z}_G=\frac{\frac{{h}_1}{2}{t}_w\left(\frac{{h}_1}{4}+\frac{t}{2}\right)+b{t}_f\left(\frac{{h}_1+t+{t}_f}{2}\right)}{{s}_et+b{t}_f+\frac{{h}_1{t}_w}{2}}$; if s E s ,   s e = s E $ {s}_E\le s,\enspace {s}_e={s}_E$, else s E s ,   s e = s $ {s}_E\ge s,\enspace {s}_e=s$; I sef = s e t z G 2 + ( h 1 2 ) 3 t w 12 + h 1 t w 2 ( h 1 4 + t 2 - z G ) 2 + b t f ( h 1 + t + t f 2 - z G ) 2 $ {I}_{{sef}}={s}_et{z}_G^2+{\left(\frac{{h}_1}{2}\right)}^3\frac{{t}_w}{12}+\frac{{h}_1{t}_w}{2}{\left(\frac{{h}_1}{4}+\frac{t}{2}-{z}_G\right)}^2+b{t}_f{\left(\frac{{h}_1+t+{t}_f}{2}-{z}_G\right)}^2$; A s = b t f + h 1 t w 2 $ {A}_s=b{t}_f+\frac{{h}_1{t}_w}{2}$, σ Epp = C pp π 2 E 10.92 ( t L r ) 2 $ {\sigma }_{{Epp}}={C}_{{pp}}\frac{{\pi }^2E}{10.92}{\left(\frac{t}{{L}_r}\right)}^2$, C pp = Ψ pp 1 + ( ρ pp ξ pp Ψ pp ) 2 $ {C}_{{pp}}={\mathrm{\Psi }}_{{pp}}\sqrt{1+{\left(\frac{{\rho }_{{pp}}{\xi }_{{pp}}}{{\mathrm{\Psi }}_{{pp}}}\right)}^2}$; ξ pp = 1.04 Z pp $ {\xi }_{{pp}}=1.04\sqrt{{Z}_{{pp}}}$, Z pp = Z ap $ {Z}_{{pp}}={Z}_{{ap}}$, ρ pp = 0.6 , Ψ pp = 2 ( 1 + 1 + γ s ) $ {\rho }_{{pp}}=0.6,{\mathrm{\Psi }}_{{pp}}=2\left(1+\sqrt{1+{\gamma }_s}\right)$.

Panel ring buckling constraints g3(X) and g4(X) are given by A Rreq A R   and   I Rreq I R $$ {A}_{{Rreq}}\le {A}_R\enspace \mathrm{and}\enspace {I}_{{Rreq}}\le {I}_R $$where, 1 δ r = 42 ε ,   ε = 235 f y $ \frac{1}{{\delta }_r}=42\epsilon,\enspace \epsilon =\sqrt{\frac{235}{{f}_y}}$; f y = 355 $ {f}_y=355$; δ r = 1 / 34 $ {\delta }_r=1/34$; A r = 3 h r t r = 3 δ r h r 2 ; $ {A}_r=3{h}_r{t}_r=3{\delta }_r{h}_r^2;$ A Rreq = ( 2 Z 2 + 0.06 ) L r t $ {A}_{{Rreq}}=\left(\frac{2}{{Z}^2}+0.06\right){L}_rt$; Z = L r 2 Rt 0.9539 $ Z=\frac{{L}_r^2}{{Rt}}0.9539$; L e = min ( L r , 2 × 1.56 Rt ) $ {L}_e=\mathrm{min}\left({L}_r,2\times 1.56\sqrt{{Rt}}\right)$; y E = L e t ( h r + t 2 ) + δ r h r 3 3 δ r h r 2 + L e t ;   I R = δ r h r 4 6 + 2 δ r h r 2 ( h r 2 - y E ) 2 + δ r h r 2 y E 2 + L e t ( h r + t 2 - y E ) 2 $ {y}_E=\frac{{L}_et\left({h}_r+\frac{t}{2}\right)+{\delta }_r{h}_r^3}{3{\delta }_r{h}_r^2+{L}_et};\enspace {I}_R=\frac{{\delta }_r{h}_r^4}{6}+2{\delta }_r{h}_r^2{\left(\frac{{h}_r}{2}-{y}_E\right)}^2+{\delta }_r{h}_r^2{y}_E^2+{L}_et{\left({h}_r+\frac{t}{2}-{y}_E\right)}^2$; I Rreq = I a + I p $ {I}_{{Rreq}}={I}_a+{I}_p$; I a = σ a t ( 1 + A s st ) R 0 4 500 E L R $ {I}_a=\frac{{\sigma }_at\left(1+\frac{{A}_s}{{st}}\right){R}_0^4}{500E{L}_R}$; R 0 = R - ( h r - y E ) $ {R}_0=R-\left({h}_r-{y}_E\right)$; I p = p F R R 0 2 L r 3 E [ 2 + 3 E y E δ 0 R 0 2 ( f y 2 - σ p ) ] $ {I}_p=\frac{{p}_FR{R}_0^2{L}_r}{3E}\left[2+\frac{3E{y}_E{\delta }_0}{{R}_0^2\left(\frac{{f}_y}{2}-{\sigma }_p\right)}\right]$; δ 0 = 0.005 R $ {\delta }_0=0.005R$.

Manufacturing constraint g5(X) is given by 2 ( R - h r 2 ) π n s - b 300   mm . $$ \frac{2\left(R-\frac{{h}_r}{2}\right)\pi }{{n}_s}-b\ge 300\mathrm{\enspace }\mathrm{mm}. $$(36)

Example 1 is solved by many researchers using different optimization techniques, including Ragsdell and Phillips [30] using geometric programming, Deb [10] using Genetic Algorithm, Coello [9] using GA based coevolution model, Coello and Montes [8] using flexibility based tournament selection scheme, Ray and Liew [31] using society and civilization algorithm and Qie and Ling [27] using effective co-evolutionary particle swarm optimization.

Jarmai et al. [21] solved Example 2 using PSO technique for the cost optimization considering the above design variables, objective functions, and the constraints. Value of cognitive learning coefficient and social learning coefficient was taken as 2 and 1.4 respectively. Population size was varied from 1 to 500 and maximum allowable function evaluations were 100000 with no improvement termination criterion as 10 iterations. The same problem was also solved using three different local optimization techniques viz. LFOPC, Dynamic – Q and ETOPC. The drawback of local optimization techniques is that it does not guarantee global optimum and also result changes with the change in starting point and thus they are time consuming with respect to global optimization techniques.

ABC algorithm used for the analysis was coded in MATLAB and all experiments were done on 1.4 GHz Intel Core 2 Duo Laptop machine. The performance of all the algorithms is compared for both the examples. To obtain the most suitable parameters a large number of experiments were conducted. For ABC the value of limit is set to N × S where N is the number of employed or onlooker bees and S is the dimension of design variable. One hundred trial runs were performed for each variant of ABC algorithm for each example. The performances of the different variants were compared based on different criteria. The first criterion for the comparison is the best solution obtained in one run. The algorithms are compared for the best value, mean value and the standard deviation for one hundred runs for the best solutions. In this comparison best value gives the global optima, and mean value serves the purpose for the average performance to search optimum result. Standard deviation gives the deviation of best result from the mean result for all the runs. The second criterion for the comparison is the success rate, which indicates the percentage success the technique has achieved in finding the global solution. Results of all the benchmark problems are shown in the Table 1.

Table 1.

Comparison of two different variants of ABC for different benchmark problems.

It can be seen from the Table 1 that modified ABC is better than classical ABC in finding the best solution for both constraint and unconstraint benchmark problems. For Griewank and pressure vessel problem classical ABC has shown better result than modified ABC for the mean of the beat solutions. These results do not indicate that classical ABC is inferior to modified ABC, but main intension of comparison is that modified ABC can find better solution with less population size and number of generations which indicates better convergence rate of the modified ABC. If population size and number of generation is increased than both the variants converge to the same solution.

All the variants were experimented by taking different population size and different number of generations to study its effect on best solution, average solution, success rate and processing time for the considered two examples of welded structures.

Tables 2 and 3 show the best solutions for both the examples It is observed from the result that ABC_2 has shown better results than other two variants for best, mean solutions and standard deviations. All the variants have given improvement in the result with the increase in population size and number of generations. Moreover ABC_2 require less function evaluation to reach optimum solution than ABC_1. For example 2 also similar kinds of results are observed. Moreover ABC_2 has given the optimum solution using less than 10 000 function evaluations, which is better than the results of PSO implemented by Jarmai (2006), which requires between 26 184 and 100 000 function evaluations considering a population size of 500. Results show the better performance of ABC over PSO in terms of solution as well as in terms of function evaluations

Table 2.

Comparison of different variants of ABC algorithm for Example 1 (welded beam) with different function evaluations for the best solutions obtained in 100 runs.

Table 3.

Comparison of different variants of ABC algorithm for Example 2 (Orthogonally stiffened cylindrical shell) with different function evaluations for the best solutions obtained in 100 runs.

Tables 4 and 5 show the success rate for both the examples. Here success is calculated if the results are within the tolerance of 1% of its global optima. It is observed that here also ABC_2 has better success rate with less population size and number of generations. With the increase in population size and number of generations success rate increases for all the variants. There is tremendous difference in the success rates of ABC_2 and ABC_1. It is clear from the above discussion than ABC_2 is a promising variant of ABC which serves all the requirement of a design engineer such as success rate and convergence of the solutions.

Table 4.

Comparison of different variants of ABC algorithm for Example 1 (welded beam) with different function evaluations for the success rate obtained in 100 runs.

Table 5.

Comparison of different variants of ABC algorithm for Example 2 (Orthogonally stiffened cylindrical shell)) with different function evaluations for the success rate obtained in 100 runs.

Example 1 has already been solved by many researchers and results from all the researchers are given in Table 6. It can be seen from the Table 5 that the best feasible solution found by ABC is better with respect to value of objective function than the best solutions found by other researchers using different optimization techniques. Results obtained by ABC are feasible because all the constraints are having negative values and thus satisfied.

Table 6.

Comparison of results obtained by different researchers using different optimization techniques for Example 1 (welded beam).

More top ten solutions for both the examples are given in Tables 7 and 8 obtained by using ABC to help other researchers to check other optimization algorithms.

Table 7.

Top ten optimum solution obtained for example 1 using ABC.

Table 8.

Top ten optimum solutions obtained for example 2 using ABC.

4 Conclusions

In the present work, cost optimization of welded structures has been presented. For orthogonally stiffened cylindrical shells cost function consist of material cost, manufacturing cost and painting cost. Shell buckling, stringer buckling, panel ring buckling and manufacturing limitations are considered as constraints. Classical ABC is modifies and it is applied on many constrained and unconstrained optimization benchmark problems and also to the optimization of welded structures. It is shown that ABC is capable of solving complicated cost optimization problems for welded structures. In terms of quality of solution, the results have shown that modified ABC has produced better results than classical ABC. It is also shown that ABC has given better results than other optimization techniques.

References

  1. AISC. 1970. Specification for the Design, Fabrication and Erection of Structural Steel for Buildings, Manual of Steel Construction, 7th edn. Chicago, Illinois: American Institute of Steel Construction. [Google Scholar]
  2. AISC. 1978. Manual of Steel Construction, 8th edn. Chicago, Illinois: American Institute of Steel Construction. [Google Scholar]
  3. AISC. 1986. Manual of Steel Construction – Load and Resistance Factor Design, 8th edn. Chicago, Illinois: American Institute of Steel Construction. [Google Scholar]
  4. Anderson KE, Chong KP. 1986. Least cost computer-aided design of steel girders. Engineering Journal. AISC, Fourth Quarter, 151–156. [Google Scholar]
  5. Basturk B, Karaboga D. 2006. An Artificial Bee Colony (ABC) Algorithm for Numeric Function Optimization, IEEE Swarm Intelligence Symposium, May 12–14, Indianapolis, Indiana, USA. [Google Scholar]
  6. Belegundu AD. 1982. A study of mathematical programming methods for structural optimization, Doctoral dissertation, Department of Civil and Environmental Engineering, University of Iowa, Iowa City, Iowa. [Google Scholar]
  7. Bhatti MA. 1996. Optimum cost design of partially composite steel beams using LRFD. Engineering Journal. AISC, First Quarter. p. 18–29. [Google Scholar]
  8. Coello CAC, Montes EM. 2002. Constraint-handling in genetic algorithms through the use of dominance-based tournament selection. Advanced Engineering Informatics, 16, 193–203. [Google Scholar]
  9. Coello CAC. 2000. Use of a self-adaptive penalty approach for engineering optimization problems. Computers in Industry, 41, 113–127. [Google Scholar]
  10. Deb K. 1991. Optimal design of a welded beam via genetic algorithms. AIAA Journal, 29(11), 2013–2015. [CrossRef] [Google Scholar]
  11. Det Norske Veritas (DNV). 1995. Buckling strength analysis, Classification Notes No. 30.1. Ed. Det Norske Veritas, Hovik, Norway. [Google Scholar]
  12. Dong Y, Tang J, Xu B, Wang D. 2005. An application of swarm optimization to nonlinear programming. Computer and Mathematics with Applications, 49, 1655–1668. [CrossRef] [Google Scholar]
  13. Douty R. 1980. Design of steel connections by math programming. Journal of the Structural Division, ASCE, 106(ST5), 1135–1154. [Google Scholar]
  14. Farkas J, Jarmai K. 2008. Minimum cost design of a square box column composed from orthogonally stiffened welded steel plates, International conference proceedings on design fabrication and economy of welded structures, Hungary. p. 70–71. [Google Scholar]
  15. Farkas J. 1991. Fabrication aspects in the optimum design of welded structures. Structural Optimization, 3(1), 51–58. [CrossRef] [Google Scholar]
  16. Farkas J, Jarmai K. 1994. Minimum cost design of laterally loaded welded rectangular cellular plates. Structural Optimization, 8(4), 262–267. [CrossRef] [Google Scholar]
  17. Farkas J, Jarmai K, Snyman J. 2009. Global minimum cost design of a welded square stiffened plate supported at four corners, Structural and Multidisciplinary Optimization, 40, 477–489. [CrossRef] [Google Scholar]
  18. Farkas J. 2002. Thickness design of axially compressed unstiffened cylindrical shells with circumferential welds. Welding in the World, 46, 26–29. [CrossRef] [Google Scholar]
  19. Hojjat A, Kamal C. 2006. Cost Optimization of Structures, Fuzzy Logic, Genetic Algorithms and Parallel Computing. John Wiley: UK. [Google Scholar]
  20. Jarmai K. 2007. Single- and multiobjective optimization of a welded stringer-stiffened cylindrical shell. Periodica Polytechnica, 51(2), 11–18. [Google Scholar]
  21. Jarmai K, Snyman JA, Farkas J. 2006. Minimum cost design of a welded orthogonally stiffened cylindrical shell. Computers and Structures, 84, 787–797. [CrossRef] [Google Scholar]
  22. Karaboga D. 2005. An idea based on honey bee swarm for numerical optimization, Technical Report-TR06, Erciyes University, Engineering Faculty, Computer Engineering Department. [Google Scholar]
  23. Karaboga D, Basturk B. 2007. Artificial bee colony (abc) optimization algorithm for solving constrained optimization problems, in IFSA 2007, LNAI 4529, Melin P et al., Editors. Springer-Verlag: Berlin, Heidelberg. p. 789–798. [Google Scholar]
  24. Karaboga D, Basturk B. 2008. On the performance of artificial bee colony (ABC) algorithm. Applied Soft Computing, 8, 687–697. [CrossRef] [Google Scholar]
  25. Lorenz RF. 1988. Understanding composite beam design methods using LRFD. Engineering Journal. AISC, First Quarter. p. 35–38. [Google Scholar]
  26. Luis MCS, Farkas J, Jarmai K. 2008. Optimization of an orthogonally stiffened plate considering fatigue constraints, International conference proceedings on design fabrication and economy of welded structures, Hungary. p. 27–34. [Google Scholar]
  27. Marcin C, Farkas J, Andrej G, Jarmai K, Katarzyna R. 2008. Optimization of steel beams and columns for variable rib configuration, International conference proceedings on design fabrication and economy of welded structures, Hungary. p. 113–118. [Google Scholar]
  28. Montes EM, Coello CAC. 2005. Useful Infeasible Solutions in Engineering Optimization with Evolutionary Algorithms. LNAI, 3789, 652–662. [Google Scholar]
  29. Qie H, Ling W. 2007. An effective co-evolutionary particle swarm optimization for constrained engineering design problems. Engineering Applications of Artificial Intelligence, 20, 89–99. [CrossRef] [Google Scholar]
  30. Ragsdell KM, Phillips DT. 1976. Optimal design of a class of welded structures using geometric programming. ASME Journal of Engineering for Industries, 98(3), 1021–1025. [CrossRef] [Google Scholar]
  31. Ray T, Liew K. 2003. Society and civilization: an optimization algorithm based on the simulation of social behavior. IEEE Transactions on Evolutionary Computation, 7, 386–396. [CrossRef] [Google Scholar]
  32. Sandgren E. 1988. Nonlinear integer and discrete programming in mechanical design, Proceedings of the ASME Design Technology Conference, Kissimine, FL. p. 95–105. [Google Scholar]

Cite this article as: Savsani V: Implementation of modified artificial bee colony (ABC) optimization technique for minimum cost design of welded structures. Int. J. Simul. Multisci. Des. Optim., 2014, 5, A11.

All Tables

Table 1.

Comparison of two different variants of ABC for different benchmark problems.

Table 2.

Comparison of different variants of ABC algorithm for Example 1 (welded beam) with different function evaluations for the best solutions obtained in 100 runs.

Table 3.

Comparison of different variants of ABC algorithm for Example 2 (Orthogonally stiffened cylindrical shell) with different function evaluations for the best solutions obtained in 100 runs.

Table 4.

Comparison of different variants of ABC algorithm for Example 1 (welded beam) with different function evaluations for the success rate obtained in 100 runs.

Table 5.

Comparison of different variants of ABC algorithm for Example 2 (Orthogonally stiffened cylindrical shell)) with different function evaluations for the success rate obtained in 100 runs.

Table 6.

Comparison of results obtained by different researchers using different optimization techniques for Example 1 (welded beam).

Table 7.

Top ten optimum solution obtained for example 1 using ABC.

Table 8.

Top ten optimum solutions obtained for example 2 using ABC.

All Figures

thumbnail Figure 1.

Normal distribution of 1e7 random numbers with mean 0 and variance 1.

In the text
thumbnail Figure 2.

Uniformly distributed random number between –1 and 1 for 1e7 random points.

In the text
thumbnail Figure 3.

Variation of solution for the sphere function for classical ABC for five different runs.

In the text
thumbnail Figure 4.

Variation of solution for the sphere function for modified ABC for five different runs.

In the text
thumbnail Figure 5.

Welded beam design (Example 1).

In the text
thumbnail Figure 6.

Orthogonally stiffened cylindrical shell with stringer and ring stiffener acted by compression and external pressure [21].

In the text
thumbnail Figure 7.

Halved rolled I – section as a stringer [21].

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.