Parametric Output Error Based Identification and Fault Detection in Structures Under Earthquake Excitation
|
|
- Louisa Maxwell
- 6 years ago
- Views:
Transcription
1 Parametric Output Error Based Identification and Fault Detection in Structures Under Earthquake Excitation J.S. Sakellariou and S.D. Fassois Department of Mechanical & Aeronautical Engr. GR 265 Patras, Greece June 2
2 Output Error Identification And Fault Detection In Civil Structures 1 TALK OUTLINE 1. Introduction 2. The Lumped-Parameter Building Model 3. Output Error Based Structural Identification 4. Fault Detection and Identification 5. Concluding Remarks
3 Output Error Identification And Fault Detection In Civil Structures 2 1. Introduction The General Problem Dynamic identification and fault detection and identification of civil structures under earthquake excitation. Problem Significance The structural identification and fault detection & identification under earthquake excitation is important for: dynamic analysis - modal parameter extraction structural health monitoring - fault diagnosis prediction & control Note that earthquake excitation is attractive for large civil structures. Overall aim: The assessment & improvement of the dynamical characteristics and safety of civil structures
4 Output Error Identification And Fault Detection In Civil Structures 3 Problem Characteristics The earthquake excitation and its characteristics introduce a number of technical difficulties due to: its transient nature its limited duration (limited data records) its limited frequency content its non-stationarity (in terms of both amplitude and frequency) & the presence of non-stationary noise its strength may lead the structure into a non-linear operating regime Coping with these characteristics requires caution & and proper selections.
5 Output Error Identification And Fault Detection In Civil Structures 4 Literature Survey A. Structural Identification Classical frequency domain (spectral) methods are not applicable to the non-stationary stochastic case. Evolutionary spectral density type methods (Staszewski, 1998). ARX & ARMAX time domain methods are not applicable. Adaptive methods (time-domain based on the EKF or alternative recursive schemes) are applicable but limited to slowly-varying nonstationarities and suboptimal accuracy (Ghanem and Shinozuka, 1995; Popescu, 1995; Loh and Tou, 1995; Kirkegaard et al., 1996).
6 Output Error Identification And Fault Detection In Civil Structures 5 B. Fault Detection & Identification Methods based on changes in the structure s stiffness matrix using FEM updating and dynamic/static test data (Hajela and Soeiro, 1989; Hearn and Testa, 1991; Capecchi and Vestroni, 1999;). (difficulties: static testing; complete eigenmode information requires extensive instrumentation/testing; a stiffness element may have contributions from several structural elements; stochastic effects often unaccounted for). Methods based on changes in the structure s modal characteristics (Yao, 1988; Salawu, 1997). (difficulties: most require complete modal information; extensive instrumentation/testing; fault magnitude estimation problematic; fault identification may be impossible from pattern changes; stochastic effects often unaccounted for). Methods based on damage indices (Dipasquale and Çakmak 199; Wiliams and Sexsmith, 1995; Garcia et al., 1999). (difficulties: measurement in every storey; extensive instrumentation; ambiguous correlation with the actual fault).
7 Output Error Identification And Fault Detection In Civil Structures 6 Study Objectives The formulation of an effective Output Error (OE) based method for the dynamic identification and fault detection of structures under earthquake excitation. Intended characteristics: Suitable for non-stationary excitation and measurement noise. Achieving high accuracy & resolution with limited data records. Effective fault detection & identification in a stochastic environment. Effectiveness with a single transfer function measurement (two sensors). Capability for fault magnitude estimation.
8 Output Error Identification And Fault Detection In Civil Structures 7 2. The Lumped-Parameter Building Model Figure 1: Schematic diagram of the 6 storey building. m i = kg, k i = KN/m, c i = Ns/m.
9 Output Error Identification And Fault Detection In Civil Structures 8 Theoretical building characteristics (Ẍ6/X transfer func- Table 1: tion). f n (Hz) ζ(%) Residue j j j j j j.824 The building longitudinal dynamics are described by the differential equation: M ẍ + C ẋ + K x = Q x x = [x 1 x 2 x 3 x 4 x 5 x 6 ] T : longitudinal storey displacement vector x = [x ẋ ] T : input vector (x : longitudinal ground displacement) Q : input shaping matrix The above model is reformulated in state space form as follows: ż = A z + B x o y = C z
10 Output Error Identification And Fault Detection In Civil Structures 9 z = [z 1... z 12 ] T : state vector, y = [ ψ 1 ψ 1... ψ 6 ψ ] T 6 A = 1 k 1+k 2 m 1 c 1+c 2 m 1 k 2 m 1 c 2 m 1 1 k 2 m 2 c 2 m 2 k 2+k 3 m 2 c 2+c 3 m 2 k 3 m 2 c 3 m 2 1 k 3 m 3 c 3 m 3 k 3+k 4 m 3 c 3+c 4 m 3 k 4 m 3 c 4 m 3 1 k 4 m 4 c 4 m 4 k 4+k 5 m 4 c 4+c 5 m 4 k 5 m 4 c 5 m 4 1 k 5 m 5 c 5 m 5 k 5+k 6 m 5 c 5+c 6 m 5 k 6 m 5 c 6 m 5 1 k 6 m 6 c 6 m 6 k 6 m 6 c 6 m 6 B = [ 1 ]T C = I(12) The expression for the 6 th floor acceleration is: ẍ 6 [t] = k 1 m 1 ψ 6 + c m 1... ψ 6
11 Output Error Identification And Fault Detection In Civil Structures 1 3. Output Error Based Structural Identification Earthquake Excitation Earthquake recorded in Patras, Greece, in Non-stationary behaviour - Energy concentration in low frequencies ( - 1 Hz) Duration of strong excitation 6.52 sec Maximum acceleration.19g Sampling rate f s = 2Hz Recording time 37.2 sec N = 748 samples Frequency (Hz) (a) Displacement (m) 6 x (b) Time (sec) Figure 2: Seismic displacement (a) and its spectrogram (b) (f s = 25Hz).
12 Output Error Identification And Fault Detection In Civil Structures 11 Modeling of the earthquake displacement Adaptive Recursive AutoRegressive Moving Average (RARMA) modeling based upon the Recursive Maximum Likelihood (RML) algorithm is used along with non-parametrically estimated innovations variance (sliding non-causal window) (Fouskitakis and Fassois, 2). A RARMA(3,8) model with forgetting factor λ =.975 is obtained. RARMA(3,8): x[t] 3 i=1 φ i [t] x[t i] = a[t] 8 i=1 θ i [t] α[t i] x[t] : ground displacement a[t] : innovations (uncorrelated) sequence E{a[t]} = E{a 2 [t]} = σα[t] 2 φ i [t] : time-varying AR parameter θ i [t] : time-varying MA parameter This RARMA(3,8) model is used for the generation of artificial ground displacement signals to be used in the Monte Carlo experiments.
13 Output Error Identification And Fault Detection In Civil Structures 12 Simulation of the Building Dynamics Integration via step invariance discretization at f s = 2Hz Low pass filtering (f c = 1Hz) of the excitation and resulting responses. Signal resampling at f s = 25Hz. The resulting signals are N = 926 samples long. The response (6 th storey acceleration) is corrupted with non-stationary uncorrelated noise at the 5% N/S ratio. Frequency (Hz) (a) Acceleration (m/s 2 ) Time (sec) Figure 3: 6 th floor acelleration (a) and its spectrogram (b) (nominal structure). (b)
14 Output Error Identification And Fault Detection In Civil Structures 13 The Output Error model Output Error (OE) identification aims at the estimation of a discrete-time structural model of the form, OE(na, nb) model : ẍ 6 [t] = B(B) A(B) x [t] + e[t] x [t]: ground displacement ẍ 6 [t]: noise-corrupted 6th storey acceleration e[t]: Output error / non-stationary disturbance B: Backshift operator (B x[t] = x[t 1]) A(B) = 1 + a 1 B a na B n a : AutoRegressive (AR) polynomial B(B) = b + b 1 B b nb B n b : EXogenous (X) polynomial n a : AR order n b : X order The main advantage of the OE model is that it is capable of accounting for noise effects without resorting on an explicit noise representation overcoming problems with non-stationary measurement noise.
15 Output Error Identification And Fault Detection In Civil Structures 14 Parameter Estimation Parameter estimation is based upon minimization of a quadratic functional of the Output Error. Thus: ˆθ N = arg min θ J(θ) = arg min θ 1 N N 1 t= e 2 [t] θ = [ a 1 a 2... a na. b b 1 b 2... b nb ] T Optimization is based upon the Levenberg-Marquardt scheme (Ljung, 87). The estimate is asymptotically Gaussian with mean equal to the true parameter vector and specific covariance: θ o : true parameter vector N(ˆθN θ o ) N (,P) (N) Model Order Selection Model order selection is based on the successive estimation of increasingly higher order models and: Examination of the achieved Output Error cost J(ˆθ) Modal Dispersion Analysis (Lee and Fassois, 1993) provides the energy contribution of each mode in the measured signal unnecessary (overfitted) modes are detected by their small dispersions
16 Output Error Identification And Fault Detection In Civil Structures 15 System Identification Results The identification procedure leads to an OE(6,6) model (b ) with J = Output Error Criterion [6,5] [6,6] [6,7] [6,8] [6,9] [7,5] [7,6] [7,7] [7,8] [7,9] [8,5] [8,5] [8,5] [8,5] Output Error model orders (na,nb) [8,5] [9,5] [9,6] [9,7] [9,8] [9,9] Figure 4: structure). Output Error for some Output error models (nominal
17 Output Error Identification And Fault Detection In Civil Structures 16 Monte Carlo Identification experiments (nominal structure) 15 2 Magnitude (db) Phase (deg) Frequency (Hz) Frequency (Hz) Figure 5: Bode diagrams of the estimated OE(6,6) models ( ) and the theoretical structure (- - -) (Ẍ6/X transfer function; 2 Monte Carlo experiments; nominal structure). Mode Natural Frequencies (Hz) Damping factors (%) Dispersions (%) (2.18 ± ) 1.21 (1.21 ± ) 5.1 (5.91 ±.23) (6.41 ± ) 3.56 (3.56 ± ) (39.12 ± 1.62) (1.26 ± ) 5.7 (4.6 ± ) 6.95 (54.97 ± 1.85) True (Estimated ± Standard Deviation)
18 Output Error Identification And Fault Detection In Civil Structures Fault Detection and Identification Structural Faults The faults considered represent reductions in the elasticity characteristics of each storey. A fault mode is defined as the continuum of faults of all possible magnitudes occurring in a particular storey. Six fault modes are presently defined: F 1 α, F 2 α,..., F 6 α α: fault magnitude (stiffness reduction) Several structural health monitoring Monte Carlo experiments are considered corresponding to: F (no fault) F 3.5 (5% reduction in K 3 ) F 1.5 (5% reduction in K 1 ) F 3.2 (2% reduction in K 3 ) F 1.2 (2% reduction in K 1 ) The effects of the faults on the considered ẍ 6 [t]/x [t] transfer function are minimal.
19 Output Error Identification And Fault Detection In Civil Structures Magnitude (db) 1 2 Phase (deg) Frequency (Hz) Frequency (Hz) Figure 6: Theoretical Bode diagram of the nominal structure ( ), the structure under F 3.5 (- - -), and the structure under F 3.2 (...) (Ẍ6/X transfer function). Acceleration (m/s 2 ) Acceleration (m/s 2 ) Acceleration (m/s 2 ) Time (sec) Figure 7: The 6 th floor acceleration in the nominal (a) and damaged (k 3 reduced by 5% and 2%) (b), (c) cases. (a) (b) (c)
20 Output Error Identification And Fault Detection In Civil Structures 19 Fault Detection Methodology Fault detection is based upon comparison of the interval estimator of the model parameter vector in the nominal and current states of the structure. Let C i represent the confidence interval (at the α =.5 level) of the i-th scalar parameter θ i : C i = (ˆθi 1.96ˆσ θi, ˆθ ) i ˆσ θi The structure is declared as healthy if all confidence intervals of the current parameter vector are not disjoint with those of the nominal parameter vector. Otherwise the structure is declared as damaged. C i C u i = { } for some i otherwise Fault detected No fault detected C i C u i : conf. interval of θ i (nominal state of the structure) : conf. interval of θ i (current state of the structure) { } : empty set C i C i C u i C u i No fault Fault
21 Output Error Identification And Fault Detection In Civil Structures 2 Fault Identification Methodology Once a fault has been detected, fault identification (localization) and magnitude estimation are accomplished via the Geometric Approach (Sadeghi & Fassois, 1998). This consists of seven stages: Stage 1: An initial feature vector θ is selected, e.g. θ = [AR parameters] T Stage 2: The feature vector is transformed into a coordinate system in which information compression may be best achieved: s = U T θ (P = U Λ U T : covariance) Stage 3: A reduced, transformed, feature vector s M is obtained by specifying the allowable information loss. θ = M s j u j + N s j u j j=1 j=m+1 }{{}}{{} reduced error feature vector = [U M.U N M ] s M... s N M Information is expressed in terms of logarithmic entropy: H R (s) = N ν j = Var[s j ]/ N j=1 j=1 ν j log 2 ( ν j ) Var[s j ] (j = 1, 2,..., N)
22 Output Error Identification And Fault Detection In Civil Structures 21 Stage 4: A stochastic feature space, spanned by the 1st & 2nd order moments of the transformed & reduced feature vector s M, is constructed. θ = [ µ T s M, Cov[s M ] ] T (ρ dimensional) Stage 5: Fault mode geometric representations (hyperplanes). A geometric representation of each fault mode Fα i (i = 1, 2,..., 6) within the stochastic feature space is constructed with the aid of experiments performed with a simulation model. The i-th fault mode is thus represented by the hyperplane equation: g i (θ K ) = θ K1 + ω1 i θ K ωρ 1 i θ Kρ ωρ i = (ω1, i ω2, i..., ωρ) i estimated via linear regression. Stage 6: Fault mode identification (localization). The vector θ u K corresponding to the current fault is estimated. Its distances from the fault mode hyperplanes are estimated by optimizing the Lagrangian: L(θ K, γ) = D(θ K, θ u K ) + 2γ gi (θ K ) D(θ K, θ u K ) : squared distance function The current fault θ u K is then identified as belonging to the fault mode with the hyperplane of which the distance is minimal.
23 Output Error Identification And Fault Detection In Civil Structures 22 Stage 7: Fault magnitude estimation. A measure of the fault magnitude is obtained by evaluating the distance between the current fault and the nominal system: α D 1/2 (θk o, θu K ) Figure 8: Principle of the geometric approach (θk u : current fault; g i (θ K ) = : i-th fault mode subspace; θk o : nominal-unfailed system)
24 Output Error Identification And Fault Detection In Civil Structures 23 Fault Detection and Identification Results Two test cases: F 3.5 and F Confidence Intervals a 1 a 2 a 3 a 4 a 5 a 6 AR Parameters Figure 9: Fault detection results for the F 3.5 (...) and F 3.2 (- - -) faults (the confidence intervals of the parameters of the nominal structure are indicated by solid line; single experiments).
25 Output Error Identification And Fault Detection In Civil Structures 24 Initial feature vector: θ = [α 1 α 2... α 6 ] T Transformed & reduced feature vector: s = [s 1 s 2 s 3 ] T θ K = [µ s1 µ s2 µ s3 σ 2 s 1 σ 2 s 2 σ 2 s 3 ] T Entropy percentage Feature Vector Dimension Figure 1: Normalized logarithmic entropy of the transformed feature vectors for every damage level, versus their dimension (F 3 case).
26 Output Error Identification And Fault Detection In Civil Structures 25 injected faults F 3 α α = Parameter θ k Parameter θ k Parameter θ k1 1.8 (a) Distance %.5 5% Magnitude (%) (b) Figure 11: (a) 3-dimensional intersection of the F 3 hyperplane along with projections (+) of various fault points; (b) the distance function versus fault magnitude for the F 3 fault mode (two examples of fault magnitude estimation for the F 3.5 and F 3.2 faults are indicated).
27 Output Error Identification And Fault Detection In Civil Structures 26 9 x 1 3 Fault F (a) Distance minimal distance F 1 F 2 F 3 F 4 F 5 F 6 Distance Fault F 3.2 minimal distance F 1 F 2 F 3 F 4 F 5 F 6 Fault Modes (b) Figure 12: Distances of the current fault point from each hyperplane: (a) Current fault F 3.5; (b) current fault F 3.2 (single experiments).
28 Output Error Identification And Fault Detection In Civil Structures 27 Monte Carlo Experiments (2 realizations each) Structure Correct Detection False Detection No Fault 2 F F F F Fault Correct Identification False Identification F F F F 1.2 2
29 Output Error Identification And Fault Detection In Civil Structures Concluding Remarks A stochastic time-domain Output Error based method for the dynamic identification and fault detection/identification of structures under earthquake excitation was introduced. The method was shown to overcome difficulties associated with: the excitation s transient nature & limited duration (N=926 samples long data records were used) the excitation s limited frequency range (-1 Hz) the excitation s strong non-stationarity & the presence of nonstationary random measurement noise (N/S = 5%) the need for multiple measurement (only the ground and 6 th storey accelerations were measured) The method was shown to achieve: Statistically consistent and highly accurate estimates of the structure s modal parameters. Effective fault detection, identification and magnitude estimation in the presence of random non-stationary noise. Effectiveness with faults as small as 5 % stiffness reduction.
NON-STATIONARY MECHANICAL VIBRATION MODELING AND ANALYSIS
NON-STATIONARY MECHANICAL VIBRATION MODELING AND ANALYSIS VIA FUNCTIONAL SERIES TARMA MODELS A.G. Poulimenos and S.D. Fassois DEPARTMENT OF MECHANICAL &AERONAUTICAL ENGINEERING GR-26500 PATRAS, GREECE
More informationA. Poulimenos, M. Spiridonakos, and S. Fassois
PARAMETRIC TIME-DOMAIN METHODS FOR NON-STATIONARY RANDOM VIBRATION IDENTIFICATION AND ANALYSIS: AN OVERVIEW AND COMPARISON A. Poulimenos, M. Spiridonakos, and S. Fassois DEPARTMENT OF MECHANICAL & AERONAUTICAL
More informationVibration Based Health Monitoring for a Thin Aluminum Plate: Experimental Assessment of Several Statistical Time Series Methods
Vibration Based Health Monitoring for a Thin Aluminum Plate: Experimental Assessment of Several Statistical Time Series Methods Fotis P. Kopsaftopoulos and Spilios D. Fassois Stochastic Mechanical Systems
More informationStationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study
Stationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study Andriana S. GEORGANTOPOULOU & Spilios D. FASSOIS Stochastic Mechanical Systems & Automation
More informationStationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study
6th International Symposium on NDT in Aerospace, 12-14th November 2014, Madrid, Spain - www.ndt.net/app.aerondt2014 More Info at Open Access Database www.ndt.net/?id=16938 Stationary or Non-Stationary
More informationVector-dependent Functionally Pooled ARX Models for the Identification of Systems Under Multiple Operating Conditions
Preprints of the 16th IFAC Symposium on System Identification The International Federation of Automatic Control Vector-dependent Functionally Pooled ARX Models for the Identification of Systems Under Multiple
More informationT.-C.J. Aravanis, J.S. Sakellariou and S.D. Fassois
Vibration based fault detection under variable non-measurable, operating conditions via a stochastic Functional Model method and application to railway vehicle suspensions T.-C.J. Aravanis, J.S. Sakellariou
More information2330. A study on Gabor frame for estimating instantaneous dynamic characteristics of structures Wei-Chih Su 1, Chiung-Shiann Huang 2 1
2330. A study on Gabor frame for estimating instantaneous dynamic characteristics of structures Wei-Chih Su 1, Chiung-Shiann Huang 2 1 National Center for High-Performance Computing, National Applied Research
More informationIdentification of Stochastic Systems Under Multiple Operating Conditions: The Vector Dependent FP ARX Parametrization
Identification of Stochastic Systems Under Multiple Operating Conditions: The Vector Dependent FP ARX Parametrization Fotis P Kopsaftopoulos and Spilios D Fassois Abstract The problem of identifying stochastic
More informationthe Functional Model Based Method
Multi-Site Damage Localization via the Functional Model Based Method Christos S. Sakaris, John S. Sakellariou and Spilios D. Fassois Stochastic Mechanical Systems & Automation (SMSA) Laboratory Department
More informationOnboard Engine FDI in Autonomous Aircraft Using Stochastic Nonlinear Modelling of Flight Signal Dependencies
Onboard Engine FDI in Autonomous Aircraft Using Stochastic Nonlinear Modelling of Flight Signal Dependencies Dimitrios G. Dimogianopoulos, John D. Hios and Spilios D. Fassois Stochastic Mechanical Systems
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 33 Probabilistic methods in earthquake engineering-2 So, we have
More informationEL1820 Modeling of Dynamical Systems
EL1820 Modeling of Dynamical Systems Lecture 9 - Parameter estimation in linear models Model structures Parameter estimation via prediction error minimization Properties of the estimate: bias and variance
More information742. Time-varying systems identification using continuous wavelet analysis of free decay response signals
74. Time-varying systems identification using continuous wavelet analysis of free decay response signals X. Xu, Z. Y. Shi, S. L. Long State Key Laboratory of Mechanics and Control of Mechanical Structures
More informationTIME SERIES ANALYSIS. Forecasting and Control. Wiley. Fifth Edition GWILYM M. JENKINS GEORGE E. P. BOX GREGORY C. REINSEL GRETA M.
TIME SERIES ANALYSIS Forecasting and Control Fifth Edition GEORGE E. P. BOX GWILYM M. JENKINS GREGORY C. REINSEL GRETA M. LJUNG Wiley CONTENTS PREFACE TO THE FIFTH EDITION PREFACE TO THE FOURTH EDITION
More informationOperational modal analysis using forced excitation and input-output autoregressive coefficients
Operational modal analysis using forced excitation and input-output autoregressive coefficients *Kyeong-Taek Park 1) and Marco Torbol 2) 1), 2) School of Urban and Environment Engineering, UNIST, Ulsan,
More informationOn Input Design for System Identification
On Input Design for System Identification Input Design Using Markov Chains CHIARA BRIGHENTI Masters Degree Project Stockholm, Sweden March 2009 XR-EE-RT 2009:002 Abstract When system identification methods
More informationElements of Multivariate Time Series Analysis
Gregory C. Reinsel Elements of Multivariate Time Series Analysis Second Edition With 14 Figures Springer Contents Preface to the Second Edition Preface to the First Edition vii ix 1. Vector Time Series
More informationOutput Only Parametric Identification of a Scale Cable Stayed Bridge Structure: a comparison of vector AR and stochastic subspace methods
Output Only Parametric Identification of a Scale Cable Stayed Bridge Structure: a comparison of vector AR and stochastic subspace methods Fotis P. Kopsaftopoulos, Panagiotis G. Apostolellis and Spilios
More informationSINGLE DEGREE OF FREEDOM SYSTEM IDENTIFICATION USING LEAST SQUARES, SUBSPACE AND ERA-OKID IDENTIFICATION ALGORITHMS
3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 24 Paper No. 278 SINGLE DEGREE OF FREEDOM SYSTEM IDENTIFICATION USING LEAST SQUARES, SUBSPACE AND ERA-OKID IDENTIFICATION
More informationStochastic Dynamics of SDOF Systems (cont.).
Outline of Stochastic Dynamics of SDOF Systems (cont.). Weakly Stationary Response Processes. Equivalent White Noise Approximations. Gaussian Response Processes as Conditional Normal Distributions. Stochastic
More informationMulti Channel Output Only Identification of an Extendable Arm Structure Under Random Excitation: A comparison of parametric methods
Multi Channel Output Only Identification of an Extendable Arm Structure Under Random Excitation: A comparison of parametric methods Minas Spiridonakos and Spilios Fassois Stochastic Mechanical Systems
More information6.435, System Identification
System Identification 6.435 SET 3 Nonparametric Identification Munther A. Dahleh 1 Nonparametric Methods for System ID Time domain methods Impulse response Step response Correlation analysis / time Frequency
More informationTemporal-Difference Q-learning in Active Fault Diagnosis
Temporal-Difference Q-learning in Active Fault Diagnosis Jan Škach 1 Ivo Punčochář 1 Frank L. Lewis 2 1 Identification and Decision Making Research Group (IDM) European Centre of Excellence - NTIS University
More informationIdentification of Time-Variant Systems Using Wavelet Analysis of Force and Acceleration Response Signals
LOGO IOMAC'11 4th International Operational Modal Analysis Conference Identification of Time-Variant Systems Using Wavelet Analysis of Force and Acceleration Response Signals X. Xu 1,, W. J. Staszewski
More informationStructural Dynamics Lecture 7. Outline of Lecture 7. Multi-Degree-of-Freedom Systems (cont.) System Reduction. Vibration due to Movable Supports.
Outline of Multi-Degree-of-Freedom Systems (cont.) System Reduction. Truncated Modal Expansion with Quasi-Static Correction. Guyan Reduction. Vibration due to Movable Supports. Earthquake Excitations.
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard c Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 Outline 1 Lecture 9: Model Parametrization 2
More informationTime series methods for fault detection and identification in vibrating structures
Time series methods for fault detection and identification in vibrating structures By Spilios D. Fassois and John S. Sakellariou Stochastic Mechanical Systems (SMS) Group Department of Mechanical & Aeronautical
More informationNon-stationary functional series modeling and analysis of hardware reliability series: a comparative study using rail vehicle interfailure times
Reliability Engineering and System Safety 68 (2000) 169 183 www.elsevier.com/locate/ress Non-stationary functional series modeling and analysis of hardware reliability series: a comparative study using
More informationOUTPUT-ONLY STATISTICAL TIME SERIES METHODS FOR STRUCTURAL HEALTH MONITORING: A COMPARATIVE STUDY
7th European Workshop on Structural Health Monitoring July 8-11, 2014. La Cité, Nantes, France More Info at Open Access Database www.ndt.net/?id=17198 OUTPUT-ONLY STATISTICAL TIME SERIES METHODS FOR STRUCTURAL
More informationA time series is called strictly stationary if the joint distribution of every collection (Y t
5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a
More informationFurther Results on Model Structure Validation for Closed Loop System Identification
Advances in Wireless Communications and etworks 7; 3(5: 57-66 http://www.sciencepublishinggroup.com/j/awcn doi:.648/j.awcn.735. Further esults on Model Structure Validation for Closed Loop System Identification
More informationStatistics 910, #5 1. Regression Methods
Statistics 910, #5 1 Overview Regression Methods 1. Idea: effects of dependence 2. Examples of estimation (in R) 3. Review of regression 4. Comparisons and relative efficiencies Idea Decomposition Well-known
More informationLecture 2: ARMA(p,q) models (part 2)
Lecture 2: ARMA(p,q) models (part 2) Florian Pelgrin University of Lausanne, École des HEC Department of mathematics (IMEA-Nice) Sept. 2011 - Jan. 2012 Florian Pelgrin (HEC) Univariate time series Sept.
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Prediction Error Methods - Torsten Söderström
PREDICTIO ERROR METHODS Torsten Söderström Department of Systems and Control, Information Technology, Uppsala University, Uppsala, Sweden Keywords: prediction error method, optimal prediction, identifiability,
More informationRozwiązanie zagadnienia odwrotnego wyznaczania sił obciąŝających konstrukcje w czasie eksploatacji
Rozwiązanie zagadnienia odwrotnego wyznaczania sił obciąŝających konstrukcje w czasie eksploatacji Tadeusz Uhl Piotr Czop Krzysztof Mendrok Faculty of Mechanical Engineering and Robotics Department of
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More information6.435, System Identification
SET 6 System Identification 6.435 Parametrized model structures One-step predictor Identifiability Munther A. Dahleh 1 Models of LTI Systems A complete model u = input y = output e = noise (with PDF).
More informationSystem Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang
System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang 1-1 Course Description Emphases Delivering concepts and Practice Programming Identification Methods using Matlab Class
More informationStructural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake
Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Seung Keun Park and Hae Sung Lee ABSTRACT This paper presents a system identification (SI) scheme
More informationHierarchical sparse Bayesian learning for structural health monitoring. with incomplete modal data
Hierarchical sparse Bayesian learning for structural health monitoring with incomplete modal data Yong Huang and James L. Beck* Division of Engineering and Applied Science, California Institute of Technology,
More informationOn the Nature of Random System Matrices in Structural Dynamics
On the Nature of Random System Matrices in Structural Dynamics S. ADHIKARI AND R. S. LANGLEY Cambridge University Engineering Department Cambridge, U.K. Nature of Random System Matrices p.1/20 Outline
More informationEL1820 Modeling of Dynamical Systems
EL1820 Modeling of Dynamical Systems Lecture 10 - System identification as a model building tool Experiment design Examination and prefiltering of data Model structure selection Model validation Lecture
More informationSEISMIC RELIABILITY ANALYSIS OF BASE-ISOLATED BUILDINGS
International Symposium on Engineering under Uncertainty: Safety Assessment and Management January 4 to 6, 2012 Paper No.: CNP 070 SEISMIC RELIABILITY ANALYSIS OF BASE-ISOLATED BUILDINGS M.C. Jacob 1,
More informationcovariance function, 174 probability structure of; Yule-Walker equations, 174 Moving average process, fluctuations, 5-6, 175 probability structure of
Index* The Statistical Analysis of Time Series by T. W. Anderson Copyright 1971 John Wiley & Sons, Inc. Aliasing, 387-388 Autoregressive {continued) Amplitude, 4, 94 case of first-order, 174 Associated
More informationControl Systems Lab - SC4070 System Identification and Linearization
Control Systems Lab - SC4070 System Identification and Linearization Dr. Manuel Mazo Jr. Delft Center for Systems and Control (TU Delft) m.mazo@tudelft.nl Tel.:015-2788131 TU Delft, February 13, 2015 (slides
More informationLessons learned from the theory and practice of. change detection. Introduction. Content. Simulated data - One change (Signal and spectral densities)
Lessons learned from the theory and practice of change detection Simulated data - One change (Signal and spectral densities) - Michèle Basseville - 4 6 8 4 6 8 IRISA / CNRS, Rennes, France michele.basseville@irisa.fr
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Basics of System Identification Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont (UBC) EECE574 - Basics of
More information1.1 OBJECTIVE AND CONTENTS OF THE BOOK
1 Introduction 1.1 OBJECTIVE AND CONTENTS OF THE BOOK Hysteresis is a nonlinear phenomenon exhibited by systems stemming from various science and engineering areas: under a low-frequency periodic excitation,
More informationResearch Article Experimental Parametric Identification of a Flexible Beam Using Piezoelectric Sensors and Actuators
Shock and Vibration, Article ID 71814, 5 pages http://dx.doi.org/1.1155/214/71814 Research Article Experimental Parametric Identification of a Flexible Beam Using Piezoelectric Sensors and Actuators Sajad
More informationSENSITIVITY ANALYSIS OF ADAPTIVE MAGNITUDE SPECTRUM ALGORITHM IDENTIFIED MODAL FREQUENCIES OF REINFORCED CONCRETE FRAME STRUCTURES
SENSITIVITY ANALYSIS OF ADAPTIVE MAGNITUDE SPECTRUM ALGORITHM IDENTIFIED MODAL FREQUENCIES OF REINFORCED CONCRETE FRAME STRUCTURES K. C. G. Ong*, National University of Singapore, Singapore M. Maalej,
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi Lecture 13-29 May, 2013 Courtesy of Prof. S. Fassois & Dr. F. Kopsaftopoulos, SMSA Group, University of Patras Statistical methods for SHM courtesy of Prof. S. Fassois & Dr. F. Kopsaftopoulos,
More informationy k = ( ) x k + v k. w q wk i 0 0 wk
Four telling examples of Kalman Filters Example : Signal plus noise Measurement of a bandpass signal, center frequency.2 rad/sec buried in highpass noise. Dig out the quadrature part of the signal while
More informationNew Developments in Tail-Equivalent Linearization method for Nonlinear Stochastic Dynamics
New Developments in Tail-Equivalent Linearization method for Nonlinear Stochastic Dynamics Armen Der Kiureghian President, American University of Armenia Taisei Professor of Civil Engineering Emeritus
More informationSystem Identification procedures for nonlinear response of Buckling Restraint Braces J. Martínez 1, R. Boroschek 1, J. Bilbao 1 (1)University of Chile
System Identification procedures for nonlinear response of Buckling Restraint Braces J. Martínez, R. Boroschek, J. Bilbao ()University of Chile. Abstract Buckling Restrained Braces (BRB) are hysteretic
More information751. System identification of rubber-bearing isolators based on experimental tests
75. System identification of rubber-bearing isolators based on experimental tests Qiang Yin,, Li Zhou, Tengfei Mu, Jann N. Yang 3 School of Mechanical Engineering, Nanjing University of Science and Technology
More informationParametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d)
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d) Electrical & Computer Engineering North Carolina State University Acknowledgment: ECE792-41 slides
More informationPredictive Control of Gyroscopic-Force Actuators for Mechanical Vibration Damping
ARC Centre of Excellence for Complex Dynamic Systems and Control, pp 1 15 Predictive Control of Gyroscopic-Force Actuators for Mechanical Vibration Damping Tristan Perez 1, 2 Joris B Termaat 3 1 School
More informationEcon 424 Time Series Concepts
Econ 424 Time Series Concepts Eric Zivot January 20 2015 Time Series Processes Stochastic (Random) Process { 1 2 +1 } = { } = sequence of random variables indexed by time Observed time series of length
More informationModeling and Identification of Dynamic Systems (vimmd312, 2018)
Modeling and Identification of Dynamic Systems (vimmd312, 2018) Textbook background of the curriculum taken. In parenthesis: material not reviewed due to time shortage, but which is suggested to be read
More informationTAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω
ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.
More informationParticle Relaxation Method of Monte Carlo Filter for Structure System Identification
Civil Structural Health Monitoring Workshop (CSHM-4) - Lecture 11 Particle Relaxation Method of Monte Carlo Filter for Structure System Identification Tadanobu SATO International Institute of Urban Engineering,
More informationPh.D student in Structural Engineering, Department of Civil Engineering, Ferdowsi University of Mashhad, Azadi Square, , Mashhad, Iran
Alireza Entezami a, Hashem Shariatmadar b* a Ph.D student in Structural Engineering, Department of Civil Engineering, Ferdowsi University of Mashhad, Azadi Square, 9177948974, Mashhad, Iran b Associate
More informationInput-Output Peak Picking Modal Identification & Output only Modal Identification and Damage Detection of Structures using
Input-Output Peak Picking Modal Identification & Output only Modal Identification and Damage Detection of Structures using Time Frequency and Wavelet Techniquesc Satish Nagarajaiah Professor of Civil and
More informationAnalysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems
Analysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems Jeremy S. Conner and Dale E. Seborg Department of Chemical Engineering University of California, Santa Barbara, CA
More information1 Estimation of Persistent Dynamic Panel Data. Motivation
1 Estimation of Persistent Dynamic Panel Data. Motivation Consider the following Dynamic Panel Data (DPD) model y it = y it 1 ρ + x it β + µ i + v it (1.1) with i = {1, 2,..., N} denoting the individual
More informationNONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM
NONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM D.G. Dimogianopoulos, J.D. Hios and S.D. Fassois DEPARTMENT OF MECHANICAL & AERONAUTICAL ENGINEERING GR-26500 PATRAS,
More informationARIMA Modelling and Forecasting
ARIMA Modelling and Forecasting Economic time series often appear nonstationary, because of trends, seasonal patterns, cycles, etc. However, the differences may appear stationary. Δx t x t x t 1 (first
More informationIdentification of ARX, OE, FIR models with the least squares method
Identification of ARX, OE, FIR models with the least squares method CHEM-E7145 Advanced Process Control Methods Lecture 2 Contents Identification of ARX model with the least squares minimizing the equation
More informationCh 5. Models for Nonstationary Time Series. Time Series Analysis
We have studied some deterministic and some stationary trend models. However, many time series data cannot be modeled in either way. Ex. The data set oil.price displays an increasing variation from the
More informationSYSTEM IDENTIFICATION & DAMAGE ASSESSMENT OF STRUCTURES USING OPTICAL TRACKER ARRAY DATA
SYSTEM IDENTIFICATION & DAMAGE ASSESSMENT OF STRUCTURES USING OPTICAL TRACKER ARRAY DATA Chin-Hsiung Loh 1,* and Chuan-Kai Chan 1 1 Department of Civil Engineering, National Taiwan University Taipei 10617,
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 32 Probabilistic Methods in Earthquake Engineering-1 (Refer Slide
More informationTime Series Models and Inference. James L. Powell Department of Economics University of California, Berkeley
Time Series Models and Inference James L. Powell Department of Economics University of California, Berkeley Overview In contrast to the classical linear regression model, in which the components of the
More information7. Forecasting with ARIMA models
7. Forecasting with ARIMA models 309 Outline: Introduction The prediction equation of an ARIMA model Interpreting the predictions Variance of the predictions Forecast updating Measuring predictability
More informationRoss Bettinger, Analytical Consultant, Seattle, WA
ABSTRACT DYNAMIC REGRESSION IN ARIMA MODELING Ross Bettinger, Analytical Consultant, Seattle, WA Box-Jenkins time series models that contain exogenous predictor variables are called dynamic regression
More informationExercises - Time series analysis
Descriptive analysis of a time series (1) Estimate the trend of the series of gasoline consumption in Spain using a straight line in the period from 1945 to 1995 and generate forecasts for 24 months. Compare
More informationIdentification of Nonlinear Mechanical Systems: State of the Art and Recent Trends
- Acceleration (m/s 2 ) Identification of Nonlinear Mechanical Systems: State of the Art and Recent Trends 5 Force level: 67 N rms Gaëtan Kerschen Space Structures and Systems Laboratory Aerospace and
More informationMultivariate Time Series: VAR(p) Processes and Models
Multivariate Time Series: VAR(p) Processes and Models A VAR(p) model, for p > 0 is X t = φ 0 + Φ 1 X t 1 + + Φ p X t p + A t, where X t, φ 0, and X t i are k-vectors, Φ 1,..., Φ p are k k matrices, with
More informationStatistics 910, #15 1. Kalman Filter
Statistics 910, #15 1 Overview 1. Summary of Kalman filter 2. Derivations 3. ARMA likelihoods 4. Recursions for the variance Kalman Filter Summary of Kalman filter Simplifications To make the derivations
More informationLecture 8: Signal Detection and Noise Assumption
ECE 830 Fall 0 Statistical Signal Processing instructor: R. Nowak Lecture 8: Signal Detection and Noise Assumption Signal Detection : X = W H : X = S + W where W N(0, σ I n n and S = [s, s,..., s n ] T
More informationMultivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8]
1 Multivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8] Insights: Price movements in one market can spread easily and instantly to another market [economic globalization and internet
More informationGiacomo Boffi. Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano
http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 21, 2017 Outline of Structural Members Elastic-plastic Idealization
More informationCodal Provisions IS 1893 (Part 1) 2002
Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893
More informationScalar and Vector Time Series Methods for Vibration Based Damage Diagnosis in a Scale Aircraft Skeleton Structure
Scalar and Vector Time Series Methods for Vibration Based Damage Diagnosis in a Scale Aircraft Skeleton Structure Fotis P. Kopsaftopoulos and Spilios D. Fassois Stochastic Mechanical Systems & Automation
More informationMatlab software tools for model identification and data analysis 11/12/2015 Prof. Marcello Farina
Matlab software tools for model identification and data analysis 11/12/2015 Prof. Marcello Farina Model Identification and Data Analysis (Academic year 2015-2016) Prof. Sergio Bittanti Outline Data generation
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationEffects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device
Effects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device A. Krishnamoorthy Professor, Department of Civil Engineering Manipal
More informationBrief Review on Estimation Theory
Brief Review on Estimation Theory K. Abed-Meraim ENST PARIS, Signal and Image Processing Dept. abed@tsi.enst.fr This presentation is essentially based on the course BASTA by E. Moulines Brief review on
More informationRESPONSE SPECTRUM METHOD FOR ESTIMATION OF PEAK FLOOR ACCELERATION DEMAND
RESPONSE SPECTRUM METHOD FOR ESTIMATION OF PEAK FLOOR ACCELERATION DEMAND Shahram Taghavi 1 and Eduardo Miranda 2 1 Senior catastrophe risk modeler, Risk Management Solutions, CA, USA 2 Associate Professor,
More informationReliability Theory of Dynamically Loaded Structures (cont.)
Outline of Reliability Theory of Dynamically Loaded Structures (cont.) Probability Density Function of Local Maxima in a Stationary Gaussian Process. Distribution of Extreme Values. Monte Carlo Simulation
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More informationTHE subject of the analysis is system composed by
MECHANICAL VIBRATION ASSIGNEMENT 1 On 3 DOF system identification Diego Zenari, 182160, M.Sc Mechatronics engineering Abstract The present investigation carries out several analyses on a 3-DOF system.
More informationTime Series: Theory and Methods
Peter J. Brockwell Richard A. Davis Time Series: Theory and Methods Second Edition With 124 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vn ix CHAPTER 1 Stationary
More informationANNEX A: ANALYSIS METHODOLOGIES
ANNEX A: ANALYSIS METHODOLOGIES A.1 Introduction Before discussing supplemental damping devices, this annex provides a brief review of the seismic analysis methods used in the optimization algorithms considered
More informationOn Identification of Cascade Systems 1
On Identification of Cascade Systems 1 Bo Wahlberg Håkan Hjalmarsson Jonas Mårtensson Automatic Control and ACCESS, School of Electrical Engineering, KTH, SE-100 44 Stockholm, Sweden. (bo.wahlberg@ee.kth.se
More informationAnalyzing the Effect of Moving Resonance on Seismic Response of Structures Using Wavelet Transforms
Analyzing the Effect of Moving Resonance on Seismic Response of Structures Using Wavelet Transforms M.R. Eatherton Virginia Tech P. Naga WSP Cantor Seinuk, New York, NY SUMMARY: When the dominant natural
More informationf-domain expression for the limit model Combine: 5.12 Approximate Modelling What can be said about H(q, θ) G(q, θ ) H(q, θ ) with
5.2 Approximate Modelling What can be said about if S / M, and even G / G? G(q, ) H(q, ) f-domain expression for the limit model Combine: with ε(t, ) =H(q, ) [y(t) G(q, )u(t)] y(t) =G (q)u(t) v(t) We know
More informationSTAT 443 Final Exam Review. 1 Basic Definitions. 2 Statistical Tests. L A TEXer: W. Kong
STAT 443 Final Exam Review L A TEXer: W Kong 1 Basic Definitions Definition 11 The time series {X t } with E[X 2 t ] < is said to be weakly stationary if: 1 µ X (t) = E[X t ] is independent of t 2 γ X
More informationOn the use of wavelet coefficient energy for structural damage diagnosis
Safety, Reliability and Risk of Structures, Infrastructures and Engineering Systems Furuta, Frangopol & Shinozuka (eds) 010 Taylor & Francis Group, London, ISBN 978-0-415-47557-0 On the use of wavelet
More informationNon-Stationary Random Vibration Parametric Modeling and its Application to Structural Health Monitoring
Non-Stationary Random Vibration Parametric Modeling and its Application to Structural Health Monitoring Luis David Avendaño-Valencia and Spilios D. Fassois Stochastic Mechanical Systems and Automation
More information