Joint parameter estimation from magnetic resonance and vertical electric soundings using a multi‐objective genetic algorithm

Magnetic resonance sounding (MRS) has increasingly become an important method in hydrogeophysics because it allows for estimations of essential hydraulic properties such as porosity and hydraulic conductivity. A resistivity model is required for magnetic resonance sounding modelling and inversion. Therefore, joint interpretation or inversion is favourable to reduce the ambiguities that arise in separate magnetic resonance sounding and vertical electrical sounding (VES) inversions. A new method is suggested for the joint inversion of magnetic resonance sounding and vertical electrical sounding data. A one‐dimensional blocky model with varying layer thicknesses is used for the subsurface discretization. Instead of conventional derivative‐based inversion schemes that are strongly dependent on initial models, a global multi‐objective optimization scheme (a genetic algorithm [GA] in this case) is preferred to examine a set of possible solutions in a predefined search space. Multi‐objective joint optimization avoids the domination of one objective over the other without applying a weighting scheme. The outcome is a group of non‐dominated optimal solutions referred to as the Pareto‐optimal set. Tests conducted using synthetic data show that the multi‐objective joint optimization approximates the joint model parameters within the experimental error level and illustrates the range of trade‐off solutions, which is useful for understanding the consistency and conflicts between two models and objectives. Overall, the Levenberg‐Marquardt inversion of field data measured during a survey on a North Sea island presents similar solutions. However, the multi‐objective genetic algorithm method presents an efficient method for exploring the search space by producing a set of non‐dominated solutions. Borehole data were used to provide a verification of the inversion outcomes and indicate that the suggested genetic algorithm method is complementary for derivative‐based inversions.     

Multi‐objective particle swarm optimization for generating optimal trade‐offs in reservoir operation

A multi‐objective particle swarm optimization (MOPSO) approach is presented for generating Pareto‐optimal solutions for reservoir operation problems. This method is developed by integrating Pareto dominance principles into particle swarm optimization (PSO) algorithm. In addition, a variable size external repository and an efficient elitist‐mutation (EM) operator are introduced. The proposed EM‐MOPSO approach is first tested for few test problems taken from the literature and evaluated with standard performance measures. It is found that the EM‐MOPSO yields efficient solutions in terms of giving a wide spread of solutions with good convergence to true Pareto optimal solutions. On achieving good results for test cases, the approach was applied to a case study of multi‐objective reservoir operation problem, namely the Bhadra reservoir system in India. The solutions of EM‐MOPSOs yield a trade‐off curve/surface, identifying a set of alternatives that define optimal solutions to the problem. Finally, to facilitate easy implementation for the reservoir operator, a simple but effective decision‐making approach was presented. The results obtained show that the proposed approach is a viable alternative to solve multi‐objective water resources and hydrology problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Obtaining refined first‐order predictive models of linear structural systems

This study presents an effective method for identifying predictive models and the underlying modal parameters of linear structural systems using only measured output and excitation time histories obtained from dynamic testing. The system under examination is modelled as a first‐order multi‐input multi‐output time‐invariant system, and the structural model is realized using the Eigensystem Realization Algorithm together with the Observer/Kalman filter IDentification algorithm. The identified state‐space model is further refined using a non‐linear optimization technique based on sequential quadratic programming. The numerical examples show that the developed methodology performs very well even in the presence of inadequate instrumentation and measurement noise, and that the methodology is highly capable of creating realistic predictive models of structural systems, as well as estimating their underlying modal parameters. Copyright © 2002 John Wiley & Sons, Ltd.     

Optimal seismic design of steel frame buildings based on life cycle cost considerations

A multi‐objective optimization procedure is presented for designing steel moment resisting frame buildings within a performance‐based seismic design framework. Life cycle costs are considered by treating the initial material costs and lifetime seismic damage costs as two separate objectives. Practical design/construction complexity, important but difficult to be included in initial cost analysis, is taken into due account by a proposed diversity index as another objective. Structural members are selected from a database of commercially available wide flange steel sections. Current seismic design criteria (AISC‐LRFD seismic provisions and 1997 NEHRP provisions) are used to check the validity of any design alternative. Seismic performance, in terms of the maximum inter‐storey drift ratio, of a code‐verified design is evaluated using an equivalent single‐degree‐of‐freedom system obtained through a static pushover analysis of the original multi‐degree‐of‐freedom frame building. A simple genetic algorithm code is used to find a Pareto optimal design set. A numerical example of designing a five‐storey perimeter steel frame building is provided using the proposed procedure. It is found that a wide range of valid design alternatives exists, from which a decision maker selects the one that balances different objectives in the most preferred way. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic design of RC structures: A critical assessment in the framework of multi‐objective optimization

The assessment of seismic design codes has been the subject of intensive research work in an effort to reveal weak points that originated from the limitations in predicting with acceptable precision the response of the structures under moderate or severe earthquakes. The objective of this work is to evaluate the European seismic design code, i.e. the Eurocode 8 (EC8), when used for the design of 3D reinforced concrete buildings, versus a performance‐based design (PBD) procedure, in the framework of a multi‐objective optimization concept. The initial construction cost and the maximum interstorey drift for the 10/50 hazard level are the two objectives considered for the formulation of the multi‐objective optimization problem. The solution of such optimization problems is represented by the Pareto front curve which is the geometric locus of all Pareto optimum solutions. Limit‐state fragility curves for selected designs, taken from the Pareto front curves of the EC8 and PBD formulations, are developed for assessing the two seismic design procedures. Through this comparison it was found that a linear analysis in conjunction with the behaviour factor q of EC8 cannot capture the nonlinear behaviour of an RC structure. Consequently the corrected EC8 Pareto front curve, using the nonlinear static procedure, differs significantly with regard to the corresponding Pareto front obtained according to EC8. Furthermore, similar designs, with respect to the initial construction cost, obtained through the EC8 and PBD formulations were found to exhibit different maximum interstorey drift and limit‐state fragility curves. Copyright © 2007 John Wiley & Sons, Ltd.     

Multi‐objective risk‐informed design of floor isolation systems

The design of floor isolation systems (FISs) for the protection of acceleration sensitive contents is examined considering multiple objectives, all quantified in terms of the probabilistic system performance. The competing objectives considered correspond to (i) maximization of the level of protection offered to the sensitive content (acceleration reduction) and (ii) minimization of the demand for the isolator displacement capacity and, more importantly, for the appropriate clearance to avoid collisions with surrounding objects (floor displacement reduction). Both of these objectives are probabilistically characterized utilizing a versatile, simulation‐based framework for quantifying seismic risk, addressing all important uncertainties related to the seismic hazard and the structural model. FIS performance is assessed through time‐history analysis, allowing for all important sources of nonlinearity to be directly addressed in the design framework. The seismic hazard is described through a stochastic ground motion model. For efficiently performing the multi‐objective optimization, an augmented surrogate modeling methodology is established, considering development of a single metamodel with respect to both the uncertain model parameters and the design variables for the FIS system. This surrogate model is then utilized to simultaneously support the probabilistic risk assessment and the design optimization to provide the Pareto front of dominant designs. Each of these designs establishes a different compromise between the considered risk‐related objectives offering a variety of potential options to the designer. Within the illustrative example, the efficiency of the established framework is exploited to compare three different FIS implementations, whereas the impact of structural uncertainties on the optimal design is also evaluated. Copyright © 2016 John Wiley & Sons, Ltd.     

Multi-objective optimal design for flood risk management with resilience objectives

In flood risk management, the divergent concept of resilience of a flood defense system cannot be fully defined quantitatively by one indicator and multiple indicators need to be considered simultaneously. In this paper, a multi-objective optimization (MOO) design framework is developed to determine the optimal protection level of a levee system based on different resilience indicators that depend on the probabilistic features of the flood damage cost arising under the uncertain nature of rainfalls. An evolutionary-based MOO algorithm is used to find a set of non-dominated solutions, known as Pareto optimal solutions for the optimal protection level. The objective functions, specifically resilience indicators of severity, variability and graduality, that account for the uncertainty of rainfall can be evaluated by stochastic sampling of rainfall amount together with the model simulations of incurred flood damage estimation for the levee system. However, these model simulations which usually require detailed flood inundation simulation are computationally demanding. This hinders the wide application of MOO in flood risk management and is circumvented here via a surrogate flood damage modeling technique that is integrated into the MOO algorithm. The proposed optimal design framework is applied to a levee system in a central basin of flood-prone Jakarta, Indonesia. The results suggest that the proposed framework enables the application of MOO with resilience objectives for flood defense system design under uncertainty and solves the decision making problems efficiently by drastically reducing the required computational time.     

Modification of stochastic ground motion models for matching target intensity measures

Stochastic ground motion models produce synthetic time‐histories by modulating a white noise sequence through functions that address spectral and temporal properties of the excitation. The resultant ground motions can be then used in simulation‐based seismic risk assessment applications. This is established by relating the parameters of the aforementioned functions to earthquake and site characteristics through predictive relationships. An important concern related to the use of these models is the fact that through current approaches in selecting these predictive relationships, compatibility to the seismic hazard is not guaranteed. This work offers a computationally efficient framework for the modification of stochastic ground motion models to match target intensity measures (IMs) for a specific site and structure of interest. This is set as an optimization problem with a dual objective. The first objective minimizes the discrepancy between the target IMs and the predictions established through the stochastic ground motion model for a chosen earthquake scenario. The second objective constraints the deviation from the model characteristics suggested by existing predictive relationships, guaranteeing that the resultant ground motions not only match the target IMs but are also compatible with regional trends. A framework leveraging kriging surrogate modeling is formulated for performing the resultant multi‐objective optimization, and different computational aspects related to this optimization are discussed in detail. The illustrative implementation shows that the proposed framework can provide ground motions with high compatibility to target IMs with small only deviation from existing predictive relationships and discusses approaches for selecting a final compromise between these two competing objectives.     

System identification of the Vincent Thomas suspension bridge using earthquake records

The Vincent Thomas Bridge in the Los Angeles metropolitan area, is a critical artery for commercial traffic flow in and out of the Los Angeles Harbor, and is at risk in the seismically active Southern California region, particularly because it straddles the Palos Verdes fault zone. A combination of linear and non‐linear system identification techniques is employed to obtain a complete reduced‐order, multi‐input–multi‐output (MIMO) dynamic model of the Vincent Thomas Bridge based on the dynamic response of the structure to the 1987 Whittier and 1994 Northridge earthquakes. Starting with the available acceleration measurements (which consists of 15 accelerometers on the bridge structure and 10 accelerometers at various locations on its base), an efficient least‐squares‐based time‐domain identification procedure is applied to the data set to develop a reduced‐order, equivalent linear, multi‐degree‐of‐freedom model. Although not the main focus of this study, the linear system identification method is also combined with a non‐parametric identification technique, to generate a reduced‐order non‐linear mathematical model suitable for use in subsequent studies to predict, with good fidelity, the total response of the bridge under arbitrary dynamic environments. Results of this study yield measurements of the equivalent linear modal properties (frequencies, mode shapes and non‐proportional damping) as well as quantitative measures of the extent and nature of non‐linear interaction forces arising from strong ground shaking. It is shown that, for the particular subset of observations used in the identification procedure, the apparent non‐linearities in the system restoring forces are quite significant, and they contribute substantially to the improved fidelity of the model. Also shown is the potential of the identification technique under discussion to detect slight changes in the structure's influence coefficients, which may be indicators of damage and degradation in the structure being monitored. Difficulties associated with accurately estimating damping for lightly damped long‐span structures from their earthquake response are discussed. The technical issues raised in this paper indicate the need for added spatial resolution in sensor instrumentation to obtain identified mathematical models of structural systems with the broadest range of validity. Copyright © 2003 John Wiley & Sons, Ltd.     

On the use of multi‐algorithm,genetically adaptive multi‐objective method for multi‐site calibration of the SWAT model

With the availability of spatially distributed data, distributed hydrologic models are increasingly used for simulation of spatially varied hydrologic processes to understand and manage natural and human activities that affect watershed systems. Multi‐objective optimization methods have been applied to calibrate distributed hydrologic models using observed data from multiple sites. As the time consumed by running these complex models is increasing substantially, selecting efficient and effective multi‐objective optimization algorithms is becoming a nontrivial issue. In this study, we evaluated a multi‐algorithm, genetically adaptive multi‐objective method (AMALGAM) for multi‐site calibration of a distributed hydrologic model—Soil and Water Assessment Tool (SWAT), and compared its performance with two widely used evolutionary multi‐objective optimization (EMO) algorithms (i.e. Strength Pareto Evolutionary Algorithm 2 (SPEA2) and Non‐dominated Sorted Genetic Algorithm II (NSGA‐II)). In order to provide insights into each method's overall performance, these three methods were tested in four watersheds with various characteristics. The test results indicate that the AMALGAM can consistently provide competitive or superior results compared with the other two methods. The multi‐method search framework of AMALGAM, which can flexibly and adaptively utilize multiple optimization algorithms, makes it a promising tool for multi‐site calibration of the distributed SWAT. For practical use of AMALGAM, it is suggested to implement this method in multiple trials with relatively small number of model runs rather than run it once with long iterations. In addition, incorporating different multi‐objective optimization algorithms and multi‐mode search operators into AMALGAM deserves further research. Copyright © 2009 John Wiley & Sons, Ltd.     

Bayesian model selection for ARX models and its application to structural health monitoring

A Bayesian framework for model order selection of auto‐regressive exogenous (ARX) models is developed and applied to actual earthquake response data obtained by the structural health monitoring system of a high‐rise building. The model orders of ARX models are selected appropriately by the Bayesian framework, and differ significantly from the optimal order estimated by AIC; in fact, in many cases AIC does not even give an optimal order. A method is also proposed for consistently selecting the same ‘genuine’ modes of interest from the whole set of modes corresponding to each of the identified models from a sequence of earthquake records. In the identification analysis based on building response records from 43 earthquakes over 9 years, the modal parameters of the first four modes in each horizontal direction are estimated appropriately in all cases, showing that the developed methods are effective and robust. As the estimates of natural frequency depend significantly on the response amplitude, they are compensated by an empirical correction so that the influence of the response amplitude is removed. The compensated natural frequencies are much more stable over the nine‐year period studied, indicating that the building had no significant change in its global dynamic characteristics during this period. Copyright © 2010 John Wiley & Sons, Ltd.     

Identification of modal parameters of non‐classically damped linear structures under multi‐component earthquake loading

A method for parametric system identification of classically damped linear system in frequency domain is adopted and extended for non‐classically damped linear systems subjected up to six components of earthquake ground motions. This method is able to work in multi‐input/multi‐output (MIMO) case. The response of a two‐degree‐of‐freedom model with non‐classical damping, excited by one‐component earthquake ground motion, is simulated and used to verify the proposed system identification method in the single‐input/multi‐output case. Also, the records of a 10 storey real building during the Northridge earthquake is used to verify the proposed system identification method in the MIMO case. In this case, at first, a single‐input/multi‐output assumption is considered for the system and modal parameters are identified, then other components of earthquake ground motions are added, respectively, and the modal parameters are identified again. This procedure is repeated until all four components of earthquake ground motions which are measured at the base level of the building are included in the identification process. The results of identification of real building show that consideration of non‐classical damping and inclusion of the multi‐components effect of earthquake ground motions can improve the least‐squares match between the finite Fourier transforms of recorded and calculated acceleration responses. Copyright © 2006 John Wiley & Sons, Ltd.     

Algorithms for time synchronization of wireless structural monitoring sensors

Dense networks of wireless structural health monitoring systems can effectively remove the disadvantages associated with current wire‐based sparse sensing systems. However, recorded data sets may have relative time‐delays due to interference in radio transmission or inherent internal sensor clock errors. For structural system identification and damage detection purposes, sensor data require that they are time synchronized. The need for time synchronization of sensor data is illustrated through a series of tests on asynchronous data sets. Results from the identification of structural modal parameters show that frequencies and damping ratios are not influenced by the asynchronous data; however, the error in identifying structural mode shapes can be significant. The results from these tests are summarized in Appendix A. The objective of this paper is to present algorithms for measurement data synchronization. Two algorithms are proposed for this purpose. The first algorithm is applicable when the input signal to a structure can be measured. The time‐delay between an output measurement and the input is identified based on an ARX (auto‐regressive model with exogenous input) model for the input–output pair recordings. The second algorithm can be used for a structure subject to ambient excitation, where the excitation cannot be measured. An ARMAV (auto‐regressive moving average vector) model is constructed from two output signals and the time‐delay between them is evaluated. The proposed algorithms are verified with simulation data and recorded seismic response data from multi‐story buildings. The influence of noise on the time‐delay estimates is also assessed. Copyright © 2004 John Wiley & Sons, Ltd.     

Joint inversion of surface wave dispersion curves and reflection travel times via multi-objective evolutionary algorithms

Due to the character of seismic energy generation and propagation, shallow high-resolution seismic-reflection surveys often fail in the identification of the shallowest horizons and, due to the limited offsets, accuracy of velocity analyses is often not very high.In recent years, Rayleigh wave dispersion analysis have proved to have good potential also for near-surface applications but dispersion curve inversion and related uncertainty evaluation pose serious problems to a completely stand-alone application.In order to overcome these problems a joint inversion scheme is proposed, which is based on the identification of the Pareto front, performed in the framework of a Multi-Objective Evolutionary Algorithm (MOEA). Seismic data considered to design the two objectives are the Rayleigh wave dispersion curve and reflection travel times.We initially analyse a set of synthetic cases and evaluate the obtained results. A significant improvement of the retrieved models is observed as long as reflection travel times are added to the dispersion curve alone.Furthermore, the proposed methodology also provides relevant indications about the consistency of the overall inversion process. In fact, the distribution of the models in the objective space, the trend of the objectives over the passing generations and the evolution of the Pareto front can provide useful information to evaluate the provisional tentative interpretation (number of strata and reflector identification) inherently adopted for the data inversion.On the basis of the results obtained from the tests on the synthetic datasets, the analyses of a field dataset are interpreted as possible evidence of lateral heterogeneities.     

Automatic calibration of HEC‐HMS using single‐objective and multi‐objective PSO algorithms

This study presents single‐objective and multi‐objective particle swarm optimization (PSO) algorithms for automatic calibration of Hydrologic Engineering Center‐ Hydrologic Modeling Systems rainfall‐runoff model of Tamar Sub‐basin of Gorganroud River Basin in north of Iran. Three flood events were used for calibration and one for verification. Four performance criteria (objective functions) were considered in multi‐objective calibration where different combinations of objective functions were examined. For comparison purposes, a fuzzy set‐based approach was used to determine the best compromise solutions from the Pareto fronts obtained by multi‐objective PSO. The candidate parameter sets determined from different single‐objective and multi‐objective calibration scenarios were tested against the fourth event in the verification stage, where the initial abstraction parameters were recalibrated. A step‐by‐step screening procedure was used in this stage while evaluating and comparing the candidate parameter sets, which resulted in a few promising sets that performed well with respect to at least three of four performance criteria. The promising sets were all from the multi‐objective calibration scenarios which revealed the outperformance of the multi‐objective calibration on the single‐objective one. However, the results indicated that an increase of the number of objective functions did not necessarily lead to a better performance as the results of bi‐objective function calibration with a proper combination of objective functions performed as satisfactorily as those of triple‐objective function calibration. This is important because handling multi‐objective optimization with an increased number of objective functions is challenging especially from a computational point of view. Copyright © 2012 John Wiley & Sons, Ltd.     

Approximate modal decomposition of inelastic dynamic responses of wall buildings

Two approximate methods for decomposing complicated inelastic dynamic responses of wall buildings into simple modal responses are presented. Both methods are based on the equivalent linear concept, where a non‐linear structure is represented by a set of equivalent linear models. One linear model is used for representing only one vibration mode of the non‐linear structure, and its equivalent linear parameters are identified from the inelastic response time histories by using a numerical optimizer. Several theoretical relations essential for the modal decomposition are derived under the framework of complex modal analysis. Various numerical examinations have been carried out to check the validity of the proposed modal decomposition methods, and the results are quite satisfactory in all cases. Fluctuating bending moment and shear at any location along the wall height contributed by each individual vibration mode can be obtained. Modal contributions to shear and flexural strength demands, as well as the corresponding modal properties, under various seismic loading conditions can also be identified and examined in detail. Furthermore, the effects of higher vibration modes on seismic demands of wall buildings are investigated by using the modal decomposition methods. Several new insights into the complicated inelastic dynamics of multi‐story wall buildings are presented. Copyright © 2004 John Wiley & Sons, Ltd.     

Assessment of the adequacy of the regional relationships between catchment attributes and catchment response dynamics,calibrated by a multi‐objective approach

This study aimed to evaluate the effectiveness of the regionalization method on the basis of a combination of a parsimonious model structure and a multi‐objective calibration technique. For this study, 12 gauged catchments in the Republic of Korea were used. The parsimonious model structure, requiring minimal input data, was used to avoid adverse effects arising from model complexity, over‐parameterization and data requirements. The IHACRES rainfall‐runoff model was applied to represent the dynamic response characteristics of catchments in Korea. A multi‐objective approach was adopted to reduce the predictive uncertainty arising from the calibration of a rainfall‐runoff model, by increasing the amount of information retrieved from the available data. The regional relationships (or models) between the model parameters and the catchment attributes were established via a multiple regression approach, incorporating correlation analysis and stepwise regression on linear and logarithmic scales. The impacts of the parameters, calibrated by the multi‐objective approach, on the adequacy of regional relationships were assessed by comparison with impacts obtained by the single‐objective approach. The regional relationships were well defined, despite limited available data. The drainage area, the effective soil depth, the mean catchment slope and the catchment gradient appeared to be the main factors for describing the hydrologic response characteristics in the areas studied. The overall model performance of the regional models based on the multi‐objective approach was good, producing reasonable results for high and low flows and for the overall water balance, simultaneously. The regional models based on the single‐objective approach yielded accurate predictions in high flows but showed limited predictive capability for low flows and the overall water balance. This was due to the optimal model parameter estimates when using a single‐objective measure. The parameters calibrated by the single‐objective approach decreased the predictability of the regional models. Copyright © 2013 John Wiley & Sons, Ltd.     

System identification of linear structures based on Hilbert–Huang spectral analysis. Part 1: normal modes

Based on the Hilbert–Huang spectral analysis, a method is proposed to identify multi‐degree‐of‐freedom (MDOF) linear systems using measured free vibration time histories. For MDOF systems, the normal modes have been assumed to exist. In this method, the measured response data, which are polluted by noises, are first decomposed into modal responses using the empirical mode decomposition (EMD) approach with intermittency criteria. Then, the Hilbert transform is applied to each modal response to obtain the instantaneous amplitude and phase angle time histories. A linear least‐square fit procedure is proposed to identify the natural frequency and damping ratio from the instantaneous amplitude and phase angle for each modal response. Based on a single measurement of the free vibration time history at one appropriate location, natural frequencies and damping ratios can be identified. When the responses at all degrees of freedom are measured, the mode shapes and the physical mass, damping and stiffness matrices of the structure can be determined. The applications of the proposed method are illustrated using three linear systems with different dynamic characteristics. Numerical simulation results demonstrate that the proposed system identification method yields quite accurate results, and it offers a new and effective tool for the system identification of linear structures in which normal modes exist. Copyright © 2003 John Wiley & Sons, Ltd.     

Multireservoir system optimization in the Han River basin using multi‐objective genetic algorithms

In this study, NSGA‐II is applied to multireservoir system optimization. Here, a four‐dimensional multireservoir system in the Han River basin was formulated. Two objective functions and three cases having different constraint conditions are used to achieve nondominated solutions. NSGA‐II effectively determines these solutions without being subject to any user‐defined penalty function, as it is applied to a multireservoir system optimization having a number of constraints (here, 246), multi‐objectives, and infeasible initial solutions. Most research by multi‐objective genetic algorithms only reveals a trade‐off in the objective function space present, and thus the decision maker must reanalyse this trade‐off relationship in order to obtain information on the decision variable. Contrastingly, this study suggests a method for identifying the best solutions among the nondominated ones by analysing the relation between objective function values and decision variables. Our conclusions demonstrated that NSGA‐II performs well in multireservoir system optimization having multi‐objectives. Copyright © 2005 John Wiley & Sons, Ltd.     

Bridge health monitoring system based on vibration measurements

A bridge health monitoring system is presented based on vibration measurements collected from a network of acceleration sensors. Sophisticated structural identification methods, combining information from the sensor network with the theoretical information built into a finite element model for simulating bridge behavior, are incorporated into the system in order to monitor structural condition, track structural changes and identify the location, type and extent of damage. This work starts with a brief overview of the modal and model identification algorithms and software incorporated into the monitoring system and then presents details on a Bayesian inference framework for the identification of the location and the severity of damage using measured modal characteristics. The methodology for damage detection combines the information contained in a set of measurement modal data with the information provided by a family of competitive, parameterized, finite element model classes simulating plausible damage scenarios in the structure. The effectiveness of the damage detection algorithm is demonstrated and validated using simulated modal data from an instrumented R/C bridge of the Egnatia Odos motorway, as well as using experimental vibration data from a laboratory small-scaled bridge section.