Underwater vibration and acoustic radiation calculation of double cylindrical shell by three-dimensional sono-elasticity of ships

Yangyang Zhang1 , Jingjun Lou2 , Xiang Yu3

1, 2, 3College of Power Engineering, Naval University of Engineering, Wuhan, Hubei Province, China

1Corresponding author

Vibroengineering PROCEDIA, Vol. 7, 2016, p. 18-24.
Received 24 June 2016; accepted 26 June 2016; published 31 August 2016

Copyright © 2016 JVE International Ltd. 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.
Creative Commons License
Table of Contents References
Cite this article
Views 182
Reads 69
Downloads 943
Abstract.

As common structure of underwater vehicle, the accurate prediction of underwater vibration and acoustic radiation of double cylindrical shell is always a difficult problem in the research of vibration and noise. Establish finite element model of double cylindrical shell in Abaqus software. Based on the three-dimensional sono-elasticity theory of ships, analysis vibration and acoustic radiation characteristics of underwater double cylindrical shell by the use of Thafts-acoustic software, and compared with experimental result. Results show that error of prediction result and experiment result at single frequency point is about 10 %~20 %, and error of total level of acceleration and sound pressure is less than 5 dB. It is accurate to predict the underwater vibration and acoustic radiation characteristics of double cylindrical shell with three-dimensional sono-elasticity theory of ships, which can provide new effective guidance for the engineering application.

Keywords: double cylindrical shell, vibration, acoustic radiation, prediction.

1. Introduction

Double cylindrical shell is one of the common structure of underwater vehicle. With the further improved requirement of modern war on acoustic stealth performance of underwater vehicle, the study of vibration and sound radiation prediction technology of underwater double cylindrical shell is becoming increasingly important. To this end, domestic and foreign scholars have carried out a lot of research work. From the shell vibration control equation, literature [1] derived finite or infinite cylindrical shell’s vibration and acoustic radiation analytic expression using modal expansion method; literature [2] derived vibration and acoustic radiation approximate analytical solution of finite stiffened cylindrical shell; literature [3] studied the acoustic radiation response of submerged finite cylindrical shell by the use of analytical and numerical methods; literature [4] established acoustic coupling equation of finite ring and longitudinal ribs stiffened shell by Flügge vibration theory, derived vibration and acoustic radiation analytical expression of stiffened shell under single-frequency multi-source excitation; literature [5] established the vibration and acoustic transfer path analysis (TPA) model of underwater structure, to effectively identify, quantify the main noise source and noise transmission path; literature [6] analysed the vibration and acoustic transfer path of double cylindrical shell under fluid - solid coupling condition by finite element method.

For double cylindrical shell with complex structure, analytical solution for solving vibration and acoustic radiation is almost impossible. In this article, three-dimensional sono-elasticity theory of ships is applied to analysis underwater vibro-acoustic characteristics of double cylindrical shell, predict the underwater vibro-acoustic characteristics with the help of software Thafts-acoustic by using numerical method and validate the accuracy of the prediction result by comparing with experimental result.

2. Model introduction

Double cylindrical shell is composed of outer shell, inner shell, ring ribs, end covers and elastic panel. Inner and outer cylindrical shell is connected via ring ribs and ring ribs are uniform along the cylindrical shell. Elastic panel locates in cylindrical shell’s horizontal middle position. The upper and lower side of the outer cylindrical shell is slotted, the water can fill the cavity between the inner and outer cylindrical shell when cylindrical shell is underwater. Model Parameters are shown in Table 1.

Establish half dry model and finite element model of double cylindrical shell using Abaqus software, which are shown in Fig. 1.

Table 1. Model parameters

Parameter
Value
Parameter
Value
Shell length L / m
2.0
Outer shell radius R / m
0.6
Ring rib spacing l / m
0.4
Inner shell radius r / m
0.5
Outer shell thickness c1 / m
0.004
Material density ρ / kg·m-3
7.8e+03
Inner shell thickness c2 / m
0.008
Elastic modulus E / n·m-2
2.1e+11
Elastic panel thickness c3 / m
0.008
Poisson’s ratio
0.3
End cover thickness c4 / m
0.012

Fig. 1. Cylindrical shell model

Cylindrical shell model

a) Geometry model

Cylindrical shell model

b) FEM model

3. Three-dimensional sono-elasticity calculation method of ships with double basin coupling

The schematic diagram of the acoustic elastic coupling decomposition of inner and outer flow field is shown in Fig. 2.

Fig. 2. Schematic diagram of acoustic elastic coupling decomposition calculation of the inner and outer flow field

Schematic diagram of acoustic elastic coupling decomposition calculation  of the inner and outer flow field

a) Coupling calculation model of outer flow field and shell

Schematic diagram of acoustic elastic coupling decomposition calculation  of the inner and outer flow field

b) Coupling calculation model of inner flow field and shell

Inner and outer flow field are not mutually connected aqueous medium. The inner and outer flow fields are regarded as two independent flow fields, and the attached water mass, attached water damping and generalized recovery coefficient caused by inner and outer flow field were calculated respectively. Fig. 2(a) shows the outer flow field and shell coupling mathematical model, where τ is the ideal sound medium outflow field, S¯ is outer flow field corresponding average wet surface, on which the unit normal vector is n¯; τi is imaginary inner flow field encircled by the wet surface Si1 which is coincident with S¯ and virtual impedance surface Si2. Fig. 2(b) shows the inner flow field and shell coupling mathematical model, where τb is inner flow field, corresponding to the average wet surface S¯b1 and Sd1; τcis imaginary outer flow field, corresponding to wet surface Sc which is coincident with S¯b1. Since the virtual outer flow field is semi-infinite flow field and there is no resonance problem, so there is no need to introduce a virtual impedance surface; τd is imaginary inner flow field encircled by S¯b2, wet surface Sd1 and a virtual impedance surface Sd2.

3.1. Coupling effect of flow field and shell

By the modal superposition method, nodal displacement of the double cylindrical shell structure discrete systems can be expressed as:

(1)
U t = D q t = r = 1 m D r q r t ,

where Dr is the modal displacement column vector corresponding to r-order modal, qt and qrt are respectively generalized primary coordinate column vector principal coordinate and r-order dry modal primary coordinate components.

Generalized kinetic equation for double cylindrical shell vibration can be expressed as:

(2)
a q ¨ t + b q ˙ t + c q t = f e t + f p t ,

where [a], [b], [c] were dry structure modal generalized mass matrix, generalized damping matrix and generalized stiffness matrix; fe, fp were generalized force column vector corresponding to outer excitation force in non-voice media flow field and generalized hydrodynamic column vector corresponding to acoustic medium flow field force on wet surfaces.

First, consider the outer flow field and shell coupling at zero speed. By adding the inner and outer flow field boundary integral equation, we can obtain:

(3)
1 4 π S ¯ σ r 0 G r , r 0 d s + 1 4 π S i 2 Z S ρ 0 u n i 2 r 0 G r , r 0 n i r 0 - i ω G r , r 0 u n i 2 r 0 d s = ϕ r ,
r τ S ¯ .

The boundary integral equation for solving source intensity can be expressed as:

(4)
ϕ r n r = - 1 2 σ r + 1 4 π S ¯ σ r 0 G r , r 0 n r d s - i ω 4 π S i 2 u n i 2 r 0 G r , r 0 n r d s
            + 1 4 π ρ 0 S i 2 Z S u n i 2 r 0 2 G r , r 0 n i r 0 n r d s ,         r S ¯ .

Boundary condition on the virtual impedance surface can be expressed as:

(5)
Z S u n i 2 r = 1 2 π S ¯ σ r 0 G r , r 0 d s - i ω 2 π S i 2 u n i 2 r 0 G r , r 0 d s
        + 1 2 π ρ 0 S i 2 Z S u n i 2 r 0 G r , r 0 n i r 0 d s ,         r S i 2 .

The r-order dry modal displacement vector of outer shell known as the boundary condition, substitute it into the boundary integral Eq. (4), and then combined with the Eq. (5) can obtain corresponding source intensity σr on the wet surface and vibration displacement uni2r on virtual impedance surface. Substitute these two quantities into the Eq. (3) can obtain velocity potential ϕrr of outer flow field corresponding to r dry modal.

Based on linear hypothesis of micro movements, omitting the second order small quantities, equivalent instantaneous wet surface integral to average wet surface integral, the attached water mass, attached water damping and generalized recovery coefficient contributed by outer flow field can be obtained in the frequency domain:

(6)
A r j ( 1 ) = ρ 0 ω 2 R e S ¯ i ω n ¯ u ¯ r ϕ j d s , B r j ( 1 ) = - ρ 0 ω I m S ¯ i ω n ¯ u ¯ r ϕ j d s C r j ( 1 ) = - ρ 0 S ¯ n ¯ u ¯ r g ω j d s . ,

In the formula above, ranges of subscripts of modal order were r= 1,..., m and j= 1,..., m.

Neglecting the influence of the inner flow field on the generalized recovery force coefficient, the generalized attached water mass and the attached water damping contributed by inner flow field can be obtained similarly:

(7)
A r j ( 2 ) = ρ 0 ω 2 R e S ¯ b 1 i ω n ¯ b u ¯ b r ϕ b j d s + S ¯ b 2 i ω n ¯ b u ¯ b r ϕ b j d s , B r j ( 2 ) = - ρ 0 ω I m S ¯ b 1 i ω n ¯ b u ¯ b r ϕ b j d s + S ¯ b 2 i ω n ¯ b u ¯ b r ϕ b j d s .

3.2. Acoustic elastic coupled solution

After taking the generalized hydrodynamic coefficients of the inner and outer flow field into account, we can obtain generalized kinetic equation of double basin acoustic elastic coupling:

(8)
- ω 2 a + A ( 1 ) + A 2 + i ω b + B ( 1 ) + B 2 + c + C ( 1 ) q
            = f e ω + T ω ,

where A1, B1 and C1 were respectively attached water mass matrix, attached water damping matrix and generalized recovery force coefficient matrix of dry modal contributed by outer flow field; A2 and B2 were respectively attached water mass matrix and damping matrix contributed of dry modal by inner flow field. According to the principle of modal superposition, the radiation acoustic velocity potential in the frequency domain can be obtained after the principal coordinate response of each dry modal is obtained:

(9)
ϕ R x , y , z = r = 1 m ϕ r x , y , z q r ,

where ϕrx,y,z can be obtained from Eq. (3). Combined with single-frequency harmonic time factor eiωt, substitute the principal coordinate response into Eq. (1), vibration response of ship structure in frequency domain can be obtained.

4. Prediction of underwater vibration and acoustic radiation of cylindrical shell

In order to verify the accuracy of three-dimensional sono-elasticity theory of ships in the prediction of underwater acoustic vibration characteristics of cylindrical shell, the result of Abaqus/Thafts-acoustic simulation and experiment were compared.

4.1. Underwater vibration prediction

The predicted vibration acceleration level of the test points on elastic panel and inner wall of the cylindrical shell are shown in Fig. 3 and Fig. 4.

Fig. 3. Vibration acceleration level and error on the elastic panel

Vibration acceleration level and error on the elastic panel

a) Acceleration level

Vibration acceleration level and error on the elastic panel

b) Error

Fig. 4. Vibration acceleration level and error on the inner wall of cylindrical shell

Vibration acceleration level and error on the inner wall of cylindrical shell

a) Acceleration level

Vibration acceleration level and error on the inner wall of cylindrical shell

b) Error

It can be seen from Fig. 3(a) that, in the 0-1000 Hz band, vibration acceleration level prediction and experimental test result on elastic panel of a cylindrical shell are close, and spectral line at some points like 42 Hz, 104 Hz and 455 Hz can corresponds good; from Fig. 3(b), we can see in the entire frequency band 0-1000 Hz, vibration acceleration level error of prediction and experimental test result are within 20 %, especially in 400-1000 Hz band prediction error is less than 10 %, the prediction result is very close to experimental result.

In 0-1000 Hz range, calculate the total acceleration level of simulation and experiment result. The total acceleration level of simulation result is 137.8 dB, while the total acceleration level of experiment result is 135.2 dB, the error is only 2.6 dB.

It can be seen from Fig. 4(a) that, in the 0-1000 Hz band, vibration acceleration level prediction and experimental test result trends of the test points on the inner wall of cylindrical shell are basically the same, prediction in the band 0-600 Hz was very close to experimental result; from Fig. 4(b), we can see in the entire frequency band 0-1000 Hz, vibration acceleration level prediction result of the test point is smaller than experimental result, positive error is small to 10 %, negative error is as greater as 20 %, but the overall error is within 20 %. The prediction result are accurate and available.

In 0-1000 Hz range, calculate the total acceleration level of simulation and experiment result. The total acceleration level of simulation result is116.1 dB, while the total acceleration level of experiment result is 120.0 dB, the error is only 3.9 dB.

4.2. Underwater acoustic radiation prediction

Test points collected from inner sound cavity and outer sound field of cylindrical shell, prediction and experimental result comparison of acoustic radiation is shown in Fig. 5 and Fig. 6.

Fig. 5. Sound pressure level and error in inner sound cavity of cylindrical shell

Sound pressure level and error in inner sound cavity of cylindrical shell

a) Sound pressure level

Sound pressure level and error in inner sound cavity of cylindrical shell

b) Error

Fig. 6. Sound pressure level and error in outer sound field of cylindrical shell

Sound pressure level and error in outer sound field of cylindrical shell

a) Sound pressure level

Sound pressure level and error in outer sound field of cylindrical shell

b) Error

It can be seen from Fig. 5(a) that, in the 0-1000 Hz frequency range, sound pressure level (SPL) prediction and experimental result trends of test points are consistent, so as resonance peak of some frequency points; from Fig. 5(b), we can see in the entire frequency band 0-1000 Hz, the SPL prediction and experimental test result error of cylindrical shell sound cavity are within 20 % and changes stably, the simulation result have high accuracy.

In 0-1000 Hz range, calculate the total sound pressure level of simulation and experiment result. The total sound pressure level of simulation result is125.6 dB, while the total sound pressure level of experiment result is 121.3 dB, the error is only 4.3 dB.

It can be seen in Fig. 6(a) that, in the 0-180 Hz band, the test point SPL of radiation outside the cylindrical shell underwater background noise described mainly, which shows that in lower frequency range radiated noise ratio is low; in 180-1000 Hz frequency, SPL prediction and experimental result trends of test points in outer sound field are consistent, but the prediction result is smaller than experimental result; from Fig. 6(b), you can see in 0-180 Hz frequency, the SPL prediction and experimental test result error of cylindrical shell in outer sound field are within 5 %, in the 180-1000 Hz band, the SPL radiation prediction error are positive error and less than 10 %, but negative error reaches 20 %. The possibly reason is background noise in outer sound field underwater is large, low signal to noise ratio result in test result is not accurate. But in general, acoustic radiation prediction of cylindrical shell underwater is of high accuracy.

In 0-1000 Hz range, calculate the total sound pressure level of simulation and experiment result. The total sound pressure level of simulation result is 119.1 dB, while the total sound pressure level of experiment result is 123.9 dB, the error is only 4.8 dB.

5. Conclusions

Establish finite element model of double cylindrical shell by software Abaqus. Based on three-dimensional sono-elasticity theory of ships, use software Thafts-acoustic to predict underwater vibration and acoustic radiation characteristics of double cylindrical shell and compared with experimental result. Results show that error of prediction result and experiment result at single frequency point is about 10 %-20 %, and error of total level of acceleration and sound pressure is less than 5 dB. Conclusions can be drawn as the following: the use of three-dimensional sono-elasticity theory of ships to predict underwater vibration and acoustic radiation characteristics of double cylindrical shell is of high accuracy, which is an effective means of underwater structure vibration and acoustic radiation prediction and can provide effective guidance for engineering applications.

References

  1. Junger M. C., Feit D. Sound Structures and Their Interaction. 2nd Ed. The MIT Press, Cambridge, 1986. [Search CrossRef]
  2. Tang Weilin, He Bingrong Approximate analytic solution of vibration and sound radiation from stiffened finite cylindrical shells in water. Acta Acoustic, Vol. 26, Issue 1, 2001, p. 1-5, (in Chinese). [Search CrossRef]
  3. Caresta M., Kessissoglou N. J. Structural and acoustic responses of a fluid-loaded cylindrical hull with structural discontinuities. Applied Acoustics, Vol. 70, 2009, p. 954-963. [Search CrossRef]
  4. Zhou Qizheng, Wang Deshi, Gao Shengyao Research on vibration and acoustic radiation of stiffened finite cylindrical shells in water under multi-excitations. Chinese Journal of Applied Mechanics, Vol. 30, Issue 5, 2013, p. 652-657, (in Chinese). [Search CrossRef]
  5. Zhang Lei, Cao Yue-yun, Yang Zi-chun Structural transfer path analysis for vibration and noise of a submerged double-shell. Journal of Vibration and Shock, Vol. 31, Issue 20, 2012, p. 12-16, (in Chinese). [Search CrossRef]
  6. Jin Guangwen, Zhang Linke, Miao Xuhong Vibration transmissibility of a submerged cylindrical double-shell based on reconstructing velocity field. Journal of Vibration and Shock, Vol. 30, Issue 5, 2011, p. 218-221, (in Chinese). [Search CrossRef]