Acoustic Pressure Waves in Vibrating 3-D Laminated Beam-Plate Enclosures

The effect of structural vibration on the propagation of acoustic pressure waves through a cantilevered 3-D laminated beam-plate enclosure is investigated analytically. For this problem, a set of well-posed partial differential equations governing the vibroacoustic wave interaction phenomenon are formulated and matched for the various vibrating boundary surfaces. By employing integral transforms, a closed form analytical expression is computed suitable for vibroacoustic modeling, design analysis, and general aerospace defensive applications. The closed-form expression takes the form of a kernel of polynomials for acoustic pressure waves showing the influence of linear interface pressure variation across the axes of vibrating boundary surfaces. Simulated results demonstrate how the mode shapes and the associated natural frequencies can be easily computed. It is shown in this paper that acoustic pressure waves propagation are dynamically stable through laminated enclosures with progressive decrement in interfacial pressure distribution under the influence of high excitation frequencies irrespective of whether the induced flow is subsonic, sonic , supersonic, or hypersonic. Hence, in practice, dynamic stability of hypersonic aircrafts or jet airplanes can be further enhanced by replacing their noise transmission systems with laminated enclosures.


Introduction
The control of vibration and noise propagation from industrial plants, aircraft engines and noise generating machines has remained an active research area for several decades.Nonetheless, limited literature exist in the area of noisestructure dynamic interaction modeling.However, within the context of analytical and experimental studies in acoustic -structure dynamics, a number of investigations have been reported in [1][2][3][4][5][6][7][8][9].For these problems, analytical techniques were employed to study active control of acoustic-structure interaction in 2-D and 3-D enclosures.
In particular, Fang et al. [10], investigated acousticstructure interaction through a 3-D rectangular enclosure to improve control design analysis.Nevertheless, reliable low-order finite model of 3-D acoustic-structure interaction has posed a formidable challenge in the field, given the fact that model density and order increase from 1-D to 3-D configurations, especially for severe multiple input and output systems.In the field of acoustic-structure interaction, the scarcity of literature in analytical modeling of 3-D enclosures having complex boundary conditions has been attributed mainly to the intractability of dynamic equations which preclude closed form solutions.
For some cases, with simplified assumptions, closedform solutions are possible usually with rigorous mathematical intrigues and manipulations.One of such methods involves the computation of a transfer function from the frequency response data.Even at that, the limit of these identification techniques is that, the system identification algorithms work satisfactorily only for low-order systems and systems with separated nodes.
However, for 3-D acoustic enclosures, with vibrating boundary surfaces, these identification techniques are fronted with difficulty, especially with two-input one-output systems.Nevertheless, investigations of vibroacoustic pressure waves propagation through a 3-D acoustic enclosure having vibrating laminated boundary surfaces has not been widely reported in literature.Within the context of transverse vibrations of laminated structures, with non-uniform interfacial pressure distribution, some interesting results have been reported recently in [11][12][13][14][15].The effect of structural vibration on the propagation of acoustic pressure waves through a cantilevered 3-D laminated beam-plate enclosure is investigated.For this problem, a set of well posed partial differential equations governing the vibroacoustic pressure waves interaction phenomenon are formulated and matched for the various vibrating boundary surfaces.This paper is organized as follows.Section 1 introduces the problem under investigation within a general context.In the next section, the essential analytical mechanics is briefly reviewed.In Section 3, these relationships are incorporated into a simplified analytical model for the mathematical analysis of the vibroacoustic pressure waves problem.Section 4 is concerned with the analysis of the transmission intensity as modulated by interfacial pressure variation and excitation frequency.
In Section 5, simulated results for the vibroacoustic pressure waves, as influenced by the linear interface pressure variations across the axes of the vibrating boundary surfaces are analyzed whilst in Section 6, the paper ends with the summary and conclusion.

Problem Formulation
The problem here is to examine analytically, the effect of the nature of load, frequency variation and the pressure gradient on the acoustic pressure waves propagation.
A general theory of the energy dissipation properties of press-fit joints in the presence of Coulomb friction as originally developed by Goodman and Klumpp provides the basis for the physics of the problem.As illustrated in Figure 1(a) below the proposed 3-D enclosures is constructed with laminated elastic structure for each boundary surface.The contact conditions between the mating layers as itemised in Damisa et al. [11] hold, namely, (i) there is continuity of stress distributions at the interface to sufficiently hold the equivalent layers together both in the pre-and post-slip conditions, (ii) a stick elastic slip with presence of Coulomb friction occurs at the interface of the sandwich vibrating boundary to dissipate energy and does not remain constant as a function of some other variables such as spatial distance, time or velocity.
The formulated vibroacoustic wave equation is given by the relation while the governing equations for each vibrating boundaries Ω 1, Ω 2 , Ω 3 , Ω 4 are derived in the appendix.In the meantime, the 2-D views of the structural loadings through Ω 3 , Ω 4 are illustrated in Figures (1(b) and 1(c)) for brevity.

Mathematical Analysis of the Vibroacoustic Pressure Waves Equation
In this investigation we will simplify the solution of (1) by assuming that the prompting vibrating membrane is through domain Ω 1 while the boundary effects on the pressure waves propagation are restricted to the following domains (Ω 2 , Ω 3 , Ω 4 ).
Advances in Acoustics and Vibration       Under this circumstance, (1) takes the form, viz, ( Equation ( 2) can be further expressed as By introducing the triple Fourier Finite Sine Transform, namely where X denotes sin(nπx/a) sin(mπ y/b) sin(kπz/c), in conjunction with the Laplace transform, viz, Advances in Acoustics and Vibration Equation (3) in the transform plane, subject to zero initial conditions takes the following form, viz, where the following have been defined So that on appropriate substitution, ( 6) is now expressed in Fourier-Laplace transform plane as where k denotes (nπ/a(−1 ). Subject to the following pressure profiles prescribed in Olunloyo et al. [12], viz, This now allows us to obtain the vibroacoustic pressure in the transform plane as ∀P av (a, 0 To complete the solution of (10), we have employed the following double Fourier sine transform relations, namely, By employing ( 11), (13), and ( 14) subject to (A.5), (A.6), (A.7), (A.8), (A.9), and (A.10), W 1 FxFy Fz (λ n , λ m , λ k , s) can be evaluated from (A.1) as Advances in Acoustics and Vibration where w denotes 6 F/EbH z 3 , and S denotes (1+(−1) k+1 )/kπ.So that on appropriate substitution into (10), the vibroacoustic pressure can be evaluated in the transform plane as where z denotes ρ 2 (s 2 + ω 2 0 ), where the following have been defined, viz, 4  , where By invoking the Fourier inversion, the solution of ( 16) in the Laplace transform plane can be written as where H denotes sin(mπx/a) sin(nπ y/b) sin(kπz/c).
For the limit case as Mach number, M → 0, the well known result in Laudau and Liftsitz [16] viz , is correctly recovered.By employing the Laplace inversion, the closed form solution for the vibroacoustic pressure wave can be computed as where r denotes (1 + (ψ 1 /2)), (32)

Mathematical Analysis of Sound Intensity as Influenced by Interfacial Pressure Distribution
In acoustics, sound intensity        energy per unit area in a direction normal to the area.By combining the fundamental equations governing sound field, mass continuity, the relation between sound pressure and density change in conjunction with Euler's equation of motion, leads to the equation relating sound intensity with the net flow of energy, viz, where (E k•e + E p•e ) is total acoustic energy.
In this paper, we shall limit analysis to intensity variation along the spatial variable x and for this purpose; the forgoing equation is modified to the following form, namely, Vibro-acoustic pressure   where P is the vibroacoustic pressure; A c is the cross sectional area of the 3-D enclosure; c m is the velocity of sound through the 3-D enclosure; ζ is a dummy variable of integration.Equation ( 31) can be normalized into to the form in conjunction with the ongoing mathematical representations, viz, Advances in Acoustics and Vibration In view of the above relations, the magnitude of the sound intensity can be evaluated from (35) as (37)

Analysis of Results
In this paper, vibration and noise propagation control emanating from complex engineering systems such as industrial power plants, aircraft engines, space propulsive devices and machine enclosures is investigated.The acoustic-structure configuration of interest is the one in which an acoustic disturbance is prompted by one of the vibrating boundaries of the enclosure such as in aircraft cabin noise transmission and systems for outer space exploration.The modeling techniques employed for this study derives from recent advances made in the mechanics of sandwich structures, with non-uniform interfacial pressure distribution.The acousticstructure dynamic interaction problem is simplified by Vibro-acoustic pressure assuming zero initial conditions prior to the excitation of the upper boundary surface.In the formulated problem, the upper layer of the sandwich elastic plates is subjected to harmonic excitation force as illustrated in Figure 1.To demonstrate the practical significant of our model problem, we have utilized some characteristic values listed in [10] as shown in the Canadian Society for Mechanical Engineers (CSME), Transaction of Applied Mechanics for simulation.Figures 2, 3, 4, 5, and 6 illustrate the characteristic acoustic natural frequencies as couched in (29).In fact, our expression allows for the incorporation of Mach number and for the case where the Mach number is approaching zero, Landau and Lifshitz [16] characteristic frequencies equation is recovered as a special case.As can be seen, the ordering of the family of curves clearly displayed the existence of subsonic, sonic, supersonic and hypersonic zones.
In each zone, the general pattern of the frequency profiles showed that the results are dependent on the geometry of the 3-D enclosure and ambient sound velocity.With respect to the effect of modal parameters, Figure 2, showed the profiles for the principal modes as the axial length increases from 0.4 to 1.2 m whilst b and c are restricted to the values listed in Table 1.
In the subsonic zone, the natural frequency profiles are inversely proportional to the axial length of the enclosure and decrease monotonically to a constant value irrespective of the axial length as we approach sonic flow.However, in the supersonic zone, the natural frequencies are further attenuated in reversed order to their respective critical values before increasing monotonically in the hypersonic zone.
With respect to the effect of higher modal parameters on the natural frequency profiles, as demonstrated in Figures 3-6, we observed similar pattern but with higher magnitude.
With respect to the profiles of vibroacoustic pressure, Figures 7,8,9,10,11,12,13,14,15, and 16 displayed the various pictures for low and very high frequency excitations.In general the pressure wave profiles are strongly influenced  by the values of the pressure gradients and the frequency of the vibrating boundary surface.Interestingly, several profiles can be simulated by playing with the values of ψ 1 and ψ 2 irrespective of the frequency of excitation, whereas the choice of any suitable combination can be arranged by progressively varying the tightening torques along the laminated surfaces.Such pressure gradients selection can have significant effect on the vibroacoustic pressure levels.
The pictures in Figures 7-10 illustrate the acoustic pressure profiles at the end of one cycle of vibration.
With respect to the effect of low excitation frequency through the boundary surface, the interfacial pressure gradients do not seem to have any strong effect on the acoustic pressure profiles.In fact we note in the lower range of Mach number, that the acoustic pressure waves rise initially before dropping at Mach 2 and increased marginally.However, as we move beyond Mach 16, the effects of the pressure gradients are noticeable.As can be observed, the pressure profiles are increasing with time.In particular, we note that pressure profile magnitude with lower boundary excitation frequency are significantly higher compare with the effect of higher boundary excitation frequency.On the other hand, Figure 13, showed the acoustic pressure profiles for an excitation frequency of 1 MHz.In respective of the pressure gradients combination, the profiles are monostatically stable.As the boundary excitation frequency increases to 1 GHz, the pressure profiles are significantly lower in the fundamental mode of vibration.
Advances in Acoustics and Vibration However, this may not be the case as we vary the values of the modal parameters by comparison of Figures 14-16.This suggests that in practice, dynamic stability of hypersonic aircrafts or jet airplanes can be further enhanced by replacing their noise transmission systems with laminated enclosures.

Summary and Conclusion
In this paper, explicit closed form solutions for the vibroacoustic characteristic frequencies, pressure waves and sound intensity or transmission quality through a 3-D enclosure with a vibrating laminated boundary surface is investigated.The vibroacoustic properties are shown to be dependent on, interfacial pressure gradients and boundary surfaces excitation frequencies.The results presented in this study can be positively exploited for the design of modern airplanes, aeroelastic structures and propulsive devices for the launching of space systems.

Appendix
For Case Ω 1 , For Case Ω 3 (mirror reflection of Ω 1 ), For Case Ω 4 (mirror reflection of Ω 2 ), Now for the 3-D cantilever enclosures under investigation, the following boundary conditions hold for each of the vibrating membranes, viz, On the other hand the effects of the excitation forces at the free end of the 3-D cantilever enclosures can be captured via the following end conditions [11].

a Ω 3 b 3 F 4 = 4 (
(b) 2-D view of the composite structure geometry through Ω F 04 e iw4t a Ω c) 2-D view of the composite structure geometry through Ω 4

I 8 Advances
, is a vector quantity and can be defined as the time average of the net flow of sound
W 3 (a, 0, t) = W 3 (a, b, t) = 0, ∀z = −h (A.15) (v) ∂W 3 (0, 0, t) ∂y = ∂W 3 (0, b, t) ∂y = 0, ∀z = −h Modulus of rigidity of the laminate F:Applied end force amplitude x:Space coordinate along the length of the laminate y:Space coordinate along the width of the laminate z:Space coordinate along the height of the laminate Dynamic response in Laplace transform plane W F : Dynamic response in Fourier transform plane W F : Dynamic response in Fourier-Laplace transform plane μ:C o e fficient of friction at the interface of sandwich composite elastic beam ρ 0 : Ambient air densityρ = ρ 0 −1 : t: Time coordinate p 0 :Clamping pressure at the interface τ (xz) 1 :Shear stress at the upper half of the laminates interface τ (xz) 2 :Shear stress at the lower half of the laminates interface

Table 1 :
Parameters for simulation of results.