1. Introduction
As human demand for marine resources increases, more and more ships and marine structures are also used to explore marine resources [
1]. This also leads to an increasing impact on marine life [
2,
3]. Noise is a special pollutant produced by ships, which will affect the living environment of marine organisms in a certain range [
4]. In addition, in the military field, there are higher requirements for the noise level of ships [
5]. Underwater vehicles are a special type of vessel. Their appearance and internal structure are relatively complex, so the noise components of this type of ship are also quite complex. In general, ship noise can be divided into three types: noise caused by mechanical equipment, noise caused by fluids, and noise caused by propellers. Among them, fluid-induced noise involves the interaction between the fluid and structure. This problem belongs to the interdisciplinary field of fluid mechanics, solid mechanics, and acoustics. Therefore, it has become a recent research hotspot. A large underwater vehicle can be approximately regarded as a cylindrical thin-walled structure. Local deformation of such structures occurs under the action of water pressure and free flow, and this deformation changes continuously during the underwater vehicle’s voyage, potentially reducing the vehicle’s service life if not effectively controlled. In order to reduce the effects of fluid-induced vibration on underwater vehicles, many scholars have used various methods to study the characteristics of fluid-induced vibration. Ying et al. [
6], taking the porous square cylinder model as the research object, studied the flow-induced vibration characteristics, influencing factors, and mechanism of the porous square cylinder via the numerical model method. Corcos et al. [
7] first measured the pulsating pressure on a cylindrical surface in a wind tunnel and analyzed the loading characteristics of the cylindrical surface in different frequency ranges. Based on the experimental results, Chase et al. [
8,
9,
10] proposed a turbulent pulsating pressure frequency wavenumber spectral model which approximately describes the pulsating pressure distribution on a flat plate. Liu et al. [
11] investigated the characteristics of pulsating pressure on the surface of subscale submarines using large vortex simulation (LES) and frequency wavenumber spectral methods and also explored the mechanism of flow-induced vibration noise. Wei et al. [
12] and Tian et al. [
13] studied the excitation load characteristics of a stern propeller under unsteady flow. The finite element and boundary element methods are also used to analyze the structural response of the submarine under the excitation load of the stern shaft. Li et al. [
14] observed the vortex structure on the surface of a ship propeller at different speeds through model tests. They also used computational fluid dynamics (CFD) methods to predict the pulsating pressure on the propeller surface, while evaluating the prediction accuracy of different frequency bands. Qin et al. [
15] studied the hydrodynamic characteristics of underwater vehicles based on LES and Ffows Williams and Hawkings (FW–H) methods and proposed a shape of underwater vehicles with low noise level.
With the continuous application of machine learning methods in engineering, some researchers have begun to apply neural network models to turbulence feature prediction. A flow-induced vibration power output prediction model was established by Li et al. [
16] through the mechanical learning method of a deep neural network (DNN). Based on this model, a power database was established, which significantly improved the prediction efficiency of the flow-induced vibration power response. Fukami et al. [
17,
18] reconstructed turbulent flow field data using a convolutional neural network model and improved the spatial resolution of turbulent flow fields. Rabault [
19] designed a reinforcement learning model for active flow control and reduced surface resistance of cylinders at moderate Reynolds numbers. Maulik [
20] applied machine learning methods to the classification problem of spatiotemporal dynamic turbulence models and verified the reliability of the model in predicting complex flow phenomena. Li et al. [
21] also used the neural network model to accurately detect the turbulent and non-turbulent interfaces in the three-dimensional space flow field. This method improves the precision of interface prediction. Raissi et al. [
22] combined nonlinear differential equations describing the laws of physics with neural networks. In solving physical problems, the convergence speed of the neural network model is accelerated, and its prediction accuracy is improved.
Under the influence of turbulent pulsating pressures, underwater structures are susceptible to flow-induced vibrations, which can significantly affect their service life. Chen et al. [
23] designed a flow-induced vibration test device for a flexible hydrofoil. Through this device, they studied the vibration characteristics of composite hydrofoils under turbulent pulsating pressure excitation and compared the accuracy of prediction via linear potential flow theory. Banafsheh [
24] and Ma [
25] et al. studied the influence of the wake flow of an upstream cylinder on vortex-induced vibration characteristics of a downstream cylinder in a tandem cylinder system through experimental tests. Xu [
26] and Ma [
27] et al. studied the vibration characteristics of a flexible cylinder under turbulent excitation, while Korkischko [
28] adjusted the boundary conditions of the cylinder to analyze the vortex-induced vibration characteristics of the cylinder under different boundary conditions. Jin et al. [
29,
30] developed an energy-based formula to describe the vibration characteristics of submarine hulls in heavy fluids. Based on this formula, they explored the coupling of stiffness with the vibration and acoustic radiation characteristics of underwater structures. Zhang et al. [
31] studied the unsteady flow characteristics of a propeller under stern wake conditions and analyzed the vibration characteristics of a propeller under unbalanced exciting force. There are two main methods for controlling typical flow-induced vibration in underwater structures: active control and passive control [
32]. A velocity feedback controller was developed by Baz [
33] to control vortex-induced vibration of a flexible cylinder, and the results were compared with theoretical predictions to verify the effectiveness of the controller. A method based on modal development techniques was proposed by Zhou et al. [
34], which can effectively predict the vibration–acoustic behavior of immersion reinforced composite plates under turbulent boundary layer excitation. A semi-analytical method is proposed by Chen [
35], which can predict the vibration and noise of annular ribbed cylindrical shells considering the internal structural coupling. A new semi-analytical method was proposed by Wang et al. [
36,
37] to predict the vibration characteristics of functionally graded composites. They also used this method to study the vibration characteristics of plate–shell structures under different boundary conditions.
From the literature review, it can be found that many scholars have studied the characteristics of turbulent pulsating pressure and flow-induced vibration of underwater structures, but there are few studies on the rapid prediction of pulsating pressure and flow-induced vibration. Machine learning methods have been applied in the field of flow field identification by some scholars, but this method has not been applied in the field of vibration prediction. The research shows that the neural network model in the machine learning method has good performance in solving the data prediction problem. Compared with other models, the BP neural network model has stronger approximation ability when solving nonlinear mapping relationships.
In this paper, the cone–cylinder–sphere combination structure is taken as the research object. An experimental model was designed to measure the pulsating pressure and vibration response of the structure surface at different flow velocities. The influence of sensor quality and water pressure on the calculation results of the model is analyzed. A neural network model is also established to predict the pulsating pressure and vibration response of the shell structure surface. The validity of the machine learning method is verified by comparing the predicted results with the experimental results.
2. Experiment Measurement
2.1. Experimental Model
The structure of the actual underwater vehicle is complicated, which is not suitable for direct theoretical and experimental research. Zou et al. [
38] simplified the underwater vehicle into a cone–cylinder–sphere combined structure when studying the flow-induced noise in underwater vehicles, which could better represent the main characteristics of underwater vehicles. The same simplification method was used in this study. The specific geometric parameters are shown in
Figure 1a, and the shape parameters of the combined structure can be represented by Equation (1).
where
L = 1.15 m represents the total length of the underwater vehicle, and
R = 0.1 m is the radius of the cylindrical section of the underwater vehicle.
According to the equation of forced vibration for objects, the vibration response of the structure under complex flow is affected not only by the fluctuating pressure of the fluid, but also by its inherent modes. In order to obtain the relationship between vibration response and excitation load in different regions of the underwater vehicle, two cone–cylinder–sphere combined structural models of the same size are designed and fabricated. In order to install the pressure sensor, several circular holes were cut into the surface of one of the models. The installation position of the pressure sensor is shown in
Figure 1b, and the parameters of the pressure sensor are consistent with those described in the literature [
38]. The surface of another model is smooth with marked locations for the installation of acceleration sensors on the internal wall. These sensors are used to measure the vibration acceleration at these locations. However, the acceleration sensor may change the mass distribution of the model itself and change the inherent vibration characteristics of the cone–pillar–ball composite structure. Considering sensitivity, mass, and size, the piezoelectric waterproof acceleration sensor was chosen. Each sensor has a mass of approximately 39 g. During the simulation, each sensor is loaded to its corresponding position on the surface of the cone–cylinder–sphere combined structure.
Figure 2 shows the influence of additional mass caused by the acceleration sensors on the natural vibration characteristics of the structure. As can be seen from the figure, the additional mass has almost no effect on the first three-order natural frequencies, but its effect gradually emerges as the frequency increases.
Figure 2c,d shows the structural modes at different frequencies. It can be seen that with the increase in frequency, the structural mode becomes more complex, and the local mode begins to dominate. This shows that the additional mass has a significant effect on the local mode. It also leads to a large error in the modal characteristics of the structure in the middle- and high-frequency range.
2.2. Experiment System
The experiment was conducted in a gravity-type water tunnel at the Key Laboratory of Underwater Acoustic Technology, Harbin Engineering University. The experimental system, detailed in the literature (Zou), met the testing requirements for this study. According to the installation requirements of the experimental system, the experimental model was processed and the installation interface was left in the upper part of the model. The basic principle of the experimental test is shown in
Figure 3. The test model was fixed to the test section through the top cover plate, and the acceleration sensor was installed to the inner wall of the model through magnetic suction. The sensor was connected to the signal acquisition system via cable. Two sets of acceleration sensors were also installed on the top cover plate of the test model to measure the background vibration of the experimental system. In order to study the fluid-induced vibration characteristics of the cone–column–sphere combined structures, the vibration responses of the structures at different flow velocities were measured.
The fluctuating pressure load and flow-induced vibration response of the cone–cylinder–sphere combined structure were separately measured. The layout of the pulsating pressure sensor requires mounting holes in the model surface. The vibration acceleration sensor is adsorbed on the inner surface of the structure via a magnetic base. On the same model, there is no way to install a vibration acceleration sensor and a pulsating pressure sensor in the same location. Two models with consistent structural parameters are designed to maintain the consistency of pulsating pressure and vibration response measurements. Model 1 was only used to measure the fluctuating pressure on the surface of the structure, while model 2 was only used to measure the vibration response on the inner wall of the structure. During the experiment, the model was lifted with a crane and waterproofed at the connection between the model and the test system. At least 3 measurements were taken for each flow velocity to reduce testing errors. After completing all the tests for all operating conditions with model 1, it was replaced with model 2, and the above process was repeated for measurement.
The uncertainty of the test results under the same working condition was compared during the test. For the measurement of pulsating pressure, it is considered that the measurement results with a dispersion of less than 3 dB are valid. In all the repeated tests of each working condition, three groups of test data that meet the requirements are selected as the basis for analysis.
The measurement points of the model are numbered 1# to 15# from the nose to the tail.
Figure 4 shows the measured results of the acceleration of the model and the top cover plate at some flow velocities. The experimental results show that there is a significant difference in vibration acceleration between the model surface and the top cover plate in the frequency range of 10 Hz–1000 Hz. This finding provides evidence that the vibration of the top cover plate has no noticeable effect on the model surface. The results show that it is reasonable to truncate the top cover plate in the subsequent analysis.
The mechanism of cavitation noise is that when the fluid moves at high speed, the local pressure decreases, and cavitation bubbles are formed. When the bubble bursts, it produces a strong pulsing noise. The premise of cavitation noise is that the flow is cavitation.
The cavitation number (
) can be used to measure whether cavitation occurs in the fluid surrounding a structure and the degree of development of cavitation. The cavitation number can be expressed as follows:
In the formula, represents the absolute pressure of the reference point, represents the free flow velocity, represents the liquid density, and represents the saturated vapor pressure of the liquid at the corresponding temperature.
In this study, the maximum flow rate is 9.81 m/s, and the flow temperature is less than 30 °C. According to the above formula, the cavitation number is about 2.2, and it is difficult to have a stable cavitation phenomenon.
3. Theoretical Methods
3.1. Turbulence Model
The large eddy simulation (LES) method is a technique that lies between direct numerical simulation (DNS) and Reynolds-averaged Navier–Stokes (RANS) methods. It provides high accuracy and computational efficiency in simulating unsteady flow. The filtered continuity equation and the incompressible Navier–Stokes equation can be expressed as follows:
where
represents the velocity component associated with
,
represents the filtered pressure of the fluid,
and
denote the averaged velocity components after filtering,
represents the density of the fluid, and
represents the dynamic viscosity coefficient of the fluid.
represents sub-grid scale Reynolds stress. This physical quantity reflects the interaction between large-scale eddies and small-scale eddies, such as energy and momentum exchange, feedback of small-scale eddies on large-scale eddies, and others. The Smagorinsky sub-grid scale model for turbulence was originally proposed in LES and remains one of the most commonly used sub-grid scale models. In this study, the adaptive wall local vortex viscosity (WALE) sub-grid scale model is selected. According to Boussinesq’s assumption, the sub-grid scale stress tensor can be expressed as follows:
where
represents the Kronecker delta function, and
represents the sub-grid scale Reynolds stress tensor.
represents the sub-grid scale eddy viscosity coefficient, and
represents the filter width.
represents the model coefficient.
represents deformation coefficient.
represents the traceless, symmetric part of the squared filtered velocity gradient tensor.
represents the identity tensor.
The WALE sub-grid scale model uses a new form of velocity gradient tensor that does not require any form of near-wall damping. It can provide accurate proportional scaling automatically at the wall. Posa and Balaras [
39] used the WALE model coupled with the standard equilibrium wall-layer model for the simulation of the DARPA Suboff at very high Reynolds number, obtaining a good agreement with the experimental data.
3.2. Modal Superposition Method
The vibration equation of a multi-degree-of-freedom elastic body under external forces can be expressed as follows:
where
represents the mass matrix of the system,
represents the damping matrix of the system,
represents the stiffness matrix of the system,
represents the generalized displacement vector of the system, and
represents the excitation load of the system.
When the system is under proportional damping, the damping matrix can be expressed as a linear combination of
and
, and
can be expressed as follows:
In the formula, and represent rayleigh damping coefficients.
Define the modal matrix of a system as
. Therefore, Equation (9) can be transformed as follows:
Therefore, Equation (10) can be further expressed as follows:
where
represents the
r-th order modal frequency of the system. When each component of the excitation load
is a determined function, each equation in Formula (15) can be decomposed into independent single-degree-of-freedom systems.
The vibration response of a structure subjected to turbulence-induced pressure fluctuations can be solved using the modal superposition method. This method assumes that the surface displacement of the structure can be represented by the linear combination of modal shapes and uses modal frequencies and structural shapes to decouple the fluid–structure coupling vibration equation. A modal coordinate system is established to calculate the response in the modal coordinate system, and the structural response in the original physical coordinate system is obtained via combination. This method can be expressed by Equation (17).
where
represents the displacement vector corresponding to each node degree of freedom,
represents the
i-th order modal shape vector, and
represents the
i-th order modal coordinate. The relationship between the cross power spectral density function (
) in generalized coordinates and the cross power spectral density function (
) of the excitation load can be expressed as follows:
where
and
represent the system transformation matrix.
By substituting Equation (18) into Equation (17), the structural vibration response can be further expressed as follows:
In the equation, represents power spectral density, and is defined as the autocorrelation function of power spectral density. By using this definition, the power spectral density function of the structural response can be calculated.
3.3. Machine Learning Framework
The internal relationships in complex dynamic problems can be effectively described with neural networks. In this study, a neural network framework is proposed. As shown in
Figure 5, the framework consists of a feature input, activation function layer, batch normalization, hidden layer, and feature output.
The feature input layer is the data collector of the neural network and performs normalization of the global data. In this study, when the surface pulsation pressure of the cone–cylinder–sphere combined structure is predicted, the input data , where f represents the frequency, u represents the distance between the measuring point and the head of the model, l represents the total length of the model, r represents the radius of the cross section where the measuring point is located, R represents the maximum radius of the model cross-section, represents the pressure at the head of the cylinder, and represents the pressure at the end of the cylinder. The output variable is the pulsating pressure at the measuring point. When the surface fluid-induced vibration response of the cone–cylindric–spherical composite structure is predicted, the input data , where represents the participation factor of the structural mode, and represents the pulsating pressure at the measuring point. The output variable is the fluid-induced vibration response at the measuring point.
The activation function layer applies the selected activation function to the feature input layer to obtain a high-dimensional representation. The activation function plays a very important role in the neural network framework, as it greatly affects the gradient features during the training process. In this study, the tanh-sigmoid function and linear unit function were chosen as the activation functions, which can be expressed as follows:
To enhance the robustness of the neural network, batch normalization was incorporated into its framework to mitigate the problem of gradient explosion. The input feature data set was divided into a training set and a testing set, which were, respectively, utilized for neural network model training and effectiveness verification. The hidden layer, which represents the core computation module of the neural network architecture, processes the input feature parameters with Equation (22) to obtain the input of the
h-th hidden layer neuron. Its output is then processed by Equation (23) to generate the input of the
j-th output neuron.
In the formula, represents the input received by the h-th neuron in the hidden layer, and represents the input received by the j-th neuron in the output layer.
If the neural network’s output for the training set
is
, then the mean squared error of the neural network on the training set can be expressed as follows:
In order to make the training result closer to the real result, the error is backpropagated, and the weight and threshold are adjusted according to the error of the hidden layer neurons. After repeated iterations, the training stops when the error meets the set requirements. The cumulative error of the entire training set can be expressed as Formula (25). To avoid overfitting of the neural network, which may result in training low set error but test set high error, a component describing the complexity of the network was added to the error function. Then Formula (25) can be represented as Formula (26).