Reliability-Based Design Optimization for Crane Metallic Structure Using ACO and AFOSM Based on China Standards

The design optimization of crane metallic structures is of great significance in reducing their weight and cost. Although it is known that uncertainties in the loads, geometry, dimensions, and materials of crane metallic structures are inherent and inevitable and that deterministic structural optimization can lead to an unreliable structure in practical applications, little amount of research on these factors has been reported. This paper considers a sensitivity analysis of uncertain variables and constructs a reliability-based design optimization model of an overhead traveling crane metallic structure. An advanced first-order second-moment method is used to calculate the reliability indices of probabilistic constraints at each design point. An effective ant colony optimization with a mutation local search is developed to achieve the global optimal solution. By applying our reliability-based design optimization to a realistic crane structure, we demonstrate that, compared with the practical design and the deterministic design optimization, the proposed method could find the lighter structure weight while satisfying the deterministic and probabilistic stress, deflection, and stiffness constraints and is therefore both feasible and effective.


Introduction
As tools for moving and transporting goods, cranes are used in various settings to aid the development of economy and undertake the heavy logistical handling tasks in factories, railways, ports, and so on.Metallic structure is mainly made out of rolled merchant steel and plate steel by welding method according to certain structural organization rules.Crane metallic structures (CMSs), also called crane bridge, bear and transfer the burden of crane and their own weights.CMSs are mechanical skeleton and form the main components of cranes.Their design qualities have direct impact on the technical and economic benefit, as well as the safety, of the whole crane.
Generally, a CMS requires a large quantity of steel and consequently weighs a considerable amount.Its cost accounts for over a third of the total cost of the crane.Thus, under the condition of satisfying the relevant design codes, improving the performances of the CMS, saving material, and reducing weight have important significance in terms of cost-savings.As a consequence, various scholars have researched on this problem, and current optimization methods mainly include finite element method [1], neural networks [2], and Lagrange multipliers [3], amongst others [4][5][6].However, these methods are all based on deterministic optimum designs and do not consider the effect of randomness in the structure and/or load.Studies have shown that the loads effected on the CMS, and the materials and geometric dimensions of the structure itself are uncertain.Deterministic optimum designs are pushed to the limits of their constraints boundaries, leaving no room for uncertainty.Optimum designs obtained without consideration of such uncertainties are therefore unreliable [7,8].Recently there have been a few reports about reliability-based design of crane structure [9,10].Literature [9] researched on the reliability-based design of structure of tower crane based on the finite element analysis (FEA) and response surface method (RSM).Meng et al. [10] analyzed the reliability and sensitivity of crane metal structure by means of BP neural network based on finite element and first-order second-moment (FOSM) method.Nevertheless, design only considering reliability could not guarantee the best work performance and the optimal design parameters.
The aim of a design is to achieve adequate safety at minimum cost under the condition of meeting specified performance requirements.Hence, the optimization based on reliability concepts appears to be a more rational design philosophy, which is why reliability-based design optimization (RBDO) has been developed.RBDO incorporates the optimization of design parameters and reliability calculations for specified limit states.At present, it is attracting increased attention, both in theoretical research and practical applications [11][12][13].Despite advances in this area, few RBDO approaches specific to CMS have appeared in the technical literature.Therefore, this paper develops an RBDO methodology for optimizing CMS that both minimizes the weight and guarantees structural reliability.The main structural behaviors are modeled by the crane design code (China Standard) [14] based on material mechanics, structural mechanics, and elasticity theory.
CMSs are engaged in busy and heavy work.It must have sufficient strength, stiffness, and stability under complex operating conditions.Their design calculations involve the hyperstatic problem of space structures.Therefore, both the calculating model and design calculation are very complicated.Furthermore, in practical production, structural dimensions are usually taken for integral multiples of millimeters and the specified thickness of the steel plates [15].Due to these manufacturing limitations the design variables cannot be considered as continuous but should be treated as discrete in a large number of practical design situations, which means that the CMS design optimization is a constrained nonlinear optimization problem with discrete variables, known as NP-complete combinatorial optimization.To solve such problems, recent studies have focused on the development of heuristic optimization techniques, such as genetic algorithms (GAs) [16,17], particle swarm optimization (PSO) [18][19][20], ant colony optimization (ACO) [21,22], big bang-big crunch (BB-BC) [23,24], imperialist competitive algorithm (ICA) [25], and charged system search (CSS) [26].These algorithms can overcome most of the limitations found in traditional methods, such as becoming trapped in local minima and impractical computational complexity [27,28].In view of the simple operation, easy implementation, and the suitability of ACO for computational problems involving discrete variables and combinatorial optimization, the optimization process of BRDO described in this paper is performed using ACO with a mutation local search (ACOM) [29].
Structural reliability can be analyzed using analytical methods, such as first-and second-order second moments (FOSM and SOSM) [30] and advanced first-order second moment (AFOSM) [31], or with simulation methods such as Monte Carlo sampling (MCS) method.The FOSM method is very simple and requires minimal computation effort but sacrifices accuracy for nonlinear limit state functions.The accuracy of the SOSM method is improved compared with that of the FOSM, but its computation effort is greatly increased and this makes it not frequently used in practices.MCS method is accurate; however, it is computationally intensive as it needs a large number of samples to evaluate small failure probabilities.The AFOSM method, a more accurate analytical approach than the FOSM method, is able to efficiently handle low-dimensional uncertainties and nonlinear limit state functions [32] and is applied in most practical cases.It is used to calculate the reliability indices of RBDO in this paper.
The paper is organized as follows.Section 2 outlines the general formulation of discrete RBDO and then Section 3 constructs the RBDO model of an overhead traveling CMS (OTCMS).Section 4 develops the ACOM algorithm used for the optimization process of the RBDO and Section 5 describes the AFOSM method applied for reliability analysis.The RBDO procedure is illustrated in Section 6.Some applied examples that demonstrate the potential of the proposed approach for solving realistic problem are presented in Section 7, followed by concluding remarks in Section 8.

Formulation of Discrete RBDO
In contrast to deterministic design optimization (DDO), RBDO assumes that quantities related to size, materials, and applied loads of a structure have a random nature to conform to the actual one.The parameters characterizing these quantities are called random variables, and these need to be taken into account in reliability analysis.These random variables may be either random design variables or random parameter variables.In optimization process, the mean values of the random design variables are treated as optimization variables.The formulation of discrete RBDO problem is generally written as follows: where

RBDO Modeling of an OTCMS
Cranes are mechanically applied to moving loads without interfering in activities on the ground.As overhead traveling cranes are the most widely used, a typical OTCMS is selected as the study object.As shown in Figure 1, it is composed of two parallel main girders that span the width of the bay between the runway girders.The OTCMS moves longitudinally, and the two end carriages located on either sides of the span house the wheel blocks.Because the main girders are the principal horizontal beams that support the trolley and are supported by the end carriages, they are the primary load-carrying components and account for more than 80% of the total weight of the OTCMS.Therefore, the RBDO of OTCMS mainly focuses on the design of these main girders.The crane's solid-web girder is usually a box section fabricated from steel plate, for the main and vice webs, top, and bottom flanges, as shown in Figure 2. Therefore, given a span, its RBDO is to obtain the minimum dead weight, that is, the minimum main girder cross-sectional area, which simultaneously satisfies the required deterministic and probabilistic constraints associated with strength, stiffness, stability, manufacturing process, and dimensional limits [14].In this paper, a practical bias-rail box main girder is considered.The mathematical model of its RBDO is given in Table 1.
A calculation diagram of section 1-1 of the main girder is shown in Figure 2. Figure 3 is a force diagram for this main girder in the vertical and horizontal planes.In Table 1 and Figures 2 and 3,   represents the uniform load;   ( = 1, 2, . ..) denotes the concentrated load; ∑  stands for the trolley wheel load applied by the sum of the lifting capacity   and trolley weight   ;   and   represent the horizontal concentrated and uniform inertial load, respectively;   is the lateral force.
In the RBDO model, the section dimensions are random variables, while their mean values are treated as design variables (see Table 1).The independent uncertain variables include section dimensions and some important parameters (Table 1), and the other parameters are considered to be deterministic.It is assumed that all random variables are normally distributed around their mean value.Effects of each random variable on the active constraints at the optimal point are illustrated in Figure 4.The inequality constraints of the RBDO problem define, in turn, the failure by stress, deflection, fatigue, and stiffness.These structural behaviors are described by the assumed limit state functions given in Table 1 [15].With respect to estimating the most critical configuration of the trolley on the main girder, there are three configurations.(i) The limit state functions   (X, Y),  = 1, 2, 3 correspond to the position of the maximum bending moment, when a fully loaded trolley is lowering and braking

Variable and parameter
is taken equal to 20 mm which facilitates welding and   is set equal to 15 so as to guarantee the local stability of flange overhanging and the rail installation.
Maximum shear stress at the middle of main web (2-2 section).
Vertical static deflection of the midspan Horizontal static displacement of the midspan where   in the midspan and at the same time the crane bridge is starting or braking.(ii) The limit state functions   (X, Y),  = 4, 5 correspond to the position of the maximum shear stress, when a fully loaded trolley is lowering and braking at the end of the span and at the same time the crane bridge is starting or braking.(iii) The limit state functions   (X, Y),  = 9, 10, 11, 12 correspond to the position of maximum deflection.This is when a fully loaded trolley is positioned in the midspan.The fatigue limit state functions   (X, Y),  = 7, 8 correspond to the normal working condition of the crane.The local stability (local buckling) of the main girder can be guaranteed by arranging transverse and longitudinal stiffeners to form the grids according to the width-to-thickness ratio of the flange and the heightto-thickness ratio of the web.Thus, only the global stability of the main girder is considered here.Therefore, the RBDO model of the OTCMS main girder is a five-dimensional optimization problem with 11 deterministic and 11 probabilistic constraints.

The ACO for the OTCMS
ACO simulates the behavior of real life ant colonies, in which individual ants deposit pheromone along a path while moving from the nest to food sources and vice versa.Thereby, the pheromone trail enables individual to smell and select the optimal routs.The paths with more pheromone are more likely to be selected by other ants, bringing on further amplification of the current pheromone trails and producing a positive feedback process.This behavior forms the shortest path from the nest to the food source and vice versa.The first ACO algorithm, called ant system (AS), was applied to solve the traveling salesman problem.Because the search parallelism of ACO is based on the components of a solution level, it is very efficient.Thus, since the introduction of AS, the ACO metaheuristic has been widely used in many fields [33,34], including for structural optimization, and has shown promising results for various applications.Therefore, ACO is selected as optimization algorithm in the present study.( The initial solution is randomly selected.We set the initial pheromone level [  (  )] 0 of all array elements in the space to be zero.The heuristic information   (  ) of array element array  [  ] is expressed by (4), so as to induce subsequent solutions to select smaller variable values as far as possible and accelerate optimization process:

Selection Probability and Construction of Solutions.
As the ants move from node to node to generate paths, they will ceaselessly select the next node from the unvisited neighbor nodes.This process forms the ant paths and, thus, in our algorithm, constructs solutions.In accordance with the transition rule of ant colony system (ACS) [35], each ant begins with the first array element array  GAs are a powerful tool for solving combinatorial optimization problems.They solve optimization problems using the idea of Darwinian evolution.Basic evolution operations, including crossover, mutation, and selection, make GAs appropriate for performing search.In this paper, the mutation operation is introduced to the proposed algorithm to perform local searches.We assume that the current global optimal solution d best = {array

Pheromone Updating.
The pheromone updating rules of ACO include global updating and local rules.When ant  has finished a path, the pheromone trails on the array elements through which the ant has passed are updated.In this process, the pheromone on the visited array elements is considered to have evaporated, thus increasing the probability that following ants will traverse the other array elements.This process is performed after each ant has found a path; it is a local pheromone update rule with the aim of obtaining more dispersed solutions.The local pheromone update rule is When all of the ants have completed their paths (which is called a cycle), a global pheromone update is applied to the array elements passed through by all ants.This process is applied in an iterative mode [29].The rule is described as follows: where (d   , P, X   , Y) represents the fitness function value of ant  on the th iteration (see the next section).array  [  ] belongs to the array elements of the path generated by ant  on the th iteration. is a phase constant ( 1 ≥  3 ≥  2 ) depending on the quality of the solutions to reinforce the pheromone of the best path and evaporate that of the worst. ∈ [0, 1] and  ∈ [0, 1] are the local and global pheromone evaporation rates, respectively.Other parameters are the same as before.

Evaluation of Solution.
The aim of OTCMS RBDO is to develop a design that minimizes the total structure weight while satisfying all deterministic and probabilistic constraints.The ACO algorithm was originally developed for unconstrained optimization problems, and hence it is necessary to somehow incorporate constraints into the ACO algorithm.Constraint-handling techniques have been explored by a number of researchers [36], and commonly employed methods are penalty functions, separation of objectives and constraints, and hybrid methods.Penalty functions are easy to implement and, in particular, are suitable for discrete RBDO.Hence, the penalty function method is selected for constraint handling.The following fitness function is used to transform a constrained RBDO problem to an unconstrained one: where (d,P,X,Y) is unconstrained objective function (the fitness function); (d,P) is the original constraint objective function (see ( 1)); , , and  are positive problem-specific constants; and   (d,P),   (X,Y) are penalty functions corresponding to the th deterministic and th probabilistic, respectively.When satisfied, these penalty functions return to a value of zero; otherwise, the values would be amplified according to the square term in (8).Other parameters are the same as before.

Termination Criterion.
Each run is allowed to continue for a maximum of 100 generations.However, a run may be terminated before this when no improvement in the best objective value is noticed.

Structural Reliability Analysis
The goal of RBDO is to find the optimal values for a design vector to achieve the target reliability level.In accordance with the RBDO model of OTCMS, randomness in the structure is expressed as the random design vector X and the random parameter vector Y.The limit state functions   (X,Y), which can represent the stress, displacement, stiffness, and so on, is defined in terms of random vectors S (define S ∈ (X,Y)).The limit states that separate the design space into "failure" and "safe" regions are   (X,Y) = 0. Accordingly, the probability of structural reliability with respect to the th limit state function in the specified mode is where  denotes the safe domain (  (S) ≤ 0).  ( 1 ,  2 , . . .,   ) is the joint probability density function (PDF) of the random vector S, and   can be calculated by integrating the PDF   ( 1 ,  2 , . . .,   ) over .Nevertheless, this integral is not a straightforward task, as   ( 1 ,  2 , . . .,   ) is not always available.To avoid this calculation, moment methods and simulation techniques can be applied to estimate the probabilistic constraints.FOSM method is broadly used for RBDO applications owing to its effectiveness, efficiency, and simplicity and it was recommended by the Joint Committee of Structural Safety.It solves structural reliability using mean value and standard derivation.At first, the performance function is expanded using the Taylor series at some point; then truncating the series to linear terms, the first-order approximate mean value and standard deviation may be obtained and the reliability index could be solved.Therefore, it is called FOSM.According to the difference of the selected linearization point, FOSM is divided into mean value first-order second moment (MFOSM) (the linearization point is mean value point) and advanced first-order second moment (AFOSM) also called Hasofer-Lind and Rackwitz-Fiessler method (the linearization point is the most probable failure point (MPP)).The advantages of AFOSM are that it is invariant with respect to different failure surface formulations in spaces having the same dimension and more accurate compared with FOSM [37,38].Consequently, the MPP-based AFOSM is used to quantify probabilistic characterization in this research.
AFOSM uses the closest point on the limit state surface to the origin in the standard normal space as a measure of the reliability.The point is called as the design point or MPP S * , and the reliability index  is defined as the distance of the design point and the origin,  = ‖S * ‖, which could be calculated by determining the MPP in random variable space.
Firstly, obtain a linear approximation of the performance function   =   (S) by using the first-order Taylor's series expansion about the MPP S * : Since the MPP S * is on the limit state surface, the limit state function equals zero ( * 1 ,  * 2 , . . .,  *  ) = 0.Here the subscript  has been dropped for the sake of simplicity of the subsequent notation.
's mean value   and standard deviation   could be expressed as follows: The reliability index  is shown as where Then Finally, combining ( 13), (15), and the limit state function (S) = 0, the MPP  *  and reliability index  could be calculated by an iterative procedure.Then, the reliability  could be approximated by  = Φ() where Φ(⋅) is the standard normal cumulative distribution function.Here, random variables   ( = 1, 2, . . ., ) are assumed to be normal distribution and are independent to each other and the same assumption will be used throughout this paper.

RBDO Procedure
Using the methods described in the previous sections, the RBDO numerical procedure illustrated in Figure 5 was developed.

Example of Reliability-Based Design Optimization
The proposed approach was coded in C++ and executed on 2.26 GHz Intel Dual Core processor and 1 GB main memory.(units in mm).In practical design, the main girder height and width are usually designed as integer multiples of 5 mm, so the number of discrete mean values available for the design variables (main girder height and width) was  1 = 65 and  5 = 100.The step size intervals were set to be 0.5 mm for web and flange thicknesses of less than 30 mm and to be 1 mm for thicknesses more than 30 mm [15].Thus, the number of discrete mean values available for the design variables (thicknesses of main web, vice web, and flange) was  2 =  3 =  4 = 60.The statistical properties of the random variables are summarized in Table 3.For the desired reliability probabilities, we refer to JCSS [39] and the current crane design code [14]; a value of 0.979 was set for

RBDO Results
. Some final points concerning the practical design and the deterministic and reliability-based optimization process are given in Table 5.Table 6 shows the performance and reliability with respect to each optimum design.Variations in the number of iterations required for the reliability indices of active constraints and cross-sectional area to converge are illustrated in Figure 6.
Tables 5 and 6 show that the deterministic design took 58 generations to find the optimum area of 30875 mm 2 , which is about 8825 mm 2 less than that of the practical design.Thus, the deterministic optimization found a better solution within just 0.1794 s.The critical constraint of the optimum design is the vertical natural vibration frequency   at the midspan point.The corresponding constraint value is 2.00067 Hz,  which is higher than the required value 2 Hz, and satisfies the performance requirement.Probabilistic analyses were also conducted for this deterministic optimum and the practical design; the results are shown in Table 6.The active constraints of the deterministic optimum have reliabilities of 0.604, 0.728, and 0.494, which are all below the desired reliabilities of 0.979, 0.968, and 0.759, respectively.The results indicate that the DDO can significantly reduce the structural area, but its ability to meet the design requirements for reliability under uncertainties is quite low.To obtain a more reliable design by considering uncertainties during the optimization process, RBDO is needed.
As shown in Tables 5 and 6, RBDO required 23 optimization iterations.The reliabilities of the active constraints at the optimum point are 0.996, 0.989, and 0.851 (corresponding to reliability indices  3 = 2.67,  8 = 2.32, and  10 = 1.04), which are all above the desired reliabilities of 0.979, 0.968, and 0.759.Compared with DDO, the final area given by the RBDO process increases from 30875 to 35756 mm 2 (an increase of 15.9%), and the CPU time is one order of magnitude more than that of DDO.However, the reliability of the RBDO results exhibits a significant increase and meets the desired levels.Therefore, considering the inherent uncertainties in material, dimensions, and loads, only the final RBDO design is both feasible and safe.

Conclusions
This paper presented an RBDO methodology that combines ACOM and AFOSM.This was applied to the design of a real-world OTCMS under uncertainties in loads, crosssectional dimensions, and materials for the first time.The design procedure directly couples structural performance calculation, numerical design optimization, and structural reliability analysis while considering different modes of failure in the OTCMS.From the results obtained, the following conclusions could be drawn.The deterministic optimization method can improve design quality and efficiency; nevertheless, it is more likely to lead to unreliable solutions once we consider uncertainty.On the contrary, RBDO can achieve a more compromised design that balances economic and safety.The inherent nature of uncertain factors in the design of CMSs means that RBDO is a more realistic design method.It is worth noting that, in such high-risk equipment, an increase in the reliability that leads to a cost decrement is financially much more beneficial rather than increasing the weight which results in the cost increments on a long view.The constructed approach is applicable and efficient for OTCMSs RBDO and might also be useful for other metallic structures with more design and random variables, as well as multiple objectives.This will be studied in a future work.

Figure 1 :
Figure 1: Metal structure for overhead travelling crane.

Figure 4 :
Figure 4: Effects of the random variables on the active constraints.

Table 1 :
Mathematical model of the RBDO for OTCMS.
Design variable mean values of the optimization problem d = [ Mathematical Problems in Engineering 4.1.Representation and Initialization of Solution.In our ACO algorithm, each solution is composed of different array elements that correspond to different design variables.One ant's search path represents a solution to the optimization problem or a set of design schemes.The path of the th ant in the -dimensional search space at iteration  could be denoted by d   ( = 1, 2, . . ., popsize; popsize denotes the population size).Continuous array elements array  [  ] are used to store the mean values of the discrete design variable   in the nondecreasing order (  denotes the array sequence number for the th design variable mean value   .This is integer with 1 ≤   ≤   . = 1, 2, . . ., , where  is the number of design variables and   is the number of discrete values available for   . lower the lower and upper bounds of   ).Thus, consider d   = {array 1 [ 1 ] , array 2 [ 2 ] , . . ., array  [  ]} are = [ 1 ,  2 , . . .,   ]   .
)] , if array  [  ] ∈   , 0 otherwise, (5b) where   is the set of feasible neighbor array elements of array  [  ] ( = +1 and  ≤ ).  [  ] and   (  ) denote the pheromone 1 [ 1 ], array 2 [ 2 ], . . ., array  [  ]} best has not been improved for a certain number of stagnation generations, sgen.One or more array elements are chosen at random from d best , and these are changed in a certain manner.Through this mutation operation, we obtain the mutated solution d  best .If d  best is better than d best , we replace d best with d  best .Otherwise, the global optimal solution remains unchanged.

Table 2 :
Main technical characteristics and mechanical properties of the metallic structure.

Table 3 :
Statistical properties of the random variables   .

Table 4 :
Parameter values used in adopted ACO.
1 (X, Y) to  7 (X, Y), 0.968 for  8 (X, Y) and  9 (X, Y), and 0.759 for  10 (X, Y) and  11 (X, Y).It should be noted that these target reliabilities serve only as examples to illustrate the proposed RBDO approach and are not recommended design values.By means of a large number of trials and experience, the parameter values of the constructed ACO algorithm were set in Table4.
i Figure 5: Numerical procedure of RBDO.

Table 6 :
The performance and reliability with respect to optimum designs.