`Mathematical Problems in EngineeringVolume 2010, Article ID 907232, 33 pageshttp://dx.doi.org/10.1155/2010/907232`
Research Article

## Solution of the Dynamic Interaction Problem between a Framed Structure and an Acoustic Cavity Using Imposed Deformation Functions at the Interface

Department of Civil Engineering, University of Brasília, 70919-970 Brasília, DF, Brazil

Received 14 August 2009; Revised 17 February 2010; Accepted 23 March 2010

Copyright © 2010 Paulo Marcelo Vieira Ribeiro and Lineu José Pedroso. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This article presents an analytical procedure for solution of the dynamic interaction problem of a vibrating framed structure connected to a bidimensional cavity, containing an acoustic fluid. Initially the pressure solution for the fluid domain is developed, using the separation of variables technique. In a next step, this solution is applied to an entirely open cavity and to a closed cavity in the transversal direction, both containing a vibrating boundary with an arbitrary deformation. The generalized parameters of the structure (mass, rigidity, and force) are obtained by means of the virtual work principle, with the generalized force represented by the dynamic pressures acting on the interface. The dynamic equilibrium equation of the system is established for an imposed deformation, making a parametric study of the involved variables possible. Finally, it is demonstrated that this procedure can be generalized, allowing the construction of practical abacuses for other boundary conditions of both the structure and the cavity, and that these results allow a reasonable interpretation of the coupling regions, including the prediction of added mass and added stiffness effects, as well as corresponding frequencies and mode shapes of the coupled problem.

#### 1. Introduction

The problems of dynamic fluid-structure interactions are of great interest to a wide range of research fields in engineering. In these cases, the movement of these two subsystems is not independent and is governed by dynamic contact conditions. According to Bathe [1], among the categories of coupled problems, the interaction between a structure and an acoustic fluid is distinguished for its simplicity. The inevitable contact of the structures with the acoustic medium (air, water, etc.) makes the acoustic-structural interaction relevant in the analysis of many physical systems [2].

The problems of fluid contained in enclosures are common in many practical applications. Great interest has been devoted to the acoustic comfort of passengers in aircraft or automobile cabins, for example. Additionally the case of water cooled nuclear reactors, rocket propeller tanks, submerged sonar domes, and reservoirs can also be cited as relevant examples of this category of problem [3].

Fluid-structure interaction problems involving acoustic cavities have already been solved by various methods, including analytical, semianalytical, and numerical approaches [49]. In the majority of these works the acoustic medium was represented by the air. Among the studies involving fluid contained in tanks or reservoirs, one can cite the studies of [1013], which evaluated the interaction effects between the fluid and a vibrating column, and [1416], which established the coupling effects between a plate and the surrounding fluid. In some of these works the problem was solved analytically, with the determination of frequencies and mode shapes of the coupled problem. However, little attention has been given to the development of abacus and graphical representations that allows the interpretation of such phenomena, including the identification of added mass and added stiffness regions, as well as physical interpretation of the corresponding mode shapes and its relation with the cavity modes.

One of the greatest advantages of the analytical treatment of coupled problems is the analysis of dimensionless parameters, which allow physical interpretation of the solutions. However, for more complex models, such set of solutions is often not possible. These are applied only to some specific cases, while numerical procedures can handle mode general problems. However, the implementation of numerical solutions is time consuming both in construction and in computer processing. According to Amabili and Kwak [15], analytical procedures show its importance in the solution of simple cases and can be used in the validation of numerical solutions.

In this paper an analytical procedure applied to a framed structure coupled to a bidimensional acoustic cavity is developed for solution of frequencies and mode shapes of the coupled problem. The separation of variables technique is employed, resulting in a general pressure solution, which is then applied to specific cases of an entirely open cavity and a closed cavity in the transversal direction, both containing a vibrating boundary. The association of this movement with the related structural vibration function, and the introduction of the dynamic fluid pressures as external forces, enables the construction of the dynamic equilibrium equation of the coupled structure, which is solved for frequencies of the equivalent system. While this procedure presents the limitation of prior knowledge of the boundary deformation, the mathematical simplicity of the solution justifies its application, resulting in equations represented with clear and well-defined parameters that allow the construction of abacuses for interpretation of frequencies and mode shapes.

#### 2. Fluid Domain General Solution

The basic assumptions of this problem are the corresponding to treatment of this medium as an acoustic fluid. With these considerations, it is assumed that the fluid transmits only pressure waves. Some applications of this theory include propagation of pressure waves in pipes and sound waves propagating through fluid-solid media. The following basic assumptions are made for the solution of this problem (based on the propositions of Chopra [17] and Rashed [18]).(i)The fluid is homogeneous, inviscid, and linearly compressible.(ii)The flow is irrotational.(iii)Displacements and their derivatives are small.(iv)Surface wave effects are neglected.(v)The movement of the fluid-structure interface is bidimensional (the same for any vertical plane perpendicular to the structural axis).(vi)The fluid-structure interface is vertical.(vii)Interface displacements are represented by an arbitrary deformation function.

The previous assumptions lead to a dynamic pressure distribution , in excess of the static pressure, given by which corresponds to the bidimensional wave equation, where is the bidimensional Laplacian operator and represents the fluid sound velocity, where indicates the fluid bulk modulus and its density. For an incompressible fluid and therefore thus (2.1) is reduced to This indicates Laplace’s equation governing the dynamic pressures in an incompressible fluid. It is evident that this expression is a particular case of (2.1).

Solution of (2.1) is achieved using the separation of variables technique. Therefore it is assumed that this expression can be separated, resulting in Substitution of (2.4) in (2.1) provides where the tiles indicate derivatives related to the corresponding variable. Division of (2.5) by results in Analysis of expression (2.6) indicates that the left-hand side of this equation depends only on , while the right-hand side depends on and . Since this equation must be satisfied for any values of , and , it is necessary that both sides are equivalent to an arbitrary constant , which can assume positive, negative, or null values. Thus, This last expression provides three ordinary differential equations: where and represent arbitrary separation constants. Equations (2.9) and (2.10) provide the following relation: The combined solutions of the differential equations (2.8) to (2.10) provide the complete solution of (2.1). However, it is important to notice that each one of these expressions requires two constants, resulting in a total of six unknown constants. A simple and time independent solution for this problem can be achieved with the hypothesis of time harmonic vibrations, with frequency Thus it is assumed that the time-related function is given by The hypothesis of time harmonic travelling waves establishes an important relation between the separation constants in the and directions. Substitution of (2.12) in (2.10) provides Substitution of (2.13) in (2.7) gives Thus From this last expression it can be concluded that Equation (2.16) indicates that the hypothesis of time harmonic travelling waves results in a separation constant in the direction that depends on the following parameters: frequency (), fluid sound velocity (), and direction separation constant ().

The solution of (2.1) can be established with the solution of the two ordinary differential equations given by (2.14) and (2.15). The separation constant can assume positive, negative, or null values. However, it should be noticed that depending on the problem, some of these solutions may result in expressions that are either trivial or without physical meaning.

Table 1 presents a resume of the possible solutions in both and directions according to the value of . Values of are also included.

Table 1: Resume of possible solutions according to .

#### 3. Solution of Rectangular Cavities with a Vibrating Boundary

##### 3.1. Entirely Open Cavity

In this case it is assumed that the cavity has the following pair of boundary conditions in the and directions, respectively: vibrating boundary-open and open-open. The vibrating boundary condition is related to the fluid-structure interaction at the interface, while the open condition implies at zero pressure at the contour. It will be shown later that this cavity, despite the lack of physical meaning, presents mathematical simplications that justify the use of such conditions, providing a useful interpretation of the phenomena. Figure 1 illustrates the analyzed domain as well as the boundary conditions.

Figure 1: Analyzed domain scheme including boundary conditions.

It is assumed that the vibrating boundary horizontal acceleration is governed by a time harmonic function, related to an arbitrary shape function and a maximum amplitude . Thus, for a given point at the interface the corresponding horizontal acceleration will be given by Equations (2.16) and (3.1) can be substituted at the interface boundary condition . Thus Therefore The remaining boundary conditions are easily established, indicating a zero pressure at the contour. Therefore, in order to avoid the trivial solution, must be null at these locations. Thus Equations (3.3) and (3.4) define the four boundary conditions needed for the solution of this problem. Table 1 will be used as guide for selection of the corresponding longitudinal and transversal solutions. Nontrivial solutions occur only when For this value the transversal solution is given by where, The corresponding longitudinal solution for this value of is given by where the separation constant is defined by The complete solution is given by the sum of every possible solution of . Thus The remaining constant is obtained with application of boundary condition (3.3), using the sine function orthogonality property. Thus Equation (3.10) can be substituted in (3.9) resulting in This last expression represents the dynamic pressure field solution for an open bidimensional cavity containing a vibrating boundary subjected to a harmonic motion, described by an arbitrary shape function .

Application of Sinusoidal Deformation Functions on Boundary
The sine function orthogonality property can be applied once again, if it is assumed that the vibrating boundary has the following shape function: which could be associated (in the case where the vibrating boundary represents deformations of a framed structure) to the corresponding mode shapes of a simple beam. In this case, (3.11) will result in Defining , where is an arbitrary value, and substituting this expression in (3.8) Therefore, applying (3.14) in (3.13) gives where This last expression can be rewritten in terms of a longitudinal function a transversal function and an amplitude factor Thus It is important to notice the presence of a cosine function in . This term will lead expression (3.17) to infinite results whenever the value of this trigonometric function approaches zero. This condition establishes critical points that result in resonance responses of the acoustic cavity. These points are defined by This last expression indicates that every value of is associated with infinite critical values that will lead the pressure response to infinity. These values will always occur for . The limit of when is given by Equations (3.18) and (3.19) are equivalent, respectively, to the frequencies and mode shapes of an acoustic cavity, closed-opened in the longitudinal direction and opened-opened in the transversal direction . Therefore, it can be concluded that critical values represent a set of frequencies that will lead to the mode shapes of the associated cavity, with the vibrating boundary condition at replaced by a rigid wall. Table 2 indicates the first seven uncoupled modes for cavities with and . The corresponding mode shapes of the vibrating boundary for a given value of are also indicated.
Analysis of Table 2 results indicates that a given mode with an arbitrary value of will produce only circular zones in the transversal direction. For example, solutions with two circular regions in are always associated to In the longitudinal direction the corresponding relation is given by . Thus, for , for example, only one and a half circular regions are expected at the direction. Cavities with a greater value of will present mode sequences with lower values of while smaller values of will provide sequences with higher values of this parameter.

Table 2: Uncoupled mode shapes sequence for the associated cavity.
##### 3.2. Closed Cavity in the Transversal Direction

In this case boundary conditions at and are replaced by Nontrivial solutions are defined by the following values of : For boundary conditions given by (3.21) provide a valid solution, given by Table 1. Thus Table 1 provides the correspondent longitudinal solution, defined by Therefore, the complete solution for is given by where indicates a remaining constant. Another valid solution is provided by In this case the transversal solution is given by where It is important to notice that in (3.25) provide a valid solution, given by (3.22). Therefore, solutions for are included in these last expressions, when The correspondent longitudinal solution is defined by where the separation constant is defined by The complete solution for is given by The remaining constant is obtained with application of boundary condition (3.3), using the cosine function orthogonality property. Thus,

Application of Sinusoidal Deformation Functions on Boundary
Assuming that the vibrating boundary is governed by (3.12) the following conditions are established by (3.30) integrals: In expression (3.32) it must be observed that even values of result in nonzero solutions only for If is an odd number, nontrivial solutions are expected only for Therefore, it can be concluded that summation, indicated on (3.29), does not vanish on this case (probably the greatest disadvantage when compared to an entirely open cavity), with remaining values of providing valid solutions. Dynamic pressure solution is defined by The complete solution is given by sum of expressions (3.33) and (3.34). For even values of expression (3.33) vanishes and solution is defined only by (3.34), which assumes nonzero values for For odd values of expressions (3.33) and (3.34) are present, with this last one assuming nontrivial solutions for
As in the previous case (entirely open cavity), critical values are also present in expressions (3.33) and (3.34), which are equivalent to the uncoupled cavity frequencies, with replaced by a rigid wall. These are defined by Limiting configurations (equivalent to the associated cavity modes) are established in expressions (3.33) and (3.34), when . Table 3 presents these results for , with corresponding values of and . It is important to notice that odd values of are related to symmetrical distributions with respect to . The opposite occurs for even values of . Thus it can be concluded that symmetrical structural modes lead to symmetrical distributions of cavity pressures with respect to the vertical mid section.

Table 3: Uncoupled mode shapes sequence for the associated cavity ().

#### 4. Fluid-Structure Coupled Solution for an Imposed Deformation

The analyzed problem is illustrated on Figure 2. It consists of a general structure with a constant cross-section subjected to an external load. The dynamic response of this system can be represented by a generalized coordinate , allowing the construction of generalized parameters (mass, stiffness, and loading) for any arbitrary mode shape, related to . This type of solution will be very useful for the introduction of fluid pressures, since the previous developed approach for the fluid domain is also dependent on the shape function.

Figure 2: Representation scheme of the structural model.

For the mathematical development of this problem the following considerations are assumed: mass per unit length , flexural rigidity , length , external distributed loading and unitary width perpendicular to the plane. The deflections are represented by , related to an arbitrary coordinate and a mode shape function , normalized at the generalized coordinate location. Therefore, The system’s dynamic equilibrium equation is obtained with the virtual work principle, equaling the work done by internal and external forces. Thus, Finally, introducing the following notations: and substitution of these last expressions in (4.2) provide which represents the system’s dynamic equilibrium equation of motion in terms of a generalized coordinate , where , and are related, respectively, to generalized mass, generalized stiffness, and generalized loading. It should be noticed that the generalized stiffness term includes only bending deformation effects. Additional contributions could be included by means of modification of this parameter. Damping could also be included, and in this case it would be more convenient to express this effect using a damping ratio (). Thus

##### 4.1. Application of Fluid Pressures at the Interface

The previous solution depends on an external loading. In the case of a coupled system this loading is represented by the dynamic pressures acting at the interface. Therefore where the function is related to the corresponding cavity (associated to the framed structure), with boundary conditions depending on the analyzed solution for the fluid domain. For a harmonic vibrating contour accelerations at the interface are given by Accelerations at the interface are equivalent for both the structure and the fluid. Therefore, Analysis of (4.9) and (4.10) provides Substitution of (4.5) and (4.8) in (4.6) gives In the above expression the generalized force is located at the left-hand side, because physically the dynamic pressure acts in the same direction of inertia and elastic forces. Equation (4.11) can be applied in (4.12) leading to a simplified expression given by Equation (4.13) represents the free vibration of the structural model, with a generalized mass produced by the interaction between fluid-solid domains. This expression can be simplified with the inclusion of a generalized added mass term, which corresponds to sum of the structural and fluid dislocated masses. Thus, Expression (4.11) provides Substitution of (4.15) in (4.14) results in For a nontrivial solution the term in brackets of (4.16) must be null. Therefore, Solution of the above expression provides frequencies of the coupled problem. It should be noticed that the generalized stiffness depends on . And the total generalized mass is composed of two parts. The first one is structure related, being dependent on . The second one is associated to the fluid dislocated mass and depends on and , which are unknown parameters of the problem, corresponding to the structural coupled mode shape and the coupled system frequency, respectively. Therefore, this type of solution establishes only one equation and two unknown variables. A simplified solution is possible with the introduction of an imposed deformation function at the interface. Thus, a frequency equation for a corresponding mode shape is constructed and the corresponding set of solutions is obtained. Latter these values can be applied in the total generalized added mass expression, resulting in the system’s dynamic equilibrium equation of motion, which can be solved for an arbitrary excitation. Or, they can be substituted in the pressure field solution, resulting in the cavity coupled mode shapes.

##### 4.2. Frequency Equation for an Open Cavity with a Sinusoidal Vibrating Boundary

In this specific case fluid pressure solution is given by (3.15). At the interface of a square cavity () this expression is reduced to Thus, the generalized added mass is given by This last equation represents the generalized fluid-added mass solution for a square open cavity with a harmonic vibrating boundary associated to . Therefore, the dynamic equilibrium equation of this system is given by substitution of (4.19) in (4.17): where the generalized parameters and are given by Substitution of (4.21) in (4.20) provides which indicates the frequency equation of the coupled problem. For the uncoupled case the second term in brackets vanishes and the corresponding solution is given by For the coupled case, solution of (4.22) is more complicated and includes infinite solutions for a given value of (as it will be demonstrated later). The second term in brackets is a function of and it is interesting to study the variation of this term with this parameter, which can be simplified to Figure 3 illustrates the variation of (4.24) along the axis for values of and . Analysis of this graphic indicates that the generalized fluid-added mass solution is hyperbolic and without critical values (resonances) in the interval. Occurrences of critical points are expected beyond this point, leading the generalized added mass to infinite values (due to the trigonometric nature of the solution). It is interesting to notice that values of produce almost constant functions, with defined values at .

Figure 3: Variation of along the axis for values of and

The limit when establishes which indicates a generalized fluid-added mass solution independent of . Values of this function can be evaluated for a given value of Table 4 illustrates these results. Reduced values of the added mass are expected with the increasing of

Table 4: Values of function for

It is also possible to rewrite (4.22) in terms of the uncoupled frequencies of both the structure and the associated cavity. Thus, where is given by (4.23) and , which corresponds to the first transversal frequency of the uncoupled cavity. Solutions of (4.26) are established when the term in curly brackets vanishes. Therefore, where the following dimensionless parameters are defined: and . This last term can be rewritten, resulting in with and . Therefore, the coupled solutions are defined in terms of four parameters: the structural frequency of the corresponding mode shape in vacuum, the first transversal frequency of the uncoupled cavity, the density relation between fluid and structure, and the thickness/height ratio of the structure.

A parametric study of (4.27) is presented on Figure 4, for . The distributed plots are associated to function , which represents the second term in this expression. The coupled values of are established at the intersection of these curves with a corresponding horizontal line, which indicates the value of . Constant values of are illustrated on this figure, providing a reference for analysis of the involved parameters.

Figure 4: Parametric abacus of the frequency equation for and

Analysis of Figure 4 indicates the presence of two distinct regions: A and B. The first one is contained on the initial interval between and the first critical value (first resonance). The behavior of function at this region is given by curves with developing amplitudes towards infinity. The second region is composed by curves with amplitudes ranging from to . In a given point there is a common intersection, where all the functions share the same amplitude. In both regions an increase of results in greater horizontal distances between function and the vertical asymptotes (corresponding to the cavity resonance frequencies). This implies greater relative differences between the coupled problem solution and the corresponding uncoupled cavity frequencies. It should also be noticed that an arbitrary horizontal line , e.g.) will provide infinite solutions, intercepting curves more than once.

Limit zones are established by the solid curve which connects the common intersection points in Region B. This condition implies leading to Thus, solutions intercepted by curve will present coupled frequencies which are equal to the corresponding in-vacuum values. Additionally, this function defines solution zones of added mass and added stiffness, which are located, respectively, above and below this curve. An added mass region implies coupled frequencies inferior to the corresponding in-vacuum values, while an added stiffness implies the opposite. It should be noted that all solutions in Region A are of added mass. Therefore, the initial solution will always be of this type. For Region B both types of solutions are possible, including the common intersection points. Solutions in Region A are unique, while for Region B an infinite set of solutions is defined. The common points in Region B are given by Solution of this last expression results in Therefore, which provides the common points solution, which is independent of The corresponding solution of when condition (4.29) is established is given by Therefore, problems with condition (4.33) satisfied will present a single coupled frequency which is equal to the corresponding in-vacuum value. As mentioned before, this type of solution is always located on Region B.

Figures 5 and 6 present, respectively, the parametric abacuses for and . In both cases Region A is extended to the proximity of . Observations made for Figure 4 are still valid.

Figure 5: Parametric abacus of the frequency equation for and
Figure 6: Parametric abacus of the frequency equation for and

It is important to notice that previous analyses were concerned with square cavities (). However, the same procedure could be extended to a rectangular geometry. In this case the generalized added mass will be dependent on the parameter. Therefore, (4.18) can be replaced by The generalized added mass resulting from this last expression is given by Thus, the corresponding frequency equation results in The effects of in the parametric abacuses are related to the resonance values given by (3.18). An increasing value of this parameter results in a greater number of critical points at a given interval along the axis, resulting in more regions of type B (Figure 7). The first resonance is also influenced by this parameter, defining the horizontal extension of Region A. Decreasing the value of will result in a larger extension of Region A, and in a smaller number of regions of type B (Figure 8).

Figure 7: Parametric abacus of the frequency equation for and
Figure 8: Parametric abacus of the frequency equation for and .

The limit of expression (4.35) when is given by Equation (4.37) can be evaluated for given values of and . Table 5 illustrates these results. Reduced values of the added mass are obtained for (short cavities). However, these effects are stronger at smaller values of exerting little influence on higher modes. Cavities with (long cavities) present the same added mass for a given value of .

Table 5: Values of function for

Figures 7 and 8 illustrate, respectively, the effects of in the parametric abacuses for . Comparison with Figure 4 indicates that regions of type B are developed earlier for , with the second region of this type appearing before . The opposite occurs for , with these regions dislocated to higher values of .

##### 4.3. Frequency Equation for a Closed Cavity in the Transversal Direction with a Sinusoidal Vibrating Boundary

In this specific case fluid pressure solution is given by (3.33) and (3.34). At the interface of a square cavity () these expressions are reduced to The generalized added mass is given by This last expression can be evaluated using orthogonality properties of sine and cosine functions. Therefore, substituting (4.38) in (4.39) It should be noticed that conditions related to dynamic pressures are still valid for these generalized added masses. Thus, expression (4.40) vanishes for even values of , with (4.41) resulting in non trivial solutions for odd values of Table 6 presents these conclusions.

Table 6: Generalized fluid-added mass conditions.

As in the previous case, a frequency equation is established by with defining the following parameter:

#### 5. Application Examples and Results

Two application examples of the previous described procedures are presented on this item. Simple beams associated to acoustic cavities entirely open and closed in the transversal direction are solved analytically and these results are compared to a finite element analysis. Figure 9 illustrates the general analysis scheme with material and geometrical properties.

Figure 9: General analysis scheme, material and geometrical properties.
##### 5.1. Analysis 1—Entirely Open Cavity

In this specific case the structural elastic modulus and density are taken, respectively, as and Parameters and are given by Therefore, the above equations indicate that parameter is constant, while depends on , which is associated to the structural vibration mode. The first three structural in-vacuum modes have the following frequencies: The first transversal frequency of the uncoupled cavity is given by Thus, for the first three modes parameter presents the following values: The above values together with (5.2) can be applied in the parametric abacuses corresponding to the associated value of or substituted in (4.27) for solution of the coupled values of in a given mode shape. Coupled frequencies are given by Table 7 illustrates hydrodynamic pressure distribution, coupled frequencies, and boundary deformation associated to the first seven vibration modes, obtained analytically and numerically using ANSYS finite element code. The last two columns on the right side indicate, respectively, frequencies and mode shapes of the uncoupled cavity (which are equal to results presented on Table 2 for ). An alternative type of representation of these results is illustrated on Table 8, which presents the relative differences between coupled analytical and numerical frequencies, as well as the relative differences between analytical coupled and uncoupled cavity frequencies.

Table 7: Analytical and numerical solutions (Analysis 1).
Table 8: Relative differences between analyzed frequencies (Analysis 01).

Analysis of Tables 7 and 8 results indicates that modes 1, 2, and 6 do not present characteristics of cavity modes (or resonant modes). The remaining coupled modes demonstrate certain proximity with the corresponding cavity frequencies and mode shapes, presenting modes slightly dislocated in the horizontal direction when compared to the cavity solution. In these cases the relative differences between coupled frequencies and cavity values are inferior to 10%. Modes 1, 2, and 6 have coupled frequencies inferior to the corresponding in-vacuum values. Thus, it can be concluded that these modes are included in the added mass region of the abacuses (above the curve).

An alternative approach could be used for solution of coupled frequencies of modes 1, 2, and 6. Application of Table 5 and (4.37) provides approximate values of the generalized fluid-added mass, for solution of coupled frequencies located in Region A of the parametric abacuses. Therefore, coupled frequencies are given by Thus, This presents an excellent agreement with corresponding values indicated on Table 7. It should be noticed that better results are expected for . Therefore, solution for  =  is almost exact, while for  =  =  an error of 7% is encountered.

##### 5.2. Analysis 2—Closed Cavity in the Transversal Direction

In this specific case the structural elastic modulus and density are taken, respectively, as and  =  Parameter is given by The first three structural in-vacuum modes have the following frequencies: The corresponding values of parameter are given by Table 9 illustrates hydrodynamic pressure distribution, coupled frequencies, and boundary deformation associated to the first seven vibration modes, obtained analytically and numerically using ANSYS finite element code. Modes 1, 2, and 5 have coupled frequencies smaller than the corresponding in-vacuum values, indicating fluid-added mass effects. Modes 2, 4, and 7 demonstrate certain proximity with cavity results, with configurations slightly dislocated in the horizontal direction. In the remaining modes, where is an odd number, this observation is not very clear. Solutions of this type include cavity modes associated to null values of . Therefore, modes 1 and 3 are related to the first cavity resonance (), while modes 5 and 6 are related to the third cavity mode ().

Table 9: Analytical and numerical solutions (Analysis 2).

It is important to notice that in this case fluid effects can modify structural in-vacuum mode shapes, as it can be noticed on modes 3, 6, and 7. The proposed procedure remains valid if a dominant configuration is still given by simple beam mode shapes, such as in modes 3 and 7. However, mode 6 presents an exception, where no dominant configuration is identified.

##### 5.3. Results Discussion for Analysis 1

The previous results indicate the division of coupled modes in two distinct categories. In the first type, coupled frequencies and mode configurations are very different from the corresponding uncoupled cavity solutions (illustrated on Table 2). In the second type, coupled modes show great resemblance with the uncoupled cavity configurations, with solutions dislocated in the horizontal direction and with a certain proximity in frequency values (which depends on both and , with smaller differences expected for solutions of located near the resonances). The first type has characteristics of added mass, with coupled frequencies smaller than the corresponding structural in-vacuum values. Solutions of this type are located on Region A in the parametric abacuses. The second type presents the possibility of coupled frequencies smaller (added mass), greater (added stiffness), or even equal to the corresponding structural in-vacuum solutions. Solutions of this type are always located on Region B.

An important conclusion can be established based on the previous results. A given mode , located on Region A, with hyperbolic solution, will present given by From this last expression it is possible to conclude that structural modes with coupled frequencies inferior to the right side of (5.12) provide typical added mass solutions. That is, coupled modes were the structure exerting higher influence in the response. It is also important to notice that this limit is given by the uncoupled cavity frequencies in the transverse direction.

Typical added mass solutions are given by the modified first cavity resonance for a given , with a hyperbolic solution in the longitudinal direction. Therefore, for these modes it is expected a decay in the dynamic pressure solution towards . Table 10 illustrates the possible configurations for typical added mass modes related to to .

Table 10: Possible configurations for the typical added mass modes.

Configurations of the coupled modes can be divided into distinct regions, according to the value of . Values of define hyperbolic variations, with the resulting solution located before the first resonance. Solutions in this region are classified as typical added mass modes (Region I). Values of define trigonometric solutions, with the corresponding mode shape resembling the first resonance configuration slightly dislocated to the left (Region II). When tends to the corresponding mode shape approaches the first resonance configuration. For the mode shape resembles the first resonance configuration dislocated to the right (Region III-a). When tends to the corresponding mode shape approaches a null pressure condition at the interface, defining the first zero point. Table 11 illustrates these domains for and , resulting in .

Table 11: Domains division for and

Table 12 illustrates the next sequence, with and , resulting in . It is important to notice that Region III is divided into two categories. The first, classified as “a”, is contained in the interval between the first resonance and the first zero. The second, classified as “b”, is limited by the first zero and the second resonance. From the first resonance and beyond only Region III is present. The zero point is defined as the value of which establishes a null pressure condition at the interface (). Analysis of expressions (3.15) and (3.17) indicates that this situation occurs whenever . Thus, Solution of the above expression results in which indicates exactly the same expression given by (4.32) when . Therefore, it is possible to conclude that the previous described common points in the parametric abacuses are equivalent to the zero points, defining conditions of null pressure at the interface and consequently coupled frequencies equal to the corresponding in-vacuum values.

Table 12: Domains division for and

It is also possible to divide the generalized fluid-added mass diagram of Figure 3 into the previously described regions. Figure 10 illustrates these results. Region II defines the beginning of the trigonometric solution and corresponds to the final interval before the first resonance. The first critical point establishes a signal change in the amplitude factor of (4.24), which becomes negative. This factor remains negative until the first zero point, defining Region III-a. From this point until the second resonance the amplitude factor is kept positive, defining Region III-b. After the second resonance the cycle restarts, alternating between Regions III-a and III-b. The negative values of have the physical meaning of an added stiffness, while the positive values define an additional mass. Therefore, solutions located on Regions I, II, and III-b provide the effects of an extra mass to the system, while solutions located in Region III-a provide an extra stiffness. Tables 11 and 12 indicate that solutions in Region III-a are given by critical configurations dislocated to the right, while solutions in Region III-b are equivalent to resonant modes dislocated to the left. Based on the previous analysis it is possible to conclude that solutions located at any of the previous described domains are equivalent to modified uncoupled cavity modes, with hyperbolic or trigonometric variations, and with corresponding configurations dislocated in the horizontal direction when compared to the uncoupled cavity solution.

Figure 10: Variation of along the axis for values of and
##### 5.4. Results Discussion for Analysis 2

This type of analysis presents an infinite summation in the generalized fluid-added mass term for even and odd values of . Therefore, even values of this parameter have a common set of resonances (antisymmetrical cavity modes). This is also true for odd values of (symmetrical cavity modes). The occurrence of this phenomena implies in composed modes, where a given structural mode shape has more than a single related to , with the moving boundary at resulting in a combination of even or odd values of . In some cases, especially in the first set of modes, a dominant mode is present and the proposed procedure is able to identify related frequencies (e.g. modes 1, 2, and 4,). However, the procedure fails if no significant dominance is present (as it can be noticed on mode 6). It is also important to notice that for an entirely open cavity (Analysis 1) there is no occurrence of common resonances, since summation is reduced to a non trivial solution only for , resulting in a unique set of critical values for each . Parametric abacuses could also be built for this type of analysis, resulting in similar conclusions to those identified in Analysis 1.

#### 6. Concluding Remarks

A closed analytical procedure for solution of frequencies and mode shapes of a simple beam connected to an open rectangular acoustic cavity was presented. The mathematical development resulted in a frequency equation defined by two dimensionless parameters, which enabled the construction of practical abacuses for solution and interpretation of the coupled problem. Analytical results are in agreement with analysis using the finite element method and indicate the presence of two distinct categories of solutions. The first type presents hyperbolic nature and a strong influence of the structure, producing typical added mass modes, with coupled frequencies inferior to the corresponding structural in-vacuum values. The second type presents a higher influence of the cavity modes, with coupled frequencies which can be smaller, greater, or even equal to the corresponding structural in-vacuum values. Mode configurations from this last category are equivalent to modified uncoupled cavity modes, with strong similarity in both modes and frequencies depending on the values of and ,

The proposed solution is useful for interpretation of the coupled phenomena, providing an understanding of the modes sequences and corresponding configurations. The resulting mathematical expressions applied at the practical abacuses allow the identification of these general conclusions, with an initial interval , resulting in added mass modes, followed by regions with , where the coupled solution can assume conditions related to the structural in-vacuum frequencies values, such as, added stiffness (higher), added mass (smaller), or a null pressure at the interface (equal). The limit curve identified by enables the identification of these regions. Solutions located above this function present characteristics of added mass modes, while solutions located under this curve present an added stiffness behavior. For a given mode only a single hyperbolic solution is possible (typical added mass mode), followed by an infinite set of trigonometric solutions (modified cavity modes).

As an extension, a simplified solution for a closed cavity in the transversal direction was also developed. The mathematical expressions indicate that summation is present in this type of solution, defining, respectively, symmetrical and antisymmetrical cavity resonances for odd and even values of . Practical consequences of this phenomena imply composed mode shapes, with represented by a set of odd or even values of . The proposed procedure is able to identify coupled frequencies whenever a dominant mode is present. However, this technique fails when solutions are related to composed modes. Properties of an entirely open cavity, such as added mass, added stiffness, or null pressures at the interface are also valid in this case.

Although results have been presented for two specific cases, solutions involving other types of boundary conditions for both the structure and the cavity are also possible, resulting in similar expressions and abacuses for solution and interpretation of this phenomena. It is important to remember that the only limitation of this procedure is the prior knowledge of the imposed deformation functions at the interface.

#### Acknowledgment

The authors are grateful for the financial support provided by CNPq scholarships.

#### References

1. K. J. Bathe, “Fluid-structure interactions,” Mechanical Engineering Magazine, vol. 120, no. 5, pp. 66–68, 1998.
2. R. B. Davis, Techniques to assess acoustic-structure interaction in liquid rocket engines, Ph.D. thesis, Department of Mechanical Engineering and Materials Science, Duke University, 2008.
3. F. Fahy, Foundations of Engineering Acoustics, Academic Press, London, UK, 1st edition, 2001.
4. E. H. Dowell and H. M. Voss, “The effect of a cavity on panel vibrations,” AIAA Journal, vol. 1, pp. 476–477, 1963.
5. A. J. Pretlove, “Free vibrations of a rectangular panel backed by a closed rectangular cavity by a closed rectangular cavity,” Journal of Sound and Vibration, vol. 2, no. 3, pp. 197–209, 1965.
6. R. W. Guy and M. C. Bhattacharya, “The transmission of sound through a cavity-backed finite plate,” Journal of Sound and Vibration, vol. 27, no. 2, pp. 207–216, 1973.
7. E. H. Dowell, G. F. Gorman III, and D. A. Smith, “Acoustoelasticity: general theory, acoustic natural modes and forced response to sinusoidal excitation, including comparisons with experiment,” Journal of Sound and Vibration, vol. 52, no. 4, pp. 519–542, 1977.
8. J. Pan and D. A. Bies, “The effect of fluid-structural coupling on sound waves in an enclosure—theoretical part,” Journal of the Acoustical Society of America, vol. 87, no. 2, pp. 691–707, 1990.
9. K. L. Hong and J. Kim, “Analysis of free vibration of structural-acoustic coupled systems—part I: development and verification of the procedure,” Journal of Sound and Vibration, vol. 188, no. 4, pp. 561–575, 1995.
10. H. Goto and K. Toki, “Vibrational characteristics and seismic design of submerged bridge piers,” Memoirs of the Faculty of Engineering (Kyoto University), vol. 27, pp. 17–30, 1965.
11. C.-Y. Liaw and A. K. Chopra, “Earthquake response of axisymmetric tower structures surrounded by water,” UCB/EERC Report 73-25, Earthquake Engineering Research Center, University of California, Berkeley, Calif, USA, 1973.
12. Y. Zhu, Z. Weng, and J. Wu, “The coupled vibration between column and water considering the effects of surface wave and compressibility of water,” Acta Mechanics Sinica, vol. 10, pp. 657–667, 1989.
13. J. T. Xing, W. G. Price, M. J. Pomfret, and L. H. Yam, “Natural vibration of a beam-water interaction system,” Journal of Sound and Vibration, vol. 199, no. 3, pp. 491–512, 1997.
14. H. Lamb, “On the vibrations of an elastic plate in contact with water,” Proceedings of the Royal Society of London, vol. A98, pp. 205–216, 1921.
15. M. Amabili and M. K. Kwak, “Free vibrations of circular plates coupled with liquids: revising the lamb problem,” Journal of Fluids and Structures, vol. 10, no. 7, pp. 743–761, 1996.
16. Y. K. Cheung and D. Zhou, “Coupled vibratory characteristics of a rectangular container bottom plate,” Journal of Fluids and Structures, vol. 14, no. 3, pp. 339–357, 2000.
17. A. K. Chopra, “Earthquake response of concrete gravity dams,” Tech. Rep. EERC 70-1, Earthquake Engineering Research Center, University of California, Berkeley, Calif, USA, January 1970.
18. A. Rashed, Dynamic analyses of fluid-structure systems, Ph.D. thesis, California Institute of Technology, Pasadena, Calif, USA, 1983.