Output-Only Modal Parameter Recursive Estimation of Time-Varying Structures via a Kernel Ridge Regression FS-TARMA Approach

Modal tracking plays a vital role in structural health monitoring since changes in modal parameters help us understand a structure’s dynamic characteristics and may reflect the potential deterioration of structural performance. Although numerous modal parameter estimation (MPE) methods exist, it is not guaranteed that an MPE process will exclude all spurious modes and not lose any physical modes every time over a long-term monitoring period. Relatively large damping of a structure, poor data quality, and significant changes in structural modal parameters may make the estimated modal parameters spurious, missing, or misclassified. It makes long-term modal tracking semiautomated or manual, which constrains timely downstream applications such as anomaly detection, condition assessment, and decision making. This research aims to propose a long-term continuous automatic modal tracking algorithm based on Bayesian inference even when the modal parameters, damping, and data quality change significantly. Bayesian inference is used to determine the physical modes from the results of existing MPE methods. Both the modes identified from the most recent response set and the modal probability model from multiple previous response sets are considered in the Bayesian model to better determine the physical modes from the results of MPE. Moreover, the proposed algorithm requires only three extra hyperparameters compared to general modal tracking algorithms, and they can be quickly determined by a grid search method. The performance of the proposed algorithm is verified by a numerical example and a real-world civil structure Z24 Bridge benchmark.

The accuracy of modal parameter estimates from estimation methods built in the framework of experimental modal analysis (EMA) can be assessed by several methods provided the information of input loadings and output responses. However, a deeper investigation is necessary for estimation methods built in the operational modal analysis (OMA) framework to establish confidence boundaries for the modal parameter estimates. This paper aims to estimate uncertainties of the modal parameter estimates by testing the sensitivity of the damped natural frequencies, damping ratios, and mode shapes to the choice of the estimation method and data acquisition settings used in this method. The uncertainty estimates of the modal parameter estimates are tested for changes in the boundary conditions. The tests provide insights into tendencies for the accuracy of the modal parameter estimates when the boundary conditions provide a different set of “true” modal parameter values. A simple cantilever beam experimental setup is used to obtain benchmark modal parameter estimates from EMA that are used to estimate the accuracy of the modal parameter estimates from OMA. The method’s robustness in estimating uncertainties and accuracies for modal parameter estimates is tested for several independent experiments. The boundary conditions for the cantilever beam are changed to test the influence of altering modal parameter estimates on the estimated uncertainties and accuracies. The boundary conditions are changed using materials with different elastic properties and the connection pressure between the clamping device and beam. A method presented in Paulsen et al. (Proceedings of the 9th International Operational Modal Analysis Conference, IOMAC, 2022) calculates a representative modal parameter estimate from a set of modal parameter estimates based on weightings obtained from the modal assurance criterion. The method’s ability to reduce uncertainties and increase accuracies of the representative estimate is demonstrated by repeating a test setup ten times. The results present simultaneously a reduction of uncertainties and an increase of accuracies for the representative modal parameter when the weighting method is used. An evaluation of dependencies between the estimated uncertainties and accuracies finds a proportional relation for the damped natural frequency when the uncertainties are described as a relative standard deviation. The analysis presents higher accuracy for the damped natural frequency and mode shape estimates compared to the damping factor, and the accuracy of the damping factor is dependent on the actual value of the damping factor. These phenomena for accuracy are well known and indicate that the uncertainties of modal parameter estimates can be assessed through a sensitivity test of estimation methods and data acquisition parameters.

A novel recursive modal parameter estimator for operational time-varying structural dynamic systems based on least squares support vector machine and time series model

ABSTRACTThe structural stiffness of a blade changes when it is defective, which affects its modal parameters. Consequently, modal parameters can be crucial to investigating defects. Therefore, developing a model that can accurately analyze the modal parameters of healthy and defective blades is essential. This study proposes a high‐fidelity finite element (FE) model development method based on multiobjective optimization. Three 10‐kW‐class wind turbine blades were manufactured, and the modal parameter changes due to debonding damage were measured. The initial FE model was constructed using structural drawings and considering specific material properties. Subsequently, the modal analysis results were compared with the experimental results. The analysis errors of normal and abnormal blades were expressed as cost functions and multiobjective optimization was performed. Machine learning–based predictive models facilitated time‐ and cost‐effective optimization, enabling modal updating that simultaneously reflects modal parameter changes for various blade models. Because accurately predicting the dynamic characteristics of blades by various defect locations and sizes is crucial for blade damage detection, the proposed method is significant for damage detection research.

This article is an experimental study of the effects of multi-site damage on the vibration response of a composite beam damaged by low-energy impact. The variation of the modal parameters with different levels of impact energy and density of impact is studied. Specimens are impacted symmetrically in order to induce a global rate of damage. A damage detection tool Damage Index is introduced in order to verify the estimation of damping ratios. Design of Experiments is used to establish the sensitivity of both energy of impact and density of damage. The DOE analysis results (using natural frequency only) indicate that impact energy for 2nd, 3rd, and 4th bending modes is the most significant factor contributing to the changes in the modal parameters for this kind of symmetrical dynamic test.

Modal analysis has developed into a major technology for the study of structural dynamics in the past several decades. Through it, complex structural dynamics phenomena can be represented in terms of structural invariants, i.e., the modal parameters: natural frequencies, damping ratios and mode shapes. Operational Modal Analysis (OMA) deals with the estimation of modal parameters on vibration data measured for operational conditions, when the excitation on the structure is not measured. In this work, OMA is performed on a wind turbine blade undergoing wind tunnel testing. These tests included different wind speeds, pitch angles and also different health conditions of the blade, where masses of different magnitude were fixed to the blade on different locations to emulate damage conditions. In order to monitor the modal parameters across multiple days of varied tests in the wind tunnel, the Polymax modal parameter estimator was implemented, coupled with an Automated Modal Analysis methodology. This methodology included an automatic modal parameter selection technique, using a Machine Learning (ML) clustering algorithm, coupled with a modal tracking procedure which applied statistical thresholds on the modal parameters’ values. The tracking procedure searches for modes similar to the ones calculated for healthy conditions. The results show how the modal parameters of the wind turbine blade vary with the different measured conditions in the wind tunnel. Moreover, a damage detection methodology is implemented to differentiate between the healthy and damaged conditions on the blade, by leveraging an anomaly detection algorithm using the Multivariate Gaussian Distribution (MGD). This algorithm takes as input the modal parameters calculated by the previous Automated Modal Analysis methodology and detects statistical deviations among them which could indicate the presence of damage. All steps of this work contribute to developing an automatic framework able to detect damages on a wind turbine blade, and therefore perform Structural Health Monitoring (SHM) for different operational conditions.

This paper presents an experimental study on the detection and localization of damage in a steel plate using vibration-based methods. Dynamic impulse tests were conducted on both undamaged and damaged plate to determine changes in modal parameters, such as eigenfrequencies and mode shapes. Several damage detection methods, including the Modal Assurance Criterion (MAC), Coordinate Modal Assurance Criterion (COMAC), Mode Shape Curvature (MSC), Mode Shape Curvature Square (MSCS), and Damage Index (DI), were applied. The results showed that while MAC and COMAC effectively detected damage, they did not provide information about the location. MSC and MSCS indices demonstrated improved damage localization, with the MSCS being particularly effective. The DI index provided the clearest and most accurate representation of the damage location. These findings highlight the potential of these methods for non-destructive testing and structural health monitoring.

Recent work with autonomous modal parameter estimation has shown great promise in the quality of the modal parameter estimation results when compared to results from traditional methods by experienced users. While autonomous modal parameter estimation means slightly different things to different researchers and practitioners, for the purpose of this discussion, autonomous will require an automated procedure which sorts and processes a large number of possible modal parameter solutions to yield one consistent estimate with no user interaction after initial thresholds are chosen. In the work discussed, this final, consistent set of modal parameters is identifiable due to the combination of temporal and spatial information in azdomain state vector of relatively high order (5–10). Since thiszdomain state vector has both complex modal frequency and modal vector information as embedded content, sorting consistent estimates from the multitude of possible solutions is relatively trivial. Because thiszdomain state vector can be developed from the results of any modal parameter estimation method, possible solutions from different traditional methods can be utilized in the autonomous procedure to yield one consistent set of modal parameters.

Modal parameters estimation plays an important role for structural health monitoring. Accurately estimating the modal parameters of structures is more challenging as the measured vibration response signals are contaminated with noise. This study develops a mathematical algorithm of solving the partially described inverse singular value problem (PDISVP) combined with the complex exponential (CE) method to estimate the modal parameters. The PDISVP solving method is to reconstruct an L2-norm optimized (filtered) data matrix from the measured (noisy) data matrix, when the prescribed data constraints are one or several sets of singular triplets of the matrix. The measured data matrix is Hankel structured, which is constructed based on the measured impulse response function (IRF). The reconstructed matrix must maintain the Hankel structure, and be lowered in rank as well. Once the filtered IRF is obtained, the CE method can be applied to extract the modal parameters. Two physical experiments, including a steel cantilever beam with 10 accelerometers mounted, and a steel plate with 30 accelerometers mounted, excited by an impulsive load, respectively, are investigated to test the applicability of the proposed scheme. In addition, the consistency diagram is proposed to exam the agreement among the modal parameters estimated from those different accelerometers. Results indicate that the PDISVP-CE method can significantly remove noise from measured signals and accurately estimate the modal frequencies and damping ratios.

Modal-based structural health monitoring (SHM) detects damage and degradation phenomena from the variations of the modal parameters over time. However, the modal parameter estimates are also influenced by environmental and operational variables (EOVs) whose effects have to be compensated. Modeling the influence of EOVs on modal parameters is very challenging, so black-box models, such as regression models, are often adopted as an alternative. However, in many applications, the set of measured EOVs is incomplete or the factors influencing the estimates cannot be identified or measured. In these conditions, output-only techniques for compensation of environmental and operational effects are an attractive alternative. Different linear as well as nonlinear methods for the compensation of the environmental and operational influence on modal parameters in the context of modal-based SHM are reviewed in the present paper. Real datasets collected from vibration-based monitoring systems are analyzed, and the results are presented and discussed to illustrate the applicative perspectives and possible drawbacks of the selected methods.KeywordsModal-based SHMEnvironmental influenceRegression modelPrincipal component analysisKernel PCASecond-order blind identification

A multivariate maximum likelihood procedure for the estimation of modal parameters is presented. The vibrating system is excited by a random force and sensors output only are used to estimate the natural frequencies and damping coefficients of the system. The method works in time domain and a vector autoregressive moving average (VARMA) process is used. The information about modal parameters is contained in the multivariate AR part, which is estimated using an iterative maximum likelihood algorithm. This algorithm uses a score technique and output data only. The order of the AR part is obtained via the Minimum Description Length associated with an instrumental variable procedure. Experimental results show the effectiveness of the method for model order selection and modal parameters estimation.

The aim of the present paper is to evaluate structural health monitoring (SHM) techniques based on modal analysis for crack detection in small wind turbine blades. A finite element (FE) model calibrated to measured modal parameters will be introduced to cracks with different sizes along one edge of the blade. Changes in modal parameters from the FE model are compared with data obtained from experimental tests. These comparisons will be used to validate the FE model and subsequently discuss the usability of SHM techniques based on modal parameters for condition monitoring of wind turbine blades.

Bladed disks are critical parts of aeronautic turbojet engines requiring numerous dynamic analysis tests to fully achieve their design process. This work presents subspace state-space identification techniques for estimation of modal parameters of bladed disks using realistic although numerically generated test case data. Indeed, modal testing of bladed disks can exhibit high modal density leading to modal estimation issues. This study focuses on subspace state-space identification framework as it is an efficient way to process vector valued time series and thus to estimate modal parameters in a context of high modal density. In order to evaluate the techniques of interest, a versatile modal model has been specifically implemented to simulate data representative of full scale rotating aeronautic fan modal tests. This model is easily adjustable allowing to evaluate the identification methods in a more or less severe estimation context. The performance of the investigated methods are discussed and compared with the prescribed parameters of the model. Moreover, different techniques to estimate the order of the associated state-space model are reviewed and tested over several simulated configurations. Finally, a method to evaluate the uncertainties over the modal parameters using covariance of estimator results and pseudospectrum computation is proposed and discussed over the investigated test cases.

In order to identify time-varying transient modal parameters only from nonstationary vibration response measurement signals for slow linear time-varying (SLTV) structures which are weakly damped, a moving window differential evolution (DE) independent component analysis- (ICA-) based operational modal analysis (OMA) method is proposed in this paper. Firstly, in order to overcome the problems in traditional ICA-based OMA, such as easy to go into local optima and difficult-to-identify high-order modal parameters, we combine DE with ICA and propose a differential evolution independent component analysis- (DEICA-) based OMA method for linear time invariant (LTI) structures. Secondly, we combine the moving widow technique with DEICA and propose a moving window differential evolution independent component analysis- (MWDEICA-) based OMA method for SLTV structures. The MWDEICA-based OMA method has high global searching ability, robustness, and complexity of time and space. The modal identification results in a three-degree-of-freedom structure with slow time-varying mass show that this MWDEICA-based OMA method can identify transient time-varying modal parameters effectively only from nonstationary vibration response measurement signals and has better performances than moving window traditional ICA-based OMA.

A two-stage, time-domain structural dynamics identification method, that circumvents some of the main drawbacks of current techniques, and combines good noise rejection properties with a low computational load and accurate modal parameter extraction, is proposed. In the first stage of the method a discrete model is obtained through a linear estimation scheme and a systematic and effective modeling strategy. In the second stage the modal parameters are evaluated according to a new procedure that is capable of detecting and correcting large errors that may be introduced into the modal parameter estimates through the discrete-to-continuous transformation. The proposed method yields consistent modal parameter estimates at a low computational load, and is well-suited for microcomputer implementation. Its effectiveness is finally demonstrated with both simulated and actual experimental data.