Underwater vibration and acoustic radiation calculation of double cylindrical shell by threedimensional sonoelasticity of ships
Yangyang Zhang^{1} , Jingjun Lou^{2} , Xiang Yu^{3}
^{1, 2, 3}College of Power Engineering, Naval University of Engineering, Wuhan, Hubei Province, China
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 7, 2016, p. 1824.
Received 24 June 2016; accepted 26 June 2016; published 31 August 2016
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 threedimensional sonoelasticity theory of ships, analysis vibration and acoustic radiation characteristics of underwater double cylindrical shell by the use of Thaftsacoustic 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 threedimensional sonoelasticity 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 singlefrequency multisource 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, threedimensional sonoelasticity theory of ships is applied to analysis underwater vibroacoustic characteristics of double cylindrical shell, predict the underwater vibroacoustic characteristics with the help of software Thaftsacoustic 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 ${c}_{1}$ / m

0.004

Material density $\rho $ / kg·m^{3}

7.8e+03

Inner shell thickness ${c}_{2}$ / m

0.008

Elastic modulus $E$/ n·m^{2}

2.1e+11

Elastic panel thickness ${c}_{3}$ / m

0.008

Poisson’s ratio

0.3

End cover thickness ${c}_{4}$ / m

0.012

Fig. 1. Cylindrical shell model
a) Geometry model
b) FEM model
3. Threedimensional sonoelasticity 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
a) Coupling calculation model of outer flow field and shell
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 $\tau $ is the ideal sound medium outflow field, $\overline{S}$ is outer flow field corresponding average wet surface, on which the unit normal vector is $\overline{n}$; ${\tau}_{i}$ is imaginary inner flow field encircled by the wet surface ${S}_{i1}$ which is coincident with $\overline{S}$ and virtual impedance surface ${S}_{i2}$. Fig. 2(b) shows the inner flow field and shell coupling mathematical model, where ${\tau}_{b}$ is inner flow field, corresponding to the average wet surface ${\overline{S}}_{b1}$ and ${S}_{d1}$; ${\tau}_{c}$is imaginary outer flow field, corresponding to wet surface ${S}_{c}$ which is coincident with ${\overline{S}}_{b1}$. Since the virtual outer flow field is semiinfinite flow field and there is no resonance problem, so there is no need to introduce a virtual impedance surface; ${\tau}_{d}$ is imaginary inner flow field encircled by ${\overline{S}}_{b2}$, wet surface ${S}_{d1}$ and a virtual impedance surface ${S}_{d2}$.
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:
where ${D}_{r}$ is the modal displacement column vector corresponding to $r$order modal, $\left\{q\left(t\right)\right\}$ and ${q}_{r}\left(t\right)$ 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:
where [$a$], [$b$], [$c$] were dry structure modal generalized mass matrix, generalized damping matrix and generalized stiffness matrix; $\left\{{f}_{e}\right\}$, $\left\{{f}_{p}\right\}$ were generalized force column vector corresponding to outer excitation force in nonvoice 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:
$\overrightarrow{r}\in \tau \cup \overline{S}.$
The boundary integral equation for solving source intensity can be expressed as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\frac{1}{4\pi {\rho}_{0}}\underset{{S}_{i2}}{\iint}{Z}_{S}{u}_{ni2}\left({\overrightarrow{r}}_{0}\right)\frac{{\partial}^{2}G\left(\overrightarrow{r},{\overrightarrow{r}}_{0}\right)}{\partial {n}_{i}\left({\overrightarrow{r}}_{0}\right)\partial n\left(\overrightarrow{r}\right)}ds,\overrightarrow{r}\in \overline{S}.$
Boundary condition on the virtual impedance surface can be expressed as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\frac{1}{2\pi {\rho}_{0}}\underset{{S}_{i2}}{\iint}{Z}_{S}{u}_{ni2}\left({\overrightarrow{r}}_{0}\right)\frac{\partial G\left(\overrightarrow{r},{\overrightarrow{r}}_{0}\right)}{\partial {n}_{i}\left({\overrightarrow{r}}_{0}\right)}ds,\overrightarrow{r}\in {S}_{i2}.$
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 ${\sigma}_{r}$ on the wet surface and vibration displacement ${u}_{ni2r}$ on virtual impedance surface. Substitute these two quantities into the Eq. (3) can obtain velocity potential ${\varphi}_{r}\left(\overrightarrow{r}\right)$ 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:
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:
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:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\left\{{f}_{e}\left(\omega \right)\right\}+\left\{T\left(\omega \right)\right\},$
where $\left[{A}^{\left(1\right)}\right]$, $\left[{B}^{\left(1\right)}\right]$ and $\left[{C}^{\left(1\right)}\right]$ were respectively attached water mass matrix, attached water damping matrix and generalized recovery force coefficient matrix of dry modal contributed by outer flow field; $\left[{A}^{\left(2\right)}\right]$ and $\left[{B}^{\left(2\right)}\right]$ 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:
where ${\varphi}_{r}\left(x,y,z\right)$ can be obtained from Eq. (3). Combined with singlefrequency harmonic time factor ${e}^{i\omega 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 threedimensional sonoelasticity theory of ships in the prediction of underwater acoustic vibration characteristics of cylindrical shell, the result of Abaqus/Thaftsacoustic 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
a) Acceleration level
b) Error
Fig. 4. Vibration acceleration level and error on the inner wall of cylindrical shell
a) Acceleration level
b) Error
It can be seen from Fig. 3(a) that, in the 01000 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 01000 Hz, vibration acceleration level error of prediction and experimental test result are within 20 %, especially in 4001000 Hz band prediction error is less than 10 %, the prediction result is very close to experimental result.
In 01000 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 01000 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 0600 Hz was very close to experimental result; from Fig. 4(b), we can see in the entire frequency band 01000 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 01000 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
a) Sound pressure level
b) Error
Fig. 6. Sound pressure level and error in outer sound field of cylindrical shell
a) Sound pressure level
b) Error
It can be seen from Fig. 5(a) that, in the 01000 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 01000 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 01000 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 0180 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 1801000 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 0180 Hz frequency, the SPL prediction and experimental test result error of cylindrical shell in outer sound field are within 5 %, in the 1801000 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 01000 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 threedimensional sonoelasticity theory of ships, use software Thaftsacoustic 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 threedimensional sonoelasticity 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
 Junger M. C., Feit D. Sound Structures and Their Interaction. 2nd Ed. The MIT Press, Cambridge, 1986. [Search CrossRef]
 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. 15, (in Chinese). [Search CrossRef]
 Caresta M., Kessissoglou N. J. Structural and acoustic responses of a fluidloaded cylindrical hull with structural discontinuities. Applied Acoustics, Vol. 70, 2009, p. 954963. [Search CrossRef]
 Zhou Qizheng, Wang Deshi, Gao Shengyao Research on vibration and acoustic radiation of stiffened finite cylindrical shells in water under multiexcitations. Chinese Journal of Applied Mechanics, Vol. 30, Issue 5, 2013, p. 652657, (in Chinese). [Search CrossRef]
 Zhang Lei, Cao Yueyun, Yang Zichun Structural transfer path analysis for vibration and noise of a submerged doubleshell. Journal of Vibration and Shock, Vol. 31, Issue 20, 2012, p. 1216, (in Chinese). [Search CrossRef]
 Jin Guangwen, Zhang Linke, Miao Xuhong Vibration transmissibility of a submerged cylindrical doubleshell based on reconstructing velocity field. Journal of Vibration and Shock, Vol. 30, Issue 5, 2011, p. 218221, (in Chinese). [Search CrossRef]