Response variance prediction in the statistical energy analysis of built-up systems

In the statistical energy analysis (SEA) of high frequency noise and vibration, a complex engineering structure is represented as an assembly of subsystems. The response of the system to external excitation is expressed in terms of the vibrational energy of each subsystem, and these energies are found by employing the principle of power balance. Strictly the computed energy is an average taken over an ensemble of random structures, and for many years there has been interest in extending the SEA prediction to the variance of the energy. A variance prediction method for a general built-up structure is presented here. Closed form expressions for the variance are obtained in terms of the standard SEA parameters and an additional set of parameters alpha(k) that describe the nature of the power input to each subsystem k, and alpha(ks) that describe the nature of the coupling between subsystems k and s. The theory is validated by comparison with Monte Carlo simulations of plate networks and structural-acoustic systems.     

Numerical and experimental validation of a hybrid finite element-statistical energy analysis method

The finite element (FE) and statistical energy analysis (SEA) methods have, respectively, high and low frequency limitations and there is therefore a broad class of "mid-frequency" vibro-acoustic problems that are not suited to either FE or SEA. A hybrid method combining FE and SEA was recently presented for predicting the steady-state response of vibro-acoustic systems with uncertain properties. The subsystems with long wavelength behavior are modeled deterministically with FE, while the subsystems with short wavelength behavior are modeled statistically with SEA. The method yields the ensemble average response of the system where the uncertainty is confined in the SEA subsystems. This paper briefly summarizes the theory behind the method and presents a number of detailed numerical and experimental validation examples for structure-borne noise transmission.     

Vibro-acoustic analysis of complex systems

A general method is presented for predicting the ensemble average steady-state response of complex vibro-acoustic systems that contain subsystems with uncertain, or random, properties. The method combines deterministic and statistical techniques to produce a non-iterative hybrid method that incorporates equations of dynamic equilibrium and power balance. The method is derived explicitly without reference to statistical energy analysis (SEA); however, it is seen that the wave approach to SEA can be viewed as a special case of the proposed method. The proposed method provides a flexible way to account for necessary deterministic details in a vibro-acoustic analysis without requiring that an entire system be modeled deterministically. The method therefore provides a potential solution to the mid-frequency problem (in which a system is neither entirely deterministic nor entirely statistical). The application of the method is illustrated with a numerical validation example.     

Extension of the statistical modal energy distribution analysis for estimating energy density in coupled subsystems

The present article deals with an extension of the Statistical modal Energy distribution Analysis (SmEdA) method to estimate kinetic and potential energy density in coupled subsystems. The SmEdA method uses the modal bases of uncoupled subsystems and focuses on the modal energies rather than the global energies of subsystems such as SEA (Statistical Energy Analysis). This method permits extending SEA to subsystems with low modal overlap or to localized excitations as it does not assume the existence of modal energy equipartition. We demonstrate that by using the modal energies of subsystems computed by SmEdA, it is possible to estimate energy distribution in subsystems. This approach has the same advantages of standard SEA, as it uses very short calculations to analyze damping effects. The estimation of energy distribution from SmEdA is applied to an academic case and an industrial example.     

STATISTICAL ENERGY ANALYSIS OF COUPLED PLATE SYSTEMS WITH LOW MODAL DENSITY AND LOW MODAL OVERLAP

Finite element methods, experimental statistical energy analysis (ESEA) and Monte Carlo methods have been used to determine coupling loss factors for use in statistical energy analysis (SEA). The aim was to use the concept of an ESEA ensemble to facilitate the use of SEA with plate subsystems that have low modal density and low modal overlap. An advantage of the ESEA ensemble approach was that when the matrix inversion failed for a single deterministic analysis, the majority of ensemble members did not encounter problems. Failure of the matrix inversion for a single deterministic analysis may incorrectly lead to the conclusion that SEA is not appropriate. However, when the majority of the ESEA ensemble members have positive coupling loss factors, this provides sufficient motivation to attempt an SEA model. The ensembles were created using the normal distribution to introduce variation into the plate dimensions. For plate systems with low modal density and low modal overlap, it was found that the resulting probability distribution function for the linear coupling loss factor could be considered as lognormal. This allowed statistical confidence limits to be determined for the coupling loss factor. The SEA permutation method was then used to calculate the expected range of the response using these confidence limits in the SEA matrix solution. For plate systems with low modal density and low modal overlap, relatively small variation/uncertainty in the physical properties caused large differences in the coupling parameters. For this reason, a single deterministic analysis is of minimal use. Therefore, the ability to determine both the ensemble average and the expected range with SEA is crucial in allowing a robust assessment of vibration transmission between plate systems with low modal density and low modal overlap.     

基于高级统计能量分析的周期加筋板振动特性研究

统计能量分析(statistical energy analysis, SEA)是复杂耦合系统中、高频动力学特性计算的有力工具. 本文以波传播理论和SEA的基本原理为基础, 研究周期加筋板中弯曲波传播特性. 分析了周期结构的频率带隙特性和加强筋对板上弯曲波的滤波特性对SEA计算结果的影响规律, 发现经典SEA由于忽视了加筋板中物理上不相邻子系统间存在的能量隧穿效应, 而导致响应预测结果产生最高近 40 dB的误差. 为了解决这一问题, 本文应用高级统计能量分析(advanced statistical energy analysis, ASEA)方法, 考虑能量在不相邻子系统间的传递、转移和转化的物理过程, 从而大幅提高子系统响应的预测精度, 将误差在大部分频段降低至小于5 dB. 设计了模拟简支边界条件的加筋板振动测试实验装置, 实验测试结果与有限元结果符合较好, 对理论模型进行了验证.     

A hybrid approach for predicting the distribution of vibro-acoustic energy in complex built-up structures

Finding the distribution of vibro-acoustic energy in complex built-up structures in the mid-to-high frequency regime is a difficult task. In particular, structures with large variation of local wavelengths and/or characteristic scales pose a challenge referred to as the mid-frequency problem. Standard numerical methods such as the finite element method (FEM) scale with the local wavelength and quickly become too large even for modern computer architectures. High frequency techniques, such as statistical energy analysis (SEA), often miss important information such as dominant resonance behavior due to stiff or small scale parts of the structure. Hybrid methods circumvent this problem by coupling FEM/BEM and SEA models in a given built-up structure. In the approach adopted here, the whole system is split into a number of subsystems that are treated by either FEM or SEA depending on the local wavelength. Subsystems with relative long wavelengths are modeled using FEM. Making a diffuse field assumption for the wave fields in the short wave length components, the coupling between subsystems can be reduced to a weighted random field correlation function. The approach presented results in an SEA-like set of linear equations that can be solved for the mean energies in the short wavelength subsystems.     

Efficient parametric uncertainty analysis within the hybrid Finite Element/Statistical Energy Analysis method

This paper is concerned with the development of efficient algorithms for propagating parametric uncertainty within the context of the hybrid Finite Element/Statistical Energy Analysis (FE/SEA) approach to the analysis of complex vibro-acoustic systems. This approach models the system as a combination of SEA subsystems and FE components; it is assumed that the FE components have fully deterministic properties, while the SEA subsystems have a high degree of randomness. The method has been recently generalised by allowing the FE components to possess parametric uncertainty, leading to two ensembles of uncertainty: a non-parametric one (SEA subsystems) and a parametric one (FE components). The SEA subsystems ensemble is dealt with analytically, while the effect of the additional FE components ensemble can be dealt with by Monte Carlo Simulations. However, this approach can be computationally intensive when applied to complex engineering systems having many uncertain parameters. Two different strategies are proposed: (i) the combination of the hybrid FE/SEA method with the First Order Reliability Method which allows the probability of the non-parametric ensemble average of a response variable exceeding a barrier to be calculated and (ii) the combination of the hybrid FE/SEA method with Laplace's method which allows the evaluation of the probability of a response variable exceeding a limit value. The proposed approaches are illustrated using two built-up plate systems with uncertain properties and the results are validated against direct integration, Monte Carlo simulations of the FE and of the hybrid FE/SEA models.     

Extension of SEA model to subsystems with non-uniform modal energy distribution

In order to widen the application of statistical energy analysis (SEA), a reformulation is proposed. Contrary to classical SEA, the model described here, statistical modal energy distribution analysis (SmEdA), does not assume equipartition of modal energies.Theoretical derivations are based on dual modal formulation described in Maxit and Guyader (Journal of Sound and Vibration 239 (2001) 907) and Maxit (Ph.D. Thesis, Institut National des Sciences Appliquées de Lyon, France 2000) for the general case of coupled continuous elastic systems. Basic SEA relations describing the power flow exchanged between two oscillators are used to obtain modal energy equations. They permit modal energies of coupled subsystems to be determined from the knowledge of modes of uncoupled subsystems. The link between SEA and SmEdA is established and make it possible to mix the two approaches: SmEdA for subsystems where equipartition is not verified and SEA for other subsystems.Three typical configurations of structural couplings are described for which SmEdA improves energy prediction compared to SEA: (a) coupling of subsystems with low modal overlap, (b) coupling of heterogeneous subsystems, and (c) case of localized excitations.The application of the proposed method is not limited to theoretical structures, but could easily be applied to complex structures by using a finite element method (FEM). In this case, FEM are used to calculate the modes of each uncoupled subsystems; these data are then used in a second step to determine the modal coupling factors necessary for SmEdA to model the coupling.     

Application of a Bayesian algorithm for the Statistical Energy model updating of a railway coach

The classical statistical energy analysis (SEA) theory is a common approach for vibroacoustic analysis of coupled complex structures, being efficient to predict high-frequency noise and vibration of engineering systems. There are however some limitations in applying the conventional SEA. The presence of possible strong coupling between subsystems and the lack of diffuseness result in a significant uncertainty. This is the main motivation for the present study, where a procedure to update SEA models is proposed. The proposed procedure is the combination of the classical SEA method and a Bayesian technique. Due to reasons such as finding a limited number of important parameters, using a limited search range, avoiding matrix inversion and taking the effect of noise into account, the proposed strategy can be considered as a proper alternative to the experimental SEA approach. To investigate the performance of the proposed strategy, the SEA model updating of a railway passenger coach is carried out. First, a sensitivity analysis is carried out to select the most sensitive parameters of the SEA model. For the selected parameters of the model, prior probability density functions are then taken into account based on published data on comparison between experimental and theoretical results, so that the variance of the theory is estimated. The Monte Carlo Metropolis Hastings algorithm is employed to estimate the modified values of the parameters. It is shown that the algorithm can be efficiently used to update the SEA models with a high number of unknown parameters.     

Ensemble energy average and energy flow relationships for nonstationary vibrating systems

This paper attempts to introduce a new point of view on energy analysis in structural dynamics with particular emphasis to its link with uncertainty and complexity. A linear, elastic system undergoing free vibrations, is considered. The system is subdivided into two subsystems and their respective energies together with the shared energy flow are analysed.First, the ensemble energy average of the two subsystems, assuming uncertain the natural frequencies, is investigated. It is shown how the energy averages follow a simple law when observing the long-term response of the system, obtained by a suitable asymptotic expansion. The second part of the analysis shows how the ensemble energy average of a set of random samples is representative even of the single case if the system is complex enough.The two previous points, combined, produce a result that applies to the energy sharing between two subsystems even independently of uncertainty: for complex systems, a simple energy sharing law is indeed stated. Moreover, in the case of absence of damping, a nonlinear relation between the energy flow and the energy (weighted) difference between the two subsystems is derived; on the other hand, when damping is present, this relationship becomes linear, including two terms: one is proportional to the energy (weighted) difference between the two subsystems, the other being proportional to its time derivative. Therefore, the approach suggests a way for deriving a general approach to energy sharing in vibration with results that, in some cases, are reminiscent of those met in Statistical Energy Analysis.Finally, computational experiments, performed on systems of increasing complexity, validate the theoretical results.     

2D quantum gravity from quantum entanglement

In quantum systems with many degrees of freedom the replica method is a useful tool to study the entanglement of arbitrary spatial regions. We apply it in a way that allows them to backreact. As a consequence, they become dynamical subsystems whose position, form, and extension are determined by their interaction with the whole system. We analyze, in particular, quantum spin chains described at criticality by a conformal field theory. Its coupling to the Gibbs' ensemble of all possible subsystems is relevant and drives the system into a new fixed point which is argued to be that of the 2D quantum gravity coupled to this system. Numerical experiments on the critical Ising model show that the new critical exponents agree with those predicted by the formula of Knizhnik, Polyakov, and Zamolodchikov.     

Sound radiation from a vibrating plate with uncertainty

Sound radiation into open space from a vibrating structure has been investigated since Rayleigh. On the other hand the sound power transferring into neighboring reverberant subsystems has also been rigorously studied using statistical energy analysis, particularly for the high frequency range. Falling between the two well-known problems, pressure and intensity fields from the sound radiation have not yet been widely studied using statistical methods. In this paper, the sound radiation from a vibrating thin plate having uncertain dynamic properties is investigated. Estimates are developed for the reverberant vibration field in the uncertain plate subjected to a point-excitation, and for the ensemble average of pressure from the direct field and from the reverberant field, leading to an estimate of the average sound intensity. The power radiated from the plate and the radiation efficiency is also derived. Monte Carlo simulations are conducted with an ensemble of plates with randomly-distributed point masses, and the simulation results compare well with the estimates.     

On the reciprocity relationship between direct field radiation and diffuse reverberant loading

This analysis is concerned with the derivation of a "diffuse field" reciprocity relationship between the diffuse field excitation of a connection to a structural or acoustic subsystem and the radiation impedance of the connection. Such a relationship has been derived previously for connections described by a single degree of freedom. In the present work it is shown that the diffuse-field reciprocity relationship also arises when describing the ensemble average response of connections to structural or acoustic subsystems with uncertain boundaries. Furthermore, it is shown that the existing diffuse-field reciprocity relationship can be extended to encompass connections that possess an arbitrary number of degrees of freedom. The present work has application to (i) the calculation of the diffuse field response of structural-acoustic systems modeled by Finite Elements, Boundary Elements, and Infinite Elements; (ii) the general calculation of the Coupling Loss Factors employed in Statistical Energy Analysis (SEA); and (iii) the derivation of an alternative analysis method for describing the dynamic interactions of coupled subsystems with uncertain boundaries (a generalized "boundary" approach to SEA).     

Prediction of airborne sound transmission across a timber–concrete composite floor using Statistical Energy Analysis

This paper concerns the development and experimental validation of prediction models using Statistical Energy Analysis (SEA) to calculate the airborne sound insulation of a timber–concrete composite floor. The complexity in modelling this floor is due to it having (1) a multilayer upper plate formed from concrete and Oriented Strand Board (OSB), (2) multiple types of rigid connector between the upper plate and the timber joists and (3) a resiliently suspended ceiling. A six-subsystem model treats the concrete–OSB plate as a single subsystem and three different five-subsystem models treat the combination of concrete, OSB and timber joists as a single orthotropic plate subsystem. For the orthotropic plate it is suggested that bending stiffnesses predicted using the theories of Huffington and Troitsky provide a more suitable and flexible approach than that of Kimura and Inoue. All SEA models are able to predict the weighted sound reduction index to within 2 dB of the measurement. The average difference (magnitude) between measurements and predictions in one-third octave bands is up to 4 dB. These results confirm that SEA can be used to model direct transmission across relatively complex floor constructions. However, this requires the inclusion of measured data in the SEA model, namely the dynamic stiffness of the resilient isolators and the cavity reverberation time.     

Direct determination of phase coexistence properties of fluids by Monte Carlo simulation in a new ensemble

A methodology is presented for Monte Carlo simulation of fluids in a new ensemble that can be used to obtain phase coexistence properties of multicomponent systems from a single computer experiment. The method is based on performing a simulation simultaneously in two distinct physical regions of generally different densities and compositions. Three types of perturbations are performed, a random displacement of molecules that ensures equilibrium within each region, an equal and opposite change in the volume of the two regions that results in equality of pressures, and random transfers of molecules that equalize the chemical potentials of each component in the two regions. The method is applied to the calculation of the liquid-gas coexistence envelope for the pure Lennard-Jones (6, 12) fluid for several reduced temperatures from the vicinity of the triple point to close to the critical point (T? = 0.75 to T? = 1.30). Good overall agreement with previously available literature results is obtained, with some deviations at the extremes of this temperature range.     

Estimating ensemble average power delivered by a piezoelectric patch actuator to a non-deterministic subsystem

Engineering systems such as aircraft, ships and automotive are considered built-up structures. Dynamically they are taught of as being fabricated from many components that are classified as ‘deterministic subsystems’ (DS) and ‘non-deterministic subsystems’ (Non-DS). Structures' response of the DS is deterministic in nature and analysed using deterministic modelling methods such as finite element (FE) method. The response of Non-DS is statistical in nature and estimated using statistical modelling technique such as statistical energy analysis (SEA). SEA method uses power balance equation, in which any external input to the subsystem must be represented in terms of power. Often, input force is taken as point force and ensemble average power delivered by point force is already well-established. However, the external input can also be applied in the form of moments exerted by a piezoelectric (PZT) patch actuator. In order to be able to apply SEA method for input moments, a mathematical representation for moment generated by PZT patch in the form of average power is needed, which is attempted in this paper. A simply-supported plate with attached PZT patch is taken as a benchmark model. Analytical solution to estimate average power is derived using mobility approach. Ensemble average of power given by the PZT patch actuator to the benchmark model when subjected to structural uncertainties is also simulated using Lagrangian method and FEA software. The analytical estimation is compared with the Lagrangian model and FE method for validation. The effects of size and location of the PZT actuators on the power delivered to the plate are later investigated.     

Vibration of L-shaped plates under a deterministic force or moment excitation: a case of statistical energy analysis application

Analytical and closed form solutions are presented in this paper for the vibration response of an L-shaped plate under a point force or a moment excitation. Inter-relationships between wave components of the source and the receiving plates are clearly defined. Explicit expressions are given for the quadratic quantities such as input power, energy flow and kinetic energy distributions of the L-shaped plate. Applications of statistical energy analysis (SEA) formulation in the prediction of the vibration response of finite coupled plate structures under a single deterministic forcing are examined and quantified. It is found that the SEA method can be employed to predict the frequency averaged vibration response and energy flow of coupled plate structures under a deterministic force or moment excitation when the structural system satisfies the following conditions: (1) the coupling loss factors of the coupled subsystems are known; (2) the source location is more than a quarter of the plate bending wavelength away from the source plate edges in the point force excitation case, or is more than a quarter wavelength away from the pair of source plate edges perpendicular to the moment axis in the moment excitation case due to the directional characteristic of moment excitations. SEA overestimates the response of the L-shaped plate when the source location is less than a quarter bending wavelength away from the respective plate edges owing to wave coherence effect at the plate boundary.     

Investigation of structural response and noise reduction of an acoustical enclosure using SEA method

Noise reduction capability of acoustical enclosures has been extensively investigated and is commonly quantified by noise reduction. Although much work has been done in predicting the noise reduction by using statistical energy analysis (SEA) method, discrepancy often exists between the predicted and measured results, especially at the low and intermediate frequencies. This is largely due to the fact that some physical mechanisms which affect the prediction of noise reduction were not addressed in the available SEA models. In this paper, the sources of the discrepancy were identified by investigating the limitation of SEA for energy transfer in the entire frequency range and the effect of structural–structural coupling, acoustical–structural coupling on the prediction of noise reduction. An improved SEA model including the non-resonant response and more accurate transmission coefficient of finite panels is presented. Finally, the predicted structural response and noise reduction of an acoustical enclosure were compared with measured results, and the improved agreement between the prediction and experiment was observed.     

Quantum Chaos,Random Matrices,and Irreversibility in Interacting Many-Body Quantum Systems

The Pauli master equation describes the statistical equilibration of a closed quantum system. Simplifying and generalizing an approach developed in two previous papers, we present a derivation of that equation using concepts developed in quantum chaos and random-matrix theory. We assume that the system consists of subsystems with strong internal mixing. We can then model the system as an ensemble of random matrices. Equilibration results from averaging over the ensemble. The direction of the arrow of time is determined by an (ever-so-small) coupling to the outside world. The master equation holds for sufficiently large times if the average level densities in all subsystems are sufficiently smooth. These conditions are quantified in the text, and leading-order correction terms are given.