Abstract

This paper studies the vibrational behavior and far-field sound radiation of a submerged stiffened conical shell at low frequencies. The solution for the dynamic response of the conical shell is presented in the form of a power series. A smeared approach is used to model the ring stiffeners. Fluid loading is taken into account by dividing the conical shell into narrow strips which are considered to be local cylindrical shells. The far-field sound pressure is solved by the Element Radiation Superposition Method. Excitations in two directions are considered to simulate the loading on the surface of the conical shell. These excitations are applied along the generator and normal to the surface of the conical shell. The contributions from the individual circumferential modes on the structural responses of the conical shell are studied. The effects of the external fluid loading and stiffeners are discussed. The results from the analytical models are validated by numerical results from a fully coupled finite element/boundary element model.

1. Introduction

The reduction of noise emitted by maritime platforms is an important topic in naval research. A truncated cone is used to support the aft-shaft-bearing of the submarine propeller. The structural and acoustic responses of a conical shell are not as widely reported in literature as in the case of cylindrical shells. This is due to the increased mathematical complexity associated with the effect of the variation of the radius along the length of the cone on the elastic waves. Much of the earlier work on the equations of motion for conical shells and forms of solution with different boundary conditions has been summarized by Leissa [1]. The approximate locations of the natural frequencies have been found using the Rayleigh-Ritz method by several authors [2, 3]. Tong presented a procedure for the free vibration analysis of isotropic and orthotropic conical shells in the form of a power series [4]. Guo studied the propagation and radiation properties of elastic waves in a cone in vacuo [5]. For fluid loaded conical shells, Guo investigated fluid loading effects of waves on conical shells in which a perturbation technique is applied by introducing a small vertex angle parameter [6]. Caresta and Kessissoglou considered the fluid loading by dividing the conical shell into narrow strips and each of these strips is treated as a locally cylindrical shell [7]. Donnell-Mushtari and Flügge equations of motion were compared by Caresta and Kessissoglou [8]. It was shown that in order to use the Flügge equations of motion with the power series solution, an approximation of the shear force is required.

Many researchers have investigated approximate methods to solve the acoustic radiation from vibrating elastic structures. The boundary integral equation formulation is a popular method. A good review on the evolution of boundary element techniques was given by Chien et al. [9]. Coupled finite element/boundary element methods have been presented by Jeans and Mathews [10]. Caresta and Kessissoglou used the Helmholtz integral to calculate the far-field radiated sound pressure of a submerged cylindrical shell closed by two truncated conical shells [11]. Since the boundary element method (BEM) requires significant computational time and memory, many researchers have been working on improving the efficiency of BEM [12] or exploring other approximate methods [13, 14]. Wang et al. proposed an approximate method named Element Radiation Superposition Method (ERSM) to predict the radiated pressure based on the surface velocity [15]. The total radiation pressure equals the superposition of the radiated pressure from each surface element, which can be obtained through weighting the vibration velocity using the transfer quantity. The proposed method was found to be superior in terms of computational efficiency. ERSM can be implemented to solve various models by adopting regular baffles such as spherical, infinite cylindrical and spherical baffles.

In this work, the structural and acoustic responses of a stiffened conical shell submerged in a heavy fluid are presented. Section 2 describes the solution for the dynamic response of a conical shell, which is presented in the form of a power series. Fluid loading is taken into account by dividing the shell into narrow strips, similar to the procedures presented by Caresta and Kessissoglou [7]. The stiffening of the structure is considered by smearing the stiffeners over the conical shell surface, which results in orthotropic material behavior of the conical shell. In Section 3, the Element Radiation Superposition Method (ERSM) is introduced to solve the far-field pressure of a conical shell, based on the theory presented in [15]. The total acoustic pressure of conical shell is assumed to be the superposition of the radiated pressure of all divided elements, which is approximated using the analytical solution of the radiated pressure from a piston with a cylindrical baffle. In Section 4, the validity and convergence of the method to solve the structural responses and the far-field acoustic pressure of a stiffened conical shell are discussed. The contributions from the individual circumferential modes for different load cases are examined. The effects of the external fluid loading and the ring stiffeners are described.

2. Dynamic Model of a Submerged Stiffened Conical Shell

2.1. Dynamic Response of a Conical Shell

For a conical shell, the coordinate system is defined in Figure 1. and are the orthogonal components of the displacement in the and directions, respectively, and is the displacement normal to the shell surface. is the semivertex angle of the conical shell. is the mean radius of the shell. is the radius of the smaller end and is the radius of the larger end. is the radius of the cone at location . is the length of the cone along its generator. According to the Flügge theory, the equations of motion are given by [1] where and are the density and the thickness of the conical shell, respectively. , , , , , and are the forces for a conical shell and are given in the Appendix. For a conical shell with stiffeners, the effect of the stiffeners can be considered by averaging their properties over the surface of the conicalshell. The smeared approach can be used when the shells are reinforced by closely spaced uniform stiffeners [16]. Figure 2 shows a schematic diagram of the ring stiffeners. is the modulus of elasticity of the stiffeners. and are the thickness and height of the ring stiffeners, respectively. is the cross-sectional area of the ring stiffeners. is the distance from the centroid of the stiffener to the shell. is the moment of inertia of the stiffeners about line of the reference (the middle line of the shell). and / are the distances between the ring stiffeners along the axial and generator direction of the conical shell, respectively. Using the smeared approach, the thickness of the conical shell, , in (2.1) is replaced by , where . and are replaced by and , where [17] and is the strain in the shell and is the variation of the curvature. Consider  , , and . The validity of the smeared approach is discussed in Section 4.1.

When the fluid loading on the surface of the stiffened conical shell is considered, the equations of motion can be expressed as With the power series solution, the displacements of the conical shell are given by [7] where the vectors , , and are is the vector of the eight unknown coefficients

The fluid loading effects on the conical shell are taken into account using an approximate method by dividing the conical shell into narrow strips and each of these strips is treated as a cylindrical shell. The low frequency range is considered here, where the variation of the real wave number that is used to calculate the fluid loading parameter is approximately a linear function of [7]. A qualitative analysis on the number of cylindrical subdivisions required to represent the conical shell is given in [7]. A tolerance is introduced due to the assumption that the wave number does not change along the length of the cone [7]. Thus, the length of each segment of the cone is determined by . depends on the geometry of the model and the frequency range, which is analyzed in Section 4.1. The parts are connected by 8 continuity equations at each interface.

The external pressure on every part of the conical shell is where is the external pressure loading on the cylindrical shell, which can be approximated using an infinite model, as shown in [18], and is given by where is the nondimensional ring frequency, is the fluid wave number, and , are the density and speed of sound in the fluid, respectively. is the hull axial wave number and is the radian frequency. is the mean radius of the conical shell segment. is the Hankel function of order . is the modified Bessel function of order . and are their derivatives with respect to the argument.

2.2. Boundary Conditions and Excitation

The dynamic response of the conical shell for each circumferential mode number is expressed in terms of for each segment of the conical shell, which can be solved by continuity conditions between the segments and boundary conditions at the ends. The various boundary conditions for the conical shell are as follows [8].

Free end:

clamped end:

shear-diaphragm (SD) end: where, for the SD boundary condition, the change from the conical coordinate system defined in Figure 1 to the cylindrical coordinate system is considered. and are expressed as The forces and moments of a conical shell are given in the Appendix. The four boundary conditions at each end of the conical shell with the eight continuity equations at each junction of two conical segments can be arranged in matrix form , where .

Two excitation cases as shown in Figure 3 are considered. is the excitation along the direction. For the free end boundary condition at the end of the conical shell, the external force results in a modification of the force equilibrium at the end given by (9), which now becomes [11] Excluding the harmonic time dependency, multiplying (16) by , and integrating from to yielded where if and if . Similarly, is the excitation normal to the shell surface. In this case, (11) becomes Point excitation in any direction can be expressed by these two excitations. The FRF can be calculated for the case of a point force located at any junction. The axial and radial displacements are obtained from the components of the displacements,

3. Far-Field Sound Pressure

In this section, the theory for the boundary element method is used to demonstrate that the far-field acoustic pressure of the shell can be solved by superposing the pressure of the individual elements. An analytical solution of the acoustic pressure radiated by a cylindrical piston in a cylindrical baffle is used to obtain the pressure of conical elements.

3.1. Acoustic Transfer Vector

According to the Helmholtz boundary integral equation, the acoustic pressure radiated by the arbitrary surface can be expressed as [19] where is the position vector on the vibrating surface and is the position of the observation field point. and are the surface pressure and normal velocity, respectively. is free-space Green’s function and .

By discretizing the surface into elements, (20) can be expressed as Equation (21) can be rearranged as [15] where and are the vectors independent of the surface velocity. The relationship between the surface pressure vector and the normal velocity can be deduced as where the matrices , are independent of the surface velocity but depend on the shape of the structure and frequency. Equation (22) can be expressed as The vector is called Acoustic Transfer Vector (ATV), which is an array of transfer functions between the surface normal velocities and the pressure at the field point. The ATV depends on frequency, geometry, acoustical parameters of the fluid, and position of the observation point.

3.2. Element Radiation Superposition Method

The conical shell is divided into elements, . Each element of the ATV, , is equal to the acoustic pressure radiated by the corresponding surface element vibrating with unity surface velocity while the other elements are at rest. For the case where only one element vibrates with the velocity , the radiated pressure of the field point at is given by Therefore, given the velocities , the radiated pressure is Hence, the total radiation pressure is the superposition of the acoustic pressures radiated by the surface elements vibrating with a rigid baffle. Since the velocities can be obtained from the solution of dynamic response of the conical shell, the main task is to find the acoustic transfer values .

According to Wang [15], the ATV of the element on a cylindrical shell is equal to the acoustic pressure radiated by a baffled piston vibrating with unity velocity. The far-field acoustic pressure radiated by a cylindrical piston in a cylindrical baffle can be obtained by solving the Helmholtz differential equation in closed form, which is given by [20] where is the axial length of the piston, is the radial angle of the piston, is the axial distance between the center of the piston and the origin, and is the azimuth angle of the center of the piston. All of these parameters are shown in Figure 4(a). is the velocity of the baffled piston. For a velocity equal to unity, the ATV of the each cylindrical shell element can be obtained from .

The radiated pressure from each element of the conical shell can be approximated by the analytical solution of the radiated pressure of a piston in a cylindrical baffle, as shown in Figure 4(b). The radius of this cylindrical baffle is equal to the mean radius of the conical strip which includes this element. The ATV of the piston in the cylindrical baffle is solved using its local coordinate. The angle between the axis of the conical shell and the axis of the cylindrical baffle is the semivertex angle of the conical shell. In this way, the ATV elements of the conical shell can be obtained and the total radiated pressure of the conical shell is solved by summing up the products of the individual transfer quantities and the corresponding velocities.

4. Results for Structural and Acoustic Responses for a Conical Shell

The structural and acoustic responses for the conical shell under harmonic excitation are presented using the following parameters:  m,  m,  m, and  m. The stiffeners have a rectangular cross-section of 0.014 m 0.04 m ( m,  m) and are evenly spaced by  m. The material properties for steel are density , Young’s modulus , and Poisson’s ratio . Structural damping is introduced using complex Young’s modulus , where is the structural loss factor. The surrounding fluid has a density of and speed of sound . The structural responses are predicted at a point equal to the drive point of the excitation, as shown in Figure 3. The distance between the far-field point and the axis of the conical shell is 100 m. The circumferential angle of the observation point is and its axial position is equal to the axial position of the drive point. The boundary conditions considered here are free at the smaller end of the conical shell and clamped at the larger end.

A finite element model is developed in ANSYS to validate the results obtained from analytical model presented here. The conical shell model and the ring stiffeners are built using Shell 63 elements and Beam 188 elements, respectively. Fluid 30 and Fluid 130 elements are used to simulate the external fluid and the absorbing boundaries around the fluid domain. ANSYS is used to obtain the structural responses. SYSNOISE is then used to calculate the far-field pressure using the direct boundary element method.

4.1. Validity and Convergence

The validity of smeared approach for the stiffened conical shell and the convergence of dividing the conical shell are discussed in what follows. A model of the stiffened conical shell in vacuo is built to validate the smeared approach. It is observed in Figure 5 that the analytical results for the structural responses due to match the finite element analysis (FEA) results in the low frequency range. However, the accuracy decreases rapidly when the frequency is above 180 Hz, where the structural wavelength is of the same order as the distance between the stiffeners. It is further observed that the increase of the spacing and width of the stiffeners affects the accuracy of this method, as shown in Figures 6(a) and 6(b), respectively. To conclude, the smeared approach is only suitable when the conical shells are reinforced by closely spaced stiffeners, that is, is small, and sufficiently small stiffeners, that is, is small, in the low frequency range, where the wavelength is larger than the stiffener spacing.

As discussed in Section 2, the conical shell is divided into segments and the length of each segment is . For the model in this paper, . The result for the structural responses due to for is compared with the results for and in Figure 7. It can be seen that the result converged for  . Two additional cases of dividing the conical shell into parts with equal lengths are also considered in Figure 7. It shows that segments with fixed yield faster convergence than fixed length segments. The FEA result is also presented in Figure 7 to validate the converged analytical result.

The validity and convergence of the Element Radiation Superposition Method (ERSM) for solving the far-field acoustic pressure of a conical shell are discussed in what follows. Each conical segment with the length of is divided into 1, 3, and 5 parts in its axial direction, corresponding to 19, 37, and 95 parts, respectively, for and 19 conical segments. The total number of elements for the conical shell is 1140, 2220, and 5700, respectively, for a circumferential angle of 6° for each element. Each element can be regarded as a piston in a cylindrical baffle. The results in Figure 8 show that the acoustic pressure due to converges when there are 3 parts per conical segment. The BEM result is also shown in Figure 8. As discussed in Section 3, the pressure is obtained from the velocity and ATV. In this part only, for both methods, the velocities are analytical solutions of the structural responses. The difference between ERSM and BEM lies in the solution of the ATV. The ATV for ERSM is obtained from the analytical solution of the Helmholtz differential equation for a piston baffled in a cylindrical shell, which is used to approximate a piston baffled in a conical shell, while the ATV for BEM is obtained by the numerical solution of the Helmholtz boundary integral formulation. The results in Figure 8 indicate that the locations of the peaks are well matched but the errors of amplitudes of the peaks are less than 5% at lower frequencies and less than 10% at higher frequencies. However, ERSM has a much quicker solution than BEM, due to its avoidance of solving the Helmholtz boundary integral formulation. Therefore, ERSM is an effective method to solve the far-field pressure of conical shell.

4.2. Effect of Excitation on a Conical Shell

The structural and acoustic responses of the conical shell due to the excitation shown in Figure 3 are presented in Figures 9 and 10. The comparison between the analytical results and the FEA results indicates that the axial and radial displacements using these two methods agree well, except for some errors at higher frequencies. These errors arise due to approximate solutions for the stiffeners and the external fluid loading. The stiffeners are considered by averaging their properties over the surface of the conical shell. The fluid loading is taken into account by dividing the conical shell into narrow strips, which are treated as cylindrical shells and the external fluid loading is approximated using an infinite model. The far-field acoustic pressures obtained using the analytical and FE/BE methods are also in good agreement. Errors are caused by the velocities, which are gained from the results of structural responses, and the ERSM, as it has a different solution of the ATV from the BEM. Thus, it can be found that some errors exist, especially at higher frequencies, due to complexity and approximations of the pressure solution from ERSM. But they are acceptable and applicable in the preliminary research, in the case that the analytical methods are able to save much time of the calculations.

The contributions of the independent circumferential modes with mode order are marked in the figures of the axial and radial displacements. It is observed that, in both cases, the modes are predominantly excited, followed by and modes. Higher order circumferential modes do not significantly contribute to the structural responses. It is also found that both structural and acoustic responses due to the excitation are significantly higher than those due to the excitation.

4.3. Effect of External Fluid Loading

The effect of the fluid loading on the structural and acoustic responses for the two excitation cases is presented for the conical shell in air in Figures 11 and 12. It is observed that the effects of the external fluid for the two cases are similar. The fluid loading reduces the peaks of the axial and radial resonant frequencies, due to damping. Moreover, the resonant peaks shift to lower frequencies because of the mass-like effect of the fluid loading. The amplitude of the far-field pressure for the conical shell in water is significantly higher than the pressure for the conical shell in air, which indicates that the sound is transmitted more efficiently in water than in air.

4.4. Effect of Ring Stiffeners

The effect of the ring stiffeners on the structural and acoustic responses of the conical shell responses is presented in Figures 13 and 14. The ring stiffeners are observed to have a stiffening effect, resulting in an increase of the resonant frequencies for both excitation cases. The shift of the resonant peaks of the far-field pressure is not as significant as that for the displacements, which shows that the added stiffness has little effect on the radiated pressure. It is further observed that the amplitudes of the pressure from conical shells with and without ring stiffeners are similar, which means that the ring stiffeners do not significantly affect the transmission of sound into the acoustic far-field.

5. Conclusion

An analytical method to predict the structural and acoustic responses of a stiffened conical shell immersed in fluid has been presented. The solution for the dynamic response of a conical shell was presented in the form of a power series. Fluid loading on the shell was taken into account by dividing the shell into narrow strips and stiffeners were modeled using a smeared approach. The Element Radiation Superposition Method was used to solve the far-field pressure of a conical shell. Good agreement between the analytical and FEM/BEM results was observed. Two cases of excitation by considering forces acting in different directions at the smaller end of the conical shell were examined, in which the forces acted along the generator and normal to the surface of the shell. In both cases, modes were predominantly excited, followed by and modes. Both the structural and acoustic responses of the conical shell under excitation in normal direction were greater than for excitation along the generator. The effects of fluid loading and ring stiffeners were examined. The results showed that external fluid loading has both damping and mass-like effects while stiffeners mainly increase the stiffness of the conical shell.

Appendix

The forces and moments for a conical shell are given by [1] where , , and , respectively, are Young’s modulus, Poisson’s ratio, and shell thickness.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

All the work in this paper obtains great support from the National Natural Science Foundation of China (51179071) and the Fundamental Research Funds for the Central Universities, HUST: 2012QN056.