Structural nonlinear vibration reliability analysis based on improved ant colony algorithm

Yancang Li1 , Chenguang Ban2 , Shujing Zhou3 , Qiuyue Cheng4 , Shuanghong Peng5

1, 2, 3, 4, 5College of Civil Engineering, Hebei University of Engineering, Handan, 056038, Hebei, China

1Corresponding author

Journal of Vibroengineering, Vol. 16, Issue 3, 2014, p. 1231-1239.
Received 31 October 2013; received in revised form 1 November 2013; accepted 3 November 2013; published 15 May 2014

Copyright © 2014 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License

To improve the efficiency of the structural nonlinear vibration reliability optimization, the solution of the structural reliability index is converted into the corresponding unconstrained optimization problem by combining with the penalty function method through the analysis of the geometric meaning. And a novel ant colony algorithm is employed to solve the optimization problem. The introduction of ants transfer strategy makes the algorithm expand the search range in the early phase, and improve the convergence rate in the late which can overcome the defects of precocity and stagnation. Finally, the improved ant colony algorithm was used to find the shortest failure paths of the structure. Engineering practice and comparison with other algorithms demonstrate the strong applicability and validity of the method.

Keywords: nonlinear vibration, reliability index, failure paths, ant colony algorithms, penalty function method.

1. Introduction

The application of reliability in structural design began in the 1940s. The calculation of structural reliability is the core content of structural design, and many approximate calculation methods have been proposed in recent years [1-6]. Calculation methods of reliability indexes commonly used in engineering include First Order Second Moment method, JC method, Gradient Optimization method, Response Surface method, Monte Carlo method and Stochastic Finite Element method, which play an important role in the reliability analysis of engineering structures. But in the actual engineering reliability analysis, the use of these methods often encounters many difficult problems, such as the limit state surface is difficult to parse and describe, effective search of solution space is difficult to achieve for finding out the shortest distance from standard normal space origin to limit state surface, and the weakest failure paths. In recent years, a number of new intelligent algorithms gradually have risen and aroused widespread interest of the majority of scholars, and ant colony optimization algorithms is one of them. The algorithm makes comprehensive use of positive feedback, distributed computing, and the greedy inspired search technology, which can be used to solve the general forms of the non-convex and nonlinear constrained optimization problems, and shows unique advantages in solving the combinatorial optimization problems. Therefore, the ant colony algorithm was employed to the structure nonlinear vibration reliability analysis.

Due to random factors, whether predetermined function of the structure can be completed in the future is unable to be prior explicitly judged, while only the completion possibility level of its predetermined function can be determined. The higher possibility level leads to more reliable structure, and vice versa. Therefore, the reliability of structure is measured by probability both at home and abroad, which is defined as the probability to complete the predetermined function within a predetermined time and a predetermined condition [7-10]. The reliability of structure is gradually developed in engineering practice based on reorganization of uncertain factors, such as loading and materials, which is indispensable in structural design.

According to this, we proposed an improved ant colony optimization algorithm and employed it to the solution of the failure paths searching and the reliability index obtaining. This study can provide an effective method for the structure nonlinear vibration reliability analysis.

The rest of the paper is organized as follows. First, the improved algorithm was proposed after the introduction of the basic knowledge of ant colony optimization algorithm. Then, it was employed to the engineering practice. Finally, the efficiency of the improved algorithm was demonstrated through the engineering practice and comparison with other algorithms.

2. Improved ant colony algorithm

2.1. Summary of basic ACO

Ant colony algorithm, proposed by M. Dorigo, V. Maniezzo and A. Colorni in the early 1990s, is a group of the intelligent optimization algorithms, which made significant effects in solving the TSP problem with the characteristics of the NP-hard in early stage, where traditional optimization methods were difficult to work, thereafter the algorithm is widely concerned by academia and industry [11-13]. The ant colony algorithm does not involve the partial derivatives of the objective function, and is able to overcome shortcoming of the other algorithm to fall into local optimum, thus possessing strong adaptability. Currently, the ant colony algorithm has been successfully applied to structural optimization, communication networks, vehicle scheduling, construction engineering and many other fields [14-19].

Here an example of the classic TSP problem is used to elaborate the ant colony system model. For other problems, the model can be used by a little change  m=i=1nbi(t), where m is the number of ants in the colony, bi(t) accounts for the number of ants in the city i at time t. dij(i,j=1,2,,n) is the distance between city i and city j, τij(t) accounts for residual amount of information on the j path at time t. In the initial time, the pheromone on each path is equal and τij(0) is assumed as a constant C. The ant k(k=1,2,,m) decides to shift the direction according to the pheromone on the path, and pijk(t) accounts for the probability of ant k transfers from the city i to the city j at time t:

p i j k = τ i j α ( t ) η i j β ( t ) s a l l o w e d k τ i s α ( t ) η i s β ( t ) , j a l l o w e d k , 0 , otherwise ,

where allowedk={0,1,,n-1}-tabuk is the cities that the ant k is allowed to select in the next step, and tabuk(k=1,2,,m) is used to record the cities the ants have traversed, which will be dynamically adjusted with the evolution process. Through the n moments, the ants complete a path traversal, and the pheromone on each path will be adjusted according to the formula:

τ i j t + n = ρ τ i j t + Δ τ i j ,
Δ τ i j = k = 1 m Δ τ i j k ,

where ρ is the pheromone evaporation coefficient, Δτijk is the pheromone left on the ij path by the kth ant in the cycle, and Δτij accounts for the increment of the pheromone on the ij path in the cycle.

2.2. Improvement of ACO

In order to overcome the deficiencies of the ant colony algorithm, a dynamic combination of the pheromone and heuristic information is considered here to control the weights of respective share by time, so as to achieve the purpose of broadening the search scope in the prophase, and improving the convergence rate in the late stage.

The improved selection strategy of ants is following:

p i j k ( t ) = 1 - ω t τ i j α t k t a b l e k τ i s α t + ω t η i j β t k t a b l e k η i s β t , j t a b l e k , 0 otherwise ,

where ωt is the dynamic weight, whose range is [0,1]. ω(t) will be adjusted and changed in real time with the search process of ants, so that the ants could maintain a balance between the “exploration” and “utilization”. The following gradient functions can be chosen:

ω ( t ) = w 1 , t T 1 , w 2 , T 1 < t T 2 , w 3 , T 2 < t ,

where wi(i=1, 2, 3) is a constant, and is different values for the corresponding step function. If the optimal solution obtained in a period of time does not change, the added pheromone should be reduced, that is, the value of ω(t) is raised to make it escape from the local minimum point. To avoid falling into local optimal solution in the early search stage, an appropriate amount of negative feedback is added into the search process and a less value of ω(t) is chosen to expand the search range of the algorithm. Fig. 1 shows the flowchart of the improved ant colony algorithm.

Fig. 1. The framework of the improved ant colony algorithm

The framework of the improved ant colony algorithm

2.3. Simulation analysis of algorithm

We verify the better search capability and stability of the improved ant colony algorithm than that of the basic ant colony algorithm by simulation experiments. Both algorithms are run for 10 times, and each time the maximum number of iterations is 200. Table 1 shows the parameter settings of the basic ant colony algorithm and the improved ant colony algorithm. The experiments solve the problems for Oliver30 and d1291 by the improved ant colony algorithm and the basic ant colony algorithm, and the simulation results are shown in Fig. 2 and Table 2.

Fig. 2. Evolutionary graph on problem Oliver30 and d1291

 Evolutionary graph on problem Oliver30 and d1291

Table 1. Parameter settings of the experiments

The basic ant colony algorithm
The improved ant colony algorithm
α 1
β 1
ρ 1
Q 1
α 2
β 2
ρ 2
Q 2
ω 1
ω 2
ω 3
T 1
T 2

Table 2. Comparative table of the simulation results

TSP problem
The basic ant colony algorithm
The improved ant colony algorithm
Optimal solution
Average solution
Average evolution algebra
Optimal solution
Average solution
Average evolution algebra

According to the experimental results, compared with the basic ant colony algorithm, the optimal solution and the average evolution algebra of the improved ant colony algorithm have been increased, indicating that the improved algorithm can search for the better solution, and global convergence has been enhanced. Moreover, the average solution of the improved ant colony algorithm here is greater than that of the basic ant colony algorithm, indicating that stability of the improved algorithm has also been enhanced.

3. Engineering practice

Here, we conduct analysis of nonlinear vibration and solving of reliability index on the simply supported concrete beam with an edge crack. The geometry size and material constant of the beam are: E=200 Gpa, ρ=7850 kg/m3,h=10 mm, b=20 mm and L=300 mm. The schematic diagram of the beam with the edge crack is shown in Fig. 3.

Fig. 3. Model of a cracked beam

 Model of a cracked beam

As shown in Fig. 3, the crack beam is under the action of bending moment M, the crack location is x0, its geometry and material characteristic value is the cross section area S, beam length is L, height is h, width is b, moment of inertia is I, beam density is ρ, Poisson ratio is γ and modulus of elasticity is E. Analyze this module by Euler Bernoulli beam with boundary surface crack, and transfer the vibration problem of the beam containing crack into the vibration problem that elastic hinge coupling two elastic beams. The formula of the bending beam is:

E I z 4 W x 4 + ρ S 2 W l 2 = 0 .

As a result of the existence of crack, establish two functions on the right and left sides of the crack x=x0, given Wi(x,l)=wi(x)Φ(l), in which Φ(l)=aexp[j(ωl+φ)], the result can be obtained:

w i n x = A i n c o s h k n x + B i n s i n h k n x + C i n c o s k n x + D i n s i n k n x .

In which: k(n)=ρSω(n)EI, n refers to the mode of vibration.

(1) The boundary condition:

when x=0, w1(0)=0, w1''(0)=0.

when x=L, w2(l)=0, w2''(l)=0.

(2) Internal continuity condition:

when x=x0, w1(x0)=w2(x0), w1''(x0)=w2''(x0), w1''(x0)=w2'''(x0).

(3) Compatibility condition:

when w1''(x0)=1βl(w2'(x0)-w1'(x0)).

We can obtain homogeneous equation: Ai(n),Bi(n),Ci(n),Di(n)(i=1, 2). Suppose the determinant equals to zero, we can obtain k(n). Finally, given the crack depth and location, we can get: Ai(n),Bi(n),Ci(n),Di(n)(i=1, 2), thus the crack beam function is obtained.

3.1. Calculation result of cracked beam

When the crack at x0/L=0.5, the crack depth are respectively: a/h=0.1, 0.3 and 0.5. The numerical analysis and calculation results are as follows:

(1) From Table 3, we can get the following conclusion: with the increase of crack depth, the natural frequency of cracked beam in case of 2r (r=1-n) is decreased gradually. Table 3 shows that the natural frequency of even order doesn’t change basically, so we choose positive symmetrical function to calculate explicit equation.

(2) Influence of various crack depths on the mode shape of cracked beam.

Table 3. Non-dimensional natural frequencies of simply supported beam with an edge crack

Crack depth a/h

Seen from Fig. 4, along with the increase of crack depth, the influence on the mode shape of cracked beam is increased gradually.

When the crack depth is not big, the mode shape of cracked beam and crack-free beam are almost the same. Fig. 5 shows the odd-order mode shape in case of a/h=0.1. And Fig. 4 indicates that in case of a/h=0.5, the influence of crack on the mode shape are particularly obvious, mainly because with the bigger rack depth, the local rigidity at the crack are bigger consequently.

Fig. 4. Comparison of the first-order mode shapes for various cracks depths

 Comparison of the first-order mode shapes for various cracks depths

Fig. 5. The first, third, fifth and seventh-order mode shapes of a crack beam in case of a/h=0.1

 The first, third, fifth and seventh-order mode shapes of a crack beam in case of a/h=0.1

3.2. Calculation result of explicit equation

The explicit equation of first-order nonlinear mode shape of cracked simply supported beam is:

w n / 1 * x , a 1 = a 1 w 1 * x + r = 2 n 3 a 1 3 b 1 11 r * 2 k 11 * + a 1 2 b 1   111 * m r r * m 11 * - k r r * .

Hence, the following analysis and calculation are conducted.

(1) Change in the non-linear frequencies with different given first function coefficient a1.

It can be seen from the Table 4 that as the given first function coefficient a1 increases, the calculated non-linear frequencies increase consequently.

(2) Comparison of explicit equation of cracked beam and that of crack-free beam mode shape.

Fig. 6. Comparison of linear and the non-linear first mode shape for a simply supported beam

 Comparison of linear and the non-linear first mode shape for a simply supported beam

As we can see in Fig. 6: (1) under the condition of unchanged given first function coefficient a1, along with the increase of crack depth, the maximum displacement increases; the mode shape of non-liner function becomes steepening, in addition, when the change of crack position is sharper near the border, the change of non-liner vibration is less than the liner vibration, while as the closer the crack position is, the bigger the non-linear mode shape is than the liner mode shape. (2) Under a certain crack depth, with the increasing of the given first function coefficient a1, both the non-linear function curvature and maximum displacement will increase.

Table 4. Comparison of the non-linear vibration frequencies

a 1
Crack depth a/h

3.3. Reliability analysis on cracked beam

To verify whether the efficiency of the algorithm proposed, we employed it to the reliability analysis of this cracked beam.

3.3.1. Calculation model of reliability index based on improved ACO

From structural reliability theory, geometric significance of reliable indicators can be defined as the shortest distance of the origin of coordinates to a limit state surface in the coordinate space of standard normal. Therefore, the problem of solving reliable indicators is changed to the problem of solving the minimum value in the following:

X = x 1 , x 2 , , x n T , β = m i n f ( x ) = m i n i = 1 n x i - μ x i σ x i 2 s.t. g ( x 1 , x 2 , . . . , x n ) = 0 . ,

In the Equitation (9), X=x1,x2,,xnT is the random variable of the structure, gx1,x2,...,xn=0 accounts for limit state equations composed of the above variables, and μxi, σxi are mean values and standard deviations of random variable xi, respectively.

The Eq. (9) is converted to the corresponding unconstrained optimization problem by Penalty function method, as following:

m i n F x , M = f x + M i = 1 n g x 2 ,

where M is corresponding penalty factor (an enough large positive number), the corresponding items are penalty items.

3.3.2. Find shortest failure paths based on improved ACO

The failure path is defined as any path between the state of the initial structure and the state of structure failure. For a structure, there are multiple failure paths, and each failure path includes a plurality of structural units. The weakest failure path is a path that the structural units successively turn failure with minimum reliability index until the entire structural system fails.

For ant colony algorithm to solve TSP problem, ηij=1dij is heuristic information, and dij is the distance between city i and city j. Thus, when through the improved ant colony algorithm to find the paths of structural failure, heuristic information is assumed as ηij=1β, where β is reliable indicator of the structural unit, m is the number of ants.

The steps to find the failure paths by the improved ant colony algorithm are as follows:

(1) Initializing parameters τij(0)=C.

(2) m ants are placed to n structure units.

(3) Calculating the probability for each ant to move to the next structure unit by the Eq. (4), and selecting moving direction of ants according to the probability.

(4) Modifying the tabu list, that is, the ant is transferred to the next city after choosing, and the city is placed at the tabu list of the individual ant.

(5) If the cities in the collection has finished traversing, that is k<m, skipping to the step (3). Otherwise, it will execute the step (6).

(6) Updating the amount of information on the paths based on the Eq. (2) and the Eq. (3).

(7) If a shutdown condition is met, the loop ends and the calculation result is output. Otherwise, emptying the tabu list and skipping to the step (2).

3.3.3. Solution of reliability indexes β

Under the premise of using improved ant colony algorithm to find the failure paths, it also is used to solve reliability index. We adopt the proposed improved ant colony algorithm to calculate the reliability indexes of the structure, among which, the set ant number is m=5 for 30 iterations. Iterative convergence process of the improved ant colony algorithm is shown in Fig. 7, which shows the improved ACO algorithm converges to achieve the optimal solution in the 6th iteration process. The result obtained by JC method is β=2.9272, the result obtained by PSO method is β=2.9237 [20], the result obtained by GA method is β=2.9241 [21], and the result obtained by the method we introduced is β=2.9232. Compared the former data, they are consistent. Moreover, iterative convergence speed of the improved ant colony algorithm is fast, and simply setting a small number of ants can achieve better computational accuracy, and the algorithm has good applicability to deal with complex and limit state equations of problems.

Fig. 7. The iterative process

 The iterative process

4. Conclusions

The structural nonlinear vibration reliability analysis is the essential part in the structural design. This paper finds the weakest failure path, and then solves the reliability index of the structure through the improved ant colony algorithm. Compared with the calculation results, the improved ant colony algorithm for the calculation of structural reliability index is feasible. The improved ant colony algorithm can overcome the restriction of traditional method which is applicable only to convex function and cumbersome derivation process. In computational accuracy and efficiency, the ant colony algorithm can meet the project requirements, having strong applicability and certain application value.


The work was supported by Natural Science Foundation of Hebei Province, China. (No. E2012402030) and the Program of Selection and Cultivating of Disciplinary Talents of Colleges and Universities in Hebei Province (BR2-206) and the Foundation of Hebei Housing and Urban Rural Construction Committee (2009-128).


  1. Stefania A. Reliability based approach for structural design and assessment: performance criteria and indicators incurrent European codes and guidelines. Int. J. Lifecycle Performance Engineering, Vol. 1, Issue 1, 2012, p. 64-91. [Search CrossRef]
  2. Frangopol D. M., Strauss A. Bridge reliability assessment based on monitoring. Journal of Bridge Engineering, Vol. 13, Issue 3, 2008, p. 258-270. [Search CrossRef]
  3. Han S. W., Kim W. T., Foutch D. A. Seismic behavior of HSS bracing members according to width-thickness ratio under symmetric cyclic loading. Journal of Structural Engineering, Vol. 133, Issue 2, 2007, p. 264-273. [Search CrossRef]
  4. Ibarra L. F., Medina R. A., Krawinkler H. Hysteretic models that incorporate strength and stiffness deterioration. Earthquake Engineering and Structural Dynamics, Vol. 34, Issue 12, 2005, p. 1489-1511. [Search CrossRef]
  5. Petrini F., Manenti S., Gkoumas K., Bontempi F. Structural design and analysis of offshore wind turbines from a system point of view. Wind Engineering, Vol. 34, Issue 1, 2010, p. 85-108. [Search CrossRef]
  6. Lehman D. E., Roeder C. W., Herman D., Johnson S., Kotulka B. Improved seismic performance of gusset plate connections. Journal of Structural Engineering, Vol. 134, 2008, p. 840-901. [Search CrossRef]
  7. Arangio S., Bontempi F., Ciampoli M. Structural integrity monitoring for dependability. Structures and Infrastructures, Vol. 7, Issue 1-2, 2011, p. 75-86. [Search CrossRef]
  8. Crosti C. Structural analysis of steel structures under fire loading. Acta Polytechnica, Vol. 49, Issue 1, 2009, p. 21-28. [Search CrossRef]
  9. Stefania A. Reliability based approach for structural design and assessment: performance criteria and indicators in current European codes and guidelines. International Journal of Lifecycle Performance Engineering, Vol. 1, Issue 1, 2012, p. 64-91. [Search CrossRef]
  10. Vrouwenvelder T., Scholten N. Assessment criteria for existing structures. Structural Engineering International, Vol. 20, Issue 1, 2010, p. 62-65. [Search CrossRef]
  11. Soheil G., Nikbakhsh J. An ant colony algorithm for solving fixed destination multi-depot multiple traveling salesmen problem. Journal Applied Soft Computing, Vol. 11, Issue 1, 2011, p. 10-16. [Search CrossRef]
  12. Dorigo M., Gambardella L. M. Ant colony system: a cooperative learning approach to the traveling salesman problem. IEEE Transactions on Evolutionary Computation, Vol. 1, Issue 1, 1997, p. 53-66. [Search CrossRef]
  13. Pablo R., Ismael R., Rubio F. Studying the application of ant colony optimization and river formation dynamics to the Steiner tree problem. Evolutionary Intelligence, Vol. 4, Issue 1, 2011, p. 51-65. [Search CrossRef]
  14. Kaveh A., Hassani B., Shojaee S., Tavakkoli S. M. Structural topology optimization using ant colony methodology. Engineering Structures, Vol. 30, Issue 9, 2008, p. 2559-2565. [Search CrossRef]
  15. Zhang Z., Li H., Huang L. The structural topology optimization based on ant colony optimization. Chinese Journal of Applied Mechanics, Vol. 28, Issue 3, 2011, p. 226-231. [Search CrossRef]
  16. Ahuja A., Pahwa A. Using ant colony optimization for loss minimization in distribution networks. Proceedings of the 37th Annual North American Power Symposium, Ames, Iowa, 2005. [Search CrossRef]
  17. Aghababa M. P., Amrollahi M. H., Borjkhani M. Application of GA, PSO, and ACO algorithms to path planning of autonomous underwater vehicles. Journal of Marine Science and Application, Vol. 11, Issue 3, 2012, p. 378-386. [Search CrossRef]
  18. Suo J. J. Study on disambiguation by statistical and logical reasoning. Advanced Materials Research, Vol. 108-111, 2010, p. 1080-1085. [Search CrossRef]
  19. Mohan B. C., Baskaran R. A survey: ant colony optimization based recent research and implementation on several engineering domain. Expert Systems with Applications, Vol. 39, Issue 4, 2011, p. 4618-4627. [Search CrossRef]
  20. Kennedy J., Eberhart R. Particle swarm optimization. Proceedings of the IEEE International Conference on Neural Networks, Perth, Australia, 1995. [Search CrossRef]
  21. Xie G., Zhang J., Li J. Slope reliability analysis based on improved genetic algorithm. Rock and Soil Mechanics, Vol. 30, Issue 6, 2009, p. 283-288. [Search CrossRef]