Predictability of Inverse Impact Force Location as Affected by Measurement Noise

1 Civil Engineering and Mechanics Laboratory, Abdelmalek Essaâdi University, 91001 Tangier, Morocco 2 Communications Systems and Detection Laboratory, Abdelmalek Essaâdi University, 93002 Tetouan, Morocco 3 Department of Physics, Faculty of Sciences at Tetouan, P.O. Box 2121, M’Hannech II, 93002 Tetouan, Morocco 4Civil and Environmental Engineering Laboratory, Institute of Applied Sciences at Lyon, 20 Albert Einstein avenue, 69621 Villeurbanne Cedex, France 5 University of Lyon, 69622 Lyon, France 6 IFSTTAR, LBMC, UMR-T9406, Université Lyon 1, Villeurbanne, France


Introduction
Identification of impact force location for impact events occurring on elastic structures can be performed by various methods that were proposed in the literature [1][2][3].To review briefly some of the important contributions in this filed, Martin and Doyle [4] have described how to find the location of an impact force using dynamic response measurements.They proposed a solution procedure using the spectral element method with a stochastic iterative search.Experimentally measured acceleration responses from two frame structures were used to achieve force localization by minimizing a fitness function.A Genetic Algorithm was used to guess iteratively the minimum through monitoring the actual error associated to a given sampling generation.The process enabled to discriminate between good and bad guesses and gave at convergence the correct impact location.An alternative technique which employs the arrival time of each frequency component of a pulse detected by means of wavelet transform was proposed by Inoue et al. [5].But this approach suffers from the lack of accuracy in measurement of small arrival times of signals.Yen and Wu [6,7] have used multiple strain responses along with a mutuality relationship based on Green's functions and measured strains to achieve identification of force location on two-dimensional plate-like structures.Choi and Chang [8] minimized the error between measured strain responses in PZT sensors and numerically evaluated impact force locations.Shin [9] proposed a technique for identifying the force location using modal displacements and transient signal measured by accelerometers.
As mentioned by Doyle [10], the implicit character of force localization requires considering the inverse problem of force identification as being rather associated to two decoupled subproblems: localization of impact force and reconstruction of the time force signal.The reason for this is that the inverse problem of force identification cannot be handled by just adding the force point coordinates as extra unknown parameters to the discrete vector of time force signal values.These positions intervene in fact as implicit parameters, in contrast with the force history values which are explicit.Iterations are then required for solution of force positions unknowns in order to determine them at first.Treating all force history values to be implicit parameters is not a good choice since it would increase dramatically the number of unknowns and would consequently penalize the computational cost.When the force localization is obtained through the solution of a nonlinear mathematical program, the force time signal reconstruction can then be performed by means of a regularization deconvolution technique.
In this work, separation of the localization and reconstruction phases of the impact force inverse problem is adopted.Elastic structures subjected to nonpunctual impacts for which the resulting force field can be assumed to be uniformaly distributed over a finite domain of the structure are considered.The localization problem is solved by using minimization procedures that are of evolutionary type.Two methods are examined in the following: Genetic Algorithm (GA) based strategy [11] and the Particle Swarm Optimization (PSO) approach [12,13].Focus will be done on the particular role related to noise affecting strain measurement as to the resulting perturbations that produce and which may impede these algorithms to converge towards the exact problem solution.
The localization problem is derived straightforwardly from the reciprocity Maxwell-Betti theorem which is valid for any elastic structure.However, as the direct expression of the fitness function formed by this theorem is ill conditioned, because it admits a lot of trivial meaningless solutions which are associated to the fixed boundary conditions, the fitness function is modified [7].The introduced modification transforms the nonlinear mathematical program to a unique solution problem and removes the trivial parasitic solutions which complicate optimization process convergence to the real solution.Some extra constraints are also introduced in order to guide the exploration of the optimal solution; these are associated to the geometric pieces of evidence that describe the domain containing the tracked unknowns.
The aim is to examine the predictability of impact location in this situation where noise is present in observation measurement.This means determining the amount of noise that can be tolerated and also selecting between the two proposed algorithms the most suited one that can be used to conduct solution of force localization problem with sufficient accuracy and minimal error.

Direct Problem Formulation.
Although the problem can be stated for any elastic structural system, a simplified model having the form of a beam with a rectangular uniform section is considered, Figure 1.The beam is assumed to be simply supported on both ends.It has length , width , and height .It is assumed to be made from a homogeneous and isotropic elastic material with Young's modulus  and density .The applied force modelling impact is assumed to result from a uniformly distributed pressure, , which is applied on a rectangular patch as shown in Figure 1.The pressure rectangle is assumed to be centred on  0 and of length 2.
The differential equation for transverse vibrations of a beam writes where  is the transverse displacement,  is the axial position,  is the time,  is the cross-section area,  is the density,  is Young's modulus,  is the moment of inertia,  is the viscous damping coefficient, and (, ) is the applied pressure over The considered boundary conditions are as follows: Using modal superposition and Duhamel's integral, the transient dynamic solution of the impacted beam in terms of axial strain  at location   is given by the following convolution problem: where (  , ,  − ) is the time response function between point  and the sensor location   .Time discretisation of (3) under the assumption that the pressure (, ) is uniform yields the following linear system: where   designates the vector containing the discrete values of axial strain at point location indicated by index ,  is the vector containing discrete time values of the applied uniform pressure, and   ( 0 , ) is the Toeplitz-like matrix.This last has the following form: where  is the number of time steps.
An explicit expression can be obtained for matrix .Let us denote by Δ the time step used in discretization and by  the order of modal truncation; the impulse response function giving the deformation of the top fiber for a section having abscissa   is given by with where   = ( 2  2 / 2 )√/ and   are, respectively, the circular eigenfrequency and damping ratio for a given eigenmode .

Localization Problem as Perturbed by the Presence of Measurement
Noise.The problem of finding the impact location for the beam considered in the present study consists in identifying the impact patch centre position  0 and parameter  defining the extent of the impacted zone.Using (4) to (6), the responses measured by strain sensors placed at points having the abscissa   and   can be expressed, respectively, under the following form: Equation ( 7) can be used to prove the following important commutativity property which does not contain the pressure vector : In this way, ( 0 , ) appears to be the solution of the following equation: where ‖ ⋅ ‖ is the Euclidian norm.
There are a lot of other trivial solutions for (9) for which the factor sin(/) sin(V/) appearing in the second hand side of ( 6) vanishes for all values of  ∈ [1, 𝑀].These are associated to the boundaries of the domain containing ( 0 , ); that is, { = 0;  = ; V = 0; V = }.To get the exact solution of impact location, the parasitic solutions should be withdrawn from (9).
Denoting by   the number of sensors used, the fitness function to be minimized in order to find the impact zone parameters ( 0 , ) is proposed under the following form: with   = ‖  (, V)‖ 2 + ‖  (, V)‖ 2 defining the weighting coefficients that are introduced in order to remove the parasitic solutions, as they vanish also for { = 0;  = ; V = 0; V = }.
In the presence of measurement noise, the perturbed fitness function takes the following form: with where   is a random number belonging to the interval [−1, 1] and ]  and ]  , respectively designate the noise level present in measurement delivered by sensor  and sensor , respectively.
To stabilize the minimization procedure and obtain the unique physical solution, the unconstrained mathematical program defined by (10) is constrained by adding the following geometrical bounds conditions: To solve the mathematical program defined by ( 6) and ( 10)-( 13), the PSO algorithm and the GA are considered.Their performance will then be assessed as function of the random noise level present in strain measurement.

PSO Algorithm.
A PSO based method was proposed initially by Eberhart and Kennedy [12,13].This approach has gained since then considerable interest as being one of the most promising optimization methods that is able to provide high speed and high accuracy.PSO mimics the social behavior that a population of individuals adapts to its environment by returning to promising regions that were previously discovered [14].This adaptation to the environment is a stochastic process that depends on both the memory of each individual, called particle, and the knowledge gained by the population, called swarm.
In the simplest numerical implementation of this method, each particle is characterised by four attributes: the position vector in the search space, the velocity vector, the best position achieved in its track, and the best position achieved by the swarm.The process steps can be outlined as follows.
Step 1. Generate the initial swarm involving  given particles placed at random.
Step 2. Calculate the new velocity vector of each particle, based on its actual attributes.
Step 3. Calculate the new position of each particle from the current position and its new velocity vector.
Step 4. If the termination condition is satisfied, stop.Otherwise, go to Step 2.
To be more specific, the new velocity vector of the ith particle at time  + 1, denoted V +1  , is calculated according to the following Shi and Eberhart [15] formula: In ( 14),   1 and   2 are random numbers between 0 and 1,    is the best position of the ith particle in its track, and    is the best position of the swarm.There are three problem dependent parameters that fix performance of this algorithm, namely, the inertia of the particle   and the two trust parameters  1 and  2 .
The new position of the ith particle at time , denoted  +1  , is then calculated as follows: where    is the current position of the ith particle at time .The ith particle actual position enables to determine the best position in its track    .When considering all the particles, the global best position of the swarm    is then obtained.PSO algorithm works such that particles concentrate on the best search position of the swarm.They cannot easily escape from the local optimal solution since the search direction vector V +1  calculated by ( 14) always includes the direction vector to the best search position of the swarm.This shows the major feature of PSO algorithm as being a robust process of continuous enhancement for optimum search.
In the presence of constraints, a particle move should be restricted in order to remain in the feasible solution space by examining the given constraints.A modified PSO version was introduced for constrained problems in order to manage this situation [16].[11].It is a probabilistic optimization method that is able to achieve global search by mimicking natural biological evolution.GA operates on a population of individuals called the set of potential solutions.Each individual is represented by an encoded string (chromosome) that contains the decision variables (genes).Traditionally, GA uses binary strings as chromosome representation.

Genetic Algorithm. Genetic Algorithm (GA) was firstly introduced by Holland
The GA has an iterative procedure structure that comprises generally the following five main steps.
Step 2. Evaluation of the performance of each individual or chromosome (  ) of the population, by means of a fitness function to be maximized.
Step 3. Selection of individuals for the reproduction of a new population.
Step 4. Application of genetic operators: Crossover and Mutation.
Step 5. Iteration of Steps 2 to 4 until a termination criterion is fulfilled.
In the localization problem considered in this work, the candidate solution is the centre position  0 and the extent of the impact zone .These variables are then coded in a chromosome using a binary coding scheme.
To start the algorithm, an initial population of individuals (chromosomes) is defined.The GA is configured, so that it creates a fixed number of initial individuals at random from the whole feasible solution space.An important parameter in initialization is the population size.In general, the population size affects both the ultimate performance and the efficiency of GA and should be determined in a case by case study.

Results and Discussion
A pinned-pinned beam having the following material and geometric properties is considered:  = 7.06 × 10 10 Pa;  = 0.5 m;  =  = 5 × 10 −3 m;  = 2660 kg⋅m −3 ;   = 2%;  0 = 0.417 m;  = 0.0417 m.The time interval considered has the duration   = 1 s.The beam is assumed to be subjected to a half sine pulse pressure having the shape and spectral content depicted in Figure 2. The maximum pressure is taken to be 10 5 Pa.
Four gauge strain sensors were used in this study.Their labeling and positions are indicated in Table 1.Previous studies conducted by means of PSO algorithm have shown that, while considering the noise free problem, localization needs at least three sensors to be achieved with adequate accuracy.Here, a fourth sensor has been added in order to enhance the performance of the localization inverse problem in the presence of measurement noise.Table 1: Positions of the gauge sensors considered for the measurement of axial strain at the upper beam fiber.
Three noise levels were considered in order to test robustness of both PSO and GA based algorithms.They correspond to the following.Case 1. Noise free case where noise amplitude is set to the value ]  = 0% for the four sensors.
Case 2. Noise level is fixed at ]  = 2% for the four sensors.
Case 3. Noise level is fixed at ]  = 5% for the four sensors.
For GA algorithm the first values of the unknown parameters impact centre and extent were initialized with  0 =  = 1.25 × 10 −2 .
Figure 3 gives evolution of impact location characteristics as a function of iterations in case of noiseless problem.One can see the difference existing between these two algorithms; PSO algorithm gives the exact solution of the impact location.On the opposite, GA has not converged to the right solution.
Figure 4 gives evolution of impact location characteristics as function of iterations in case of noisy measurement conditions with a noise level ]  = 2%.One can observe that GA has not converged at all, as the calculated solution is too far from the exact solution.The relative error reached 92% for the impact zone extent.Meanwhile, PSO based algorithm has continued to give the exact solution with only a moderate error, representing 13% of the impact zone extent.
Figure 5 gives the evolution of impact location characteristics as function of iterations in case of noise level value given by ]  = 5%.Unexpectedly, GA is better in this condition than for noise level ]  = 2% as the relative error for impact zone extent has decreased from 92% to 68%.Moreover, one can see that both GA and PSO algorithms fail to give the exact solution as the error on the impact zone extent is too large, 61% and 68%, respectively.
Tables 2 and 3 recall the obtained relative error as the function of the algorithm used and noise level.One can see that GA fails to predict the exact solution, with a behavior that is not monotonous, while PSO based algorithm performs well for small measurement noise levels not exceeding 2%.The performance decreases after that when identifying the impact zone extent, while it remains sufficiently good for the impact zone centre as the relative error is smaller than 5%.
To emphasize the irregular behavior of GA, a comparison has been made between the results of two runs of GA. Figure 6 shows a comparison between the obtained evolutions.They are not the same because of the stochastic nature of GA.GA is found thus to be inappropriate for the inverse force location problem as it is too sensitive to noise level and its repeatability is too poor.
On the opposite, PSO algorithm runs always in the same manner with perfect repeatability.

Conclusion
Based on the separation approach that decouples force location from force signal reconstruction in an inverse     as a function of measurement noise level.It was found that PSO algorithm continues to achieve exact prediction of the solution even in the presence of 2% of noise intensity, while GA algorithm fails at that noise as the associated error was too large.Both of these algorithms fail to predict correct impact extent when noise level reaches 5%, even if they continue to predict the position of the impact centre with rather acceptable accuracy.
To assess predictability in a large sense, system model noise should also be integrated in the future in conjunction with measurement noise.

Figure 1 :
Figure 1: An elastic beam with a uniform rectangular crosssection and loaded with a distributed uniform pressure; four sensors positions are indicated.

Figure 2 :
Figure 2: Characteristics of the impulse impact pressure; (a) time history and (b) frequency content.

Figure 3 :
Figure 3: Noise free case; evolution of impact zone characteristics.(a) Impact centre and (b) impact zone extent.

Figure 6 :
Figure 6: Noise level ]  = 5% for the four sensors; evolution of impact centre and extent as function of iterations (a) first test GA 1 and (b) second test GA 2.

Table 2 :
Relative error affecting impact zone extent  as compared with the exact solution.
impact problem occurring on an elastic beam, robustness of a particular localization procedure was analyzed.This uses a modified fitness function derived from Maxwell-Betti theorem by applying some filtering coefficients that

Table 3 :
Relative error affecting impact zone centre  0 as compared with the exact solution.
enable to remove parasitic solutions.Solution of the obtained constrained nonlinear mathematical program that provides the impact zone location was performed by GA and PSO based algorithms.Predictability of force location was studied