Understanding Power Spectrum Density Transmissibility
Power spectrum density transmissibility (PSDT) is a type of complex frequency domain function proposed recently. It describes the relation between cross-spectra of system outputs. Since PSDTs with same local-reference degree of freedom (DOF) combination but with different transferring output DOFs cross each other at the system’s poles under certain load condition, the functions have been used as the primary data in operational modal analysis (OMA) to extract modal parameters, and such technique is named as PSDT-based OMA (PSDTOMA). Because PSDT is a concept that appears recently, researches on which, especially, in-depth discussions aim at the essence and properties of PSDT are very rare. For appropriate application of PSDTOMA, it is necessary to perform further study on such problem obviously. In this paper, two paths to get PSDT, which, respectively, are referred to as PSDT estimator and PSDT syntheticism, are given firstly; some properties about PSDT are explored based on the PSDT syntheticism, and the relation between PSDT and single reference transmissibility function (STF) is analyzed. Finally, the above conclusions are verified with numerical values and experimental data.
Similar to more common-used transfer functions, transmissibility is a type of (complex) frequency domain functions, which describe the relation between physical quantities of the same type at different locations in a dynamical system. According to the type of physical quantities described, the transmissibility about system dynamic characteristics can be divided into force transmissibility and motion transmissibility.
Power spectrum density transmissibility (PSDT), which this paper studies, was first proposed by Yan and Ren  in 2012 as the third kind of motion transmissibility following the traditional single reference transmissibility function (STF) [2–4] and multi-/poly-reference transmissibility function (MTF/PTF) [5–7], and it is defined as the ratio of cross-spectra between responses at different locations in a system. Similar to STF, PSDT is a scalar (complex) frequency domain function as well, representing the relation between responses of local degree of freedom (DOF) and reference DOF. Because PSDT introduces the concept of transferring output DOF additionally and the PSDT, depending on different transferring output DOFs, intersects at the system poles under certain loading condition [1, 8], in theory, the STF, depending on different loading conditions, can be replaced with the PSDT, depending on different transferring output DOFs, which is used as the primary data of operational modal analysis (OMA) to lower identification method’s demand on the number of loading conditions. This OMA process is called PSDT-based OMA (PSDTOMA). Same as STF-based OMA (sTOMA), PSDTOMA has the characteristic that identification results are unaffected by harmonic component in excitation, compared to the traditional OMA methods.
Because PSDT was proposed recently, it has not been studied thoroughly and deeply. To the authors’ knowledge, all the current PSDT-related studies focus on the application layer [1, 8], and no discussion on its essence and property has been made. In this paper, two paths to get PSDT, which, respectively, are referred to as PSDT estimator and PSDT syntheticism, are given firstly; some properties about PSDT are explored based on the PSDT syntheticism, and the relation between PSDT and single reference transmissibility function (STF) which is similar to PSDT is analyzed. Finally, the above conclusions are verified with numerical values and experimental data.
2. Power Spectrum Density Transmissibility
Under a certain load condition, PSDT is defined as the ratio of cross-spectral densities between response at DOF and responses at DOFs , , respectively [1, 8]. It is expressed aswhere and are cross-spectral densities between responses at DOFs , and responses at DOFs , , respectively; DOFs , , and are referred to as local DOF, reference DOF, and transferring output DOF, respectively.
In the frequency domain, the relation between inputs and outputs of a linear time invariant (LTI) system can be expressed aswhere is output power spectrum matrix; is input power spectrum matrix; is transfer function matrix; denote the complex conjugate and transpose of . According to (2), the cross-spectral density between response at two DOFs, such as and , which also is the entry in the th row and th column of , can be expressed aswhere and are row vectors consisting of entries in the th and th row of , respectively. Combining (1) and (3), we get another expression of PSDT:Equations (1) and (4) are two opposite paths to get the value of PSDT. Equation (1) is referred to as PSDT estimator, which is used for getting PSDT estimation from outputs, while (4) is referred to as PSDT syntheticism, which is used for getting PSDT synthesis from inputs and transfer function. The relation among PSDT syntheticism, PSDT estimator, PSDT synthesis, and PSDT estimation is shown in Figure 1. In theory, PSDT estimation is a consistent approach of PSDT synthesis, and this makes PSDT synthesis be regarded as an evaluation criterion for PSDT estimation.
Actually, the result obtained by putting random excitations into (4) directly is an approximation of PSDT. It is because there exists interference caused by spectral leakage. In some cases, however, (4) can be further simplified to get the exact value of PSDT. In theory , the energy of the stationary random process, , is distributed uniformly in the frequency domain, so the amplitude of its autospectrum is unrelated to (circle) frequent and proportional to variance of the process , and it can be expressed aswhere is a real nonzero proportion factor, related to the sampling length. Therefore, if two stationary random processes, and , are fully correlative, the amplitude of their cross-spectrum is proportional to covariance of these two processes , and it can be expressed asNevertheless, if and are mutually uncorrelated zero-mean random processes, the amplitude of cross-spectrum will approach zero after averaging many times; that is, , as shown in Figure 2.
According to (5) and (6), the amplitude of matrix can be uniformly expressed asIn the case that the system’s excitations are correlated completely, the expression of matrix is as follows:where are the standard deviations of random excitations. In the case that all the excitations are fully independent, the value of is the diagonal matrix as follows:Finally, if excitations are partially correlative, the value of can be obtained by appropriately choosing off-diagonal entries. Based on the above discussion, (4) can be rewritten asThus, PSDT’s amplitude can be obtained exactly by using (10).
Let be the following vector: in (8) can be expressed asBy putting (12) into (10), we getObviously, in the case that the system’s excitations are fully independent, the value of PSDT’s amplitude is unrelated to the way of selecting transferring output DOF . More particularly, in the case that the system is subject to single concentrated excitation, for instance, the force is acting on DOF ; (10) degrades into the following form:which means that the amplitude of PSDT depends on the location, rather than the strength of excitation.
3. Power Spectrum Density Transmissibility and Single Reference Transmissibility
STF is the other kind of scalar motion transmissibility for MDOF system. In frequency domain, STF is defined as the ratio of Fourier transforms of responses at local DOF and DOF [2–4]; that is,Besides (2), the relation between inputs and outputs of LTI system can also be expressed aswhere is input spectral vector and is output spectral vector. According to (16), the response at some DOF, such as , can be expressed asBy combining (15) and (17), we can getSimilar to (4), (18) can be used to synthesize STF from input information and transfer function. It is important to note that because there exists infinite error introduced by end effects , nothing but useless result would be obtained by putting random excitations directly. Hence, (18) can only be used to get STF synthesis in deterministic case.
Consider that the autospectrum of random process , , can be expressed aswhere is the length of time history; is the Fourier transform of ; denotes the complex conjugate of . By combining (5) and (19), we can get the following equation:Therefore, the amplitude of input vector can be expressed aswhere are the standard deviations of random excitations, respectively. Putting (21) into (18), we getEquation (22) can be used to synthesize the amplitude of STF under random excitation, since it is unaffected by end effects and spectral leakage. Again, in the case that the system is subject to single concentrated excitation, (22) degrades into the following form:By comparing (22) with (13), as well as (23) with (14), we getTherefore, in the cases that the system is subject to single concentrated excitation and all the inputs of the system are correlated completely, PSDT is exactly identical with the STF with the same local-reference DOF combination.
4. Numerical Study
Next, the above conclusions were verified by using the simply supported slab shown in Figure 3. The slab’s finite element model (FEM) was established by using ANSYS software. The model, whose size is 3.6 m × 2 m × 0.1 m, was uniformly divided into 720 elements. Its modulus of elasticity is 2 × 1011 N/m2; material density is 7800 kg/m3; Rayleigh damping coefficients α and β are 3.0 s−1 and 1 × 10−4 s, respectively. 21 nodes on the upper surface of the model were selected at equal spacing as the control points used to apply excitations and obtain responses. The positions and numbers of these control points are shown in Figure 3.
13 kinds of random excitation loading conditions, listed in Table 1, were applied to the model successively and the responses at all the 21 control points were extracted. After that, amplitudes of PSDT measurements and syntheses were obtained by using (1) and (10), respectively. For comparison, amplitudes of STFs were synthesized by using (13).
By using the loading conditions in case 1 and case 2, the PSDT syntheses, depending on different transferring output DOFs, are shown in Figure 4. It can be seen that the PSDT syntheses coincide exactly with the corresponding STF synthesis under the two loading conditions, no matter what the transferring output DOF is. Figure 5 shows STF and PSDT syntheses in case 13, from which it can be seen that STF and PSDTs are no longer identical, and the value of PSDT changes with transferring output DOF’s changing.
The PSDT measurements, respectively, in single input condition (case 1) and single source condition (case 2), with corresponding STF syntheses under the two loading conditions, are shown in Figure 6. With reference to Figure 4, it can be seen that PSDT measurements coincide exactly with PSDT syntheses under the two loading conditions. The PSDT measurements and PSDT syntheses in case 13 are shown in Figure 7. Despite the uncertain components caused by fewer averaging times (about 100 times), spectral leakage in PSDT measurements, and other factors, it still can be seen clearly that PSDT syntheses and PSDT measurements are identical. Hence the correctness of PSDT syntheticism ((10), (13) and (14)) is proved.
Comparing with loading condition in case 2, the input applied to point 21 was replaced with a random process uncorrelated with all the other inputs in case 3. From Figure 8 it can be seen that the PSDT measurements, depending on , were no longer identical with STF synthesis.
Furthermore, the function is defined to measure the similarity between STF syntheses and PSDT measurements, and it is expressed aswhere is within 0-1; the bigger the value is, the closer the PSDT measurements and STF syntheses are, with “1” indicating the two are identical. In case 2–case 13, the relation between , depending on different transferring output DOFs, and loading conditions is shown in Figure 9. In the case that excitations are fully correlated, the values of between all the PSDT measurements and STF syntheses were close to 1, while, under other loading conditions, no obvious correlation between values and the number of independent excitations was found.
5. Experimental Study
Single source excitation is a commonly used mode of excitation, so the above conclusions can be verified easily in practical applications. Figure 10 shows an experimental setup of submarine pipeline vibration-table model. 10 acceleration sensors were installed on the model from south to north, respectively, at the measure points 1–10, to measure the model’s responses in the east-west direction, as shown in Figure 11. When the submarine pipeline is not immersed in water (left graph in Figure 10), the input of the model is totally provided by the vibration-table surface, and thus it can be regarded as a single source excitation system. When the submarine pipeline is immersed in water (right graph in Figure 10), the model is subject to the vibration table, water flow, and waves simultaneously, and thus it can be seen as a random multiple excitation system.
This study takes the loading condition that the model, which was not immersed in water, was only subject to the 0–80 Hz random east-west excitation provided by the vibration table as case 1 and takes the loading condition that the model, which was immersed in water, was subject to both the excitation in case 1 and the east-west action of water flow and waves as case 2. Figures 12 and 13 show the PSDT measurements obtained based on model responses and, respectively, taking the measure points 2, 6, 5, and 9 as local DOFs and reference DOFs. From them it can be seen that, in case 1, the PSDT obtained by taking the response of any one of the measure points as transferring output was basically identical to each other (Figures 12(a) and 13(a)). Though the excitation provided by water flow and waves was much weaker than that provided by the vibration table in the experiment, significant difference still could be seen among the PSDTs depending on different transferring output DOFs (Figures 12(b) and 13(b)).
This paper presented the syntheticism of PSDT (4) and STF (18), which can be used to calculate the theoretical value of transmissibility with known transfer functions and excitations. In particular, because the theoretical power spectrum’s amplitude of a random process is unrelated to frequency domain variables , the amplitude of PSDT synthesis can be obtained by using (10) after appropriately choosing off-diagonal entries and establishing matrix according to the correlation of excitation. In addition, the amplitude of excitation can be replaced with its standard deviation to calculate STF syntheses by using (18) in the case of random excitation and to obtain better results. By comparison, it can be seen that the results obtained by the above calculation are consistent with system response-based PSDT measurements obtained by using (1).
Based on the above syntheticism, this paper further explored some property of PSDT and reached the following conclusions:(1)PSDT and STF are two physical quantities that are closely related to each other and different from each other. Under the loading condition of completely correlated excitations (including single excitation), the amplitude of PSDT, depending on same local-reference DOF combination, coincides exactly with the amplitude of STF.(2)If inputs are fully correlative to each other (including the situation of single input), PSDT’s amplitude does not change with transferring output DOF ; or else, PSDT’s amplitude changes with different modes of selection and is generally different from that of STF.(3)In the case that excitations are not fully correlated, it seems that there is no obvious correlation between the similarity of PSDTs with the same local-reference DOF combination and the number of independent excitations.The correctness of the above conclusions has been proven by the numerical experiment and underwater vibration-table model experiment.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This study is supported by the Ministry of Science and Technology of China (Grants nos. 2011CB013702 and 2013CB035906) and National Natural Science Foundation of China (Grant no. 51378088).
W.-J. Yan and W.-X. Ren, “Operational modal parameter identification from power spectrum density transmissibility,” Computer-Aided Civil and Infrastructure Engineering, vol. 27, no. 3, pp. 202–217, 2012.View at: Publisher Site | Google Scholar
C. Devriendt, P. Guillaume, and G. De Sitter, “Identification of modal parameters from transmissibility measurements with multiple loads,” in Proceedings of the 12th International Congress on Sound and Vibration 2005 (ICSV '05), pp. 2919–2927, IEEE, Lisbon, Portugal, July 2005.View at: Google Scholar
C. Devriendt and P. Guillaume, “The use of transmissibility measurements in output-only modal analysis,” Mechanical Systems and Signal Processing, vol. 21, no. 7, pp. 2689–2696, 2007.View at: Publisher Site | Google Scholar
C. Devriendt and P. Guillaume, “Identification of modal parameters from transmissibility measurements,” Journal of Sound and Vibration, vol. 314, no. 1-2, pp. 343–356, 2008.View at: Publisher Site | Google Scholar
A. M. R. Ribeiro, J. M. M. Silva, and N. M. M. Maia, “On the generalisation of the transmissibility concept,” Mechanical Systems and Signal Processing, vol. 14, no. 1, pp. 29–35, 2000.View at: Publisher Site | Google Scholar
M. Fontul, A. M. R. Ribeiro, J. M. M. Silva, and N. M. M. Maia, “Transmissibility matrix in harmonic and random processes,” Shock and Vibration, vol. 11, no. 5-6, pp. 563–571, 2004.View at: Publisher Site | Google Scholar
N. M. M. Maia, R. A. B. Almeida, and A. P. V. Urgueira, “Understanding transmissibility properties,” in Proceedings of the Conference & Exposition on Structural Dynamics—Technologies for Civil Structures (IMAC XXVI '08), Orlando, Fla, USA, February 2008.View at: Google Scholar
I. G. Araújo and J. E. Laier, “Operational modal analysis using SVD of power spectral density transmissibility matrices,” Mechanical Systems and Signal Processing, vol. 46, no. 1, pp. 129–145, 2014.View at: Publisher Site | Google Scholar
J. S. Bendat and A. G. Piersol, Random Data: Analysis and Measurement Procedures, John Wiley & Sons, 3rd edition, 2000.View at: MathSciNet
P. M. T. Broersen, “A comparison of transfer function estimators,” IEEE Transactions on Instrumentation and Measurement, vol. 44, no. 3, pp. 657–661, 1995.View at: Publisher Site | Google Scholar