Research Article  Open Access
Scaling Mode Shapes in OutputOnly Structure by a MassChangeBased Method
Abstract
A masschangebased method based on outputonly data for the rescaling of mode shapes in operational modal analysis (OMA) is introduced. The mass distribution matrix, which is defined as a diagonal matrix whose diagonal elements represent the ratios among the diagonal elements of the mass matrix, is calculated using the unscaled mode shapes. Based on the theory of null space, the mass distribution vector or mass distribution matrix is obtained. A small mass with calibrated weight is added to a certain location of the structure, and then the mass distribution vector of the modified structure is estimated. The mass matrix is identified according to the difference of the mass distribution vectors between the original and modified structures. Additionally, the universal set of modes is unnecessary when calculating the mass distribution matrix, indicating that modal truncation is allowed in the proposed method. The massscaled mode shapes estimated in OMA according to the proposed method are compared with those obtained by experimental modal analysis. A simulation is employed to validate the feasibility of the method. Finally, the method is tested on outputonly data from an experiment on a fivestorey structure, and the results confirm the effectiveness of the method.
1. Introduction
Operational modal analysis (OMA) is an engineering field in which the modal properties of structures under ambient vibrations are studied [1]. Its development can be traced to 1965, when the method using the crosscorrelation function to estimate the frequency response function (FRF) was proposed by Clarkson and Mercer [2], which was also the origin of using the crosscorrelation function instead of the impulse response function (IRF) to identify the mode parameters in cases where the excitation is unknown. After approximately 50 years of development, especially due to studies of autoregressive moving average (ARMA) [3], random decrement technique (RDT) [4], operational deflection shapes (ODS) [5], natural excitation technique (NExT) [6], and stochastic subspace identification (SSI) [7], OMA has matured and has been extended to a wider range of applications in engineering.
Furthermore, if the identified mode parameters are further used to predict structural response, structural modification, or health monitoring, the FRF or IRF should be identified. In fact, the FRF is obtained using experimental modal analysis (EMA) in the laboratory. Since the boundary conditions of the structure during a laboratory test are not entirely identical with real inoperation working conditions of the structure, the FRF of the system in the laboratory is different from that in an operating state. Therefore, reconstructing the FRF based on the mode parameters (natural frequencies, damping ratios, and massscaled mode shapes) identified using the OMA method is necessary. Because the excitation cannot be measured in OMA, the mode shapes cannot be normalized or scaled by mass. As a consequence, an additional procedure to obtain the massscaled modal shapes is needed. In the past few years, many methods have been proposed to solve this problem. Doebling and Farrar [8] used a finiteelement model (FEM) mass matrix that is reduced to the measurement degree of freedom (DOF) to normalize the mode shapes, according to the work of Guyan [9]. However, their method is only valid for the lowerfrequency modes. Another method used the mode parameters identified using the OMA technology to upgrade the FEM, from which the mass matrix that scales the mode shapes was obtained [10–15]. Randall et al. [16] strictly limited the form of excitation and developed cepstral methods for updating mode parameters identified using the OMA method when the logarithmic spectrum of the excitation is smooth and flat enough. Using the mass perturbation as a mode shape scaling method originated from Parloo et al. [17, 18]. By adding one or more masses whose weights are calibrated to the structure, the operational mode shapes can be massscaled by the shift in mode frequencies between the original and modified structures. On the basis of this idea, much work has been done. The approximate formula for determination of the scaling factors based on the frequency shift when introducing mass perturbation on the structure was derived directly from the governing equation of motion [19]. It explained how testing should be performed in order to significantly reduce approximation errors. The uncertainties of the estimated scaling factors were also investigated [20]. An exact formulation computing the square of the scaling factors was given by Bernal [21, 22]. This formulation can be used in structures with both nonclosely spaced modes and closely spaced modes. LópezAenlle et al. [23] provided two new exact formulations for scaling factors, illustrated the large difference in accuracy between the approximate and exact formulations, and confirmed that the results are similar only when the mass perturbation is relatively small. Taking into account a drawback that the first mode of a structure is less sensitive to the mass perturbation, a new massstiffness perturbation method was suggested to scale the mode shapes [24]. A procedure to optimize the mass perturbation strategy was proposed by LópezAenlle et al. [25, 26], which used the mode parameters (natural frequencies and shapes) of the original structure as basic information. Different mass perturbation methods were applied to estimate the scaling factors of a 15T concrete slab by Fernández et al. [27]. The results showed that all of the methods estimate the scaling factors with good accuracy when an appropriate mass change strategy is adopted. Poozesh et al. employed the drive point scaling method to scale optically measured operating deflection shapes (ODS) [28]. Additionally, they demonstrated the benefits and drawbacks of the mass sensitivity technique. However, the accuracy of the mass matrix obtained by the finiteelement method commonly depends on the accuracy of the modeling of the structure. Meanwhile, the mass change in the aforementioned method should not be too high in order to minimize the difference of the mode shapes between the original and modified structures.
In this paper, the mass distribution matrix is defined and calculated using unscaled mode shapes given that the structure is discrete and the mass matrix is diagonal [8, 9, 29]. The mass matrix is the product of the mass distribution matrix and a coefficient. Then, a small mass with calibrated weight is added to a certain location of the structure, and the mass distribution vector of the structure is estimated. It is noteworthy that removing a small mass with calibrated weight from a certain location of the structure is also feasible. Based on the difference of the mass distribution vectors between the original and modified structures, the mass matrix is identified. Mode shapes are scaled by the mass matrix. It is also worth noting that a universal set of modes is unnecessary when calculating the mass distribution matrix. This means that modal truncation is allowed in the proposed method.
The rest of this paper is organized as follows. In Section 2, the mass distribution matrix of the structure is defined and calculated. The mass matrix is obtained based on the method that a small mass with calibrated weight is added to a certain DOF of the structure. Then, the unscaled mode shapes are scaled by the mass matrix. Finally, a simulation and an experiment are presented in Sections 3 and 4, respectively.
2. Scaling Mode Shapes
In this section, the mass distribution matrix is defined and calculated by unscaled mode shapes. Then, the mass matrix is obtained by the mass change method and the mode shapes are scaled by the mass matrix. Finally, modal truncation is discussed.
2.1. Mass Distribution Matrix
The equation of motion of a linear timeinvariant discrete system with DOFs can be written aswhere , , and represent the mass, damping, and stiffness matrices, respectively. is a real, symmetric, positive definite matrix. and are real, symmetric, and positive or semipositive definite matrices. denotes time. is the excitation, and represents the displacement response vector.
According to the property that mode shapes are weighted orthogonal with respect to the mass matrix, one haswhere is the massscaled mode shape matrix and is the identity matrix.
The relation between the unscaled or arbitrary scaled mode shape matrix and the massscaled mode shape matrix can be expressed aswhere is the unscaled or arbitrarily scaled mode shape matrix. The scaling factor is a diagonal matrix with unknown but positive elements.
Substituting (3) into (2) yields
On the basis of the assumption that the structure is a discrete system whose mass matrix can be assumed to be diagonal [8, 9, 29], the expanded form of (4) can be written aswhereand returns a square diagonal matrix with the elements of vector on the main diagonal.
Considering the symmetry on the lefthand side of (5), equations can be obtained from (5), in which is the independent variable.
When , equations can be extracted from (5) as
Equation (7) cannot be solved since the scaling factor matrix is unknown.
When , equations can be extracted from (5) aswhere
Owing to the symmetry on the lefthand side of (5), the case when is similar to that when and is omitted here. In (8), there are equations. is the null space of . Given that exists and is unique, the dimension of the null space of is 1. Therefore, the rank of is . The mass distribution vector , which denotes the ratios among the elements of , is the unique basis of the null space of . Since singular value decomposition (SVD) is an efficient method of calculating the null space of a matrix, it is adopted to obtain :
Because the rank of is and the singular values in are sorted in descending order, can be obtained aswhere is scaled by its 2norm because is a unit orthogonal matrix. Consequently, the mass distribution matrix is expressed as
2.2. Mass Matrix
A small mass with calibrated weight is added to the DOF of the structure. According to the method similar to that used in Section 2.1, the mass distribution vector of the structure modified by , denoted as , can be obtained.
Owing to the fact that is scaled by its 2norm, can be rewritten aswhere is the 2norm of vector .
Similarly, can be rewritten as
Combining the kth equation in (13) and the kth equation in (14), the following equations can be obtained:
Solving the equations in (15), one obtains
Substituting (16) into (6) and (13) yields
Then, the mass matrix of the structure is obtained as
2.3. Scaling Factor
Since the mass matrix has been obtained in Section 2.2, the mode shape matrix can be scaled. Substituting (18) into (4), one has
Since the matrix is diagonal and all diagonal elements are positive, can thus be written aswhere represents the square root of the diagonal element of matrix .
Substituting (20) into (3), one can scale the mode shapes by the mass matrix:
2.4. Modal Truncation
Generally, a universal set of modes cannot be obtained due to modal truncation. When the order of the modes is truncated as (), (8) becomes
Considering the fact that the rank of is , both and should satisfy the following inequality:
Solving the inequality in (23), one haswhere means rounding to the nearest integer greater than or equal to and is the smallest positive integer that satisfies the inequality in (24). The details of the relationship between and are presented in Figure 1.
When , as is shown in Figure 1, . This means that the universal set of modes is unnecessary when calculating the mass distribution matrix. In other words, the modal truncation is allowed. The growth rate of is much smaller than the growth rate of . When is sufficiently large, it follows that .
It is worth noting that (24) is the premise of calculating the mass distribution matrix in this paper. If the total number of natural modes identified is smaller than , the mass matrix or mass distribution matrix cannot be obtained by the aforementioned method.
Considering the modal truncation, when arbitrary mode shapes are taken from all mode shapes, (22) can be written aswhere ( is the total number of combinations of taking items at a time from ). According to (11), mass distribution vectors can be obtained.
Similarly, when arbitrary mode shapes are taken from all ones, (22) can be written aswhere . Thus, mass distribution vectors can be extracted.
Inductively, when mode shapes are taken, mass distribution vector can be obtained.
Using the enumerating method, mass distribution vectors can be obtained from identified modes; namely,
Theoretically, the mass distribution vectors calculated by different mode shape combinations should be identical.
3. Simulation
In this section, a 10DOF structure is designed to verify the method of scaling mode shapes.
The mass matrix is
The stiffness matrix is
Consequently, the natural mode frequencies are 8.34, 18.92, 26.67, 33.74, 42.02, 48.69, 54.57, 70.05, 116.68, and 122.89 Hz. Taking the damping ratio for every order of the structure to be , , , , , , , , , and , the damping matrix of the system is computed bywhere , , and .
3.1. Mass Distribution Matrix
The loads on the structure are uncorrelated random signals subjected to Gaussian distribution. In this simulation, they fit , , , , , , , , , and , where denotes a Gaussian distribution with being a mean value and the variance. In order to illustrate the advantage of modal truncation, all of the random signals are filtered by Chebyshev Type I filters, and the passband edge frequency Hz. The sampling frequency Hz. All responses of the structure are recorded.
The modal parameters , , and are identified by the SSI method. Owing to the bandlimited white noise, only the first seven modes can be obtained. According to (11) and (27), the 29 mass distribution vectors or the mass distribution matrices of the structure can be calculated by the unscaled mode shapes . In practice, according to (24), . Therefore, at least five mode shapes of the structure are required to calculate the mass distribution vector or mass distribution matrix. A different total number of identified mode shapes can be used to calculate the mass distribution vectors, the details of which are shown in Table 1.

The mass distribution vectors calculated by the first seven modes are illustrated in Figure 2. Apparently, the given mass distribution vector is
In Figure 2, the dark, thin solid lines coincide well with the thick dashed lines. Calculating the average of the 29 results, the identified mass distribution vector is
The consistency of the lines means that the method used to calculate the mass distribution matrix is reliable.
3.2. Mass Matrix and MassScaled Mode Shapes
To identify the mass matrix of the structure, a small mass, kg, is added to DOF 9 of the structure. According to a method similar to that used in Section 3.1, can be obtained as
According to (17) and (18), and can be acquired. The relative error between the given mass matrix and the identified mass matrix is shown in Table 2.

As can be seen from Table 2, the maximum relative error of the identified mass is 4.69% (DOF 5). Therefore, the relative error between the given mass matrix and the identified mass matrix is acceptable.
Since the mass matrix was obtained in Section 3.2, the mode scaling factor can be computed according to (20) and is given by
According to (21) or (3), the massscaled mode shapes can be obtained.
4. Experimental Verification
To verify the aforementioned theory, a fivestorey prototype structure was constructed for experimentation and is shown in Figure 3.
As can be seen from Figure 3, in this experiment, is white noise. The acceleration responses are measured and stored by the SO Analyze© system (m+p international, Germany). The sampling frequency Hz.
The mode parameters are identified by the SSI method, and the stability diagram is reproduced in Figure 4. The natural frequency and the damping ratio are shown in Table 3.

Five modes can be identified. However, in Figure 4, the poles at the first modes are not as stable as the others. Meanwhile, the coherence functions at the first modes are not as satisfied as at the other modes. One reason for this is the fact that the responses of the first mode are larger than that of the others. This is also because the shaker in our laboratory is not good enough in the low frequency band. According to (24), . Therefore, at least four modes shapes of the structure are required to calculate the mass distribution vector or mass distribution matrix. According to the stability diagram in Figure 4, the second, third, fourth, and fifth mode shapes are employed to calculate the mass distribution vector that is shown in Figure 5(a).
(a) Original structure
(b) Modified structure
In Figure 5(a), the mass distribution vector can be identified as
To identify the mass matrix of the structure, a small mass, kg, is added to DOF 3 of the structure. Repeated modal testing after changing the mass is undertaken. Similarly, the mass distribution vector of the modified structure can be computed, with the results shown in Figure 5(b). Then, can be obtained as
According to (17) and (18), and can be acquired, where
According to (21) or (3), massscaled mode shapes can be obtained, which are shown in Figure 6. For comparison, the massscaled mode shapes estimated by EMA serve as the control group.
(a) Hz
(b) Hz
(c) Hz
(d) Hz
(e) Hz
In Figure 6, the massscaled mode shapes obtained using EMA are represented by the light, thick, and dashed lines, while those identified using the proposed method are represented by the dark, thin, and solid lines. The dark, thin, and solid lines coincide with the light, thick, and dashed lines. The consistency of the dashed and solid lines means that the relative error between the massscaled mode shapes identified using EMA and those identified using the method proposed in this paper is acceptable.
5. Conclusions
In this paper, a masschangebased method for the rescaling of mode shapes in OMA is introduced, which only requires response data. The features of the proposed method are listed as follows.
The mass distribution matrix, which is defined as a diagonal matrix whose diagonal elements represent the ratios among the diagonal elements of the mass matrix, was calculated using the unscaled mode shapes given that the structure is discrete and the mass matrix is diagonal. Based on the theory of null space, the mass distribution vector or matrix was obtained.
A small mass with calibrated weight was added to a certain DOF of the structure, and then the mass distribution vector of the modified structure was estimated. Based on the difference of the mass distribution vectors between the original and modified structures, the mass matrix was identified. In this method, the mass with calibrated weight was not necessary to be of light weight in order to minimize the difference between the original and modified structures.
A universal set of modes was unnecessary when calculating the mass distribution matrix, indicating that modal truncation is allowed in the proposed method.
A simulation was employed to validate the feasibility of the method. Furthermore, the method was tested on the outputonly data from a fivestorey structure under laboratory conditions. The massscaled mode shapes estimated by OMA according to the method proposed in this paper were compared with those obtained using EMA. The effectiveness of the method was confirmed and validated.
This method is not only enforceable in the simulation and experimental cases proposed in this paper but is also applicable for other engineering problems.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the National Natural Science Foundation of China under Grant no. 11272235.
References
 R. Brincker, “Some elements of operational modal analysis,” Shock and Vibration, vol. 2014, Article ID 325839, 11 pages, 2014. View at: Publisher Site  Google Scholar
 B. L. Clarkson and C. A. Mercer, “Use of cross correlation in studying the response of lightly damped structures to random forces,” AIAA Journal, vol. 3, no. 12, pp. 2287–2291, 1965. View at: Publisher Site  Google Scholar
 H. Akaike, “Power spectrum estimation through autoregressive model fitting,” Annals of the Institute of Statistical Mathematics, vol. 21, pp. 407–419, 1969. View at: Publisher Site  Google Scholar  MathSciNet
 S. R. Ibrahim, “Random decrement technique for modal identification of structures,” Journal of Spacecraft and Rockets, vol. 14, no. 11, pp. 696–700, 1977. View at: Publisher Site  Google Scholar
 O. Døssing and C. H. Staker, “Operational deflection shapes: background measurement and application,” in Proceedings of the 5th International Modal Analysis Conference, pp. 1372–1378, London, UK, 1987. View at: Google Scholar
 G. H. James III, T. G. Carne, and J. P. Lauffer, “The natural excitation technique (NExT) for modal parameter extraction from operating structures,” Modal Analysis—The International Journal of Analytical and Experimental Modal Analysis, vol. 4, pp. 260–277, 1995. View at: Google Scholar
 B. Peeters, G. De Roeck, and T. Pollet, “Stochastic subspace techniques applied to parameter identification of civil engineering structures,” in Proceedings of the “New Advances in Modal Synthesis of Large Structures: Nonlinear Damped and Nondeterministic Cases”, pp. 151–162, Lyon, France, 1995. View at: Google Scholar
 S. W. Doebling and C. R. Farrar, “Computation of structural flexibility for bridge health monitoring using ambient modal data,” in Proceedings of the 11th Conference on Engineering Mechanics, pp. 1114–1117, Fort Lauderdale, Fla, USA, May 1996. View at: Google Scholar
 R. J. Guyan, “Reduction of stiffness and mass matrices,” AIAA Journal, vol. 3, no. 2, p. 380, 1965. View at: Publisher Site  Google Scholar
 C. E. Ventura, J. F. Lord, M. Turek, R. Brincker, P. Andersen, and E. Dascotte, “FEM updating of tall buildings using ambient vibration data,” in Proceedings of the 6th International Conference on Structural Dynamics (EURODYN '05), pp. 237–242, Paris, France, 2005. View at: Google Scholar
 H. Shahverdi, C. Mares, and J. E. Mottershead, “Model structure correction and updating of aeroengine casings using fictitious mass modifications,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 219, no. 1, pp. 19–30, 2005. View at: Publisher Site  Google Scholar
 B. Schwarz and M. Richardson, “Using FEA modes to scale experimental mode shapes,” in Proceedings of the 24th Conference and Exposition on Structural Dynamics, Saint Louis, Mo, USA, February 2006. View at: Google Scholar
 B. Jaishi, H.J. Kim, M. K. Kim, W.X. Ren, and S.H. Lee, “Finite element model updating of concretefilled steel tubular arch bridge under operational condition using modal flexibility,” Mechanical Systems and Signal Processing, vol. 21, no. 6, pp. 2406–2426, 2007. View at: Publisher Site  Google Scholar
 D. Hanson, R. B. Randall, J. Antoni, T. P. Waters, D. J. Thompson, and R. A. J. Ford, “Cyclostationarity and the cepstrum for operational modal analysis of MIMO systems—part II: obtaining scaled mode shapes through finite element model updating,” Mechanical Systems and Signal Processing, vol. 21, no. 6, pp. 2459–2473, 2007. View at: Publisher Site  Google Scholar
 M. L. Aenlle and R. Brincker, “Modal scaling in operational modal analysis using a finite element model,” International Journal of Mechanical Sciences, vol. 76, pp. 86–101, 2013. View at: Publisher Site  Google Scholar
 R. B. Randall, Y. Gao, and J. Swevers, “Updating modal models from response measurements,” in Proceedings of the 23rd International Conference on Noise and Vibration Engineering, pp. 1153–1160, Leuven, Belgium, 1999. View at: Google Scholar
 E. Parloo, P. Verboven, P. Guillaume, and M. Van Overmeire, “Sensitivitybased operational mode shape normalisation,” Mechanical Systems and Signal Processing, vol. 16, no. 5, pp. 757–767, 2002. View at: Publisher Site  Google Scholar
 E. Parloo, P. Verboven, P. Guillaume, and M. Van Overmeire, “Force identification by means of inoperation modal models,” Journal of Sound and Vibration, vol. 262, no. 1, pp. 161–173, 2003. View at: Publisher Site  Google Scholar
 R. Brincker and P. Andersen, “A way of getting scaled mode shapes in output only modal testing,” in Proceedings of the 21th International Modal Analysis Conference, Orlando, Fla, USA, 2003. View at: Google Scholar
 R. Brincker, J. Rodrigues, P. Andersen, and S. V. S. Aps, “Scaling the mode shapes of a building model by mass changes,” in Proceedings of the 22th International Modal Analysis Conference, Dearborn, Mich, USA, 2004. View at: Google Scholar
 D. Bernal, “Modal scaling from known mass perturbations,” Journal of Engineering Mechanics, vol. 130, no. 9, pp. 1083–1088, 2004. View at: Publisher Site  Google Scholar
 D. Bernal, “A receptance based formulation for modal scaling using mass perturbations,” Mechanical Systems and Signal Processing, vol. 25, no. 2, pp. 621–629, 2011. View at: Publisher Site  Google Scholar
 M. LópezAenlle, R. Brincker, F. Pelayo, and A. F. Canteli, “On exact and approximated formulations for scalingmode shapes in operational modal analysis by mass and stiffness change,” Journal of Sound and Vibration, vol. 331, no. 3, pp. 622–637, 2012. View at: Publisher Site  Google Scholar
 M. M. Khatibi, M. R. Ashory, and A. Malekjafarian, “Scaling of mode shapes using massstiffness change method,” in Proceedings of the 3rd International Operational Modal Analysis Conference, pp. 699–706, Portonovo, Italy, 2009. View at: Google Scholar
 M. LópezAenlle, P. Fernández, R. Brincker, and A. FernándezCanteli, “Erratum to ‘Scalingfactor estimation using an optimized masschange strategy’ [Mech. Syst. Signal Process. 24 (5)(2010) 1260–1273],” Mechanical Systems and Signal Processing, vol. 24, pp. 3061–3074, 2010. View at: Google Scholar
 M. L. Aenlle, R. Brincker, A. F. Canteli, and L. M. Villa García, “Scaling factor estimation by the mass change method,” in Proceedings of the 1st International Operational Modal Analysis Conference (IOMAC '05), pp. 53–64, Copenhagen, Danmark, April 2005. View at: Google Scholar
 P. Fernández, P. Reynolds, and M. López Aenlle, “Experimental evaluation of mass change approaches for scaling factors estimation,” in Proceedings of the 28th IMAC, A Conference on Structural Dynamics, pp. 109–118, Springer, February 2010. View at: Google Scholar
 P. Poozesh, D. D. Sabino, J. Baqersad, P. Avitabile, and C. Niezrecki, Practical Techniques for Scaling of Optically Measured Operating Deflection Shapes, Springer, Berlin, Germany, 2016.
 D. Bernal and B. Gunes, “Damage localization in outputonly systems: a flexibility based approach,” in Proceedings of the 20th International Modal Analysis Conference (IMAC '02), pp. 1185–1191, Los Angeles, Calif, USA, February 2002. View at: Google Scholar
Copyright
Copyright © 2017 Liangliang Yu and Hanwen Song. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.