Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 607418, 11 pages

doi:10.1155/2012/607418

Research Article

A Hybrid Particle Swarm Optimization (PSO)-Simplex Algorithm for Damage Identification of Delaminated Beams

Xiangdong Qian,¹ Maosen Cao,1,2 Zhongqing Su,³ and Jiangang Chen³

¹Department of Engineering Mechanics, Hohai University, Nanjing 210098, China

²Institute of Fluid Flow Machinery, Polish Academy of Sciences, 80-952 Gdansk, Poland

³Department of Mechanical Engineering, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong Special Administrative Region, Hong Kong

Received 5 May 2012; Revised 24 July 2012; Accepted 29 July 2012

Academic Editor: Alessandro Marzani

Copyright © 2012 Xiangdong Qian et al. 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.

Abstract.

Delamination is a type of representative damage in composite structures, severely degrading structural integrity and reliability. The identification of delamination is commonly treated as an issue of nondestructive testing. Differing from existing studies, a hybrid optimization algorithm (HOA), combining particle swarm optimization (PSO) with simplex method (SM), is proposed to identify delamination in laminated beams. The objective function of the optimization problem is created using delamination variables (optimization parameters) together with actually measured modal frequencies. The HOA adopts a hierarchical and cooperative regime of global search and local search to optimize the objective function. The PSO performs global search for objective function space to achieve a preliminary solution specifying a local potential space. Initialized by this preliminary solution, the SM executes local search for the local potential space to explore the optimal solution. The HOA is validated by a series of simulated delamination scenarios, and the results show that it can identify delamination in laminated beams with decent accuracy, reliability and efficiency. The method proposed holds promise for establishing online damage detection system beneficial for health monitoring of laminated composite structures.

1. Introduction

Composite structures have been extensively used in aerospace, civil, and mechanical fields due to their superior fatigue characteristics, damage tolerance, and stiffness to weight ratio compared to that of metals [1]. In laminated composite structures, delamination is a type of representative damage, commonly occurring as crack in matrix materials, separation of adjoining piles, imperfect bonding, and broken fibers [2]. The presence of delamination may severely influence the mechanical properties of composite structures. Therefore, damage detection of local and tiny delamination in the framework of structural health monitoring (SHM) has become a technological and scientific issue of significant importance [3–7].

Traditionally, delamination identification is performed by model-free nondestructive testing (NDT) methods such as acoustic, ultrasonic, radiographs, magnetic, thermal, and eddy current [8–12]. In spite of distinct capabilities of these methods, their limitations for use in SHM are noteworthy and can be described as (1) a priori knowledge about the vicinity of the delamination is required; (2) ready accessibility to the structural portion of interest is necessary; (3) commonly inapplicable for determining the details of delamination, except its occurrence; and (4) the procedure to carry out these methods is usually time consuming and costly. These deficiencies limit the use of these methods to develop an online system of delamination identification for laminated structures.

As an alternative of model-free NDT methods, model-dependent methods of delamination identification have received increasing interest in the SHM community [3, 13–16]. This type of methods uses analytical or numerical model together with experimentally acquired dynamic responses as damage indictor to identify delamination in a composite structure. Damage indictor is used to calibrate the analytical or numerical model that changes itself with varying delamination variables. A practical damage indictor is modal frequency due to its effectiveness in reflecting delamination [17–19] and facilitation in experimental acquisition.

In general, model-dependent delamination identification can be tactically treated as an optimization problem, for which the objective function is created by relating measured modal frequencies to delamination variables [14–16]. An appropriate optimization method is utilized to optimize the objective function, giving the estimate of delamination parameters. Existing optimization methods can be categorized into local- and global-search methods [20, 21]. The local-search methods, typically Nelder-Mead simplex method (SM) [22], feature great efficiency in local search due to small cost of function evaluation, but they are somewhat susceptible to local minima [23]. The global-search algorithms, particularly particle swarm optimization (PSO) [24–26], have the merit of searching an optimal solution by casting off local minima along with the inadequacy of low efficiency in convergence.

Differing from existing studies, this study addresses a hybrid optimization algorithm (HOA) featuring cooperative PSO and SM for the identification of delamination in laminated composite beams. The HOA adopts a hierarchical regime of global search using the PSO incorporating local search using the SM to explore an optimal solution indicating delamination. The HOA is demonstrated by a series of simulated delamination cases, and the results show that it can identify delamination in composite beams with great accuracy, reliability, and efficiency.

2. Optimization Problem Describing Delamination

A particular optimization problem for delamination identification in laminated composite beams is set up based on Euler-Bernoulli beam theory [27] and delamination indicators referring to modal frequencies. The procedure to formulate the optimization problem, stated below, is applicable to any type of delaminated beams considering different boundary conditions. Without loss of generality, the formulation is exemplified by a double-layered clamped-clamped laminated beam bearing a delaminated segment, shown in Figure 1. The beam has a rectangular cross-section of height , width , and length , and the delaminated segment of span consists of Subbeams 2 and 3, specified by elasticity moduli and , density of materials and , and heights and , respectively. For vibration analysis, the entire beam is treated as a combination of four subbeams, of which each one is modeled as an Euler-Bernoulli beam that is governed by the following equation of motion:

where , and denote the transverse displacement, modulus of elasticity, density of material, cross-sectional moment of inertia, and cross-sectional area for the th subbeam, respectively.

(a)

(b)

(c)

Figure 1: Model for double-layered clamped-clamped laminated beam bearing a delaminated segment. (a) geometry, (b) free bodies, (c) continuity conditions: compatibility of axial deformations (left) and equilibrium of axial forces (right).

Treating (2.1) using variable separation method, mode shapes of each subbeam can be resulted in the form

where is the nondimensional -coordinate along beam length, and denotes the nondimensional modal frequencies. is considered to specify a constrained mode model; that is, subbeams 2 and 3 are stipulated to have identical transverse deformations while they can slide along the delaminated surface. The selection of constrained mode model rather than the free mode model is attributed to the insignificant relative transverse deformation between subbeams 2 and 3 for first several mode shapes, adopted for delamination identification.

An eigenvalue problem can be derived from (2.2) by taking into account the continuity and boundary conditions detailed in Figures 1(b) and 1(c), given as

where and signify the left and right ends of the delamination segment, is a column vector of twelve arbitrary constants, and is a matrix containing entries made up of elastic and geometrical properties. For a nontrivial solution the determinant of has to vanish to yield:

where denotes the determinant of matrix . In (2.4), elastic and geometrical properties are assumed to be known quantities arising from regular measurement, and thus is simply and solely related to the delamination parameters, and . This relation sets up a paradigm of delamination identification, briefed as follows. (1) Modal frequencies are acquired by actual measurement, (2)  , are individually inputted into (2.4) to establish a simultaneous system of characteristic equations, and (3) the simultaneous system is solved to produce the estimate of delamination parameters. In spite of clear technology route for the paradigm, the treatment of the simultaneous system is extremely complicated due to tremendous complexity of transcendental equation (2.4). Herein, to cope with the simultaneous system, an optimization strategy is adopted, of which the objective function is expressed as

where is specified by measured , in which the known quantities of elasticity and geometry are omitted for clarity. To estimate and , a particular optimization problem is defined by

Two modal frequencies, matching two unknowns and , are ideally suited to set up the optimization problem as defined in (2.6). Considering the uncertainty in solving optimization problem, three or even more modal frequencies are needed to increase robustness in problem description, leading to an intricate optimization problem.

3. Hybrid Optimization Algorithm

A hybrid optimization algorithm (HOA) is proposed to solve the optimization problem for delamination identification. The HOA is established based on combining two well-known optimization methods: particle swarm optimization (PSO) [24–26] and simplex method (SM) [22]. The PSO has conspicuous capability in global search while the SM owns distinctive merit in local search [22, 24]. As a combinative offspring, the HOA features a hierarchical and cooperative regime of global search using PSO and local search relying on SM. The global and local searches refer to “global screening” and “local canvassing,” respectively. The embodiment of the HOA is presented in Figure 2, and the procedure to implement the algorithm is briefly described as follows.

Figure 2: Illustrations of “global screening” and “local canvassing” of HOA. The dashed arrows symbolize evolution of interim global best solutions, filled polygon signifies preliminary solution, and star denotes optimal solution. (PS: particle swarm; SM: simplex method).

3.1. Global Screening

A sequence of particle swarms (PSs) performs “global screening” for objective function space to achieve a preliminary solution. The intention of using sequential PSs rather than a single PS for global screening is to enhance the diversity of particles, ensuring a thorough global search. The main steps involved in the operation of global screening are briefed as the following.

(i)PS₁ goes through the objective function space in a straight way, that is, with no iteration operations, to get the first interim global best solution, gbs_1. “Interim” is used to distinguish from the final global best solution of optimization, that is, optimal solution.(ii)PSi, , examines the objective function space with possible iterations until it obtains an interim global best solution, gbs_i, satisfying that gbs_i just exceeds gbs_ (i−1).(iii) Similar to that in (ii), PSn brings forth gbs_n, termed preliminary solution.

The sequence, gbs_1, gbs_2, …, gbs_i, …, gbs_n, manifests the evolution of interim global best solutions. In this process, the method of PSO fully exerts its advantage in global search to realize a preliminary solution by avoiding local minima in all probability.

3.2. Local Canvassing

The diversity and evolution of the global screening makes the preliminary solution potentially specify a local potential space probably bearing an optimal solution. The local potential space provides the SM a platform to explore optimal solution without the interference of dropping into local minima. The SM, launched by the preliminary solution, canvasses the local potential space to reach an optimal solution in a smooth and efficient manner.

As a hybrid algorithm, the HOA tactically consolidates the good character of PSO and that of SM while evading their respective deficiencies, forming an enhanced strategy for optimization identification of delamination in composite beams.

3.3. Algorithm Verification

The HOA is demonstrated on the delamination scenario VI in Table 1. As a preliminary work, the objective function involving first two modal frequencies is calculated as a forward analysis on sparse grid over domain as described in (2.6). The negative objective function is visualized in Figure 3, where the dominant peak signifying the global minimum is caused by and , that is, the actual parameters of delamination. The striking protrusion of the peak in objective function space indicates the effectiveness of the objective function in reflecting delamination.

Table 1: Delamination scenarios and identification results from HOA.

Case no.
,
Delamination type
, ,  

(Hz)
,

I
0.2, 0.4
, ,
26.135, 74.354, 129.702
0.200, 0.400

II
, ,
26.454, 74.504, 135.005
0.199, 0.401

III
, ,
26.886, 74.696, 142.414
0.200, 0.401

IV
, ,
21.757, 61.692, 108.963
0.201, 0.400

V
, ,
15.980, 44.757, 82.814
0.200, 0.401

VI
0.4, 0.6
, ,
27.114, 67.920, 144.182
0.400, 0.600

VII
, ,
27.122, 70.045, 145.022
0.399, 0.599

VIII
, ,
27.131, 73.024, 146.074
0.400, 0.601

IX
, ,
22.484, 56.889, 119.783
0.399, 0.600

X
, ,
16.278, 42.751, 87.291
0.399, 0.599

Figure 3: Interim global best solution, gbs_1 (red ball), along with particle swarm (asterisks) distribution for PS₁.

As an inverse analysis, the objective of the HOA is to optimize the objective function by adjusting variables and to determine delamination parameters, that is, to approach the peak in Figure 3. On the stage of global screening, PS₁, PS₂, , PS₁₅, each comprising thirty particles, are employed to screen the objective function space as per the steps stated in Section 3.1. The actions of the PS₁, PS₃, and PS₁₅ are exemplified in Figures 3, 4, and 5, respectively, in which three red balls mark the interim global best solutions, gbs_1, gbs_3, and gbs_15, respectively, and the black asterisks represent particles. gbs_15, with fitness function value −3.092, gives the preliminary solution. From Figure 3 to Figure 5, the interim global best solutions (red balls) gradually approach the optimal solution (striking peak), accounting for the evolution of interim global best solutions. On the stage of local canvassing, the SM initialized by the preliminary solution swiftly converges at an optimal solution, and by canvassing the local potential space for fifty iterations. The convergence process of the SM is depicted by the solid curve in Figure 6, in which the stable segment of the solid curve labels a fitness function value −0.005 for the optimal solution. The distinction of fitness function values between the preliminary solution and the optimal solution indicates that the precision improvement induced by the SM.

Figure 4: Interim global best solution, gbs_3 (red ball), along with particle swarm (asterisks) distribution for PS3.

Figure 5: Preliminary solution, gbs_15 (red ball), along with particle swarm (asterisks) distribution for PS₁₅.

Figure 6: Convergence processes for SM (solid curve) and PS₁₅ (dotted curve) starting with preliminary solution.

In addition, the behavior of PSO to pursue an optimal solution is also exhibited. The PS₁₅ is arranged to proceed from the preliminary solution, gbs_15. After 300 iterations, the PS₁₅ converges at gbs_300, with fitness function value being −0.01, much cruder than −0.005 for the SM. The convergence process is depicted by the dotted line in Figure 6. Clearly, the dotted curve is more zigzag than the solid one. The final particle swarm distribution of the PS₁₅ in the plane is shown in Figure 7, in which the coordinates of gbs_300 remarkably deviate from the point A, the actual delamination parameters. These results imply the weaker capability of PSO in advancing the preliminary solution to an optimal solution.

Figure 7: Particle swarm distribution in plane at the 300th iteration for PS₁₅ (Figure 6).

The investigation above demonstrates the distinctive superiority of the HOA to individual PSO and SM. A module of software for implementing HOA has been developed in Matlab program language.

4. Delamination Identification

The proposed HOA is testified on a series of delamination scenarios (Table 1) concerning a double-layered clamped-clamped delaminated beam, shown in Figure 1. These scenarios are created by using the analytical mode of delaminated beam and changing the contrast ratio of densities, elastic moduli, or heights of Subbeams 2 and 3. The elastic and geometrical properties of the delaminated beam are given as, s, s, in,  psi, and  lb/in³, from which the other data can be generated in terms of the contrast ratio.

The first three modal frequencies, , generated by analytical model using (2.1)–(2.4), with slight random noise, are used as experiment data acquired by a scanning laser vibrometer (SLV) [28] to set up an optimization problem (2.6). The resulting objective function is of considerable complexity, bringing difficulties for existing optimization methods. The HOA is utilized to tackle this optimization problem. As per the procedure described in Section 3, a sequence of fifteen PSs, each comprising thirty particles, are successively performed to get a preliminary solution. Afterwards, the SM, launched by this preliminary solution, is carried out to achieve an optimal solution. The identification results for a group of delamination cases are presented in Table 1. Comparing the estimated delamination parameters with the original ones, it can be concluded that the HOA has strong capability of identifying delamination in laminated composite beams.

5. Conclusions

A hybrid optimization algorithm (HOA) is proposed to deal with delamination identification in laminated composite structures. The HOA is based on combined particle swarm optimization (PSO) and simplex method (SM) and features several distinct characteristics: (1) the algorithm provides an innovative regime of hierarchical optimization by fusing typical global-search and local-search optimization methods; (2) the algorithm features consolidated merits of both PSO and SM while casting off their respective deficiencies; (3) the algorithm has strong adaptability for discrete and continuous optimization problems because of the property of free-of-derivative for PSO and SM; (4) the objective function of the HOA is underpinned by structural modal frequencies, allowing for low cost of measurement. The proposed method, aided by advanced sensing technology, holds promise for developing online systems of health monitoring for laminated composite beams.

Acknowledgments

M. Cao is grateful to a Marie Curie Industry Academia Partnership and Pathways Grant (Grant no. 251309 STA-DY-WI-CO) within the 7th European community Framework Programme, and acknowledges the partial support provided by a Foundation for the Author of National Excellent Doctoral Dissertation of PR China (Grant no. 201050) and a National Natural Science Foundation of China (Grant no. 50978084).

References

P. M. Pawar and R. Ganguli, “Genetic fuzzy system for damage detection in beams and helicopter rotor blades,” Computer Methods in Applied Mechanics and Engineering, vol. 192, no. 16–18, pp. 2031–2057, 2003.

G. Z. Voyiadjis, Damage in Composite Materials, Elsevier, New York, NY, USA, 1993.

Y. Zou, L. Tong, and G. P. Steven, “Vibration-based model-dependent damage (delamination) identification and health monitoring for composite structures—a review,” Journal of Sound and Vibration, vol. 230, no. 2, pp. 357–378, 2000.

P. Qiao, K. Lu, W. Lestari, and J. Wang, “Curvature mode shape-based damage detection in composite laminated plates,” Composite Structures, vol. 80, no. 3, pp. 409–428, 2007.

M. S. Cao and P. Z. Qiao, “Novel Laplacian scheme and multiresolution modal curvatures for structural damage identification,” Mechanical Systems and Signal Processing, vol. 23, no. 4, pp. 1223–1242, 2009.

Z. Wei, L. H. Yam, and L. Cheng, “Detection of internal delamination in multi-layer composites using wavelet packets combined with modal parameter analysis,” Composite Structures, vol. 64, no. 3-4, pp. 377–387, 2004.

N. Hu, T. Shimomukai, C. Yan, and H. Fukunaga, “Identification of delamination position in cross-ply laminated composite beams using S0 Lamb mode,” Composites Science and Technology, vol. 68, no. 6, pp. 1548–1554, 2008.

Y. J. Yan, L. Cheng, Z. Y. Wu, and L. H. Yam, “Development in vibration-based structural damage detection technique,” Mechanical Systems and Signal Processing, vol. 21, no. 5, pp. 2198–2211, 2007.

C. R. Farrar, S. W. Doebling, and D. A. Nix, “Vibration-based structural damage identification,” Philosophical Transactions of the Royal Society A, vol. 359, no. 1778, pp. 131–149, 2001.

Q. Wang and N. Wu, “Detecting the delamination location of beam with a wavelet transform: an experimental study,” Smart Materials and Structures, vol. 20, no. 1, Article ID 012002, 7 pages, 2011.

Z. Su and L. Ye, “Lamb wave-based quantitative identification of delamination in CF/EP composite structures using artificial neural algorithm,” Composite Structures, vol. 66, no. 1–4, pp. 627–637, 2004.

Y. J. Yan and L. H. Yam, “Detection of delamination damage in composite plates using energy spectrum of structural dynamic responses decomposed by wavelet analysis,” Computers & Structures, vol. 82, no. 4-5, pp. 347–358, 2004.

S. J. Wildy, A. G. Kotousov, B. S. Cazzolato, and J. D. Codrington, “New damage detection technique based on governing differential equations of continuum mechanics. Part I: out-of-plane loading,” in Proceedings of the 6th Australasian Congress on Applied Mechanics (ACAM 6 '10), Perth, Australia, 2010.

J. Lee, R. T. Haftka, O. H. Griffin, L. T. Watson, and M. D. Sensmeier, “Detecting delaminations in a composite beam using anti-optimization,” Structural and Multidisciplinary Optimization, vol. 8, no. 2-3, pp. 93–100, 1994.

A. Nag, D. R. Mahapatra, and S. Gopalakrishnan, “Identification of delamination in composite beams using spectral estimation and a genetic algorithm,” Smart Materials and Structures, vol. 11, no. 6, pp. 899–908, 2002.

A. Maranon, P. D. Ruiz, A. D. Nurse, J. M. Huntley, L. Rivera, and G. Zhou, “Identification of subsurface delaminations in composite laminates,” Composites Science and Technology, vol. 67, no. 13, pp. 2817–2826, 2007.

S. Dey and A. Karmakar, “Natural frequencies of delaminated composite rotating conical shells—a finite element approach,” Finite Elements in Analysis and Design, vol. 56, pp. 41–51, 2012.

P. Tan and L. Tong, “Identification of delamination in a composite beam using integrated piezoelectric sensor/actuator layer,” Composite Structures, vol. 66, no. 1–4, pp. 391–398, 2004.

R. Matsuzaki and A. Todoroki, “Wireless detection of internal delamination cracks in CFRP laminates using oscillating frequency changes,” Composites Science and Technology, vol. 66, no. 3-4, pp. 407–416, 2006.

G. M. L. Gladwell, “Inverse problems in vibration,” Applied Mechanics Reviews, vol. 49, no. 10, pp. S25–S34, 1996.

R. Fletcher, Practical Methods of Optimization, A Wiley-Interscience Publication, John Wiley & Sons, New York, NY, USA, 2000.

J. A. Nelder and R. Mead, “A simplex method for function minimization,” The Computer Journal, vol. 7, pp. 308–320, 1965.

C. T. Kelley, Iterative Methods for Optimization, SIAM, Philadelphia, Pa, USA, 1999.

K. E. Parsopoulos and M. N. Vrahatis, “Recent approaches to global optimization problems through particle swarm optimization,” Natural Computing, vol. 1, no. 2-3, pp. 235–306, 2002.

M. Clerc and J. Kennedy, “The particle swarm-explosion, stability, and convergence in a multidimensional complex space,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 1, pp. 58–73, 2002.

O. Begambre and J. E. Laier, “A hybrid Particle Swarm Optimization—Simplex algorithm (PSOS) for structural damage identification,” Advances in Engineering Software, vol. 40, no. 9, pp. 883–891, 2009.

S. Timoshenko, Vibration Problems in Engineering, John Wiley & Sons, New York, NY, USA, 1974.

R. B. Bai, M. S. Cao, Z. Q. Su, W. Ostachowicz, and H. Xu, “Fractal dimension analysis of higher-order mode shapes for damage identification of Beam structures,” Mathematical Problems in Engineering, vol. 2012, Article ID 454568, 16 pages, 2012.