An example of correlation matrix based mode shape expansion in OMA
|
|
- Edwina Campbell
- 6 years ago
- Views:
Transcription
1 An example of correlation matrix based mode shape expansion in OMA Rune Brincker 1 Edilson Alexandre Camargo 2 Anders Skafte 1 1 : Department of Engineering, Aarhus University, Aarhus, Denmark 2 : Institute of Aeronautics and Space - IAE, São José dos Campos, Brazil Abstract In cases of reducing the number of modes in an operating response for instance when using band pass filtering often time domain identification techniques suffer from a tendency to over fitting. As a consequence, it might be useful to decrease the number of measured DOF s to moderate the number of modes in the modal model. In these cases the dimension of the identified mode shapes is diminished accordingly and therefore it is preferable to have a way to expand back to the full set of DOF s so that the estimated mode shapes can be animated in detail. In the present paper it is considered to use the correlation matrix for the filtered response including all measured DOF s as a basis for the mode shape expansion. This matrix contains the normal modes that can easily be extracted from the column space of the correlation matrix. In this investigation the main focus has not been on the reduction problem, so engineering judgment has been used to secure a reasonable choice of reduced channels. The technique is illustrated on an OMA case where the modes of the tail part of a Panther helicopter is estimated during different flight conditions. Keywords: OMA, DOF reduction, Expansion, Helicopter, Time domain 1. Introduction Nearly all classical modal analysis techniques in the time domain can be used for OMA using the correlation functions as input for the techniques instead of impulse response functions. This is true for instance for techniques like the Ibrahim time domain (ITD) technique, [2-5], the poly reference (PR) technique [6-7] and the Eigen realization algorithm (ERA) technique [8-10] that has been used in the present analysis. These techniques typically have a lower limit to (or a fixed size like the classical ITD) to the number of modes included in the model. This lower limit is normally equal to the number of measured degrees of freedom (DOF s). For instance in the PR technique, normally the lowest model order to be considered corresponds to 2 auto regressive matrices, defining the lowest number of modes in the model to be equal to the number of measured DOF s. This is also the case for the classical ITD technique. It is well known that it is normally an advantage to use an oversized model including a reasonable number of noise modes in order to obtain good modal estimates. However, it is also well known, that if too many noise modes are being used, problems with over fitting might arise, and the modal results can be affected in a negative way. The over fitting problem often becomes severe in cases with a high amount of noise in the signals. In such case it might be useful to consider a reduced number of measured DOF s in the identification, so that the number of noise modes can be reduced accordingly. If the modes are identified in a reduced number of DOF s, then the estimated mode shapes are only known in these DOF s, and therefore, after the identification in the reduced set, it is useful to expand the mode shapes to the full set of DOF s. In the following we will consider the case where the random response to be used for OMA is band pass filtered in order to reduce the number of modes present in the random response. This means that the number of modes is often quite small such as 2-4 modes. Further we will consider the possibility to perform the mode shape expansion based on the correlation matrix of the band pass filtered response.
2 First we will present the basic theory for why it makes sense to use the correlation matrix for expansion, the main reasons being that in case of band pass filtered signals containing only a limited number of modes and in case of OMA tests with a reasonable channel count, the correlation matrix provides a very good way of checking the number of modes in the band pass filtered signal and for expansion. The resulting modes will be normal modes (real mode shapes) so that direct comparison can be made with a finite element model. It is clear that going from a high number of measured DOF s to a small number of DOF s for identification might also be dangerous, because in such case there is a risk of excluding important modal information if the number of identification DOF s are not chosen reasonably. In this investigation however the main focus is not been on the reduction problem, so engineering judgment has been used to secure a reasonable choice of reduced channels. In the following we will illustrate how the classical techniques mentioned above combined with the proposed reduction/expansion technique perform on a difficult OMA case with a relative high channel count and high amount of noise. The considered case is OMA of the tail part of a Panther helicopter estimated during different flight conditions. 2. Correlation matrix theory Let us consider a zero mean random response y (t) with N d number of measured DOF s, i.e. the response vector y (t) is a column vector with N d number of components. In this case the response correlation matrix is Nd Nd and is defined by T (1) C E y ( t) y ( t) y For a measured time series, the matrix will normally be calculated in terms of time averaging (2) 1 T C y y( t) y T ( t) dt T0 0 0 where T 0 is the length of the considered random response. However, if the response is given as a spectral density function, we can use that the correlation function matrix is defined as T (3) R ( ) E y ( t) y ( t ) So that the correlation matrix can be found as (4) Cy R y (0) y Since the correlation function matrix is the inverse Fourier transform of the of spectral density matrix we can obtain the correlation matrix from (5) C y G y ( ) d Thus, the correlation matrix can be found as the integral of the spectral density matrix. It should be noted that in case we are dealing with sampled data, the integration should be performed only over the Nyquist band. Further it should be noted that since the imaginary part of a spectral density function is an odd function of frequency, the integral will always be real as it is given by the definition in Eq. (1). Assuming that the measured response is a linear combination of N normal modes we have
3 (6) y( t) a1q1 ( t) a2q2( t) an qn ( t) Aq( t) where A is the mode shape matrix A a1, a2,, a N q t ) q ( t), q ( t),, q T N ( t). Using this is Eq. (1) we have ( 1 2 T T T (7) Cy AE q( t) q ( t) A and q (t) is the modal coordinate vector ACqA where C q is the correlation matrix of the modal coordinates. 3. Expansion from a smaller set of DOF s We shall assume that the external loading on the considered structure is not of a kind where C q has a rank that is smaller than the number of modes, in fact we shall consider a situation where the modal coordinates are approximately stochastically independent. Further we shall assume that the number of modes N in the response is significantly smaller than the number N d of measured DOF s. In this case, when taking a singular value decomposition of the response correlation matrix (8) T Cy USU any mode shape a participating in the response given by Eq. (6) is in the subspace U and the rank of C y is equal to the number of modes N in the response. We then have the linear relation between the modal matrix A and the subspace U (9) A UT where the transformation matrix T is DOF s, N N r N N. If we have reduced the response to consider only N r number of measured d, then we have the similar equation in the reduced set of DOF s (10) Ar UrT If we have obtained the experimental modes in the reduced set of DOF s, we can find the transformation matrix from Eq. (10). If the reduced number of DOF s Nr is larger than the number of modes, then the equation can be solved approximately (11) Tˆ U r Ar where U r is the pseudo inverse of the matrix U r. In this case we solve the equation in an over determined way and the estimate Tˆ is then less sensitive to the noise on the mode shape estimates. In cases where the reduced set of DOF s is equal to the number modes (usually a model without noise modes), Eq. (10) is a normal equation that can be solved taking the inverse of the matrix U r that is now square (and full rank) (12) Tˆ 1 U r Ar In this case it should be considered if the noise on the mode shape estimate might significantly influence the transformation matrix estimate Tˆ.
4 When we have an estimate of the transformation matrix we can expand the mode shapes to full size using Eq. (9) to get an estimate of the mode shapes in all the measured DOF s (13) Aˆ UTˆ 4. Panther helicopter test case The test case is an investigation of the tail part of the Panther helicopter shown in Figure 1.a. The helicopter was tested in different flight conditions where some results for two of the flight conditions are reported here. A more comprehensive reporting of the results of all investigated flight conditions can be found in Camargo and Brincker [11]. The measurement of the flight condition analyzed in this article, were performed with the aircraft flying in a fixed altitude of 5000 feet with cruising speed under constant weather conditions and payload (Figure 1a). The altitude of the flying, at cruising speed, is taking place according to the simple principle that the pitch of the rotor determines the forward speed, and power determines the altitude of the aircraft. As a consequence, in the cruising speed condition considered here, the flight condition is characterized by constant rotor pitch angle and constant power consumption of the engine. The experimental setup was configured with 40 measurement channels, divided into three measurement units with 17 channels each, using 16 unidirectional accelerometers 4515 type (8 in X axis and 8 in Y axis) and 10 tri-axial accelerometer 4520 type. Of the 10 tri-axial accelerometers four where using all three channels, and 2 where using only two of the three channels in the OMA. The placement of the sensors is shown in Figure 1b. The data acquisition was performed using the software called LabShop from Brüel & Kjaer (B&K) and the analysis were done, in time domain, using the OMA MatLab toolbox related to the OMA book by Brincker and Ventura [1]. (a) (b) Figure 1: (a) The Panther helicopter at cruiser speed, (b) The sensor locations, left: all 40 measured DOF s, right: The four DOF s used in the OMA (reduced set). The measurements of the full set of 40 DOFs were performed for 60 seconds at an acquisition rate of 2,000 samples per second. Analysis were performed using just 4 DOFs, located at the edge of the horizontal stabilizers (HS), see Figure 1b right, decimating the data to a Nyquist frequency of 80 Hz with sampling interval of seconds. This resulted in a spectral density estimate with a resolution of 1024 lines giving a frequency step of 0.04 Hz as shown in Figure 2. A main goal of the test was to identify the first four modes of the tail part of the helicopter, that is the XX mode, that is the first tail torsion mode, the first HS bending mode, the first tail bending mode, and finally the first HS torsional mode. The natural frequencies of these modes were known to be around 12,8 Hz, 14,8 Hz, 17,7 Hz and 19,5 Hz respectively.
5 In order to isolate these four modes, the data were band-pass filtered using the fftfilt.m function in the OMA toolbox for each mode with the center frequency of 12,8 Hz, 14,8 Hz, 17,7 Hz and 19,5 Hz respectively for each mode, a flat band around the center frequency with a band width of 1 Hz and roll-off bands on each side with a width of 1 Hz. Figure 2 shows the highlighted peaks of the four modes with the small cross marking the central frequency. The first harmonic of the main rotor of the helicopter at 23,6 Hz is also indicated in Figure 2. Figure 2. FDD plot showing the singular values of the spectral density matrix of the cruiser speed measurements. The highlighted peaks indicate the first four modes of the tail part and the first harmonic of the rotor is also indicated in the plot. The numbers of block rows on the AR/PR and ITD were chosen to na = 2. For the ERA, the numbers of block rows of the Hankel matrices used was equal to s = 2, and the realization was not reduced so that the returned number of modes were the same for the three techniques. In all three techniques, the correlation matrix was used as multiple input, and 100 discrete time lags were used in the identification. Using all 40 measurement channels the techniques return in principle 40 modes of which some might be non-physical due to non-positive damping, so typically modes were returned. Since we have band pass filtered the signal to contain only one single mode, among the returned modes we need to identify only the single mode that we are looking for. In the following we will illustrate the identification principle on the first of the four modes indicated in Figure 2, this mode has a natural frequency around 12.8 Hz. The modal participation for each mode was estimated as described in Brincker and Ventura [1], and used to identify the physical mode following the principle that noise modes should have a small modal participation and a physical mode should have a strong modal participation. Some of the returned modes are shown in Table 1 for the case where we are using all 40 DOF s. From the results in Table 1 we see that the mode with the highest participation using the PR technique is mode 13 with a modal participation of 22.9 %. Table 1: OMA of the first mode (12.8 Hz) using 40 DOF s PR Mode 1 Mode 2 Mode 3 Mode 4 Mode 11 Mode 13 Mode 15 Mode 27 Frequency, [Hz] Damping, [%] Contribution, [%] ITD Mode 1 Mode 2 Mode 3 Mode 4 Mode 7 Mode 9 Mode 11 Mode 19 Frequency, [Hz] Damping, [%] Contribution, [%] ERA Mode 1 Mode 2 Mode 3 Mode 4 Mode 11 Mode 13 Mode 15 Mode 27 Frequency, [Hz] Damping, [%] Contribution, [%]
6 It seems like a reasonable result since the corresponding natural frequency of Hz is also close to the expected value. However, we see that the modal participation is also high for adjacent modes with frequencies quite close to the target frequency. A similar picture is seen to be present for identification using ITD and ERA identification. As a result in this case it is difficult to use the modal participation clearly to distinguish between noise modes and physical modes Table 2: OMA of the first mode (12.8 Hz) using 4 DOF s PR Mode 1 Mode 2 Mode 3 Mode 4 Frequency, [Hz] Damping, [%] Contribution, [%] ITD Mode 1 Mode 2 Mode 3 Mode 4 Frequency, [Hz] Damping, [%] Contribution, [%] ERA Mode 1 Mode 2 Mode 3 Mode 4 Frequency, [Hz] Damping, [%] Contribution, [%] Reducing the number of DOFs to the 4 sensor signals shown in Figure 1.b, the identification techniques return only 4 modes. Table 2 shows the similar results for the first mode. We see that in this case mode 2 has a modal participation that is significantly larger than all other modal participation factors, and thus, using the reduced set of DOF s we can much easier make a decision of which mode to choose as the physical one. The distribution of the modal participating factors of the modes is shown in Figure 3 and as it appears, all three techniques have a high modal participation of 2-3 modes making it difficult to distinguish noise modes from the physical mode. The plot to the left shows the similar distribution of the modal participation factors in the case of using the reduced set of DOF s and as it appears, in this case only one mode has a modal participation factor that is significantly larger than the other ones. Figure 3. The plots show the modal participation factors of the three different techniques PR (blue), ITD (red) and ERA (gray). 5. OMA results using expansion The top plots of Figure 4 shows the results when only the four DOF s are being used, and as it appears from the plots, using only 4 DOFs the mode shapes are not well defined in a graphical representation. The mode shapes were then expanded using Equation (11/13) and the results shown the bottom plots of Figure 4 clearly shows the advantage of using the mode shape in full resolution.
7 4 DOFs 40 DOFs First tail torsion mode (12.8 Hz) First HS bending mode (15Hz) First tail bending mode (17.8Hz) First HS torsional mode (19.5Hz) Figure 4: Mode shapes of 4 DOFs and after expand the modes for 40 DOFs. Top plot: Mode shape plots using only 4 DOFs. Middle plot: Full mode shapes for flight condition cruising speed, bottom plot: Full mode shapes for the flight condition maximum speed. 6. Conclusions It is shown that the correlation matrix can be used of expansion of mode shapes and it has been illustrated how the expansion can be applied in cases with high channels counts and where over fitting might occur when performing identification using classical time domain techniques. The use of some classical time domain techniques has been illustrated using OMA on a Panther helicopter under different flight conditions. Acknowledgements The financial support received from CAPES, through the Science without Borders Program, for this research, is gratefully acknowledged. References [1] Brincker, R. and Ventura, C.E.: Introduction to operational modal anlaysis. Wiley [2] Ibrahim, S.R. and Milkulcik, E.C.: A method for direct identification of vibration parameters from the free response. The Shock and Vibration Bulletin, 47, p , [3] Ibrahim, S.R.: Random decrement technique for modal identification of structures. J. of Spacecraft and Rockets, 14, p , [4] Ibrahim, S.R.: Modal confidence factor in vibration testing. J. of Spacecraft and Rockets, 15, p , [5] Fukuzono, K.: Investigation of multiple-reference Ibrahim time domain modal parameter estimation technique. M.Sc. Thesis, Department of Mechanical and Industrial Engineering, University of Cincinnati, [6] Vold, H., Kundrat, J., Rocklin, G.T. and Russell, R.:A multi-input modal estimation algorithm for mini-computers. SAE Paper NUmber , 1082.
8 [7] Vold, H. and Rocklin, G.T.: The numerical implementation of a multi-input modal estimation method for mini-computers. In Proc. of the International Modal Analysis Conference (IMAC), p , [8] Juang, J.N. and Pappa, R.S.: An eigen system realization algorithm for modal parameter identification and modal reduction. J. Guidance, V. 8, No. 5, p , [9] Pappa, R.S., Elliott, K.B. and Schenk, A.: Consistent-mode indicator for the eigensystem realization algorithm. J. Guidance, Control and Dynamics, V. 16, No. 5, p , [10] Pappa, R.S.: Eigensystem realization algorithm, user s guide for VAX/VMS computers. NASA technical memorandum , [11] Camargo, E.A. and Brincker, R.: Operational modal analysis of a Panther helicopter during different flight conditions. To be submitted to XX.
COMPARISON OF MODE SHAPE VECTORS IN OPERATIONAL MODAL ANALYSIS DEALING WITH CLOSELY SPACED MODES.
IOMAC'5 6 th International Operational Modal Analysis Conference 5 May-4 Gijón - Spain COMPARISON OF MODE SHAPE VECTORS IN OPERATIONAL MODAL ANALYSIS DEALING WITH CLOSELY SPACED MODES. Olsen P., and Brincker
More informationModal identification of output-only systems using frequency domain decomposition
INSTITUTE OF PHYSICS PUBLISHING SMART MATERIALS AND STRUCTURES Smart Mater. Struct. 10 (2001) 441 445 www.iop.org/journals/sm PII: S0964-1726(01)22812-2 Modal identification of output-only systems using
More informationIdentification Techniques for Operational Modal Analysis An Overview and Practical Experiences
Identification Techniques for Operational Modal Analysis An Overview and Practical Experiences Henrik Herlufsen, Svend Gade, Nis Møller Brüel & Kjær Sound and Vibration Measurements A/S, Skodsborgvej 307,
More informationImprovement of Frequency Domain Output-Only Modal Identification from the Application of the Random Decrement Technique
Improvement of Frequency Domain Output-Only Modal Identification from the Application of the Random Decrement Technique Jorge Rodrigues LNEC - National Laboratory for Civil Engineering, Structures Department
More informationModal Based Fatigue Monitoring of Steel Structures
Modal Based Fatigue Monitoring of Steel Structures Jesper Graugaard-Jensen Structural Vibration Solutions A/S, Denmark Rune Brincker Department of Building Technology and Structural Engineering Aalborg
More informationIN-FLIGHT MODAL IDENTIFICATION OF AN AIRCRAFT USING OPERATIONAL MODAL ANALYSIS ABSTRACT
9 th ANKARA INTERNATIONAL AEROSPACE CONFERENCE AIAC-2017-056 20-22 September 2017 - METU, Ankara TURKEY IN-FLIGHT MODAL IDENTIFICATION OF AN AIRCRAFT USING OPERATIONAL MODAL ANALYSIS Çağrı KOÇAN 1 and
More informationAutomated Modal Parameter Estimation For Operational Modal Analysis of Large Systems
Automated Modal Parameter Estimation For Operational Modal Analysis of Large Systems Palle Andersen Structural Vibration Solutions A/S Niels Jernes Vej 10, DK-9220 Aalborg East, Denmark, pa@svibs.com Rune
More informationRANDOM DECREMENT AND REGRESSION ANALYSIS OF TRAFFIC RESPONSES OF BRIDGES
RANDOM DECREMENT AND REGRESSION ANALYSIS OF TRAFFIC RESPONSES OF BRIDGES J.C. Asmussen, S.R. Ibrahid & R. Brincker Department of Building Technology and Structural Engineering Aalborg University, Sohngaardsholmsvej
More informationEliminating the Influence of Harmonic Components in Operational Modal Analysis
Eliminating the Influence of Harmonic Components in Operational Modal Analysis Niels-Jørgen Jacobsen Brüel & Kjær Sound & Vibration Measurement A/S Skodsborgvej 307, DK-2850 Nærum, Denmark Palle Andersen
More informationOperational Modal Analysis of Rotating Machinery
Operational Modal Analysis of Rotating Machinery S. Gres 2, P. Andersen 1, and L. Damkilde 2 1 Structural Vibration Solutions A/S, NOVI Science Park, Niels Jernes Vej 10, Aalborg, DK 9220, 2 Department
More informationOutput-Only Modal Analysis by Frequency Domain Decomposition Brincker, Rune; Zhang, Lingmi; Andersen, Palle
Aalborg Universitet Output-Only Modal Analysis by Frequency Domain Decomposition Brincker, Rune; Zhang, Lingmi; Andersen, Palle Published in: Proceedings of ISMA25 Publication date: 2000 Document Version
More informationAn Indicator for Separation of Structural and Harmonic Modes in Output-Only Modal Testing Brincker, Rune; Andersen, P.; Møller, N.
Aalborg Universitet An Indicator for Separation of Structural and Harmonic Modes in Output-Only Modal Testing Brincker, Rune; Andersen, P.; Møller, N. Published in: Proceedings of the European COST F3
More informationApplication of Classical and Output-Only Modal Analysis to a Laser Cutting Machine
Proc. of ISMA2002; Leuven, Belgium; 2002 Application of Classical and Output-Only Modal Analysis to a Laser Cutting Machine Carsten Schedlinski 1), Marcel Lüscher 2) 1) ICS Dr.-Ing. Carsten Schedlinski
More informationAalborg Universitet. Published in: Proceedings of ISMA2006. Publication date: Document Version Publisher's PDF, also known as Version of record
Aalborg Universitet Using Enhanced Frequency Domain Decomposition as a Robust Technique to Harmonic Excitation in Operational Modal Analysis Jacobsen, Niels-Jørgen; Andersen, Palle; Brincker, Rune Published
More informationMODAL IDENTIFICATION AND DAMAGE DETECTION ON A CONCRETE HIGHWAY BRIDGE BY FREQUENCY DOMAIN DECOMPOSITION
T1-1-a-4 SEWC2002, Yokohama, Japan MODAL IDENTIFICATION AND DAMAGE DETECTION ON A CONCRETE HIGHWAY BRIDGE BY FREQUENCY DOMAIN DECOMPOSITION Rune BRINCKER 1, Palle ANDERSEN 2, Lingmi ZHANG 3 1 Dept. of
More informationSmooth Orthogonal Decomposition for Modal Analysis of Randomly Excited Systems
Smooth Orthogonal Decomposition for Modal Analysis of Randomly Excited Systems U. Farooq and B. F. Feeny Department of Mechanical Engineering, Michigan State University, East Lansing, Michigan 48824 Abstract
More informationEstimation of Unsteady Loading for Sting Mounted Wind Tunnel Models
52nd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference 19th 4-7 April 2011, Denver, Colorado AIAA 2011-1941 Estimation of Unsteady Loading for Sting Mounted Wind Tunnel
More informationIndex 319. G Gaussian noise, Ground-loop current, 61
Index A Absolute accuracy, 72 ADCs. See Analog-to-digital converters (ADCs) AFDD-T. See Automated frequency domain decomposition-tracking (AFDD-T) Aliasing, 81, 90 Analog-to-digital converters (ADCs) conversion
More informationMASS, STIFFNESS AND DAMPING IDENTIFICATION OF A TWO-STORY BUILDING MODEL
COMPDYN 2 3 rd ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 25-28 May 2 MASS,
More informationAalborg Universitet. Modal Indicators for Operational Modal Identification Zhang, L.; Brincker, Rune; Andersen, P.
Aalborg Universitet Modal Indicators for Operational Modal Identification Zhang, L.; Brincker, Rune; Andersen, P. Published in: Proceedings of IMAC 19 Publication date: 2001 Document Version Publisher's
More informationIOMAC' May Guimarães - Portugal IMPACT-SYNCHRONOUS MODAL ANALYSIS (ISMA) AN ATTEMPT TO FIND AN ALTERNATIVE
IOMAC'13 5 th International Operational Modal Analysis Conference 2013 May 13-15 Guimarães - Portugal IMPACT-SYNCHRONOUS MODAL ANALYSIS (ISMA) AN ATTEMPT TO FIND AN ALTERNATIVE Abdul Ghaffar Abdul Rahman
More informationAssessment of the Frequency Domain Decomposition Method: Comparison of Operational and Classical Modal Analysis Results
Assessment of the Frequency Domain Decomposition Method: Comparison of Operational and Classical Modal Analysis Results Ales KUYUMCUOGLU Arceli A. S., Research & Development Center, Istanbul, Turey Prof.
More informationCONTRIBUTION TO THE IDENTIFICATION OF THE DYNAMIC BEHAVIOUR OF FLOATING HARBOUR SYSTEMS USING FREQUENCY DOMAIN DECOMPOSITION
CONTRIBUTION TO THE IDENTIFICATION OF THE DYNAMIC BEHAVIOUR OF FLOATING HARBOUR SYSTEMS USING FREQUENCY DOMAIN DECOMPOSITION S. Uhlenbrock, University of Rostock, Germany G. Schlottmann, University of
More informationVARIANCE COMPUTATION OF MODAL PARAMETER ES- TIMATES FROM UPC SUBSPACE IDENTIFICATION
VARIANCE COMPUTATION OF MODAL PARAMETER ES- TIMATES FROM UPC SUBSPACE IDENTIFICATION Michael Döhler 1, Palle Andersen 2, Laurent Mevel 1 1 Inria/IFSTTAR, I4S, Rennes, France, {michaeldoehler, laurentmevel}@inriafr
More informationIOMAC'15 6 th International Operational Modal Analysis Conference
IOMAC'15 6 th International Operational Modal Analysis Conference 2015 May12-14 Gijón - Spain PARAMETER ESTIMATION ALGORITMS IN OPERATIONAL MODAL ANALYSIS: A REVIEW Shashan Chauhan 1 1 Bruel & Kjær Sound
More informationIOMAC'15 6 th International Operational Modal Analysis Conference
IOMAC'15 6 th International Operational Modal Analysis Conference 2015 May12-14 Gijón - Spain COMPARISON OF DIFFERENT TECHNIQUES TO SCALE MODE SHAPES IN OPERATIONAL MODAL ANALYSIS. Luis Borja Peral Mtnez
More informationInvestigation of Operational Modal Analysis Damping Estimates MASTER OF SCIENCE
Investigation of Operational Modal Analysis Damping Estimates A thesis submitted to the Division of Research and Advanced Studies of the University of Cincinnati in partial fulfillment of the requirements
More informationDamping Estimation Using Free Decays and Ambient Vibration Tests Magalhães, Filipe; Brincker, Rune; Cunha, Álvaro
Aalborg Universitet Damping Estimation Using Free Decays and Ambient Vibration Tests Magalhães, Filipe; Brincker, Rune; Cunha, Álvaro Published in: Proceedings of the 2nd International Operational Modal
More information2671. Detection and removal of harmonic components in operational modal analysis
2671. Detection and removal of harmonic components in operational modal analysis Zunping Xia 1, Tong Wang 2, Lingmi Zhang 3 State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing
More informationOperational Modal Analysis of the Braga Sports Stadium Suspended Roof
Operational Modal Analysis of the Braga Sports Stadium Suspended Roof Filipe Magalhães 1, Elsa Caetano 2 & Álvaro Cunha 3 1 Assistant, 2 Assistant Professor, 3 Associate Aggregate Professor Faculty of
More informationIdentification of modal parameters from ambient vibration data using eigensystem realization algorithm with correlation technique
Journal of Mechanical Science and Technology 4 (1) (010) 377~38 www.springerlink.com/content/1738-494x DOI 107/s106-010-1005-0 Identification of modal parameters from ambient vibration data using eigensystem
More informationComparison of the Results Inferred from OMA and IEMA
Comparison of the Results Inferred from OMA and IEMA Kemal Beyen, Kocaeli University, Kocaeli, Turkey, kbeyen@kocaeli.edu.tr Mustafa Kutanis Sakarya University, Sakarya, Turkey, mkutanis@gmail.com.tr ABSTRACT:
More informationavailable online at
Acta Polytechnica CTU Proceedings 7:6 68, 7 Czech Technical University in Prague, 7 doi:./app.7.7.6 available online at http://ojs.cvut.cz/ojs/index.php/app EXPERIMENTAL IDENTIFICATION OF AN ELASTO-MECHANICAL
More informationMeasurement and Prediction of the Dynamic Behaviour of Laminated Glass
Paper 173 Measurement and Prediction of the Dynamic Behaviour of Laminated Glass Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures Technology, B.H.V.
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 informationEXPERIMENTAL AND THEORETICAL SYSTEM IDENTIFICATION OF FLEXIBLE STRUCTURES WITH PIEZOELECTRIC ACTUATORS
24 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES EXPERIMENTAL AND THEORETICAL SYSTEM IDENTIFICATION OF FLEXIBLE STRUCTURES WITH PIEZOELECTRIC ACTUATORS Aghil Yousefi-Koma*, David Zimcik* and Andrei
More informationDEVELOPMENT AND USE OF OPERATIONAL MODAL ANALYSIS
Università degli Studi di Napoli Federico II DEVELOPMENT AND USE OF OPERATIONAL MODAL ANALYSIS Examples of Civil, Aeronautical and Acoustical Applications Francesco Marulo Tiziano Polito francesco.marulo@unina.it
More informationAn Overview of Major Developments and Issues in Modal Identification
An Overview of Major Developments and Issues in Modal Identification Lingmi Zhang Institute of Vibration Engineering Nanjing University of Aeronautics & Astronautics Nanjing, 210016, CHINA, P.R. ABSTRACT
More informationGROUND VIBRATION TEST ON ULTRALIGHT PLANES
GOUND VBATON TEST ON ULTALGHT PLANES Dr. Eric Groß, Technische Universität Hamburg-Harburg, Germany Burkhard Martyn, Technische Universität Hamburg-Harburg, Germany e.gross@tuhh.de Abstract Modal parameters
More informationA comparison between modal damping ratios identified by NExT-ERA and frequency domain impact test
A comparison between modal damping ratios identified by NEx-ERA and frequency domain impact test Ali Nouri a, * and Sajjad Hajirezayi b a Department of Aerospace Engineering, Shahid Sattari Aeronautical
More informationSystem Theory- Based Iden2fica2on of Dynamical Models and Applica2ons
System Theory- Based Iden2fica2on of Dynamical Models and Applica2ons K. C. Park Center for Aerospace Structures Department of Aerospace Engineering Sciences University of Colorado at Boulder, CO, USA
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 informationA vector to a system output vector is required for many
VOL. 9, NO. 4 J. GUIDANCE JUL-AUGUST 986 Identifying Approximate Linear Models for Simple Nonlinear Systems Lucas G. Horta* and Jer-Nan Juangt NASA Langley Research Center, Hampton, Virginia This paper
More informationGeneralized Eigenvalue Decomposition in Time Domain Modal Parameter Identification
Wenliang Zhou 1 Lang Mekra Nth America, LLC, 101 Tillessen Bloulevard, Ridgeway, SC 9130 e-mail: zhouwl@gmail.com David Chelidze Department of Mechanical Engineering & Applied Mechanics, University of
More informationModal Testing and System identification of a three story steel frame ArdalanSabamehr 1, Ashutosh Bagchi 2, Lucia Trica 3
8th European Workshop On Structural Health Monitoring (EWSHM 2016), 5-8 July 2016, Spain, Bilbao www.ndt.net/app.ewshm2016 Modal Testing and System identification of a three story steel frame ArdalanSabamehr
More informationModal Parameter Identification from Responses of General Unknown Random Inputs Ibrahim, S. R.; Asmussen, J. C.; Brincker, Rune
Aalborg Universitet Modal Parameter Identification from Responses of General Unknown Random Inputs Ibrahim, S. R.; Asmussen, J. C.; Brincker, Rune Published in: Proceedings of the 14th International Modal
More informationModal parameter identification from output data only
MATEC Web of Conferences 2, 2 (215) DOI: 1.151/matecconf/21522 c Owned by the authors, published by EDP Sciences, 215 Modal parameter identification from output data only Joseph Lardiès a Institut FEMTO-ST,
More informationWind tunnel sectional tests for the identification of flutter derivatives and vortex shedding in long span bridges
Fluid Structure Interaction VII 51 Wind tunnel sectional tests for the identification of flutter derivatives and vortex shedding in long span bridges J. Á. Jurado, R. Sánchez & S. Hernández School of Civil
More informationNon-stationary Ambient Response Data Analysis for Modal Identification Using Improved Random Decrement Technique
9th International Conference on Advances in Experimental Mechanics Non-stationary Ambient Response Data Analysis for Modal Identification Using Improved Random Decrement Technique Chang-Sheng Lin and Tse-Chuan
More informationOperational modal analysis based prediction of actual stress in an offshore structural model
Downloaded from orbit.dtu.dk on: ug 24, 28 Operational modal analysis based prediction of actual stress in an offshore structural model Glindtvad Tarpø, Marius; Silva Nabuco, runa; Skafte, nders; Kristoffersen,
More informationStructural System Identification (KAIST, Summer 2017) Lecture Coverage:
Structural System Identification (KAIST, Summer 2017) Lecture Coverage: Lecture 1: System Theory-Based Structural Identification Lecture 2: System Elements and Identification Process Lecture 3: Time and
More informationEstimation of Rotational Degrees of Freedom by EMA and FEM Mode Shapes
Estimation of Rotational Degrees of Freedom by EMA and FEM Mode Shapes A. Sestieri, W. D Ambrogio, R. Brincker, A. Skafte, A. Culla Dipartimento di Ingegneria Meccanica e Aerospaziale, Università di Roma
More informationNonparametric identification of added masses in frequency domain: a numerical study
IPPT Reports 4f/2013 1 Nonparametric identification of added masses in frequency domain: a numerical study Grzegorz Suwała 1, Jan Biczyk 2, Łukasz Jankowski 1,* 1 Institute of Fundamental Technological
More informationDamage Localization under Ambient Vibration Using Changes in Flexibility
Damage Localization under Ambient Vibration Using Changes in Flexibility Yong Gao and B.F. Spencer, Jr. Ph.D. student, Department of Civil Engineering & Geological Science, University of Notre Dame, Notre
More informationApplied Modal Analysis of Wind Turbine Blades
Risø-R-1388(EN) Applied Modal Analysis of Wind Turbine Blades Henrik Broen Pedersen, Ole Jesper Dahl Kristensen Risø National Laboratory, Roskilde, Denmark February 2003 Abstract In this project modal
More informationResponse Spectrum Analysis Shock and Seismic. FEMAP & NX Nastran
Response Spectrum Analysis Shock and Seismic FEMAP & NX Nastran Table of Contents 1. INTRODUCTION... 3 2. THE ACCELEROGRAM... 4 3. CREATING A RESPONSE SPECTRUM... 5 4. NX NASTRAN METHOD... 8 5. RESPONSE
More informationABSTRACT Modal parameters obtained from modal testing (such as modal vectors, natural frequencies, and damping ratios) have been used extensively in s
ABSTRACT Modal parameters obtained from modal testing (such as modal vectors, natural frequencies, and damping ratios) have been used extensively in system identification, finite element model updating,
More informationDamage Identification in Wind Turbine Blades
Damage Identification in Wind Turbine Blades 2 nd Annual Blade Inspection, Damage and Repair Forum, 2014 Martin Dalgaard Ulriksen Research Assistant, Aalborg University, Denmark Presentation outline Research
More informationAutomated Estimation of an Aircraft s Center of Gravity Using Static and Dynamic Measurements
Proceedings of the IMAC-XXVII February 9-, 009 Orlando, Florida USA 009 Society for Experimental Mechanics Inc. Automated Estimation of an Aircraft s Center of Gravity Using Static and Dynamic Measurements
More informationModal Parameter Estimation from Operating Data
Modal Parameter Estimation from Operating Data Mark Richardson & Brian Schwarz Vibrant Technology, Inc. Jamestown, California Modal testing, also referred to as Experimental Modal Analysis (EMA), underwent
More informationAN INVESTIGATION INTO THE LINEARITY OF THE SYDNEY OLYMPIC STADIUM. David Hanson, Graham Brown, Ross Emslie and Gary Caldarola
ICSV14 Cairns Australia 9-12 July, 2007 AN INVESTIGATION INTO THE LINEARITY OF THE SYDNEY OLYMPIC STADIUM David Hanson, Graham Brown, Ross Emslie and Gary Caldarola Structural Dynamics Group, Sinclair
More information1330. Comparative study of model updating methods using frequency response function data
1330. Comparative study of model updating methods using frequency response function data Dong Jiang 1, Peng Zhang 2, Qingguo Fei 3, Shaoqing Wu 4 Jiangsu Key Laboratory of Engineering Mechanics, Nanjing,
More informationBLIND SOURCE SEPARATION TECHNIQUES ANOTHER WAY OF DOING OPERATIONAL MODAL ANALYSIS
BLIND SOURCE SEPARATION TECHNIQUES ANOTHER WAY OF DOING OPERATIONAL MODAL ANALYSIS F. Poncelet, Aerospace and Mech. Eng. Dept., University of Liege, Belgium G. Kerschen, Aerospace and Mech. Eng. Dept.,
More information2.0 Theory. 2.1 Ground Vibration Test
2.0 Theory The following section provides a comprehensive overview of the theory behind the concepts and requirements of a GVT (Ground Vibration Test), as well as some of the background knowledge required
More informationESTIMATION OF FREQUENCY RESPONSE FUNCTIONS BY RANDOM DECREMENT
ESTIMATION OF FREQUENCY RESPONSE FUNCTIONS BY RANDOM DECREMENT J.C. Asmussen & F Brincker Department of Building Technology and Structural Engineering Aalborg University, Sohngaardsholmsvej 57, 9000 Aalborg,
More informationMeasurement of wind-induced response of buildings using RTK-GPS and integrity monitoring
Measurement of wind-induced response of buildings using RTK-GPS and integrity monitoring Akihito Yoshida a,1, Yukio Tamura a, Masahiro Matsui a, Sotoshi Ishibashi b, Luisa Carlotta Pagnini c a Tokyo Institute
More informationViscoelastic vibration damping identification methods. Application to laminated glass.
Viscoelastic vibration damping identification methods. Application to laminated glass. J. Barredo a *, M. Soriano b, M. S. Gómez b, L. Hermanns' 5 "Centre for Modelling in Mechanical Engineering (CEMIMF2I2),
More informationA priori verification of local FE model based force identification.
A priori verification of local FE model based force identification. M. Corus, E. Balmès École Centrale Paris,MSSMat Grande voie des Vignes, 92295 Châtenay Malabry, France e-mail: corus@mssmat.ecp.fr balmes@ecp.fr
More informationMeasurement 44 (2011) Contents lists available at SciVerse ScienceDirect. Measurement
Measurement 44 (2011) 2112 2123 Contents lists available at SciVerse ScienceDirect Measurement journal homepage: www.elsevier.com/locate/measurement Enhancement of coherence functions using time signals
More informationAEROSPACE ENGINEERING
AEROSPACE ENGINEERING Subject Code: AE Course Structure Sections/Units Topics Section A Engineering Mathematics Topics (Core) 1 Linear Algebra 2 Calculus 3 Differential Equations 1 Fourier Series Topics
More informationFREQUENCY DOMAIN IDENTIFICATION TOOLBOX
NASA Technical Memorandum 109039 FREQUENCY DOMAIN IDENTIFICATION TOOLBOX Lucas G. Horta and Jer-Nan Juang NASA Langley Research Center Hampton, Virginia 23681-0001 Chung-Wen Chen North Carolina State University
More informationThe effect of environmental and operational variabilities on damage detection in wind turbine blades
The effect of environmental and operational variabilities on damage detection in wind turbine blades More info about this article: http://www.ndt.net/?id=23273 Thomas Bull 1, Martin D. Ulriksen 1 and Dmitri
More informationIDENTIFICATION OF THE MODAL MASSES OF AN UAV STRUCTURE IN OPERATIONAL ENVIRONMENT
IDENTIFICATION OF THE MODAL MASSES OF AN UAV STRUCTURE IN OPERATIONAL ENVIRONMENT M.S. Cardinale 1, M. Arras 2, G. Coppotelli 3 1 Graduated Student, University of Rome La Sapienza, mariosalvatore.cardinale@gmail.com.
More informationInvestigation of traffic-induced floor vibrations in a building
Investigation of traffic-induced floor vibrations in a building Bo Li, Tuo Zou, Piotr Omenzetter Department of Civil and Environmental Engineering, The University of Auckland, Auckland, New Zealand. 2009
More informationMODAL IDENTIFICATION OF STRUCTURES USING ARMAV MODEL FOR AMBIENT VIBRATION MEASUREMENT
MODAL IDENTIFICATION OF STRUCTURES USING MODEL FOR AMBIENT VIBRATION MEASUREMENT 72 C S HUANG SUMMARY A procedure is presented for evaluating the dynamic characteristics of structures from ambient vibration
More informationUsing SDM to Train Neural Networks for Solving Modal Sensitivity Problems
Using SDM to Train Neural Networks for Solving Modal Sensitivity Problems Brian J. Schwarz, Patrick L. McHargue, & Mark H. Richardson Vibrant Technology, Inc. 18141 Main Street Jamestown, California 95327
More informationSHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 1B. Damping
SHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 1B. Damping By Tom Irvine Introduction Recall the homework assignment from Unit 1A. The data.txt time history represented a rocket vehicle dropped from
More information8/13/2010. Applied System Identification for Constructed Civil Structures. Outline. 1. Introduction : Definition. 1.Introduction : Objectives
Outline Applied System Identification for Constructed Civil Structures Dionysius Siringoringo Bridge & Structure Lab, Univ of Tokyo July 22, 2010 Introduction Definition Objectives Scope Experimental Methods
More informationEMD-BASED STOCHASTIC SUBSPACE IDENTIFICATION OF CIVIL ENGINEERING STRUCTURES UNDER OPERATIONAL CONDITIONS
EMD-BASED STOCHASTIC SUBSPACE IDENTIFICATION OF CIVIL ENGINEERING STRUCTURES UNDER OPERATIONAL CONDITIONS Wei-Xin Ren, Department of Civil Engineering, Fuzhou University, P. R. China Dan-Jiang Yu Department
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 informationOperational Modal Parameter Estimation from Short Time-Data Series
Operational Modal Parameter Estimation from Short Time-Data Series A thesis submitted to the Division of Research and Advanced Studies of the University of Cincinnati in partial fulfillment of the requirements
More informationComparison of In-Flight Measured and Computed Aeroelastic Damping: Modal Identification Procedures and Modeling Approaches
doi: 1.528/jatm.v8i2.558 Comparison of In-Flight Measured and Computed Aeroelastic Damping: Modal Identification Procedures and Modeling Approaches Roberto da Cunha Follador 1, Carlos Eduardo de Souza
More informationAalborg Universitet. Damping Estimation by Frequency Domain Decomposition Brincker, Rune; Ventura, C. E.; Andersen, P.
Aalborg Universitet Damping Estimation by Frequency Domain Decomposition Brincker, Rune; Ventura, C. E.; Andersen, P. Published in: Proceedings of MAC 19 Publication date: 2001 Document Version Publisher's
More informationCOMPARATIVE STUDY OF OMA APPLIED TO EXPERIMENTAL AND SIMULATED DATA FROM AN OPERATING VESTAS V27 WIND TURBINE
COMPARATIVE STUDY OF OMA APPLIED TO EXPERIMENTAL AND SIMULATED DATA FROM AN OPERATING VESTAS V27 WIND TURBINE Oscar Ramírez 1, Dmitri Tcherniak 2, and Gunner Chr. Larsen 3 1 Research Assistant, DTU Wind
More informationOptimal System Realization in Frequency Domain
Optimal System Realization in Frequency Domain Jer-Nan Juang and Peiman G. Maghami NASA Langley Research Center Hampton, VA 23681 Abstract Several approaches are presented to identify an experimental system
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 informationOBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL
OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL Dionisio Bernal, Burcu Gunes Associate Proessor, Graduate Student Department o Civil and Environmental Engineering, 7 Snell
More informationGrandstand Terraces. Experimental and Computational Modal Analysis. John N Karadelis
Grandstand Terraces. Experimental and Computational Modal Analysis. John N Karadelis INTRODUCTION Structural vibrations caused by human activities are not known to be particularly damaging or catastrophic.
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 informationThe use of transmissibility properties to estimate FRFs on modified structures
Shock and Vibration 7 (00) 56 577 56 DOI 0./SAV-00-058 IOS Press The use of transmissibility properties to estimate FRFs on modified structures R.A.B. Almeida a,, A.P.V. Urgueira a and N.M.M. Maia b a
More informationAN EVALUATION METHOD OF MODAL DAMPLING AND ITS APPLICATION BASED ON A FVT OF CABLE STAYED BRIDGE
AN EVALUAION MEHOD OF MODAL DAMPLING AND IS APPLICAION BASED ON A FV OF CABLE SAYED BRIDGE Masataka NAKAMURA 1, Shuji UEDA 2, Jun-ichi SUZUMURA 3, Kazufumi HANADA 4, Makoto NAKAI 5 And Mitsuhide YOSHIDA
More informationDesign of Structures for Earthquake Resistance
NATIONAL TECHNICAL UNIVERSITY OF ATHENS Design of Structures for Earthquake Resistance Basic principles Ioannis N. Psycharis Lecture 3 MDOF systems Equation of motion M u + C u + K u = M r x g(t) where:
More informationThe Receptance Approach to Complex Operational Mode Normalization
he Receptance Approach to Complex Operational Mode Normalization Dionisio Bernal Northeastern University, Civil and Environmental Engineering Department, Center for Digital Signal Processing, Boston, MA
More informationExperimental Aerodynamics. Experimental Aerodynamics
Lecture 3: Vortex shedding and buffeting G. Dimitriadis Buffeting! All structures exposed to a wind have the tendency to vibrate.! These vibrations are normally of small amplitude and have stochastic character!
More informationUsing Operating Deflection Shapes to Detect Misalignment in Rotating Equipment
Using Operating Deflection Shapes to Detect Misalignment in Rotating Equipment Surendra N. Ganeriwala (Suri) & Zhuang Li Mark H. Richardson Spectra Quest, Inc Vibrant Technology, Inc 8205 Hermitage Road
More informationAutomated Frequency Domain Decomposition for Operational Modal Analysis
Autoated Frequency Doain Decoposition for Operational Modal Analysis Rune Brincker Departent of Civil Engineering, University of Aalborg, Sohngaardsholsvej 57, DK-9000 Aalborg, Denark Palle Andersen Structural
More informationAN CHARACTERIZATION OF JACKET-TYPE OFFSHORE STRUCTURE BY OP- ERATIONAL MODAL ANALYSIS
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm AN CHARACTERIZATION OF JACKET-TYPE OFFSHORE STRUCTURE BY OP- ERATIONAL MODAL ANALYSIS A.
More informationA hybrid system identification methodology for wireless structural health monitoring systems based on dynamic substructuring
A hybrid system identification methodology for wireless structural health monitoring systems based on dynamic substructuring Kosmas Dragos and Kay Smarsly Bauhaus University Weimar, Chair of Computing
More informationJ. A. Jurado, S. Hernández, A. Baldomir, F. Nieto, School of Civil Engineering of the University of La Coruña.
EACWE 5 Florence, Italy 19 th 23 rd July 2009 Flying Sphere image Museo Ideale L. Da Vinci Aeroelastic Analysis of Miradoiros Bridge in La Coruña (Spain) J. A. Jurado, S. Hernández, A. Baldomir, F. Nieto,
More informationOutput only modal analysis -Scaled mode shape by adding small masses on the structure
Master's Degree Thesis ISRN: BTH-AMT-EX--2008/D-07--SE Output only modal analysis -Scaled mode shape by adding small masses on the structure Sun Wei Department of Mechanical Engineering Blekinge Institute
More information