1. Introduction
In the past two decades, as a new type of passive vibration reduction and noise attenuation technology, the acoustic black hole (ABH) was studied a lot in terms of theoretical development, modeling and simulation, and analysis of vibration and acoustic characteristics, and it is expected to provide a new solution to the problems of vibration and noise in engineering [
1,
2,
3,
4,
5]. Mironov found in a thin wedge-shaped structure that if the thickness of the structure is reduced in a certain power exponential form (the power exponent is not less than 2), the velocity of the flexural wave in the structure will gradually decrease as the thickness of the structure decreases [
6]. Ideally, the velocity of the flexural wave can be reduced to zero to achieve zero reflection. Krylov first proposed the concept of an “acoustic black hole” to describe the above phenomenon and applied the ABH to beam and plate structures. Krylov used the geometric acoustic method [
7,
8] and experimental method [
9] to study the characteristics of ABH structure, such as flexural wave propagation, vibration suppression, and sound radiation reduction. The ABH structure can achieve the concentration and manipulation effects of flexural waves in solid structures in a wide frequency range, thereby achieving vibration reduction and noise attenuation; in addition, it also contributes to the light-weighting of the structures [
10].
There is much research on the vibration characteristics of the ABH structures. Tang et al. established a semi-analytical model by using Mexican wavelet function to fit the sharply changing wave group at the edge of the ABH to analyze the vibration characteristics of one-dimensional ABH structure and studied the influences of damping layer thickness and damping layer stiffness on the dynamic response of the ABH structure [
11,
12]. Deng et al. established a dynamic model of the truncated one-dimensional ABH beam containing passive constrained viscoelastic layers (PCVL) using the semi-analytical method, and both the analysis and the experimental results show that the truncated one-dimensional ABH beam with PCVLs can achieve better vibration attenuation effects than the one covered with a viscoelastic damping layer [
13]. Deng et al. also analyzed the influence of ABH parameters on the vibration characteristics of the two-dimensional ABH circular plate based on the semi-analytical method [
14]. Additionally, they studied the vibration response of the new ring-shaped two-dimensional ABH, and the results showed that the ring-shaped ABH exhibited excellent vibration isolation performance over the entire frequency range [
15]. Ma et al. established a semi-analytical model to analyze the vibration of a two-dimensional ABH plate and used the Daubechies wavelet function to express the lateral displacement of the plate in enhancing the effectiveness of the Rayleigh–Ritz method [
16]. Based on the approximate analytical method, O’boy et al. analyzed the frequency response characteristics of a square plate containing a two-dimensional square ABH [
17] and a circular plate containing a two-dimensional circular ABH [
18] and investigated the influence of the damping layer. The result shows that the cross mobility reaches 17 dB attenuation at the resonance peaks. Conlon arranged multiple two-dimensional ABHs in a thin plate with a certain periodicity and analyzed the vibration response of the plate by establishing a finite element model [
19]. Zhao et al. evaluated the vibration performance of the thin-walled ABH structure under transient and steady-state excitation through numerical simulation and experimental methods [
20], and both the numerical and experimental results show that the proposed structure can achieve an efficient and broadband vibration attenuation performance. Gao et al. proposed a V-shaped folded ABH beam, and the numerical simulation and experimental results show that the V-shaped folded ABH beam has an ultra-wide and complete band gap below 1 kHz, which is expected to be used for low-frequency vibration reduction and wave attenuation [
21]. Feurtado et al. experimentally measured the average mobility and radiated sound power of the surface of the aluminum plate embedded with 20 two-dimensional ABH arrays attached with damping layers [
22]. The results showed that compared with the reference plate of uniform thickness, the average surface mobility of the fully damped ABH plate is reduced by up to 14 dB. Bowyer et al. experimentally studied the influence of the power-law contour grooves in the plate on the damped flexural vibration of the one-dimensional ABH structure [
23] and carried out experiments to study the vibration response of the two-dimensional ABH structures, mainly including rectangular and circular ABH plates which are made of steel, aluminum or composite materials and contain circular indentations or dimples wedge [
24,
25]. Besides, traditional ABH structures cannot be changed to different frequencies, and Ning et al. found that the active control can be applied to black hole characteristics of the flexural wave, which can provide a much wider frequency region [
26,
27].
The acoustic black hole can concentrate flexural waves in the structure at or near the ABH region, and with the help of the damping layer or other additional components to dissipate flexural wave energy, it can achieve a good effect of suppressing structural vibration with a very small amount of damping layer. However, most of the current research on ABH still stays in the basic theory, and expanding the application research of ABH structures in different engineering fields to achieve better vibration reduction effects in practice is one of the main issues needed to be studied for the ABH [
28]. The compressor is the core component of the refrigerator, and when the compressor works, the internal motor drives the piston to reciprocate. At this time, the working torque of the piston, the unbalanced moment of inertia of the structure, the friction between the moving pairs, the electromagnetic force inside the motor, etc., can all cause vibration [
29], and the generated vibration is transmitted to the compressor shell and further to the support plate and the body of the refrigerator. The vibration of the compressor is the main source of vibration for all types of refrigerators. Since the refrigerator compressor is generally a fully enclosed structure, it is often difficult to optimize internal structures to reduce compressor vibration. Therefore, blocking the transmission of compressor vibration to the refrigerator body as much as possible is an important method to solve such problems. The support plate is an important component connecting the compressor and the refrigerator body, and almost all vibrations generated by the compressor are transmitted to the refrigerator body through the support plate. Utilizing the concentration effect of the ABH on flexural waves, it is expected to concentrate the vibration energy of the compressor in the ABH region by embedding ABHs on the support plate. At this time, the gathered vibration energy cannot be transferred to the refrigerator body through the support plate. Therefore, it is feasible to apply the ABH to suppress the propagation of refrigerator compressor vibration.
In this paper, the vibration and noise measurement and analysis of key components such as the compressor and the support plate were carried out, and then the damping ABHs were embedded in the support plate, and the effect of the damping ABHs in suppressing the vibration transmission of the compressor was studied through numerical simulation. It is expected to provide a certain reference and guidance for the application of ABHs in the field of home appliances for vibration reduction. To our best knowledge, this paper is the first to study the vibration suppression effect of ABH on refrigerators.
2. Measurement of Acoustic and Vibration Characteristics of Key Components of Frost-Free Refrigerator
Figure 1 shows the physical image of the back of the frost-free refrigerator. The dimensions of the refrigerator are: 1722 mm (height) × 910 mm (width) × 610 mm (length). The compressor, condenser, and condenser fan are all fixed on the support plate at the bottom, and the support plate and the refrigerator body are connected by bolts. The refrigerant absorbs the heat in the refrigerator at the evaporator, and then the temperature of the refrigerant rises and flows back to the compressor [
30]. The compressor provides power to transport the refrigerant to the condenser for heat dissipation and cooling and then flows to the evaporator again for heat exchange and reciprocating circulation.
To investigate the vibration of the refrigerator compressor and its influence on the vibration and radiated sound of the refrigerator body, the vibration and acoustic measurement experiment of the compressor, the support plate, and the refrigerator body was carried out. The physical image of the experimental platform is shown in
Figure 2. The measured vibration and acoustic signals can be used to analyze the transmission of the compressor vibration and the influence of the compressor vibration on the radiated sound of the refrigerator body, and they will also be used as the excitation input for subsequent numerical simulations. Six unidirectional acceleration sensors are placed on the surface of the compressor, the four positions on the support plate connected to the compressor, and the center position on the back of the refrigerator, and an acoustic sensor is placed 0.8 m behind the back of the refrigerator (the vertical distance between the acoustic sensor and the ground is also 0.8 m). The sensor types and the location of the measuring points are listed in
Table 1. The m + p measurement system is used to measure the vibration of the refrigerator compressor, support plate, and the sound and vibration of the back wall of the refrigerator. Both the acceleration sensors and the acoustic sensor are connected to the m + p VibPilot. The m + p VibPilot can collect the time-domain signals of each sensor and transmit them to the laptop computer, and then use the SO Analyser software to process and further analyze the time-domain signals.
Figure 3,
Figure 4 and
Figure 5 show the vibration acceleration spectrum of the compressor surface, the vibration acceleration spectrum of the outer wall on the back of the refrigerator, and the 1/3-octave frequency spectrum of the radiated sound from the outer wall on the back of the refrigerator. The acoustic measurement experiment was carried out in a semi-anechoic chamber from Landtop Faist, with a background noise of 12 dB (A). From these spectrograms, the main vibration frequencies and corresponding amplitudes of the compressor can be obtained, and the transmission of compressor vibration to the refrigerator body can also be studied. The 1/3-octave frequency spectrum of the radiated sound from the outer wall of the refrigerator is also used to analyze the coherence between the vibration of the compressor and the sound. The vibration acceleration frequency spectra of the four positions on the support plate connected to the compressor can be used as an excitation input in numerical modeling to study the effect of suppressing the vibration of the support plate before and after the ABHs embedded in the support plate.
3. Coherence Analysis of Vibration and Noise
The coherence theory can be used to analyze the radiated sound of the back wall of the refrigerator and the vibration acceleration signals of other measuring points to determine the main vibration source of the radiated noise. For the measurement in this paper, the radiated sound from the outer wall on the back of the refrigerator belongs to a multiple-input single-output system, which can be analyzed using a coherent power spectrum. The coherence function of the vibration acceleration signal
X(
t) and the sound pressure signal
Y(
t) is defined as [
31]:
where
and
are the self-power spectra of the vibration acceleration signal and the sound pressure signal, respectively,
is the cross power spectrum of the vibration acceleration signal and the sound pressure signal.
The calculation equation of self-power spectral density
is:
where
is the self-correlation function
.
The calculation equation of cross power spectral density
is:
where
is the cross-correlation function
.
Since there is always some noise interference in the actual signals, the value of the coherence function between the measured signals is always less than 1, and the closer the value of the coherence function is to 1, the greater the degree of coherence between the sound radiated from the outer wall of the refrigerator and the vibration acceleration signal at that point; on the contrary, the closer the value of the coherence function is to 0, the smaller the degree of coherence between the sound radiated from the outer wall of the refrigerator and the vibration acceleration signal.
Figure 6 shows the coherence function curve between the radiated sound of the back wall of the refrigerator and the vibration acceleration of the compressor surface. It can be seen from
Figure 6 that the frequencies of the greater coherence between the radiated sound and the vibration signal are mainly concentrated at 57.813 Hz, 68.750 Hz, 117.19 Hz, 637.50 Hz, 695.31 Hz, 1390.6 Hz, 2087.5 Hz, and 2260.9 Hz. At these frequencies, the coherence function values between the radiated sound of the back wall of the refrigerator and the vibration acceleration of the compressor surface are greater than 0.8, which can be considered to have high coherence, that is, the structurally radiated sound at these frequencies is mainly caused by the vibration at the measuring point 5.
Figure 7 shows the coherence function curve between the radiated sound of the back wall of the refrigerator and the vibration acceleration on the back wall of the refrigerator. It can be seen from
Figure 7 that the frequencies of the greater coherence between the radiated sound and the vibration signal are mainly concentrated at 57.813 Hz, 68.750 Hz, 100 Hz, 115.63 Hz, 251.56 Hz, 1390.6 Hz, and 2087.5 Hz. At these frequencies, the coherence function values between the radiated sound of the back wall of the refrigerator and the vibration signal of the back wall of the refrigerator are also greater than 0.8, which can be considered to have high coherence, that is, the structurally radiated sound at these frequencies is mainly caused by the vibration at the measuring point 6.
By combining the 1/3-octave frequency spectrum of the radiated sound from the back wall of the refrigerator in
Figure 5, it can be concluded that the radiated sound from the back wall of the refrigerator at 63 Hz, 125 Hz, 630 Hz, 1250 Hz, 2000 Hz, and 2500 Hz is mainly caused by the vibration of the compressor, and the sound radiated from the back wall of the refrigerator at other center frequencies is less coherent with the vibration of the compressor. Therefore, after the above analysis, it can be concluded that the vibration of the compressor is transmitted to the refrigerator body through the support plate and causes the back wall of the refrigerator to radiate sound. Suppressing the transmission of the support plate to the vibration of the refrigerator compressor can significantly reduce the vibration of the refrigerator body and the radiation sound of the back wall of the refrigerator.
4. Acoustic Black Hole Theory and Verification of the ABH Effect
A two-dimensional (2-D) ABH was embedded on the rectangular elastic plate; the schematic diagram of the structure is shown in
Figure 8. The thickness of the ABH plate can be expressed by the following equation:
where
ε is a constant, and
m is the power exponent of the thickness change in the ABH region,
,
b is the distance between the center of the ABH region and the coordinate origin in the
z-axis direction, and the thickness of the uniform region of the ABH plate is
hp. However, due to the limitation of the processing technology, the ideal ABH structure cannot be produced in the actual manufacturing process, so there must be a truncation existing in the actual ABH structure. For the model in
Figure 8, the thickness of the truncated region can be expressed as:
The existence of the truncation will significantly increase the reflection coefficient of the flexural wave. For an isotropic thin plate, the propagation speed of the flexural wave in the plate is:
where
B is the bending stiffness of the ABH plate material,
,
E is the Young’s modulus of the ABH plate material,
ν is the Poisson’s ratio;
ω is the angular frequency,
m is the areal density of the ABH plate,
, and
is the density of the ABH plate material. When the flexural wave propagates from the uniform region of the ABH plate to the ABH region, as the thickness of the plate gradually decreases, the velocity of the flexural wave also gradually decreases [
32]. At this time, the wavelength is compressed, and the amplitude of the flexural wave gradually increases so as to achieve the concentration effect of the flexural wave.
In order to verify the concentration effect of the ABH plate on flexural waves, a square plate embedded with a 2-D ABH was designed. The geometric size of the ABH plate is 200 × 200 mm, the thickness of the uniform region is 4 mm, and the coordinate origin is established at the geometric center of the middle plane in the thickness direction of the ABH plate. The thickness equation of the ABH region is . According to Equation (5), the thickness of the truncation of the ABH region is 0.5 mm, and the outer diameter of the ABH region is about 48.3 mm. The multi-physics simulation software COMSOL is used to carry out FEM modeling and numerical calculation of flexural wave propagation in the ABH plate. To ensure the calculation accuracy, the extremely fine level is used to mesh the 2-D ABH plate, and the ABH regions are encrypted to ensure that the number of grids in each wavelength is more than six.
The FEM model of the ABH plate is shown in
Figure 9. The material of the 2-D ABH plate is 1060 aluminum, the Young’s modulus is 69 GPa, the Poisson’s ratio is 0.33, and the density is 2710 kg/m
3. The left side face of the ABH plate is clamped, the remaining faces and edges are set as free boundaries, and the excitation point is located at the midpoint of the right side of the upper surface of the ABH plate. A simple harmonic excitation force perpendicular to the ABH plate is applied to the excitation point, the amplitude of the excitation force is 100 N, the frequency is 80 MHz, and the action time is 0~2 × 10
−5 s. The transient calculation of the flexural wave propagation in the ABH plate was carried out. The time length of the transient calculation is 0~5 × 10
−4 s, and the time step is set to 1.5 × 10
−7 s.
Figure 10 shows the stress distribution of the ABH plate at different times to characterize the propagation process of flexural waves in the ABH plate. It can be seen from
Figure 10 that at
t = 7.5 × 10
−5 s, the flexural waves, excited by the force, begin to propagate in the plate, and there are obvious wave crests and troughs. At
t = 1.2 × 10
−4 s, the flexural waves propagate to the ABH region, the wavelength of the flexural wave is compressed, and an obvious concentration effect is produced. When
t = 1.95 × 10
−4 s, the flexural waves reached the upper and lower boundaries of the ABH plate and were reflected. At that time, the ABH region still had a significant aggregation effect on the flexural waves so that the flexural waves near the ABH region did not significantly continue to propagate to the left. When
t = 2.31 × 10
−4 s, the flexural waves transmitted to the two corners on the left side of the ABH plate, but due to the concentration effect of the ABH region on the flexural wave, the flexural waves in the area on the left side of the ABH region did not fluctuate strongly, and most of the flexural waves were still manipulated by the ABH region. The concentration effect of the 2-D ABH plate on flexural waves is clearly and intuitively verified by the above numerical calculation results.
5. Modeling of the Support Plate Embedded with Acoustic Black Holes
From the analysis in
Section 3, it can be seen that the support plate is a key component to suppress the vibration of the refrigerator body and the radiation sound of the back wall of the refrigerator. To facilitate the subsequent meshing, the three-dimensional model of the support plate was simplified, mainly by removing chamfers, smaller openings, and more complicated small curved faces. These small geometric features have no obvious influence on the overall modes and dynamic response of the support plate, and the simplified geometric model of the support plate is shown in
Figure 11. The length of the support plate is 901 mm, the width is 203.5 mm, and the thickness is 1.8 mm. The four circular holes on the support plate are used to fix the compressor. The second-order tetrahedral unit C3D10 is used to mesh the support plate, and the grid size is 4 mm. The mesh in the small curved area of the support plate is encrypted to ensure a smooth transition from small cells to large cells and the minimum cell size is 0.8 mm during encryption.
The 2 × 4 two-dimensional ABHs are embedded on the support plate, and the geometric model of the ABH support plate is shown in
Figure 12. The geometric parameters of the two-dimensional ABH structure are:
,
, and the thickness of the truncation is 0.3 mm. It can be calculated that the outer diameter of the 2-D ABHs is 34.641 mm according to Equations (4) and (5). The centers of the two rows of 2-D ABHs are symmetric about the centerline of the support plate, and the distance is 120 mm.
The eight two-dimensional ABH regions are completely attached with damping layers with a thickness of 1.5 mm, and the parameters of the damping material are shown in
Table 2. Dynamic analysis, which can use structural damping assumptions, includes steady-state response analysis and random response analysis. The structural damping entered in the material definition is actually the loss factor of the material, and the relationship between loss factor and damping ratio is
ξ = η/2 for the Kelvin–Voigt model. When the loss factor of the viscoelastic material is obtained, the loss factor of the component made of this viscoelastic material can be evaluated. Then the damping ratio of the component and the damping coefficient of the component is obtained, which can be further used for structural analysis. The second-order tetrahedral unit C3D10 is used to mesh the damping ABH support plate. The grid size of the ABH regions and damping layers is 1 mm, and the mesh in the small curved area of the ABH regions and damping layers is encrypted to ensure a smooth transition from small cells to large cells, and the minimum cell size is 0.2 mm during encryption. The grid size of the other regions of the damping ABH support plate is the same as that of the support plate without ABHs. The entire FEM model of the damping ABH support plate has a total of 821,000 grids and 236,000 nodes. Since the grids in the ABH regions and damping layers are very dense, for ease of display, a partial view of the upper left part of the FEM model of the damping ABH support plate is shown in
Figure 13.
6. Influence of ABHs on the Vibration Characteristics of the Support Plate
Based on the FEM models established in
Section 5, the modal superposition method is used to carry out the modal analysis and the vibration response calculation of the original support plate of the refrigerator and the damping ABH support plate. The left and right sides of these two support plates are clamped, and the remaining surfaces are free boundaries. The modal calculation frequency is 0~5 kHz, and
Figure 14 shows the typical modes of the original support plate and the damping ABH support plate in the range of 0~5 kHz. It can be seen from
Figure 14 that in the frequency range below 1 kHz, there is no obvious difference between the mode shapes of the damping ABH support plate and the mode shapes of the original support plate; only the displacement in the ABH regions is increased in the 13th-order mode. However, for the frequency range above 1000 Hz, the modal distribution of the support plate was changed significantly after the damping ABHs were embedded. At the 21st, 43rd, 69th order (corresponding to the 68th order of the original support plate), due to the combined action of the 2-D ABHs and the damping layers, the modal displacement outside the ABH regions of the support plate was significantly reduced, and the vibration of the support plate was mainly concentrated in the ABH regions at these frequencies.
Taking the vibration acceleration signals of the measured points 1, 2, 3, and 4 in
Section 2 as the excitation input and applying the above excitation to the grid nodes corresponding to these four measuring points in the FEM models of the support plate, the vibration response of the original support plate and the damping ABH support plate can be obtained through numerical calculation in the commercial software LMS Virtual Lab.
Figure 15 shows the vibration acceleration cloud diagram of the original support plate and the damping ABH support plate. It can be seen from
Figure 15 that in the frequency band below 1000 Hz, embedding damping ABHs has no significant effect on the vibration acceleration amplitude or distribution of the support plate. In the range of 2 kHz to 4 kHz, most of the vibration energy is concentrated by the ABH regions of the damping ABH support plate, and the vibration acceleration of the rest of the damping ABH support plate is significantly reduced. In addition, the peak values of the vibration acceleration of the entire damping ABH support plate is much smaller than the peak values of the original support plate at the same frequencies, indicating that the 2-D damping ABHs have an obvious effect on the vibration suppression of the refrigerator support plate at these frequencies.
To better analyze the influence of the damping ABHs on the vibration acceleration of the support plate in the whole frequency band, three monitoring points were established at the left, middle and right positions of the support plate, as shown in
Figure 16.
Figure 17,
Figure 18 and
Figure 19 are the vibration acceleration response spectrograms of the three monitoring points. It can be seen from
Figure 17 that for the left monitoring point, the 2-D damping ABHs can significantly reduce the vibration acceleration of the refrigerator support plate in the range of 3 kHz to 5 kHz, and the maximum attenuation value can reach 6.69 dB (with a reference of 0.1 m/s
2); in the 0~3 kHz frequency band, the acceleration peak values of the left monitoring point are reduced at some frequencies but increased at other frequencies. It can be concluded from
Figure 18 that for the middle monitoring point, the 2-D damping ABHs have no obvious suppression effect on the acceleration peak values of the support plate; however, the acceleration peak values are increased at some frequencies. It can be seen from
Figure 19 that for the right monitoring point, the 2-D damping ABHs show a good vibration suppression effect in almost the whole frequency range, and the maximum acceleration attenuation can reach more than 10 dB, indicating that the 2-D damping ABHs have an obvious vibration suppression effect on the right side of the support plate.
Considering that the left and right sides of the support plate are connected to the refrigerator body, and the mass and volume of the refrigerator body are much larger than the mass and volume of the support plate, it can be considered that the true boundary conditions of the support plate are basically the same as in the simulation. Combined with the sound and vibration measurement in
Section 2, it can be concluded that the 2-D damping ABHs can suppress the vibrations generated by the compressor at 1390.6 Hz, 2087.5 Hz, and 2260.9 Hz from being transmitted to the refrigerator body through the support plate, thereby suppressing the radiated sound of the back wall of the refrigerator at frequencies of 1250 Hz, 2000 Hz, and 2500 Hz.