This book outlines many aspects of vibroacoustic simulation. We made the journey from lumped elements to distributed systems, from deterministic and random over to hybrid systems models. The idea was to present the full frequency portfolio of simulation methods but based on simple cases that are comprehensive but of practical relevance. It is clear that we had to stay with manageable systems and examples for clarity. Even the industrial cases in Chapter 12 can only provide some aspects of complex models.
The goal was to provide a comprising description of deterministic and random methods as far as the link between them, the hybrid FEM/SEA theory. To the author’s experience, most acousticians are exclusively focussing on one method. Or even worse, they claim that either deterministic or random methods are the best.
Hopefully after this lecture it becomes clear that there is no good or bad method. There is only an appropriate frequency range and sometimes – when the conditions of hybrid FEM/SEA methods are predominant – both approaches are required within one model.
However, for simple systems constructed by plates, cavities, and noise control treatment a complete toolbox is provided to perform simulation on such systems over the full audible frequency range. The code and libraries of the examples can be found on the author’s website: www.docpeiffer.com.
This textbook can cover only a part of the extensive subject of vibroacoustic simulation, and the author would like to propose to the reader some literature and topics for further reading and studies. The list is by no means complete but should give suggestions to either develop an interest for the ongoing research in the field or for more details of more precise material and system descriptions.
Experiments and tests are mandatory for every acoustician. They provide the link to the real world and validate the models that are used. There is sometimes a ridiculous competition between people doing simulation or experiments that does not make any sense. Experimental people claim to be more practical than theorists, and the simulation people may feel superior due to the fact that they can manage complex simulation tools of mathematical equations.
It is the author’s conviction that both groups require a deep knowledge of the other’s competence to become a true acoustic expert. An experimentalist who is just comparing curves without any idea about the theoretical reason for the differences may make the wrong design decisions. A numerical model helps a lot to understand the vibroacoustic system well and to draw the right conclusions. On the other hand, a theorist that is never visiting a lab will not develop any feeling for practical problems and limitations. She or he may also miss important source and transmission phenomena in the model. In addition modern test analysis methods are so advanced and complicated that experimental people also require a deep theoretical knowledge.
This section provides some literature proposals of advanced experimental methods that are required to provide test results that can be used in model updates.
The transfer path analysis (TPA) is based on the measurement of the frequency transfer function of vibroacoustic systems. Hence, the dynamics of the systems is reduced to specific points of interest: mount connection, passenger ear positions, or excited locations. This transfer function is measured for multiple inputs (forces, moments, and acoustic volume strength) and responses (acceleration, displacement, and pressure). A good overview about the history and state of the art can be found, for example, in (van der Seijs et al., 2016).
The system is subdivided into subsystems that are connected via specific point shaped connections so that the transfer functions of the main paths are included. Those transfer functions are similar to the matrix (1.198) as introduced in section 1.7.
Both response and excitation are measured. For structural excitation this is done with piezoelectric transducers that measure the force generated by impact hammers or shakers and accelerometers for the response. Instead of accelerometers a scanning laser vibrometer can be used. Acoustic sources with given source strength are realized by special loudspeaker systems sometimes connected to long tubes in order to establish a point source. The pressure response is determined by microphones.
The determination of the transfer functions is performed using broadband signals (white noise, sweeps, or the given hammer impact history) and by applying equation (1.194). Hammer measurements are usually averaged over several impacts.
The transfer function of the subsystem is then inverted to get the stiffness or mobility matrix of the subsystems of interest. Especially this inversion is the critical point in TPA because noisy signals and test errors may lead to badly conditioned matrices and erroneous results. Many papers deal with methods to increase the robustness of this inversion by using pseudo inverses and additional reference sensors to create overdetermined systems.
In combination with an operational analysis, the TPA determines the forces occurring at the interface connections. From the forces and the local stiffness, the power input to the system can be calculated, and transfer paths are rated in terms of significance. The precision of these and power forces depends strongly on the quality of the matrix inversion and the test result in general.
When the systems become random, phase and magnitude of the system responses become uncertain. In this case the TPA can no longer be applied.
The mode shapes of structural or fluid systems can be measured directly. The tests and the related analysis methods are called experimental modal analysis(EMA). The system of interest is either equipped with a sufficient number of sensors to capture the modal shape and excited at one or a few excitation points. Instead, few response sensors and a large number of excitation positions can also be used. Due to reciprocity both methods provide the same result.
The are several methods to analyze measured frequency response functions and to extract the modal properties. A comprising overview about the existing methods is given by He and Fu (2001). The measured mode shapes can be used in modal frequency analysis similar to numerics, or the shapes are used as experimental reference. In this case the simulated modes are compared to the experimental shapes.
The major challenge is to estimate clean transfer functions for precise mode shapes, in particular for higher order modes. One advanced method is to estimate these functions by a best fit to a fraction of polynomial functions named PolyMAX parameter estimation (Peeters et al., 2004).
It is common practice that test and simulation results are compared based on point responses, e.g. pressure at driver’s ear from simulation and test. This might work for low frequencies and deterministic systems but may be misleading for complex systems. It is much better to compare the system responses for many positions distributed over this full system.
It does not make much sense to compare hundreds of response curves in one diagram. Thus, specific criteria are introduced to compare shapes, for example shapes of modes or frequency responses. Here, the most used criteria for mode comparison is given, the so called modal assurance criteria (MAC):
Under the condition of normal modes and perfect agreement the MAC is supposed to be orthogonal, hence for equal mode shapes and for different shapes. For the global agreement of the dynamic behavior the modal frequencies must also be equal.
The application of this criterion is state of the art in model updating methods and included in many test software packages. Nevertheless, there are many practical pitfalls. Note, that structural modes are mass normalized as given by equation (1.119) and that the mass matrix is usually not available in experimental modal analysis. So, for systems with very irregular system and mass distribution, this must be considered. This affects the orthogonality of the modes; thus, different modes don’t have zero MAC. See Allemang (2001) for more details. In addition the shape might be undersampled, meaning that the selected sensors are not dense enough to catch all details of the high frequency shapes.
There are several approaches besides engineering judgement to overcome such issues. One option is to check the so called AUTOMAC first. This is the MAC when the modes are correlated to themselves. If the result is the unit matrix, the selected test mesh and the mass distribution are supposed to be correct. A further issue is the case when modes have the same frequency and thus cannot be identified as distinct modes in the test. Obviously, only deterministic systems can be compared by this criterion.
When it comes to correlation of an SEA simulation with tests, things are getting more complicated. The SEA result is the average result of an ensemble of systems. Thus, curve comparison of single tests are reasonable when the systems are clearly random, as it is the case, for example, in the academic test cases of Chapter 10.
For the correlation of SEA models to tests, the system responses must be measured at many locations to get a good spatial average. Theoretically, even an ensemble of systems should be tested, which is not possible in most cases. Especially, single point response comparison are definitely no solution. However, when automatic test systems such as laser scanning systems or robots are used, a certain average can be achieved, and the comparison of physical units or the energy is possible. The investigation of the structure borne wave propagation along the aircraft fuselage in Chapter 12 is one example of this practice (Weineisen, 2014). A detailed correlation of the fuselage sidewall to SEA models is dealt with in Wang (2015).
The experimental SEA can be considered as a transfer path analysis of random systems. The concept of virtual SEA is to measure the input power (from force and acceleration) at several excitation points on each subsystem and the velocity or pressure response on several locations of the connected and responding subsystems. The energy of such systems is estimated from (6.116) and (6.118). For the latter equation assumptions for the mass distribution must be made.
We start from the nonsymmetric version of the SEA equation (6.102), so the the modal density is not required. When we assume that the power introduced into the th subsystem is given by and the energy response of the th subsystem due to this excitation is denoted by we get the following equation system for three systems
The subsystem energy is derived with (6.116) or (6.118). This can become a complicated task for irregular shapes and mass distribution, e.g. a ribbed plate of a fuselage panel (Bouhaj et al., 2017).
In equation (13.2) there are nine equations for nine parameters when we write the power load and energy matrices in one column. This system of equations is solved for . There are several investigations on experimental SEA, for example Lalor (1987), (1990), (1997), continued by the investigations on practical systems from Borello (2009), (2015)and Borello and Gagliardini (2007). Similar to the transfer path analysis, the determination of the coupling loss factors require an invertible matrix that is not always provided by tests. One key measure to overcome this is the method of averaging. Lalor proposed calculating the coupling loss factor from every measurement and average the coupling loss factors. Other authors average all excitation positions over one subsystem. Bouhaj et al. (2017) propose a statistical method to bring more robustness into the matrix inversion. In any case, virtual experimental SEA can be a promising tool to create SEA models from test data.
In addition the same method can be applied using FE simulations. In this case the introduced power and the energy estimation is determined from simulations. The biggest advantage here is, that the amount of sensors is nearly unlimited. The global analysis procedure is exactly the same. Such studies can be used to investigate, for example, very complicated junctions where no simplified equations for the radiation stiffness are available.
Note that the general decision how the substructuring is made is the duty of the experimentalist. Hence, also experimental SEA requires decisions about the SEA model set-up namely on the separation of the system into subsystems. This is further evidence for the fact that experimental people need theoretical knowledge.
Further reading on SEA is provided by Le Bot (2015) who introduces many aspects of the statistical concept of SEA and constitutes the link to statistical physics. Langley (2007), (2016)developed an extended version of the diffuse field reciprocity (7.19) based on a modal approach and for subsystems that don’t carry a perfectly diffuse wave field. An interesting field of research is the automatic recognition of components of a hybrid FEM/SEA model as developed by Langley and Kovalewski (2012).
Some research projects aimed at developing or collecting methods for the “mid frequency” range. One example is the Marie Curie project “MID-Frequency”– CAE-Methodologies for Mid-Frequency-Analysis in Vibration and Acoustics (Atak et al., 2012). A text book on current research in vibroacoustic simulation and many practical examples was edited by Hambric et al. (2016).
The work from Shorter (1998) and Mace (2005) provides a more sophisticated approach to determine the coupling loss factor using modal analysis of the full system and subsystems. In this case no power injection is required. Many software tools have implemented this method. Similar to virtual SEA the energy flow method (EFM) is a powerful tool to investigate complicated subsystems and junctions.
In this process of subsystem selection, it is also very important to know what is not an SEA subsystem. Tewes and Peiffer (2007) showed that the inclusion of deterministic subsystems into the SEA power balance leads to incorrect coupling loss factors in the experimental SEA process.
In the initial chapters of this book we stayed with canonical and simple systems, such as homogenous plates and cavities. They were perfect for the exemplification of the method. Unfortunately, reality is neither canonical nor simple. Thus, a next step is to include more realistic but usually also complicated subsystems into the determination of modal density, wavenumber, and damping. Typical candidates are specific lay-ups, for example composite or sandwich plates. Integrated structures can be corrugated or beaded plates, rib-stiffened plates, or extruded profiles as shown in Chapter 12.
In many cases an analytical formulation of such materials or shapes does not make sense. The required equations may become more complicated and computationally even more demanding than numerical models. However, some more advanced descriptions of subsystems are available that can be managed with reasonable effort. Note that the thin shell is not correct for high frequencies, too. Here, thick plate theory must be applied, introduced, for example, in Cremer et al. (2005).
All examples in this book are using isotropic materials, meaning that the mechanical or acoustical properties are independent from the orientation. The mechanical strength and superiority of carbon fiber or glass fiber reinforced plastics comes from the fact that highest strength is designed into the direction of the highest load. Hence, those materials must be anisotropic to fulfill this specification. Anisotropic materials are used in plates, shells, or beams. So, they are mostly used for example in plate theory.
In section 9.3.4 the limp fiber model was introduced, which is a useful description for porous material without frame dynamics. One reason for using this model was that the theory is not too complex. For some foam materials that are used in the automotive industry, for example as firewall mass-spring noise control treatment, the Biot (1962) model must be used. The consideration of such layers makes special layer formulations necessary plus interface conditions. The topic of porous materials is a research subject of its own and still an ongoing field especially in the experimental field. The book from Allard and Atalla (2009) gives a comprehensive treatment of this topic. There are many recent formulations (Horoshenkov et al., 2019; Lafarge et al., 1997) of the thermal and viscous dissipation during the wave propagation in the porous frame showing that the topic is still a field of basic research.
The simulation of coupled fluid structure models in combination with porous finite elements is still in the process of being implemented as standard procedure.
In section 3.7 we dealt with the in-plane forces (3.148) and moments (3.185). The linear laminate theory is well established in structural lightweight design and specific equations for bending and in-plane waves are elaborated, for example, in Baker and Scott (2016). Basically the behavior of a plate wave doesn’t change; the main difference to isotropic plates is that stiffness ratios between bending and in-plane strain can be different. In addition, the line moments and stresses can be coupled in case of unsymmetrical lay-ups.
More optimised plate structures in terms of lightweight design are sandwich plates that put most of the strength in the skin layer. The task of the very lightweight core material is to keep the distance between both skins. The dynamics of sandwich panels is special because the skin layers decouple from a certain frequency, so sandwich panels may be statically stiff and dynamically soft. This was shown by Kurtze and Watters (1959) and further detailed in (Moore and Lyon, 1991). This effect of decoupling depends mainly on the shear modulus of the core material. Note that due to the high stiffness and low weight, sandwich structures have the potential to become a vibroacoustic disaster because of a low and extended coincidence frequency range. However, by careful selection of the shear stiffness, this can be mitigated by dynamical decoupling.
In addition sandwich plates with heavy layers can constitute double walls due to the symmetric mode of wave propagation. In (Moore and Lyon, 1991) the transfer impedance for symmetric and nonsymmetric modes is given and can be used in the transfer matrix method to calculate the transmission loss.
In this book we dealt with flat plates. In most engineering applications shells with curvature are used. There are several detailed textbooks dealing exclusively with shells of isotropic material, for example the books from Leissa (1993) and Ventsel and Krauthammer (2001). Special shapes often used in analytical SEA are cylindrical shells, also called singly curved, or spherical shells called doubly curved. A treatment of natural modes of cylindrical (and ribbed) shells is given by Szechenyi (1971). The work of Koval deals with the transmission loss and extends this treatment from cylinders of isotropic material (Koval, 1976) including pre-stress from internal pressure (aircraft fuselage) to composite cylinders.
The wave finite element method can be considered as a bridge between FEM and SEA in that sense that small FE models are used to calculate the parameters of SEA. Please note that this doesn’t mean the EFM or virtual SEA methods. The idea of the WFE is to take a representative unit cell and apply periodic boundary conditions. With these assumptions quantities such as speed of sound, wavenumber, and damping can be calculated from this unit cell.
This is a perfect method for simple systems as plates (Manconi et al., 2013) but especially for periodic systems such as rib-stiffened shells or extruded profiles (Cotoni et al., 2008; Orrenius et al., 2014).
In addition the WFE can be extended to the calculation of coupling loss factors as shown by Mitrou et al. (2017). This further opens the path to a systematic derivation of the coupling loss factor.
In section 6.3.1.5 it was shown that there is an upper limit for SEA. This is the case when the dissipation along the wave path is so large that a diffuse field cannot exist, because the waves are damped out and reverberation is not possible. The frequency limit for systems with the maximum length with this condition is
This is the domain of methods that include the energy distribution in the systems. One option to take care of this is ray tracing, also called geometrical acoustics, as presented in (Pierce, 1991; Kuttruff, 2014). Examples for the application of ray tracing in the SEA context are presented by Gardner and Macarios (2014).
A different approach is introduced by Tanner. He developed an element based ray tracing method called dynamic energy analysis (DEA). In Tanner et al. (2012) the method is applied to plate systems.
All examples in this book are coded in python except those from Chapter 12. This led to an open source simulation software prototype for SEA and hybrid simulation that is presented at the author’s website www.docpeiffer.com. The author would highly appreciate further contributions from readers by creating their own test cases based on the example scripts. Further extensions to the above mentioned more advanced methods are intended in the near future. Contributions from the acoustic community are very welcome.