# Modal Participation Estimated from the Response Correlation Matrix

Authors: Rune Brincker; Sandro D. R. Amador; Martin Juul; Manuel Lopez-Aenelle
Journal: Shock and Vibration
 (2019)

---

## Abstract

In this paper, we are considering the case of estimating the modal participation vectors from the operating response of a structure. Normally, this is done using a fitting technique either in the time domain using the correlation function matrix or in the frequency domain using the spectral density matrix. In this paper, a more simple approach is proposed based on estimating the modal participation from the correlation matrix of the operating responses. For the case of general damping, it is shown how the response correlation matrix is formed by the mode shape matrix and two transformation matricesT1 and T1 that contain information about the modal parameters, the generalized modal masses, and the input load spectral density matrix Gx. For the case of real mode shapes, it is shown how the response correlation matrix can be given a simple analytical form where the corresponding real modal participation vectors can be obtained in a simple way. Finally, it is shown how the real version of the modal participation vectors can be used to synthesize empirical spectral density functions.


---

## Body

## 1. Introduction

Modal participation is a classic subject in traditional experimental modal analysis where the modal parameters are obtained from an experimentally estimated frequency response function (FRF) matrix. Recent examples of FRF-based modal participation are found in [1–3].In the case of operational modal analysis, where the modal parameters cannot be estimated based on the FRF because the excitation forces are unknown, the modal participation has to be estimated solely from the operating responses. This means using either the correlation function (CF) matrix if the modal identification is done in the time domain or using the spectral density (SD) matrix if the modal identification is done in the frequency domain.It is a common practice in modal analysis to estimate solely the natural frequencies, damping ratios, and the corresponding mode shape vectors from the vibration measurements. Though the modal mass and modal participation factors are normally not estimated, they also contribute to describe the modal behavior of the tested structures. In [4], it is described how the model participation factor can be used to quantify the participation and, therefore, the importance of each identified mode to the total response of the tested structures.Another fact that justifies the necessity of estimating the modal participation factors is that it is always a good idea to synthesize the CF and/or the SD function from the estimated modal parameters in order to check the accuracy of these estimates. Furthermore, the state-of-art modal parameter estimation techniques such as those based on maximum likelihood [5] can only be used to estimate/optimize the modal parameters if a good guess of the poles, mode shapes, and modal participation factors are known a priori as discussed, for instance, in [6].Unlike in classic modal analysis, it is not possible to estimate the modal mass in operational modal analysis. On the other hand, it is perfectly possible to estimate the modal participation factors from output-only vibration measurements. Both in classic and in operational modal analysis, the estimation of the modal participation factors is carried in a linear least squares sense when the poles and mode shape vectors are estimated from the vibration data, as described, for instance, in [7, 8].In this paper, another approach is proposed to estimate the modal participation factors from the response correlation matrix. The main advantage of such approach with regard to the traditional least squares based ones is that the former estimates the modal participation factors by means of a closed expression. In order to illustrate the performance of the proposed approach from a practical point of view, it is applied to two application examples in the final part of the paper.


## 2. Theory

Before proceeding to the new idea that we like to propose in this paper, we will shortly introduce the classic ways of obtaining the modal participation from the operating responses. Considering the response vectoryt∈ℝNo containing the No responses of a linear system at time t, the correlation function (CF) matrix, Rτ∈ℝNo×No, can be defined in the classic way [9] as(1)Rτ=EytyTt+τ,with τ∈ℝ denoting a time lag. Assuming M modes with mode shape vectors bn∈ℂNo, poles λn∈ℂ, and a constant input load spectral density (SD) matrix Gx∈ℂNo×No of the load vector xt∈ℝNo containing No inputs, it can be shown (refer to, for instance, [4, 10]) that the analytical solution for the CF matrix is given by(2)Rτ=2π∑n=1MAnGxBne−λnτ+An∗GxBn∗e−λn∗τ,τ≤0,2π∑n=1MBnGxAneλnτ+Bn∗GxAn∗eλn∗τ,τ≥0,where ⋅∗ denotes the conjugate of a complex quantity and An∈ℂNo×No are the rank one residues defined by the mode shapes and the generalized modal mass an:(3)An=bnbnTan.The matricesBn∈ℂNo×No are a weighted sum of the residual matrices As and As∗ given by(4)Bn=∑s=1MAs−λn−λs+As∗−λn−λs∗.In [4, 10], the modal participation vector, γn∈ℂNo, is defined as(5)γn=BnGxbnan.This eventually leads to the following simplified expressions for the theoretical CF matrix solution:(6)Rτ=2π∑n=1MbnγnTe−λnτ+bn∗γnHe−λn∗τ,τ≤0,2π∑n=1MγnbnTeλnτ+γn∗bnHeλn∗τ,τ≥0,where ⋅H denotes the conjugate transpose (Hermitian) of a complex matrix. In [4], it is proposed to estimate the modal participation vectors, γn, either in the time domain fitting the empirical CF function matrices to the theoretical solutions (2) and (6) or in the frequency domain by fitting the empirical SD matrices to the similar theoretical solution for the SD matrix of the response, Gyω∈ℂNo×No, where Gyω is the Fourier transform of Rτ. This normally would include either the positive or the negative part of the CF matrix and thus leads to fitting on the half spectral density function matrix.When the modal parameters including the modal participation vectors are known, the corresponding response SD matrix can easily be found from equation (6) by taking the Fourier transform of both the negative and the positive time part of the CF matrix solution [4]:(7)Gyω=∑n=1MbnγnT−iω−λn+bn∗γnH−iω−λn∗+γnbnTiω−λn+γn∗bnHiω−λn∗.However, since this solution contains both terms of the typebnγnT and their respective transpose of the type γnbnT, there is no simple known way to fit the classic SD function matrix to find the modal participation vectors. Therefore, in modal identification with modern parametric frequency-domain identification techniques such as the polyreference least squares frequency-domain (pLSCF) [11, 12], the fitting is carried out by making use of the half spectrum (i.e, the SD computed solely from the CF function matrix with positive time lags τ) which consists only of the last two terms in equation (7).However, in this paper a more simple approach is proposed based on the response correlation matrix. Since the response correlation matrix is the value of the CF matrixRτ at time lag zero, but also equal to the integral of the SD matrix Gyω as expressed by the Parseval equation(8)C=R0=∫−∞+∞Gyωdω,the approach can be considered as a common approach for both time domain and frequency domain. In the following section, we will show that the solution given by equations (2) and (6) leads to an expression for the response correlation matrix that is real and symmetric. Furthermore, we will formulate a simple way to estimate the modal participation vectors from the correlation matrix, and finally we will illustrate how this information can be used to synthesize the SD matrix in case of random response data.


### 2.1. General Form of the Response Correlation Matrix

It is obvious that since the correlation matrixC=EytyTt can now be calculated in two different ways using equations (2) and (6), C+=R0+ and C_=R0_, we will start showing that the two different ways of calculation leads to the same result.From the bottom part of equation (2), we have(9)C+=2π∑n=1MBnGxAn+Bn∗GxAn∗.Substituting equation (4) in equation (9) and carrying out the multiplications, we get(10)C+=2π∑n=1M∑s=1MAs−λn−λs+As∗−λn−λs∗GxAn+∑s=1MAs∗−λn∗−λs∗+As−λn∗−λsGxAn∗=2π∑n=1M∑s=1MAs−λn−λsGxAn+As−λn∗−λsGxAn∗+As∗−λn∗−λs∗GxAn∗+As∗−λn−λs∗GxAn.Swapping the summation indexes in equation (10) and contracting the equation back to short form yields(11)C+=2π∑s=1M∑n=1MAsGxAn−λs−λn+AsGxAn∗−λs−λn∗+As∗GxAn∗−λs∗−λn∗+As∗GxAn−λs∗−λn=2π∑s=1MAsGx∑n=1MAn−λs−λn+An∗−λs−λn∗+As∗Gx∑n=1MAn∗−λs∗−λn∗+As−λs∗−λn=2π∑s=1MAsGxBs+As∗GxBs∗,which actually is equal to what we get from the top part of equation (2) for R0−, meaning that the two formulations C+=R0+ and C−=R0− are identical.In order to get a little closer to the structure of this theoretical solution for the response correlation matrix, we continue manipulating equations (9) and (10). By substituting equation (3) in these equations and renaming the summation indexes (s⟶r and n⟶c), we obtain(12)C+=2π∑r,c=1MbrbrTGxbcar−λr−λcacbcT+br∗brHGxbcar∗−λr∗−λcacbcT+brbrTGxbc∗ar−λr−λc∗ac∗bcH+br∗brHGxbc∗ar∗−λr∗−λc∗ac∗bcH,or in a more compact form(13)C+=2π∑r,c=1Mbrt1rcbcT+br∗t2rcbcT+brt3rcbcH+br∗t4rcbcH,with tirc=tir,c∈ℂ∀i∈1,2,3,4 given by(14)t1rc=brTGxbcar−λr−λcac,t2rc=brHGxbcar∗−λr∗−λcac,t3rc=brTGxbc∗ar−λr−λc∗ac∗,t4rc=brHGxbc∗ar∗−λr∗−λc∗ac∗.As a side remark, it is worth noticing that we without any loss of generality can remove the generalized modal massesan from the equation using the mass scaled mode shapes bn/an. However, since we later will use mode shapes with arbitrary scaling, we will keep the explicit modal masses in the equations.Considering the first term in the summation of equation (12), we see that the scalar t1rc=brTGxbc/ar−λr−λcac can be seen as a matrix element with row index r and column index c and that we can write the summation as the outer product BT1BT where B is the mode shape matrix B=b1b2⋯bM∈ℂNo×M and(15)T1=t1rc=t1r,c=t11,1⋯t11,M⋮⋱⋮t1M,1⋯t1M,M=brTGxbcar−λr−λcac∈ℂM×Mis the matrix containing the elements t1rc∀r,c∈1,2,…,M. Using this notation on all four terms of equation (12), we get(16)C+=2πBT1BT+B∗T2BT+BT3BH+B∗T4BH,where(17)T2=brHGxbcar∗−λr∗−λcac∈ℂM×M,T3=brTGxbc∗ar−λr−λc∗ac∗∈ℂM×M,T4=brHGxbc∗ar∗−λr∗−λc∗ac∗∈ℂM×M.Evaluating the bottom part of equation (6) for τ=0, we now obtain a similar relation expressed in the modal participation vectors:(18)C+=2π∑n=1MγnbnT+γn∗bnH=2πΓBT+Γ∗BH,where Γ=γ1γ2⋯γM∈ℂNo×M is the modal participation vector matrix. Since the two first terms of equation (16) come from the first term of equation (9) and thus equal ΓBT, we can conclude that the modal participation matrix can be expressed by(19)Γ=BT1+B∗T2,and we see that the last term in equation (18) Γ∗BH equals the last two terms in equation (16).Taking origin in the top part of equation (2) instead of equation (6) and using the definition of the modal participation matrix given by equation (19) leads to satisfaction of the equation C−=2πBΓT+B∗ΓH, which is the version of equation (18) corresponding to the top part of equation (6).Finally, we can conclude that equation (19) is a general result for the modal participation matrix and that our final and general result for the response correlation matrix is(20)Γ=2πBT1+B∗T2BT+B∗T1∗+BT2∗BH.


### 2.2. Correlation Matrix for Proportional Damping

In case of real modes, equation (20) becomes(21)C=4πBReTBT,where Re⋅ denotes the real value of a complex quantity and(22)T=T1+T2.For proportional damping, the generalized modal mass becomes pure imaginary and is given by(23)an=2iωdnmn,where ωdn is the damped angular natural frequency and mn is the classic modal mass. Furthermore, it is known (as seen, for instance, in [7]) that from the orthogonality equations, we have in the case of general damping(24)λr+λcmrc+crc=0,r≠c,2λnmn+cn,r=c.In the case of proportional damping, the full modal mass and damping matricesmrc and crc in equation (24) reduce to the corresponding diagonal forms mn and cn. Combining T1 and T2 from equations (15) and (17) with equation (23), we get(25)T1=brTGxbc4ωdrmrωdcmc1λr+λc,T2=−brTGxbc4ωdrmrωdcmc1λr∗+λc,which leads to the following expression for the matrix T:(26)T=brTGxbc4ωdrmrωdcmc1λr+λc−1λr∗+λc.If we want to obtain the real version in equation (21),(27)ReT=brTGxbc4ωdrmrωdcmcRe1λr+λc−1λr∗+λc.Equation (27) can also be obtained from equation (16). By combining equation (16) with equation (23) and assuming real mode shapes, it follows that(28)C=2πBbrTGxbc4ωdrmrωdcmc1λr+λc−1λr∗+λc−1λr+λc∗−1λr∗+λc∗BT.Now, combining equation (28) with equation (21) leads directly to equation (27). It is worth noticing that having an estimate of the matrix T might be helpful in estimating the modal masses if the input SD matrix Gx is known or it might be helpful estimating the input matrix if the modal masses are known.It is also worth checking the result against the classic solution for a 1-degree-of-freedom (1-DOF) system with the poleλ. In this case, the mode shape reduces to 1, the input SD matrix Gx becomes a scalar, and the correlation matrix C reduces to the variance σy2 of the response of the 1-DOF system. Rewriting equations (21) and (27) for 1-DOF system and combining the resulting equations, we get(29)σy2=πGxωd2m2Re12λ+12Reλ.After some reduction using that the classic definition of the poleλ=−ζω0+iωd, we obtain(30)σy2=πGx2ζm2ω03=πGxck,where c is the damping coefficient and k is the stiffness coefficient in the classic 1-DOF differential equation. This is the solution known from the literature (see example 1 in chapter 8 of [13]).At this point, it is worth highlighting that, similar to equation (20), equation (21) is only a general definition of the modal participation Γ. In reality, since the load spectral density, Gx, is not known in an output-only vibration test, the modal participation cannot be computed with this expression. Instead, a new expression needs to be formulated to estimate the modal participation by making use of the covariance matrix, C. This is elaborated in the following section.


## 2.1. General Form of the Response Correlation Matrix

It is obvious that since the correlation matrixC=EytyTt can now be calculated in two different ways using equations (2) and (6), C+=R0+ and C_=R0_, we will start showing that the two different ways of calculation leads to the same result.From the bottom part of equation (2), we have(9)C+=2π∑n=1MBnGxAn+Bn∗GxAn∗.Substituting equation (4) in equation (9) and carrying out the multiplications, we get(10)C+=2π∑n=1M∑s=1MAs−λn−λs+As∗−λn−λs∗GxAn+∑s=1MAs∗−λn∗−λs∗+As−λn∗−λsGxAn∗=2π∑n=1M∑s=1MAs−λn−λsGxAn+As−λn∗−λsGxAn∗+As∗−λn∗−λs∗GxAn∗+As∗−λn−λs∗GxAn.Swapping the summation indexes in equation (10) and contracting the equation back to short form yields(11)C+=2π∑s=1M∑n=1MAsGxAn−λs−λn+AsGxAn∗−λs−λn∗+As∗GxAn∗−λs∗−λn∗+As∗GxAn−λs∗−λn=2π∑s=1MAsGx∑n=1MAn−λs−λn+An∗−λs−λn∗+As∗Gx∑n=1MAn∗−λs∗−λn∗+As−λs∗−λn=2π∑s=1MAsGxBs+As∗GxBs∗,which actually is equal to what we get from the top part of equation (2) for R0−, meaning that the two formulations C+=R0+ and C−=R0− are identical.In order to get a little closer to the structure of this theoretical solution for the response correlation matrix, we continue manipulating equations (9) and (10). By substituting equation (3) in these equations and renaming the summation indexes (s⟶r and n⟶c), we obtain(12)C+=2π∑r,c=1MbrbrTGxbcar−λr−λcacbcT+br∗brHGxbcar∗−λr∗−λcacbcT+brbrTGxbc∗ar−λr−λc∗ac∗bcH+br∗brHGxbc∗ar∗−λr∗−λc∗ac∗bcH,or in a more compact form(13)C+=2π∑r,c=1Mbrt1rcbcT+br∗t2rcbcT+brt3rcbcH+br∗t4rcbcH,with tirc=tir,c∈ℂ∀i∈1,2,3,4 given by(14)t1rc=brTGxbcar−λr−λcac,t2rc=brHGxbcar∗−λr∗−λcac,t3rc=brTGxbc∗ar−λr−λc∗ac∗,t4rc=brHGxbc∗ar∗−λr∗−λc∗ac∗.As a side remark, it is worth noticing that we without any loss of generality can remove the generalized modal massesan from the equation using the mass scaled mode shapes bn/an. However, since we later will use mode shapes with arbitrary scaling, we will keep the explicit modal masses in the equations.Considering the first term in the summation of equation (12), we see that the scalar t1rc=brTGxbc/ar−λr−λcac can be seen as a matrix element with row index r and column index c and that we can write the summation as the outer product BT1BT where B is the mode shape matrix B=b1b2⋯bM∈ℂNo×M and(15)T1=t1rc=t1r,c=t11,1⋯t11,M⋮⋱⋮t1M,1⋯t1M,M=brTGxbcar−λr−λcac∈ℂM×Mis the matrix containing the elements t1rc∀r,c∈1,2,…,M. Using this notation on all four terms of equation (12), we get(16)C+=2πBT1BT+B∗T2BT+BT3BH+B∗T4BH,where(17)T2=brHGxbcar∗−λr∗−λcac∈ℂM×M,T3=brTGxbc∗ar−λr−λc∗ac∗∈ℂM×M,T4=brHGxbc∗ar∗−λr∗−λc∗ac∗∈ℂM×M.Evaluating the bottom part of equation (6) for τ=0, we now obtain a similar relation expressed in the modal participation vectors:(18)C+=2π∑n=1MγnbnT+γn∗bnH=2πΓBT+Γ∗BH,where Γ=γ1γ2⋯γM∈ℂNo×M is the modal participation vector matrix. Since the two first terms of equation (16) come from the first term of equation (9) and thus equal ΓBT, we can conclude that the modal participation matrix can be expressed by(19)Γ=BT1+B∗T2,and we see that the last term in equation (18) Γ∗BH equals the last two terms in equation (16).Taking origin in the top part of equation (2) instead of equation (6) and using the definition of the modal participation matrix given by equation (19) leads to satisfaction of the equation C−=2πBΓT+B∗ΓH, which is the version of equation (18) corresponding to the top part of equation (6).Finally, we can conclude that equation (19) is a general result for the modal participation matrix and that our final and general result for the response correlation matrix is(20)Γ=2πBT1+B∗T2BT+B∗T1∗+BT2∗BH.


## 2.2. Correlation Matrix for Proportional Damping

In case of real modes, equation (20) becomes(21)C=4πBReTBT,where Re⋅ denotes the real value of a complex quantity and(22)T=T1+T2.For proportional damping, the generalized modal mass becomes pure imaginary and is given by(23)an=2iωdnmn,where ωdn is the damped angular natural frequency and mn is the classic modal mass. Furthermore, it is known (as seen, for instance, in [7]) that from the orthogonality equations, we have in the case of general damping(24)λr+λcmrc+crc=0,r≠c,2λnmn+cn,r=c.In the case of proportional damping, the full modal mass and damping matricesmrc and crc in equation (24) reduce to the corresponding diagonal forms mn and cn. Combining T1 and T2 from equations (15) and (17) with equation (23), we get(25)T1=brTGxbc4ωdrmrωdcmc1λr+λc,T2=−brTGxbc4ωdrmrωdcmc1λr∗+λc,which leads to the following expression for the matrix T:(26)T=brTGxbc4ωdrmrωdcmc1λr+λc−1λr∗+λc.If we want to obtain the real version in equation (21),(27)ReT=brTGxbc4ωdrmrωdcmcRe1λr+λc−1λr∗+λc.Equation (27) can also be obtained from equation (16). By combining equation (16) with equation (23) and assuming real mode shapes, it follows that(28)C=2πBbrTGxbc4ωdrmrωdcmc1λr+λc−1λr∗+λc−1λr+λc∗−1λr∗+λc∗BT.Now, combining equation (28) with equation (21) leads directly to equation (27). It is worth noticing that having an estimate of the matrix T might be helpful in estimating the modal masses if the input SD matrix Gx is known or it might be helpful estimating the input matrix if the modal masses are known.It is also worth checking the result against the classic solution for a 1-degree-of-freedom (1-DOF) system with the poleλ. In this case, the mode shape reduces to 1, the input SD matrix Gx becomes a scalar, and the correlation matrix C reduces to the variance σy2 of the response of the 1-DOF system. Rewriting equations (21) and (27) for 1-DOF system and combining the resulting equations, we get(29)σy2=πGxωd2m2Re12λ+12Reλ.After some reduction using that the classic definition of the poleλ=−ζω0+iωd, we obtain(30)σy2=πGx2ζm2ω03=πGxck,where c is the damping coefficient and k is the stiffness coefficient in the classic 1-DOF differential equation. This is the solution known from the literature (see example 1 in chapter 8 of [13]).At this point, it is worth highlighting that, similar to equation (20), equation (21) is only a general definition of the modal participation Γ. In reality, since the load spectral density, Gx, is not known in an output-only vibration test, the modal participation cannot be computed with this expression. Instead, a new expression needs to be formulated to estimate the modal participation by making use of the covariance matrix, C. This is elaborated in the following section.


## 3. Estimating the Modal Participation from the Response CF Matrix (Proposed Approach)

In this section, we will consider several issues of practical interest using the expressions of the previous section. First let us simplify equation (21) by writing(31)C=4πBTrBT,where Tr=ReT. Similarly, we can now write equation (18) for real modes as(32)C=2πΓBT+Γ∗BH=4πΓrBT,where Γr=ReΓ, so that(33)Γr=BTr.In case we have estimated all modal parameters from the operating responses of a structure and now want to estimate the modal participation from the response correlation matrix, the approach is simple. For a full set of modes, where we know the square mode shape matrixB and the response correlation matrix C, the real version of the modal participation matrix is simply found by solving equation (32) for Γr, which gives(34)Γr=14πCB−T,where B−T is the inverse of the transpose mode shape matrix. However, when using experimental data, both the mode shape matrix and the response correlation matrix are estimated quantities and thus influenced by noise. This means that equations (31)–(33) can only be approximately satisfied, and we have to solve the equations in a least square sense.If the number of modes is smaller than the number of degrees of freedom (DOFs), which is often the case in experimental dynamics, then we have the preferred situation of an over determinate problem, and we can then use a least square solution. If we start with equation (31) and solve it for Tr, we obtain(35)T^r=14πB+CB+T,where T^r is the least square estimate Tr and B+ and B+T are the pseudo inverse of the mode shape matrix and its transpose, respectively. Substituting this solution into equation (31) defines the smoothed version of the experimentally obtained response correlation matrix:(36)C^=4πBT^rBT.Furthermore, using equation (33) defines the estimated modal participation matrix as(37)Γ^r=BT^r.We could also start with equation (32), and we now get the following estimate:(38)Γ^r=14πCB+T,which is the least square equivalent to the solution given by equation (34). The two estimates given by equations (37) and (38) are only the same if we use the smoothed version of the response correlation matrix in equation (38). In this case, substituting equation (36) into equation (38), we obtain(39)Γ^r=BT^rBTB+T=BT^r.If there is some noise of significance on the response correlation matrix estimate, it should be considered to use the smoothed estimate.In the cases where the number of modes (M) is equal or greater than measured DOFs (No), we might need to look at ways of reducing the number of modes to secure an overdetermined problem so that we can use the abovementioned least square solutions. We can do that using a (frequency) band-by-band approach.Another reason for using such approach might be that we like to estimate the input SD matrixGx or its modal load counterpart:(40)BTGxB=brTGxbc.Also, we like to be able to obtain different estimates ofGx or BTGxB in different frequency bands. We solve this problem by generalizing equation (5) to the case where the modal participation vectors are frequency dependent γn=γnω:(41)γnω=BnGxbnan.We can now think about modelling a general response SD matrix similar to that given in equation (7) using the frequency dependent modal participation vectors instead in the constant ones. The band-by-band approach now materializes by using the Parseval equation given by equation (8) but now formulated for a certain frequency band(42)ΔB=ω;ω+Δω.Also, we have(43)ΔC=2Re∫ΔBGyωdω.Here, we have used that we need to integrate over the frequency band and its negative frequency counterpart, which provides the complex conjugate of the positive frequency integration, so by adding the two terms, we get the real result as given by equation (43).Thus, if we have a frequency band where a number of modes are dominating and we integrate over this frequency band as given by equation (43), then ΔC is approximately the contribution to the response correlation matrix from these modes. ΔC is exactly the contribution to the response correlation matrix from the frequency band ΔB but only approximately the contribution from the considered modes because the response from these modes covers the whole frequency band, but we will assume that the response from the modes is negligible outside of ΔB.Finally, let us discuss how noise can be reduced using the Parseval equation integrating the response SD matrix as given by equations (8) and (43). It is well known (as seen, for instance, in [14]) that the physics in the SD matrix can be decoupled from noise using a singular value decomposition like(44)Gyω=UωSωUHω,where the singular values in the diagonal matrixSω contain information about the modal coordinates in the frequency domain, and the singular vectors of the matrix Uω contain information about the mode shapes. Furthermore, it is known that the large singular values (and the corresponding singular vectors) represent the physics and lower singular values represent the noise.So, if we are consideringNm number of modes (in the abovementioned frequency band), then we can reduce the noise by keeping only the first n singular values of Sω and setting the rest of the singular values to zero. In this case, n≥Nm because if n<Nm, then we will exclude information from one of the mode shapes when estimating ΔC using equation (43).In practice, we will be using a sampled version of the SD matrixGyn=GynΔω. If we are using a sampled version of the theoretical SD matrix, then we obtain(45)ΔC=2Re∑nΔω∈ΔBGynΔω=2Δω∑nΔω∈ΔBReGyn,where Δω is the frequency resolution measured in rad/s, given by Δω=2πΔf, where Δf is the traditional frequency resolution measured in Hertz.


## 4. Simulation Case

We are considering a case with 5 modes and modal parameters as shown in Table1; all modal masses are assumed equal to unity. We are assuming a time step Δt=0.103s corresponding to the frequency resolution Δf=0.0095Hz or Δω=0.0596rad/s for the assumed time domain segment size equal to Ns=1024 data points. This corresponds to the Nyquist frequency fν=4.8544Hz. In the following, the approach introduced in Section 3 is applied to estimate the modal participation both from the exact and empirical SD matrices simulated with the modal parameters shown in Table 1.Table 1
Modal participation vector estimates for the first mode.
Parameter
Mode 1
Mode 2
Mode 3
Mode 4
Mode 5
Natural frequency (Hz)
1.0
1.2
3.0
3.2
4.0
Damping ratio (%)
2.0
1.0
1.5
1.0
0.5
Mode shape (mass scaled)
−0.3073
−0.4975
−0.4028
0.2276
−17.240
−11.535
0.6350
−0.2009
−0.9302
−10.723
0.5554
−0.0830
−0.580
−0.6404
0.2468
0.2995
−0.4454
−0.766
−0.1306
0.2896
−0.0115
0.8423
17.276
0.5856
0.3342
−0.0200
−0.1401
11.322
−0.9856
−0.0964


### 4.1. Estimation Using the Exact SD Matrix

The exact complex modal participation vectors are then calculated according to equation (5), and the exact real modal participation vectors are calculated according to equations (27) and (33). A full band estimate of C is found according to equation (45), i.e, by using numerical integration of the response SD matrix over the Nyquist band and subsequently by estimating the modal participation for all five modes using equations (35) and (37). Similarly, a band-by-band estimate is found by using partial integration according to equation (45) over the three bands(i)
ΔB1, containing the 2 first modes: f∈0:2.08 Hz(ii)
ΔB2, containing modes 3 and 4: f∈2.08:3.59 Hz(iii)
ΔB3, containing mode 5: f∈3.59:4.85 Hzto obtain the contributions to the correlation matrix from each of the bands and then by using equations (35) and (37) for each of the correlation matrix contributions. When the modal participation vectors have been estimated, the resulting synthesized SD matrix is obtained from equation (7).The resulting SD plots are shown in Figure1, where the top and middle plots are based on the full band estimate of C and the bottom plot is based on the band-by-band estimates. In the top plot, the singular values of the SD matrix synthesized with the estimated modal participation is compared to those of the exact SD matrix, whereas in the two bottom plots, this comparison is illustrated in terms of the cross spectral density g12f=Gy1,2f.Figure 1
Exact density functions of the simulation case (red) and estimated density functions (blue). Top plot shows the singular values of the spectral density matrix,Sω, whereas the two bottom plots show the cross spectral density function, g12f=Gy1,2f, between the first and second DOFs. The middle plot shows the full band estimate, and the bottom plot shows the band-by-band estimate.In general, we can see that the real modal participation vector estimates of the spectral densities are quite good. We see from the single density plots (2 bottom plots of the Figure1) that the estimates mainly influence the zeros of the SD function, but since it is known that the zeros do not have much physical meaning for SD functions, this should not be of any concern. This influence is mainly due to the fact that the zero estimates are more sensitive to the loss of mode shape phase and, in turn, to the assumption of real mode shapes. Finally, we see that the difference between the full band and band-by-band estimates is quite small—only the zeros are moving a bit.Table2 shows the values of the modal participation for the first mode. We see that the real modal participation vector is just the real value of the vector, and the full band and band-by-band estimates show only small deviation from the exact real vector.Table 2
Modal participation vector estimates for the first mode.
Exact complex
Exact real
Full band estimate
First band estimate
−0.3670 – 0.0557i
−0.3670
−0.3668
−0.3728
−1.4270 + 0.1262i
−14.270
−14.266
−14.283
0.6827 – 0.0264i
0.6827
0.6826
0.6805
0.3762 – 0.0754i
0.3762
0.3760
0.3736
−0.0309 + 0.1238i
−0.0309
−0.0308
−0.0234
−0.0221 + 0.0005i
−0.0221
−0.0219
−0.0232


### 4.2. Estimation Using the SD Matrix Estimated from the Vibration Responses

In order to illustrate the application of the approach introduced in Section3 from an experimental perspective, the modal participation was also computed using the empirical SD matrix estimated from the simulated responses.These responses were obtained using Gaussian white noise uncorrelated univariate inputs in all DOFs with 32 k data points. The empirical SD matrix is found using Welch averaging with 50% overlap with 1024 data in each segment. In this case, the correlation matrix is estimated as a full band estimate from the SD matrix estimate just like in the previous application example. The results are shown in Figure2. In this case, where the average of the empirical (i.e., estimated) spectral density is equal to the correlation matrix, the theoretical SD matrix obtained from equation (7) has to be divided by the total band width measured in rad/s NsΔω=2πfs in order to provide a fit to the empirical spectral densities.Figure 2
Estimated empirical density functions of the simulation case (red) and synthesized density functions (blue). Top plot shows the singular values,Sω, of the spectral density matrix, whereas the bottom plot shows the cross spectral density function, g12f=Gy1,2f, between the first and second DOFs.The obtained result is shown in Figure2. Also in this figure, the singular values of the SD matrix estimated from the simulated responses are compared to those from SD matrix synthesized with the estimated modal participation in the top plot, whereas the bottom plot shows this comparison in terms of the cross spectral density g12f=Gy1,2f. Again, we see that the fit is quite good and only the zeros of the single SD function show some deviation.Many different simulation cases have been carried out, and the results shown here are typical. Cases with more closely spaced modes do not show results that are different in nature.


## 4.1. Estimation Using the Exact SD Matrix

The exact complex modal participation vectors are then calculated according to equation (5), and the exact real modal participation vectors are calculated according to equations (27) and (33). A full band estimate of C is found according to equation (45), i.e, by using numerical integration of the response SD matrix over the Nyquist band and subsequently by estimating the modal participation for all five modes using equations (35) and (37). Similarly, a band-by-band estimate is found by using partial integration according to equation (45) over the three bands(i)
ΔB1, containing the 2 first modes: f∈0:2.08 Hz(ii)
ΔB2, containing modes 3 and 4: f∈2.08:3.59 Hz(iii)
ΔB3, containing mode 5: f∈3.59:4.85 Hzto obtain the contributions to the correlation matrix from each of the bands and then by using equations (35) and (37) for each of the correlation matrix contributions. When the modal participation vectors have been estimated, the resulting synthesized SD matrix is obtained from equation (7).The resulting SD plots are shown in Figure1, where the top and middle plots are based on the full band estimate of C and the bottom plot is based on the band-by-band estimates. In the top plot, the singular values of the SD matrix synthesized with the estimated modal participation is compared to those of the exact SD matrix, whereas in the two bottom plots, this comparison is illustrated in terms of the cross spectral density g12f=Gy1,2f.Figure 1
Exact density functions of the simulation case (red) and estimated density functions (blue). Top plot shows the singular values of the spectral density matrix,Sω, whereas the two bottom plots show the cross spectral density function, g12f=Gy1,2f, between the first and second DOFs. The middle plot shows the full band estimate, and the bottom plot shows the band-by-band estimate.In general, we can see that the real modal participation vector estimates of the spectral densities are quite good. We see from the single density plots (2 bottom plots of the Figure1) that the estimates mainly influence the zeros of the SD function, but since it is known that the zeros do not have much physical meaning for SD functions, this should not be of any concern. This influence is mainly due to the fact that the zero estimates are more sensitive to the loss of mode shape phase and, in turn, to the assumption of real mode shapes. Finally, we see that the difference between the full band and band-by-band estimates is quite small—only the zeros are moving a bit.Table2 shows the values of the modal participation for the first mode. We see that the real modal participation vector is just the real value of the vector, and the full band and band-by-band estimates show only small deviation from the exact real vector.Table 2
Modal participation vector estimates for the first mode.
Exact complex
Exact real
Full band estimate
First band estimate
−0.3670 – 0.0557i
−0.3670
−0.3668
−0.3728
−1.4270 + 0.1262i
−14.270
−14.266
−14.283
0.6827 – 0.0264i
0.6827
0.6826
0.6805
0.3762 – 0.0754i
0.3762
0.3760
0.3736
−0.0309 + 0.1238i
−0.0309
−0.0308
−0.0234
−0.0221 + 0.0005i
−0.0221
−0.0219
−0.0232


## 4.2. Estimation Using the SD Matrix Estimated from the Vibration Responses

In order to illustrate the application of the approach introduced in Section3 from an experimental perspective, the modal participation was also computed using the empirical SD matrix estimated from the simulated responses.These responses were obtained using Gaussian white noise uncorrelated univariate inputs in all DOFs with 32 k data points. The empirical SD matrix is found using Welch averaging with 50% overlap with 1024 data in each segment. In this case, the correlation matrix is estimated as a full band estimate from the SD matrix estimate just like in the previous application example. The results are shown in Figure2. In this case, where the average of the empirical (i.e., estimated) spectral density is equal to the correlation matrix, the theoretical SD matrix obtained from equation (7) has to be divided by the total band width measured in rad/s NsΔω=2πfs in order to provide a fit to the empirical spectral densities.Figure 2
Estimated empirical density functions of the simulation case (red) and synthesized density functions (blue). Top plot shows the singular values,Sω, of the spectral density matrix, whereas the bottom plot shows the cross spectral density function, g12f=Gy1,2f, between the first and second DOFs.The obtained result is shown in Figure2. Also in this figure, the singular values of the SD matrix estimated from the simulated responses are compared to those from SD matrix synthesized with the estimated modal participation in the top plot, whereas the bottom plot shows this comparison in terms of the cross spectral density g12f=Gy1,2f. Again, we see that the fit is quite good and only the zeros of the single SD function show some deviation.Many different simulation cases have been carried out, and the results shown here are typical. Cases with more closely spaced modes do not show results that are different in nature.


## 5. Application on the Heritage Court Tower Data

In order to illustrate the approach proposed in Section3 to estimate the modal participation factors, the vibration responses of the Heritage Court Tower (HCT) depicted in Figure 3 are used as primary data. The HCT vibration responses are a well-known set of data in the literature of operational modal analysis.Figure 3
North (a) and east (b) building faces of the Heritage Court Tower in Vancouver, Canada.
(a)
(b)The original publication of the ambient vibration test of the HCT is found in [15], and a full operational modal analysis was presented at the same year of this publication by Brincker and Andersen [16]. The case is also one of the examples in the text book by Brincker and Ventura [4], and the data are a part of the associated Matlab toolbox [17].The HCT data consist of four datasets, and the first datasets to be considered here have six measurement channels (remaining datasets have eight measurements channels). We will use the modal parameters estimated by example 9.1 in [4], where identification is carried out using the time-domain polyreference technique [18, 19].In this section, we will just synthesize the first three modes of dataset 1 using equation (45), where we have estimated the partial response correlation matrix by numerical integration from 0.5 Hz up to 3.0 Hz. The result is shown in Figure 4. The sampling time step was 0.025 s corresponding to a Nyquist frequency of 50 Hz. Spectral density was estimated by Welch averaging using 2048 data points in the data segments. This gives only a few averages and thus a relatively high estimation error.Figure 4
Estimated empirical density functions of the Heritage Court Tower case (red) and synthesized density functions (blue). Top plot shows the singular values of the spectral density matrix,Sω, and the bottom plot shows the cross spectral density function, g12f=Gy1,2f, between the first and second DOFs.Again, we see that the fit is quite good, both for the singular value plot and for single SD function plot. Even though the experimental data might look quite noisy, the obvious noise content is not due to measurement noise, but nearly solely due to the abovementioned estimation error. Therefore, in this case, noise cannot be reduced by cancelling of lowest singular values as described in equation (44), and thus, the raw SD matrix has been used as the basis of the this analysis.


## 6. Conclusions

We have shown that two matrices can be defined, a symmetric matrixT1 and a Hermitian matrix T2, so that in the case of general damping, the modal participation matrix Γ can be obtained as Γ=BT1+B∗T2, where B is the mode shape matrix of the considered system.A general expression for the response correlation matrixC has been formulated based on the abovementioned matrices T1 and T2. In case of real mode shapes, this expressions turn into the simple relation C=4πBTrBT, where Tr is the real part of the sum of the two matrices T1 and T2. This also naturally leads to a corresponding real version of the participation matrix given by Γr=BTr, which has the right physical interpretation since it relates the response correlation matrix with the mode shape matrix in a simple way as C=4πΓrBT.Finally, it has been shown how the real version of modal participation matrix can be estimated from the response correlation matrix and how the so-estimated modal participation vectors can be used to synthesize empirical spectral density functions. This has been illustrated both on simulated as well as real data.



---

*Source: 9347075-2019-08-28.xml*